Analyses of genetic population structure of two ecologically important mangrove tree species, Bruguiera gymnorrhiza and Kandelia obovata from different river basins of Iriomote Island of the Ryukyu Archipelago, Japan

We analyzed nuclear and chloroplast microsatellite makers to assess genetic diversity and examine genetic structure of two mangrove tree species, Bruguiera gymnorrhiza and Kandelia obovata recovered from nine major river basins of Iriomote Island of the Ryukyu Archipelago, Japan. The average number of alleles per nuclear locus per population was 2.6 in B. gymnorrhiza and 1.7 in K. obovata. Bayesian clustering analysis using InStruct identified two genetic clusters in B. gymnorrhiza and three clusters in K. obovata. Chloroplast microsatellites revealed two dominant haplotypes from B. gymnorrhiza and three haplotypes from K. obovata. The overall result suggests low genetic diversity for both species. AMOVA for nuclear microsatellites showed 9.3?% of population variation in B. gymnorrhiza. Although genetic differentiation between several populations was significant in this species, F _ST suggested low to moderate level of differentiations between most of the populations. Distribution of genetic clusters and chloroplast haplotypes also suggested weak differentiations among B. gymnorrhiza populations. In K. obovata, population variation was, however, relatively high (27.8?%). The high differentiation between K. obovata populations was also suggested from the F _ST and genetic clusters from nuclear microsatellites, and chloroplast haplotypes. A significant correlation between chloroplast genetic distances and coastline distances as well as haplotype distributions for both species suggest that propagules from each species mostly disperse to the neighboring river basins. While significant F _IS and higher percentage of admixed clusters from nuclear microsatellites suggested inbreeding, continual gene exchange by propagule dispersal among populations especially among neighboring populations was suggested for both species from maternally inherited chloroplast microsatellites analyses.     

Chloroplast microsatellite markers for the mangrove tree species Bruguiera gymnorrhiza, Kandelia candel, and Rhizophora stylosa, and cross-amplification in other mangrove species

Chloroplast microsatellite (cpSSR) markers were developed for three ecologically and economically important tree species in the mangrove family, Rhizophoraceae: Bruguiera gymnorrhiza, Kandelia candel, and Rhizophora stylosa. Noncoding regions of chloroplast DNA (cpDNA) from each species were separately amplified using universal chloroplast primers. Six, two, and three polymorphic cpSSR loci in B. gymnorrhiza, K. candel, and R. stylosa, respectively, were developed from amplified noncoding cpDNA regions. Characterization of 216, 156, and 253 individuals of B. gymnorrhiza, K. candel, and R. stylosa, respectively, collected from different natural mangrove populations (B. gymnorrhiza, 9; K. candel, 7; R. stylosa, 9) on Iriomote Island in Japan showed that these loci provide cpSSR markers with polymorphisms ranging from two to four alleles per locus and gene diversity between 0.027 and 0.480. These cpSSR markers will be useful for analyzing the maternal lineage distributions and population genetic structures of the three species. Several of these markers may also be useful in similar studies of other mangrove species.     

Patterns of genetic variation in the Chinese endemic Psilopeganum sinense (Rutaceae) as revealed by nuclear microsatellites and chloroplast microsatellites

Psilopeganum sinense is a perennial herb endemic to Three-Gorges Reservoir Area (TGRA) in China. Genetic diversity of this endangered species was assessed by using 11 nuclear microsatellites and three chloroplast microsatellite (cpSSR) markers. A total of 8 haplotypes were identified in a survey of 212 individuals sampled from nine populations encompassing most of the natural range of the species. A low level of genetic diversity was detected: H_E = 0.301 for SSR and H_E = 0.28 for cpSSR. Populations were highly differentiated from one another: an AMOVA analysis that showed that 56.3% and 68.2% genetic variation resided between populations based on SSR and cpSSR analysis, respectively, and F_ST and F_STc (0.467 for SSR and 0.644 for cpSSR, respectively) were high. Significant differences were found between estimates of haplotypic differentiation calculated by using unordered alleles (G_STc = 0.857) and ordered alleles (N_STc = 0.728), which indicated the existence of phylogeographical structure in P. sinense. The indirect ratio of pollen flow/seed flow derived from estimates of haplotypic and nuclear DNA differentiation indicated that gene flow via pollen is less efficient than via seed. Two distinct evolutionary lineages (evolutionary significant units, ESUs) were recognized for P. sinense on the basis of both the PCoA and NCA analysis. Sampling strategies for conserving this endangered plant were discussed.     

