<figure><img src="image_1.png" /></figure>

Article

# Distinguishing the Effects of Water Volumes versus Stocking Densities on the Skeletal Quality during the Pre-Ongrowing Phase of Gilthead Seabream (*Sparus aurata*)

Zachary Dellacqua<sup>1,2,*</sup>, Claudia Di Biagio<sup>1,3 ID</sup>, Corrado Costa<sup>4 ID</sup>, Pedro Pousão-Ferreira<sup>5 ID</sup>, Laura Ribeiro<sup>5 ID</sup>, Marisa Barata<sup>5</sup>, Paulo J. Gavaia<sup>6 ID</sup>, Francesco Mattei<sup>7</sup>, Andrea Fabris<sup>8</sup>, Marisol Izquierdo<sup>2 ID</sup> and Clara Boglione<sup>1 ID</sup>

1. Department of Biology, University of Rome 'Tor Vergata', 00133 Rome, Italy

2. Ecoaqua Institute, University of Las Palmas de Gran Canaria, 35214 Telde, Gran Canaria, Spain

3. Laboratory of Evolutionary Developmental Biology, University of Ghent, 9000 Ghent, Belgium

4. CREA—Consiglio per la Ricerca in Agricoltura e L'analisi Dell'economia Agraria (CREA)—Centro di Ricerca Ingegneria e Trasformazioni Agroalimentari, 00015 Rome, Italy

5. IPMA—Instituto Portugues do Mar e Atmosfera—Research Station, 8700-305 Olhão, Portugal

6. CCMAR—Centre of Marine Sciences, University of the Algarve, 8005-139 Faro, Portugal

7. UMR 7093, Laboratoire d'Oceanographie de Villefranche (LOV), Sorbonne University, 06230 Villefranche-sur-Mer, France

8. Associazione Piscicoltori Italiani, 37135 Verona, Italy

* Correspondence: dellacquaz95@gmail.com; Tel.: +39-351-857-0196

**Simple Summary:** The development of skeletal anomalies in the early life stage of gilthead seabream (*Sparus aurata*) poses a significant challenge for farmers, affecting their profit margins, animal welfare, and the consumers' perception of the aquaculture industry. Although many factors have been considered to be causative in the development of skeletal anomalies in marine finfish, the stocking density and available swimming space represent two key parameters which can be easily manipulated by the farmers during the critical phase of pre-ongrowing (prior to being placed in sea cages). This work aims at distinguishing which among tank volume and stocking density is the driving factor eliciting the development of skeletal anomalies during the pre-ongrowing phase in gilthead seabream, a productive cycle in which many skeletal anomalies can arise, particularly those affecting the vertebral axis. The results from this work indicated that particular cranial and axis deformities affected fish in higher incidences when they were reared at higher densities. The results are discussed through an eco-evo-devo approach in relation to the potential mechanisms at play affecting the increased prevalence of skeletal anomalies found. This research represents an intriguing contribution to aquaculture with results that can be applied directly to the production methods used by fish farmers to ameliorate the skeletal and morphological quality of farmed gilthead seabream.

**Abstract:** Gilthead seabream (*Sparus aurata*) production is a highly valued aquaculture industry in Europe. The presence of skeletal deformities in farmed gilthead seabream represents a major bottleneck for the industry leading to economic losses, negative impacts on the consumers' perception of aquaculture, and animal welfare issues for the fish. Although past work has primarily focused on the hatchery phase to reduce the incidence of skeletal anomalies, this work targets the successive pre-ongrowing phase in which more severe anomalies affecting the external shape often arise. This work aimed to test the effects of: (i) larger and smaller tank volumes, stocked at the same density; and (ii) higher and lower stocking densities maintained in the same water volume, on the skeleton of gilthead seabream fingerlings reared for ~63 days at a pilot scale. Experimental rearing was conducted with gilthead seabream juveniles (~6.7 ± 2.5 g), which were selected as 'non-deformed' based on external inspection, stocked at three different densities (Low Density (LD): 5 kg/m³; Medium Density (MD): 10 kg/m³; High Density (HD): 20 kg/m³) in both 500 L and 1000 L tanks. Gilthead seabream were sampled for growth performance and radiographed to assess the skeletal elements at the beginning and end of the experimental trial. Results revealed that (i) LD fish were significantly longer than HD fish, although there were no differences in final weights, regardless of the water volume; (ii) an

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Dellacqua, Z.; Di Biagio, C.; Costa, C.; Pousão-Ferreira, P.; Ribeiro, L.; Barata, M.; Gavaia, P.J.; Mattei, F.; Fabris, A.; Izquierdo, M.; et al.

Distinguishing the Effects of Water Volumes versus Stocking Densities on the Skeletal Quality during the Pre-Ongrowing Phase of Gilthead Seabream (*Sparus aurata*). *Animals* **2023**, *13*, 557. https://doi.org/10.3390/ani13040557

Academic Editor: Constanze Pietsch

Received: 11 January 2023

Revised: 31 January 2023

Accepted: 1 February 2023

Published: 5 February 2023

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2023 by the authors.
Licensee MDPI, Basel, Switzerland.
This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).
increase in the prevalence of seabream exhibiting cranial and vertebral axis anomalies was found to be associated with increased density. These results suggest that farmers can significantly reduce the presence of some cranial and axis anomalies affecting pre-ongrown gilthead seabream by reducing the stocking density.

**Keywords:** morphometric quality; *Sparus aurata*; skeletal anomalies; stocking density; swimming space; tank volume

# 1. Introduction

Gilthead seabream (*Sparus aurata*) production is one of the main aquaculture industries in the Mediterranean, producing 258,754 tons of fish in 2019 [1]. Modern ongrowing farms (farms in which reared fish are fattened up to reach the commercial size of 200–350 g) often employ intensive strategies with high energy demands, tailored feeds, and high stocking densities, demanding elevated water exchange rates and increased levels of dissolved oxygen [2].

Generally, weaned juveniles (5 g) are pre-ongrown in land-based tanks for 8–10 months up to reaching 50–60 g, prior to being transferred to sea cages for the final grow out. In recent years, there has been a great deal of uncertainty with regards to the economic profitability for many gilthead seabream farmers. Since production costs often outweigh the profits due to a rapid market expansion in the 1980s, with consequent oversupply, a lower market value has been established in the last two decades [3]. The presence of skeletal anomalies in gilthead seabream production significantly hinders profitability [4]. Considering that gilthead seabream hatcheries produce varying amounts of suboptimal quality juveniles (i.e., with externally detectable skeletal anomalies), ongrowing farms are usually willing to pay premium prices for high quality fry. Unfortunately, this extra cost is not always effective: firstly, certain deformities are undetectable in early life stages, only becoming evident later on [5]; secondly, skeletal anomalies can arise at varying points throughout the life cycle [6].

Several factors may augment the economic uncertainty for farmers. The presence of deformed fish in farm products provokes higher management expenditures (energy, oxygenation, workforce, feeding for rearing deformed fish that are destined to die or be thrown out prior to being marketed), lower market value for suboptimal products, and/or the need for post-harvest processing to other less-marketable forms (fillet or pet food) [7–9].

The rearing cycle of pre-ongrowing is typically studied as an integrated system coupled with the ongrowing phase for economic, growth performance, or qualitative studies. The pre-ongrowing phase is rarely investigated disjointedly from ongrowing, although such studies could offer unique insights and opportunities to improve rearing strategies during an extremely pertinent phase. In ongrowing farms, the primary production costs come from the feed (46%) and fry (14%) [10]. These data highlight the importance of pre-ongrowing fry quality as they represent a significant cost for producers. For example, the occurrence of dental prognathism (underbite), which was found to affect 29% of juvenile gilthead seabream, had increased to 57% in the same lot once they had reached market-size [11]. Consequently, the high feed cost mentioned for grow-out farms could be effectively compensated for by identifying rearing procedures capable of reducing the percentages of deformed fish, thereby increasing the final marketable products.

Taking into consideration that the pre-ongrowing phase is mainly carried out in land-based facilities, it represents the last opportunity for farmers to apply quality controls and cull out deformed fish before their transfer into sea cages. Furthermore, the pre-ongrowing phase enables a more effective sorting/quality guarantee than culling conducted during or at the end of the hatchery phase. In fact, while fish with severe cranial anomalies could easily be culled out at the end of larval rearing due to the precocity of cranium osteogenesis, vertebral axis deformities are rare (or lethal if due to precocious notochordal defects) at this
stage, often arising in more advanced developmental stages. Currently, routine practices
during pre-ongrowing entails regular culling to both reduce the stocking density and to
eliminate deformed fish from the production cycle [12]. This strategy can inflict additional
stress on the animals [13] and requires a great amount of technical labor and know-how by
the technicians.

Pre-ongrowing in tanks is generally carried out in intensive conditions that require frequent tank-cleaning as well as high water exchange rates and oxygen saturation [2]. While the pros and cons of choosing intensive versus semi-intensive rearing are weighed out by the farmers based primarily on the space availability, farm dimensions, and willingness to wait for economic returns (i.e., the quantity and size of fish that can be produced in the shortest time while minimizing the production costs), other incentives are gaining in popularity and economic competitiveness in European aquaculture. In particular, improved animal welfare, reduction of feed, water consumption, labor and energetic costs, organic labeling, and added-value marketing (gained by labeling schemes for high-quality products) all represent modernized approaches to animal production that can offer farmers a competitive edge and niche-market opportunities [14]. Obviously, the principal goal of producing both a high quantity and quality of fish remains paramount for the farmers. In order to efficiently optimize the desired quantity with the ideal quality, the best rearing parameters must be identified for each farmed species and life stage.

One particular biological issue that presents a significant loss in market value and
overall profits is the presence of skeletal deformities in reared fish. In fact, in the European
Union economic losses due to skeletal anomalies in fish farming were estimated to be
around 50 million euros annually [4]. This value is likely the same or higher nowadays,
as the problem has persisted [15] and oversupply/low demand in recent years is forcing
farmers to optimize internal production efficiency to improve their profit margins [3]. The
development of skeletal anomalies is undeniably complex due to the many etiologies
and synergistic factors eliciting the development of skeletal anomalies [15-17]. Presently,
many studies have been conducted to improve skeletal quality and reduce the presence
of skeletal anomalies in gilthead seabream by using tailored rearing strategies to target
life-stage specific requirements. Previous works among the main branches of research have
focused their efforts on investigating the influence of dietary/nutritional factors [18-24],
genetic components [25-27], and environmental dynamics [28-32] on the development
of anomalies.

