Oceanic dispersal barriers, adaptation and larval retention: an interdisciplinary assessment of potential factors maintaining a phylogeographic break between sister lineages of an African prawn
BMC Evolutionary Biology volume 8, Article number: 341 (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.
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 .
The South African coastline provides several examples of phylogeographic breaks associated with discontinuities in water temperature. These include breaks at Cape Agulhas and Cape Point, both of which separate lineages associated with cool-temperate and warm-temperate biogeographic regions [12–14], a break in the Wild Coast region on the southeast coast that is associated with the boundary of the warm-temperate and subtropical provinces [13, 15–17], and a break in the northeast of the country near St Lucia that separates subtropical and tropical lineages [18, 19]. While phylogeographic breaks associated with these biogeographic discontinuities are mostly subtle and are based on significant genetic structure or monophyletic mtDNA phylogroups, the Wild Coast case includes examples of lineages for which morphological [15, 16] and physiological (Zardi & Nicastro, in prep.) differences have been identified. The oceanography of the Wild Coast is strongly influenced by the warm, southward-flowing Agulhas Current, which lies about 10 km offshore along most of South Africa's east and southeast coast as it follows the narrow continental shelf . Inshore of the current, there are frequent reversals of current direction, with surface water flow being driven predominantly by longshore wind . From about Port Alfred southwards (Fig. 1), the shelf starts to widen and the current is deflected away from the coast (a region hereafter referred to as the Agulhas Current Deflection Zone, or ACDZ). This results in a decreasing influence of the current that manifests itself in sea surface temperatures being cooler south of this point . A persistent, localised upwelling cell is located inshore of the current where the shelf widens, although it does not always affect surface water temperature . Many tropical and subtropical species have been reported in temperate waters beyond the ACDZ, especially during the summer months, but they do not establish themselves permanently .
The objective of this study was to assess the relative importance of various oceanographic and biological factors in maintaining the separation of genetic lineages of coastal invertebrates at the Wild Coast phylogeographic break. As the study organism, we selected the estuarine mudprawn Upogebia africana (Ortmann, 1894), because of all the species that have phylogeographic breaks in this region, its biology has been comparatively well studied. Based on mtDNA sequence data , the species can be divided into two major lineages, one associated with temperate biogeographic provinces (with some differentiation between cool-temperate and warm-temperate provinces) and the other subtropical. In previous surveys [13, 14], the distribution of the temperate lineage was found to range from the Olifants Estuary in south-western Africa to the Mthatha Estuary on the southeast coast, and that of the subtropical lineage from Haga Haga to Inhambane in Mozambique, with some overlap in the Wild Coast region (Fig. 1). The two lineages can be distinguished on the basis of a single morphological character: the subtropical lineage has a subdistal spine on the first pereiopod that is absent in the temperate lineage (Fig. 2). Dubula and Lasiak  further reported that individuals from the central Wild Coast were significantly smaller than those from the southern Wild Coast. The size difference observed may have a genetic basis rather than being merely a function of water temperature, because individuals of the temperate lineage grow larger in warmer water .
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.
Haplotypes belonging to the subtropical lineage were exclusively found in the two northernmost populations (Mthatha and Mngazana), and the southernmost population (Gqunube) had only haplotypes from the temperate lineage. The geographically intermediate Wavecrest population comprised a mixture of haplotypes from both genetic lineages (Fig. 3). While haplotype diversity indices were of similar magnitude for all populations (Gqunube 0.94 ± 0.03 [mean ± SD]; Wavecrest: 0.95 ± 0.02; Mthatha: 0.95 ± 0.02; Mngazana: 0.90 ± 0.05), the nucleotide diversity index of the Wavecrest population was considerably higher than those of the others, indicating that the haplotypes present in this population represent a combination of haplotypes from different lineages that are distantly related to each other (Gqunube: 0.009 ± 0.005; Wavecrest: 0.032 ± 0.016; Mthatha: 0.007 ± 0.004; Mngazana: 0.006 ± 0.004).
To test whether gene flow in the Wild Coast mostly takes place in a south-westerly direction (which would indicate that most dispersal of mudprawn larvae takes place while these are temporarily entrained in the Agulhas Current), we used the program LAMARC 2.02  to simultaneously calculate θ (effective female population size multiplied by mutation rate) and M (migration rate). As this program performs best when populations have significantly different genetic structure, we performed exact tests of sample differentiation based on haplotype frequencies. This test identified significant differentiation (P < 0.05) in all pairwise population comparisons, with the exception of the Mthatha vs. Mngazana comparison (P = 0.62). As the Mthatha Estuary is geographically closer to the remaining two estuaries, we excluded the Mngazana population from the LAMARC analyses. Estimates of θ and M for runs using the same three sub-sampled data-sets were consistent, trace plots generated by LAMARC using TRACER indicated that log-likelihood values had reached stationarity following burnin, and the lowest effective sample size found during any run was 388. This indicates that the LAMARC analyses were run for long enough to obtain reliable estimates. Migration rates close to zero were found between the Gqunube and Mthatha populations, in which only temperate and subtropical haplotypes were present, respectively (Table 1). Estimates for the remaining population pairs ranged from 67 for gene flow from Mthatha to Wavecrest to 407 from Wavecrest to Gqunube. Mean estimates of θ were 0.02 for the Gqunube and Mthatha populations, and 0.03 for the Wavecrest population. Given the large standard deviations and confidence intervals, as well as the fact that estimates were based on a single marker, we consider our results sufficient to indicate that gene flow along the Wild Coast is not unidirectional, but insufficient to obtain any more accurate estimates.
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.
Upogebia africana larvae of the warm-temperate lineage completed development to the Decapodid at all temperatures (Fig. 4). In contrast, larvae of the subtropical lineage were unable to complete development to the Decapodid at 12 and 14°C (Fig. 4), all Zoea 1 larvae dying at 12°C and those that survived at 14°C subsequently dying during the Zoea 2 stage. Apart from this difference in temperature tolerance, absolute survival at various temperatures was comparable. Larval survival for both lineages was highest at 23°C, decreasing progressively at either side of this temperature. Survival to the Decapodid for the warm-temperate lineage was significantly higher at 26°C, but there were no differences at other temperatures in which larvae from both lineages survived to the Decapodid.
For both lineages the mean duration of development for each zoeal stage and, consequently, also the mean cumulative duration of development decreased sharply with increasing temperature until a threshold temperature was reached, after which duration remained constant or more or less constant (Fig. 5). The threshold temperature for the Zoea 1 was 20°C and for the Zoea 2 and 3 it was 23°C. Each incremental decrease in temperature below these threshold temperatures resulted in a highly significant increase in duration of development. The temperature dependence of duration of development for the Zoea 1 of both lineages at temperatures below the threshold could be described by statistically highly significant linear regressions and for the Zoea 2 and 3 and cumulative duration by power-model regressions (Fig. 5). Slopes of regressions for each zoeal stage and for cumulative duration were significantly steeper for the subtropical lineage (in all cases p < 0.001).
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 184.108.40.206 . 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.
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.
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.
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.
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.
Kinlan BP, Gaines SD: Propagule dispersal in marine and terrestrial environments: a community perspective. Ecology. 2003, 84: 2007-2020. 10.1890/01-0622.
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.
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.
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.
Waters JM: Marine biogeographic disjunction in temperate Australia: historic landbridge, contemporary currents, or both?. Divers Distrib. 2008
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.
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.
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, 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.
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.
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.
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.
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.
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.
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.
Tripp RT: An atlas of coastal surface drifts: Cape Town to Durban. 1967, South African Oceanographic Data Centre, Department of Oceanography, University of Cape Town
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.
Harris TFW: Sources of the Agulhas Current in the spring of 1964. Deep-Sea Res. 1972, 19: 633-650.
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.
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.
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.
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.
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.
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.
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.
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.
Strathman RR: Evolution of local recruitment and its consequences for marine populations. Bull Mar Sci. 2002, 70: 377-396.
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.
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.
Kuhner MK: LAMARC 2.0: Maximum likelihood and Bayesian estimation of population parameters. Bioinformatics. 2006, 22: 768-770. 10.1093/bioinformatics/btk051.
Rambaut A, Drummond AJ: TRACER v1.4. 2007, [http://beast.bio.ed.ac.uk/Tracer]
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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.
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 Elizabeth
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.
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.
Rocha LA, Robertson DS, Roman J, Bowen BW: Ecological speciation in tropical reef fishes. Proc Biol Sci. 2005, 272 (1563): 573-579.
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.
Raymond M, Rousset F: An exact test for population differentiation. Evolution. 1995, 49: 1280-1283. 10.2307/2410454.
Bandelt H-J, Forster P, Röhl A: Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999, 16: 37-48.
Posada D, Crandall KA: MODELTEST: testing the model of DNA substitution. Bioinformatics. 1998, 14: 817-818. 10.1093/bioinformatics/14.9.817.
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.
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.
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.
Lutjeharms JRE: Three decades of research on the greater Agulhas Current. Ocean Sci Discuss. 2006, 3: 939-995.
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.
PRT conceived of the study and coordinated the research, collected samples, generated and analysed the genetic data, helped with the larval rearing experiments, and drafted the manuscript. IP collected samples, did larval rearing experiments and helped analyse the genetic data. BKN analysed the larval rearing data and helped draft the manuscript. PCD identified the subdistal spine used to distinguish the two Upogebia africana lineages and prepared the drawings of the pereiopods. NPB and CDM provided conceptual guidance and logistical support. All authors contributed to the preparation of the manuscript, and read and approved the final version.
About this article
Cite this article
Teske, P.R., Papadopoulos, I., Newman, B.K. et al. Oceanic dispersal barriers, adaptation and larval retention: an interdisciplinary assessment of potential factors maintaining a phylogeographic break between sister lineages of an African prawn. BMC Evol Biol 8, 341 (2008). https://doi.org/10.1186/1471-2148-8-341