Abstract

Phenotypic plasticity has been demonstrated in fishes but rarely addressed with respect to skeletogenesis. The influence of water velocity on the sequence of chondrification and ossification is studied for the median fins of Arctic charr, Salvelinus alpinus, during a period of 90 days post hatching. Time of appearance, relative position within sequences, and direction of development among serially repeated elements are compared between two velocity treatments. Water velocity has induced changes in the timing of events and to a lesser extent on the relative sequence events in the locomotor system. Ossification is more responsive to water velocity than chondrification, and early-forming elements are less responding than late-forming elements. Directions of development are fairly conservative. It is suggested that a faster sustained swimming (behavioural adaptation to a higher water velocity) could induce differential mechanical stresses on developing skeletal elements involved in locomotion and therefore induce changes primarily in the timing of the ossification.

1. Introduction

Developmental sequences are known to be controlled genetically as well as constrained environmentally [1, 2]. A specific type of developmental sequence, the sequences of chondrification and ossification, has been investigated in a large diversity of fishes [313], amphibians [1419], reptiles (including birds) [2028], and mammals [2932]; both sequences refer to the specific ontogenetic order in which anatomical cartilaginous and bony structures appear during the early development of vertebrates. Information derived from such sequences could be used to validate homology [12, 13, 33] and phylogenetic position [14, 20], and to infer evolutionary developmental patterns and processes (e.g., heterochrony) [11, 29, 30, 32, 34, 35].

It has been suggested that ossification sequences conform to functional needs [2931, 35, 36]. In a general way, structures that are required functionally earlier during development will ossify earlier in the sequence although inconsistencies may occur [3, 37, 38]. Among the functional requirements most influential on the ossification, locomotion has been shown to induce differential mechanical stress, which may change the shape, size, or toughness of the bones [3942]. In order to understand if the observed patterns of ossification are the result of processes controlled genetically or developmentally, as well as a response to environmental constraints, experimental studies on living fish morphology are of primary interest.

The effects of many environmental parameters (e.g., pH, temperature, salinity, and water viscosity) have been tested experimentally on the early development of larval and juvenile fishes from hatcheries and wild stocks. It has long been known that experimental changes in temperature, photoperiod and salinity affect the meristic count of serial (or iterated) bones such as fin rays and vertebrae [43, 44]. However, studies that have shown an experimental effect on the timing of ossification have mainly focused on temperature [35, 4548]. Pavlov and Moksness [46] demonstrated that wolffish Anarhichas lupus incubated at higher temperature exhibited morphological modifications of the axial and appendicular skeleton (e.g., fused hypurals, fused epurals, and absence of distal radials and fin rays in median fins). Mabee et al. [35] also tested the effect of temperature during incubation and rearing periods of the zebrafish Danio rerio and the Siamese fighting fish Betta splendens. They found out that there was no significant temperature effect on the sequence of ossification, or only minor variation in the sequence. However, they observed that fishes incubated and reared at higher temperature were characterised with higher levels of intraspecific variation in iterated cranial bones (e.g., infraorbitals, branchial arch bones, and branchiostegal rays) and in cartilage bones (preformed in a cartilaginous template sensu Patterson [49]). Campinho et al. [48] studied the sensitivity of chondrification and ossification to temperature in the Mozambique tilapia Oreochromis mossambicus. They found that the sequence of ossification and rate at which it occurred is more sensitive to temperature than is chondrogenesis. Sequences of chondrification and ossification in O. mossambicus were not altered by temperature. All previous studies investigated different responses of the ossification (primarily) to a temperature treatment during the early development. However, no study evaluating the effect of water velocity on the early chondrification and ossification of fishes has been published.

According to Davison’s [50] review on the effect of exercise training on fish, most of the literature deals with food conversion ratio, social behaviour, growth rate, and physiology. Azuma [51] and Azuma et al. [52] have demonstrated that juvenile masu salmon Oncorhynchus masou masou reared in higher flow speed had a better growth profile after c. 100 days. Furthermore, the effects of water velocity (in natural and laboratory conditions) on the plasticity of salmonid body size and shape have been investigated in a few salmonid species [36, 5359]. All of these studies have shown significant effect of water velocity on size and differentially on shape. Water velocity is one of the most limiting environmental factors regulating the distribution and viability of larval fishes. Different types of larval behaviour (limited to high mobility) are associated with different environments subjected to flow variation [60]. Developmental patterns are known to reflect adaptive escape behaviour in fishes in order to optimize survival [61] and are thus hypothesized to be maintained in the course of evolution. Changes of velocity are hypothesized to modify mechanical stress imposed on the locomotor system, pointing to the usefulness of studying the development of the postcranial system under differential water velocities in the laboratory.

