# Spatial genetic structure in Beta vulgaris subsp. maritima and Beta macrocarpa reveals the effect of contrasting mating system, influence of marine currents, and footprints of postglacial recolonization routes.

Marie Leys, Eric Petit, Yasmina El-Bahloul, Camille Liso, Sylvain Fournet, Jean-François Arnaud

► To cite this version:

Marie Leys, Eric Petit, Yasmina El-Bahloul, Camille Liso, Sylvain Fournet, et al. Spatial genetic structure in Beta vulgaris subsp. maritima and Beta macrocarpa reveals the effect of contrasting mating system, influence of marine currents, and footprints of postglacial recolonization routes. Ecology and Society, 2014, 4 (10), pp.1828-1852. 10.1002/ece3.1061. hal-01060474

HAL Id: hal-01060474

https://univ-rennes.hal.science/hal-01060474

Submitted on 3 Sep 2014

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Spatial genetic structure in *Beta vulgaris* subsp. *maritima* and *Beta macrocarpa* reveals the effect of contrasting mating system, influence of marine currents, and footprints of postglacial recolonization routes

Marie Leys<fnref n="1" />, Eric J. Petit<fnref n="2" />, Yasmina El-Bahloul<fnref n="3" />, Camille Liso<fnref n="1" />, Sylvain Fournet<fnref n="4" /> & Jean-François Arnaud<fnref n="1" />

**Keywords**

Cytonuclear genetic diversity, glacial refugia, marine currents, postglacial expansion, spatial genetics.

**Correspondence**

Jean-François Arnaud, Laboratoire de Géné-
tique et Évolution des Populations Végétales,
UMR CNRS 8198, Bâtiment SN2, Université des
Sciences et Technologies de Lille - Lille 1,
Villeneuve d'Ascq Cedex, F-59655, France.

Tel: +33 (0)3 20 33 63 97;

Fax: +33 (0) 3 20 43 69 79;

E-mail: jean-francois.arnaud@univ-lille1.fr

**Present address**

Marie Leys, Swiss Federal Institute of Aquatic
Science and Technology, Department Aquatic
Ecology, Ueberlandstrasse 133, Duebendorf
8600, Switzerland

**Funding information**

This work was funded by grants from the
"PRAD-PHC Maroc (09-01)," the "FRB/
Région Nord-Pas de Calais" (GENEFRAG
project), and the "INRA SPE."

Received: 3 February 2014; Revised: 17
March 2014; Accepted: 18 March 2014

Ecology and Evolution 2014; 4(10): 1828–1852

doi: 10.1002/ece3.1061

**Abstract**

Understanding the factors that contribute to population genetic divergence across a species' range is a long-standing goal in evolutionary biology and ecological genetics. We examined the relative importance of historical and ecological features in shaping the present-day spatial patterns of genetic structure in two related plant species, *Beta vulgaris* subsp. *maritima* and *Beta macrocarpa*. Using nuclear and mitochondrial markers, we surveyed 93 populations from Brittany (France) to Morocco – the southern limit of their species' range distribution. Whereas *B. macrocarpa* showed a genotypic structure and a high level of genetic differentiation indicative of selfing, the population genetic structure of *B. vulgaris* subsp. *maritima* was consistent with an outcrossing mating system. We further showed (1) a strong geographic clustering in coastal *B. vulgaris* subsp. *maritima* populations that highlighted the influence of marine currents in shaping different lineages and (2) a peculiar genetic structure of inland *B. vulgaris* subsp. *maritima* populations that could indicate the admixture of distinct evolutionary lineages and recent expansions associated with anthropogenic disturbances. Spatial patterns of nuclear diversity and differentiation also supported a stepwise recolonization of Europe from Atlantic-Mediterranean refugia after the last glacial period, with leading-edge expansions. However, cytoplasmic diversity was not impacted by postglacial recolonization: stochastic long-distance seed dispersal mediated by major oceanic currents may mitigate the common patterns of reduced cytoplasmic diversity observed for edge populations. Overall, the patterns we documented here challenge the general view of reduced genetic diversity at the edge of a species' range distribution and provide clues for understanding how life-history and major geographic features interact to shape the distribution of genetic diversity.

**Introduction**

Plant populations are often composed of distinct breeding units, and patterns of genetic variation within and among

populations are determined primarily by mating system and the neutral microevolutionary process of gene flow (Hamrick and Godt 1996; Duminil et al. 2009). Historical events such as Pleistocene climatic changes have also

<footnote n="1">¹Laboratoire de Génétique et Évolution des Populations Végétales, UMR CNRS 8198, Bâtiment SN2, Université des Sciences et Technologies de Lille - Lille 1, Villeneuve d'Ascq Cedex F-59655, France</footnote>

<footnote n="2">²UMR CNRS 6553 ECOBIO, Station biologique, Université de Rennes 1, Paimpont F-35380, France</footnote>

<footnote n="3">³Unité d'Amélioration des Plantes Conservation et Valorisation des Ressources Phytogénétiques, Centre Régional de la Recherche Agronomique de Rabat, INRA-Maroc, Rabat-Instituts 10101, Morocco</footnote>

<footnote n="4">⁴UMR 1349 IGEPP, INRA - Agrocampus Ouest-Université de Rennes 1, Bât 320, BP35327, Le Rheu Cedex 35653, France</footnote>
played a fundamental role in generating phylogeographic structure within species' ranges (e.g., Hewitt 2004). In the northern hemisphere, the alternation of Quaternary glacial and interglacial episodes resulted in substantial changes in species' distribution. These events have left footprints in their genetic structure and may explain the common poleward decrease in genetic diversity exhibited by terrestrial taxa (Hewitt 1999).

From the extensive body of phylogeographic work on European temperate species, it is acknowledged that the Mediterranean Basin acted as a major refugium during these glacial cycles. In particular, the Iberian, Italian, and Balkan Peninsulas have been considered as differentiation centers and origins of postglacial northward expansion (Hewitt 2000). However, the literature remains largely biased with respect to geographic coverage and the type of organism studied. In contrast to Western Europe, North Africa is largely underrepresented. This imbalance is unfortunate, because North Africa harbors hot spots of plant biodiversity and endemism (Médail and Quézel 1997). Furthermore, most of the studies conducted in North Africa have focused on inland woody species (e.g., Lumaret et al. 2002; Migliore et al. 2012) rather than on coastal plants (Ortiz et al. 2007). Coastal plant species offer, however, ideal case studies for interpreting phylogeographic patterns. They have a very wide latitudinal and longitudinal range, but are restricted to a narrow belt along the shoreline, suggesting limited migration routes and/or postglacial recolonization pathways. In addition to climatic oscillations, the regional distribution of coastal plant species may also be influenced by marine currents *via* seed drift (dispersal), leading to sharp genetic discontinuities or nonintuitive genetic affinities among populations (Kadereit et al. 2005).

Beyond the influence of historical events, patterns of genetic structure across a species range are also informative about contemporary processes. Current spatiotemporal dynamics of plant species' distribution may indeed be influenced by anthropogenic disturbance, acting as dispersal agents through exchanges or interfering with natural recolonization by modifying the landscape (D'Hondt et al. 2012). Human impact on the environment is a long-standing feature of the Mediterranean Basin (Blondel 2006). For example, the long history of agriculture and other ancient and recent human activities have resulted in an evolving landscape mosaic, with land uses changing in space and time. Therefore, population-level processes are highly complex, and spatial genetic structure is thus tightly coupled to ongoing processes that superimpose new signatures on existing natural phylogeographic patterns. The interplay between mating system, geographic distribution, oceanographic features/sea circulation, past climatic changes, and anthropogenic pressures is likely

responsible for the current genetic diversity and species distribution of coastal plants along the southern Atlantic coast and the Mediterranean Basin (Rodríguez-Sánchez et al. 2008).

The present study set out to depict large-scale patterns of cytoplasmic and nuclear genetic structure in two related species thought to have contrasting mating systems: the wild sea beet (*Beta vulgaris* subsp. *maritima*) and its close relative, *Beta macrocarpa*. Sea beet is a typical coastal taxon, mainly outcrossing, with a wide distribution from the Canary Islands in the west, northward along Europe's Atlantic coast and Baltic Seas (Biancardi et al. 2012). Sea beet also extends eastward through the Mediterranean Basin along the coastline but also in inland man-made habitats (Desplanque et al. 1999; Arnaud et al. 2009). *B. macrocarpa*, a self-compatible species, has a more restricted geographic distribution located in inland or coastal habitats in western and eastern Mediterranean areas (Buttler 1977; Letschert 1993; Castro et al. 2013). Elucidating the underlying mechanisms of range shifts and the determinants of population establishment are fundamental questions, especially in the context of climate change. Based on a comprehensive sampling in Western Europe and in North Africa and using nuclear and cytoplasmic data, we thus applied classical *F*-statistics, Bayesian clustering, and spatial multivariate methods to depict the phylogeographic and genetic structure of both species. Specifically, we asked the following questions:

1. Do *B. vulgaris* subsp. *maritima* and *B. macrocarpa* show significant genetic differentiation and different genotypic structure indicative of contrasting mating system, as suggested by previous studies based on greenhouse experiments and allozyme survey?

2. What are the roles of alternative evolutionary processes in shaping different lineages along the southern European/northern African distribution of these two species? We hypothesized that range shift driven by last Quaternary glaciation, ocean circulation patterns, and anthropogenic disturbances could be fundamental processes explaining large-scale patterns of genetic diversity and divergence we observed.

## Materials and Methods

### The species

The genus *Beta* is divided into two sections, that is, the *Beta* section and the *Corollinae* section (Biancardi et al. 2012). The *Beta* section includes five taxa: (1) the sea beet (*Beta vulgaris* subsp. *maritima*, see Fig. 1), which is considered as the wild ancestor of all cultivated beets; (2) the different forms of cultivated beets (*Beta vulgaris* subsp.
*vulgaris*), including fodder, leaf, or sugar beets; (3) *Beta macrocarpa*, which is an annual self-compatible plant thought to reproduce predominantly by autogamy; (4) *Beta patula*, which is endemic to an islet of the Madeira Archipelago; and (5) *Beta vulgaris* subsp. *adanensis*, a self-compatible species that is only found in Eastern Mediterranean areas, Greece and Turkey (Frese et al. 1990; Kade-reit et al. 2006; Biancardi et al. 2012). Therefore, in our extensive sampling area in the western Mediterranean, only two species of the *Beta* section can be found in coastal and inland ruderal habitats: *Beta vulgaris* subsp. *maritima* and *Beta macrocarpa*.

Sea beet is a short-lived perennial species growing in the supralittoral fringe (Fig. 1), typically on beaches of sand or pebbles, or associated with rocky substrates such as cliffs overhanging the sea (De Cauwer et al. 2010). These ecosystems are characterized by linear physical boundaries, increased salinity, periodic strong winds, and soil instability due to water movements during the highest tides. In addition, a particular wild form of *B. vulgaris* subsp. *maritima* has colonized human-disturbed inland habitats in the western part of the Mediterranean area (Desplanque et al. 1999). These inland ruderal beets have founded populations associated with past and present road works (roadsides, car parks), wastelands, rubble deposits, or gardens (Arnaud et al. 2009). Sea beet is a diploid (2n = 18) gynodioecious species, has an outcrossing mating system – enforced by a gametophytic self-incompatibility system with up to four gametophytic S loci – and relies on wind pollination (Owen 1942). No vegetative reproduction occurs, and dispersal is only mediated by seed and/or pollen movements. Seeds are clustered in an irregular dry fruit, forming a multigerm

seedball composed of one to eight seeds (De Cauwer et al. 2011). Although this seedball is mainly dispersed by gravity with a large proportion of seeds falling in the vicinity of the mother plant (Arnaud et al. 2011; De Cauwer et al. 2012), it has a corklike layer around the seeds, which may confer buoyancy and make it amenable to dispersal in seawater (see Fievet et al. 2007).

Only flowering characteristics such as tepal morphology allow distinguishing between *B. v.* subsp. *maritima* and *B. macrocarpa*. *B. macrocarpa* is a self-compatible species thought to be predominantly self-fertilizing (Jung et al. 1993; Bruun et al. 1995). Nonmorphologically distinguishable diploid (2n = 18) and tetraploid (2n = 36) cytotypes occur in *B. macrocarpa* (Castro et al. 2013). Using diverse accessions of western Mediterranean areas, both cytotypes were found to be annual and fully self-compatible (Buttler 1977; Lange and Debock 1989). Letschert (1993), Bartsch and Ellstrand (1999), and Villain (2007) provided some evidence suggesting nearly complete autogamy in both diploid and tetraploid accessions of *B. macrocarpa*.

## Study area and sampling

We sampled 93 localities distributed along the eastern Atlantic coast, from Brittany, France, to southern Morocco, but also in the Mediterranean region on both sides of the Strait of Gibraltar (Fig. 2A). In Western Europe, samples were collected at 36 sites and labelled from 1 to 36. In Morocco, 52 coastal and inland localities were sampled, among which 32 sites were representative of typical coastal populations collected along the shoreline (labelled from 37 to 68, Fig. 2A) and 20 were from inland ruderal sites (labelled from *a* to *t*, Fig. 2A). Five sites in the Madeira archipelago were also sampled and numbered from 69 to 73. In total, 1816 individuals were collected with an average sample size per site of 19.53 (SE ± 6.07). The population labelled 67 was made up of only four individuals and was not included in analyses dealing with mean levels of genetic diversity, such as the estimation of allelic richness. Sampling was carried out during five successive surveys between 2007 and 2011. Fresh leaves from the collected individuals were dried and stored in silica gel at room temperature until DNA extraction. Appendix 1 provides full details on the sampling localities.

## Genetic data collection

Extraction and purification of total DNA were performed using a DNeasy 96 Plant Kit following the standard protocol for isolation of DNA from plant leaf tissue outlined in the DNeasy 96 Plant protocol handbook (QIAGEN Inc., Courtaboeuf, France).

<figure><img src="image_2.png" /><figcaption>**Figure 1.** Wild sea beet (*Beta vulgaris* subsp. *maritima*), the wild ancestor of all cultivated beets. Sea beet is a typical coastal taxon growing in the supralittoral fringe, typically on beaches of sand or pebbles or associated with rocky substrates.</figcaption></figure>
<figure><img src="image_3.png" /><figcaption>**Figure 2.** (A) Map of the geographic location of the 93 sampled stations for *Beta vulgaris* subsp. *maritima* and *Beta macrocarpa*. Each *B. vulgaris* subsp. *maritima* station is represented by a black diamond. *B. macrocarpa* stations are symbolized by blue diamonds and stations composed of a mixture of both species by red diamonds. (B) Bayesian clustering following Pritchard et al. (2000). Assignment probabilities of individual memberships of *B. vulgaris* subsp. *maritima* and *B. macrocarpa* into three inferred clusters for the first modal value (K = 3). Each individual is represented by a thin horizontal line (y axis) partitioned into three colored segments that represent the individual's estimated membership fractions to the K = 3 clusters (x-axis). (C) Map of mean population membership probabilities for the three clusters.</figcaption></figure>
## Nuclear microsatellites

All individuals were examined for nuclear genetic variation using eight microsatellite loci, named *Bvm3* (Mörchen et al. 1996), *Gtt1*, *Caa1* (Viard et al. 2002), *Sb04*, *Sb06*, *Sb07*, *Sb15* (Richards et al. 2004), and *FDSB1027* (McGrath et al. 2007). Polymerase chain reaction (PCR) conditions and allele-sizing procedures can be found in Fénart et al. (2007). Electrophoresis and genotyping were performed on an ABI automated sequencer model 3110 (Applied Biosystems Inc.) as described in De Cauwer et al. (2010). To minimize the rate of genotyping errors, a second round of PCR and electrophoresis were performed for individuals with dubious multilocus genotypes, that is, with missing data or displaying rare alleles.

