Next Article in Journal
Habitat Suitability Analysis and Future Distribution Prediction of Giant Panda (Ailuropoda melanoleuca) in the Qinling Mountains, China
Next Article in Special Issue
Diversity and Distribution of Monocot Understory Herbs during Tropical Forest Succession in Northeastern Costa Rica
Previous Article in Journal
Root-Zone Bacterial Diversity in Field-Grown Individual Plants from Alfalfa Lines with Wild Relatives in Their Genetic Backgrounds
Previous Article in Special Issue
Hilly, Semi-Mountainous and Mountainous Areas Harbor Landraces Diversity: The Case of Messinia (Peloponnese-Greece)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Genetic Characterization of the Canarian Endemic Palm (Phoenix canariensis) by Simple Sequence Repeats and Chloroplast Markers: A Tool for the Molecular Traceability of Phoenix Hybridization

by
Isabel Saro
1,
Priscila Rodríguez-Rodríguez
1,
Diego Rivera
2,
Concepción Obón
3,
Fredérique Aberlenc
4,
Antonio Díaz-Pérez
5,6,
Salwa Zehdi-Azouzi
7,
Leticia Curbelo
1 and
Pedro A. Sosa
1,*
1
Instituto Universitario de Estudios Ambientales y Recursos Naturales (IUNAT), Universidad de Las Palmas de Gran Canaria, Campus de Tafira, 35017 Las Palmas, Spain
2
Departamento de Biología Vegetal (Botánica), Faculty of Biología, Universidad de Murcia, 30100 Murcia, Spain
3
Departmento de Biología Aplicada, Escuela Técnica Superior de Orihuela, Universidad Miguel Hernández de Elche, Ctra. Beniel Km, 3.2, 03312 Orihuela, Spain
4
Plant Diversity, Adaptation and Development, Université de Montpellier, Institut de Recherche pour Développement, Centre de Coopération Internationale en Recherche Agronomique pour le Développement, 911 Av. Agropolis, BP 64501, 34394 Montpellier CEDEX 5, France
5
Gestión y Planeamiento Territorial y Medioambiental S.A. C/León y Castillo 54, 35003 Las Palmas de Gran Canaria, Spain
6
Facultad de Agronomía, Instituto de Genética, Universidad Central de Venezuela, 1050 Maracay, Venezuela
7
Laboratoire de Génétique Moléculaire, Faculté des Sciences de Tunis, Immunologie et Biotechnologie (LR99ES12), Université de Tunis El Manar, Campus Universitaire Farhat Hached, Tunis 1068, Tunisia
*
Author to whom correspondence should be addressed.
Diversity 2024, 16(7), 411; https://doi.org/10.3390/d16070411
Submission received: 25 June 2024 / Revised: 9 July 2024 / Accepted: 10 July 2024 / Published: 16 July 2024
(This article belongs to the Special Issue Plant Diversity Hotspots in the 2020s)

Abstract

:
The endemic palm from the Canary Islands, Phoenix canariensis, is one of the most distinctive elements of the Canarian vegetation landscape, contributing to cultural, economic and environmental aspects. One of the main conservation problems facing this iconic palm is anthropogenic hybridization with other Phoenix species, particularly Phoenix dactylifera, which has been introduced extensively throughout its geographical range. Therefore, it is important to obtain a genetic tool that addresses different issues that may have an impact on the protection of P. canariensis, including ornamental applications and wild population conservation purposes. Our main goals were to detect a molecular tracer that could reliably distinguish between Phoenix canariensis and P. dactylifera in the Canary archipelago and to characterize the presence and extent of genetic hybridization events between the two species. We used 19 nuclear microsatellites and 1 chloroplast minisatellite set and analysed a large sample size (N = 433) of plants using both Bayesian methods and ordination techniques. Our data showed that a set of 13 nuclear markers revealed diagnostic alleles for P. canariensis, which were defined as the Canarian nuclear genotype (CNG). Moreover, P. canariensis exhibited an exclusive chlorotype of 266 bp that together with the GNC serve as an indicator of genetic purity in the Canarian palm. These markers are sufficient to detect any hybrid, even if it is not related to morphological differences. The occurrence of a considerable number of specimens with different degrees of hybridization is discussed in terms of the existence of different generations of hybrids and different types of crosses. Thus, the genetic tracers represent an invaluable tool to address any proposal for the genetic conservation of Phoenix canariensis.

1. Introduction