In the present study, we test the influence of water velocity during rearing in the Arctic charr, Salvelinus alpinus. The effect of water velocity on the sequence of chondrification and ossification (i.e., timing, order, and direction of formation among elements) is investigated for the dorsal, anal, and caudal fins, which are three major morphological systems associated with fish locomotion. Patterns of chondrification and ossification will be analyzed in terms of absolute and relative timing, rank shift variation, and direction of formation (i.e., uni- or bidirectionally) among serially repeated elements (e.g., radials, hypurals, and caudal vertebrae).

2. Materials and Methods

Fertilized Arctic charr, Salvelinus alpinus, eggs used for this experiment originated from the Centre Aquacole Marin de Grande-Rivière from their captive Nauyuk strain broodstock popolation (1 female × 1 male).

2.1. Rearing Conditions

Incubation temperature was fixed at 5°-6°C (with 85% survival during incubation). At hatching, fry were initially kept in the incubation trays for part of the yolk absorption period, and, afterwards, 500 fish were transferred to the rearing units where the water temperature was raised gradually from 6°C to 9°C (0.5°C per day) and kept constant for the duration of the experiment (90 days). When the absorption phase was completed, first feeding was initiated using enriched artemia and formulated feed to ensure proper food intake of both experimental treatments despite different water flow through the experimental units. Enriched artemia (500 artemia L−1 min−1) supplementation was stopped after 14 days, at which time inert food ingestion was generally acquired. Fish were fed manually to satiety four times a day. The experimental units were low volume raceways (20 cm wide × 60 cm long × 15 cm high) with individual oxygen supply. Water inflow was dechlorinated and UV treated. Mortality, temperature, and dissolved oxygen were monitored twice a day. Overall mortality was low at less than 1%. Photoperiod condition followed a constant cycle 12 : 12 (L : D).

2.2. Experimental Design

The original design of the experiment took into account water flow rate rather than water velocity. The water flow, kept constant throughout the experiment, was fixed at 2 L min−1 in treatment S (slow) and at 6 L min−1 in treatment F (fast); the optimal rearing condition of S. alpinus being 6 L min−1 [62]. Velocities were calculated with a flowmeter SonTek ATV. Three depth measurements were taken for each treatment: (1) half the water depth in the centre of the unit, (2) half the water depth where the alevins used to gather, and (3) at the location where the alevins used to gather. The average velocity for the current in treatment S (2 L min−1) was 0.37733 cm s−1, whereas for treatment F (6 L min−1) it was 0.50629 cm s−1.

2.3. Specimens Examined

For each treatment, two individuals were sampled every other day from 2 to 90 days post hatching (dph). The fixation of alevins follows Mabee's [63] protocol. A total of 33 individuals per treatment ranging from 2 dph [c. 13 mm total length ( 𝐿 𝑇 )] to 90 dph (c. 28 mm 𝐿 𝑇 ) were cleared and double stained for cartilage and bone [64]. The total length ( 𝐿 𝑇 ) and standard length ( 𝐿 𝑆 ) were measured with an ocular micrometer.

2.4. Anatomical Nomenclature

Nomenclature for the appendicular and axial skeletons follows that of Arratia and Schultze [44], Arratia et al. [65], and Grünbaum and Cloutier [33]. The caudal fin involves all the elements posterior to neural spine of the preural centrum 2 (i.e., corresponding to the neural spine associated with the second preural caudal vertebra) inclusively [44]. The dorso-caudal lepidotrichia are located dorsally to the notochord articulating with the epurals, in contrast to the ventro-caudal lepidotrichia located ventrally to the notochord and articulating with the haemal spines, parhypural, and hypurals.

2.5. Determination of Developmental Sequences

Individuals were examined and scored under a stereodissecting Wild MZ9.5 microscope. A structure was considered to be forming or ossified when it uptakes the stain (Alcian blue for cartilage and Alizarin red S for bone). Chondrification and ossification sequences are defined in terms of successive events. A chondrification event corresponds to the appearance of a single cartilaginous element. An ossification event corresponds to the ossification of one cartilaginous precursor as well as the direct ossification of an element.

Patterns of chondrification and ossification for 70 endoskeletal elements including the proximal and distal radials of the dorsal and anal fins, and the hypurals, epurals, uroneurals, and vertebrae of the caudal fin are determined plus the ossification patterns of the fin rays of all fins. Chondrification and ossification sequences were described in terms of three major components: (1) the absolute timing (in terms of dph) which provides the temporal sequence of appearance of the elements (cartilages and bones), (2) the temporal information converted to relative timing, or the sequential order of formation of these structures, and (3) when dealing with serially repeated elements, direction of chondrification or ossification (i.e., antero-posteriorly, postero-anteriorly, and bidirectional with respect to the body axis) is inferred. Sequences of chondrification and ossification are compared between treatments in terms of (1) the time of first appearance of an element, (2) the relative timing within morphological units (i.e., dorsal, anal, and caudal fins defined by their anatomical boundaries and respective elements), and (3) the direction of chondrification and ossification among serial elements within morphological units.

