Evidence of lack of population substructuring in the Brazilian freshwater fish Prochilodus costatus

FISH MOLECULAR GENETICS

SHORT COMMUNICATION

Evidence of lack of population substructuring in the Brazilian freshwater fish Prochilodus costatus

Luis F. Carvalho-Costa; Terumi Hatanaka; Pedro M. Galetti Jr.

Departamento de Genética e Evolução, Universidade Federal de São Carlos, São Carlos, SP, Brazil

^{Send correspondence to}

ABSTRACT

Curimbatá-pióa (Prochilodus costatus Valenciennes, 1850) is an endemic species from the São Francisco River Basin with migratory habits, important ecological role in the ecosystem and relevant fishing importance. The present study aimed to assess the genetic variation in P. costatus in order to investigate its population genetic structure. Genetic variation was studied at three sites downstream the Três Marias dam (Upper-middle São Francisco River, State of Minas Gerais) through six specific microsatellite loci. Fish from the three sites had quite similar genetic diversity levels and no genetic differentiation was detected, suggesting that P. costatus might represent a single reproductive unit in the studied area. Alternatively, the present study was not able to detect putative coexistent and comigrating populations along the main channel river. Our findings could be very helpful for the management and conservation of this fish.

Key words: Prochilodontidae, Curimbatá, microsatellites, genetic variation.

The São Francisco River is one of the main Brazilian rivers located in an area inhabited by approximately 13 million people (Cappio et al., 1995). The river was effectively fragmented into two parts since the Três Marias dam was built in the upper-middle region of the basin (Sato and Godinho, 2004). Although fish from 152 species inhabits the river basin, and some of them migrate in the reproduction season, no mechanism for fish passage exists at the Três Marias dam (Sato and Godinho, 2004).

One of these species is the endemic Prochilodus costatus (= P. affinis, Valenciennes, 1850), the "curimbatá-pióa", a long distance migratory species (Froese and Pauly, 2006) and an important resource for artisanal fishing communities, representing one of the commonly landed species of the region (Camargo and Petrere Jr, 2001). Moreover, this species can play an important ecological role considering its detritivorous habit that involves the processing of sediments (Flecker, 1996), representing a key-species for conservation programs. Little is known on population genetics of this fish. Previous molecular studies in some other migratory fish species inhabiting the upper-middle São Francisco River revealed the occurrence of population substructuring, and posed the question whether this is a common pattern in all migratory Neotropical freshwater fish (Wasko and Galetti Jr., 2002; Hatanaka and Galetti Jr., 2003; Hatanaka et al., 2006). Confronted with this scenario, the goals of this study were to assess the genetic diversity of Prochilodus costatus, and to investigate its population genetic structure downstream of the Três Marias dam.

The specimens were collected during the reproductive seasons of 2003/2004 downstream of the Três Marias dam near the municipality of Três Marias (State of Minas Gerais, southeastern Brazil). Fish were sampled at three sampling sites (Figure 1): A) the main channel of the São Francisco River immediately downstream of the dam; B) the Abaeté River; and C) the main channel of the São Francisco River downstream of the junction between the Abaeté and São Francisco Rivers (approximately 40 km downstream of the Três Marias dam).

DNA was extracted from liver tissues stored in 95% ethanol following the method of Aljanabi and Martinez (1997). Six previously isolated species-specific microsatellite loci were used to assess genetic variation following the protocol of Carvalho-Costa et al. (2006). The number of alleles per locus was assessed by genotyping the fragments in an automated sequencer using the protocol developed by Schuelke (2000). Allele sizes were scored against the internal size standard GeneScan-350-Rox (PE Applied Biosystems) through the Genescan 3.1 and Genotyper 2.5 softwares (ABI).

Allele frequencies, expected (He) and observed (Ho) heterozygosities, departures from Hardy-Weinberg expectations (HWE) and linkage disequilibrium were obtained in the program Genepop 3.3 (Raymond and Rousset, 1995). Population differentiation was estimated through Wright's fixation indices (Wright, 1978) according to the variance method of Weir and Cockerham (1984) implemented in the program Fstat 2.9.3.2 (Goudet, 2001). Allelic richness was also calculated in the program Fstat 2.9.3.2 and tested nonparametrically: Kruskal-Wallis and Mann-Whitney (for independent samples) and Wilcoxon (for related samples). These tests were also used to compare the observed heterozygosities among the sampling sites. All these tests were performed in Bioestat 3.0 (Ayres et al., 2003).