However, these strategies are not always applied to the actual production due to various reasons, such as expensive live or inert feed particularities (i.e., floatability, detectability, digestibility, etc.), tank/cage shape and water flow diversity, and novel environmental or technical challenges.

Nevertheless, two environmental parameters can be easily adapted and optimized
by gilthead seabream farmers: the stocking density and the available swimming space.
Previous studies have looked at stress parameter responses of different rearing densities in
gilthead seabream [33-38], and others have compared the effects on skeletal anomalies of
semi-intensive ('Large Volumes' and 'Mesocosms') to intensive rearing strategies [39-41].
Two recent studies on model fish species have also revealed that higher stocking densities
provoke an increase in the presence of skeletal anomalies [42,43]. In a previous experiment
carried out by Dellacqua and colleagues [44], it was found that reducing the stocking
density during the hatchery phase improved the survival, increased final juvenile size, and
reduced the presence of skeletal anomalies, while the increased water volume (swimming
space) had a secondary positive effect by reducing the presence of jaw anomalies.

In this scenario, this study represents the first study aimed at delineating which,
between tank volume and stocking density, is the main driver in the formation and prolif-
eration of skeletal anomalies during the pre-ongrowing phase of gilthead seabream. The
experimentation was designed to mimic industrial rearing management in order to be
directly applied at a commercial scale. The expected results could greatly contribute to the
optimization of gilthead seabream farming within the scenario of improving sustainability,
profitability, and animal welfare. In order to facilitate the implementation of improved rearing strategies presented in this work, this study was conducted (i) at a pilot scale, following the standardized management practices typically applied during the pre-ongrowing phase; (ii) by testing three stocking densities (5, 10, and 20 kg/m³) suggested by the Italian Fish Farmers Association (API) and (iii) two tank volumes (500 L and 1000 L), which would typically be available for on-site use in any commercial gilthead seabream farm.

## 2. Materials and Methods

Experimental rearing was conducted at the Portuguese Institute of the Sea and Atmosphere (IPMA) facilities in Olhão (Portugal) during the late summertime for 63 days from August to October using gilthead seabream produced at the facility. Animal handling and experimentation was directed by trained scientists (following FELASA category C recommendations) and conducted according to the European guidelines on protection of animals used for scientific purposes (Directive 2010/63/UE of European Parliament and of the European Union Council), and the related guidelines and Portuguese legislation (Decreto- Lei 113/2013) for animal experimentation and welfare. IPMA's Aquaculture Research Station was certified by the Direção Geral Direcção de Veterinária to conduct animal experimentation under the authorization 2018/12/17—025516.

### 2.1. Rearing Setup

Gilthead seabream juveniles (average weight of 6.7 ± 2.5 g and length of 7.8 ± 1.1 cm), obtained from different spawns and parents, were previously selected for skeletal anomalies based on an external inspection carried out by experienced technicians. The 'non-deformed' fish (hereafter referred to as T₀) were stocked at 3 different densities LD: 5 kg/L; MD: 10 kg/L; HD: 20 kg/L in both 500 L and 1000 L water volumes, with a replacement rate of 100% tank water volume/hour, using sand-filtered natural seawater. The tanks for the 6 different experimental conditions were labeled LD1000, MD1000, HD1000, LD500, MD500, and HD500. The cylindrical tanks were all of the same diameter (1 m); however, they had different heights to hold two different water volumes (1000 L vs. 500 L). The stocking densities were periodically controlled during the experiment by sampling and weighing fish from each tank and removing excess fish to maintain the original biomass/L for each of the respective densities in the lots (Supplementary Table S1). The light regimen was based on natural light and photoperiod in the summertime in Olhão. The water temperature in the tanks ranged from 19–29 °C throughout the ~63 days of rearing (Supplementary Figure S1a). Oxygen diffusers were placed at the bottom of the tanks in order to maintain an O₂ saturation level above 70% and DO > 5 mg/L, which was checked several times daily (Supplementary Figure S1b). Fish were fed 3% of their body weight/day with AquaGold 5 pellets (AquaSoja, Ovar, Portugal).

### 2.2. Sampling and Analyses

A sample of 161 T₀ specimens were euthanized with 700 ppm of 2-phenoxyethanol and fixed in 1.5% PFA and 1.5% glutaraldehyde in 0.1 M sodium cacodylate buffer (pH 7.4), then transferred to 70% ethanol and successively radiographed (Gilardoni CPX 160/4 System, Unleaded film Kodak Mx 125, 55 Kv, 4 mA) for anatomical analyses.

During the trial, a total of 4769 individuals were sampled along 3 different time points from all of the experimental lots and weighed (wet weight—WW, weighed immediately once fish were anesthetized with 200 ppm of 2-phenoxyethanol) to calculate the growth in biomass and subsequently the density. After ~63 days of experimental rearing (T<sub>F</sub>), a total of 853 fish were euthanized following the same protocol that was applied to the samples from T₀. T<sub>F</sub> fish were measured for fork length, total length (TL), and wet weight (WW) (Supplementary Table S1). The Fulton's condition factor (K) was determined on T<sub>F</sub> samples based on the formula: K= 100 × W<sub>t</sub>/L<sup>3/2</sup>, in which W<sub>t</sub> represents the individual's WW, while L represents the TL.
The Food Conversion Ratio (FCR), Specific Growth Rate (SGR), and % Weight Gain were calculated using the T₀ and T_f weights according to formulas described in [45].

Specimens were radiographed with a digital DXS 4000 Pro X-ray (Carestream) (0.7 mm double top filter, 17 mm focal plane, with a single acquisition 120 s exposure time), and digitally visualized with the Carestream package at the University of Algarve. The anatomical investigation was carried out on the radiographs of T₀ and T_f samples.

Dissections of the vertebral column were made on a total of 12 lordotic and 12 normal T<sub>F</sub> fish individuated from the digital radiograph images, to assess potential differences in the calcium-phosphorous ratio between lordotic and non-lordotic vertebral centra (3 vertebral centra/fish). After removing attached soft tissues, the vertebrae were digested in nitric acid, and then the calcium and phosphorous mineral content was quantified using Microwave Plasma-Atomic Emission Spectroscopy (MP-AES; Agilent) at the University of Algarve.

Another 72 fish were also selected for specific anomaly types to be used for histology and sections of both deformed and non-deformed bones were dissected and were fixed as described above. Samples were then decalcified for 4 weeks by submersion in 10% EDTA and 0.5% PFA at 4 °C before being embedded in a GMA (glycol methacrylate) resin, according to the protocols described by Witten et al. [46]. Sections (5 µm) were cut on a standard rotary microtome (Microm HM360, Marshall Scientific (Hampton, NH, USA)) and stained with toluidine blue and Verhoeff's elastin stain.

Datasets from each experimental group were first tested for the normality of the distribution (Shapiro-Wilk, Anderson-Darling, Lilliefors L, and Jarque-Bera JB tests) and then subsequently tested either with a parametric ANOVA or a nonparametric Kruskal-Wallis test for equal medians and a Dunn's post hoc with Bonferroni correction.

## 2.3. Skeletal Anomalies Survey

Meristic counts and skeletal anomaly data (types and frequencies) were obtained from all of the 853 gilthead seabream individuals (T<sub>F</sub>) and 161 individuals (T<sub>0</sub>) by detailed examination of the radiographic images using FIJI [47]. Data was recorded using an alphanumeric code expressing the body region and the anomaly type (Supplementary Table S2). The examination was carried out based on the following assumptions: (1) incomplete fused bone elements were counted as discrete elements; (2) supernumerary bones with normal morphology were not considered to be anomalies but were included in meristic count variations; conversely, anomalously-shaped supernumerary elements were included among anomalies; (3) only clearly and unquestionably identifiable variations in shape were considered as skeletal anomalies: i.e., only the axis deviation associated to deformation of the vertebrae involved.

Skeletal anomaly data were then expressed in a raw matrix (RM) and used to calculate the frequencies (%) of each type of anomaly over the total number of anomalies, in each group (Supplementary Table S3). The RM was subsequently transformed into a binary matrix (BM), which was used to calculate the prevalence of individuals affected by each anomaly type (Supplementary Table S4). The following metrics were calculated for each group: (1) relative frequency (%) of individuals with at least one anomaly; (2) number of anomaly typologies observed; (3) average charge of anomalies (total amount of anomalies/number of malformed individuals); (4) relative frequency (%) of individuals with at least one anomaly; (5) frequency (%) of each anomaly type with respect to the total number of anomalies. The representative data from the RM and BM are expressed in tables and histograms. A Principal Component Analysis (PCA) was performed on the frequency of individuals affected by selected grouped anomaly types (Supplementary Table S5). The selection of the grouped anomaly types was based on the results obtained by testing the significance of differences in the frequencies of individuals affected by certain anomalies among the rearing conditions using a χ<sup>2</sup> test and pairwise post hoc [48].

Axial deviations were also described by measuring the angle of bending from the apex, located at the point of greatest angular flexion with vectors drawn to either the dorsal
or ventral (dorsal in the case of lordosis and ventral in the case of kyphosis) tip of the two vertebral centra preceding and the two successive centra following the apex angle (Figure 1). A modified categorization of angular class grouping made by Sfakianakis et al. [49] was applied to characterize 7 classes of axis deviations based on severity, according to the rankings shown in Table 1. A Kolmogorov-Smirnov (KS) test was applied to test for differences between the angular class frequency distributions among the lots.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Gilthead seabream individual exhibiting a hemal lordosis (classified as 'quite severe' with an angle represented by θ of 142.5°) and a hemal vertebral body fusion marked by the asterisk (*). Bar = 1 cm.</figcaption></figure>

<table><caption>**Table 1.** Ranking of different angular classes of axis deviation.</caption><tr><td>Most severe</td><td>&lt;137.9°</td></tr><tr><td>Quite severe</td><td>138.0°–144.9°</td></tr><tr><td>Severe</td><td>145.0°–151.9°</td></tr><tr><td>Moderate</td><td>152.0°–158.9°</td></tr><tr><td>Slight</td><td>159.0°–165.9°</td></tr><tr><td>Relatively mild</td><td>166.0°–172.9°</td></tr><tr><td>Normal</td><td>173.0°–180°</td></tr></table>

Statistical analyses and graphs were performed using Past 4.02 [50] and Python 3.8.

