Genetic diversity and population structure of bocachico Prochilodus magdalenae (Pisces, Prochilodontidae) in the Magdalena River basin and its tributaries, Colombia

RESEARCH ARTICLE

Genetic diversity and population structure of bocachico Prochilodus magdalenae (Pisces, Prochilodontidae) in the Magdalena River basin and its tributaries, Colombia

Gilberto Orozco Berdugo; Juan C. Narváez Barandica

Universidad del Magdalena, Grupo de Biodiversidad y Ecología Aplicada, Laboratorio de Genética Molecular, Santa Marta, Magdalena, Colombia.

^{Send correspondence to}

ABSTRACT

Prochilodus magdalenae is an endemic freshwater fish that occurs in the Magdalena, Sinú and Atrato hydrographic basins. It has an important economic role and is a food resource for the artisanal fishing communities. Its socioeconomic importance contrasts with the current status of its fisheries, where stocks are being depleted. Considering its importance and lack of information on its genetic structure, we used seven microsatellite markers to assess the genetic structure of wild populations of P. magdalenae. The genetic diversity was assessed and the population genetic structure was estimated through Fst, analysis of molecular variance and Bayesian analysis. A total of 290 alleles were found in all loci throughout all population. The high polymorphism contrasts with the levels of observed heterozygosity (Ho = 0.276), which are the lowest values recorded for the family. We found three populations of bocachico coexisting throughout the studied system, contradicting the hypothesis that freshwater migratory fish form panmictic populations. These results on the genetic structure of P. magdalenae constitute tools for a better understanding of the behavior and biology of this species, contributing to fish management and conservation programs.

Key words: neotropical fish, microsatellite, conservation genetics, population genetics.

Introduction

The Magdalena River is the principal hydrographic system in Colombia and the major axis of economical development in the country (Jimenez-Segura et al., 2010). Its basin is intensely populated and nearly 80% of the Colombian human population lives on it (Galvis and Mojica, 2007). The basin has an area of 257,438 km2, occupying approximately 24% of the Colombian territory (Galvis and Mojica, 2007) and has a large biomass and diversity of freshwater fish, harboring 167 species (Galvis and Mojica, 2007), some of which, like bocachico (Prochilodus magdalenae), migrate in the reproductive season. The"bocachico" is an endemic species occurring in the Magdalena, Sinú and Atrato hydrographic basins (Mojica et al., 2002; Maldonado-Ocampo et al., 2005). This is an important fish resource for the artisanal fishing communities, representing one of the commonly captured species. The bocachico is a detrivorous species that feeds on the fine bottom sediment containing organic particles. Is considered a key element in the ecosystem and may play an ecological role due to its sediments processing habit (Flecker, 1996).

Like many other prochilodontids, P. magdalenae has high fecundity and spawns all eggs at once in the open waters of the main river channel (Agostinho et al., 1993). Larvae drift passively towards flooded areas where they feed and remain until maturation (Agostinho et al., 1993). Their life cycle is reliant on the hydrological patterns of the Magdalena River basin, where the fish migrate to the main river in twice-yearly hydrological periods. In the first onewhen water levels begin to decrease (December-February), spawning migration to the upstream starts and fishes remain there during the dry season. Spawning begins with the onset of the first high-water period (March-June). The second one takes place when water levels start to decrease (July-September) and a second spawning migration occurs. This spawning has a lower magnitude and the process is not very clear, but apparently relies on fishes that failed to reproduce in the previous period (Jimenéz-Segura et al., 2010).

Over the last years, genetic investigations have revealed that freshwater migratory fish can display population structuring within a single hydrographic system (Wasko and Galetti Jr, 2002; Hatanaka et al., 2006; Sanches and Galetti Jr, 2007). During the spawning season, fish schooling may exhibit a behavior that enables the maintenance of the genetic integrity of the populations (Hatanaka et al., 2006; Sanches and Galetti Jr, 2007). This means that the reproductive behavior constitutes a decisive factor in the species population subdivision (Sanches and Galetti Jr, 2007).