The results for intra-population genetic variation are shown in Table 1. We found allele numbers ranging from 43 to 48, including 28 private alleles. Site A presented 18.75%, site B 18.6% and site C 21.75% of private alleles. All sites presented departures from Hardy-Weinberg expectations for some loci (A: Pcos04 and Pcos18; B: Pcos03 and Pcos18; C: Pcos17 and Pcos18) after Bonferroni correction, also shown by positive F_IS values (Table 1). There was no significant non-random association among the genotypes of the six loci, suggesting that the analyses could be performed assuming statistical independence of the loci.

The observed and expected heterozygosities varied from 0.160 (Pcos20, site A) to 1.000 (Pcos14, site C) and from 0.150 (Pcos20, site C) to 0.950 (Pcos14, site C), respectively. Fish from the dam (A) were the ones that presented the largest mean Ho (0.480). Altogether, the average Ho was 0.45, a value comparable to that of other freshwater fish (Ho = 0.46, DeWoody and Avise, 2000). No significant Ho differences between sampling sites were observed by Kruskal-Wallis (H = 0.390, p = 0.823), Mann-Whitney (AxB: p = 0.689; BxC: p = 0.522; AxC: p = 0.749) and Wilcoxon tests (AxB: p = 0.249; BxC: p = 0.345; AxC: p = 0.753). Similarly, the allelic richness was not significantly different (Kruskal-Wallis: H = 0.246, p = 0.884; Mann-Whitney: AxB: p = 0.936, BxC: p = 0.631 and AxC: p = 0.522; Wilcoxon: AxB: p = 0.655, BxC: p = 0.5 and AxC: p = 0.5).

These results show similar levels of genetic variability for the three sites, differing from the results previously reported for the congeneric P. argenteus (Hatanaka and Galetti Jr., 2003; Hatanaka et al., 2006) and Brycon orthotaenia (Wasko and Galetti Jr., 2002). In both of the latter species, significant heterozygosity differences were observed when the downstream population nearest the dam (comparable to sample site A) was compared with the population further downstream (equal to site C), revealing population divergence. In all cases, the lower heterozygosity levels of the populations nearest the dam were correlated to the degraded environmental condition of this area (Wasko and Galetti Jr., 2002; Hatanaka and Galetti Jr., 2003).

The estimated pairwise F_ST values show no significant differences between the three sites (Dam x Abaeté River = -0.009, p = 0.6; Dam x Downstream junction of Abaeté/São Francisco Rivers = 0.004, p = 0.6; Abaeté River x Downstream junction of Abaeté/São Francisco Rivers = 0.006, p = 0.52). This pattern was also observed when loci Pcos17 and Pcos18 (loci showing departure from HWE) were eliminated from the analyses (data not presented).

Although sampling error could affect the accuracy of the estimated genetic population parameters, we raise two hypotheses to explain our results. The first one states that this species may constitute a single panmictic population in the Três Marias region, as was already reported for some fishes such as Prochilodus lineatus (Revaldaves et al., 1997; Sivasundar et al., 2001). This pattern may arise from the annual reproductive migration acting as a homogenizing agent of different gene pools. Combined with large effective population sizes (Turner et al., 2004) and strong gene flow, it could minimize the erosive effect of genetic drift on within-population genetic variation, thus homogenizing inter-population diversity.

On the other hand, we might not be able to detect a significant signal for population substructuring, considering previous evidence indicating that different genetic populations of other reophilic fishes may coexist and comigrate along the main river channel (Wasko and Galetti Jr., 2002; Hatanaka and Galetti Jr., 2003; Hatanaka et al., 2006). If this is true for P. costatus, hidden substructuring may explain the heterozygote deficit found in some loci (Wahlund effect).

With the use of microsatellites we were able to perform the first diagnosis on the genetic structure of this important endemic fish, which could be valuable for management and conservation purposes. Further analyses with a larger sample size, more sampling sites and temporal replicates are encouraged. Comparing the genetic diversity between fish from both up and downstream sites of the Três Marias dam would help to evaluate the impact of the dam on the genetic structure and the diversity of this fish.