## Mitochondrial minisatellites

To gain information on cytoplasmic genetic diversity, mitochondrial polymorphism was characterized by genotyping all individuals at four mitochondrial minisatellite loci: *TR1*, *TR2*, *TR3*, and *TR4* (Nishizawa et al. 2000). PCR amplifications, cycling conditions, and genotyping were performed according to protocols previously described in Fénart et al. (2008). As the mitochondrial genome is maternally inherited as a single linkage unit, each genotype combination for these four minisatellite loci was analyzed as a single haplotype for the levels of cytoplasmic diversity and population differentiation, as in De Cauwer et al. (2012).

# Statistical analyses of cytoplasmic and nuclear genetic variation

## Distinguishing between *B. vulgaris* subsp. *maritima* and *B. macrocarpa*

As pointed out by Buttler (1977), only tepal characters are relevant criteria to distinguish between *B. v.* subsp. *maritima* and *B. macrocarpa*. However, sampling was carried out at the beginning of the flowering season, when *B. vulgaris* subsp. *maritima* shows large phenotypic plasticity and is morphologically similar to *B. macrocarpa*. However, as *B. macrocarpa* is a self-fertilizing species, the detection of multilocus fixed-homozygote genotypes can help to determine a posteriori which individuals belong to *B. macrocarpa*. Because (1) the two species may possibly hybridize (Kishima et al. 1987; Lange and Debock 1989; Bartsch and Ellstrand 1999) and (2) *B. macrocarpa* tetraploid forms with fixed heterozygosity can occur (Letschert 1993), we also discriminated individuals on the basis of their probability of belonging to one or the other species due to allele frequency differences in multilocus genotypes. To do so, we used the model-based Bayesian algorithm

implemented in STRUCTURE and described in Pritchard et al. (2000). Bayesian clustering assumes that each individual has admixed ancestral origins in gene pools that have correlated allele frequencies because of migration or shared ancestry (Falush et al. 2003). Using the whole data set, we assessed the number of potential clusters (*K*) from 30 runs along a range of *K* varying from 1 to 50 with 2.10⁶ Markov chain Monte Carlo replications following a burn-in period of 100,000 iterations without any prior information on the putative population affiliation of individuals. Among the 30 replicates, 15 with the highest likelihood were retained for subsequent analyses. The ad hoc statistic Δ*K* was then calculated to determine the most accurate number of *K* clusters (Evanno et al. 2005).

The Bayesian clustering described in Pritchard et al. (2000) assigns individuals by creating groups within which Hardy-Weinberg (HW) disequilibrium and linkage disequilibrium (LD) are minimized. Nonetheless, inbreeding may induce departures from HW expectations and LD among loci, which can lead to bias or erroneous inferences. To circumvent this, we ran two parallel extensive runs by keeping or removing *B. macrocarpa* individuals to verify the consistency of clustering results.

Additional analyses were conducted using spatially explicit Bayesian clustering, that is, incorporating spatial trends and autocorrelation in the prior distribution on individual admixture coefficients, as described in Chen et al. (2007) and Durand et al. (2009). Using TESS version 2.3 (Grenoble INP -TIMC-IMAG, Faculty of Medicine, La Tronche, France), we ran 1.10⁶ (100,000 burn-in) MCMC iterations from *K* = 1 to *K* = 20 (25 replicates per *K* value) using the BYM admixture model described in Durand et al. (2009) and the geographic distances between individuals. Of the 25 replicates per *K* tested, 20 with the lowest deviance information criterion (DIC) were retained. We set the Dirichlet parameter of the allele frequency model, the trend degree to 1.0, and the admixture and spatial interaction parameters to default values. To reveal which *K* may provide the best fit to the genetic data, average DIC values were plotted against *K*.

Similarity coefficients between runs and the average matrices of individual membership proportions were estimated using CLUMPP version 1.1.2 (Jakobsson and Rosenberg 2007). Clusters were displayed using DISTRUCT version 1.1 (Rosenberg 2004).

## Basic parameters of genetic diversity

Genetic variation was examined per locus and per population with descriptive statistics: the total and mean number of alleles (*A*<sub>t</sub> and *A*<sub>m</sub>, respectively), allelic richness (*A*<sub>r</sub>), and level of gene diversity (*H*<sub>E</sub>), all using FSTAT version 2.9.3.2 (Goudet 1995). Allelic richness was
estimated using the rarefaction approach proposed by El Mousadik and Petit (1996) and based on a minimal population size of eight diploid individuals. For additional information on population genetic uniqueness, estimates of private allelic richness ($A_rP$) were computed following a rarefaction procedure ($n = 8$) using ADZE software (Szpiech et al. 2008).

Genotypic linkage disequilibrium among all locus-pair combinations was assessed prior to other analyses using a Markov chain approximation of the Fisher's exact test, based on the contingency tables for all pairs of loci in each population, as implemented in GENEPOP version 4.0.1.1 (Raymond and Rousset 1995). Departures from HW equilibrium for each microsatellite locus and overall loci were quantified using the intrapopulation fixation index ($F_{IS}$). Statistical significance of $F_{IS}$ values was subsequently assessed using the randomization procedure provided by FSTAT. Sequential Bonferroni corrections for multiple comparisons were applied following Rice (1989).

Levels of genetic variation ($A_r, A_rP$, and $H_E$) were compared between the genetically distinct groups depicted through Bayesian clustering and multivariate analyses described later. Statistical significance of differences was tested using a one-way analysis of variance followed by Tukey's multiple-comparisons test. To assess whether a northward-declining diversity gradient occurred in *B. vulgaris* subsp. *maritima*, linear regressions were applied on multilocus genetic diversity estimates and spatial attributes of populations: latitude and geographic distance measured along the coastline from the most southern sampling site. Linear regressions were performed on separate population data sets to take into account the presence of genetic discontinuities. A potential bias could arise because levels of allelic diversity often varied among the loci. Therefore, we also performed linear mixed models in analyzing genetic diversity as a function of geographic locations, with loci as random intercepts.

## Extent of genetic differentiation

We investigated population structure using classical F-statistics according to Weir and Cockerham (1984) to test for population differentiation across all populations, using FSTAT. Similarly to $F_{IS}$ values, we used a permutation test to determine whether observed values of $F_{ST}$ were significantly different from 0 for each locus and overall loci (10,000 randomizations of multilocus genotypes among populations). Furthermore, we compared the mean pairwise genetic differentiation between genetic clusters following the statistical procedure accounting for nonindependent data described in Coulon et al. (2006). The population grouping corresponded to the genetic clusters depicted through Bayesian clustering and multivariate

analyses described later. The apportionment of genetic variation was also assessed following the hierarchical procedure proposed by Yang (1998). Hierarchical F-statistics were carried out within and among genetically distinct groups depicted by clustering and multivariate analyses described later. Using the R package HIERFSTAT (Goudet 2005), the significance of genetic differentiation among populations was assessed by randomly permuting individuals (10,000 permutations) among populations within groups. Genetic differentiation among genetic groups was tested by carrying out 10,000 permutations of populations among genetic groups.

## Spatial multivariate analysis

We also analyzed our data set using a spatially explicit multivariate method: the spatial principal component analysis (sPCA) described in Jombart et al. (2008). The advantage of sPCA is that it imposes no genetic assumptions on mating system, population structure, or allele frequency models (Jombart et al. 2009). Basically, the method aims to find independent synthetic variables that maximize the product of the spatial autocorrelation calculated on a set of allelic frequencies using Moran's index and the genetic variance among individuals or populations. Spatial information used for the computation of Moran's *I* is stored in a spatial weighting binary matrix determined through a neighborhood graph. The Gabriel graph was selected to perform the sPCA because it showed the best fit with our genetic data (see Table S1). To evaluate the consistency of spatial patterns observed, the significance of global and local spatial structures was assessed using permutation tests as described in Jombart et al. (2008). Main results provide a geographic map of entity scores allowing a visual assessment of spatial genetic structure. To draw a comprehensive synthetic representation, each of the first three principal component scores was simultaneously represented into a channel of colors as in Menozzi et al. (1978). All the analyses were performed using the ADEGENET package implemented in R (Jombart 2008; R Development Core Team 2011).

## Analysis of isolation-by-distance patterns

We studied the relationships between genetic divergence and spatial isolation of populations using unidirectional Mantel correlograms (Oden and Sokal 1986; Arnaud et al. 2009). Genetic divergence was quantitated using the $D_{CE}$ chord distance (Cavalli-Sforza and Edwards 1967). These analyses were used to compare the strength of spatial structuring as depicted by nuclear and cytoplasmic data. Mantel correlograms were designed using the normalized Mantel statistics (rz, Smouse et al. 1986) following the
procedure described in Oden and Sokal (1986). For each
distance class, designed with a minimum of 15 pairwise
comparisons, the significance of *rz* values was tested with
a Mantel test (10,000 permutations). All analyses were
performed using PASSAGE version 1.1.2.3 (Rosenberg
and Anderson 2011).

**Results**

Distinguishing between *B. vulgaris* subsp.
*maritima* and *B. macrocarpa*