2.6. Statistical Analyses

Simple linear regressions were used to describe the relationship between log-transformed length (Log10 𝐿 𝑇 and Log10 𝐿 𝑆 ) and age (dph). Analysis of covariance (ANCOVA) was used to compare growth between treatments.

Water velocity treatment might have an effect on the timing of a developmental event and on the order of an event in a sequence. Changes in the order might also imply changes in the direction of formation of serial elements. Differences in the timing of first appearance (in dph) between treatments were calculated for each skeletal element composing the three median fins. Parametric one-sample 𝑡 -tests were used to compare the mean difference in absolute timing of appearance of skeletal element in each fin between treatments with a theoretical value of zero (no treatment effect). The relative timing of events was investigated by attributing a rank value to each event in each sequence. For instance, the five hypothetical elements A–E are forming in the following sequence of events [E – B – (A and C) – D]: E appears first, followed by B, then A and C form simultaneously and finally D is the last element to form (i.e., D is the last event of the sequence). The sequence is then composed of five events with two events being simultaneous. The rank value of each event is then scored as follows: E is the first element to form defining an event with a rank value of 1. Element B is the second event in the sequence; its rank value is then 2. As elements A and C have simultaneously formed, their order as events within the sequence cannot be discriminate. Indeed, A and C are the 3rd or 4th and vice versa. Then, both events are each ranked as 3.5 and finally the last element D is the 5th event within the sequence and is consequently ranked as 5. Equal rank was readjusted according to the procedure of Dagnelie [66]. The almost synchronous appearance of elements might also be an artifact of coarse sampling. Spearman rank correlation coefficients are used in order to compare the concordance between sequences [7]. For the purpose of chondrification and ossification sequences, rank comparisons between only those elements with a cartilaginous precursor and a corresponding ossified bone were integrated into the analysis. Statistical analyses were performed using SYSTAT version 10 and SAS 9.1.3.

3. Results

3.1. Growth

A linear model best described the relationship between the age and L o g 1 0 𝐿 𝑇 ( 𝑛 = 2 6 ; Treatment S 𝑅 2 = 0 . 7 9 0 5 ; Treatment F 𝑅 2 = 0 . 9 0 9 0 ) and L o g 1 0 𝐿 𝑆 ( 𝑛 = 2 6 ; Treatment S 𝑅 2 = 0 . 7 9 4 4 ; Treatment F 𝑅 2 = 0 . 8 8 7 7 ) for the 90 dph period. Regression lines do not differ significantly between treatments (ANCOVA: Log10 𝐿 𝑇 : 𝐹 1 , 5 1 = 1 . 1 2 ; 𝑃 = . 2 9 5 8 ; L o g 1 0 𝐿 𝑆 𝐹 1 , 5 1 = 0 . 2 8 ; 𝑃 = . 6 0 1 4 ).

3.2. Median Fin Anatomy

The dorsal fin is composed of 13 proximal radials and 12 distal radials (Figures 1(a) and 1(d)); occasionally a small 14th proximal radial is formed. The first proximal radial is not associated with a distal radial; rarely, the second proximal radial lacks a distal radial. The first and second proximal radials are short and relatively straight. The distal end of proximal radials 4–13 is curved posteriorly. The first and second proximal radials are shorter than the remaining radials. The distal radials are small and grossly rounded. In some individuals, proximal radials 1 and 2, and proximal radials 12 and 13 show fusion of their proximal parts (i.e., innermost part within the body). The dorsal fin includes c. 14–16 lepidotrichia.

The general anatomy of the anal fin is similar to that of the dorsal fin. The supports of the anal fin include 12-13 proximal radials and 11-12 distal radials (Figures 1(b) and 1(d)). In some individuals, proximal radial 6 and 12 are poorly developed. The first and second proximal radials are not associated with a distal radial. Similarly to the distal radials of the dorsal fin, the anal distal radials are small and rounded. Proximal radial 12 is occasionally associated with two or three distal radials. Fusions of proximal parts are relatively frequent between proximal radials 1 and 2, 8 and 9, 9 and 10, 10 and 11, and 11 and 12. Approximately 13–15 lepidotrichia are associated with the anal fin.

