Genetic variation in the endemic and endangered rosmarinus tomentosus huber-morath & maire (labiatae) using rapd markers

ABSTRACT _Rosmarinus tomentosus_ Huber-Morath & Maire, an endemic species of southern Spain, is critically endangered as a consequence of habitat destruction by anthropogenic activities.

Random amplified polymorphic DNA (RAPD) markers were used for initial evaluation of genetic variation in this species; among zones, among populations (within zones and independently of

zones), and among individuals (within populations and zones). The eight primers used to study five populations provided 126 reliable bands, of which 109 (86%) were polymorphic. Analyses of

the RAPD markers conducted using different approaches (cluster analysis, Shannon’s diversity measure and analysis of molecular variance-AMOVA) showed similar results concerning genetic

structure of populations of _R. tomentosus_. The three zones studied, and their corrresponding populations, were clearly defined through UPGMA analysis. AMOVA analyses indicated that 18.0%

of the total genetic variation was attributable to differences among zones. When total variance was partitioned without considering the zonal distribution of the populations, 34.1% was

attributable to differences among populations, and 65.9% to differences among individuals within populations. Similar results were obtained using Shannon’s information measure. Significant

genetic distances (Phist values) were obtained between any two populations. This strong genetic differentiation among zones, populations and subpopulations may be explained by geographical

distances and limited gene flow. This study will be of considerable help during subsequent formulations of conservation strategies for the endangered species _R. tomentosus_. SIMILAR CONTENT

BEING VIEWED BY OTHERS LOW GENETIC DIVERSITY INDICATING THE THREATENED STATUS OF _RHIZOPHORA APICULATA_ (RHIZOPHORACEAE) IN MALAYSIA: DECLINED EVOLUTION MEETS HABITAT DESTRUCTION Article

Open access 05 November 2020 GENETIC DIVERSITY AND CONSERVATION OF SIBERIAN APRICOT (_PRUNUS SIBIRICA_ L.) BASED ON MICROSATELLITE MARKERS Article Open access 11 July 2023 INTER SIMPLE

SEQUENCE REPEAT MARKERS TO ASSESS GENETIC DIVERSITY OF THE DESERT DATE (_BALANITES AEGYPTIACA_ DEL.) FOR SAHELIAN ECOSYSTEM RESTORATION Article Open access 11 September 2020 INTRODUCTION

_Rosmarinus tomentosus_ Huber-Morath & Maire is a Spanish endemic species found in a very small coastal region within the province of Granada and barely extending into the province of

Málaga (both provinces in southern Spain). It is a diploid (2_n_=24) canescent, chasmophyte shrub with abundant branching, and has stems and small leaves covered with thick, greyish-white

tomentum (Maire, 1940; Rosúa, 1981). This species is predominantly outcrossing, as are most species of _Labiatae_ (Owens & Ubera-Jiménez, 1992), with entomophilous pollination mainly by

Hymenoptera, such as bees and bumblebees (J. E. Hernández Bermejo, pers. obs.). _R. tomentosus_ plants grow on coastal cliffs at altitudes between 0 and 250 m a.s.l.; some of these cliffs

have formed inland, as a result of fluvial sedimentation from the Guadalfeo river. According to recent estimates, fewer than 40 000 individuals survive distributed in the three zones. (i)

The western zone is made up of two coastal populations, CG (the cliffs of Cerro Gordo– Cantarriján) has fewer than 20 000 plants and PM (Punta de la Mona), separated from the first location,

has between 400 and 500 plants. (ii) The central and interior zone has one population (G) of about 300 individuals that grow in the Guadalfeo gorge, a former estuary that is now 10 km

inland from the coast. (iii) The eastern zone, on the coast between localities of Calahonda–Punta Cerrón (C-PC) and Castell de Ferro– Punta Gorda (CF-PG), has no more than 15 000 plants.

These plants are well adapted to this hostile ecological niche, characterized by extreme sun exposure, scarce precipitation, strong winds and saline environment (Hernández Bermejo &

Clemente Muñoz, 1994). Although the species is well adapted in nature, it is categorized as ‘Endangered’ (Hernández Bermejo & Clemente Muñoz, 1994). This is principally a result of the

anthropogenic pressure in its distribution area (e.g. tourism, housing developments, pollution, widening of highways), leading to its rapid demographic and geographical decline. Moreover,

