Acessibilidade / Reportar erro

Evidence of polygamy in the socially monogamous Amazonian fish Arapaima gigas (Schinz, 1822) (Osteoglossiformes, Arapaimidae)

Abstracts

Arapaima gigas is one of the largest freshwater fishes of the world. It is socially monogamous, forming pairs, constructing a nest and providing parental care. We performed a paternity analysis under three scenarios in captive, semi-natural and natural areas using 10 microsatellite markers. As a positive control, we analyzed three pairs and their offspring isolated individually in artificial breeding ponds (a priori very high probability of monogamy). We then analyzed two samples of offspring from large artificial ponds with multiple adults but only one reproductive pair (a priori high probability of monogamy), two samples from semi-natural breeding station with multiple adults but only one reproductive pair (a priori high probability of monogamy), and a sample from a natural lake with multiple adults, some potentially breeding (a priori medium probability of monogamy). Analysis of patterns of Mendelian heredity suggested an extra-pair contribution for all broods except the positive controls. Similarly, results based on multilocus analysis estimated at least two sib-groups per nest. These results reject monogamy as a system of breeding in Arapaima gigas. From a management perspective, this behavior may be exploited to maintain genetic diversity in captive and as well in wild populations of Arapaima gigas.

Mate choice; Microsatellites; Multiple paternity; Pirarucu; Polygamy


O pirarucu Arapaima gigas é um dos maiores peixes de água doce do mundo. É socialmente monogâmico, forma casais, constrói ninhos e fornece cuidado parental. Com o objetivo de acessar o sistema de acasalamento do pirarucu, analisamos três cenários: em áreas de cativeiro, semi-naturais e naturais, utilizando 10 marcadores microssatélites. Como controle positivo, analisamos três casais e suas ninhadas isolados em açudes individuais (probabilidade a priori muito alta de monogamia). A seguir, analisamos duas amostras de ninhadas de um açude com vários adultos, mas somente um casal reprodutivo (probabilidade a priori alta de monogamia), duas amostras de estação de criação semi-natural com vários adultos mas somente um casal reprodutivo (probabilidade a priori alta de monogamia), e uma amostra de lago natural com vários adultos alguns potencialmente em fase de reprodução (probabilidade a priori média de monogamia). Análises de padrões mendelianos de hereditariedade sugerem contribuição extra-par para todas as ninhadas, exceto as do controle positivo. Similarmente, resultados baseados em análises multilocus realizadas no programa KINALYZER estimaram pelo menos dois grupos-irmãos por ninhada. Nossos resultados rejeitam a monogamia como sistema de acasalamento em Arapaima gigas. Da perspectiva de manejo, esse comportamento pode ser explorado para manter a diversidade genética em cativeiro assim como em populações naturais de Arapaima gigas.


Introduction

Arapaima gigas (the 'pirarucu' or 'arapaima') is one of the largest freshwater fishes of the world reaching up to three meters, and weighing up to 200 kg. It is also one of the economically most important fish species in the Amazon basin, forming an integral part of the Amazon people's traditional diet, constituting an important food source, and thus making it subject to intense commercial and subsistence exploitation. Indeed, this species became scarce in the 70's and 80's, and in some cases even extinct near large Amazon cities (Goulding, 1980Goulding, M. 1980. Fishes and the Forest: Explorations in Amazonian Natural History. Los Angeles, CA, University of California Press.). In 1975, the 'arapaima' was listed in Appendix II of CITES (Convention on International Trade in Endangered Species of Wild Fauna and Flora). Currently on the red list of the IUCN (International Union for Conservation of Nature), Arapaima gigas is classified as "Data Deficient", a category that characterizes a species for which there is not enough knowledge for it to fit into one of the categories. Even though it has been exploited commercially for many years, much basic life history data are still lacking. However, several recent publications have made significant advances of our understanding of growth and reproductive behavior (Arantes et al., 2010Arantes, C. C., L. Castello, D. J. Stewart, M. Certa & H. L. Queiroz. 2010. Population density, growth and reproduction of arapaima in an Amazonian river-floodplain. Ecology of Freshwater Fish, 19: 455-465.), population dynamics (Coutinho et al., 2010Coutinho, S. dos S. de S., L. Bevilacqua & H. L. de Queiroz. 2010. Population dynamics modeling of Arapaima gigas. Acta Amazonica, 40: 333-346.; Castello et al., 2011Castello, L., D. J. Stewart & C. C. Arantes. 2011. Modeling population dynamics and conservation of arapaima in the Amazon. Reviews in Fish Biology and Fisheries, 21: 623-640.), spatial and temporal distribution (Arantes et al., 2011) and trophic ecology (Watson et al., 2013Watson, L. C., D. J. Stewart & M. A. Teece. 2013. Trophic ecology of Arapaima in Guyana: giant omnivores in Neotropical floodplains. Neotropical Ichthyology, 11: 341-349.) as well as population structuring (Hrbek et al., 2005Hrbek, T., I. P. Farias, M. Crossa, I. Sampaio, J. I. R. Porto & A. Meyer. 2005. Population genetic analysis of Arapaima gigas, one of the largest freshwater fishes of the Amazon basin: implications for its conservation. Animal Conservation, 8: 297-308.; Hrbek et al., 2007Hrbek, T., M. Crossa & I. P. Farias. 2007. Conservation strategies for Arapaima gigas (Schinz, 1822) and the Amazonian várzea ecosystem. Brazilian Journal of Biology, 67: 909-917.; Araripe et al., 2013Araripe, J., P. S. do Rêgo, H. Queiroz, I. Sampaio & H. Schneider. 2013. Dispersal capacity and genetic structure of Arapaima gigas on different geographic scales using microsatellite markers. PloS One, 8: e54470.).