The anatomy of the caudal fin includes, dorsal to the notochord, the neural arches and spines of preural centra 1 and 2, two epurals, two uroneurals, and the stegural (Figures 1(c) and 1(d)). Ventrally, the elements are the haemal arches and spines, the parhypural (= hypural 1 of Vladykov [67]), and a series of six hypurals (= hypurals 2–7 of Vladykov [67]). Finally, surrounding the notochord, we coded for one preural centrum (i.e., preural centrum 1, PU1) and two ural chordacentra (chcU2, 4–6). As mentioned by Arratia and Schultze [44], for Salvelinus spp., ural centrum 2 articulates mainly with hypural 2 whereas ural centrum 4 articulates with hypural 4. Ural chordacentra 5 and 6 are not present in all individuals but are nevertheless included in the sequence of ossification. When present ural chordacentra 5 and 6 are poorly developed and articulate with their respective hypurals. Two sets of lepidotrichia compose the caudal fin: 5–7 dorso-caudal lepidotrichia and 26-27 ventro-caudal lepidotrichia.

3.3. Time of First Appearance for an Element

The time of first chondrification and ossification is provided for each element of the dorsal and anal fins (Figure 2) and the caudal fin (Figure 3) with the exception of their fin rays for both treatments during the 90 days of experiment. The general pattern of skeletogenesis for the dorsal and anal fins is similar (Figure 2). The chondrification of the dorsal and anal fins is primarily taken into account during the time of the experiment whereas it is primarily the ossification pattern of the caudal fin that occurs during the same period (Figure 3).

Cartilaginous elements are present in the three median fins as early as 2 dph. With respect to the caudal and anal fins, more cartilaginous elements are present at day 2 in treatment S than in treatment F. All cartilaginous precursors form between 2 and 86 dph in treatment S whereas they form between 2 and 80 dph in treatment F. The completion of the cartilaginous skeleton is achieved earlier for the three fins in treatment F.

Most elements ossified between 26 and 90 dph, with the exception of the dorsal distal radials and anal distal radials that ossify later than 90 days (Figure 2). The caudal fin is the first of the three fins to ossify some elements (e.g., 26 dph, haemal spine 2, hypurals 1–4 and parhypural) followed by the dorsal (58 and 52 dph for treatments S and F, resp.) and anal fin (62 and 54 dph for treatments S and F, resp.) (Figures 2 and 3).

Figure 4 compares the time of first appearance between the two treatments for all cartilaginous ( 𝑛 = 6 5 ; Figure 4(a)) and ossified elements ( 𝑛 = 4 5 ; Figure 4(b)) of the three median fins simultaneously. Elements figuring above the synchrony line are those forming in treatment S before treatment F (21 cartilaginous and 6 ossified elements), the elements on the line are those appearing simultaneously in both treatments (30 cartilaginous and 10 ossified elements), whereas the elements below the line form earlier in treatment F than in treatment S (14 cartilaginous and 28 ossified elements). Independently of the three median fins, the timing of chondrification for c. 46% of the elements remains unchanged in both treatments whereas c. 64% of the elements ossify earlier in treatment F (Figures 24). Certain bones, such as neural spine 2, display a temporal shift of 16 days later in treatment F (this shift is even more pronounced for this element in chondrification, 24 days, Figures 3 and 4). The sequence of chondrification ( 𝑛 = 6 5 ) and ossification ( 𝑛 = 4 4 ) combining all fins elements (only the elements for which a specific day of appearance is known) are fairly congruent between treatments (chondrification: 𝑟 𝑠 = 0 . 8 9 5 2 2 , 𝑃 < . 0 0 0 1 ; ossification: 𝑟 𝑠 = 0 . 9 2 5 0 3 , 𝑃 < . 0 0 0 1 ), revealing similar timing of first appearance of skeletal elements between treatments.

On average, the cartilaginous elements of the caudal and anal fins form approximately at the same time with respect to both treatments (0.67 day and 0.27 day later in treatment F; Table 1), whereas dorsal fin elements appear on average 2.00 days earlier in treatment F. Although not significant, there is a weak tendency for the chondrification of the caudal and anal fins to be slower in treatment F (Table 1). The treatment effect is highly significant for the chondrification of the dorsal fin (Table 1). The chondrification of the caudal and anal fins displays greater variability than the dorsal fin whether it is expressed in terms of magnitude of difference inter-treatment or standard deviation (Table 1).

On average, caudal fin elements ossify 1.43 days earlier in treatment F; this difference is not significant (Figure 3; Table 1). In treatment F, the elements of the dorsal fin ossify on average 4.83 days earlier whereas for the anal fin they ossify on average 5.4 days earlier. The treatment effect is highly significant for the ossification of the anal and dorsal fins (Table 1).