The ΔK statistic showed a strong mode at K = 3: Bayes-
ian clustering analysis separated the whole data set into
three genetically distinct groups (Figs. 2B,C, S1). The
depicted clusters corresponded to (1) European sea beet
populations also including some Moroccan populations
located near Gibraltar; (2) the remaining Moroccan sea
beet populations; and (3) all B. macrocarpa populations.
We are very confident that the three depicted genetic
clusters match B. vulgaris subsp. maritima and B. macro-
carpa delimitations. Indeed, along with particular multilo-
cus genotypes found in the third cluster for which we
hypothesize to be exclusively composed of B. macrocarpa
individuals, the hierarchical nuclear genetic differentiation
between the two species based on this clustering pattern
was the strongest and highly significant (F<sub>CT</sub> = 0.288,
P < 0.001). Among the two groups of sea beets, a signifi-
cant hierarchical differentiation supported this clustering:
Genetic differentiation among populations within these
two clusters accounted for most of the genetic variability
(F<sub>PR</sub> = 0.131; P < 0.001), and the differentiation between
the two clusters, although weaker, was highly significant
(F<sub>RT</sub> = 0.039; P < 0.001). The K versus ΔK distribution
was multimodal with another mode at K = 6 suggesting
further hierarchical structure within the sampled area
(Figs. S1, S2). The similarity coefficients across indepen-
dent runs were all larger than 0.998, which suggests the
absence of genuine multimodality across runs. Inbreeding
did not induce any bias because an identical clustering pat-
tern was found when only analyzing B. vulgaris subsp.
maritima populations (Fig. S2). Assisting the Bayesian
clustering using spatial distribution of individuals yielded
identical results but with lower statistical support: five
B. vulgaris subsp. maritima genetic clusters were also
inferred in TESS analyses, but similarity coefficients among
runs were considerably low, indicating genuine multimo-
dality (H' = 0.599, from 0.457 to 0.737; see Fig. S3).

Overall, based on genotypic structure and on the exten-
sive runs of Bayesian clustering analyses, eight Moroccan
sites and two localities on Santo Porto island were
composed of *B. macrocarpa* individuals and of an admix-
ture of *B. macrocarpa* and *B. vulgaris* subsp. *maritima*

individuals (sites labelled 41, 65, 66, b, h, j, 69, and 70, see Fig. 2 and Appendix). Although in sympatry, no admixture indicative of hybridization was observed between the two species.

Genotypic linkage, Hardy-Weinberg
disequilibrium, and levels of nuclear and
cytoplasmic polymorphism

Within each *B. vulgaris* subsp. *maritima* population, exact tests for LD between microsatellite loci showed only 31 significant *P*-values of 2380 comparisons (1.30%) after Bonferroni correction (119 expected from type I error at $\alpha$ = 0.05). LD was only found in 13 of the 85 *B. vulgaris* subsp. *maritima* populations and involved different pairs of loci. This suggests spurious LD due to genetic drift and/or population substructuring effects (De Cauwer et al. 2012).

Using multilocus probability tests, we detected signifi-
cant departures from HW equilibrium in 32 of the 85
populations of *B. vulgaris* subsp. *maritima*. However, only
a few of these populations displayed consistent F<sub>IS</sub> values
across loci, suggesting departures from panmixia (Appen-
dix). Over all *B. vulgaris* subsp. *maritima* populations,
five of eight loci showed significant deviation from HW
genotypic proportions, although these F<sub>IS</sub> estimates were
very low, ranging from 0.007 to 0.087 (Table 1). Com-
bined probability tests over all loci and across all popula-
tions resulted in an overall departure from HW due to
only few significant single-locus departures within popu-
lations or significant multilocus departures in only a few
populations. In contrast, all but two populations of the
third cluster defined as corresponding to *B. macrocarpa*
exhibited very high, significantly positive, single- and
multilocus F<sub>IS</sub> values (Table 1 and Appendix). The only
exceptions were the *B. macrocarpa* populations located on
Santo Porto Island in Madeira archipelago, one of which
(population labelled 70) had a highly negative mean F<sub>IS</sub>
value due to fixed heterozygosity at multilocus genotypes
(Appendix). The other one was composed of individuals
characterized by either fixed homozygosity or fixed het-
erozygosity.

Nuclear microsatellite loci for *B. vulgaris* subsp. *maritima* populations displayed moderate to high numbers of alleles, ranging from 10 (*Gtt1*) to 56 (*Bvm3*), with an average of 32.25 alleles and a mean gene diversity (*H*<sub>E</sub>) of 0.679 (Table 1). Mitochondrial genetic diversity was much lower, with a number of alleles varying from 3 (*TR2*) to 15 (*TR1*). Based on a single haplotype determined by combining the alleles at four loci, we found a total of 45 haplotypes, of which four represented 77.31% of the total number of haplotypes. In sharp contrast, *B. macrocarpa* displayed much lower levels of genetic
<table><caption>**Table 1.** Genetic diversity measures and *F*-statistics following Weir and Cockerham (1984) estimated over all *Beta vulgaris* subsp. *maritima* (*n* = 85) and *B. macrocarpa* (*n* = 8) populations for eight nuclear microsatellite and four cytoplasmic minisatellite loci.</caption><thead><tr><th>Locus</th><th>A<sub>t</sub></th><th>A<sub>r</sub></th><th>H<sub>E</sub></th><th>Allelic size range (bp)</th><th>F<sub>IT</sub></th><th>F<sub>ST</sub></th><th>F<sub>IS</sub></th></tr></thead><tbody><tr><td colspan="8"><strong>B. maritima</strong></td></tr><tr><td colspan="8">Nuclear</td></tr><tr><td>FDBS1027</td><td>32</td><td>6.426</td><td>0.614</td><td>170-228</td><td>0.184*</td><td>0.170*</td><td>0.018</td></tr><tr><td>Bvm3</td><td>56</td><td>12.043</td><td>0.831</td><td>093-157</td><td>0.195*</td><td>0.134*</td><td>0.070*</td></tr><tr><td>Caa1</td><td>52</td><td>10.059</td><td>0.809</td><td>128-230</td><td>0.205*</td><td>0.130*</td><td>0.087*</td></tr><tr><td>Gtt1</td><td>10</td><td>3.650</td><td>0.418</td><td>104-127</td><td>0.200*</td><td>0.177*</td><td>0.028</td></tr><tr><td>Sb04</td><td>22</td><td>5.320</td><td>0.608</td><td>167-205</td><td>0.223*</td><td>0.172*</td><td>0.062*</td></tr><tr><td>Sb06</td><td>14</td><td>6.675</td><td>0.699</td><td>145-183</td><td>0.198*</td><td>0.168*</td><td>0.037*</td></tr><tr><td>Sb07</td><td>44</td><td>11.505</td><td>0.827</td><td>240-286</td><td>0.181*</td><td>0.127*</td><td>0.062*</td></tr><tr><td>Sb15</td><td>28</td><td>6.622</td><td>0.624</td><td>131-189</td><td>0.132*</td><td>0.126*</td><td>0.007</td></tr><tr><td>Average</td><td>32.250</td><td>7.788</td><td>0.679</td><td>-</td><td>0.190*</td><td>0.148*</td><td>0.050*</td></tr><tr><td colspan="8">Cytoplasmic</td></tr><tr><td>TR1</td><td>15</td><td>7.167</td><td>-</td><td>396-780</td><td>-</td><td>0.361*</td><td>-</td></tr><tr><td>TR2</td><td>3</td><td>1.109</td><td>-</td><td>359-456</td><td>-</td><td>0.175*</td><td>-</td></tr><tr><td>TR3</td><td>4</td><td>2.589</td><td>-</td><td>309-440</td><td>-</td><td>0.389*</td><td>-</td></tr><tr><td>TR4</td><td>7</td><td>2.648</td><td>-</td><td>304-485</td><td>-</td><td>0.392*</td><td>-</td></tr><tr><td>Average</td><td>7.250</td><td>3.37825</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Haplotypes</td><td>45</td><td>3.573</td><td>-</td><td>-</td><td>-</td><td>0.350*</td><td>-</td></tr><tr><td colspan="8"><strong>B. macrocarpa</strong></td></tr><tr><td colspan="8">Nuclear</td></tr><tr><td>1027</td><td>8</td><td>3.511</td><td>0.264</td><td>171-201</td><td>1.000*</td><td>0.688*</td><td>1.000*</td></tr><tr><td>Bvm3</td><td>10</td><td>3.374</td><td>0.219</td><td>096-117</td><td>0.949*</td><td>0.759*</td><td>0.786*</td></tr><tr><td>Caa1</td><td>9</td><td>3.069</td><td>0.180</td><td>128-179</td><td>0.901*</td><td>0.714*</td><td>0.655*</td></tr><tr><td>Gtt1</td><td>5</td><td>2.794</td><td>0.199</td><td>107-121</td><td>0.601*</td><td>0.051</td><td>0.580*</td></tr><tr><td>Sb04</td><td>6</td><td>2.826</td><td>0.192</td><td>167-183</td><td>0.858*</td><td>0.676*</td><td>0.561*</td></tr><tr><td>Sb06</td><td>5</td><td>2.279</td><td>0.167</td><td>146-168</td><td>0.886*</td><td>0.841*</td><td>0.279</td></tr><tr><td>Sb07</td><td>13</td><td>4.063</td><td>0.360</td><td>245-269</td><td>0.910*</td><td>0.300*</td><td>0.870*</td></tr><tr><td>Sb15</td><td>6</td><td>2.469</td><td>0.174</td><td>131-169</td><td>0.886*</td><td>0.842*</td><td>0.279</td></tr><tr><td>Average</td><td>7.750</td><td>3.048</td><td>0.219</td><td>-</td><td>0.905*</td><td>0.588*</td><td>0.769*</td></tr><tr><td colspan="8">Cytoplasmic</td></tr><tr><td>TR1</td><td>3</td><td>2.207</td><td>-</td><td>494-527</td><td>-</td><td>0.748*</td><td>-</td></tr><tr><td>TR2</td><td>1</td><td>1.000</td><td>-</td><td>359</td><td>-</td><td>-</td><td>-</td></tr><tr><td>TR3</td><td>3</td><td>1.210</td><td>-</td><td>374-440</td><td>-</td><td>0.873*</td><td>-</td></tr><tr><td>TR4</td><td>1</td><td>1.000</td><td>-</td><td>364</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Average</td><td>2</td><td>1.354</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Haplotypes</td><td>5</td><td>2.073</td><td>-</td><td>-</td><td>-</td><td>0.743*</td><td>-</td></tr></tbody></table>

A$_{t}$, total number of alleles; A$_{r}$, allelic richness; H$_{E}$, gene diversity; F$_{IT}$, F$_{ST}$, F$_{IS}$: fixation indices following Weir and Cockerham (1984); -, not applicable (haploid data).

Significance of departure from Hardy-Weinberg equilibrium within (F$_{IS}$), over all populations (F$_{IT}$), and genetic differentiation (F$_{ST}$) was tested by 10,000 randomizations (*P* < 0.0001). Bonferroni corrections were applied to the mean nuclear F$_{ST}$ estimates and cytoplasmic loci. Allelic size ranges in base pairs (bp) are also indicated for each locus.

diversity, with an average of 7.750 alleles for nuclear variation and a total number of 5 haplotypes for cytoplasmic polymorphism (see Table 1). The geographic distribution of mtDNA haplotypes within populations can be visualized in Fig. S4.

## Geographic trends in the levels of genetic diversity

Estimates of genetic diversity (A$_{t}$, A$_{r}$P, and H$_{E}$) varied greatly among *B. vulgaris* subsp. *maritima* populations

(Appendix). A striking northward decrease in private nuclear alleles within populations was observed over the sampled area, with the highest private allelic richness being found on both sides of the Strait of Gibraltar and in Moroccan inland ruderal populations (Fig. 3A). Nuclear gene diversity and allelic richness significantly decreased with latitude from the Gibraltar region (H$_{E}$: r² = 0.283, P < 10⁻⁴; A$_{t}$: r² = 0.323, P < 10⁻⁴; Fig. 3B). Similar trends were described using the geographic distance calculated along the coastline (H$_{E}$: r² = 0.211, P < 0.003; A$_{t}$: r² = 0.242, P < 0.001). In contrast, no
<figure><img src="image_4.png" /></figure>

**Figure 3.** Spatial patterns of genetic diversity. (A) Nuclear private allelic richness ($A_rP$) of *B. vulgaris* subsp. *maritima* populations is proportional to the size of the dot mapped onto the population location. (B) Linear regressions of gene diversity ($H_E$) and allelic richness ($A_r$) based on nuclear polymorphism with respect to latitude. Black dots represent genetic diversity values for northern populations, from Brittany to northern Moroccan populations at the Strait of Gibraltar (populations labelled 1–40 in Fig. 1). Gray dots represent values for Moroccan populations at lower latitudes (populations labelled 43 to t). Dashed lines indicate the 95% confidence intervals.

significant trends were observed for Moroccan *B. vulgaris* subsp. *maritima* populations at lower latitudes (Fig. 3B). Finally, regarding cytoplasmic polymorphism, no significant relationships with latitude or coastline distance were observed for any parameter of genetic diversity based on multilocus data (data not shown). Very similar trends of decreasing nuclear gene diversity with latitude were depicted using linear mixed models with loci as random intercept (see Table S2 and Fig. S5).

We then focused on Morocco, a region considered as a hot spot of genetic diversity, to compare coastal (C) and inland ruderal (R) populations. Inland ruderal populations exhibited higher levels of gene diversity ($H_E$, C: 0.675, R: 0.734; $P < 10^{-3}$), allelic richness ($A_r$, C: 5.538, R: 6.400; $P < 10^{-3}$), and mean number of alleles ($A_n$, C: 7.569, R: 8.835; $P < 10^{-2}$) for nuclear data. Inland ruderal populations also showed a higher level of allelic richness for cytoplasmic data ($A_r$, C: 1.830, R: 2.209; $P < 10^{-3}$; $A_r$ *Haplo*, C: 3.349, R: 4.532; $P < 10^{-2}$). However, no significant differences were observed for private allelic richness ($A_rP$, C: 0.066, R: 0.093; $P = 0.255$). Regarding genotypic structure, mean multilocus $F_{IS}$ values marginally differed between inland ruderal populations and coastal

populations (C: 0.030, R: 0.073; $P = 0.057$), and no difference was observed for nuclear multilocus $F_{ST}$ (C: 0.115, R: 0.074; $P = 0.112$).

## Mean levels of genetic differentiation

A strong genetic differentiation among *B. vulgaris* subsp. *maritima* populations was found using mitochondrial minisatellite markers: The mean $F_{ST}$ estimate, based on haplotypes, was 0.350 ($P < 0.01$, Table 1). Highly significant spatial genetic differentiation also occurred for nuclear polymorphism, with consistent single-locus $F_{ST}$ estimates ranging from 0.126 to 0.177 (all at $P < 0.0001$) with a mean multilocus $F_{ST}$ of 0.148 ($P < 0.01$, see Table 1). Among the 3570 comparisons, mean pairwise $F_{ST}$ between populations ranged from 0 to 0.461 (94.64% at $P < 0.05$) and from 0 to 0.973 (44.14% at $P < 0.05$) for nuclear and cytoplasmic data, respectively (Fig. S6). Altogether, these results clearly indicated genetically structured populations for both kinds of markers over the whole sampled area. Among the eight sampled populations of *B. macrocarpa*, nuclear and cytoplasmic differentiation among populations was extremely high, with a
mean significant $F_{ST}$ value of 0.588 and of 0.743 for nuclear and cytoplasmic variation, respectively (Table 1 and Fig. S6).

## Population genetic structure of *B. vulgaris* subsp. *maritima* inferred from spatial multivariate ordination method

The sPCA performed on the cytoplasmic data revealed no significant structure, either global ($P = 0.681$) or local ($P = 0.184$). When performing sPCA on nuclear data, a significant global structure was observed ($P < 0.001$), indicating that neighboring entities are genetically more similar than would be expected at random. The three first axes in the sPCA unambiguously captured all global spatial patterns (see Fig. S7). The resulting synthetic map defined four clearly different geographic entities: (1) the Bay of Biscay; (2) the Portuguese Atlantic coast, with a genetic cline between northern and southern populations; (3) the Strait of Gibraltar region on the Spanish and Moroccan sides; and (4) Moroccan southern Atlantic coast (Fig. 4A). Moreover, Moroccan inland ruderal populations clustered strongly with the group of populations located on the Moroccan side of the Strait of Gibraltar and southern Spain. Interestingly, populations located across the Strait of Gibraltar showed high autocorrelation values, suggesting no major genetic discontinuities in this area (Figs. 4A, S7). This analysis further traced back the evolutionary history of populations sampled on Madeira Island, which appear to originate from southern Portugal and Spain or northern Morocco near Gibraltar. Strongly concordant results were obtained with Bayesian clustering for the second modal value at $K = 6$ (see Fig. S2). Levels of genetic diversity for these genetically distinct groups of populations are depicted by both sPCA and Bayesian clustering in Figure 4B. Significant differences were found among groups for all variables (all at $P < 10^{-2}$). Overall, the general trend was that populations located in Morocco and along the Strait of Gibraltar displayed the highest levels of cytoplasmic and nuclear genetic diversity. Figure 5 illustrates how pairwise genetic differentiation among populations differed between genetic clusters and shows a general increase in nuclear differentiation northward along the Atlantic coastline. Interestingly, cytoplasmic and nuclear data showed similar patterns with higher

pairwise $F_{ST}$ values observed along the Portuguese coastline (Fig. 5).

Lastly, including *B. macrocarpa* individuals in the sPCA clearly supported the high genetic distinctiveness between the two species observed with Bayesian analyses (see Fig. S8).

## Isolation-by-distance patterns

To design Mantel correlograms based on (1) previous results on population memberships and (2) a sufficient number of pairwise comparisons within distance classes to ensure robust statistical support, we only focused on three groups of *B. vulgaris* subsp. *maritima* populations that can be visualized in Figure 6. Using nuclear data, Mantel correlograms for northern populations up (north of the Strait of Gibraltar) revealed a clear continuous decrease in genetic similarity with increasing geographic distance (Fig. 6A). Focusing on Moroccan Atlantic populations, we observed moderately positive *rz* values stabilizing up to 280 km. From this distance, *rz* values substantially decreased (Fig. 6C). Regarding inland ruderal populations for nuclear data, only a positive short-distance autocorrelation and a long-distance genetic differentiation were observed (Fig. 6E). Finally, for cytoplasmic data, no significant spatial trends were observed: The correlograms showed fluctuations indicative of random genetic drift effects (Fig. 6B, D and F).

## Discussion

### Nuclear and cytoplasmic polymorphisms in *B. vulgaris* subsp. *maritima* and *B. macrocarpa*

One aim of this study was to elucidate spatial patterns of population genetic structure and draw indirect inferences on mating system based on genotypic structure, in two related species, *B. vulgaris* subsp. *maritima* and *B. macrocarpa*. One may wonder whether the modal *K* depicting three genetically distinct clusters clearly discriminated the two species. Nonetheless, we are very confident of the species assignation for the following reasons: (1) Concordant patterns of genetic distinctiveness were obtained using both sPCA and Bayesian clustering, (2) the hierarchical genetic

**Figure 4.** Genetic distinctiveness and mean levels of genetic diversity estimates for *B. vulgaris* subsp. *maritima* populations. (A) Genetic structure of *B. vulgaris* subsp. *maritima* populations revealed by a multivariate ordination method (sPCA) as described in Jombart et al. (2008). Color plot synthesizes the three first sPCA principal components at their respective spatial population coordinates. Each of the three axes was plotted on a red, green, and blue color scale, and the intensity of the color is proportional to the spatial principal component scores. Arrows indicate major surface currents redrawn from Koutsikopoulos and LeCann (1996), Millot (1999), Pelegrí et al. (2005), and Peliz et al. (2005). (B) Bar graph of the mean genetic diversity estimates ($H_e$, $A_r$, and $A_r_{Haplo}$) for each genetically distinct group of populations revealed by sPCA and for *B. macrocarpa* populations. Bars with different letters are significantly different at $P < 0.05$.
<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>
<figure><img src="image_7.png" /><figcaption>**Figure 5.** Box-plots of nuclear (A) and cytoplasmic (B) pairwise differentiation of populations ($F_{ST}$) within genetic clusters depicted by sPCA. Box-plots indicate the median (horizontal line), the 25th and 75th percentiles (bottom and the top of the box), and the minimum/ maximum values (vertical dashed lines). Plots with different letters are significantly different at $P < 0.05$.</figcaption></figure>

differentiation was the highest between the two clusters of sea beets and the putative cluster composed of *Beta macrocarpa* individuals, and (3) contrasting patterns of genetic differentiation among populations and genotypic structure within populations were found between the two clusters of sea beets and the cluster composed of *Beta macrocarpa* individuals, the observed patterns being concordant with what is known about the mating system in these two species (Letschert 1993; Bruun et al. 1995) and more generally in plant species displaying different reproductive traits (Hamrick and Godt 1996). Therefore, the third cluster we depicted is clearly unlikely to correspond to an additional cluster of sea beets that would be characterized by complete autogamy. Overall, significant cytoplasmic and nuclear genetic differentiation was found among 85 *B. vulgaris* subsp. *maritima* and 8 *B. macrocarpa* populations. Mitochondrial DNA diversity revealed considerably more population structuring than that found for nuclear DNA, as shown by the strong asymmetry in $F_{ST}$ values. This finding can be explained by (1) the larger effective population size for nuclear genes as compared to organellar genes, indicating a stronger effect of genetic drift on cytoplasmic DNA,

and (2) higher gene flow for nuclear DNA, which occurs via seed and pollen dispersal, than for cytoplasmic DNA, which takes place via seed movements only (Ennos 1994). In *B. vulgaris* subsp. *maritima*, the amount of genetic differentiation was fairly high compared with those found in previous studies involving a metapopulation structure at microgeographic scale less than a few hundred meters (De Cauwer et al. 2010) or involving a regional structure (Fievet et al. 2007). This is in accordance with the fact that genetic differentiation generally increases with geographic distance (reviewed in Petit et al. 2005). However, the contrast between nuclear and cytoplasmic differentiation observed over a wide species-range scale appeared reduced compared with those found in studies performed at finer geographic scales. The large differences between seed and pollen dispersal may indeed be more difficult to detect due to the greater action of genetic drift and mutation, with differentiation being random at large geographic scales (Hutchison and Templeton 1999; Fievet et al. 2007). Furthermore, founder events involving rare propagules can undergo strong genetic drift, as suggested by the lack of large-scale structuring for cytoplasmic diversity (Whitlock and McCauley 1990). This type of metapopulation dynamics is plausible given that sea beets occur in coastal environments where frequent disturbance may lead to extinctions and recolonizations (De Cauwer et al. 2012). This explanation is consistent with hierarchical F-statistics analyses, revealing that most differentiation occurs among populations within geographic groups.

*Beta macrocarpa* populations exhibited very low genetic variation, large heterozygote deficiency, and a high degree of genetic differentiation compared with *B. vulgaris* subsp. *maritima* populations. Herbaceous annuals with seed dispersal by gravity and mixed-mating or pure selfing mating systems usually exhibit higher $F_{IS}$ values and stronger genetic differentiation than outcrossing species (Hamrick and Godt 1996; Duminil et al. 2009). Accordingly, this suggested that selfing could be the main reproductive mode of *B. macrocarpa*. This is therefore consistent with a self-compatible mating system documented by Buttler (1977), Lange and Debock (1989), Letschert (1993), and Bruun et al. (1995). Interestingly, two populations located in Santo Porto showed mean single-locus $F_{IS}$ values and multilocus genotypes indicative of individual fixed heterozygosity, as also documented by Letschert (1993) on tetraploid ($2n = 36$) accessions from Canary Islands. Along with the recent discovery of diploid and tetraploid forms occurring in sympatry in south Portugal (Castro et al. 2013), our findings also suggest the mixing of diploid and tetraploid variants in *B. macrocarpa* located in Madeira Islands. Inbreeding depression is weaker in tetraploids, which may be a determining factor in the evolution of polyploidy (Husband and Schemske 1997). Indeed, poly-
<figure><img src="image_8.png" /><figcaption>**Figure 6.** Mantel correlograms based on nuclear microsatellite (A, C, E) and mitochondrial minisatellite (B, D, F) data. Three groupings were investigated: northern populations (north of the Strait of Gibraltar) (A, B), Moroccan Atlantic populations (C, D), and inland ruderal populations (E, F). Black circles indicate a significant *rz* value at *P* < 0.05.</figcaption></figure>
ploidy appears to be especially prevalent in self-fertilizing annual colonizing plants, and it is often argued that heterozygosity is a key element in their colonizing ability (reviewed in Soltis and Soltis 2000). In *B. macrocarpa*, tetraploidy may thus be adaptively significant because this species often occurs in disturbed environments. Nonethe- less, further experiments involving chromosome counting and flow cytometry are necessary to accurately document the occurrence of diploid and tetraploid variants in sympatry.

**B. vulgaris subsp. maritima: do inland populations have peculiar evolutionary histories?**

To date, there are only few examples of typical coastal species that also occur inland – albeit sparsely – especially halophytic species. Comes and Abbott (1998) observed lower genetic diversity in inland populations of *Senecio gallicus* and suggested that they originated from coastal populations. However, lower genetic diversity for inland populations is not a common feature (e.g., Krüger et al. 2002; Prinz et al. 2009). In our study, Moroccan inland ruderal populations exhibited higher levels of genetic diversity than coastal populations. In addition, no clear isolation-by-distance (IBD) pattern was detected, contrary to spatial genetic structure commonly observed in coastal populations (Fievet et al. 2007). Only positive autocor- relation of short-distance classes was observed, likely reflect- ing genetic clustering and family structuring caused by limited short-range dispersal of the multigerm seedball (De Cauwer et al. 2012). The lack of IBD patterns at lar- ger spatial scales can be attributed – at least partially – by recent spreads along with human-caused habitat distur- bances, which blurs the spatial patterns of genetic related- ness by modifying the movement of dispersal vectors and significantly changing the number and distribution of recruitment sites (González-Martínez et al. 2010). Pollen flow may also be very localized and constrained by inland topography, as opposed to the propitious pollen dispersal in open habitats such as shorelines (Arnaud et al. 2009; De Cauwer et al. 2010).

Finally, Moroccan inland ruderal populations clustered with coastal populations in the Gibraltar region, a result consistent with both sPCA and Bayesian analyses (Figs. 3, S2). Again, these results could highlight a recent expan- sion of *B. vulgaris* subsp. *maritima* from the Mediterra- nean coast, possibly mediated by human-induced long- dispersal events as suggested by Desplanque et al. (1999) and Arnaud et al. (2009) for southern European inland wild beet populations. Ruderal habitats are common ele- ments in the Mediterranean landscape due to land-use changes (Naveh and Vernet 1991). Therefore, opportunis-

tic growth and rapid reproduction in unstable habitats appear to have played an important role in the inland expansion of weedy beet relatives (Arnaud et al. 2011). Inland ruderal Moroccan populations further appeared to have admixed origin from coastal European and Moroc- can lineages. The admixture of different lineages and sub- sequent multiple introductions in a new habitat can provide large amounts of genetic diversity, novel sources of genetic combinations in plant species complexes, and may thus offer opportunities to intensify the adaptability of populations to a new niche on microevolutionary timescales (Ellstrand and Schierenbeck 2000; Barrett et al. 2008). A striking example of recurring pattern of adaptive evolution in *B. vulgaris* is the evolution of invasive weed beets from crop-wild hybrids that entails numerous changes in life history, morphology, and ecology. It has indeed been shown that the invasiveness of weedy lineages arises through the introgression of a wild genetic back- ground into the crop gene pool (Arnaud et al. 2010). This wild genetic background comes from inland ruderal beets found in disturbed habitats in southern France, which have very similar ecology as inland Moroccan beets (Des- planque et al. 1999; Arnaud et al. 2009). Indeed, in southern France, these noncoastal wild beet populations have colonized disturbed, human-influenced habitats, thanks to specific life-history traits that enhance coloniz- ing abilities: the timing of flowering, the seed dormancy associated with a long-lived seed bank, and the break- down of self-incompatibility (Desplanque et al. 1999; Hautekèete et al. 2002; Arnaud et al. 2011). These traits may also be crucial in inland Moroccan populations shar- ing very similar ecological conditions, that is, unpredict- able and man-made environments at low latitude. Overall, in such disturbed-like environments, the combination of reproductive traits and life-history traits due to admix- ture of different lineages may facilitate the adaptive evolu- tion of inland beet lineages and could facilitate their colonization capabilities, as they do for weedy invasive lin- eages found in cultivated fields (Arnaud et al. 2010). In the future, this could be of particular concern in the major areas of cultivation of sugar beets in central west Morocco.

Lastly, founder effects during colonization are likely to involve ample mixing of seeds due to pastoralism and grazing activities along field edges (Pereira et al. 2009). Multiple sources of introduction and/or different succes- sive independent colonization events from the Mediterra- nean coast may therefore explain the admixture, the lack of IBD, and the higher levels of genetic diversity found in inland populations. Newly developed ABC methods could be useful to elucidate the origins of inland populations (e.g., Lombaert et al. 2011). However, our data set does not conform to current ABC computations in two aspects. First, we dealt with a very large number of popu-
lations not amenable to current ABC analyses without size reduction. Second, the among-population genetic structure observed in each compartment (coastal and ruderal) is high, which can lead to the inference of erroneous scenarios. Potential solutions exist to circumvent these difficulties, by modeling unsampled populations and/or using different sets of sampled populations (E. Lombaert, pers. comm.). Testing the appropriateness of these different solutions is, however, beyond the scope of this paper. As such, we acknowledge that concluding on a particular evolutionary history of inland populations remains speculative as long as these scenarios have not been thoroughly tested.

## Phylogeographic patterns in *B. vulgaris* subsp. *maritima*

Understanding the underlying mechanisms of range shifts and the determinants of population establishment are fundamental questions in evolutionary biology and ecological genetics and more importantly in the context of changing climate. The patterns we document in *B. vulgaris* subsp. *maritima* provide clues for understanding how life-history and major geographic features interact to shape the distribution of genetic diversity.

## Does the Strait of Gibraltar impede gene flow?

The Mediterranean-Atlantic transition has been reported to be an efficient barrier to the expansion of many coastal plants and marine species along the coast (Kadereit et al. 2005; Westberg and Kadereit 2009). Nonetheless, this barrier effect remains uncertain in *B. vulgaris* subsp. *maritima* because we found a clear genetic similarity among populations located on both sides of the Strait of Gibraltar. This lack of genetic structuring suggested that there is no major biogeographic barrier in this area, which can be explained by: i) colonization before barrier formation and/or ii) effective dispersal after barrier formation. During the last glacial maxima, the relative sea level dropped 120 m relative to its present level and the distance between the two continents was thereby reduced to 10 km (Yokoyama et al. 2000). In addition to the proximity between shores, small islands and islets emerged and may have provided pathways of migration through a stepping-stone-like network between the two continents (Collina-Girard 2001). Alternatively, it may involve long-distance seed dispersal events favored by sea currents and efficient pollen flow. As found in *B. vulgaris* subsp. *maritima*, genetic continuities on both sides of the strait have been reported for several herbaceous plant species characterized by ruderal life-history traits (e.g., Caujapé-Castells and Jansen 2003; Ortiz et al. 2007). These find-