Chloroplast microsatellite diversity of Clintonia udensis (Liliaceae) populations in East Asia

Nineteen populations of Clintonia udensis Trautv. & Mey. were examined to quantify genetic diversity and genetic structure by chloroplast DNA microsatellites (cpSSR). Significant cpSSR genetic diversity (PPB = 63.64%) was detected in this species. Tetraploid populations demonstrated approximately the same level of genetic diversity as the diploid ones. A significant differentiation, however, was found between tetraploids and diploids. Most of the sixteen chloroplast haplotypes were limited to a single population. The level of haplotype diversity was high (Hd = 0.915). AMOVA, PCA and Bayesian clustering analysis revealed that there were significant genetic differences among populations. Inter-population genetic distances among population sites correlated significantly with geographic distances. These results indicate that the mixed-mating – breeding system, limited gene flow, environmental stress, and historical factors may be the main factors causing geographical differentiation in the genetic structure of C. udensis.     

Geographic patterns of genetic variation in native pecans

A structured collection of 80 seedling pecan trees [Carya illinoinensis (Wangenh.) K. Koch], representing 19 putatively native pecan populations across the species range, was evaluated at three plastid and 14 nuclear microsatellite (simple sequence repeat, SSR) loci. Data were analyzed using a priori population designations and also within a Bayesian framework, in which individuals were assigned to clusters regardless of population of origin. Population genetic analyses using a priori populations, clusters based on chloroplast microsatellite data (cpSSR), and clusters based on nuclear microsatellite data (nSSR) yielded consistent results. For all groupings, cpSSR variation exhibited more geographic structure than the nSSR data. Furthermore, cpSSR microsatellite diversity decreased with increasing latitude, but this pattern was not observed with the nuclear data. Contrasting patterns in plastid and nuclear genetic diversity demonstrate unique aspects of postglacial recolonization reflected in the movement of seeds versus pollen. These data suggest that plastid SSRs are useful tools for identifying population structure in pecan and hold promise for ongoing efforts to identify and conserve representative germplasm in ex situ collections.     

Population genetic structure of Vitex negundo (Verbenaceae) in Three-Gorge Area of the Yangtze River: The riverine barrier to seed dispersal in plants

Previous studies have demonstrated that large rivers can influence inter- and intra-specific gene flow for many animals. The effects of large rivers on the genetics of plant populations have focused on either hydrochoric impacts of water current on gene flow or genetic differentiation among populations from different watersheds. Few studies have explicitly tested the barrier effects on plant gene flow across banks of large rivers, especially their relative effects on pollen and seed dispersals. The Yangtze River (Changjiang River), one of the major rivers of the world, provides an excellent model to evaluate the impacts of rivers on gene flow in plants. Using RAPD (random amplified polymorphic DNA) and cpDNA (chloroplast DNA) markers, we investigated the genetic structure of 10 populations of Vitex negundo in two regions of Three-Gorge Area along the Yangtze River. Each region contained two populations on the north bank, two on the south bank and one island population along the river. The analyses indicated low RAPD between banks, and similar or a little higher differentiation between populations within the same bank. In contrast, a large proportion of chloroplast polymorphism was ascribed to among-bank variation but much lower cpDNA differentiation was among populations within the same bank. These results indicate that the Yangtze River represents a general barrier to the dispersal of seeds but not to the movement of pollen in V. negundo. The cpDNA genetic distances or differentiations between the island populations and those on either bank of the river are intermediate to those between the banks across the river, implying that the islands in the Yangtze River may serve as a stepping-stone for seed dispersal. Our results suggest that large rivers may serve as a general barrier, not only for the movement of animals, but also for the dispersal of plants, which should be of great significance for the conservation of biodiversity around the rivers.     

Consequences of multiple mating‐system shifts for population and range‐wide genetic structure in a coastal dune plant

