1 Introduction
The Brazilian rivers have been sorely segmented by dams. Fish stocking is one of the management actions used to minimize the impacts of dams on native fish fauna (Agostinho et al., 2008). It consists of the breeding and cultivation of fish in captivity for later release into natural water bodies. It is normally applied when the management goal is to supplement the native populations (Agostinho et al., 2007). However, despite being historically applied, the efficiency of stocking, as a management measure, has been not satisfactory in Brazil, since the target stocks have not been recovered (Agostinho et al., 2007). This can be explained by poor habitat quality of the stocked rivers (Pompeu, et al., 2012) and by differences between hatchery and wild environment.
Hatchery tanks are characterized by a high degree of environmental homogeneity (Pakkasmaa et al., 1998). They are normally made of concrete, with no substrate, obstacles or shelter and have constant food availability, high population density and species monoculture. As a result, farmed fish differ from wild fish in several aspects, especially in the morphological ones, exhibiting distinct colorations (Maynard et al., 1995); different body, head and fin proportions (Swain et al., 1991; Fleming et al., 1994; Hard et al., 2000; Solem et al., 2006; Belk et al., 2008) and less morphological variability (Taylor, 1986). Moreover, density-dependent processes probably influence early growth and survival (Winemiller and Dailey, 2002). On the other hand, natural environments are marked by many environmental variations, whether physical, chemical or biological ones. Marginal lagoons inhabited by many neotropical migratory fish species in the early stages of development, are rich in food and shelter provided by aquatic plants, wood debris, and other structures (Agostinho and Júlio-Junior, 1999), and have high species richness with diverse social interactions (Pompeu and Godinho, 2003). As a result of these environmental differences, the survival rates of released fish are low, because they are adapted to the rearing environment and not to the wild (Maynard et al., 1995). The selection of phenotypes adapted to the captivity is called domestication, which can affect the physiology, behavior and morphology of fish (Araki et al., 2008).
Specifically in relation to morphology, polymorphisms between alternative environments may occur in a same species, in response to different selective pressures (Robinson and Wilson, 1996; Langerhans et al., 2003). In this respect, a range of studies conducted, mainly in temperate countries, has demonstrated the existence of wide morphological differences between wild and farmed fish (Taylor, 1986 ; Swain et al., 1991; Pakkasmaa et al., 1998; Hard et al., 2000; Ojanguren and Brana, 2003; Cramon-Taubadel et al., 2005; Berejikian and Tezak, 2005; Solem et al., 2006; Belk et al., 2008).
Fisheries management in South America has been based on stocking, fisheries regulation and construction of fish passes (Pompeu et al., 2012). However, although fish stocking have been performing for years, we lack basic studies evaluating survival rates of released fish, and the efficiency of this strategy in recovering populations (Vieira and Pompeu, 2001). Until now, there are no studies comparing wild and hatchery fish in respect to morphological aspects. Based on this, the primary goal of this study was to verify if hatchery fish are morphologically different from wild individuals of the same species and size class.
2 Material and Methods
2.1 Species studied
Prochilodus lineatus (Achille Valenciennes, 1837) (Actinopterygii: Prochilodontidae), popularly known as curimba, curimbatá, curimatã or papa-terra, is a native fish species of the Paraná River Basin. Its feeding habits are detritivorous, eating mostly decaying organic matter (Hahn et al., 2004). Its reproductive habits are migratory, swimming upriver to spawn in the breeding season. The adults inhabit the main channel of the river, showing preference for river stretches surrounded by wetlands (Pesoa and Schulz, 2010). The young inhabit the marginal lagoons of the floodplain, where they remain during the first years of life (Agostinho and Júlio-Júnior, 1999; Graça and Pavanelli, 2007; Pesoa and Schulz, 2010). Prochilodus lineatus was chosen due to its ecological relevance, importance for fishing, and because is one of the most used native fish in stocking and repopulation. It is also used in Brazilian hatchering programs as live food for other native species, mainly piscivores, because of its easy management in pisciculture (Viveiros et al., 2009).
2.2 Collections
Hatchery individuals of Prochilodus lineatus were obtained at the hatchery station of Furnas Hidroeletric Plant, located on the Grande river, in Passos county, Minas Gerais state. A total of 62 individuals with standard length ranging from 10.80 to 16.70 cm were obtained. In this station, fish are produced from wild breeders collected around Furnas reservoir, since its purpose is to perform the repopulation of the areas affected by the dam, including the main tributaries of the reservoir. Fish tanks at the hatchery station have homogeneous conditions, because are made up of concrete walls, characterized by the absence of substrates, barriers or shelters and by high population density, food availability and species monoculture.
Since the aim of this study is compare wild populations with juveniles from hatcheries, wild specimens were captured using gillnets in four marginal lagoons of the Mortes river, a tributary of the Grande river (see Table 1). Floodplain lagoons are considered the main nurseries areas for this species, and the Mortes river is one of the main lotic remnant of the region (Suzuki et al., 2011). The marginal lagoons of this region are rich in food and shelter provided by aquatic plants, twigs and other structures. Exhibits considerable seasonality, provided by variations in hydrological conditions and physicochemical parameters. The lagoons still have high species richness with diverse social interactions (pers. obs.). Collections were conducted during the rainy season of 2011, and a total of 70 individuals were captured, with standard length ranging from 10.99 to 16.57 cm. All fish were sacrificed by overdose of eugenol, fixed in 10% formalin and preserved in 70% alcohol. The voucher specimens were deposited in the Ichthyological Collection of The Federal University of Lavras (CI-UFLA 0450 and 0859).
2.3 Morphometry
Seventeen morphometric measurements (see Table 2) were taken from each specimen. The linear measurements were obtained directly from the fish body, on the left side, with a digital caliper (0.01 mm precision), and drawings of the outline of the body and fins of each fish were made on paper to obtain the body and fin areas. The drawings were then scanned and the areas were calculated using the Image J software (Abramoff et al., 2004). To minimize differences in the measurements due to the specimen size, 12 ecomorphological attributes were calculated (see Table 3), which were chosen according to their ecological significance defined in the literature (Hora, 1930; Gatz Junior, 1979; Watson and Balon, 1984; Pouilly et al., 2003).
Table 2 Morphometric measurements and their descriptions.
Morphometric measurements | |||
---|---|---|---|
Measurement | Acronym | Description | |
1. | Standard Length | SL | Distance from tip of the snout to the posterior margin of caudal peduncle |
2. | Body Height | BH | Greatest dorsal-ventral distance perpendicular to the longest body axis |
3. | Body Width | BW | Greatest side to side body width |
4. | Average Body Height | ABH | Distance from the abdomen to the line that cuts the body between the mouth and the tail |
5. | Body Area | BA | Body area plus the area of the caudal fin |
6. | Head Length | HL | Distance between the tip of the snout to the end of the operculum |
7. | Head Height | HH | Distance between the dorsal and ventral part of the head in the region of the eyes |
8. | Eye Height | EH | Distance from the center of the eye to the lower jaw |
9. | Caudal Peduncle Length | CPL | Distance between the end of the anal fin to the top of the caudal |
10. | Caudal Peduncle Height | CPH | Height of peduncle measured at the same point as the width |
11. | Caudal Peduncle Width | CPW | Peduncle width measured at its midpoint |
12. | Caudal Fin Width | CFW | Maximum distance between the two ends of the fin fully extended |
13. | Caudal Fin Area | CFA | Fully extended Caudal fin area |
14. | Pelvic Fin Area | PlFA | Fully extended pelvic fin area |
15. | Pectoral Fin Length | PtFL | Distance between the base of the fin and its end |
16. | Pectoral Fin Width | PtFW | Largest width of the fin in a plane perpendicular to the length axis of the fin fully open |
17. | Pectoral Fin Area | PtFA | Fully extended Pectoral fin area |
Table 3 Descriptions of ecomorphological attributes and their ecological significance.
Ecomorphological attributes | |||
---|---|---|---|
Attribute | Formula | Significance | |
1. | Compression Index | CI = BH/BW | Higher values indicate laterally compressed fish that inhabit areas of low water velocity (Watson and Balon, 1984). |
2. | Ventral flattening Index | VFI = ABH/BH | Low values are typical of fish living in waters with high speed, enabling it to maintain its position without swimming (Hora, 1930). |
3. | Relative height of the body | RH = BH/SL | Inversely related to the speed of the water and directly related to the ability to develop vertical displacements (Gatz Junior, 1979). |
4. | Relative position of the eye | RPE = EH/HH | This index is related to food detection and provides information about the use of vision in predation activities (Pouilly et al., 2003). It is indicative of the preferred position in the water column (Gatz Junior, 1979). |
5. | Relative length of the head | RLH = HL/SL | Directly related to the size of the prey, high values suggest species which are predators of relatively large prey (Gatz Junior, 1979). |
6. | Relative length of the caudal peduncle | RLCP = CPL/SL | Relatively long peduncles indicate fish that inhabit turbulent waters and have good swimming skills (Gatz Junior, 1979; Watson and Balon, 1984). |
7. | Compression Index of the caudal peduncle | CICP= CPH/CPW | It is inversely related to the amplitude of the swimming movements, higher values indicate compressed peduncles, typical of less active swimmers (Gatz Junior, 1979). |
8. | Relative area of the caudal fin | RACdF = CdFA/BA | High values indicate fins capable of producing rapid pulses, typical for swimming of many benthic fish (Watson and Balon, 1984); |
9. | Aspect ratio of the caudal fin | ARCdF = CdFW/CdFA | Directly proportional to the amount of swimming that the fish conducts (Gatz Junior, 1979). |
10. | Relative area of the pelvic fin | RAPlF = PlFA/BA | Benthic fish have relatively large areas (Gatz Junior, 1979). |
11. | Relative area of the pectoral fin | RAPtF = PtFA/BA | Higher values indicate slower swimmers who use these fins to perform maneuvers. However, fish inhabitants of fast water that use them as water deflection surfaces, keeping their body close to the substrate, also have high values (Watson and Balon, 1984). |
12. | Aspect ratio of the pectoral fin | ARPtF = PtFL/PtFW | High values indicate long fins and suggest fish with great swimming ability (Watson and Balon, 1984). |
Ecomorphological attributes are ratios resulting from the division of a morphometric measurement by another, being used from the assumption that this measurements are proportional to each other (Watson and Balon, 1984). To verify the proportionality of the morphometric measurements used, we performed linear regressions between them. All the regressions were positive and significant (p = 0.0001), confirming the proportionality between the measurements and allowing to use the ecomorphological attributes resulting from its division.
2.4 Statistical analysis
Initially, through the Statistica software, the similarity between the mean standard length of both population (wild and hatchery fish) was confirmed using a T test. Similar lengths are necessary to enable comparison between the two populations, since there is a dependency between shape and body size (Peres-Neto, 1999). Then, Principal Component Analysis (PCA) was used to evaluate the distribution of individuals in morphological space. To investigate the differences between the populations in the first two PCA axes a T test was conducted and also a Discriminant Analysis (DA), as an alternative test.
Since the tested hypotheses tested are strongly associated with the response of the individual, resulting in a population response, the individuals accounted for the replicates. For each population a normalized Euclidean distance matrix was constructed, in the Past software (Hammer et al., 2001), between the pairs of individuals from the morphological attributes, according to the Equation 1.
where Djk is the Euclidean distance between individuals j and k, n is the number of attributes, xij and xik are the values of the attribute i for the pair of individuals jk.
In the case of normalized distance, the data were adjusted to the Gaussian distribution with average zero and standard deviation one, according to Equation 2.
where xik is the average of attribute k of population i, Xk the average for all xik and DPK, the standard deviation for the attribute k.
From each distance matrix were calculated: the distance of the closest neighbor (DCN), as a dispersion index of individuals in ecomorphological space (Ricklefs et al., 1981); and the distance to the centroid of each individual of the population (DC), as an estimate of the relative size of the morphological hypervolume occupied by the population (Winemiller, 1991).
To evaluate possible differences between the distance from the centroid and the closest neighbor we applied an averages test, through the Statistica software. The data of the distance from the centroid were normalized by logarithmization in base ten and tested by the T test. The data of the distance from the closest neighbor could not be normalized and thus the nonparametric Kruskal-Wallis test was used.
3 Results
The analyzed fish of wild and hatchery populations showed a similar standard length (t(1;131) = 0.1441; p = 0.7048), allowing comparison between the two groups.
The Principal Component Analysis (PCA) indicated that each population occupies completely distinct morphological spaces (see Figure 1). The first two PCA axes were the most important in explaining the total variation of the sample (PCA1 = 34.73% and PCA 2 = 12.08%), jointly explaining 46.81% of variance existing between the two populations and due to the ecomorphological attributes analyzed.

