Identification of Astyanax altiparanae (Teleostei, Characidae) in the Iguaçu River, Brazil, based on mitochondrial DNA and RAPD markers

Identification of Astyanax altiparanae (Teleostei, Characidae) in the Iguaçu River, Brazil, based on mitochondrial DNA and RAPD markers

Sônia M.A.P. Prioli^1,2, Alberto J. Prioli^1,2, Horácio F. Júlio Jr.^1,2, Carla S. Pavanelli², Alessandra V. de Oliveira^2,3 , Helaine Carrer⁴, Dirce M. Carraro⁵ and Laudenir M. Prioli2,3

¹Departamento de Biologia Celular e Genética, Universidade Estadual de Maringá, Maringá, PR, Brasil.

2Núcleo de Pesquisas em Limnologia Ictiologia e Aqüicultura (NUPELIA), Universidade Estadual de Maringá, Maringá, PR, Brasil.

3Departamento de Biologia, Universidade Estadual de Maringá, Maringá, PR, Brasil.

4Departamento de Ciências Biológicas, Escola Superior de Agricultura Luiz de Queiroz, Universidade de São Paulo, Piracicaba, SP, Brasil.

5Ludwig Institute for Cancer Research, São Paulo, SP, Brasil.

ABSTRACT

Astyanax fishes are among the most important food-web components of South America rivers. In the Iguaçu River basin, the Astyanax genus is represented mainly by endemic species. For millions of years, that hydrographic basin has been geographically isolated from the Paraná River basin by the Iguaçu Falls. Recently, a species from the Upper Paraná River basin identified as Astyanax bimaculatus was revised and described as a new species named Astyanax altiparanae Garutti & Britski, 2000. Fauna endemism and geographic isolation triggered interest in investigations to evaluate the identification and genetic relatedness among two A. altiparanae populations from the Upper Paraná River basin and the population identified as A. bimaculatus in the Iguaçu River, upstream from the Iguaçu Falls. Mitochondrial DNA sequences and RAPD markers revealed high genetic diversity within each population, as well as low genetic distance, high gene flow, and high mitochondrial DNA similarity among all three populations. In conjunction with morphological similarities, these results demonstrated that the population presently known as Astyanax bimaculatus in the Iguaçu River should actually be stated as Astyanax altiparanae. Furthermore, it could be inferred that the A. altiparanae population is not endemic and most likely it was recently introduced in the Iguaçu River basin, maintaining the ancestral genetic identity.

Key words:Astyanax, fish, D-loop, control region, mitochondrial DNA, RAPD.

Received: July 31, 2002; accepted: November 4, 2002.

INTRODUCTION

Species of the genus Astyanax Baird & Girard, 1854 are among the most important components of the food-web in South America rivers, with a significant participation in the diet of big fishes. There are about one hundred nominal species and subspecies described in this genus, but many aspects of their taxonomy are still unknown. In the Iguaçu River basin, Paraná State, Brazil, Astyanax is represented mainly by seven or eight endemic species. The Iguaçu ichthyofauna has a peculiar evolutionary history, with a high rate of endemism in a cataract-ridden fluvial environment. Approximately 80% of the species are exclusive to this hydrographic basin and non-existent elsewhere on the planet (Agostinho and Gomes, 1997). Among the many falls in the river’s 1,060 km course, the Iguaçu Falls are the most outstanding. They are located about 22 km from the mouth of the Iguaçu River where it discharges into the Paraná River. Dated from the Oligocene, about 22 million years ago, and dropping 72 m, the Iguaçu Falls isolate the ichthyofauna of the Iguaçu River basin from that of the Paraná River basin (Haseman, 1911; Sampaio, 1988).

All Astyanax species existing in the Paraná River basin are absent in the hidrographic basin of Iguaçu, with exception of a fish population popularly known as tambiú and lambari-de-rabo-amarelo. In the Paraná River basin, this fish population had been identified as Astyanax bimaculatus, but recently it was revised and described by Garutti and Britski (2000) as a new species named Astyanax altiparanae Garutti & Britski, 2000. It is indicated that only the species from Surinam, apparently the holotype of this species, should a priori be called A. bimaculatus. The geographical isolation of the Iguaçu River basin and the special characteristics of its ichthyofauna endemism instigated an investigation as to whether the tambiú population from the Iguaçu River, above the falls, referred to as A. bimaculatus (Garavello et al., 1997) belongs to the A. altiparanae species found in tributaries of the Upper Paraná River basin. Current work in systematics of fishes has counted on the foundation and starting point provided by the remarkable success achieved by morphological studies. In the last decade, new suites of characters provided by molecular approaches have helped to invigorate studies of fish systematics. Interspecific and intraspecific genetic variation analyzed by mitochondrial DNA (mtDNA) sequences and Random Amplified Polymorphic DNA (RAPD) are among the molecular tools most widely used in these studies.

