Evolutionary and Biogeographic Insights on the Macaronesian Beta-Patellifolia Species (Amaranthaceae) from a Time-Scaled Molecular Phylogeny

The Western Mediterranean Region and Macaronesian Islands are one of the top biodiversity hotspots of Europe, containing a significant native genetic diversity of global value among the Crop Wild Relatives (CWR). Sugar beet is the primary crop of the genus Beta (subfamily Betoideae, Amaranthaceae) and despite the great economic importance of this genus, and of the close relative Patellifolia species, a reconstruction of their evolutionary history is still lacking. We analyzed nrDNA (ITS) and cpDNA gene (matK, trnH-psbA, trnL intron, rbcL) sequences to: (i) investigate the phylogenetic relationships within the Betoideae subfamily, and (ii) elucidate the historical biogeography of wild beet species in the Western Mediterranean Region, including the Macaronesian Islands. The results support the Betoideae as a monophyletic group (excluding the Acroglochin genus) and provide a detailed inference of relationships within this subfamily, revealing: (i) a deep genetic differentiation between Beta and Patellifolia species, which may have occurred in Late Oligocene; and (ii) the occurrence of a West-East genetic divergence within Beta, indicating that the Mediterranean species probably differentiated by the end of the Miocene. This was interpreted as a signature of species radiation induced by dramatic habitat changes during the Messinian Salinity Crisis (MSC, 5.96–5.33 Mya). Moreover, colonization events during the Pleistocene also played a role in shaping the current diversity patterns among and within the Macaronesian Islands. The origin and number of these events could not be revealed due to insufficient phylogenetic resolution, suggesting that the diversification was quite recent in these archipelagos, and unravelling potential complex biogeographic patterns with hybridization and gene flow playing an important role. Finally, three evolutionary lineages were identified corresponding to major gene pools of sugar beet wild relatives, which provide useful information for establishing in situ and ex situ conservation priorities in the hotspot area of the Macaronesian Islands.