The first studies of reproductive biology of the Arapaima gigas in captivity were conducted by Oliveira (1944)Oliveira, C. E. de. 1944. Piscicultura amazônica. A Voz do Mar, 23: 104-106. and Fontenele (1948)Fontenele, O. 1948. Contribuição para o conhecimento da biologia do pirarucu, Arapaima gigas (Cuvier), em cativeiro (Actinopterygii, Osteoglossidae). Revista Brasileira de Biologia, 8: 445-459., who reported important biological data such as that it lays its eggs in lentic water, and belongs to the group of species with partitioned spawning and parental care. These data were confirmed by Imbiriba (1994)Imbiriba, E. P. 1994. Reprodução, larva e alevinagem de pirarucu (Arapaima gigas), recomendações básicas. Coleção CRIAR: 94., Imbiriba (2001)Imbiriba, E. P. 2001. Potencial de criação de pirarucu em cativeiro. Acta Amazonica, 31: 299-316., and Campos Baca (2001)Campos Baca, L. 2001. Historia Biológica del paiche o pirarucu Arapaima gigas (Cuvier) y bases para su cultivo em la Amazônia, Iquitos, Peru. Programa de Biodiversidad. Iquitos, Peru, Instituto de Investigaciones de la Amazonia Peruana.. Furthermore, these authors state that the spawning coincides with the rainy season, which, in the Amazon, extends from December to June. Imbiriba (2001)Imbiriba, E. P. 2001. Potencial de criação de pirarucu em cativeiro. Acta Amazonica, 31: 299-316. reports that the arapaima's reproduction is probably related with the water-level variation on the floodplains of the Amazon River. Sexual maturity occurs around the third to fifth year of life, when the fish reach an average length of 1.6 m and weigh 40 - 50 Kg (Lüling, 1964Lüling, V. 1964. Zur Biologie und Ökologie von Arapaima gigas (Pisces, Osteoglossidae). Zeitschrift für Morphologie und Ökologie der Tiere, 54: 436-530.; Godinho et al., 2005Godinho, H. P., J. E. Santos, P. S. Formagio & R. J. Guimarães-Cruz. 2005. Gonadal morphology and reproductive traits of the Amazonian fish Arapaima gigas (Schinz, 1822). Acta Zoologica, 86: 289-294.; Arantes et al., 2010Arantes, C. C., L. Castello, D. J. Stewart, M. Certa & H. L. Queiroz. 2010. Population density, growth and reproduction of arapaima in an Amazonian river-floodplain. Ecology of Freshwater Fish, 19: 455-465.). Fontenele (1953)Fontenele, O. 1953. Hábitos de desova do pirarucu, Arapaima gigas (Cuvier) (PISCES: Isospondyli, Arapaimidae), e a evolução de sua larva: 22. observed that during the spawning season secondary extra-genital sexual characteristics develop in the breeders, such as that males begin to present a dark coloring on the upper part of the head and flanks, and the abdomen and caudal peduncle turn reddish. However, as observed by Campos Baca (2001)Campos Baca, L. 2001. Historia Biológica del paiche o pirarucu Arapaima gigas (Cuvier) y bases para su cultivo em la Amazônia, Iquitos, Peru. Programa de Biodiversidad. Iquitos, Peru, Instituto de Investigaciones de la Amazonia Peruana., these features may not always be observed in all specimens. In both sexes only the left testicle/ovary is functional, while the right testicle/ovary is atrophied (Imbiriba et al., 1996Imbiriba, E. P., J. B. Lourenço Júnior, L. O. D. Moura Carvalho, L. B. Góes, D. Uliana & L. Brito Filho. 1996. Criação de pirarucu. Coleção CRIAR: 93.; Campos Baca, 2001Campos Baca, L. 2001. Historia Biológica del paiche o pirarucu Arapaima gigas (Cuvier) y bases para su cultivo em la Amazônia, Iquitos, Peru. Programa de Biodiversidad. Iquitos, Peru, Instituto de Investigaciones de la Amazonia Peruana.).

In accordance with Imbiriba (2001)Imbiriba, E. P. 2001. Potencial de criação de pirarucu em cativeiro. Acta Amazonica, 31: 299-316., the reproduction of the 'arapaima' occurs in a similar manner whether in nature or in captivity, including courtship, nest construction, mating and care of the offspring. In captivity, the sequence begins with a lack of interest in food on the part of the adults entering reproductive mode. In the following mating phase fights sometimes take place, although it not known if these are for the establishment of reproductive territories of they represent disputes for females. A four to five day period of tranquility follows, disturbed only at long intervals by the slow ascent of one or the other individual to breathe. Thereafter, the reproductive adults can be seen, one at a time, in a vertical position, head pointing down. During this period, the water becomes murky due to the excavation of the nest, which is circular, and commonly known as the "pan". The favorite place for preparation of the nest is the bottom of lakes, with clay soil and without vegetation, and where the water is shallow and still. After preparation of the nest, the female lays her eggs in it, and the male fertilizes them. During the nesting phase, the couple jealously guards the nest from up to a 10 m distance. After eclosion of the eggs, the larvae remain in the nest for five days, closely guarded by the male until absorption of the vitella/yolk. Brauner et al. (2004)Brauner, C. J., V. Matey, J. M. Wilson, N. J. Bernier & A. L. Val. 2004. Transition in organ function during the evolution of air-breathing; insights from Arapaima gigas, an obligate air-breathing teleost from the Amazon. Journal of Experimental Biology, 207: 1433-1438. also suggest that the males mouth brood or at least carry eggs in their mouths, a behavior observed by fishermen in distressed males. The larvae are black and swim over the head and dorsal region of the male protecting them, and after reaching one week of life, they become visible to the naked eye. In this period, they already surface to breathe (Imbiriba et al., 1996Imbiriba, E. P., J. B. Lourenço Júnior, L. O. D. Moura Carvalho, L. B. Góes, D. Uliana & L. Brito Filho. 1996. Criação de pirarucu. Coleção CRIAR: 93.). The juveniles are cared for exclusively the male, while the female circles in the vicinity of the group, defending against possible predators. The fry are cared for until their independence at between 35-40 cm in length and 4-5 months of age. During the free-swimming period, the couple swim slowly but continuously, guiding and protecting their offspring in what can be though of as a dynamic territory that can reach 1 ha/day. The size of this territory depends on the availability of food resources and refuge, including vegetation cover and shade. At the end of the parental care phase males are normally very thin, whereas females may have mature gonads and be able to spawn again in the very near future.