ings are in sharp contrast with other studies that examined long-lived plants such as trees, which tend to show genetic breaks associated with the strait (e.g., Lumaret et al. 2002). Overall, it appears that the Strait of Gibraltar filtered plant migration according to their dispersal and life-history traits (Lavergne et al. 2013).

## Influence of marine currents

Seed dispersal is a major determinant of the spatial genetic structure of populations (Ennos 1994). In most cases, *B. vulgaris* subsp. *maritima* seeds are dispersed by gravity in the vicinity of the mother plant (De Cauwer et al. 2012). However, secondary seed dispersal through water transport has been suggested because buoyant seeds can be carried away during the highest tides and subsequently drifted in marine waters (Fievet et al. 2007). Interestingly, the main phylogeographic breaks depicted by both sPCA and Bayesian clustering analysis matched the main surface current discontinuities (see Fig. 3). The sharp genetic change observed in Galicia, at the western tip of the Iberian Peninsula, coincides with a water-mass front located off Cape Finisterre and has been reported as a major dispersal barrier for several marine organisms such as algae and marine snails with planktonic larvae being released in seawater (Perez et al. 1993; Alberto et al. 1999; Piñeira et al. 2008). Furthermore, an additional biogeographic barrier has been identified in Portugal at 41°N where the eastern North Atlantic surface currents bifurcate (e.g., Diekmann et al. 2005; Coyer et al. 2011). The resulting south-flowing Portuguese current and the northeast-flowing Portuguese counter-current may thus account for the genetic cline and increased levels of pairwise genetic differentiation observed in *B. vulgaris* subsp. *maritima* populations along the Portuguese Atlantic coast (see Martins et al. 2002). Finally, the substantial genetic distinctiveness between southern Moroccan Atlantic populations and populations found in the Strait of Gibraltar region on the Spanish and Moroccan sides corresponds to a physical discontinuity of marine currents, as shown in Figure 3. Thus, large-scale oceanographic features along the coastline may affect seed dispersal and population connectivity in wild sea beet, leading to genetically differentiated lineages similar to what can be observed in marine organisms (e.g., Alberto et al. 1999; Coyer et al. 2011). Lastly, drifting through seawater can facilitate seed dispersal, which is very limited along the terrestrial shoreline. Indeed, sea beet seeds have no particular dispersal mechanism except a gravity-driven pattern of dispersal. As a consequence, most seeds are dispersed at short distances near the mother plants in terrestrial areas (Arnaud et al. 2011; De Cauwer et al. 2012). A high genetic differentiation among a mosaic of genetically distinct demes is thus found
within populations, with striking genetic discontinuities tak-
ing place within a few tens of meters (De Cauwer et al.
2012). Moreover, although being a wind-pollinated species,
pollen flow is spatially localized and constrained by land-
scape topography (Fievet et al. 2007; De Cauwer et al. 2010,
2012). This explains the strong IBD pattern found for coastal
populations at nuclear microsatellite loci.

Geographic trends in the levels of genetic
diversity and structure: footprints of postglacial
recolonization?

It seems unlikely that *B. vulgaris* subsp. *maritima* survived
the Würm glaciation: Due to ice and permafrost, it would
not have been able to maintain a continuous distribution
along large parts of the Atlantic coast of Europe. Today,
the northern range limit of *B. vulgaris* subsp. *maritima* is
situated at approximately 58°N near the 8°C January and
14°C July isotherms. Therefore, populations probably had
to retreat to at least northern Portugal, where the 14°C
July isotherm would have been located during the last
glacial maximum (see Fig. S9). In accordance, *Eryngium
maritimum*, a coastal species that shares the same habi-
tats, temperature requirements, and distribution range as
*B. vulgaris* subsp. *maritima*, was estimated to have
retreated to as far south as North Africa (Clausing et al.
2000). Because rare alleles are expected to be lost with
bottlenecks occurring during subsequent recolonization
and range expansion, patterns of central-marginal and lat-
itudinal gradients in genetic diversity can point to refuge
sources (Austerlitz et al. 1997). Consistently, a clear and
significant decrease in allelic richness and gene diversity
with increasing latitude was observed from the Gibraltar
region, which exhibited a high level of genetic diversity
and many private alleles at nuclear loci. This pattern of
"southern richness" is characteristic of temperate terres-
trial taxa that persisted through past glacial conditions in
southern refugia (Hampe and Petit 2005). This is also
indicative of a central-marginal trend in decreasing diver-
sity (see Guo 2012) and may corroborate the Gibraltar
region as a historical refuge for *B. vulgaris* subsp. *mariti-
ma*. Indeed, the Gibraltar region is considered as a buffer
zone that was climatically stable during geological times
due to its latitude and the influence of the Atlantic Ocean
and the Mediterranean Sea (Rodríguez-Sánchez et al.
2008). However, we cannot rule out secondary contacts
between two divergent lineages originating from the east-
ern Mediterranean and southern Moroccan coast. The
high levels of admixture found in northern Moroccan
coastal populations are consistent with the occurrence of
genetically diverse suture zones similar to glacial refugia
(Sakaguchi et al. 2011). The high genetic diversity
observed in this area may also suggest past abrupt increase

in migration among populations. Indeed, such transient
dynamics has been shown to leave genetic diversity peaks
that can be conserved over thousands of generations, as
theoretically demonstrated by Alcala et al. (2013).

Subsequently, the significant decrease in nuclear gene
diversity and allelic richness northward from the Strait of
Gibraltar region can be attributed to a series of founder
effects during range expansion along the shoreline. In this
type of narrow expanding population front, drift effects
may be accentuated by growth of genetically uniform
populations ahead of the main colonization front that
would inhibit subsequent gene flow (Waters et al. 2013).
The strong IBD pattern observed for nuclear diversity also
points to a scenario analogous to the propagule-pool
model developed for landscapes with extinctions-recolon-
ization dynamics, that is, founders are drawn from a lim-
ited number of source populations consequently reducing
population genetic diversity and increasing the genetic
differentiation among populations (Wade and McCauley
1988; Whitlock and McCauley 1990; Giles and Goudet
1997; Haag et al. 2005). Indeed, in Europe, sea beets are
generally patchily located along the sea shore, where high
tides, seasonal storms, and anthropogenic disturbances
cause frequent extinctions of populations.