Introduction
. Native geographical distribution, ecology, and IUCN conservation status of taxa from subfamily Betoideae; taxonomy according to Kadereit et al. [16]. Patellifolia. The latter differs from Beta by having short tepals that do not overtop the fruit vs. long tepals that overtop the fruit [12]. Previous studies based on morphological features (e.g. [10,[13][14]) failed to recognize Patellifolia as a separate genus but rather as part of the Beta section Procumbentes. Recent molecular phylogenetic studies (e.g. [15,16]) modified the subfamily classification previously proposed. It was also suggested that Acroglochin should be excluded from this subfamily and that the other five genera (i.e. Beta, Aphanisma, Hablitzia, Oreobliton, and Patellifolia) should fall into two clades, i.e. Beteae comprising Beta only, and Hablitzieae with the remaining four genera. These studies have been hampered by the undersampling of species from the Western Mediterranean Region, including the hotspot area of the Macaronesian Islands, where some endemic species are found (i.e. B. patula in the Madeira archipelago, P. webbiana in the Canary Islands, and P. procumbens in all the Macaronesian archipelagos except the Azores). Two of these Macaronesian endemics (i.e. B. patula and P. webbiana) were recently classified as Critically Endangered in the European Red List of Vascular Plants [17]. Though the importance of conservation of these wild taxa has been widely recognized [18], it is also important to understand the relationships within the Beta s.l. gene pools, which will offer an effective approach for utilization of the wild-beet germplasm. For instance it is pointed out that Patellifolia species can transmit traits providing resistance to the most serious diseases of sugar beets worldwide, such as sugar beet cyst nematode (Heterodera schachtii Schmidt), leaf spot disease caused by Cercospora beticola Sacc., curly top virus, rhizomania, and powdery mildew (Erysiphe polygoni DC.) [19][20][21][22].
Despite the great economic importance of the Beta and Patellifolia species [18], a reconstruction of the evolutionary history with a dated molecular phylogeny for the subfamily Betoideae is still lacking. The aims of this study are to: (1) present a hypothesis of the phylogenetic relationships within the subfamily Betoideae, and (2) gain a better understanding of the spatiotemporal history of the wild beet species which occur in the hotspot area of Western Mediterranean Region, including for the first time the endemic species from the Macaronesian Islands and samples from all the five archipelagos (i.e. the Azores, Canaries, Cape Verde, Madeira including the Desertas, and Savage Islands). Taxa marked with an * were included in the phylogenetic analyses; those marked with ** were collected for this study. a Data from fieldwork and bibliography [23][24][25][26][27][28][29][30]. Permissions to collect protected species from protected areas were issued by Portuguese authorities  GenBank accession numbers are provided in S1 Table for all the studied samples. Additionally, data on sampling sites of the samples collected in this study, including their geographical coordinates, details about vouchers and their respective herbaria, are also included in S1 Table. Molecular data DNA was extracted using DNeasy Plant Mini Kit (QIAGEN, Valencia, California, USA) and purified, using QIAquick columns (QIAGEN, Valencia, California, USA) or Silica Bead DNA Gel Extraction kit (Fermentas), according to the manufacturer's protocols. Polymerase chain reaction (PCR) amplifications using 20-30ηg of genomic DNA were performed to amplify the complete internal transcribed spacer (ITS) region, using the primers ITS4 and ITS5 [31]. Two coding regions of the chloroplast genome were amplified using primer pairs, matK [32] and rbcL [33], plus two non-coding regions using trnL intron [34] and trnH-psbA [35].
Sixty-eight samples were sequenced for the ITS region (S1 Table). From a preliminary study at population level of the four chloroplast regions, all individuals of each species, of Beta and Patellifolia, and those from a given sampling location exhibited the same cpDNA sequences. Therefore, cpDNA sequences were generated by selecting a subset of 1-5 individuals/island, resulting in 27 representative accessions sequenced. This sub-sample was generated for matK, trnH-psbA, trnL intron and rbcL (S2 Table), but three specimens for the rbcL gene and two specimens for the trnL intron and trnH-psbA spacer could not be sequenced due to PCR amplification problems. Amplified products were purified with Sureclean Plus (Bioline, London, UK) and sent to STAB Vida, Lda (Monte da Caparica, Portugal) for Sanger sequencing. For all the markers, amplicons were sequenced using both directions in ABI 3730 XL DNA Analyzer (Applied Biosystems). Raw sequences were edited and cleaned by hand in SEQUENCHER v4.0.5 (Gene Codes Corporation).

Phylogenetic analyses
Multiple sequence alignments were built for each locus dataset in MAFFT v6.717b [36], using the L-INS-i method as recommended in the manual for difficult alignments. Each dataset was concatenated into a combined matrix using ElConcatenero [37]. Maximum Likelihood (ML) and Bayesian Inference (BI) methods were used to reconstruct phylogenies from the separate (i.e. ITS, matK, trnH-psbA, trnL intron, and rbcL) and combined datasets. ML searches were performed in RAxML v8.0.9 under the GTRGAMMA model with 1000 bootstrap replicates. The best fit model for each locus data set was selected under the AIC, as implemented in MRMODELTEST v2.3 [38] and used in the Bayesian analysis performed in MRBAYES v3.1.2 [39]. Each locus was allowed to have partition-specific substitution parameters. Analyses were generated for 3x10 7 generations, sampled every 3000 th generation and using the default chain heating temperature. The analysis was run three times with one cold and three incrementally heated Metropolis-coupled Monte Carlo Markov chains, starting from random trees. Output files were analyzed and the convergence and mixing of the independent runs were assessed for all parameters using TRACER v1.4 [40]. Trees from the different runs and their associated posterior probabilities (PP) were then combined and summarized in a 50% consensus tree. All computational analyses were performed using the CIPRES Gateway cloud servers [41]. A clade with a PP value > 0.95 or a BS value > 85% was considered well supported. Additionally, for the subfamily Betoideae and using the concatenation of all loci, the NeighborNet algorithm [42] as implemented in SplitsTree v4.0 [43], was used with the default settings to visualize possible incongruences in the dataset. This method relaxes the assumption that evolution follows a strictly bifurcating path and allows for the identification of reticulated evolution or incomplete lineage sorting among the dataset.
Statistics for the alignments and phylogenetic analyses, as well as the model of evolution for the datasets, are presented in S2 Table. Divergence time analyses Divergence times within the subfamily Betoideae were estimated using the Bayesian MCMC algorithm implemented in BEAST v1.7.2 [44]. For this analysis, we used the combination of the ITS and two cpDNA markers (matK and rbcL), for which outgroup sequences could be obtained. The GTR model of sequence substitution was used for all partitions, except for rbcL, which used the GTR+G model. The fossil of Chenopodipollis multiplex, from a pollen record found in the United States and dated to the early Paleocene , was used to calibrate the root of our phylogenetic tree, which was previously suggested as the best constraint location for the fossil [15]. Therefore, a normal prior was applied to the root of the phylogenetic tree of this study with a mean of 60.5 and a standard deviation of 2, in order to accommodate the fossil age uncertainty . A relaxed lognormal molecular clock was used for all partitions, the Yule process was implemented for the tree prior with a constant speciation rate per lineage, and a random tree was used as the starting tree. The Bayesian MCMC was run for 5x10 7 generations, sampling parameters at every 5000 generations. This analysis was conducted three independent times. Tracer v1.4 [40] was used to assess convergence and correct mixing of all parameters by visually inspecting the log traces and estimating the Effective Sample Size (ESS) of each parameter. Results from the three runs were combined with LogCombiner v1.7.2 [44], after discarding the 10 first % of each analysis as burn-in. The remaining trees were summarized using a Maximum Clade Credibility target tree in TreeAnnotator v1.7.2 [44], as well as Bayesian posterior probability (PP), MEDIAN/ MEAN height and the 95% highest posterior density heights interval (95% HPD) of each node. All computational analyses were performed in the CIPRES Gateway cloud servers [41].

Phylogenetic analyses of Betoideae
Here a new phylogeny of the Betoideae is presented based on nuclear (ITS) and cpDNA markers (matK, trnH-psbA, trnL intron, rbcL), covering a widespread sampling within this subfamily, and an outgroup information from other plants from Amaranthaceae family (S2 and S3 Figs). The results provide support for the (i) monophyly of this subfamily, with the exclusion of the genus Acroglochin; (ii) a deep genetic differentiation between Beta and Patellifolia, which are monophyletic groups; and (iii) the identification of three monophyletic lineages, corresponding to major gene pools of sugar beet CWR (i.e. GP1, GP2 and GP3).
Maximum Likelihood (ML) and Bayesian Inference (BI) were used to test phylogenetic hypotheses within the Betoideae subfamily. Topology of the ML tree (S2 Fig) obtained using the concatenated ITS and cpDNA markers was similar to that obtained from Bayesian analysis (tree not shown). Both ML and BI found essentially identical tree topologies, revealing the same major clades. The monophyly of the Betoideae is well-supported by the data (BS = 100%; PP = 1), but the monotypic genus Acroglochin is excluded from this subfamily. Instead, Acroglochin persicarioides constitutes a robust clade (BS = 90%; PP = 1) with Corispermum chinganicum (Corispermoideae subfamily), and both are closely related to Atriplex prostrata (Chenopodioideae subfamily).
Relationships within the Betoideae subfamily remain somewhat uncertain, since the most basal branches are poorly supported. Therefore, the basal relationships shown among the five genera (i.e. Aphanisma, Beta, Hablitzia, Oreobliton, and Patellifolia) could be interpreted as polytomic, according to our results. Nevertheless, the close relationship between Oreobliton and Aphanisma is well-supported (BS = 100%; PP = 1).
Both Beta and Patellifolia appear to be well-supported monophyletic groups: clade I that includes all samples of the genus Beta (BS = 100%; PP = 1), and clade II (BS = 100%; PP = 1), gathering all the Patellifolia representatives. Within clade I the Beta species found in coastal areas of the Western Mediterranean Region and in the Macaronesian Islands (i.e. B. vulgaris subsp. maritima and subsp. vulgaris, B. macrocarpa, and B. patula) form a well-supported monophyletic group (BS = 99%; PP = 1), which is sister to the remaining members of Beta (i.e. B. corolliflora, B. nana, and B. trigyna) from the Eastern Mediterranean Region (BS = 80%; PP = 1). Moreover, the Macaronesian endemic species, B. patula (from Madeira) and P. webbiana (from the Canary Islands) were placed within Beta (clade I-together with the other species from the GP1) or Patellifolia (clade II-together with the other species from the GP3) clades, respectively. However, our analyses failed to resolve with confidence the relationships among these endemics and the rest of the species resulting in a polytomy (S1 and S2 Figs).
Applying the NeighborNet algorithm to the concatenated dataset reveals a substantial degree of conflicting phylogenetic signal at the divergences of Beta, Patellifolia, Hablitzia, Aphanisma and Oreobliton. This is evidenced by the substantial number of loops found in these points of the phylogenetic network (S3 Fig). Further loops are found within Beta and Patellifolia genera, albeit in a smaller number.

Divergence time analyses
Date estimates for nodes within the subfamily Betoideae are presented in Fig 1 (see C1 to C6). Our analysis indicates that the Betoideae must have diverged around 32.5 million years ago (Mya), representing the split between Hablitzia, and the other four Betoideae genera. Within this group, the split between Beta and Patellifolia (C2) was estimated to have occurred at 25.  Table 2)