One approach that has been successfully applied in genetics studies on fish is the use of microsatellite markers (Piorski et al., 2008), which are short tandem repeats of 1-6 nucleotides widely distributed in the genome (Litt and Luty, 1989). Due to their abundant distribution, high degree of polymorphism, Mendelian inheritance (Selkoe and Toonen, 2006) and simple evolutionary mechanisms (Piorski et al., 2008), microsatellites have been largely used to access genetic diversity (Barroso et al., 2005; Matsumoto and Hilsdorf, 2009; Barroca et al., 2012) and population structure in fish (Barroso et al., 2005; Hatanaka et al., 2006; Carvalho-Costa et al., 2008; Calcagnotto and DeSalle, 2009; Matsumoto and Hilsdorf, 2009; Sanches et al., 2012; Barroca et al., 2012). However, few studies using microsatellite markers were carried out to detect the genetic structure of bocachico in the Magdalena River (Santacruz, 2003).

The study of genetic diversity is essential for fish conservation, which depends on the knowledge of the amount of variation existing in a local reproductive unit (Carvalho, 1993). The importance of this approach lies on its potential for delimiting priority areas for species conservation and sustainable use (Sanches and Galetti Jr, 2007). The purpose of this study was to employ seven loci of microsatellite markers to gather information on the genetic structure of Prochilodus magdalenae in the Magdalena River and its tributaries, providing data for the conservation and management planning of this fish.

Material and Methods

Sample collection

Specimens were captured with the help of artisanal fishermen at 25 sites from upstream (Neiva) to downstream (Ciénaga Grande de Santa Marta) of the Magdalena River basin and main tributaries, like Sogamoso River, Cauca River and San Jorge River, from April to December 2010 (Figure 1). Muscle tissue was removed from 759 specimens sampled immediately after capture and stored in 96% ethanol until use.

DNA extraction and amplification of loci

The genomic DNA was extracted from P. magdalenae muscle using the MasterPure kit (Epicentre Biotechnologies®). Genetic diversity was analyzed using seven microsatellite loci described for Prochilodus lineatus (PL3, PL14, PL23, PL28, PL34, PL64 and PL119) that cross-amplify in P. magdalenae (Rueda et al., 2011). PCRs were carried out in a final volume of 10 µL containing 100 ng of DNA, 1Xreaction buffer (20mMTris-HCl pH 8; 50 mM KCl), 2 mM MgCl2, 200 µM dNTPs, 0.2 µM of each primer (forward and reverse) and 0.25 U of Taq polymerase (Bioline Meridian Life Science). PCR conditions were as follows: 5 min at 95 °C, 30 cycles of 30 s at 94 °C, 30 s at the annealing temperature, 30 s at 72 °C, and a final extension at 72 °C for 10 min. The reactions were performed with the ESCO-SWIF MaxPro gradient thermocycler. The amplification products were analyzed by capillary electrophoresis QIAxcel Advance (QIAGEN), using a high resolution kit (High Resolution Kit QIAGEN) and a molecular weight ladder of known concentration (DNA Size Marker 50-800 bp v2.0 QIAGEN). The size of each amplified product was determined with the ScreenGel QIAxcel v1.0 QIAGEN program, which allows quantifying the weight of each band thus distinguishing heterozygous from homozygous individuals.

Microsatellite statistical analysis

The genetic diversity was estimated for each population analyzing the number of alleles per locus (NA), observed heterozygosity (Ho), and expected heterozygosity (He) computed using the GenAlex 6.0 software (Peakall and Smouse, 2006). The inbreeding coefficients (Fis) perloci were calculated in FSTAT (Goudet, 1995). Departure from Hardy-Weinberg expectations was calculated using a test analogous to Fisher's exact test (Guo and Thompson, 1992) that has as null hypothesis the random union of gametes. This test was estimated with a Markov Chain Monte Carlo (MCMC) series of permutations (10,000 batches/1000 iterations), implemented in GENEPOP (Raymond and Rousset, 1995). The independent segregation of genotypes (linkage disequilibrium) was also tested using the exact test with a MCMC (10,000 batches/1000 iterations) (Guo and Thompson, 1992) facilitated also by GENEPOP. The presence of genotyping errors arising from technical artifacts, namely null alleles or large allele dropout due to poor DNA quality was assessed with the software MICRO-CHECKER (van Oosterhout et al., 2004).