The molecular characterization of the genetic patterns of the Phoenix species has been subjected to numerous studies in the last decades, mainly focused on varieties of Phoenix dactylifera due to their great economic relevance [1]. Frequent hybridization between species and the scarcity of suitable taxonomic and morphological characters to distinguish different Phoenix species have motivated the search for a useful molecular marker [2,3,4,5]. Hypervariable simple sequence repeats (SSRs) and minisatellite markers have been well adapted to generate barcoding markers given the lack of sufficient sequence variation in both the nuclear and chloroplastic genomes to discriminate between Phoenix species [3]. Currently, a battery of highly informative markers has been developed to assess genetic processes and to identify hybrids among very closely related taxa such as Phoenix species [3,6,7,8]. In addition, these molecular markers have been proved to be effective in several Phoenix studies to reconstruct phylogenetic trees [3,9,10], to clarify genetic structure across biogeographic ranges of species [5,8,11,12,13,14,15], to determine paternal linages through sex-specific genetic markers [16] and even to identify wild relatives of P. dactylifera [17].
Phoenix canariensis, the Canarian date palm is one of the most representative endemic plant species from the Canarian archipelago, being designated as its plant symbol by the Canarian Government [18]. It is one of the most distinctive elements of the Canarian vegetation landscape, contributing to cultural, economic and environmental aspects. Also, the Canarian palm has an important economic relevance as ornamental palm, and it is widely distributed across all the archipelago, but it is also planted across the subtropical and tropical world for its ability to grow on a wide range of site types [19]. Different Phoenix species have been introduced in the Canary Islands as ornamental species. P. dactylifera is by far the most widespread species. It was possibly introduced across different Canarian islands since the aboriginal occupation from Northern Africa 2000 years ago [20]. Later, it was widely dispersed across the archipelago considering its agricultural and ornamental importance. This scenario has promoted anthropogenic P. dactylifera geneflow into P. canariensis, generating hybridization processes in geographic areas in which both species overlap, threatening the genetic integrity of the island endemism. There is also hybridization with other Phoenix species, but it is not so relevant and is of a much lesser extent than that with P. dactylifera [18,21,22]. In fact, in the Canarian archipelago, Phoenix dactylifera is included in the list of alien and invasive species following the Spanish Catalogue of Invasive Alien Species regulation (BOE No. 185, 3 August 2013). Therefore, conservation management of P. canariensis must preserve and expand genetically pure populations since the genetic value of hybrid populations is considerably reduced due to genetic erosion of the endemic species [23]. In fact, the Canarian Island Government published a specific decree (Decree 62/2006 of 16 May 2006) establishing measures to promote the protection and conservation of the genetic identity of the P. canariensis but did not specify the methods and the requirements to be met. Moreover, the morphological similarity and the great ecological and morphological plasticity between both species of the genus Phoenix make their classification and the morphological identification of pure and hybrid populations a very subjective task currently. In this sense, an assessment of interspecific hybridization and introgression between both species must be implemented to develop an appropriate genetic conservation strategy for P. canariensis, compelling researchers to generate genetic tracers that unambiguously identify the existence of genetic hybridization events in the Canarian date palm germplasm.
Previous genetic studies have demonstrated the existence of hybrids between the endemic Phoenix canariensis and the widespread P. dactylifera in the Canary Islands using RAPD markers [18,21]. Other attempts to assess genetic hybridization between both species using isozyme markers supported considerable genotype differentiation between both species [22] but with limited discriminatory power due to low levels of detectable genetic variability [24]. Nevertheless, previously detected SSR marker polymorphism in Phoenix species could be used to fix RAPD and isoenzyme drawbacks, being especially useful to identify genetic purity tracers of P. canariensis [3]. Additionally, [25] have identified a polymorphic cpDNA minisatellite locus composed of 12 bp repeats located between the trnG(GCC)-trnfM(CAU) intergenic spacer, whose haplotypes were strongly associated with the taxonomy of Phoenix species [3,26,27,28]. Thus, P. canariensis seems to harbour a specific haplotype characterized by a 5X-tandem-repeat motif resulting in a fragment of 266 pb [3,26]. On the other hand, the date palm always exhibits two chlorotypes associated with an Occidental (242 pb fragment size) and an Oriental (254 pb fragment size) geographic distribution, respectively [1,3,8,17]. Overall, the combined implementation of nuclear and chloroplast markers represents a good starting point to characterize the genetic identity of the P. canariensis for genetic conservation of the species.
Nevertheless, although the date palm has been extensively sampled and analysed with the chloroplast minisatellite [1,8,14], only a small sample of P. canariensis individuals has been studied together with other Phoenix species [3,26,27], mostly from cultivated palms of France, Italy, United Kingdom and Spain. In this sense, this sample represents a minimal and biased fraction of the wild genetic pool of P. canariensis populations, which, in addition, could have originated from the wild restricted areas of La Palma and Northern Tenerife islands [13,18]. Therefore, it makes sense to consider that other haplotypes could have not been sampled across the archipelago.
In this study, we utilized nuclear SSR and chloroplast minisatellite markers to characterize genetic hybridization events between the Canarian endemic species and its Continental relative P. dactylifera. Specifically, we addressed the following aims: (i) to evaluate a group of nuclear and plastid markers to discriminate among Phoenix canariensis, P. dactylifera and their hybrids; (ii) to analyse the utility of this set of genetic markers to characterize the extent and the grade of the genetic introgression in the wild populations of the Canarian palm and (iii) to identify the taxonomic specificity of the Phoenix canariensis chlorotype employing extensive sampling from natural populations. The use of these molecular markers will serve as a powerful tool to be incorporated into genetic conservation programs of Phoenix canariensis.

2. Materials and Methods

2.1. Plant Material

We used a total of 433 samples from Phoenix: 271 of P. canariensis, 88 of wild putative hybrids distributed in natural and seminatural palm groves of Tenerife and Gran Canaria islands (Canary Islands), 26 of P. dactylifera specimens from Gran Canaria and 48 non-Canary Islands P. dactylifera from Tunisia and the Iberian Peninsula. The taxonomic classification of samples was based on carefully scored morphological traits such as form and colour of the leaf, trunk shape, spines sizes, presence of offshoots and shape of the crown and trunk [2,4,18] (Table 1 and Figure S1).
Samples were collected from two different islands to provide a wider range of morphological and genetic variation in P. canariensis and putative hybrids. Specimens with intermediate morphological features between P. canariensis and P. dactylifera were considered as putative hybrids. Moreover, we expected that P. dactylifera specimens would exhibit the Oriental (254 bp fragment size) or Occidental chlorotype (242 bp), as these had been previously detected in this species across its distribution [1,8,9,17]. Overall, we assumed six morpho-geographic types according to a combination of morphological features and geographic sampling origin: Phoenix canariensis morphotype from Gran Canaria (PCGC, N = 169) and Tenerife (PCTF, N = 102); putative hybrids from Gran Canaria (HBGC, N = 51) and Tenerife (HBTF, N = 37); and Phoenix dactylifera morphotype with eastern chlorotype (PD254, N = 31) and western chlorotype (PD242, N = 43).

2.2. DNA Isolation and Molecular Analysis

Total genomic DNA was extracted from approximately 50 mg of silica-dried leaf using the Invisorb® DNA Plant HTS 96 Kit, following the manufacturer’s protocol. Subsequently, the DNA templates were amplified with 19 nuclear microsatellite loci previously developed for P. dactylifera, including dinucleotide repeats from an SSR genomic library [6], introns coding for transcription factors [29], the cup-shaped cotyledon 3-related protein (CUC3) gene and the Pd-AP3-SSR1 sequence [8]. We also included tri-/hexanucleotide repeats from coding sequences identified by in silico mining [7] in the whole-date-palm genome sequence [30]. The chlorotypes of all the palms were characterized according to the number of scored repeats of the dodecanucleotide minisatellite located in the intergenic spacer psbZ-trnfM [25]. Amplification reactions and alleles size scoring were performed as in [8], using a combination of all the loci to define an exclusive multilocus genotype for each sample.

2.3. The Genetic Characterization of the Hybridization Process between Phoenix canariensis and Phoenix dactylifera