Phylogenetic relationships
Our study provides new insights into the phylogenetic relationships within the Amaranthaceae family, and our major findings should help in further refinement of the taxonomy of the subfamily Betoideae. The monophyly of Betoideae was resolved with confidence in our results, but the monotypic genus Acroglochin was excluded from this subfamily, supporting earlier phylogenetic studies [15,16]. This genus occurs in the remote areas of the Himalayas and forms a strongly supported clade with Corispermum chinganicum (subfamily Corispermoideae), which is also distributed in Asian regions. Our molecular data provides evidence for the inclusion of this monotypic genus within the subfamily Corispermoideae, and is consistent with previous works (e.g. [45]). Nevertheless, we considered that further investigation is necessary to effectively test this hypothesis since there is limited taxonomic information currently available for these two genera [29].
The five extant genera of the Betoideae subfamily seem to have a relatively old origin, but their sister relationships within this clade remain unknown. The difficulty in determining the phylogenetic relationships among members of the Betoideae was evidenced by the low support for the basal nodes of this group. Indeed, our phylogenetic network reveals a rather large number of loops at the base of divergence between the five Betoideae genera (see S3 Fig). There are several reasons that may cause such an unresolved phylogenetic pattern, for example a rapid radiation of most genera that leaves little time for the accumulation of mutations and creates a substantial signal of incomplete lineage sorting [46]. Alternatively, ancient hybridization events may have occurred, which created a mosaic pattern of sequence variation. This remains to be seen when more detailed phylogenetic and population data become available, but the data presented in this study is more consistent with Kühn [14] classification, which placed the five genera in one tribe. Our results contradict those of previous works [15,16] which suggested the inclusion of the Beta species within the tribe Beteae, while the Patellifolia species were included in the tribe Hablitzieae, together with three other monotypic genera (i.e. Aphanisma, Hablitzia, and Oreobliton).
Even though our results cannot confidently place Beta and Patellifolia genera relative to each other, their ancient divergence reinforces their recognition as different genera, and this is supported by former morphological studies (e.g. [47]). Therefore, based on our results and previous morphological studies (M.C. Duarte, et al. unpublished data), our study indicated that Patellifolia species, formerly included in Beta section Procumbentes (i.e. B. patellaris, B. webbiana and B. procumbens), should be regarded as a separate genus. Together with previous molecular phylogenetic studies, our results represent a starting point for a thorough taxonomic revision of the subfamily Betoideae.