this species has serious reproductive limitations, both asexually (difficulties in rooting of cuttings) and sexually (scant flowering, low seed productivity, very low germination capacity

and prolonged latency periods) (J.E. Hernández Bermejo, pers. obs.). Therefore, very few young plants are found in the wild. Another risk for the species is its propensity to hybridize with

_R. officinalis_ L. The hybrid, _R. x mendizabalii_ Sagredo ex Rosúa (Rosúa, 1981, 1986), is especially frequent in the central and interior population in the Guadalfeo valley, although it

is much less frequent in the eastern and western coastal populations, where the hostile climatic conditions put _R. tomentosus_ at an advantage over _R. officinalis_ and the hybrids. For all

these reasons _R. tomentosus_ should be re-evaluated as a ‘Critically Endangered’ taxon (IUCN, 1994). Different DNA fingerprinting methods are useful tools with a wide range of applications

in plant population studies, such as detection of genetic variation within and among populations, analysis of breeding systems, and studies of ecogeographical variation (Weising et al.,

1995). The development of random amplified polymorphic DNA (RAPD) markers (Williams et al., 1990) has proved useful for initial assessment of genetic variation in population and ecological

contexts (Hadrys et al., 1992; Gillies et al., 1997; Bussell, 1999). The RAPD procedure often yields large numbers of markers, it is technically simple to perform and does not need either

large amounts of DNA or previous information of DNA sequences, which are often unknown in wild plants. The utilization of these markers in genetic diversity and population structure studies

of plant species has increased in recent years (e.g. reviewed in Bussell, 1999), and especially for endemic, rare and/or endangered species, e.g. _Grevillea scapigera_ (Rossetto et al.,

1995), _Erodium paularense_ (Martín et al., 1997), _Caesalpinia echinata_ (Cardoso et al., 1998), and _Cerastium fischerianum_ var. _molle_ (Maki & Horie, 1999). Analysis of population

genetic structure using RAPD markers can be carried out using different approaches. Shannon’s diversity measure (King & Schaal, 1989) and Analysis of Molecular Variance (AMOVA; Excoffier

et al., 1992) overcome some of the difficulties of the other methods, such as those given by Clark & Lanigan (1993) and Lynch & Milligan (1994), which have limitations in some of

their assumptions (Bussell, 1999). The Shannon’s information measure has been applied for obtaining a partition of genetic variability within and between populations using RAPD markers in

several plant species (Martín et al., 1997; Wolff et al., 1997; Cardoso et al., 1998; Bussell, 1999). The AMOVA program was designed for haploid data but it has been used successfully for

RAPD data (Huff et al., 1993; Rossetto et al., 1995; Gillies et al., 1997; Martín et al., 1997; Cardoso et al., 1998). These reports have shown the strength of AMOVA as an analytical tool

that serves to discern regional and population differences. Accurate estimates of diversity are of great importance for formulating appropriate conservation and management strategies of the

genetic resources of a species. Keeping this in view, in this paper we have assessed the genetic variation among zones, within and among populations of the endemic and endangered species

_Rosmarinus tomentosus_ using RAPD markers. MATERIALS AND METHODS PLANT MATERIAL One hundred and ten individuals of _Rosmarinus tomentosus_ (22 individuals in each of the five populations

studied) were labelled in the field, and leaf material was collected for assessment of genetic variability. There were two coastal populations, which represented the eastern limit of the

distribution area of the species: Castell de Ferro–Punta Gorda (CF-PG) and Calahonda–Punta del Cerrón (C-PC). Two other coastal populations, Punta de la Mona (PM) and Cerro Gordo (CG),

represented the western limit. The last population (G) was located in a central zone, more interior towards the land where the Guadalfeo river reaches its final course, next to the edge of

Motril (Granada) (Fig. 1, Table 1). In the populations CF-PG and G, the 22 individuals were sampled as two subpopulations (CF-PG I and II; GI and II) and were separated by a distance of 200

m and 500 m, respectively. Individuals were sampled randomly and separated by at least 5–10 m. DNA EXTRACTION, AMPLIFICATION AND ELECTROPHORESIS DNA was extracted from frozen young leaves of

single individuals following the protocol of Torres et al. (1993). Extracted DNA was quantified using a spectrophotometer, and also by comparing band intensities with known standards of

lambda DNA on 1% (w/v) agarose gels. Working solutions of DNA (3 ng μL–1) were made using sterile, double-distilled water. In a preliminary study, 20 decamer primers of arbitrary sequence