For population structure analysis the ARLEQUIN 2000 (Schneider et al., 2000) program was employed to calculate pairwise FST value (Weir and Cockerham, 1984) estimates for each population pair. We examined the distribution of genetic variability into hierarchical levels at different geographical scales in the Magdalena River and tributaries through AMOVA (Analysis of Molecular Variance) performed in ARLEQUIN 2000. In this case, the hierarchical levels were: 1) all populations formed a single group and 2) a group from all localities of the principal channel in the Magdalena River without the tributaries like Sogamoso, Cauca and San Jorge River. Thus, AMOVA provided Φ ST values analogous to FST of Wright (1978), and through a non-parametric procedure of random permutations (10,000), the initial hypothesis of no genetic structuring between populations was tested.

The Bayesian clustering methodology of STRUCTURE version 2.3.3 (Hubisz et al., 2009) was applied to identify clusters of genetically similar individuals and determine the level of genetic substructure in the data set independently of sampling areas. To estimate the number of subpopulations (K), three independent runs of K = 1-15 were carried out at 100,000 Markov Chain Monte Carlo repetitions with a 100,000 burn-in period using no prior information and assuming correlated allele frequencies and the admixture model. The migration pattern of P. magdalenae involves the existence of high gene flow between populations; so we used an admixture model where each individual is assumed to have inherited some proportion of its ancestry from each population. To determine the number of populations (K) present in the Magdalena River basin and its tributaries, we used the method proposed by Evanno et al. (2005). This value was obtained using the software STRUCTURE HARVESTER 0.56.3 (Earl, 2009).

Results

Genetic diversity

All seven microsatellites were polymorphic (100%) in the wild populations of P. magdalenae within the Magdalena River basin and its tributaries. A total of 290 alleles were found in all loci throughout the whole population ranging from 33 to 59 by PL64 and PL119. Average number of alleles per locus in the population was 41.4± 9.6 and the mean in each sampling site ranged from 9.57 (Ciénaga de Canta Gallo) to 18.43 alleles (Ciénaga de Pijiño) (Table 1).

From the 290 alleles found in the population of bocachico of the Magdalena River basin, 61 were private alleles with a frequency below 11%. Such alleles were detected in 21 of the 25 sampling sites and ranged from one in the Barranco de Loba, Sogamoso River, Neiva and Gambote to eight in Gamarra (Table 1). The observed heterozygosity (Ho) was low in the wild population of bocachico and ranged between 0.19 (Cauca River) to 0.33 (Neiva) (Table 1). The highest expected heterozygosity (He) was found in Puerto Boyacá (0.9188), and the lowest in Ciénaga de Canta Gallo (0.8289). The average FIS (inbreeding coefficient) was highest in all sampled sites in the Magdalena River basin and ranged from 0.624 (Ciénaga San Silvestre) to 0.777 (Ciénaga de Pajarales).

All sites presented departures from the Hardy-Weinberg expectations (HWE) for all loci (p < 0.01). The U test (Raymond and Rousset, 1995) showed that the deviation from HWE is due to the heterozygote deficit (p < 0.0001). This deficit was also revealed by the positive Fis values. The occurrence of null alleles was identified in all loci and in all localities in the Magdalena River basin through tests performed on MICRO-CHECKER.

No significant association among the genotypes of seven loci presented linkage disequilibrium, suggesting that the analyses could be performed assuming statistical independence of the loci.

Population differentiation and structure

The Fst test revealed presence and absence of differentiation among all pairwise population comparisons, ranging from -0.00038 between the localities of Mompós and Puerto Boyacá to 0.15 in the localities of Ciénaga de Canta Gallo and Sogamoso River (Table 2). Variation among and within populations was assessed by AMOVA using ARLEQUIN 2000. To test the genetic structure, the sampling sites were artificially divided into the different hierarchical groups described above. The first analysis revealed a significant genetic differentiation between all populations in the Magdalena River basin and its tributaries (Φ_ST = 0.06511; p < 0.00) showing that nearly 6.5% of the total microsatellite DNA diversity was explained by variability among populations and 93.49% was found within population (Table 3). When the tributaries were excluded from the analysis, the AMOVA continued to show a significant genetic differentiation (Φ_ST = 0.0638; p < 0.00) revealing that nearly 6.38% of the variation was explained by the variability between populations and not within them.