The mitochondrial genome co-evolves, with a certain degree of independence, at its own faster rate together with the host organism. Analyses of nucleotide variations in mtDNA sequences have been extremely helpful in polymorphism studies and phylogenetic analysis of natural fish populations (Kocher and White 1989; Stepien and Kocher 1997; Sivasundar et al., 2001). The mitochondrial control region is the main noncoding region and the most variable of metazoan mtDNA, including fishes. It is characterized by the displacement loop (D-loop), a DNA segment that is complementary to the L-strand, and that displaces the H-strand. In addition to conserved segments, the control region contains interposed sequences in the 5’ and 3’ domains, which have high levels of polymorphism. The rate of evolution of the control region is two to five times higher than that of the mitochondrial protein coding genes. Because mutations accumulate fastest in the control region, mainly in the 5’ domain, this is the main mtDNA sequence of choice for population and phylogenetic studies among closely related species, including fish species (Meyer, 1994; Faber and Stepien, 1977; Sivasundar et al., 2001).

The RAPD technique consists of DNA analysis of fragments amplified by Polymerase Chain Reaction (PCR) with single primers of randomized sequences (Williams et al., 1990). These primers, consisting of randomized sequences, do not discriminate between coding and non-coding regions. It may be thus admitted that the RAPD technique samples the genome in a highly randomized way (Williams et al., 1990). In fishes, RAPD has been successfully used to supplement systematic and phylogenetic studies of natural populations, including species and sub-species (Stepien and Kocher 1997; Almeida et al. 2001; Oliveira et al., 2002).

In the present work, mitochondrial DNA sequences and RAPD markers were used to study the genetic relatedness between the non-endemic Astyanax population from the Iguaçu River, upstream from the Iguaçu Falls, and two A. altiparanae populations from tributaries of the Paraná River basin, attempting to study the identification and possible origin of this Astyanax population existent in the Iguaçu River basin.

MATERIALS AND METHODS

Fish sampling and DNA extraction

Astyanax altiparanae¹ was sampled in the Keller and Pirapó tributaries of the Upper Paraná River basin, near Maringá, Northern Paraná State, Brazil. The Pirapó belongs to the Paranapanema River sub-basin and the Keller belongs to the Ivaí River sub-basin, both in the hydrographic basin of the Upper Paraná River. In the Iguaçu River, Astyanax² specimens of interest were collected at the Salto Caxias Hydroelectric reservoir, soon after impounding, located about 200 km from Iguaçu Falls. The related species, Astyanax sp B² was chosen as an outgroup for genetic analyses. This species is endemic to the Iguaçu River basin, and it was sampled at the Salto Caxias Hydroelectric reservoir. Specimens were captured by gill nets and immediately frozen in liquid nitrogen. Species identification was undertaken on the basis of morphological characteristics, in agreement with Sampaio (1988), Garavello et al. (1997), and Garutti and Britski (2000). Samples of muscle tissues were taken for DNA extraction. Sampled fishes were numbered, preserved in 10% formaldehyde for a month, and then they were transferred to 70 °GL ethanol.

Extraction of total DNA was performed according to Whitmore et al. (1992), with a few modifications. Approximately 1 g of muscle tissue was macerated in the presence of liquid nitrogen. Extraction buffer (500 mL) was added, consisting of 100 mM Tris-HCl pH 8.0, 10 mM EDTA (ethylenediaminetetraacetic acid), 0.1% SDS (sodium dodecyl sulphate), 50 mM DTT (dithiothreitol), and 0.5 mg/mL proteinase K. After shaking for one hour at 60 Srpm, at 42 °C, DNA was extracted with an equal volume of phenol, and then with chloroform/isoamylic alcohol (24:1 v:v). DNA was then precipitated by addition of two volumes of cold ethanol 100%, followed by incubation at -20 °C for 12-24 h. After centrifugation at 10,000 x g for 10 min, the precipitate was washed with ethanol 80% and resuspended in diluted TE buffer (0.1 mM Tris-HCl pH 8.0, 0.01 mM EDTA), containing RNAase (20 mg/mL). The DNA suspension was stored at -21 °C. Aliquots of each sample were used for visual estimation of DNA quantity by comparison with known quantities of phage l DNA in agarose gel (0.8%), and stained with ethydium bromide (20 mg/100 mL).