(‘Kit O’ provided by Operon Technologies Inc., Alameda, CA) were tested for PCR amplification. Eight of these RAPD primers were chosen (based on number and polymorphism of amplified

fragments) to analyse the genetic variability of the samples (Table 2). DNA amplifications were performed in a 20-μL reaction volume containing about 30 ng template DNA, 0.2 μM of a single

primer, 200 μM of each dNTP and 1.5 units of BIOTAQTM DNA polymerase in the buffer provided by the enzyme manufacturers (Bioprobe systems). Amplification was performed in a PTC-100 thermal

cycler (MJ Research, Inc.) with a heated lid, and programmed for an initial step of 2 min at 94°C, followed by 35 cycles of 30 s at 94°C, 1 min at 37°C and 2 min at 72°C, and finally a 5

minute extension at 72°C. The amplification products were separated on 1.5% (w/v) agarose gels, in 1× TBE buffer, at 3 V cm–1 for 4 h, and then stained with ethidium bromide (0.5 μg mL–1)

for 12 min. Gels with amplification fragments were visualized and photographed under UV light. Molecular weights were estimated by reference to a 100-bp ladder (Pharmacia). All the

amplifications were repeated at least twice and only reproducible bands were considered for analysis. DATA ANALYSIS The RAPD bands were scored as 1 (present) or 0 (absent). A matrix of

similarities between every pair of individuals was created using the Dice’s similarity index (Dice, 1945), _S_=2_n__xy_/(_n__x_ + _n__y_), where _n__x_ and _n__y_ are the total number of

fragments analysed in individuals _x_ and _y_, respectively, and _n__xy_ is the number of fragments shared by the two individuals. The similarity matrix was employed to construct a UPGMA

dendrogram, using the SAHN-clustering and TREE programs from NTSYS-pc, v. 1.6 (Rohlf, 1992). A cophenetic matrix was derived from the similarity matrix to test goodness-of-fit of the

clusters, by comparing the two matrices using the Mantel matrix correspondence test (Mantel, 1967) in the MXCOMP program of the NTSYS-pc package. Genetic diversity was estimated using

Shannon’s information measure (King & Schaal, 1989), _H_o= − ∑_p__i_ log2 _p__i_, where _p__i_ is the frequency of a given RAPD fragment. _H_o was calculated at two levels: the average

diversity within populations (_H_pop), and the total diversity (_H_sp). Then the proportion of diversity within populations was estimated as _H_pop_/H_sp, and the proportion of diversity

among populations as (_H_sp − _H_pop)/_H_sp. In addition, an analysis of molecular variance (AMOVA; Excoffier et al., 1992) was applied to estimate variance components for RAPD phenotypes,

partitioning the variation among zones, among populations (within zones and independently of zones), and among individuals (within populations and within zones). Separate analyses were

carried out within eastern and western zones, and within CF-PG and G populations. For these analyses, the distance metric used was _D_=100(1 − _S_) for all pairs of individuals, where _S_ is

Dice’s estimator of similarity (Dice, 1945). Levels of significance of variance component estimates were computed by nonparametric permutational procedures (i.e. 2000 random permutations).

Homogeneity of molecular variance among zones, populations and subpopulations were tested with Bartlett’s statistics. Pairwise Phist distances (analogous to _F_-statistics at the molecular

level; Excoffier et al., 1992) were calculated among populations, and their significance were also tested by a permutation procedure, using 2000 permutations. All analyses were carried out

using WINAMOVA v. 1.55, provided by Laurent Excoffier (Genetics and Biometry Laboratory, University of Geneva, Switzerland). A Mantel matrix correspondence test (Mantel, 1967) was used to

analyse correlation between genetic (Phist values) and geographical distances among populations. RESULTS RAPD PROFILE Of the 126 reliable bands analysed in the amplifications with the eight

chosen primers, 109 bands (86%) were polymorphic in the five populations (Table 2). Considering each population separately, 48 polymorphic markers (58%) were found in population CF-PG, 72

(72%) in population C-PC, 71 (69%) in population G, 51 (54%) in population PM, and 47 (55%) in population CG. The eight primers revealed polymorphic bands, ranging between 61% (OPO-8) and

100% (OPO-1 and OPO-2); and the size range of the analysed fragments varied between 270 and 2050 bp (Table 2). Ten individual-specific bands were observed in four populations: five

individuals in population C-PC, two in population G, two in population PM, and one in population CG. Also, 14 population-specific bands were found. Nine of these bands were present in less