Evolutionary transitions from outcrossing to selfing can strongly affect the genetic diversity and structure of species at multiple spatial scales. We investigated the genetic consequences of mating‐system shifts in the North American, Pacific coast dune endemic plant Camissoniopsis cheiranthifolia (Onagraceae) by assaying variation at 13 nuclear (n) and six chloroplast (cp) microsatellite (SSR) loci for 38 populations across the species range. As predicted from the expected reduction in effective population size (N_e) caused by selfing, small‐flowered, predominantly selfing (SF) populations had much lower nSSR diversity (but not cpSSR) than large‐flowered, predominantly outcrossing (LF) populations. The reduction in nSSR diversity was greater than expected from the effects of selfing on N_e alone, but could not be accounted for by indirect effects of selfing on population density. Although selfing should reduce gene flow, SF populations were not more genetically differentiated than LF populations. We detected five clusters of nSSR genotypes and three groups of cpSSR haplotypes across the species range consisting of parapatric groups of populations that usually (but not always) differed in mating system, suggesting that selfing may often initiate ecogeographic isolation. However, lineage‐wide genetic variation was not lower for selfing clusters, failing to support the hypothesis that selection for reproductive assurance spurred the evolution of selfing in this species. Within three populations where LF and SF plants coexist, we detected genetic differentiation among diverged floral phenotypes suggesting that reproductive isolation (probably postzygotic) may help maintain the striking mating‐system differentiation observed across the range of this species.     

Genetic diversity,structure, and demography of Pandanus boninensis (Pandanaceae) with sea drifted seeds,endemic to the Ogasawara Islands of Japan: Comparison between young and old islands

Pandanus boninensis, endemic to the Ogasawara Islands, Japan, is distributed on both the older Bonin and younger Volcano Islands. In this study, we conducted population genetic analyses of P. boninensis on these islands to examine the population diversity and structure across old and young islands, to assess potential differences in population demography with island age, and to collect any evidence of migration between old and young islands. We found that the genetic diversity of expressed sequence tag (EST)–based microsatellite (SSR) markers, the nucleotide diversity of nuclear DNA sequences, and the haplotype diversity of chloroplast DNA on young islands were lower than those on old islands. Clustering analyses of EST‐SSR indicated that populations on old islands were strongly diverged from those on young islands. Approximate Bayesian computation analysis of EST‐SSR suggested that population expansion occurred on old islands while population reduction occurred on young islands. We also found evidence of migration among old islands (mostly from south to north), while it appears that there have been very few migration events between old and young islands. These differences could be due to the fact that young islands tend to be geographically isolated and support smaller populations that began a shorter time ago from limited founders. The P. boninensis populations on the Volcano Islands are interesting from an evolutionary perspective as they constitute a classic example of the early stages of progressive colonization on oceanic islands with small effective population sizes and low genetic diversity.     

Genetic structure of disjunct Argentinean populations of the subtropical tree Anadenanthera colubrina var. cebil (Fabaceae)

Anadenanthera colubrina var. cebil is a native South American tree species inhabiting seasonally dry tropical forests (SDTFs). Its current disjunct distribution presumably represents fragments of a historical much larger area of this forest type, which has also been highly impacted by human activities. In this way the hypothesis of this study is that the natural populations of A. colubrina var. cebil from Northern Argentina represent vestiges of ancient fragmentation, but they are additionally influenced by a certain degree of gene flow among them. We aimed to analyze the genetic structure of both nuclear and chloroplast DNA to evaluate the relative role of ancient and recent fragmentation on intraspecific diversity patterns. Sixty-nine individuals of four natural populations were analyzed using eight nuclear microsatellites (ncSSR) and four chloroplast microsatellite loci (cpSSR). The level and distribution of genetic variation were estimated by standard population genetic parameters and Neighbor Joining as well as Bayesian analyses. The eight ncSSR loci were highly polymorphic, while genetic diversity of cpSSRs was low. Nuclear SSRs displayed lower genetic differentiation among populations than cpSSR haplotypes (F _ST 0.11 and 0.95, respectively). However, high differentiation between phytogeographic provinces was observed in both genomes. The high genetic differentiation detected emphasizes the role of ancient fragmentation. However, the Paranaense province also shows the effects of recent fragmentation on genetic structure, whereas gene flow by pollen preserves the effects of genetic drift in the Yungas province.     