PCR and sequencing of mitochondrial DNA

A mtDNA fragment (approx. 800 bp) was amplified by PCR from DNA samples of six to seven fishes of each Astyanax population. For each fish specimen, the mtDNA fragment was amplified in three independent PCR reactions and they were then sequenced and analyzed, as replicates. The primers H16498 5’-CACCTGAAGTAGGAACCA GATG-3’ (Meyer et al., 1990), and L15774M 5’-ACATGAATTGGAGGTATACCAGT-3’, which was modified from the L15774 primer described by Irwin et al. (1991), were used for PCR amplifications. PCR reaction mix consisted of Tris-KCl (20 mM Tris-HCl pH 8.4 plus 50 mM KCl), 1.5 mM MgCl₂, 2.5 mM of each primer, 0.1 mM of each dNTP, 2.5 U Taq-polymerase, 15 ng total DNA and water to a total volume of 25 mL. PCR amplification was performed in a first cycle of 94 °C for 4 min, 50 °C for 30 s, and 72 °C for 2 min, followed by 40 cycles of 94 °C for 15 s, 56 °C for 30 s, 72 °C for 2 min, and a final extension step of 72 °C for 10 min. After amplification, an aliquot was taken from each PCR reaction mix, and the amplified DNA fragment was quantified in agarose gel by comparison with known quantities of phage l DNA.

The final DNA product obtained from each PCR reaction was used directly as template for sequencing on an ABI-3100 (Perkin Elmer, Norwalk, CT, USA) automated sequencer. Sequencing was performed with BigDye^TM Terminator Cycle Sequencing v. 2.0 kit (Perkin Elmer Biosystems), according to the manufacturers instructions. Approximately 50 ng of template DNA and 20 pmol of either L15774M or H16498 primer were added to each sequencing reaction. Immediately after heating the reaction mix at 94 °C for 4 min, amplifications were performed in 35 cycles of 30 s at 94 °C, 30 s at 55 °C, and 1 min 30 s at 60 °C, followed by 5 min at 60 °C. Sequencing data, under chromatogram format, were directly collected by the Sequencing Analysis program (Perkin Elmer) and submitted to quality check, assembly and alignment on the Vector NTI Suite version 6.0 software (Informax, Inc.). The nucleotide sequence of each specimen was confirmed in at least three replicates, consisting of three independent PCR reactions.

A 569-nucleotide portion of the 5’ side of the mitochondrial control region of 24 individuals was used in genetic relatedness analyses. In addition, a corresponding mtDNA segment of the control region of the closely related species Astyanax sp B was sequenced and included as an outgroup in the analysis. Genetic analyses were performed with MEGA version 2.1 (Kumar et al., 2001). Genetic distances of Tamura and Nei (1993) with gamma distribution were calculated and clustered into a distance neighbor-joining tree (Saitou and Nei, 1987). Bootstrap analyses were based on 1000 replications. The Gamma coefficient was estimated by using PAUP* version 4.0 (Swofford, 2002).

RAPD electrophoresis profile analysis

