- Research article
- Open Access
Evidence for panmixia despite barriers to gene flow in the southern African endemic, Caffrogobius caffer (Teleostei: Gobiidae)
BMC Evolutionary Biologyvolume 8, Article number: 325 (2008)
Oceanography and life-history characteristics are known to influence the genetic structure of marine species, however the relative role that these factors play in shaping phylogeographic patterns remains unresolved. The population genetic structure of the endemic, rocky shore dwelling Caffrogobius caffer was investigated across a known major oceanographic barrier, Cape Agulhas, which has previously been shown to strongly influence genetic structuring of South African rocky shore and intertidal marine organisms. Given the variable and dynamic oceanographical features of the region, we further sought to test how the pattern of gene flow between C. caffer populations is affected by the dominant Agulhas and Benguela current systems of the southern oceans.
The variable 5' region of the mtDNA control region was amplified for 242 individuals from ten localities spanning the distributional range of C. caffer. Fifty-five haplotypes were recovered and in stark contrast to previous phylogeographic studies of South African marine species, C. caffer showed no significant population genetic structuring along 1300 km of coastline. The parsimony haplotype network, AMOVA and SAMOVA analyses revealed panmixia. Coalescent analyses reveal that gene flow in C. caffer is strongly asymmetrical and predominantly affected by the Agulhas Current. Notably, there was no gene flow between the east coast and all other populations, although all other analyses detect no significant population structure, suggesting a recent divergence. The mismatch distribution suggests that C. caffer underwent a population expansion at least 14 500 years ago.
We propose several possible life-history adaptations that could have enabled C. caffer to maintain gene flow across its distributional range, including a long pelagic larval stage. We have shown that life-history characteristics can be an important contributing factor to the phylogeography of marine species and that the effects of oceanography do not necessarily suppress its influence on effective dispersal.
With no obvious physical barriers to gene flow, marine species with large ranges are expected to show genetic homogeneity over long stretches of ocean – a panmixia paradigm . Many studies have, however provided evidence to the contrary, finding strong genetic structure and detecting oceanic barriers to gene flow that often correspond to biogeographic regions [2, 3]. It is now known that several factors, including oceanography [4, 5], life history traits  and contemporary factors such as local adaptation [ and references therein] may strongly influence the genetic structuring of marine organisms despite the passive dispersal opportunities that the marine environment offers.
The South African coastline offers a unique setting to study the relative importance of oceanography on the genetic structuring of marine organisms. There are two contrasting current systems that meet as the continental shelf widens between Cape Agulhas and the Cape Peninsula (Figure 1). These are the warm Agulhas Current, which flows southwards from Mozambique along the eastern coast of the country, and the cold Benguela Current that runs northwards along the western coast . The Agulhas Current deflects away from the coastline near Algoa Bay, resulting in mixing of warm and cool water around the South Coast [[7, 8]; Figure 1]. Based on the distribution patterns of currents, topographical features and species compositions , three biogeographical regions are recognized along the South African coastline. These are the subtropical East Coast (from KwaZulu-Natal to east of Algoa Bay), the warm-temperate South Coast (from Algoa Bay to Cape Point), and the cold-temperate West Coast (from Cape Point to the Namibian border) [; Figure 1].
It has traditionally been accepted that species with good dispersal abilities should show little or no genetic structure over large stretches of ocean, whereas live-bearers and species with short larval durations (or larval retention strategies) should be more strongly structured genetically . However, many closely related species with the same potential dispersal abilities have been found to show contrasting phylogeographic patterns over the same geographic range [11, 12]. In addition, data derived from 70 studies that investigated the population genetic structure of species across the Atlantic-Mediterranean break have recently been re-analyzed and it has been suggested that in this region there may not be a relationship between biological traits and phylogeographic structure of marine species . In South Africa, all intertidal and shallow subtidal species studied to date have exhibited a genetic break over Cape Agulhas (the southernmost point of Africa), irrespective of their dispersal abilities [13–17]. Furthermore, larval duration, swimming abilities and retention strategies (e.g. brooding of eggs) only seem to influence the phylogeography of South African invertebrates on a regional scale (within biogeographic regions), as even the mudprawn (Upogebia africana) and caridean shrimp (Palaemon peringueyi) with long-lived, actively swimming larvae exhibit at least one genetic break .
Caffrogobius caffer is a goby endemic to South Africa. These fish are abundant in rock pools and primarily restricted to the high intertidal zone. As such they are well adapted to survive extreme temperature and salinity ranges [; authors pers. obs]. They have been cited as occurring from False Bay (Cape Peninsula) on the south-west Coast to Delagoa Bay in southern Mozambique , although Hoese  states that these gobies do not occur in Mozambique and we, despite intensive sampling, did not catch any east of Haga Haga on the east coast. Importantly, mark-and-recapture studies strongly suggest that adult movement is limited to nearby pools only , making adult dispersal over long distances highly unlikely. The reproductive biology of C. caffer is not well documented, but insights from the closely related Caffrogobius gilchristi, an estuarine inhabitant, could be used as an indicator. However, the duration of their pelagic larval stage can only be speculated on, compared with findings for closely related species (Strydom, pers. comm.). It is also likely that male gobies guard their eggs until they hatch.
Our objective was to determine whether life history influences the population structure of a widely distributed fish species that is primarily restricted to life in the high intertidal zones. Given the length of coastline investigated (~1300 km) and that C. caffer probably has limited dispersal abilities due to an abbreviated life history where male fish guard the eggs, we expect to find evidence of population genetic structuring along the distribution range, particularly over the regions strongest barrier to gene flow, Cape Agulhas. Further, we hypothesised that gene flow patterns are strongly influenced by the dominant current systems in the range of C. caffer, in particular by the Agulhas Current, and accordingly, we should find asymmetric gene flow patterns between populations. Further, we aim to elucidate whether life-history characteristics are indeed unimportant predictors of genetic structure of South Africa's coastal biota [13–16], i.e. that oceanography is the dominating force in shaping populations in the region.
Sampling protocol and laboratory procedures
Fishes were collected from ten localities covering a wide distributional range of the species using ethylenglycolmonophenylether (phenoxyethanol; Merck) added to rockpools and handnets (Figure 1). Sampling was carried out within the laws of South Africa. Despite intensive sampling, we could not obtain samples from the central and northern Kwazulu-Natal coastline where C. caffer seems to be replaced by C. natalensis. All fishes have been deposited in the JLB Smith collection at the South African Institute for Aquatic Biodiversity, Grahamstown. A small piece of muscle tissue was digested overnight in an extraction buffer containing 20 μl ProteinaseK at 55°C. DNA was extracted using the DNeasy™ Tissue Kit (Qiagen), following the protocol for animal tissues. Primers and PCR cycling conditions of Lee et al.  were used for amplification of the 5'-end of the mitochondrial DNA control region; a marker that has proved useful in numerous population genetic studies, especially of fishes. PCR products were gel purified using the Illustra™ DNA and Gel Purification Kit (GE Healthcare). Cycle sequencing was carried out using BigDye Terminator chemistry (Applied Biosystems) and the products were analysed on an automated sequencer (AB 3100, Applied Biosystems).
Sequences were edited and aligned by eye with BioEdit ver220.127.116.11 . Sequences were checked against goby control region sequences on GenBank, with a tblastx search, to ensure that the correct region was amplified. Haplotypes were identified in Collapse1.2 http://darwin.uvigo.es. Haplotype (h) and nucleotide diversity (π) indices were calculated for each population individually and for the whole dataset, in Arlequin3.1 . All haplotypes have been deposited in [GenBank: EU335086–EU335140].
Analysis of Molecular Variance (AMOVA) was carried out in Arlequin3.1 to investigate the geographical distribution of haplotypes and provide the fixation indices for between- and within population variation. Φ ST values were calculated to take the number of mutational changes between the haplotypes into consideration . SAMOVA was used to calculate the geographic distribution of genetic structure. This maximises the proportion of the total genetic variation between groups of populations, as well as identifying possible genetic barriers between them, without pre-defining populations as is necessary for AMOVA . The geographical coordinates of each sampling location were obtained from Google Earth®. To detect possible isolation by distance, a Mantel test was performed in Arlequin3.1 using 10 000 permutations and the measured geographical distances between localities determined from a scaled map as the shortest land-sea interface distance, but excluding bays.
Fu's F S test  was used to test for possible population expansion, for the whole dataset, as well as for each sampling location separately. A mismatch distribution  was also calculated to further investigate the demographic history of C. caffer populations. The Sum of Squared deviations and Harpending's Raggedness index were calculated to test for the goodness of fit of the data in the mismatch analyses. A parsimony minimum spanning network was constructed using TCS1.21 , which is a more appropriate way than constructing bifurcating phylogenetic trees to view the evolutionary relationship between closely related haplotypes where ancestry may not be strictly bifurcating .
The time of expansion was calculated with the formula T = τ/2 μ, where μ = generation time × number of base pairs per sequence × mutation rate for the marker used , and τ was calculated in Arlequin3.1. Although the generation time of C. caffer is not known, estuarine fishes (including C. gilchristi) are said to disperse to the marine environment to join spawning adults after one to three years residence time in the estuary . Estuarine spawning does; however, occur from time to time , and therefore it is reasonable to assume a generation time of between one and three years for calculations of divergence times. A mutation rate of 3.6% per million years is generally used for teleosts , although the fast evolving control region of the mtDNA molecule has been found to have a mutation rate of 11–13% per million years in the white sturgeon . Based on this, two separate calculations were done – one using a 3.6% per million years mutational rate and one based on a 11% per million years mutational rate.
Migration between sampling localities
To estimate past migration rates, as well as the directionality of gene flow between populations, Migrate-n version 2.4 was used [34, 35]. In order to maximize the statistical power of our gene flow analyses we sought to minimize the number of parameters [see e.g. ] and constructed a stepping-stone model with asymmetrical gene flow. Our study system is particularly suited to this model, given the linear nature of the southern African coastline. For the stepping-stone migration model, two analytical runs were conducted, an initial short run, followed by a second longer run. For both runs, the starting values of the population mutation parameter and the ratio between the immigration rate and the respective population and mutation rate per generation were estimated from F ST values . For the long-run 10 short-chains, each with a total of 25,000 generations and a sampling increment of 20 generations, and two long-chains each with a total of 250,000 generations and a sampling increment of 50 generations were run twice. A total of 50,000 and 12,500,000 genealogies (recorded steps multiplied by the sampling increment) were visited by the short and long chains, respectively. For both, the short and long chains, the first 10,000 genealogies were discarded (the burnin). An adaptive heating scheme with four chains (starting values of 5.00, 2.50, 1.50, 1.00) and a swapping interval of one was used to ensure that efficient mixing occurred. For the other settings, default values were implemented.
Between 21 and 26 fishes were obtained from each location; a total of 242 individuals (Table 1). After editing and alignment by eye in BioEdit, 390 base pairs remained for analyses. There were 26 polymorphic nucleotide positions (6.67% variable sites), with 20 observed transitions and eight transversions. Fifty-five haplotypes were recovered, of which 28 were restricted to single individuals. Haplotype diversity (h) was found to be very high (h = 0.956, ± 0.004) and nucleotide diversity (π), was low (π = 0.010, ± 0.005). Interestingly, theta, as estimated by MIGRATE [34, 35] was lowest for the two peripheral populations at Wooley's Pool and Haga Haga.
Pairwise Φ ST values among sampled localities were low and all but one of these comparisons were not significant (Jongensfontein and Port Alfred, Φ ST = 0.046, P < 0.05). SAMOVA analyses maximized at two groups with Jongensfontein grouped separately from the rest of the sampling sites, although this grouping was not biologically meaningful as Jongensfontein is in the middle of the distribution range (F CT = 0.022, P < 0.05). No significant relationship between genetic and geographical distance was detected by the Mantel test (0.072, P > 0.05).
Fu's F S test for all sites was highly negative and significant (Table 1), as would be expected for populations that have undergone recent demographic change . The mismatch distribution did not differ significantly from the null model of sudden expansion, with the observed and expected curves both indicative of population expansion (data not shown). The parsimony haplotype network constructed recovered no clear phylogeographic structure, indicating panmixia (Figure 2). There are two dominant haplotypes, each being present in 24 individuals (~10% of individuals), with representatives from all three biogeographic regions. Although all populations were characterized by between one and five private haplotypes (Table 1), extensive haplotype sharing occurred between the West Coast, South Coast and East Coast of the southern African coastline.
When using τ = 3.859 and the 11% per million years mutation rate for the mtDNA control region , the time of population expansion is estimated to be between 15 000 (3 year generation time) and 45 000 (1 year generation time) years ago, or between 45 800 and 138 000 years ago when estimated with a more conserved mutation rate of 3.6% per million years .
Migration between sampling localities
Interestingly, coalescent analyses of gene flow revealed strong asymmetry along the distributional range of C. caffer (Figure 1, Table 2). There was limited gene flow from western to eastern populations, with strong gene flow from east to west, suggesting an important influence of the Agulhas Current. Notably, there is limited gene flow between Port Alfred and Haga Haga, the two most eastern populations and all other sampling sites, suggesting that they are becoming effectively isolated.
The mitochondrial control region proved to have sufficient variability for population genetic analysis. The combination of high haplotypic diversity and comparatively low nucleotide diversity is typical of the signature of population expansion after a bottleneck or a founder event [37–40], which is generally associated with little or no population structure [3, 40]. Both the mismatch distribution and the significant Fu's F S (-25.6, P < 0.001) support that for C. caffer both as a species and at the individual population level has undergone recent population expansion (Table 1). The time of expansion is indicative of a growing population since the late Pleistocene (275 000 to 15 000 years ago, depending on the mutation rate and generation time used). Interestingly, similar estimates for expansion time have been found in other South African marine organisms, for example the lobsters Palinurus delagoae Palinurus gilchristi  and Jasus tristani , as well as a hake species, Merluccius capensis , and it is likely that changing sea levels and temperatures during this time contributed to demographic changes.
Interestingly, although the two populations at the edge of C. caffer distribution also show signals of population expansions, their theta (4N e mμ) values are lower than when compared to all other sampling sites. This may be a signal of an expanding population into these areas or alternatively may reflect source-sink dynamics at the putative range edge of this species.
Population structure and gene flow patterns
SAMOVA failed to identify any biologically meaningful groups and AMOVA revealed panmixia between populations of all three biogeographic regions. Further, isolation by distance  should not be expected in species with planktonic larvae that can successfully disperse over great distances and that are able to exchange individuals between the populations at the edges of their distribution ranges . In concordance with this theory, significant isolation by distance was only found for South African organisms with possible larval retention or direct development and not for those with pelagic larvae [14–16] including Caffrogobius caffer. As it is highly unlikely that adult C. caffer disperse over large distances, the only possible dispersal option is that by larvae.
Interestingly, despite the signatures of a mtDNA panmictic population of C. caffer in South Africa, coalescent analyses of gene flow show that there is very limited gene flow (west to east N m = 0 (0,0); east to west N m = 1.5 (0.7, 2.5) between east coast (Port Alfred and Haga Haga) and all other populations (Figure 1, Table 2), suggesting that the east coast is effectively isolated. This is congruent with gene flow patterns in the clinid fish Clinus cottoides , where analyses also recovered no gene flow to or from the east coast populations. This is probably a function of distance between suitable rocky shore habitats. As one moves towards the east coast of South Africa, rocky shore habitats become progressively more sparse, with large interspersing tracts of sandy beaches and it probably becomes much less likely that larvae are able to migrate between populations. One possible reason for the absence of significant population structure between the east coast and all other populations is that the reduction in gene flow may be a geologically recent event, and that not enough time has passed to fix any subsequent population specific mutations at the mtDNA level.
Notably, the Agulhas Current plays a large role in the migration patterns of these fishes. There is little gene flow from west to east, even between closely situated populations. In contrast, extensive gene flow occurs from the east to the west coast, which is in the direction of the flow of the Agulhas Current. Some fishes have been shown to utilise the inshore counter-current that runs up the east coast (sardines;  and Clinus cottoides, ), yet this is not evident for C. caffer. This may be as a result of a longer pelagic larval stage, which allows C. caffer to disperse with the Agulhas Current, unlike fishes which either actively seek to use the counter-current or that have a very short pelagic larval stage.
Remarkably, these results contrast sharply to what has been found for rocky shore and shallow intertidal organisms in South Africa to date, which have shown population genetic structuring irrespective of their life-history patterns and pelagic larval duration [13–17]. Our most surprising finding is that gene flow is taking place across Cape Agulhas, which has been identified as major oceanographic barrier to gene flow in all intertidal and shallow water inhabitants studied to date [13–15, 17]. Based on these findings it has been stated that oceanography strongly influences the dispersal potential of marine organisms in South Africa, possibly more so than having a pelagic larval stage [15, 16]. Our results suggest that C. caffer has better dispersal abilities than expected (based on mark-recapture studies of adults) and that this may overcome the effects of such major oceanographic barriers.
Dispersal capabilities of Caffrogobius caffer
Having the capability to release a large amount of eggs into the open ocean at the best possible time should greatly increase the chance of effective dispersal between populations, as in theory it only takes one (practically 5–10) successful migrant per generation to maintain some gene flow between populations . In South Africa, larvae of members of the Gobiidae are numerically dominant in the nocturnal ebb tide – a synchronized hatching strategy used by many estuarine fishes to promote dispersal into the ocean [31, 46]. Caffrogobius gilchristi in particular can potentially release millions of larvae into the open ocean during a single tidal cycle, thereby increasing their dispersal success . As larvae of the genus Caffrogobius have not to date been identified to species level for all six species in the genus, it is not possible to know how numerically dominant C. caffer is in the surfwater (Strydom, personal communication). Judging from its panmictic mtDNA pattern though, we hypothesise that the species most probably is as successful as C. gilchristi in terms of reproductive output and dispersal. Our findings are suggestive of C. caffer having a longer larval duration than expected, for it is highly unlikely that individuals could be exchanged between the East Coast and South-west Coast within a matter of days if utilizing the large-scale Agulhas current system.
Together with larval duration, larval behaviour may also play an important role in the dispersal abilities of C. caffer, like it does in many other species [e.g. [4, 47]]. Larvae of cryptobenthic fishes, including individuals from the family Gobiidae, usually possess good swimming and orientation abilities at time of settlement [48, 49], which may greatly alter the expected population connectivity of the species . Since post-flexion larvae in South African waters are known to be excellent swimmers , the wide dispersal of C. caffer larvae in local current systems seems probable.
Our results show that successful dispersal across major oceanographic barriers such as Cape Agulhas is possible for some marine species and that availability of habitat is another key determinant in population structuring of marine organisms along the South African coast. We have shown that life-history characteristics could in fact be an important contributing factor to the population genetic structuring of marine species in southern Africa and illustrate that the influence on the effective dispersal of marine organisms is not necessarily suppressed by the effects of oceanography as has recently been advocated [3, 15, 16].
Vermeij G: The dispersal barrier in the tropical Pacific: implications for molluscan speciation and extinction. Evolution. 1987, 41: 1046-1058. 10.2307/2409190.
Barber PH, Palumbi SR, Erdmann MV, Moosa KM: A marine Wallace's Line?. Nature. 2000, 406: 692-693. 10.1038/35021135.
Patarnello T, Volckaert FAMJ, Castilho R: Is the Atlantic-Mediterranean trasition a phylogeographical break?. Mol Ecol. 2007, 16: 4426-4444. 10.1111/j.1365-294X.2007.03477.x.
Chambers MD, van Blaricom GR, Hauser L, Utter F, Friedman CS: Genetic structure of black abalone (Haliotis cracherodii) populations in the California islands and central California: impacts of larval dispersal and decimation from withering syndrome. J Exp Mar Biol Ecol. 2006, 331: 173-185. 10.1016/j.jembe.2005.10.016.
Hemmer-Hansen J, Nielsen EEG, Grønkjær P, Loeschcke V: Evolutionary mechanisms shaping the genetic population structure of marine fishes; lessons from the European flounder (Platichthys flesus L.). Mol Ecol. 2007, 16: 3104-3118. 10.1111/j.1365-294X.2007.03367.x.
Riginos C, Nachman MW: Population subdivision in the marine environments: the contribution of biogeography, geographical distance and discontinuous habitat to genetic differentiation in a blennioid fish, Axoclinus nigricaudus. Mol Ecol. 2001, 10: 1439-1453. 10.1046/j.1365-294X.2001.01294.x.
Branch GM, Branch ML: The Living Shores of Southern Africa. 1995, Cape Town: Struik Publishers
Shillington FA: Oceanography of the southern African region. Smiths' Sea Fishes. Edited by: Smith MM, Heemstra PC. 1995, Cape Town: Struik Publishers, 22-28. 2
Stephenson TA, Stephenson A: Life between the tidemarks on rocky shores. 1972, San Fransisco: W.H Freeman and Company
Palumbi SR: Genetic divergence, reproductive isolation and marine speciation. Annu Rev Ecol Syst. 1994, 25: 547-572. 10.1146/annurev.es.25.110194.002555.
Bargelloni L, Alarcon JA, Alvarez MC, Penzo E, Magoulas A, Reis C, Patarnello T: Discord in the family Sparidae (Teleostei): divergent phylogeographical patterns across the Atlantic-Mediterranean divide. J Evol Biol. 2003, 16: 1149-1158. 10.1046/j.1420-9101.2003.00620.x.
Bird CE, Holland BS, Bowen BW, Toonen RJ: Contrasting phylogeography in three endemic Hawaiian limpets (Cellana spp.) with similar life histories. Molecular Ecology. 2007, 16: 3173-3186. 10.1111/j.1365-294X.2007.03385.x.
Evans BS, Sweijd NA, Bowie RCK, Cook PA, Elliott NG: Population genetic structure of the perlemoen Haliotis midae in South Africa: evidence of range expansion and founder events. Mar Ecol Prog Ser. 2004, 270: 163-172. 10.3354/meps270163.
Teske PR, Papadopoulos I, Zardi GI, McQuaid CD, Griffiths CL, Edkins MT, Barker NP: Implications of life history for genetic structure and migration rates of five southern African coastal invertebrates: planktonic, abbreviated and direct development. Mar Biol. 2007, 152: 697-711. 10.1007/s00227-007-0724-y.
Teske PR, Froneman PW, McQuaid CD, Barker NP: Phylogeographic structure of the caridean shrimp Palaemon peringueyi in South Africa: further evidence for intraspecific genetic units associated with marine biogeographic provinces. Afr J Mar Sci. 2007, 29: 253-258. 10.2989/AJMS.2007.29.2.9.192.
Zardi GI, McQuaid CD, Teske PR, Barker NP: Unexpected genetic structure of mussel populations in South Africa: indigenous Perna perna and invasive Mytilus galloprovincialis. Mar Ecol Prog Ser. 2007, 337: 135-144. 10.3354/meps337135.
von der Heyden S, Prochazka K, Bowie RCK: Significant population structure and asymmetric gene flow patterns amidst expanding populations of Clinus cottoides (Perciformes, Clinidae): application of molecular data to marine conservation planning in South Africa. Mol Ecol. 2008, 17: 4812-4826. 10.1111/j.1365-294X.2008.03959.x.
Butler GS: Aspects of the biology of Caffrogobius caffer (Günther) (Pisces: Teleostei: Gobiidae) in the Eastern Cape. MSc Thesis. 1980, Rhodes University, Grahamstown
Goren M: A review of the southern African Gobiid fish genus Caffrogobius Smitt, 1990. 1996, J.L.B. Smith Institute of Ichthyology, Special Publication, 57:
Hoese DF: Family No. 240: Gobiidae. Smiths' Sea Fishes. Edited by: Smith MM, Heemstra PC. 1995, Cape Town: Struik Publishers, 774-807. 2
Lee W-J, Conroy J, Howell WH, Kocher TD: Structure and evolution of teleost mitochondrial control regions. J Mol Evol. 1995, 41: 54-66. 10.1007/BF00174041.
Hall TA: BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucleic Acids Symp Ser. 1999, 41: 95-98.
Excoffier L, Laval G, Schneider S: Arlequin version 3: an integrated software package for population genetics data analysis. Evol Biol Online. 2005, 1:
Allendorf FW, Luikart G: Conservation and the genetics of populations. 2007, Oxford: Blackwell Publishing Ltd
Dupanloup I, Schneider S, Excoffier L: A simulated annealing approach to define the genetic structure of populations. Mol Ecol. 2002, 11: 2571-2581. 10.1046/j.1365-294X.2002.01650.x.
Fu Y-X: Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics. 1997, 147: 915-925.
Harpending HC: Signature of ancient population growth in a low-resolution mitochondrial DNA mismatch distribution. Hum Biol. 1994, 66 (4): 591-600.
Clement M, Posada D, Crandall KA: TCS: a computer program to estimate gene genealogies. Mol Ecol. 2000, 9: 1657-1660. 10.1046/j.1365-294x.2000.01020.x.
Posada D, Crandall KA: Intraspecific gene genealogies: tree grafting into networks. TREE. 2001, 16: 37-45.
Rogers AR, Harpending H: Population growth makes waves in the distribution of pairwise generic differences. Mol Biol Evol. 1992, 9: 552-569.
Whitfield AK: Biology and ecology of fishes in South African estuaries. 1998, Ichthyological Monographs of the J. L. B. Smith Institute of Ichthyology, No. 2
Donaldson KA, Wilson RW: Amphi-panamaic germinates of snook (Percoidei: Centropomidae) provide a calibration of the divergence rate in the mitochondrial DNA control region of fishes. Mol Phyl Evol. 1999, 13: 208-213. 10.1006/mpev.1999.0625.
Brown JR, Beckenbach AT, Smith MJ: Intraspecific DNA sequence variation of the mitochondrial control region of white sturgeon (Acipenser transmontanus). Mol Biol Evol. 1993, 10: 326-341.
Beerli P, Felsenstein J: Maximum likelihood estimation of a migration matrix and effective population sizes in two populations using a coalescent approach. Genetics. 1999, 152: 763-773.
Beerli P, Felsenstein J: Maximum likelihood estimation of a migration matrix and effective population sizes in n subpopulations by using a coalescent approach. Proc Nat Acad Sci USA. 2001, 98: 4563-4568. 10.1073/pnas.081068098.
Bowie RCK, Fjeldså J, Hackett SJ, Bates JM, Crowe TM: Coalescent models reveal the relative roles of ancestral polymorphism, vicariance and dispersal in shaping phylogeographical structure of an African montane forest robin. Mol Phylogen Evol. 2006, 38: 171-188. 10.1016/j.ympev.2005.06.001.
Grant WS, Bowen BW: Shallow population histories in deep evolutionary lineages of marine fishes: insights from sardines and anchovies and lessons for conservation. J Hered. 1998, 89: 415-126. 10.1093/jhered/89.5.415.
Waples RS: Separating the wheat from the chaff: patterns of genetic differentiation in high gene flow species. J Hered. 1998, 89: 438-450. 10.1093/jhered/89.5.438.
Zhang J, Cai Z, Huang L: Population genetic structure of crimson snapper Lutjanus erythropterus in East Asia, revealed by analysis of the mitochondrial control region. J Mar Sci. 2006, 63: 693-704.
von der Heyden S, Lipinski MR, Matthee CA: Mitochondrial DNA analyses of the Cape hakes reveal an expanding population for Merluccius capensis and population structuring for mature fish in Merluccius paradoxus. Mol Phyl Evol. 2007, 42: 517-527. 10.1016/j.ympev.2006.08.004.
Gopal K, Tolley KA, Groeneveld JC, Matthee CA: Mitochondrial DNA variation in spiny lobster Palinurus delagoae suggests genetically structured populations in the south-western Indian Ocean. Mar Ecol Prog Ser. 2006, 319: 191-198. 10.3354/meps319191.
Tolley KA, Groeneveld JC, Gopal K, Matthee CA: Mitochondrial DNA panmixia in spiny lobster Palinurus gilchristi suggests a population expansion. Mar Ecol Prog Ser. 2005, 297: 225-231. 10.3354/meps297225.
von der Heyden S, Groeneveld JC, Matthee CA: Long current to nowhere?' – genetic connectivity of Jasus tristani in the southern Atlantic Ocean. Afr J Mar Sci. 2007, 29: 491-497. 10.2989/AJMS.2007.29.3.15.345.
Wright S: Isolation by distance. Genetics. 1943, 28: 114-138.
Planes S, Fauvelot C: Isolation by distance and vicariance drive genetic structure of a coral reef fish in the Pacific ocean. Evolution. 2002, 56: 378-399.
Strydom NA, Wooldridge TH: Diel and tidal variations in larval fish exchange in the mouth region of the Gamtoos Estuary, South Africa. Afr J Aquat Sci. 2005, 30: 131-140.
Matthee CA, Cockcroft AC, Gopal K, von der Heyden S: Mitochondrial DNA variation of the west coast rock lobster, Jasus lalandii: marked genetic diversity differences among sampling sites. Mar Freshwater Res. 2007, 58: 1130-1135. 10.1071/MF07138.
Hickford MJH, Schiel D: Comparative dispersal of larvae from demersal versus pelagic spawning fishes. Mar Ecol Prog Ser. 2003, 252: 255-271. 10.3354/meps252255.
Leis JM, Hay AC, Trnski T: In situ ontogeny behaviour in pelagic larvae of three temperate, marine, demeral fishes. Mar Biol. 2005, 148: 655-669. 10.1007/s00227-005-0108-0.
Paris CB, Chérubin LM, Cowen RK: Surfing, spinning, or diving from reef to reef: effects on population connectivity. Mar Ecol Prog Ser. 2007, 347: 285-300. 10.3354/meps06985.
Thanks to V. Bowie and C. Muller for help with the collection of fishes and Dr N. Strydom of the Nelson Mandela Metropolitan University for helpful discussion on the the larval biology of members of the genus Caffrogobius. NM was supported by a NRF Freestanding Honours bursary and SvdH by a NRF Prestigious Free-standing postdoctoral fellowship. This project was funded by the National Research Foundation (NRF) of South Africa (FA2005040400057).
MN carried out the majority of the field and laboratory work and wrote the draft version of the manuscript. CAM participated in field work, helped with the interpretation of data and finalisation of the manuscript. RCKB collected samples, helped with statistical analyses and contributed to writing the manuscript. SvdH conceived the study, helped with sampling and data analysis and contributed to writing the final manuscript. All authors have read and approved the final manuscript.