Figure 1 Projection of two populations of Prochilodus lineatus [H = hatchery; W = wild] in the first two axes of the Principal Component Analysis.
The T test conducted between the scores of the first two PCA axes showed significant differences only for the first axis (PCA 1: p = 0.0001; MS = 453.7531; t = 640.2777; PCA 2: p = 0.9113; MS = 0.0181; t = 0.0124). Thus, the ecomorphological attributes that more contributed to the differentiation between populations in the first axis of the PCA, make up the main source of morphological differentiation among populations. The most important attributes in PCA 1 (considering contributions higher than 70%) were compression index (CI), relative body height (RH), relative area of the pelvic fin (RAPlF) and aspect ratio of the caudal fin (ARCdF).
The discriminant analysis was also significant between the two analyzed populations (Wilks’ Lambda: 0.10259; F(12,119) = 86.748 p = 0.0001). The ecomorphological attributes that contributed significantly to the differentiation between them were relative height (RH) (F (3,176) = 72.03539, p = 0.00000), compression index of the caudal peduncle (CICP) (F (3,176) = 17.33353, p = 0.00006), relative area of the pectoral fin (RAPtF) (F (3,176) = 17.22258, p = 0.00006), aspect ratio of the caudal fin (ARCdF) (F (3,176) = 6.37739, p = 0.01287) and relative length of the head (RLH) (F (3,176) = 4.76763, p = 0.03096) (see Table 4).
Table 4 Discriminant analysis between wild and hatchery populations.
Prochilodus lineatus | ||
---|---|---|
Attributes | F (3,176) | p Value |
CI | 0.64579 | 0.423226 |
RH | 72.03539 | 0.000000 |
RLCP | 0.00026 | 0.987277 |
CICP | 17.33353 | 0.000060 |
VFI | 0.54897 | 0.460200 |
RAPtF | 17.22258 | 0.000063 |
ARPtF | 3.36385 | 0.069141 |
RAPlF | 0.55145 | 0.459190 |
RACdF | 1.96947 | 0.163108 |
ARCdF | 6.37739 | 0.012873 |
RPE | 0.64820 | 0.422364 |
RLH | 4.76763 | 0.030962 |
(Wilks’ Lambda: 0.10259, F (12,119) = 86.748 p = 0.0001). The attributes in bold were those that best contributed to discrimination.
Considering the attributes individually, all the important ones to the differentiation in PCA, as well as in the DA analysis performed, were significantly different between wild and hatchery fish (Wilks’ Lambda: 0.10259; F (12,119) = 86.748; p = 0.0001). Two attributes, RH and ARCdF, were pointed out by the two analyses as the most important. These two attributes showed higher values for the wild population, indicating that the wild fish had a higher body (RH: KW(1; 132) = 97.71758, p = 0.0001) and larger caudal fin (ARCdF: t(1; 130) = 90.4828, p = 0.0001) (see Figure 2).