Spatial and temporal population genetic variation and structure of Nothotsuga longibracteata (Pinaceae), a relic conifer species endemic to subtropical China

Nothotsuga longibracteata, a relic and endangered conifer species endemic to subtropical China, was studied for examining the spatial-temporal population genetic variation and structure to understand the historical biogeographical processes underlying the present geographical distribution. Ten populations were sampled over the entire natural range of the species for spatial analysis, while three key populations with large population sizes and varied age structure were selected for temporal analyses using both nuclear microsatellites (nSSR) and chloroplast microsatellites (cpSSR). A recent bottleneck was detected in the natural populations of N. longibracteata. The spatial genetic analysis showed significant population genetic differentiation across its total geographical range. Notwithstanding, the temporal genetic analysis revealed that the level of genetic diversity between different age class subpopulations remained constant over time. Eleven refugia of the Last Glacial Maximum were identified, which deserve particular attention for conservation management.     

A range wide geographic pattern of genetic diversity and population structure of Castanea mollissima populations inferred from nuclear and chloroplast microsatellites

Understanding geographical pattern of genetic diversity and population structure is of great importance for formulating conservation and utilization strategies. In this study, we investigated the genetic diversity and population structure of 28 natural populations of Castanea mollissima in China using eight nuclear and six chloroplast microsatellite makers (nSSRs and cpSSRs). Populations from central China harbored the highest genetic diversity at both nSSR and cpSSR markers (nSSR: H _E?=?0.705; cpSSR: H?=?0.461). The standardized measure of genetic differentiation estimated as G′ _ST was 0.447 for nSSR and 0.803 for cpSSR, respectively. The G′_ST-based pollen to seed flow ratio is 3.043, indicating that pollen flow is not extensive among C. mollissima populations. No obvious population genetic structure by geographical locations was found by STRUCTURE analysis based on nSSR data, and similarly, no signal of phylogeographic structure was detected for cpSSR analysis. Five boundaries defining zones of maximum genetic differences within the network of the C. mollissima populations were found, and the locations of those barriers were consistent with those of four mountains, i.e., Daloushan Mountain, Dabashan Mountain, Wushan Mountain, and Qingliangfeng Mountain, indicating that those mountains might act as genetic barriers obstructing the genetic exchange among natural C. mollissima populations. These results provide valuable baseline data for conservation and utilization of this species.     

Dissimilar patterns of Pinus heldreichii Christ. populations in Bulgaria revealed by chloroplast microsatellites and terpenes analysis

In the present study we investigated the genetic structure and genetic diversity of Pinus heldreichii populations in Bulgaria using chloroplast microsatellite markers and terpene analysis. We were interested in addressing the following questions: (1) can population structuring in Bosnian pine be detected via chloroplast microsatellite markers; (2) are there differences in population differentiation as determined by terpenes and microsatellites; and (3) how are the patterns of size variant frequencies and geographical distances related. Four provenances were chosen throughout the species' range in Bulgaria. Following DNA extraction, chloroplast microsatellite (cpSSR) loci were surveyed using 6 primer pairs. Between 2 and 5 size variants were identified at each locus. A total of 16 size variants at the 6 loci were identified, 4 occurring at low frequencies. They were combined in 21 different haplotypes including 11 that were unique. AMOVA analysis revealed that 18.25% of the variation was found among populations, while 81.75% was expressed within populations. The cpSSR analysis divided Bosnian pine populations into two groups, the first represented by populations B, C and D located in the south and north-western part of the Pirin and Slavianka mountains, while the second group, represented by population A, is located in the north-eastern Pirin mountain. Terpene analysis revealed that on average, 59% of the monoterpene pool in P. heldreichii is accounted for by limonene (range 36–48%) followed by α-pinene (range 16–17%). The presence of two distinct groups (Pop-A, Pop-D and Pop-B, Pop-C) is more consistent with physical distances between populations. No significant correlation between genetic distance determined by chloroplast microsatellites analysis and chemotype distance determined by terpenes was observed. Our results suggest that the structural pattern of genetic diversity of cpDNA in Bosnian pine populations is the consequence of historical and biogeographical processes.     