## 2.4. Geometric Morphometrics

The total of 853 (T<sub>F</sub>) digitally radiographed gilthead seabream were analyzed with geometric morphometrics [51] and used to investigate differences in body shapes between the T<sub>F</sub> individuals from different rearing conditions. Exactly 19 landmarks were manually selected on each digital radiograph, as displayed in Figure 2.
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Radiograph of gilthead seabream individual with selected landmarks superimposed. Morphological features associated with each landmark: 1: Frontal tip of premaxillary; 2: rostral head point in line with the center of the eye; 3: dorsal head point in line with the center of the eye; 4: distal edge of the 1st (anterior-most) predorsal bone; 5: insertion point of the 1st dorsal hard ray; 6: insertion of the 1st dorsal soft ray pterygophore; 7: insertion of the last dorsal soft ray pterygophore; 8: dorsal concave inflection-point of caudal peduncle; 9: middle point between the bases of hypurals 2 and 3 (fork); 10: ventral concave inflection-point of caudal peduncle; 11: insertion of the last anal pterygophore; 12: insertion of the 1st anal ray pterygophore; 13: insertion of the pelvic fin on the body profile; 14: preopercle ventral insertion on body profile; 15: frontal tip of dentary; 16: neural arch insertion on the 1st abdominal vertebral body (carrying the first ribs and following the 2 cephalic vertebrae); 17: neural arch insertion on the 1st hemal vertebral body; 18: neural arch insertion on the 6th hemal vertebral body; 19: middle point between the pre- and post-zygapophysis of the 1st and 2nd caudal vertebral bodies. Scale bar = 1 cm.</figcaption></figure>

A Procrustes transformation [52] was applied to adjust landmark configurations for centroid size, eliminating potential effects that are irrelevant to shape. A two-dimensional between-group PCA (bgPCA), which is analogous to a Relative Warp Analysis [53], and a Linear Discriminant Analysis (LDA) were applied to the transformed x and y coordinates, calculating between-group variance to visualize the effects of the rearing condition on body shape. An Eigenvalue scale was applied to the results of the PCA in order to enhance visualization [54]. Thin Plate Spline (TPS) interpolation diagrams were obtained from the residuals of the Procrustes transformation [55] and visualized as spline deviations from the consensus along the weight matrix (the components of the bgPCA). The LDA differences between experimental groups were tested with a MANOVA (Wilks' λ) [56]; the Mahalanobis squared distances were calculated for each pairwise group (i.e., LD1000 vs. LD500) and the differences between lots were tested for significance using a post hoc test with Bonferroni corrections.

**3. Results**

**3.1. Performance Indicators**

There were no sound differences in the survival, as mortality occurred with a total of 6 individuals (1 in LD1000, 3 in HD1000, 1 in MD500, and 1 in HD500).

There were no significant differences found in the final weight, FCR, SGR, nor % of Weight Gain (54.9 ± 1.9 g, 1.43 ± 0.20, 2.82 ± 0.07, and 7.16 ± 0.28%, respectively) among the tanks.

HD lots resulted in fish which were significantly (*p* < 0.05) smaller (TL) than the LD ones, regardless of the tank volumes (Figure 3a). Significant differences were also found in the Fulton's condition factor (K) between the three densities tested in 1000 L, and between
LD and the other two tested densities (MD and HD) in 500 L (Figure 3b). However, there were no significant differences in the final weight of fish between the conditions (Figure 3c).

<figure><img src="image_6.png" /><figcaption>**Figure 3.** (a) Box and whisker plots of TL of T<sub>F</sub> samples from the different experimental lots. (b) Differences in Fulton's K factor between conditions. (c) Final wet weights among the lots. Box indicates the central percentile, the line inside the box indicates the median while the whiskers represent the minimum and maximum values. Different letters indicate significant differences (*p* < 0.05, Kruskal-Wallis, Dunn's post hoc with Bonferroni corrections) among experimental lots.</figcaption></figure>

## 3.2. *Anatomical Analyses*

Due to the disparity in the sharpness of the images between the film (T<sub>0</sub>) and digital (T<sub>F</sub>) radiographs, in this study only vertebral bodies and specific splanchnocranial anomalies were considered.

The observed anomalies affecting the vertebral region were deformities of vertebral centra (fusions, partial fusions, and hemivertebrae) and axis deviations. The anomalies observed in the head were opercular plate anomalies (Figure 4b) and anomalies affecting the maxillary, premaxillary, and dentary, defined as underbite (Figure 4c) or overbite (Figure 4d). The frequencies of individuals which resulted to be affected by each type of anomaly are shown in Figure 5.

The T<sub>0</sub> lot resulted to display the highest anomalies' charge: an average of three anomalies were detected in each deformed fish (Table 2). Elevated incidences of individuals with anomalies, particularly lordosis (38% of individuals; Figure 5) in the hemal region (32% of individuals), were observed in this lot. Furthermore, vertebral centra fusions affected a greater frequency (17%) of T<sub>0</sub> than the T<sub>F</sub> individuals (Figure 5).

<table><caption>**Table 2.** Metrics calculated at T<sub>0</sub> and in the experimental lots.</caption><thead><tr><th>Lots</th><th>T<sub>0</sub></th><th>LD1000</th><th>MD1000</th><th>HD1000</th><th>LD500</th><th>MD500</th><th>HD500</th></tr></thead><tbody><tr><td>Number of observed specimens</td><td>161</td><td>140</td><td>115</td><td>197</td><td>124</td><td>96</td><td>181</td></tr><tr><td>% of individuals with at least one anomaly</td><td>79</td><td>63</td><td>69</td><td>80</td><td>69</td><td>75</td><td>88</td></tr><tr><td>Anomalies charge (n. of anomalies/affected fish)</td><td>3</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td><td>2</td></tr><tr><td>Observed types of anomalies</td><td>21</td><td>18</td><td>23</td><td>23</td><td>17</td><td>21</td><td>23</td></tr></tbody></table>
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Splanchnocranium anomalies affecting gilthead seabream: (a) normal jaws and opercular plate; (b) reduced opercular plate, the arrow indicates a gap exposing the gill opening; (c) underbite: contemporary deformation of premaxillary and maxillary coupled with a protrusion of the dentary; (d) overbite: contemporary deformation (twisting, in this case) of the premaxillary and maxillary coupled with a reduction in length of the dentary resulting in an upper-jaw overlap. Bar = 5 mm.</figcaption></figure>

Concerning the individuals in T<sub>F</sub> lots, the majority had at least one anomaly, with occurrences ranging from 63 to 88% (Table 2). A clear density effect is discernable by the fact that both of the HD conditions displayed the highest percentages of individuals with at least one anomaly. However, for this metric, only HD500 resulted to be significantly different than the other conditions in a χ<sup>2</sup> pairwise post hoc test (*p* < 0.001). A similar effect of density was recognizable in the number of observed anomaly types.

A secondary volume effect was also detectable: the 500 L lots displayed higher frequencies of deformed individuals than the 1000 L lots of equivalent densities even though these values were not significantly different (HD500 > HD1000, MD500 > MD1000, and LD500 > LD1000; Table 2).
<figure><img src="image_8.png" /><figcaption>**Figure 5.** Frequency of individuals affected by the considered skeletal anomalies in the vertebral column and the cranium of gilthead seabream individuals evaluated from radiograph images in T₀ and in each experimental density and volume condition (LD1000, MD1000, HD1000, LD500, MD500, HD500) at T_F.</figcaption></figure>

The following sections report detailed descriptions of the observed anomalies in the different body regions.

### 3.3. Cephalic Skeleton

The T<sub>F</sub> samples revealed higher frequencies of individuals exhibiting underbites, particularly in the HD500 lot (36.5%). The incidences of underbites among the conditions demonstrated a clear density effect as the frequency of affected individuals linearly increased with increasing density for both volume treatments (LD < MD < HD; Figure 6a). Additionally, a volume-effect was revealed by an augmented occurrence of affected fish in the 500 L volume conditions compared with their 1000 L equivalent density counterparts (LD1000 < LD500, MD1000 < MD500, HD1000 < HD500; Figure 6a). Although rarer overall, even overbites tended to be more frequent when the density was greater in the 500 L conditions (LD < MD < HD; Figure 6b), whereas the frequency in the LD1000 lot was less than those observed in the other two 1000 L lots (HD and MD), which displayed nearly equal frequencies (LD < MD ≅ HD; Figure 6b).