However, the AMOVA and FST test require a priori definition of subpopulations or localities sampled, and this grouping may not correspond to biological groups or populations. So we used the STRUCTURE 2.3.3 program to verify if individuals were clustered into two or more populations when their genotypes indicated that they were admixed. The Bayesian analysis provided the most probable number of subpopulations of P. magdalenae in the Magdalena basin through the method proposed by Evanno et al. (2005) and indicated at least three clusters or populations (K = 3) genetically differentiated within the data set (Figure 2). These populations are distributed along the basin, where at least one fish of each sampled site was assigned to one cluster (red, blue or green).

Discussion

Genetic diversity in wild populations

Despite the large biodiversity of freshwater fishes in South America, there are few studies using microsatellite markers to assess the genetic characteristics of migratory fish populations (Oliveira et al., 2009). The present study demonstrates that P. magdalenae has a high genetic diversity, He = 0.877, compared to that observed for other species of the genus (Santacruz, 2003; Hatanaka et al., 2006; Galzerani, 2007; Carvalho-Costa et al., 2008; Silva, 2011; Rueda et al., 2011) and to other commercially important species that have the same reproductive pattern and migrate long distances (Sanches, 2007; Batista, 2010; Dantas, 2010) (Table 4). Likewise, the microsatellite genetic diversity was slightly lower than that reported by Aguirre et al. (2013) for mitochondrial DNA (mtDNA control region) in P. magdalenae (HD = 0.997). These results demonstrate the vagility of the P. magdalenae, in which the strong migratory behavior facilitates the maintenance of high levels of genetic variability, as in species of the genus Prochilodus (Lassala and Renesto, 2007; Santos et al., 2007).

On the other hand, the observed heterozygosity of the P. magdalenae populations are not comparable to those registered for others species of fishes (Santacruz, 2003; Hatanaka et al., 2006; Galzerani, 2007; Sanches, 2007; Carvalho-Costa et al., 2008; Silva, 2011; Rueda et al., 2011), because these values are among the lowest recorded for the whole family (Table 4). The present study demonstrates that all loci exhibited heterozygosity deficiency throughout the populations surveyed and the presence of null alleles is suggested as the cause of the deficiency. This has also been the explanation for the heterozygote deficiency proposed by Barroso et al. (2005) and Matsumoto and Hilsdorf (2009). However, several factors can originate heterozygote deficiency, such as: selection on a specific locus, inbreeding and cryptic population structure (Garcia De Leon et al., 1997). The first factor was not considered since we did not observe loci under selection (Slatkin, 1995). In the freshwater fish Prochilodus argenteus, a heterozygote deficit was attributed to a combination of random sampling effects and null alleles (Hatanaka et al., 2006). Thus, the heterozygote deficit can not be explained by a single factor, since the interaction of several factors may be contributing to it (Sanches, 2007). The microsatellite used in this study allowed the estimation of a high genetic diversity in the P. magdalenae population, despite the presence of some unique alleles. Due to the great number of individuals sampled in this study (759), it is likely that these unique alleles are rare. If this is so, they can be used as population markers, powerful to show genetic flow (Slatkin, 1985).

Population structure

The microsatellite data suggest the presence of a significant population structuring in the migratory fish P. magdalenae from the Magdalena River basin and its tributaries revealed by the FST and AMOVA statistics. This suggests that these fish organize themselves during the spawning period in a way that maintains the integrity of each subunit of the system (Sanches and Galetti Jr, 2007). Pereira et al. (2009) studied the population structure of Pseudoplatystoma corruscans and argued that, due to the magnitude and geographical scale of their study, the genetic differentiation found could only be explained by homing behavior. However, it is important to note that their samples were collected in tributaries during their spawning migration movements to avoid collecting mixed populations in their feeding areas. This prevents comparisons between our work and that of Pereira et al. (2009) as we included samples from the main canals of the Magdalena River and their feeding areas (like lagoon and marshes), making the structuring pattern of P. magdalenae unclear. Further studies with samples from headwaters of the tributary rivers could clarify whether P. magdalenae presents homing behavior during spawning migration. Thus, the genetic differentiation found in the population of P. magdalenae is not consistent with the currently accepted idea of panmictic populations of neotropical migratory freshwater fish (Sivasundar et al., 2001; Castro and Vari, 2004; Aguirre et al., 2013), despite the enormous vagility and long migrations of these fish.