Based on what is known of the reproductive behavior of Arapaima gigas, it is assumed that the species is monogamous due to formation of couples and investment in parental care, however, monogamy has not been tested explicitly. There are many cases where apparently monogamous species in fact show high incidences of extra-pair mating, and thus levels of polygamy are high across a broad taxonomic range (Goossens et al., 1998Goossens, B., L. Graziani, L. P. Waits, S. Magnolon, J. Coulon, M. C. Bel, P. Taberlet & D. Allainé. 1998. Extra-pair paternity in the monogamous Alpine marmot revealed by nuclear DNA microsatellite analysis. Behavioral Ecology and Sociobiology, 43: 281-288.; Liebgold et al., 2006Liebgold, E. B., P. R. Cabe, R. G. Jaeger & P. L. Leberg. 2006. Multiple paternity in a salamander with socially monogamous behaviour. Molecular Ecology, 15: 4153-4160.; Ophir et al., 2008Ophir, A. G., S. M. Phelps, A. B. Sorin & J. O. Wolff. 2008. Social but not genetic monogamy is associated with greater breeding success in prairie voles. Animal Behaviour, 75: 1143-1154.; Sefc et al., 2008Sefc, K. M., K. Mattersdorfer, C. Sturmbauer & S. Koblmüller. 2008. High frequency of multiple paternity in broods of a socially monogamous cichlid fish with biparental nest defence. Molecular Ecology, 17: 2531-2543.).

The reproductive pattern and reproductive success of males is largely unknown in the majority of fish species. The physical nature of habitat normally hinders direct observation of reproductive behavior and kinship relations. Besides this, most fish species possess external fertilization, increasing the possibility of there occurring sperm competition (Stockley et al., 1997Stockley, P., M. J. G. Gage, G. A. Parker & A. P. Møller. 1997. Sperm competition in fishes: the evolution of testis size and ejaculate characteristics. American Naturalist, 149: 933-954.). The study of reproductive behavior in fish species with manifestation of polyandry (a female copulating with various males), polygyny (a male copulating with various females) or the observation of a monogamous reproductive system, can be efficiently evaluated and tested through the use of molecular methods, such as the analysis of DNA microsatellite data. The use of molecular approaches has allowed addressing many questions that would otherwise have been difficult using only field observations. Additionally, molecular markers can be used to quantify the incidence of "hidden" reproductive behaviors, that otherwise is unlikely to be observed in the field.

Considering the importance for the management and conservation of the species, in the present study we assess the mating system of the 'arapaima' under three experimental conditions, which provide different levels of a priori confidence in assuming monogamous system of mating. In the first instance, breeding pairs were isolated from all other individuals of their species (a priori very high probability of monogamy). In the second instance, breeding pairs were restricted to an area with other individuals of their own species, but no other breeding adults within the restricted area (a priori high probability of monogamy). In the third instance, offspring were collected from breeding pairs in nature (a priori medium probability of monogamy). To analyze mating system patterns under these three situations, we perform a parentage analysis of eight clutches of Arapaima gigas using 10 highly variable microsatellite markers developed for the 'arapaima' (Farias et al., 2003Farias, I. P., T. Hrbek, H. Brinkmann, I. Sampaio & A. Meyer. 2003. Characterization and isolation of DNA microsatellite primers for Arapaima gigas, an economically important but severely over-exploited fish species of the Amazon basin. Molecular Ecology Notes, 3: 128-130.).

Material and Methods

The samples utilized for this study come from an enclosed system, a breeding station, Fazenda Santo Antônio, in the municipality of Presidente Figueiredo near Manaus, Amazonas, a semi-natural system, a dammed lake in the municipality of Itacoatiara near Manaus, Amazonas, and a natural system, from the São Miguel Island, near the municipality of Santarém, Pará.

At the Fazenda Santo Antônio (FSA), 117 adult 'arapaima' individuals were maintained in breeding stations for approximately five years. Fourteen matched pairs were kept isolated from other adults and potential mates in individual 20 x 20 x 2 m tanks. From these fourteen pairs we sampled three groups of parents and offspring: Couple 1 (both parents plus 16 offspring), Couple 2 (both parents plus 19 offspring), and Couple 3 (both parents plus 23 offspring). With these samples we carried out a control analysis testing for the presence of any additional alleles not compatible with monogamy.

Two additional samples of offspring were obtained from Fazenda Santo Antônio Station, collected in different time periods, which we call brood FSA1 (25 offspring collected), and brood FSA2 (20 offspring collected). These samples were obtained from a communal 200 x 40 x 1.5 m tank containing 117 adults. In the collection period, only one brood under the care of one adult couple was observed at the breeding station.

In the Itacoatiara dammed lake, the breeders were kept for over seven years and there were 68 adults, but the proportion of males and females is not known. Two samples of offspring were obtained from the Itacoatiara dammed lake, collected in different periods, which we call brood Ita1 (20 offspring collected) and brood Ita2 (24 offspring collected). In the period of the two collections, only one brood under the care of one adult couple was observed in the breeding station.

In São Miguel Island (SMI), the offspring were collected from one couple and it was possible to collect only five individuals. Collection of additional individuals would have led to a major disturbance, resulting in the modification of behavior of the adults, and potential loss of all offspring. No other broods or spawning pairs were observed in the vicinity, however, given that this is a natural open system, other breeding adults may have been present in the area.