Spatio-temporal history of Betoideae
The results of the dated molecular phylogeny suggest a relatively old origin for Betoideae, which may have taken place during the Early Oligocene Glacial Maximum (EOGM). The transition . Bayesian tree obtained from the BEAST analysis based on the concatenated dataset of ITS and cpDNA markers (matK and rbcL), illustrating the estimated divergence ages at selected calibrated nodes. Posterior probabilities (PP) are given above each branch. The geographic origin of each specimen is provided (right side) with a color code for continental areas and several Macaronesian archipelagos and grey bars differentiate the three gene pools previously described (for further details see Frese [10]) and concordant with the present phylogenetic analysis. C1 to C6 as described in Table 2 from Eocene to Oligocene was characterized by major climatic changes, which triggered extinctions in plant and animals [48]. Although the relationships among the five genera remain weakly resolved, their early diversification (ca. 32 Mya) tends to support a model of allopatric speciation within this subfamily. This could be the result of past range contractions of the most recent common ancestor of Betoideae, confirmed for instance by long branches of each genus, and reflecting possible speciation by isolation and/or extinction events during the EOGM. At least in line with such a scenario are the narrow distributions, in distant geographic regions, presented by four of the five genera, this being linked to their different ecology. Specifically, Aphanisma occurs in coastal habitats of California; Hablitzia is native in the deciduous forests of the Caucasus Region; Oreobliton is distributed on the chalk rocks of the Atlas Mountains in North-Africa; and Patellifolia is found on coastal vegetation, in maritime rocks, sea-cliffs and seashore habitats in Southern-Western Europe, with its center of diversity in the Macaronesian archipelagos. Conversely, Beta is the only genus of the subfamily Betoideae that presents more species diversity with a broader distribution, mainly throughout the circum-Mediterranean Region [8], and the large range observed in this genus seems to be the result of more recent climatic and geological events. Moreover, a deep genetic differentiation between Beta and Patellifolia species, which may have occurred in the Late Oligocene was correlated with the major gene pools, reflecting an ancient divergence between Beta (GP1, GP2) and Patellifolia (GP3) species. Between the divergence of these two genera and their own diversification, there was around 15-20 Mya of uncertain evolution, considering their respective long branches, and comparing stem and crown ages of these two clades (see Fig 1). Although, our molecular data did not allow us to provide a more reliable evolutionary scenario, Beta and Patellifolia species occur in very constraining living conditions (e.g. aridity, high salinity levels), providing additional support that both lineages may have had more ability to survive the past dramatic aridity events that have occurred within the Mediterranean Region, compared to other more vulnerable plant lineages.
The second biogeographical pattern that was revealed was the occurrence of two well-differentiated clades on each side of the Mediterranean, in the western coastal areas (GP1: B. vulgaris subsp. maritima and subsp. vulgaris, B. macrocarpa, and B. patula) and in the easternmost part of the species' distribution (GP2: B. corolliflora, B. nana and B. trigyna). The Mediterranean Beta species were probably beginning to differentiate around seven million years ago, which matches the Messinian Age of the Late Miocene. This coincides with the Messinian Salinity Crisis (MSC, 5.96-5.33 Mya) [49,50], a period where the connection between the Mediterranean Sea and the Atlantic Ocean closed, causing the Mediterranean Sea to desiccate and probably generating widespread salt marshes or coastal and halophytic habitats across the Mediterranean coast [50]. Such dramatic changes would have promoted the differentiation between GP1 and GP2. These two groups currently occur in different geographical areas and quite differentiated habitat types. GP1 occur in coastal cliffs, salt marshes and ruderal places of the Western Mediterranean Region and Macaronesia Islands, while GP2 is mainly present in continental mountainous zones of the Eastern Mediterranean.
This West-East disjunction pattern has also been found in other plants currently occupying the Mediterranean Basin, both in tree genera such as Laurus L. [51] and Juniperus L. [52], with representatives in Macaronesia, and in herbaceous genera such as Erophaca Boiss. [53]. In these studies current patterns are explained by the contraction of favorable areas, mainly due to an increase of aridity (see [53]) or by the distribution of tectonic microplates and the appearance of water barriers during the Neogene (see [51,52]). Both in the cases of Juniperus and Erophaca the authors suggest a western to eastern speciation sequence, while in Laurus [51], the opposite is hypothesized, with westward expansion of a single haplotype, which colonized across the Western Mediterranean, reaching the Macaronesia Islands.
The subsequent end of the MSC could additionally have promoted further differentiation by vicariance, with Western Mediterranean Beta populations, previously adapted to prevailing salt conditions (e.g. salt marshes habitats), being isolated by loss and fragmentation of this habitat due to post-MSC conditions. Later influential events occurred during the Plio-Pleistocene, with sea level and climate oscillations [54], leading to repeated isolation and connection of taxa, and possible subsequent speciation within the Western Mediterranean Beta. Some of the western wild beets would have later expanded and colonized the Macaronesian Islands. The diaspore adaptations of Beta and Patellifolia species towards sea dispersal (thalassochory) would have promoted their long-distance dispersal and have been clearly advantageous in the colonization of these archipelagos (see [55]). A key role of marine currents in dispersal was also suggested in a recent population study of B. macrocarpa and B. vulgaris subsp. maritima, which encompass the shoreline from France to Morocco [56]. This study suggested that B. vulgaris subsp. maritima went through a postglacial recolonization scenario from the Mediterranean-Atlantic region, with southern Iberia and Morocco, including the Strait of Gibraltar, acting as a long-term refuge.
Altogether, our results support the hypothesis that the Messinian Salinity Crisis and subsequent climatic changes in the Mediterranean Region during the Plio-Pleistocene were probably the major drivers of diversification in the genus Beta, thus explaining the current geographical ranges.