Similar trends have been reported in Prochilodus argenteus in the São Francisco River (Brazil) with microsatellite markers, (Hatanaka et al., 2006). A similar situation was also reported for P. costatus and P. argenteus studied in the same hydrographic basin (Barroca et al., 2012). In other species with similar reproductive and migratory behaviors, like Brycon orthotaenia, two defined population were identified in the São Francisco River (Sanches et al., 2012). It has also been found that B. insignis are currently structured possibly due to anthropogenic actions (Matsumoto and Hilsdorf, 2009). A similar situation was reported in P. magdalenae in the Sinú River (Colombia) (Santacruz, 2003), in which a significant genetic structure was identified with heterologous microsatellite markers. However, population studies of P. magdalenae from the Magdalena River basin with the mtDNA control region indicated a single panmictic population (Aguirre et al., 2013). The difference observed with the two markers may be explained by the recent impact of human activities that probably has not allowed enough time to leave traces in the mtDNA control region.

The structuring pattern was also clearly demonstrated in the Bayesian analysis, in which three populations were identified. These populations are distributed along the Magdalena River basin, assuming that at least one individual collected from each sampled location was assigned to one of three clusters, where each location represent a mixture of different populations. This result demonstrates the dispersal capabilities of bocachico due to its reproductive behavior and the geographic proximity of some localities, where many of the sampled sites (like lagoon and marshes) are interconnected during periods of flooding. Thus, our hypothesis is that the population genetic structure of P. magdalenae may be maintained by an event of"reproductive waves" (Jorgensen et al., 2005), represented by genetically differentiated groups that breed in the same place at different time periods with some overlap. This idea is supported by the fact that P. magdalenae has two reproductive peaks during the year in the Magdalena River basin (Valderrama, 1972; Valderrama-Barco and Petrere Jr, 1994; Jimenez-Segura et al., 2010). These authors suggested that a fraction of the population reproduces during the first hydrological cycle and the remaining in the second one. This second migration pulse may involve individuals that could not migrate during the first hydrological cycle.

However, this "reproductive waves" behavior would imply the existence of two populations in the Magdalena River basin. The third population detected by the Bayesian analysis could be a consequence of the repeated restocking programs conducted in this basin. Machado-Schiaffino et al. (2007) provided evidence of significant genetic variation losses in Atlantic salmon stocks (Salmo salar) created for supportive breeding in which the juveniles released in the rivers possessed significantly lower allelic richness than the wild stocks. The cultivated population of Brycon opalinus also presented heterozygosity reduction, indicating a loss of genetic variability in the reproductive supply currently kept in the hatchery (Barroso et al., 2005).Asimilar trend was reported by Matsumoto and Hilsdorf (2009), in which the broodstock of B. insignis kept at the hatchery has likely maintained the genetic diversity formerly present in some rivers and no longer existing in natural populations. Further studies may clarify whether the broodstock used in stocking programs has influenced the genetic structure and diversity of the P. magdalenae population in the Magdalena basin.

Implications for conservation of P. magdalenae

The IUCN recognizes the need to conserve genetic diversity as one of three global priorities for biodiversity conservation focused on those species vulnerableand at risk of over exploitation like P. magdalenae, which is considered vulnerable (Mojica et al., 2012). Therefore, the main goal of this study was to improve the genetic information of P. magdalenae population structure for fishery management and conservation purposes. Despite the high levels of genetic diversity, the current distribution of P. magdalenae is discontinued or heterogeneous, with at least three populations distributed along the Magdalena River basin. This suggests that management and conservation strategies for P. magdalenae should aim at preserving the diversity of each population.

In this scenario, the recovery of P. magdalenae depends on the strategies implemented for conservation or even on reversing the habitat degradation of the river and riparian environments. On the other hand, other measures could be implemented, such as restocking programs in areas where the population has a lower genetic variability. In this case, the genetic divergence found in the P. magdalenae populations should be taken into account and appropriate breeding management aimed at reducing the risks of genetic drift, inbreeding, and the bottleneck effect should be implemented (Matsumoto and Hilsdorf, 2009). The genetic monitoring of the broodstock and juveniles used for supportive breeding is also essential as it will result in more efficient management and conservation strategies (Barroso et al., 2005; Machado-Schiaffino et al., 2007; Lopera Barrero et al., 2009; Matsumoto and Hilsdorf, 2009).