than 20% of the individuals of the corresponding populations, four were represented in 20–40% of the individuals of the populations, and one unique band (OPO15-1500) was present exclusively

in all the individuals of population C-PC (Fig. 2). Five additional bands that characterized the individuals of eastern and western zones were observed. Bands OPO8-700, OPO8-560 and

OPO16-1200 were present in 35% to 90% of the individuals of the eastern (CF-PG, C-PC) and central (G) zones, but were absent in all individuals from the western (PM, CG) zone. Conversely,

bands OPO7-1140 and OPO16-620 were present in 55% to 80% of the individuals of the western (PM, CG) and central (G) zones, but were absent in all individuals from the eastern (CF-PG, C-PC)

zone. Comparing all the populations, band OPO15–2050 was present in all the individuals of one population from each of the three zones: CF-PG (eastern zone), PM (western zone) and G (central

zone), and was absent from the populations C-PC and CG, belonging to eastern and western zones, respectively (Fig. 2). CLUSTER ANALYSIS The UPGMA dendrogram obtained using the 126 RAPD

markers scored in the 110 individuals showed five well-defined groups, which correspond to the five populations. Furthermore, population clusters corresponding to eastern zone, central zone

and the western zone were also observed in the phenogram (Fig. 3). The 22 individuals from each of populations CF-PG, PM and CG formed distinct clusters. In population C-PC, 20 out of 22

individuals formed a group, whilst the remaining two individuals were outside of this group but close to it. Similarly, in the Guadalfeo population two individuals were not in the same

cluster as the other individuals from the population. However, these two individuals were outside of the principal cluster of the remaining 108 individuals of _R. tomentosus_ (Fig. 3). Also,

clustering of individuals within subpopulations was observed in populations CF-PG and G. In population CF-PG, 64% (7/11) of individuals of subpopulation I and 55% (6/11) of subpopulation II

formed two clusters, respectively. Similarly, 75% (9/12) of individuals of subpopulation G I and 50% (5/10) of G II formed two groups (Fig. 3). The Mantel test revealed a significant

cophenetic correlation (_r_=0.76; _P_=0.001); thus, the dendrogram provides a good fit to the Dice’s similarity matrix. The information from the cluster analysis was complemented by an

analysis of variance of the coefficients of the similarity matrix calculated for each population. No significant differences were found between the mean similarity values for pairwise

comparisons of individuals within the five populations. GENETIC DIVERSITY ESTIMATES Genetic diversity measures within populations (_H_o) for each primer and mean values with respect to all

primers were calculated (Table 3). The highest genetic diversity values were obtained with primer OPO-1, in populations C-PC, G and CG, whilst the lowest were obtained with primers OPO-7 and

OPO-8; the latter did not detect diversity in population CG. The highest and the lowest mean genetic diversity values were obtained in the populations G (2.915) and CF-PG (1.812),

respectively. The mean diversity within the five populations of _R. tomentosus_ was _H_pop=2.335, and the total diversity was _H_sp=3.588 (Table 3). The proportion of average diversity

within populations (_H_pop_/H_sp) was 65.9%, and ranged between 52.6% (primer OPO-15) and 78.0% (primer OPO-2). Of the total diversity, 34.1% was among the populations (Table 3); all primers

detected more variability within than among populations. AMOVA PARTITION Highly significant (_P_<0.0005) genetic differences were detected in the analysis of molecular variance (AMOVA)

among zones, among populations (within zones and independently of zones), and among individuals (within populations and zones) (Table 4). Of the total molecular variance, 18.0% was

attributable to zonal divergence, 18.5% to populational differences within zones, and 63.5% to individual differences within populations. When the total variance was partitioned without

considering the zonal distribution of the populations, 34.1% was attributable to population divergence and 65.9% to individual differences within populations. Likewise, the partition

considering only the zonal distribution of the individuals resulted in 29.2% of the total variance due to diversity among zones and 70.8% to individual differences within zones (Table 4).

Bartlett’s test for homogeneity of the RAPD variance among zones did not indicate significant differences (_P_=0.092). On the other hand, the same test for homogeneity of variance for all

pairwise population comparisons were highly significant (_P_<0.0005). Similar variance values were found among populations for eastern and western zones (22.4% and 24.9%, respectively)

and among individuals within populations (77.6% and 75.1%, respectively) (Table 4). Bartlett’s test for homogeneity of the molecular variance indicated a greater significance between the two