DNA extraction and genotyping. Total genomic DNA was extracted from scales and muscle tissues of the offspring preserved in 95% ethanol following the Qiagen Kit protocol. The microsatellite loci used were those developed by Farias et al. (2003)Farias, I. P., T. Hrbek, H. Brinkmann, I. Sampaio & A. Meyer. 2003. Characterization and isolation of DNA microsatellite primers for Arapaima gigas, an economically important but severely over-exploited fish species of the Amazon basin. Molecular Ecology Notes, 3: 128-130.: (CAm2, CAm13, CAm15, CAm16, CAm20, CTm3, CTm4, CTm5, CTm7 and Ctm8).

Polymerase Chain Reaction (PCR) was performed in 10 μl reaction volume containing 4.1 μl of MilliQ water, 0.8 µL of MgCl2 (25mM), 0.8 μl of dNTP (2,5 mM for each dNTP), 1.0 μl of 10x Buffer (100 mM of Tris-HCl, 500 mM of KCl), 1.0 μl of each primer (2.0 mM), 0.2 μl of Taq Polymerase (5U/μl) and 1.0 μl of DNA (~10 ng). The PCR reactions were performed in a thermal cycler Thermo (PXE 0.2) with the following steps: initial denaturation at 92º C for 2 minutes, followed by 35 cycles at 92º C for 40 seconds, seconds at the locus specific annealing temperature (according to Table 1 of Farias et al., 2003) and 72º C for 1 minute at 30 seconds. Last step was for 30 seconds a final extension at 72º C. PCR products generated with labeled primers were visualized on a MegaBACE 1000 Fragment Profiler v1.2 software (GE-Healthcare). Allele sizes were scored against the size standard ET-400 ROX (GE-Healthcare).

Table 1.
Main genetic pattern for each loci analysed in each brood of Arapaima gigas. A = number of alleles, Ho = observed heterozygosity, I = probability of genetic identity, Q = probability of paternity exclusion, IC = Joint probability of genetic identity at all loci, QC = Joint probability of paternal exclusion at all loci. Ita1 = Itacoatiara 1, Ita2 = Itacoatiara 2, FSA1 = Fazenda Santo Antônio 1, FSA2 = Fazenda Santo Antônio 2, SMI = São Miguel Island.

Allele size variation, observed and expected heterozygosity and Hardy-Weinberg equilibrium were implemented in the program Arlequin 3.5 (Excoffier & Lischer, 2010Excoffier, L. & H. E. L. Lischer. 2010. Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. Molecular Ecology Resources, 10: 564-567.). These analyses allowed inferring the genetic variability in each family group. Additional genetic parameters like distribution and frequency of each microsatellite loci were also implemented in Genetix v.4 (Belkhir et al., 2004Belkhir, K., P. Borsa, L. Chikhi, N. Raufaste & F. Bonhomme. 2004. GENETIX 4.05, logiciel sous Windows TM pour la génétique des populations. Laboratoire Génome, Populations, Interactions, CNRS UMR 5171, Université de Montpellier II, Montpellier (France).).

Parentage Analysis. We used two tests to estimate the efficiency of the markers in detecting multiple paternities in 'arapaima': the probability of genetic identity for each locus (I) and the probability of paternity exclusion for each locus (Q). Additionally, we estimated the joint probability of genetic identity â€" IC (Paetkau et al., 1995Paetkau, D., W. Calvert, I. Stirling & C. Strobeck. 1995. Microsatellite analysis of population structure in Canadian polar bears. Molecular Ecology, 4: 347-354.) and the joint probability of paternity exclusion method â€" QC (Weir, 1996Weir, B. S. 1996. Genetic Data Analysis II: Methods for Discrete Population Genetic Data. Sunderland, MA, Sinauer Associates.) for each test, respectively.

For the parentage analysis we implemented three methods. First, we used a simple method based in the Mendelian genetics of counting the number of observed alleles, and comparing it to the expected number of alleles in each nest (Myers & Zamudio, 2004Myers, E. M. & K. Zamudio. 2004. Multiple paternity in an aggregate breeding amphibian: the effect of reproductive skew on estimates of male reproductive success. Molecular Ecology, 13: 1951-1963.). The maximum number of expected alleles in each nest under the assumption of monogamy is four, unless one observes homozygous progeny (Fitzsimmons, 1998Fitzsimmons, N. N. 1998. Single paternity of clutches and sperm storage in the promiscuous green turtle (Chelonia mydas). Molecular Ecology, 7: 575-584.; Valenzuela, 2000Valenzuela, N. 2000. Multiple paternity in side-neck turtles Podocnemis expansa: evidence from microsatellite DNA data. Molecular Ecology, 9: 99-105.). Each homozygous genotype indicates that alleles are shared between parents, and therefore the number of expected alleles decrease to three with one homozygous genotype, and to two with two homozygous genotypes.

Considering that the counting method does not consider the combination of all alleles across the different loci, and is therefore likely to underestimate the number of sires contributing to the brood, we carried out a second analysis using the program KINALYZER (Berger-Wolf et al., 2007Berger-Wolf, T. Y., S. I. Sheikh, B. DasGupta, M. V Ashley, I. C. Caballero, W. Chaovalitwongse & S. L. Putrevu. 2007. Reconstructing sibling relationships in wild populations. Bioinformatics, 23: i49-i56.) which analyzes all loci simultaneously making use of a minimum 2-allele set cover approach based on Mendelian properties and parsimoniously finds the smallest number of sibling groups (Berger-Wolf et al., 2007Berger-Wolf, T. Y., S. I. Sheikh, B. DasGupta, M. V Ashley, I. C. Caballero, W. Chaovalitwongse & S. L. Putrevu. 2007. Reconstructing sibling relationships in wild populations. Bioinformatics, 23: i49-i56.). The inference of full sibling groups (groups of individuals sharing the same two parents) is possible even with no a priori knowledge of parental genotypes, as this program performs a heuristic analysis based on allelic inheritance in diploid organisms â€" every offspring inherits one allele at each locus from each of its parents.

