- Research article
- Open Access
Contrasting population genetic structure in three aggregating groupers (Percoidei: Epinephelidae) in the Indo-West Pacific: the importance of reproductive mode
BMC Evolutionary Biology volume 18, Article number: 180 (2018)
Understanding the factors shaping population genetic structure is important for evolutionary considerations as well as for management and conservation. While studies have revealed the importance of palaeogeographic changes in shaping phylogeographic patterns in multiple marine fauna, the role of reproductive behaviour is rarely considered in reef fishes. We investigated the population genetics of three commercially important aggregating grouper species in the Indo-West Pacific, namely the camouflage grouper Epinephelus polyphekadion, the squaretail coral grouper Plectropomus areolatus, and the common coral trout P. leopardus, with similar life histories but distinct spatio-temporal characteristics in their patterns of forming spawning aggregations.
By examining their mitochondrial control region and 9–11 microsatellite markers, we found an overarching influence of palaeogeographic events in the population structure of all species, with genetic breaks largely coinciding with major biogeographic barriers. The divergence time of major lineages in these species coincide with the Pleistocene glaciations. Higher connectivity is evident in E. polyphekadion and P. areolatus that assemble in larger numbers at fewer spawning aggregations and in distinctive offshore locations than in P. leopardus which has multiple small, shelf platform aggregations.
While palaeogeographic events played an important role in shaping the population structure of the target species, the disparity in population connectivity detected may be partly attributable to differences in their reproductive behaviour, highlighting the need for more investigations on this characteristic and the need to consider reproductive mode in studies of connectivity and population genetics.
A primary objective of phylogeography is to identify the historical events that have shaped the current population genetic structure within species . In the Indo-Pacific, the Pleistocene glaciations have left pronounced evolutionary footprints on a wide range of marine fauna. These glacial cycles resulted in significant reduction of coastal marine habitats and the emergence of land barriers, such as the Indo-Pacific Barrier (IPB) between the Indian and Pacific Oceans , and the Red Sea Barrier (RSB) at the shallow (137 m) and narrow (18 km) strait at Bab al Mandab, the sole connection between the Red Sea and the Indian Ocean . The associated habitat loss [4, 5] and population isolation  led to dramatic population bottlenecks , which might have facilitated genetic differentiation among allopatric populations as genetic drift operates more efficiently in smaller populations. Moreover, environmental differences such as salinity or temperature may generate selective pressures that fast track divergence between isolated populations in different environments . Nonetheless, although some reef fishes exhibit strong genetic partitioning at the IPB and RSB with very limited genetic exchange, such as damselfish , parrotfish , angelfish and squirrelfish , some species, particularly those with high dispersal potential, are genetically homogeneous throughout the Indo-Pacific, such as unicornfishes [11, 12], and moray eels . Hence, while phylogeographic patterns of organisms are often interpreted through past environmental disturbances, mediated by climate changes and geographic barriers, they may also be strongly influenced by species-specific traits . The latter has not received much consideration.
The evolution of groupers (family Epinephelidae), a commercially important and speciose taxon of considerable biomass in reef ecosystems, was significantly impacted by environmental change associated with glacial cycles, as revealed by a recent historical biogeographic study . These changes resulted in increased allopatric divergences across major biogeographic barriers (e.g. IPB and RSB) during the Pliocene and Pleistocene when sea levels dropped (5.3–0.01 million years ago (MYA)). However, little is known about the effect of Pleistocene climate change on intraspecific diversification of groupers, particularly whether or not sea level changes over the last four glacial periods has led to finer-scale population differentiation within biogeographic regions, as has been documented for other species of coral reef fishes such as the ocellaris clownfish  and the common coral trout . Moreover, there are limited studies on whether co-distributed species that differ only in one or a few major trait(s) exhibit marked differences in their population genetic structure.
Here, we examined the population genetic structure of three common and commercially important groupers in the Indo-West Pacific: the camouflage grouper Epinephelus polyphekadion (Bleeker, 1849), the square-tail coral grouper Plectropomus areolatus (Rüppell, 1830), and the common coral trout, P. leopardus (Lacepède, 1802). These groupers are top predators on coral reefs that feed predominantly on fish  and hence have an important ecological role as part of the predator biomass that helps to shape the reef ecosystem. These piscivores are large (some exceeding 1 m in total length), roaming in shallow (usually < 50 m, but some occur at 300 m) tropical and subtropical waters throughout the Indo-Pacific and largely co-occur in the Pacific , although P. leopardus has a restricted distributional range relative to the other two species (Fig. 1). All of them are common enough throughout their distribution ranges to support fisheries . Their pelagic larval durations (PLDs) are thought to range from 3 to 6 weeks [20, 21] and home ranges in non-reproductive periods are generally less than 5 km2 . The three species all aggregate to spawn whereby gamete release takes place over a relatively short period (a few days) during a few to multiple spawning months. Individuals can display high spawning site fidelity, as determined from tagging studies [23, 24]. However, they differ substantially in the degree to which their spawning aggregation is spatially, numerically and temporally concentrated (see paragraphs below).
The formation of spawning aggregations is a reproductive strategy exhibited by many fish families, including coral reef species such as groupers, snappers, jacks (trevallies), surgeonfishes, damselfishes and parrotfishes [25,26,27,28,29]. Fish spawning aggregations (FSA) of groupers that exhibit this reproductive mode are typically brief, spatially restricted, gatherings of tens to tens of thousands of reproductive individuals drawn from a wider catchment area, which can include fish that have travelled tens to hundreds of kilometres to reach the aggregation site . This is in marked contrast to many reef fishes which spawn within or close to their home range over long periods of time each year [31, 32]. FSAs may couple with short-term oceanographic conditions such as slack currents which enhance fertilization success or they may coincide with rapid currents that could strongly influence either larval retention or dispersal [33,34,35]. Some aggregations may occur on the reef platform, others only at the shelf edge exposed to extreme oceanic conditions. Such extreme temporal and spatial concentrations (point sources and locations) of gamete release during annual spawning events may confer significant impacts on a species’ population connectivity, but this factor has rarely been considered.
Among our target species, P. leopardus forms multiple small (a few hundreds of individuals) aggregations for 1–2 months each year and can migrate up to 10 km between spawning sites and home ranges , while P. areolatus forms relatively fewer, yet much larger (up to thousands of individuals), aggregations for multiple months (up to 8 months each year, although fewer months in some locations), with spawning migrations of up to 30 km [22, 36]. By contrast, E. polyphekadion forms few, large (thousands of individuals), spawning aggregations for one or sometimes two months a year, with adults known to travel up to 40 km to spawning sites from their home reefs (K. Rhodes – personal communication). The latter two species aggregate at shelf edge/drop-off areas, often associated with outer reef channels and close to offshore waters, for a few days around the full or new moon (depending on location) , whereas P. leopardus characteristically spawns on the shelf platform on the full moon in Australia . The distinctive temporal and spatial differences in the spawning aggregations of the three species may have differentially impacted their phylogeographic patterns. Specifically, a more concentrated (i.e. in time and space) release of gametes and/or more frequent spawning at offshore sites, as evident in P. areolatus, may have evolved to favour long distance larval dispersal, which may have weakened the phylogeographic signal in this species.
This study aims to test whether the phylogeographic patterns of the three grouper species (1) are shaped by past environmental changes, and (2) whether disparity in population genetic structure among the species, if any, could be attributable to differences in their reproductive traits. In this regard, we would predict P. areolatus to exhibit the lowest level of genetic population structure across its geographic range because of the concentrated release of gametes and more frequent spawning that could have evolved to favour long distance dispersal. By contrast, we expected P. leopardus to have the strongest genetic population structure across its range, being a poor disperser relative to the other two species, due to smaller and spatially scattered aggregations that are not close to offshore waters. To test these hypotheses, we conducted an extensive genetic examination over the entire range of the three species by using the mitochondrial control region and nuclear microsatellite markers. Moreover, since our target species are among the most commercially important reef fish in the Indo-Pacific, and some are already threatened by overfishing, connectivity information would be valuable to improve their management and conservation. We evaluate our findings against comparable population genetics work in other groupers and consider possible management and conservation implications.
Sample collection and DNA extraction
Fin clips or muscle tissues from the three species were obtained from 28 locations throughout the Indo-Pacific (Fig. 1). Sampling was undertaken by selective spearing or line fishing, or were acquired from local fish markets. All tissues were preserved immediately in 90% ethanol or NetStar solution. Total DNA was extracted using a QIAamp Tissue Kit (QIAGEN) or 5% Chelex solution .
Mitochondrial DNA (mtDNA) analyses
A 580-base pair (bp) fragment of mtDNA control region was amplified using primers developed for this study (see Additional file 1: Table S1). Polymerase chain reaction (PCR) mixtures contained 1–5 μl of template DNA, 1X PCR reaction buffer, consisting of 3 mM MgCl2, 200 μM dNTPs, 200 nM of each primer, 1.5 U of Taq polymerase (Amersham) and ddH2O to a total volume of 50 μl. The thermal cycle consisted of an initial denaturation at 94 °C for 3 min followed by 38 cycles at 94 °C for 30 s, 54 °C for 50 s, 72 °C for 1 min 30 s, with a final extension at 72 °C for 10 min. PCR products were purified using Millipore Montage PCR96 Cleanup Kit, following manufacturer’s instructions. Sequencing was performed in the forward direction using an Applied Biosystems 3100 sequencer and followed standard cycle sequencing protocols (BGI, China). Sequence data were edited using MEGA 5.0  and aligned using MUSCLE  with default settings, as implemented in MEGA 5.0 and adjusted manually. The haplotype diversity (h) and nucleotide diversity (π)  were assessed using Arlequin 3.5 . Neutrality tests (Tajima’s D , and Fu’s FS ) were conducted in Arlequin with 1000 permutations.
Population genetic structure was examined by three means: (1) SAMOVA 2.0  was used to analyse the spatial genetic structure with 1000 permutations and the number of initial conditions was set to 100. The maximum proportion of missing data was set as 0.5%. If the number of individuals analysed at a site (e.g. Bali and Borneo) was less than 16, the site(s) would either be removed or grouped with nearby sites (as noted in the footnote to Additional file 1: Tables S2–4). We tested the number (K) of groups of populations from 2 to the number of populations minus 1; the K value with the highest and significant FCT was regarded as the best representation of the number of distinctive groups for that dataset. SAMOVA was run with and without geographic information, which gave consistent results. To reduce errors in significance calculations due to missing data, selected grouping schemes were subjected to locus-by-locus Analysis of Molecular Variance (AMOVA) in Arelquin 3.5; (2) Pairwise ΦST statistics were computed using Arlequin with 99,999 permutations; and (3) The evolutionary relationships among haplotypes and genetic structure were evaluated and visualized by a minimum spanning network (MSN) analysed using Arlequin and drawn using Gelphi  and Adobe Illustrator.
Based on the average nucleotide substitution rate of the control region of 10% Myr− 1 between lineages , the divergence time among genetic clades was approximately estimated.
Between nine and eleven microsatellite loci were analysed in the three species examined (Additional file 1: Table S5). All PCR mixes contained 1.25–1.5 μl of template DNA, 1X KAPA2G Buffer A, 3 mM MgCl2, 150 μM dNTPs, 60–240 nM of each fluorescent-labelled (FAM, VIC, NED or PET) forward primer, 60–240 nM of each reverse primer, 4% Bovine Serum Albumin (BSA) and 1.5 U KAPA2G Fast DNA polymerase in a total volume of 6.25 μl. A touch-down thermal programme was used: initial denaturation at 95 °C for 2 min, followed by 30 cycles of denaturation at 95 °C for 10 s, 12 s at annealing temperature stepping down 0.2 °C per cycle from 62 °C to 56 °C and an extension at 68 °C for 3 s. Subsequently, another 10 cycles of denaturation at 95 °C for 10 s, 12 s annealing at 56 °C and extension at 68 °C for 3 s, were run. PCR products were genotyped using an Applied Biosystems 3100 DNA Analyzer along with GeneScan LIZ-500 (Applied Biosystems) as internal size standard. Allele sizes were analysed using GENEMARKER 2.4.0 (Softgenetics, LLC). MICROCHECKER 2.2.3  was used to test for null alleles, large allele dropout and scoring errors. We used FreeNA  to test for the influence of null alleles by estimating FST before and after ENA (excluding null alleles). Correction for null alleles and bootstrapping (10,000 pseudoreplicates) was used to determine significance levels. The observed and expected heterozygosities (Ho and He) were estimated for each locus using GenALEx 6.5 . Significant deviations from Hardy-Weinberg equilibrium (HWE) at each locus and linkage disequilibrium (LD) between all loci  were tested by running 10,000 Markov chain iterations in the online version of GENEPOP 4.0 .
Genotypic population structure was assessed using four approaches: (1) GenALEx was used to estimate pairwise F-statistics; (2) SAMOVA was used to analyse the spatial genetic structure as previously described for mtDNA. The maximum proportion of missing data was set at 5%; (3) A discriminant analysis of principal components (DAPC)  was employed to produce scatterplots of discriminant functions derived from the spatial distribution of microsatellite genotypes. We retained the number of principal components (PCs) after which little information was gained by adding PCs, and chose the optimal number of groups according to Bayesian information criterion (BIC); and (4) Population genetic structure was examined using a Bayesian approach implemented in STRUCTURE 2.3.3  under an admixture model with correlated allele frequencies and sampling locations as priors. Runs of 1 million steps and 20% burn-in for ten replicates of each value of K (from K = 1 to K = 12) were performed. STRUCTURE HARVESTER web 0.6.93  was used to determine the best K for the dataset using the method of Evanno, Regnaut & Goudet . We then divided the data by major groups to find hierarchical structure [57, 58].
This study analysed 723-bp, 580-bp and 556-bp segments of control region from 270 E. polyphekadion, 341 P. areolatus and 340 P. leopardus, respectively (Additional file 1: Tables S2–4, GenBank accession numbers KM656497-KM656787, MH853841-MH854500). Overall haplotype diversity (h) was high, ranging from 0.920 ± 0.01 in P. areolatus to 0.993 ± 0.002 in E. polyphekadion, while nucleotide diversity (π) ranged from 0.025 ± 0.013 in P. areolatus to 0.131 ± 0.063 in P. leopardus.
A total of 261, 397 and 365 individuals of E. polyphekadion, P. areolatus and P. leopardus, respectively, were genotyped for the microsatellite loci (Additional file 1: Tables S6–8). Large allele dropout and scoring errors were not detected in our datasets, but null alleles were present in 4 to 8 loci in 1 to 7 populations (Additional file 1: Table S5). Null allele frequency estimated by FreeNA ranged from 0 to 0.387 (mean 0.05) in E. polyphekadion, from 0 to 0.315 (mean 0.02) in P. areolatus, and from 0 to 0.317 (mean 0.076) in P. leopardus (data not shown). No significant disparity in global FST values computed with and without ENA correction was found (Additional file 1: Table S5). No pair of loci exhibited significant linkage disequilibrium in more than two populations (with α = 0.01, Additional file 1: Table S16). Therefore, we used all microsatellite loci in subsequent analyses. Locus-specific Ho and He ranged from 0.000 to 0.938 and from 0.000 to 0.909 in E. polyphekadion, from 0.000 to 1.000 and from 0.000 to 0.898 in P. areolatus, from 0.054 to 0.929 and from 0.053 to 0.913 in P. leopardus (Additional file 1: Tables S6–8). Significant deviation from HWE (Benjamini-Hochberg adjusted P < 0.01) was detected in 15/99, 8/126 and 28/108 cases E. polyphekadion, P. areolatus and P. leopardus, respectively, most of which were due to heterozygote deficiency (Additional file 1: Tables S6–8).
Population genetic structure
The MSN of the control region (Fig. 2a) revealed two distinct groups separated by 107 mutational steps: one in the Western Indo-Pacific (WIP); the other in the Central Indo-Pacific (CIP) combined with the Eastern Indo-Pacific (EIP) (following Spalding et al. ). The average genetic distance amongst these two groups was 14% (± 1.2%), and the molecular clock suggested a divergence time of 1.41 ± 0.12 MYA (Table 1). Fine-scale population genetic structure within groups was not evident in the MSN, with the exception that all the peripheral French Polynesian haplotypes were closely related to each other, unlike any other populations. Results from microsatellite DAPC analyses were consistent with the mtDNA data, with the Maldives partitioned substantially along the first principal component, and the others only partitioned along the second principal component, which accounted for about a third of the total genetic variation (Fig. 3a). A separate DAPC analysis of the CIP populations (i.e. excluding the Maldives data) revealed five clusters, one containing mostly the Cocos (Keeling) Islands (CK) samples, two containing most of the Western Australian, Philippines, Palau and Pohnpei samples, and the remaining two containing mostly the New Caledonia, Fiji and French Polynesia samples (Additional file 1: Figure S1). Consistent with this, STRUCTURE analysis of the microsatellite data identified a WIP and a CIP-EIP group (Fig. 4a). A further STRUCTURE assignment analysis within the latter group revealed finer scale population genetic structure as noted before, with almost all CK samples assigned to group 1, most samples from Western Australia, Palau, Philippines and Pohnpei assigned to group 2, and an increasing trend of assignments to group 3 from New Caledonia eastward to French Polynesia (Fig. 4a).
SAMOVA computed the highest FCT when K = 2 in both the control region (standard FCT = 0.852, P = 0.105, locus by locus FCT = 0.852, P < 0.0001) and in microsatellites (standard FCT = 0.265, P = 0.129, locus by locus FCT = 0.197, P < 0.0001, Additional file 1: Table S9), suggesting that E. polyphekadion is comprised of two distinctive groups: (1) WIP, including the Maldives and Phuket, and (2) the CIP-EIP, including samples from CK to French Polynesia (see Fig. 1A). SAMOVA within the latter group revealed considerable genetic partitioning between CIP and EIP in the control region (standard FCT = 0.219, P = 0.123, locus by locus FCT = 0.218, P < 0.0001) and in microsatellites (standard FCT = 0.066, P = 0.129, locus by locus FCT = 0.038, P < 0.05, Additional file 1: Table S9). SAMOVA of the CIP revealed very weak, but significant population structuring (control region standard FCT = 0.046, P < 0.01, locus by locus FCT = 0.046, P < 0.0001; microsatellite standard FCT = 0.051, P < 0.05, locus by locus FCT = 0.033, P < 0.01, Additional file 1: Table S9). F-statistics indicated the existence of a strong genetic break located west of CK (ΦST > 0.842 in mtDNA and FST > 0.127 in microsatellites, with all P < 0.001, Additional file 1: Tables S10, 11). Another lower magnitude genetic break was identified east of Fiji (ΦST > 0.290 and FST > 0.033) in mtDNA and microsatellites with all P < 0.001). Population differentiation within the CIP was generally weak (ΦST - 0.108 to 0.132, P < 0.00001) in four of 36 cases and FST - 0.035 to 0.092, P < 0.00001 in six of 15 cases), but CK was a genetically distinct population from all other populations based on microsatellites (FST > 0.036, all P < 0.001). We detected no sign of recent population expansion in the WIP group (Tajima’s D -1.173, P = 0.105; Fu’s FS -5.045, P = 0.038), while in the CIP-EIP group, Fu’s FS, which is more sensitive to sudden demographic expansion than Tajima’s D, was significantly negative (- 23.705, P < 0.01), while Tajima’s D was not significant (- 1.618, P = 0.017).
The MSN of the control region (Fig. 2b) partitioned P. areolatus into three groups: Red Sea (RS), Maldives-CK, and the Western Australia–Pacific Ocean (WAPO), with the WAPO population genetically less distinct from the Maldives-CK population (separated by eight mutational steps) than from the RS clade (separated by 28 mutational steps). The average genetic distance amongst the WAPO and Maldives-CK group was 1.6% (± 0.6%), and amongst the RS and the rest was 6.2% (± 1.2%) and the molecular clock estimates of divergence time were about 0.16 ± 0.06 MYA and 0.62 ± 0.12 MYA, respectively (Table 1). Fine-scale population genetic structure within both the RS and the WAPO groups was not evident. Indeed, a dominant haplotype was shared in each geographic clade, by all samples. In the Maldives-CK group, however, the two sites mostly harboured private haplotypes and shared a single haplotype (Fig. 2b). DAPC analyses of the microsatellites suggested 10 clusters, but the scatter plot depicted significant dissimilarity only among three cluster-groups corresponding to the RS, Central Indian Ocean Islands (CIOI sensu Spalding et al. 2007, i.e. Maldives), and CIP samples, with the CK cluster showing some distinction from the remaining CIP samples along the second principal component only, accounting for ~ 25% of the genetic variation captured by the first principal component (Fig. 3b). DAPC analysis for the RS populations found no fine-scale population structuring among sites (Additional file 1: Figure S2). By contrast, DAPC analysis within Western and Central Indo-Pacific populations found distinction among the Maldives, CK and the WAPO populations, while a further DAPC analysis of the WAPO populations revealed subtle differences between the Western Australian populations and the others (Additional file 1: Figure S4). Similarly, the best-fitting K was 2 using STRUCTURE for the entire microsatellite dataset, partitioning it into a RS and a CIOI-CIP group (Fig. 4b), but an additional STRUCTURE assignment analysis of the CIOI-CIP group alone identified two groups corresponding to Maldives-CK and WAPO populations. No further partitioning was detected in STRUCTURE analyses of the RS and WAPO datasets.
The optimal grouping inferred by SAMOVA with the highest FCT was congruent amongst both datasets in partitioning P. areolatus into three clusters, largely corresponding to the aforementioned pattern (control region standard and locus by locus FCT = 0.898; microsatellite standard FCT = 0.292, locus by locus FCT = 0.290), all with P < 0.0001 (Additional file 1: Table S9). Similar to E. polypekadion, only very weak population structure was revealed within WAPO (control region standard and locus by locus FCT = 0.045, P < 0.05; microsatellite standard and locus by locus FCT = 0.043, P < 0.05, Additional file 1: Table S9). However, CK changed group membership between the two datasets. Specifically, CK grouped with either the Maldives in CIOI (control region, ΦST > 0.750, P < 0.001) or the WAPO region (microsatellites, FST > 0.228, P < 0.001). Both datasets identified the genetic break between the RS and Maldives (Additional file 1: Table S9). Similarly, F-statistics (Additional file 1: Tables S12, 13) revealed that the RS population was strongly isolated (control region pairwise ΦST > 0.926, microsatellite pairwise FST > 0.184, all P < 0.001) and that the CK population was genetically distinct from all others (control region pairwise ΦST > 0.161, P < 0.001, microsatellite pairwise FST > 0.049, all P < 0.001). Population differentiation within the WAPO was generally low and insignificant (ΦST - 0.922 to 0.146), with significant ΦST (P < 0.00001) only between Fiji and Rowley Shoals, based on the control region (Additional file 1: Table S12). However, microsatellite data suggested significant, albeit subtle differentiation (FST 0.002 to 0.075), specifically between West Australian sites and the rest, with P < 0.00001 in 22 of 36 comparisons (Additional file 1: Table S13). Likewise, within the RS population, there was no significant differentiation between locations for either dataset (Additional file 1: Tables S12, 13). All three groups showed signs of recent demographic expansion with their significantly negative Tajima’s D and Fu’s FS: WAPO Tajima’s D -2.033, P < 0.01; Fu’s FS -26.984, P < 0.01; Maldives-CK: Tajima’s D -1.507, P < 0.05; Fu’s FS -7.311, P < 0.01; RS: Tajima’s D -2.622, P < 0.01, Fu’s FS -27.606, P < 0.01).
The control region MSN (Fig. 2c) revealed two major groups that were separated by 101 mutational steps partitioned along Sahul Shelf. We herein termed the meta-population east of this break (from Torres Strait to Fiji) the Eastern population, and the meta-population west of it (from Western Australia to East Indonesia) the Western population. The average genetic distance amongst the two groups was 16.9% (± 1.8%), and the molecular clock suggested a divergence time of about 1.69 ± 0.18 MYA (Table 1). The Western population consisted of two groups separated by 15 mutational steps. One group contained all the Western Australian and Palauan P. leopardus, as well as some from East Indonesia and a Taiwanese individual. The other group was divided into two subgroups, ten mutational steps apart. One subgroup harboured mostly individuals from Taiwan and Cebu (central Philippines), while individuals from Palawan (west Philippines) and Indonesia dominated the other subgroup. The Eastern population contained three major groups: one group consisted mostly of Solomon Islands individuals, one harboured exclusively Fijian P. leopardus, and the remaining group contained individuals from the Torres Strait, Great Barrier Reef and New Caledonia. The DAPC analysis of microsatellites (Fig. 3c) also revealed an East-West division in this species, but the subpopulations were more mixed than revealed by the control region MSN. A further DAPC analysis within the Western population revealed four clusters, with one cluster containing almost exclusively all the Palauan samples, one containing predominantly Taiwanese samples, one containing mostly the Western Australian samples, and one containing mostly Philippine samples (Additional file 1: Figure S5). The DAPC analysis of the Eastern population identified six clusters, and the scatter plot showed strong genetic division between Fiji and the remaining populations (Additional file 1: Figure S6). Population assignment by STRUCTURE also identified two distinct P. leopardus clusters - an Eastern and a Western group, each of which could be further divided into three subgroups (Fig. 4c). The Western population contained three subgroups corresponding to Western Australian, Palauan, and Philippine-Taiwanese populations. The Eastern population contained three subgroups corresponding to a Great Barrier Reef, a New Caledonian and a Fijian population.
The K value for P. leopardus estimated by SAMOVA for the control region was nine (standard FCT = 0.879, locus by locus FCT = 0.795, all P < 0.0001) and almost all sites were partitioned, except Torres Strait grouped with the Capricorn Bunker Reefs, Cebu grouped with Taiwan, and Scott Reef grouped with Palau (Additional file 1: Table S9). In the microsatellite dataset, K was ten (standard FCT = 0.250, P = 0.105, locus by locus FCT = 0.133, P < 0.0001), with almost all sites being partitioned as before, except that three northern Great Barrier Reef sites (Lizard Island, Hick’s Reef, and Townsville) formed a single group (Additional file 1: Table S9). The strong population genetic structure was also supported by the fact that > 80% of control region pairwise ΦST and > 96% of microsatellite pairwise FST values were high (averaged 0.696 for ΦST and 0.116 for FST) and were statistically significant (P < 0.001) (Additional file 1: Tables S14, 15). The F-statistics revealed a very strong genetic break between East Indonesia and Torres Strait (control region) and between Palau and the Great Barrier Reef (microsatellites), as pairwise F-statistics between sites on opposite side of this break were higher (mean ΦST = 0.935; mean FST = 0.102) than those on the same side (mean ΦST = 0.463; mean FST = 0.057). Both Tajima’s D (- 0.228, P = 0.466) and Fu’s FS (- 10.640, P = 0.030) were negative but not significant in the Eastern group, while Fu’s FS was significantly negative in the Western group (- 23.690, P < 0.01), Tajima’s D was not (- 0.782, P = 0.234).
Influence of past environmental changes and geographic factors
This study identifies phylogeographic signatures of Pleistocene glacial cycles in all three species, with genetic breaks generally coinciding with major biogeographic boundaries - the most prominent one being the IPB (see Fig. 1). As shallow reef inhabitants, groupers could have experienced marked population bottlenecks due to coastal habitat loss associated with marine regression during Pleistocene glaciation [4, 5, 60], followed by recent demographic expansion, as suggested by the significantly negative Fu’s FS in many of the grouper lineages from this study. The emergence of multiple land barriers including the IPB at the Sahul and Sunda Shelves  was one of the most significant barriers that impacted multiple fauna and could have physically segregated the grouper populations in CIOI and Indo-Polynesia (IP) , leading to genetic differentiation and/or speciation amongst these biogeographic compartments. Such patterns have frequently been observed in marine fauna, such as parrotfishes [9, 62], gobies , angelfish , gastropods , prawns  and multiple reef fishes [10, 67]. The notion that Pleistocene glaciation caused the genetic differentiation across IPB in groupers is supported in E. polyphekadion and P. areolatus based on the mid-Pleistocene divergence estimated by conventional molecular clock methods (1.41–0.16 MYA). Moreover, the IPB appears to present a dispersal boundary for P. leopardus such that the barrier might have resulted in the divergence of P. leopardus from its sister species P. pessuliferus (the roving coral grouper), which largely inhabited the west to central Indian Ocean, in the Pleistocene .
At the western periphery of the Indo-Pacific, the RSB exerted a profound effect on P. areolatus, causing the divergence of the RS lineage at about 0.63 MYA (mid-Pleistocene) based on the conventional molecular clock. Previous comparative phylogeographic work that focused on examining this barrier found genetic differentiation between the RS and Western Indian Ocean populations in five out of seven reef fish species that exhibit a wide spectrum of biological traits . Similar genetic partitioning across the RSB was also detected in damselfish , crab , and starfish , but not in lionfish . The RS was repeatedly isolated during the Pleistocene glacial cycle owing to the marked sea level drop – either through emergence of a land barrier or reduction of water flow associated with increased salinity and temperature [72, 73]. Such physical isolation and the strong selection pressure associated with the divergent environmental conditions at the peripheral habitats of a species’ distribution range are likely to have contributed to the genetic diversification observed in P. areolatus (this study) and in other fauna .
Intriguingly, despite the marine transgression that allowed for secondary contact of previously isolated populations, as evidenced in the peacock grouper (Cephalopholis argus ), we observed no contemporary connectivity across the IPB in either E. polyphekadion or P. areolatus. One likely explanation is that the extensive separation (> 2000 km) of suitable coral habitats and the relative absence of stepping stones in the deep ocean (> 200 m) for dispersal between the WIOI and IP, and between RS and WIOI (see Fig. 1) present a challenge to dispersal of these groupers even without land barriers . Habitat discontinuity can result in significant population differentiation over short distances, even for marine fish with a pelagic larval phase, such as the corkwing wrasse  and triplefin blenny . A combination of ecological, biological and environmental factors, rather than oceanography per se, may be more important in determining population genetic structure across major biogeographic barriers for the grouper species examined here.
At the eastern boundary of the species’ range (Tuamotus of French Polynesia), we detected subtle but significant isolation of the French Polynesian E. polyphekadion, which has lower genetic diversity (in both mitochondrial and microsatellite loci) than other populations. Multiple factors, including geographic isolation, oceanographic currents that influence gene flow, and selection pressure in distinct environmental conditions, may play a role in creating the genetic differentiation in this marginal population of E. polyphekadion [79, 80]. In fact, such patterns of extensive population connectivity across the wide range of the WIP with isolation at the eastern margin of the Pacific range is frequently detected in reef fishes, such as the peacock grouper , the brown surgeonfish , the bluestriped snapper , and the flame angelfish .
Our results revealed multiple genetically differentiated populations of P. leopardus in the WIP, with the deepest divergence at the Sahul Shelf (Fig. 1c). While the other populations were likely isolated by regional oceanographic features , the divergence of the East and West populations of P. leopardus was estimated at about 1.69 MYA (early Pleistocene), suggesting the possible contribution of the Pleistocene glaciation that resulted in the recurrent emergence of the Sahul shelf land bridge. The Sahul shelf land bridge has been documented as a barrier associated with genetic differentiation for several other tropical marine species, from starfish  to sea snakes . Though the land bridge has been submerged for the past 9000 years, the region remained uninhabited by P. leopardus (see Fig. 1), and thus, given the lower dispersal potential of this species (see discussion below), the two meta-populations may maintain low connectivity and genetic differentiation. Another possible explanation for the pattern observed is that each isolated region has developed region-specific traits, though this notion remains speculative until the ecology of P. leopardus across its range is studied in more detail. Theoretically, region-specific traits may prevent realized dispersal during transgressions by selection against immigrants from other environments, as documented for Atlantic herring using genomic approaches . This was also demonstrated in the Dascyllus trimaculatus (Rüppell, 1829) threespot dascyllus species complex, where historical allopatry developed during periods of isolation and subsequent ecological factors facilitated divergence and speciation of these species .
Notably, the affinity of the CK population, which is located at the WIP-CIP boundary, is different between E. polyphekadion and P. areolatus. In E. polyphekadion, the CK population groups with the IP population in all analyses (mitochondrial and microsatellite data). However, in P. areolatus, the CK population has a close affinity with the Maldives population in the WIOI, based on the mitochondrial control region, but it is more closely related, albeit also genetically distinct from, the IP populations based on FST, SAMOVA and DAPC analyses of microsatellite data. However, STRUCTURE results were congruent with mitochondrial data. While we cannot fully account for the phylogeographic pattern observed, the genetic differentiation inferred from microsatellite analyses of both species essentially reflects the relative geographic distances between the Maldives, CK and Western Australia, and, more importantly, implies that the CK population of both species was in effect isolated from the majority of CIP and WIP stocks at an ecological time scale.
Influence of species-specific features
A number of features could have, perhaps in concert, resulted in the drastically different phylogeographic patterns observed in the three grouper species. A longer evolutionary history could have allowed population differentiation to build up and intensify the genetic signal. The divergence of the P. areolatus WAPO lineage from the rest was much more recent (late Pleistocene) than divergence of the P. leopardus East and West lineages (early Pleistocene), which might partly explain the weak population structure observed in the former species. Nonetheless, length of evolutionary history alone cannot explain the similarly weak population structure observed in E. polyphekadion, whose WIP lineage also emerged in the early Pleistocene. The much lower dispersal ability of P. leopardus in comparison to the other species is the most probable explanation for the disparity observed. In fact, the smaller distribution range of P. leopardus might be an indication that this species is less capable of long distance dispersal.
The patterns of reproductive output per spawning event at each FSA differ substantially among the three species studied, leading to possible differences in connectivity linked to larval dispersal. Spawning aggregations of E. polyphekadion and P. areolatus in unexploited or lightly exploited areas can involve thousands of reproductive adults at fewer sites, with spawning occurring in reef passages adjacent to the open ocean around the full or new moon periods. In contrast, P. leopardus FSAs usually consist of a few hundred adults  at multiple sites, spawning down current during medium to strong current flows associated with the new moon (at least in Australia) . Hence, massive concentrated pulses of gametes are generated for the former two (reef passage-spawning) species, while P. leopardus aggregations generate an order of magnitude smaller and much more spatially distributed series of gamete source pulses within an area. Given that the three species have similar PLDs, and that monthly spawning events for all three species are similarly brief, often during times when tidal flows are greatest (full and/or new moon phases), one might expect the immense number of gametes released by E. polyphekadion and P. areolatus FSAs to have a greater probability of long-distance dispersal due to oceanography compared to the lower number of propagules released in dispersed smaller pulses (smaller groups of adults). Moreover, E. polyphekadion and P. areolatus aggregate next to the open ocean where currents can potentially carry gametes long distances from the aggregation area, which contrasts with the inshore/reef platform spawning locations of P. leopardus. The dispersal of a few concentrated larger cohorts could potentially be more strongly influenced by meso-scale oceanographic events than multiple, scattered, smaller cohorts, which would probably be more affected by short-term localized events.
On the other hand, given that P. areolatus spawn over multiple months while E. polyphekadion and P. leopardus only spawn for one or two months each year, we might predict wider overall dispersal in P. areolatus than in the other two species if oceanographic conditions vary across seasons and thereby spread offspring more widely. More/less extensive larval dispersal may be linked to weaker/stronger population genetic structure, respectively as has been documented in many terrestrial and aquatic (including marine) species [84, 85]. However, dispersal ability more broadly, like PLD per se, is not a good predictor of population genetic structure  or range size  in the marine realm. Maybe the sheer number and concentration of gametes released into strong currents, and possibly the timing and location of these releases, contribute to the increased range of dispersal between populations and the distributional range size of both E. polyphekadion and P. areolatus.
Few other studies have examined population connectivity in groupers and considered the possible effect of reproductive mode with similar outcomes to our study. Portnoy et al. tested the hypothesis that aggregative spawning behaviour affects gene flow by examining population genetic structure and connectivity of two groupers: the aggregate spawning Epinephelus guttatus and the non-aggregating Cephalopholis fulva . Over a relatively small geographical scale in the Caribbean Sea, they found weaker population genetic structure in the aggregate-spawner, contrary to the expectation of the authors that such a reproductive mode would lower the genetic exchange between catchment areas of spawning aggregations. This result was partly attributed by the authors to the fact that C. fulva is sedentary, spawning over large areas of the shallow reef platform, and may maintain the same territory, forming multiple small mating groups over multiple months in the year . Additionally, Portnoy et al. proposed that E. guttatus, which aggregates over just a few months each year and at the shelf drop off close to open ocean, may join other aggregations in subsequent (unexamined) years, thereby homogenizing the gene pool, given that E. guttatus can migrate more than 30 km to spawning sites. This proposal, however, does not consider pelagic larval distribution and recruitment as an alternative avenue of connectivity. Consistent with this finding, E. striatus, which forms large spatially and temporally restricted spawning aggregations like E. guttatus in the Caribbean, was also found to exhibit weak population genetic structure across a long distance (~ 1600 km) . Results from parentage analysis of this species suggested significant external recruitment caused by high population connectivity .
In summary, stronger genetic structuring appears to be closely associated with species that form multiple small spawning groups away from the shelf edge (e.g. P. leopardus and C. fulva). Reduced population structuring is evident among species that assemble in large numbers at relatively few spawning sites at the shelf edge close to the open ocean (e.g. E. polyphekadion, E. guttatus, E. striatus and P. areolatus) (Additional file 1: Figure S7). We have assumed that the population sizes (both historical and present) of the three species are comparable, and that the aggregation type for each species remains invariant throughout their distributional ranges. Given their similar habitat and dietary requirements, and considering that all are considered to be common enough to exploit throughout their ranges, the first assumption is reasonable. While the latter assumption has not been comprehensively assessed because data on aggregations on all three species are not available throughout their ranges, all available publications show that this assumption holds true (including where they co-distribute) and where comparisons have been possible (e.g. ; database from www.scrfa.org), fully supporting our approach. Nonetheless, further studies on reproductive biology and spawning behaviour of these species are encouraged to advance our understanding of factors shaping population connectivity. Such information could also provide important information for fishery management frameworks for these economically important reef fishes . Although other regionally relevant factors, such as retention of larvae at natal grounds in some locations, are known to occur (e.g. P. maculatus , and P. areolatus ), such studies of self- or local recruitment were performed on single cohorts of recruits from a limited part of the distributional range, using parentage analyses. Hence, these findings are unlikely to contradict the connectivity measured here over evolutionary time scales by sampling mixed age adult cohorts from across the range. While our study suggests an association between reproductive mode and population structure, to what extent this factor impacts connectivity at ecological time scales is largely unknown. Nonetheless, our results highlight the need to consider the possible influence of reproductive mode (spatial and temporal aspects of reproduction) in future studies of population genetics and connectivity.
Conservation and management implications
Groupers are highly sought after and many are also often intrinsically vulnerable to fishing pressure because of their longevity, late sexual maturation, high site fidelity, complex social structure and, for some species, protogynous hermaphroditism [95,96,97]. As many as 25% of all exploited epinepheline species may be at risk from fishing activities or coastal development . Species that form spawning aggregations heavily exploited by aggregation-fisheries may be particularly at risk because such reproductive gatherings of spawning fish tend to be temporally and spatially predictable and are easily and quickly fished out once discovered and if not suitably managed . There are declines or disappearances of E. polyphekadion and P. areolatus spawning aggregations in Palau, Fiji, Pohnpei, French Polynesia and Indonesia for example, (www.scrfa.org). Our study species are listed as threatened or near-threatened on the IUCN Red List, largely as a result of uncontrolled aggregation-fishing . Better fishery management and conservation (e.g. by marine protected areas (MPAs), seasonal protection from fishing during spawning, or, ideally, some combination) are urgently needed .
The results of this study highlight the importance of population structuring in determining possible management units and the spatial scale that management would need to consider. For example, species with large catchment areas like E. polyphekadion and P. areolatus would need larger MPAs and regional scale management than a species like P. leopardus. Conversely, other management methods may be more practical for the former two species such as seasonal protection during spawning complemented by seasonal sales bans because enforcement may be particularly challenging for offshore sites . However, connectivity analyses at smaller scales have also revealed high self- or local recruitment within a single cohort sometimes justifying the benefit of more local scale management strategies (e.g. P. maculatus , and P. areolatus, ). Taken together, both regional- and local-scale management may need to be considered to ensure sustainable fisheries of these economically important reef fishes.
This study identified phylogeographic signatures of Pleistocene glacial cycles in all three species, with genetic breaks generally coinciding with major biogeographic boundaries, while species-specific reproductive traits determine the magnitude of these phylogeographic signals. As hypothesized, P. leopardus exhibited markedly stronger population genetic structure at various geographic scales within the CIP than the other two groupers, while E. polyphekadion displayed only slightly stronger population structuring than P. areolatus. The populations in the CIP were generally highly connected in the latter two species, and a significant genetic break only occurred among ocean basins. Results from this study and a comparable study in the Caribbean both suggested that stronger genetic structuring appeared to be closely associated with species that form multiple small spawning groups on the shelf platform and away from the open ocean. Our results highlighted the need for more investigations on this characteristic and the need to consider reproductive mode in studies of connectivity and population genetics as well as in more sustainable fisheries management.
Central Indian Ocean Islands
Discriminant analysis of principal components
Fish spawning aggregation
Marine Protected Area
Minimum spanning network
Pelagic larval durations
Red Sea Barrier
Western Australia-Pacific Ocean
Avise JC. Phylogeography: the history and formation of species. Cambridge: Harvard University Press; 2000.
Randall JE. Zoogeography of shore fishes of the Indo-Pacific region. Zool Stud. 1998;37:227–68.
Klausewitz W. Evolutionary history and zoogeography of the Red Sea ichthyofauna. Fauna Saudi Arab. 1989;10:310–37.
Klicka J, Zink RM. The importance of recent ice ages in speciation: a failed paradigm. Science. 1997;277:1666–9.
Fauvelot C, Bernardi G, Planes S. Reductions in the mitochondrial DNA diversity of coral reef fish provide evidence of population bottlenecks resulting from Holocene sea-level change. Evolution. 2003;57:1571–83.
Briggs JC. Modes of speciation: marine Indo-West Pacific. Bull Mar Sci. 1999;65:645–56.
Lamichhaney S, Barrio AM, Rafati N, Sundström G, Rubin C-J, Gilbert ER, et al. Population-scale sequencing reveals genetic differentiation due to local adaptation in Atlantic herring. Proc Natl Acad Sci. 2012;109:19345–50.
Leray M, Beldade R, Holbrook SJ, Schmitt RJ, Planes S, Bernardi G. Allopatric divergence and speciation in coral reef fish: the three-spot dascyllus, Dascyllus trimaculatus, species complex. Evolution. 2010;64:1218–30.
Bay LK, Choat JH, van Herwerden L, Robertson DR. High genetic diversities and complex genetic structure in an Indo-Pacific tropical reef fish (Chlorurus sordidus): evidence of an unstable evolutionary past? Mar Biol. 2004;144:757–67.
DiBattista JD, Berumen ML, Gaither MR, Rocha LA, Eble JA, Choat JH, et al. After continents divide: comparative phylogeography of reef fishes from the Red Sea and Indian Ocean. J Biogeogr. 2013;40:1170–81.
Horne JB, van Herwerden L, Choat JH, Robertson DR. High population connectivity across the Indo-Pacific: congruent lack of phylogeographic structure in three reef fish congeners. Mol Phylogenet Evol. 2008;49:629–38.
Klanten OS, Choat JH, van Herwerden L. Extreme genetic diversity and temporal rather than spatial partitioning in a widely distributed coral reef fish. Mar Biol. 2007;150:659–70.
Reece JS, Bowen BW, Joshi K, Goz V, Larson A. Phylogeography of two moray eels indicates high dispersal throughout the Indo-Pacific. J Hered. 2010;101:391–402.
Papadopoulou A, Knowles LL. Toward a paradigm shift in comparative phylogeography driven by trait-based hypotheses. Proc Natl Acad Sci. 2016;113:8018–24.
Ma KY, Craig MT, Choat JH, van Herwerden L. The historical biogeography of groupers: clade diversification patterns and processes. Mol Phylogenet Evol. 2016;100:21–30.
Timm J, Figiel M, Kochzius M. Contrasting patterns in species boundaries and evolution of anemonefishes (Amphiprioninae, Pomacentridae) in the centre of marine biodiversity. Mol Phylogenet Evol. 2008;49:268–76.
van Herwerden L, Howard Choat J, Newman SJ, Leray M, Hillersøy G. Complex patterns of population structure and recruitment of Plectropomus leopardus (Pisces: Epinephelidae) in the Indo-West Pacific: implications for fisheries management. Mar Biol. 2009;156:1595–607.
Heemstra PC, Randall JE. Groupers of the world: FAO species catalogue 16. Rome: food and agriculture organization; 1993.
Craig MT, Sadovy de Mitcheson YJ, Heemstra PC. Groupers of the world. South Africa: NISC; 2011.
Doherty PJ, Fowler AJ, Samoilys MA, Harris DA. Monitoring the replenishment of coral trout (Pisces: Serranidae) populations. Bull Mar Sci. 1994;54:343–55.
Leis JM. Larval development in four species of Indo-Pacific coral trout Plectropomus (Pisces: Serranidae: Epinephelinae) with an analysis of the relationships of the genus. Bull Mar Sci. 1986;38:525–52.
Green AL, Maypa AP, Almany GR, Rhodes KL, Weeks R, Abesamis RA, et al. Larval dispersal and movement patterns of coral reef fishes, and implications for marine reserve network design. Biol Rev. 2015;90:1215–47.
Johannes RE, Squire L, Graham T, Sadovy Y, Renguul H. Spawning aggregations of groupers (Serranidae) in Palau. Nat Conserv Mar Res Ser Publ. 1999;1:1–144.
Zeller DC. Spawning aggregations: patterns of movement of the coral trout Plectropomus leopardus (Serranidae) as determined by ultrasonic telemetry. Mar Ecol Prog Ser. 1998;162:253–63.
Erisman BE, Buckhorn ML, Hastings PA. Spawning patterns in the leopard grouper, Mycteroperca rosacea, in comparison with other aggregating groupers. Mar Biol. 2007;151:1849–61.
Gerhardinger LC, Hostim-Silva M, Medeiros RP, Matarezi J, Bertoncini ÁA, Freitas MO, et al. Fishers’ resource mapping and goliath grouper Epinephelus itajara (Serranidae) conservation in Brazil. Neotrop Ichthyol. 2009;7:93–102.
Gladstone W. Temporal patterns of spawning and hatching in a spawning aggregation of the temperate reef fish Chromis hypsilepis (Pomacentridae). Mar Biol. 2007;151:1143–52.
Sadovy de Mitcheson Y, Cornish A, Domeier M, Colin PL, Russell M, Lindeman KC. A global baseline for spawning aggregations of reef fishes. Conserv Biol. 2008;22:1233–44.
Sala E, Aburto-Oropeza O, Paredes G, Thompson G. Spawning aggregations and reproductive behavior of reef fishes in the Gulf of California. Bull Mar Sci. 2003;72:103–21.
Bolden SK. Long-distance movement of a Nassau grouper (Epinephelus striatus) to a spawning aggregation in the Central Bahamas. Fish Bull Ocean Atmos Adm. 2000;98:642–5.
Nemeth RS, Blondeau J, Herzlieb S, Kadison E. Spatial and temporal patterns of movement and migration at spawning aggregations of red hind, Epinephelus guttatus, in the US Virgin Islands. Environ Biol Fish. 2007;78:365–81.
Rhodes KL, McIlwain J, Joseph E, Nemeth RS. Reproductive movement, residency and fisheries vulnerability of brown-marbled grouper, Epinephelus fuscoguttatus (Forsskål, 1775). Coral Reefs. 2012;31:443–53.
Cherubin LM, Nemeth RS, Idrisi N. Flow and transport characteristics at an Epinephelus guttatus (red hind grouper) spawning aggregation site in St. Thomas (US Virgin Islands). Ecol Model. 2011;222:3132–48.
Choat JH. Spawning aggregations in reef fishes; ecological and evolutionary processes. In: Sadovy de Mitcheson Y, Colin PL, editors. Reef fish spawning aggregations: biology, research and management. Dordrecht: Springer; 2012. p. 85–116.
Colin PL. Timing and location of aggregation and spawning in reef fishes. In: Sadovy de Mitcheson Y, Colin PL, editors. Reef fish spawning aggregations: biology, research and management. Dordrecht: Springer; 2012. p. 117–58.
Tupper M. Identification of nursery habitats for commercially valuable humphead wrasse Cheilinus undulatus and large groupers (Pisces: Serranidae) in Palau. Mar Ecol Prog Ser. 2007;332:189–99.
Rhodes KL, Sadovy Y. Temporal and spatial trends in spawning aggregations of camouflage grouper, Epinephelus Polyphekadion, in Pohnpei, Micronesia. Environ Biol Fishes. 2002;63:27–39.
Samoilys MA. Periodicity of spawning aggregations of coral trout Plectropomus leopardus (Pisces: Serranidae) on the northern Great Barrier Reef. Mar Ecol Prog Ser. 1997;160:149–59.
Walsh PS, Metzger DA, Higuchi R. Chelex 100 as a medium for simple extraction of DNA for PCR-based typing from forensic material. BioTechniques. 1991;10:506–13.
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S. MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011;28:2731–9.
Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32:1792–7.
Nei M. Molecular evolutionary genetics. New York: Columbia University Press; 1987.
Excoffier L, Lischer HEL. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour. 2010;10:564–7.
Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989;123:585–95.
Fu YX. Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics. 1997;147:915–25.
Dupanloup I, Schneider S, Excoffier L. A simulated annealing approach to define the genetic structure of populations. Mol Ecol. 2002;11:2571–81.
Bastian M, Heymann S, Jacomy M. Gephi: an open source software for exploring and manipulating networks. Icwsm. 2009;8:361–2.
Bowen BW, Muss A, Rocha LA, Grant WS. Shallow mtDNA coalescence in Atlantic pygmy angelfishes (genus Centropyge) indicates a recent invasion from the Indian Ocean. J Hered. 2006;97:1–12.
van Oosterhout C, Hutchinson WF, Wills DPM, Shipley P. MICRO-CHECKER: software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes. 2004;4:535–8.
Chapuis MP, Estoup A. Microsatellite null alleles and estimation of population differentiation. Mol Biol Evol. 2007;24:621–31.
Peakall R, Smouse PE. GenAlEx 6.5: genetic analysis in excel. Population genetic software for teaching and research--an update. Bioinformatics. 2012;28:2537–9.
Guo SW, Thompson EA. Performing the exact test of Hardy-Weinberg proportion for multiple alleles. Biometrics. 1992:361–72.
Rousset F. genepop’007: a complete re-implementation of the genepop software for Windows and Linux. Mol Ecol Resour. 2008;8:103–6.
Jombart T, Devillard S, Balloux F. Discriminant analysis of principal components: a new method for the analysis of genetically structured populations. BMC Genet. 2010;11:94.
Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics Genetics Soc America. 2000;155:945–59.
Earl DA, vonHoldt BM. STRUCTURE HARVESTER: a website and program for visualizing STRUCTURE output and implementing the Evanno method. Conserv Genet Resour. 2012;4:359–61.
Evanno G, Regnaut S, Goudet J. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol Ecol. 2005;14:2611–20.
Willis SC, Macrander J, Farias IP, Ortí G. Simultaneous delimitation of species and quantification of interspecific hybridization in Amazonian peacock cichlids (genus Cichla) using multi-locus data. BMC Evol Biol. 2012;12:96.
Spalding MD, Fox HE, Allen GR, Davidson N, Ferdaña ZA, Finlayson M, et al. Marine ecoregions of the world: a bioregionalization of coastal and shelf areas. Bioscience. 2007;57:573–83.
Dibattista JD, Howard Choat J, Gaither MR, Hobbs JPA, Lozano-Cortés DF, Myers RF, et al. On the origin of endemic species in the Red Sea. J Biogeogr. 2016;43:13–30.
Briggs JC. Marine zoogeography. New York: McGraw-Hill; 1974.
Fitzpatrick JM, Carlon DB, Lippe C, Robertson DR. The West Pacific diversity hotspot as a source or sink for new species? Population genetic insights from the Indo-Pacific parrotfish Scarus rubroviolaceus. Mol Ecol. 2011;20:219–34.
Hoareau TB, Boissin E, Berrebi P. Evolutionary history of a widespread Indo-Pacific goby: the role of Pleistocene sea-level changes on demographic contraction/expansion dynamics. Mol Phylogenet Evol. 2012;62:566–72.
Coleman RR, Eble JA, DiBattista JD, Rocha LA, Randall JE, Berumen ML, et al. Regal phylogeography: range-wide survey of the marine angelfish Pygoplites diacanthus reveals evolutionary partitions between the Red Sea, Indian Ocean, and Pacific Ocean. Mol Phylogenet Evol. 2016;100:243–53.
Crandall ED, Frey MA, Grosberg RK, Barber PH. Contrasting demographic history and phylogeographical patterns in two Indo-Pacific gastropods. Mol Ecol. 2008;17:611–26.
Tsoi KH, Ma KY, Wu TH, Fennessy ST, Chu KH, Chan TY. Verification of the cryptic species Penaeus pulchricaudatus in the commercially important kuruma shrimp P. japonicus (Decapoda:Penaeidae) using molecular taxonomy. Invertebr Syst. 2014;28:476–90.
Fernandez-Silva I, Randall JE, Coleman RR, DiBattista JD, Rocha LA, Reimer JD, et al. Yellow tails in the Red Sea: phylogeography of the Indo-Pacific goatfish Mulloidichthys flavolineatus reveals isolation in peripheral provinces and cryptic evolutionary lineages. J Biogeogr. 2015;42:2402–13.
Froukh T, Kochzius M. Species boundaries and evolutionary lineages in the blue green damselfishes Chromis viridis and Chromis atripectoralis (Pomacentridae). J Fish Biol. 2008;72:451–7.
Fratini S, Vannini M. Genetic differentiation in the mud crab Scylla serrata (Decapoda: Portunidae) within the Indian Ocean. J Exp Mar Bio Ecol. 2002;272:103–16.
Benzie JAH. Major genetic differences between crown-of-thorns starfish (Acanthaster planci) populations in the Indian and Pacific Oceans. Evolution. 1999;53:1782–95.
Kochzius M, Söller R, Khalaf MA, Blohm D. Molecular phylogeny of the lionfish genera Dendrochirus and Pterois (Scorpaenidae, Pteroinae) based on mitochondrial DNA sequences. Mol Phylogenet Evol. 2003;28:396–403.
Siddall M, Rohling EJ, Almogi-Labin A, Hemleben C, Meischner D, Schmelzer I, et al. Sea-level fluctuations during the last glacial cycle. Nature. 2003;423:853–8.
Bailey G. The Red Sea, coastal landscapes, and hominin dispersals. In: Petraglia M, Rose L, editors. The Evolution of Human Populations in Arabia. Dordrecht: Springer; 2009. pp. 15–37.
Suzuki G, Keshavmurthy S, Hayashibara T, Wallace CC, Shirayama Y, Chen CA, et al. Genetic evidence of peripheral isolation and low diversity in marginal populations of the Acropora hyacinthus complex. Coral Reefs. 2016;35:1419–32.
Gaither MR, Bowen BW, Bordenave T-R, Rocha LA, Newman SJ, Gomez JA, et al. Phylogeography of the reef fish Cephalopholis argus (Epinephelidae) indicates Pleistocene isolation across the Indo-Pacific barrier with contemporary overlap in the coral triangle. BMC Evol Biol. 2011;11:189.
Kemp J. Zoogeography of the coral reef fishes of the Socotra archipelago. J Biogeogr. 1998;25:919–33.
Gonzalez EB, Knutsen H, Jorde PE. Habitat discontinuities separate genetically divergent populations of a rocky shore marine fish. PLoS One. 2016;11:e0163052.
Riginos C, Nachman MW. Population subdivision in marine environments: the contributions of biogeography, geographical distance and discontinuous habitat to genetic differentiation in a blennioid fish, Axoclinus nigricaudus. Mol Ecol. 2001;10:1439–53.
Gaither MR, Bowen BW, Toonen RJ, Planes S, Messmer V, Earle J, et al. Genetic consequences of introducing allopatric lineages of Bluestriped Snapper (Lutjanus kasmira) to Hawaii. Mol Ecol. 2010;19:1107–21.
Eble JA, Rocha LA, Craig MT, Bowen BW. Not all larvae stay close to home: insights into marine population connectivity with a focus on the Brown Surgeonfish (Acanthurus nigrofuscus). J Mar Biol. 2011;2011:1–12.
Schultz JK, Pyle RL, DeMartini E, Bowen BW. Genetic connectivity among color morphs and Pacific archipelagos for the flame angelfish, Centropyge loriculus. Mar Biol. 2007;151:167–75.
Williams ST, Benzie JAH. Evidence of a biogeographic break between populations of a high dispersal starfish: congruent regions within the Indo-West Pacific defined by color morphs, mtDNA, and allozyme data. Evolution. 1998;52:87–99.
Lukoschek V, Waycott M, Marsh H. Phylogeography of the olive sea snake, Aipysurus laevis (Hydrophiinae) indicates Pleistocene range expansion around northern Australia but low contemporary gene flow. Mol Ecol. 2007;16:3406–22.
Bohonak AJ. Dispersal, gene flow, and population structure. Q Rev Biol. 1999;74:21–45.
Kinlan BP, Gaines SD, Lester SE. Propagule dispersal and the scales of marine community process. Divers Distrib. 2005;11:139–48.
Weersing K, Toonen RJ. Population genetics, larval dispersal, and connectivity in marine systems. Mar Ecol Prog Ser. 2009;393:1–12.
Lester SE, Ruttenberg BI, Gaines SD, Kinlan BP. The relationship between dispersal ability and geographic range size. Ecol Lett. 2007;10:745–58.
Portnoy DS, Hollenbeck CM, Renshaw MA, Cummings NJ, Gold JR. Does mating behaviour affect connectivity in marine fishes? Comparative population genetics of two protogynous groupers (family Serranidae). Mol Ecol. 2013;22:301–13.
De Araujo JN, Martins AS. Age and growth of coney (Cephalopholis fulva), from the central coast of Brazil. J Mar Biol Assoc United Kingdom. 2006;86:187–91.
Bernard AM, Feldheim KA, Nemeth R, Kadison E, Blondeau J, Semmens BX, et al. The ups and downs of coral reef fishes: the genetic characteristics of a formerly severely overfished but currently recovering Nassau grouper fish spawning aggregation. Coral Reefs. 2016;35:273–84.
Sadovy de Mitcheson Y, Colin PL. Reef fish spawning aggregations: biology, research and management. Dordrecht: Springer; 2011.
Russell MW, Luckhurst BE, Lindeman KC. Management of spawning aggregations. In: Sadovy de Mitcheson Y, Colin PL, editors. Reef fish spawning aggregations: biology, research and management. Dordrecht: Springer; 2012. p. 371–404.
Harrison HB, Williamson DH, Evans RD, Almany GR, Thorrold SR, Russ GR, et al. Larval export from marine reserves and the recruitment benefit for fish and fisheries. Curr Biol. 2012;22:1023–8.
Almany GR, Hamilton RJ, Bode M, Matawai M, Potuku T, Saenz-Agudelo P, et al. Dispersal of grouper larvae drives local resource sharing in a coral reef fishery. Curr Biol. 2013;23:626–30.
Coleman FC, Koenig CC, Eklund A-M, Grimes CB. Management and conservation of temperate reef fishes in the grouper-snapper complex of the southeastern United States. Am Fish Soc Symp. 1999:233–42.
Musick JA. Life in the slow lane: ecology and conservation of long-lived marine animals. In: American fisheries society symposium, vol. 23. Bethesda, Maryland: American Fisheries Society; 1999.
Sadovy YJ. Reproduction of reef fishery species. In: Polunin NVC, Roberts CM, editors. Reef Fish. Dordrecht: Springer; 1996. p. 15–59.
Sadovy de Mitcheson Y, Craig MT, Bertoncini AA, Carpenter KE, Cheung WWL, Choat JH, et al. Fishing groupers towards extinction: a global assessment of threats and extinction risks in a billion dollar fishery. Fish Fish. 2013;14:119–36.
Thierry C, Sadovy Y, Yeeting B. Plectropomus areolatus. IUCN Red List Threat Species. 2008;2008:e.T64411A12779260.
We would like thank Alec Hughes and David Duchene (James Cook University), Ling Ming Tsang (The Chinese University of Hong Kong), Min Liu (Xiamen University), Noah DesRosiers (King Abdullah University of Science and Technology), and researchers in the Marine Research Centre in Maldives for their assistance in sample collections and/or molecular lab work, as well as inspiring discussion on this study. We also thank Wayne Mallett (James Cook University High Performance Computing Unit) for help with computer analyses. We thank two anonymous reviewers for their constructive comments that improved the manuscript.
This study was supported by a research fund from the Research Grants Council of the Hong Kong Special Administrative Region awarded to YSM (University of Hong Kong) and KHC (The Chinese University of Hong Kong) [project no. HKU 773908 M], funding from KAUST Red Sea Research Center awarded to MLB (King Abdullah University of Science and Technology), and the James Cook University IRA fund to KYM.
Availability of data and materials
Control region sequences will be deposited in GenBank with the accession numbers KM656497-KM656787, MH853841-MH854500. The datasets supporting the conclusions of this article are included within the article and in the Additional file 1.
This research presented and reported in this manuscript was conducted in compliance with the National Health and Medical Research Council (NHMRC) Australian Code of Practice for the Care and Use of Animals for Scientific Purposes, 7th Edition, 2004 and the QLD Animal Care and Protection Act, 2001. The proposed research study received animal ethics approval from the Animal Ethics Committee of James Cook University, with Ethic Approval Number A1659.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional details of sampling, methodology, and results: Table S1. Primers used for PCR amplification of control region sequences. Table S2-S8. Basic information and genetic diversity of control region and microsatellite datasets. Table S9. SAMOVA results. Table S10-S15. Pairwise ΦST and FST of control region and microsatellite datasets. Table S16. Results of linkage disequilibrium test. Figure S1-S6. Results of DAPC analyses. Figure S7. Summary of population genetic structure and reproductive characteristics of five grouper species analysed in this study and Portnoy et al. (2013). (PDF 4110 kb)