eastern populations (CF-PG and C-PC; _P_<0.0005) compared to the two western populations (PM and CG; _P_<0.05). Variance values among subpopulations of populations CF-PG and G were

low (approx. 10%) but genetic differences were highly significant (_P_<0.0005) (Table 4). Bartlett’s test for homogeneity of variance indicated highly significant differences among

subpopulations of population G (_P_ < 0.0005), although the test was not significant in population CF-PG (_P_=0.034). Genetic distances among populations of _R. tomentosus_ obtained from

AMOVA (Phist values between pairs of populations) revealed a larger separation among populations from different zones than a pair of populations located in the same zone (Table 5). All

distances between pairs of populations were significantly different from zero. The Mantel matrix correspondence test indicated a significant correlation between the matrix of genetic

distances among populations with the corresponding matrix of geographical distances (_r_=0.926; _P_=0.038). DISCUSSION Analyses of the RAPD markers using different approaches (cluster

analysis, Shannon’s diversity measure and AMOVA) demonstrated similar interpretations of the genetic structure of populations of _R. tomentosus_. The present study revealed a high level of

polymorphism (86%) in the five populations of _R. tomentosus;_ considering each population separately, this value ranged between 54% and 72%. The percentage of polymorphic bands within

populations of _R. tomentosus_ is relevant for a rare endemic species, and is comparatively higher than in the other allogamous rare endemic species, such as _Allium aaseae_ (40–63%; Smith

& Pham, 1996) and _Erodium paularense_ (44–51%; Martín et al., 1997). Populations C-PC and G showed the highest number of polymorphic markers (about 70%), and, hence, the highest

estimates of genetic diversity (Table 3). The dendrogram (Fig. 3) shows that two of the 22 analysed individuals from each of these two populations were displaced from their corresponding

population clusters, and may be one of the reasons for these populations showing a greater diversity. The separation of the two variant individuals in the Guadalfeo valley population (G)

from the principal cluster could be the result of hybridization between _R. tomentosus_ and _R. officinalis_ (Rosúa, 1981, 1986), which is especially frequent in this interior central zone.

The AMOVA analysis shows that 18.0% of the total variance is attributable to zonal differences, and indicates significant genetic differentiation among zones. However, Bartlett’s test point

out no significant differences for the homogeneity of the variance among zones; hence, the genetic divergence was similar in the three zones. The results of the partition of variation among

