Barilius bendelisis, commonly known as Indian Hill Trout is an upland water fish of South East Asia. It belongs to the family Cyprinidae and dwells in shallow, clear and cold water. In this study, the intraspecific variation of Barilius bendelisis, on the basis of morphometric characters, was investigated. Altogether, 402 specimens were collected from four rivers in the Central Indian Himalaya. A truss network was constructed by interconnecting 12 landmarks to yield 30 distance variables that were extracted from digital images of specimens using tpsDig2 and PAST software. Allometric transformed truss measurements were subjected to univariate analysis of variance, factor analysis and discriminant analysis. All variables exhibited significant differences between the populations. Altogether 88% of the specimens were classified into their original populations (81.98% under a ‘leave-one-out’ procedure). With factor analysis measurements of the head region, the middle portion and the caudal region had high loadings on the first and second axis. The results indicated that B. bendelisis has significant phenotypic heterogeneity between the geographically isolated regions of Central Indian Himalaya. We hypothesize that the marked interspecific variation in B. bendelisis is the result of local ecological conditions.
Key words: discriminant function analysis, shape, truss box, fishes, Cyprinidae.
Resumen
Barilius bendelisis, comúnmente ]]>
Barilius bendelisis basado en caracteres morfométricos. En total se obtuvieron 402 muestras de cuatro ríos en el centro del Himalaya hindú. Se construyó ]]>
B. bendelisis tiene heterogeneidad fenotípica significativa entre las regiones geográficamente aisladas del centro del Himalaya hindú. Nuestra hipótesis es que hay una marcada variación interespecífica de B. bendelisis como resultado de las condiciones ecológicas locales. Palabras clave: función de análisis discriminante, forma, caja truss, peces, Cyprinidae.
Morphological characters are commonly used in fisheries and systematics to measure discreteness and relationships among various taxa and have long been used to delineate stocks of fish. Morphometry are continuous characters describing the body shape, which have provided evidence for stock discreteness and for species identification (Murta, 2000; Costa, Almeidam, & Costa, ]]>
However, there is a major limitation in using morphological characters at the intra specific ]]>
Barilius bendelisis (Hamilton), commonly known as Indian Hill Trout is an upland water fish of South East Asia. It belongs to the family Cyprinidae and dwells in shallow, clear and cold water. In Northeastern India, it is commonly distributed in hilly streams and rivers of Himalayan region. The fish plays significant role in the capture fishery in several parts of the Himalayan region of Uttarakhand, inhabiting shallow lotic and seasonal lentic water bodies where Indian major carps and exotic carps cannot be raised successfully (Sahoo, Saikia, & Das, 2009). This species is ]]>
B. bendelisis is a popular ornamental species among aquarists as well as a food fish for local population, but is poorly studied; as a result this species has suffered a drastic reduction in the recent past (Sah, Barat, Pande, Sati, & Goel, 2011). Therefore, the present study is a first step to determine the stock structure of this species based on morphometry, ]]>
Material and Methods
Study area: The Indian Himalayas are drained by 19 main rivers, including three major river systems; the Indus, Ganga and Brahmaputra. The Indus Basin system is the longest river system which originates from Western Indian Himalayas (160 000km2) and consists of five rivers. The Ganga basin system consists of nine rivers and originates from Central Indian Himalayas (150 000km2) and the Brahmaputra Basin is the second longest river system which starts in Eastern Himalayas (150 000km2) having five rivers (Hora, 1954). This study includes four major and lesser drainages of Ganga basin from Central Indian Himalayas (Gaula, Kosi, Alaknanda and Mandakini). River Gaula and Kosi are the main water sources of Kumaon region of Uttarakhand and River Alaknanda and Mandakini flow in Garhwal region of Uttarakhand, both regions are geographically isolated (Hora, 1954). The Gaula River is approximately 500km long. It originates in the Sattal lakes of Uttarakhand state and joins the Ramganga River (Major tributary of Ganga basin) about 15km. The river Kosi is one of the major tributaries of Ramganga River. It originates from village Budha Peenath of Kausani region of Uttarakhand (Kumar, & Bahadur, 2009). River Alaknanda is a major source of water in Garhwal region of Uttarakhand; it originates from Satopanth Glacier in Garhwal Himalaya (Menon, 1954). It is about 190km long. River Mandakini is a major tributary of river Alaknanda and originates from the Chorabari Glacier near Kedarnath in Garhwal region of Uttarakhand. It joins river Alaknanda at Rudraprayag (Mir et al., 2013). These are the four least explored rivers of Central Indian Himalaya and are an important fishery resource for the daily living of local people (Mir et al., 2013).
Sampling: A total of 402 specimens of Barilius bendelisis were collected from four different rivers across Central Indian Himalaya viz. Gaula, Kosi, Alaknanda and Mandakini, by using different fishing gear (cast nets and gill nets) between November 2012 and February 2013. The specimens of B. bendelisis were obtained before the breeding season and after the spawning period to avoid a bias toward size difference; every sampling station was visited six times during this period. The brooders and other small sized specimens obtained from the natural habitat during fishing were kept for breeding purpose and voucher specimens were kept in the museum at the Directorate of Coldwater Fisheries Research, India for possible future comparison studies. The mesh size of the fishing gear (cast nets: 9m length, 9m breadth and 1/2cm mesh size and drag nets: 100m length, 20m breadth, 1⁄2cm mesh size) was designed for the large sized specimens to avoid any fingerling and fry capture. The specimens were measured by using digital caliper and the total length range of the fish was 7.8-18.6cm. GPS coordinates, altitude (m.a.s.l), number of samples, min-max length and weight of B. bendelisis across the Central Indian Himalaya are included in Table 1.
Specimen digitization and truss distance measurements: The sampled specimens were first cleaned in running water, drained and placed on a flat platform with graph paper as a background to equilibrate scaling information to the digital images. The fins were erected to make the origin and insertion points visible. Each individual was labeled with a specific code to identify it. A Sony Cybershot DSC-W300 digital camera was used for capturing the digital images. To avoid errors in image capture, all photos were taken by a single person from same angle and height at every shot. After image capture, some specimens were dissected to identify the sex by macroscopic examination of the gonads. Gender was used as the class variable in ANOVA to test for significant differences in the morphometric characters, if any, between males and females of B. bendelisis.
The extraction of truss distances from the digital images of specimens was conducted using a linear combination of three software platforms, tpsUtil, tpsDig2 v2.1 (Rohlf, 2006) and Paleontological Statistics (Hammer, Harper, & Ryan, 2001). The truss protocol used for the B. bendelisis was based on 12 landmarks (Fig. 1A, Fig. 1B). A box truss of 30 lines connecting these landmarks was generated for each fish to represent the basic shape of the fish (Strauss & Bookstein, 1982; Mir et al., 2013). All the measurements were transferred to a spreadsheet file (Excel 2007), and the X-Y coordinate data transformed into linear distances by computer (using the Pythagorean Theorem) for subsequent analysis.
Multivariate data analysis: Size dependent variation was corrected by adopting an allometric method as suggested by Elliott, Haskard, & Koslow (1995):
Madj = M (Ls/L0)b
where M is the original measurement, Madj is the size adjusted measurement, L0 is the standard length of the fish, Ls the overall mean standard length, and b was estimated for each character from the observed data as the slope of the regression of log M on log L0 using all fish from every group. The results derived from the allometric method were confirmed by testing significance of the correlation between transformed variables and standard length (Turan, 1999). Univariate analysis of variance (ANOVA) was performed for the 30 morphometric characters to evaluate the significant difference between the four locations. These 30 transformed truss measurements were subjected to FACTOR analysis, to explain these variables in terms of their common underlying dimensions. A maximum likelihood method was used to extract the factors. With the assistance of Scree plot, the cumulative variance explained by the factors and the meaningful biological groupings of the traits loading on each factor were taken into consideration to retain the number of factors for a rotation procedure. The retained factors were subjected to a Varimax rotation procedure and to identify the variables demonstrating high loadings for a given factor, the rotated factors were subjected to a scratching procedure, as described by Hatcher (2003). The Wilks’ λ was used to compare the difference between all groups. The discriminant function analysis (DFA) was used to calculate the percentage of correctly classified (PCC) fish and a cross-validation using PCC was done to estimate the expected actual error rates of the classification functions. Statistical analyses for morphometric data were performed using the SPSS vers. 16.1.0 and Microsoft Excel 2007.
Results
None of the size-adjusted truss measurements showed a significant correlation with the standard length of the fish, indicating that the variation in body length had been successfully removed by the allometric transformation. Among four selected rivers, means of all the truss measurements of B. bendelisis were found to be significantly (p<0.001) different in univariate analysis of variance. The morphometric characters between two sexes did not differ significantly (p>0.05) hence; the data for both sexes were pooled for all subsequent analysis.
Wilks’ lambda tests of DFA indicated significant differences in morphometric characters of all populations (p<0.001). In discriminant function analyses, three discriminant functions (DFs) were formed. The first discriminant function (DF I) accounted for 98.9% of the total variation. The second discriminant function (DF II) and third discriminant function (DF III) accounted for 0.6 and 0.5%, respectively, between group variability among the populations. The morphometric measurements, 1-2, 1-12, 1-10, 4-11, 3-11, 6-8, 4-6, 4-7, 3-8, 5-8, 4-5 and 4-8 contributed to DF I, 5-6, 9-11, 2-3, 8-11 and 2-12 contributed to DF II and 6-7, 3-4, 10-12, 11-12, 5-7, 4-9, 2-11, 8-9, 10-11, 2-10, 7-8 and 3-9 contributed to DF III (Table 2), showing that these characters were the most important in distinguishing the population.
The DF I vs. DF II plot explained 99.5% of total variance among the samples and formed two separate ]]>
Fig. 2). All the samples from the rivers of Central Indian Himalayan regions were clearly separated from each other in discriminant space. However, Alaknanda River showed overlapping with Mandakini River and river Gaula showed intermingling with Kosi River.
]]>
The classification of individuals into their original population varied between 83.8% and 72.8% by discriminant analysis and 88.0% of individuals could be classified into their original a priori grouping (Table 3). The proportion of Gaula River samples into their original population was the highest (83.8%). A cross-validation test using leave one out procedure was also performed by which 81.98% of the samples were classified into their original ]]>
Table 3).
Discussion
]]>
The results of present study indicated that the B. bendelisis showed significant phenotypic heterogeneity between the two geographically isolated Central Himalayan regions. Discriminant function analysis (DFA) could be a useful method to distinguish different stocks of the same species (Karakousis, Triantaphyllidis, & Economidis, 1991). In the present paper, 78% of ]]>
Schizothorax richardsonii from different regions of Indian Himalaya, finding significant morphometric heterogeneity among different populations by applying DFA and it to migration of fish and environmental factors. The DFA confirmed that the variation in morphological measurements was evident in the whole body, contributed mainly by middle portion, head region, body depth and caudal peduncle, between these morphologically distinct populations of
B. bendelisis. Hossain, Nahiduzzaman, Habiba, Mst & Alam (2010) applied DFA and PCA on three populations of Labeo calbasu from river Jamuna, Halda and hatchery, and reported morphological discrimination among them due to the environmental factors and local migration of the fish. Similar observations were noticed by Khan, Miyan & Khan (2012) in case of Channa punctatus from three ]]>
Schizothorax curvifrons from Indus river basin, and attributed to changing physical and ecological conditions of water bodies.
These morphological ]]>
]]>
The causes of morphological differences among different populations are often quite difficult to explain (Poulet, Berrebi, Crivelli, Lek, & Argillier, 2004). It has been suggested that the morphometric differentiation of a fish can be influenced by genetic, environment and the interaction between them (Mir et al., 2013). Apparently, the fragmentation of river impoundments can lead to an enhancement of pre-existing genetic differences, providing a high interpopulational structuring (Esguicero, & ]]>
The water fragmentation obstructs the migration of fish populations resulting in an ecological trap (Jager, Chandler, Lepla, & Winkle, 2001). A series of barrages have been commissioned in the upper segment of river Ganges from Rishikesh ]]>
The inferences of this study leads to the conclusion that B. bendelisis has distinct stocks in the two geographically isolated regions of Central Indian Himalaya, while the ]]>
]]>
Acknowledgments
The Department of Biotechnology, New Delhi is acknowledged for funding this study. Thanks are due to Director, Directorate of Coldwater Fisheries Research, Bhimtal, Uttarakhand for providing necessary facilities for the study.
References
Akbarzadeh, A., Farahmand, H., Shabani, A. A., Karami, M., Kaboli, M., Abbasi, K., & Rafiee, G. R. (2009). Morphological variation of the pikeperch Sander lucioperca (L.) in the southern Caspian Sea, using a truss system. Journal of Applied Ichthyology, 25, 576-82. [ Links ]
AnvariFar, H., Khyabani, A., Farahmand, H., Vatandoust, S., AnvariFar, H., & Jahageerdar S. (2011). Detection of morphometric differentiation between isolated up- and downstream populations of Siah Mahi (Capoeta capoeta gracilis) (Pisces: Cyprinidae) in the Tajan River (Iran). Hydrobiologia, 673, 41-52. [ Links ]
Costa, J. L., Almeidam P. R., & Costa, M. J. (2003). A morphometric and meristic investigation of Lusitanian toadfish Halobatrachus didactylus (Bloch and Schneider, 1081): evidence of population fragmentation on Portuguese coast. Scientia Marina, 67, 219-231. [ Links ]
Daud, S. K., Mohammad, M., Siraj, S. S., & Zakaria, M. P. (2005). Morphometric analysis of Malaysian Oxudercine Goby, Boleophthalmus boddarti (Pallas, 1770). Pertanika Journal of Tropical Agricultural Sciences, 28(2), 121-134. [ Links ]
Eschmeyer, W. N., & Fricke, R. (Eds.). (2011). Catalog of Fishes electronic version (29 March 2011). Retrieved from http://research.calacademy.org/ichthyology/catalog/fishcatmain.asp.Accessed on 30th March 2011. [ Links ]
Esguicero, A. L. H., & Arcifa, S. A. (2010). Fragmentation of a Neotropical migratory fish population by a century old dam. Hydrobiologia, 638, 41-53. [ Links ]
Hammer, O., Harper, D. A. T., & Ryan, P. D. (2001). PAST: paleontological statistics software package for education and data analysis. Palaeontologia Electronica, 4, 9. [ Links ]
Hatcher, L. (2003). A step by step approach to using SAS for factor analysis and structural equational modeling. SAS Institute Inc., 57, 125. [ Links ]
Hora, S. L. (1954). The evolution of the Indian torrential environment and its fishes. Bulletin of National Institute of Sciences of India, 17, 437-444. [ Links ]
Hossain, M. A. R., Nahiduzzaman, M., Habiba, K. D. S., Mst, U., & Alam, M. S. (2010). Landmark-based morphometric and meristic variations of the endangered Carp, Kalibaus Labeo calbasu, from stocks of two isolated rivers, the Jamuna and Halda, and a hatchery. Zoological Studies, 49, 556-563. [ Links ]
IUCN (2012). Red List of Threatened Species, vers. 2012.1. [ Links ]
Jager, H. I., Chandler, J. A., Lepla, K. B., & Winkle, W. V. (2001). A theoretical study of river fragmentation by dams and its effect on white sturgeon populations. Environmental Biology of Fishes, 60, 347-361. [ Links ]
Karakousis, Y., Triantaphyllidis, C., & Economidis, P. S. (1991). Morphological variability among seven populations of brown trout, Salmon trutta L., in Greece. Journal of Fish Biology, 38, 807-817. [ Links ]
Khan, M. A., Miyan, K., & Khan, S. (2012). Morphometric variation of snakehead fish, Channa punctatus, populations from three rivers. Journal of Applied Ichthyology, 28, 154-155. [ Links ]
Kumar, A., & Bahadur, Y. (2009). Physico-chemical studies on the pollution potential of river Kosi at Rampur (India). World Journal of Agricultural Sciences, 5, 01-04. [ Links ]
McAllister, D. E., Craig, J. F., Davidson, N., Delany, S., & Seddon, M. (2001). Biodiversity Impacts of Large Dams: Dams, ecosystem functions and environmental restoration. On behalf of IUCN-The World Conservation Union. [ Links ]
Menon, A. G. K. (1954). Fish geography of the Himalayas. Zoological Survey of India, 15, 467-493. [ Links ]
Menon, A. G. K. (1999). Check list - fresh water fishes of India. Zoology Survey of India, 175, 366. [ Links ]
Mir, F. A., Mir, J. I., & Chandra, S. (2013). Phenotypic variation in the Snowtrout Schizothorax richardsonii (Gray, 1832) (Actinopterygii: Cypriniformes: Cyprinidae) from the Indian Himalayas. Contributions to Zoology, 82(3), 115-122. [ Links ]
Mir, F. A., Mir, J. I., & Chandra, S. (2014). Detection of morphometric differentiation in Sattar snowtrout, Schizothorax curvifrons (Cypriniformes: Cyprinidae) from Kashmir Himalaya using a truss network system. Revista de Biologia Tropical, 62(1), 119-127. [ Links ]
Murta, A. G. (2000). Morphological variation of horse mackerel (Trachurus trachurus) in the Iberian and North African Atlantic: Implications for stock identification. ICES Journal of Marine Science, 57, 1240-1248. [ Links ]
Poulet, N., Berrebi, P., Crivelli, A. J., Lek, S., & Argillier, C. (2004). Genetic and morphometric variations in the pikeperch (Sander lucioperca L.) of a fragmented delta. Archiv für Hydrobiologie, 159, 531-554. [ Links ]
Rahman, A. K. A. (1989). Freshwater fishes of Bangladesh. Zoological Society of Bangladesh, 364. [ Links ]
Rao, R. J. (2001). Biological resources of the Ganga River. Hydrobiologia, 458, 159-168. [ Links ]
Rohfritsch, A., & Borsa, P. (2005). Genetic structure of Indian scad mackerel Decapterus russelli: Pleistocene variance and secondary contact in the Central Indo-West Pacific seas. Heredity, 95, 315-322. [ Links ]
Rohlf, F. J. (2006). tpsDig2, vers. 2.1. State University of New York, Stony Brook. Retrieved from: http://life.bio.sunysb.edu/morph. [ Links ]
Sah, S., Barat, A., Pande, V., Sati, J., & Goel, C. (2011). Population Structure of Indian Hill Trout (Barilius bendelisis) inferred from variation in Mitochondrial DNA sequences. Advances in Biological Research, 5, 93-98. [ Links ]
Sahoo, P. K., Saikia, S. K., & Das, D. N. (2009). Natural food resources and niche breadth of Barilius bendelisis (Hamilton) (Pisces, Cyprinidae) in river Dikrong, an upland riverine ecosystem in India. Pan-American Journal of Aquatic Sciences, 4(1), 12-16. [ Links ]
Samaee, S. M., Mojazi-Amiri, B., & Hosseini-Mazinani, S. M. (2006). Comparison of Capoeta capoeta gracilis (Cyprinidae, Teleostei) populations in the south Caspian Sea River basin, using morphometric ratios and genetic markers. Folia Zoologica, 55, 323-335. [ Links ]
Strauss, R. E., & Bookstein, F. L. (1982). The truss: body form reconstruction in morphometrics. Systematic Zoology, 31, 113-135. [ Links ]
Swain, D. P., & Foote, C. J. (1999). Stocks and chameleons: the use of phenotypic variation in stock identification. Fisheries Research, 43, 113-128. [ Links ]
Talwar, P. K., & Jhingran, A. G. (1991). Inland fishes of India and adjacent countries. Vol. 1. Balkema AA, Rotterdam. [ Links ]
Turan, C. (1999). A note on the examination of morphometric differentiation among fish populations: the truss system. Turkish Journal of Zoology, 23, 259-263. [ Links ]
Tzeng, T. D. (2004). Morphological variation between populations of spotted mackerel (Scomber australasicus) off Taiwan. Fisheries Research, 68, 45-55. [ Links ]
Valdiya, K. S. (1991). “Hydrogeological Studies of springs in the catchment of the Gaula River”. Mountain Research and Development, 11, 239-258. [ Links ]
Wimberger, P. H. (1992). Plasticity of fish body shape. The effects of diet, development, family and age in two species of Geophagus (Pisces, Cichlidae). Biological Journal of Linnaean Society, 45, 197-218. [ Links ]