The nuclear allele and chlorotype frequencies and sizes were calculated with GENALEX v6.5 [31]. To explore the molecular relationships among the samples and infer pure and putative hybrid specimens according to their group or intermediate membership in the main components axis, a principal component analysis (PCA) was carried out on the allele’s frequency according to the dudi.pca function implemented in the ade4 package for R software v2.15.3 (R Development Core Team, Vienna, Austria, 2011).
To infer the genetic ancestry of multilocus genotypic data, we applied two Bayesian clustering methods. First, the STRUCTURE 2.3.4 software [32,33] was used to detect admixture between P. canariensis and P. dactylifera, considering that both species were well differentiated into two clusters under an admixture and uncorrelated alleles frequency ancestral model. Then, we assumed that pure samples showed a priori ancestry (q-value) greater than 95% for one of the two clusters according to the Popinfo setting of [32]. Admixed ancestry is modelled by assuming that an individual, i, has inherited some fraction (q) of its genome from any of the two ancestors. Thus, hybrids could be inferred when observing individuals with q-values less than 95% for any of the two clusters. We used Markov chain Monte Carlo (MCMC) iterations as implemented in the software’s algorithm to explore a parameter space fixing individual memberships to a K = 2 cluster (pre-defined pure species class), which was replicated 10 independent times for each K, each with a burn-in and chain length of 105 and 106 replicates, respectively. Second, to further validate our STRUCTURE analysis, we characterized specimens according to the posterior distribution of individual assignment into pure parental (P. canariensis, pure I class; P. dactylifera, pure II class) and hybrid categories (F1, F2 or backcross with pure species I and II) as implemented in NEWHYBRIDS [34]. We used default genotypic classes with no prior information on allelic frequencies and included uniform and Jeffereys priors for θ and π following the default setting [34]. Five runs were conducted with 105 sweeps of burn-in and 107 sweeps of data collection. We calculated Pearson’s correlation between STRUCTURE’s proportion of membership (qi) and NEWHYBRIDS’ posterior probabilities (pps) to visualize the relationship between these two analyses, using R software v2.15.3.

2.4. The Feasibility of Molecular Markers for the Detection of a Hybrid Signal in Phoenix canariensis

The molecular genotypes were checked for diagnostic alleles to assess the purity and hybrid status of the samples of Phoenix according to three strategies. First, we identified diagnostic nuclear alleles that had frequencies greater than 0.9 in P. canariensis PCGC and PCTF and were not detected in P. dactylifera PD254 and PD242, which also contributed the most to the PCA load plot specified in the previous section. Such alleles were considered as a tracer of absolute genetic purity in the Canarian endemism, defining altogether the Canarian nuclear genotype (hereafter, CNG). Thus, a measure of the hybridization extent was established as the expected loss of purity in natural populations of Phoenix, assuming different molecular types according to the proportion of CNG alleles: (i) pure P. canariensis: higher than 90% of CNG alleles; (ii) hybrid: between 10% and 90% of CNG alleles; (iii) pure P. dactylifera: lower than 10% of CNG alleles. Second, we analysed the frequency and distribution of different expected chlorotypes along the molecular types. Third, we used a q-value threshold to make a distinction between hybrid and pure individuals in STRUCTURE analysis, as it has been recommended for conservation purpose [35]. So, individuals with a q-value between 0 and <0.10 or >0.9 and 1 were classified as pure, and any individuals with a q-value between 0.1 and 0.9 were classified as hybrid. Also, we computed the proportion of times we observed a correspondence between the morpho-geographic types and the molecular types; that is, if the following condition is met: PCGC + PCTF = pure P. canariensis; HBGC + HBTF = hybrids; PD254 + PD242 = pure P. dactylifera.

3. Results

3.1. The Genetic Characterization of Phoenix canariensis and P. dactylifera and Their Hybrids

According to the 20 SSR loci, the PCA separated the samples of P. canariensis (PCGC + PCTF) and P. dactylifera (PD254 + PD242) morpho-geographic types into different groups and differentiated the within-island genetic structure in Canarian palm samples (Figure 1a).
Both species were clearly separated by the combination of the first two axes of the PCA (44.4% of total variance), indicating a strong genetic differentiation between them. Samples located at intermediate positions of the PCA, suggesting a hybrid genetic composition, were mostly related to specimens with the putative hybrid morpho-geographic types HBGC and HBTF. The PC2 discriminated P. canariensis samples according to island origins, separating specimens from Gran Canaria (PCGC) and Tenerife (PCTF) (Figure 1a). On the other hand, both Bayesian inference analyses showed strong agreement with each other, supported by a significant correlation between the ancestry coefficients generated by STRUCTURE and NEWHYBRIDS (r = 0.952; p > 0.0001). Both software revealed two distinct clusters associated with the P. canariensis and P. dactylifera morpho-geographic types, respectively, whose samples showed assignation scores of q (STRUCTURE) or pp (NEWHYBRIDS) higher than 0.99. Also, substantial genetic introgression levels were detected in some samples, which were assigned to a greater extent to F2 genotype classes according to STRUCTURE, or backcrosses to pure P. canariensis as indicated by NEWHYBRIDS’ outputs (Figure 2).

3.2. The Feasibility of SSR Markers for the Detection of a Hybrid Signal in Phoenix canariensis