populations and among individuals within populations were identical when comparing Shannon’s information measure and AMOVA (Tables 3 and 4. A close association exists between breeding system

and the distribution of genetic diversity (Schoen & Brown, 1991; Hamrick & Godt, 1996). Inbreeding species are generally characterized by high levels of genetic differentiation

among populations, whilst outbreeding plants tend to retain considerable variability within populations (Hamrick, 1990; Hamrick & Godt, 1996). Genetic differentiation estimates among

populations for outbreeding species when RAPD markers have been analysed using Shannon’s index ranged between 15% and 38%, and using the analysis of molecular variance this value was usually

<40%. For inbred species, interpopulation genetic variation was usually >50% (reviewed in Bussell, 1999). The genetic variation among populations of _R. tomentosus_ (34.1%) is towards

the higher limit and is in agreement with the outbreeding nature of the species. There are other outbreeding plant species that present similar values, e.g. _Alkanna orientalis_ (32–38%;

Wolff et al., 1997), _Cedrela odorata_ (35–37%; Gillies et al., 1997), and slightly higher values, as in the other endangered plant, _Caesalpinia echinata_ (41%; Cardoso et al., 1998). The

analysis of molecular variance shows highly significant differences among populations, and Bartlett’s test also indicates highly significant heterogeneity of the variance among populations.

Also, all Phist distances between pairs of _R. tomentosus_ populations are significantly different from zero; thus, the five populations may be considered genetically to be significantly

different from each other. These results indicate high genetic isolation among populations. However, a significant loss of diversity through genetic drift has not been detected in the

present study, probably because the size of the populations is not critically low as shown by Ellstrand & Elam (1993). The significant correlation between the matrix of genetic distances

(Phist values) and the matrix of geographical distances among populations (_r_=0.926; _P_=0.038) indicates that geographical isolation is an important factor influencing the genetic

structure of _R. tomentosus_ populations. Similar observations were made in _Caesalpina echinata_, which also has fragmented populations (Cardoso et al., 1998). Both the eastern and western

regions show similar levels of variation among populations (about 25%) and the genetic differentiation was highly significant. Moreover, Bartlett’s test for the homogeneity of variance

indicates a significantly greater genetic isolation between the two eastern populations (CF-PG and C-PC; _P_ < 0.0005) than between the two western populations (PM and CG; _P_ < 0.05).

This significantly higher isolation between the two eastern populations may be related to the increase of the climatic aridity from west to east in this coastal region, which may

consequently slow pollinator activity. The partitioning of variance among subpopulations of CF-PG and G (approx. 10%) are significant. Nevertheless, Bartlett’s test for homogeneity of

variance indicates a significant isolation between the two subpopulations of G, and not for the two subpopulations of CF-PG. In other words, it appears that a separation of 200 m (in

subpopulations of CF-PG) between groups of _R. tomentosus_ individuals does not impede gene flow, although a separation of about 500 m (in subpopulations of G) produces significant

isolation. The UPGMA dendrogram (Fig. 3) and the AMOVA analysis (Table 4) show a clear genetic differentiation among zones, among populations (within zones and independently of zones), and a

considerable divergence between subpopulations. The movement of genes among populations has a significant influence on the distribution of genetic variation. Species with restricted gene

flow should exhibit more genetic differentiation than species with widely dispersed pollen and seeds (Hamrick, 1990). The strong genetic differentiation in _R. tomentosus_ would suggest that

gene flow among different groups of individuals is very low. Genetic differentiation among populations attributed to a limited or absent interpopulation gene flow has been reported for

other endemic or endangered plant species (Brauner et al., 1992; Travis et al., 1996; Martín et al., 1997; Cardoso et al., 1998). Although the majority of _Labiatae_ species favour

outcrossing, their lack of self-incompatibility forces them to activate mechanisms that help them avoid excessive self-pollination (Owens & Ubera-Jiménez, 1992). _R. tomentosus_ staggers

its flowering throughout the year, resulting in just a few open flowers on any given plant at one time (Hernández Bermejo & Clemente Muñoz, 1994). The hymenopteran visitors therefore

have to go from plant to plant searching for the few _R. tomentosus_ flowers that are open; this is especially true in spring and summer. Thus, it is possible that staggering flowering may

slightly increase gene flow in _R. tomentosus_, but it is still quite low as indicated by this analysis using RAPD markers. Moreover, _R. tomentosus_ is a species with mainly

gravity-dispersed seeds, and most seedlings establish close to the mother plant (J.E. Hernández Bermejo, pers. obs.). The limited seed dispersal also contributes to reduced gene flow,

increasing genetic variation among different groups of individuals (Hamrick, 1990). In short, the limited gene flow among different populations and subpopulations of _R. tomentosus_,

together with the geographical separation, may explain the high level of isolation and the high genetic differentiation observed among populations and even among subpopulations from the

three zones studied. The use of RAPD markers to characterize and determine the genetic structure of the populations has proved to be a useful tool in the study of _Rosmarinus tomentosus_.

Considering the fragmented distribution of this endemic and endangered species, and the high genetic differentiation among populations as detected in the present study, we suggest that _in

situ_ conservation measures should include all the existing clusters of individuals (from subpopulation level to zonal level). The _ex situ_ conservation may be carried out by collecting

seeds from each of the zones, populations, and include different subpopulations for long-term storage in germplasm banks. REFERENCES * Brauner, S., Crawford, D. J. and Stuessy, T. F. (1992).

Ribosomal DNA and RAPD variation in the rare plant family Lactoridaceae. _Am J Bot_, 79: 1436–1439. Article  CAS  Google Scholar  * Bussell, J. D. (1999). The distribution of random

amplified polymorphic DNA (RAPD) diversity amongst populations of _Isotoma petraea_ (Lobeliaceae). _Mol Ecol_, 8: 775–789, 10.1046/j.1365-294X.1999.00627.x. Article  CAS  Google Scholar  *

Cardoso, M. A., Provan, J., Powell, W., Ferreiras, P. C. G. and de Oliveira, D. E. (1998). High genetic differentiation among remnant populations of the endangered _Caesalpinia echinata_

Lam. (Leguminosae-Caesalpinioideae). _Mol Ecol_, 7: 601–608. Article  Google Scholar  * Clark, A. G. and Lanigan, C. M. S. (1993). Prospects for estimating nucleotide divergence with RAPDs.

_Mol Biol Evol_, 10: 1096–1111. CAS  PubMed  Google Scholar  * Dice, L. R. (1945). Measures of the amount of ecologic association between species. _Ecology_, 26: 297–302. Article  Google

Scholar  * Ellstrand, N. C. and Elam, D. R. (1993). Population genetic consequences of small population size: implications for plant conservation. _Ann Rev Ecol Syst_, 24: 217–242. Article 

Google Scholar  * Excoffier, N. C., Smouse, P. E. and Quattro, J. M. (1992). Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human

mitochondrial DNA restriction data. _Genetics_, 131: 479–491. CAS  PubMed  PubMed Central  Google Scholar  * Gillies, A. C. M., Cornelius, J. P., Newton, A. C., Navarro, C., Hernández, M.

and Wilson, J. (1997). Genetic variation in Costa Rican populations of the tropical timber species _Cedrela odorata_ L., assessed using RAPDs. _Mol Ecol_, 6: 1133–1145. Article  CAS  Google

Scholar  * Hadrys, H., Balick, M. and Schierwater, B. (1992). Application of random amplified polymorphic DNA (RAPD) in molecular ecology. _Mol Ecol_, 1: 55–63. Article  CAS  PubMed  Google

Scholar  * Hamrick, J. L. (1990). Isozymes and the analysis of genetic structure in plant populations. In: Soltis, D. E. and Soltis, P. S. (eds) _Isozymes in Plant Biology_. pp. 87–105.

Chapman & Hall, London. Google Scholar  * Hamrick, J. L. and Godt, M. J. W. (1996). Effects of the life history traits on genetic diversity in plant species. _Phil Trans R Soc B_, 351:

1291–1298. Article  Google Scholar  * Hernández Bermejo, J. E. and Clemente Muñoz, M. (1994). _Protección de la Flora En Andalucía_. Consejería de Cultura y Medio Ambiente, Junta de

Andalucía, Seville. * Huff, D. R., Peakall, R. and Smouse, P. E. (1993). RAPD variation within and among natural populations of outcrossing buffalograss [_Buchloë dactyloides_ (Nutt.)

Engelm.]. _Theor Appl Genet_, 86: 927–934. Article  CAS  PubMed  Google Scholar  * IUCN (1994). _IUCN Red List Categories_. IUCN, Gland, Switzerland. * King, L. M. and Schaal, B. A. (1989).

Ribosomal DNA variation and distribution in _Rudbeckia missouriensis_. _Evolution_, 43: 1117–1119. Article  PubMed  Google Scholar  * Lynch, M. and Milligan, B. G. (1994). Analysis of

population genetic structure with RAPD markers. _Mol Ecol_, 3: 91–99. Article  CAS  PubMed  Google Scholar  * Maire, R. (1940). Sur un Romarin nouveau de l’Espagne meridionale. _Bull Soc

Hist Nat Afr Nord_, 31: 79–79. Google Scholar  * Maki, M. and Horie, S. (1999). Random amplified polymorphic DNA (RAPD) markers reveal less genetic variation in the endangered plant

_Cerastium fischerianum_ var. _molle_ than in the widespread conspecific _C. fischerianum_ var. _fischerianum_ (Caryophyllaceae). _Mol Ecol_, 8: 145–150, 10.1046/j.1365-294X.1999.00517.x.

Article  CAS  Google Scholar  * Mantel, N. (1967). The detection of disease clustering and a generalized regression approach. _Cancer Res_, 27: 209–220. CAS  PubMed  Google Scholar  *

Martín, C., González-Benito, M. E. and Iriondo, J. M. (1997). Genetic diversity within and among populations of a threatened species: _Erodium paularense_ Fern. Gonz. & Izco. _Mol Ecol_,

6: 813–820. Article  Google Scholar  * Owens, S. J. and Ubera-Jiménez, J. L. (1992). Breeding systems in _Labiatae_. In: Harley, R. M. and Reynolds, T. (eds) _Advances in Labiatae Science_.

pp. 257–280. Royal Botanic Gardens, Kew. Google Scholar  * Rohlf, F. J. (1992). NTSYS-PC: _Numerical taxonomy and multivariate analysis system_ v. 1.60. Exeter Software, New York. *

Rossetto, M., Weaver, P. K. and Dixon, K. W. (1995). Use of RAPD analysis in devising conservation strategies for the rare and endangered _Grevillea scapigera_ (Proteaceae). _Mol Ecol_, 4:

321–329. Article  CAS  PubMed  Google Scholar  * Rosúa, J. L. (1981). El complejo _Rosmarinus eriocalyx–tomentosus_ en la Península Ibérica. _Anales Jard Bot Madrid_, 37: 587–595. Google

Scholar  * Rosúa, J. L. (1986). Contribución al estudio del género _Rosmarinus_ L. en el Mediterráneo occidental. _Lagascalia_, 14: 179–187. Google Scholar  * Schoen, D. J. and Brown, A. H.

D. (1991). Intraspecific variation in population gene diversity and effective population size correlates with the mating system in plants. _Proc Natl Acad Sci USA_, 88: 4494–4497. Article 

CAS  PubMed  PubMed Central  Google Scholar  * Smith, J. F. and Pham, T. V. (1996). Genetic diversity of the narrow endemic _Allium aaseae_ (Alliaceae). _Am J Bot_, 83: 717–726. Article 

Google Scholar  * Torres, A. M., Weeden, N. F. and Martín, A. (1993). Linkage among isozyme, RFLP and RAPD markers in _Vicia faba_. _Theor Appl Genet_, 85: 937–945. Article  CAS  PubMed 

Google Scholar  * Travis, S. E., Maschinski, J. and Keim, P. (1996). An analysis of genetic variation in _Astragalus cremnophylax_ var. _cremnophylax_ a critically endangered plant, using

AFLP markers. _Mol Ecol_, 5: 735–745. Article  CAS  PubMed  Google Scholar  * Weising, K., Nybom, H., Wolff, K. and Meyer, W. (1995). _DNA Fingerprinting in Plants and Fungi_. CRC Press,

Boca Raton, FL. * Williams, J. G. K., Kubelik, A. R., Livak, K. J., Rafalski, J. A. and Tingey, S. V. (1990). DNA polymorphism amplified by arbitrary primers are useful as genetic markers.

_Nucl Acids Res_, 18: 6531–6535. Article  CAS  PubMed  PubMed Central  Google Scholar  * Wolff, K., El-Akkad, S. and Abbott, R. J. (1997). Population substructure in _Alkanna orientalis_

(Boraginaceae) in the Sinai Desert, in relation to its pollinator behaviour. _Mol Ecol_, 6: 365–372. Article  Google Scholar  Download references ACKNOWLEDGEMENTS This work was supported by

the Spanish Government, CICYT project no. AMB96–1021-C02–01. We personally thank our colleagues M. Clemente-Muñoz, A. Escudero, A. Lora, J.M. Iriondo, J. Prados and J.L. Vivero for their

collaboration with the field work in the _R. tomentosus_ habitat. AUTHOR INFORMATION AUTHORS AND AFFILIATIONS * Departamento de Biología Vegetal, Escuela Técnica Superior de Ingenieros

Agrónomos, Universidad Politécnica de Madrid, Ciudad Universitaria, Madrid, 28040, Spain Juan Pedro Martín * Jardín Botánico de Córdoba (Departamento de Ciencias y Recursos Agrícolas y

Forestales), Universidad de Córdoba. Avda Linneo s/n, Córdoba, 14004, Spain J Esteban Hernández Bermejo Authors * Juan Pedro Martín View author publications You can also search for this

author inPubMed Google Scholar * J Esteban Hernández Bermejo View author publications You can also search for this author inPubMed Google Scholar CORRESPONDING AUTHOR Correspondence to Juan

Pedro Martín. RIGHTS AND PERMISSIONS Reprints and permissions ABOUT THIS ARTICLE CITE THIS ARTICLE Martín, J., Hernández Bermejo, J. Genetic variation in the endemic and endangered

_Rosmarinus tomentosus_ Huber-Morath & Maire (Labiatae) using RAPD markers. _Heredity_ 85, 434–443 (2000). https://doi.org/10.1046/j.1365-2540.2000.00774.x Download citation * Received:

03 March 2000 * Accepted: 02 August 2000 * Published: 01 November 2000 * Issue Date: 01 November 2000 * DOI: https://doi.org/10.1046/j.1365-2540.2000.00774.x SHARE THIS ARTICLE Anyone you

share the following link with will be able to read this content: Get shareable link Sorry, a shareable link is not currently available for this article. Copy to clipboard Provided by the

Springer Nature SharedIt content-sharing initiative KEYWORDS * conservation * endangered species * population genetic structure * RAPD * _Rosmarinus tomentosus_