君，已阅读到文档的结尾了呢~~
分子标记技术及其研究进展,分子标记技术,分子标记,分子标记辅助选择,分子标记辅助选择育种,issr分子标记,分子标记辅助育种,ssr分子标记,分子标记方法,srap分子标记
扫扫二维码，随身浏览文档
手机或平板扫扫即可继续访问
分子标记技术及其研究进展
举报该文档为侵权文档。
举报该文档含有违规或不良信息。
反馈该文档无法正常浏览。
举报该文档为重复文档。
推荐理由：
将文档分享至：
分享完整地址
文档地址：
粘贴到BBS或博客
flash地址：
支持嵌入FLASH地址的网站使用
html代码：
&embed src='/DocinViewer--144.swf' width='100%' height='600' type=application/x-shockwave-flash ALLOWFULLSCREEN='true' ALLOWSCRIPTACCESS='always'&&/embed&
450px*300px480px*400px650px*490px
支持嵌入HTML代码的网站使用
您的内容已经提交成功
您所提交的内容需要审核后才能发布，请您等待！
3秒自动关闭窗口PMCID: PMC4768523Phylogeography of a widespread species: pre-glacial vicariance, refugia, occasional blocking straits and long-distance migrations,1,* ,2 ,1 and
22Área de Ecolog , Facultad de Ciencias, Campus A Zapateira, Universidad de A Coruña, 15071 A Coruña, Spain*Corresponding author's e-mail address: Contributed by Associate Editor: F. Xavier Picó
This article has been
other articles in PMC.Phylogeographic studies give us the opportunity to reconstruct the historical migrations of species and link them with climatic and geographic variation. They are, therefore, a key tool to understanding the relationships among biology, geology and history. One of the most interesting biogeographical areas of the world is the Mediterranean region. However, in this area, the description of concordant phylogeographic patterns is quite scarce, which limits the understanding of evolutionary patterns related to climate. Species with one-dimensional distribution ranges, such as the strawberry tree (Arbutus unedo), are particularly useful to unravel these patterns. Here, we describe its phylogeographic structure and check for concordance with patterns seen in other Mediterranean plants: longitudinal/latitudinal clines of diversity, evidence for glacial refugia and the role of sea straits in dispersal. We also identify the most likely source for the disjunct Irish population. With this aim, we sequenced four chloroplast non-coding fragments of A. unedo from 23 populations covering its whole distribution. We determined the genetic diversity, population structure, haplotype genealogy and time to the most recent common ancestor. The genealogy revealed two clades that separated during the last 700 ky but before the last glacial maximum. One clade occupies Atlantic Iberia and North Africa, while the other occurs in the Western Mediterranean. The Eastern Mediterranean is inhabited by newer haplotypes derived from both clades, while the Irish population is closely related to Iberian demes. The straits of Sicily and Gibraltar partially restricted the gene flow. We concluded that a vicariance event during the Late Quaternary in the western end of the species' range followed by eastward migration seems a likely explanation for the observed phylogeographic pattern. The role of straits indicates an occasional communication between Europe and North Africa, suggesting that the latter was a novel refugia. The East–West genetic split in Iberia is consistent with the refugia-within-refugia model. Finally, the strawberry tree possibly reached Ireland from Iberia instead of throughout the maritime fringe of France as previously thought.Keywords: Arbutus unedo, clades, cpDNA, Lusitanian, Mediterranean, refugia, strait, vicariancePhylogeography uses molecular information to infer the historical forces that have shaped the genealogical architecture of populations and closely related species (). The typical way to do so is by reference to the framework of genealogical concordance (). In particular, multispecies concordance has stimulated a number of comparative phylogeographic surveys for faunas and floras with mixed success (). One region where the description of multispecies concordant patterns is particularly difficult is the Mediterranean basin, one of the earth's 25 biodiversity hotspots, with 4.3 % of the world's 300 000 plant species (). The basin serves as a contact zone for three continents separated by a tortuous sea, resulting in a very heterogeneous region in terms of climate and geography (), which is mirrored by an intricate phylogeography. Thus, a recent review of the phylogeographic studies conducted in the region concluded that, compared with the Alps or North America, the Mediterranean lacks largely common patterns across plant species (). Moreover, the phylogeographic complexity of the region might also be a consequence of a blurring of genetic footprints over time. Thus, the less drastic effects of the Pleistocene glaciations in the Mediterranean area are thought to have facilitated an accumulation of species responses to the successive palaeoenvironmental changes ().Despite the above, commonalities still exist, even though they frequently show some inconsistencies. Similarly to other areas of the world, there is a south–north decrease in genetic diversity (). This latitudinal cline must have been caused by a leading-edge expansion from southern refugia, and many studies have corroborated the role of the three large Mediterranean peninsulas (Iberia, Italy and Balkans) as refugia for the survival of species and engines of speciation (, ). Moreover, it is now widely recognized that the ‘refugia-within-refugia’ model () explains the phylogeographic breaks identified within these peninsulas, particularly in Iberia (; ). Nonetheless, it is growingly acknowledged that refugia were not confined to major peninsulas and many of them have been identified in areas that were previously attributed a lesser role (large islands, North Africa, Turkey and Catalonia-Provence). Again, this complex arrangement of refugia possibly evidences the cumulative effects of historical and environmental factors that occurred since the Tertiary ().Another interesting pattern is the finding that genetic diversity increases from west to east along the Mediterranean in many species (). Several Mediterranean trees even show a clear division between the eastern and the western ends of the Mediterranean, with eastern lineages commonly pre-dating western ones (; ). These longitudinal patterns have been attributed to relatively recent processes such as stronger demographic impacts under local last glacial maximum (LGM) climate in west populations and an east–west recolonization during the Holocene (). Alternatively, the east–west vicariance has also been interpreted as a genetic footprint of much older, pre-Quaternary range dynamics (; ; , ). Finally, the effectiveness of the sea straits (like Gibraltar or Sicily) as phylogeographic barriers has varied over time given the eustatic sea-level shifts during the Pleistocene () and even the same strait (e.g. Gibraltar) served as barrier or not depending on the species (; ).Given its phylogeographic complexity, the study of simplified systems such as organisms with one-dimensional distribution range may help in the search for congruence in the Mediterranean basin (; ). This is the case of Arbutus unedo (strawberry tree), a small tree with a neat circum-Mediterranean distribution (; ) that occupies a narrow coastal fringe from Tunisia to Morocco and from Spain to Turkey (Fig. ). From a historical perspective,
cited the genus Arbutus as a component of the Madrean–Tethyan flora that existed between western North America and Eurasia until the end of the Oligocene (25 Mya). The divergence between Mediterranean and North American Arbutus was estimated in at least 21.2 Mya (). The ancestral A. unedo remained in the area currently occupied by Hungary–Poland–Bulgaria until the Messinian (), A. andrachne being the closest relative. Interestingly, the strawberry tree also lives outside the Mediterranean bioclimatic region along the Atlantic coasts of Morocco, Iberian Peninsula, France and Southwest Ireland. In the Atlantic, the strawberry tree shows the discontinuous distribution typically found in members of the ‘Lusitanian flora’, a group of species that have puzzled biogeographers because they occur in Ireland and northern Iberia while they are largely absent from intervening countries (; ). Early interpretations viewed Lusitanian species as relicts of the Tertiary flora that survived through the glacial period (; ), but a post-glacial colonization of Ireland from southern refugia is currently considered a more likely alternative. Accordingly, Iberia has been described as the likely origin for Irish oaks (; ), whereas the precise location of the putative southern refugia remains elusive for other species (, ). Some Lusitanian species may even have recolonized Ireland with human help rather than naturally (; ). In the particular case of the strawberry tree, its presence in Ireland since at least 4000 years ago has been confirmed by both pollen () and wood remains (). Accordingly, it might have reached Ireland in the mid-Holocene thermal optimum, 9000&#x years ago (; ).
Population frequencies of the cpDNA lineages across the range (greyed) of A. unedo.In this study, we sequenced four chloroplastidial non-coding fragments to address a number of specific issues: (i) to estimate the time to the most recent common ancestor (TMRCA) of A. unedo and to infer possible colonization routes along the Mediterranean and the Atlantic façade, (ii) to test whether some patterns detected in other Mediterranean plants also apply to the strawberry tree (longitudinal/latitudinal clines of diversity, placement of glacial refugia and role of the straits as barriers to dispersal) and (iii) to identify the likely source for the disjunct Irish population. To our knowledge, this is the first range-wide study of the phylogeographic structure of A. unedo.The speciesArbutus unedo is a small–medium tree from the family Ericaceae, subfamily Arbutoideae. Genus Arbutus includes several members in USA and four in Europe: A. andrachne, A. pavarii, A. canariensis and A. unedo (). Flowering spans from September to December and fecundation is entomophilus (). Fruits take 12 months to ripen and seeds are dispersed mainly by birds (). Regarding soil conditions, the strawberry tree prefers siliceous or decarbonated substrata, although it also occurs on formations of carbonated substrata, on heavy clay soils, sandstone and limestone, with pH ranging from 4 to 7 (; ). The distribution of the species is limited by a mean January temperature of 40 °F (4.4 °C) ().Sample collection and DNA extractionThe sampling plan covered the whole range of A. unedo and included leaf samples from 23 wild populations in nine countries (Table ). Additionally, we collected samples from two individuals of A. andrachne, the closest relative of A. unedo, from Çanakkale (Turkey). At each population, leaves were collected from 12 adult trees separated at least 15 m from each other to minimize the risk of sampling relatives. Leaves were wrapped in Kimtech Science wipes (Kimberly-Clark Europe Ltd, UK) and stored in silica gel until DNA extraction. DNA was extracted with the Realpure Genomic DNA extraction kit (REAL) according to manufacturer’s instructions except that we extended the incubation times for cell lysis (2 h at 55 °C) and protein precipitation (20 min at −20 °C). Extractions were conducted in batches of 23 samples plu each batch combined individuals from at least five populations. The quality of extracted DNA and negative controls was checked on 1.5 % agarose gels.
Arbutus unedo sampling localities, UTM coordinates, sample size (n), number of segregating sites (S), number of total (Ht) and private (Hp) haplotypes, haplotype (hd) and nucleotide (π × 100) diversities (with SD) obtained based on the ...SequencingAfter reviewing previous phylogeographic studies, initial trials were conducted with five non-coding cpDNA regions commonly used with trees and, specifically, in studies of closely related species: trnS(GCU)-trnG(UUC), trnH(GUG)-psbA, trnL(UAG)-rpl32, trnT(UGU)-trnL(UAA) and Intrón L (; , ). These initial trials revealed a widespread occurrence of illegible sequences due to long poly-A fragments in non-coding regions (polymerase slippage). Sequencing problems were solved and clear sequences were obtained for four out of the five non-coding regions, all except trnS(GCU)-trnG(UUC), by using a high-fidelity polymerase (Phusion High-Fidelity ADN Polymerase, Thermo Scientific) (). Polymerase chain reactions (PCRs) were performed using 10 µL of 2× Phusion Master Mix (with 1 U of Phusion High-Fidelity ADN Polymerase), 0.5 µmol L−1 of each primer and 1 µL of DNA (diluted 1 : 10). Amplification conditions were 30 s at 98 °C for DNA 30 cycles of 10 s at 98 °C, 20 s at the specific annealing temperature of each primer and 30 s at 72 °C; and a final extension of 7 min at 72 °C. Annealing temperatures were 60.8 °C (trnT-trnL), 59.7 °C (trnL-rpl32), 61.8 °C (Intron L) and 66.8 °C (trnH-psbA). Before sequencing, PCR products were checked on 1.5 % agarose gels and purified with 1 µL of Exonuclease I (20 U µL−1) and 2 µL of FastAP (10 U µL−1) (Fermentas, Waltham, MA, USA). Purified PCR products were sequenced under BigDye Terminator cycling conditions on an Automatic Sequencer 3730XL (Applied Biosystems, USA). The program Geneious R6 v.6.1.4 was used to check the quality of the chromatograms and to perform alignments. Singleton polymorphisms and unique haplotypes (i.e. detected in one single individual) were corroborated with reverse sequencing reactions to discard sequencing artefacts. As the non-recombinant nature of cpDNA makes it equivalent to a single-locus marker, sequences from the four fragments were combined into a single haplotype for every individual. These sequences can be obtained from GenBank using the accession numbers: BankIt 1868954: KU194976–KU195222 (Intron L fragment); BankIt 1874447: KU205359–KU205584 (trnT-tnrL fragment); BankIt 1874452: KU205585–KU205820 (trnL-rpl32 fragment); BankIt 1874463: KU205821–KU206025 (trnH-psbA fragment). It is noteworthy that the identification of each population is represented by the last two letters and two numbers of the sequence definition name.
shows the equivalence between the population codes used in this article and the population codes that identify the populations in GenBank. Still, the identity of each fragment was retained when we calculated the TMRCA to allow for differences in the mutation rate of each region.Data analysisHaplotypes were identified with DnaSP v.5 (), which also provided estimates for haplotype diversity (Hd) (). Genetic diversity was evaluated as nucleotide diversity (π) using Arlequin 3.5 (). Arlequin was also used to carry out an analysis of molecular variance (AMOVA) to assess population structure and to calculate the Fu's FS test to detect any evidence of demographic expansion (). Population structure was further investigated by defining clusters of genetically similar populations with a spatial analysis of molecular variance (SAMOVA) (). Spatial analysis of molecular variance group populations into a user-defined number of groups (K) using a simulated annealing procedure to maximize the variance (FCT) among those groups. The analysis was performed in SAMOVA 2.0 () running 10 000 iterations from 100 initial conditions, testing K = 2–9 with and without constraint for the geographic composition of the groups. The final number of groups was chosen for which FCT began to plateau. In addition, Permut 1.0 was used to estimate population differentiation by calculating GST and NST which are two measures of populational differentiation, GST use only of the allelic frequencies and NST also take into account the similarities between the haplotypes. The analysis was performed under the assumption that a significantly higher NST value suggests the existence of phylogeographic structure ().The genealogy of the haplotypes was inferred with the median-joining network algorithm and maximum parsimony calculation implemented in Network 4.6 (; ). We ran the program with ε = 0 and mutations weighted following recommended guidelines: 10 as a default value, 5 for very common mutations (where the various alternative stages occurred in similar proportions) and 20 for indels. To improve the interpretation of the resulting tree, we removed the resulting loops, but always without modifying the interpretation. Complementarily, we estimated the genealogy using TCS 1.21 (), setting the connection limit at 95 % and considering the gaps (missing value) as fifth state.Estimates of the TMRCA were conducted with the 207 A. unedo sequences in BEAST 2.1 () by combining two searches with 50 million Markov chain Monte Carlo (MCMC) each and a sample frequency of 1000. The first 5 million generations (10 %) were discarded as a burn-in. We used jModelTest 2.1 () to determine with Akaike and the Bayesian information criteria the simplest model of sequence evolution that best fitted our data. The four cpDNA fragments were tested both separately and concatenated into a single sequence. As F81 and HKY were consistently identified as best models, we used the model HKY available in BEAST. A Coalescent Constant Population was employed as tree prior and we implemented the strict molecular clock recommended for intra-species analysis (no expected rate variation among branches). Given the wide variation of substitution rates seen among flowering plants (; ), and considering that the mutation rate in the literature for the cpDNA fragments used in our study was 1.0−9 (), we intentionally specified a broad prior using an uniform distribution with bounds 1.0−10 and 1.0−8 (substitutions site−1 year−1), which allowed BEAST to obtain the precise mutation rate. The four cpDNA fragments were designed as unlinked partitions to allow for variable mutation rates. Chain convergence was assessed with Tracer 1.6 (available in the BEAST2 package at ) and by checking the effective sampling size (ESS) values. Trees were summarized in a maximum clade credibility tree obtained in TreeAnotator 2.1 and visualized in FigTree 1.4.1 (both also available in the BEAST2 package).The calculation of the TMRCAs for the groups later detected was based on independent running of BEAST, by combining two searches with 10 million MCMC, using only the haplotypes assigned to each group.Climatic reconstructionThe climatic reconstruction observed in Fig.  took as reference a couple of maps created thanks to the ‘stage 3 project’ (). The first was the map ‘tmin_21k’, which represents the seasonal minimum air temperature 21 000 years ago. The second was the map ‘tmin_mod’, which describes the current seasonal minimum air temperature. Both maps are available at .Areas where the minimum monthly temperatures suitable for A. unedo (>4 °C) occurred: during the LGM (21 ky, blue) and in modern times (green).Genetic diversity and population structureWe produced sequences for 207 individuals of A. unedo (Table ) plus the two individuals of A. andrachne. The four non-coding fragments revealed polymorphisms along the range of A. unedo and were retained for phylogenetic analyses. Once concatenated, and depending on the occurrence of indels, the four fragments resulted in sequences 2401&#x bp long, equivalent to 1.6 % of the chloroplast genome of A. unedo (). The data set of A. unedo contained 16 polymorphic sites that included 7 point mutations and 9 indels. Three indels were 1 bp long, two indels were 9 bp long and the remaining four indels were 10, 11, 35 and 36 bp long, respectively.A total of 28 haplotypes were found [see
], and half of them (14) were private (Hp); interestingly, five Hp were detected in North Africa (Morocco). The two individuals of A. andrachne produced a single haplotype that differed by 16 new point mutations from the consensus sequence of A. unedo. Haplotype diversity in A. unedo was 0.901 (±0.01) (standard error), while nucleotide diversity (π) (×102) was 1.238 (±0.60). Nonetheless, genetic diversity was unevenly partitioned among populations as Hd and π estimates for the 23 populations ranged 0.000&#x and 0.000&#x, respectively (Table ). The highest π values were consistently recorded in populations that contained a mixture of haplotypes separated by long indels (35 or 36 bp).Total gene diversity (ht) was 0.919 (±0.02), while the average within-population gene diversity (hs) was 0.585 (±0.04). We also calculated the diversity at a regional scale by partitioning the populations into three groups suggested by the phylogeographic analyses (see below): Eastern Mediterranean (GR1, GR2, TU1 and TU2), Western Mediterranean (FR1, FR2, SP7, SP8, IT1-IT3 and TUN) and Atlantic coast (IRL, SP1-SP6, PO1, PO2, MO1 and MO2). Diversity at regional scale increased from the Eastern (Ht = 7, Hd = 0.693 ± 0.076, Π × 102 = 0.58 ± 0.30) to the Western Mediterranean (Ht = 15, Hd = 0.762 ± 0.044, Π × 102 = 0.63 ± 0.32), and to the Atlantic coast (Ht = 13, Hd = 0.787 ± 0.029, Π × 102 = 0.75 ± 0.38).The AMOVA revealed that 38.8 % of the genetic variation was due to differences between populations, and the resulting FST value was highly significant (P < 0.001). NST (0.645 ± 0.067) was significantly higher than GST (0.363 ± 0.040) (P < 0.001) indicating the existence of a strong phylogeographic structure (). Furthermore, SAMOVA showed that FCT began to plateau at 0.697 for three groups of populations regardless of the alternative used (with or without constraint for the geographic composition of the groups) (Table ). One group was geographically homogeneous and included the populations sampled in the Eastern Mediterranean plus Bordeaux (FR1) in the Atlantic France. Another group was split into two sets of non-adjacent populations because populations from Atlantic Iberia–North Africa grouped with populations from Western Mediterranean. Finally, a third, smaller group clustered the Irish population with two locations in Atlantic Iberia. Less than 8 % of the variation was among populations within groups and nearly 23 % was within populations.
Results of the SAMOVA for the number of groups (K = 3) for which FCT reached a plateau. ***P < 0.001.Haplotype genealogy and distributionHaplotypes H2 and H4 dominated the data set (37 %), whereas 12 haplotypes were very rare, as they were detected in just one (10 unique haplotypes) or two individuals (2 haplotypes) [see
]. The median-joining network (Fig. ) yielded a genealogy where many haplotypes were just one mutational step away from their closest relative. The topology of the network also revealed two clades separated from the ancestral node inferred with the help of A. andrachne by two mutational steps, and from each other by four mutational steps. One clade was dominated by H2 at the centre of 17 haplotypes arranged in a star-like pattern (henceforth H2-derived), while the other clade was dominated by H4, again surrounded by 9 less-common haplotypes (H4-derived).
Median-joining haplotype network for 18 A. unedo non-coding cpDNA haplotypes with A. andrachne as external rooting. Each circle represents a haplotype (name outside) while colours indicate sampling regions (SW Iberia–Morocco = PO1, PO2, SP6, MO1 and ...The partition of haplotypes in two large lineages was remarkably consistent with their geographical location as well as with the clusters of populations identified by SAMOVA (Fig. ). Accordingly, the dominant haplotypes H2 and H4 were regionally segregated and their ranges barely overlapped because H4 was detected only in the Mediterranean, while plants assigned to H2 came from sites along the Atlantic and North Africa. Kroumerie (TUN) in North Africa was the only site where H2 and H4 were found living in sympatry. Like H2, most of the H2-derived haplotypes occurred along the Atlantic from Ireland to Morocco. In fact, H2 and H2-derived sequences (henceforth the Atlantic clade) were the only haplotypes detected in populations from the Atlantic shoreline other than FR1 (see below). Interestingly, however, another subset of five closely related H2-derived sequences had a non-Atlantic provenance as they belonged to plants sampled at Sardinia (IT3), Rome (IT2) and, particularly, the Eastern Mediterranean (Greece and Turkey). In comparison, no haplotype from the Atlantic clade was ever found in the NW Mediterranean (NE Spain, W France and N Italy). The H4-derived haplotypes (Mediterranean clade) also resembled the geographical distribution seen in their ancestral haplotype H4 as they were largely confined to the Mediterranean basin. The exceptions were two H4-derived haplotypes (H6 and H12) that eventually reached a small section of the Atlantic species' range in W France (FR1), where they lived with congeners from the Atlantic phylogroup. Interestingly, the ancestral haplotypes H2 and H4 seemed largely absent from the populations sampled at the northern (Ireland) and eastern (Greece and Turkey) edges of the species' only H4 was detected in three individuals from Turkey (population TU2). Additionally, we analysed the trees obtained with TCS [see
] and BEAST [see
]. The previous description in two sublineages was corroborated and also the distribution of the main haplotypes. Even if the general patterns were concordant, some exceptions were found. TCS assigned H11 and H10 to the sublineage H2, while H20 was placed into sublineage H4 and it placed H7 and H9 as a different group, although closest to H2 than to H4 (that was checked with complementary analysis) (not shown). Likewise, BEAST analysis located H15 into sublineage H2, while H18 and H20 were placed into sublineage H4. Despite these minor modifications, the previous interpretation of the Network tree is still consistent.Timing of the diversificationInitial trials in BEAST were run using exponential growth as the coalescent tree prior. Regardless of the data set tested (complete data set, Atlantic populations, Mediterranean populations, individuals from the Atlantic clade or individuals from the Mediterranean clade), the posterior distribution of the growth rate [95 % highest posterior density (HPD)] always included zero indicating no support for the existence of exponential growth and establishes the plausibility of the constant population growth model, based on our data. The absence of demographic expansion was confirmed by Fu's FS tests run for the complete data set and for each genetic clade treated separately (P-value always &#x). Therefore, the final estimates of TMRCA relied on a constant growth model. However, the presence of 11 out of the 14 private haplotypes in the Atlantic clade might indicate a higher ability to diverge and spread. The latter is further supported by the higher number of haplotypes and more extensive geographical distribution of this clade. The Bayesian tree inferred by BEAST reproduced the neat division in two clades (Atlantic and Mediterranean) found in the network analysis [see
]. The BEAST inference suggested that the split between Atlantic and Mediterranean lineages must have occurred before the LGM, around the Middle/Late Pleistocene boundary. Thus, the TMRCA for A. unedo was 365 ky (95 % HPD: 48&#x effective sample size, ESS = 1295). The calculation of the TMRCAs for the Atlantic and the Mediterranean lineages was 258 ky (95 % HPD: 43&#x ky, ESS = 804) for the Atlantic and 194 ky (95 % HPD: 9&#x ky, ESS = 596) for the Mediterranean.East–West patternThe most conspicuous attribute in our data set is the separation of the strawberry tree into two neat matrilineal lineages with longitude. East–west phylogeographic breaks seem common in the Mediterranean, and extreme cases include plants with disjunct populations on both ends of the Mediterranean (; ; ). In plants, however, most examples involve species with continuous ranges that still show a phylogeographic split between the Eastern and the Western Mediterranean (; ; ; ). The peculiarity in the strawberry tree is that both lineages appear to have originated in the western end of Mediterranean region: one in the Atlantic coasts of Iberia and North Africa, and the other in the Mediterranean basin itself. The presence of several genetic clades in the Western Mediterranean has been observed in other plants and was linked to range contractions during the Quaternary glaciations (; ). In the case of A. unedo, we infer that the diversification in two lineages also occurred in the late Pleistocene but before the LGM, suggesting that it may have coincided with the hardest glaciations recorded in the Quaternary (; ). Moreover, the strawberry tree possibly survived the late Quaternary in the Western Mediterranean given the observation that the ancestral haplotypes H2 and H4 are largely restricted to this region. In this regard, several areas that had previously been attributed a lesser role are now considered to have played an important role as glacial refugia (). Interestingly, the latter include several areas located in the Western Mediterranean (large Mediterranean islands, North Africa and Catalonia) that may also have played a role in the case of A. unedo.The more recent colonization of the Eastern Mediterranean by the strawberry tree contrasts with most studies of plants where eastern lineages typically are more ancient (; ; ; ; ). In trees, the predominance of westward migrations has been attributed to divergence that predates the Pleistocene (, ; ). Alternatively, however, genetic and fossil evidence linked the East–West break to range fragmentations and to more severe climate in the west during the Pleistocene (; ). In comparison, there are few examples of colonization in the opposite direction. Still, North African populations were likely candidates for the ancestral pool of Europhaca (), while the submediterranean alpine Anthyllis montana also migrated eastward from Iberia up to the Balkans along the northern edge of the Mediterranean (). Similarly, Pinus pinaster () and Erica arborea () seem to have migrated eastward along North Africa to Europe. Nonetheless, the closer example to the case of the strawberry tree is Myrtus communis, another circum-Mediterranean shrub with fruits dispersed by birds. Following differentiation events during the Pleistocene, some matrilineal lineages of M. communis spread from west to east, in particular, to the Balkan Peninsula where they met with older lineages (). In A. unedo, an eastward colonization would seem consistent with (i) a predominance of younger haplotypes in this area (H2-derived and H4-derived) and (ii) the likely unsuitability of the West Mediterranean during glacial maxima. Accordingly, paleoclimatic reconstructions of the LGM indicate that minimum monthly temperatures in the Aegean–Anatolia area were below the 4 °C limit required by the strawberry tree for survival () (Fig. ).
The role of the straitsThe likely route for the movement of the Atlantic lineage towards the Eastern Mediterranean was through North Africa to mainland Italy and the Balkans by crossing the Strait of Sicily. The latter has been involved in biogeographical connections for many taxa. In trees, a crossing by matrilineal lineages from North Africa has been inferred for some species (; ; ) but not for others (; ). Nonetheless, the strait still had an impact on the migration of the Atlantic clade of A. unedo because only a handful of H2-derived haplotypes managed to reach the Eastern Mediterranean, contrarily to other trees where several haplotypes were common to both sides of the strait (; ). Our estimate of the TMRCA for the Atlantic clade suggests that the strawberry tree crossed the Strait of Sicily in the last 200 000 years but still before the LGM. Interestingly, the impact of the strait on the migration of the Mediterranean lineage was different because the same set of haplotypes occurs on both sides of the strait. However, the Mediterranean clade seemed unable to spread beyond Tunisia, a confinement that would be consistent with the ‘high density blocking hypothesis’ () if North Africa was already colonized by the Atlantic lineage before the arrival of the Mediterranean clade.Gibraltar, the other large strait in the region, also had partial success in interrupting the gene flow of A. unedo. Only plants with Atlantic ancestry live on both sides of the strait, and the ancestral haplotype H2 is equally common to the north and to the south of Gibraltar. However, none of the H2-derived haplotypes recorded in Iberia occurs in North Africa and the other way around. Hence, the isolation of the populations on both sides of the strait must have lasted long enough to allow the diversification of distinct sets of haplotypes on each area, such as many of the private haplotypes, supporting the role of North Africa as a refugium for A. unedo. Interestingly,
noted that colonization ability, rather than dispersal potential, is a better predictor of the genetic structure across the strait.Genetic diversity patterns and refugia within refugiaThe phylogeographic arrangement of A. unedo does not show the west to east increase in genetic diversity seen in many other plants and animals (). Instead, genetic diversity within populations is notably variable and populations with high or low diversity are interspersed all along the species' range. If anything, genetic diversity increases from East to West at a regional level due to the dominance of plants with an Atlantic ancestry in the West and the small number of haplotypes detected in this region.The predominance of H4-derived haplotypes in FR1 reveals that the Mediterranean lineage managed to reach the Atlantic coast. Interestingly, FR1 shares a private haplotype (H6) with north-east Iberia (SP7) and the Balearics (SP8). Indeed, north-east Iberia (SP7) is more closely related to French demes than to other Iberian populations. Thus, the strawberry tree adds to the many cases of trees with a clear-cut East–West division in the Iberian Peninsula (), providing another example of the pattern typically attributed to the ‘refugia-within-refugia’ model (). Since western and eastern Iberia are characterized by different soil types, it has been speculated that soil type might have promoted the isolation of lineages in the peninsula (). However, the potential of ecological factors to determine the genetic structure of A. unedo is uncertain since common garden experiments show the absence of local adaptation for morphological and physiological traits in this species ().Finally, our results discard a Tertiary origin for A. unedo in Ireland (). Instead, a relatively recent colonization of Ireland seems a more plausible alternative given (i) the absence of private haplotypes in Ireland and (ii) the genetic proximity between Irish and Iberian population (SAMOVA results). The arrival of the strawberry tree to Ireland has been tentatively attributed to movements along northern Iberia and the maritime fringe of France (; ). Our results do not necessarily support this hypothesis because none of the haplotypes recorded in Ireland was detected in FR1, at the northern edge of the continuous range of A. unedo in the Atlantic. As noted above, FR1 is dominated by haplotypes with a Mediterranean ancestry that are conspicuously absent from Ireland. Moreover, the few H2-derived haplotypes found in FR1 are different from the only H2-derived haplotype detected in Ireland. Alternatively, France could still appear as a plausible source for the Irish populations if the set of haplotypes that currently inhabit FR1 arrived there after the colonization of Ireland. However, this alternative explanation would require a total displacement/disappearance of any previous haplotype in FR1, and we cannot find any likely explanation for such replacement. Hence, our data seem more easily consistent with an Iberian origin for the Irish populations of the strawberry tree. A similar Iberian origin has been proposed for oaks (; ), and the arrival to Ireland by crossing the Celtic Sea has been suggested for other tree species ().The view that species respond on an individual basis to climate change and create new community patterns has gained increasing support in the last decades (). Here, we have shown concordances and discrepancies between the phylogeographic structure of the strawberry tree and that of other plants investigated in the Mediterranean basin. The partition of the strawberry tree into two clades could be attributed to vicariance events during the Quaternary glaciations of the last 700 ky, but before the LGM. These vicariance events possibly occurred in glacial refugia located in the western end of the Mediterranean region. Thus, the Atlantic clade possibly survived in North Africa–Atlantic Iberia, whereas the Mediterranean clade may have resisted in the large islands and the European rim of the Western Mediterranean sub-basin. We also found evidence of occasional connections between Europe and North Africa through the straits of Gibraltar and Sicily, but North Africa still retains a distinctive genetic composition that provides further support to its role as a glacial refugium. From these refugia, the Eastern Mediterranean possibly was colonized more recently, and its lower genetic diversity contrasts with the pattern commonly seen in most plants. The East–West genetic split shown by the strawberry tree in the Iberian Peninsula resembles the pattern found in many other trees and provides a new example of genetic complexity that seems consistent with the refugia-within-refugia model (). Our results discard a Tertiary origin for A. unedo in Ireland and do not support a continuous range along the maritime fringe of France either. Instead, the strawberry tree may have arrived to Ireland in post-glacial times from N Iberia by crossing the Celtic Sea. Altogether, our results reveal a considerable ability to disperse for the strawberry tree, migrating over thousands of kilometres and even crossing stretches of sea (Mediterranean straits). This ability for long-distance migration could be useful in a scenario of future changes, allowing the species to migrate to new areas such as Northern France, where a drier and warmer climate, with a greater number of extreme weather events, it is expected (; ).These sequences can be obtained from GenBank using the accession numbers: BankIt 1868954: KU 194976–KU195222 (Intron L fragment); BankIt 1874447: KU205359–KU205584 (trnT-tnrL fragment); BankIt 1874452: KU205585–KU205820 (trnL-rpl32 fragment); BankIt 1874463: KU205821–KU206025 (trnH-psbA fragment).This work was supported by the Spanish Ministry of Science and Innovation (research grant CGL), the European Regional Development's Fund (ERDF) and also by the Spanish Ministry of Education (FPU fellowship AP- to X.S.).R.R. and R.B. X.S. and R.R. X.S. and L.L. performed th X.S., L.L. and R.B. analysed the data and X.S., L.L., R.B. and R.R. led the writing.None declared.The following additional information is available in the online version of this article –Figure S1. The variable positions in the 28 haplotypes detected in A. unedo. The first 18 haplotypes belong to the Atlantic clade and the last 10 to the Mediterranean one.Figure S2. The genealogy of the haplotypes obtained with FigTree, which visualized the analysis made with BEAST and TreeAnotator. We observed two main clades, the upper one (Atlantic) with 17 haplotypes and the lower one (Mediterranean) with 11 haplotypes.Figure S3. The genealogy of the haplotypes obtained with TCS. They are grouped into three clades, being the ancestral haplotype of each clade represented inside the rectangle.Table S1. The number of individuals per haplotype in 23 populations of A. unedo.
Additional Information:
We thank Dr Fatima Ain-Lhout and her father who generously sampled and supplied plant material. We also thank Dr R. BermÞz, S. Louro, Chris Barron, Dr X. Sánchez, Dr S. Roiloa and Dr Murat Seker for their helpful contribution in sampling.Avise JC.
What is the field of biogeography, and where is it going?
53:893&#x. Avise JC.
Phylogeography: retrospect and prospect. Journal of Biogeography
36:3–15. Avise JC, Ball RM
Principles of genealogical concordance in species concepts and biological taxonomy. Oxford Surveys in Evolutionary Biology
7:45–67.Axelrod DI.
Evolution and biogeography of madrean-tethyan sclerophyll vegetation. Annals of the Missouri Botanical Garden
62:280&#x. Bandelt HJ, Forster P, Röhl A
Median-joining networks for inferring intraspecific phylogenies. Molecular Biology and Evolution
16:37–48.
[]Beatty GE, Provan J
Post-glacial dispersal, rather than in situ glacial survival, best explains the disjunct distribution of the Lusitanian plant species Daboecia cantabrica (Ericaceae). Journal of Biogeography
40:335&#x. Beatty GE, Provan J
Phylogeographical analysis of two cold-tolerant plants with disjunct Lusitanian distributions does not support in situ survival during the last glaciation. Journal of Biogeography
41:2185&#x. Bouckaert R, Heled J, K�hnert D, Vaughan T, Wu CH, Xie D, Suchard MA, Rambaut A, Drummond AJ
BEAST 2: a software platform for Bayesian evolutionary analysis. PLoS Computational Biology
10:e1003537
[]Bucci G, González-Mart�nez SC, Le Provost G, Plomion C, Ribeiro MM, Sebastiani F, Al R, Vendramin GG
Range-wide phylogeography and gene zones in Pinus pinaster Ait. revealed by chloroplast microsatellite markers. Molecular Ecology
16:2137&#x.
[]Casimiro-Soriguer R, Talavera M, Balao F, Terrab A, Herrera J, Talavera S
Phylogeny and genetic structure of Erophaca (Leguminosae), a East–West Mediterranean disjunct genus from the Tertiary. Molecular Phylogenetics and Evolution
56:441&#x.
[]Clausing G, Vickers K, Kadereit JW
Historical biogeography in a linear system: genetic variation of Sea Rocket (Cakile maritima) and Sea Holly (Eryngium maritimum) along European coasts. Molecular Ecology
9:1823&#x.
[]Clement M, Posada D, Crandall KA
TCS: a computer program to estimate gene genealogies. Molecular Ecology
9:1657&#x.
[]Conord C, Gurevitch J, Fady B
Large-scale longitudinal gradients of genetic diversity: a meta-analysis across six phyla in the Mediterranean basin. Ecology and Evolution
2:2600&#x.
[]Cox C, Moore PD
Biogeography: an ecological and evolutionary approach, 7th edn
Oxford, UK: Wiley-Blackwell.Darriba D, Taboada GL, Doallo R, Posada D
jModelTest 2: more models, new heuristics and parallel computing. Nature Methods
[]Davis B, Brewer S, Stevenson A, Guiot J
The temperature of Europe during the Holocene reconstructed from pollen data. Quaternary Science Reviews
22:1701&#x. Désamoré A, Laenen B, Devos N, Popp M, González-Mancebo JM, Carine MA, Vanderpoorten A
Out of Africa: north-westwards Pleistocene expansions of the heather Erica arborea. Journal of Biogeography
38:164&#x. De Vries H, Haarsma RJ, Hazeleger W
Western European cold spells in current and future climate. Geophysical Research Letters
Dupanloup I, Schneider S, Excoffier L
A simulated annealing approach to define the genetic structure of populations. Molecular Ecology
11:2571&#x.
[]Escudero M, Vargas P, Arens P, Ouborg NJ, Luceño M
The east-west-north colonization history of the Mediterranean and Europe by the coastal plant Carex extensa (Cyperaceae). Molecular Ecology
19:352&#x.
[]Excoffier L, Lischer HEL
Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Molecular Ecology Resources
10:564&#x.
[]Fazekas AJ, Steeves R, Newmaster SG
Improving sequencing quality from PCR products containing long mononucleotide repeats. Biotechniques
48:277&#x.
[]Fernández-Mazuecos M, Vargas P
Ecological rather than geographical isolation dominates Quaternary formation of Mediterranean Cistus species. Molecular Ecology
19:1381&#x.
[]Foss P, Doyle G, Nelson E
The distribution of Erica erigena R. Ross in Ireland. Watsonia
16:311&#x.Fu Y-X.
Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics
147:915&#x.
[]Gómez A, Lunt DH
Refugia within refugia: patterns of phylogeographic concordance in the Iberian Peninsula. In: Weiss S, Ferrand N, eds. Phylogeography of southern European refugia. The Netherlands: Springer.Habets F, Boé J, Déqué M, Ducharne A, Gascoin S, Hachour A, Martin E, Pagé C, Sauquet E, Terray L, Thiéry D, Oudin L, Viennot P
Impact of climate change on the hydrogeology of two basins in northern France. Climatic Change
121:771&#x. Hampe A, Arroyo J, Jordano P, Petit RJ
Rangewide phylogeography of a bird-dispersed Eurasian shrub: contrasting Mediterranean and temperate glacial refugia. Molecular Ecology
12:3415&#x.
[]Hewitt G.
Post-glacial re-colonization of European biota. Biological Journal of the Linnean Society
68:87&#x. Hewitt G.
The genetic legacy of the Quaternary ice ages. Nature
405:907&#x.
[]Hewitt G.
Mediterranean peninsulas: the evolution of hotspots. In: Zachos FE, Habel JC, eds. Biodiversity hotspots. Berlin: Springer.Hileman LC, Vasey MC, Parker VT
Phylogeny and biogeography of the Arbutoideae (Ericaceae): implications for the Madrean-Tethyan hypothesis. Systematic Botany
26:131&#x.Huang C-C, Hung K-H, Wang WK, Ho CW, Huang CL, Hsu TW, Osada N, Hwang CC, Chiang TY
Evolutionary rates of commonly used nuclear and organelle markers of Arabidopsis relatives (Brassicaceae). Gene
499:194&#x.
[]Kadereit G, Yaprak AE
Microcnemum coralloides (Chenopodiaceae-Salicornioideae): an example of intraspecific East-West disjunctions in the Mediterranean region. Anales del Jard�n Botánico de Madrid
65:415&#x.Kadereit JW, Arafeh R, Somogyi G, Westberg E
Terrestrial growth and marine dispersal? Comparative phylogeography of five coastal plant species at a European scale. Taxon
54:861&#x. Kelleher C, Hodkinson T, Kelly D, Douglas G
Characterisation of chloroplast DNA haplotypes to reveal the provenance and genetic structure of oaks in Ireland. Forest Ecology and Management
189:123&#x. Kropf M, Kadereit JW, Comes HP
Late Quaternary distributional stasis in the submediterranean mountain plant Anthyllis montana L. (Fabaceae) inferred from ITS sequences and amplified fragment length polymorphism markers. Molecular Ecology
11:447&#x.
[]Librado P, Rozas J
DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics
25:1451&#x.
[]López de Heredia U, Carrión JS, Jiménez P, Collada C, Gil L
Molecular and palaeoecological evidence for multiple glacial refugia for evergreen oaks on the Iberian Peninsula. Journal of Biogeography
34:1505&#x. Lowe A, Unsworth C, Gerber S, Davies S, Munro R, Kelleher C, King A, Brewer S, White A, Cottrell J
Route, speed and mode of oak postglacial colonisation across the British Isles: integrating molecular ecology, palaeoecology and modelling approaches. Botanical Journal of Scotland
57:59–81. Lumaret R, Mir C, Michaud H, Raynal V
Phylogeographical variation of chloroplast DNA in holm oak (Quercus ilex L.). Molecular Ecology
11:2327&#x.
[]Lumaret R, Tryphon-Dionnet M, Michaud H, Sanuy A, Ipotesi E, Born C, Mir C
Phylogeographical variation of chloroplast DNA in cork oak (Quercus suber). Annals of Botany
96:853&#x.
[]Mart�nez-Alberola F, Del Campo EM, Lázaro-Gimeno D, Mezquita-Claramonte S, Molins A, Mateu-Andrés I, Pedrola-Monfort J, Casano LM, Barreno E
Balanced gene losses, duplications and intensive rearrangements led to an unusual regularly sized genome in Arbutus unedo chloroplasts. PLoS ONE
[]Médail F, Diadema K
Glacial refugia influence plant diversity patterns in the Mediterranean Basin. Journal of Biogeography
36:1333&#x. Migliore J, Baumel A, Juin M, Médail F
From Mediterranean shores to central Saharan mountains: key phylogeographical insights from the genus Myrtus. Journal of Biogeography
39:942&#x. Mitchell FJG.
The biogeographical implications of the distribution and history of the strawberry tree, Arbutus unedo in Ireland. In: Costello MJ, Kelly KS, eds. Biogeography of Ireland: past, present, and future. Dublin: Irish Biogeographical Society.Mitchell FJG.
Where did Ireland's trees come from?
Biology & Environment: Proceedings of the Royal Irish Academy
106:251&#x. Myers N, Mittermeier RA, Mittermeier CG, Da Fonseca GAB, Kent J
Biodiversity hotspots for conservation priorities. Nature
403:853&#x.
Molecular evolutionarygenetics. New York: Columbia University Press.Nieto Feliner G.
Patterns and processes in plant phylogeography in the Mediterranean Basin. A review. Perspectives in Plant Ecology, Evolution and Systematics
16:265&#x. Palamarev E.
Paleobotanical evidences of the Tertiary history and origin of the Mediterranean sclerophyll dendroflora. Plant Systematics and Evolution
162:93&#x. Petit RJ, Hampe A, Cheddadi R
Climate changes and tree phylogeography in the Mediterranean. Taxon
54:877&#x. Peyron O, Magny M, Goring S, Joannin S, Beaulieu JD
Pollen-inferred quantitative reconstruction of the Holocene climate in the central Mediterranean area (Italy). Geophysical Research Abstracts
14:7377.Polzin T, Daneshmand SV
On Steiner trees and minimum spanning trees in hypergraphs. Operations Research Letters
31:12–20. Pons O, Petit RJ
Measuring and testing genetic differentiation with ordered versus unordered alleles. Genetics
144:1237&#x.
[]Quézel P, Médail F
Ecologie et biogéographie des for�ts du bassin méditerranéen. Paris: Elsevier.Rodr�guez-Sánchez F, Pérez-Barrales R, Ojeda F, Vargas P, Arroyo J
The Strait of Gibraltar as a melting pot for plant biodiversity. Quaternary Science Reviews
27:2100&#x. Rodr�guez-Sánchez F, Guzmán B, Valido A, Vargas P, Arroyo J
Late Neogene history of the laurel tree (Laurus L., Lauraceae) based on phylogeographical analyses of Mediterranean and Macaronesian populations. Journal of Biogeography
36:1270&#x. Rodr�guez-Sánchez F, Hampe A, Jordano P, Arroyo J
Past tree range dynamics in the Iberian Peninsula inferred through phylogeography and palaeodistribution modelling: a review. Review of Palaeobotany and Palynology
162:507&#x. Santiso X, López L, Gilbert KJ, Barreiro R, Whitlock MC, Retuerto R
Patterns of genetic variation within and among populations in Arbutus unedo and its relation with selection and evolvability. Perspectives in Plant Ecology, Evolution and Systematics
17:185&#x. Sealy JR.
Arbutus unedo. The Journal of Ecology
37:365&#x. Shaw J, Lickey EB, Beck JT, Farmer SB, Liu W, Miller J, Siripun KC, Winder CT, Schilling EE, Small RL
The tortoise and the hare II: relative utility of 21 noncoding chloroplast DNA sequences for phylogenetic analysis. American Journal of Botany
92:142&#x.
[]Shaw J, Lickey EB, Schilling EE, Small RL
Comparison of whole chloroplast genome sequences to choose noncoding regions for phylogenetic studies in angiosperms: the tortoise and the hare III. American Journal of Botany
94:275&#x.
[]Smith RJ, Waldren S
Patterns of genetic variation in Colchicum autumnale L. and its conservation status in Ireland: a broader perspective on local plant conservation. Conservation Genetics
11:1351&#x. Smith SA, Donoghue MJ
Rates of molecular evolution are linked to life history in flowering plants. Science
322:86–89.
[]Stewart JR, Lister AM, Barnes I, Dalén L
Refugia revisited: individualistic responses of species in space and time. Proceedings of the Royal Society B: Biological Sciences
277:661&#x.
[]Taberlet P, Gielly L, Pautou G, Bouvet J
Universal primers for amplification of three non-coding regions of chloroplast DNA. Plant Molecular Biology
17:1105&#x.
[]Torres JA, Valle F, Pinto C, Garc -Fuentes A, Salazar C, Cano E
Arbutus unedo L. communities in southern Iberian Peninsula mountains. Plant Ecology
160:207&#x. Van Andel TH.
The climate and landscape of the middle part of the Weichselian glaciation in Europe: the stage 3 project. Quaternary Research
57:2–8. Van Rijn P.
The analysis of charcoal from Ross Island. In: O'brien W, editor. , ed. Ross Island. Mining, metal and society in early Ireland. Bronze age studies 6
Galway: National University of Ireland.Webb DA.
The flora of Ireland in its European context. Journal of Life Sciences, Royal Dublin Society
4:143&#x.Articles from AoB Plants are provided here courtesy of Oxford University Press&