Eighteen fishes from the Iguaçu River, and eighteen from both tributaries of the Upper Paraná River basin were evaluated by RAPD markers. Initial tests of PCR primers were undertaken by using Astyanax DNA in amplifications with two sets of 20 decamers of Kits X and W from Operon Technologies (Alameda, CA, USA). Primers amplifying the highest number of fragments and producing accurate bands were chosen for RAPD analysis. PCR reaction mixture consisted of buffer Tris-KCl (20 mM Tris-HCl pH 8.4 and 50 mM KCl), 2 mM MgCl₂, 0.46 mM primer, 0.2 mM dNTP, 1U/Reaction of Taq-polymerase (Gibco BRL), 10 ng DNA, and sterile deionized water to a total volume of 13 mL. The reaction mixture was heated at 92 °C for 4 min, followed by 40 cycles of 1 min at 92 °C, 1 min 30 s at 40 °C and 2 min at 72 °C. Immediately after the last amplification cycle, the reaction mixture was maintained at 72 °C for 5 min and then cooled at 4 °C for 20 min. PCR amplified DNA fragments were separated by electrophoresis in agarose gel (1.4%), stained with ethydium bromide (20 mg/100 mL), at 3 V.cm^-1 for 5-6 h. A sample was included as a negative control in each experiment, consisting of the PCR reaction mixture without template DNA. In addition, 1 to 3 samples of PCR products, containing Astyanax DNA fragments previously amplified and analyzed with the same primer, were included in each agarose gel. The electrophoretic profile was visualized under UV radiation and photographed with Kodak EDAS-290. Sizes of DNA fragments were estimated by comparison with the standard Ladder 100 pb (Gibco BRL).

Unbiased genetic distances based on RAPD markers were estimated according to Nei (1978), with Lynch and Milligan’s (1994) correction, using the software Rapddist (Black and Antolin, 1997). The Shannon gene diversity indexes were calculated with the software Popgen 1.31 (Yeh et al., 1999). The Shannon diversity index was used to quantify levels of genetic diversity in isolated populations and for the partition of components within and between populations. Shannon’s diversity index has been useful for RAPD data analysis due to its relative non-sensitivity to bias, which could be introduced in data because of the impossibility of detecting heterozygous individuals. Fixation index F_st, and gene flow, equivalent to the number of migrants per generation (Nm), were estimated according to Lynch and Milligan (1994), with the computer program RAPDFST (Black, 1997). The arithmetic complement of the coefficient of similarity among individuals, according to Nei and Li (1979) was estimated by using RAPDPLOT 3.0 (Black, 1997). Clustering was performed by the neighbor-joining method (Saitou and Nei, 1987), with MEGA 2.1 (Kumar et al., 2001).

RESULTS

Mitochondrial DNA

A mtDNA fragment of approximately 800 bp was amplified by PCR from all Astyanax specimens. This mtDNA fragment consisted of approximately 90 bp of the 3’ end of the mitochondrial cytochrome b gene (GenBank Accession Numbers AF531481 to AF531504), immediately followed by the complete nucleotide sequences of tRNA^Thr and tRNA^Pro genes (GenBank Accession Numbers AF531874 to AF531897) and by a partial sequence of the 5’ side of the control region (GenBank Accession Numbers AY125820 to AY125843). Nucleotide sequences of the control region and cytochrome b were identified by similarity to GenBank accesses via BLASTN (Altschul et al., 1997) from NCBI. The tRNA genes were identified with the tRNAScan-SE program (Lowe and Eddy, 1997). Altogether, twenty mtDNA sequences, amplified from twenty specimens, were analyzed. The cytochrome b, tRNA^Thr, and tRNA^Pro nucleotide sequences were identical among individuals of all three Astyanax populations, thus these genes were excluded from analyses. Nucleotide polymorphism was restricted to the amplified control region sequence, from which 569 bp were analyzed in this work. As shown in Figure 1, multiple alignment of these Astyanax sequences revealed 27 polymorphic nucleotides, approximately 4.7%, in the Keller, Pirapó, and Iguaçu populations. Altogether, a total of eleven different control-region haplotypes were identified among twenty individuals of these populations. Changes were mostly single nucleotide substitutions. On the other hand, in comparisons of these three populations with the outgroup Astyanax sp B, the percentage of polymorphic sites was increased to approximately 30%. A high polymorphism level was also detected when the nucleotide sequence of the mitochondrial control region from another characid species, Oligosarcus paranensis, was tested as the outgroup, replacing Astyanax sp B (data not shown).

