BACKGROUND
the early recognition of functionally significant levels of biodiversity within a species is a central objective of evolutionary biology and biodiversity conservation  <cit> .
this is particularly important at times when biodiversity loss is one of the planet´s major global problems  <cit> .
biodiversity can be investigated, and its importance evaluated, at different levels including species, their genes and functional traits  <cit> .
given the recurrent detection of distinct phylogeographic lineages, it is clearly necessary to consider the adaptive significance of intraspecific differences in order to understand the evolutionary potential of a species  <cit> .
complementary approaches based on neutral genetic markers, commonly employed to estimate recent demographic connectivity, and on methods that can reveal ecologically relevant adaptive traits across heterogeneous habitats are pivotal for the assessment and management of intraspecific diversity .
the identification of evolutionarily significant units uniquely based on neutral molecular markers, while overlooking intraspecific patterns of phenotypic divergence, can lead to underrepresentation of locally adapted populations and thus underestimate biodiversity  <cit> .
local environmental conditions resulting in spatially divergent selective regimes can maintain local adaptation across populations of the same species independently of varying and/or contrasting neutral genetic structure  <cit> .
species with wide distributions over variable environments and large effective population sizes often display marked intraspecific biodiversity with distinct genotypes/populations within a species, resulting from adaptation to local environmental conditions but capable of interbreeding with other ecotypes of the same species .
marine broadcast spawners with broad geographic ranges extending over abrupt environmental clines are ideal model organisms for studying the interplay between local adaptation and genetic structure due to neutral processes, because divergent selection produced by distinct environments is coupled with a high potential for gene flow among geographically separated subpopulations through larval transport.
the southern african coastline covers a wide range of climatic and oceanic conditions and can be divided into three major biogeographic regions: cool-temperate , warm-temperate  and sub-tropical .
the environmental discontinuities characteristic of these different parts of the coastline are mirrored by a degree of endemism among the different biogeographic provinces and phylogeographic breaks in species that occur in more than one province .figure  <dig> study area and genetic results.
a) map of the study area and sampling sites.
black dots: localities either mentioned in the text or used for ecologic data; blue dots: localities used for genetic data abbreviated as in table  <dig>  continuous, black line: distribution of the western lineage ; dashed, black line: distribution of the eastern lineage .
the three biogeographic regions of southern africa  are delimited by grey dotted lines.
b) results of the bayesian approach using the discriminant analysis of principal components  scatter plot of individuals with a priori-defined locations, to investigate genetic structure.
labels were placed at the centre of dispersion of each group, which are delineated by inertia ellipses.
c) dapc a posteriori genetic cluster assignment for k =  <dig>  the size of each square is proportional to the percentage of individuals assigned to each cluster at each location.
d) stacked bar graph of assignment probabilities per individual within p. perna locations across the southern tip of south africa.
assignment proportions for each individual were obtained from structure, based on genotypes from seven nuclear microsatellite loci, for k =  <dig>
the brown mussel perna perna is a habitat-forming species distributed widely around the world in warm waters.
in southern africa, it dominates intertidal shores in the sub-tropical and warm-temperate bioregions .
it is absent from there to central namibia due to the cold, upwelled waters of the benguela system .
it then extends northwards along the west coast of africa to the southern iberian peninsula and into the mediterranean sea as far as the gulf of tunis .
analyses of mitochondrial dna  sequences indicate a sharp phylogeographic break on the southeast coast of south africa  forming a western and an eastern lineage  <cit> .
the western lineage includes animals from both the namibian coast and the south coast of south africa .
the eastern lineage includes individuals from the southeast and east coasts of south africa.
the two lineages overlap in their distributions on the southeast coast over a distance of approximately  <dig> km.
the presence of the two lineages was recently confirmed by nuclear  sequence data  <cit> .
most importantly, mitochondrial  and nuclear  markers clearly show a non-sister relationship for the two south african lineages  <cit> .
the most likely biogeographic scenario explaining the independent origins of the two south african lineages involves an indo-pacific origin for p. perna, dispersal into the mediterranean and atlantic through the tethys seaway, and recent secondary contact after southward expansion of the western and eastern south african lineages  <cit> .
manipulative experiments and the release of oceanographic drifters suggest that phylogeographic patterns of p. perna may be maintained by a combination of local conditions experienced in the different biogeographic regions, and by the isolating effect of regional oceanographic dynamics  <cit> .
here, we examine the hypotheses that contemporary demographic connectivity and functional traits of p. perna populations are influenced by the tropical/subtropical biogeographic transition.
first, we adopted a population genetic approach using microsatellite markers, which are unaffected by natural selection, to investigate differentiation and migration patterns of p. perna between populations inhabiting different bioregions.
the re-examination of genetic structure in south african p. perna with microsatellites and comparison with previously reported mitochondrial data is critical to assess potential mito-nuclear discordance and contemporary biogeographic conditions.
second, we used a common-garden design to compare fundamental behavioural traits known to affect survival along large  and small  scale environmental gradients.
in particular, we assessed shell gaping and attachment strength, key attributes affecting tolerance of the two major stresses in intertidal habitats, desiccation and wave action.
a common-garden design is necessary to distinguish between adaptive and plastic responses of population subjected to distinct selective pressures.
finally, we compared intertidal elevation  of the two lineages by using in situ temperature profiles to compare intertidal limits of populations from the two bioregions.
materials and methods
ethics statement
mussel collection were performed under permit number res2014/ <dig> issued by the department of agriculture forestry and fisheries to the department of zoology and entomology at rhodes university.
genetic data
dna extraction, amplification and genotyping
mussels  were collected between february  <dig> and may  <dig> at six locations in southern africa  in the middle portion of the species’ vertical distribution.
mantle tissue  was dissected from each individual, preserved in 96% ethanol and stored at -20°c.
total genomic dna extraction was performed using a standard proteinase k protocol adapted from  <cit> .
molecular genotyping of a total of  <dig> samples was carried out using a set of  <dig> microsatellite loci in three multiplex reactions, with subsequent separation of the pcr products using an abi prism 3130xl dna analyser  with genescan liz  <dig> as size standard  and following the procedure of coelho et al.
<cit> .
data were scored using peak-scanner v.  <dig>  software .table  <dig>
genetic diversity of each population
locations
codes
coordinates
n
h
e
h
o
f
is
Â
locations are depicted in the map of figure 1a and are ordered from west to east.
decimal coordinates ; n: sample sizes; he and ho: expected and observed heterozygosity; fis: inbreeding coefficient with significant values after correction for multiple test indicated with *; Â : allelic richness  ± standard deviation.
location-based analyses
scoring errors, large allele dropout and null alleles were checked employing the program microchecker  <cit>  and frequency of null alleles estimated based on the algorithm presented in brookfield  <cit> .
additionally, tests for hardy–weinberg equilibrium were conducted using genetix  <dig>  software  <cit>  for all locus-location combinations by assessing the inbreeding estimator wright’s fis and estimating significant values after  <dig>  permutations.
tests for linkage disequilibrium between all pairs of loci were performed according to the method of black and kraftsur  <cit>  implemented in genetix.
for each location, allele frequencies, observed heterozygosity  and expected heterozygosity , plus fis  were computed with genetix.
standardized allelic richness  was estimated by resampling  <dig> times to standardize to the smallest sample size  to account for differences between sampled locations, using standarich r package  <cit> .
pairwise location estimates of genetic differentiation plus 95% bootstrapped confidence intervals  were estimated using weir and cockerham’s  <cit>  fst estimator  and jost’s  <cit>  dst.
calculations were executed in the diversity r package  <cit> .
estimating admixture
to evaluate the extent of admixture between the two mitochondrial lineages using microsatellite data, we used structure  <dig>   <cit> .
clusters  varied from one to a maximum of seven, corresponding to the number of locations in the study data plus one.
twenty replicate runs  were performed under the admixture and the correlated allele frequency model.
the most probable value of k was inferred based on the ad hoc criteria l proposed by pritchard et al.
<cit>  and delta k by evanno et al.
<cit> .
similar replicate runs were grouped based on the symmetric similarity coefficient of > <dig>  using the fullseach algorithm in clumpp  <dig>   <cit>  and visualized with bar_plotter.rb  <dig>  implemented in ruby.
to complement the results of structure, we used discriminant analysis of principal components , which can be performed when uncovering population admixture in the absence of assumptions.
this method first transforms the data using principal components analysis, which ensures that the variables are not correlated and that the number of variables is smaller than the number of individuals.
we ran two dapc analyses in adegenet r package  <cit> , following jonker et al.
<cit> : a location assignment and a genetic cluster assignment.
we assigned each individual a priori to its location of origin and obtained for each one the probability of assignment to the correct sampling site.
this procedure maximizes the among-location variation and minimizes the within-location variation  <cit> .
using the find.clusters function , we ran successive k-means clustering of the individuals from k =  <dig> to k =  <dig>  and identified the best supported number of clusters through comparison with the bayesian information criterion  for the different values of k. we determined the number of clusters and assigned each individual to a genetic cluster without providing any a priori population assignment.
in that way, the grouping factor is the genetic cluster and not the sampling location, offering an unbiased interpretation of population structure.
to avoid unstable assignments of individuals to clusters, we retained  <dig> pcs , but used all discriminant functions, in a preliminary dapc run.
the results were then reiterated by the optim.a.score function with  <dig> simulations to determine the optimal number of pcs, and a final dapc was subsequently carried out with the optimal number of pcs.
estimating the directionality of gene flow
to compare different biogeographic hypotheses on migration rates of p. perna, we used the coalescent-based program migrate-n migrate version  <dig> .
<dig>  <cit> .
this approach assumes the wright-fisher model, where locations have a constant effective size through time, the rate of mutation is constant, and locations exchange migrants with constant rates per generation, but those rates can vary among locations.
we conducted the analyses on a random subsample of  <dig> individuals, with the dataset structured into three groups according to the genetic clusters recovered : west , centre  and east .
we tested four variations of the three-group  migration model: all directional migration , strict west to east migration , strict east to west migration  and from the centre to the margins .
testing the directionality of gene flow is justified because the dominant ocean current between the ocean basins, the agulhas current, runs westerly from the indian towards the atlantic ocean  and is thought to play a role in limiting marine dispersal in the opposite direction  <cit>  though coastal counter-currents also occur  <cit>  .
models  <dig> and  <dig> consider indistinct migration among groups and migration directionality towards the transition area between the two clusters respectively.
initial values were calculated using fst.
mutation rates were set to be constant among loci.
the migrate-n run parameters were calibrated on the full model for convergence of the markov chain monte carlo sampling method.
the prior distributions were uniform for mutation-scaled population size parameters theta , that are four times the product of the effective population size and the mutation rate, and mutation-scaled migration rates m, that is, migration rate scaled by the mutation rate, over the range of θ = 0– <dig> and m = 0– <dig>  these settings resulted in converged posterior distributions with a clear maximum for each estimate.
the three replicate bayesian runs consisted of one long chain with a total of  <dig> million states and genealogy changes after discarding the first  <dig>  genealogies as “burn-in” for each locus.
for all the analyses we used an adaptive heating scheme with  <dig> simultaneous chains using different acceptance ratios .
overall loci information was combined into a single estimate by bézier approximation of the thermodynamic scores as described by beerli and palczewski  <cit>  and we averaged the bézier score over three different runs and used this as input to evaluate multiple models using bayes factors  <cit> .
ecological data
attachment strength and gaping behaviour
all mussels used for the assessment of behavioural traits were collected from locations known to support pure mtdna lineages .
specimens  were collected in february  <dig> .
common garden experiments involve the acclimation and comparison of individuals from distinct geographical locations under identical environmental conditions.
such experimental set-ups are widely used to disentangle the effects of genetic and environmental conditions on observed phenotype differences .
mussels were divided in three separate tanks , where they were placed at least 10cm apart so that they maintained a solitary position  and attachment strengths were not influenced by nearby conspecifics.
after three days, the attachment strength of each mussel was measured by drilling a hole into one of the valves and using a spring balance to record the force required to dislodge the mussel perpendicularly to the substratum  <cit> .
mussels  were exposed to air for 3h in a controlled environment chamber at 35°c and 60-80% humidity.
for each one hour interval, the number of valve movements  were noted by visual observation, without recording the width of the gape  <cit> .
intertidal distributional limits
temperature dataloggers  were used to relate the effective shore height of p. perna at two open coast sites on the east coast  and two within the western range of the species .
at each site, two transects were established, approximately  <dig> m apart and running perpendicular to the shoreline.
along each transect, two logger were deployed at the upper and two at the lower limit of p. perna distribution.
loggers were programmed to record data at  <dig> min intervals between february 1st to 14th of  <dig>  and march  <dig> to  <dig> of  <dig> , covering the two-week tidal cycle in both regions .
temperature logger profiles clearly reveal when loggers are first inundated by the returning tide .
this was used to calculate the effective shore height using the online forecasts service of shom database .
because tide cycles are not in phase across large geographical scales, durban and mossel bay tidal elevations were used for the east and south coast respectively.
where both loggers were recovered at a transect, the values were averaged and used to estimate populations’ intertidal limits at each transect.
data analyses
each trial was analysed separately.
attachment strength data were analysed using a 2-way anova with lineage  and block  as fixed and random factors respectively.
gaping behaviour data were analysed using a 2-way anova with lineage  and time  as fixed factors.
for each intertidal limit , data were analysed under a nested design with range  as a fixed factor, site  nested within range and transect nested within range and site.
heteroscedasticity was tested using levene’s test and post-hoc separation of significant means by snk tests.
analyses were performed in spss .
RESULTS
genetic data
null alleles at high frequencies  were indicated by microchecker for three loci  for several localities .
for these three loci, departures from hwe were consistently reported across all localities.
to avoid overestimation of fis and fst estimators  <cit>  by null alleles, these three loci were removed from all subsequent analyses, including descriptive analyses reported in table  <dig>  particularly high null allele frequencies have often been reported in several invertebrate taxa, including corals and bivalves, due to high levels of dna sequence variation  <cit> .
specifically, nucleotide variations in microsatellite flanking regions can cause a high frequency of null-alleles, even within populations used for characterization of microsatellites  <cit> .
all microsatellite loci were highly polymorphic among the six populations.
the  <dig> individuals contained  <dig> alleles in these loci, with an average number of alleles per locus of  <dig>  ranging from  <dig>  to  <dig> .
allelic richness and gene diversity he varied between  <dig> - <dig>  and  <dig> - <dig>  respectively and the highest allelic richness was recorded at du .
mean heterozygote deficiency was detected at all locations with the exception of po .
genetic differentiation among populations was detected as indicated by pairwise fst and dst values .
while dst pairwise differentiation was lower within than among two groups , fst values did not clearly differentiate pa as belonging to either group.table  <dig>
genetic differentiation between pairs of populations
location
hm
pl
pa
pe
du
po
<dig>
<dig>
<dig>
<dig>
<dig>
<dig>
<dig>
<dig>
<dig>
<dig>
<dig>
<dig>
<dig>
<dig>
<dig>
<dig>
<dig>
<dig>
<dig>
<dig>
<dig>
<dig>
<dig>
<dig>
<dig>
<dig>
codes correspond to locations in figure 1a and are ordered from west to east.
pairwise genetic differentiation .
95% confidence intervals are shown in brackets and confidence intervals not enclosing  <dig> are represented in bold.
discriminant analysis of principal components , run with geographic location defined a priori, clearly differentiated the ‘western’ and ‘eastern’ groups on the scatter plot x-axis, confirming the presence of the genetic structure indicated by the fst analysis .
individuals from pa contained an admixture of alleles from both groups and indistinctly clustered with each of them.
when the a posteriori genetic cluster assignment was run for k =  <dig> , most individuals of the ‘western’ group were assigned to the same cluster .
the individuals of the ‘eastern’ group were also mostly assigned to one cluster , and in a higher rate than the populations of the first cluster.
this is particularly true for pe, where all individuals were assigned to cluster  <dig>  individuals from pa were assigned to both clusters with a slightly higher number of individuals assigned to cluster  <dig>
this pattern of subdivisions was largely consistent with the results obtained with structure, which revealed the most significant increase of Δk at two clusters, , thereafter Δk remained unchanged.
the l plot rejected the option of k =  <dig>  we concluded that k =  <dig> is the most likely number of genetic clusters .
while at most locations individuals displayed an admixture higher than  <dig>  for either group or cluster, pa showed a split of  <dig>  and  <dig>  for the western and eastern groups respectively.
when comparing the four biogeographic hypotheses with migrate-n, model iii, strictly east to west migration, was the only model supported .table  <dig>
bayes factors model comparison of migration models for
p. perna
models
model parameters
bézier
dbézier
probability
set 1
set 2
set 3
model parameters code as follows: the first set of three signs indicate estimation of theta for west, migration from centre to west, migration from east to west; the second set indicate migration from west to centre, estimation of theta for centre, migration from east to centre; the third set indicates migration from west to east, migration from centre to east and estimation of theta for east.
an asterisk indicates that that particular parameter was estimated by the model, and a  <dig> indicates that no migration was allowed.
ecological data
in each trial, lineage  was a highly significant factor influencing attachment strength , with the eastern lineage having a higher attachment strength than the western.figure  <dig> ecological results.
a) gaping behaviour of the eastern and western lineages, each trial is plotted separately, b) strength of attachment , each trial plotted separately, c) high and low intertidal limits of p. perna western and eastern lineages.
dots depict mild outliers  and stars are extreme outliers .
in each trail, individuals from the eastern lineage always gaped significantly more than those from the western lineage .
post-hoc snk tests indicated that both lineages gaped significantly more at hour one than at hours two or three during the first trial , however, during the second trail, there was a time × lineage interaction  because at hour two eastern lineage individuals gaped more that at hour three.
p. perna in the western range reached a significantly higher intertidal elevation than on the east coast .
however, mussels in the western and eastern ranges displayed the same lower vertical limit.
discussion
the microsatellite data revealed strong genetic division between the two geographically defined groups of populations, supporting previous findings based on mitochondrial and nuclear markers  <cit> .
populations from cluster one inhabit the warm-temperate, south-western coast of south africa, while cluster two comprises populations distributed along the sub-tropical eastern shores.
most importantly, genetic admixture was only moderate, asymmetric and not restricted to the previously described  <dig> km contact zone, which lies mostly in the warm-temperate transition zone  <cit> ).
the admixture detected with microsatellite markers is probably the result of secondary contact between ancient lineages; phylogeographic patterns based on mtdna  and ndna  sequence data indicate an indo-pacific origin of the species and a dispersal into the mediterranean and atlantic through the tethys seaway.
after the closing of the seaway, the two non-sister lineages diverged  and came into secondary contact after their respective eastern, or western, southward expansion  <cit> .
however, genetic differentiation between lineages can be primarily attributed to the magnitude of the differences in allele frequencies rather than to unique alleles, as would be expected in ancient isolated lineages.
the detected gene exchange was primarily unidirectional, from east to west.
additionally, the sampled location at the geographic border of the two lineages  shows the strongest gene intermediacy.
gene exchange and the maintenance of this pattern could be attributed to several non-mutually exclusive effects.
eastward gene flow could be impeded by dispersal barriers , or by selection  causing outbreeding depression if lineages are locally adapted.
alternatively, poor mating success between lineages could also play a role, either due to divergent reproductive ecology, or to some level of genetic incompatibility.
barriers to dispersal
mussel larval dispersal is largely influenced by local hydrodynamics.
the main oceanographic influence on the east and south coasts of south africa is the agulhas current.
this is a major ocean current that flows towards the southeast, bringing warm, oligotrophic water from the mozambique channel  <cit> .
the inshore thermal front of this current varies geographically and can bring its warm waters very close inshore  <cit> .
from near pa southwards the current is deflected away from the coast.
the trajectories of nearshore drogues released off the south coast revealed no mixing of southern waters into the eastern water mass, while drifters released off the east coast moved southward, caught up in the agulhas current  <cit> .
although there are wind-driven surface currents that could allow limited eastward larval dispersal  <cit> , water exchange is predominantly unidirectional, promoting westerly larval dispersal of the eastern lineage and explaining the asymmetrical gene flow described here.
previous mtdna data on several species with planktonic larvae  support our results obtained with microsatellite data, by showing asymmetric east-to-west migration along the east coast .
interestingly, one of our most eastern locations  displayed a fair amount of genetic admixture.
du is the largest and busiest shipping terminal on the african continent, and globally one of the seven key epicentres for inter-regional exchange of species  <cit> .
shipping activities are known vectors of larval dispersal  <cit> .
the fact that du is the only location on the east coast with these characteristics suggests that the introduction of non-indigenous genes from the western range detected only here is a consequence of human-mediated larval transport.
departures from random mating
out of the six populations sampled, five displayed significant heterozygote deficiency.
null alleles are unlikely to be the cause of high fis values in this case, since microsatellite loci with a high frequency of null alleles were excluded from the analyses, while the remaining loci compiled in the final dataset displayed similar homozygote excess simultaneously across all loci.
inbreeding is also unlikely in a broadcast spawner with long a larval pelagic phase  as recruitment is expected to produce spatially random kinship.
one possible explanation for heterozygote deficiency is a wahlund effect due to the presence of genetically distinct subpopulations  in each population sampled.
in this context, it is interesting that pa, which is at the contact zone between the two lineages, displayed the lowest mean observed heterozygosity, although this was not true across all loci .
such effects can be expected when pooling discrete subpopulations with different allele frequencies that do not interbreed as a single randomly mating unit.
in contrast, non-significant levels of heterozygote deficiency were detected in po, which is the farthest population from the contact zone and located on the east side of the coast.
reproductive isolation of the lineages, leading to heterozygote deficiency, could result from differences in the timing and frequency of spawning events  <cit> .
it is not known when each lineage spawns nor what triggers the event, but lack of synchrony in spawning by co-existing mussel species has been reported, and analogous intraspecific differences may occur between lineages of p. perna  <cit> .
selection
behaviour can radically moderate an organism’s experience of the environment and so dictate physiological reactions, and consequently, tolerance limits to the stress imposed by those environmental conditions.
we show that the two p. perna mtdna lineages exhibit clear differences in behavioural traits, indicating an intraspecific behavioural divergence that potentially influences vertical zonation and habitat segregation over large and small spatial scales.
by acclimating and comparing individuals from different locations in a common laboratory environment  we could differentiate between adaptive and plastic responses to distinct environmental conditions.
the maintenance of behavioural phenotypic divergence in experimental populations indicated that differences are due to underlying genetic divergence.
mussel gaping, the intermittent opening and closure of the valves, allows the maintenance of aerobic respiration by sustaining an oxygen gradient across the mantle wall and the gills during emersion.
respiration rates of organisms usually increase with temperature and several studies have shown that oxygen consumption by mussels increases exponentially with temperature .
in p. perna, gaping significantly increased with temperature, supporting the important role of this behaviour in aerobic respiration  <cit>  and the correlation between gaping and heat stress.
although gaping has advantageous respiratory effects, it results in water loss caused by both evaporation and incidental discharge of water during valve closure, thus increasing the risk of desiccation  <cit> .
although valve opening suggests possible evaporative cooling, several studies have excluded any link between this behaviour and the body temperature of isolated individuals  <cit>  .
the adaptive significance of gaping behaviour helps explain the vertical intertidal zonation of the western lineage of p. perna and the coexisting invasive mytilus galloprovincialis, a non-gaping mussel species  <cit>  highlighting how even very simple behaviour can influence the outcome of interactions between coexisting species.
our results and previous studies  <cit>  suggest a determinant effect of gaping on the small and large scale spatial distributions of p. perna lineages.
more frequent gaping in the eastern lineage would promote efficient respiration under high, subtropical air temperatures, but more gaping would simultaneously increase the risk of desiccation, restricting the eastern lineage to lower on the shore than western conspecifics .
in addition to differences in gaping rates, the eastern lineage showed much stronger attachment to the substratum.
firm byssal attachment is a prerequisite for survival in the intertidal  but this difference is unlikely to explain large scale spatial segregation of the lineages as the temperate and subtropical regions are not known to differ significantly in terms of wave action.
temperature offers a more likely explanation as it has direct or indirect effects on attachment.
higher temperatures promote byssal thread production while weakening attachment by accelerating thread deterioration and negatively affecting the curing and moulding of byssal threads  <cit> .
additionally, collagenases are common in marine bacteria  <cit>  and their activity is enhanced by temperature  <cit> , potentially contributing to the weakening of the byssus.
thus, the greater attachment strength of mussels from the east coast could be an adaptive response to the thermally more challenging subtropical environment.
CONCLUSIONS
this study shows how distinct genetic lineages, partially segregated along a sharp geographic environmental gradient display adaptive responses to different selective pressures.
p. perna inhabits shores that are experiencing rapid and geographically uneven environmental changes  <cit> .
our results clearly indicate that, rather than the species exhibiting a common overall response, the two lineages are likely to be affected by, and react quite differently to, climate change.
our findings highlight how critical it is to understand intraspecific diversity if we wish to have a grasp of the real magnitude of biodiversity and how within-species diversity will shape ecosystem resilience/resistance in a changing environment.
availability of supporting data
the data set supporting the results of this article is included within the article .
additional files
additional file 1: figure a <dig>  frequencies of null alleles per locus per population.
estimates based on the algorithm presented in brookfield  <cit>  by microchecker.
additional file 2: figure a <dig>  allele frequencies of each location.
each box represents one locus, according to the following order: p <dig>  p <dig>  p <dig>  p <dig>  p <dig>  p <dig>  p <dig>  p <dig>  p <dig>  p <dig>
additional file 3: table a <dig>  observed heterozygosity calculated for each locations and locus.
estimates made in genetix program  <cit> .
additional file 4: figure a <dig>  identifying the best number of k cluster.
plot of mean posterior probability ) values  per cluster  based on  <dig> replicates per k, generated by the structure program  <cit>  and delta k analysis  of lnp  <cit> .
additional file 5: table a <dig>  genotyping table.
scored alleles per location and per locus.
gerardo i zardi and katy r nicastro contributed equally to this work.
competing interests
there are no financial or non-financial competing interests.
authors’ contributions
giz, krn: conceived and designed the experiment.
jc, giz, krn, cdm: performed the experiments.
krn, rc, jc: analyzed the data.
giz, eas, gap: contributed reagents/materials/analysis tools.
giz, krn, cdm, gap, eas, rc: contributed to the writing of the manuscript.
all authors read and approved the final manuscript.
this research was funded by project ptdc/bia-bec/103916/ <dig> from the fundação para a ciência e a tecnologia  and supported by an award from the south africa research chairs initiative  of the department of science and technology .
the funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
we are grateful to r. jacinto for technical support and p. teske for kindly supplying specimens from ponta do ouro.

 the early recognition of functionally significant levels of biodiversity within a species is a central objective of evolutionary biology and biodiversity conservation. this is particularly important at times when biodiversity loss is one of the planet´s major global problems. biodiversity can be investigated, and its importance evaluated, at different levels including species, their genes and functional traits. given the recurrent detection of distinct phylogeographic lineages, it is clearly necessary to consider the adaptive significance of intraspecific differences in order to understand the evolutionary potential of a species. complementary approaches based on neutral genetic markers, commonly employed to estimate recent demographic connectivity, and on methods that can reveal ecologically relevant adaptive traits across heterogeneous habitats are pivotal for the assessment and management of intraspecific diversity. the identification of evolutionarily significant units uniquely based on neutral molecular markers, while overlooking intraspecific patterns of phenotypic divergence, can lead to underrepresentation of locally adapted populations and thus underestimate biodiversity. local environmental conditions resulting in spatially divergent selective regimes can maintain local adaptation across populations of the same species independently of varying and/or contrasting neutral genetic structure. species with wide distributions over variable environments and large effective population sizes often display marked intraspecific biodiversity with distinct genotypes/populations within a species, resulting from adaptation to local environmental conditions but capable of interbreeding with other ecotypes of the same species. marine broadcast spawners with broad geographic ranges extending over abrupt environmental clines are ideal model organisms for studying the interplay between local adaptation and genetic structure due to neutral processes, because divergent selection produced by distinct environments is coupled with a high potential for gene flow among geographically separated subpopulations through larval transport. the southern african coastline covers a wide range of climatic and oceanic conditions and can be divided into three major biogeographic regions: cool-temperate, warm-temperate and sub-tropical. the environmental discontinuities characteristic of these different parts of the coastline are mirrored by a degree of endemism among the different biogeographic provinces and phylogeographic breaks in species that occur in more than one province.figure study area and genetic results. a) map of the study area and sampling sites. black dots: localities either mentioned in the text or used for ecologic data; blue dots: localities used for genetic data abbreviated as in table  continuous, black line: distribution of the western lineage; dashed, black line: distribution of the eastern lineage. the three biogeographic regions of southern africa are delimited by grey dotted lines. b) results of the bayesian approach using the discriminant analysis of principal components scatter plot of individuals with a priori-defined locations, to investigate genetic structure. labels were placed at the centre of dispersion of each group, which are delineated by inertia ellipses. c) dapc a posteriori genetic cluster assignment for k =  the size of each square is proportional to the percentage of individuals assigned to each cluster at each location. d) stacked bar graph of assignment probabilities per individual within p. perna locations across the southern tip of south africa. assignment proportions for each individual were obtained from structure, based on genotypes from seven nuclear microsatellite loci, for k =  the brown mussel perna perna is a habitat-forming species distributed widely around the world in warm waters. in southern africa, it dominates intertidal shores in the sub-tropical and warm-temperate bioregions. it is absent from there to central namibia due to the cold, upwelled waters of the benguela system. it then extends northwards along the west coast of africa to the southern iberian peninsula and into the mediterranean sea as far as the gulf of tunis. analyses of mitochondrial dna sequences indicate a sharp phylogeographic break on the southeast coast of south africa forming a western and an eastern lineage. the western lineage includes animals from both the namibian coast and the south coast of south africa. the eastern lineage includes individuals from the southeast and east coasts of south africa. the two lineages overlap in their distributions on the southeast coast over a distance of approximately km. the presence of the two lineages was recently confirmed by nuclear sequence data. most importantly, mitochondrial and nuclear markers clearly show a non-sister relationship for the two south african lineages. the most likely biogeographic scenario explaining the independent origins of the two south african lineages involves an indo-pacific origin for p. perna, dispersal into the mediterranean and atlantic through the tethys seaway, and recent secondary contact after southward expansion of the western and eastern south african lineages. manipulative experiments and the release of oceanographic drifters suggest that phylogeographic patterns of p. perna may be maintained by a combination of local conditions experienced in the different biogeographic regions, and by the isolating effect of regional oceanographic dynamics. here, we examine the hypotheses that contemporary demographic connectivity and functional traits of p. perna populations are influenced by the tropical/subtropical biogeographic transition. first, we adopted a population genetic approach using microsatellite markers, which are unaffected by natural selection, to investigate differentiation and migration patterns of p. perna between populations inhabiting different bioregions. the re-examination of genetic structure in south african p. perna with microsatellites and comparison with previously reported mitochondrial data is critical to assess potential mito-nuclear discordance and contemporary biogeographic conditions. second, we used a common-garden design to compare fundamental behavioural traits known to affect survival along large and small scale environmental gradients. in particular, we assessed shell gaping and attachment strength, key attributes affecting tolerance of the two major stresses in intertidal habitats, desiccation and wave action. a common-garden design is necessary to distinguish between adaptive and plastic responses of population subjected to distinct selective pressures. finally, we compared intertidal elevation of the two lineages by using in situ temperature profiles to compare intertidal limits of populations from the two bioregions. materials and methods ethics statement mussel collection were performed under permit number res2014/ issued by the department of agriculture forestry and fisheries to the department of zoology and entomology at rhodes university. genetic data dna extraction, amplification and genotyping mussels were collected between february and may at six locations in southern africa in the middle portion of the species’ vertical distribution. mantle tissue was dissected from each individual, preserved in 96% ethanol and stored at -20°c. total genomic dna extraction was performed using a standard proteinase k protocol adapted from. molecular genotyping of a total of samples was carried out using a set of microsatellite loci in three multiplex reactions, with subsequent separation of the pcr products using an abi prism 3130xl dna analyser with genescan liz as size standard and following the procedure of coelho et al.. data were scored using peak-scanner v. software.table genetic diversity of each population locations codes coordinates n h e h o f is Â locations are depicted in the map of figure 1a and are ordered from west to east. decimal coordinates; n: sample sizes; he and ho: expected and observed heterozygosity; fis: inbreeding coefficient with significant values after correction for multiple test indicated with *; Â: allelic richness  ± standard deviation. location-based analyses scoring errors, large allele dropout and null alleles were checked employing the program microchecker and frequency of null alleles estimated based on the algorithm presented in brookfield. additionally, tests for hardy–weinberg equilibrium were conducted using genetix software for all locus-location combinations by assessing the inbreeding estimator wright’s fis and estimating significant values after permutations. tests for linkage disequilibrium between all pairs of loci were performed according to the method of black and kraftsur implemented in genetix. for each location, allele frequencies, observed heterozygosity and expected heterozygosity, plus fis were computed with genetix. standardized allelic richness was estimated by resampling times to standardize to the smallest sample size to account for differences between sampled locations, using standarich r package. pairwise location estimates of genetic differentiation plus 95% bootstrapped confidence intervals were estimated using weir and cockerham’s fst estimator and jost’s dst. calculations were executed in the diversity r package. estimating admixture to evaluate the extent of admixture between the two mitochondrial lineages using microsatellite data, we used structure. clusters varied from one to a maximum of seven, corresponding to the number of locations in the study data plus one. twenty replicate runs were performed under the admixture and the correlated allele frequency model. the most probable value of k was inferred based on the ad hoc criteria l proposed by pritchard et al. and delta k by evanno et al.. similar replicate runs were grouped based on the symmetric similarity coefficient of > using the fullseach algorithm in clumpp and visualized with bar_plotter.rb implemented in ruby. to complement the results of structure, we used discriminant analysis of principal components, which can be performed when uncovering population admixture in the absence of assumptions. this method first transforms the data using principal components analysis, which ensures that the variables are not correlated and that the number of variables is smaller than the number of individuals. we ran two dapc analyses in adegenet r package, following jonker et al.: a location assignment and a genetic cluster assignment. we assigned each individual a priori to its location of origin and obtained for each one the probability of assignment to the correct sampling site. this procedure maximizes the among-location variation and minimizes the within-location variation. using the find.clusters function, we ran successive k-means clustering of the individuals from k =  to k =  and identified the best supported number of clusters through comparison with the bayesian information criterion for the different values of k. we determined the number of clusters and assigned each individual to a genetic cluster without providing any a priori population assignment. in that way, the grouping factor is the genetic cluster and not the sampling location, offering an unbiased interpretation of population structure. to avoid unstable assignments of individuals to clusters, we retained pcs, but used all discriminant functions, in a preliminary dapc run. the results were then reiterated by the optim.a.score function with simulations to determine the optimal number of pcs, and a final dapc was subsequently carried out with the optimal number of pcs. estimating the directionality of gene flow to compare different biogeographic hypotheses on migration rates of p. perna, we used the coalescent-based program migrate-n migrate version.. this approach assumes the wright-fisher model, where locations have a constant effective size through time, the rate of mutation is constant, and locations exchange migrants with constant rates per generation, but those rates can vary among locations. we conducted the analyses on a random subsample of individuals, with the dataset structured into three groups according to the genetic clusters recovered: west, centre and east. we tested four variations of the three-group migration model: all directional migration, strict west to east migration, strict east to west migration and from the centre to the margins. testing the directionality of gene flow is justified because the dominant ocean current between the ocean basins, the agulhas current, runs westerly from the indian towards the atlantic ocean and is thought to play a role in limiting marine dispersal in the opposite direction though coastal counter-currents also occur. models and consider indistinct migration among groups and migration directionality towards the transition area between the two clusters respectively. initial values were calculated using fst. mutation rates were set to be constant among loci. the migrate-n run parameters were calibrated on the full model for convergence of the markov chain monte carlo sampling method. the prior distributions were uniform for mutation-scaled population size parameters theta, that are four times the product of the effective population size and the mutation rate, and mutation-scaled migration rates m, that is, migration rate scaled by the mutation rate, over the range of θ = 0– and m = 0– these settings resulted in converged posterior distributions with a clear maximum for each estimate. the three replicate bayesian runs consisted of one long chain with a total of million states and genealogy changes after discarding the first genealogies as “burn-in” for each locus. for all the analyses we used an adaptive heating scheme with simultaneous chains using different acceptance ratios. overall loci information was combined into a single estimate by bézier approximation of the thermodynamic scores as described by beerli and palczewski and we averaged the bézier score over three different runs and used this as input to evaluate multiple models using bayes factors. ecological data attachment strength and gaping behaviour all mussels used for the assessment of behavioural traits were collected from locations known to support pure mtdna lineages. specimens were collected in february. common garden experiments involve the acclimation and comparison of individuals from distinct geographical locations under identical environmental conditions. such experimental set-ups are widely used to disentangle the effects of genetic and environmental conditions on observed phenotype differences. mussels were divided in three separate tanks, where they were placed at least 10cm apart so that they maintained a solitary position and attachment strengths were not influenced by nearby conspecifics. after three days, the attachment strength of each mussel was measured by drilling a hole into one of the valves and using a spring balance to record the force required to dislodge the mussel perpendicularly to the substratum. mussels were exposed to air for 3h in a controlled environment chamber at 35°c and 60-80% humidity. for each one hour interval, the number of valve movements were noted by visual observation, without recording the width of the gape. intertidal distributional limits temperature dataloggers were used to relate the effective shore height of p. perna at two open coast sites on the east coast and two within the western range of the species. at each site, two transects were established, approximately m apart and running perpendicular to the shoreline. along each transect, two logger were deployed at the upper and two at the lower limit of p. perna distribution. loggers were programmed to record data at min intervals between february 1st to 14th of and march to of, covering the two-week tidal cycle in both regions. temperature logger profiles clearly reveal when loggers are first inundated by the returning tide. this was used to calculate the effective shore height using the online forecasts service of shom database. because tide cycles are not in phase across large geographical scales, durban and mossel bay tidal elevations were used for the east and south coast respectively. where both loggers were recovered at a transect, the values were averaged and used to estimate populations’ intertidal limits at each transect. data analyses each trial was analysed separately. attachment strength data were analysed using a 2-way anova with lineage and block as fixed and random factors respectively. gaping behaviour data were analysed using a 2-way anova with lineage and time as fixed factors. for each intertidal limit, data were analysed under a nested design with range as a fixed factor, site nested within range and transect nested within range and site. heteroscedasticity was tested using levene’s test and post-hoc separation of significant means by snk tests. analyses were performed in spss.  genetic data null alleles at high frequencies were indicated by microchecker for three loci for several localities. for these three loci, departures from hwe were consistently reported across all localities. to avoid overestimation of fis and fst estimators by null alleles, these three loci were removed from all subsequent analyses, including descriptive analyses reported in table  particularly high null allele frequencies have often been reported in several invertebrate taxa, including corals and bivalves, due to high levels of dna sequence variation. specifically, nucleotide variations in microsatellite flanking regions can cause a high frequency of null-alleles, even within populations used for characterization of microsatellites. all microsatellite loci were highly polymorphic among the six populations. the individuals contained alleles in these loci, with an average number of alleles per locus of ranging from to. allelic richness and gene diversity he varied between - and - respectively and the highest allelic richness was recorded at du. mean heterozygote deficiency was detected at all locations with the exception of po. genetic differentiation among populations was detected as indicated by pairwise fst and dst values. while dst pairwise differentiation was lower within than among two groups, fst values did not clearly differentiate pa as belonging to either group.table genetic differentiation between pairs of populations location hm pl pa pe du po codes correspond to locations in figure 1a and are ordered from west to east. pairwise genetic differentiation. 95% confidence intervals are shown in brackets and confidence intervals not enclosing are represented in bold. discriminant analysis of principal components, run with geographic location defined a priori, clearly differentiated the ‘western’ and ‘eastern’ groups on the scatter plot x-axis, confirming the presence of the genetic structure indicated by the fst analysis. individuals from pa contained an admixture of alleles from both groups and indistinctly clustered with each of them. when the a posteriori genetic cluster assignment was run for k = , most individuals of the ‘western’ group were assigned to the same cluster. the individuals of the ‘eastern’ group were also mostly assigned to one cluster, and in a higher rate than the populations of the first cluster. this is particularly true for pe, where all individuals were assigned to cluster individuals from pa were assigned to both clusters with a slightly higher number of individuals assigned to cluster this pattern of subdivisions was largely consistent with the results obtained with structure, which revealed the most significant increase of Δk at two clusters,, thereafter Δk remained unchanged. the l plot rejected the option of k =  we concluded that k =  is the most likely number of genetic clusters. while at most locations individuals displayed an admixture higher than for either group or cluster, pa showed a split of and for the western and eastern groups respectively. when comparing the four biogeographic hypotheses with migrate-n, model iii, strictly east to west migration, was the only model supported.table bayes factors model comparison of migration models for p. perna models model parameters bézier dbézier probability set 1 set 2 set 3 model parameters code as follows: the first set of three signs indicate estimation of theta for west, migration from centre to west, migration from east to west; the second set indicate migration from west to centre, estimation of theta for centre, migration from east to centre; the third set indicates migration from west to east, migration from centre to east and estimation of theta for east. an asterisk indicates that that particular parameter was estimated by the model, and a indicates that no migration was allowed. ecological data in each trial, lineage was a highly significant factor influencing attachment strength, with the eastern lineage having a higher attachment strength than the western.figure ecological results. a) gaping behaviour of the eastern and western lineages, each trial is plotted separately, b) strength of attachment, each trial plotted separately, c) high and low intertidal limits of p. perna western and eastern lineages. dots depict mild outliers and stars are extreme outliers. in each trail, individuals from the eastern lineage always gaped significantly more than those from the western lineage. post-hoc snk tests indicated that both lineages gaped significantly more at hour one than at hours two or three during the first trial, however, during the second trail, there was a time × lineage interaction because at hour two eastern lineage individuals gaped more that at hour three. p. perna in the western range reached a significantly higher intertidal elevation than on the east coast. however, mussels in the western and eastern ranges displayed the same lower vertical limit. discussion the microsatellite data revealed strong genetic division between the two geographically defined groups of populations, supporting previous findings based on mitochondrial and nuclear markers. populations from cluster one inhabit the warm-temperate, south-western coast of south africa, while cluster two comprises populations distributed along the sub-tropical eastern shores. most importantly, genetic admixture was only moderate, asymmetric and not restricted to the previously described km contact zone, which lies mostly in the warm-temperate transition zone ). the admixture detected with microsatellite markers is probably the result of secondary contact between ancient lineages; phylogeographic patterns based on mtdna and ndna sequence data indicate an indo-pacific origin of the species and a dispersal into the mediterranean and atlantic through the tethys seaway. after the closing of the seaway, the two non-sister lineages diverged and came into secondary contact after their respective eastern, or western, southward expansion. however, genetic differentiation between lineages can be primarily attributed to the magnitude of the differences in allele frequencies rather than to unique alleles, as would be expected in ancient isolated lineages. the detected gene exchange was primarily unidirectional, from east to west. additionally, the sampled location at the geographic border of the two lineages shows the strongest gene intermediacy. gene exchange and the maintenance of this pattern could be attributed to several non-mutually exclusive effects. eastward gene flow could be impeded by dispersal barriers, or by selection causing outbreeding depression if lineages are locally adapted. alternatively, poor mating success between lineages could also play a role, either due to divergent reproductive ecology, or to some level of genetic incompatibility. barriers to dispersal mussel larval dispersal is largely influenced by local hydrodynamics. the main oceanographic influence on the east and south coasts of south africa is the agulhas current. this is a major ocean current that flows towards the southeast, bringing warm, oligotrophic water from the mozambique channel. the inshore thermal front of this current varies geographically and can bring its warm waters very close inshore. from near pa southwards the current is deflected away from the coast. the trajectories of nearshore drogues released off the south coast revealed no mixing of southern waters into the eastern water mass, while drifters released off the east coast moved southward, caught up in the agulhas current. although there are wind-driven surface currents that could allow limited eastward larval dispersal, water exchange is predominantly unidirectional, promoting westerly larval dispersal of the eastern lineage and explaining the asymmetrical gene flow described here. previous mtdna data on several species with planktonic larvae support our results obtained with microsatellite data, by showing asymmetric east-to-west migration along the east coast. interestingly, one of our most eastern locations displayed a fair amount of genetic admixture. du is the largest and busiest shipping terminal on the african continent, and globally one of the seven key epicentres for inter-regional exchange of species. shipping activities are known vectors of larval dispersal. the fact that du is the only location on the east coast with these characteristics suggests that the introduction of non-indigenous genes from the western range detected only here is a consequence of human-mediated larval transport. departures from random mating out of the six populations sampled, five displayed significant heterozygote deficiency. null alleles are unlikely to be the cause of high fis values in this case, since microsatellite loci with a high frequency of null alleles were excluded from the analyses, while the remaining loci compiled in the final dataset displayed similar homozygote excess simultaneously across all loci. inbreeding is also unlikely in a broadcast spawner with long a larval pelagic phase as recruitment is expected to produce spatially random kinship. one possible explanation for heterozygote deficiency is a wahlund effect due to the presence of genetically distinct subpopulations in each population sampled. in this context, it is interesting that pa, which is at the contact zone between the two lineages, displayed the lowest mean observed heterozygosity, although this was not true across all loci. such effects can be expected when pooling discrete subpopulations with different allele frequencies that do not interbreed as a single randomly mating unit. in contrast, non-significant levels of heterozygote deficiency were detected in po, which is the farthest population from the contact zone and located on the east side of the coast. reproductive isolation of the lineages, leading to heterozygote deficiency, could result from differences in the timing and frequency of spawning events. it is not known when each lineage spawns nor what triggers the event, but lack of synchrony in spawning by co-existing mussel species has been reported, and analogous intraspecific differences may occur between lineages of p. perna. selection behaviour can radically moderate an organism’s experience of the environment and so dictate physiological reactions, and consequently, tolerance limits to the stress imposed by those environmental conditions. we show that the two p. perna mtdna lineages exhibit clear differences in behavioural traits, indicating an intraspecific behavioural divergence that potentially influences vertical zonation and habitat segregation over large and small spatial scales. by acclimating and comparing individuals from different locations in a common laboratory environment we could differentiate between adaptive and plastic responses to distinct environmental conditions. the maintenance of behavioural phenotypic divergence in experimental populations indicated that differences are due to underlying genetic divergence. mussel gaping, the intermittent opening and closure of the valves, allows the maintenance of aerobic respiration by sustaining an oxygen gradient across the mantle wall and the gills during emersion. respiration rates of organisms usually increase with temperature and several studies have shown that oxygen consumption by mussels increases exponentially with temperature. in p. perna, gaping significantly increased with temperature, supporting the important role of this behaviour in aerobic respiration and the correlation between gaping and heat stress. although gaping has advantageous respiratory effects, it results in water loss caused by both evaporation and incidental discharge of water during valve closure, thus increasing the risk of desiccation. although valve opening suggests possible evaporative cooling, several studies have excluded any link between this behaviour and the body temperature of isolated individuals. the adaptive significance of gaping behaviour helps explain the vertical intertidal zonation of the western lineage of p. perna and the coexisting invasive mytilus galloprovincialis, a non-gaping mussel species highlighting how even very simple behaviour can influence the outcome of interactions between coexisting species. our results and previous studies suggest a determinant effect of gaping on the small and large scale spatial distributions of p. perna lineages. more frequent gaping in the eastern lineage would promote efficient respiration under high, subtropical air temperatures, but more gaping would simultaneously increase the risk of desiccation, restricting the eastern lineage to lower on the shore than western conspecifics. in addition to differences in gaping rates, the eastern lineage showed much stronger attachment to the substratum. firm byssal attachment is a prerequisite for survival in the intertidal but this difference is unlikely to explain large scale spatial segregation of the lineages as the temperate and subtropical regions are not known to differ significantly in terms of wave action. temperature offers a more likely explanation as it has direct or indirect effects on attachment. higher temperatures promote byssal thread production while weakening attachment by accelerating thread deterioration and negatively affecting the curing and moulding of byssal threads. additionally, collagenases are common in marine bacteria and their activity is enhanced by temperature, potentially contributing to the weakening of the byssus. thus, the greater attachment strength of mussels from the east coast could be an adaptive response to the thermally more challenging subtropical environment.  this study shows how distinct genetic lineages, partially segregated along a sharp geographic environmental gradient display adaptive responses to different selective pressures. p. perna inhabits shores that are experiencing rapid and geographically uneven environmental changes. our results clearly indicate that, rather than the species exhibiting a common overall response, the two lineages are likely to be affected by, and react quite differently to, climate change. our findings highlight how critical it is to understand intraspecific diversity if we wish to have a grasp of the real magnitude of biodiversity and how within-species diversity will shape ecosystem resilience/resistance in a changing environment. availability of supporting data the data set supporting the results of this article is included within the article. additional files additional file 1: figure a frequencies of null alleles per locus per population. estimates based on the algorithm presented in brookfield by microchecker. additional file 2: figure a allele frequencies of each location. each box represents one locus, according to the following order: p p p p p p p p p p additional file 3: table a observed heterozygosity calculated for each locations and locus. estimates made in genetix program. additional file 4: figure a identifying the best number of k cluster. plot of mean posterior probability ) values per cluster based on replicates per k, generated by the structure program and delta k analysis of lnp. additional file 5: table a genotyping table. scored alleles per location and per locus. gerardo i zardi and katy r nicastro contributed equally to this work. competing interests there are no financial or non-financial competing interests. authors’ contributions giz, krn: conceived and designed the experiment. jc, giz, krn, cdm: performed the experiments. krn, rc, jc: analyzed the data. giz, eas, gap: contributed reagents/materials/analysis tools. giz, krn, cdm, gap, eas, rc: contributed to the writing of the manuscript. all authors read and approved the final manuscript. this research was funded by project ptdc/bia-bec/103916/ from the fundação para a ciência e a tecnologia and supported by an award from the south africa research chairs initiative of the department of science and technology. the funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. we are grateful to r. jacinto for technical support and p. teske for kindly supplying specimens from ponta do ouro.