Figure 2 Most important ecomorphological attributes in morphological differentiation between populations of Prochilodus lineatus in PCA1 and DA concomitantly.
The morphological variability of the fish was analyzed by comparing the distance from the centroid and closest neighbor (see Figure 3). The morphological variability of the wild population was higher than that of the hatchery population, because the wild fish presented higher dispersion in morphological space (ANOVA log(DC+1): F(1; 131) = 11.9207, p = 0.0007) and thus occupied a larger morphological hypervolume (Kruskal-wallis’ test CDN: KW(1; 133) = 18.7552; p = 0.0001).

Figure 3 Log of the distance from the centroid plus one [log (DC +1)] and the distance from the closest neighbor. The ANOVA performed for the distance from the centroid was significantly higher for the population of wild fish, as was the Kruskal-Wallis test conducted for the distance from the nearest neighbor.
4 Discussion
The results show that the fish hatcheries creates a less diverse population and individuals morphologically different from its wild counterparties of the same size class. Selective pressures present in the fish farming tank are quite different from those existing in the natural environment, which may result in very different morphologies generated by phenotypic plasticity (Taylor, 1986; Swain et al., 1991; Fleming et al., 1994; Hard et al., 2000; Belk et al., 2008). Therefore, in the present study, the morphological differences between wild and hatchery fish could have been generated by phenotypic plasticity caused by adaptation to the growth environment, or may have been caused by genetic differences between populations.
The absence of genetic data does not allow separating the effects of these two sources of variation. But as the hatchery population also has a wild origin, it is concluded that the morphological differences detected were mostly due to the environmental component. One evidence of this, it is that the differentiation between the two studied populations can be easily understood considering the environmental differences between the lagoons and hatchery tanks. For example, the relative body height (RH) was higher for wild fish. This attribute is directly related to the ability of the fish to develop vertical displacements (Gatz Junior, 1979). Therefore, the higher RH of wild fish is consistent with the need to perform maneuvers to avoid obstacles, considering the structural richness of the lagoons. Wild fish also had higher aspect ratio of the caudal fin (ARCdF), attribute that is directly proportional to the amount of swimming that the fish does (Gatz Junior, 1979). This can be an indication that wild fish swim more than the hatchery, probably because they have a larger and more complex environment to explore. Morphological differences between populations of wild and hatchery fish have been previously highlighted, and environmental factors have been suggested as the principal source of variation (Taylor, 1986; Fleming et al., 1994; Ellis et al., 1997; Cramon-Taubadel et al., 2005; Belk et al., 2008). However, the role of intraspecific interactions in wild environments as a morphological differentiation mechanism cannot be excluded. During interspecific competition, for instance, individuals from one species suffer a reduction in fertility, growth or survival as a result of resource exploitation or interference by another species (Begon et al. 2007).
In relation to the morphological variability, the explanation for the lower value found for the hatchery population may lie in the lower genetic variability of this population (Vuorinen, 1982; Agostinho et al., 2007), which may consequently result in minor morphological variability. However, recent studies (Gomes et al., 2008; Chittenden et al., 2010) have shown the presence of minor genetic differences when comparing wild specimens with hatchery, since currently genetic diversity of stocked fish has been expanded from the use of a high number of breeder individuals of wild origin. Thus, the low morphological variability found for the hatchery population likely is due to the high homogeneity of the tanks (Hynes et al., 1981; Taylor, 1986), where environmental pressures are constant, due to low disturbance frequency. The intraspecific polymorphism is typically caused by different selective pressures (Langerhans et al., 2003). Thus, in heterogeneous environments such as lagoons, subject to numerous pressures and disturbances, various morphological types can coexist, as each may be favored at different times, depending on the disorder or pressure acting on the population. On the other hand, in homogeneous environments, such as farming tanks, the constant selective pressure favors a basic type of morphology, so that individuals tend to be morphologically closer.
Although stocking is intended to replenish stocks that have experienced genetic and/or demographic losses, the can carry environmental risks, even when performed with native species. These potential impacts include dissemination of pathogens and parasites and deleterious effects related to the genetic quality of matrices and fingerlings (bottleneck effects, loss of genetic variability and fitness, domestication) (Agostinho et al., 2010). The release of a less variable and morphologically distinct population does not represent itself an impact. However, these aspects must represent lower probability to survive, hindering the use of this strategy as a conservation tool.
In summary, our results demonstrate that hatchery populations of Prochilodus lineatus are morphologically different and less diverse when compared with wild ones. However, for the actual consequences for hatchery management it is necessary to carry out specific experiments and studies that aim to elucidate how much and for how long these morphological differences could affect performance and survival of fish released in the natural environment. Additionally, considering the lower morphological variability of hatchery populations and the concept of environmental heterogeneity, the enrichment of hatchery tanks can be proposed as a way of increasing the morphological variability of the produced fingerlings.