We also use the program ML-Relate (Kalinowski et al., 2006Kalinowski, S. T., A. P. Wagner & M. L. Taper. 2006. ML-RELATE: a computer program for maximum likelihood estimation of relatedness and relationship. Molecular Ecology Notes, 6: 576-579.) to estimate the relationships among individuals of the same clutch. This program uses maximum likelihood estimates of relatedness to classify pairs of individuals as unrelated, half-sibs, full-sibs, or parent-offspring relationships (Jones et al., 2010Jones, A. G., C. M. Small, K. A. Paczolt & N. L. Ratterman. 2010. A practical guide to methods of parentage analysis. Molecular Ecology Resources, 10: 6-30.).

Results

Statistical power and genetic diversity parameters. The microsatellite loci had high statistical power to discriminate between individuals and exclude paternity. The probability of genetic identity (I) varied from 0.09657709 (CTm7 in Ita2) to 0.50743744 (Cam13 in FSA) while the joint probability of genetic identity (IC) ranged from 3.5 x 10-6 (FSA2) to less than 1 x 10-9 (FSA2), demonstrating the high discriminatory power of these loci. On the other hand, the probability of paternity exclusion (Q), varied from 0.15687672 (CTm8 in FSA) to 0.56049347 (CTm7 in Ita2) while the joint probability of exclusion (QC) ranged from 0.975539087 (FSA1) to 0.99958039 (SMI) again demonstrating the high discriminatory power of these loci (Table 1). In other words, the set of markers used for this analysis of the relationship among the 'arapaima' offspring proved to be very robust for the analyses performed.

The average genetic diversity for all the loci in each locality varied from 0.486 (FSA) to 0.583 (Itacoatiara 2) (Table 2). Wright's endogamy index (FIS) was low and not significant in any of the broods. Of the ten microsatellite loci genotyped, the loci CTm4 and CTm8 had the lowest number of alleles in all the broods (Table 1). Locus Cam20 was monomorphic in the FSA brood, and CTm4 and CTm8 were monomorphic in the SMI brood. In most cases, the remaining loci showed a considerable degree of polymorphism.

Table 2.
Main genetic pattern for each brood analysed in Arapaima gigas. Note: HO= observed heterozygosity, HE= expected heterozygosity, FIS = fixation index.

Parentage Analysis. The control analyzes of the three pairs and their offspring sampled from FSA showed that all the juveniles were offspring of their hypothetical parents. There were no new alleles detected, and all loci and alleles followed correct Mendelian inheritance patterns (results not shown). Monogamy was, therefore, detected in the experimental set-up where monogamy was the only expected system of mating, given that the adult couple was isolated from all other adults.

Considering the set of alleles found in the groups of young in each nest, and assuming absence of information on the parents, we began with a simple analysis based on the patterns of Mendelian inheritance. Considering that each offspring possesses 50% of the maternal and 50% of the paternal genome, and each parent can be heterozygous at any given locus, the total number of alleles at a locus observed in each brood can be at most four. The presence of five or more alleles in a brood suggests a possible extra contribution in addition to the parents (two maternal alleles, two paternal alleles, and the additional allele(s) from an extra-pair individual). Out of the loci analyzed in this study, five presented five or more alleles in the broods analyzed: Ita1, locus CAm20 (7 alleles); Ita2, loci CAm20 (7 alleles), CTm5 (5 alleles), CTm7 (5 alleles); FSA2, loci CTm5 (7 alleles), CTm7 (5 alleles); SMI, loci CAm2 (7 alleles), CAm13 (5 alleles) (Table 3). The FSA1 brood had a maximum of four alleles at each locus.

Table 3.
Alleles frequencies (allele size/frequency) and number of homozygotes (allele size/number of homozygote individuals) observed in each brood of Arapaima gigas. The number of sib-groups inferred in Kinalyzer is also listed for each brood. N = number of offspring. NHm = number of homozygote individuals.