Some SSR loci revealed very similar allelic frequencies between the three morpho-geographic types (e.g., PdCUC3-ssr1, mPdIRD013 or mPdIRD057). However, seven loci showed alleles with contrasting frequencies between P. canariensis (frequency > 0.9) and P. dactylifera (frequency = 0). These loci were mPdCIR015 (P15A and P15B), mPdCIR016 (P16), mPdCIR035 (P35), PdAP3-ssr-F4 (F4), mPdIRD033 (P33) and PdAG1-ssr (PAG1) (Table 2).
Some of these alleles were found in varying proportions in the putative hybrids. The locus mPdCIR015 (P15) presented four peaks in some individuals of P. canariensis and in some of the putative hybrids; therefore, it was considered as a putative duplication, being treated as two independent loci (P15A and P15B). Interestingly, none of the specimens considered as molecularly pure P. dactylifera had duplication of the mPdCIR015 (P15) locus. That is, they all had only two alleles (which were included in the P15A locus; Table 2). However, all P. canariensis specimens and all molecular hybrid specimens had four alleles for this P15 locus (distributed in loci P15A and P15B; Table 2). Other loci, mPdCIR010 (P10), mPdCIR063 (P63), mPdCIR085 (P85), mPdIRD031 (P31) and mPdIRD040 (P40), harboured more than one predominant allele in P. canariensis, but all these alleles were exclusives of Canarian species. Moreover, the sum of these allele frequencies ranged between 0.838 for mPdCIR010 (P10) and 0.991 for mPdIRD040 (P40) in P. canariensis, and they were not detected in P. dactylifera (Table 2). The loci and alleles of this group were considered as a tracer of genetic purity in the Canarian endemism, defining altogether the CNG. Therefore, we classified the samples into three different molecular types according to the proportion of CNG alleles: Pure P. canariensis specimen, showing a total of between 22 and 24 alleles of the CNG (i.e., 90–100% of the CNG alleles). On the contrary, a pure P. dactylifera specimen harboured a range of alleles between zero and two CNG alleles (i.e., 0–10%). Between these two values (9 to 21 CNG alleles), we classified samples as hybrid specimens but with different degrees of hybridization or genetic introgression. No hybrids were found harbouring between three and eight CNG alleles (Table 3).
The vast majority of specimens of pure P. canariensis with 22–24 CNG alleles also harboured chlorotype 266. Only 11 of 275 specimens with 22–24 CNG alleles (4%) showed chlorotype 242, but they had morphological traits of P. canariensis. Then, among the samples with 10–22 CNG alleles that were classified as hybrids, 47 samples (61.0%) harboured chlorotype 266, 35.1% harboured chlorotype 242 and only 3 hybrid specimens (3.9%) had the eastern chlorotype (C254). In addition, the Bayesian inference largely agreed with the molecular types established on the samples. Furthermore, the samples attributed to hybrids with the date palm chlorotype (Clp242 and Clp254) were more consistently assigned to Bayesian hybrid clusters than those with the Canarian palm chlorotype (Clp266). This seems to indicate that introgression from date palms in Canarian palms with Clp266 is less detectable in appearance than that in hybrids and introgression from mother date palm linages (e.g., Clp242 and Clp254).
Lastly, when comparing the established morpho-geographic and molecular types, in general, a high proportion (85.1% to 92.7%) of pure molecular types agree with the morpho-geographic type, suggesting a scenario in which it is likely to detect the correct species by the morphological traits (Table 3). No individual morphologically classified as P. dactylifera was included as a pure P. canariensis molecular type and vice versa; only a few proportions (7.3% to 14.9%) of pure molecular type individuals were wrongly attributed to hybrid morpho-geographic type individuals. However, the matching success significantly decreases when we consider the hybrid molecular type, in which only 51.9 were correctly assigned to putative hybrids based on their morphological traits (Table 3).

4. Discussion

In this study, the molecular markers used generated substantial genetic differentiation to reliably characterize P. canariensis and P. dactylifera and infer their putative hybrids. Of the 217 alleles detected, more than half were specific to one of the two species, as suggested by the allelic loadings of the PCA (Figure 1b). Overall, our data supported previous expectations in the Phoenix genus about the ability of the SSR markers to discriminate Phoenix species [3,7,8,10]. Out of the 20 loci, 12 of them, named the Canarian nuclear genotype (CNG), exhibited 30 specific alleles for P. canariensis and, consequently, were considered as genetic tracers of the purity of this species, which appeared at different frequencies in the hybrid individuals (Table 2). Nevertheless, this set of diagnostic alleles should be treated as sample-specific markers to the P. canariensisP. dactylifera complex, since some alleles not found in P. dactylifera could exist in other Phoenix species.
We observed that the most frequent chlorotype within the P. canariensis pure type was the 266 bp fragment (Cl266), confirming its specificity in the Canarian palm (Figure S2), which agrees with previous studies with smaller sample sizes representing the Canarian endemism [3,26,27,28]. Nevertheless, some exceptions to this trend were observed in a very small sample of 11 specimens with the chlorotype 242, showing at the same time the CNG and morphological traits of P. canariensis (Table 3). These specimens were detected in three localities of Gran Canaria, Acusa (7), Telde (1) and Arucas (3), which coincidentally harbour several of the hybrids detected in this study and wild date palms, coinciding with the fact that specimens of uncertain origin have been planted in the surrounding area. We did not find this molecular type in any other population from Tenerife or the La Gomera islands, despite the large number of Phoenix specimens analysed so far by our research group. Considering that we cannot completely rule out that these specimens could be pure Canary palm individuals with the chlorotype 242, we are inclined to think that they could be in fact ancestral hybrids, whose primary parents were a female P. dactylifera (Clp242) and a male P. canariensis (Clp266). Consequently, current individuals were then generated by continuous crossbreeding by several generations of backcrossing to the more abundant pure P. canariensis palms, diluting the original hybrid signal in the nuclear genome but maintaining the 242 chlorotype by maternal chloroplast inheritance. The discordance among the nuclear and plastid genes is very common in hybrid lineages, a process referred to as chloroplast capture, which would perfectly explain these results [36,37,38].
The rest of the hybrids found fit into a very varied casuistry, although all 77 specimens detected as true molecular hybrids seem to come from the crosses between specimens of P. canariensis and P. dactylifera, since we only detected the exclusive chlorotype of these two species (Clp242, Clp254 or Clp266), but none from other Phoenix species [26]. The Bayesian analyses showed that the detection of backcross hybrids and evidence of historic introgression between P. dactylifera and P. canariensis are complex, and not limited to the F1. It is very difficult to determine the exact origin of the P. dactylifera specimens introduced into the Canary Islands. However, it is important to mention that the aboriginal Canary Islanders were Berbers from North Africa, who had relevant knowledge of agriculture. They arrived on the islands some almost 1800 years ago, with goats, cultivated various vegetable “crop packages”, especially hulled barley (Hordeum vulgare), durum wheat (Triticum durum) or lentil (Lens culinaris) among others [39]; therefore, we cannot rule out the possibility that they brought seeds of Phoenix dactylifera from their places of origin (North of Africa) and cultivated them. Under this scenario, it is likely to think that hybrid palms were generated more than 1500 years ago. More recently, we do know, however, that most of the palms imported into the Canary Islands as ornamental elements are mainly from the Mediterranean coast, where the western date palm type prevails. This is possibly part of the reason why most of the date palms in the Canary Islands are Clp242, with a scarce representation of eastern individuals. However, it is interesting to note that the western chlorotype of date palm (Clp242) was not observed in samples collected in Tenerife (Figure S2). Furthermore, the eastern specimens found in Tenerife are in the same locality (Rambla de Castro), a cultivated population that was created in a very recent period. Anyway, taking to account the great amount of F2 hybrids or backcrosses with P. canariensis detected in this study, it would be interesting to study the evolutionary implications of late-generation hybrids such as a possible loss of fertility in P. canariensis, as well as the environmental conditions that favour the anthropogenic hybridization [40].
Hybrids and specimens of P. dactylifera are widely distributed throughout the Canary Islands, either for historical reasons, in which the aboriginal Canary Islanders brought seeds of P. dactylifera, or due to their very recent introduction in tourist and urban areas as ornamentals. Most of these specimens are found in urban and peri-urban areas, where massive planting of these specimens has been very important. The effect of these urban specimens on the genetic conservation of the natural populations of P. canariensis is nil, due to the remoteness and, therefore, certain isolation from the natural palm groves. However, many other hybrid specimens (and date palms) grow in areas close to pure P. canariensis palm groves, so the threat of hybridization could progressively increase in these areas. The lack of knowledge and the difficulty of morphologically identifying the species to which the young individuals belong has led to the plantations established in semi-natural areas becoming areas of strong hybridization (for example, Rambla de Castro in Tenerife or Acusa in Gran Canaria).
In conclusion, we detected a molecular marker set that unambiguously serves as a tool to identify and separate both species and their hybrids, especially when morphological identification is in doubt. We also recommend the use of this molecular set to check the levels of genetic introgression in the natural palm groves of the Canary Islands. The use of this group of molecular markers could be implemented permanently in the official identification of the natural palm groves of Phoenix canariensis as seed sources for the conservation of the genetic resources of Phoenix canariensis in the Canary Islands archipelago.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/d16070411/s1: Figure S1: Morphological differences between Phoenix dactylifera and P. canariensis [41]; Figure S2: Number of individuals of different chlorotypes for each population; Figure S3: Number of individuals in each molecular category, according to the number of CNG alleles.