Between day 20 and 90, the formation of each fin ray (lepidotrichium) includes two distinct phases (Figure 5): (1) the appearance of the fin rays support, the actinotrichia, that are composed of elastoidin fibers that are laid down by the ectoderm [68] and (2) the dermal ossification of the lepidotrichia. The number of fin rays for each fin increases slowly during the 90 days period. The number of fin rays per fin in treatment S is more variable than in treatment F; however, the low sample size might have introduced this variation. The number of fin rays per fin does not differ between treatments; it is slightly higher than the meristic counts provided by Vladykov [67] for adult S. alpinus. The lepidotrichia of the ventral part of the caudal fin are the first to ossify in both treatments, followed by a simultaneous ossification of the lepidotrichia in the dorsal and anal fins as well as the dorso-caudal lepidotrichia. The lepidotrichia ossify earlier in treatment F than in treatment S: 12 days earlier for the ventro-caudal rays and 10 days earlier for the dorsal and anal fins as well as the dorso-caudal rays.

3.4. Relative Timing within Morphological Systems

In the following section, the effect of treatment is evaluated in terms of magnitude of rank shift for element formation or ossification (i.e., relative order of event) rather than the difference in terms of the number of days. Rank shifts and associated statistics are represented for the chondrification and ossification of the three fins in Figure 6. Out of 65 elements, 27 elements (41.5%) show an acceleration of chondrification in the treatment S, 15 elements (23.1%) do not change their rank in the sequence, and 23 elements (35.4%) show acceleration in treatment F. With respect to the ossification of the 44 elements, 20 elements (45.4%) show an acceleration of ossification in treatment S, eight elements (18.2%) do not change their rank between treatments, and 16 elements (36.4%) show an acceleration of ossification in treatment F. Globally, the relative proportions of rank shift are similar during chondrification and ossification. In contrast to the absolute changes, more changes occur earlier in treatment S than treatment F in terms of rank shift.

During chondrification, the dorsal fin has the lowest variation of rank shift (SE = 0.4000). In terms of chondrification, 22% of the 18 caudal elements form earlier in treatment S than in treatment F. Neural spine 2 shows the greatest rank shifts in terms of both chondrification and ossification. For the dorsal fin, approximately half of the 25 elements do not change ranks between treatments; 14 elements display rank shifts with seven elements shifting ranks either to occur faster in favour of treatments S or F. In the dorsal fin, Drp13 is the cartilaginous element displaying the greatest magnitude of rank shift (6) forming earlier in treatment F. The chondrification of the anal fin is characterized by the greatest proportion of elements forming in treatment S before treatment F (73%), the greatest magnitude (10.5 ranks for Arp12), and the greatest variation (SE = 0.7458). For the dorsal and anal fins, the greatest magnitudes are associated with some of the last elements to chondrify in their respective fins and they are located at the posterior anatomical limit of these fins.

During ossification, the caudal fin has the lowest variation of rank shift (SE = 0.3686). Approximately 36% of the caudal elements (8) occupied the same rank in both treatments. Nine elements shifted earlier in the sequence in treatment F whereas five shifted later. Most elements of the dorsal fin show very low rank shift with the exception of Drp11. However, 75% of the elements ossify earlier in treatment S; this high proportion is associated with the relative importance of rank shift for Drp11. In contrast to the chondrification, proportions of rank shift are similar for the ossification of the anal and dorsal fins.

The general patterns of chondrification and ossification are congruent between treatments as can be validated by the Spearman rank correlation coefficients (Table 2). Coefficients (chondrification in treatment S versus chondrification in treatment F, and ossification in treatment S versus ossification in treatment F) are highly significant, meaning that the patterns are congruent independently of the treatment. Most correlations between chondrification and ossification of different treatments are also significant but with lower rank correlation and lower level of significance; non-significant patterns observed for the anal and dorsal fins might be a result of the low number of ossified elements compared ( 𝑛 = 1 0 (anal) and 12 (dorsal)). The pattern of relative timing is more conservative in the caudal fin. The relative timing of ossification among the bones within each one of the fins is fairly consistent in the two treatments. Faster water velocity has induced minor inversions of events in the sequences for each fin (e.g., anterior and posterior proximal radials of the dorsal and anal fins). However, such inversions do not affect radically the cohesion of the developmental sequence of skeletogenesis.

3.5. Direction of Chondrification and Ossification

In the caudal fin, directions are provided for the two epurals, two uroneurals, six hypurals, and the four ural chordocentra whereas numerous proximal and distal radials are used for the dorsal and anal fins (Table 3). With respect to the seven elements for which direction of chondrification is available, two of them have the potential to have been affected by the treatment whereas only one out of eight has been clearly affected with regard to the ossification (Table 3).