Diversification of wild beets on Macaronesia
The estimation of divergence times provides information on the genetic distance among wild beet species, and facilitates understanding the process and timing of evolution within Beta and Patellifolia species, revealing that the diversification was quite recent, during the Pleistocene. Although this pattern was also reported for other Macaronesian native plant lineages (e.g. [57,58]), the available data does not allow us to recognize this due to insufficient phylogenetic information, and consequently we could not discard the existence of a soft polytomy, meaning that some of the Beta and Patellifolia species may have diverged at different times.
Within the Patellifolia clade there are unresolved polytomies and our study could not infer the monophyly of each species as well as the relationships among P. procumbens (from Madeira, Canary Islands and Cape Verde), P. patellaris (from mainland, Madeira, Canary Islands and Cape Verde) and P. webbiana (from Canary Islands) (see S2 Fig). This pattern can be the consequence of recent island colonization and differentiation, recurrent gene flow with the ancestral mainland populations or congeneric species, or even incomplete lineage sorting. Regarding the latter, the DNA regions sequenced in our study cannot provide a clear resolution for this shallow evolutionary event. Likewise, within the Beta clade, encompassing B. vulgaris subsp. maritima, the cultivated forms and all the Macaronesian species group, the phylogenetic relationships also remain unresolved. Sequences from B. vulgaris subsp. maritima from Madeira Island clustered with B. macrocarpa from Canary Islands, could be explained by introgression or hybridization processes, which are in accordance with the loops observed in our networks analyses (see S3 Fig). A recent study using a flow cytometry analysis, revealed the existence of mixed-ploid populations of B. vulgaris subsp. maritima and B. macrocarpa, in the South of Portugal [59]. Consequently our results suggest that these clustered sequences could reflect an ancient hybridization between the diploids, B. vulgaris subsp. maritima and B. macrocarpa, as was previously suggested by Villain [60]. This author suggested ranking the tetraploid B. macrocarpa from Canary Islands as a separate taxon, and it was proposed that these tetraploid populations result from at least two independent colonization/hybridization events in that archipelago [61]. The increased number of polyploids among these island species can be attributed to the higher adaptive potential of the polyploids [62], which might have been particularly successful in periods of ecological upheaval when new ecological niches were occupied by vigorous polyploids and less competitive diploids were outcompeted [63].
Within the Beta and Patellifolia genera, potential hybridization and the risk of demographic or genetic assimilation of rare endemics (i.e. B. patula in Madeira and P. webbiana in Canary Islands) by other native congener may occur. One possible reason is that the weakness of genetic barriers to hybridization in many islands groups is a by-product of a small genetic differentiation in recently radiated species [64]. As B. patula is classified as Critically Endangered (CR) and is one of the closest wild relatives (GP1) of domestic B. vulgaris subsp. vulgaris, it should be afforded higher conservation priority over the more distantly related species [2]. Thus prioritizing threatened species and conserving the entire extent of their natural ranges was recently recognized as a crucial step towards a better strategy to conserve the endemic flora in the Macaronesia archipelagos [65]. Beyond the actual or potential socio-economic value of these wild relatives as a genetic resource for crop improvement, their extinction would entail the loss of genetic resources that could help such plants overcome the future climatic shifts [66].