gene flow between regions was predominantly east-to-west and was only moderate, with higher heterozygote deficiency where the two lineages co-occur. we investigated genetic differentiation, contemporary patterns of demographic connectivity and intraspecific variation of adaptive behavioural traits in two lineages of an intertidal mussel across a tropical/subtropical biogeographic transition. microsatellite analyses revealed clear genetic differentiation between western and eastern populations, confirming divergence previously detected with mitochondrial and nuclear markers. common garden experiments showed that gaping and attachment to the substratum differed significantly between the two lineages. this can be explained by differential selection and/or oceanographic dynamics acting as a barrier to larval dispersal.
 intraspecific variability is seen as a central component of biodiversity. we investigated genetic differentiation, contemporary patterns of demographic connectivity and intraspecific variation of adaptive behavioural traits in two lineages of an intertidal mussel across a tropical/subtropical biogeographic transition.  microsatellite analyses revealed clear genetic differentiation between western and eastern populations, confirming divergence previously detected with mitochondrial and nuclear markers. gene flow between regions was predominantly east-to-west and was only moderate, with higher heterozygote deficiency where the two lineages co-occur. this can be explained by differential selection and/or oceanographic dynamics acting as a barrier to larval dispersal. common garden experiments showed that gaping and attachment to the substratum differed significantly between the two lineages. western individuals gaped more and attached less strongly to the substratum than eastern ones.  these behavioural differences are consistent with the geographic and intertidal distributions of each lineage along sharp environmental clines, indicating their strong adaptive significance. we highlight the functional role of diversity below the species level in evolutionary trends and the need to understand this when predicting biodiversity responses to environmental change. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords genetic diversitymarine connectivityperna pernaissue-copyright-statement© the author 2015 