In contrast, we did not find any IBD patterns for cyto-
plasmic variation, nor any evidence for decreasing levels
of cytoplasmic diversity or increase in differentiation
among newly established northern Atlantic populations
following the postglacial recolonization. A migrant-pool
model – where colonization occurs from multiple source
populations (Slatkin 1977) – would be more likely to fit
the spatial patterns of cytoplasmic variation. At first
glance, these results appear puzzling as several sea beet
seeds are aggregated in a single fruit, which provides the
opportunity for kin-structured founder effects (Whitlock
and McCauley 1990; De Cauwer et al. 2012). Similar con-
trasting results between nuclear and cytoplasmic diversity
have been observed in plants (e.g., García-Verdugo et al.
2010) and could be interpreted by the occurrence of mul-
tiple stochastic long-distance seed migration mediated by
major marine currents, as suggested in our study by the
depicted phylogeographic breaks that matched some
major marine current discontinuities. Seed dispersal by
marine currents may indeed possibly compensate the
founder effects by gene mixing despite a narrow coastline
habitat (Fayard et al. 2009).

**Conclusion**

Nearly complete autogamy was verified in *B. macrocarpa*.
Gradients in allelic richness and gene diversity in
*B. v.* subsp. *maritima* suggested a scenario of postglacial
recolonization from the Mediterranean-Atlantic region and
supported the status of southern Iberia and Morocco as a long-term refuge as documented for a wide range of organisms. Further, specific life-history traits along with the admixture from genetically differentiated lineages are likely to increase invasiveness for ruderal lineages in anthropogenically disturbed habitats. Our results also indicate that environmental factors, such as marine currents, may play a role in shaping genetic diversity in coastal plant species concordantly of what can be observed in marine organisms releasing larvae in water flows. Range expansions do not always have straightforward consequences: in our case study, long-distance seed dispersal through marine currents may drive to a mixing of cytoplasmic diversity and buffer the northward decrease in nuclear allelic richness observed for edge populations recently established. This challenges the usual view that populations located at the edge of a species' geographic distribution are genetically depauperate.

## Acknowledgments