The diversity of chloroplast microsatellite loci in Siberian fir (Abies sibirica Ledeb.) and two Far East fir species A. nephrolepis (Trautv.) Maxim. and A. sachalinensis Fr. Schmidt

The genetic variability in 29 populations of Abies sibirica, three of A. nephrolepis, and seven of A. sachalinensis was studied using SSR markers of chloroplast DNA. Stable amplification and polymorphic products were obtained using primer pairs Pt71936 and Pt30204 (with nine and forteen alleles, respectively) of 10 pairs. Totally, 70 haplotypes were found, 43 in A. sibirica, 49 in A. sachalinensis, and 31 in A. nephrolepis. The highest values of genetic diversity parameters were observed in A. sachalinensis, and the lowest in A. nephrolepis. The Siberian fir differs from Far East species by the uneven multimodal frequency distributions of allele length in both cpSSR loci, which is explained by the presence of few separated from each other dominating haplotypes. This fact indicates that A. sibirica and the Far East species have different demographic histories. In A. sibirica, the proportion of diversity between populations in the total genetic diversity, calculated taking into account the differences between haplotypes (R(ST)) was 8.34% and 4.42% without accounting for haplotypes differences (R(ST) > G(ST), P= 0.01). The pairwise G(ST) correlate significantly with geographic distances between the populations A. sibirica and with genetic distances D calculated from allozyme data. No such correlations were found with the R(ST) parameter. The results of cpSSR variability analysis strongly support the conclusions inferred from allozyme data: several geographic groups of comparatively genetically close populations are identified, which may be explained by the invasion of colonization of the present-day Siberian fir range.     

Plio-Pleistocene history and phylogeography of Acacia senegal in dry woodlands and savannahs of sub-Saharan tropical Africa: evidence of early colonisation and recent range expansion

Drylands are extensive across sub-Saharan Africa, socio-economically and ecologically important yet highly sensitive to environmental changes. Evolutionary history, as revealed by contemporary intraspecific genetic variation, can provide valuable insight into how species have responded to past environmental and population changes and guide strategies to promote resilience to future changes. The gum arabic tree (Acacia senegal) is an arid-adapted, morphologically diverse species native to the sub-Saharan drylands. We used variation in nuclear sequences (internal transcribed spacer (ITS)) and two types of chloroplast DNA (cpDNA) markers (PCR-RFLP, cpSSR) to study the phylogeography of the species with 293 individuals from 66 populations sampled across its natural range. cpDNA data showed high regional and rangewide haplotypic diversity (h_T(cpSSR)=0.903–0.948) and population differentiation (G_ST(RFLP)=0.700–0.782) with a phylogeographic pattern that indicated extensive historical gene flow via seed dispersal. Haplotypes were not restricted to any of the four varieties, but showed significant geographic structure (G_ST(cpSSR)=0.392; R_ST=0.673; R_ST>R_ST (permuted)), with the major division separating East and Southern Africa populations from those in West and Central Africa. Phylogenetic analysis of ITS data indicated a more recent origin for the clade including West and Central African haplotypes, suggesting range expansion in this region, possibly during the Holocene humid period. In conjunction with paleobotanical evidence, our data suggest dispersal to West Africa, and across to the Arabian Peninsula and Indian subcontinent, from source populations located in the East African region during climate oscillations of the Plio-Pleistocene.     

Structure of Pinus nigra Arn. populations in Bulgaria revealed by chloroplast microsatellites and terpenes analysis: Provenance tests

Pinus nigra is a forest and low elevation mountain species found around the Mediterranean Sea that has had its distribution reduced and fragmented by anthropogenic disturbance. Due to commercial interest it is currently being replanted, however, the genetic structure of populations is little known and current planting strategies could threaten its genetic diversity. In the present study we investigated the genetic structure and genetic diversity of P. nigra populations in Bulgaria using chloroplast microsatellite markers and terpene analysis. Nine provenances were chosen throughout the species' range in Bulgaria. Following DNA extraction, chloroplast microsatellite (cpSSR) loci were surveyed using three primer pairs. Between 5 and 9 size variants were identified at each locus. A total of 22 size variants at the 3 loci were identified, that were combined in 68 different haplotypes, of which 7 represent 39.81% of the genetic structure. AMOVA analysis revealed that 6.06% of the variation was found among populations, while 93.94% was expressed within populations. The cpSSR analysis divided European Black pine populations into four groups, the first represented by populations located the eastern Rhodopes, Sr. Gora and St. Planina mountains, while the second group is primarily located in the Phodopes and Slavianca mountains. The populations from Pirin and Osogovo mountains show different genetic patterns. Terpene analysis revealed that most of the monoterpene pool in P. nigra was accounted for by α-pinene followed by β-pinene. The presence of four distinct terpene groups is not consistent with physical distances between populations, and a similar non-significant correlation between genetic distance determined by chloroplast microsatellites analysis and chemotype distance (determined by terpenes) was observed. Our results suggest that the structural pattern of genetic diversity of cpDNA in European Black pine populations is the consequence of historical biogeographic processes.     