Acknowledgments

The authors would like to thank to COLCIENCIAS and the Universidad Del Magdalena for financial support through the projects Establecimiento de Sistemas de Parentales con Criterios Genéticos Para Producir Semillas de Bocachico Prochilodus magdalenae en el norte de Colombia. (Number: 1117-489-25459) and Evaluación de la Ecología Molecular de los Bocachicos (Prochilodus spp) Asociados a los Ríos que Drenan al Caribe Colombiano (Number: 1117-521-28352). We also thank J. Aguirre, E. Muñoz, Y. Julio, D. Castañeda and A. Torregroza for their support in the collections and in the laboratory.

References

Aguirre JC, Narváez JC and Castro LR (2013) Mitochondrial DNA variation of the bocachico Prochilodus magdalenae (Characiformes, Prochilodontidae) in the Magdalena River Basin, Colombia. Aquat Conserv 23:594-605.

Agostinho AA, Vazzoler AEA, Gomes LC and Okada EK (1993) Estratificação e comportamento de Prochilodus scrofa em distintas fases de ciclo de vida, em lá planície de inundação del alto rio Paraná y embalse de Itaipu, Paraná, Brasil. Rev Hydrobiol Trop 26:79-90.

Barroca T, Arantes F, Magalhaes B, Siqueira F, Horta C, Pena I, Dergam J and Kalapothakis E (2012) Genetic diversity and population structure of Prochilodus costatus and Prochilodus argenteus preceding dam construction in the Paraopeba River, São Francisco River Basin, Minas Gerais, Brazil. Open J Genet 2:121-130.

Barroso RM, Hilsdorf AWS, Moreira HLM, Cabello pH and Traub-Cseko YM (2005) Genetic diversity of wild and cultured populations of Brycon opalinus (Cuvier, 1819) (Characiforme, Characidae, Bryconiae) using microsatellites. Aquaculture 247:51-65.

Batista JS (2010) Caracterização genética da dourada - Brachyplatystoma rousseauxii, Castelnau, 1855 (Siluriformes, Pimelodidae) na Amazônia por meio de marcadores moleculares mitocondriais e microssatélites: Subsídios para conservação e manejo. Doctoral Thesis, Instituto Nacional de Pesquisas da Amazônia, Manaus, Brazil.

Calcagnotto D and DeSalle R (2009) Population genetics structuring in Pacu (Piaractus mesopotamicus) across the Parana- Paraguay Basin: Evidence from microsatellite. Neotrop Ichthyol 7:607-616.

Carvalho GR (1993) Evolutionary aspects of fish distribution: Genetic variability and adaptation. J Fish Biol 43 (Suppl A):53-73.

Carvalho-Costa LF, Hatanaka T and Galetti Jr PM (2008) Evidence of lack of population substructuring in the Brazilian freshwater fish Prochilodus costatus. Genet Mol Biol 31 (Suppl 1):377-380.

Castro R and Vari R (2004) Detritivores of South American fish family Prochilodontidae (Teleostei, Ostariophysi, Characiformes) a phylogenetic and revisionary study. Smithson. Contrib Zool 622:1-189.

Dantas HL (2010). Avaliação da estrutura genética do surubim, Pseudoplatystoma corruscans (Actinopterigy, Siluriformes) como subsidio para o repovoamento do submédio São Francisco. Masters Thesis Universidade Federal Rural do Pernambuco. Recife, Brazil.

Evanno G, Regnaut S and Goudet J (2005) Detecting the number of cluster of individuals using the software structure: A simulation study. Mol Ecol 14:2611-2620.

Flecker AS (1996). Ecosystem engineering by a dominant detritivore in a diverse tropical stream. Ecology 77:845-854.

Galvis G and Mojica J (2007) The Magdalena river fresh water fishes and fisheries. Aquat Ecosyst Health 10:127-139.

Galzerani F (2007) Análise da variabilidade genética de Prochilodus argenteus (Pisces, Prochilodontidae) do rio São Francisco, região de Três Marias, através de marcadores microsatélites. Monography, Dept. Genetics and Evolution, Universidade Federal de São Carlos, São Carlos, Brazil.