The presence of homozygous genotypes in the offspring is also important in inferring of how many adult individuals could be contributing genetically to each group of offspring (Table 3), mainly in cases where there occur three or four alleles per locus. For example, a shared parental allele may be inferred when an offspring is homozygous for a given locus (e.g., offspring genotype = AA; parental genotype = Aâ€") which leads to the expectation of only three alleles being observed in a brood. Similarly, the complete parental genotype may be inferred, when the offspring are homozygous for two different alleles (e.g., offspring genotypes = AA and BB; parental genotype = AB) which leads to the expectation of only two alleles being observed in a brood. Additional alleles will have to have been contributed via extra-pair mating events. For example, upon analysis of the locus CAm2 in the brood Ita1, we observe three alleles (298, 312, 322) and two types of homozygotes: five offspring 312/312 and one offspring 298/298. Based on the presence of these homozygotes, we can infer that the genotype of the parents is 298/312, and that the appearance of allele 322 observed in the brood took place by extra-pair contribution of at least one individual 322/-. This individual may have been a homozygote (322/322), a heterozygote sharing the second allele with the other same sex parent (322/312 or 322/298), or a heterozygote where the second allele was not passed on or not detected in our sample. Using this analysis, extra-pair contribution was also detected in the broods Ita2, FSA1 and FSA2 (Table 3).

In summary, extra-pair contribution was observed in all broods. In the broods Ita1, Ita2, FSA2 and SMI the presence of five or more alleles was inferred, while extra alleles were inferred in the broods Ita1, Ita2, FSA1, and FSA2.

The presence of multiple sib-groups in all five broods was inferred in the program KINALYZER. At least two sib-groups per brood were inferred (Table 3). Analyses in ML-Relate also indicate the presence of multiple sib-groups in each group resulting from the contribution of multiple unrelated individuals (results not shown). The ML-Relate analysis thus indicated that in each group more than one female and more than one male contributed to each brood. All analyses therefore suggest that Arapaima gigas does not have a monogamous system of mating.

Discussion

The analysis of microsatellite data provides powerful information on the relatedness and parentage of individuals, including fishes (Avise et al., 2002Avise, J. C., A. G. Jones, D. Walker, J. A. DeWoody, B. Dakin, A. C. Fiumera, D. E. Fletcher, M. Mackiewicz, D. E. Pearse, B. Porter & S. D. Wilkins. 2002. Genetic mating systems and reproductive natural histories of fishes: Lessons for ecology and evolution. Annual Review of Genetics, 36: 19-45.). Various authors have discussed the diverse reproductive behaviors and tactics found in fishes through genetic assessments of parentage (Taborsky, 2001Taborsky, M. 2001. The evolution of bourgeois, parasitic, and cooperative reproductive behaviors in fishes. Journal of Heredity, 92: 100-110., 2008Taborsky, M. 2008. Alternative reproductive tactics in fish. Pp. 251-299 In: Oliveira, R. F., M. Taborsky & H. J. Brockmann (Eds.). Alternative Reproductive Tactics. Cambridge, UK, Cambridge University Press.; Wilson & Ferguson, 2002Wilson, A. B. & M. M. Ferguson. 2002. Molecular pedigree analysis in natural populations of fishes: approaches, applications, and practical considerations. Canadian Journal of Fisheries and Aquatic Sciences, 59: 1696-1707.; Hain & Neff, 2007Hain, T. J. A. & B. D. Neff. 2007. Multiple paternity and kin recognition mechanisms in a guppy population. Molecular Ecology, 16: 3938-3946.), and their findings have revolutionized the study of reproductive behavior, revealing that individuals of many species engage in extra-pair copulations, including parasitic spawning behavior between satellite and territorial males such as sneaking, egg piracy and female mimicry of cooperative breeding between satellite and territorial males.

The results obtained in this study suggest that Arapaima gigas broods have contributions from multiple unrelated males and females. Therefore, we can reject the hypothesis that Arapaima gigas is a monogamous species in natural situations. Given that multiple male and female individuals contribute to each brood, the spawning behavior of Arapaima gigas is likely to be complex. Since males are the primary caregivers and guardians of the brood, males may chose to mate with multiple females, and the contribution of additional males to the brood would then likely to be result of parasitic behavior of satellite males (Roldán & Soler, 2011Roldán, M. & M. Soler. 2011. Parental-care parasitism: how do unrelated offspring attain acceptance by foster parents? Behavioral Ecology, 22: 679-691.). Another, not necessarily exclusive possibility, is the existence of alloparental care, where other pairs spawn in the nest of the caregiver male (Wisenden, 1999Wisenden, B. D. 1999. Alloparental care in fishes. Reviews in Fish Biology and Fisheries, 9: 45-70.). Alloparental care is thus defined as care of unrelated young.

There are many reasons for why a male would choose to mate with multiple females (Ridley, 1978Ridley, M. 1978. Paternal care. Animal Behaviour, 26: 904-932.). A male may choose to mate with an additional female or females, if they are likely to increase his fitness. The male does not need to give up the current brood to mate again, and actually may even become more attractive for females if already with offspring (Ridley, 1978Ridley, M. 1978. Paternal care. Animal Behaviour, 26: 904-932.). Given that parental care in majority of fishes including Arapaima gigas takes form of protecting against predators rather provisioning of offspring (Wisenden, 1999Wisenden, B. D. 1999. Alloparental care in fishes. Reviews in Fish Biology and Fisheries, 9: 45-70.), the increase in size in the clutch has minimal effect on parental investment. Thus, it may be advantageous for the male to solicit multiple matings, thus increasing the total number of offspring, and hedging his bets, and consequently increasing his fitness through increasing his chances of a larger portion of his offspring surviving. Alloparental care may also increase reproductive fitness. While alloparental care appears to be disadvantageous (Wisenden, 1999Wisenden, B. D. 1999. Alloparental care in fishes. Reviews in Fish Biology and Fisheries, 9: 45-70.), this behavior may have some advantages for the male caregiver. One of the possible advantages of alloparental care is the dilution of predation (Wisenden, 1999Wisenden, B. D. 1999. Alloparental care in fishes. Reviews in Fish Biology and Fisheries, 9: 45-70.): if predation pressure is high, having unrelated fry mixed with the male's offspring lessens the chance of his young being preyed upon. The percentage of the offspring not related to the male caregiver will exist as an equilibrium between investment into unrelated offspring and probability of predation of one owns offspring. However, since in fishes parental care is usually only a protection against predators, and involves no investment in provisioning of offspring (McKaye, 1981McKaye, K. R. 1981. Natural selection and the evolution of interespecific brood care in fishes. Pp. 173-183. In: Alexander, R. D. & D. W. Tinkle (Eds.). Natural Selection and Social Behavior: Recent Research and New Theory. New York, NY, USA, Chiron.), an increase in the size of a brood may not be much of a disadvantage for caregivers, since increase in brood size increases energetic investment only marginally (McKaye, 1981McKaye, K. R. 1981. Natural selection and the evolution of interespecific brood care in fishes. Pp. 173-183. In: Alexander, R. D. & D. W. Tinkle (Eds.). Natural Selection and Social Behavior: Recent Research and New Theory. New York, NY, USA, Chiron.). Unrelated fry in the male's clutch may also be the result of parental-care parasitism, but may have similar indirect benefits as alloparental care. Parental-care parasitism is defined as an interaction in which an individual (extra-pair) obtains reproductive benefits while reducing or completely eliminating its own costs of parenting by exploiting any type of offspring care provided by other individuals (Roldán & Soler, 2011Roldán, M. & M. Soler. 2011. Parental-care parasitism: how do unrelated offspring attain acceptance by foster parents? Behavioral Ecology, 22: 679-691.), such as the protection to offspring provided by the male caregiver in Arapaima gigas.

The results observed here correspond mainly to an artificial environment (breeding station) or semi-artificial (dammed reservoir) that do not possess many of the same characteristics as the natural environment. However, contributions of multiple males and females to a brood were also observed in a natural situation in the São Miguel Island (SMI) brood. Therefore multiple paternity is likely a general phenomenon of the reproductive system of Arapaima gigas.

We believe these results will stimulate more complex studies of the mating system of the 'arapaima' and will encourage fish farmers to consider possible alternative reproductive tactics for the management in captivity. From a management perspective, this behavior may be exploited to maintain genetic diversity in captive and as well in wild populations of Arapaima gigas. Polygamous mating systems maintain genetic diversity and heterozygosity of populations, and thus may be exploited in the genetic management of captive populations and breeding stocks. The discovery of polygamous mating system may also eliminate the perceived need and the associated cost of isolating reproducing couples from the rest of the adults.

Acknowledgments

This research was supported by grants from CNPq/CT-Amazônia (n.o 553036/2005-0 and n.o 554057/2006-9) and FINEP/DARPA (Convênio n.o 01.09.0472.00) to IPF and CNPq/MPA No. 406893/2012-8 to TH. IPF and TH are supported by a Bolsa de Pesquisa scholarship from CNPq and YA was supported by a fellowship from SUFRAMA. This study formed a portion of YA's Masters Project in the Biotechnology graduate program at UFAM.

References

  • Arantes, C. C., L. Castello, D. J. Stewart, M. Certa & H. L. Queiroz. 2010. Population density, growth and reproduction of arapaima in an Amazonian river-floodplain. Ecology of Freshwater Fish, 19: 455-465.
  • Arantes, C. C., L. Castello, M. Cetra & A. Schilling. 2011. Environmental influences on the distribution of arapaima in Amazon floodplains. Environmental Biology of Fishes, 96: 1257-1267.
  • Araripe, J., P. S. do Rêgo, H. Queiroz, I. Sampaio & H. Schneider. 2013. Dispersal capacity and genetic structure of Arapaima gigas on different geographic scales using microsatellite markers. PloS One, 8: e54470.
  • Avise, J. C., A. G. Jones, D. Walker, J. A. DeWoody, B. Dakin, A. C. Fiumera, D. E. Fletcher, M. Mackiewicz, D. E. Pearse, B. Porter & S. D. Wilkins. 2002. Genetic mating systems and reproductive natural histories of fishes: Lessons for ecology and evolution. Annual Review of Genetics, 36: 19-45.
  • Belkhir, K., P. Borsa, L. Chikhi, N. Raufaste & F. Bonhomme. 2004. GENETIX 4.05, logiciel sous Windows TM pour la génétique des populations. Laboratoire Génome, Populations, Interactions, CNRS UMR 5171, Université de Montpellier II, Montpellier (France).
  • Berger-Wolf, T. Y., S. I. Sheikh, B. DasGupta, M. V Ashley, I. C. Caballero, W. Chaovalitwongse & S. L. Putrevu. 2007. Reconstructing sibling relationships in wild populations. Bioinformatics, 23: i49-i56.
  • Brauner, C. J., V. Matey, J. M. Wilson, N. J. Bernier & A. L. Val. 2004. Transition in organ function during the evolution of air-breathing; insights from Arapaima gigas, an obligate air-breathing teleost from the Amazon. Journal of Experimental Biology, 207: 1433-1438.
  • Campos Baca, L. 2001. Historia Biológica del paiche o pirarucu Arapaima gigas (Cuvier) y bases para su cultivo em la Amazônia, Iquitos, Peru. Programa de Biodiversidad. Iquitos, Peru, Instituto de Investigaciones de la Amazonia Peruana.
  • Castello, L., D. J. Stewart & C. C. Arantes. 2011. Modeling population dynamics and conservation of arapaima in the Amazon. Reviews in Fish Biology and Fisheries, 21: 623-640.
  • Coutinho, S. dos S. de S., L. Bevilacqua & H. L. de Queiroz. 2010. Population dynamics modeling of Arapaima gigas. Acta Amazonica, 40: 333-346.
  • Excoffier, L. & H. E. L. Lischer. 2010. Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. Molecular Ecology Resources, 10: 564-567.
  • Farias, I. P., T. Hrbek, H. Brinkmann, I. Sampaio & A. Meyer. 2003. Characterization and isolation of DNA microsatellite primers for Arapaima gigas, an economically important but severely over-exploited fish species of the Amazon basin. Molecular Ecology Notes, 3: 128-130.
  • Fitzsimmons, N. N. 1998. Single paternity of clutches and sperm storage in the promiscuous green turtle (Chelonia mydas). Molecular Ecology, 7: 575-584.
  • Fontenele, O. 1948. Contribuição para o conhecimento da biologia do pirarucu, Arapaima gigas (Cuvier), em cativeiro (Actinopterygii, Osteoglossidae). Revista Brasileira de Biologia, 8: 445-459.
  • Fontenele, O. 1953. Hábitos de desova do pirarucu, Arapaima gigas (Cuvier) (PISCES: Isospondyli, Arapaimidae), e a evolução de sua larva: 22.
  • Godinho, H. P., J. E. Santos, P. S. Formagio & R. J. Guimarães-Cruz. 2005. Gonadal morphology and reproductive traits of the Amazonian fish Arapaima gigas (Schinz, 1822). Acta Zoologica, 86: 289-294.
  • Goossens, B., L. Graziani, L. P. Waits, S. Magnolon, J. Coulon, M. C. Bel, P. Taberlet & D. Allainé. 1998. Extra-pair paternity in the monogamous Alpine marmot revealed by nuclear DNA microsatellite analysis. Behavioral Ecology and Sociobiology, 43: 281-288.
  • Goulding, M. 1980. Fishes and the Forest: Explorations in Amazonian Natural History. Los Angeles, CA, University of California Press.
  • Hain, T. J. A. & B. D. Neff. 2007. Multiple paternity and kin recognition mechanisms in a guppy population. Molecular Ecology, 16: 3938-3946.
  • Hrbek, T., M. Crossa & I. P. Farias. 2007. Conservation strategies for Arapaima gigas (Schinz, 1822) and the Amazonian várzea ecosystem. Brazilian Journal of Biology, 67: 909-917.
  • Hrbek, T., I. P. Farias, M. Crossa, I. Sampaio, J. I. R. Porto & A. Meyer. 2005. Population genetic analysis of Arapaima gigas, one of the largest freshwater fishes of the Amazon basin: implications for its conservation. Animal Conservation, 8: 297-308.
  • Imbiriba, E. P. 1994. Reprodução, larva e alevinagem de pirarucu (Arapaima gigas), recomendações básicas. Coleção CRIAR: 94.
  • Imbiriba, E. P. 2001. Potencial de criação de pirarucu em cativeiro. Acta Amazonica, 31: 299-316.
  • Imbiriba, E. P., J. B. Lourenço Júnior, L. O. D. Moura Carvalho, L. B. Góes, D. Uliana & L. Brito Filho. 1996. Criação de pirarucu. Coleção CRIAR: 93.
  • Jones, A. G., C. M. Small, K. A. Paczolt & N. L. Ratterman. 2010. A practical guide to methods of parentage analysis. Molecular Ecology Resources, 10: 6-30.
  • Kalinowski, S. T., A. P. Wagner & M. L. Taper. 2006. ML-RELATE: a computer program for maximum likelihood estimation of relatedness and relationship. Molecular Ecology Notes, 6: 576-579.
  • Liebgold, E. B., P. R. Cabe, R. G. Jaeger & P. L. Leberg. 2006. Multiple paternity in a salamander with socially monogamous behaviour. Molecular Ecology, 15: 4153-4160.
  • Lüling, V. 1964. Zur Biologie und Ökologie von Arapaima gigas (Pisces, Osteoglossidae). Zeitschrift für Morphologie und Ökologie der Tiere, 54: 436-530.
  • McKaye, K. R. 1981. Natural selection and the evolution of interespecific brood care in fishes. Pp. 173-183. In: Alexander, R. D. & D. W. Tinkle (Eds.). Natural Selection and Social Behavior: Recent Research and New Theory. New York, NY, USA, Chiron.
  • Myers, E. M. & K. Zamudio. 2004. Multiple paternity in an aggregate breeding amphibian: the effect of reproductive skew on estimates of male reproductive success. Molecular Ecology, 13: 1951-1963.
  • Oliveira, C. E. de. 1944. Piscicultura amazônica. A Voz do Mar, 23: 104-106.
  • Ophir, A. G., S. M. Phelps, A. B. Sorin & J. O. Wolff. 2008. Social but not genetic monogamy is associated with greater breeding success in prairie voles. Animal Behaviour, 75: 1143-1154.
  • Paetkau, D., W. Calvert, I. Stirling & C. Strobeck. 1995. Microsatellite analysis of population structure in Canadian polar bears. Molecular Ecology, 4: 347-354.
  • Ridley, M. 1978. Paternal care. Animal Behaviour, 26: 904-932.
  • Roldán, M. & M. Soler. 2011. Parental-care parasitism: how do unrelated offspring attain acceptance by foster parents? Behavioral Ecology, 22: 679-691.
  • Sefc, K. M., K. Mattersdorfer, C. Sturmbauer & S. Koblmüller. 2008. High frequency of multiple paternity in broods of a socially monogamous cichlid fish with biparental nest defence. Molecular Ecology, 17: 2531-2543.
  • Stockley, P., M. J. G. Gage, G. A. Parker & A. P. Møller. 1997. Sperm competition in fishes: the evolution of testis size and ejaculate characteristics. American Naturalist, 149: 933-954.
  • Taborsky, M. 2001. The evolution of bourgeois, parasitic, and cooperative reproductive behaviors in fishes. Journal of Heredity, 92: 100-110.
  • Taborsky, M. 2008. Alternative reproductive tactics in fish. Pp. 251-299 In: Oliveira, R. F., M. Taborsky & H. J. Brockmann (Eds.). Alternative Reproductive Tactics. Cambridge, UK, Cambridge University Press.
  • Valenzuela, N. 2000. Multiple paternity in side-neck turtles Podocnemis expansa: evidence from microsatellite DNA data. Molecular Ecology, 9: 99-105.
  • Watson, L. C., D. J. Stewart & M. A. Teece. 2013. Trophic ecology of Arapaima in Guyana: giant omnivores in Neotropical floodplains. Neotropical Ichthyology, 11: 341-349.
  • Weir, B. S. 1996. Genetic Data Analysis II: Methods for Discrete Population Genetic Data. Sunderland, MA, Sinauer Associates.
  • Wilson, A. B. & M. M. Ferguson. 2002. Molecular pedigree analysis in natural populations of fishes: approaches, applications, and practical considerations. Canadian Journal of Fisheries and Aquatic Sciences, 59: 1696-1707.
  • Wisenden, B. D. 1999. Alloparental care in fishes. Reviews in Fish Biology and Fisheries, 9: 45-70.

Publication Dates

  • Publication in this collection
    31 Mar 2015
  • Date of issue
    Jan-Mar 2015

History

  • Received
    04 Feb 2014
  • Accepted
    06 Oct 2014
Sociedade Brasileira de Ictiologia Neotropical Ichthyology, Núcleo de Pesquisas em Limnologia, Ictiologia e Aquicultura, Universidade Estadual de Maringá., Av. Colombo, 5790, 87020-900, Phone number: +55 44-3011-4632 - Maringá - PR - Brazil
E-mail: neoichth@nupelia.uem.br