Genetic differentiation and delimitation between ecologically diverged Populus euphratica and P. pruinosa

Background

The fixed genetic differences between ecologically divergent species were found to change greatly depending on the markers examined. With such species it is difficult to differentiate between shared ancestral polymorphisms and past introgressions between the diverging species. In order to disentangle these possibilities and provide a further case for DNA barcoding of plants, we examine genetic differentiation between two ecologically divergent poplar species, Populus euphratica Oliver and P. pruinosa Schrenk using three different types of genetic marker.

Methodology/Principal Findings

We genotyped 290 individuals from 29 allopatric and sympatric populations, using chloroplast (cp) DNA, nuclear (nr) ITS sequences and eight simple sequence repeat (SSR) loci. Three major cpDNA haplotypes were widely shared between the two species and between-species cpDNA differentiation (F_CT) was very low, even lower than among single species populations. The average SSR F_CT values were higher. Bayesian clustering analysis of all loci allowed a clear delineation of the two species. Gene flow, determined by examining all SSR loci, was obvious but only slightly asymmetrical. However, the two species were almost fixed for two different nrITS genotypes that had the highest F_CT, although a few introgressed individuals were detected both in allopatric and sympatric populations.

Conclusions

The two species shared numerous ancestral polymorphisms at cpDNA and a few SSR loci. Both ITS and a combination of nuclear SSR data could be used to differentiate between the two species. Introgressions and gene flow were obvious between the two species either during or after their divergence. Our findings underscore the complex genetic differentiations between ecologically diverged species and highlight the importance of nuclear DNA (especially ITS) differentiation for delimiting closely related plant species.     

An extreme cytoplasmic bottleneck in the modern European cultivated potato (Solanum tuberosum) is not reflected in decreased levels of nuclear diversity

We have used the polymorphic chloroplast (cp) and nuclear simple sequence repeats (SSRs) to analyse levels of cytoplasmic and nuclear diversity in the gene pool of the European cultivated potato (Solanum tuberosum ssp. tuberosum). Primers designed from the complete chloroplast sequence of tobacco (Nicotiana tabacum) were used to amplify polymorphic products in a range of potato cultivars. Combining the data from seven polymorphic cpSSR loci gave 26 haplotypes, one of which (haplotype A) accounted for 151 out of the 178 individuals studied and corresponded to the T-type cytoplasm previously identified in cultivated potatoes using chloroplast restriction fragment length polymorphism analysis. Phylogenetic and diversity analyses of the relationships between cpSSR haplotypes confirmed much higher levels of cytoplasmic diversity outwith the T-type group. Diversity levels at eight nuclear SSR loci, however, were not significantly different between cytoplasmic groups, suggesting a severe maternal bottleneck in the evolution of the modern cultivated potato. These results highlight the importance in quantifying levels of cytoplasmic as well as nuclear diversity and confirm the need for a change in breeding practices to increase levels of non-T-type cytoplasm in the cultivated gene pool, thus helping reduce problems associated with pollen sterility. This may be facilitated by germplasm analysis using cpSSRs, which will allow efficient selection of diverse cytoplasm donors.     

Applications of Multiple Nuclear Genes to the Molecular Phylogeny,Population Genetics and Hybrid Identification in the Mangrove Genus Rhizophora