We wish to express our gratitude to C. Godé for technical assistance in molecular analyses; S. Gallina for computing assistance with Bayesian clustering; C. Roux for graphical help in R; and S. Villain, F. Roux, N. Saidi, and J. Cuguen for help in sampling. Numerical results presented in this paper were carried out using the European Grid Infrastructure (http://www.egi.eu) with the Biomed virtual organization (http://lsgc.org/en/Biomed:home) via DIRAC portal (http://diracgrid.org) supported by France Grille (http://www.france-grilles.fr/). We thank the technical staff of the European Grid Infrastructure and the supporting National Grid Initiatives for providing technical support and infrastructure. We are also grateful to C.R. Engel, the associate editor, and three anonymous referees for very helpful comments on the manuscript. This work was funded by a grant from the "PRAD-PHC Maroc (09-01)," the "FRB/Région Nord-Pas de Calais" (GENEFRAG project), and a grant from "INRA SPE."

## Data Accessibility

Coordinates of sample locations, cytoplasmic, and nuclear genetic data are available in DRYAD: doi: 10.5061/dryad.m1n38.

## Conflict of Interest

None declared.

## References

Alberto, F., R. Santos, and J. M. Leitão. 1999. Assessing patterns of geographic dispersal of *Gelidium sesquipedale*

(Rhodophyta) through RAPD differentiation of populations. Mar. Ecol. Prog. Ser. 191:101–108.

Alcala, N., D. Streit, J. Goudet, and S. Vuilleumier. 2013. Peak and persistent excess of genetic diversity following an abrupt migration increase. Genetics 193:953–971.

Arnaud, J.-F., S. Fénart, C. Godé, S. Deledicque, P. Touzet, and J. Cuguen. 2009. Fine-scale geographical structure of genetic diversity in inland wild beet populations. Mol. Ecol. 18:3201–3215.

Arnaud, J.-F., S. Fénart, M. Cordellier, and J. Cuguen. 2010. Populations of weedy crop-wild hybrid beets show contrasting variation in mating system and population genetic structure. Evol. Appl. 3:305–318.

Arnaud, J.-F., J. Cuguen, and S. Fénart. 2011. Metapopulation structure and fine-scaled genetic structuring in crop-wild hybrid weed beets. Heredity 107:395–404.

Austerlitz, F., B. Jung-Muller, B. Godelle, and P.-H. Gouyon. 1997. Evolution of coalescence times, genetic diversity and structure during colonization. Theor. Popul. Biol. 51:148–164.

Barrett, S. C. H., R. I. Colautti, and C. G. Eckert. 2008. Plant reproductive systems and evolution during biological invasion. Mol. Ecol. 17:373–383.

Bartsch, D., and N. C. Ellstrand. 1999. Genetic evidence for the origin of Californian wild beets (genus *Beta*). Theor. Appl. Genet. 99:1120–1130.

Biancardi, E., L. W. Panella, and R. T. Lewellen. 2012. *Beta maritima*: the origin of beets. Springer, New York, NY.

Blondel, J. 2006. The 'design' of Mediterranean landscapes: a millennial story of humans and ecological systems during the historic period. Hum. Ecol. 34:713–729.

Bruun, L., A. Haldrup, S. G. Petersen, L. Frese, T. S. M. de Bock, and W. Lange. 1995. Self-incompatibility reactions in wild species of the genus *Beta* and their relation to taxonomical classification and geographical origin. Genet. Resour. Crop Evol. 42:293–301.

Buttler, K. P. 1977. Variation in wild populations of annual beet (*Beta*, Chenopodiaceae). Plant Syst. Evol. 128:123–136.

Castro, S., M. M. Romeiras, M. Castro, M. C. Duarte, and J. Loureiro. 2013. Hidden diversity in wild *Beta* taxa from Portugal: insights from genome size and ploidy level estimations using flow cytometry. Plant Sci. 207:72–78.

Caujapé-Castells, J., and R. K. Jansen. 2003. The influence of the Miocene Mediterranean desiccation on the geographical expansion and genetic variation of *Androcymbium gramineum* (Cav.) McBride (Colchicaceae). Mol. Ecol. 12:1515–1525.

Cavalli-Sforza, L. L., and A. W. F. Edwards. 1967. Phylogenetic analysis: models and estimation procedures. Evolution 21:550–570.

Chen, C., E. Durand, F. Forbes, and O. François. 2007. Bayesian clustering algorithms ascertaining spatial population structure: a new computer program and a comparison study. Mol. Ecol. Notes 7:747–756.
Clausing, G., K. Vickers, and J. W. Kadereit. 2000. Historical biogeography in a linear system: genetic variation of Sea Rocket (*Cakile maritima*) and Sea Holly (*Eryngium maritimum*) along European coasts. Mol. Ecol. 9:1823–1833.

Collina-Girard, J. 2001. L'Atlantide devant le détroit de Gibraltar ? Mythe et géologie. Earth Planet. Sci. 333:233–240.

Comes, H. P., and R. J. Abbott. 1998. The relative importance of historical events and gene flow on the population structure of a Mediterranean ragwort, *Senecio gallicus* (Asteraceae). Evolution 52:355–367.

Coulon, A., J. F. Cosson, N. Morellet, J.-M. Angibault, B. Cargnelutti, M. Galan, et al. 2006. Dispersal is not female biased in a resource-defence mating ungulate, the European roe deer. Proc. R. Soc. B Biol. Sci. 273:341–348.

Coyer, J. A., G. Hoarau, J. F. Costa, B. Hogerdijk, E. A. Serrão, E. Billard, et al. 2011. Evolution and diversification within the intertidal brown macroalgae *Fucus spiralis/F. vesiculosus* species complex in the North Atlantic. Mol. Phylogenet. Evol. 58:283–296.

De Cauwer, I., M. Dufay, J. Cuguen, and J.-F. Arnaud. 2010. Effects of fine-scale genetic structure on male mating success in gynodioecious *Beta vulgaris* ssp. *maritima*. Mol. Ecol. 19:1540–1558.

De Cauwer, I., J.-F. Arnaud, A. Courseaux, and M. Dufay. 2011. Sex-specific fitness variation in gynodioecious *Beta vulgaris* ssp. *maritima*: do empirical observations fit theoretical predictions? J. Evol. Biol. 24:2456–2472.

De Cauwer, I., M. Dufay, B. Hornoy, A. Courseaux, and J.-F. Arnaud. 2012. Gynodioecy in structured populations: understanding fine-scale sex ratio variation in *Beta vulgaris* ssp. *maritima*. Mol. Ecol. 21:834–850.

Desplanque, B., P. Boudry, K. Broomberg, P. Saumitou-Laprade, J. Cuguen, and H. Van Dijk. 1999. Genetic diversity and gene flow between wild, cultivated and weedy forms of *Beta vulgaris* L (Chenopodiaceae), assessed by RFLP and microsatellite markers. Theor. Appl. Genet. 98:1194–1201.

D'Hondt, B., P. Breyne, W. Van Landuyt, and M. Hoffmann. 2012. Genetic analysis reveals human-mediated long-distance dispersal among war cemeteries in *Trifolium micranthum*. Plant Ecol. 213:1241–1250.

Diekmann, O. E., J. A. Coyer, J. Ferreira, J. L. Olsen, W. T. Stam, G. A. Pearson, et al. 2005. Population genetics of *Zostera noltii* along the west Iberian coast: consequences of small population size, habitat discontinuity and near-shore currents. Mar. Ecol. Prog. Ser. 290:89–96.

Duminil, J., O. J. Hardy, and R. J. Petit. 2009. Plant traits correlated with generation time directly affect inbreeding depression and mating system and indirectly genetic structure. BMC Evol. Biol. 9:177.

Durand, E., F. Jay, O. Gaggiotti, and O. François. 2009. Spatial inference of admixture proportions and secondary contact zones. Mol. Biol. Evol. 26:1963–1973.

El Mousadik, A., and R. J. Petit. 1996. High level of genetic differentiation for allelic richness among populations of the

argan tree *Argania spinosa* (L) Skeels endemic to Morocco. Theor. Appl. Genet. 92:832–839.

Ellstrand, N. C., and K. A. Schierenbeck. 2000. Hybridization as a stimulus for the evolution of invasiveness in plants? Proc. Natl Acad. Sci. USA 97:7043–7050.

Ennos, R. A. 1994. Estimating the relative rates of pollen and seed migration among plant-populations. Heredity 72:250–259.

Evanno, G., S. Regnaut, and J. Goudet. 2005. Detecting the number of clusters of individuals using the software STRUCTURE: a simulation study. Mol. Ecol. 14:2611–2620.

Falush, D., M. Stephens, and J. K. Pritchard. 2003. Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies. Genetics 164:1567–1587.

Fayard, J., E. K. Klein, and F. Lefevre. 2009. Long distance dispersal and the fate of a gene from the colonization front. J. Evol. Biol. 22:2171–2182.

Fénart, S., F. Austerlitz, J. Cuguen, and J.-F. Arnaud. 2007. Long distance pollen-mediated gene flow at a landscape level: the weed beet as a case study. Mol. Ecol. 16:3801–3813.

Fénart, S., J.-F. Arnaud, I. De Cauwer, and J. Cuguen. 2008. Nuclear and cytoplasmic genetic diversity in weed beet and sugar beet accessions compared to wild relatives: new insights into the genetic relationships within the *Beta vulgaris* complex species. Theor. Appl. Genet. 116:1063–1077.

Fievet, V., P. Touzet, J.-F. Arnaud, and J. Cuguen. 2007. Spatial analysis of nuclear and cytoplasmic DNA diversity in wild sea beet (*Beta vulgaris* ssp. *maritima*) populations: do marine currents shape the genetic structure? Mol. Ecol. 16:1847–1864.

Frese, L., E. Demeijer, and J. Letschert. 1990. New wild beet genetic resources from Portugal and Spain. Zuckerindustrie 115:950–955.

García-Verdugo, C., A. D. Forrest, M. F. Fay, and P. Vargas. 2010. The relevance of gene flow in metapopulation dynamics of an oceanic island endemic, *Olea europaea* subsp. *guanchica*. Evolution 64:3525–3536.

Giles, B. E., and J. Goudet. 1997. A case study of genetic structure in a plant metapopulation. Pp. 429–454 in I. Hanski and M. E. Gilpin, eds. Metapopulation biology: ecology, genetics, and evolution. Academic Press, San Diego, CA.

González-Martínez, S. C., M. Dubreuil, M. Riba, G. G. Vendramin, F. Sebastiani, and M. Mayol. 2010. Spatial genetic structure of *Taxus baccata* L. in the western Mediterranean Basin: past and present limits to gene movement over a broad geographic scale. Mol. Phylogenet. Evol. 55:805–815.

Goudet, J. 1995. FSTAT (Version 1.2): a computer program to calculate *F*-statistics. J. Hered. 86:485–486.

Goudet, J. 2005. HIERFSTAT, a package for R to compute and test hierarchical *F*-statistics. Mol. Ecol. Notes 5:184–186.
Guo, Q. 2012. Incorporating latitudinal and central-marginal trends in assessing genetic variation across species ranges. Mol. Ecol. 21:5396-5403.

Haag, C. R., M. Riek, J. W. Hottinger, V. I. Pajunen, and D. Ebert. 2005. Genetic diversity and genetic differentiation in *Daphnia* metapopulations with subpopulations of known age. Genetics 170:1809-1820.

Hampe, A., and R. J. Petit. 2005. Conserving biodiversity under climate change: the rear edge matters. Ecol. Lett. 8:461-467.

Hamrick, J. L., and M. J. W. Godt. 1996. Effects of life history traits on genetic diversity in plant species. Philos. Trans. R. Soc. Lond. B Biol. Sci. 351:1291-1298.

Hautekèete, N. C., Y. Piquot, and H. Van Dijk. 2002. Life span in *Beta vulgaris* ssp. *maritima*: the effects of age at first reproduction and disturbance. J. Ecol. 90:508-516.

Hewitt, G. M. 1999. Post-glacial recolonization of European biota. Biol. J. Linn. Soc. 68:87-112.

Hewitt, G. M. 2000. The genetic legacy of the Quaternary ice ages. Nature 405:907-913.

Hewitt, G. M. 2004. Genetic consequences of climatic oscillations in the Quaternary. Philos. Trans. R. Soc. Lond. B Biol. Sci. 359:183-195.

Husband, B. C., and D. W. Schemske. 1997. The effect of inbreeding in diploid and tetraploid populations of *Epilobium angustifolium* (Onagraceae): implications for the genetic basis of inbreeding depression. Evolution 51:737-746.

Hutchison, D. W., and A. R. Templeton. 1999. Correlation of pairwise genetic and geographic distance measures: inferring the relative influences of gene flow and drift on the distribution of genetic variability. Evolution 53:1898-1914.

Jakobsson, M., and N. A. Rosenberg. 2007. CLUMPP: a cluster matching and permutation program for dealing with label switching and multimodality in analysis of population structure. Bioinformatics 23:1801-1806.

Jombart, T. 2008. ADEGENET: a R package for the multivariate analysis of genetic markers. Bioinformatics 24:1403-1405.

Jombart, T., S. Devillard, A. B. Dufour, and D. Pontier. 2008. Revealing cryptic spatial patterns in genetic variability by a new multivariate method. Heredity 101:92-103.

Jombart, T., D. Pontier, and A. B. Dufour. 2009. Genetic markers in the playground of multivariate analysis. Heredity 102:330-341.

Jung, C., K. Pillen, L. Frese, S. Fahr, and A. E. Melchinger. 1993. Phylogenetic relationships between cultivated and wild species of the genus *Beta* revealed by DNA fingerprinting. Theor. Appl. Genet. 86:449-457.

Kadereit, J. W., R. Arafeh, G. Somogyi, and E. Westberg. 2005. Terrestrial growth and marine dispersal? Comparative phylogeography of five coastal plant species at a European scale. Taxon 54:861-876.

Kadereit, G., S. Hohmann, and J. W. Kadereit. 2006. A synopsis of Chenopodiaceae subfam. Betoideae and notes on the taxonomy of *Beta*. Willdenowia 36:9-19.

Kishima, Y., T. Mikami, A. Hirai, M. Sugiura, and T. Kinoshita. 1987. *Beta* chloroplast genomes: analysis of fraction I protein and chloroplast DNA variation. Theor. Appl. Genet. 73:330-336.

Koutsikopoulos, C., and B. LeCann. 1996. Physical processes and hydrological structures related to the Bay of Biscay anchovy. Scientia Marina 60:9-19.

Krüger, A. M., F. H. Hellwig, and C. Oberprieler. 2002. Genetic diversity in natural and anthropogenic inland populations of salt-tolerant plants: random amplified polymorphic DNA analyses of *Aster tripolium* L. (Compositae) and *Salicornia ramosissima* Woods (Chenopodiaceae). Mol. Ecol. 11:1647-1655.

Lange, W., and T. S. M. Debock. 1989. The diploidised meiosis of tetraploid *Beta macrocarpa* and its possible application in breeding sugar-beet. Plant Breeding 103:196-206.

Lavergne, S., A. Hampe, and J. Arroyo. 2013. In and out of Africa: how did the Strait of Gibraltar affect plant species migration and local diversification? J. Biogeogr. 40:24-36.

Letschert, J. P. W. 1993. *Beta* section *Beta*: biogeographical patterns of variation and taxonomy. Wageningen Agricultural University papers 93:1-137.

Lombaert, E., T. Guillemaud, C. E. Thomas, L. J. L. Handley, J. Li, S. Wang, et al. 2011. Inferring the origin of populations introduced from a genetically structured native range by approximate Bayesian computation: case study of the invasive ladybird *Harmonia axyridis*. Mol. Ecol. 20:4654-4670.

Lumaret, R., C. Mir, H. Michaud, and V. Raynal. 2002. Phylogeographical variation of chloroplast DNA in holm oak (*Quercus ilex* L.). Mol. Ecol. 11:2327-2336.

Martins, C. S., M. Hamann, and A. F. G. Fiúza. 2002. Surface circulation in the eastern North Atlantic, from drifters and altimetry. J. Geophys. Res., 107:1-22.

McGrath, J. M., D. Trebbi, A. Fenwick, L. Panella, B. Schulz, V. Laurent, et al. 2007. An open-source first-generation molecular genetic map from a sugarbeet x table beet cross and its extension to physical mapping. Crop Sci. 47:S27-S44.

Médail, F., and P. Quézel. 1997. Hot-spots analysis for conservation of plant biodiversity in the Mediterranean basin. Ann. Mo. Bot. Gard. 84:112-127.

Menozzi, P., A. Piazza, and L. Cavallisforza. 1978. Synthetic maps of human gene-frequencies in Europeans. Science 201:786-792.

Migliore, J., A. Baumel, M. Juin, and F. Médail. 2012. From Mediterranean shores to central Saharan mountains: key phylogeographical insights from the genus *Myrtus*. J. Biogeogr. 39:942-956.

Millot, C. 1999. Circulation in the Western Mediterranean Sea. J. Mar. Syst. 20:423-442.
Mörchen, M., J. Cuguen, G. Michaelis, C. Hanni, and P. Saumitou-Laprade. 1996. Abundance and length polymorphism of microsatellite repeats in *Beta vulgaris* L. Theor. Appl. Genet. 92:326-333.

Naveh, Z., and J. L. Vernet. 1991. The palaeohistory of the Mediterranean biota. Pp. 19-32 *in* R. H. Groves, F. Di Castri, eds. Biogeography of Mediterranean invasions, Cambridge University Press, Cambridge, U.K.

Nishizawa, S., T. Kubo, and T. Mikami. 2000. Variable number of tandem repeat loci in the mitochondrial genomes of beets. Curr. Genet. 37:34-38.

Oden, N. L., and R. R. Sokal. 1986. Directional autocorrelation – an extension of spatial correlograms to two dimensions. Syst. Zool. 35:608–617.

Ortiz, M. Á., K. Tremetsberger, S. Talavera, T. Stuessy, and J. L. García-Castaño. 2007. Population structure of *Hypochaeris salzmanniana* DC. (Asteraceae), an endemic species to the Atlantic coast on both sides of the Strait of Gibraltar, in relation to Quaternary sea level changes. Mol. Ecol. 16:541–552.

Owen, F. V. 1942. Inheritance of cross- and self-sterility and self-fertility in *Beta vulgaris*. J. Agric. Res. 64:0679-0698.

Pelegrí, J. L., J. Arístegui, L. Cana, M. González-Dávila, A. Hernández-Guerra, S. Hernández-León, et al. 2005. Coupling between the open ocean and the coastal upwelling region off northwest Africa: water recirculation and offshore pumping of organic matter. J. Mar. Syst. 54:3-37.

Peliz, Á., J. Dubert, A. M. P. Santos, P. B. Oliveira, and B. Le Cann. 2005. Winter upper ocean circulation in the Western Iberian basin - fronts, eddies and poleward flows: an overview. Deep Sea Res. I 52:621-646.

Pereira, F., S. Queirós, L. Gusmão, I. J. Nijman, E. Cuppen, J. A. Lenstra, et al. 2009. Tracing the History of Goat Pastoralism: new clues from Mitochondrial and Y chromosome DNA in North Africa. Mol. Biol. Evol. 26:2765-2773.

Perez, F. F., C. Mouriño, F. Fraga, and A. F. Rios. 1993. Displacement of water masses and remineralization rates off the Iberian peninsula by nutrients anomalies. J. Mar. Res. 51:869-892.

Petit, R. J., J. Duminil, S. Fineschi, A. Hampe, D. Salvini, and G. G. Vendramin. 2005. Comparative organization of chloroplast, mitochondrial and nuclear diversity in plant populations. Mol. Ecol. 14:689-701.

Piñeira, J., H. Quesada, E. Rolán-Alvarez, and A. Caballero. 2008. Genetic discontinuity associated with an environmentally induced barrier to gene exchange in the marine snail *Littorina saxatilis*. Mar. Ecol. Prog. Ser. 357:175-184.

Prinz, K., K. Weising, and I. Hensen. 2009. Genetic structure of coastal and inland populations of the annual halophyte *Suaeda maritima* (L.) dumort. in Central Europe, inferred from amplified fragment length polymorphism markers. Plant Biol. 11:812-820.

Pritchard, J. K., M. Stephens, and P. Donnelly. 2000. Inference of population structure using multilocus genotype data. Genetics 155:945-959.

R Development Core Team. 2011. R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. http://www.R-project.org; ISBN: 3-900051-07-0.

Raymond, M., and F. Rousset. 1995. GENEPOP (Version-1.2) - Population-genetics software for exact tests and ecumenicism. J. Hered. 86:248-249.

Rice, W. R. 1989. Analyzing tables of statistical tests. Evolution 43:223-225.

Richards, C. M., M. Brownson, S. E. Mitchell, S. Kresovich, and L. Panella. 2004. Polymorphic microsatellite markers for inferring diversity in wild and domesticated sugar beet (*Beta vulgaris*). Mol. Ecol. Notes, 4:243-245.

Rodríguez-Sánchez, F., R. Pérez-Barrales, F. Ojeda, P. Vargas, and J. Arroyo. 2008. The Strait of Gibraltar as a melting pot for plant biodiversity. Quatern. Sci. Rev. 27:2100-2117.

Rosenberg, N. A. 2004. DISTRUCT: a program for the graphical display of population structure. Mol. Ecol. Notes 4:137-138.

Rosenberg, M. S., and C. D. Anderson. 2011. PASSaGE: pattern analysis, spatial statistics and geographic exegesis. version 2. Methods Ecol. Evol. 2:229-232.

Sakaguchi, S., Y. Takeuchi, M. Yamasaki, S. Sakurai, and Y. Isagi. 2011. Lineage admixture during postglacial range expansion is responsible for the increased gene diversity of *Kalopanax septemlobus* in a recently colonised territory. Heredity 107:338-348.

Slatkin, M. 1977. Gene flow and genetic drift in a species subject to frequent local extinctions. Theor. Popul. Biol. 12:253-262.

Smouse, P. E., J. C. Long, and R. R. Sokal. 1986. Multiple regressions and correlation extensions of the Mantel test of matrix correspondence. Syst. Zool. 35:627-632.

Soltis, P. S., and D. E. Soltis. 2000. The role of genetic and genomic attributes in the success of polyploids. Proc. Natl Acad. Sci. USA 97:7051-7057.

Szpiech, Z. A., M. Jakobsson, and N. A. Rosenberg. 2008. ADZE: a rarefaction approach for counting alleles private to combinations of populations. Bioinformatics 24:2498-2504.

Viard, F., J. Bernard, and B. Desplanque. 2002. Crop-weed interactions in the *Beta vulgaris* complex at a local scale: allelic diversity and gene flow within sugar beet fields. Theor. Appl. Genet. 104:688-697.

Villain, S. 2007. *Histoire évolutives de la section Beta: mise en évidence des phénomènes d'hybridation et de spéciation au sein de la Section dans le bassin Méditerranéen* PhD Thesis, Université Lille 1-Sciences et Technologies.

Wade, M. J., and D. E. McCauley. 1988. Extinction and recolonization: their effects on the genetic differentiation of local populations. Evolution 42:995-1005.
Waters, J. M., C. I. Fraser, and G. M. Hewitt. 2013. Founder takes all: density-dependent processes structure biodiversity. TREE 28:78–85.

Weir, B. S., and C. C. Cockerham. 1984. Estimating *F*-statistics for the analysis of population structure. Evolution 38:1358–1370.

Westberg, E., and J. W. Kadereit. 2009. The influence of sea currents, past disruption of gene flow and species biology on the phylogeographical structure of coastal flowering plants. J. Biogeogr. 36:1398–1410.

Whitlock, M. C., and D. E. McCauley. 1990. Some population genetic consequences of colony formation and extinction - genetic correlations within founding groups. Evolution 44:1717–1724.

Yang, R. C. 1998. Estimating hierarchical *F*-statistics. Evolution 52:950–956.

Yokoyama, Y., K. Lambeck, P. De Deckker, P. Johnston, and I. K. Fifield. 2000. Timing of the Last Glacial Maximum from observed sea-level minima. Nature 406:713–716.

## Supporting Information

Additional Supporting Information may be found in the online version of this article:

**Figure S1.** Assignment results from Bayesian clustering following Pritchard et al. (2000) performed on (A) *B. vulgaris* subsp. *maritima* and *B. macrocarpa* individuals (B) *B. vulgaris* subsp. *maritima* individuals only; and, (C) following Durand et al. (2009) on both species. STRUC-TURE analyses (A,B): mean (±SD) probabilities of the data Ln Pr(X|K) over 15 replicated runs plotted as a function of the putative number of clusters *K* (blue dots) and the standardized second-order rate of change of Ln Pr(X|K), Δ*K*, as a function of *K* (red dots). TESS analyses (C): average Deviation Index Criterion (DIC) over 20 replicated runs plotted against *K*.

**Figure S2.** Assignment probabilities of membership of (A) *B. vulgaris* subsp. *maritima* individuals into the five inferred clusters for the second modal value, and (B) *B. vulgaris* subsp. *maritima* and *B. macrocarpa* individuals into the six inferred clusters for the second modal value. Each individual is represented by a thin horizontal line (y axis) partitioned into coloured segments that represented the individual's estimated membership coefficients (x axis) (C) Map of mean population membership probabilities for the six clusters.

**Figure S3.** Population genetic structure of *B. vulgaris* subsp. *maritima* inferred from TESS analyses assuming

*K* = 5. The individual estimates membership coefficients for each cluster are shown for 20 independent runs.

**Figure S4.** Geographical distribution of minisatellites haplotypes within populations.

**Figure S5.** Linear regressions of genetic diversity indices based on nuclear (*A*<sub>r</sub>, *H*<sub>E</sub>, *A*<sub>r</sub>*P*) and based on cytoplasmic (*A*<sub>r</sub>) polymorphism with respect to latitude and coastline distance for each locus. Regressions were performed on (A) northern populations (labelled 1 to 40 in Fig. 1), and (B) Moroccan populations at lower latitudes (labelled 43 to *t*).

**Figure S6.** Matrices of population-pairwise differentiation values (*F*<sub>ST</sub>) of *B. vulgaris* subsp. *maritima* and *B. macrocarpa* individuals represented for (A) nuclear and (B) cytoplasmic data.

**Figure S7.** Individual representation of the three first global axes of the spatial principal component analysis (sPCA) performed on *B. vulgaris* subsp. *maritima* populations. According to the first three global eigenvalues, First (A), second (B), and third (C) sPCA scores, are represented on each plot as squares whose size is proportional to the value of the score, so that the maximum differentiation is between large white squares and large black squares.

**Figure S8.** Individual representation of the three first global axes of the spatial principal component analysis (sPCA) performed on *B. vulgaris* subsp. *maritima* and *Beta macrocarpa* populations.

**Figure S9.** Modern (12°C, 14°C and 16°C mean July, northern Europe) and reconstructed (June/July/August isotherms, southern Europe) Last Glacial Maximum (LGM) isotherms from Kadereit et al. 2005; Bold lines indicate the ice shield during the LGM and dotted lines indicate the coastline during the LGM.

**Table S1.** Mantel tests carried out on genetic distance matrices based on nuclear and cytoplasmic data for *B. vulgaris* subsp. *maritima* populations and either (1) Euclidian geographical distance matrices, (2) matrices of geographical distance measured along the coastline or (3) geographical distance matrices determined through a neighbourhood graph.

**Table S2.** Linear mixed models with loci as random intercept testing the relationship between genetic diversity parameters (based on nuclear and cytoplasmic polymorphism) and explanatory variables (latitude and coastline distance). Models were performed on two datasets to take into account the genetic discontinuities indicated as G1 for northern populations (labelled 1 to 40 in Fig. 1), and G2 for Moroccan populations at lower latitudes (labelled 43 to *t*).
Appendix

Population localities, label, geographic coordinates, sample sizes, and summary of genetic diversity for nuclear polymorphism at eight microsatellite loci and for cytoplasmic polymorphism at four minisatellite loci among 93 wild *Beta vulgaris* subsp. *maritima* (Mari) and *Beta macrocarpa* (Macro) populations. *B. macrocarpa* populations are shaded in gray. *H*<sub>E</sub>, gene diversity within populations; *A*<sub>r</sub>, allelic richness; *A*<sub>r</sub>*P*, private allelic richness; *A*<sub>n</sub>, mean number of alleles; *F*<sub>IS</sub>, the intrapopulation fixation index (*P* < 0.0001); *Nb Haplo*, number of haplotypes; *A*<sub>r</sub>*Haplo*, haplotypic richness. *A*<sub>r</sub> is the allelic richness corrected to *n* = 8 using the rarefaction method of El Mousadik and Petit (1996) implemented in FSTAT version 2.9.3.2.

<table><thead><tr><th rowspan="2">Region</th><th rowspan="2">Label</th><th rowspan="2">Sampling locality</th><th rowspan="2">Latitude</th><th rowspan="2">Longitude</th><th rowspan="2">n</th><th rowspan="2">Species</th><th colspan="5">Nuclear</th><th colspan="3">Cytoplasmic</th></tr><tr><th>H<sub>E</sub></th><th>A<sub>r</sub></th><th>A<sub>r</sub>P</th><th>A<sub>n</sub></th><th>F<sub>IS</sub></th><th>Nb Haplo</th><th>A<sub>r</sub> Haplo</th><th>A<sub>r</sub></th></tr></thead><tbody><tr><td rowspan="6">France</td><td>1</td><td>Bretagne Aber Ildut</td><td>48.46672</td><td>-4.75677</td><td>15</td><td>Mari</td><td>0.537</td><td>4.072</td><td>0.001</td><td>5.000</td><td>0.158</td><td>2</td><td>2.000</td><td>1.500</td></tr><tr><td>2</td><td>Bretagne Pors Gored</td><td>48.46492</td><td>-4.76597</td><td>20</td><td>Mari</td><td>0.516</td><td>3.711</td><td>0.000</td><td>4.500</td><td>0.079</td><td>3</td><td>2.870</td><td>1.913</td></tr><tr><td>3</td><td>St Gilles-Croix-de-Vie</td><td>46.70918</td><td>-1.98033</td><td>11</td><td>Mari</td><td>0.715</td><td>4.790</td><td>0.091</td><td>5.125</td><td>-0.048</td><td>5</td><td>4.803</td><td>2.244</td></tr><tr><td>4</td><td>St Gilles-Croix-de-Vie</td><td>46.69930</td><td>-1.97313</td><td>11</td><td>Mari</td><td>0.660</td><td>4.142</td><td>0.053</td><td>4.375</td><td>-0.034</td><td>1</td><td>1.000</td><td>1.000</td></tr><tr><td>5</td><td>Yves</td><td>46.03008</td><td>-1.05638</td><td>20</td><td>Mari</td><td>0.761</td><td>5.938</td><td>0.009</td><td>8.000</td><td>0.085</td><td>6</td><td>4.813</td><td>2.008</td></tr><tr><td>6</td><td>Fouras</td><td>45.95538</td><td>-1.06833</td><td>20</td><td>Mari</td><td>0.692</td><td>5.350</td><td>0.003</td><td>7.875</td><td>-0.011</td><td>9</td><td>6.882</td><td>2.805</td></tr><tr><td rowspan="6">Iberian peninsula</td><td>7</td><td>Île Madame</td><td>45.95347</td><td>-1.11120</td><td>20</td><td>Mari</td><td>0.729</td><td>5.384</td><td>0.010</td><td>7.250</td><td>0.045</td><td>5</td><td>3.929</td><td>1.954</td></tr><tr><td>8</td><td>Aquitaine Contis plage</td><td>44.08869</td><td>-1.32183</td><td>14</td><td>Mari</td><td>0.604</td><td>4.158</td><td>0.004</td><td>4.750</td><td>0.069</td><td>3</td><td>2.976</td><td>1.498</td></tr><tr><td>9</td><td>Lekeitio</td><td>43.36063</td><td>-2.49835</td><td>24</td><td>Mari</td><td>0.637</td><td>4.629</td><td>0.006</td><td>6.375</td><td>0.149</td><td>4</td><td>3.375</td><td>1.619</td></tr><tr><td>10</td><td>Oyembre</td><td>43.38925</td><td>-4.32932</td><td>39</td><td>Mari</td><td>0.578</td><td>4.444</td><td>0.001</td><td>8.250</td><td>0.148</td><td>4</td><td>2.892</td><td>1.609</td></tr><tr><td>11</td><td>Verdicio</td><td>43.62828</td><td>-5.87637</td><td>26</td><td>Mari</td><td>0.556</td><td>3.668</td><td>0.001</td><td>5.250</td><td>0.082</td><td>2</td><td>1.994</td><td>1.249</td></tr><tr><td>12</td><td>Verdicio</td><td>43.62478</td><td>-5.87912</td><td>15</td><td>Mari</td><td>0.386</td><td>2.430</td><td>0.003</td><td>2.750</td><td>0.118</td><td>2</td><td>1.963</td><td>1.246</td></tr><tr><td>13</td><td>Xago</td><td>43.60868</td><td>-5.91265</td><td>15</td><td>Mari</td><td>0.643</td><td>4.512</td><td>0.002</td><td>5.500</td><td>0.129</td><td>2</td><td>2.000</td><td>1.250</td></tr><tr><td>14</td><td>Celtigos</td><td>43.72675</td><td>-7.79440</td><td>16</td><td>Mari</td><td>0.619</td><td>4.126</td><td>0.001</td><td>4.750</td><td>0.056</td><td>3</td><td>2.516</td><td>1.613</td></tr><tr><td>15</td><td>Celtigos</td><td>43.72542</td><td>-7.79452</td><td>15</td><td>Mari</td><td>0.583</td><td>3.979</td><td>0.033</td><td>4.750</td><td>0.139</td><td>3</td><td>2.790</td><td>1.674</td></tr><tr><td>16</td><td>Espasante</td><td>43.72255</td><td>-7.81105</td><td>20</td><td>Mari</td><td>0.678</td><td>4.763</td><td>0.006</td><td>6.625</td><td>-0.033</td><td>3</td><td>2.646</td><td>1.926</td></tr><tr><td>17</td><td>Texueiro</td><td>42.80988</td><td>-9.11333</td><td>18</td><td>Mari</td><td>0.517</td><td>3.974</td><td>0.067</td><td>5.000</td><td>0.015</td><td>2</td><td>1.698</td><td>1.189</td></tr><tr><td>18</td><td>San Francisco</td><td>42.75742</td><td>-9.07678</td><td>14</td><td>Mari</td><td>0.598</td><td>4.327</td><td>0.008</td><td>6.125</td><td>0.111*</td><td>3</td><td>2.853</td><td>2.198</td></tr><tr><td>19</td><td>San Francisco</td><td>42.76038</td><td>-9.06135</td><td>21</td><td>Mari</td><td>0.628</td><td>4.591</td><td>0.000</td><td>5.500</td><td>-0.081</td><td>4</td><td>3.861</td><td>2.250</td></tr><tr><td>20</td><td>Paço</td><td>41.75697</td><td>-8.87668</td><td>18</td><td>Mari</td><td>0.543</td><td>3.290</td><td>0.003</td><td>3.750</td><td>-0.137*</td><td>2</td><td>1.996</td><td>1.250</td></tr><tr><td>21</td><td>Areosa</td><td>41.72123</td><td>-8.87067</td><td>17</td><td>Mari</td><td>0.432</td><td>3.104</td><td>0.003</td><td>3.750</td><td>-0.072</td><td>3</td><td>2.758</td><td>2.347</td></tr><tr><td>22</td><td>Areosa</td><td>41.72217</td><td>-8.86568</td><td>17</td><td>Mari</td><td>0.407</td><td>2.800</td><td>0.050</td><td>3.500</td><td>-0.04</td><td>2</td><td>2.000</td><td>1.500</td></tr><tr><td>23</td><td>Praia da Tocha</td><td>40.33647</td><td>-8.83982</td><td>13</td><td>Mari</td><td>0.734</td><td>5.030</td><td>0.000</td><td>6.000</td><td>0.069</td><td>4</td><td>3.723</td><td>2.186</td></tr><tr><td>24</td><td>Ericeira</td><td>38.96197</td><td>-9.41917</td><td>10</td><td>Mari</td><td>0.675</td><td>4.184</td><td>0.000</td><td>4.750</td><td>0.053</td><td>5</td><td>4.833</td><td>2.676</td></tr><tr><td>25</td><td>Ericeira</td><td>38.96460</td><td>-9.43267</td><td>15</td><td>Mari</td><td>0.735</td><td>5.097</td><td>0.000</td><td>5.375</td><td>0.064</td><td>3</td><td>2.968</td><td>1.499</td></tr><tr><td>26</td><td>Lizandro</td><td>38.94468</td><td>-9.41478</td><td>15</td><td>Mari</td><td>0.713</td><td>5.382</td><td>0.000</td><td>7.000</td><td>-0.028</td><td>5</td><td>4.509</td><td>2.415</td></tr></tbody></table>
Appendix. Continued.

<table><thead><tr><th rowspan="2">Region</th><th rowspan="2">Label</th><th rowspan="2">Sampling locality</th><th rowspan="2">Latitude</th><th rowspan="2">Longitude</th><th rowspan="2">n</th><th rowspan="2">Species</th><th colspan="5">Nuclear</th><th colspan="3">Cytoplasmic</th></tr><tr><th>H<sub>E</sub></th><th>A<sub>r</sub></th><th>A<sub>r</sub>P</th><th>A<sub>n</sub></th><th>F<sub>IS</sub></th><th>Nb</th><th>Haplo</th><th>A<sub>r</sub> Haplo</th><th>A<sub>r</sub></th></tr></thead><tbody><tr><td rowspan="19">Morocco</td><td>27</td><td>Samoqueira</td><td>37.86950</td><td>-8.79295</td><td>15</td><td>Mari</td><td>0.715</td><td>4.988</td><td>0.025</td><td>6.125</td><td>-0.095*</td><td>1</td><td>1.000</td><td>1.212</td></tr><tr><td>28</td><td>Porto Covo</td><td>37.84960</td><td>-8.79315</td><td>15</td><td>Mari</td><td>0.728</td><td>6.045</td><td>0.067</td><td>7.875</td><td>0.01</td><td>6</td><td>5.450</td><td>2.815</td></tr><tr><td>29</td><td>Pessegueiro</td><td>37.82962</td><td>-8.79122</td><td>10</td><td>Mari</td><td>0.756</td><td>6.400</td><td>0.009</td><td>6.875</td><td>0.103*</td><td>5</td><td>5.000</td><td>2.995</td></tr><tr><td>30</td><td>Carvoeiro</td><td>37.09582</td><td>-8.47155</td><td>26</td><td>Mari</td><td>0.684</td><td>5.511</td><td>0.074</td><td>8.500</td><td>0.161*</td><td>4</td><td>3.919</td><td>1.988</td></tr><tr><td>31</td><td>Chiclana de la Frontera</td><td>36.42903</td><td>-6.15942</td><td>15</td><td>Mari</td><td>0.821</td><td>7.081</td><td>0.099</td><td>9.250</td><td>0.107*</td><td>5</td><td>4.576</td><td>2.348</td></tr><tr><td>32</td><td>Barbate</td><td>36.19350</td><td>-5.90813</td><td>21</td><td>Mari</td><td>0.748</td><td>5.940</td><td>0.010</td><td>8.750</td><td>-0.02</td><td>2</td><td>1.646</td><td>1.510</td></tr><tr><td>33</td><td>Zahara</td><td>36.13983</td><td>-5.85178</td><td>15</td><td>Mari</td><td>0.752</td><td>6.394</td><td>0.046</td><td>8.250</td><td>0.047</td><td>5</td><td>4.713</td><td>1.953</td></tr><tr><td>34</td><td>Tarifa</td><td>36.0095</td><td>-5.60633</td><td>22</td><td>Mari</td><td>0.833</td><td>7.329</td><td>0.097</td><td>11.500</td><td>-0.03</td><td>6</td><td>5.298</td><td>2.844</td></tr><tr><td>35</td><td>Estepona</td><td>36.41867</td><td>-5.17850</td><td>32</td><td>Mari</td><td>0.699</td><td>5.042</td><td>0.061</td><td>7.625</td><td>0.071*</td><td>2</td><td>2.000</td><td>1.250</td></tr><tr><td>36</td><td>Almayate</td><td>36.72283</td><td>-4.14383</td><td>32</td><td>Mari</td><td>0.816</td><td>6.325</td><td>0.005</td><td>9.625</td><td>0.071*</td><td>4</td><td>3.775</td><td>1.961</td></tr><tr><td>37</td><td>Cap Mazarí</td><td>35.53830</td><td>-5.21467</td><td>24</td><td>Mari</td><td>0.744</td><td>6.099</td><td>0.076</td><td>9.500</td><td>0.019</td><td>4</td><td>3.120</td><td>1.864</td></tr><tr><td>38</td><td>Cabo Negro</td><td>35.66600</td><td>-5.28470</td><td>24</td><td>Mari</td><td>0.653</td><td>4.824</td><td>0.049</td><td>6.875</td><td>0.198*</td><td>3</td><td>2.866</td><td>1.476</td></tr><tr><td>39</td><td>Ksar-as-Seghir</td><td>35.84418</td><td>-5.55490</td><td>24</td><td>Mari</td><td>0.753</td><td>5.651</td><td>0.077</td><td>7.750</td><td>-0.016</td><td>2</td><td>2.000</td><td>1.250</td></tr><tr><td>40</td><td>Ksar-as-Seghir</td><td>35.82791</td><td>-5.65704</td><td>33</td><td>Mari</td><td>0.692</td><td>4.421</td><td>0.226</td><td>6.500</td><td>-0.033</td><td>3</td><td>1.858</td><td>1.602</td></tr><tr><td>41</td><td>Cap Malabata</td><td>35.79852</td><td>-5.75094</td><td>26</td><td>Macro</td><td>0.045</td><td>1.234</td><td>0.000</td><td>1.250</td><td>1.000*</td><td>1</td><td>1.000</td><td>1.000</td></tr><tr><td>42</td><td>North of Asilah</td><td>35.53147</td><td>-6.00637</td><td>25</td><td>Mari</td><td>0.724</td><td>6.146</td><td>0.017</td><td>9.000</td><td>0.012</td><td>3</td><td>2.786</td><td>1.460</td></tr><tr><td>43</td><td>Asilah</td><td>35.47471</td><td>-6.02639</td><td>15</td><td>Mari</td><td>0.683</td><td>5.332</td><td>0.037</td><td>6.875</td><td>0.127*</td><td>4</td><td>3.801</td><td>2.158</td></tr><tr><td>44</td><td>Asilah</td><td>35.46424</td><td>-6.04214</td><td>19</td><td>Mari</td><td>0.790</td><td>5.841</td><td>0.279</td><td>7.375</td><td>0.01</td><td>1</td><td>1.000</td><td>1.000</td></tr><tr><td>45</td><td>Larache</td><td>35.19856</td><td>-6.15395</td><td>34</td><td>Mari</td><td>0.713</td><td>5.431</td><td>0.120</td><td>8.750</td><td>0.006</td><td>4</td><td>3.254</td><td>1.585</td></tr><tr><td>46</td><td>Mehdiya harbor</td><td>34.26198</td><td>-6.66486</td><td>10</td><td>Mari</td><td>0.573</td><td>3.609</td><td>0.004</td><td>3.625</td><td>-0.025</td><td>1</td><td>1.000</td><td>1.000</td></tr><tr><td>47</td><td>Mehdiya beach</td><td>34.24847</td><td>-6.68020</td><td>15</td><td>Mari</td><td>0.522</td><td>3.925</td><td>0.000</td><td>4.750</td><td>0.043</td><td>3</td><td>2.963</td><td>1.741</td></tr><tr><td>48</td><td>Bouknadel</td><td>34.14674</td><td>-6.73757</td><td>31</td><td>Mari</td><td>0.636</td><td>5.399</td><td>0.008</td><td>8.125</td><td>0.096</td><td>5</td><td>3.652</td><td>1.699</td></tr><tr><td>49</td><td>Rabat</td><td>34.03328</td><td>-6.83532</td><td>24</td><td>Mari</td><td>0.679</td><td>5.541</td><td>0.090</td><td>8.500</td><td>0.041</td><td>4</td><td>3.914</td><td>2.234</td></tr><tr><td>50</td><td>Rabat</td><td>33.95541</td><td>-6.92240</td><td>14</td><td>Mari</td><td>0.692</td><td>6.664</td><td>0.027</td><td>8.625</td><td>-0.021</td><td>7</td><td>6.425</td><td>3.153</td></tr><tr><td>51</td><td>Temara</td><td>33.92765</td><td>-6.96061</td><td>17</td><td>Mari</td><td>0.701</td><td>6.361</td><td>0.069</td><td>8.750</td><td>0.082*</td><td>5</td><td>4.168</td><td>2.035</td></tr><tr><td>52</td><td>Skhirat beach</td><td>33.89583</td><td>-7.00148</td><td>12</td><td>Mari</td><td>0.643</td><td>5.697</td><td>0.091</td><td>6.500</td><td>0.087</td><td>3</td><td>2.933</td><td>1.737</td></tr><tr><td>53</td><td>Casablanca</td><td>33.57103</td><td>-7.70927</td><td>15</td><td>Mari</td><td>0.692</td><td>6.305</td><td>0.032</td><td>8.125</td><td>0.036</td><td>5</td><td>4.372</td><td>1.920</td></tr><tr><td>54</td><td>Casablanca</td><td>33.53502</td><td>-7.78873</td><td>16</td><td>Mari</td><td>0.654</td><td>6.102</td><td>0.090</td><td>8.250</td><td>0.104*</td><td>4</td><td>3.515</td><td>1.863</td></tr><tr><td>55</td><td>Tamaris Beach</td><td>33.52242</td><td>-7.84612</td><td>17</td><td>Mari</td><td>0.684</td><td>6.064</td><td>0.010</td><td>8.625</td><td>0.054</td><td>3</td><td>2.991</td><td>1.499</td></tr><tr><td>56</td><td>Rit Reima</td><td>33.47979</td><td>-7.94681</td><td>13</td><td>Mari</td><td>0.702</td><td>6.361</td><td>0.159</td><td>7.750</td><td>0.083*</td><td>4</td><td>3.723</td><td>1.956</td></tr><tr><td>57</td><td>El Jadida</td><td>33.25606</td><td>-8.50058</td><td>23</td><td>Mari</td><td>0.653</td><td>5.782</td><td>0.032</td><td>8.750</td><td>0.019</td><td>5</td><td>4.096</td><td>2.055</td></tr></tbody></table>
Appendix. Continued.

<table><thead><tr><th rowspan="2">Region</th><th rowspan="2">Label</th><th rowspan="2">Sampling locality</th><th rowspan="2">Latitude</th><th rowspan="2">Longitude</th><th rowspan="2">n</th><th rowspan="2">Species</th><th colspan="5">Nuclear</th><th colspan="3">Cytoplasmic</th></tr><tr><th>H<sub>E</sub></th><th>A<sub>r</sub></th><th>A<sub>r</sub>P</th><th>A<sub>n</sub></th><th>F<sub>IS</sub></th><th>Nb</th><th>Haplo</th><th>A<sub>r</sub> Haplo</th><th>A<sub>r</sub></th></tr></thead><tbody><tr><td>58</td><td>Sidi Moussa</td><td></td><td>32.54697</td><td>-9.27294</td><td>15</td><td>Mari</td><td>0.617</td><td>5.413</td><td>0.060</td><td>6.875</td><td>-0.013</td><td>3</td><td>2.952</td><td>1.992</td></tr><tr><td>59</td><td>Safi</td><td></td><td>32.32353</td><td>-9.26135</td><td>25</td><td>Mari</td><td>0.660</td><td>5.797</td><td>0.003</td><td>8.250</td><td>0.09</td><td>6</td><td>5.590</td><td>2.413</td></tr><tr><td>60</td><td>El Mehattat</td><td></td><td>31.82972</td><td>-9.54357</td><td>24</td><td>Mari</td><td>0.700</td><td>6.497</td><td>0.030</td><td>9.500</td><td>0.018</td><td>4</td><td>3.802</td><td>1.996</td></tr><tr><td>61</td><td>Moulay-Bouzerktoun</td><td></td><td>31.64687</td><td>-9.66897</td><td>21</td><td>Mari</td><td>0.689</td><td>6.279</td><td>0.053</td><td>9.500</td><td>0.069*</td><td>7</td><td>5.951</td><td>2.462</td></tr><tr><td>62</td><td>Essaouira</td><td></td><td>31.49162</td><td>-9.76452</td><td>26</td><td>Mari</td><td>0.716</td><td>6.665</td><td>0.081</td><td>11.250</td><td>0.122*</td><td>6</td><td>5.396</td><td>2.126</td></tr><tr><td>63</td><td>Immesouanne Beach</td><td></td><td>30.84024</td><td>-9.82156</td><td>20</td><td>Mari</td><td>0.524</td><td>3.337</td><td>0.021</td><td>4.250</td><td>-0.113</td><td>1</td><td>1.000</td><td>1.000</td></tr><tr><td>64</td><td>Immesouanne Beach</td><td></td><td>30.83945</td><td>-9.82058</td><td>20</td><td>Mari</td><td>0.647</td><td>4.310</td><td>0.004</td><td>5.125</td><td>-0.067*</td><td>4</td><td>3.993</td><td>2.150</td></tr><tr><td>65</td><td>Tarhazoute</td><td></td><td>30.56927</td><td>-9.74449</td><td>18</td><td>Macro</td><td>0.100</td><td>1.735</td><td>0.000</td><td>1.750</td><td>0.926*</td><td>1</td><td>1.000</td><td>1.000</td></tr><tr><td>66</td><td>North of Tarhazoute</td><td></td><td>30.56915</td><td>-9.74450</td><td>30</td><td>Macro</td><td>0.062</td><td>1.470</td><td>0.235</td><td>1.750</td><td>0.723*</td><td>1</td><td>1.000</td><td>1.000</td></tr><tr><td>67</td><td>Imourane Beach</td><td></td><td>30.50862</td><td>-9.68645</td><td>4</td><td>Mari</td><td>0.766</td><td>-</td><td>-</td><td>4.375</td><td>-0.184*</td><td>-</td><td>-</td><td>-</td></tr><tr><td>68</td><td>Aglou Beach</td><td></td><td>29.80457</td><td>-9.83363</td><td>20</td><td>Mari</td><td>0.687</td><td>5.217</td><td>0.101</td><td>7.375</td><td>-0.146</td><td>-</td><td>-</td><td>-</td></tr><tr><td>A</td><td>Tétouan</td><td></td><td>35.56475</td><td>-5.59577</td><td>24</td><td>Mari</td><td>0.766</td><td>5.808</td><td>0.090</td><td>8.375</td><td>0.109*</td><td>4</td><td>3.374</td><td>1.772</td></tr><tr><td>B</td><td>Khemisset</td><td></td><td>33.86032</td><td>-5.86353</td><td>20</td><td>Mari</td><td>0.727</td><td>5.998</td><td>0.019</td><td>8.375</td><td>-0.049</td><td>7</td><td>5.586</td><td>2.792</td></tr><tr><td>B</td><td>Khemisset</td><td></td><td>33.86032</td><td>-5.86353</td><td>4</td><td>Macro</td><td>-</td><td>-</td><td>-</td><td>2.875</td><td>0.657*</td><td>1</td><td>-</td><td>-</td></tr><tr><td>C</td><td>Tiflét</td><td></td><td>33.83300</td><td>-6.13092</td><td>24</td><td>Mari</td><td>0.762</td><td>6.843</td><td>0.163</td><td>11.500</td><td>0.071*</td><td>5</td><td>4.458</td><td>2.136</td></tr><tr><td>D</td><td>Tiflét</td><td></td><td>33.88387</td><td>-6.2798</td><td>24</td><td>Mari</td><td>0.673</td><td>5.289</td><td>0.000</td><td>7.250</td><td>0.080*</td><td>3</td><td>2.245</td><td>1.547</td></tr><tr><td>E</td><td>Rabat Center</td><td></td><td>34.01871</td><td>-6.82016</td><td>13</td><td>Mari</td><td>0.726</td><td>6.351</td><td>0.120</td><td>7.750</td><td>0.192*</td><td>6</td><td>5.570</td><td>2.455</td></tr><tr><td>F</td><td>Bouznika</td><td></td><td>33.77795</td><td>-7.23266</td><td>21</td><td>Mari</td><td>0.729</td><td>6.489</td><td>0.031</td><td>9.750</td><td>0.118*</td><td>2</td><td>1.998</td><td>1.476</td></tr><tr><td>G</td><td>Azemmour</td><td></td><td>33.28874</td><td>-8.33113</td><td>19</td><td>Mari</td><td>0.691</td><td>5.734</td><td>0.021</td><td>8.000</td><td>0.076*</td><td>5</td><td>4.316</td><td>1.864</td></tr><tr><td>H</td><td>Sebt-des-Oulad-Hassine</td><td></td><td>33.11637</td><td>-8.44159</td><td>2</td><td>Mari</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td><td>-</td></tr><tr><td>H</td><td>Sebt-des-Oulad-Hassine</td><td></td><td>33.11637</td><td>-8.44159</td><td>6</td><td>Macro</td><td>0</td><td>-</td><td>-</td><td>1.000</td><td>-</td><td>1</td><td>-</td><td>-</td></tr><tr><td>I</td><td>Settat</td><td></td><td>33.01755</td><td>-7.62048</td><td>33</td><td>Mari</td><td>0.787</td><td>6.826</td><td>0.021</td><td>11.750</td><td>0.151*</td><td>9</td><td>6.15</td><td>2.880</td></tr><tr><td>J</td><td>El Guelb</td><td></td><td>32.59028</td><td>-7.81497</td><td>4</td><td>Mari</td><td>0.620</td><td>-</td><td>-</td><td>3.375</td><td>-0.109</td><td>1</td><td>-</td><td>-</td></tr><tr><td>J</td><td>El Guelb</td><td></td><td>32.59028</td><td>-7.81497</td><td>37</td><td>Macro</td><td>0.073</td><td>1.770</td><td>0.062</td><td>2.375</td><td>0.716*</td><td>1</td><td>1.000</td><td>1.000</td></tr><tr><td>K</td><td>Benguerir</td><td></td><td>32.09525</td><td>-7.95204</td><td>15</td><td>Mari</td><td>0.722</td><td>6.887</td><td>0.182</td><td>9.250</td><td>-0.022</td><td>7</td><td>6.466</td><td>2.452</td></tr><tr><td>L</td><td>Marrakech</td><td></td><td>31.62862</td><td>-8.11320</td><td>18</td><td>Mari</td><td>0.806</td><td>7.253</td><td>0.087</td><td>10.375</td><td>0.070*</td><td>5</td><td>4.954</td><td>2.187</td></tr></tbody></table>
Appendix. Continued.

<table><thead><tr><th rowspan="2">Region</th><th rowspan="2">Label</th><th rowspan="2">Sampling locality</th><th rowspan="2">Latitude</th><th rowspan="2">Longitude</th><th rowspan="2">n</th><th rowspan="2">Species</th><th colspan="5">Nuclear</th><th colspan="3">Cytoplasmic</th></tr><tr><th>H<sub>E</sub></th><th>A<sub>r</sub></th><th>A<sub>r</sub>P</th><th>A<sub>n</sub></th><th>F<sub>IS</sub></th><th>Nb Haplo</th><th>A<sub>r</sub> Haplo</th><th>A<sub>r</sub></th></tr></thead><tbody><tr><td></td><td>M</td><td>Zoudia</td><td>31.62891</td><td>-8.26613</td><td>19</td><td>Mari</td><td>0.744</td><td>6.629</td><td>0.164</td><td>9.750</td><td>0.115*</td><td>7</td><td>6.051</td><td>2.611</td></tr><tr><td></td><td>N</td><td>Oudaya</td><td>31.59095</td><td>-8.43793</td><td>11</td><td>Mari</td><td>0.764</td><td>6.294</td><td>0.010</td><td>6.625</td><td>0.091</td><td>4</td><td>3.993</td><td>2.250</td></tr><tr><td></td><td>O</td><td>Chichaoua East</td><td>31.55745</td><td>-8.72626</td><td>14</td><td>Mari</td><td>0.767</td><td>6.654</td><td>0.072</td><td>8.625</td><td>-0.013</td><td>4</td><td>3.933</td><td>2.213</td></tr><tr><td></td><td>P</td><td>Chichaoua Center</td><td>31.54526</td><td>-8.75929</td><td>18</td><td>Mari</td><td>0.779</td><td>7.559</td><td>0.100</td><td>11.500</td><td>0.101*</td><td>8</td><td>6.890</td><td>2.671</td></tr><tr><td></td><td>Q</td><td>Chichaoua</td><td>31.53588</td><td>-8.76460</td><td>29</td><td>Mari</td><td>0.748</td><td>7.170</td><td>0.090</td><td>12.625</td><td>0.081*</td><td>6</td><td>4.696</td><td>2.231</td></tr><tr><td></td><td>R</td><td>South of Essaouira</td><td>31.39083</td><td>-9.71322</td><td>14</td><td>Mari</td><td>0.634</td><td>4.357</td><td>0.125</td><td>5.000</td><td>-0.029</td><td>3</td><td>2.823</td><td>1.470</td></tr><tr><td></td><td>S</td><td>Immesouanne</td><td>30.84051</td><td>-9.78695</td><td>19</td><td>Mari</td><td>0.754</td><td>6.664</td><td>0.021</td><td>9.250</td><td>0.048</td><td>5</td><td>4.559</td><td>2.597</td></tr><tr><td></td><td>T</td><td>North of Guelmin</td><td>29.26903</td><td>-9.73838</td><td>16</td><td>Mari</td><td>0.739</td><td>6.406</td><td>0.355</td><td>8.750</td><td>0.006</td><td>4</td><td>3.516</td><td>2.159</td></tr><tr><td rowspan="2">Santo Porto Island</td><td>69</td><td>Serra de Dentra</td><td>33.08833</td><td>-16.3102</td><td>16</td><td>Macro</td><td>0.535</td><td>3.375</td><td>0.001</td><td>3.375</td><td>-0.461*</td><td>2</td><td>2.000</td><td>1.250</td></tr><tr><td>70</td><td>Adega das Levadas</td><td>33.03647</td><td>-16.37742</td><td>17</td><td>Macro</td><td>0.528</td><td>4.051</td><td>0.000</td><td>4.125</td><td>0.234*</td><td>4</td><td>3.996</td><td>1.749</td></tr><tr><td rowspan="3">Madeira Island</td><td>71</td><td>Porto Moniz</td><td>32.86643</td><td>-17.1663</td><td>21</td><td>Mari</td><td>0.605</td><td>3.656</td><td>0.017</td><td>4.375</td><td>0.115*</td><td>1</td><td>1.000</td><td>1.000</td></tr><tr><td>72</td><td>Ponta do Pargo</td><td>32.81427</td><td>-17.26288</td><td>24</td><td>Mari</td><td>0.702</td><td>4.121</td><td>0.000</td><td>4.750</td><td>-0.01</td><td>2</td><td>1.815</td><td>1.215</td></tr><tr><td>73</td><td>Caniço de Baixa</td><td>32.64597</td><td>-16.82443</td><td>18</td><td>Mari</td><td>0.745</td><td>5.242</td><td>0.006</td><td>7.000</td><td>0.194*</td><td>5</td><td>4.596</td><td>1.924</td></tr><tr><td colspan="2">Mean multilocus estimates</td><td></td><td></td><td></td><td></td><td>Mari</td><td>0.678</td><td>5.367</td><td>0.052</td><td>7.239</td><td>0.050*</td><td>3.965</td><td>3.573</td><td>1.915</td></tr><tr><td colspan="2"></td><td></td><td></td><td></td><td></td><td>Macro</td><td>0.223</td><td>2.272</td><td>0.050</td><td>2.313</td><td>0.769*</td><td>1.500</td><td>1.666</td><td>1.167</td></tr></tbody></table>