In the caudal fin, the chondrification and ossification of the epurals, uroneurals and hypurals in both treatments proceed primarily from anterior to posterior. The minor shift in the timing of chondrification and ossification observed in epurals suggests that both epurals could form simultaneously; the simultaneity results from opposing effect of water velocity on the two forming epurals (see Figure 3) or the low number of elements from which a direction is determined. Because hypurals 1–4 form as a block very early in ontogeny (or separate prehatching stages not sampled for this study), directionality is provided by the development of hypurals 5 and 6. The direction of ossification of the ural chordacentra is problematical owing to the sporadic presence of chondacentra 5 and 6 among the observed specimens. In both treatments, ural chordacentrum 4 ossifies before chordacentrum 2. This suggests a direction of formation from posterior to anterior. However, the ossification of ural chordacentra 5 and 6 precedes the formation of ural chordacentrum 2 in treatment F whereas it is formed after ural chordacentrum 2 in treatment S. Thus, a bidirectional formation is observed in both treatments with a variation in the alternation of the relative order of ural chordacentra 2 and 5-6 with reference to ural chordacentrum 4.

Although most proximal radials of the dorsal fin are present at 2 dph, chondrification proceeds bidirectionally as indicated by the late appearance of proximal radials 1 and 13 (Figure 2). In both treatments, the ossification of the dorsal fin is bidirectional. This bidirectionality is clearer in treatment S whereas a slight tendency to shift posteriorly is observed in treatment F. A bidirectional formation for the distal radials of the dorsal fin is interpreted solely because of the chondrification of distal radial 2; however, the general tendency is from anterior to posterior. The direction of ossification for the distal radials remains unknown because they were still cartilaginous at 90 dph in both treatments.

The proximal radials of the anal fin show indication of bidirectional chondrification in treatment S although only proximal radial 1 suggests an anterior component to the bidirectionality. However, in treatment F, the chondrification of the proximal radials of the anal fin is clearly bidirectional. In contrast to the proximal radials, chondrification of the distal radials of the anal fin proceeds clearly from anterior to posterior in treatment F; however, in treatment S, the formation of distal radial 3 suggests a bidirectional formation (Figure 2).

In both the dorsal and anal fins, proximal radials form before the distal radials; thus a proximo-distal direction of chondrification and ossification is similar in both fins.

Directions of chondrification and ossification among serially repeated elements for the median fins are highly conservative. In general, the formation of the different caudal elements proceeds from anterior to posterior. The development of the dorsal and anal fins (Table 3) proceeds (1) from proximal to distal within each fin and (2) primarily bidirectionally for the radials. However, there might be a tendency to shift from a bidirectional pattern towards an antero-posterior pattern in the formation of the dorsal and anal fins. With the exception of the epurals in the caudal fin and the distal radials of the anal fin, the treatment has no effect on the direction of chondrification and ossification.

4. Discussion

4.1. Developmental Plasticity

The morphology of salmonids is known to be plastic. Numerous studies [5159, 69] have demonstrated the effect of water velocity on body size and shape for different species of salmonids. None of these studies addressed the potential effect of water velocity on the skeletal anatomy of these fishes. However, Fischer-Rousseau et al. [36], in comparing the development of brook charrs Salvelinus fontinalis captured in stream and lake, not only demonstrated differences in shape, level of integration, and structure of variation between the habitats but also observed that ossification events (for the complete postcranial skeleton) in stream fishes occur at a significantly smaller size. It was suggested that a muscular functional demand was potentially responsible for this difference in ossification timing. In the present study, we show experimentally for the first time the effect of water velocity on the timing of first appearance, and to a lesser extent on the relative timing of appearance and the direction of chondrification and ossification. Water velocity has a differential effect on the absolute timing of chondrification versus ossification, as well as on the different median fins. However, the observed timing changes have only a low impact on the order of events within the sequences themselves.

Developmental plasticity on the timing and rate of development in fishes has already been suggested. Sfakianakis et al. [70] reported the absence of temperature effect on the pattern of osteological development of the postcranial skeleton of the common pandora Pagellus erythrinus but demonstrated an effect on the rate of development. Fuiman et al. [71] mentioned that temperature might accelerate the rate of development more than the growth rate. In the present study, similar patterns have been found for S. alpinus although water velocity rather than temperature was tested. The treatment had no significant effect on growth although significant changes in the absolute timing were found where elements ossify earlier in the higher velocity treatment. The absence of effect on growth during the 90-day period is most likely owing to the weak difference of water velocity between our two treatments; significant differences in growth rate of S. alpinus were found associated with greater differential water velocities [58].

