- Research article
- Open Access
Oceanic dispersal barriers, adaptation and larval retention: an interdisciplinary assessment of potential factors maintaining a phylogeographic break between sister lineages of an African prawn
© Teske et al; licensee BioMed Central Ltd. 2008
- Received: 23 June 2008
- Accepted: 24 December 2008
- Published: 24 December 2008
Genetic breaks separating regional lineages of marine organisms with potentially high broadcasting abilities are generally attributed either to dispersal barriers such as currents or upwelling, or to behavioural strategies promoting self-recruitment. We investigated whether such patterns could potentially also be explained by adaptations to different environmental conditions by studying two morphologically distinguishable genetic lineages of the estuarine mudprawn Upogebia africana across a biogeographic disjunction in south-eastern Africa. The study area encompasses a transition between temperate and subtropical biotas, where the warm, southward-flowing Agulhas Current is deflected away from the coast, and its inshore edge is characterised by intermittent upwelling. To determine how this phylogeographic break is maintained, we estimated gene flow among populations in the region, tested for isolation by distance as an indication of larval retention, and reared larvae of the temperate and subtropical lineages at a range of different temperatures.
Of four populations sampled, the two northernmost exclusively included the subtropical lineage, a central population had a mixture of both lineages, and the southernmost estuary had only haplotypes of the temperate lineage. No evidence was found for isolation by distance, and gene flow was bidirectional and of similar magnitude among adjacent populations. In both lineages, the optimum temperature for larval development was at about 23°C, but a clear difference was found at lower temperatures. While larvae of the temperate lineage could complete development at temperatures as low as 12°C, those of the subtropical lineage did not complete development below 17°C.
The results indicate that both southward dispersal of the subtropical lineage inshore of the Agulhas Current, and its establishment in the temperate province, may be limited primarily by low water temperatures. There is no evidence that the larvae of the temperate lineage would survive less well in the subtropical province than in their native habitat, and their exclusion from this region may be due to a combination of upwelling, short larval duration with limited dispersal potential near the coast, plus transport away from the coast of larvae that become entrained in the Agulhas Current. This study shows how methods from different fields of research (genetics, physiology, oceanography and morphology) can be combined to study phylogeographic patterns.
- Zoeal Stage
- Biogeographic Province
- Agulhas Current
- Larval Retention
- Phylogeographic Break
The marine habitat is considered to contain few physical barriers to dispersal, and marine organisms can therefore potentially disperse over large distances [1, 2]. The collection of larvae of coastal marine invertebrates at often great distances from the coast , and the often limited genetic structuring of widely separated populations of high-dispersal species [4, 5] reinforce this notion. Despite this, sudden breaks in the distributions of marine organisms with potentially high dispersal abilities have been identified throughout the world, and genetic methods have often identified the presence of closely related sister lineages on either side of these breaks [6–8]. Such phylogeographic patterns are often attributed to allopatric divergence resulting from the formation of past land bridges, with subsequent maintenance because propagules of populations on either side of the former barrier primarily recruited to the parent population [6, 9]. There are also a number of examples of phylogeographic breaks that cannot be directly linked to geological vicariance events, but for which incomplete isolation of populations in the form of sudden discontinuities in environmental parameters has been invoked. Examples include the freshwater plume of the Amazon-Orinoco outflow that inhibits gene flow of marine organisms between the Caribbean and Brazil  and coldwater upwelling in south-western Africa that limits dispersal of tropical and subtropical marine organisms from the Indian Ocean to the Atlantic Ocean .
We investigated how the phylogeographic break between the two mudprawn lineages is likely to be maintained by using a combination of genetic and physiological methods. As temperature is the most obvious environmental parameter that changes along the Wild Coast, we studied the development of larvae of the two lineages at a range of different temperatures. Adult mudprawns are restricted to estuaries and are thus physiologically adapted to tolerate variations and extremes in temperature and salinity, suggesting that changes in sea surface temperatures along the Wild Coast are unlikely to impose strong selection pressure on them. Larval development, however, must obligatorily be completed in the marine environment. We therefore hypothesised that if the lineages are adapted to the temperature regimes characteristic of their biogeographic province, then this would be most readily detectable in the more sensitive larvae.
The influence of the following three factors on maintaining the phylogeographic break was investigated. In each case, a number of expectations in terms of genetic and physiological patterns are listed.
The core of the southeast African coastal upwelling region is located in the vicinity of Port Alfred, but upwelling may occur as far north as Port St Johns (Fig. 1; ). If upwelling represented a significant dispersal barrier, then one would expect the phylogeographic break separating the two mudprawn lineages to be near Port Alfred. The development and survival of larvae should be severely affected at the low water temperatures that may be experienced during upwelling events (≤ 17°C ), but would not necessarily differ for the two lineages. In addition, it is possible that the southward-flowing Agulhas Current may limit gene flow from the temperate to the subtropical province. In contrast, larvae of the subtropical lineage can become temporarily entrained in the Agulhas Current may eventually settle in estuaries in the warm-temperate province. Gene flow would then be expected to be mostly unidirectional in a south-westerly direction.
Larvae of the two mudprawn lineages may be adapted to temperature regimes characteristic of each biogeographic province, which in nearshore regions range from around 14 to 22°C in the warm-temperate province  and from around 20 to >26°C in the subtropical province . A few degrees of difference in mean sea surface temperature between northern and southern Wild Coast may be sufficient to prevent each lineage from establishing itself in the habitat of the other. In terms of genetic structure, one would expect a gradual changeover of lineages along the Wild Coast. In terms of larval development, differences in temperature tolerance ranges and/or optimum conditions for growth and development should be identified.
Recent studies have questioned the notion that the benefits traditionally associated with long-distance dispersal are indeed advantageous in terms of recruitment success [29, 30]. Since larvae may disperse away from habitat suitable for settlement, it has been argued that retention near the parent population may be more beneficial since conditions there are obviously suitable for survival (so-called self-recruitment ). The larvae of several crustacean species have been shown to employ vertical movement to avoid being swept away from suitable habitats [32, 33]. If such mechanisms are employed by mudprawn larvae and contribute significantly towards maintaining genetic structure, one would expect each population to have a unique combination of mtDNA haplotypes. This should manifest itself in a pattern of isolation by distance along the coast.
Four estuaries spanning the transition zone between the warm-temperate and subtropical biogeographic provinces were selected on the basis of having permanent mudprawn populations at high densities, as established during earlier surveys [13, 14]. These were the Gqunube Estuary, a site at the confluence of the Nxaxo and Ngqusi rivers (hereafter referred to as Wavecrest), the Mthatha Estuary and the Mngazana Estuary (Fig. 1). In each of these systems, tissue samples from 35 mudprawns were collected, and a portion of the cytochrome oxidase subunit I gene (COI) was amplified and sequenced.
Migration rates among Upogebia africana populations
241 ± 148
(38 – 698)
1 ± 0
(0 – 46)
407 ± 390
(55 – 746)
155 ± 115
(26 – 415)
0 ± 1
(0 – 35)
67 ± 26
(1 – 291)
To determine whether larval retention in the region may have resulted in genetic differentiation among populations that is correlated with geographic distance, we tested for isolation by distance (IBD) using a Mantel test. As tests for IBD are only meaningful when performed on single regional genetic lineages  and require a minimum of three populations (see Methods), only samples from the subtropical lineage were included. No significant correlation was found between genetic and geographic distance among the subtropical lineage haplotypes of the Mngazana, Mthatha and Wavecrest populations (r = 0.0007; P = 0.33), indicating that if larval retention occurs in mudprawns, then it occurs too rarely in the study area to result in genetic structure that can be adequately described by a model of isolation by distance.
Relationships between genetic and morphological data
Samples were not originally collected with the intention of identifying the subdistal spine on the merus of the first pereiopods, so the relevant portion of the pereiopods was absent in 22 specimens of the subtropical lineage and seven specimens of the temperate lineage. Of the remaining 61 individuals belonging to the subtropical lineage and 50 specimens of the temperate lineage, a single temperate lineage individual from the Gqunube population had the spine, one subtropical individual from the Mthatha population did not have it, and it was not well developed in a small juvenile from the Wavecrest population belonging to the subtropical lineage. An examination of samples collected previously [13, 14] revealed that none of the mudprawns collected west of the Gqunube Estuary had the spine (n = 85), and that it was present in all of the specimens collected north of the Mngazana Estuary (n = 8). These results indicate that the subdistal spine is generally a useful diagnostic feature to distinguish the mtDNA lineages of U. africana, but it is likely that there is some interbreeding along the Wild Coast (i.e. the region where the ranges of the two lineages overlap).
Larval rearing experiments
Larvae of the two mudprawn lineages were reared at various constant temperatures, and their survival and duration of development were compared. Larvae used in these experiments came from two estuaries that are each located well within one of the two major biogeographic provinces inhabited by Upogebia africana, the temperate Swartkops Estuary and the subtropical Mngazana Estuary. This choice of populations ensured that in each case, only representatives from a single lineage were studied. In both lineages, the larval sequence comprised three zoeal stages followed by a metamorphic moult to the Decapodid.
The possibility that regional genetic lineages of a species present on either side of a phylogeographic break may have undergone some adaptive divergence is not usually investigated in studies of this nature, although it has previously been suggested [17, 36]. Phylogeographic studies are often exclusively conducted by geneticists, but the most suitable genetic methods available for studying adaptation require complete genome analyses and are only feasible for economically important or charistmatic species . The present study shows that the coupling of methods from different fields of research (genetics, physiology, oceanography and morphology) presents an alternative approach to determine whether genetic patterns identified using presumably neutral markers are linked to environmental discontinuities rather than merely being the result of stochastic processes. The results indicate that the distribution patterns of the two Upogebia africana lineages are probably maintained by a combination of environmental and adaptive factors. We therefore suggest that in cases where genetic structure is identified and where environmental differences are evident, the assumption that some local adaptations may exist seems valid and warrants further research into how each lineage is adapted to its environment.
Below, we discuss how each of the three factors hypothesised to be involved in maintaining the Wild Coast phylogeographic break is likely to impact on the two mudprawn lineages.
Upwelling was identified as a potentially important factor limiting dispersal of mudprawns along the Wild Coast, as cold sea surface temperatures typical during upwelling events negatively affected the development of the larvae of both lineages. The larvae of the subtropical lineage were unable to complete development at the low temperatures experienced during upwelling events, suggesting that few larvae that disperse inshore of the Agulhas Current survive transport towards the warm-temperate biogeographic province during upwelling events. Although the effects of cold water were not as drastic in the case of the warm-temperate lineage, a decreased survival rate was recorded, and larval development took approximately twice as long at 17°C, and even longer at colder temperatures. As prolonged larval duration increases the chances of predation or displacement from suitable habitat, the upwelling region thus represents an incomplete dispersal barrier that primarily influences recruitment success.
The core of the upwelling cell is located south of the southernmost estuary sampled, whereas the point at which the transition zone of mudprawn lineages occurred was identified near Wavecrest. This apparent discrepancy should perhaps not be overemphasized, because upwelling may occur over a large area that stretches from Algoa Bay in the south to Port St Johns in the north (Fig. 1), and changes in the genetic composition of mudprawn populations (haplotypes from the warm-temperate lineage were previously found in the Mthatha Estuary ) indicate that the exact position of the transition zone may vary over time. In other coastal invertebrates that are represented in the Wild Coast region by distinct temperate and subtropical mtDNA lineages, the exact positions of phylogeographic breaks or clines differed considerably. In the limpet Patella granularis, the ranges of the two lineages overlapped at Coffee Bay (near the Mthatha Estuary ), those of the prawn Palaemon peringueyi at the Mbhanyana Estuary (south of the Mthatha Estuary ), those of the snail Nassarius kraussianus at Port Alfred , those of the mussel Perna perna between approximately Port Alfred and the Haga Haga Estuary (Fig. 1), and a sharp phylogeographic break just south-west of Port Alfred was identified in the cumacean Iphinoe truncata [13, 14].
Although the genetic data-set used has limitations because a single locus (mtDNA) was sequenced, the fact that bidirectional gene flow was identified indicates that at the relatively small scale at which sampling was conducted, the effect of the southward-flowing Agulhas Current did not seem to impact strongly on gene flow. Oceanographic data indicate that currents in near-shore regions are mostly wind-driven, and that the direction of the wind changes frequently [21, 40], suggesting that gene flow among mudprawn populations in this region may mostly occur close to the coast. Genetic data indicate that dispersal by the Agulhas Current becomes important at a larger scale (e.g. over the entire distribution of the subtropical mudprawn lineage), as the migrants that become entrained in the Agulhas Current in the north of the species' distribution and are moved over greater distances have an impact on the genetic diversity of populations downstream . It is also probable that a high proportion of larvae that approach the ACDZ from both south and east and that become entrained in the Agulhas Current are transported away from the coast and are unable reach suitable habitat in time to complete development. This would indicate that the Agulhas Current itself represents a dispersal barrier to both lineages, but the scale of sampling employed in this study is unsuitable to test this hypothesis. Nicastro, Zardi, Roberts and McQuaid (in prep.) released drifters from the eastern south coast during both summer and winter and found that none moved in a north-western direction beyond the Haga Haga Estuary (Fig. 1). All eventually became entrained in the Agulhas Current and were transported away from the coast. Drifters were released ≥ 10 km offshore, and the paths taken may thus not fully indicate how planktonic larvae disperse in the nearshore region. Nonetheless, these results might indicate that larval loss is high in the ACDZ and could further explain why mudprawn larvae of the warm-temperate lineage are unable to establish themselves in the subtropical province.
Adaptation of the two mudprawn lineages to environmental conditions prevalent in their biogeographic provinces is likely to constitute an important additional factor limiting their mixing. The lower tolerance of larvae from the subtropical lineage to cold water temperature thus not only results in upwelling presenting a dispersal barrier, but it is also likely to be important in limiting recruitment success of larvae that have reached the warm-temperate province either during periods when there is no inshore upwelling (about 55% of the time at Port Alfred and about 95% of the time at Port St Johns ) or that have bypassed the upwelling region while temporarily entrained in the Agulhas Current. Temperatures below 20°C are not uncommon on the south coast even during the summer months [27, 28], and the longer duration of development for subtropical mudprawn larvae increases the chance of recruitment failure.
Other adaptive features having a genetic basis that were not studied here may also play a role in limiting the establishment of the subtropical lineage in the warm-temperate province. Firstly, it is possible that adult mudprawns are also less able to tolerate low water temperatures, and individuals that have established themselves in the warm-temperate province may not be able to survive for long. This concept is illustrated by fiddler crab (Uca sp.) populations that have their southern limit around the middle of the Wild Coast, establishing populations in estuaries in the eastern portion of the warm-temperate biogeographic province. These populations eventually become extinct during the cold winter months (Newman, pers. obs.). Secondly, the embryos of the subtropical mudprawns that have settled in the warm-temperate province may be affected by the cooler, more saline conditions of those estuaries .
Barber et al.  suggested that antagonistic behaviour among potential cryptic species of mantis shrimps, Haptosquilla pulchella, may contribute to maintaining sharp genetic breaks between regional lineages. Mudprawns are comparatively more limited in terms of social interactions because they are rarely encountered outside their burrows. Moreover, a few exceptions to the trend that individuals from the subtropical mtDNA lineage have the subdistal spine on the first pereiopod, while the temperate lineage individuals do not, may indicate that the two lineages are able to interbreed. This should be tested further using nuclear genetic markers, and if hybrids are identified, then this could be followed by experiments on hybrid fitness. Despite the likely absence of antagonistic behaviour among the lineages, it is nonetheless possible that one lineage competitively excludes the other from its habitat, as mudbanks having the ideal consistency for mudprawn settlement are often extremely crowded. For example, a mean density of 218 ± 166 individuals m-2 has been estimated for the warm-temperate Swartkops Estuary . If recruitment success of the one mudprawn lineage in its region of occupancy is considerably greater than that of migrants from the other lineage, then this will limit the frequency of occurrence of sympatric distribution patterns.
Although no evidence was found for larval retention, it must be acknowledged that mtDNA sequence data have limitations. Recent studies of high-dispersal coastal invertebrates using microsatellites, for example, have identified subtle genetic structure or high autocorrelation values at local scales where none were identified with mtDNA [43, 44]. Nonetheless, mtDNA data often provide sufficient signal to detect genetic structure, even in planktonic dispersers within biogeographic regions [45, 46]. The lack of genetic structure and isolation by distance found among populations of the subtropical mudprawn lineage indicates that it is unlikely that larval retention is important in maintaining genetic breaks when compared to environmental and adaptive factors.
While there is no clear evidence that larvae of U. africana have behavioural strategies that allow them to avoid displacement from suitable habitat, a certain amount of larval retention may nonetheless be possible because of the regional oceanography coupled with a fairly short planktonic development phase. This may contribute to the number of migrants successfully crossing the ACDZ being too low to establish viable populations. Both hypotheses are discussed below.
Inshore current flow
The finding that gene flow is bidirectional along the Wild Coast suggests that a large proportion of larvae remain inshore of the Agulhas Current, where wind-driven longshore currents are generally weak and frequently change direction. This suggests that even in the absence of significant isolation by distance (because levels of gene flow are high at this relatively small scale), much larval settlement may nonetheless occur in the vicinity of the parent habitat. McQuaid and Phillips  found that 90% of mussel larvae that were dispersed by wind-driven currents on the south coast settled within less than 5 km of the parent habitat, which they attributed to frequent wind reversals. Similarly, in a study on dispersal potential of mangrove propagules by ocean currents using plastic drift cards, Steinke & Ward  found that twice as many cards released near Wavecrest were washed ashore in the vicinity than cards released from two more northerly sites (Richards Bay and Durban). This suggests that the oceanography of this region is not conducive to long-distance dispersal. Because the optimum temperature for larval development is the same for the two lineages, adaptive differences are unlikely to play a major role in excluding the warm-temperate lineage from the east coast, and oceanography is a more probable explanation.
Short development phase
The larval ontogeny of both lineages comprised three zoeal stages prior to the metamorphic moult. For upogebiid prawns, this is a slightly abbreviated mode of development, since many species pass through four zoeal stages prior to metamorphosis [48, 49]. The primary implication of abbreviated development is that larvae are able to complete development within a shorter period compared to sibling species, reducing their exposure to predators and other stressors and, in the context of this study, limiting wide-scale dispersal.
This study illustrates the value of incorporating information on regional adaptation and morphological differentiation as indicators of evolutionary divergence when studying phylogeographic breaks, in addition to the more commonly applied approach of investigating correlations between phylogeographic patterns and oceanographic features. Warm-temperate and subtropical mtDNA lineages of the mudprawn Upogebia africana can not only be distinguished morphologically, but the former was found to be able to complete larval development at low temperatures at which the latter can not. This indicates that the structuring effects of environmental discontinuities and currents not only limit gene flow between the two biogeographic regions, but that the establishment of each lineage in the habitat of the other may be limited by physiological adaptations to regional environmental conditions (possibly including additional ones not studied here). While the work presented here indicates incipient speciation of the two lineages, much work still remains to be done to determine whether the lineages are reproductively isolated, whether there is competitive exclusion, and how embryos and adults are affected by environmental differences between the biogeographic provinces.
Multi-locus and multi-species studies may eventually elucidate how phylogeographic breaks between evolutionary lineages of marine invertebrates in south-eastern Africa became established in the first place. Adaptations to environmental conditions can only arise when units of populations of a species are geographically isolated from other such units, but there are no obvious geological features along the shallow continental shelf in this region that could have acted as long-term vicariant boundaries. A rough divergence time estimate based on mtDNA data indicates that the two mudprawn lineages diverged around 3 million years ago , i.e. during a time when the present flow path of the Agulhas Current had been in place for over 2 million years . This suggests that evolutionary divergence of regional lineages may have been driven by oceanographic features similar to the ones observed today, and that the lineages were never completely isolated from each other. Support for parapatric speciation scenarios among populations of marine organisms is still comparatively rare, but may be more common than is generally acknowledged [51, 52].
Tissue samples from 35 mudprawns each were collected in four Wild Coast estuaries by either breaking off a specimen's right first pereipod (large specimens) or preserving complete individuals (small specimens) in preservation medium containing 70% ethanol and 30% TE buffer. DNA extraction, amplification of a portion of the COI gene and sequencing were performed as described previously using crustacean forward primer 5'-TCA ACA AAT CAY AAA GAY ATT GG-3' and decapod reverse primer 5'-AAT TAA AAT RTA WAC TTC TGG-3' . The sequences generated in this study were submitted to GenBank (accession numbers FJ416903–FJ417042). The program ARLEQUIN 3.1  was used to calculate each population's haplotype and nucleotide diversity based on pairwise differences, and an exact test of sample differentiation based on haplotype frequencies  was performed for pairs of populations. To visualise which lineages were represented in each of the estuaries sampled, a median-joining haplotype network was reconstructed using NETWORK 18.104.22.168 . We used the program LAMARC 2.02  to estimate migration rates (M) among populations in conjunction with each population's parameter θ (effective female population size multiplied by mutation rate). A Bayesian Metropolis Hastings Markov chain Monte Carlo search strategy was specified that consisted of a single final chain with 60 000 samples, a sampling interval of 1000, and a burnin of 6 000 000 steps. As MODELTEST 3.7  identified the HKY model  as being most suitable for the COI data, we specified a transition/transversion ratio of 5.54 as determined in MEGA 4 , as well as empirical base frequencies. Coalescent-based programs perform best when sample sizes are low , so we limited our analyses to 20 individuals per population. However, in order to incorporate as much information as possible from the original data-sets of 35 individuals per population, we randomly sub-sampled each population three times to create different data-sets, each consisting of a different sample of 20 individuals per population. Values of θ and M, as well as their standard deviations and 95% confidence intervals, were then reported as the mean values of the three runs. Each run was repeated three times to check for consistency of results, and TRACER v. 1.4 was used to assess whether the program had run for sufficiently long, as indicated by effective sample sizes of no less than 200, and trendlines with stable likelihood values. A Mantel test in ARLEQUIN was used to test for isolation by distance (IBD). Most southern African coastal invertebrates with planktonic dispersal phases comprise mtDNA phylogroups whose distributions are limited by marine biogeographic disjunctions . Pooling sequence data from different regional lineages will therefore incorrectly result in the detection of significant correlations between genetic and geographic distance [e.g., ], even in species whose phylogroups exhibit panmixia within their respective biogeographic region. As a minimum of three populations are required to test for IBD, we used only haplotypes of the subtropical mudprawn lineage. One matrix comprised genetic distances (ΦST values based on pairwise differences) and the other geographic distances (shortest long-shore distances among estuarine populations). To test for significant correlation between matrices, we specified 10 000 random permutations.
Larval rearing experiments
Ovigerous female mudprawns with embryos in an advanced state of development (i.e. with well-developed eyespots) were collected from the muddy intertidal of the lower reaches of the Swartkops and Mngazana estuaries (Fig. 1). In the laboratory, females were held in filtered (0.5 μm), UV-irradiated seawater (salinity of 35, hereafter referred to as culture water) at 20°C under an artificial 12:12 h light:dark (06h00 to 18h00) photoperiod until larvae hatched, which always occurred at night. Actively swimming larvae were collected between 06h00–08h00 the morning after hatching using a wide bore pipette and individually transferred to glass rearing vials filled with 15 ml of culture water at 20°C and fed newly hatched (<8 h old) Artemia sp. nauplii. Thereafter, larvae were acclimated to rearing temperatures (12, 14, 17, 23, 26 and 29°C) at a rate of 3°C per hour, by transferring rearing vials in trays between water baths set at the appropriate temperature. After acclimation, three broods of 20 individually reared larvae from the warm-temperate lineage, and three broods of 10 larvae from the subtropical lineage, were reared at each temperature (females from the subtropical lineage were generally smaller and carried fewer eggs).
Daily, larvae were removed from the rearing vials, and examined under a dissecting microscope to determine survival and development stage. Larvae were then transferred to new rearing vials filled with culture water at the appropriate temperature, and fed. Daily analysis continued until all larvae had either metamorphosed to the Decapodid stage or died. Larvae were considered dead when opaque and/or no movement of internal or external structures and appendages could be detected under moderate magnification.
Plots of survival and duration of development data were presented as the arithmetic mean ± standard deviation (Figs 4 and 5). Survival and duration of development for larvae from females of each lineage at each temperature treatment were compared using t-tests (after survival proportions were arc-sine transformed). Quantitative relationships between stage-specific and cumulative duration of development and temperature were described by means of linear and non-linear (power function) regressions. Slopes of regressions were tested for significant differences from zero with t-tests, after linearisation of power models by log transformation of both independent and dependent variables. Equality of regression slopes and intercepts was analysed using a procedure equivalent to ANCOVA.
We are grateful to Magda Hawkins for drawing the map and to Jo Palmer for surveying the subdistal spine on the dorsal margin of the merus. PRT was supported by a National Research Foundation (NRF) Postdoctoral Research Fellowship and a study grant from the Ernest Oppenheimer Memorial Trust. This study was funded by a National Research Foundation Grant (GUN 2069119) awarded to NPB.
- Lessios HA, Kessing BD, Robertson DR: Massive gene flow across the world's most potent marine biogeographic barrier. Proc Roy Soc B. 1998, 265: 583-588. 10.1098/rspb.1998.0334.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Scheltema RS: Initial evidence for the transport of teleplanic larvae of benthic invertebrates across the East Pacific Barrier. Biol Bull. 1988, 174: 145-152. 10.2307/1541781.View ArticleGoogle Scholar
- Palumbi SR: Marine reserves and ocean neighbourhoods: the spatial scale of marine populations and their management. Annu Rev Environ Resour. 2004, 29: 31-68. 10.1146/annurev.energy.29.062403.102254.View ArticleGoogle Scholar
- Kinlan BP, Gaines SD: Propagule dispersal in marine and terrestrial environments: a community perspective. Ecology. 2003, 84: 2007-2020. 10.1890/01-0622.View ArticleGoogle Scholar
- Barber PH, Palumbi SR, Erdmann VM, Moosa MK: Sharp genetic breaks among populations of Haptosquilla pulchella (Stomatopoda) indicate limits to larval transport: patterns, causes and consequences. Mol Ecol. 2002, 11: 659-674. 10.1046/j.1365-294X.2002.01468.x.View ArticlePubMedGoogle Scholar
- Teske PR, Hamilton H, Palsbøll PJ, Choo CK, Gabr H, Lourie SA, Santos M, Sreepada A, Cherry MI, Matthee CA: Molecular evidence for long-distance colonization in an Indo-Pacific seahorse lineage. Mar Ecol Prog Ser. 2005, 286: 249-260. 10.3354/meps286249.View ArticleGoogle Scholar
- Waters JM, King TM, O'Loughlin PM, Spencer HG: Phylogeographic disjunction in abundant high-dispersal littoral gastropods. Mol Ecol. 2005, 14: 2789-2802. 10.1111/j.1365-294X.2005.02635.x.View ArticlePubMedGoogle Scholar
- Waters JM: Marine biogeographic disjunction in temperate Australia: historic landbridge, contemporary currents, or both?. Divers Distrib. 2008Google Scholar
- Lessios HA, Kane J, Robertson DR: Phylogeography of the pantropical sea urchin Tripneustes: contrasting patterns of population structure between oceans. Evolution. 2003, 57: 2026-2036.View ArticlePubMedGoogle Scholar
- Rocha LA, Robertson DR, Rocha CR, van Tassell JL, Craig MT, Bowen BW: Recent invasion of the tropical Atlantic by an Indo-Pacific coral reef fish. Mol Ecol. 2005, 14: 3921-3928. 10.1111/j.1365-294X.2005.02698.x.View ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- Teske PR, McQuaid CD, Froneman PW, Barker NP: Impacts of marine biogeographic boundaries on phylogeographic patterns of three South African estuarine crustaceans. Mar Ecol Prog Ser. 2006, 314: 283-293. 10.3354/meps314283.View ArticleGoogle Scholar
- Teske PR, Papadopulos 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.View ArticleGoogle Scholar
- Ridgeway TM, Stewart BA, Branch GM, Hodgson AN: Morphological and genetic differentiation of Patella granularis (Gastropoda: Patellidae): recognition of two sibling species along the coast of southern Africa. J Zool. 1998, 245: 317-333. 10.1111/j.1469-7998.1998.tb00107.x.View ArticleGoogle Scholar
- Edkins MT, Teske PR, Papadopoulos I, Griffiths CL: Morphological and genetic data suggest that southern African crown crabs, Hymenosoma orbiculare, represent five distinct species. Crustaceana. 2007, 80: 667-683. 10.1163/156854007781360694.View ArticleGoogle Scholar
- 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.View ArticleGoogle Scholar
- Gopal K, Tolley KA, Groeneveld JC, Matthee CA: Mitochondrial DNA variation in spiny lobster Palinurus delagoae suggests genetically structured populations in the southwestern Indian Ocean. Mar Ecol Prog Ser. 2006, 319: 191-198. 10.3354/meps319191.View ArticleGoogle Scholar
- Teske PR, Barker NP, McQuaid CD: Lack of genetic differentiation among four sympatric southeast African intertidal limpets (Siphonariidae): phenotypic plasticity in a single species?. J Moll Stud. 2007, 73: 223-228. 10.1093/mollus/eym012.View ArticleGoogle Scholar
- Tripp RT: An atlas of coastal surface drifts: Cape Town to Durban. 1967, South African Oceanographic Data Centre, Department of Oceanography, University of Cape TownGoogle Scholar
- Anderson FP, Grundlingh ML, Stavropoulos CC: Kinematics of the southern Natal coastal circulation: some historic measurements 1962–63. S Afr J Sci. 1988, 84: 857-860.Google Scholar
- Harris TFW: Sources of the Agulhas Current in the spring of 1964. Deep-Sea Res. 1972, 19: 633-650.Google Scholar
- Lutjeharms JRE, Cooper J, Roberts M: Upwelling at the inshore edge of the Agulhas Current. Cont Shelf Res. 2000, 20: 737-761. 10.1016/S0278-4343(99)00092-8.View ArticleGoogle Scholar
- Beckley LE: The ichthyoplankton assemblage of the Algoa Bay nearshore region in relation to coastal zone utilization by juvenile fish. S Afr J Zool. 1985, 21: 244-252.View ArticleGoogle Scholar
- Dubula O, Lasiak TA: Spatial variability in the mudprawn Upogebia africana on the south-east coast of South Africa. Afr J Mar Sci. 2003, 25: 275-282.View ArticleGoogle Scholar
- Hill BJ: The effect of heated effluent on egg production in the estuarine prawn Upogebia africana (Ortmann). J Exp Mar Biol Ecol. 1977, 29: 291-302. 10.1016/0022-0981(77)90072-7.View ArticleGoogle Scholar
- Schumann EH, Cohen AL, Jury MR: Coastal sea surface temperature variability along the south coast of South Africa and the relationship to regional and global climate. J Mar Res. 1995, 53: 231-248. 10.1357/0022240953213205.View ArticleGoogle Scholar
- Thibault-Botha D, Lutjeharms JRE, Gibbons MJ: Siphonophore assemblages along the east coast of South Africa; mesoscale distribution and temporal variations. J Plankton Res. 2004, 26: 1115-1128. 10.1093/plankt/fbh104.View ArticleGoogle Scholar
- Knowlton N, Keller BD: Larvae which fall far short of their potential: highly localized recruitment in an alpheid shrimp with extended larval development. Bull Mar Sci. 1986, 39: 213-223.Google Scholar
- Taylor MS, Hellberg ME: Genetic evidence for local retention of pelagic larvae in a Caribbean Reef Fish. Science. 2003, 299: 107-109. 10.1126/science.1079365.View ArticlePubMedGoogle Scholar
- Strathman RR: Evolution of local recruitment and its consequences for marine populations. Bull Mar Sci. 2002, 70: 377-396.Google Scholar
- Queiroga H, Blanton J: Interactions between behaviour and physical forcing in the control of horizontal transport of decapod crustacean larvae. Adv Mar Biol. 2004, 47: 107-214. 10.1016/S0065-2881(04)47002-3.View ArticleGoogle Scholar
- Marta-Almeida M, Dubert J, Peliz Á, Queiroga H: Influence of vertical migration pattern on retention of crab larvae in a seasonal upwelling system. Mar Ecol Prog Ser. 2006, 307: 1-19. 10.3354/meps307001.View ArticleGoogle Scholar
- Kuhner MK: LAMARC 2.0: Maximum likelihood and Bayesian estimation of population parameters. Bioinformatics. 2006, 22: 768-770. 10.1093/bioinformatics/btk051.View ArticlePubMedGoogle Scholar
- Rambaut A, Drummond AJ: TRACER v1.4. 2007, [http://beast.bio.ed.ac.uk/Tracer]Google Scholar
- Sotka EE, Wares JP, Barth JA, Grosberg RD, Palumbi SR: Strong genetic clines and geographical variation in gene flow in the rocky intertidal barnacle Balanus glandula. Mol Ecol. 2004, 13: 2143-2156. 10.1111/j.1365-294X.2004.02225.x.View ArticlePubMedGoogle Scholar
- Latta RG: Conservation genetics as applied evolution: from genetic pattern to evolutionary process. Evol Appl. 2008, 1: 84-94. 10.1111/j.1752-4571.2007.00008.x.PubMed CentralView ArticlePubMedGoogle Scholar
- 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.View ArticleGoogle Scholar
- Teske PR, Papadopoulos I, Newman BK, McQuaid CD, Barker NP: Climate change, genetics or human choice: why were the shells of mankind's earliest ornament larger in the Pleistocene than in the Holocene?. PLoS ONE. 2007, 2: e614-10.1371/journal.pone.0000614.PubMed CentralView ArticlePubMedGoogle Scholar
- McQuaid CD, Phillips TE: Limited wind-driven dispersal of mussel larvae: in situ evidence from the plankton and the spread of the invasive species Mytilus galloprovincialis in South Africa. Mar Ecol Prog Ser. 2000, 201: 211-220. 10.3354/meps201211.View ArticleGoogle Scholar
- Harrison TD: Physico-chemical characteristics of South African estuaries in relation to the zoogeography of the region. Estuar Coast Shelf Sci. 2004, 61: 73-87. 10.1016/j.ecss.2004.04.005.View ArticleGoogle Scholar
- Hanekom N, Baird D, Erasmus T: A quantitative study to assess standing biomasses of macrobenthos in soft substrata of the Zwartkops estuary, South Africa. S Afr J Mar Sci. 1988, 6: 163-174.View ArticleGoogle Scholar
- Banks SC, Piggott MP, Williamson JE, Bové U, Holbrook NJ, Beheregaray LB: Oceanic variability and coastal topography shape genetic structure in a long-dispersing sea urchin. Ecology. 2007, 88: 3055-3064. 10.1890/07-0091.1.View ArticlePubMedGoogle Scholar
- York KL, Blacket MJ, Appleton BR: The Bassian Isthmus and the major ocean currents of southeast Australia influence the phylogeography and population structure of a southern Australian intertidal barnacle Catomerus polymerus (Darwin). Mol Ecol. 2008, 17: 1948-1961. 10.1111/j.1365-294X.2008.03735.x.View ArticlePubMedGoogle Scholar
- Nicastro KR, Zardi GI, McQuaid CD, Teske PR, Barker NP: Mesoscale topographically driven genetic structure in marine mussels. Mar Ecol Prog Ser. 2008, 368: 189-195. 10.3354/meps07607.View ArticleGoogle Scholar
- Matthee CA, Cockroft AC, Gopal K, Heyden von der 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.View ArticleGoogle Scholar
- Steinke TD, Ward CJ: Use of plastic drift cards as indicators of possible dispersal potential of propagules of the mangrove Avicennia marina by ocean currents. Afr J Mar Sci. 2003, 25: 169-176.View ArticleGoogle Scholar
- Konishi K: Larval development of the mud shrimp Upogebia (Upogebia) major (De Haan) under laboratory conditions, with comment on larval characters of thalassinid families. Bull Nat Res Inst Aquacult. 1989, 15: 1-17.Google Scholar
- Newman BK: Life-history strategies of the mudprawn, Upogebia africana Ortmann 1894 (Thalassinidea: Decapoda: Crustacea), during the marine phase of development. 2000, PhD thesis, University of Port ElizabethGoogle Scholar
- Martin AK, Goodlad SW, Salmon DA: Sedimentary basin in-fill in the northernmost Natal Valley, hiatus development and Agulhas Current palaeo-oceanography. J Geol Soc. 1982, 139: 183-201. 10.1144/gsjgs.139.2.0183.View ArticleGoogle Scholar
- Beheregaray LB, Sunnucks P: Fine-scale genetic structure, estuarine colonization and incipient speciation in the marine silverside fish Odontesthes argentiniensis. Mol Ecol. 2001, 10: 2849-2866. 10.1046/j.1365-294X.2001.t01-1-01406.x.View ArticlePubMedGoogle Scholar
- Rocha LA, Robertson DS, Roman J, Bowen BW: Ecological speciation in tropical reef fishes. Proc Biol Sci. 2005, 272 (1563): 573-579.PubMed CentralView ArticlePubMedGoogle Scholar
- Excoffier L, Laval G, Schneider S: ARLEQUIN (version 3.0): An integrated software package for population genetics data analysis. Evol Bioinf Online. 2005, 1: 47-50.Google Scholar
- Raymond M, Rousset F: An exact test for population differentiation. Evolution. 1995, 49: 1280-1283. 10.2307/2410454.View ArticleGoogle Scholar
- Bandelt H-J, Forster P, Röhl A: Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999, 16: 37-48.View ArticlePubMedGoogle Scholar
- Posada D, Crandall KA: MODELTEST: testing the model of DNA substitution. Bioinformatics. 1998, 14: 817-818. 10.1093/bioinformatics/14.9.817.View ArticlePubMedGoogle Scholar
- Hasegawa M, Kishino K, Yano T: Dating the human-ape splitting by a molecular clock of mitochondrial DNA. J Mol Evol. 1985, 22: 160-174. 10.1007/BF02101694.View ArticlePubMedGoogle Scholar
- Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24: 1596-1599. 10.1093/molbev/msm092.View ArticlePubMedGoogle Scholar
- Felsenstein J: Accuracy of coalescent likelihood estimates: do we need more sites, more sequences or more loci?. Mol Biol Evol. 2006, 23: 691-700. 10.1093/molbev/msj079.View ArticlePubMedGoogle Scholar
- Lutjeharms JRE: Three decades of research on the greater Agulhas Current. Ocean Sci Discuss. 2006, 3: 939-995.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.