Body shape variation between farms of tilapia (Oreochromis sp.) in Colombian Andes using landmark-based geometric morphometrics.
Variation in the shape of the body has been extensively studied in fish and particularly in Cichlids (Clabaut et al., 2007; Kassam et al., 2007; Kerschbaumer & Sturmbauer, 2011). An essential element of these studies is their focus on natural populations in the context of evolutionary biology. By contrast in aquaculture, the differences in morphology between farmed populations are of great concern for selective breeding because the shape is a commercially important trait that contributes to the market value of the product (Colihueque & Araneda, 2014; de Oliveira et al, 2016).
One way to describe the shape quantitatively is by geometric morphometrics (GM) instead of traditional methods based on linear measurements between reference points. Landmark-based GM involves summarizing shape regarding a constellation of discrete anatomical loci, each described by Cartesian coordinates (Webster & Sheets, 2010). Key advantages of GM include: a) emphasis on the complete retention of geometric information throughout the research process, b) much higher statistical power to detect shape differences with sufficient sample sizes, c) localization of the spatial morphological variation, and d) visualization of the shape differences directly as illustrations (Zelditch et al., 2004; Slice, 2007; Klingenberg, 2013).
The morphological plasticity of tilapia has been characterized previously by GM. For example, Lorenz et al. (2014) determined that after an eradication attempt with rotenone, tilapia were deeper in body and head shape than pre-management individuals. Similarly, Ndiwa et al. (2016) found variations in the head, caudal peduncle and anal fin base of Nile tilapia from extreme environmental conditions compared to populations experiencing less extreme conditions. These authors also registered morphological differences between populations with similar genetic background. Also, Firmat et al. (2012) described that invasive populations of Mozambique tilapia exhibited a more elongated body shape, a shorter caudal peduncle and a more expanded anterior region relative to native populations. Concerning shape changes during growth, Fujimura & Okada (2008) assessed the developmental trajectory that leads to the adult lower jaw shape in Nile tilapia and concluded that differences in adult shapes might be due to differences arising early in development. Differences in shape between lines, farms or rearing conditions have been found previously using GM on aquaculture species such as: European sea bass (Dicentrarchus labrax), seabream (Sparus aurata), brown trout (Salmo trutta) and rainbow trout (Oncorhynchus mykiss) (Costa et al., 2010; Vehanen & Huusko, 2011; Pulcini et al., 2013; Fragkoulis et al., 2016). Therefore, we hypothesized that a phenotypic variable species as tilapia would show differences in shape between separate rearing sites.
In this study, morphological variation between 265 fishes from one Nile (Oreochromis niloticus, n = 87) and two red farms (mainly Oreochromis mossambicus x Oreochromis aureus, red 1 n = 89, red 2 n = 89) was analyzed using landmark-based GM (Fig. 1).
All three farms are land-based aquaculture systems. Red 1 is geographically separated from red 2 by 5.45 km, red 1 from the Nile by 120.65 km and, finally, red 2 from the Nile by 116.5 km. The Nile farm is a raceway system of export-oriented production of fillets. This intensive farm uses a high water flow rate and high stocking densities. Red farms are semi-intensive productions oriented to the national market. In the three farms, the water demand is calculated from a mass oxygen balance, and no supplemental oxygen is used. The farms were chosen according to the following criteria: a) their relative high contribution to tilapia production in Antioquia region (Colombia), b) they had their breeders and hatcheries, and c) they had their processing facilities. Although the broodstocks have been kept in these farms for more than five generations, specific sources and pedigree information was not registered by the farmers. In the three farms, fry sex reversal was accomplished by oral administration of 17a-methyltestosterone.
General methods followed those of Kavembe et al. (2016), and the handling procedures followed the section seven of the Aquatic Animal Health Code about the welfare of farmed fish (OIE, 2016). Images of the left side of fish with a scale included were taken after harvesting and before stunning at each one of the processing facilities, using an 18-megapixel EOS 7D digital camera with a 50 mm 1:2.5 lens (Canon USA, Inc.) mounted on a tripod stand. The coordinates of 11 landmarks (Fig. 2) were digitalized in the same order on each image after setting the scale factor using TPSDIG2 v2.30 (Rohlf, 2015). For shape analysis, the data set containing the x-y coordinates was then imported into MORPHOJ v1.06d (Klingenberg, 2011).
In order to translate, rotate and uniformly scale the specimens relative to each other so as to minimize a total sum of squares, a full Procrustes fit and a projection of the data to the tangent space (Dryden & Mardia, 1998) was conducted. Next, an inspection for outliers of the new dataset was performed. A regression analysis with Procrustes coordinates as the dependent variable and log-transformed centroid size as the independent variable with a permutation test against the null hypothesis of independence including 10,000 randomization rounds (Klingenberg, 2016) was carried out to examine the statistical association between size and shape. A canonical variate analysis (CVA) of the covariance matrix of the shape coordinates (Mitteroecker & Bookstein, 2011) was used to assess body shape differences between farmed populations of tilapia.
Shape changes were visualized using a wireframe graph as well as a transformation grid superimposed with their warped outline drawing for each canonical variate. To test the significance of the shape differences between farmed populations, a multivariate analysis of variance (MANOVA) was performed using PAST v3.15 (Hammer et al., 2001). A random permutation of individuals testing the significance of each pair-wise Mahalanobis distance among groups with crossvalidation was performed to assess the accuracy of the morphometric classification, using the module PAD of the CLIC package (Dujardin, 2008).
The regression analysis with the group centered scores of Procrustes coordinates as the dependent variable, and log-transformed centroid size as an independent variable (Fig. 3) showed: a) superimposition of the values of the majority of the three groups of tilapia examined, b) a high range of values of the regression score for a small range of log-transformed centroid sizes for each group of tilapia, c) only 2.28% of shape variation in the tilapias studied covaried with size (P-value < 0.0001 at 10,000 permutations), and d) in the present study, a small range of sizes were examined (no fingerlings or juveniles were measured and the range of values of centroid size in Figure 3 was small).
Therefore, only a short section of the growth trajectory was covered by our data. Consequently, no further size correction was applied. The canonical variate analysis displayed a distinct separation between the three tilapia farms. In this analysis, the first and second axes accounted for 70.7 and 29.3% of the total shape variation respectively (Fig. 4). Also, the CVA indicated significant differences among farms of tilapia (P-value < 0.0001 in all cases). The Nile group in this defined morphospace was entirely separated by the first CV from the red groups, while red 1 and red 2 were mainly separated along the second CV.
As illustrated by the wireframe graphs and transformation grids, individuals from the Nile farm were more elongated and had a more ventral position of the posterior extreme of the orbit and the insertion of pectoral fin than red tilapias. Moreover, the Nile group showed a shorter space between the mouth profile (defined by the landmarks one and two) and the posterior extreme of the orbit, compared to red groups. On the other hand, individuals from farm red 2 were deep bodied and had a smaller head compared to tilapias from farm red 1. Finally, the MANOVA detected significant differences in shape variables among farms of tilapia (Wilk's lambda: 0.03775, F: 56,45; P < 0.001) and the cross-validated classification, correctly reassigned 92 and 100% of the individuals from the farms red 1, red 2 and Nile, respectively.
We found that individuals from different farms of tilapia from Antioquia showed significant differences in body shape. Prior research has identified that both environmental and genetic factors influence the body shape of fishes. For example, Crichigno et al. (2012) achieved plastic induction of body shape of Odontesthes hatchery by manipulation of incubation temperature and diet. In the same fashion, Staszny et al. (2013) reared under identical environmental conditions or different diets two inbred lines of Danio rerio. These authors found that genetic and environmental factors markedly determined the shape of scales. Divanach & Koumoundouros (2014) concluded that developmental temperature significantly affected the position of the bases of some head bones and fins of Sparus aurata juveniles. Environmental factors such as water velocity, depth, rearing density, diet, farming method (i.e., pond or cage) and temperature as well as genetics would affect the shape of fishes under culture conditions (Pakkasmaa & Piironen, 2001; Kause et al., 2003; Ramler et al., 2014). In this way, individuals from the same gene pool but reared in different conditions show differences in body shape (Costa et al., 2010; Vehanen & Huusko, 2011; Fragkoulis et al., 2016). Recently, Montoya-Lopez (unpubl. data) characterized the genetic diversity and population structure of the broodstocks from the three farms analyzed in this study using short tandem repeats. This author found that broodstocks from both red farms belonged to a single genetic cluster.
In contrast, the Nile broodstocks formed a separate cluster. However, an important difference between the two red farms was the presence and number of private alleles, particularly in the farm red 2. The Tables 1 and 2 show the water quality of the three tilapia farms evaluated from Betancur et al. (2016). Dissolved oxygen, oxygen saturation, pH, phosphate, ammoniacal nitrogen, and total solids, show no significant difference between farms. Conversely, the temperature was significantly lower in the Nile farm than in red 1, alkalinity and nitrate was significantly lower in the Nile than in red 2 and dissolved solids were significantly higher in red 1 and red 2 than in the Nile. However, the contribution of both genetic and environmental effects on tilapia shape remains to be experimentally determined.
We found that Nile individuals were more elongated and had a more ventral position of the posterior extreme of the orbit while individuals from the two red farms were separated by differences in body depth and head shape. These findings are similar to previous studies in tilapia by GM, which identified changes in body depth and head shape as the main variable characteristics (Firmat et al., 2012; Lorenz et al., 2014). In like manner, Clabaut et al. (2007) compared specimens from 45 species of Lake Tanganyika cichlids and concluded that the most important differences in body shape between species were related to body length as well as the proportion of sizes of head and caudal peduncle.
Our results provide clear evidence that tilapias from different farms in Colombia display differences in body shape. This fact can be applied to selective breeding programs after establishing the preferences of consumers for the body shape of tilapia because consumer perceptions and public attitudes toward specific characteristics of shape in this species remain unclear in Colombia. Therefore, future work should include experiments such as progeny tests to clarify the influence of genetics, environment and their interaction in the body shape of this species.
This research was funded by Secretaria de Agricultura y Desarrollo Rural de Antioquia, Fondo de Ciencia Tecnologia e Innovacion, Sistema General de Regalias, Agreement 4600000970 SADRA-ASOACUICOLA.
Associacao Cultural e Educacional Brasil (ACEB). 2014. 1o Anuario brasileiro da pesca e aquicultura. Itaipu Binacional, Ministerio da Pesca e Aquicultura, Florianopolis.
Betancur, J.J., Correa-Agudelo, L. & Montoya-Lopez, A.F. 2016. Programa de monitoreo sanitario-ambiental y diagnostico molecular de infecciones bacterianas como herramienta estrategica para la certificacion sanitaria de centros de produccion de trucha arcoiris y tilapia. Informe Final Convenio 4600000970. Asociacion Colombiana de Acuicoltores--ASOACUICOLA, Secretaria de Agricultura y Desarrollo Rural de Antioquia, Fondo de Ciencia Tecnologia e Innovacion, Sistema General de Regalias, Medellin, 59 pp.
Clabaut, C., Bunje, P.M., Salzburger, W. & Meyer, A. 2007. Geometric morphometric analyses provide evidence for the adaptive character of the Tanganyikan cichlid fish radiations. Evolution, 61: 560-578.
Colihueque, N. & Araneda, C. 2014. Appearance traits in fish farming: Progress from classical genetics to genomics, providing insight into current and potential genetic improvement. Frontiers in Genetics, 5: 1-8.
Comision Nacional de Acuacultura y Pesca (CONAPESCA). 2013. Anuario estadistico de acuacultura y pesca 2013. Comision Nacional de Acuacultura y Pesca, Mazatlan, 296 pp.
Costa, C., Vandeputte, M., Antonucci, F., Boglione, C., Menesatti, P., Cenadelli, S., Parati, K., Chavanne, H. & Chatain, B. 2010. Genetic and environmental influences on shape variation in the European sea bass (Dicentrarchus labrax). Biological Journal of the Linnean Society, 101: 427-436.
Crichigno, S.A., Battini, M.A. & Cussac, V.E. 2012. Early morphological variation and induction of phenotypic plasticity in Patagonian pejerrey. Neotropical Ichthyology, 10: 341-348.
de Oliveira, C.A., Ribeiro, R.P., Yoshida, G.M., Kunita, N.M., Rizzato, G.S., de Oliveira, S.N., Dos Santos, A.L. & Nguyen, N.H. 2016. Correlated changes in body shape after five generations of selection to improve growth rate in a breeding program for Nile tilapia Oreochromis niloticus in Brazil. Journal of Applied Genetics, 57: 487-493.
Divanach, P. & Koumoundouros, G. 2014. Thermallyinduced phenotypic plasticity in gilthead seabream Sparus aurata L. (Perciformes, Sparidae). Aquaculture, 432: 383-388.
Dryden, I.L. & Mardia, K.V. 1998. Statistical shape analysis. Wiley, Chichester, 347 pp.
Dujardin, J. 2008. Morphometrics applied to medical entomology. Infection, Genetics and Evolution, 8: 875-890.
Food and Agriculture Organization of the United Nations (FAO). 2014. Contribucion de la pesca y la acuicultura a la seguridad alimentaria y el ingreso familiar en Centroamerica. FAO, Panama, 91 pp.
Food and Agriculture Organization of the United Nations (FAO). 2016. The state of world fisheries and aquaculture 2016, Contributing to food security and nutrition for all. FAO, Rome, 200 pp.
Firmat, C., Schliewen, U.K., Losseau, M. & Alibert, P. 2012. Body shape differentiation at global and local geographic scales in the invasive cichlid Oreochromis mossambicus. Biological Journal of the Linnean Society, 105: 369-381.
Fragkoulis, S., Christou, M., Karo, R., Ritas, C., Tzokas, C., Batargias, C. & Koumoundouros, G. 2016. Scaling of body-shape quality in reared gilthead seabream Sparus aurata L. Consumer preference assessment, wild standard, and variability in reared phenotype. Aquaculture Research, 48: 2402-2410.
Fujimura, K. & Okada, N. 2008. Shaping of the lower jaw bone during growth of Nile tilapia Oreochromis niloticus and a Lake Victoria cichlid Haplochromis chilotes: a geometric morphometric approach. Development Growth & Differentiation, 50: 653-663.
Hammer, D., Harper, T. & Ryan, P. 2001. PAST: Paleontological statistics software package for education and data analysis. Palaeontologia Electronica, 4: 1-9.
Kassam, D., Yamaoka, K., Rusuwa, B. & Hori, M. 2007. The robustness of geometric morphometrics in testing the morphological equivalence hypothesis among cichlid species from East African Great Lakes. Biological Journal of the Linnean Society, 91: 1-9.
Kause, A., Ritola, O., Paananen, T., Eskelinen, U. & Mantysaari, E. 2003. Big and beautiful? Quantitative genetic parameters for the appearance of large rainbow trout. Journal of Fish Biology, 62: 610-622.
Kavembe, G.D., Kautt, A.F., Machado-Schiaffino, G. & Meyer, A. 2016. Eco-morphological differentiation in Lake Magadi tilapia, an extremophile cichlid fish living in hot, alkaline and hypersaline lakes in east Africa. Molecular Ecology, 25: 1610-1625.
Kerschbaumer, M. & Sturmbauer, C. 2011. The utility of geometric morphometrics to elucidate pathways of cichlid fish evolution. International Journal of Evolutionary Biology, 2011: 290-245.
Klingenberg, C.P. 2011. MorphoJ: an integrated software package for geometric morphometrics. Molecular Ecology Resources, 11: 353-357.
Klingenberg, C.P. 2013. Visualizations in geometric morphometrics: how to read and how to make graphs showing shape changes. Hystrix, 24: 15-24.
Klingenberg, C.P. 2016. Size, shape, and form: concepts of allometry in geometric morphometrics. Development Genes and Evolution, 226: 113-137.
Lorenz, O., Smith, P. & Coghill, L. 2014. Condition and morphometric changes in tilapia (Oreochromis sp.) after an eradication attempt in Southern Louisiana. NeoBiota, 20: 49-59.
Mitteroecker, P. & Bookstein, F. 2011. Linear discrimination, ordination, and the visualization of selection gradients in modern morphometrics. Evolutionary Biology, 38: 100-114.
National Marine Fisheries Service (NMFS). 2016. Annual trade data by product, country/association 2016. [https://www.st.nmfs.noaa.gov/pls/webpls/trade_prdc t_cntry_ind.results?qtype=IMP&qyearfrom=2015&q yearto=2016&qprod_name=TILAPIA&qcountry=%2 5&qsort=COUNTRY&qoutput=TABLE]. Reviewed: 20 November 2017.
Ndiwa, T.C., Nyingi, D.W., Claude, J. & Agnese, J.F. 2016. Morphological variations of wild populations of Nile tilapia (Oreochromis niloticus) living in extreme environmental conditions in the Kenyan Rift-Valley. Environmental Biology of Fishes, 99: 473-485.
Pakkasmaa, S. & Piironen, J. 2001. Water velocity shapes juvenile salmonids. Evolutionary Ecology, 14: 721-730.
Pulcini, D., Wheeler, P.A., Cataudella, S., Russo, T. & Thorgaard, G.H. 2013. Domestication shapes morphology in rainbow trout Oncorhynchus mykiss. Journal of Fish Biology, 82: 390-407.
Ramler, D., Mitteroecker, P., Shama, L.N., Wegner, K.M. & Ahnelt, H. 2014. Nonlinear effects of temperature on body form and developmental canalization in the threespine stickleback. Journal of Evolutionary Biology, 27: 497-507.
Rohlf, F.J. 2015. The tps series of software. Hystrix, 26: 1-4.
Slice, D.E. 2007. Geometric morphometrics. Annual Review of Anthropology, 36: 261-281.
Staszny, A., Havas, E., Kovacs, R., Urbanyi, B., Paulovits, G., Bencsik, D. & Csenki, Z. 2013. Impact of environmental and genetic factors on the scale shape of zebrafish, Danio rerio (Hamilton 1822): a geometric morphometric study. Acta Biologica Hungarica, 64: 462-475.
Vehanen, T. & Huusko, A. 2011. Brown trout Salmo trutta expresses different morphometrics due to divergence in the rearing environment. Journal of Fish Biology, 79: 1167-1181.
Webster, M. & Sheets, H.D. 2010. A practical introduction to landmark-based geometric morphometrics. The Paleontological Society Papers, 16: 163-188.
World Organisation for Animal Health (OIE). 2016. Aquatic Animal Health Code. World Organisation for Animal Health, Paris, 285 pp.
Zelditch, M.L., Swiderski, D.L., Sheets, H.D. & Fink, W.L. 2004. Geometric morphometries for biologists : a primer. Elsevier Academic Press, New York, 437 pp.
Corresponding editor: Fernando Vega-Villasante
Received: 28 February 2018; Accepted: 25 June 2018
Andres Montoya-Lopez (1), Cintia Moreno-Arias (2), Ariel Tarazona-Morales (3) Martha Olivera-Angel (1) & James Betancur (4)
(1) Grupo Biogenesis, Facultad de Ciencias Agrarias, Universidad de Antioquia, Medellin, Colombia
(2) Laboratorio de Ictiologia, Unidad de Ecologia y Sistematica (UNESIS), Departamento de Biologia Facultad de Ciencias, Pontificia Universidad Javeriana, Bogota, Colombia
(3) Grupo Biogenesis, Departamento de Produccion Animal, Facultad de Ciencias Agrarias Universidad Nacional de Colombia, Medellin, Colombia
(4) Asociacion Colombiana de Acuicultores (ASOACUICOLA), Medellin, Colombia
Corresponding author: Andres Montoya-Lopez (firstname.lastname@example.org)
Caption: Figure 1. Representative images for each group of tilapia examined. a) Nile, b) red 1, and, c) red 2. Scale bar = 1 cm.
Caption: Figure 2. The position of landmarks used in the present study. 1) the intersection between the upper lip and body outline in the nasal-palatine anterior area, 2) most posterior corner of the maxilla when the mouth is closed, 3) the posterior extreme of the orbit, 4) anterior insertion of the first dorsal spine, 5) posterior insertion of the last dorsal ray, 6) last pore of the lateral line, 7) posterior insertion of the last anal ray, 8) anterior insertion of the first anal spine, 9) anterior insertion of the first pelvic spine, 10) upper insertion of pectoral fin, and 11) most ventral corner of interoperculum-suboperculum joint. Image from farm red 2.
Caption: Figure 3. Group-centered scores of the regression between Procrustes coordinates and log-transformed centroid size.
Caption: Figure 4. Canonical variate analysis (CVA) of the covariance matrix of the shape coordinates. The axis represents the canonical variates 1 and 2 (CV1 and CV2) with their respective graphs for visualizations of shape changes.
Table 1. Water quality of the tilapia farms. Data are shown as mean [+ or -] SD; different letters indicate significant (P < 0.05) differences among the farms. Farm Dissolved oxygen (mg [L.sup.-1]) Inlet Outlet Nile 7.16 [+ or -] 1.0 5.59 [+ or -] 1.49 Red 1 7.5 [+ or -] 1.78 6.02 [+ or -] 2.12 Red 2 5.81 [+ or -] 1.42 4.8 [+ or -] 1.95 Farm Oxygen saturation (%) Inlet Outlet Nile 93.61 [+ or -] 10.17 70.75 [+ or -] 18.50 Red 1 94.83 [+ or -] 25.59 78.41 [+ or -] 28.79 Red 2 72.71 [+ or -] 16.77 60.43 [+ or -] 24.16 Farm Temperature ([degrees]C) pH Nile 24.88 [+ or -] 1.35b 7.39 [+ or -] 0.60 Red 1 27.21 [+ or -] 4.46a 7.54 [+ or -] 0.43 Red 2 26.84 [+ or -] 1.12ab 7.71 [+ or -] 0.71 Table 2. Water quality of the tilapia farms. Data are shown as mean [+ or -] SD; different letters indicate significant (P < 0.05) differences among the farms. Alkalinity Phosphate Farm (mg [L.sup.-1]) (mg [L.sup.-1]) Nile 14.88 [+ or -] 9.17a 0.08 [+ or -] 0.04 Red 1 54.38 [+ or -] 20.76a 0.10 [+ or -] 0 Red 2 101.63 [+ or -] 9.23b 0.10 [+ or -] 0 Ammoniacal Nitrate Farm nitrogen (mg (mg [L.sup.-1]) [L.sup.-1]) Nile 0.70 [+ or -] 0.84 1.50 [+ or -] 0.58a Red 1 1.14 [+ or -] 1.92 3.05 [+ or -] 1.32ab Red 2 1.17 [+ or -] 2.02 4.10 [+ or -] 1.82b Total solids Dissolved Farm (mg [L.sup.-1]) solids (mg [L.sup.-1]) Nile 91.02 [+ or -] 46.38 55.13 [+ or -] 27.07a Red 1 187.63 [+ or -] 59.81 132.88 [+ or -] 35.38b Red 2 213.63 [+ or -] 13.46 145.01 [+ or -] 22.80b
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||Short Communication|
|Author:||Montoya-Lopez, Andres; Moreno-Arias, Cintia; Tarazona-Morales, Ariel; Olivera-Angel, Martha; Betancu|
|Publication:||Latin American Journal of Aquatic Research|
|Date:||Mar 1, 2019|
|Previous Article:||Southernmost record of the white shark Carcharodon carcharias (Chondrichthyes: Lamnidae) in the Mexican Pacific.|
|Next Article:||Breaking strength evaluation of biodegradable twines to reduce ghost fishing in the pot and trap fisheries of Chile.|