As shown in Table I, low genetic distances of Tamura and Nei (1993) with gamma distribution, based on control region nucleotide sequences, were estimated among the Keller, Pirapó, and Iguaçu populations, at a maximum value of 0.019. On the other hand, high genetic distances, ranging from 3.176 to 3.527, were calculated between these three populations and Astyanax sp B. As expected, in the neighbor-joining tree, the Keller, Pirapó, and Iguaçu specimens were combined as a single cluster, distant from the outgroup (Figure 2). Within the subtree, an intermix of individuals from the Keller, Pirapó and Iguaçu rivers clearly demonstrated the high molecular similarity among these Astyanax specimens. The distances among individuals of the three populations were in a similar magnitude observed among the individuals from the outgroup population, as represented in the tree branches (Figure 2). Low bootstrap values can be attributed to the low molecular differences among individuals within the cluster. As a contrast, this subtree was clearly separated in a high distance, with 100% bootstrap value, from the related Astyanax sp B species used as outgroup.

RAPD markers

In each Astyanax population, a total of 10 selected primers were used for DNA amplifications by PCR. Amplified DNA fragments counted for 79 polymorphic bands, seven 100% monomorphic bands and one monomorphic band at 5% probability level. RAPD electrophoresis profiles are represented in Figure 3, for one of the selected primers. Sizes of polymorphic bands ranged from approximately 0.4 to 3.0 kb (Table II). The number of bands per primer varied from five to ten. Polymorphic bands were taken as RAPD markers, and they were combined in a simple matrix of binary data, with registers of presence or absence of bands. The Shannon genetic diversity index within each population was 0.58 (±0.15) for Iguaçu, 0.54 (±0.20), for Keller, and 0.50 (±0.25) for Pirapó. Estimates of Nei’s unbiased genetic distances ranged from 0.050 to 0.070 (Table III). As shown in Table IV, the estimated value of fixation index F_st, was not significant between the Iguaçu and Keller populations, and it was significant, only at the 5% level, between the Pirapó and Keller, as well as Pirapó and Iguaçu. The migration model among the three populations was adopted to estimate gene flow (Nm) from F_st (Lynch and Milligan, 1994). Each population was considered sufficiently large to present negligible genetic drift. Allele frequency of migrants among the three populations was considered to be equal to the mean of the allele frequency in the entire population. Estimates of gene flow on all loci among populations varied from 2.4 to 4.0 migrants per generation (Table IV). The low RAPD divergence among populations, as derived from the neighbor-joining tree (Figure 4), indicated that the three studied Astyanax populations form a homogeneous group; therefore, a result congruent with the high genetic similarity estimated by all analyses performed for both mtDNA and RAPD markers.

DISCUSSION

As inferred by the mitochondrial DNA sequences and RAPD markers, the genetic similarity among the three studied Astyanax populations is so high that it is not consistent with their current taxonomic identification as two different species. Genetic relatedness estimated by molecular markers indicated that the population from the Iguaçu River and the Astyanax altiparanae population from the Keller tributary of theUpper Paraná River basin are close to genetic homogeneity and they did not differentiate, not even by genetic drift. Pirapó was the most differentiated population, as compared either to Keller or to Iguaçu. The molecular genetic evidence, in conjunction with morphological similarities (Graça and Pavanelli, 2002), demonstrated that the species status of Astyanax bimaculatus from the Iguaçu River should actually be stated as Astyanax altiparanae Garutti & Britski, 2000. Furthermore, high genetic diversity within and high similarities between the populations, as indicated by of RAPD markers, suggested that this species was most probably introduced recently in the Iguaçu River basin, and that this introduction involved a large number of individuals.

Molecular tools associated to morphological analyses have been effectively applied in defining relationships of fishes in different levels, including among genus and species populations. The mitochondrial control region has been used extensively as a molecular marker because of its rapid rate of evolution compared to nuclear DNA and the significant literature available on the theory and application of mtDNA sequence analysis. In the present study, low genetic distances, according to Tamura and Nei (1993), were found among the Iguaçu, Keller, and Pirapó populations, as inferred from control region sequences (Table I; Figure 2). The high genetic similarity data estimated on the basis of the control region were reinforced by the cytochrome b partial sequence and the complete sequences of tRNA^Pro and tRNA^Thr, which were identical among all Astyanax individuals, from both populations of the Upper Paraná basin and the population from the Iguaçu River. These results are consistent with findings in other vertebrate species, showing that although the mitochondrial protein encoding gene and the tRNAs do not evolve as rapidly as the control region, their nucleotide sequences are highly conserved at the intraspecific level, while they have been found polymorphic at the interspecific level (Meyer, 1994; Saccone et al., 2002). Furthermore, the genetic relatedness of all three Astyanax populations was obvious upon the contrast of their mtDNA sequences with those of the related species Astyanax sp B. As shown in Figure 2, individuals from the Iguaçu and Upper Paraná rivers were not discriminated by the neighbor-joining method, but they were clearly separated, with 100% bootstrap value, from Astyanax sp B. The genetic distances among specimens of all three studied Astyanax populations were in similar magnitudes to those observed among the individuals of the outgroup Astyanax sp B (Figure 1). The molecular genetic relatedness inferred from mtDNA sequences indicated that the studied Astyanax population sampled in the Iguaçu River should be joined together with the Keller and Pirapó populations, within the Astyanax altiparanae species.