Author Contributions

I.S.: Conceptualization, Methodology, Software, Data Curation, Writing— Original Draft Preparation. P.R.-R.: Investigation, Writing—Review and Editing, Visualization. D.R.: Methodology, Data Curation, Writing—Review and Editing. C.O.: Methodology, Data Curation, Writing—Review and Editing. F.A.: Resources. A.D.-P.: Data curation, Writing—Review and Editing. S.Z.-A.: Resources. L.C.: Resources, Investigation. P.A.S.: Conceptualization Writing—Original draft Preparation, Supervision, Funding Acquisition. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Autonomous Government of the Canary Islands; Agencia Canaria de Investigación, Innovación y Sociedad de la Información (Project PALMHYBRID), grant number ProID2020010029; Programa Operativo FEDER Canarias 2014–2020 and Consejería de Medio Ambiente y Emergencias of Cabildo de Gran Canaria (Realización de estudio genético que caracterice las fuentes semilleras de los palmerales silvestres de Gran Canaria; grant number 267).

Institutional Review Board Statement

Not applicable.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors on request.

Acknowledgments

We thank all the staff of the Consejerías de Medio Ambiente de los Cabildos de Gran Canaria and Tenerife for their collaboration in the collection of the samples. We also thank Agustín Naranjo Cigala for the geographical location of the samples.

Conflicts of Interest

