versión On-line ISSN 0718-560X
Lat. Am. J. Aquat. Res. vol.40 no.2 Valparaíso jun. 2012
Lat. Am. J. Aquat. Res., 40(2): 408-417, 2012
Population genetic structure revealed by a school of the freshwater migratory fish, Brycon hilarii
Estructura genética poblacional revelada por un cardumen del pez migratorio de agua dulce, Brycon hilarii
Alexandra Sanches1 & Pedro Manoel Galetti Jr.2
1Departamento de Ecologia, Universidade Estadual Júlio de Mesquita Filho P.O. Box 13506-900, Rio Claro, SP, Brazil
2Departamento de Genética e Evolução, Universidade Federal de São Carlos Via Washington Luis, P.O. Box 676, 13565-905, São Carlos, SP, Brazil
ABSTRACT. It is believed that different genetic populations of migratory fishes can co-exist in a single hydrographic system. Although different populations may occupy and explore the river together, they segregate during the spawning season and consequently the population genetic structuring is maintained. Genetic variation of a Brycon hilarii spawning school and samples from different collection sites in the Miranda River basin were surveyed using seven microsatellites. Population structuring was revealed by a significant differentiation of the spawning school related to the supposed local populations. The genetic differentiation detected may be supported by behavior during the reproductive season that promotes the maintenance of the genetic integrity of different populations. These data may contribute toward the understanding of the behavior and biology of these fish as well as fishery management and species conservation programs.
Keywords: Characidae, Bryconinae, neotropical fish, population structuring, microsatellite, Brasil.
RESUMEN. Se ha documentado que diferentes poblaciones genéticas de peces migratorios pueden coexistir en un único sistema hidrográfico. Diferentes poblaciones pueden ocupar y explorar el río juntas, pero se segregan durante la temporada de desove y consecuentemente la estructuración genética poblacional se mantiene. La diversidad genética de un cardumen reproductivo de Brycon hilarii y muestras de diferentes sitios en la cuenca del Río Miranda fueron analizadas mediante siete microsatélites. La estructura poblacional fue revelada por una diferenciación genética significativa del cardumen reproductivo con las muestras de las poblaciones locales. La diferenciación genética detectada puede ser resultado de un probable comportamiento durante la temporada reproductiva, que promueve el mantenimiento de la integridad genética de las diferentes poblaciones. Estos datos pueden contribuir a la comprensión del comportamiento y biología de estos peces, así como amparar programas de gestión de la pesca y conservación de las especies.
Palabras chave: Characidae, Bryconinae, peces neotropicales, estructuración poblacional, microsatélites, Brasil.
During flood seasons, Neotropical freshwater migratory fish form large schools that swim upstream at reproduction time (Lowe-McConnell, 1987; Carolsfeld et al., 2003). Studies on radio-telemetry have shown the ability of these fishes to migrate several hundred kilometers for spawning or feeding (Lowe-McConnell, 1987; Godinho & Kynard, 2006; Godinho et al., 2007). Therefore, a single large panmitic population within a continuous hydrographic system would be an expected situation for these migratory fishes.
However, a number of studies have indicated that even vagile fish species may contain populations genetically different within a single hydrographic basin. It is believed that two or more genetic populations of migratory fish can co-exist within a continuous hydrographic system (Hatanaka et al., 2006; Sanches & Galetti, 2007; Piorski et al., 2008).
Molecular analysis of Brycon orthotaenia (= Brycon lundii) enable to detect a significant population structuring in the São Francisco River basin, Brazil (Wasko & Galetti, 2002). These authors hypothesized that this migratory fish constitute different coexisting and comigrating populations along the main river channel. Similar results were reported for another migratory fish that also inhabits the São Francisco River - Prochilodus argenteus (Hatanaka & Galetti, 2003; Hatanaka et al., 2006).
Although different populations may occupy and explore the river together, they segregate during the spawning season and consequently the population genetic structuring is maintained. Homing behavior (Gerking, 1959; Hatanaka et al., 2006) and spawning waves (Jorgensen et al., 2005) have been hypothesized as mechanisms responsible for maintaining the genetic cohesiveness of a fish population. Homing is a behavior mainly observed in salmonids and carps (Stähl, 1987), in which mature individuals return to the birthplace for reproduction rather than going to other equally viable places (Gerking, 1959). Mark-recapture and radio tracking studies have evidenced that migratory fish species, such as Prochilodus lineatus (Godoy, 1975), P. argenteus (Godinho & Kynard, 2006) and Pseudoplatystoma corruscans (Godinho et al., 2007) show homing behaviour. Spawning waves constitute different populations that spawn at a same location but at different times, even if with some overlap, as observed in some herring species (McPherson et al., 2003; Jorgensen et al., 2005).
Fish have diverse reproductive behaviors with a rich natural-history literature on mating systems ranging from pelagic group spawning, cooperative breeding to social monogamy (for a review, see Avise et al., 2002). Several reports have documented the inbreeding avoidance through the kin recognition among fish species (Pusey & Wolf, 1996; Gerlach & Lysiak, 2006; Ala-Honkola et al., 2010). On the other hand, recent studies have demonstrated that a number of fish species are able to recognize relatives or non-relatives individuals, and that they can show preference for schooling or mating with the relative ones (Odling-Smee & Braithwaite, 2003; Ward & Hart, 2003; Fraser et al., 2005). These mechanisms could provide a limited exchange of reproductive individuals between populations, resulting in genetic differences. However, the life history movements and spawning of freshwater Neotropical fish are poorly known (Godinho et al., 2007).
Significant genetic differentiation was recently revealed between populations of Brycon hilarii inhabiting small streams of the Miranda River basin in the Bonito region of western Brazil (Sanches & Galetti, 2007). These authors assessed the genetic variation of this species during different times of the year and detected a significant population differentiation. According to the authors, since significant differences were detected between population subunits, spawning schools organize themselves in such way as to maintain the integrity of each subunit residing in the system.
B. hilarii is a medium-sized migratory species, widely distributed throughout the Paraguay River basin (Lima, 2003). It is appreciated for meat quality and sport fishing, as it is very voracious at the time of capture. Due to the beauty of its color, it is also the main tourist attraction for underwater observation activities in the clear river waters of the Bonito region (Sabino & Andrade, 2003). Although this species is not vulnerable, the impact of visitors on the ichthyofauna has increased in the Bonito region, leading to a reduction in species richness (Sabino & Andrade, 2003).
Brazil has about 21% of all freshwater fish species in the world with around two thousand catalogued species (Agostinho et al., 2005; Buckup & Menezes, 2003). However, according to a Brazilian environmental agency (MMA, 2004), 134 of these species were considered endangered, including members of the genus Brycon. This loss of biodiversity is directly related to the impact on aquatic continental ecosystems, including pollution, eutrophication, silting, the construction of dams, fisheries and the introduction of exotic species (Agostinho et al., 2005).
According to the previous study (Sanches & Galetti, 2007), in which significant differences were detected between population subunits of B. hilarii, and assuming that this structuring reflects the reproductive moment of these fish, the present study surveyed the genetic variation of a B. hilarii spawning school and local populations from different collection sites in the Miranda River basin using microsatellites. The goals of this study were to identify the genetic composition of a spawning school and local populations and examine the genetic differentiation between these samples, providing useful genetic data for knowledge on biology, ecology and conservation of this fish.
MATERIALS AND METHODS
The collection was performed between 2000 and 2004. A total of 128 individuals of B. hilarii were obtained at three small streams of the Miranda River basin pertaining to the larger Paraguay River basin (MS, Brazil): Formoso River (Fo; n = 30), Peixe River (Pe; n = 19) and the main channel of the Miranda River (Mi; n = 19) (Fig. 1). A spawning school (Sch) was sampled (n = 62) in the main channel of the Miranda River. Fo and Pe are independent tributaries and both flow into the Miranda River.
Figure 1. a) Location of the Paraná-Paraguay River Basin in South America. The small rectangle delimits the Miranda River basin, b) study areas in the Miranda River basin: Formoso River-Fo (21°4'37.2"S, 56°21' 43.9"W); Peixe River-Pe (20°51'45"S, 56°30'51"W); Miranda River-Mi (20°43'26.1"S, 56°8'3.5"W). The reproductive school (Sch) was collected at the Miranda River (Mi).
Figura 1. a) Ubicación de la cuenca del río Paraná-Paraguai en América del Sur. El pequeño rectángulo delimita la cuenca del río Miranda, b) áreas de estudio en la cuenca del río Miranda: río Formoso-Fo (21°4'37,2"S, 56°21'43,9"W); río Peixe-Pe (20°51' 45"S, 56°30' 51"W); río Miranda-Mi (20°43'26,1"S, 56°8'3,5"W). El cardumen reproductivo (Sch) fue obtenido en el río Miranda (Mi).
Blood samples were collected from the caudal peduncle of the fish using hypodermic syringes containing 0.5 M EDTA (0.5 mL blood: 0.1 mL EDTA). Before releasing the fish back into the capture sites, all specimens were marked with nylon strings tied to the dorsal fin base in order to prevent recapture.
Total genomic DNA was extracted from blood samples using a saline buffer method and precipitation with ethanol (Lahiri & Nurnberger, 1991). Genetic variability was analyzed at seven microsatellite loci, Bh5, Bh6, Bh8, Bh13, Bh15, Bh16, Bh17, previously reported (Sanches & Galetti, 2006). PCRs were carried out using the method described by Schuelke (2000) in a final volume of 15 μL, containing 50 ng of DNA, 0.2 mM dNTPs, 1 x PCR buffer (20 mM Tris-HCl, pH 8.4 and 50 mM KCl; LGC Biotecnologia), 8 pmol of each reverse primer and 6-FAM, HEX or NED M13 (-21) fluorescent labeled primer as well as 2 pmol of the forward primer, 1.5 mM MgCl2 and 0.5 U of Taq DNA Polymerase (LGC Biotecnologia). PCR conditions for all loci were as follows: 1 cycle at 95°C (5 min), 35 cycles at 94°C (30 s), 56°C (45 s) and 72° C (45 s), followed by 12 cycles at 94°C (30 s), 53°C (45 s), 72°C (45 s) and a final extension at 72°C for 10 min. The microsatellite loci were analyzed on a MegaBACE 1000 automated sequencer and the alleles were scored with the Genetic Profiler (Amersham Biosciences).
Within population genetic variability, as measured by the number of alleles (NA), inbreeding coefficient (Weir & Cockerham, 1984) (FIS), observed hetero-zygosity (ZZ0) and expected (He) heterozygosity were obtained using the GENEPOP software (Raymond & Rousset, 1995). Allelic richness (Petit et al., 1998) was calculated using the FSTAT software (Goudet, 1995). Significant differences in heterozygosity and allelic richness values were evaluated between sites using the non-parametric Kruskal-Wallis test through BIOESTAT software (Ayres et al., 2003).
Exact-significance testing methods (Guo & Thompson, 1992), implemented in GENEPOP, were used to analyze deviations from the Hardy-Weinberg equilibrium (10,000 batches/10,000 iterations). Evidence for linkage disequilibrium for all pairs of loci was investigated through the Markov Chain Monte Carlo method in GENEPOP. The presence of null allele was investigated for each locus using Micro-Checker (Van Oosterhout et al., 2004).
A hierarchical partioning of the genetic diversity was carried out using an analysis of variance framework as implemented in software Arlequin 2.0 (Schneider et al., 2000). Population differentiation also was estimated using FST (Weir & Cockerham, 1984) obtained with the Arlequin. Pairwise population differentiation in genotype and gene frequencies were assessed using the exact test implemented in GENEPOP. The critical values for multiple comparisons were adjusted using the Benjamini and Yekutieli correction (Benjamini & Yekutieli, 2001). An analysis of isolation by distance was performed with the Mantel test using GENEPOP considering only the local populations.
Population structure was assessed using a model-based Bayesian procedure implemented on the Structure 2.3.3 program (Pritchard et al., 2000, 2007). This analysis was carried out assuming the admixture model and non correlated allele frequencies. Three individual repetitions of each K estimate (1-8) were run (500,000 iterations, burn-in of 200,000 iterations).
Among the relatedness estimators there is not one estimator that overpowers the others (Van de Casteele et al., 2001; Csillery et al., 2006). We used the most widely used relatedness index (RQG), (Queller & Goodnight 1989; Van Horn et al., 2008) obtained for all possible pair of individuals of each population by means of GeneAlEx 6 (Peatkall & Smouse, 2006). This relatedness index measures the extent to which two individuals share alleles that are identical by descent (Queller & Goodnight, 1989). We used Kruskal-Wallis test implemented in SPSS 13.0 to compare average relatedness obtained for each identified genetic population.
A total of 57 alleles (mean of 40.75 per population) were identified over all loci, ranging from 4 (Bh16) to 11 (Bh8 and Bh17) alleles per locus (Table 1). No evidence of linkage disequilibrium (P > 0.05) was detected for any pair of loci. Sixteen alleles were exclusive to a single population, four of which were found in Fo, one in Pe, two in Mi and nine in the Sch population. It is important to state that one of the Sch-exclusive alleles was the second most frequent allele among these fish (Bh8/ allele 227 bp/ frequency = 0.200).
Table 1. Analysis on a locus-by-locus basis of four populations of B. hilarii using seven microsatellite loci. Number of alleles per locus (Na), allelic richness (Ra), observed heterozygosity (HO), expected heterozygosity (HE), inbreeding coefficient (Fis) and P-value for departures from Hardy-Weinberg expectations. (*) Significant P-value <0.05.
Tabla 1. Análisis por locus de las cuatro poblaciones de B. hilarii utilizando siete loci de microsatélites. Número de alelos por locus (NA), riqueza alélica (Ra), heterocigosidad observada (HO), heterocigosidad esperada (HE), coeficiente de endogamia (Fis) y valores de P para el análisis de las expectativas de Hardy-Weinberg. (*) Valor significativo de P <0,05.
Two populations (Fo and Sch) exhibited a departure from Hardy-Weinberg expectations with a deficit of heterozigosity (FIS significantly higher than zero) when the data was analyzed across all loci (Table 2). The analysis on a locus-by-locus basis demonstrated that the loci Bh6 and Bh13 exhibited departures from HWE in two populations (Fo and Mi, Fo and Sch, respectively) (Appendix). According to the Micro-Checker analysis, both these loci (Bh6 and Bh13) demonstrated evidence for null alleles.
Table 2. Number of alleles (NA), allelic richness (Ra), observed heterozygosity (H0), expected heterozygosity (HE), inbreeding coefficient (Fis), P-values for departures from Hardy-Weinberg expectations (Phw) and average relatedness (Rqg) found in four populations of B. hilarii using seven microsatellite loci.
Tabla 2. Número de alelos (Na), riqueza alélica (Ra), heterocigosidad observada (HO), heterocigosidad esperada (HE), coeficiente de endogamia (Fis), valores de P para el análisis de las expectativas de Hardy-Weinberg (Phw) y índice de parentesco promedio (Rqg) encontrados para B. hilarii utilizando siete loci de microsatélites.
The populations studied showed similar genetic variation, as measured by heterozygosity, gene diversity and allelic richness values (Table 2), and no significant differences were found among them (Kruskal-Wallis, P = 1). The allele richness ranged from 5.20 to 5.66 and the gene diversity from 0.66 to 0.68. The observed and expected heterozygosities ranged of 0.59-0.61 and 0.66-0.67, respectively.
The analysis of molecular variance (AMOVA) revealed a significant genetic differentiation between populations (Fst = 0.0184; P = 0.00032). Of the total genetic variation, 1.84% was attributed to interpopulation divergence and 98.16% to the individual differences within populations (Table 3). Pair-wise analysis identified significant differentiation as revealed by the distribution of gene and genotypes frequencies (Table 4). Fst values were low, but demonstrated significant divergence between populations (Table 4). These analyses pointed to a same result, in which the spawning school exhibited significant difference when compared to any supposed local population. Considering the local populations, Fo and Pe presented significant genetic differentiation. On the other hand, when the 128 genotypes were submitted to Bayesian clustering analysis, one cluster were identified with a higher probability (K = 1; P(K/X) = 0.99). The Mantel test identified no correlation between the genetic and geographic distance of the local populations.
Table 3. Analysis of molecular variance (AMOVA) between four populations of Brycon hilarii. The data present the degrees of freedom (df), sum of squared deviation (SSD), variance component estimates and partioning of the total genetic diversity.
Tabla 3. Análisis de varianza molecular (AMOVA) entre cuatro poblaciones de Brycon hilarii. Los datos muestran los grados de libertad (df), suma de desviaciones al cuadrado (SSD), estimaciones de la varianza y repartición de la diversidad genética total.
Table 4. Genetic differentiation between four populations of Brycon hilarii. P-values from exact test of genotype and gene differentiation; FST and respective P-values.
Tabla 4. Diferenciación genética entre cuatro poblaciones de Brycon hilarii. Valores de P de la prueba exacta de la diferenciación génica y genotípica; FST y valores correspondientes de P.
The average intrapopulation relatedness (Table 2) were significantly different among the populations (Kruskal-Wallis test H = 19.6404; P = 0.002). Fish from the school were more closely related as demonstrated by the higher average relatedness value (Dunn post-test, P < 0.001).
A number of studies have reported genetic structuring of Neotropical freshwater migratory fish (Wasko & Galetti, 2002; Hatanaka et al., 2006), but the factors that facilitate this structuring are poorly known. The microsatellites used here revealed a significant genetic differentiation between the spawning school and the supposed local populations sampled in the Miranda River Basin. Moreover the two local populations Fo and Pe also presented significant genetic divergence. This same result was not achieved with the Bayesian analysis, but in cases where genetic differentiation is more subtle, testing for allelic frequency differences between groups (FST) can be more powerful than applying Structure (Pritchard et al., 2007).
In the previous study with B. hilarii, RAPD markers also revealed genetic differences between the same sites of the present study for fish collected during the feeding season (Sanches & Galetti, 2007), corroborating with the differentiation detected between Fo and Pe. These populations can comprise resident fish group into these two independent tributaries of the Miranda River Basin. Miranda River with 43.787 km2 of drainage area is one of the main rivers of the Pantanal, where there are numerous marginal and nutrient-rich lagoons that are ideal for feeding (SEMA/MS, 1994). In this way, its large size and the availability of food resources may facilitate the accumulation of several populations minimizing the chance of detecting genetic structuring through F statistic.
The genetic divergence found in the spawning school may be the result of behavior at reproduction time that reflects in the genetic structure of the fish populations (Laikre et al., 2005). A number of studies have demonstrated the relationship between reproductive behavior, social organization and genetic structure of fish populations, such as the importance of homing behavior in the genetic structuring of salmon populations that spawn in different sites (Stähl, 1987). Similarly, genetic distinctiveness has been found between resident and anadromous Oncorhynchus mykiss fish (Narum et al., 2004) as well as temporal differentiations in a herring species (Clupea harengus), which were attributed to sampling in different spawning waves (Jorgensen et al., 2005). Recent studies have demonstrated familiarity and kin association in fish (Ward & Hart, 2003; Fraser et al., 2005). It has been suggested that schooling with relatives could promote benefits (Hamilton, 1964; Quinn & Busack, 1985) by performing co-operative behaviours such as predator inspection (Wisenden & Smith, 1998), reduced levels of intraspecific aggression (Brown & Brown, 1993), enhancement of social learning in groups (Swaney et al., 2001), improvement of natal homing to breeding areas in migratory species (Olsen, 1989). The high relatedness value found for fish from school may be evidence of kin association or reproductive behaviors, such as homing and spawning waves.
However, life history and spawning are insufficiently known in Neotropical migratory freshwater species. In the São Francisco Basin, recent studies using telemetry were carried out with two commercially important species, P. argenteus and P. corruscans, in which migrations were investigated. Both species exhibited a dualistic migration pattern (Godinho & Kynard, 2006; Godinho et al., 2007). In P. argenteus, some fish were resident presenting a small total linear home range (<1 km) and used the same grounds for spawning and non-spawning activities. Though, most fish were migratory, with a large total linear home range (>53 km) and different reaches for spawning and foraging (Godinho & Kynard, 2006). Moreover, radio-tagged fish exhibited spawning-site homing behavior over two successive seasons and some fish also demonstrated very precise homing to non-spawning sites. In P. corruscans, the spatial distribution of the fish was greater during the non-spawning season than the spawning season (Godinho et al., 2007). These studies corroborate the genetic data when homing was addressed to explain the maintenance of population structuring observed in P. argenteus (Hatanaka & Galetti, 2003).
Studies on fish populations commonly encounter significant genetic differentiations between populations, but with low Fst values (Wirth & Bernatchez, 2001; Hatanaka et al., 2006). Weak differentiation has been especially detected in vagile species of great abundance, with wide distribution and no visible barriers to gene flow (Knutsen et al., 2003; Jorgensen et al., 2005). The degree and pattern of differentiation is related to the amount of gene flow between population units, and even when genetic structuring is identified, different populations present some exchange of individuals (Pettersson et al., 2001; Narum et al., 2004; Jorgensen et al., 2005), as was observed here in the spawning school.
Although the magnitude of Fst is still an issue of intense debate (Jorgensen et al., 2005), certainly it is advisable to consider the existence of subdivisions within a population than treating it as a single population, which could result in the depletion of genetic variation (Laikre et al., 2005).
The maintenance of genetic diversity depends on the preservation of gene flow among the populations of a system (Narum et al., 2004; Laikre et al., 2005). Therefore, the population structuring could maintain the existence of distinct genetic populations and avoid panmixia. On the other hand, gene flow ensures the introduction of variation into the genetic populations, avoiding inbreeding and the loss of genetic variation in the species as a whole.
In the present study, the population structuring of B. hilarii was demonstrated by the genetic differentiation of a spawning school in relation to the local populations and between Fo and Pe. This population structuring may be supported by behavior that promotes the maintenance of the genetic integrity of different subpopulations during the reproductive season. Before and after reproductive seasons, fish from different population units may be resident (as Fo and Pe) or mixed and spread throughout the hydrographic system (as Mi). But, the spawning school also may consist of a population unit from a different locality than the examined sites of the present study.
The understanding of the genetic structure of natural populations as well as the identification of population units and the knowledge of how these units are distributed throughout the environment are essential to management decisions and the conservation of the genetic resources (Moritz, 1994; Paetkau, 1999). Thus, data on the genetic structure of B. hilarii constitute tools of great importance contributing to the understanding of the behavior and biology of these fish as well as fishery management and species conservation programs.
The authors are grateful to the Fundação O Boticário de Proteção à Natureza (FBPN), Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES), Fundação de Amparo a Pesquisa do Estado de São Paulo (FAPESP) and Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq) for financial support.
Agostinho, A.A., S.M. Thomaz & L.C. Gomes. 2005. Conservação da biodiversidade em águas continentais do Brasil. Megadiversidade, 1: 70-78. [ Links ]
Ala-Honkola, A., A. Uddström, B. Diaz-Pauli & K. Lindström. 2010. Strong inbreeding depression in male mating behaviour in a poeciliid fish. J. Evol. Biol., 22: 1396-1406. [ Links ]
Avise, J.C., A.G. Jones, D. Walker & J.A. DeWoody. 2002. Genetic mating systems and reproductive natural histories of fishes: lessons for ecology and evolution. Annu. Rev. Genet., 36: 19-45. [ Links ]
Ayres, M., M. Ayres Jr., D.L. Ayres & A.S.S. Santos. 2003. BioEstat 3.0: aplicaçöes estatísticas nas áreas das ciências bio-médicas. CNPq/MCT, Belém, 259 pp. [ Links ]
Benjamini, Y. & D. Yekutieli. 2001. The control of false discovery rate under dependency. Ann. Statist., 29: 1165-1188. [ Links ]
Brown, G.E. & J.A. Brown. 1993. Social dynamics in salmonids fishes: do kin make better neighbours? Anim. Behav., 45: 863-871. [ Links ]
Carolsfeld, J., B. Harvey, C. Ross & A. Baer. 2003. Migratory fish of South America. World Fisheries Trust, Victoria, British Columbia, 372 pp. [ Links ]
Csillery, K., T. Johnson, D. Beraldi, T. Clutton-Brock, D. Coltman, B. Hansson, G. Spong & J.M. Pemberton. 2006. Performance of marker-based relatedness estimators in natural populations of outbred vertebrates. Genetics, 173: 2091-2101. [ Links ]
Fraser, D.J., P. Duchesne & L. Bernatchez. 2005. Migratory charr schools exhibit population and kin associations beyond juvenile stages. Mol. Ecol., 14: 3133-3146. [ Links ]
Gerking, S.D. 1959. The restricted movement of fish populations. Biol. Rev. Camb. Philos. Soc., 34: 221-242. [ Links ]
Gerlach, G. & N. Lysiak. 2006. Kin recognition and inbreeding avoidance in zebrafish, Danio rerio, is based on phenotype matching. Anim. Behav., 71: 1371-1377. [ Links ]
Godinho, A.L. & B. Kynard. 2006. Migration and spawning of radio-tagged zulega Prochilodus argenteus in a dammed Brazilian River. Trans. Am. Fish. Soc., 135: 811-824. [ Links ]
Godinho, A.L., B. Kynard & H.P. Godinho. 2007. Migration and spawning of female surubim (Pseudoplatystoma corruscans, Pimelodidae) in the São Francisco River, Brazil. Environ. Biol. Fishes, 80: 421-433. [ Links ]
Godoy, M.P. 1975. Peixes do Brasil: Ordem Characoidei. Editora Franciscana, Piracicaba, 847 pp. [ Links ]
Goudet, J. 1995. FSTAT (version 1.2): a computer program to calculate F-statistics. J. Hered., 86: 485-486. [ Links ]
Guo, S.W. & W.A. Thompson. 1992. Performing the exact test of Hardy-Weinberg proportions for multiple alleles. Biometrics, 48: 361-372. [ Links ]
Hamilton, W.D. 1964. The genetical theory of social behavior. J. Theor. Biol., 7: 1-52. [ Links ]
Hatanaka, T. & P.M. Galetti Jr. 2003. RAPD markers indicate the occurrence of structured populations in a migratory fish species. Genet. Mol. Biol., 26: 19-25. [ Links ]
Hatanaka, T., F.H. Silva & P.M. Galetti Jr. 2006. Population substructuring in a migratory freshwater fish Prochilodus argenteus (Characiformes, Prochilodontidae) from the São Francisco River. Genetica, 126: 153-159. [ Links ]
Jorgensen, H.B.H., M.M. Hansen, D. Bekkevold, D.E. Ruzzante & V. Loeschcke. 2005. Marine landscapes and population genetic structure of herring (Clupea harengus L.) in the Baltic Sea. Mol. Ecol., 14: 3219-3234. [ Links ]
Knutsen, H., P.E. Jorde, C. André & N.C. Stenseth. 2003. Fine-scaled geographical population structuring in a highly mobile marine species: the Atlantic cod. Mol. Ecol., 12: 385-394. [ Links ]
Lahiri, D.K. & J.I. Nurnberger. 1991. A rapid non-enzymatic method for the preparation of HMW DNA from blood for RFLP studies. Nucleic Acids Res., 19: 5444. [ Links ]
Laikre, L., S. Palm & N. Ryman. 2005. Genetic population structure of fishes: implications for coastal zone management. Ambio, 34: 111-119. [ Links ]
Lima, F.C.T. 2003. Subfamily Bryconinae (Characins, tetras). In: R.E. Reis, S.O. Kullander & C.J. Ferraris Jr. (eds.). Check List of the freshwater fish of south and Central America. EDIPUCRS, Porto Alegre, 742 pp. [ Links ]
Lowe-McConnell, R.H. 1987. Ecological studies in tropical fish communities. Cambridge University Press, New York, 382 pp. [ Links ]
McPherson, A.A., R.L. Stephenson & C.T. Taggart. 2003. Genetically different Atlantic herring Clupea harengus spawning waves. Mar. Ecol. Prog. Ser., 247: 303-309. [ Links ]
Ministério do Meio Ambiente (MMA). 2004. Lista nacional de espécies de invertebrados aquáticos e peixes ameaçadas de extinção. Ministério do Meio Ambiente. Diário Oficial da União, 1: 136-142. [ Links ]
Moritz, C. 1994. Defining evolutionary significant units for conservation. Trends Ecol. Evol., 9: 373-375. [ Links ]
Narum, S.R., C. Contor, A. Talbot & M.S. Powell. 2004. Genetic divergence of sympatric resident and anadromous forms of Oncorhynchus mykiss in the Walla Walla River, USA. J. Fish Biol., 65: 471-488. [ Links ]
Odling-Smee, L. & V.A. Braithwaite. 2003. The role of learning in fish orientation. Fish Fisheries, 4: 235-246. [ Links ]
Olsen, K.H. 1989. Sibling recognition in juvenile Arctic charr, Salvelinus alpinus. J. Fish Biol., 34: 571-581. [ Links ]
Peatkall, R. & P.E. Smouse. 2006. GENALEX 6: genetic analysis in Excel. Population genetic software for teaching and research. Mol. Ecol. Notes., 6: 288-295. [ Links ]
Paetkau, D. 1999. Using genetics to identify intraspecific conservation units: a critique of current methods. Conserv. Biol., 13: 1507-1509. [ Links ]
Petit, R.J., A. El Mousadik & O. Pons. 1998. Identifying populations for conservation on the basis of genetic markers. Conserv. Biol., 12: 844-855. [ Links ]
Pettersson, J.C.E., M.M. Hansen & T. Bohlin. 2001. Does dispersal from landlocked trout explain the coexistence of resident and migratory trout females in a small stream? J. Fish Biol., 58: 487-495. [ Links ]
Piorski, N.M., A. Sanches, L.F. Carvalho-Costa, T. Hatanaka, M. Carrillo-Avila, P.D. Freitas & P.M. Galetti Jr. 2008. Contribution of conservation genetics in assessing neotropical freshwater fish biodiversity. Brazil. J. Biol., 68 (Suppl.): 1039-1050. [ Links ]
Pritchard, J.K., M. Stephens & P. Donnelly. 2000. Inference of population structure using multilocus genotypic data. Genetics, 155: 945-959. [ Links ]
Pusey, A. & M. Wolf. 1996. Inbreeding avoidance in animals. Tree, 11: 201-205. [ Links ]
Queller, D.C. & K.F. Goodnight. 1989. Estimating relatedness using genetic markers. Evolution, 43: 258-275. [ Links ]
Quinn, T.P. & C.A. Busack. 1985. Chemosensory recognition of siblings in juvenile coho salmon (Oncorhynchus kisutch). Anim. Behav., 33: 51-56. [ Links ]
Raymond, M. & F. Rousset. 1995. GENEPOP (version 1.2): a population genetic software for exact test and ecumenism. J. Hered., 86: 248-249. [ Links ]
Sabino, J. & L.P. Andrade. 2003. Uso e conservação da ictiofauna no ecoturismo da região de Bonito, Mato Grosso do Sul: o mito da sustentabilidade ecológica no rio Baía Bonita (Aquário Natural de Bonito). Biota Neotrop., 3: 1-9. [ Links ]
Sanches, A. & P.M. Galetti Jr. 2006. Microsatellites loci isolated in the freshwater fish Brycon hilarii. Mol. Ecol. Notes, 6: 1045-1046. [ Links ]
Sanches, A. & P.M. Galetti Jr. 2007. Genetic evidence of population substructuring in the neotropical freshwater fish Brycon hilarii. Brazil. J. Biol., 67: 889 -895. [ Links ]
Schneider, S., D. Roessli & L. Excoffier. 2000. Arlequin ver. 2.0: a software for population genetics data analysis. Genetics and Biometry Laboratory, University of Geneva, Switzeland, 174 pp. [ Links ]
Schuelke, M. 2000. An economic method for the fluorescent labeling of PCR fragments. Nat. Biotechnol., 18: 233-234. [ Links ]
SEMA/MS. 1994. Relatório sobre a implantação da rede básica de monitoramento da qualidade das águas da bacia do Alto Paraguai. Secretaria do Meio Ambiente do Estado do Mato Grosso do Sul, Campo Grande, 127 pp. [ Links ]
Stähl, G. 1987. Genetic population structure of Atlantic salmon. In: N. Ryman & F. Utter (eds.). Population genetics and fishery management. Washington Sea Grant Program. University of Washington Press, Seattle, pp. 121-140. [ Links ]
Swaney, W., J. Kendal, H. Capon, C. Brown & K.N. Laland. 2001. Familiarity facilitates social learning of foraging behaviour in the guppy. Anim. Behav., 62: 591-598. [ Links ]
Van de Casteele, T., P. Galbusera & E. Matthysen. 2001. A comparison of microsatellite-based pairwise relatedness estimators. Mol. Ecol., 10: 1539-1549. [ Links ]
Van Horn, R.C., J. Altmann & S.C. Alberts. 2008. Can't get there from here: inferring kinship from pairwise genetic relatedness. Anim. Behav., 75: 1173-1180. [ Links ]
Van Oosterhout, C., W.F. Hutchinson, D.P.M. Wills & P. Shipley. 2004. Micro-Checker: software for identifying and correcting genotyping errors in microsatellite data. Mol. Ecol. Notes, 4: 535-538. [ Links ]
Ward, A.J.W. & P.J.B. Hart. 2003. The effects of kin and familiarity on interactions between fish. Fish Fisheries, 4: 348-358. [ Links ]
Wasko, A.P. & P.M. Galetti Jr. 2002. RAPD analysis in the Neotropical fish Brycon lundii: genetic diversity and its implications for the conservation of the species. Hydrobiologia, 474: 131-137. [ Links ]
Weir, B.S. & C.C. Cockerham. 1984. Estimating F-statistics for the analysis of population structure. Evolution, 38: 1358-1370. [ Links ]
Wisenden, B.D. & R.J.F. Smith. 1998. A re-evaluation of the effect of shoalmate familiarity on the proliferation of alarm substance cells in ostariophysan fishes. J. Fish Biol., 53: 841-846. [ Links ]
Wirth, T. & L. Bernatchez. 2001. Genetic evidence against panmixia in the European eel. Nature, 409: 1037-1040. [ Links ]
Received: 15 September 2011; Accepted: 9 June 2012
Corresponding author: Alexandra Sanches (firstname.lastname@example.org)