The genus Rhizophora is one of the most important components of mangrove forests. It is an ideal system for studying biogeography, molecular evolution, population genetics, hybridization and conservation genetics of mangroves. However, there are no sufficient molecular markers to address these topics. Here, we developed 77 pairs of nuclear gene primers, which showed successful PCR amplifications across all five Rhizophora species and sequencing in R. apiculata. Here, we present three tentative applications using a subset of the developed nuclear genes to (I) reconstruct the phylogeny, (II) examine the genetic structure and (III) identify natural hybridization in Rhizophora. Phylogenetic analyses support the hypothesis that Rhizophora had disappeared in the Atlantic-East Pacific (AEP) region and was re-colonized from the IWP region approximately 12.7 Mya. Population genetics analyses in four natural populations of R. apiculata in Hainan, China, revealed extremely low genetic diversity, strong population differentiation and extensive admixture, suggesting that the Pleistocene glaciations, particularly the last glacial maximum, greatly influenced the population dynamics of R. apiculata in Hainan. We also verified the hybrid status of a morphologically intermediate individual between R. apiculata and R. stylosa in Hainan. Based on the sequences of five nuclear genes and one chloroplast intergenic spacer, this individual is likely to be an F1 hybrid, with R. stylosa as its maternal parent. The nuclear gene markers developed in this study should be of great value for characterizing the hybridization and introgression patterns in other cases of this genus and testing the role of natural selection using population genomics approaches.     

Chloroplast microsatellites reveal that metallicolous populations of the Mediterranean shrub Cistus ladanifer L have multiple origins

Cistus ladanifer L. (Cistaceae) is a Mediterranean shrub covering different kinds of soils in the Western Mediterranean area. This species has colonised several metalliferous areas (serpentine outcrops as well as human-polluted sites) throughout its distribution range, and is therefore an interesting species to study the possible effects on genetic diversity and differentiation produced by the colonisation of areas polluted with heavy metals. The genetic structure of 33 natural populations distributed across its entire natural distribution range (Morocco, Portugal and Spain) and growing on either metalliferous or non-metalliferous soils was investigated using chloroplast microsatellites. Population genetic parameters were estimated and genetic groups were identified using Bayesian inference. In addition, we compared the genetic diversity and differentiation among metallicolous and non-metallicolous populations within each Bayesian-defined group. The cpSSR data suggested that metallicolous populations of Cistus ladanifer have arisen through multiple independent evolutionary origins within two different chloroplast lineages. Evidence that the soil type provoked genetic bottlenecks in metallicolous populations or genetic differentiation among metallicolous and non-metallicolous populations was not observed. Historical factors are the main cause of the present genetic structure of C. ladanifer. The nature of tolerance to heavy metals as a species-wide trait in this shrub is discussed.     

Complex population genetic structure in the endemic Canary Island pine revealed using chloroplast microsatellite markers

The Canary archipelago, located on the northwestern Atlantic coast of Africa, is comprised of seven islands aligned from east to west, plus seven minor islets. All the islands were formed by volcanic eruptions and their geological history is well documented providing a historical framework to study colonization events. The Canary Island pine (Pinus canariensis C. Sm.), nowadays restricted to the westernmost Canary Islands (Gran Canaria, Tenerife, La Gomera, La Palma and El Hierro), is considered an old (Lower Cretaceous) relic from an ancient Mediterranean evolutionary centre. Twenty seven chloroplast haplotypes were found in Canary Island pine but only one of them was common to all populations. The distribution of haplotypic variation in P. canariensis suggested the colonization of western Canary Islands from a single continental source located close to the Mediterranean Basin. Present-day populations of Canary Island pine retain levels of genetic diversity equivalent to those found in Mediterranean continental pine species, Pinus pinaster and Pinus halepensis. A hierarchical analysis of variance (AMOVA) showed high differentiation among populations within islands (approximately 19%) but no differentiation among islands. Simple differentiation models such as isolation by distance or stepping-stone colonization from older to younger islands were rejected based on product-moment correlations between pairwise genetic distances and both geographic distances and population-age divergences. However, the distribution of cpSSR diversity within the islands of Tenerife and Gran Canaria pointed towards the importance of the role played by regional Pliocene and Quaternary volcanic activity and long-distance gene flow in shaping the population genetic structure of the Canary Island pine. Therefore, conservation strategies at the population level are strongly recommended for this species.Communicated by D.B. NealeA. Gómez and S.C. González-Martínez as joint authors