Garcia De Leon FJ, Chikhi L and Bonhomme F (1997) Microsatellite polymorphism and population subdivision in natural populations of European sea bass Dicentrarchus labrax (Linneo 1758) Mol Ecol 6:51-62.

Goudet J (1995). FSTAT: A computer program to calculate F-statistics. J Hered 86:485-486.

Guo SW and Thompson EA (1992) Performing the exact test of Hardy-Weinberg proportion for multiple alleles. Biometrics 48:361-372.

Hatanaka T, Silva FH and Galetti Jr PM (2006) Population substructuring in a migratory freshwater fish Prochilodus argenteus (Characiformes, Prochilodontidae) from the São Francisco River. Genetica 126:153-159.

Hubisz MJ, Falush D, Stephens M and Pritchard JK (2009) Inferring weak population structure with the assistance of sample group information. Mol Ecol Resourc 9:1322-1332.

Jimenez-Segura L, Palacio J and Leite R (2010) River flooding and reproduction of migratory fish species in the Magdalena river basin, Colombia. Ecol Freswather Fish 19:178-186.

Jorgensen H, Hansen M, Bekkevold D, Ruzzante DE and Loeschcke V (2005) Marine landscape and population genetic structure of herring (Clupea harengus L.) in the Baltic Sea. Mol Ecol 14:3219- 3234.

Lassala MDP and Renesto E (2007) Reproductive strategies and genetic variability in tropical freshwater fish. Genet Mol Biol 30:690-697.

Litt M and Luty JA (1989) A hypervariable microsatellite revelead by in vitro amplification of a dinucleotide repeat within the cardiac muscle actin gene. Am J Hum Genet 44:397-401.

Lopera Barrero N, Ribeiro R, Povh J, Vargas L, Jacometo C and Gomes P (2009) Genetic diversity in Piaractus mesopotamicus stocks used in stock enhancement programs and implications for conservation. Agrociencia 43:249-256.

Machado-Schiaffino G, Dopico E and Garcia-Vazquez E (2007) Genetic variation losses in Atlantic salmon stocks created for supportive breeding. Aquaculture 264:59-65.

Maldonado-Ocampo JA, Ortega-Lara A, Usma O, Galvis V, Villa- Navarro FA, Vasquez G, Prada-Pedreros S and Ardila R (2005) Peces de los Andes de Colombia. Instituto de Investigaciones y Recursos Biológicos Alexander von Humboldt, Bogotá DC, 346 pp.

Matsumoto CK and Hilsdorf AS (2009) Microsatellite variation and population genetic structure of a Neotropical endangered Bryconinae species Brycon insignis Steindachner, 1877: Implications for its conservation sustainable management. Neotrop Ichthyol 7:395-402.

Mojica J, Castellanos C, Usma J and Álvarez R (2002) Libro rojo de peces dulceacuícolas de Colombia. Serie Libros Rojos de Especies Amenazadas de Colombia. Universidad Nacional de Colombia y Ministerio del Medio Ambiente, Bogotá, 288 pp.

Mojica J, Castellanos C, Usma J, Álvarez R and Lasso C (2012) Libro rojo de peces dulceacuícolas de Colombia. Serie Libros Rojos de Especies Amenazadas de Colombia. Universidad Nacional de Colombia, WWF Colombia y Universidad de Manizales, Bogotá, 319 pp.

Oliveira C, Foresti F and Hilsdorf AWS (2009) Genetics of neotropical fish: From chromosomes to populations. Fish Physiol Biochem 35:81-100.

Peakall R and Smouse P (2006) GenAlEx 6: Genetic analysis in Excel. Population genetic software for teaching and research. Mol Ecol Notes 6:288-295.

Pereira LHG, Foresti F and Oliveira C (2009) Genetic structure of the migratory catfish Pseudoplatystoma corruscans (Siluriformes, Pimelodidae) suggests homing behaviour. Ecol Freshwater Fish 18:215-225.

Piorski NM, Sanches A, Carvalho-Costa LF, Hatanaka T, Carrillo-Avila M, Freitas PD and Galetti PM (2008) Contribution of conservation genetics in assessing neotropical freshwater fish biodiversity. Braz J Biol 68(Suppl):1039- 1050.