The authors declare no conflicts of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. Gros-Balthazard, M.; Baker, W.J.; Leitch, I.J.; Pellicer, J.; Powell, R.F.; Bellot, S. Systematics and Evolution of the Genus Phoenix: Towards Understanding Date Palm Origins. In The Date Palm Genome, Vol. 1. Compendium of Plant Genomes; Springer: Cham, Switzerland, 2021; Volume 1, pp. 29–54. ISBN 978-3-030-73746-7. [Google Scholar]
  2. Barrow, S.C. A Monograph of Phoenix L. (Palmae: Coryphoideae). Kew Bull. 1998, 53, 513–575. [Google Scholar] [CrossRef]
  3. Pintaud, J.-C.; Zehdi, S.; Couvreur, T.; Barrow, S.; Henderson, S.; Aberlenc-Bertossi, F.; Tregear, J.; Billotte, N. Species delimitation in the Genus Phoenix (Arecaceae) based on SSR Markers, with emphasis on the identity of the Date Palm (Phoenix dactylifera). In Diversity, Phylogeny and Evolution in the Monocotyledons; Seberg, O., Petersen, G., Barford, A., Davis, J., Eds.; Aarhus University Press: Aarhus, Denmark, 2010; pp. 267–285. [Google Scholar]
  4. Gros-Balthazard, M. Hybridization in the genus Phoenix: A review. Emir. J. Food Agric. 2013, 25, 831–842. [Google Scholar] [CrossRef]
  5. Saro, I.; García-Verdugo, C.; González-Pérez, M.A.; Naranjo, A.; Santana, A.; Sosa, P.A. Genetic Structure of the Canarian Palm Tree (Phoenix canariensis) at the Island Scale: Does the ‘Island within Islands’ Concept Apply to Species with High Colonisation Ability? Plant Biol. 2019, 21, 101–109. [Google Scholar] [CrossRef] [PubMed]
  6. Billotte, N.; Marseillac, N.; Brottier, P.; Noyer, J.L.; Jacquemoud-Collet, J.P.; Moreau, C.; Couvreur, T.; Chevallier, M.H.; Pintaud, J.C.; Risterucci, A.M. Nuclear Microsatellite Markers for the Date Palm (Phoenix dactylifera L.): Characterization and Utility across the Genus Phoenix and in Other Palm Genera. Mol. Ecol. Notes 2004, 4, 256–258. [Google Scholar] [CrossRef]
  7. Aberlenc-Bertossi, F.; Castillo, K.; Tranchant-Dubreuil, C.; Chérif, E.; Ballardini, M.; Abdoulkader, S.; Gros-Balthazard, M.; Chabrillange, N.; Santoni, S.; Mercuri, A.; et al. In Silico Mining of Microsatellites in Coding Sequences of the Date Palm (Arecaceae) Genome, Characterization, and Transferability. Appl. Plant Sci. 2014, 2, 1300058. [Google Scholar] [CrossRef] [PubMed]
  8. Zehdi-Azouzi, S.; Cherif, E.; Moussouni, S.; Gros-Balthazard, M.; Naqvi, S.A.; Ludeña, B.; Castillo, K.; Chabrillange, N.; Bouguedoura, N.; Bennaceur, M.; et al. Genetic Structure of the Date Palm (Phoenix dactylifera) in the Old World Reveals a Strong Differentiation between Eastern and Western Populations. Ann. Bot. 2015, 116, 101–112. [Google Scholar] [CrossRef] [PubMed]
  9. Pintaud, J.C.; Ludeña, B.; Aberlenc-Bertossi, F.; Gros-Balthazard, M.; Ivorra, S.; Terral, J.F.; Tengberg, M.; Hernández, I.S.; González-Pérez, M.A.; Sosa, P.; et al. Biogeography of the Date Palm (Phoenix dactylifera L., Arecaceae): Insights on the Origin and on the Structure of Modern Diversity. Acta Hortic. 2013, 994, 19–38. [Google Scholar] [CrossRef]
  10. Zehdi, S.; Cherif, E.; Rhouma, S.; Santoni, S.; Hannachi, A.S.; Pintaud, J.C. Molecular Polymorphism and Genetic Relationships in Date Palm (Phoenix dactylifera L.): The Utility of Nuclear Microsatellite Markers. Sci. Hortic. 2012, 148, 255–263. [Google Scholar] [CrossRef]
  11. Elshibli, S.; Korpelainen, H. Microsatellite Markers Reveal High Genetic Diversity in Date Palm (Phoenix dactylifera L.). Germplasm from Sudan. Genetica 2008, 134, 251–260. [Google Scholar] [CrossRef]
  12. Khanam, S.; Sham, A.; Bennetzen, J.L.; Aly, M.A.M. Analysis of Molecular Marker-Based Characterization and Genetic Variation in Date Palm (Phoenix dactylifera L.). Aust. J. Crop Sci. 2012, 6, 1236–1244. [Google Scholar]
  13. Saro, I.; González-Pérez, M.A.; García-Verdugo, C.; Sosa, P.A. Patterns of Genetic Diversity in Phoenix canariensis, a Widespread Oceanic Palm (Species) Endemic from the Canarian Archipelago. Tree Genet. Genomes 2015, 11, 815. [Google Scholar] [CrossRef]
  14. Zango, O.; Cherif, E.; Chabrillange, N.; Zehdi-Azouzi, S.; Gros-Balthazard, M.; Naqvi, S.A.; Lemansour, A.; Rey, H.; Bakasso, Y.; Aberlenc, F. Genetic Diversity of Southeastern Nigerien Date Palms Reveals a Secondary Structure within Western Populations. Tree Genet. Genomes 2017, 13, 75. [Google Scholar] [CrossRef]
  15. Bourguiba, H.; Hachef, A.; Saffar, W.; Chabrillange, N.; Aberlenc, F.; Cherif, E.; Zehdi-Azouzi, S. Male Diversity Matters: Genetic Structuring of Insular Male Date Palm (Phoenix dactylifera L.) Revealed Valuable Breeding and Conservation Resources. J. Hortic. Sci. Biotechnol. 2024, 99, 212–222. [Google Scholar] [CrossRef]
  16. Cherif, E.; Zehdi, S.; Castillo, K.; Chabrillange, N.; Abdoulkader, S.; Pintaud, J.C.; Santoni, S.; Salhi-Hannachi, A.; Glémin, S.; Aberlenc-Bertossi, F. Male-Specific DNA Markers Provide Genetic Evidence of an XY Chromosome System, a Recombination Arrest and Allow the Tracing of Paternal Lineages in Date Palm. New Phytol. 2013, 197, 409–415. [Google Scholar] [CrossRef] [PubMed]
  17. Saffar, W.; Bourguiba, H.; Cherif, E.; Chabrillange, N.; Hachef, A.; Ivorra, S.; Moussini, S.; Naqvi, S.A.; Littardi, C.; Zango, O.; et al. Towards the Identification of Wild African Date Palms in Kerkennah Islands from Tunisia by Morphometric and Genetic Analyses. Sci. Hortic. 2022, 297, 110973. [Google Scholar] [CrossRef]
  18. Sosa, P.A.; Saro, I.; Johnson, D.; Obón, C.; Alcaraz, F.; Rivera, D. Biodiversity and Conservation of Phoenix canariensis: A Review. Biodivers. Conserv. 2021, 30, 275–293. [Google Scholar] [CrossRef]
  19. Spennemann, D.H.R. Canary Island Date Palms (Phoenix canariensis) in Australia: Introduction and Early Dispersal. Palms 2018, 62, 185–201. [Google Scholar]
  20. Zona, S. The Horticultural History of the Canary Islands Date Palm (Phoenix canariensis). Gard. Hist. 2008, 36, 301–309. [Google Scholar]
  21. González-Pérez, M.A.; Caujapé-Castells, J.; Sosa, P.A. Allozyme Variation and Structure of the Canarian Endemic Palm Tree Phoenix canariensis (Arecaceae): Implications for Conservation. Heredity 2004, 93, 307–315. [Google Scholar] [CrossRef]
  22. González-Pérez, M.A.; Sosa, P.A. Hybridization and Introgression Between the Endemic Phoenix canariensis and the Introduced P. dactylifera in the Canary Islands. Open For. Sci. J. 2009, 2, 78–85. [Google Scholar] [CrossRef]
  23. Allendorf, F.W.; Leary, R.F.; Spruell, P.; Wenburg, J.K. The Problems with hybrids: Setting conservation guidelines. Trends Ecol. Evol. 2001, 16, 613–622. [Google Scholar] [CrossRef]
  24. González-Pérez, M.A.; Caujapé-Castells, J.; Sosa, P.A. Molecular Evidence of Hybridisation between the Endemic Phoenix canariensis and the Widespread P. dactylifera with Random Amplified Polymorphic DNA (RAPD) Markers. Plant Syst. Evol. 2004, 247, 165–175. [Google Scholar] [CrossRef]
  25. Henderson, S.A.; Billotte, N.; Pintaud, J.C. Genetic Isolation of Cape Verde Island Phoenix Atlantica (Arecaceae) Revealed by Microsatellite Markers. Conserv. Genet. 2006, 7, 213–223. [Google Scholar] [CrossRef]
  26. Ballardini, M.; Mercuri, A.; Littardi, C.; Abbas, S.; Couderc, M.; Ludeña, B.; Pintaud, J.C. The Chloroplast DNA Locus PsbZ-TrnfM as a Potential Barcode Marker in Phoenix L. (Arecaceae). Zookeys 2013, 365, 71–82. [Google Scholar] [CrossRef] [PubMed]
  27. Carreño, E. Diversidad Genética en especies del Género Phoenix. Ph.D. Dissertation, Universidad Miguel Hernández, Orihuela, Spain, 2017. [Google Scholar]
  28. Carreño, E.; Rivera, D.; Obón, C.; Alcaraz, F.; Johnson, D.; Bartual, J. What Are Candits? Study of a Date Palm Landrace in Spain Belonging to the Western Cluster of Phoenix dactylifera L. Genet. Resour. Crop Evol. 2021, 68, 135–149. [Google Scholar] [CrossRef]
  29. Ludeña, B.; Chabrillange, N.; Aberlenc-Bertossi, F.; Adam, H.; Tregear, J.W.; Pintaud, J.C. Phylogenetic Utility of the Nuclear Genes AGAMOUS 1 and PHYTOCHROME B in Palms (Arecaceae): An Example within Bactridinae. Ann. Bot. 2011, 108, 1433–1444. [Google Scholar] [CrossRef] [PubMed]
  30. Al-Dous, E.K.; George, B.; Al-Mahmoud, M.E.; Al-Jaber, M.Y.; Wang, H.; Salameh, Y.M.; Al-Azwani, E.K.; Chaluvadi, S.; Pontaroli, A.C.; Debarry, J.; et al. De Novo Genome Sequencing and Comparative Genomics of Date Palm (Phoenix dactylifera). Nat. Biotechnol. 2011, 29, 521–527. [Google Scholar] [CrossRef] [PubMed]
  31. Peakall, R.; Smouse, P.E. GenALEx 6.5: Genetic Analysis in Excel. Population Genetic Software for Teaching and Research-an Update. Bioinformatics 2012, 28, 2537–2539. [Google Scholar] [CrossRef]
  32. Pritchard, J.K.; Stephens, M.; Donnelly, P. Inference of Population Structure Using Multilocus Genotype Data. Genetics 2000, 155, 945–959. [Google Scholar] [CrossRef]
  33. Falush, D.; Stephens, M.; Pritchard, J.K. Inference of Population Structure Using Multilocus Genotype Data: Linked Loci and Correlated Allele Frequencies. Genetics 2003, 164, 1567–1587. [Google Scholar] [CrossRef]
  34. Anderson, E.C.; Thompson, E.A. A Model-Based Method for Identifying Species Hybrids Using Multilocus Genetic Data. Genetics 2002, 160, 1217–1229. [Google Scholar] [CrossRef] [PubMed]
  35. Vähä, J.P.; Primmer, C.R. Efficiency of Model-Based Bayesian Methods for Detecting Hybrid Individuals under Different Hybridization Scenarios and with Different Numbers of Loci. Mol. Ecol. 2006, 15, 63–72. [Google Scholar] [CrossRef] [PubMed]
  36. Arnold, M.L. Natural Hybridization as an Evolutionary Process. Annu. Rev. Ecol. Syst. 1992, 23, 237–261. [Google Scholar] [CrossRef]
  37. Tsitrone, A.; Kirkpatrick, M.; Levin, D.A. A Model for Chloroplast Capture. Evolution 2003, 57, 1776–1782. [Google Scholar] [CrossRef] [PubMed]
  38. Baldwin, E.; McNair, M.; Leebens-Mack, J. Rampant Chloroplast Capture in Sarracenia Revealed by Plastome Phylogeny. Front. Plant Sci. 2023, 14, 1237749. [Google Scholar] [CrossRef] [PubMed]
  39. Morales, J.; Speciale, C.; Rodríguez-Rodríguez, A.; Henríquez-Valido, P.; Marrero-Salas, E.; Hernández-Marrero, J.C.; López, R.; Delgado-Darias, T.; Hagenblad, J.; Fregel, R.; et al. Agriculture and Crop Dispersal in the Western Periphery of the Old World: The Amazigh/Berber Settling of the Canary Islands (ca. 2nd–15th Centuries Ce). Veg. Hist. Archaeobot 2023, 1, 1–15. [Google Scholar] [CrossRef]
  40. Abbott, R.J.; Barton, N.H.; Good, J.M. Genomics of Hybridization and Its Evolutionary Consequences. Mol. Ecol. 2016, 25, 2325–2332. [Google Scholar] [CrossRef]
  41. Kunkel, G.; Kunkel, M.A. Flora y Vegetación de Gran Canaria, 2nd ed.; Cabildo de Gran Canaria: Las Palmas de Gran Canaria, Spain, 2015; Volume 1. [Google Scholar]