Certain environmental or experimental constraints induce differential mechanical stress and loading by muscular forces, which may involve plastic response of an organism or more specifically of its bones [7275]. Such mechanical stress may originate from functional and behavioural needs [35, 41, 7678]. Rearing fishes in different water velocities most likely affected their locomotion and therefore the mechanical stress exerted on the musculoskeletal system. Different lines of evidence have demonstrated the correlative link between muscles and cartilages-bones ((1) induced paralysis in developing chick embryo [79, 80], (2) comparison between normal and muscleless mouse embryo mutants [74, 81], and (3) osteological comparison in asymmetrical training in humans [82, 83]). Responses correlated to higher muscular stress or activities include increases in numerous structural parameters (total bone volume, bone mineral density, bone mineral content, bone mass, and osteon density [74, 75, 7886]). In addition to these structural differences, Hosseini and Hogg [79] observed a delay in the timing of first appearance for the ossification for 15 elements in the paralyzed chick embryo. Furthermore, earlier ossifying bones were less affected than later developing bones [79].

Overall, the ossification of fin supports and fin rays is faster in the higher water velocity, and the responses are more important in the dorsal and anal fins whereas the pattern is globally more conservative for the caudal fin. Although the formation of the ventral elements of the caudal fin (haemal spine 2, parhypural, and hypurals 1–4) is similar in both treatments, the sequences of chondrification and ossification of the dorsal elements differ. The absolute timing of epural 1, uroneural 1 and neural arch 2 display a greater plasticity than the parhypural and hypurals 1–4. These dorsal structures of the caudal fin are located in proximity of the stegural, where mechanical stresses are important due to the thrust propagation along the posterior part of the notochord [87]. Most of the mechanical stresses associated with the subcarangiform locomotion might be taken by the stegural that strengthens the caudal fin. Given that skeletogenesis and to a greater extent ossification are sensitive to environmental inductions and mechanical stimulations [41], the timing of dorsal elements might be more plastic. Changes in timing of ossification of certain elements provide adaptive responses of skeletal systems involved in functional needs [35, 75, 87].

Most of the elements forming the skeleton of the median fins are cartilage bones, either endochondral or perichondral, with the exception of the ural chordacentra associated with the caudal fin that do not have cartilaginous precursors. Differential responses between chondrification and ossification to an experimental treatment have already been reported in fishes. Mabee et al. [35] noted differential responses to temperature treatment between cranial cartilage bones and dermal bones. In zebrafish, cartilage and dermal bones did not differ in their average timing of development among temperatures, but cartilage bones were more variable than dermal bones [35]. Campinho et al. [48] reported the ossification to be more sensitive than chondrification to a temperature treatment. A difference in velocity, as minimal as the one tested experimentally in the present study, has an effect on the timing of the ossification events and to a lesser extent on the chondrification events. Since the effects differ in the two sets of developmental events, we can assume that the pattern of ossification is partially independent from the pattern of chondrification. Partial independence between these sequences has already been shown with regards to the development of the caudal fin of some catostomids [7]. Most of the chondrification of the median fins of S. alpinus is completed during the first 50 days of the experiment whereas most of the ossification is not initiated before day 50. Because of the duration of exposure to the experimental constraint, it is suggested that the chondrification pattern is more under genetic control with a minimal epigenetic component whereas the ossification pattern, although also controlled genetically, might have a greater epigenetic component than the chondrification pattern. Evidently, the process of endochondral bone formation is strongly dependent of the coordination between chondrocyte and osteoblast differentiation associated with a cascade of several factors at the molecular and cellular levels [88]. However, skeletogenesis might be regulated by epigenetic effects such as mechanical stimulation and environmental modifications [41, 74, 85]. Thus, one may argue that because of the potential epigenetic effect of environmental modifications on skeletal structure development, the cartilaginous-bony dependence at the molecular level is not necessarily a warranty of a complete dependence of cartilages and bones at the structural level.

Magnitude rank shift is very low in average for each fin, the majority of elements shifting only from one or two ranks. This low variability is responsible for the highly conservative sequence of events whether comparing chondrification sequences or ossification sequences. Such low level of variability had been reported for the cranial elements of Danio rerio and Betta splendens reared at different temperatures [35]. Our results suggest that median fin structures are relatively well integrated given the conservative sequence and the low level of variability.

4.2. Anatomical Variation

We provided the first description of the chondrification and ossification sequences for the median fins of Salvelinus alpinus. This description completes previous descriptions of the median fin supports of S. alpinus [33, 44, 67]. The characterization of this developmental sequence also provided the opportunity to report anatomical variation among the components of the median fins.

