Services on Demand
On-line version ISSN 1678-4685
Genet. Mol. Biol. vol.30 no.3 São Paulo 2007
Maria Dolores Peres Lassala; Erasmo Renesto
Núcleo de Pesquisas em Limnologia, Ictiologia e Aqüicultura, Universidade Estadual de Maringá, Maringá, PR, Brazil
We estimated the genetic variability of nine fish species from the Brazilian upper Paraná River floodplain (Astyanax altiparanae, Hoplias malabaricus, Leporinus lacustris, Loricariichthys platymetopon, Parauchenipterus galeatus, Pimelodus maculatus, Rhaphiodon vulpinus, Roeboides paranensis and Serrasalmus marginatus) based on data for 36 putative allozyme loci obtained using corn starch gel electrophoresis of 13 enzymatic systems: aspartate aminotransferase (EC 184.108.40.206), acid phosphatase (EC 220.127.116.11), esterase (EC 18.104.22.168), glycerol-3-phosphate dehydrogenase (EC 22.214.171.124), glucose-6-phosphate dehydrogenase (EC 126.96.36.199), glucose-6-phosphate isomerase (EC 188.8.131.52), Iditol dehydrogenase (EC 184.108.40.206), isocitrate dehydrogenase NADP+ (EC 220.127.116.11), L-lactate dehydrogenase (EC 18.104.22.168), malate dehydrogenase (EC 22.214.171.124), malate dehydrogenase-NADP+ (EC 126.96.36.199), phosphoglucomutase (EC 188.8.131.52) and superoxide dismutase, (EC 184.108.40.206). The mean expected heterozygosity varied from zero to 0.147. When data from the literature for 75 species of tropical fish were added to the nine species of this study, the heterozygosity values differed significantly among the groups of different reproductive strategies. The highest mean heterozygosity was for the non-migratory without parental care, followed by the long-distance migratory, and the lowest mean was for the non-migratory with parental care or internal fecundation.
Key words: genetic variability, heterozygosity, Paraná River, polymorphism, reproductive strategies, tropical fish.
The reactions that species offer to environmental selective pressures is reflected in their physiological strategies, which can reduce maintenance energy, increase the efficiency of energy acquisition and, consequently, optimize reproductive efficiency. The success of such strategies may be estimated by the individual ability of being genetically present in the offspring (Agostinho and Júlio Jr., 1999).
In fish, reproductive strategies are extremely diversified, encompassing variations in the way they meet and attract partners, in spawning locations and parental care (Wootton, 1990). External fertilization and embryonic development are the dominant pattern for the majority of teleost fish, although internal fertilization and external embryonic development occurs in some Brazilian freshwater species such as Parauchenipterus galeatus. Blumer (1982) described parental care as any type of investment that enhances the survival probability of offspring (reproductive success) and which vary from pre-fertilization activities such as nest building (e.g. in Hoplias malabaricus) to mouth brooding and live-bearing.
Adaptation to heterogeneous environments depends upon the genetic variability of species, reduction in genetic variability resulting in reduced feeding and reproductive performance leading to inefficient exploration of the habitats (Kirpichnikov, 1992). Several authors have tried to explain the mechanisms which maintain genetic variability in natural populations. Nevo (1988) analyzed the relationship between heterozygosity and ecological and biological features such as climate (arctic, temperate and tropical), geographic range (wide, narrow, endemic), habitat type (underground, terrestrial, aquatic), habitat range (specialist, generalist) aridity (arid, sub-arid, sub-humid, mesic with a moderate or well-balanced moisture supply) and territoriality (territorial, non-territorial) in 1,111 species (vertebrates, invertebrates, and plants) and concluded that environmental heterogeneity is the major factor in maintaining and structuring genetic diversity in natural populations. Ward et al. (1994) compared the heterozygosity of marine, freshwater and anadromous (fish which live mostly in the sea but breed in fresh water) fish and concluded that marine fish have greater genetic variability than the other types. Heithaus and Laushman (1997) investigated the effects of ecology, life history and water quality on genetic variation of three stream-dwelling fish species and revealed that the genetic variability tends to decrease as species become more ecologically specialized.
Allozyme electrophoresis has been extensively used to estimate genetic variability in natural populations (Nevo, 1988; Van Der Bank et al., 1989; Ward et al., 1994) and our study employed this technique to estimate the genetic variability of nine species of fish from the Brazilian upper Paraná River floodplain and, jointly with literature data, verify the relationships between genetic variability and reproductive strategies.
Material and Methods
From March to December, 2002 in the Brazilian upper Paraná River floodplain, nine species from the most abundant fish with different reproductive strategies were sampled from Baía River (22°44' S, 53°17' W) and Paraná River (22°46' S, 53° 15'22" W) (Figure 1). The species sampled were Astyanax altiparanae Garutti & Britski 2000, Hoplias malabaricus (Bloch 1794), Leporinus lacustris Campos 1945, Loricariichthys platymetopon Isbrücker & Nijssen 1979, Parauchenipterus galeatus (Linnaeus 1766), Pimelodus maculatus Lacépède 1803, Rhaphiodon vulpinus Agassiz 1829, Roeboides paranensis Pignalberi 1975 and Serrasalmus marginatus Valenciennes 1847. A total of 267 specimens were sampled.
Immediately after capture, white skeletal muscle, gill, heart and liver tissue were removed from each specimen and frozen in liquid nitrogen. Tissues were homogenized with plastic sticks in 1.5 mL microcentrifuge tubes containing Tris/0.02 M HCl (1:1 w:v) buffer, pH 7.5. Carbon tetrachloride (CCl4) was added to the homogenized liver samples (1:2 v:v) due to the large amounts of fat present in the tissues (Pasteur et al., 1988). The homogenized samples were centrifuged at 45,114 x g for 30 min at temperatures between 1° and 5 °C and the supernatants submitted to horizontal electrophoresis in 15% corn starch gel (Val et al., 1981).
We evaluated 16 enzymatic systems (Table 1), enzyme nomenclature following the proposals of Murphy et al. (1996). Electrophoreses conditions were according to the following authors: Boyer et al., (1963) for SOD; Ruvolo-Takasusuki et al., (2002) for ACP and EST; Shaw and Prasad (1970) for AAT, G3PDH, G6PDH, GPI, IDDH, IDHP, LDH, MDH, MDHP and PGM. Standard histochemical staining procedures were used to visualize specific enzymes (Aebersold et al., 1987). Genetic interpretation of the gels was based on the quaternary structure of the enzymes (Ward et al., 1992). Data were analyzed using the POPGENE program version 1.31 (Yeh et al., 1997). Genetic variability was estimated using Nei's unbiased heterozygosity (He) or gene diversity (Nei, 1978). The observed (Ho) and expected (He) heterozygosities for each putative loci and the overall loci means were also calculated. Genotypic frequencies were tested for Hardy-Weinberg equilibrium using the chi squared (c2) test.
The species were organized in three distinct groups of reproductive strategy according to Vazzoler and Menezes (1992): group 1 (G1), containing the long-distance (> 100 km) migratory species P. maculatus and R. vulpinus; group 2 (G2), containing the sedentary or short-distance (< 100 km) migratory species A. altiparanae, L. lacustris, R. paranensis which provide no parental care; and group 3, containing the sedentary or short-distance migratory species H. malabaricus, L. platymetopon, P. galeatus, S. marginatus with internal fertilization or which provide parental care. Since it is unlikely that the species analyzed by us were representative of all tropical fish, we added data from the literature for a further 75 tropical fish species (61 from Brazil and 14 from Africa) for which the heterozygosity values had been estimated using allozyme data (Table 2), hence the total number of species analyzed was 84.
The significance of differences between groups was evaluated using a null models analysis of variance and the EcoSim 7 program (Gotelli & Entsminger, 2006).
The calculated genetic variability values for each species are presented in Table 3, from which it can be seen that the greatest frequency of polymorphic loci, average number of alleles per locus and heterozygosity were obtained for A. altiparanae, followed by H. malabaricus. No polymorphic locus was detected in R. vulpinus for the 29 loci analyzed. Except for P. maculatus the expected heterozygosity were higher than obtained heterozygosity, which indicates, on average, an excess of homozygotes for all loci.
All loci were in Hardy-Weinberg equilibrium (HWE) for L. lacustris, L. platymetopon, P. maculatus, R. vulpinus and R. paranensis. The loci not in HWE were G6pdh-1, Iddh-1 and Mdhp-1 for A. altiparanae, Gdh-1 for S. marginatus, Gdh-2 for P. galeatus, Adh-1, G6pdh-1, Gpi-A, Gpi-b, sMdh-B and Sod-1 for H. malabaricus.
Heterozygosity and reproductive strategies
The lowest heterozygosity values in our sample of nine species were for the G1 long-distance migratory species (R. vulpinus = 0 and P. maculatus = 0.011), which therefore also presented the lowest group average (He = 0.005). The highest average heterozygosity (He = 0.110) occurred in the G2 sedentary or short-distance migratory species without parental care, with A. altiparane having the highest heterozygosity (0.152) of all the species analyzed. The G3 sedentary or short-distance migratory species with internal fertilization or parental care exhibited intermediate heterozygosity (0.067).
When we considered the nine tropical fish species analyzed by us plus the 75 species from the other studies we found that the number of species in each reproductive strategy group and the average heterozygosity (He) ± the standard error for each group was as follows: G1 = 7 (Ho = 0.064 ± 0.021), G2 = 8 (He = 0.081 ± 0.016) and G3 = 69 (He = 0.046 ± 0.002) (Table 3, Figure 2), with He varying from zero to 0.142 for G1, from 0.006 to 0.152 for G2, and from zero to 0.143 for G3 (Tables 2 and 3). These data show that the highest average heterozygosity was presented by the G2 sedentary or short-distance migratory species without parental care, the G1 long-distance migratory species had and intermediate average heterozygosity and the G3 species with parental care or internal fertilization presented the lowest average heterozygosity.
A Shapiro-Wilk test for normality showed that the He values were not normally distributed (W = 0.9199, p < 0.0001) and the Levene test for homogeneity of variances showed a significant value (F = 4.843, p = 0.0103), indicating that the variances were heterogeneous. We conducted ANOVA but since the assumptions of normality and homogeneity of variances were not met by the data the significance of the ANOVA was tested by null models using the EcoSim 7 programs. The ANOVA detected a significant difference (F2, 81 = 4.02; p < 0.05) between the three reproductive strategy groups.
The genetic variability estimated in nine fish species from the upper Paraná River floodplain using Nei's gene diversity (He) varied from zero to 0.1518 with an average of 0.066, very near to the average of 0.051 described for 195 piscine species from several world-wide localities reviewed by Ward et al. (1992).
It has been generally accepted that populations of long-distance migratory species have higher genetic variability than non-migratory species because of their higher gene flow. Although this is true for several taxa (Ward et al., 1992), our data are not in agreement with this hypothesis. In our study, no polymorphic locus was found in a sample of 29 analyzed loci of R. vulpinus, a long-distance migratory species. This is surprising, since other cases of absence of genetic variability have been found in two sedentary species, Neoplecostomus paranensis from the upper Paraná River (Zawadzki et al., 2004b) and Hypostomus sp. 2 from the third-order Paraná River tributary Ribeirão Maringá (Paiva et al., 2005). Furthermore, P. maculatus, the other long-distance migratory species analyzed in our study, showed low genetic variability (He = 0.011), contrasting with the value of He = 0.132 estimated for other migratory species such as Prochilodus lineatus (Revaldaves et al., 1997) and Leporinus friderici (Chiari and Sodré, 1999).
The genetic variability estimated by us for S. marginatus (He = 0.045) and P. galeatus (He = 0.062) were intermediate compared to the average of 49 freshwater fish species (He = 0,046) estimated by Ward et al. (1994). On the other hand, our estimated values for A. altiparanae (He = 0,152), H. malabaricus (He = 0,137), R. paranensis (He = 0,084) and L. lacustris (He = 0,081) were high compared to the average found by Ward et al. (1994).
In the nine species analyzed by us the average heterozygosity of species with different reproductive strategies was not statistically different when compared amongst each other, which could lead to the conclusion that reproductive strategies are not related to genetic variability. However, our results could have been biased because the group of migratory species analyzed by us had only two species, so to compensate for this we carried out an analysis of our set of nine species combined with data on 75 tropical fish species from the literature. When the data for the 84 tropical fish species were analyzed, the relationship between heterozygosity and reproductive strategies were changed and there was a significant difference between the reproductive strategies groups (Figure 2). With this analysis, the average heterozygosity values which differed most were those between the G3 parental care group with and the G2 group without parental care. The heterogeneity of the heterozygosity variances between groups may have been due to the number of species analyzed, since the minimum and maximum heterozygosity values were similar among the groups.
The number of surviving offspring is likely to be greater in species with parental care than in species without parental care (Vazzoler, 1996), and species exhibiting parental care tend to be less variable. This can be observed in L. platymetopon, which brood their fry inside the male's mouth and is the most abundant species of the Paraná River floodplain. In species without parental care, few offspring are likely to survive, and a greater genetic variability is important to face environmental challenges.
Nei (1987) pointed out three main factors related to heterozygosity level are the quaternary structure of proteins, the molecular weight of the protein subunit and the species population size. In addition to these factors, we found that parental care is another factor that should be taken into account to maintain gene diversity. However, since a large number of factors influence the amount of genetic variability, we believe that further studies will be necessary to elucidate the relationship between genetic variability and ecological features of fish species.
The authors would like to thank Dr. Horácio Ferreira Júlio Jr. and Gilmar Lassala for helping with the capture and identification of fish, Jaime L. Pereira for drawing the map, Gislaine Iachstel Manetta for personal help and Dr. Luiz Carlos Gomes for statiscal help. We also thank NUPELIA-UEM for material support, and the Brazilian National Counsel of Technological and Scientific Development (Conselho Nacional de Pesquisa CNPq) and Coordenação de Aperfeiçoamento de Pessoal de Ensino Superior (CAPES) for financial support.
Aebersold PB, Winans GA, Tell DJ, Milner GB and Utter M (1987) Manual for starch gel electroforesis: A method for the detection of genetic variation. NOAA Technical Report NMFS 61:1-17. [ Links ]
Agostinho AA and Júlio HF Jr (1999) Peixes da bacia do alto rio Paraná. In: Lowe-McConnell RH (ed) Estudos Ecológicos de Comunidades de Peixes Tropicais. Editora da Universidade de São Paulo, São Paulo, pp 374-400. [ Links ]
Almeida FS and Sodré LMK (1998) Analysis of genetic variability in three species of Pimelodidae (Ostariophysi, Siluriformes). Genet Mol Biol 21:487-492. [ Links ]
Boyer SH, Fainer DC and Naughton MA (1963) Myoglobin: Inherited structural variation in man. Science 140:1228-1231. [ Links ]
Blumer LS (1982) A bibliography and categorization of bony fishes exhibiting parental care. Zool J Lin Soc 75:1-22. [ Links ]
Chiari L and Sodré LMK (1999) Genetic variability in five species of anostomidae (Ostariophysi, Characiformes). Genet Mol Biol 4:517-523. [ Links ]
Heithaus MR and Laushman RH (1997) Genetic variation and conservation of stream fishes: Influence of ecology, life history, and water quality. Can J Fish Aquat Sci 54:1822-1836. [ Links ]
Kirpichnikov VS (1992) Adaptative nature of intrapopulational biochemical polymorphism in fish. J Fish Biol 40:1-16. [ Links ]
Murphy RW, Sites JW, Buth DG Jr and Haufler CH (1996) Proteins: Isozyme electrophoresis. In: Hillis DM, Moritz C and Mable BK (eds) Molecular Systematics. 2nd edition. Sinnauer Associates, Sunderland, pp 51-120. [ Links ]
Nei M (1987) Molecular Evolutionary Genetics. Columbia University Press, New York, 512 pp. [ Links ]
Nei M (1978) Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics 89:583-590. [ Links ]
Nevo E (1988) Genetic diversity in nature. Evol Biol 23:217-246. [ Links ]
Paiva S, Renesto E and Zawadzki CH (2005) Genetic variability of Hypostomus (Teleostei, Loricariidae) from the Ribeirão Maringá, a stream of the Upper Rio Paraná basin, Brazil. Genet Mol Biol 28:370-375. [ Links ]
Pasteur N, Pasteur G, Bonhomme F, Catalan J and Britton-Davidian J (1988) Practical Isozyme Genetics. Ellis Horwood Limited, Chichester, 215 pp. [ Links ]
Renesto E, Zawadzki CH and Revaldaves E (2001) Biochemical taxonomy of Crenicichla (Pisces, Perciformes, Cichlidae) of the Iguaçu river, Brazil. Brazil Arch Biol Technol 44:15-22. [ Links ]
Renesto E, Zawadzki CH and Revaldaves E (2000) Genetic evidence for two species of genus Pimelodus Lacépède, 1803 (Siluriformes, Pimelodidae) in the Iguaçu River (Brazil). Genet Mol Biol 23:809-813. [ Links ]
Revaldaves E, Renesto E and Machado MFPS (1997) Genetic variability of Prochilodus lineatus (Characiformes, Prochilodontidae) in the Upper Paraná River. Rev Brazil Genet 20:381-388. [ Links ]
Ruvolo-Takasusuki MCC, Machado MFPS and Conte H (2002) Esterase-3 polymorphism in the sugarcane borer Diatraea saccharalis. Genet Mol Biol 25:61-64. [ Links ]
Shaw CR and Prasad R (1970) Starch gel electrophoresis of enzymes: A compilation of recipes. Biochem Genet 4:297-320. [ Links ]
Val AL, Schwantes AR, Schwantes MLB and Luca PH (1981) Amido hidrolisado de milho como suporte eletroforético. Ciência e Cultura 33:737-741 (Abstract in English). [ Links ]
Van Der Bank FH, Grobler JP and Du Preez HH (1992) A comparative biochemical genetic study of three populations of domesticated and wild African catfish, Clarias gariepinus. Comp Biochem Physiol 101:387-390. [ Links ]
Van Der Bank FH, Grant WS and Ferreira JT (1989) Electrophoretically detectable genetic data for fifteen southern African cichlids. J Fish Biol 34:465-483. [ Links ]
Vazzoler AEA and Menezes NA (1992) Síntese de conhecimentos sobre o comportamento reprodutivo dos Characiformes da América do Sul (Teleostei, Ostariophysi). Rev Bras Biol 52:627-640. [ Links ]
Vazzoler AEA (1996) Biologia da Reprodução de Peixes Teleósteos: Teoria e Prática. Editora da Universidade Estadual de Maringá, Maringá, 169 pp. [ Links ]
Ward RD, Skibinski DOF and Woodward M (1992) Protein heterozygosity, protein structure and taxonomic differentiation. Evol Biol 26:73-59. [ Links ]
Ward RD, Woodward M and Skibinski DOF (1994) A comparison of genetic diversity levels in marine, freshwater, and anadromous fishes. J Fish Biol 44:213-232. [ Links ]
Wootton RJ (1990) Ecology of Teleost Fishes. Chapman & Hall, London, 404 pp. [ Links ]
Yeh FC and Boyle TJB (1997) Population genetic analysis of co-dominant and dominant markers and quantitative traits. Belg J Bot 129:156-157. [ Links ]
Zawadzki CH, Renesto E and Bini LM (1999) Allozyme discrimination of three species of the genus Hypostomus Lacépède, 1803 (Osteichthyes, Loricariidae) from the rio Iguaçu basin (Brazil). Rev Suis Zool 106:91-105. [ Links ]
Zawadzki CH, Reis RE and Renesto E (2000) Allozyme discrimination of three species of Loricariichthys (Siluriformes, Loricariidae) from Southern Brazil. Rev Suis Zool 107:1-12. [ Links ]
Zawadzki CH, Weber C, Pavanelli CS and Renesto E (2002) Morphological and biochemical comparison of two allopatric populations of Hypostomus margaritifer (Regan, 1907) (Osteichthyes, Loricariidae) from the upper Paraná River basin, Brazil. Acta Scientiarum 24:499-505. [ Links ]
Zawadzki CH, Renesto E, Paiva S and Lara-Kamei MCS (2004a) Allozyme differentiation of four populations of Hypostomus (Teleostei, Loricariidae) from Ribeirão Keller, a small stream in the upper Rio Paraná basin, Brazil. Genetica 121:251-257. [ Links ]
Zawadzki CH, Alves AL, Renesto E and Oliveira C (2004b) Biochemical evidence of a possible new species of Neoplecostumus (Teleostei, Loricariidae) from the upper Rio Paraná basin, Brazil. Biochem Syst Ecol 32:573-582. [ Links ]
Zawadzki CH, Renesto E, Reis RE, Moura MO and Mateus RP (2005) Allozyme relationships in hypostomines (Teleostei, Loricariidae) from the Itaipu Reservoir, Upper Rio Paraná basin, Brazil. Genetica 123:271-283. [ Links ]
Gotelli NJ and Entsminger GL (2006) EcoSim: Null models software for ecology. Version 7. Acquired Intelligence Inc. & Kesey-Bear. Jericho, VT 05465. http://garyentsminger.com/ ecosim/index.html.
Send correspondence to:
Núcleo de Pesquisas em Limnologia, Ictiologia e Aqüicultura
Universidade Estadual de Maringá
Av. Colombo 5790
87020-900 Maringá, PR, Brazil
Received: July 19, 2006; Accepted: March 16, 2007.
Associate Editor: João S. Morgante