Acknowledgments

We thank the Instituto Florestal de Minas Gerais for authorizing the capture of specimens. We also thank Yoshimi Sato and CEMIG/CODEVASF for collecting the samples, Flávio Henrique da Silva for enabling the use of the automated sequencer and Anna Carolina Dal Ri Barbosa for technical support. This research was supported by Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq).

Internet Resources

Received: August 24, 2006; Accepted: May 4, 2007.

Associate Editor: Maria Iracilda Sampaio

License information: This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

• Aljanabi SM and Martinez I (1997) Universal and rapid salt-extraction of high quality genomic DNA for PCR- based techniques. Nucleic Acids Res 25:4692-4693.
• Ayres M, Ayres Jr M, Ayres DL and Santos AS (2003) Bioestat: Versão 3.0. Sociedade Civil Mamirauá. MCT-CNPQ, Belém.
• Camargo SAF and Petrere Jr M (2001) Social and financial aspects of the artisanal fisheries of Middle São Francisco River, Minas Gerais, Brazil. Fish Manag Ecol 8:163-171.
• Cappio LF, Martins A and Kirchner R (1995) O Rio São Francisco. Editora Vozes, Petrópolis 110 pp.
• Carvalho-Costa LF, Hatanaka T and Galetti Jr PM (2006) Isolation and characterization of polymorphic microsatellite markers in the migratory freshwater fish Prochilodus costatus Mol Ecol Notes 6:818-819.
• DeWoody JA and Avise, JC (2000) Microsatellite variation in marine, freshwater and anadromous fishes compared with other animals. J Fish Biol 56:461-473.
• Flecker AS (1996) Ecosystem engineering by a dominant detritivore in a diverse tropical stream. Ecology 77:1845-1854.
• Hatanaka T and Galetti Jr PM (2003) RAPD markers indicate the occurrence of structured populations in a migratory freshwater fish species. Genet Mol Biol 26:19-25.
• 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.
• Raymond M and Rousset M (1995) Genepop (version 1.2): Population genetics software for exact tests and ecumenicism. J Heredity 86:248-249.
• Revaldaves E, Renesto E and Machado MFPS (1997) Genetic variability of Prochilodus lineatus (Characiformes, Prochilodontidae) in the Upper Paraná River. Braz J Genet 20:381-388.
• Sato Y and Godinho HP (2004) Migratory fishes of the São Francisco River. In: Carolsfeld J, Harvey B, Ross C and Baer A (eds) Migratory Fishes of South America: Biology, Fisheries and Conservation Status. IDRC, Victoria, 380 pp.
• Schuelke M (2000) An economic method for the fluorescent labeling of PCR fragments. Nature Biotech 18:233-234.
• Sivasundar A, Bermingham E and Ortí G (2001) Population structure and biogeography of migratory freshwater fishes (Prochilodus, Characiformes) in major South American rivers. Mol Ecol 10:407-417.
• Turner TF, McPhee MV, Campbell P and Winemiller KO (2004) Phylogeography and intraspecific genetic variation of prochilodontid fishes endemic to rivers of northern South America. J Fish Biol 64:186-201.
• Wasko AP and Galetti Jr PM (2002) RAPD analysis in the Neotropical fish Brycon lundii: Genetic diversity and its implications for the conservation of the species. Hydrobiol 474:131-137.
• Weir BS and Cockerham CC (1984) Estimating F-statistics for the analysis of population structure. Evolution 38:1358-1370.
• Wright S (1978) Evolution and the Genetics of Populations, v. 4: Variability Within and Among Natural Populations. The University of Chicago Press, Chicago, 590 pp.
• Goudet J (2001) Fstat, a program to estimate and test gene diversities and fixation indices, version 2.9.3. Available in http://www.unil.ch/izea/softwares/fstat.html
• Froese R and Pauly D (2006) Fishbase, versão 21. World wide web publication. Available in www.fishbase.org (access in March 2006).

Publication Dates

• Publication in this collection
  03 June 2008
• Date of issue
  2008

History

• Accepted
  04 May 2007
• Received
  24 Aug 2006