RAPD markers corroborated the inference of high genetic similarity among the three Astyanax populations, as estimated from mtDNA sequences, and they indicated high nuclear genetic diversity within each population. The Shannon indexes based on RAPD markers revealed high genetic diversity within the Iguaçu population, at very similar levels to those of the Keller and Pirapó populations. This result indicated that the genetic variability was not lost during the establishment of these three populations. The highest Nei’s unbiased genetic distance, estimated as 0.070 from RAPD markers (Table III), is considered low and consistent with intraspecific distances. Nei’s genetic distances between species from the same genus have been found no lower than 0.162 and over 0.3 in most species, as estimated from enzyme studies (Thorpe and Solé-Cava, 1994). Analyses of RAPD markers of fish species from the floodplain of the Upper Paraná River basin revealed relatively high Nei’s unbiased genetic distances, as 0.504 between Pseudoplatystoma corruscans (‘pintado’) and Pseudoplatystoma fasciatum (‘cachara’), and 0.729 between Steindachnerina brevipinna and Steindachnerina insculpta (Sekine et al., 2000; Oliveira et al., 2002). Even though the studied species within each of these two genus are so closely related that they crossbreed generating hybrids, their genetic distances were obviously much higher than those among the three Astyanax populations. Therefore, the low genetic distances among the two Astyanax altiparanae populations from the Upper Paraná river basin and the studied Astyanax population from Iguaçu river are coherent with expected intraspecific distances. Nei’s unbiased genetic distances are considered suitable for long evolutionary processes, with populational divergences due to genetic drift and mutational events (Weir, 1990), and its value is proportional to the time of divergence and the rate of gene substitution per locus and per generation. Nei’s unbiased estimates of low genetic distances based on RAPD did not discriminate the Astyanax populations from the Iguaçu and Upper Paraná river basins, indicating that they maintain the genetic identity of a single recent common ancestor. Intraspecific differentiation was significant only at 5% level for the Pirapó population as compared to both the Iguaçu and Keller populations, by the fixation index F_ST. This result seems to indicate an expected level of intraspecific variation. RAPD estimates of gene flow on all loci among all populations was high, ranging from 2.4 to 4.0 migrants per generation (Table IV), indicating that the all three populations are very similar, and did not differentiate, not even by genetic drift.

The Iguaçu Falls have been considered the principal factor for fish isolation and speciation in the Iguaçu River basin, as initially suggested by Haseman (1911). In fact, most fish species from the Iguaçu River basin have been characterized as endemic, including species of the genus Astyanax (Sampaio, 1988; Severi and Cordeiro, 1994; Garavello et al., 1997). Nevertheless, introduced species, admittedly not very abundant, have also been reported in this hydrographic basin (Agostinho and Gomes, 1997). The high genetic similarity and high gene flow estimated among A. altiparanae populations from tributaries of the Upper Paraná River basin and from the Iguaçu population, as well as the high genetic diversity within each population, at very similar levels in the two hydrographic basins, demonstrated that the Astyanax altiparanae population from the Iguaçu River is not endemic. Such high genetic similarity agrees with the hypothesis of introduction of this species in the Iguaçu River from the Upper Paraná River basin (Garavello et al., 1997). Furthermore, the molecular data in the present work suggest a recent separation of the Iguaçu and Upper Paraná populations, maintaning the intrapopulational genetic diversity representative of the ancestral population. In addition, the time or number of generations would not be yet sufficient for mutation accumulation of mutations, as detected by the methodology used. Another possible hypothesis to explain these results, apparently less likely, would be an old introduction of A. altiparanae from the Upper Paraná River basin into the Iguaçu River, followed by recurrent migration or transposition of a high number of individuals, to such a degree that the genetic differentiation would be mitigated or eliminated.

