The Goliath grouper ( Epinephelus itajara ) is one of the most endangered species of fish of the subfamily Epinephelinae. Slow to develop and mature, and dependent on mangrove habitats for breeding, the species also suffers intense harvesting, which has reduced drastically in numbers in many areas. To contribute to the understanding of the characteristics of E. itajara populations, we conducted a molecular genetics study of the species, focusing on populations from the Northern Brazilian coast. The mtDNA control region (D-loop) of 116 individuals from five localities (Bragança, Ajuruteua, Parnaíba, Fortaleza and Natal) was analysed, and a sequence of 499 base pairs identified. Analyses of the sequences indicated that genetic variability was generally lower in E. itajara than in other endangered species of the genus. AMOVA found no significant grouping structure among the populations. Nested Clade Analysis revealed a significant association between genetic variability and geographic distribution among only three populations (Ajuruteua, Parnaíba and Natal). Genetic diversity was higher in populations from the Amazon region, which may be related to the better conservation of mangrove habitats in this area. Therefore, the present study could be used for the implementation of conservation and management measures in order to protect and consolidate these populations.
Goliath grouper; Epinephelus itajara; control region; mtDNA; population genetics
Genetic characterisation of populations of the critically endangered Goliath grouper (Epinephelus itajara, Serranidae) from the Northern Brazilian coast through analyses of mtDNA
Gláucia C. Silva-Oliveira1
Péricles S. do Rêgo1
1Laboratório de Genética e Biologia Molecular, Instituto de Estudos Costeiros, Universidade Federal do Pará, Bragança, PA, Brazil
2Centro de Investigação em Biodiversidade e Recursos Genéticos, Campus Agrário de Vairão, Universidade do Porto, Vairão, Portugal
The Goliath grouper (Epinephelus itajara) is one of the most endangered species of fish of the subfamily Epinephelinae. Slow to develop and mature, and dependent on mangrove habitats for breeding, the species also suffers intense harvesting, which has reduced drastically in numbers in many areas. To contribute to the understanding of the characteristics of E. itajara populations, we conducted a molecular genetics study of the species, focusing on populations from the Northern Brazilian coast. The mtDNA control region (D-loop) of 116 individuals from five localities (Bragança, Ajuruteua, Parnaíba, Fortaleza and Natal) was analysed, and a sequence of 499 base pairs identified. Analyses of the sequences indicated that genetic variability was generally lower in E. itajara than in other endangered species of the genus. AMOVA found no significant grouping structure among the populations. Nested Clade Analysis revealed a significant association between genetic variability and geographic distribution among only three populations (Ajuruteua, Parnaíba and Natal). Genetic diversity was higher in populations from the Amazon region, which may be related to the better conservation of mangrove habitats in this area. Therefore, the present study could be used for the implementation of conservation and management measures in order to protect and consolidate these populations.
The Goliath grouper is much sought after by both sport and commercial fishermen (Szpilman, 2000). In general, the species is targeted during the shoaling phase, which permits the harvesting of large numbers of individuals. This was the main cause of the decline E. itajara populations in Central America in the 1950s (Randall, 1983; Coleman and Koening Laboratory Research). The formation of reproductive shoals in predictable periods, the dependence of the species on mangrove habitats and its slow maturation rates all contribute to its vulnerability (Sadovy and Eklund, 1999; Frias-Torres, 2006). According to Morris et al. (2000), around 40% of the epinephelids are currently listed at some level of extinction risk by the International Union for Conservation of Nature and Natural Resources (IUCN, 2007). At the present time, the Goliath grouper is listed as critically endangered, and is potentially vulnerable to drastic reductions in numbers in the near future (Tak-Chuen and Ferrera, 2006, in IUCN, 2007). Given its conservation status, E. itajara was the first species of marine fish to be protected under Brazilian law, starting in 2002 (IBAMA, portaria nº 121 20/09/2002 and portaria n. 42, 19/09/2007). Similar measures have already been taken since 1990 in the United States and 1993 in the Caribbean (Tak-Chuen and Ferrera, 2006, in IUCN, 2007).
In the light of the present status of E. itajara, there is an urgent need for population studies, in particular on genetic variability. However, no data are available on the genetic characteristics of this species, but some information is available for closely-related species. Three studies of population genetics are available for other members of the genus Epinephelus (Rivera et al., 2004; Zatcoff et al., 2004; Maggio et al., 2006), although they are based on different methods and molecular markers. Given this situation, the aims of the present study were to estimate the genetic variability of E. itajara populations from northern Brazilian coast based on the analyses of the control region (D-loop) of the mitochondrial genome (mtDNA), investigate migratory patterns and demographic and structural traits of the different populations.
In the present study, 116 tissue samples of western Atlantic E. itajara were collected from five sites on the northern Brazilian coast, through donations from fishermen and fishmongers. Amazonian populations are represented by samples from the Ajuruteua peninsula (00°56'33” S and 47°06'58” W; n = 20) and the municipal fish market in Bragança (1°3'46” S and 46°46'22” W; n = 73), both in the state of Pará (Figure 1). In addition, samples were also obtained from the cities of Fortaleza, in Ceará state (3°43'1” S and 38°32'34” W; n = 2), Parnaíba, Piauí (2°54'18” S and 41°46'37” W; n = 12) and Natal, in Rio Grande do Norte (Rio Potengí, 5°45' 22" S and 35°12'5" W; n = 9). Samples were stored in 95% ethanol until processing for the isolation of genetic material.
Total DNA was extracted from the samples following a phenol-chloroform method (Sambrook et al., 1989). Taxon identification was confirmed with the assistance of a species-specific molecular marker, in this case, a segment of 16S rRNA from the mitochondrial DNA (Palumbi et al., 1991). The sequences of all samples were compared with the complete sequence of a juvenile E. itajara deposited in the tissue collection of the Genetic and Molecular Biology Laboratory of the Bragança campus of the Universidade Federal do Pará (UFPA). Polymerase chain reactions (PCR) were conducted in order to isolate the D-loop fragment, using primers A-1 (Lee et al., 1995) and Perc12S1R (Santa Brígida et al., 2007). Approximately 100 ng of total DNA were used for each reaction, together with 4 μL dNTP (1.25 mM), 2.5 μL of 10X buffer (Invitrogen - Tris-HCl and KCl, pH 7.8), 1.5 μL MgCl2 (50 mM), 0.25 μL of each primer (200 ng/μL), 0.25 μL of Taq polymerase (5 U/μL - Invitrogen) and extra-purified water to complete a final volume of 25 μL. The conditions required for amplification were 94 °C for 3 min, followed by 30 cycles at 94 °C for 1 min, 47.5 °C for 1 min, and 72 °C for 1 min, with an additional final extension cycle at 72 °C for 5 min. The PCR products were purified using an ExoSAP-IT kit (Amersham Pharmacia Biotech, Inc., UK), following the manufacturer's protocol. The samples were sequenced separately on the basis of the primers used for PCR, as well as an internal primer (EIT-INT 5'-GAATATTCCTTCAAC ATTAC-3') designed for the present study. The dideoxyterminal method was used for sequencing, with the commercial DYEnamicTM ET dye Terminator kit (Amersham Pharmacia Biotech, Inc., UK) in a MegaBACE 1000 sequencer (Amersham Pharmacia Biotech, Inc., UK).
Sequences were aligned using the BIOEDIT program (Hall, 1999). For analysis, the sequences were grouped in two distinct files. The first file was used for the calculation of genetic variability, and contained all the samples from the five localities (n = 116), whereas the second (n = 41) was reserved for the analyses of population structuring, and excluded the samples from Bragança (73 individuals), the origin of which was unknown, and Fortaleza, which included only two specimens.
The first group was used to calculate the percentage of variable and informative sites for parsimony analysis, obtained with the MEGA program, version 3.1 (Kumar et al., 2004). The number and respective frequencies of the different haplotypes were estimated using DnaSP, version 4.0 (Rozas et al., 2003). Haplotypic (h; Nei, 1987) and nucleotidic (π; Nei, 1987) diversities were also calculated in Arlequin, version 3.1 (Excoffier et al., 2005).
Population structuring and demographic history were analysed in the DnaSP and Arlequin programs. The genetic differentiation among the three selected populations (Ajuruteua - AJU, Parnaíba - PAR and Natal - NAT) was analysed through pairwise estimates of Fst, the significance of which was tested in 1000 permutations.
Among-population partitioning of genetic diversity was tested using the hierarchical analysis of molecular variance, AMOVA (Excoffier et al., 1992). Deviations from theoretical neutral evolution were verified through Fs test (Fu 1997), based on 1000 permutations. A variety of tests were used to detect possible processes of demographic fluctuation: pairwise mismatch distribution of the haplotypes (Rogers and Harpending, 1992); the sum of the squared deviations (SSD) and its respective p value, which represents the probability of observing a match between the model of expansion and the random distribution; estimate of the expansion time of the population (t), following the Rogers and Harpending (1992) model; and the parameters τ (unit of population growth), θ0 and θ1 (θ values before and after expansion, respectively) based on the mismatch distribution, according to the equation t = τ/2u, where u is the mutation rate per sequence per generation (Rogers, 1995; Schneider and Excoffier, 1999). The D-loop mutation rate adopted here was 3.6 × 10-8 (Donaldson and Wilson, 1999), which was the same one used for other epinephilids (Rivera et al., 2004) and generation time for E. itajara was considered to be six and a half years (Bullock et al., 1992). Effective population size (Ne) was inferred given the estimate of diversity index (θ = 2Neμ) based on Watterson (1975) and Tajima (1983) estimates of θ (θS and θ, respectively) in Arlequin 3.1.
The haplotype network was defined using the TCS program, version 1.17 (Clement et al., 2000), in which it was nested manually into increasingly inclusive clades (or nestings) following the rules described by Templeton et al. (1987) and Templeton and Sing (1993). This nested haplotype network was used in a nested clade analysis (NCA: Templeton et al., 1987; Templeton and Sing, 1993). Interpretation of the results was based on the inference key available at the site http://darwin.uvigo.es/download/geodisKey_14Jul04.pdf. An analysis based on the correlation between the genetic and geographic distances among populations was conducted using the Mantel test run in Arlequin 3.1.
A total of 499 base pairs were obtained starting from the 3' region of the mtDNA control region in 116 E. itajara specimens from the northern Brazilian coast. This fragment of the mtDNA presented only 17 variable sites with just eight informative. The simple deletions were used in the analyses, although a deletion of eight base pairs was observed in one individual (position 254-261), and was considered as a single substitution. A total of 27 haplotypes were identified on the basis of these variable sites and gaps (Table 1). Four of these (H1, H2, H3, and H11) were well represented, whereas 15 were singletons. The population from Parnaíba presented only two haplotypes (H1 and H2), both of which were shared with Ajuruteua, Natal and Bragança.
As expected due to the largest sample size among the studied areas, Bragança had the largest number of haplotypes (22), 17 of which were exclusive to this locality, and 14 were singletons. However, it is important to note here that the samples from this site were obtained from a fish market which receives catches from a wide area surrounding the Amazon estuary, from the state of Amapá in the north, to Maranhão in the east (Braga CF Msc Thesis, Universidade Federal do Pará, Belém-PA, 2002). Because of its legally protected status, fishermen are reluctant to provide information on their fishing sites of this species, but it seems likely that the samples collected in this study were originated from a relatively wide geographic area. This also supports the idea that the Amazon region is an important refuge for E. itajara, given its extensive areas of well-preserved mangrove habitat, and the high productivity of local fisheries in Pará state (Braga CF Msc Thesis, Universidade Federal do Pará, Belém-PA, 2002; Tak-Chuen and Ferrera, 2006).
Haplotypic diversity varied from 0.53 in Parnaíba to 0.86 in Natal (Table 2). The values recorded here for the Goliath grouper are generally lower than those recorded in studies of other grouper species - Epinephelus marginatus (Maggio et al., 2006), and Epinephelus labriformis and Epinephelus clippertonensis (Craig et al., 2006) - however, based on a marker (the mitochondrial cyt b gene) which evolves more slowly than the control region. Nucleotide diversity is low (0.1% to 0.5%) in all E. itajara populations, which may reflect anthropogenic pressures - not just harvesting, but also the loss of mangrove habitat - or other evolutionary process, and will be an extremely important information to plan the management and conservation of the species. Conversely, the high diversity found in Natal could also be due to gene flow from populations from further south, which can be confirmed, in the future, by enlarging the sampling area.
The pairwise Fst values were significant when comparing Parnaíba with Ajuruteua (0.24971, p < 0.01) and Natal (0.28711, p < 0.01), while between Ajuruteua and Natal low and not significant value was observed (0.04891, p > 0.05). These three populations were used in AMOVA, which indicated that 19% of the observed variation was the result of variation among populations. Despite the significant Φst value observed in this analysis (0.185, p < 0.01), no significant structuring was found in the three populations (data not shown).
All populations except Parnaíba showed negative values of Fs (Table 2), although the values were only significant for the sample from Bragança, and the general population, which is probably due to the excess of uncommon haplotypes, a characteristic of populations in expansion (Fu, 1997). The unimodal curve of the mismatch distributions of these populations (data not shown), the Raggedness (0.03 in both cases) and SSD (p > 0.05) values, and the different values of θ are all typical of populations in expansion. The τ values recorded for Bragança (3.55) and the general population (3.11) indicate that they have been expanding over the past 197 thousand and 173 thousand years, respectively. The estimates of the effective population size for the general population were approximately 27,300 and 17,660 using the Watterson (1975) and Tajima (1983) estimates of θ ± SD, respectively θS ± SD = 3.192 ± 1.059 and θΠ ± SD = 2.062 ± 1.289.
The nested clade analysis (Figure 2) only revealed a significant association between genetic variability and geographic distribution. The category at the 1-1 nested clade level showed significant differences (p < 0.001) for clade (Dc) and nested clade (Dn) distances, which indicates restricted gene flow with isolation by distance. However, the Mantel test found no association between the Fst values and geographic distances (p > 0.05). These results may be explained by the small number of haplotypes (two) recorded in the population from Parnaíba, and the fact that the one of these haplotypes is probably the ancestral form. While haplotypic diversity is moderate in this case (h = 0.53), it is important to remember that where a population is declining, the number of alleles will be affected before any changes in genetic diversity (Frankham et al., 2002). These results indicate that E. itajara populations may suffer greater pressures in northeastern Brazil (Fortaleza, Parnaíba and Natal), where both fishery activity and habitat loss are more intense. However, a much wider range of populations must be studied throughout the geographic distribution of the species, as well as the use of nuclear markers, such as microsatellites, before it will be possible to analyse recent gene flow and population structure and dynamics more reliably. The sub-structuring of populations observed in the present study may be favoured by life history traits of the species, such as its territoriality, habitat preferences, reproductive shoaling and fidelity to shoaling sites (Heemstra and Randall, 1993).
The present study provides the first data on the genetic variability and population structure of the Goliath grouper, Epinephelus itajara, a critically endangered New World fish species considered to be an important fishery resource. The analyses of populations from the northern Brazilian coast indicated overfishing, especially in the Northeast. It is important to remember that a reduction in the effective size of populations of species such as E. itajara may leave the species extremely vulnerable to evolutionary processes. While samples of a larger number of individuals and localities will be necessary for a more reliable characterisation of the population dynamics of this species, the levels of genetic variability and structuring of stocks recorded here indicate an urgent need for the implementation of conservation and management measures in order to protect and consolidate these populations.
Many different fish species present clear signs of decline in abundance, which includes many different sharks (Baum et al., 2003). Therefore, more effective controls of fishery activity are essential here, but the preservation of mangrove habitats and shoaling areas will be equally important for the long-term conservation of the species.
Figure 1 Location of the collecting localities for the Goliath grouper on the Northern Brazilian coast (base map was created at Online Map Creations, www.aquarius.geomar.de/omc).
Figure 2 Minimum-spanning network of Epinephelus itajara haplotypes and the nested clade design on which the nested clade analysis was based. The two specimens from Fortaleza (H7) were added to this analysis.
This work was supported by CNPq (Programa Taxonomia n. 563967/2005-6) for a grant to G. Silva-Oliveira. We also thank Stephen Ferrari for revision of the manuscript.
Baum J.K., Myers R.A., Kehler D.G., Worm B., Harley S.J., Doherty P.A. (2003) Collapse and conservation of shark populations in the northwest Atlantic. Science 299:389-392.
Bullock L.H., Murphy M.D., Godcharles M.F., Mitchell M.E. (1992) Age, growth, and reproduction of jewfish Epinephelus itjara in the eastern Gulf of Mexico. Fish Bull (Wash D C) 90:243-249.
Clement M., Posada D., Crandall K.A. (2000) TCS: A computer program to estimate gene genealogies. Mol Ecol 9:1657-1659.
Craig M.T., Hastings P.A., Pondela D.J., Robertson D.R., Rosales-Casian J.A. (2006) Phylogeography of the flag cabrilla Epinephelus labriformis (Serranidae): Implications for the biogeography of the Tropical Eastern Pacific and the early stages of speciation in a marine shore fish. J Biogeogr 33:969-979.
Donaldson K.A., Wilson R.R.J. (1999) Amphi-panamic geminates of snook (Percoidei, Centropomidae) provide a calibration of the divergence rate in the mitochondrial DNA control region of fishes. Mol Phylogenet Evol 13:208-213.
Excoffier L., Laval G., Schneider S. (2005) Arlequin v. 3.0: An integrated software package for population genetics data analysis. Evol Bioinform Online 1:47-50.
Excoffier L., Smouse P.E., Quattro J.M. (1992) Analysis of molecular variance inferred from metric distances among DNA haplotypes: Application to human mitochondrial DNA restriction data. Genetics 131:479-491.
Francesconi J.J., Schwartz F.J. (2000) Jewfish, Epinephelus itajara, from North Carolina, with range correction and body comparisons. J Elisha Mitchell Sci Soc 116:167-170.
Frankham R., Balou J.D., Briscoe D.A. (2002) Introduction to Conservation GeneticsNew YorkCambridge University Press
Frias-Torres S. (2006) Habitat use of juvenile Goliath grouper Epinephelus itajara in the Florida Keys, USA. Endang Species Res 2:1-6.
Fu Y.X. (1997) Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics 147:915-925.
Hall T.A. (1999) BioEdit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser 41:95-98.
Heemstra P.C., Randall J.E. (1993) FAO species catalogue: Groupers of the world (Family Serranidae, Subfamily Epinephelinae). An annotated and illustrated catalogue of the grouper. Rockcod, Hind, Coral grouper and Lyretail species known to date. FAO Fish Synops n. 125.RomeFAO
Kumar S., Tamura K., Nei M. (2004) Mega 3: Integrated software for molecular evolutionary genetics analysis and sequence alignment. Bioinformatics 5:150-163.
Lee W.J., Conroy J., Howell W.H., Kocher T.D. (1995) Structure and evolution of teleost mitochondrial control regions. J Mol Evol 41:54-66.
Maggio T., Andaloro F., Arculeo M. (2006) Genetic population structure of Epinephelus marginatus (Pisces, Serranidae) revealed by two molecular markers. Ital J Zool 73:275-283.
Morris A.V., Roberts C.M., Hawkings J.P. (2000) The threatened status of groupers (Epinephelinae). Biodivers Conserv 9:919-942.
Nei M. (1987) Molecular Evolutionary GeneticsNew YorkColumbia University Press
Palumbi S.R., Martin A., Romano S., McMillan W.O., Stice L., Grabowski G. (1991) The Simple Fool's Guide to PCRHonoluluUniversity of Hawaii Press
Randall J.E. (1983) Caribbean Reef Fishes2nd editionNeptune CityTFH Publications
Rivera M.A.J., Kelley C.D., Roderick G.K. (2004) Subtle population genetic structure in the Hawaiian Grouper, Epinephelus quernus (Serranidae) as revealed by mitochondrial DNA analyses. Biol J Linn Soc 81:449-468.
Rogers A.R. (1995) Genetic evidence for a Pleistocene population explosion. Evolution 49:608-615.
Rogers A.R., Harpending H. (1992) Population growth makes waves in the distribution of pairwise genetic differences. Mol Biol Evol 9:552-569.
Rozas J., Sánchez-DelBarrio J.C., Messeguer X., Rozas R. (2003) DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics 19:2496-2497.
Sadovy Y., Eklund A.-.M. (1999) Synopsis of biological data on the Nassau grouper, Epinephelus striatus (Bloch, 1792), and the jewfish, E. itajara (Lichtenstein, 1822). NOAA Technical Report 146SeattleNational Marine Fisheries Service
Sambrook J., Fritsch E.F., Maniatis T. (1989) Molecular Cloning: A Laboratory Manual2nd editionNew YorkCould Spring Harbor Laboratory Press
Santa Brígida E.L., Cunha D.B., Rego P.S., Sampaio I., Schneider H., Vallinoto M. (2007) Population analysis of Scomberomorus cavalla (Cuvier, 1829) (Perciformes, Scombridae) from the Northern and Northeastern coast of Brazil. Braz J Biol 67:919-924.
Schneider S., Excoffier L. (1999) Estimation of past demographic parameters from the distribution of pairwise differences when the mutation rates vary among sites: Application to human mitochondrial DNA. Genetics 152:1079-1089.
Szpilman M. (2000) Peixes Marinhos do Brasil: Guia Prático de IdentificaçãoRio de JaneiroInstituto ecológico Aqualung e MAUAD Editora
Tajima F. (1983) Evolutionary relationship of DNA sequences in finite populations. Genetics 105:437-460.
Templeton A.R., Sing C.F. (1993) A cladistic analysis of phenotypic associations with haplotypes inferred from restriction endonuclease mapping. IV. Nested analyses with cladogram uncertainty and recombination. Genetics 134:659-669.
Templeton A.R., Boerwinkle E., Sing C.F. (1987) A cladistic analysis of phenotypic associations with haplotypes inferred from restriction endonuclease mapping. I. Basic theory and an analysis of alcohol dehydrogenase activity in Drosophila. Genetics 117:343-351.
Watterson G. (1975) On the number of segregating sites in genetical models without recombination. Theor Popul Biol 7:256-276.
Zatcoff M.S., Ball A.O., Sedberry G.R. (2004) Population genetic analysis of red grouper, Epinephelus morio, and scamp, Mycteroperca phenax, from the southeastern U. S. Atlantic and Gulf of Mexico. Mar Biol 144:769-777.
Coleman and Koening Laboratory Research (2008) Available from: http://www.bio.fsu.edu/coleman_lab/goliath_grouper.php (July30,2008)
Tak-Chuen TC and Ferrera BPEpinephelus itajara. In: IUCN 2007. 2007 IUCN Red List of Threatened Species. (2006) Available from: http://www.iucnredlist.org (August2,2008)
- Haplotypes of the control region of the mtDNA of the five populations of Epinephelus itajara analysed in the present study.
1AJU = Ajuruteua; BRA = Bragança; FOR = Fortaleza; NAT = Natal; PAR = Parnaíba.
- Indexes of genetic variability and results of the neutrality test for the different populations of E. itajara analysed in the present study.
*Significant value (p < 0,05).
Received: June 13, 2008; Accepted: August 19, 2008
Marcelo Vallinoto. Laboratório de Genética e Biologia Molecular, Instituto de Estudos Costeiros, Universidade Federal do Pará, Al. Leandro Ribeiro s/n, 68.600-000 Bragança, PA, Brazil. E-mail: firstname.lastname@example.org.
Publication in this collection
31 Oct 2008
Date of issue