bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
1 Monsters with a shortened vertebral column: A population phenomenon in
2 radiating fish Labeobarbus (Cyprinidae)
3
4 Alexander S. Golubtsov1, Nikolai B. Korostelev1, Boris A. Levin2,3*
5 1 Severtsov Institute of Ecology and Evolution, Russian Academy of Sciences, Moscow, Russia
6 2 Papanin Institute for Biology of Inland Waters, Russian Academy of Sciences, Borok, Russia
7 3 Cherepovets State University, Cherepovets, Russia
8
9 ☯ These authors contributed equally to this work.
10 * [email protected]
11
12
13 Short title: Monstrosity in fish as novelty?
14
15
16
17
18
19
20
21
22
23
24
25
26 1 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
27 ABSTRACT
28 The phenomenon of a massive vertebral deformity was recorded in the radiating Labeobarbus
29 assemblage from the middle reaches of the Genale River (south-eastern Ethiopia, East Africa).
30 Within this sympatric assemblage, five trophic morphs – generalized, lipped, piscivorous and two
31 scraping feeders – were reported between 1993 and 2019. In 2009, a new morph with prevalence
32 of ~10% was discovered. The new morph, termed ‘short’, had an abnormally shortened vertebral
33 column and a significantly heightened body. This type of deformity is common in farmed Atlantic
34 salmon and other artificially reared fish, but is rare in nature. In the Genale Labeobarbus
35 assemblage, the deformity was present exclusively within the generalized and lipped morphs. The
36 short morph had between seven and 36 deformed (compressed and/or fused) vertebrae. Their body
37 height was positively correlated with number of deformed vertebrae. In another collection in 2019,
38 the short morph was still present at a frequency of 11%. Various environmental and genetic factors
39 could contribute to the development of this deformity in the Genale Labeobarbus, but based on
40 the available data, it is impossible to confidently identify the key factor(s). Whether the result of
41 genetics, the environment, or both, this high-bodied phenotype is assumed to be an anti-predator
42 adaptation, as there is evidence of its selective advantage in the generalized morph. The Genale
43 “monstrosity” is the first reported case of a massive deformity of the vertebral column in a natural
44 population of African fishes.
45
46
47
48
49
50
51
52
2 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
53 “We have also what are called monstrosities; but they
54 graduate into varieties. By a monstrosity I presume is
55 meant some considerable deviation of structure in one
56 part, either injurious to or not useful to the species,
57 and not generally propagated. If it could be shown
58 that monstrosities were even propagated for a
59 succession of generations in a state of nature,
60 modifications might be effected (with the aid of
61 natural selection) more abruptly than I am inclined to
62 believe they are.”
63 Darwin (1860, pp. 46, 426).
64
65 INTRODUCTION
66 The emergence and establishment of morphological novelties in a population is central to
67 morphological evolution, although this process remains poorly understood. As initially highlighted
68 by Darwin [1], there is no clear discrimination between morphological abnormalities
69 (monstrosities) and regular variation. Morphological abnormalities that are caused by genetic
70 factors – and hence may be shaped by natural selection – are particularly interesting to evolutionary
71 biologists as potential novelties.
72 Ray-finned fishes form one of the major vertebrate groups, comprised of more than 30,000
73 species with extreme variations in morphology and body plans [2]. For example, the ocean
74 sunfishes (Molidae, Tetraodontiformes) exhibit one of the most impressive evolutionary
75 transformations of the axial skeleton among ray-finned fishes [3,4]. It has been suggested that the
76 enigmatic loss of caudal fin and caudal part of the vertebral column is related to “the existence in
77 collections of a number of malformed adult tetraodontiforms without caudal fins” [5,6].
3 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
78 Various body and skeletal deformities in fish have been reviewed in many works [7-13];
79 earlier reports can also be found in Dawson’s bibliographies [13-17]. Many of these deformities
80 apparently can be treated as monstrosities, particularly when the vertebral column is shortened
81 without pronounced curvature, resulting in an altered body form that is extremely short and high.
82 Such deformities are known in at least 26 species of 15 families (S1 Table). Remarkably, all
83 individuals of Aischgrunder Karpfen, the South German carp breed (extinct since 1956), exhibited
84 such deformity [18,19].
85 The riverine adaptive radiations in the large African barbs of the genus Labeobarbus
86 Rüppell 1835 (Cyprinidae) – the dominating fish group in the waters of the Ethiopian Highlands
87 – have been studied by us for almost 30 years. Labeobarbus belongs to the African Torini, a lineage
88 of hexaploids [20-24] that originated in the Middle East via hybridization of tetraploids (maternal
89 Tor lineage) and diploids (paternal Cyprinion lineage), and then dispersed throughout Africa [25].
90 Although the continental-scale phylogeny of Labeobarbus is still poorly resolved [26], the mt-
91 DNA based phylogeny of this group in Ethiopian waters is relatively well-studied [27-32]. It is
92 important to note that the barbs inhabiting the southeastern part of Ethiopia (drained by the Wabi-
93 Shebele and Juba river systems into the Indian Ocean) form a monophyletic group, Labeobarbus
94 gananensis (Vinciguerra, 1895) complex. This is a sister group to Labeobarbus inhabiting
95 enclosed basins of the Ethiopian Rift Valley, as well as all Labeobarbus from the waters of western
96 and northern Ethiopia belonging to the Omo-Turkana and Nile systems, and additional some
97 Kenyan barbs [27,31,32].
98 In addition to the famed Labeobarbus radiation in Lake Tana [33-63], this group has also
99 experienced parallel riverine radiations in four Ethiopian basins. These riverine radiations are
100 similar in terms of their morphological and ecological differentiation, but differ in the ways of
101 genetic divergence [32]. Outside of Ethiopia, the riverine Labeobarbus radiation seems to occur
102 in the Inkisi River basin (Lower Congo) [64].
4 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
103 The Labeobarbus assemblage in the Genale River (the main tributary of the Juba River) is
104 the only radiation occurring in southeastern Ethiopia. It is geographically separated from the other
105 Ethiopian Labeobarbus radiations by the Ethiopian Rift valley. It exhibits the deepest divergence
106 in terms of morphology, ecology and genetics compared to the other riverine radiations [32].
107 During an intensive sampling of the middle Genale River in 2009, we observed a substantial
108 portion of barbs with extremely short and high bodies [31].
109 In 1993, an assemblage of the sympatric morphologically distinct morphs of the L.
110 gananensis complex was discovered in the middle reaches of the Genale River [65]. This area was
111 re-sampled twice in the 1990s [66,67], then in 2009 [27,31], and again in 2019. Together with the
112 omnivorous generalized morph of L. gananensis and the scraping-feeder L. jubae (Banister, 1984),
113 three additional morphs specialized in their morphology and feeding habits were always present
114 in these samples from the middle Genale. These include the lipped morph of L. gananensis (having
115 hypertrophied lips), the scraping morph Labeobarbus sp. 1 (provisionally called ‘smiling’) and the
116 large-mouthed morph Labeobarbus sp. 2 (called ‘piscivorous’), as well as a hybrid morph called
117 ‘smiling hybrids’. In 2009, a new morph, ‘short’, was discovered, which has an abnormally short
118 and high body (Fig. 1). It composed approximately 10% of the total barb catch in both 2009 and
119 2019 samplings [31].
120
121 Fig. 1. External appearance (left) and x-ray images (right) of the generalized and short
122 Labeobarbus morphs from the middle Genale assemblage. A: generalized (SL = 242 mm), B:
123 short, with eight deformed vertebrae, marked by red asterisks (SL = 245 mm), C: short, with
124 28 deformed vertebrae (SL = 216 mm), and D: short, with 31 deformed vertebrae (SL = 201 mm).
125
126 Trophic resource partitioning has been demonstrated among most morphs from the middle
127 Genale Labeobarbus assemblage [31]. In contrast to the other morphs, the short morph does not
128 diverge from the generalized morph in its diet [31]. Moreover, based on mtDNA markers, the short
5 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
129 morph shares a common gene pool with the generalized and lipped morphs [31]. Preliminary
130 investigation of the vertebral column revealed that each short individual had a substantial number
131 of deformed (compressed and/or fused) vertebrae (Fig. 1).
132 The objectives of this study were: (1) to investigate the structure of the vertebral column
133 in the short morph compared to other sympatric Labeobarbus morphs from the Genale River
134 assemblage, (2) to analyze morphological differences between the short and related (generalized
135 and lipped) morphs, with clarification of the contribution of the vertebral column shortening to
136 these differences, and (3) to discuss the possible causes of the emergence of the aberrant morph in
137 light of its life history characteristics and decade-long presence in the local Labeobarbus
138 population.
139
140 MATERIALS AND METHODS
141
142 Ethics statement
143 Fish were collected in southeastern Ethiopia under the umbrella of the Joint Ethiopian-Russian
144 Biological Expedition (JERBE), with permission from the National Fishery and Aquatic Life
145 Research Center under the Ethiopian Institute of Agricultural Research and the Ethiopian Ministry
146 of Innovation and Technology. Fish were sacrificed humanely using an anesthetic overdose
147 (American Veterinary Medical Association). The experiments were carried out in accordance with
148 the rules of the Papanin Institute of Biology of Inland Waters (IBIW), Russian Academy of
149 Sciences, and approved by IBIW’s Ethics Committee.
150
151 Study sites and sampling
152 The Ethiopian Highlands are divided by the Rift Valley into the western and eastern
153 plateaus [68]. The southern part of the eastern plateau is drained by the Wabi Shebeli (Webi
154 Shabeelle or Uebi Scebeli) and Juba (Jubba) river systems into the Indian Ocean via Somalia.
6 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
155 Three large tributaries, the Genale (Ganale Doria), the Dawa (Daua) and the Weyb (Gestro), meet
156 near the Ethiopia-Somalia border to form what is known as the Juba River within Somalia. The
157 Genale, Weyb and Wabi Shebeli originate in the Bale Mountains, the highest part of the eastern
158 plateau (Fig. 2).
159 Since 1990, we have sampled a total of 28 sites across the Wabi-Shebele and Juba river
160 systems (dots in Fig. 2). Among these, Labeobarbus were sampled from 15 sites, and not found in
161 the other 13 sites. An additional four sites important for the present work are numbered in Fig. 2:
162 (1) main sampling site, the middle reaches of the Genale River (5º42' N 39º32' E, altitude 1125 m
163 above sea level, asl); (2) the upper reaches the Genale River north of Kibre Mengist (6°02'40" N
164 39°07'19" E, 1324 m asl); (3) the middle reaches of the Welemele River, a northern tributary of
165 the Genale (6°14'30" N 39°47'22" E, altitude c. 1040 m asl); and (4) the upper reaches the Awata
166 River, northern tributary of the Dawa (5º47'07"N 38º55'46"E, 1630 m asl). All sampling sites are
167 listed in S2 Table. Most samples were taken in the different years from the end of January to mid-
168 May during the main rainy season in southeastern Ethiopia [69], but before substantial rise of water
169 level in the rivers.
170
171 Fig. 2. Map of sampling sites in the Juba and Wabi Shebele drainages. Black and red circles
172 indicate sites where Labeobarbus was detected; white circles indicate sites where barbs were not
173 found. Numbered red circles designate the most important sites for the current study – 1: Genale
174 River, middle reaches, 2: Genale River, upper reaches, 3: Welmel River, and 4: Awata River.
175 Geographical coordinates for all sites are given in S2 Table.
176
177 The main sampling site (no. 1, Fig. 2) included two habitats: the continuous pool starting
178 from the shallows around the island, characterized by moderate current and sandy/silty bottom,
179 and the upper part of continuous rapids with small waterfalls, rocky bottom and very fast current
180 (S3 Fig.). Samples were mostly taken from the pool by gill netting (at dawn and overnight), but an
7 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
181 additional sample was obtained from the rapids with cast nets in 2009. The cast net catches differed
182 substantially from the gill net catches, in terms of fish size and frequency of Labeobarbus morphs
183 (S4 Table), therefore these catches were analyzed separately. The remaining sites (nos. 2-4) were
184 sampled by both cast and gill netting during the day, and the catches were analyzed together.
185 The selected fish were maintained for several hours in large barrels submerged into the
186 river, then killed with an overdose of MS-222 and preliminarily investigated for morphology.
187 Tissue samples were taken for genetic and stable isotope analyses because the genetic and trophic
188 specialization studies have been conducted on the same individuals [27,31,32]. Most specimens
189 were preserved in 10% formalin in the field, and subsequently transferred to 70% ethanol in the
190 laboratory. Of the 291 specimens collected in 2009, 132 were preserved with salt to make the dried
191 bone preparations. These specimens were dissected for determination of sex and gonad maturity
192 stages.
193
194 Examination of vertebral column
195 Most formalin-preserved specimens were x-rayed. Based on the film or digital radiograms,
196 the structure of vertebrae were examined. The total vertebral number and numbers of pre-dorsal,
197 pre-anal, trunk, transitional and caudal vertebrae were counted according to Naseka [70].
198 Examination of vertebral structure and counts of vertebrae mentioned above were made by means
199 of visual inspection of the dried bone preparations. The head magnifier (×4) was used for film and
200 preparation inspection if necessary. In total, the vertebral column structure was studied in 364
201 Genale barbs representing all morphs. This number includes the radiograms of 34 and 45
202 specimens collected before 2009 and in 2019, respectively, as well as 285 radiograms and dried
203 bone preparations of fish collected in 2009.
204 Witten et al. [11] propose a classification of vertebral column deformities that are
205 repetitively observed in salmon (Salmon salar L.) under farming conditions. In our material, we
206 observed all nine types of deformities united by these authors into the category of ‘compressed
8 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
207 and/or fused vertebral bodies’. Taking into account the relatively small size of our samples, we
208 did not consider the different types of such deformities separately. All vertebrae with compressed
209 and/or fused bodies were considered herein as deformed vertebrae. The deformities of other
210 categories [11] were not recorded in our material, except in one specimen of the smiling morph
211 that had vertebral column curvature (kyphosis). Other deformities such as a presence of abnormal
212 additional ribs and neural spines, and splitting of neural spine [71] were recorded in our material,
213 however they were not analyzed because they were not related to the shortening of the vertebral
214 column.
215 Examination of other morphologic characters
216 Prior to preservation, the following measurements were taken (with a ruler, to the nearest
217 0.5 mm) on all specimens: standard length (SL), head length (HL), dorsal spine length (DL),
218 maximum (H) and minimum (h) height of body, and length of caudal peduncle (cpd). For further
219 analysis, the body length (BL) was calculated as SL minus HL.
220 Of all morphs sampled in 2009 (N = 129, SL = 98-308 mm), we randomly selected a sub-
221 set of formalin-preserved specimens and took 23 measurements in the laboratory (with a caliper,
222 to the nearest 0.1 mm): standard length (SL), head length (HL), snout length (R), orbit horizontal
223 diameter (O), interorbital width (IO), postorbital length of head (PO), height of head at mid-orbit
224 (Ch), height of head at occiput (CH), mouth width (MW), dorsal spine length (DL), pre-dorsal
225 distance (PrD), post-dorsal distance (PD), caudal peduncle length (cpd), length of dorsal fin base
226 (lD), maximum (H) and minimum (h) height of body, length of pectoral (lP) and ventral (lV) fins,
227 pectoral-ventral distance (PV), ventral-anal distance (VA), height of anal fin (hA), length of anterior
228 barbel (Ab), and length of posterior barbel (Pb). In the present work, only data on these extended
229 measurements were included for 22 and 28 specimens of the generalized (SL = 108-257 mm) and
230 short (SL = 119-244 mm) morphs, respectively. All measurements were taken by the same person
231 for consistency [72].
9 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
232 In most sampled specimens, the following counts were taken: number of scales in the
233 lateral line, number of scale rows above and below lateral line, number of predorsal scales, number
234 of scales around the caudal peduncle, number of gill rakers on the first gill arch, number of
235 branched rays in the dorsal fin.
236
237 Examination of age and growth
238 In most salt-preserved specimens of the different morphs, age was determined from
239 vertebrae. Annual rings were counted on the anterior and posterior sides of the fifth vertebra after
240 removing the dried intervertebral discs. Counts were made with a binocular microscope at
241 magnification ×16, with a drop of glycerin for contrasting. Age estimates were calibrated with
242 vertebrae from the artificially reared Labeobarbus from Lake Tana [63].
243 Growth rate estimates for the different Labeobarbus morphs were based on comparisons
244 of individual size variation in the different year classes. We did not calculate values of growth rate
245 or parameters of the von-Bertalanffy equation because of the small sample sizes.
246 Statistical analysis
247 Various R packages [73] realized in R-studio v.1.2.5033 were used for statistical analyses
248 and plot construction: summarytools library was used for obtaining basal descriptive statistics,
249 ggpubr library was used to analyze Pearson correlation and plot the result,
250 posthoc.kruskal.dunn.test function was applied for posthoc Dunn’s test, FSA library was applied
251 for Mann-Whitney U test, prcomp function was used for the principal component analysis (PCA)
252 and plotting the results, ggplot2 library was applied for plotting boxplots and histogram bars. The
253 proportions of head and body were used for analyses of single characters, as well as for PCA – all
254 measurements were divided by head length. Data was standardized for PCA.
255
256 Deposition of material
10 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
257 All samples were deposited to the Severtsov Institute of Ecology and Evolution Russian
258 Academy of Sciences and Papanin Institute for Biology of Inland Waters Russian Academy of
259 Sciences under provisional labels of the Joint Ethiopian-Russian Biological Expedition.
260
261 Abbreviation of morph names
262 The names of the morphs from the Genale Labeobarbus assemblage have been assigned
263 the following abbreviations: short SH, generalized GN, lipped LP, L. jubae JB, smiling SM,
264 smiling hybrids HB, and piscivorous PS.
265
266 RESULTS
267 Short definition
268 Initially, in the field, the short morph was identified based on altered body proportions:
269 these individuals had relatively higher and shorter bodies compared to the normal individuals (Fig.
270 1). Preliminary investigation of the vertebral column revealed that each short individual had a
271 substantial number of deformed vertebrae. The data on the structure of the vertebral column for
272 all morphs are presented in Table 1.
273
274 Table 1. Occurrence of individuals with various numbers of deformed vertebrae in the different
275 morphs in the total sample (1993-2019) (N, number of individuals; %dv, percent of individuals
276 with deformed vertebrae in a particular morph).
277
Number of deformed vertebrae per individual
Morph N %dv 0 1 2 3 4 5 6 7 8 ≥ 9
Shorts (SH) 52 100 1 2 49
Generalized 78 16.7 65 1 8 2 2
(GN)
11 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
Lip (LP) 26 11.5 23 2 1
Jubae (JB) 78 12.8 68 1 7 1 1
Smiling (SM) 48 14.6 41 1 4 1 1
Hybrid (HY) 36 5.6 34 2
Piscivorous 46 9.8 41 3 1 1
(PS)
278
279 We defined the SH morph as individuals with seven or more deformed vertebrae. It is
280 important to note that the mouth structure in the SH morph was similar to that in GN and LP
281 morphs that had normal vertebral columns. These morphs – with either shortened or normal
282 vertebral columns – differed in the degree of lip development (S5 Fig.).
283 There were almost no individuals with markedly shortened bodies among the trophically
284 specialized morphs (JB, SM, HY, PS). The only PS individual with eight deformed vertebrae
285 (Table 1) displayed noticeable body shortening (S6 Fig.).
286 Temporal and spatial distribution of the short morph
287 In the samples taken before 2009, SH individuals were not discovered among the 34
288 radiographed individuals of the different morphs, or among the hundreds of superficially examined
289 barbs in catches of 1993, 1997 and 1998. To estimate SH prevalence in the samples of 2009 and
290 2019, we analyzed the total gill net catches. The catch composition appeared to be quite similar in
291 2009 and 2019 (S7 Table). The ratios of morphs remained roughly stable. It is important to note
292 that SH prevalence in the total Labeobarbus catches was almost unchanged: they represented 10%
293 of the 400 barbs sampled in 2009, and 11% of the 179 barbs sampled in 2019.
294 We did not sample the SH morph from the upper reaches of the Genale River (sampling
295 site no. 2, Fig. 2), where most other morphs common to the main sampling site were found.
296 Moreover, the SH morph was not recorded from other sampling sites in the Wabi-Shebele and
297 Juba river systems (Fig. 2, S2 Table) or from other parts of Ethiopia (sampling sites may be found
12 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
298 in [32,74,75]). However, at sampling site no. 3, we found barbs that had shortened bodies but
299 without deformed vertebrae (discussed below).
300 Short vs generalized and lipped morphs: plastic characters
301 As expected, body height of the SH morph was positively correlated (Fig. 3A) with the
302 number of deformed vertebrae, which in turn was negatively correlated with body length (Fig.
303 3B). In other morphs of the Genale Labeobarbus assemblage, no such correlations were found.
304
305 Fig. 3. Pearson correlation of (A) body height and (B) body length with number of
306 deformed vertebrae in the short (SH) morph.
307
308 The changes of body proportions in the SH morph (caused by the deformity of the vertebral
309 column) resulted in a pronounced difference in appearance between this and the other morphs.
310 Based on PCA, the SH and GN morphs were well differentiated from each other in PC1 (Fig. 4).
311 Indices of 21 measurements relative to head length (HL) rather than SL were used in the PCA
312 because of the variable influence of vertebral deformity on SL in the short individuals (Fig. 3B).
313 PC1 explained > 45% of the variance, while PC2 was less than 24%. Eigenvectors of the 10 most
314 loaded characters for PC1 and PC2 are given in S7 Table.
315
316 Fig. 4. PCA of short (SH) and generalized (GN) morphs based on 21 indices of head and
317 body proportions.
318
319 Using the a larger data set (which included the gill net catch of 2009, S5 Table), we
320 compared the characters that had the strongest influence on body form (head length, heights of
321 body and caudal peduncle) between the SH, GN and LP morphs (Fig. 5). The LP and SH morphs
322 were similarly characterized by relatively longer heads (Fig. 5A). In the LP morph, the longer head
323 length arose from elongation of the snout (S5 Fig.), while in the SH morph the head length
13 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
324 appeared longer relative to the shortened body. The LP, SH and lipped SH morphs all had
325 relatively longer heads than the GN; the lipped SH morph had a notably longer head than the SH
326 morph, but the difference was not significant (Fig. 5A). The relative heights of body and caudal
327 peduncle in the SH morph differed significantly from both GN and LP, however there were no
328 such differences between the lipped SH and any of the other morphs (Fig. 5B-C). Notably, lipped
329 SH exhibited the most variation in relative heights of body and caudal peduncle compared to the
330 other morphs. This variation was caused by different degrees of snout elongation and variation
331 caused by body shortening (determined by varying number of deformed vertebrae), which seemed
332 to act additively.
333
334 Fig. 5. Indices of (A) head length, (B) body height, and (C) caudal peduncle height in
335 generalized, lipped, lipped short, and short morphs sampled in 2009. Median is shown as the
336 horizontal black line inside the box. The box represents 1st and 3rd quartiles of variation.
337 Lowercase letters above the boxplots indicate significant differences between morphs (p < 0.05,
338 Dunn’s test with Holm adjustment of p-value).
339
340 The comparison of barbs sampled in 2009 and 2019 yielded unexpected results. Although
341 the later sample was rather small, significant differences in relative head length and height of
342 caudal peduncle, as well as a nearly significant (p = 0.054) difference in body height, were revealed
343 between GN sampled in 2009 and 2019 (Fig. 6). Remarkably, in only 10 years the GN morph
344 became significantly more similar to the SH morph in all three of the indices most important for
345 their discrimination. No differences in SL and absolute values of the corresponding measurements
346 were found between samples of 2009 and 2019 (Mann-Whitney U test).
347
348 Fig. 6. Indices of (A) head length, (B) body height, and (C) caudal peduncle height of
349 generalized and short morphs sampled in 2009 and 2019. Median is shown as the horizontal black
14 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
350 line inside the box. The box represents 1st and 3rd quartiles of variation. Lowercase letters above
351 the boxplots indicate significant differences between morphs (p < 0.05, Mann-Whitney U test; * -
352 nearly significant, p = 0.054).
353
354 Short vs generalized and lipped morphs: structure of vertebral column and other counts
355 Comparison of vertebral counts between the SH morph and the GN and LP morphs did not
356 reveal any differences in total number of vertebrae or in numbers of vertebrae in the three segments
357 of the vertebral column (trunk, transitional and caudal). Moreover, there were no differences in
358 the numbers of pre-dorsal and pre-anal vertebrae (S9 Table).
359 Distribution of the deformed vertebrae along the vertebral column was not random, and
360 was different in the SH morph compared to the other morphs (Fig. 7). When the data on the
361 positions of the deformed vertebra was merged for individuals with the different total number of
362 vertebrae, the frequency estimates for each specific vertebra in the caudal region were obscured.
363 For this reason, only the data for individuals with 41 total vertebrae are presented in Fig. 7. The
364 respective data for individuals with 40, 42 and 43 vertebrae are shown in S10 Figure.
365
366 Fig. 7. Distribution of deformed vertebrae along the vertebral column in (A) short morph
367 (n = 17), and (B) other morphs (n = 47). Only data for individuals with 41 vertebrae are shown
368 here; the deformities were not analyzed for the first four vertebrae (Weberian apparatus) or the 41st
369 (pre-ural 1) vertebra for technical reasons.
370
371 As shown in Fig. 7A , the lowest frequencies of each specific vertebra deformity in the SH
372 morph were found in 5th vertebra, and also in the caudal region, where the minimum frequencies
373 were found in the 36-40th vertebrae. In the other morphs, however, the deformed vertebrae
374 occurred mostly in the posterior half of the vertebral column, with maximum frequencies in the
375 37-38th vertebrae (Fig. 7B).
15 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
376 No differences were found in the counts of scales, gill rakers and fin rays between the GN
377 and LP morphs, on the one hand, and the SH and lipped SH morphs, on the other hand. No changes
378 in the counts for samples 2009 and 2019 were recorded.
379
380 Short vs generalized and lipped morphs: size, age, growth rate, sex ratio and gonad
381 conditions
382 The range of size variation in the SH morph (SL 107-291 mm) was less than that in any
383 other morph from the middle Genale assemblage in 2009 and 2019. The ranges for the GN and LP
384 morphs were 51-375 mm and 106-459 mm, respectively (S4 Table). Most of the salt-preserved
385 barbs were aged from vertebrae. A maximum age of six years was recorded for the SH and HB
386 morphs, seven years for the SM morph, eight years for the GN morph, nine years for the JB morph,
387 12 years for the LP morph and 14 years for the PS morph. There was no correlation between the
388 percentage of deformed vertebrae and age.
389 Analysis of size variation within each year class of the SH and GN morphs revealed that
390 when SL was used as a measure of linear growth, SH did not differ in size from GN during the
391 first four years (S11 Fig.). Modification of the vertebral column in the SH morph, however,
392 apparently resulted in shortening of their body and respective reduction of SL. Hence, it is
393 reasonable to consider the head length as a measure of linear growth in comparisons of the SH
394 morph with the other morphs. In this case, the SH morph demonstrated a markedly fast growth
395 rate during the first four years, which then leveled off during the fifth and sixth years (Fig. 8).
396
397 Fig. 8. Age vs head length in generalized (GN) and short (SH) morphs in 2009. Median is
398 shown as the horizontal black line inside the box. The box represents 1st and 3rd quartiles of
399 variation. Letters above each boxplot designate significant (red S), or non-significant (black NS)
400 differences between morphs within each year class, based on Mann-Whitney U tests (p < 0.05).
401
16 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
402 Gonads were examined in 30 SH individuals: five juveniles, 11 females and 14 males were
403 recorded. There were 16 fish at maturity stage I, 11 fish at stage II, three fish at stages II-III and
404 III. Among the 46 examined individuals of GN and LP morphs, 16 juveniles, 14 females and 16
405 males were recorded. There were 24 fish at maturity stage I, 20 fish at stage II, one fish at stages
406 II-III and one female at stage IV. Thus, neither sex ratio nor gonad condition differed in the SH
407 morph compared to GN and LP morphs.
408
409 Discussion
410 It is impossible to clearly distinguish the SH morph from the GN and LP morphs using
411 only the plastic characters reflecting their body form, which aligns with Darwin’s observation that
412 “monstrosities … graduate into varieties” [1]. Discrimination is possible, however, by using the
413 number of deformed vertebrae per individual: the GN and LP morphs have ≤ 4 such vertebrae,
414 while the SH morph has ≥ 7 (Table 1). Considering the emergence of the SH morph in the
415 Labeobarbus population predominantly composed of GN and LP morphs, one may suggest that it
416 is determined by some factor(s) acting qualitatively or the threshold reaction of the Labeobarbus
417 individuals to some factor(s) with continuously varying parameters.
418 In regards to Darwin’s idea about the propagation of monstrosities “for a succession of
419 generations in a state of nature” [1], we are motivated to estimate the period of SH existence in a
420 population. Taking into account the absence of SH individuals in our sample from 1998 and their
421 maximum age of six years in the sample from 2009, it is reasonable to suggest that the SH morph
422 appeared in the population between 1999 and 2005. The minimum size of SH individuals in the
423 sample from 2019 (SL = 177 mm) corresponds to an age of 2-3 years according to the SH sample
424 from 2009 (S4 Table). Thus, the SH individuals regularly appeared in the population from 2005
425 (or slightly earlier) until at least 2016. This situation could be explained by the influence of some
426 external factor(s) that abruptly appeared and then became permanent in the environment, or by the
17 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
427 reproduction of SH individuals during the ten year period, which obviously includes the life-span
428 of several generations.
429 A central question has arisen from the discovery of the SH morph: does the emergence and
430 10-year persistence of this morph have a genetic background? Based on mtDNA variation, there
431 is no significant genetic difference between the SH, GN and LP morphs [31], but it is still possible
432 that these differences are present but subtle. In the absence of direct genetic evidence, the factors
433 causing vertebral deformities in fish are considered more generally.
434 Factors responsible for abnormal body shortening in fish
435 There is extensive literature devoted to the question of vertebral deformities in fish. The
436 most comprehensive data available pertains to artificially propagated fish species, especially
437 farmed salmonids. Fish with abnormally shortened bodies displaying multiple compressed and/or
438 fused vertebrae are quite frequent in farmed stocks of Atlantic salmon Salmo salar L. 1758 [11,76-
439 86]. Externally (imaged in [76]), they are similar to the short morph from the Genale Labeobarbus
440 assemblage. These salmons, called ‘short tails’ or ‘short-spined’ individuals [78,81], pose a serious
441 problem for the farming industry because of their reduced commercial value [87]. Therefore,
442 special efforts have been undertaken to elucidate the factors determining the appearance of short
443 tails, however the results are vague, at least in terms of the genetic factors. Initially, results
444 suggested that the deformities were heritable [76,79], but this was questioned later [84,85]. As
445 described by Witten et al. [81, p. 244]: “one cannot rule out the possibility that some fish are
446 genetically predisposed and develop vertebral compression as a reaction to external cues”.
447 The factors that increase the risk of vertebral deformities in farmed Atlantic salmon and
448 other salmonids are given by Witten et al. [11]: bacterial and parasitic infections, vitamin C
449 deficiency, phosphorous deficiency, elevated egg incubation temperature, fast growth in under-
450 yearling smolts, inappropriate light regimes, vaccination, inappropriate water current and quality,
451 as well as environmental pollution. All of these factors (except vaccination), along with radiation
452 [71] may be present in the natural population of the middle Genale Labeobarbus. Notably, there
18 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
453 is a relationship between the increased frequency of vertebral deformities and high growth rate in
454 Atlantic salmon smolts [11,82,85]. We also found an increased growth rate in the SH morph
455 compared to the GN and LP morphs at the age of 2-4 years. Thus, the deformations of the vertebral
456 column in the SH morph could be caused by the accelerated early individual growth and
457 mechanical muscle overload of bones, as suggested for Atlantic salmon [11,82,85]. For the
458 Atlantic salmon post-smolts, a retardation of growth is reported in fish with high numbers of
459 deformed vertebrae [86], however the authors used fork length as the measure of linear growth,
460 which may substantially bias the growth estimates.
461 As for other salmonids, a short tailed brown trout, Salmo trutta L. 1758 was reported from
462 Scotland over a century ago [89]. Much later, >5% of the brown trout from a fish farm in
463 Hampshire (England) exhibited the short tail phenotype and extensive compression of the vertebral
464 column [90]. The author of the study indicated that inbreeding depression was a probable cause of
465 the abnormalities. However, this explanation is hardly applicable to the Genale barbs for two
466 reasons. First, the generalized morph, which is the ancestor of most individuals of the short morph,
467 is the most common morph among the middle Genale barbs, and hence is unlikely to suffer from
468 inbreeding depression. Second, spatial subdivision of the Genale barb populations can occur only
469 in the spawning grounds; in this case, an increase of inbreeding would only result from almost
470 perfect homing, which is highly unlikely in Labeobarbus.
471 In a study of the vertebral fusion patterns in coho salmon, Oncorhynchus kisutch (Walbaum
472 1792) [91], the author demonstrated that the distribution of vertebral fusions along the spinal
473 column differed significantly among crosses from two hatchery stocks, indicating a genetic basis
474 for this character. Recently, differences in the distribution of vertebral fusions along the spinal
475 column were found between the genetically distinct year-classes of Atlantic salmon and its hybrids
476 with brown trout and Arctic char Salvelinus alpinus (L. 1758) [92]. This result is interesting for
477 understanding the difference in distribution of deformed vertebrae along the spinal column
478 between the SH morph and other Genale Labeobarbus (Fig. 7).
19 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
479 In non-salmonid fishes, it has been reported that body shortening from the compression of
480 the vertebral column was apparently heritable in the south German common carp Aischgrunder
481 Karpfen [18], as well as in laboratory strains of the banded topminnow, Fundulus cingulatus
482 Valenciennes [93] and two poeciliids, blackstripe livebearer Poeciliopsis prolifica Miller 1960
483 [94] and guppy Poecilia reticulata Peters 1859 [95]. In zebrafish, Danio rerio (Hamilton 1822),
484 column shortening and vertebral fusions are exhibited by type I collagen mutants, as well as by
485 individuals with knockout alleles in two genes involved in type I collagen processing [96]. At the
486 same time, the increased frequency of vertebral deformities in wild-type zebrafish is induced by
487 high rearing densities [97]. The stumpbody phenotypes with the shortened vertebral column in
488 channel catfish Ictalurus punctatus (Rafinesque 1818) [98] and blue tilapia Oreochromis aureus
489 (Steindachner 1864) [12,99] are described as not heritable.
490 The situation in the Japanese rice fish, Oryzias latipes (Temminck & Schlegel 1846) is
491 especially informative [100-102]. The shortened vertebral column and vertebral fusions are found
492 in (1) fused mutants obtained from the laboratory strain (‘a simple recessive Mendelian character’
493 with expression modified by temperature [100,102]), (2) ‘wild-fused’ phenotype from ‘a certain
494 region on the eastern outskirts of the City of Nagoya’ that were not heritable, and (3) individuals
495 of the normal strain treated by phenylthiourea at the early embryonic stage [101]. In general, it is
496 obvious that in different species and even in different lineages of the same species, the emergence
497 of fish with shortened columns and fused vertebrae may be determined by various external and
498 genetic cues.
499 Considering the possible role of external cues, the following points must be mentioned.
500 First, among the middle Genale Labeobarbus assemblage, the deformation of the vertebral column
501 is found almost exclusively within the isolated genetic pool including the GN, LP, and SH morphs
502 but excluding the trophically specialized morphs. If the deformation has no genetic background
503 and is triggered by external cues such as mechanic, thermal, chemical or radioactive environmental
504 stress, these impacts are only acting selectively. For example, they would have temporal or spatial
20 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
505 effects on only some of the GN and LP barbs, but not the rest of the assemblage that is possible,
506 in our understanding, mostly on the spawning grounds.
507 Second, environmental stress usually results in multiple malformations of different
508 morphological structures [71,103]. We did not observe such multiple morphological abnormalities
509 in SH. Third, external cues like bacterial [104] or parasitic [105,106] infection, and nutrient
510 deficiency (e.g., vitamin C or phosphorous [107]) are usually manifested as a decrease in the
511 individuals’ conditions, especially in the growth rate, which is often retarded in poor health.
512 However, the growth rate was exceptionally high in the abnormal individuals (SH) observed here.
513 Taking into account all of the above, we suggest that genetic factors (possibly together with some
514 environmental ones) are the most parsimonious explanation of the emergence and 10-year
515 persistence of the SH morph in the middle Genale Labeobarbus assemblage.
516 Nevertheless, to definitively prove this suggestion, we must seek answers to the following
517 questions: do SH individuals reproduce in nature? Would be the progeny of artificial crossing be
518 different for the SH breeders vs. normal GN and LP breeders? Are there any genetic/genomic
519 differences between the SH morph and other morphs? Hopefully, these will be addressed in further
520 studies.
521
522 Shortened vertebral column as a target of natural selection
523 Among several streams in the State of Washington, Patten [108] noted that ‘the presence
524 of many severely deformed fish (some of a mature size) shows that some factor or factors permitted
525 their survival’. Möller [109] also stated that ‘anatomical abnormalities … may become particularly
526 evident in populations that do not suffer from a high selection pressure by predators or food
527 competitors’. Essentially these authors are referring to a relaxation of stabilizing selection; an
528 influence of natural selection is often reasonable to expect, but very difficult to prove
529 unambiguously. There are many records on the increased frequency and severity of skeletal
530 abnormalities in artificially reared fish stocks compared to natural populations [110-117].
21 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
531 Interestingly, the opposite situations are also reported [117,118]. Demonstrating how the intensity
532 of rearing conditions plays a role in increasing the frequency and severity of skeletal abnormalities
533 for several fish species [117,119-122] can lend support to the role of selection relaxation, both in
534 cultivated stocks and natural populations exhibiting the mass abnormalities.
535 Some researchers [123-125] have indicated the relaxation of natural selection as an
536 important mechanism in increasing variation of morphological, ecological and other traits in the
537 course of adaptive radiation. Such an increase often allows the radiating species to occupy new
538 adaptive zones. This phenomenon has been called ‘extralimital specialization’ by Myers [126]. In
539 this perspective, the occurrence of the vertebral deformations at a relatively high frequency only
540 in the radiating Labeobarbus assemblage in the middle Genale River may be evidence for relaxed
541 selection in this particular site.
542 Concerning the possibility of directional selection – positive or negative – in shaping the
543 abnormal barbs (SH morph) in the middle Genale River, the following points must be mentioned.
544 First, some retardation of growth in the 5-6 year classes of the abnormal barbs (Fig. 8), along with
545 their relatively short life span (< 7 years), could be treated as evidence for negative selection at
546 later stages of ontogenesis. For Atlantic cod, Gadus morhua L. 1758 from the German Wadden
547 Sea, a few individuals with spinal compression were found in the larger size classes [127]. This
548 reduced frequency in later life stages was explained by negative selection for this phenotype. At
549 the same time, seasonal variation of spinal compression frequencies in the smaller size classes was
550 explained by the migratory activity that differed between the abnormal and normal individuals
551 [127], where the former are considered to be less active.
552 Second, one may expect that the barbs with deformed vertebral columns should exhibit
553 reduction in swimming performance, but this is not straightforward. Studying the triploid lines of
554 Atlantic salmon, Powell et al. [128] did not find a difference in swimming performance between
555 normal fish and those with visible spinal shortening. They did note a poor recovery from
556 exhaustive swimming exercises in the deformed fish compared to the normal fish. Based on this
22 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
557 result, we cannot exclude the possibility of variable swimming performances between the
558 deformed and normal Genale barbs. At the same time, it is clear that the deformed barbs can avoid
559 habitats with fast water current in order to reduce the energy cost of their abnormality.
560 Third, there could be various reasons for the difference in distribution of the deformed
561 vertebrae along the vertebral column between barbs with few and many compressed and/or fused
562 vertebrae (Fig. 7). As described earlier, most of these vertebrae in the individuals with slightly
563 deformed vertebral columns occur in the posterior end of the fish [9,83,86,97,121,129,130]. This
564 phenomenon is explained by the morphogenetic influence of the caudal complex, whose normal
565 development includes several vertebral fusions [130]. The increased frequency of abnormalities in
566 the caudal region was also detected as a result of thyroid hormone disruption [131,132]. However,
567 in the deformed (SH) barbs, the frequency of deformed vertebrae was lower in the pre-caudal
568 region (Fig. 7). This is unlike the farmed stock of Atlantic salmon [86]. In the wild Genale
569 Labeobarbus, this frequency distribution of deformed vertebrae may be caused by the genetic
570 distinctiveness [91,92] of these deformed individuals, or by selection against the individuals with
571 deformed vertebrae in the caudal region, for example because of their potentially decreased
572 swimming capability.
573 The data obtained from individuals with normal vertebral columns allows us to shed light
574 on the susceptibility of column deformations to directional natural selection. First, comparison of
575 the samples from 2009 and 2019 unexpectedly revealed changes in the body form of the GN morph
576 (Fig. 6). The GN individuals became more similar to the SH morph, particularly in the relative
577 body height. At the same time, the SH morph also became more high-bodied than in 2009. The
578 most parsimonious explanation for this phenomenon is a selective pressure for the increase in
579 relative body height in the population studied. Second, in our small sample from the Welemele
580 River (the northern tributary of the Genale; Fig. 2, site no. 3), together with the normal generalized
581 Labeobarbus morph, we found a few individuals with body proportions very similar to that in the
582 SH morph from the middle Genale assemblage. However, their vertebral columns were not
23 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
583 deformed (Fig. 9). This finding might provide evidence for a selective pressure for an increase in
584 the relative body height in another Labeobarbus population.
585 The increased body height is considered as an anti-predator adaptation in the crucian carp
586 Carassius carassius (L. 1758) [133,134]. Indeed, the differentiation among Labeobarbus in the
587 Chamo-Abaya lake basin (Ethiopian Rift valley) corroborates this idea. The lakes are inhabited by
588 the high-bodied barbs classified as L. bynni (Fabricius 1775) [135], while the tributaries are
589 inhabited by the low-bodied barbs classified as L. intermedius (Rüppell 1835) [135]. However,
590 phylogeny based on mt-DNA [30] demonstrates a close relationship between these two forms and
591 their distant relation to L. bynni from the Nile basin. Thus, high-bodied forms evolved
592 independently in the Chamo-Abaya and Nile basins most probably as an anti-predator adaptation,
593 as the faunal data [74,136,137] indicate the greater predator diversity in these basins compared to
594 the Chamo-Abaya tributaries and most other rivers of the Ethiopian Highlands that are populated
595 by the low-bodied Labeobarbus. In our understanding, the importance of the predator-induced
596 high-bodied phenotype is evident in some Labeobarbus populations. In these cases, natural
597 selection could favor the Genale deformed barbs (SH morph) because of their anti-predator high-
598 bodied phenotype.
599
600 Fig. 9. X-ray images of two Labeobarbus from the Welemele River with (A) short-like
601 exterior (SL = 160 mm) and (B) normal exterior (SL = 177 mm). No deformed vertebrae were
602 detected in either individual.
603
604 Shortened vertebral column as a population phenomenon
605 In natural fish populations, the frequency of individuals with few compressed and/or fused
606 vertebrae usually ranges from less than one percent to several percents (usually <10%) [109-112,
607 114,118,138-143], but sometimes it appears substantially higher (up to or over 50%)
608 [115,117,143-145]. Considering these facts, the frequencies of individuals with 1 to 6 compressed
24 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
609 and/or fused vertebrae (ranging from 5% to 15% in the morphs other than SH; Table 1) do not
610 look outstanding. However, among the 34 individuals sampled from the Awata River (Fig. 2, site
611 no. 4), such deformed vertebrae were not found at all. Hence, in the absence of data for other
612 Labeobarbus populations, we may consider the frequency of individuals with deformed vertebrae
613 in the middle Genale populations somewhat increased. Unfortunately, we cannot say with certainty
614 how this is related to an emergence of the SH morph there.
615 There are many reports of individuals with the markedly shortened bodies and deformed
616 vertebral columns in natural populations of different fish groups (S1 Table). Among wild
617 cyprinids, there are reports of a Mesopotamian barb, Mesopotamichthys sharpeyi (Günther 1874),
618 with nine deformed vertebrae, and a yellowfin barbell, Luciobarbus xanthopterus Heckel 1843 ,
619 with 22 deformed vertebrae [146,147]. The situations, as in the Genale Labeobarbus, where the
620 abnormality becomes a population phenomenon occurring with substantial frequency among
621 several successive generations are rare. It is important to note that in the Genale barb assemblage,
622 only two morphs (GN and JB) were always present in catches with frequencies higher than 10%
623 recorded for the abnormal (SH) morph (S7 Table).
624 To the best of our knowledge, apart from the Genale Labeobarbus, the similar situations
625 are described for only two other fish species, Atlantic cod Gadus morhua from the Elbe estuary
626 and German Wadden Sea [127,148-150] and Arctic char Salvelinus alpinus from Transbaikalian
627 Lake Dzhelo [151]. The high prevalence of vertebral column deformities was recently reported in
628 a wild population of lumpfish Cyclopterus lumpus L. 1758 from Masfjorden, Norway [117], but
629 changes of the body form and temporal stability of the abnormal phenotypes are not evident in the
630 latter species.
631 In 1971, Wunder [148] reported the high prevalence (10-15%, sometimes up to 20%) of
632 cod individuals with shortened bodies and deformed vertebrae in catches from the Elbe estuary
633 and adjacent regions of the Northern Sea. In each deformed individual, 40-80% of vertebrae were
634 compressed and/or fused. This abnormality was most frequent in the young cods (length < 20 cm).
25 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
635 The situation persisted at least until the late 1980s [127,149,150]. The geographical range of this
636 cod deformity expanded to some regions of the Wadden Sea; the mean frequency of the deformed
637 individuals in the catches of all inshore regions of the German Wadden Sea was 9.2% in the size
638 class of 11-23 cm [127].
639 Dzhelo is a small lake (1.4×0.3 km, maximum depth 25 m) in the upper reaches of the Lena
640 River system in Transbaikalia. According to Alekseyev [151], the lake is inhabited by dwarf char
641 (maximum length 24 cm). There are numerous deformed individuals with the shortened body,
642 which also show multiple fusions of vertebrae. Their prevalence is 12%, and a change of the body
643 form in the individuals with ≤ 6 deformed vertebrae is not noticeable. The lake was sampled once
644 in 2003 [151].
645 The main question about these deformities in cod and char is whether they still persist with
646 the substantial frequencies. Comparative etiology of vertebral deformities in cyprinid, salmonid
647 and gadid species could be of great interest. Moreover, deeper understanding of the action of
648 natural selection on these abnormalities can shed the light on how morphological novelties are
649 established in a population, especially if the abnormalities in question have the genetic
650 background.
651
652 Conclusions
653 The striking emergence of individuals with deformed (shortened) vertebral columns
654 resulting in a high-bodied phenotype is a phenomenon rarely detected among wild fish
655 populations. For the particular case of the radiating Labeobarbus assemblage from the middle
656 Genale River, the following circumstances must be highlighted: i) the recent emergence (~ 15
657 years ago) of the high-bodied phenotype with the shortened vertebral column, ii) the persistence
658 of this phenotype with a substantial frequency (~ 10%) in several generations (>10 years) in a
659 population, iii) the fact that this is the only type of morphological deformity; virtually all other
660 morphological abnormalities are absent, iv) the presence of the deformity in the isolated gene pool
26 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
661 within the radiating Labeobarbus assemblage, and v) the rapid growth of the abnormal individuals
662 at 2-4 years of life compared to other sympatric morphs. Taking into account all the above, as well
663 as the evidence for a genetic contribution to such abnormality in some other fish species,
664 particularly the common carp, we assume that this phenomenon in this Labeobarbus population
665 most likely has a genetic background. Moreover, the high-bodied phenotype in Labeobarbus may
666 have a selective advantage as an anti-predator defense.
667
668
669 Acknowledgments
670
671 Material for this study was collected within the scope of the Joint Ethiopian-Russian Biological
672 Expedition (JERBE). We gratefully acknowledge the JERBE coordinators A.A. Darkov
673 (Severtsov Institute of Ecology and Evolution, Moscow, Russia - IEE) and Simenew Keskes
674 Melaku (Ministry of Innovation & Technology, Addis Ababa, Ethiopia) for administrative
675 assistance. We express our gratitude to S.E. Cherenkov and Yu.Yu. Dgebuadze (both from IEE),
676 Fekadu Tefera and Genanaw Tesfaye (both from National Fishery and Other Aquatic Life
677 Research Center of the Ethiopian Institute of Agricultural Research - EIAR, Sebeta), M.V. Mina
678 (Institute of Developmental Biology, Moscow, Russia) for sharing field operations and assistance
679 in collecting material. We thank S.E. Cherenkov for his help with photography, F.N. Skil (IEE)
680 for aging materials from experimentally reared barbs, O.N. Artaev for creating the map, M.V.
681 Mina for critical comments on the manuscript, and Jacquelin DeFaveri for linguistic corrections.
682 This work was financially supported by the Russian Science Foundation, project no. 19-14-00218.
683
684
685
686
27 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
687 References
688 1. Darwin C. On the origin of species by means of natural selection, or the preservation of
689 favoured races in the struggle for life. A new ed. New York: D. Appleton & Company;
690 1860. http://darwin-online.org.uk/converted/pdf/1861_OriginNY_F382.pdf
691 2. Nelson JS, Grande T, Wilson MVH. Fishes of the world. 5th ed. Hoboken (NJ): John
692 Wiley & Sons; 2016.
693 https://batrachos.com/sites/default/files/pictures/Books/Nelson_ea_2016_Fishes%20of%
694 20the%20World.pdf
695 3. Fraser-Brunner A. The ocean sunfishes (Family Molidae). Bull Br Mus (Nat Hist) Zool.
696 1951;1:87-121.
697 4. Pope EC, Hays GC, Thys TM, Doyle TK, Sims D, Queiroz N, et al. The biology and
698 ecology of the ocean sunfish Mola mola: a review of current knowledge and future
699 research perspectives. Rev Fish Biol Fisher. 2010;20(4):471-87.
700 https://doi.org/10.1007/s11160-009-9155-9
701 5. Tyler JC. Abnormal fin and vertebral growth structures in plectognath fishes. Proc Acad
702 Nat Sci Phila. 1970;122:249-71.
703 6. Johnson GD, Britz R. Leis’ conundrum: Homology of the clavus of the ocean sunfishes.
704 2. Ontogeny of the median fins and axial skeleton of Ranzania laevis (Teleostei,
705 Tetraodontiformes, Molidae). J Morphol. 2005;266:11–21.
706 https://doi.org/10.1002/jmor.10242
707 7. Barahona-Fernades MH. Body deformities in hatchery reared European sea bass
708 (Dicentrarchus labrax L). Types, prevalence, and effect on fish survival. J Fish Biol.
709 1982;21:239-49.
710 8. Divanach P, Boglione C, Menu B, Koumoundouros G, Kentouri M, Cataudella S. (1996).
711 Abnormalities in finfish mariculture: an overview of the problem, causes and solutions.
712 In: European Aquaculture Society (EAS) International Workshop on Sea Bass and Sea
28 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
713 Bream Culture: Problems and Prospects: Handbook of contributions and short
714 communications; 1996 16–18 Oct; Verona, Italy. 1996. p. 45-66.
715 https://core.ac.uk/download/pdf/80177647.pdf
716 9. Gavaia PJ, Dinis MT, Cancela ML. Osteological development and abnormalities of the
717 vertebral column and caudal skeleton in larval and juvenile stages of hatchery-reared
718 Senegal sole (Solea senegalensis). Aquaculture 2002;211:305-23.
719 https://doi.org/10.1016/S0044-8486(02)00167-9
720 10. vertebral stability Gorman KF, Breden F. Teleosts as models for human and deformity.
721 Comp Biochem Physiology C. 2007;145:28-38.
722 https://doi.org/10.1016/j.cbpc.2006.10.004 PMID: 17240199
723 11. Witten PE, Gil Martens L, Huysseune A, Takle H, Hjelde K. Classification and
724 developmental relationships of vertebral body malformations in farmed Atlantic salmon
725 (Salmo salar). Aquaculture. 2009;295:6-14.
726 https://doi.org/10.1016/j.aquaculture.2009.06.037
727 12. Tave D, Jo J-Y, Kim DS. Gross abnormalities in tilapia. Fish Aquatic Sci. 2011;14:148-
728 60. https://doi.org/10.5657/FAS.2011.0148
729 13. Boglione C, Gisbert E, Gavaia P, Witten PE, Moren M, Fontagné S, et al. (2013).
730 Skeletal anomalies in reared European fish larvae and juveniles. Part 2: main typologies,
731 occurrences and causative factors. Rev Aquacult. 2013;5:121-67.
732 https://doi.org/10.1111/raq.12016
733 14. Dawson CE. A bibliography of anomalies of fishes. Gulf Res Rep. 1964;1(6):308-99.
734 Retrieved from https://aquila.usm.edu/gcr/vol1/iss6/1
735 15. Dawson CE. A bibliography of anomalies of fishes. Supplement 1. Gulf Res. Rep.
736 1966;2(2):169-76. Retrieved from https://aquila.usm.edu/gcr/vol2/iss2/3
737 16. Dawson CE. A bibliography of anomalies of fishes. Supplement 2. Gulf Res. Rep.
738 1971;3(2):215-39. Retrieved from https://aquila.usm.edu/gcr/vol3/iss2/5
29 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
739 17. Dawson CE, Heal E. A bibliography of anomalies of fishes. Supplement 3. Gulf Res.
740 Rep. 1976;5(2):35-41. Retrieved from https://aquila.usm.edu/gcr/vol5/iss2/5
741 18. Wunder W. Wirbelsäulenverkürzung als rassenbildendes Merkmal beim Aischgründer
742 Karpfen. Wilhelm Roux' Arch Entwicklung. 1949;144;1-24.
743 19. Steffens W. Der Karpfen. 5 Auflage. Wittenberg, Lutherstadt: A. Ziemsen Verlag; 1980.
744 20. Oellerman LK, Skelton PH. Hexaploidy in yellowfish species (Barbus, Pisces, Cyprinidae)
745 from southern Africa. J Fish Biol. 1990;37:105-15. https://doi.org/10.1111/j.1095-
746 8649.1990.tb05932.x
747 21. Golubtsov AS, Krysanov EY. Karyological study of some cyprinid species from
748 Ethiopia. The ploidy differences between large and small Barbus of Africa. J Fish Biol.
749 1993;42:445-55. https://doi.org/10.1111/j.1095-8649.1993.tb00347.x
750 22. Guégan J-F, Rab P, Machordom A, Doadrio I. New evidence of hexaploidy in 'large'
751 African Barbus with some considerations on the origin of hexaploidy. J Fish Biol.
752 1995;47:192-8. https://doi.org/10.1111/j.1095-8649.1995.tb01888.x
753 23. Krysanov EYu, Golubtsov AS. Karyotypes of some Ethiopian Barbus and Varicorhinus
754 from the Nile Basin including Lake Tana morphotypes. Folia Zool. 1996;45(Suppl.
755 1):67-75.
756 24. Naran D, Skelton PH, Villet MH. Karyology of three evolutionarily hexaploid southern
757 African species of yellowfish, Labeobarbus Rüppel, 1836 (Cyprinidae). Afr. Zool.
758 2007;42:254-60. https://doi.org/10.3377/1562-7020(2007)42[254:KOTEHS]2.0.CO;2
759 25. Yang L, Sado T, Hirt MV, Pasco‐Viel E, Arunachalam M, Li J, et al. Phylogeny and
760 polyploidy: resolving the classification of cyprinine fishes (Teleostei: Cypriniformes).
761 Mol. Phylogenet. Evol. 2015;85:97-116. https://doi.org/10.1016/j.ympev.2015.01.014
762 PMID: 25698355
30 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
763 26. Vreven E.J, Musschoot T, Snoeks J, Schliewen UK. The African hexaploid Torini
764 (Cypriniformes: Cyprinidae): Review of a tumultuous history. Zool J Linnean Soc.
765 2016;177:231-305. https://doi.org/10.1111/zoj.12366
766 27. Levin BA, Golubtsov AS, Dgebuadze YY, Mugue NS. New evidence of homoplasy
767 within the African genus Varicorhinus (Cyprinidae): an independent origin of specialized
768 scraping forms in the adjacent drainage systems of Ethiopia inferred from mtDNA
769 analysis. Afr. Zool. 2013;48:400-6. https://doi.org/10.1080/15627020.2013.11407609
770 28. Beshera KA, Harris PM. Mitochondrial DNA phylogeography of the Labeobarbus
771 intermedius complex (Pisces, Cyprinidae) from Ethiopia. J Fish Biol. 2014;85:228-45.
772 https://doi.org/10.1111/jfb.12408
773 29. Beshera KA, Harris PM, Mayden RL. Novel evolutionary lineages in Labeobarbus
774 (Cypriniformes; Cyprinidae) based on phylogenetic analyses of mtDNA sequences.
775 Zootaxa. 2016;409(3):363-81. http://dx.doi.org/10.11646/zootaxa.4093.3.4
776 30. Levin BA, Golubtsov AS. An evidence of past introgressive hybridization between
777 Labeobarbus ethiopicus and L. intermedius in the Ethiopian Rift Valley, East Africa.
778 Ethiop J Biol Sci. 2017;16(Suppl.):45-60.
779 31. Levin BA, Casal-López M, Simonov E, Dgebuadze YY, Mugue NS, Tiunov AV, et al.
780 (2019). Adaptive radiation of barbs of the genus Labeobarbus (Cyprinidae) in an East
781 African river. Freshw Biol. 2019;64:1721-36. https://doi.org/10.1111/fwb.13364
782 32. Levin BA, Simonov E, Dgebuadze YY, Levina M, Golubtsov AS. (2020). In the rivers:
783 Multiple adaptive radiations of cyprinid fishes (Labeobarbus) in Ethiopian Highlands.
784 Sci Rep. 2020;10(1):7192. https://doi.org/10.1038/s41598-020-64350-4
785 33. Bini G. I pesci del Lago Tana. Missione di studio al Lago Tana richerge Limnologiche.
786 B. Chimica e Biologia. Reale Accademia d'Italia. 1940;3:137-79.
31 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
787 34. Brunelli G. Le mutazioni del genere Barbus del Lago Tana. Missione di studio al Lago
788 Tana richerge Limnologiche. B. Chimica e Biologia. Reale Accademia d'Italia. 1940;3:
789 207-13.
790 35. Nagelkerke LAJ, Sibbing FA, van den Boogaart JGM, Lammens EHRR, Osse JWM. The
791 barbs (Barbus spp.) of Lake Tana: a forgotten species flock? Env Biol Fish. 1994;39:1-
792 21. https://doi.org/10.1007/BF00004751
793 36. Nagelkerke LAJ, Mina MV, Tesfaye Wudneh, Sibbing FA, Osse J.M. In Lake Tana, a
794 unique fish fauna needs protection. BioScience. 1995;45:772-775.
795 https://doi.org/10.1007/BF00004751
796 37. Nagelkerke LAJ, Sibbing FA, Osse JWM. 1995, Morphological divergence during
797 growth in the large barbs (Barbus spp.) of Lake Tana, Ethiopia. Neth J Zool.
798 1995;45:431-54. https://doi.org/10.1163/156854295X00401
799 38. Alekseyev SS, Dgebuadze YY, Mina MV, Mironovsky AN. Small ‘large barbs’
800 spawning in tributaries of Lake Tana: what are they? Folia Zool. 1996;45(suppl.1):85-96.
801 39. Mina M.V., Mironovsky AN, Dgebuadze YY. Lake Tana large barbs: phenetics, growth
802 and diversification. J Fish Biol. 1996;48:383–404. https://doi.org/10.1111/j.1095-
803 8649.1996.tb01435.x
804 40. Mina M.V., Mironovsky AN, Dgebuadze YY. Morphometry of barbel of Lake Tana,
805 Ethiopia: multivariate ontogenetic channels. Folia Zool. 1996;45(suppl. 1):109-16.
806 41. Nagelkerke LAJ, Sibbing FA. Reproductive segregation among the large barbs (Barbus
807 intermedius complex) of Lake Tana, Ethiopia. An example of intralacustrine speciation?
808 J Fish Biol. 1996;49:1244-66. https://doi.org/10.1111/j.1095-8649.1996.tb01793.x
809 42. Nagelkerke LAJ, Sibbing FA. A revision of the large barbs (Barbus spp., Cyprinidae,
810 Teleostei) of Lake Tana, Ethiopia, with a description of seven new species. pp. 105–170.
811 In: Nagelkerke LAJ. The barbs of Lake Tana, Ethiopia: Morphological diversity and its
812 implications for taxonomy, trophic resource partitioning and fisheries [Ph.D. Thesis].
32 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
813 Wageningen (The Netherlands): Wageningen Agricultural University; 1997. p. 105-70.
814 https://edepot.wur.nl/121222
815 43. Berrebi P, Valushok D. Genetic divergence among morphotypes of Lake Tana (Ethiopia)
816 barbs. Biol J Linn Soc. 1998;64:369-84. https://doi.org/10.1111/j.1095-
817 8312.1998.tb00338.x
818 44. Sibbing FA, Nagelkerke LAJ, Stet RJM, Osse JWM. Speciation of endemic Lake Tana
819 barbs (Cyprinidae, Ethiopia) driven by trophic resource partitioning; a molecular and
820 ecomorphological approach. Aquatic Ecol. 1998;32:217-27.
821 https://doi.org/10.1023/A:1009920522235
822 45. Dgebuadze YY, Mina MV, Alekseyev SS, Golubtsov AS. Observations on reproduction
823 of the Lake Tana barbs. J Fish Biol. 1999;54:417-23. https://doi.org/10.1111/j.1095-
824 8649.1999.tb00840.x
825 46. De Graaf M, Eshete Dejen, Sibbing FA, Osse JWM. The piscivorous barbs of Lake Tana
826 (Ethiopia): major questions on their evolution and exploitation. Neth J Zool.
827 2000;50:215-23. https://doi.org/10.1163/156854200X00081
828 47. Nagelkerke LAJ, Sibbing FA. 2000. The large barbs (Barbus spp., Cyprinidae, Teleostei)
829 of Lake Tana (Ethiopia), with a description of a new species, Barbus osseensis. Neth J
830 Zool. 2000;50:179-214. https://doi.org/10.1163/156854200X00072
831 48. Mina MV, Mironovsky AN, Golani D. Consequences and modes of morphological
832 diversification of East African and Eurasian barbins (genera Barbus, Varicorhinus and
833 Capoeta) with particular reference to Barbus intermedius complex. Env Biol Fish.
834 2001;61: 241-52. https://doi.org/10.1023/A:1010952108081
835 49. Sibbing FA, Nagelkerke LAJ. Resource partitioning by Lake Tana barbs predicted from
836 fish morphometrics and prey characteristics. Rev Fish Biol Fisher. 2001;10:393-437.
837 https://doi.org/10.1023/A:1012270422092
33 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
838 50. Mina MV, Golani D. Changes of skull proportions in ontogeny of Barbus canis from
839 Lake Kinneret and of the Lake Tana large-mouthed barbs, Barbus intermedius complex.
840 Env Biol Fish. 2004;70:315-20. https://doi.org/10.1023/B:EBFI.0000033346.86973.17
841 51. de Graaf M, Nentwich ED, Osse JWM, Sibbing FA. (2005). Lacustrine spawning : is this
842 a new reproductive strategy among ‘large’ African cyprinid fishes? J Fish Biol.
843 2005;66:1214-36. https://doi.org/10.1111/j.0022-1112.2005.00671.x
844 52. Smirnov SV, Dzerzhinskii KF, Levin BA. (2006). On the relationship between scale
845 number in the lateral line in the African barbell Barbus intermedius and the rate of
846 ontogeny (by experimental data). J Appl Ichthyol. 2006;46:129-32.
847 https://doi.org/10.1134/S0032945206010176
848 53. Zworykin DD, Budaev SV, Darkov AA, Dzerzhinskii KF, Lyovin BA, Mina MV.
849 Assessment of the role of chemoreception in the mate choice in barbs of the Barbus
850 intermedius complex from Lake Tana, Ethiopia. J Ichthyol. 2006;46:661-7.
851 https://doi.org/10.1134/S0032945206080133
852 54. de Graaf M, Megens H-J, Samallo J, Sibbing FA. Evolutionary origin of Lake Tana’s
853 (Ethiopia) small Barbus species: indications of rapid ecological divergence and
854 speciation. Anim Biol. 2007;57:39-48. https://doi.org/10.1163/157075607780002069
855 55. Dzerzhinskii KF, Shkil FN, Belay Abdissa, Wondie Zelalem, Mina MV. Spawning of
856 large Barbus (Barbus intermedins complex) in a small river of the Lake Tana Basin
857 (Ethiopia) and relationships of some putative species. J Ichthyol. 2007;47:639-46.
858 https://doi.org/10.1134/S0032945207080103
859 56. de Graaf M, Eshete Dejen, Osse JWM, Sibbing FA. 2008. Adaptive radiation of Lake
860 Tana’s (Ethiopia) Labeobarbus species flock (Pisces, Cyprinidae). Mar Freshwater Res.
861 2008;59:391-407. https://doi.org/10.1071/MF07123
862 57. de Graaf M, Megens H-J, Samallo J, Sibbing F. Preliminary insight into the age and
863 origin of the Labeobarbus fish species flock from Lake Tana (Ethiopia) using the
34 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
864 mtDNA cytochrome b gene. Mol Phylogenet Evol. 2010;54:336-43.
865 https://doi.org/10.1016/j.ympev.2009.10.029
866 58. de Graaf M, Nagelkerke LAJ, Palstra AP, Sibbing FA. 2012. Experimental evidence for
867 the biological species status in Lake Tana's Labeobarbus flock (Cyprinidae). Anim Biol.
868 2012;60:183-93. https://doi.org/10.1163/157075610X491725
869 59. Shkil FN, Kapitanova DV, Borisov VB, Abdissa B, Smirnov SV. Thyroid hormone in
870 skeletal development of cyprinids: Effects and morphological consequences. J Appl
871 Ichthyol. 2012;28:398-405. https://doi.org/10.1111/j.1439-0426.2012.01992.x
872 60. Smirnov SV, Borisov VB, Kapitanova DV, Abdissa B, Shkil FN. Heterochronies in skull
873 development of Lake Tana large African barbs (Labeobarbus; Cyprinidae; Teleostei). J
874 Appl Ichthyol. 2012;28:406-10. https://doi.org/10.1111/j.1439-0426.2012.01993.x
875 61. Nagelkerke LAJ, Leon-Kloosterziel KM, Megens H-J, de Graaf M, Diekmann OE,
876 Sibbing FA. Shallow genetic divergence and species delineations in the endemic
877 Labeobarbus species flock of Lake Tana, Ethiopia. J Fish Biol. 2015;87:1191-208.
878 62. Shkil FN, Lazebnyi OE, Kapitanova DV, Abdissa B, Borisov VB, Smirnov SV.
879 Ontogenetic mechanisms of explosive morphological divergence in the Lake Tana
880 (Ethiopia) species flock of large African barbs (Labeobarbus; Cyprinidae; Teleostei).
881 Russ J Dev Biol. 2015;46:294-306. https://doi.org/10.1134/S1062360415050069
882 63. Shkil FN, Smirnov SV, Kapitanova DV. Internal sources of infraorbital series and lateral
883 line scale count variability in six Labeobarbus species composing the Lake Tana
884 (Ethiopia) species flock. J Appl Ichthyol. 2018;34:465-71.
885 https://doi.org/10.1111/jai.13665
886 64. Vreven EJ, Musschoot T, Decru E, Wamuini Lunkayilakio S, Obiero K, Cerwenka AF, et
887 al. (2019). The complex origins of mouth polymorphism in the Labeobarbus
888 (Cypriniformes: Cyprinidae) of the Inkisi River basin (Lower Congo, DRC, Africa):
35 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
889 Insights from an integrative approach. Zool J Linnean Soc. 2019;186:414-82.
890 https://doi.org/10.1163/157075610X491725
891 65. Golubtsov AS. Biogéographie des ‘grands Barbus’ d'Éthiopie avec référence spéciale à
892 des formes à statuts taxinomiques incertains. Cah Ethol. 1993;13:227-30.
893 66. Mina MV, Mironovsky AN, Golubtsov AS, Dgebuadze YuYu. 'Barbus' intermedius
894 species flock in Lake Tana (Ethiopia). II. Morphological diversity of "large barbs" from
895 Lake Tana and neighbouring areas: homoplasies or synapomorphies? Ital J Zool.
896 1998;65(Suppl.):9-14. https://doi.org/10.1080/11250009809386789
897 67. Dimmick WW, Berendzen PB, Golubtsov AS. Genetic comparison of three Barbus
898 (Cyprinidae) morphotypes from the Genale River, Ethiopia. Copeia. 2001(4):1123-29.
899 https://doi.org/10.1643/0045-8511(2001)001[1123:GCOTBC]2.0.CO;2
900 68. Paugy G. The Ethiopian subregion fish fauna: an original patchwork with various
901 origins. Hydrobiologia. 2010;649:301–15. https://doi.org/10.1007/s10750-010-0273-0
902 69. Berhanu B, Seleshi Y, Melesse AM. (2014). Chapter 6. Surface water and groundwater
903 resources of Ethiopia: Potentials and challenges of water resources development. In:
904 Melesse AM, Abtew W, Setegn SG, editors. Nile River Basin: Ecohydrological
905 Challenges, Climate Change and Hydropolitics. Switzerland: Springer International
906 Publishing; 2014. p. 97-213. https://doi.org/10.1007/978-3-319-02720-3_6
907 70. Naseka AM. Comparative study on the vertebral column in the Gobioninae (Cyprinidae,
908 Pisces) with special reference to its systematics. Publ Espec Inst Esp Oceanogr.
909 1996;21:149-67.
910 71. Bogutskaya, NG, Zuykov MA, Naseka AM, Anderson B. Normal axial skeleton structure
911 in common roach Rutilus rutilus (Actinopterygii: Cyprinidae) and malformations due to
912 radiation contamination in the area of the Mayak (Chelyabinsk Province, Russia) nuclear
913 plant. J Fish Biol. 2011;79:991-1016. https://doi.org/10.1111/j.1095-8649.2011.03078.x
36 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
914 72. Mina MV, Levin BA, Mironovsky AN. On the possibility of using character estimates
915 obtained by different operators in morphometric studies of fish. J Ichthyol. 2005;45:284-
916 94. https://doi.org/10.1134/S0032945206020056
917 73. RStudio Team (2020). RStudio: Integrated Development for R. RStudio, PBC, Boston,
918 MA URL http://www.rstudio.com/.
919 74. Golubtsov AS, Dgebuadze YY, Mina MV. Fishes of the Ethiopian Rift Valley. In:
920 Tudorancea C, Taylor WD, editors. Ethiopian Rift Valley Lakes (pp. 167–258). Leiden,
921 The Netherlands: Backhuys Publishers; 2002. p. 167-258.
922 75. Levin BA. New data on morphology of the African scraping feeder Varicorhinus beso
923 (Osteichthyes: Cyprinidae) with the special reference to specialized traits. J Ichthyol.
924 2012;52:908-23. https://doi.org/10.1134/S0032945212110069
925 76. McKay LR, Gjerde B. Genetic variation for a spinal deformity in Atlantic salmon, Salmo
926 salar. Aquaculture. 1986;52:263-72. https://doi.org/10.1016/0044-8486(86)90369-8
927 77. Vågsholm I, Djupvik HO. Risk factors for spinal deformities in Atlantic salmon, Salmo
928 salar L. J Fish Dis. 1998;21:47-53. https://doi.org/10.1046/j.1365-2761.1998.00069.x
929 PMID: 29739172
930 78. Kvellestad A, Høie S, Thorud K, Tørud B, Lyngøy A. Plathyspondyly and shortness of
931 vertebral column in farmed Atlantic salmon Salmo salar in Norway—description and
932 interpretation of pathologic changes. Dis Aquat Organ. 2000;39:97-108.
933 https://doi.org/10.3354/dao039097 PMID: 10715815
934 79. Gjerde B, Pante MJR, Baeverfjord G. Genetic variation for a vertebral deformity in
935 Atlantic salmon (Salmo salar). Aquaculture.2005;244:77-87.
936 https://doi.org/10.1016/j.aquaculture.2004.12.002
937 80. Wargelius A, Fjelldal PG, Hansen T. Heat shock during early somitogenesis induces
938 caudal vertebral column defects in Atlantic salmon (Salmo salar). Dev Genes Evol.
939 2005;215:350-57. https://doi.org/10.1007/s00427-005-0482-0
37 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
940 81. Witten PE, Gil Martens L, Hall BK, Huysseune A, Obach A. Compressed vertebrae in
941 Atlantic salmon (Salmo salar): evidence for metaplastic chondrogenesis as a skeletogenic
942 response late in ontogeny. Dis Aquat Organ. 2005;64:237-46.
943 http://dx.doi.org/10.3354/dao064237 PMID: 15997822
944 82. Fjelldal PG, Lock E-J, Grotmol S, Totland GK, Nordgarden U, Flik G, et al. Impact of
945 smolt production strategy on vertebral growth and mineralisation during smoltification
946 and the early seawater phase in Atlantic salmon (Salmo salar L.). Aquaculture.
947 2006;261:715-28. https://doi.org/10.1016/j.aquaculture.2006.08.008
948 83. Witten PE, Obach A, Huysseune A, Baeverfjord G. Vertebrae fusion in Atlantic salmon
949 (Salmo salar): development, aggravation and pathways of containment. Aquaculture.
950 2006;258:164-72. https://doi.org/10.1016/j.aquaculture.2006.05.005
951 84. Sullivan M, Guy DR, Roberts RJ, Manchester NJ. The aetiology of spinal deformity in
952 Atlantic salmon, Salmo salar L.: influence of genetic factors on the frequency and
953 severity in freshwater stages. J Fish Dis. 2007;30:753-58. https://doi.org/10.1111/j.1365-
954 2761.2007.00888.x
955 85. Fjelldal PG, Glover KA, Skaala O, Imsland A, Hansen TJ. Vertebral body mineralization
956 and deformities in cultured Atlantic salmon (Salmo salar L.): effects of genetics and off-
957 season smolt production. Aquaculture. 2009;296:36–44.
958 https://doi.org/10.1016/j.aquaculture.2009.08.020
959 86. Hansen T, Fjelldal PG, Yurtseva A, Berg A. A possible relation between growth and
960 number of deformed vertebrae in Atlantic salmon (Salmo salar L.). J Appl Ichthyol.
961 2010;26:355-59. https://doi.org/10.1111/j.1439-0426.2010.01434.x
962 87. Fjelldal PG, Hansen T, Breck O, Ørnsrud R, Lock E-J, Waagbø R, et al. Vertebral
963 deformities in farmed Atlantic salmon (Salmo salar L.) — etiology and pathology. J Appl
964 Ichthyol. 2012;28:433-40. https://doi.org/10.1111/j.1439-0426.2012.01980.x
38 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
965 88. Fjelldal PG, Hansen T, Breck O, Sandvik R, Waagbø R, Berg A, et al. Supplementation
966 of dietary minerals during the early seawater phase increase vertebral strength and reduce
967 the prevalence of vertebral deformities in fast-growing under-yearling Atlantic salmon
968 (Salmo salar L.) smolt. Aquacult Nut. 2009;15:366-78. https://doi.org/10.1111/j.1365-
969 2095.2008.00601.x
970 89. Ritchie J. A hump-backed trout from Stranraer. Ann Scott Nat Hist. 1908;68:223-7.
971 90. Poynton S. (1987). Vertebral column abnormalities in brown trout, Salmo trutta L. J
972 Fish Dis. 1987:10:53-7. https://doi.org/10.1111/j.1365-2761.1987.tb00718.x
973 91. Campbell WB. (1995). Genetic variation of vertebral fusion patterns in coho salmon. J
974 Fish Biol. 1995;46:717-20. https://doi.org/10.1111/j.1095-8649.1995.tb01106.x
975 92. Fraser TWK, Hansen TJ, Sambraus F, Fjelldal PG. Vertebral deformities in interspecific
976 diploid and triploid salmonid hybrids. J Fish Biol. 2020.
977 https://doi.org/10.1111/jfb.14353
978 93. Thomerson JE. Tbe stubby vertebral deformity in aquarium-raised Fundulus cingulatus.
979 Copeia. 1966;(4):878-79. https://www.jstor.org/stable/i262904
980 94. Schultz RJ. Stubby, a hereditary vertebral deformity in the viviparous fish Poeciliopsis
981 prolifica. Copeia. 1963;(2):325-30. https://www.jstor.org/stable/1441351
982 95. Lodi E. the guppy Palla, A hereditary vertebral deformity in, Poecilia reticulata Peters,
983 Pisces, Osteichthyes. Genetica. 1978;48:197-200. https://doi.org/10.1007/BF00155570
984 96. Gistelinck C, Kwon RY, Malfait F, Symoens S, Harris MP, Henke K, et al. Zebrafish
985 type I collagen mutants faithfully recapitulate human type I collagenopathies. Proc Natl
986 Acad Sci USA. 2018;115(34);E8037-46. https://doi.org/10.1073/pnas.1722200115
987 PMID: 2197991
988 97. Martini A, Huysseune A, Witten PE, Boglione C. Plasticity of the skeleton and skeletal
989 deformities in zebrafish (Danio rerio) linked to rearing density. J Fish Biol. 2020.
990 https://doi.org/10.1111/jfb.14272
39 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
991 98. Dunhan RA, Smitherman RO, Bondari K. Lack of inheritance of stumpbody and
992 taillessness in channel catfish. Prog Fish-Cult. 1991;53:101-5.
993 https://doi.org/10.1577/1548-8640(1991)053<0101:LOIOSA>2.3.CO;2
994 99. Tave D, Bartels JE, Smitherman RO. (1982). Stumpbody Sarotherodon aureus
995 (Steindachner) (Tilapia aurea) and tail-less S. niloticus (L.) (T. nilotica): two vertebral
996 anomalies and their effects on body length. J Fish Dis. 1982;5:487-94.
997 100. Aida T. Further genetical studies of Aplocheilus latipes. Genetics. 1930;15:1-16.
998 Retrieved from https://www.genetics.org/content/15/1/1
999 101. Yamamoto T, Tomita H, Matsuda N. (1963). Japan J Genet. 1963;38:36-47.
1000 https://doi.org/10.1266/jjg.38.36
1001 102. Ogawa NM. A phylogenetic study of the vertebral fused (f) in the medaka,
1002 Oryzias latipes. Embryologia. 1965;9:13-33. https://doi.org/10.1111/j.1440-
1003 169X.1965.tb00212.x
1004 103. Hubbs C. (1959). High incidence of vertebral deformities in two natural
1005 populations of fishes inhabiting warm springs. Ecology. 1959;40:154-55. DOI:
1006 10.2307/1929941 https://www.jstor.org/stable/1929941
1007 104. Madsen L, Arnbjerg J, Dalsgaard I. Radiological examination of the spinal
1008 column in farmed rainbow trout Oncorhynchus mykiss (Walbaum): experiments with
1009 Flavobacterium psychrophilum and oxytetracycline. Aquac Res. 2001;32:235-41.
1010 https://doi.org/10.1046/j.1365-2109.2001.00552.x
1011 105. Cunningham ME, Markle DF, Watral VG, Kent ML, Curtis LR. Patterns of fish
1012 deformities and their association with trematode cysts in the Willamette River, Oregon.
1013 Env Biol Fish. 2005;73:9-19. https://doi.org/10.1007/s10641-004-3153-5
1014 106. Villeneuve DL, Curtis LR, Jenkins JJ, Warner KE, Tilton F, Kent ML, et al.
1015 Environmental stresses and skeletal deformities in fish from the Willamette River,
1016 Oregon. Environ Sci. 2005;39:3495-506. https://doi.org/10.1021/es048570c
40 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
1017 107. Witten PE, Owen MAG, Fontanillas R, Soenens M, McGurk C, Obach A. A
1018 primary phosphorus-deficient skeletal phenotype in juvenile Atlantic salmon Salmo
1019 salar: the uncoupling of bone formation and mineralization. J Fish Biol. 2015:88:690-
1020 708. https://doi.org/10.1111/jfb.12870
1021 108. Patten BG. Abnormal freshwater fishes in Washington stream. Copeia.
1022 1968;(2):399-401. https://doi.org/10.2307/1441768
1023 109. Möller H. Geographical distribution of fish diseases in the NE Atlantic.
1024 Meeresforschung. 1979;27:217-35.
1025 110. Gabriel, M. L. Factors affeeting the number and form of vertebrae in Fundulus
1026 heteroclitus. J Exp Zool. 1944;95:105-47. https://doi.org/10.1002/jez.1400950108
1027 111. Komada N. Incidence of gross malformations and vertebral anomalies of natural
1028 and hatchery Plecoglossus altivelis. Copeia. 1980;(1):29-35.
1029 https://doi.org/10.2307/1444131
1030 112. Ferreri F, Nicolais C, Boglione C, Bmertoline B. Skeletal characterization of wild
1031 and reared zebrafish: anomalies and meristic characters. J Fish Biol. 2000;56:1115-28.
1032 https://doi.org/10.1111/j.1095-8649.2000.tb02127.x
1033 113. Boglione C, Costa C, Di Dato P, Ferzini G, Scardi M, Cataudella S. Skeletal
1034 quality assessment of reared and wild sharpsnout sea bream and pandora juveniles.
1035 Aquaculture. 2003;227:373-94. https://doi.org/10.1016/s0044-8486(03)00515-5
1036 114. Fjelldal PG, van der Meeren T, Jørstad KE, Hansen TJ. A radiological study on
1037 vertebral deformities in cultured and wild Atlantic cod (Gadus morhua, L.). Aquaculture.
1038 2009;289:6-12. https://doi.org/10.1016/j.aquaculture.2008.12.025
1039 115. Sambraus F, Glover KA, Hansen T, Fraser TWK, Solberg MF, Fjelldal PG.
1040 Vertebra deformities in wild Atlantic salmon caught in the Figgjo River, Southwest
1041 Norway. J Appl Ichthyol. 2014;30:777-82. https://doi.org/10.1111/jai.12517
41 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
1042 116. Lü H, Zang X, Fu M, Dan X, Shenggi S, Weizhi Y. Vertebral deformities in
1043 hatchery-reared and wild-caught juvenile Japanese flounder, Paralichthys olivaceus.
1044 Chin J Oceanol Limn. 2015;33:84-91. http://dx.doi.org/10.1007/s00343-015-4041-x
1045 117. Fjelldal PG, Madaro A, Hvas M., Stien HL, Oppedal F, Fraser TWK. Skeletal
1046 deformities in wild and farmed cleaner fish species used in Atlantic salmon Salmo salar
1047 aquaculture. J Fish Biol. 2020. https://doi.org/10.1111/jfb.14337
1048 118. McCrimmon HR, Bidgood B. Abnormal vertebrae in the rainbow trout with
1049 particular reference to electrofishing. Trans Am Fish Soc. 1965;94:84-8.
1050 https://doi.org/10.1577/1548-8659(1965)94[84:AVITRT]2.0.CO;2
1051 119. Boglione C, Marino G, Giganti M, Longobardi A, De Marzi P, Cataudella S.
1052 Skeletal anomalies in dusky grouper Epinephelus marginatus (Lowe 1834) juveniles
1053 reared with different methodologies and larval densities. Aquaculture. 2009;291:48-60.
1054 https://doi.org/10.1016/j.aquaculture.2009.02.041
1055 120. Roo J, Socorro J, Izquierdo MS. Effect of rearing techniques on skeletal
1056 deformities and osteological development in red porgy Pagrus pagrus (Linnaeus, 1758)
1057 larvae. J Appl Ichthyol. 2010;26:372-76. https://doi.org/10.1111/j.1439-
1058 0426.2010.01437.x
1059 121. Prestinicola L, Boglione C, Makridis P, Spanò A, Rimatori V., Palamara E, et al.
1060 (2013). Environmental conditioning of skeletal anomalies typology and frequency in
1061 gilthead seabream (Sparus aurata L., 1758) juveniles. PLoS ONE. 2013;8(2),e55736.
1062 https://doi.org/10.1371/journal.pone.0055736
1063 122. Boglione C, Pulcini D, Scardi M, Palamara E, Russo T, Cataudella S. (2014).
1064 Skeletal anomaly monitoring in rainbow trout (Oncorhynchus mykiss, Walbaum 1792)
1065 reared under different conditions. PLoS ONE. 2014;9(5),e96983.
1066 https://doi.org/10.1371/journal.pone.0096983
42 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
1067 123. Golubtsov AS. [Species flocks of fishes in rivers and lakes: Sympatric divergence
1068 in faunistically depauperate fish communities as especial modus of evolution]. In: Pavlov
1069 DS, Dgebuadze YY, Shatunovskiy MI, editors. Current problems of ichthyology. To 100
1070 ‐jubelee of G.V. Nikolsky. Moscow, Russia: KMK; 2010. p. 96-123. Russian.
1071 124. Yoder JB, Clancey E, Des Roches S, Eastman JM, Gentry L, et al. Ecological
1072 opportunity and the origin of adaptive radiations. J Evol Biol. 2010;23:1581-96.
1073 https://doi.org/10.1111/j.1420-9101.2010.02029.x
1074 125. Stroud JT, Losos JB. Ecological opportunity and adaptive radiation. Annu Rev
1075 Ecol Evol Syst. 2016;47:507-32. https://doi.org/10.1146/annurev-ecolsys-121415-032254
1076 126. Myers GS. The endemic fish fauna of Lake Lanao and evolution of higher
1077 taxonomic categories. Evolution. 1960;14:323-333. https://doi.org/10.1111/j.1558-
1078 5646.1960.tb03094.x
1079 127. Hilger I. Spinal compression of Atlantic cod Gadus morhua from the German
1080 Wadden Sea. Dis Aquat Organ. 1992;13:83-8. https://doi.org/10.3354/dao013083
1081 128. Powell M, Jones M, Lijalad M. Effects of skeletal deformities on swimming
1082 performance and recovery from exhaustive exercise in triploid Atlantic salmon. Dis
1083 Aquat Organ. 2009;85:59-66. https://doi.org/10.3354/dao02056
1084 129. Bensimon-Brito A, Cancela ML, Huysseune A, Witten PE. The zebrafish (Danio
1085 rerio) caudal complex - a model to study vertebral body fusion. J Appl Ichthyol.
1086 2010;26:235-38. https://doi.org/10.1111/j.1439-0426.2010.01412.x
1087 130. Bensimon-Brito A, Cancela ML, Huysseune A, Witten PE. (2012). Vestiges,
1088 rudiments and fusion events: the zebrafish caudal fin endoskeleton in an evo-devo
1089 perspective. Evol Dev. 2012;14:116-27. https://doi.org/10.1111/j.1525-
1090 142x.2011.00526.x
43 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
1091 131. Shkil FN, Kapitanova DV, Borisov VB, Abdissa B, Smirnov SV. Thyroid
1092 hormone in skeletal development of cyprinids: effects and morphological consequences. J
1093 Appl Ichthyol. 2012;28:398-405. https://doi.org/10.1111/j.1439-0426.2012.01992.x
1094 132. Bolotovskiy AA, Levin BA. Effects of thyroid hormones on vertebral numbers in
1095 two cyprinid fish species: Rutilus rutilus (Linnaeus, 1758) and Abramis brama (Linnaeus,
1096 1758). J Appl Ichthyol. 2018;34:449-54. https://doi.org/10.1111/jai.13659
1097 133. Nilsson PA, Brönmark C. Prey vulnerability to a gape‐size limited predator:
1098 behavioural and morphological impacts on northern pike piscivory. Oikos. 2000;88:539-
1099 46. https://doi.org/10.1034/j.1600-0706.2000.880310.x
1100 134. Domenici P, Turesson H, Brodersen J, Brönmark C. Predator-induced
1101 morphology enhances escape locomotion in crucian carp. P Roy Soc B-Biol Sci.
1102 2007;275:195-201. https://doi.org/10.1098/rspb.2007.1088
1103 135. Banister KE. A revision of the large Barbus (Pisces, Cyprinidae) of East and
1104 Central Africa. Part II. Bull Br Mus Nat Hist (Zool). 1973;26:1-148.
1105 https://doi.org/10.5962/bhl.part.204
1106 136. Golubtsov AS, Darkov AA, Dgebuadze YY, Mina MV. An artificial key to fish
1107 species of the Gambela region (the White Nile basin in the limits of Ethiopia). Addis
1108 Abeba: Joint Ethio-Russian Biological Expedition; 1995.
1109 137. Moritz T, El Dayem ZM, Abdallah MA, Neumann D. New and rare records of
1110 fishes from the White Nile in the Republic of the Sudan. Cybium. 2019;43:137-51.
1111 https://doi.org/10.26028/cybium/2019-423-011
1112 138. Ford E, Bull HO. Abnormal vertebrae in herring. J Mar Biol Ass UK.
1113 1926;14:509-17. https://doi.org/10.1017/S0025315400007979
1114 139. Orska J. Anomalies in the vertebral columns of the pike (Esox lucius L.). Acta
1115 Biol Crac Ser Zool. 1962;5:327-45.
44 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
1116 140. Gill CD, Fisk DM. Vertebral abnormalities in sockeye, pink and chum salmon.
1117 Trans Am Fish Soc. 1966;95:177-82. https://doi.org/10.1577/1548-
1118 8659(1966)95[177:VAISPA]2.0.CO;2
1119 141. Komada N. Occurrence and formation of vertebral anomalies in the cyprinid fish,
1120 Zacco platypus. Japan J Ichthyol. 1983;30:150-7. https://doi.org/10.11369/jji1950.30.150
1121 142. Pohl C. Skeletal deformities and trace metal contents of European smelt, Osmerus
1122 eperlanus, in the Elbe estuary. Meeresforschung. 1990;33:76–89.
1123 143. Ershov PN. The vertebral abnormalities in eelpout Zoarces viviparus (Linnaeus,
1124 1758) (Pisces, Zoarcidae). Proceedings ZIN. 2008;312:74-82.
1125 144. Fraser TWK, Hansen T, Mayer I, Skjæraasen JE, Glover KA, Sambraus F, et al.
1126 The effect of triploidy on vaccine side-effects in Atlantic salmon. Aquaculture.
1127 2014;433:481-90. https://doi.org/10.1016/j.aquaculture.2014.07.009
1128 145. Jawad, L. A., Hansen, T. J., Meier, S., Skjaeraasen, J. E., & Fjelldal, P. G.
1129 Vertebral column deformity in adult wild haddock (Melanogrammus aeglefinus). J Fish
1130 Dis. 2018. https://doi.org/10.1111/jfd.12903
1131 146. Al-Hassan LAJ, Naama AK. More vertebral anomalies in fishes from Euphrates
1132 River, Iraq. Indian J Fish. 1986;33:237-41.
1133 147. Jawad LA, Al-Faisal AJ, Al-Mukhtar M. A case of vertebral coalescence in
1134 Luciobarbus xanthopterus (Heckel, 1843) (Pisces: Cyprinidae) obtained from the lower
1135 reaches of Mesopotamia. Trav Mus Natl Hist Nat Grigore Antipa. 2015;57:127-32.
1136 https://doi.org/10.1515/travmu-2015-0006
1137 148. Wunder W. Missbildungen beim Kabeljau (Gadus morhua) verursacht durch
1138 Wirbelsaulenverkiirzung. Helgol Meeresunters. 1971;22:201-12.
1139 149. Möller, H. High skeletal deformation rates of cod in the Elbe estuary. B Eur Ass
1140 Fish Pat. 1983;3:7-8.
45 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
1141 150. Möller, H. Dynamics of fish diseases in the lower Elbe river. Helgol
1142 Meeresunters. 1984;37:389-413.
1143 151. Alekseyev SS. [Distribution, diversity and diversification of the Arctic char
1144 Salvelinus alpinus (L.) complex (Salmoniformes, Salmonidae) in Siberia] [dissertation].
1145 Moscow (Russia): Moscow State University, Faculty of Biology; 2016. Russian.
1146 Retrieved from
1147 http://www.bio.msu.ru/res/Dissertation/757/DISSERTATION_FILENAME/Alekseyev_d
1148 isser.pdf
1149
1150
1151
1152
1153
1154
1155
1156
1157
1158
1159
1160
1161 Supporting information
1162
1163 S1 Table. Occurrence of aberrant individuals or stocks with shortened vertebral column and
1164 abnormally short and high body in the different groups of fish.
1165 (DOCX)
1166 S2 Table. Complete list of sampling sites.
46 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.
1167 (XLSL)
1168 S3 Figure. Photographs of two habitats at the main sampling (no. 1) of the middle Genale River
1169 – (A) the continuous pool and (B) the upper section of rapids downstream of the pool.
1170 (DOCX)
1171 S4 Table. Size characteristics of the cast and gill net catches for different Labeobarbus forms
1172 from the middle Genale assemblage sampled in 2009 and 2019. Size characteristics of the aged
1173 individuals from the gill net catches of 2019.
1174 (DOCX)
1175 S5 Figure. Variation in lower lip development in short individuals.
1176 (DOCX)
1177 S6 Figure. Two piscivorous barbs with (A) normal and (B) deformed vertebral column.
1178 (DOCX)
1179 S7 Table. Gill net catch composition and percentage of different Labeobarbus morphs in the
1180 middle Genale assemblage in 2009 and 2019.
1181 (DOCX)
1182 S8 Table. Eigenvectors of the 10 most loaded characters from Fig. 4.
1183 (DOCX)
1184 S9 Table. Data on vertebral counts.
1185 (DOCX)
1186 S10 Figure. Supplement to Fig. 7. Distribution of deformed vertebrae along the vertebral column
1187 in individuals with 40 and 42 total vertebrae.
1188 (DOCX)
1189 S11 Figure. Growth of short (SH), generalized (GN) and lipped (LP) forms in the middle Genale
1190 Labeobarbus assemblage sampled in 2009, estimated with standard length (SL).
1191 (DOCX)
1192
47 bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license. bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license. bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license. bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license. bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license. bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license. bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license. bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license. bioRxiv preprint doi: https://doi.org/10.1101/2020.09.11.292763; this version posted September 11, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY 4.0 International license.