In addition to the plasticity in the timing and the position of events in the sequence, certain elements display a greater level of anatomical variation. Although we have only qualified these anatomical variations, they are worth mentioning and interpreting; further studies would have to quantify this variation and test the effect of the water velocity treatment on the pattern of variation. The greatest variation in terms of anomalies (e.g., fusion, additional elements) and magnitude of rank-shift is present in the most anterior and most posterior proximal and distal radials of the dorsal and anal fins; these elements correspond to the peripheral serially repeated elements of these series. Two hypotheses are suggested to explain the topographic position of this variation. First, the variability could be the result of suboptimal experimental conditions [46, 89]; however, rearing conditions were selected to be acceptable or optimal for S. alpinus [58, 62]. Significantly higher occurrence of skeletal malformations of the vertebral column were reported in sea bass Dicentrarchus labrax [90] and red sea bream Pagrus major [91] as a result of a higher water velocity. None of the anatomical variation observed in the dorsal and anal fins of S. alpinus have been identified in the sea bass and the red sea bream. Thus, the median fin variation of S. alpinus does not seem to be related to suboptimal rearing condition or the water velocity treatment. We are suggesting a second hypothesis in which the source of variation might be in relation to the peripheral topographic position of serially repeated elements at the edge of the fins. In a bidirectional mode of formation (such as the proximal radials of the dorsal and anal fins), peripheral serially repeated elements are the latest elements to appear in the sequence. One might expect that developmental constraints, such as integration or canalization, might be relaxed in peripheral serially repeated elements, and thus epigenetic constraints be more influential.

4.3. Modularity

Phenotypic developmental modules in fishes have been suggested based on direction of ossification. Mabee et al. [92] have identified various fin patterning as well as one positioning module and three patterning modules with respect to the median fins of actinopterygians. In the present study, we can only address two of the three patterning modules. Faustino and Power [93] and Mabee et al. [92] mentioned that the overall pattern of development of the dorsal and anal fin supports in teleosts is similar with a few exceptions. The dorsal and anal fins are characterized by a patterning module in which the skeletal elements of both fins differentiate bidirectionally. In S. alpinus, the proximal radials of the anal and dorsal fins form and ossify bidirectionally whereas no clear pattern was identified for the distal radials of both fins. Similar patterns of chondrification and ossification for these two fins have been described for sea bream Sparus aurata [93] and brook charr Salvelinus fontinalis [36]. Not only the patterns of formation are similar, but the initiation of ossification for both fins is similar in both treatments, the absolute timing is similar, and the response in terms of rank shift is also similar. Our results document for the first time the presence of the dorsal-anal fin patterning module in the Salmoniformes.

The second patterning module concerns the relation between exo- and endoskeleton in the development of all median fins [92]. In the dorsal and anal fins of S. alpinus, the proximal radials are the first structures to develop as cartilages, followed by the distal radials and their respective fin rays. The formation of the lepidotrichia requires the presence of the radials but not their ossification. Similar sequences have been reported previously for different teleosts [9295]. The conservation of this patterning is expected since an induction mechanism has been shown between the radial and the corresponding fin rays in the bastard halibut, Paralichthys olivaceus [96].

Patterning modules have been investigated using developmental sequences at the species level without consideration of environmental constraints. In the present study, we showed that an experimental constraint, namely the water velocity, has the potential to modify direction of both chondrification and ossification. Thus, one has to be aware that one of the defining criteria of patterning modules could be plastic. However, the plastic response has a tendency to be coordinated between modular components. On the other hand, this coordination of phenotypic plastic response could be used to test the robustness of developmental patterning modules. In further studies, it would be worth testing the cohesion and variation within modules on larger sample size as well as on a wider range of water velocity.

Abbreviations

Ard 3–12:Anal distal radials 3–12
Arp 1–12:Anal proximal radials 1–12
chcU2, 4–6:Ural chordacentra 2, 4–6
c.pl:Cartilaginous plate
Drd 2–13:Dorsal distal radials 2–13
Drp 1–13:Dorsal proximal radials 1–13
E1-2:Epural 1-2
H1–H6:Hypural 1–6
HA:Haemal arch
HS:Haemal spine
HS2:Haemal spine 2
NA1-2:Neural arches 1-2
nc:Notochord
NS2:Neural spine 2
PH:Parhypural
PU1:Preural centrum 1
ST:Stegural
UN1-2:Uroneural 1-2.

Acknowledgments

The authors specially thank M. Bélanger and V. Canévet (UQAR) for providing substantial assistance with the clearing and staining of the specimens. They thank B. Archer who was responsible for the fish rearing and the Centre Aquacole Marin de Grande-Rivière (Québec, Canada) for providing facilities during the experiment. They thank G. Arratia, K. Chu, L. Fischer-Rousseau, and B. Vincent for stimulating discussions. This research was supported by Natural Sciences and Engineering Research Council of Canada operating grant 238612 to the first author. Additional support was granted by a FIR from the Université du Québec à Rimouski to R. Cloutier and N. R. Le François.