Figure 1. (a) PCA analyses based on allele frequency obtained from 19 nuclear SSR loci on Phoenix samples. The percentage of explained variance of two first components is given in parenthesis. The genetic diversity is represented by a gradient of colour, the PC1 separated Phoenix canariensis (in green) from Phoenix dactylifera (in orange). The PC2 discriminated Phoenix canariensis samples according to island origins, from Gran Canaria (in green) and Tenerife (in grey). (b) Loading plot of PCA representing the alleles that exceeds a threshold established by 75% of the values plotted, indicating the weight of these alleles on the multidimensional ordination of Phoenix simples. Only the alleles above a 0.05 threshold in the loadings axis are labelled.
Figure 1. (a) PCA analyses based on allele frequency obtained from 19 nuclear SSR loci on Phoenix samples. The percentage of explained variance of two first components is given in parenthesis. The genetic diversity is represented by a gradient of colour, the PC1 separated Phoenix canariensis (in green) from Phoenix dactylifera (in orange). The PC2 discriminated Phoenix canariensis samples according to island origins, from Gran Canaria (in green) and Tenerife (in grey). (b) Loading plot of PCA representing the alleles that exceeds a threshold established by 75% of the values plotted, indicating the weight of these alleles on the multidimensional ordination of Phoenix simples. Only the alleles above a 0.05 threshold in the loadings axis are labelled.
Diversity 16 00411 g001
Figure 2. Identification and classification of hybrids: (a) STRUCTURE results showing the proportion of the genome of every individual originating from each of the two inferred clusters; orange is cluster 1 (Phoenix canariensis), and blue is cluster 2 (P. dactylifera). The bars between the two vertical black lines show the introgressed individuals with Q-values between 0.1 and 0.9. (b) NEWHYBRIDS classification for every individual with posterior probabilities for a given class; Cl266; Cl242 and Cl254: chlorotype of P. canariensis, P. dactylifera occidental and P. dactylifera oriental, respectively.
Figure 2. Identification and classification of hybrids: (a) STRUCTURE results showing the proportion of the genome of every individual originating from each of the two inferred clusters; orange is cluster 1 (Phoenix canariensis), and blue is cluster 2 (P. dactylifera). The bars between the two vertical black lines show the introgressed individuals with Q-values between 0.1 and 0.9. (b) NEWHYBRIDS classification for every individual with posterior probabilities for a given class; Cl266; Cl242 and Cl254: chlorotype of P. canariensis, P. dactylifera occidental and P. dactylifera oriental, respectively.
Diversity 16 00411 g002
Table 1. List of traits considered for the morphological assignation of the Phoenix species that samples belonged to. Hybrids were assigned when the samples exhibited intermediate characteristics between both species of Phoenix. See Figure S1.
Table 1. List of traits considered for the morphological assignation of the Phoenix species that samples belonged to. Hybrids were assigned when the samples exhibited intermediate characteristics between both species of Phoenix. See Figure S1.
AttributesPhoenix canariensisPhoenix dactylifera
StemCompletely straight. Columnar appearanceUsually, crooked. Not excessively thick
Presence of offshootsAbsenceOccasional presence
Size of leavesLong, broad leaves with elongated, dense leafletsSlightly shorter leaves, leaflets not very dense
Colour of the leavesDark greenGlaucus
Spines Long and thick, sometimes yellowish Short and thin
Shape of the crownSpherical Open leaf crown
Table 2. Allele frequency of loci that characterize the Canarian nuclear genotype. PC: Phoenix canariensis; HB: hybrids; PD242: P. dactylifera Occidental; PD254: P. dactylifera Oriental. In bold, allele frequencies higher than 0.900.
Table 2. Allele frequency of loci that characterize the Canarian nuclear genotype. PC: Phoenix canariensis; HB: hybrids; PD242: P. dactylifera Occidental; PD254: P. dactylifera Oriental. In bold, allele frequencies higher than 0.900.
LocusAbbreviated NameAllelePCHBPD242
242
PD254
254
mPdCIR010P101280.4780.4550.0000.000
1290.2920.1880.0000.000
1310.0680.0340.0000.000
mPdCIR015P15A1190.9780.7950.0000.000
P15B1510.9800.8300.0000.000
mPdCIR016P161110.9700.8010.0000.000
mPdCIR035P351750.9940.8690.0000.000
mPdCIR063P631390.1790.1820.0000.000
1410.4820.4030.0000.000
1430.2270.1190.0000.000
mPdCIR085P851470.0020.0000.0000.000
1490.6940.6080.0000.000
1510.1250.0740.0000.000
1530.0280.0000.0000.000
1550.0480.0110.0000.000
1610.0330.0000.0000.000
PdAP3-ssr-F4F42290.9940.8470.0000.000
mPdIRD031P313430.0020.0000.0000.000
3550.0090.0060.0000.000
3600.0450.0340.0000.000
3630.0170.0110.0000.000
3660.6640.6080.0000.000
3680.1410.1420.0000.000
3720.0990.0060.0000.000
mPdIRD033P331950.9760.8580.0000.000
PdAG1-ssrPAG12200.9110.5740.0810.000
mPdIRD040P401990.1240.1420.0000.000
2050.3540.2500.0000.000
2110.5090.3580.0000.000
2170.0040.0000.0000.000
Table 3. Distribution of the Canarian nuclear genotype (CNG) in molecular assignation and morphologic type groups. No. CNG: number of alleles of Canarian nuclear genotype. NT: total number of individuals. Morphology group: PD: Phoenix dactylifera (both types, Oriental and Occidental). PC_GC: Phoenix canariensis from Gran Canaria. PC_TF: Phoenix canariensis from Tenerife Island. HB_GC: hybrids from Gran Canaria. HB_TF: hybrids from Tenerife. PC: Phoenix canariensis (both islands). HB: hybrids (both islands).
Table 3. Distribution of the Canarian nuclear genotype (CNG) in molecular assignation and morphologic type groups. No. CNG: number of alleles of Canarian nuclear genotype. NT: total number of individuals. Morphology group: PD: Phoenix dactylifera (both types, Oriental and Occidental). PC_GC: Phoenix canariensis from Gran Canaria. PC_TF: Phoenix canariensis from Tenerife Island. HB_GC: hybrids from Gran Canaria. HB_TF: hybrids from Tenerife. PC: Phoenix canariensis (both islands). HB: hybrids (both islands).
Molecular
Assignation
No. CNGNTChlorotypeMorphology Group% Morphotype Match
C266C242C254PDPC GCPC TFHB GCHB TFPCHBPD
Phoenix canariensis pure22–24275264110015183202185.114.90
Molecular hybrids9–21774727301819291148.151.91
Phoenix dactylifera pure0–2810443774002507.392.7
Total-4333118240741691025137
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Saro, I.; Rodríguez-Rodríguez, P.; Rivera, D.; Obón, C.; Aberlenc, F.; Díaz-Pérez, A.; Zehdi-Azouzi, S.; Curbelo, L.; Sosa, P.A. The Genetic Characterization of the Canarian Endemic Palm (Phoenix canariensis) by Simple Sequence Repeats and Chloroplast Markers: A Tool for the Molecular Traceability of Phoenix Hybridization. Diversity 2024, 16, 411. https://doi.org/10.3390/d16070411