Conclusions
This study uncovered the phylogenetic relationships between sugar beet (Beta vulgaris subsp. vulgaris) and the wild species, with particular emphasis on the Beta and Patellifolia species that are commonly found in coastal areas of the Western Mediterranean Region and Macaronesian Islands. The phylogeny recovered on a time-calibrated Bayesian-tree revealed a deep genetic differentiation between Beta and Patellifolia species, which may have occurred in the Late Oligocene. Furthermore, we hypothesized that ecological divergence of Beta in the Mediterranean Basin may have occurred during the Messinian Salinity Crisis (MSC, 5.96-5.33 Mya). Western and Eastern Beta species inhabit very contrasting ecological areas, from salt marshes to mountainous zones respectively. The MSC with its deep habitat modifications and extension could have provided an extraordinary period for Western Mediterranean Beta adaptation to these extreme ecological conditions. The subsequent end of the MSC could additionally have promoted further differentiation by vicariance due to fragmentation and isolation of a previously extended habitat. Some of the western wild beets later expanded and colonized the Macaronesian Islands during the Pleistocene. Moreover, the two endemic taxa (i.e. B. patula and P. webbiana), classified as threatened according to IUCN criteria, are associated with short phylogenetic branches and polytomic groups revealing that the diversification was quite recent in these archipelagos, and unraveling a potentially complex biogeographic pattern with hybridization and gene flow playing an important role. Finally, our phylogenetic analysis of the Betoideae sheds light on the genetic differentiation among the major gene pools of sugar beet wild relatives which are of high evolutionary, ecological, and economic relevance, providing useful data for establishing conservation priorities in the hotspot area of the Macaronesian Islands. We considered that only the conservation of populations in their natural habitats ensures renewal of gene pools and the continued supply of novel genetic material potentially critical for future crop improvement, which is recognized as an asset in maintaining global food security.