Conversely to underbites, the incidences of gilthead seabream with overbites did not display the same volume effect, as the frequencies varied haphazardly and the differences were relatively small (LD1000 < LD500, MD1000 > MD500, and HD1000 < HD500; Figure 6b). Nevertheless, the differences witnessed in the frequencies of individuals affected both by underbites and overbites did not result to be significant among the experimental lots (Pearson's χ² test: *p* values = 0.2 and 0.09, respectively).

In the case of opercular anomalies, which were inspected only from the left side of the fish, greater frequencies of fish affected by opercular anomalies were associated with greater densities and discernable in both volume groups (LD < MD < HD; Figure 6c). Analogous to what was found in the case of underbites, a secondary volume effect was also detectable as the smaller volume lots displayed higher incidences than the large volume density equivalents (LD1000 < LD500, MD1000 < MD500, and HD1000 < HD500; Figure 6c). Conversely to the jaw anomalies, these differences in frequencies of fish affected by opercular anomalies resulted to be significant between LD1000 vs. HD1000, LD1000 vs. MD500, LD1000 vs. HD500, and MD1000 vs. HD500 (Figure 6c).
<figure><img src="image_9.png" /><figcaption>**Figure 6.** Frequencies of gilthead seabream individuals with (a) an underbite (%); (b) an overbite (%); and (c) a left-opercular plate anomaly (%), among the experimental conditions. Neither (a) nor (b) resulted to have any significant differences, while in (c) asterisks indicate pairwise significant differences (χ² pairwise test with Bonferroni corrections, * = *p* < 0.05; ** = *p* < 0.005).</figcaption></figure>

### 3.4. Axial Skeleton

Vertebral body counts ranged from 22 to 25 in the LD1000 and MD1000 lots, while in the other lots the counts ranged from 22 to 24 (Supplementary Figure S2). When analyzing increasing densities, a tendency towards lower vertebra counts appears (Supplementary Figure S2) mainly due to higher occurrences of vertebral centra fusions, recognizable by the presence of a fully remodeled (amphicoelous) centrum and the presence of more than one neural or hemal arch. Specimens with vertebrae fusions and lower vertebral counts in the LD and MD varied between 0–16.7%, conversely to the HD lots displaying higher frequencies (28% and 52.6% for HD1000 and HD500, respectively). However, the
overall frequency of individuals exhibiting fused vertebrae regardless of the vertebral body count did not vary greatly between the lots (8.6–11.7%). Although lower percentages were observed in the LD lots with respect to MD and HD ones (Figure 7a), these differences were not significant among the lots (χ² test).

<figure><img src="image_10.png" /><figcaption>**Figure 7.** Frequencies of gilthead seabream individuals with (a) fusions of vertebral body centra (%); (b) lordosis in the vertebral axis (%); and (c) kyphosis in the vertebral axis (%), among the experimental conditions. (a) resulted to not have any significant differences, while in (b,c) asterisks indicate pairwise significant differences (χ² pairwise test with Bonferroni corrections, * = *p* < 0.05; ** = *p* < 0.005).</figcaption></figure>

Concerning axis deformations, 42% of all individuals presented lordosis and 6% kyphosis. The histogram of the frequencies shown in Figure 7b displays a progressive increment of individuals with lordosis associated with increasing density. The highest frequency of lordotic individuals was found in HD1000 (52.3%) followed by HD500 (50.3%), while the other conditions displayed lower frequencies which were in the same order of magnitude for each density group, indicating a clear density effect, but no secondary volume effect (HD1000 (52.3%) > HD500 (50.3%) > MD1000 (37.4%) ≅ MD500 (37.4%) > LD1000 (31.4%) ≅ LD500 (32.3%); Figure 7b). Differences between the frequencies of LD and HD individuals affected by lordosis were statistically significant based on pairwise comparisons (Figure 7b).
The region most commonly affected by lordosis among all of the conditions was the hemal region with the HD1000 lot, followed by the H500 lot, displaying the greatest frequency of affected individuals (49.2% and 44.2%, respectively; Figure 8). Furthermore, the severity of lordosis, expressed by angular classes (Table 1), resulted to be significantly different among the experimental lots (Figure 9; Table 3), although no clear trend was discernable.

<figure><img src="image_11.png" /><figcaption>**Figure 8.** Subdivisions of vertebral regions in which individual gilthead seabream were affected by lordosis among the different experimental lots.</figcaption></figure>

<figure><img src="image_12.png" /></figure>

<figure><img src="image_13.png" /><figcaption>**Figure 9.** Distribution of different angular classes of lordotic curvatures based on the classes described in Table 1: (a) in the lot T₀ at the beginning and (b) the different T<sub>F</sub> lots of different densities and volumes at the end of the experimental rearing. Significant differences from the KS pairwise post hoc are described in Table 3.</figcaption></figure>
<table><caption>**Table 3.** *p*-values obtained from a Kolmogorov–Smirnov (KS) pairwise post hoc test with distribution of lordosis severity (excluding non-deformed individuals). Only significant (*p* < 0.05) comparisons are reported.</caption><thead><tr><td>Pairwise Post Hoc</td><td><i>p</i>-Values</td></tr></thead><tbody><tr><td>LD1000 vs. HD1000</td><td>0.002</td></tr><tr><td>MD1000 vs. HD1000</td><td>0.002</td></tr><tr><td>HD1000 vs. LD500</td><td>0.002</td></tr><tr><td>HD1000 vs. MD500</td><td>0.02</td></tr><tr><td>LD500 vs. HD500</td><td>0.02</td></tr></tbody></table>

Kyphotic gilthead seabream were less frequent than the lordotic ones; however, an effect of density was detected between all of the 1000 L lots (LD < MD < HD), while in the 500 L conditions the trend is upheld only between HD and the other two conditions (LD500 (2.4%) ≅ MD500 (2.1%) < HD500: 12.2%; Figure 7c). Overall, the highest prevalence of individuals displaying kyphosis was found in the HD500 (12.2%) followed by the HD1000 (8.6%), and the χ² pairwise post hoc test with Bonferroni corrections confirmed that significant differences were found between LD1000 vs. HD1000, LD1000 vs. HD500, HD500 vs. MD1000, HD500 vs. LD500, and HD500 vs. MD500 (Figure 7c). As far as the effect of water volume (1000 L vs. 500 L conditions with the same density) on the frequencies of kyphotic individuals is concerned, none of the differences observed resulted to be statistically significant (Figure 7c). Neither were there any significant differences in the distributions of severity classes (Supplementary Figure S3), contrary to lordosis.

### 3.5. Histological Evaluations

Histological analyses highlighted interesting features of some vertebral anomalies. Partial (Figure 10) and complete vertebral body fusions display a pronounced dorsal-ventral asymmetry. Presumably, Figure 10a is representative of an early stage hemal fusion, in which the hemal arches begin to first fuse, with a narrowing of the ventral intervertebral space, which still persists (Figure 10b,c). Lateral bending in the orientation of the collagen fibers constituting the trabeculae was evident in the fusing centra in Figure 10b as well as in Figure 11b. In a complete fusion of hemal centra (Figures 11 and 12c,d), the components of the intervertebral space (the collagen type-II based notochord sheath and its outer elastin layer) are still present ventrally, whereas modelling has already progressed dorsally (Figure 12c,d). In the case of hemal lordosis, the histological investigations reveal morphological changes in the shape of the centra, confirming the radiographic observations. All the components of the intervertebral spaces, neural, and hemal arches are maintained and healthy (Figure 13).

### 3.6. Principal Components Analysis

A PCA was applied to six anomaly types in order to enhance the visualization of differences between the conditions via an ordination plot. The anomalies considered were underbites, overbites, opercular anomalies, fusions, lordosis, and kyphosis (Figure 14). This PCA explained a total of 91.31% of the variance in the first two axes (76.62% and 14.69%, respectively for the first and second axes). In this ordination, a clear effect of density is distinguishable: the HD conditions were positioned towards the negative direction of component 1, in the same direction of the skeletal category vectors; MD conditions were positioned towards the center of the plane and LD groups plotted in the positive direction of component 1. Although component 2 explains less of the variance (14.69% vs. 76.62%), it hints at the secondary volume effect, detectable by the positioning of 500 L lots (green circles) in the negative hemiplane of component 2 (in proximity of the opercular and underbite anomaly vectors), in opposition to the 1000 L lots.
<figure><img src="image_14.png" /><figcaption>**Figure 10.** Compressed/partially fusing vertebral centra with an already fused dorsal zygopophysis. (a) Radiographic image of a T₀ seabream exhibiting two compressed vertebral centra (rectangle). Bar = 5 mm; (b,c) Toluidine blue-stained sections of the deformed vertebral bodies shown in (a). The black arrowheads in (b) point to bending trabeculae associated with remodeling of the misshapen centra. Note the dorsal-ventral asymmetry, where the ventral regions of the two centra appear enlarged. Bar = 1 mm; (c) higher magnification of the region indicated with the black dashed rectangle in (b): the asterisk marks the insertion point of the hemal arches in the ventral intravertebral space between the two centra. The black arrowheads point to the bone of the hemal arches. Bar = 0.1 mm.</figcaption></figure>
<figure><img src="image_15.png" /><figcaption>**Figure 11.** Complete fusion of vertebral bodies. (a) Radiograph image of a T₀ seabream exhibiting a multi-vertebra fusion (rectangle). Bar = 5mm; (b,c) Toluidine blue-stained sagittal section of the multiple fusion event shown in (a), photographed with polarized light. The black arrowheads in (b) point to laterally bending trabeculae that surround the fusing vertebral centra. The presence of intravertebral space between the trabeculae of the fusing centra (asterisk) indicates a lateral deviation of the fusing vertebral bodies as in a localized scoliosis. The white arrowhead points to a normal vertebral body endplate. Bar = 1mm; (c) the black arrowheads point to putative reactive scar tissue of the notochord strand (white arrowheads), occurring proximal to the point of the fusing centra (asterisk). Note the cartilage (in pink) of the hemal and neural arches, located both dorsally and ventrally to the fusing centra above and below the asterisk in (c). Bar = 1mm.</figcaption></figure>
<figure><img src="image_16.png" /><figcaption>**Figure 12.** (a) Non-fused vertebral centra and, (b) multiple fusions among preural centra a, (c) Toluidine blue staining; (b,d) Verhoeff's elastin-staining. The black arrowhead in (c) indicating that all of the components of the intravertebral ligaments are maintained; and in (d) the arrowhead indicates the persisting elastin layer in the fused intervertebral space. Bar = 0.5 mm.</figcaption></figure>

<figure><img src="image_17.png" /><figcaption>**Figure 13.** (a) Radiograph of a lordotic seabream (the white rectangle highlights the localized affected area). Bar = 5 mm; (b) Toluidine blue parasagittal section photographed with polarized light of the lordotic vertebral body centra shown in a. The double arrows highlight the different widths of the dorsal profile of lordotic vertebrae. Bar = 2 mm.</figcaption></figure>
<figure><img src="image_18.png" /><figcaption>**Figure 14.** PCA ordination model applied to the frequencies of individuals affected by underbites, overbites, opercular anomalies, fusions, lordosis, and kyphosis in each experimental condition. The resulting ordination model indicates that the high-density conditions were highly associated with these skeletal anomalies.</figcaption></figure>

## 3.7. Geometric Morphometrics

A bgPCA was applied to the residuals from the Procrustes transformation derived from landmarks annotated on the radiograph images. In this PCA, the ordination model revealed that inter-group differences can be explained by 62.98% of the variance on the first axis and 24.01% of the variance on the second axis (Figure 15a). The LD lot centroids are both located on the positive plane of component 1, while the other lots are positioned on the opposite space in the negative hemi-plane of component 1. The component 2 enhances differences between the LD500 and LD1000 (quadrant 1 and 4, respectively) and between the HD (quadrant 2) and MD (quadrant 3) conditions (Figure 15a). Furthermore, lots with the same density and different volumes are positioned proximally to each other, suggesting a sound and strong effect of rearing density on fish shape. Thin plate spline (TPS) deformation grids, which were superimposed on the PCA in Figure 15a, assisted in further discerning shape differences among the conditions. The TPS deformation grid in the positive hemi-plane of component 1 (proximal to the LD lot centroids) shows that LD individuals exhibited a slenderer shape than the roundish shape of the MD and HD individuals located on the negative hemiplane of component 1 (Figure 15a). In particular, the differences reside mainly from landmark convergence within the caudal peduncle and divergence in the frontal-rostral region.

The LDA (Linear Discriminant Analysis) ordination explained 54.5% and 27.1% of the total variance in the first two axes, respectively (Figure 15b). Additionally, the LDA further enhanced the detection of shape differences between the groups displaying a strong effect of density as the LD condition centroids are positioned in the second and the third quadrants, in the negative hemi-plane of LD1, while the MD and HD centroids are located in the positive hemi-plane of LD1, in the first and the fourth quadrants, respectively. The LDA confirms that the centroids of the lots reared with the same density are located proximally to each other, independently from the tank volumes.
<figure><img src="image_19.png" /><figcaption>**Figure 15.** (a) bgPCA scatter with deformation grids plotted on the extremities of the first component. TPS deformations are color-coded in which cool colors (blue/green) represent convergence with respect to the consensus and warm colors (yellow, orange, and red) represent divergence. (b) Linear Discriminant Analysis (LDA) scatter results with 95% ellipticals and lot centroids, highlighting group shape differences among the T<sub>F</sub> experimental lots of different densities and volumes.</figcaption></figure>

Lastly, a MANOVA (Wilks' λ) applied to the uniform and non-uniform scores derived
from the weight matrix revealed significant inter-group differences, while a pairwise
test between the Mahalanobis squared distances validates differences between particular
conditions (Table 4, Bonferroni corrected *p*-values, *p* < 0.00001). In particular, the greatest Mahalanobis squared distance was 14.9, which was found between LD500 and HD500.

<table><caption>**Table 4.** Mahalanobis squared distances inter-groups.</caption><thead><tr><th></th><th>MD1000</th><th>HD1000</th><th>LD500</th><th>MD500</th><th>HD500</th></tr></thead><tbody><tr><th>LD1000</th><td>6.7184</td><td>6.5169</td><td>4.0499</td><td>11.3400</td><td>8.7823</td></tr><tr><th>MD1000</th><td></td><td>2.5247</td><td>9.1422</td><td>3.6365</td><td>3.9081</td></tr><tr><th>HD1000</th><td></td><td></td><td>11.321</td><td>6.9757</td><td>2.9881</td></tr><tr><th>LD500</th><td></td><td></td><td></td><td>11.575</td><td>14.898</td></tr><tr><th>MD500</th><td></td><td></td><td></td><td></td><td>5.905</td></tr></tbody></table>

Significantly different distances were found between all groups (Bonferroni correction, *p* < 0.0001).

### 3.8. Calcium-Phosphorous Mineral Content

Results from the mineral quantification revealed the absence of significant differences in the Ca/P ratios between malformed lordotic vertebral bodies and normal vertebral bodies, although a greater variability of Ca/P ratios from lordotic vertebral centra than normal centra is discernable (Figure 16).

<figure><img src="image_20.png" /><figcaption>**Figure 16.** Calcium-phosphorous ratio in ppm (µg/mL) of Ca to ppm (µg/mL) P between lordotic and normal centra (three vertebral bodies). Each dot is indicative of a single individual while the central line denotes the median. (ANOVA, NS).</figcaption></figure>

### 4. Discussion

This study has demonstrated that pre-ongrowing gilthead seabream stocked in low densities (LD) for two months were significantly longer than those reared in higher densities, regardless of the water volume. The fact that individuals among the conditions were not heavier but presented significant differences in the condition factor K indicated that LD gilthead seabream acquired a slenderer profile than the fish reared in MD and HD, which appeared more roundish, as confirmed by the geometric morphometrics analyses. Both of these features could be due to the higher prevalence of lordotic fish found in high
density conditions since the presence of axis deviations affecting the length and the body
height could have given the deformed fish a stumpy body aspect. A previous study
investigating growth and stress parameters on adult (commercial size: 217 ± 1.9 g) gilthead
seabream reared for 6 weeks at 5, 10, and 20 kg/m³ found fish with greater weights
in the LD conditions than in the HD conditions [33], confirming that the density can
influence the size profile even in larger adults. The results from geometric morphometrics
demonstrated that after only 63 days of rearing gilthead seabream sub-adults in the six
different experimental conditions, significant differences in shape between the lots were
already established. This seemed to be a response to the stocking density, even though the
fish belonged to the same initial starting group of gilthead seabream from different origins.
This is striking since morphometric analyses typically enhance shape differences between
different populations and/or species from different ecological (i.e., feeding habits; ref. [56]),
evolutive, or phylogenetic contexts [57,58], or in the same species, between differences
in size classes [59], origins [60], life-stages, and/or long-term environmentally distinct
factors (i.e., between wild and reared gilthead seabream; ref. [61]). It is likely that the
differences we found after only 63 days of experimental rearing, in a relatively advanced
life-stage dominated by isometric growth in gilthead seabream [62], were possibly due
to the greater frequency of deformed individuals determined by HD conditions. The
graphical representation in Supplementary Figure S4, in which the spline of the centroid of
a non-lordotic LD1000 individual is compared to the TPS of a lordotic HD500 individual,
epitomizes this postulation.

The present study highlighted that when the effects of the rearing density were applied
and analyzed separately from those of the tank volume, both may have exerted some effects
on the skeletal phenotype in gilthead seabream. However, these effects varied according
to the skeletal region. For example, while an effect of stocking density was evident in the
frequency of opercular anomalies and axial deviations, no clear effects were identifiable
regarding the incidence of vertebral fusions.

Overall, the most frequently observed anomaly found in all of the lots was lordosis.
This anomaly is particularly problematic for gilthead seabream farms as its presence greatly
reduces the products' potential value. In this study, we found that the most common
region of the body affected by lordosis was the hemal region. This finding is in accordance
with other studies on the presence of hemal lordosis in gilthead seabream and has already
been recognized to be frequently present in reared fishes [63]. Furthermore, the severity
of lordosis can range from light deformations, which do not present any clear effects on
the external morphology, to severe, with apparent impacts on external shape [64]. The
increased severity of lordosis based on the angle of curvature was found to be significantly
related with the number of affected vertebral centra in European seabass [49], therefore
angular measurements can represent a reliable proxy for overall shape changes in the
vertebral column. The classification of lordosis severity could also have useful applications
for studies investigating recovery, due to recent findings that, in certain cases, affected fish
can recover from lordosis [65-67]. In this current study, the classification of lordosis severity
demonstrated that, at this life stage, gilthead seabream show a wide array of severity and
that differences among the distributions of severity classes significantly differed in response
to the density. This is the first study that has statistically confirmed differences between
severity classes, laying the foundations for future applications and studies of axial deviation
severity. The classification of different degrees of severity could also be useful for future
applications on the recovery potential. For example, farmers could screen for slight axis
deviations and apply different rearing strategies (such as reduced stocking density) to
effectively recover those fish before the severity increases beyond a threshold of no-return,
saving both time and money.

Lordosis had been initially proposed to be present only in association with improper
inflation of the swim bladder [68–70]. However, Boglione et al. [71] found lordotic seabass
juveniles with functional swim bladders, and Andrades et al. [72] described seabream
larvae displaying a lordotic notochord even prior to the normal inflation of the swim
bladder, suggesting that other etiologies may be responsible for eliciting this malformation. Some studies, in fact, have suggested that the primary causes for lordosis in reared fish are linked to the swimming activity [49,65,73,74] and that lordosis is an adapted response to increased mechanical loading exerted by trunk muscles [75]. This hypothesis seems to be corroborated by the healthy appearance of lordotic vertebrae from our histological analyses. Additionally, the finding of no significant differences in the calcium-phosphorous ratio (Ca:P) between lordotic centra and normal centra suggests that lordosis could be either triggered by other metabolic factors or by biomechanical responses. Similarly, to our findings, Boursiaki et al. [76] found no differences in the Ca:P ratio between scoliotic and non-scoliotic centra in seabream individuals. In our study, we observed a variability in the values of Ca:P ratio in the lordotic vertebrae. Previous studies have found that increased flow improved the growth rates in seabream [32,77]. Enhanced growth due to exercise can occur when fish are subjected to higher water replacement regimes, with the exercised fish experiencing a replacement rate of 700 L/h which establishes hyperplasia and hypertrophy in the white muscles [78,79]. Although the growth was improved in higher flow regimes in previous studies, a recent study by Palstra et al. [32] found an association with high flow regimes (~0.78 m/s) and increased frequencies of lordotic gilthead seabream. While in this current study an increased inlet flow velocity in the larger volume tanks was necessary to maintain the same water exchange rate as the small volumes, there were no differences in the frequencies of lordotic seabream between the different volume (flow) conditions. The water replacement rate and flow velocities tested in this current work, however, were lower than the replacement rate and inlet velocities from these previous studies (current study: maximum water exchange rate: 100%/h; flow velocities of ~0.73 and ~0.31 m/s for 1000 L and 500 L tanks at the inlets, respectively). Regardless, this potential effect of flow did not reveal itself in our experimental rearing since in cases of lordosis, no differences were present based on the different volumes (nor, consequently, flow velocity).

Contrary to water volumes, the stocking density did play a driving role in the prevalence of lordosis. The HD conditions could have influenced the fish behavior, based on the increased interactions among more individuals within the tanks (avoidance, escape, competition, aggression, hierarchies, etc.). Interestingly, in a previous work on the behavioral response of gilthead seabream to low stocking density, it was found that individual fish exhibited less overall movement within the tanks but a greater diversity in their responses to behavioral and social tests [80], while tighter interactions with consequent rapid changes of directions and burst swimming were more frequent in high density conditions. Burst swimming is characterized by high speeds maintained for a few seconds (<20 s) [81], rapid changes of direction by angular bending in the hemal region, and energy primarily supplied to myotomal white muscle through anaerobic processes [82]. All of these features yield potentially dire consequences for the (hemal) musculoskeletal system [34], with an onset of axis deviations in the hemal region induced by higher biomechanical pressure of muscle on underlying vertebrae.

Nevertheless, what has been discussed above should be considered in the context that the initial lot, T₀, displayed a high presence of anomalies, whose occurrence in some cases were found to be reduced, and in others augmented in the final T<sub>F</sub> seabream. The consequence is that the quality of the fish used to start the ongrowing phase define the quality threshold from which amelioration or worsening can be obtained by the predominant modulating factor: rearing density.

**5. Conclusions**

In the pre-ongrowing phase, the results of a primarily strong effect of density followed by a secondary benefit of increased water volume confirms previously found results regarding the incidences of skeletal anomalies in response to the rearing density during the hatchery phase [44]. Despite the fact that the life stage, production phase, genetic origin, and experimental location were different, some effects of the same physical driver (density) were identifiable in the same skeletal elements. Bear in mind that the differences in the life stage
of fish under experimental trials were not negligible: the larval stage represents a critical and unique life phase, characterized by ontogenetic phases during which the differentiation of larval organs into the definitive (juvenile/adult) organs develop in succession, with consequences affecting the allometric/isometric growth, changes in habitat, trophic ecology, feeding modality, type of swimming, social responses, and behavior. The life stage under examination in the present study, the pre-ongrowing phase, is characterized by a more canalized [83] and exclusively isometric growth since differentiation has already been completed [62]. The genetic differences which characterize the gilthead seabream collected from different lots of different ages that only share the same weight-range are dramatic when compared to the same-batch origin of the specimens analyzed in the hatchery study.

The gilthead seabream skeleton seems to react to challenging environmental conditions by developing similar homeorhetic trajectories for modeling processes affecting the same bones. In particular, in both of the production/life phases, it was found that the use of higher densities determines higher occurrences of opercular, jaw, and vertebral axis anomalies. Additionally, it was found that the use of larger water volume mitigates the incidences of jaw anomalies, in both of the stages. However, considering that differences in occurrences of these anomalies were not significant, this study clearly individuates the stocking density as the main factor modulating the development of specific skeletal anomalies in gilthead seabream.

Future research carried out at different levels (biochemical, behavioral, and epigenetic) and on different species could actively help in unveiling the underlying mechanisms leading to the augmentation of anomalies in high density conditions.

**Supplementary Materials:** The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ani13040557/s1, Figure S1: (a) Average daily temperature throughout rearing. (b) Average daily dissolved oxygen throughout rearing; Figure S2: Frequency classes of vertebrae counts in the experimental lots.; Figure S3: Distribution of different angular classes of kyphotic curvatures among the density and volume conditions.; Figure S4: Normal LD1000 gilthead seabream TPS on left and lordotic HD500 gilthead seabream TPS on right.; Table S1: Pre-ongrowing rearing established tank densities and performative parameters; Table S2: List of considered anomaly types; Table S3: Relative frequency of specific skeletal anomalies among the conditions; Table S4; Frequency of individuals exhibiting specific skeletal anomalies among the conditions; Table S5: Descriptors used to preform PCA ordination model (see Sup. Table S2 for the codes of the anomalies).

**Author Contributions:** Z.D.: experimental setup, running experiment, sampling, data collection, data curation, and statistical analyses, formal analysis, writing of manuscript. C.D.B.: histological processing and interpretation. C.C.: geometric morphometrics data processing, statistical analyses, and interpretation. P.P.-F.: experimental setup, equipment and facilities, supply of fish. L.R.: experimental setup, facilities, sampling, data collection. M.B.: experimental setup, technical training, sampling, data collection. P.J.G.: Facilities, data collection, technical training. F.M.: statistical analyses and interpretation. A.F.: experimental design. M.I.: interpretation of results and editing manuscript. C.B.: Experimental design, data collection, training on data curation, interpretation of results, writing and editing of manuscript. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was funded by European Union's Horizon 2020 research and innovation program 'BIOMEDAQU' under the Marie Skłodowska-Curie grant agreement No.766347. CCMAR receives funding from the Portuguese Foundation for Science and Technology (FCT) through projects UIDB/04326/2020, UIDP/04326/2020, and LA/P/0101/2020. IPMA also receives funding from the MAR2020 program with the DIVERSIAQUA project (Mar2020-P02M01-0656P).

**Institutional Review Board Statement:** Animal handling and experimentation was directed by trained scientists (following FELASA category C recommendations) and conducted according to the European guidelines on protection of animals used for scientific purposes (Directive 2010/63/UE of European Parliament and of the European Union Council), and the related guidelines and Portuguese legislation (Decreto-Lei 113/2013) for animal experimentation and welfare. IPMA's Aquaculture Research Station was certified by the Direção Geral Direcção de Veterinária to conduct animal experimentation under the authorization 2018/12/17-025516.
**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The data presented in this study are available in the supplementary materials. Additional data presented in this study are available on request from the corresponding author. The data are not yet publicly available due to unfinalized database selection.

**Acknowledgments:** The authors would like to sincerely thank the entire team at IPMA in Olhão for their support. In particular, we would also like to acknowledge Leticia Luján and Cátia Marques for their help with the final sampling. Additionally, from UALG, we would like to express our gratitude to Nádia Silva for her technical and analytical support in conducting MP-EAS and to Alessio Carletti for his help in collecting/storing radiograph images. We would also like to express our gratitude to Eckhard Witten at the University of Ghent, as well as the research team there, for hosting and training us in the particularities and interpretations of histology applied to skeletal tissues. We would also like to express our gratitude to Arianna Martini for her help in annotating images. The authors would like to express our thanks to all of the reviewers that greatly helped in improving the manuscript.

**Conflicts of Interest:** The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

**Abbreviations**

<table><tr><td>bgPCA</td><td>Between group principal components analysis</td></tr><tr><td>BM</td><td>Binary Matrix</td></tr><tr><td>FCR</td><td>Food conversion ratio</td></tr><tr><td>HD</td><td>High Density</td></tr><tr><td>K</td><td>Fulton's condition factor</td></tr><tr><td>KS</td><td>Kolmogrov-Smirnov test for equal distribution</td></tr><tr><td>LDA</td><td>Linear discriminant analysis</td></tr><tr><td>LD</td><td>Low Density</td></tr><tr><td>MD</td><td>Medium Density</td></tr><tr><td>PCA</td><td>Principal components analysis</td></tr><tr><td>RM</td><td>Raw Matrix</td></tr><tr><td>SGR</td><td>Specific Growth Rate</td></tr><tr><td>T<sub>0</sub></td><td>Initial batch of gilthead seabream before being placed in experimental conditions</td></tr><tr><td>T<sub>F</sub></td><td>Final samples of gilthead seabream after being reared in experimental conditions</td></tr><tr><td>TL</td><td>Total Length</td></tr><tr><td>TPS</td><td>Thin plate splines</td></tr><tr><td>WW</td><td>Wet Weight</td></tr></table>

**References**

1. FAO Fisheries & Aquaculture—FishStatJ. Available online: http://www.fao.org/fishery/statistics/software/fishstatj/en (accessed on 22 October 2022).

2. Pavlidis, M.A.; Mylonas, C.C. *Sparidae*; John Wiley & Sons: Chichester, UK, 2011.

3. Llorente, I.; Fernández-Polanco, J.; Baraibar-Diez, E.; Odriozola, M.D.; Bjørndal, T.; Asche, F.; Guillen, J.; Avdelas, L.; Nielsen, R.; Cozzolino, M.; et al. Assessment of the economic performance of the seabream and seabass aquaculture industry in the European Union. *Mar. Policy* **2020**, *117*, 103876. [CrossRef]

4. Hough, C. *Manual of Control of Malformations in Fish Aquaculture*. *Science and Practice*; Federation of European Aquaculture Producers, RapidPRess: Luxembourg, 2009; p. 150.

5. Bardon, A.; Vandeputte, M.; Dupont-Nivet, M.; Chavanne, H.; Haffray, P.; Vergnet, A.; Chatain, B. What is the heritable component of spinal deformities in the European sea bass (*Dicentrarchus labrax*)? *Aquaculture* **2009**, *294*, 194–201. [CrossRef]

6. Daoulas, C.; Economou, A.; Bantavas, I. Osteological abnormalities in laboratory reared sea-bass (*Dicentrarchus labrax*) fingerlings. *Aquaculture* **1991**, *97*, 169–180. [CrossRef]

7. Le Vay, L.; Carvalho, G.; Quinitio, E.; Lebata, J.; Ut, V.; Fushimi, H. Quality of hatchery-reared juveniles for marine fisheries stock enhancement. *Aquaculture* **2007**, *268*, 169–180. [CrossRef]

8. Lijalad, M.; Powell, M.D. Effects of lower jaw deformity on swimming performance and recovery from exhaustive exercise in triploid and diploid Atlantic salmon *Salmo salar* L. *Aquaculture* **2009**, *290*, 145–154. [CrossRef]
9. Boglione, C.; Gavaia, P.; Koumoundouros, G.; Gisbert, E.; Moren, M.; Fontagné, S.; Witten, P.E. Skeletal anomalies in reared European fish larvae and juveniles. Part 1: Normal and anomalous skeletogenic processes. *Rev. Aquac.* **2013**, *5*, S99-S120.
[CrossRef]

10. Nielsen, R.; Ankamah-Yeboah, I.; Llorente, I. Technical efficiency and environmental impact of seabream and seabass farms. *Aquac. Econ. Manag.* **2020**, *25*, 106–125. [CrossRef]

11. Boglione, C.; Costa, C. Skeletal Deformities and Juvenile Quality. *Sparidae* **2011**, *233–294*. [CrossRef]

12. Koumoundouros, G.; Gagliardi, F.; Divanach, P.; Boglione, C.; Cataudella, S.; Kentouri, M. Normal and abnormal osteological development of caudal fin in *Sparus aurata* L. fry. *Aquaculture* **1997**, *149*, 215–226. [CrossRef]

13. Costa, C.; Antonucci, F.; Boglione, C.; Menesatti, P.; Vandeputte, M.; Chatain, B. Automated sorting for size, sex and skeletal anomalies of cultured seabass using external shape analysis. *Aquac. Eng.* **2013**, *52*, 58–64. [CrossRef]

14. Puszkarski, J.; Śniadach, O. Instruments to implement sustainable aquaculture in the European Union. *Mar. Policy* **2022**, *144*, 105215. [CrossRef]

15. Boglione, C.; Gisbert, E.; Gavaia, P.; Witten, P.E.; Moren, M.; Fontagné, S.; Koumoundouros, G. Skeletal anomalies in reared European fish larvae and juveniles. Part 2: Main typologies, occurrences and causative factors. *Rev. Aquac.* **2013**, *5*, S121-S167.
[CrossRef]

16. Eissa, A.E.; Abu-Seida, A.M.; Ismail, M.M.; Abu-Elala, N.M.; Abdelsalam, M. A comprehensive overview of the most common skeletal deformities in fish. *Aquac. Res.* **2021**, *52*, 2391–2402. [CrossRef]

17. Kourkouta, C.; Tsipourlianos, A.; Power, D.M.; Moutou, K.A.; Koumoundouros, G. Variability of key-performance-indicators in commercial gilthead seabream hatcheries. *Sci. Rep.* **2022**, *12*, 17896. [CrossRef]

18. Izquierdo, M.S.; Ghrab, W.; Roo, J.; Hamre, K.; Hernández-Cruz, C.M.; Bernardini, G.; Terova, G.; Saleh, R. Organic, inorganic and nanoparticles of Se, Zn and Mn in early weaning diets for gilthead seabream (*Sparus aurata*; Linnaeus, 1758). *Aquac. Res.* **2016**, *48*, 2852–2867. [CrossRef]

19. Izquierdo, M.; Domínguez, D.; Jiménez, J.I.; Saleh, R.; Hernández-Cruz, C.M.; Zamorano, M.J.; Hamre, K. Interaction between taurine, vitamin E and vitamin C in microdiets for gilthead seabream (*Sparus aurata*) larvae. *Aquaculture* **2018**, *498*, 246–253.
[CrossRef]

20. Dominguez, D.; Montero, D.; Zamorano, M.J.; Castro, P.; Fontanillas, R.; Prabhu, P.A.J.; Izquierdo, M. Effects of vitamin D3 supplementation in gilthead seabream (*Sparus aurata*) juveniles fed diets high in plant based feedstuffs. *Aquaculture* **2021**, *543*, 736991. [CrossRef]

21. Tseng, Y.; Dominguez, D.; Bravo, J.; Acosta, F.; Robaina, L.; Geraert, P.-A.; Kaushik, S.; Izquierdo, M. Organic Selenium (OH-MetSe) Effect on Whole Body Fatty Acids and *Mx* Gene Expression against Viral Infection in Gilthead Seabream (*Sparus aurata*) Juveniles. *Animals* **2021**, *11*, 2877. [CrossRef]

22. Sivagurunathan, U.; Dominguez, D.; Tseng, Y.; Eryalçın, K.M.; Roo, J.; Boglione, C.; Prabhu, P.A.J.; Izquierdo, M. Effects of dietary vitamin D3 levels on survival, mineralization, and skeletal development of gilthead seabream (*Sparus aurata*) larvae. *Aquaculture* **2022**, *560*, 738505. [CrossRef]

23. Dominguez, D.; Castro, P.; Lall, S.; Montero, D.; Zamorano, M.J.; Fontanillas, R.; Izquierdo, M. Effects of Menadione Sodium Bisulphite (Vitamin K3) Supplementation of the Diets Based on Plant Feed Ingredients on Growth and Bone Health of Gilthead Seabream (*Sparus aurata*) Fingerlings. *Aquac. Nutr.* **2022**, *2022*, 1613030. [CrossRef]

24. Ferosekhan, S.; Sarih, S.; Afonso, J.M.; Zamorano, M.J.; Fontanillas, R.; Izquierdo, M.; Kaushik, S.; Montero, D. Selection for high growth improves reproductive performance of gilthead seabream Sparus aurata under mass spawning conditions, regardless of the dietary lipid source. *Anim. Reprod. Sci.* **2022**, *241*, 106989. [CrossRef] [PubMed]

25. Berillis, P. Skeletal Deformities in Seabreams. Understanding the Genetic Origin Can Improve Production? *J. Fish. Sci.* **2017**, *11*, 57. [CrossRef]

26. Fragkoulis, S.; Batargias, C.; Kolios, P.; Koumoundouros, G. Genetic parameters of the upper-jaw abnormalities in Gilthead seabream *Sparus aurata*. *Aquaculture* **2018**, *497*, 226–233. [CrossRef]

27. Fragkoulis, S.; Economou, I.; Moukas, G.; Koumoundouros, G.; Batargias, C. Caudal fin abnormalities in Gilthead seabream (*Sparus aurata* L.) have a strong genetic variance component. *J. Fish Dis.* **2020**, *43*, 825–828. [CrossRef] [PubMed]

28. Balbuena-Pecino, S.; Riera-Heredia, N.; Vélez, E.J.; Gutiérrez, J.; Navarro, I.; Riera-Codina, M.; Capilla, E. Temperature Affects Musculoskeletal Development and Muscle Lipid Metabolism of Gilthead Sea Bream (*Sparus aurata*). *Front. Endocrinol.* **2019**, *10*, 173. [CrossRef]

29. Suniaga, S.; Rolvien, T.; Scheidt, A.V.; Fiedler, I.A.K.; Bale, H.A.; Huysseune, A.; Witten, P.E.; Amling, M.; Busse, B. Increased mechanical loading through controlled swimming exercise induces bone formation and mineralization in adult zebrafish. *Sci. Rep.* **2018**, *8*, 3646. [CrossRef]

30. Hall, B.K.; Witten, P.E. Plasticity and variation of skeletal cells and tissues and the evolutionary development of actinopterygian fishes. In *Evolution and Development of Fishes*; Johanson, Z., Underwood, C., Richter, M., Eds.; Cambridge University Press: Cambridge, UK; New York, NY, USA, 2018; pp. 126–143.

31. Boglione, C. Skeletal abnormalities. In *Climate Change and Non-Infectious Fish Disorders*; Woo, P.T.K., Iwama, G.K., Eds.; Cabi: Wallingford, UK; Boston, MA, USA, 2020; pp. 54–79.

32. Palstra, A.P.; Roque, A.; Kruijt, L.; Jéhannet, P.; Pérez-Sánchez, J.; Dirks, R.P. Physiological Effects of Water Flow Induced Swimming Exercise in Seabream *Sparus aurata*. *Front. Physiol.* **2020**, *11*, 610049. [CrossRef]
33. Araújo-Luna, R.; Ribeiro, L.; Bergheim, A.; Pousão-Ferreira, P. The impact of different rearing condition on gilthead seabream welfare: Dissolved oxygen levels and stocking densities. *Aquac. Res.* **2018**, *49*, 3845–3855. [CrossRef]

34. Arechavala-Lopez, P.; Nazzaro-Alvarez, J.; Jardí-Pons, A.; Reig, L.; Carella, F.; Carrassón, M.; Roque, A. Linking stocking densities and feeding strategies with social and individual stress responses on gilthead seabream (*Sparus aurata*). *Physiol. Behav.* **2020**, *213*, 112723. [CrossRef]

35. Canario, A.V.M.; Condeça, J.; Power, D.M.; Ingleton, P.M. The effect of stocking density on growth in the gilthead sea-bream, *Sparus aurata* (L.). *Aquac. Res.* **1998**, *29*, 177–181. [CrossRef]

36. Montero, D.; Blazer, V.; Socorro, J.; Izquierdo, M.; Tort, L. Dietary and culture influences on macrophage aggregate parameters in gilthead seabream (*Sparus aurata*) juveniles. *Aquaculture* **1999**, *179*, 523–534. [CrossRef]

37. Montero, D.; Robaina, L.; Socorro, J.; Vergara, J.; Tort, L.; Izquierdo, M. Alteration of liver and muscle fatty acid composition in gilthead seabream (*Sparus aurata*) juveniles held at high stocking density and fed an essential fatty acid deficient diet. *Fish Physiol. Biochem.* **2001**, *24*, 63–72. [CrossRef]

38. Tort, L.; Sunyer, J.; Gómez, E.; Molinero, A. Crowding stress induces changes in serum haemolytic and agglutinating activity in the gilthead sea bream *Sparus aurata*. *Vet. Immunol. Immunopathol.* **1996**, *51*, 179–188. [CrossRef]

39. Prestinicola, L.; Boglione, C.; Makridis, P.; Spanò, A.; Rimatori, V.; Palamara, E.; Scardi, M.; Cataudella, S. Environmental Conditioning of Skeletal Anomalies Typology and Frequency in Gilthead Seabream (*Sparus aurata* L., 1758) Juveniles. *PLoS ONE* **2013**, *8*, e55736. [CrossRef]

40. Boglione, C.; Gagliardi, F.; Scardi, M.; Cataudella, S. Skeletal descriptors and quality assessment in larvae and post-larvae of wild-caught and hatchery-reared gilthead sea bream (*Sparus aurata* L. 1758). *Aquaculture* **2001**, *192*, 1–22. [CrossRef]

41. Divanach, P.; Kentouri, M. Hatchery techniques for specific diversification in Mediterranean finfish larviculture. *Cah. Opt. Medit.* **2000**, *47*, 75–87.

42. Martini, A.; Huysseune, A.; Witten, P.E.; Boglione, C. Plasticity of the skeleton and skeletal deformities in zebrafish (*Danio rerio*) linked to rearing density. *J. Fish Biol.* **2020**, *98*, 971–986. [CrossRef]

43. Di Biagio, C.; Dellacqua, Z.; Martini, A.; Huysseune, A.; Scardi, M.; Witten, P.E.; Boglione, C. A Baseline for Skeletal Investigations in Medaka (*Oryzias latipes*): The Effects of Rearing Density on the Postcranial Phenotype. *Front. Endocrinol.* **2022**, *13*. [CrossRef]

44. Dellacqua, Z.; Di Biagio, C.; Martini, A.; Mattei, F.; Rakaj, A.; Williams, J.; Fabris, A.; Izquierdo, M.; Boglione, C. Distinguishing the Impacts of Rearing Density versus Tank Volume on the Skeletal Quality and Development of Gilthead Seabream (*Sparus aurata*) during the Hatchery Phase. *Animals* **2023**, submitted.

45. Meena, D.; Sahoo, A.; Jayant, M.; Sahu, N.; Srivastava, P.; Swain, H.; Behera, B.; Satvik, K.; Das, B. Bioconversion of Terminalia arjuna bark powder into a herbal feed for Labeo rohita: Can it be a sustainability paradigm for Green Fish production? *Anim. Feed. Sci. Technol.* **2021**, *284*, 115132. [CrossRef]

46. Witten, P.E.; Hansen, A.; Hall, B.K. 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. [CrossRef] [PubMed]

47. Schindelin, J.; Arganda-Carreras, I.; Frise, E.; Kaynig, V.; Longair, M.; Pietzsch, T.; Preibisch, S.; Rueden, C.; Saalfeld, S.; Schmid, B.; et al. Fiji: An open-source platform for biological-image analysis. *Nat. Methods* **2012**, *9*, 676–682. [CrossRef] [PubMed]

48. Pearson, K.X. On the criterion that a given system of deviations from the probable in the case of a correlated system of variables is such that it can be reasonably supposed to have arisen from random sampling. *Lond. Edinb. Dublin Philos. Mag.* **1900**, *50*, 157–175. [CrossRef]

49. Sfakianakis, D.; Georgakopoulou, E.; Papadakis, I.; Divanach, P.; Kentouri, M.; Koumoundouros, G. Environmental determinants of haemal lordosis in European sea bass, *Dicentrarchus labrax* (Linnaeus, 1758). *Aquaculture* **2006**, *254*, 54–64. [CrossRef]

50. Hammer, Ø.; Harper, D.A.; Ryan, P.D. PAST: Paleontological statistics software package for education and data analysis. *Palaeontol. Electron.* **2001**, *4*, 9.

51. Loy, B.A.; Boglione, C.; Cataudella, S. Geometric morphometrics and morpho-anatomy: A combined tool in the study of sea bream (*Sparus aurata*, Sparidae) shape. *J. Appl. Ichthyol.* **1999**, *15*, 104–110. [CrossRef]

52. Rohlf, F.J.; Slice, D. Extensions of the Procrustes Method for the Optimal Superimposition of Landmarks. *Syst. Zool.* **1990**, *39*, 40–59. [CrossRef]

53. Costa, C.; Tibaldi, E.; Pasqualetto, L.; Loy, A. Morphometric comparison of the cephalic region of cultured *Acipenser baerii* (Brandt, 1869), *Acipenser naccarii* (Bonaparte, 1836) and their hybrid*. *J. Appl. Ichthyol.* **2006**, *22*, 8–14. [CrossRef]

54. Legendre, P.; Legendre, L. *Numerical Ecology*, 2nd ed.; Elsevier Sci.: Amsterdam, The Netherlands; 1998; 853p.

55. Bookstein, F. Principal warps: Thin-plate splines and the decomposition of deformations. *IEEE Trans. Pattern Anal. Mach. Intell.* **1989**, *11*, 567–585. [CrossRef]

56. Loy, A.; Mariani, L.; Bertelletti, M.; Tunesi, L. Visualizing Allometry: Geometric Morphometrics in the Study of Shape Changes in the Early Stages of the Two-Banded Sea Bream, *Diplodus vulgaris* (*Perciformes, Sparidae*). *J. Morphol.* **1998**, *237*, 137–146. [CrossRef]

57. Colangelo, P.; Ventura, D.; Piras, P.; Bonaiuti, J.P.G.; Ardizzone, G. Are developmental shifts the main driver of phenotypic evolution in *Diplodus* spp. (*Perciformes: Sparidae*)? *BMC Evol. Biol.* **2019**, *19*, 106. [CrossRef]

58. Kerschbaumer, M.; Sturmbauer, C. The Utility of Geometric Morphometrics to Elucidate Pathways of Cichlid Fish Evolution. *Int. J. Evol. Biol.* **2011**, *2011*, 290245. [CrossRef]

59. Costa, C.; Pasqualetto, L.; Tibaldi, E.; Loy, A. Shape changes and differences in the cranial region in *Acipenser baerii* and *A. naccarii* of different size. *Ital. J. Zool.* **2004**, *71*, 57–61. [CrossRef]
60. Loy, A.; Boglione, C.; Gagliardi, F.; Ferrucci, L.; Cataudella, S. Geometric morphometrics and internal anatomy in sea bass shape analysis (Dicentrarchus labrax L., Moronidae). *Aquaculture* **2000**, *186*, 33–44. [CrossRef]

61. Yiannis, G.R.; Katerina, K.K.; Alexios, R.; Evagelos, D.; George, N.K. Comparative morphology of wild, farmed and hatchery-released gilthead sea bream (*Sparus aurata*) in western Greece. *Int. J. Fish. Aquac.* **2011**, *3*, 1–9.

62. Russo, T.; Costa, C.; Cataudella, S. Correspondence between shape and feeding habit changes throughout ontogeny of gilthead sea bream *Sparus aurata* L., 1758. *J. Fish Biol.* **2007**, *71*, 629–656. [CrossRef]

63. Koumoundouros, G. Morpho-anatomical abnormalities in Mediterranean marine aquaculture. In *Recent Advances in Aquaculture Research*; Transworld Research Network: Trivandrum, India, 2010; Volume 661, pp. 125–148.

64. Fragkoulis, S.; Koumoundouros, G. Simple morphometrics for predicting lordosis-induced deviations of body shape in reared Gilthead seabream (*Sparus aurata* L.). *J. Fish Dis.* **2021**, *44*, 1265–1267. [CrossRef]

65. Printzi, A.; Fragkoulis, S.; Dimitriadi, A.; Keklikoglou, K.; Arvanitidis, C.; Witten, P.E.; Koumoundouros, G. Exercise-induced lordosis in zebrafish *Danio rerio* (Hamilton, 1822). *J. Fish Biol.* **2021**, *98*, 987–994. [CrossRef]

66. Fragkoulis, S.; Printzi, A.; Geladakis, G.; Katribouzas, N.; Koumoundouros, G. Recovery of haemal lordosis in Gilthead seabream (*Sparus aurata* L.). *Sci. Rep.* **2019**, *9*, 9832. [CrossRef]

67. Fragkoulis, S.; Kourkouta, C.; Geladakis, G.; Printzi, A.; Glaropoulos, A.; Koumoundouros, G. Recovery of Haemal Lordosis in European Seabass *Dicentrarchus labrax* (Linnaeus 1758). *Aquac. J.* **2022**, *2*, 1–12. [CrossRef]

68. Chatain, B.; Dewavrin, G. Influence des anomalies de développement de la vesse natatoire sur la mortalité de Dicentrarchus labrax au cours du sevrage. *Aquaculture* **1989**, *78*, 55–61. [CrossRef]

69. Chatain, B. Abnormal swimbladder development and lordosis in sea bass (*Dicentrarchus labrax*) and sea bream (*Sparus auratus*). *Aquaculture* **1994**, *119*, 371–379. [CrossRef]

70. Kitajima, C.; Watanabe, T.; Tsukashima, Y.; Fujita, S. Lordotic Deformation and Abnormal Development of Swim Bladders in Some Hatchery-Bred Marine Physoclistous Fish in Japan. *J. World Aquac. Soc.* **1994**, *25*, 64–77. [CrossRef]

71. Boglione, C.; Marino, G.; Fusari, A.; Ferreri, A.; Finoia, M.G.; Cataudella, S. Skeletal anomalies in *Dicentrarchus labrax* juveniles selected for functional swimbladder. *ICES Mar. Sci. Symp.* **1995**, *201*, 163–169.

72. Andrades, J.; Becerra, J.; Fernández-Llebrez, P. Skeletal deformities in larval, juvenile and adult stages of cultured gilthead sea bream (*Sparus aurata* L.). *Aquaculture* **1996**, *141*, 1–11. [CrossRef]

73. Divanach, P.; Papandroulakis, N.; Anastasiadis, P.; Koumoundouros, G.; Kentouri, M. Effect of water currents on the development of skeletal deformities in sea bass (*Dicentrarchus labrax* L.) with functional swimbladder during postlarval and nursery phase. *Aquaculture* **1997**, *156*, 145–155. [CrossRef]

74. Kihara, M.; Ogata, S.; Kawano, N.; Kubota, I.; Yamaguchi, R. Lordosis induction in juvenile red sea bream, Pagrus major, by high swimming activity. *Aquaculture* **2002**, *212*, 149–158. [CrossRef]

75. Kranenbarg, S.; Waarsing, J.H.; Muller, M.; Weinans, H.; van Leeuwen, J.L. Lordotic vertebrae in sea bass (*Dicentrarchus labrax* L.) are adapted to increased loads. *J. Biomech.* **2005**, *38*, 1239–1246. [CrossRef]

76. Boursiaki, V.; Theochari, C.; Zaoutsos, S.P.; Mente, E.; Vafidis, D.; Apostologamvrou, C.; Berillis, P. Skeletal Deformity of Scoliosis in Gilthead Seabreams (*Sparus aurata*): Association with Changes to Calcium-Phosphor Hydroxyapatite Salts and Collagen Fibers. *Water* **2019**, *11*, 257. [CrossRef]

77. Blasco, J.; Moya, A.; Millán-Cubillo, A.; Vélez, E.J.; Capilla, E.; Pérez-Sánchez, J.; Gutiérrez, J.; Borrás, J.F. Growth-promoting effects of sustained swimming in fingerlings of gilthead sea bream (*Sparus aurata* L.). *J. Comp. Physiol. B* **2015**, *185*, 859–868. [CrossRef]

78. Moya, A.; Torrella, J.; Fernández-Borràs, J.; Rizo-Roca, D.; Millán-Cubillo, A.; Vélez, E.; Arcas, A.; Gutiérrez, J.; Blasco, J. Sustained swimming enhances white muscle capillarisation and growth by hyperplasia in gilthead sea bream (*Sparus aurata*) fingerlings. *Aquaculture* **2019**, *501*, 397–403. [CrossRef]

79. Ibarz, A.; Felip, O.; Fernández-Borràs, J.; Martín-Pérez, M.; Blasco, J.; Torrella, J.R. Sustained swimming improves muscle growth and cellularity in gilthead sea bream. *J. Comp. Physiol. B* **2011**, *181*, 209–217. [CrossRef]

80. Sánchez-Muros, M.J.; Sánchez, B.; Barroso, F.G.; Toniolo, M.; Trenzado, C.E.; Rus, A.S. Effects of rearing conditions on behavioural responses, social kinetics and physiological parameters in gilthead sea bream *Sparus aurata*. *Appl. Anim. Behav. Sci.* **2017**, *197*, 120–128. [CrossRef]

81. Beamish, F.W.H. Swimming Capacity. *Fish Physiol.* **1978**, *7*, 101–187.

82. Gui, F.; Wang, P.; Wu, C. Evaluation approaches of fish swimming performance. *Agric. Sci.* **2014**, *05*, 106–113. [CrossRef]

83. Kozhara, A.V. Phenotypic Variance of Bilateral Characters as an Indicator of Genetic and Environmental Conditions in Bream *Abramis brama* (L.) (Pisces, Cyprinidae) Populations. *J. Appl. Ichthyol.* **1994**, *10*, 167–181. [CrossRef]