ACKNOWLEDGMENTS

The authors gratefully acknowledge Wilson M. Furuya for helping with fish sampling, Weferson J. Graça for assisting with morphological comparisons of studied populations, Reimar Padovani for helping in DNA sequencing, and COPEL and NUPELIA-UEM for logistic support. Part of this research was supported by grants from CNPq, CAPES, and COPEL.

Oliveira AV, Prioli AJ, Prioli SMAP, Pavanelli CS, Julio Jr HF and Panarari RS (2002) Diversity and genetic distance in populations of Steindachnerina in the Upper Paraná river floodplain. Genetica (in press).

Send correspondence to A.J.P. E-mail: [e-mail].

² Voucher Specimens from Iguaçu River basin: Astyanax altiparanae - NUP1623 (100ex.), and Astyanax sp B – NUP2457 (44 ex.), Iguaçu River, Caxias Reservoir, Capitão Leônidas Marques and Nova Prata do Iguaçu townships, 25° 30’ S/53° 25’ W, PR, Brazil, collected by Nupelia staff, January 25 ^th, 1999.

References

• Agostinho AA and Gomes LC (1997) Manejo e monitoramento de recursos pesqueiros: perspectivas para o reservatório de Segredo In: Agostinho AA and Gomes LC (eds) Reservatório de Segredo: Bases Ecológicas para o Manejo. Eduem, Universidade Estadual de Maringá, PR, Brazil, pp 319-364.
• Almeida FS, Pelegrinelli MH and Sodré LMK (2001) RAPD and isoenzyme analysis of genetic variability in three allied species of catfish (Siluriformes: Pimelodidae) from the Tibagi River, Brazil. J Zool Lond 253:113-120.
• Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W and Lipman DJ (1997) Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 25:3389-3402.
• Black, B (1997) RAPDFST 401 - A FORTRAN Program to estimate F(ST) and effective migration rates among subpopulations using RAPD-PCR files. Colorado State University.
• Black B and Antolin M (1997) Explanation of RAPDDIST 10. Colorado State University.
• Faber JE and Stepien CA (1977) The utility of mitochondrial DNA control region sequences for analyzing phylogenetic relationships among populatios, species, and genera of the Percidae In: Kocher TD and Stepien (eds) Molecular Systematics of Fishes. Academic Press, New York, pp 129-143.
• Garavello JC, Pavanelli CS and Suzuki HI (1997) Caracterização da ictiofauna do rio Iguaçu In: Agostinho AA and Gomes, LC (eds) Reservatório de Segredo: Bases Ecológicas para o Manejo. Eduem Universidade Estadual de Maringá, pp 61-84.
• Garutti V and Britski HA (2000) Descrição de uma espécie nova de Astyanax (Teleostei: Characidae) da bacia do alto Rio Paraná e considerações gerais sobre as demais espécies do gênero na bacia. Comum Mus Ciênc Tecnol PUCRS Sér Zool 13:65-88.
• Graça WJ and Pavanelli CS (2002) Comments on the geographic distribution of Astyanax altiparanae Garutti & Britski, 2000 (Pisces, Characiformes) in the Rio Iguaçu basin. Acta Scientarium 24(2):(In Press).
• Haseman, JD (1911) Some new species of fishes from the rio Iguaçu. Ann Carnegie Mus 7:374-387.
• Irwin DM, Kocher TD and Wilson AC (1991) Evolution of cytochrome b of mammals. J Mol Evol 32:128-144.
• Kocher TD and White TJ (1989) Evolutionary analysis via PCR. In: Erlich HA (ed) PCR Technology: Principles and Applications for DNA Amplification. Stockton Press. New York. pp 137-147.
• Kumar S, Tamura K, Jakobsen IB and Nei M (2001) MEGA2: Molecular Evolutionary Genetics Analysis software. Bioinformatics. Arizona State University, Tempe, Arizona, USA.
• Lowe TM and Eddy SR (1997) tRNAscan-SE: A program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 25:995-964.
• Lynch M and Milligan BG (1994) Analysis of population structure with RAPD markers. Mol Ecol 3:91-99.
• Meyer A (1994) DNA technology and phylogeny of fish In: Beaumont AR (ed) Genetics and Evolution of Aquatic Organisms. Chapman & Hall, London, pp 219-249.
• Meyer A, Kocher TD, Basasibwaki P et al. (1990) Monophyletic origin of Victoria cichlid fish suggested by mitochondrial DNA sequences. Nature 347:550-553.
• Nei M (1978) Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics 89:583-590.
• Nei M and Li W-H (1979) Mathematical model for studying genetic variation in terms of restriction endonucleases. Proc Natl Acad Sci USA 76:5269-5273.
• Saccone C, Gissi C, Reyes A, Larizza A, Sbisà E and Pesole G (2002) Mitochondrial DNA in Metazoa: degree of freedom in a frozen event. Gene 286:3-12.
• Saitou N and Nei M (1987) The Neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol Biol Evol 4:406-425.
• Sampaio FAA (1988) Estudos taxonômicos preliminares dos ChaNeighbor-joining tree, based on the matrix of the arithmetic complement of the similarity coefficient of Nei and Li (1979), obtained from RAPD markers of Astyanax altiparanae individuals from the Iguaçu River, and from the Keller and Pirapó tributaries of the Upper Paraná River basin.raciformes (Teleostei, Ostariophysi) da bacia do rio Iguaçu, com comentários sobre o endemismo dessa fauna. MSc Thesis, Universidade Federal de São Carlos, São Carlos.
• Sekine ES, Prioli AJ, Julio Jr HF, Prioli SMAP and Castro LL (2000) Análise da divergência genética entre populações de Pseudoplatystoma fasciatum e Pseudoplatystoma corruscans (Pisces-Siluriformes) utilizando marcadores RAPD. Genet Mol Biol (Suppl) 23(3):101
• Severi W and Cordeiro AAM (1994) Catálogo de Peixes da Bacia do Rio Iguaçu. IAP/GTZ, Curitiba, 112 pp.
• Sivasundar A, Bermingham E and Orti G (2001) Population structure and biogeography of migratory freshwater fishes (Prochilodus: Characiformes) in major South American rivers. Mol Ecol 10(2):407-417.
• Stepien CA and Kocher TD (1997) Molecules and morphology in studies of fish evolution In: Kocher TD and Stepien (eds) Molecular Systematics of Fishes. Academic Press, New York, pp 1-11.
• Swofford DL (2002) PAUP*: Phylogenetic Analysis using Parsimony (*and other methods), vers 4. Software, Sinauer Associates, Sunderland.
• Tamura K and Nei M (1993) Estimation of the number of nucleotide substitutions in the control region of mitochondrial DNA in humans and chimpanzees. Mol Biol Evol 10:512-526.
• Thorpe JP and Solé-Cava AM (1994) The use of allozyme electrophoresis in invertebrate systematics. Zool Scr 23:3-18.
• Weir BS (1990) Genetic Data Analysis. Sinauer Associates Inc., Sunderland, 376 pp.
• Whitmore DH, Thai TH and Craft CM (1992) Gene amplification permits minimally invasive analysis of fish mitochondrial DNA. Trans Am Fish Soc 121:170-177.
• Williams JGK, Kubelik AR, Livak KJ, Rafalski JA and Tingey SV (1990) DNA polymorphisms amplified by arbitrary primers are useful as genetic markers. Nucl Acids Res 18:6531-6535.
• Yeh FC, Boyle T, Ye Z and Xiyan JM (1999) POPGENE Version 131: Microsoft Window-based freeware for population genetic analysis. University of Alberta and Center for International Forestry Research.

¹

Voucher Specimens from Upper Paraná River basin: Astyanax altiparanae- NUP2700 (20 ex.), Keller River, tributary to Ivaí River, Marialva, Paraná State (PR), Brazil, collected by Prioli AJ & Júlio Jr HF, October, 2000, and NUP2701 (12 ex.), Cidade Canção fisher-ground, Águas Maringá stream, affluent to Pirapó River, Maringá, PR Brazil, collected by Furuya WM, Prioli JA & Prioli SMAP, October 10

^th, 2000.

Publication Dates

• Publication in this collection
  12 Mar 2003
• Date of issue
  2002

History

• Received
  31 July 2002
• Accepted
  04 Nov 2002