Rueda E, Sommer J, Scarabotti P, Markariani R and Ortí G (2011) Isolation and characterization of polymorphic microsatellite loci in the migratory freshwater fish Prochilodus lineatus (Characiformes, Prochilodontidae). Conserv Genet Resource 3:681-684.

Raymond M and Rousset F (1995) GENEPOP v. 1.2.: Population genetics software for exact test and ecumenicism. J Hered 86:248-249.

Sanches A (2007) Estrutura genética populacional de Brycon hilarii na sub-bacia do Rio Miranda, e seu significado para programas de conservação. Doctoral Thesis, UFSCar, São Carlos, Brazil. Sanches A and Galetti Jr PM (2007) Genetic evidence of population substructuring in the freshwater fish Brycon hilarii. Braz J Biol 67:889-895.

Sanches A, Galetti Jr P, Galzerani F, Derazo J, Cutilak-Bianchi B and Hatanaka T (2012) Genetic population structure of two migratory freshwater fish species (Brycon orthotaenia and Prochilodus argenteus) from the São Francisco River in Brazil and its significance for conservation. Latin Am J Aquat Res 40:177-186.

Santacruz BDH (2003) Evaluación de la variabilidad genética con marcadores microsatélites del bocachico Prochilodus magdalenae (Steindachner 1878) en el Río Sinú, Colombia. Thesis, Universidad Nacional de Colombia, Bogotá, 126 pp.

Santos MCF, Ruffino ML and Farias IP (2007) High levels of genetic variability and panmixia of the tambaqui Colossoma macropomum (Cuvier, 1816) in the main channel of theAmazon River. J Fish Biol 71(A):33-44.

Schneider S, Roessli D and Excoffier L (2000). Arlequin: A software for population genetics data analysis. Ver 2.000. Genetics and Biometry Lab, Dept. of Anthropology, University of Geneva, Switzerland. Selkoe KA and Toonen RJ (2006) Microsatellites for ecologists: A practical guide to using and evaluating microsatellite markers. Ecol Lett 9:615-629

Silva A (2011) Estrutura genética populacional de Prochilodus costatus Valenciennes 1850 (Characiformes, Prochilodontidae) no alto São Francisco. Masters Thesis, Universidade Federal de São Carlos. São Carlos, Brazil.

Sivasundar A, Eldredge B and Orti G (2001) Population structure and biogeography of migratory freshwater fishes (Prochilodus, Characiformes) in major South American rivers. Mol Ecol 10:407-417.

Slatkin, M (1985) Rare alleles as indicators of gene flow. Evolution 39:53-65. SlatkinM(1995) Hitchhiking and associative overdominance at a microsatellite locus. Mol Biol Evol 12:473- 480.

Valderrama M (1972) Operación subienda. Investigación pesquera. Instituto de Desarrollo de los Recursos Naturales Renovables, Bogotá, 94 pp. Valderrama-Barco MV and Petrere Jr M (1994) Crecimiento del bocachico Prochilodus magdalenae Steindachner 1878 (Prochilodontidae) y su relación con el régimen hidrológico en la parte baja de la cuenca del Río Magdalena (Colombia). Bol Cientif INPA 2:136-151.

Van Oosterhout C, Hutchinson WF, Wills DPM and Shipley P (2004) MICRO-CHECKER: Software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes 4:535-538.

Wasko AP and Galetti Jr P (2002) RAPD analysis in the Neotropical fish Brycon lundii: Genetic diversity and its implications for the conservation of the species. Hydrobiology 474:131- 137.

Weir B and Cockerham C (1984) Estimating F-statistics for the analysis of population structure. Evolution 38:1358-1370.

Wright E (1978) Evolution and Genetics of Populations. Vol. 2: The Theory of Gene Frequencies. University of Chicago Press, London. 519 pp.

Internet Resources

Earl D (2009) Structure Harvester v. 0.3, http://taylor0.biology.ucla.edu/structureHarvester/# (September 6, 2012).

Received: March 11, 2013;

Accepted: October 28, 2013.

Associate Editor: Fausto Foresti

All the content of the journal, except where otherwise noted, is licensed under a Creative Commons License CC BY-NC.

Publication Dates

• Publication in this collection
  24 Mar 2014
• Date of issue
  2014

History

• Accepted
  28 Oct 2013
• Received
  11 Mar 2013