AMA Style

Saro I, Rodríguez-Rodríguez P, Rivera D, Obón C, Aberlenc F, Díaz-Pérez A, Zehdi-Azouzi S, Curbelo L, Sosa PA. The Genetic Characterization of the Canarian Endemic Palm (Phoenix canariensis) by Simple Sequence Repeats and Chloroplast Markers: A Tool for the Molecular Traceability of Phoenix Hybridization. Diversity. 2024; 16(7):411. https://doi.org/10.3390/d16070411

Chicago/Turabian Style

Saro, Isabel, Priscila Rodríguez-Rodríguez, Diego Rivera, Concepción Obón, Fredérique Aberlenc, Antonio Díaz-Pérez, Salwa Zehdi-Azouzi, Leticia Curbelo, and Pedro A. Sosa. 2024. "The Genetic Characterization of the Canarian Endemic Palm (Phoenix canariensis) by Simple Sequence Repeats and Chloroplast Markers: A Tool for the Molecular Traceability of Phoenix Hybridization" Diversity 16, no. 7: 411. https://doi.org/10.3390/d16070411

APA Style

Saro, I., Rodríguez-Rodríguez, P., Rivera, D., Obón, C., Aberlenc, F., Díaz-Pérez, A., Zehdi-Azouzi, S., Curbelo, L., & Sosa, P. A. (2024). The Genetic Characterization of the Canarian Endemic Palm (Phoenix canariensis) by Simple Sequence Repeats and Chloroplast Markers: A Tool for the Molecular Traceability of Phoenix Hybridization. Diversity, 16(7), 411. https://doi.org/10.3390/d16070411

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop