Spatial genetic structure of Taxus baccata L. in the western Mediterranean Basin: past and present limits to gene movement over a broad geographic scale

Assessment of genetic conservation units of an endangered glacial relict insular species, Amentotaxus formosana, based on fine-scale genetic structures of multiple fragmented mountainous populations in Taiwan

Insular species are usually endemic and prone to long-term population reduction, low genetic diversity, and inbreeding depression, which results in difficulties in species conservation. The situation is even more challenging for the glacial relict species whose habitats are usually fragmented in the mountainous regions. Amentotaxus formosana is an endangered and endemic relict tree species in Taiwan. It is typically found scattered across different mountain regions, leading to fragmented and discontinuous populations in various habitats. Consequently, reproductive isolation may lead to deep and cryptic population structures between and within populations. To test this hypothesis and assess the most effective conservation units for the species, comprehensive genetic analyses were conducted using multiple microsatellite loci to uncover fine-scale population structures in three mountainous regions. In this study, a substantial sample of over 600 individuals, predominantly middle-aged, was collected, representing a majority of the existing individuals in the Chachayalaishan and Dawu habitats. Population genetic structure analyses were conducted using three complementary Bayesian clustering approaches (STRUCTURE, InStruct, and GENELAND) to ensure robust identification of distinct genetic clusters in three key habitats (CHA, DAWU, and DL). Results indicated low genetic diversity, distinct genetic differentiation, and severe inbreeding within fragmented populations. Additionally, demographic analysis suggested a substantial downfall in effective population sizes and limited gene flow within and between populations. Based on these findings, we recommend several management strategies to ensure the effective preservation of A. formosana: (1) establishing genetic conservation units corresponding to distinct genetic clusters in the CHA (CG-8-1 to CG-8-8), DAWU (DG-15-1 to DG-15-15), and DL (LG-6-1 to LG-6-6) populations; (2) implementing interpopulation cross-breeding programs to enhance genetic diversity; (3) promoting habitat restoration efforts to establish connectivity among fragmented populations; and (4) implementing vegetative propagation of selected trees for ex-situ conservation, along with establishing conservation nurseries and seed production areas. These localized conservation approaches, combined with the comprehensive genetic insights provided by this study, serve as crucial directives for the genetic monitoring, policy formulation, and sustainable conservation of A. formosana.

The Genetic Differences and Structure of Selected Important Populations of the Endangered Taxus baccata in the Czech Republic

Taxus baccata L. (common yew) is an endangered tree species in the Czech Republic. However, its natural occurrence has not been adequately protected in all areas of the country. The aim of this study is to determine whether the yew population in the newly established Mařeničky seed orchard (TS_L) enables mixing with other Czech yew populations. Using a set of nuclear microsatellites, the genetic diversity in the Lužické Mountains (TS_L) and in selected wild-provenance populations from the Czech Republic (Jílovské yews, TS_J; Březinské yews, TS_B, and yews from Moravský Karst, TS_M) was studied, as they could be donor sources for potential translocation activities. We observed that the level of genetic diversity within the four Czech yew units that were investigated was high. An analysis of the molecular variance (AMOVA) showed 7% variation among populations, and the genetic differentiation values were low to moderate (FST = 0.042–0.108). According to a STRUCTURE analysis, high genetic similarity was observed between the TS_L and TS_B units. Our results provide important genetic suggestions on how conservation management can be designed to maximize its success.

High Level of Phenotypic Differentiation of Common Yew (Taxus baccata L.) Populations in the North-Western Part of the Balkan Peninsula

Common or English yew (Taxus baccata L., Taxaceae) is a conifer species, native to Europe, northern Africa, Asia Minor and Caucasus. It is a dioecious, wind-pollinated and animal-dispersed tree, known for its high-quality wood and medicinal properties, albeit poisonous. The species is rare and has been legally protected at the European and national levels. In addition, its low population density and disjunct character of distribution have reinforced the need for its protection as regeneration is mostly lacking. The aim of this study was to phenotypically characterise the north-western Balkan yew populations. Phenotypic diversity was examined for seven populations, using morphometric analysis of nine phenotypic traits of needles. Descriptive and multivariate statistical analyses were conducted to evaluate the inter- and intrapopulation variability. In addition, to test correlations between geographic, climatic and phenotypic data, Mantel test was used. We identified a geographic structure across studied populations that exhibited high levels of variability on intra- and interpopulation levels. Two groups of populations have been defined and are consistent with previously described genetic divergent lineages from separate refugia. In addition, a significant correlation between phenotypic and geographic data were revealed, i.e., isolation by distance (IBD). However, the Mantel test revealed no significant correlation between morphometric and environmental data. In conclusion, our data reveal that the historical events and persistent IBD acted in combination to produce the morphological patterns observed in common yew populations in the north-western part of the Balkan Peninsula. Finally, we suggested conservation measures to be implemented on a stand level, with habitat preservation as the main goal. In addition, ex situ conservation should be considered, both in the form of collections and urban planting, as both provide additional gene pool reserves.

Distinct Genetic Structure Reflects Ploidy Level Differentiation in Newly Discovered, Extremely Small Populations of Xanthocyparis vietnamensis from Southwestern China

Population genetic assessment is crucial for the conservation and management of threatened species. Xanthocyparis vietnamensis is an endangered species that is currently restricted to karst mountains in southwestern China and Vietnam. This rare conifer was first recorded in 2002 from northern Vietnam and then in 2013 from Guangxi, China, yet nothing is known about its genetic diversity nor ploidy level variation, although previous cytological study suggest that Vietnamese populations are tetraploids. There have been about 45 individuals found to date in Guangxi, China. Here, we genotyped 33 X. vietnamensis individuals using 20 newly developed, polymorphic microsatellite loci, to assess the genetic variability of its extremely small populations. The genetic diversity of X. vietnamensis (H_E = 0.511) was lower than that of two other heliophile species, Calocedrus macrolepis and Fokienia hodginsii, which have similar distribution ranges. This is consistent with the signature of a genetic bottleneck detected in X. vietnamensis. Although the population genetic differentiation coefficient across loci is moderate (F_ST = 0.125), STRUCTURE analysis revealed two distinct genetic clusters, namely the northern and southern population groups; DAPC analysis grouped the southern populations together in one cluster separate from the northern populations; AMOVA analysis detected a significant genetic differentiation between the two population groups (F_RT = 0.089, p < 0.05), and BARRIER analysis detected a genetic barrier between them. Moreover, we detected differentiation in ploidy level between northern and southern populations, sampled individuals from the former and the later are all diploid and tetraploid cytotypes with mean genome sizes of 26.08 and 48.02 pg/2C, respectively. We deduced that heterogeneous geomorphology and historical events (e.g., human deforestation, Quaternary climate oscillations) may have contributed to population fragmentation and small population size in X. vietnamensis. Considering both genetic and ploidy level differentiation, we propose that two different management units (northern and southern) should be considered and a combination of in situ and ex situ conservation measures should be employed to preserve populations of this endangered species in southwestern China in the light of our findings.

EST-SSR-based landscape genetics of Pseudotaxus chienii, a tertiary relict conifer endemic to China

Pseudotaxus chienii, belonging to the monotypic genus Pseudotaxus (Taxaceae), is a relict conifer endemic to China. Its populations are usually small and patchily distributed, having a low capacity of natural regeneration. To gain a clearer understanding of how landscape variables affect the local adaptation of P. chienii, we applied EST-SSR markers in conjunction with landscape genetics methods: (a) to examine the population genetic pattern and spatial genetic structure; (b) to perform genome scan and selection scan to identify outlier loci and the associated landscape variables; and (c) to model the ecological niche under climate change. As a result, P. chienii was found to have a moderate level of genetic variation and a high level of genetic differentiation. Its populations displayed a significant positive relationship between the genetic and geographical distance (i.e., "isolation by distance" pattern) and a strong fine-scale spatial genetic structure within 2 km. A putatively adaptive locus EMS6 (functionally annotated to cellulose synthase A catalytic subunit 7) was identified, which was found significantly associated with soil Cu, K, and Pb content and the combined effects of temperature and precipitation. Moreover, P. chienii was predicted to experience significant range contractions in future climate change scenarios. Our results highlight the potential of specific soil metal content and climate variables as the driving force of adaptive genetic differentiation in P. chienii. The data would also be useful to develop a conservation action plan for P. chienii.

As old as the hills: Pliocene palaeogeographical processes influence patterns of genetic structure in the widespread, common shrub Banksia sessilis

The impact of Quaternary glaciation on the development of phylogeographic structure in plant species is well documented. In unglaciated landscapes, phylogeographic patterns tend to reflect processes relating to persistence and stochasticity, yet other factors, associated with the palaeogeographical history of the landscape, including geomorphological events, can also have a significant influence. The unglaciated landscape of south-western Western Australia is an ideal location to observe these ancient drivers of lineage diversification, with tectonic activity associated with the Darling Fault in the late Pliocene attributed to patterns of deep phylogeographic divergence in a widespread tree from this region. Interestingly, other species within this region have not shown this pattern and this palaeogeographical boundary therefore presents an opportunity to examine age and historical distribution of plant species endemic to this region. In this study, we assess patterns of genetic diversity and structure across 28 populations of the widespread shrub Banksia sessilis using three cpDNA markers and nine nuclear microsatellite markers. Sixteen cpDNA haplotypes were identified, comprising two major chloroplast DNA lineages that are estimated to have diverged in the Pliocene, approximately 3.3 million years ago. This timing coincides with major geomorphological processes in the landscape, including the separation of the Darling Plateau from the adjacent Swan Coastal Plain, as well as eustatic changes on the Swan Coastal Plain that are likely to have resulted in the physical isolation of historical plant lineages. Chloroplast lineages were broadly aligned with populations associated with older lateritic soils of the Darling Plateau and Geraldton sandplains or the younger sandy soils associated with the Swan Coastal Plain and Southern Coastline. This structural pattern of lateritic versus non-lateritic division was not observed in the nuclear microsatellite data that identified three genetic clades that roughly corresponded to populations in the North, South, and Central portions of the distributions.

A multiscale approach to detect selection in nonmodel tree species: Widespread adaptation despite population decline in Taxus baccata L

Detecting the molecular basis of local adaptation and identifying selective drivers is still challenging in nonmodel species. The use of purely population genetic approaches is limited by some characteristics of genetic systems, such as pleiotropy and polygenic control, and parallel evidence from phenotypic-based experimental comparisons is required. In long-lived organisms, the detection of selective pressures might also be precluded by evolutionary lag times in response to the environment. Here, we used the English yew to showcase an example of a multiscale integrative approach in a nonmodel species with limited plant and genomic resources. We combined information from two independent sources, phenotypes in a common environment and genomic data in natural populations, to investigate the signature of selection. Growth differences among populations in a common environment, and phenological patterns of both shoot elongation and male strobili maturation, were associated with climate clines, providing evidence for local adaptation and guiding us in the selection of populations for genomic analyses. We used information on over 25,000 SNPs from c. 1,200 genes to infer the demographic history and to test for molecular signatures of selection at different levels: SNP, gene, and biological pathway. Our results confirmed an overall demographic history of population decline, but we also found evidence for putative local adaptation at the molecular level. We identified or confirmed several candidate genes for positive and negative selection in forest trees, including the pseudo-response regulator 7 (PRR7), an essential component of the circadian clock in plants. In addition, we successfully tested an approach to detect polygenic adaptation in biological pathways, allowing us to identify the flavonoid biosynthesis pathway as a candidate stress-response pathway that deserves further attention in other plants. Finally, our study contributes to the emerging view that explaining contemporary standing genetic variation requires considering adaptation to past climates, especially for long-lived trees.

Spatial genetic structure and diversity of natural populations of Aesculus hippocastanum L. in Greece

Horse-chestnut (Aesculus hippocastanum L.) is an endemic and relict species from the Mediterranean biodiversity hotspot and a popular ornamental tree. Knowledge about the evolutionary history of this species remains scarce. Here, we ask what historical and ecological factors shaped the pattern of genetic diversity and differentiation of this species. We genotyped 717 individuals from nine natural populations using microsatellite markers. The influence of distance, topography and habitat variables on spatial genetic structure was tested within the approaches of isolation-by-distance and isolation-by-ecology. Species niche modeling was used to project the species theoretical range through time and space. The species showed high genetic diversity and moderate differentiation for which topography, progressive range contraction through the species’ history and long-term persistence in stable climatic refugia are likely responsible. A strong geographic component was revealed among five genetic clusters that are connected with very limited gene flow. The environmental variables were a significant factor in the spatial genetic structure. Modeling results indicated that future reduction of the species range may affect its survival. The possible impact of climate changes and high need of in situ conservation are discussed.

Iso-Seq analysis of the Taxus cuspidata transcriptome reveals the complexity of Taxol biosynthesis

BACKGROUND

Taxus cuspidata is well known worldwide for its ability to produce Taxol, one of the top-selling natural anticancer drugs. However, current Taxol production cannot match the increasing needs of the market, and novel strategies should be considered to increase the supply of Taxol. Since the biosynthetic mechanism of Taxol remains largely unknown, elucidating this pathway in detail will be very helpful in exploring alternative methods for Taxol production.

RESULTS

Here, we sequenced Taxus cuspidata transcriptomes with next-generation sequencing (NGS) and third-generation sequencing (TGS) platforms. After correction with Illumina reads and removal of redundant reads, more than 180,000 nonredundant transcripts were generated from the raw Iso-Seq data. Using Cogent software and an alignment-based method, we identified a total of 139 cytochrome P450s (CYP450s), 31 BAHD acyltransferases (ACTs) and 1940 transcription factors (TFs). Based on phylogenetic and coexpression analysis, we identified 9 CYP450s and 7 BAHD ACTs as potential lead candidates for Taxol biosynthesis and 6 TFs that are possibly involved in the regulation of this process. Using coexpression analysis of genes known to be involved in Taxol biosynthesis, we elucidated the stem biosynthetic pathway. In addition, we analyzed the expression patterns of 12 characterized genes in the Taxol pathway and speculated that the isoprene precursors for Taxol biosynthesis were mainly synthesized via the MEP pathway. In addition, we found and confirmed that the alternative splicing patterns of some genes varied in different tissues, which may be an important tissue-specific method of posttranscriptional regulation.

CONCLUSIONS

A strategy was developed to generate corrected full-length or nearly full-length transcripts without assembly to ensure sequence accuracy, thus greatly improving the reliability of coexpression and phylogenetic analysis and greatly facilitating gene cloning and characterization. This strategy was successfully utilized to elucidate the Taxol biosynthetic pathway, which will greatly contribute to the goals of improving the Taxol content in Taxus spp. using molecular breeding or plant management strategies and synthesizing Taxol in microorganisms using synthetic biological technology.

Seed and pollen gene dispersal in Taxus baccata, a dioecious conifer in the face of strong population fragmentation

Background and Aims

Dispersal is crucial due to its direct impact on dynamics of a species' distribution as well as having a role in shaping adaptive potential through gene flow. In plants forming scarce and small populations, knowledge about the dispersal process is required to assess the potential for colonizing new habitats and connectivity of present and future populations. This study aimed to assess dispersal potential in Taxus baccata, a dioecious gymnosperm tree with a wide but highly fragmented distribution.

Methods

Seed and pollen dispersal kernels were estimated directly in the framework of the spatially explicit mating model, where genealogies of naturally established seedlings were reconstructed with the help of microsatellite markers. In this way, six differently shaped dispersal functions were compared.

Key Results

Seed dispersal followed a leptokurtic distribution, with the Exponential-Power, the Power-law and Weibull being almost equally best-fitting models. The pollen dispersal kernel appeared to be more fat-tailed than the seed dispersal kernel, and the Lognormal and the Exponential-Power function showed the best fit. The rate of seed immigration from the background sources was not significantly different from the rate of pollen immigration (13.1 % vs. 19.7 %) and immigration rates were in agreement with or below maximum predictions based on the estimated dispersal kernels. Based on the multimodel approach, 95 % of seeds travel <109 m, while 95 % of pollen travels <704 m from the source.

Conclusions

The results showed that, at a local spatial scale, yew seeds travel shorter distances than pollen, facilitating a rapid development of a kinship structure. At the landscape level, however, although yew exhibits some potential to colonize new habitats through seed dispersal, genetic connectivity between different yew remnants is strongly limited. Taking into account strong population fragmentation, the study suggests that gene dispersal may be a limiting factor of the adaptability of the species.

De novo assembly of English yew (Taxus baccata) transcriptome and its applications for intra- and inter-specific analyses

We provide novel genomic resources for Taxus baccata in the form of a reference transcriptome, SSR and SNP markers, and orthologous single-copy genes, useful for phylogenomic and population genomic applications. English yew (T. baccata) is the only European representative of the Taxaceae family, a conifer group originated in the Jurassic period. The wide extent of environmental heterogeneity within the species' range, together with its long presence in Europe, make English yew an ideal species to investigate adaptive evolution in conifers. To enlarge the genomic resources available for this species, we used Illumina short read sequencing followed by de novo assembly to build the transcriptome of English yew. In addition to a fully annotated transcriptome as well as large sets of new potential SSR and SNP markers for T. baccata, we provide a data set of orthologous single-copy genes across three Taxus species using Picea sitchensis as outgroup, and discuss ortholog uses and limitations for phylogenomic and population genomic applications.

Is the diversification of Mediterranean Basin plant lineages coupled to karyotypic changes?

The Mediterranean Basin region, home to 25,000 plant species, is included in the worldwide list of hotspots of biodiversity. Despite the indisputably important role of chromosome transitions in plant evolution and diversification, no reference study to date has dealt with the possible relationship between chromosome evolution and lineage diversification in the Mediterranean Basin. Here we study patterns of diversification, patterns of chromosome number transition (either polyploidy or dysploidy) and the relationship between the two for 14 Mediterranean Basin angiosperm lineages using previously published phylogenies. We found a mixed pattern, with half of the lineages displaying a change in chromosome transition rates after the onset of the Mediterranean climate (six increases, one decrease) and the other half (six) experiencing constant rates of chromosome transitions through time. We have also found a heterogeneous pattern regarding diversification rates, with lineages exhibiting moderate (five phylogenies) or low (six) initial diversification rates that either increased (six) or declined (five) through time. Our results reveal no clear link between diversification rates and chromosome number transition rates. By promoting the formation of new habitats and driving the extinction of many species, the Mediterranean onset and the posterior Quaternary climatic oscillations could have been key for the establishment of new chromosomal variants in some plant phylogenies but not in others. While the biodiversity of the Mediterranean Basin may be partly influenced by the chromosomal diversity of its lineages, this study concludes that lineage diversification in the region is largely decoupled from karyotypic evolution.

How does the functional diversity of frugivorous birds shape the spatial pattern of seed dispersal? A case study in a relict plant species

Genetic markers used in combination with network analysis can characterize the fine spatial pattern of seed dispersal and assess the differential contribution of dispersers. As a case study, we focus on the seed dispersal service provided by a small guild of frugivorous birds to the common yew, Taxus baccata L., in southern Spain. We build the spatial networks of seed dispersal events between trees and seed-plots within the studied population-local network-and the spatial network that includes all dispersal events-regional network. Such networks are structured in well-defined modules, i.e. groups of tightly connected mother trees and seed-plots. Neither geographical distance, nor microhabitat type explained this modular structure, but when long-distance dispersal events are incorporated in the network it shows a relative increase in overall modularity. Independent field observations suggested the co-occurrence of two complementary groups, short- and long-distance dispersers, mostly contributing to the local and regional seed rain, respectively. The main long-distance disperser at our site, Turdus viscivorus, preferentially visits the most productive trees, thus shaping the seed rain at the landscape scale and affecting the local modular organization. We end by discussing how DNA barcoding could serve to better quantify the role of functional diversity.

Genetic Structure in the Northern Range Margins of Common Ash, Fraxinus excelsior L

During post glacial colonization, loss of genetic diversity due to leading edge effects may be attenuated in forest trees because of their prolonged juvenile phase, allowing many migrants to reach the colonizing front before populations become reproductive. The northern range margins of temperate tree taxa in Europe are particularly suitable to study the genetic processes that follow colonization because they have been little affected by northern refugia. Here we examined how post glacial range dynamics have shaped the genetic structure of common ash (Fraxinus excelsior L.) in its northern range compared to its central range in Europe. We used four chloroplast and six nuclear microsatellites to screen 42 populations (1099 trees), half of which corresponded to newly sampled populations in the northern range and half of which represented reference populations from the central range obtained from previously studies. We found that northern range populations of common ash have the same chloroplast haplotypes as south-eastern European populations, suggesting that colonization of the northern range took place along a single migration route, a result confirmed by the structure at the nuclear microsatellites. Along this route, diversity strongly decreased only in the northern range, concomitantly with increasing population differentiation and complex population substructures, a pattern consistent with a leading edge colonization model. Our study highlights that while diversity is maintained in the central range of common ash due to broad colonizing fronts and high levels of gene flow, it profoundly decreases in the northern range, where colonization was unidirectional and probably involved repeated founder events and population fluctuations. Currently, common ash is threatened by ash dieback, and our results on northern populations will be valuable for developing gene conservation strategies.

Chemical Composition and Interpopulation Variability of Essential Oils of Taxus baccata L. from Serbia

The composition of the essential oil of the twigs and needles of Taxus baccata L. from three natural populations in Serbia has been determined by GC/MS analysis. Of the 91 detected compounds, 87 were identified. The most abundant compound classes were aliphatic alcohols, terpenes, aliphatic hydrocarbons, and aliphatic aldehydes, which together comprised ca. 86.92% of the total oil composition. The dominant constituents were oct-1-en-3-ol (23.48%), (3Z)-hex-3-en-1-ol (11.46%; aliphatic alcohols), and myrtenol (11.38%; oxygenated monoterpene). The PCA of 22 selected compounds revealed differentiations of populations based on geographic distribution. The CA showed that Populations I and II from the Dinaric Alps were similar, and that Population III from the Balkan mountain system was distinct. This was the first investigation of interpopulation variability of T. baccata populations based on essential oil composition. The results of this study were compared with those of previous studies concerning volatile compounds produced by Taxus species. The results indicate that the essential oil content of T. baccata populations from this study is unique, mostly resembling the population from southeast Serbia.

Adapting through glacial cycles: insights from a long-lived tree (Taxus baccata)

Despite the large body of research devoted to understanding the role of Quaternary glacial cycles in the genetic divergence of European trees, the differential contribution of geographic isolation and/or environmental adaptation in creating population genetic divergence remains unexplored. In this study, we used a long-lived tree (Taxus baccata) as a model species to investigate the impact of Quaternary climatic changes on genetic diversity via neutral (isolation-by-distance) and selective (isolation-by-adaptation) processes. We applied approximate Bayesian computation to genetic data to infer its demographic history, and combined this information with past and present climatic data to assess the role of environment and geography in the observed patterns of genetic structure. We found evidence that yew colonized Europe from the East, and that European samples diverged into two groups (Western, Eastern) at the beginning of the Quaternary glaciations, c. 2.2 Myr before present. Apart from the expected effects of geographical isolation during glacials, we discovered a significant role of environmental adaptation during interglacials at the origin of genetic divergence between both groups. This process may be common in other organisms, providing new research lines to explore the effect of Quaternary climatic factors on present-day patterns of genetic diversity.

The threatened epiphytic lichen Lobaria pulmonaria in the Iberian Peninsula: genetic diversity and structure across a latitudinal gradient

The current genetic diversity and structure of a species plays a marked role in the species' future response to environmental changes. Identification of the factors that might ensure the long-term viability of populations along its distribution area is therefore important for conserving biodiversity. In this work, infraspecific genetic diversity and structure of the threatened lichen Lobaria pulmonaria was investigated along a latitudinal gradient, spanning the Spanish latitudinal range of L. pulmonaria. Eighteen populations in Northern, Central, and Southern Spain were analysed using six specific fungal microsatellites of L. pulmonaria. Genetic diversity indices were calculated and compared among populations. Genetic differentiation was assessed using AMOVA and Bayesian methods. Additionally, a redundancy analysis was used to estimate the relative importance of environmental factors on the genetic variation among populations. Annual precipitation was the only factor affecting the genetic diversity probably through its influence on population and thallus size of L. pulmonaria, and significantly higher levels of genetic diversity were detected in southern populations. Isolation by distance was not significant, being environmental variables most important factors controlling genetic variation in L. pulmonaria populations.

Genetic diversity, genetic structure and demographic history of Cycas simplicipinna (Cycadaceae) assessed by DNA sequences and SSR markers

BACKGROUND

Cycas simplicipinna (T. Smitinand) K. Hill. (Cycadaceae) is an endangered species in China. There were seven populations and 118 individuals that we could collect were genotyped in this study. Here, we assessed the genetic diversity, genetic structure and demographic history of this species.

RESULTS

Analyses of data of DNA sequences (two maternally inherited intergenic spacers of chloroplast, cpDNA and one biparentally inherited internal transcribed spacer region ITS4-ITS5, nrDNA) and sixteen microsatellite loci (SSR) were conducted in the species. Of the 118 samples, 86 individuals from the seven populations were used for DNA sequencing and 115 individuals from six populations were used for the microsatellite study. We found high genetic diversity at the species level, low genetic diversity within each of the seven populations and high genetic differentiation among the populations. There was a clear genetic structure within populations of C. simplicipinna. A demographic history inferred from DNA sequencing data indicates that C. simplicipinna experienced a recent population contraction without retreating to a common refugium during the last glacial period. The results derived from SSR data also showed that C. simplicipinna underwent past effective population contraction, likely during the Pleistocene.

CONCLUSIONS

Some genetic features of C. simplicipinna such as having high genetic differentiation among the populations, a clear genetic structure and a recent population contraction could provide guidelines for protecting this endangered species from extinction. Furthermore, the genetic features with population dynamics of the species in our study would help provide insights and guidelines for protecting other endangered species effectively.

Spatial genetic structure in Beta vulgaris subsp. maritima and Beta macrocarpa reveals the effect of contrasting mating system, influence of marine currents, and footprints of postglacial recolonization routes

Understanding the factors that contribute to population genetic divergence across a species' range is a long-standing goal in evolutionary biology and ecological genetics. We examined the relative importance of historical and ecological features in shaping the present-day spatial patterns of genetic structure in two related plant species, Beta vulgaris subsp. maritima and Beta macrocarpa. Using nuclear and mitochondrial markers, we surveyed 93 populations from Brittany (France) to Morocco – the southern limit of their species' range distribution. Whereas B. macrocarpa showed a genotypic structure and a high level of genetic differentiation indicative of selfing, the population genetic structure of B. vulgaris subsp. maritima was consistent with an outcrossing mating system. We further showed (1) a strong geographic clustering in coastal B. vulgaris subsp. maritima populations that highlighted the influence of marine currents in shaping different lineages and (2) a peculiar genetic structure of inland B. vulgaris subsp. maritima populations that could indicate the admixture of distinct evolutionary lineages and recent expansions associated with anthropogenic disturbances. Spatial patterns of nuclear diversity and differentiation also supported a stepwise recolonization of Europe from Atlantic-Mediterranean refugia after the last glacial period, with leading-edge expansions. However, cytoplasmic diversity was not impacted by postglacial recolonization: stochastic long-distance seed dispersal mediated by major oceanic currents may mitigate the common patterns of reduced cytoplasmic diversity observed for edge populations. Overall, the patterns we documented here challenge the general view of reduced genetic diversity at the edge of a species' range distribution and provide clues for understanding how life-history and major geographic features interact to shape the distribution of genetic diversity.

Phylogeography and genetic effects of habitat fragmentation on endangered Taxus yunnanensis in southwest China as revealed by microsatellite data

It is not known how the profoundly complex topography and habitat heterogeneity generated by the uplift of the Qinghai-Tibetan Plateau (QTP) during the late Tertiary affected population genetic structure of endangered Taxus yunnanensis. In addition, the effects of habitat fragmentation due to anthropogenic disturbance on genetic diversity and population differentiation of this species have not been studied. T. yunnanensis is an ancient tree/shrub mainly distributed in southwest China. Recently, the species has suffered a sharp decline due to excessive logging for its famous anticancer metabolite taxol, resulting in smaller and more isolated populations. To understand the phylogeography and genetic consequences of habitat fragmentation of this endangered species, using 11 polymorphic microsatellites, we genotyped 288 individuals from 14 populations from a range-wide sampling in China. Our results suggest that two different population groups that were once isolated have persisted in situ during glacial periods in both areas, and have not merged since. Habitat fragmentation has led to significant genetic bottlenecks, high inbreeding and population divergence in this species. The two different population groups of T. yunnanensis could be attributed to restricted gene flow caused through isolation by geographical barriers and by habitat heterogeneity during uplift of the QTP, or the existence of two separate glacial refugia during the Pleistocene. In situ and ex situ conservation of the two Evolutionarily Significant Units (ESUs), artificial gene flow between populations and a comprehensive understanding of the pollination system in this endangered species are suggested from this study.

Climate warming and the decline of Taxus airborne pollen in urban pollen rain (Emilia Romagna, northern Italy)

Woody plant performance in a changing global environment has always been at the centre of palaeoenvironmental and long-term climate reconstructions carried out by means of pollen analysis. In Mediterranean regions, Taxus constitutes the highest percentage in past pollen diagrams from cold or cool periods, and therefore it is generally considered a good index to infer climate features from past records. However, a comparison of these inferences with the true current trends in pollen production has not been attemped until now. This study reports the decline of airborne pollen of Taxus observed in Emilia Romagna, a region of northern Italy, during the period 1990-2007. Phenological observations on four male specimens and microscopic examination of fresh pollen were made in order to check Taxus flowering time and pollen morphology. Airborne pollen was monitored through continuous sampling with a Hirst volumetric sampler. In the 18-year long period of investigation, Taxus pollen production has decreased, while total woody pollen abundance in air has increased. The trend of the Taxus pollen season shows a delay at the beginning, a shortening of the pollen period, and an advance of the end of the pollen season. This was interpreted as a response to climate warming. In particular, Taxus follows the behaviour of winter-flowering plants, and therefore earlier pollination is favoured at low autumn temperatures, while late pollination occurs more often, most likely after warm autumn temperatures.

Environmental heterogeneity explains the genetic structure of Continental and Mediterranean populations of Fraxinus angustifolia Vahl

Tree species with wide distributions often exhibit different levels of genetic structuring correlated to their environment. However, understanding how environmental heterogeneity influences genetic variation is difficult because the effects of gene flow, drift and selection are confounded. We investigated the genetic variation and its ecological correlates in a wind-pollinated Mediterranean tree species, Fraxinus angustifolia Vahl, within a recognised glacial refugium in Croatia. We sampled 11 populations from environmentally divergent habitats within the Continental and Mediterranean biogeographical regions. We combined genetic data analyses based on nuclear microsatellite loci, multivariate statistics on environmental data and ecological niche modelling (ENM). We identified a geographic structure with a high genetic diversity and low differentiation in the Continental region, which contrasted with the significantly lower genetic diversity and higher population divergence in the Mediterranean region. The positive and significant correlation between environmental and genetic distances after controlling for geographic distance suggests an important influence of ecological divergence of the sites in shaping genetic variation. The ENM provided support for niche differentiation between the populations from the Continental and Mediterranean regions, suggesting that contemporary populations may represent two divergent ecotypes. Ecotype differentiation was also supported by multivariate environmental and genetic distance analyses. Our results suggest that despite extensive gene flow in continental areas, long-term stability of heterogeneous environments have likely promoted genetic divergence of ashes in this region and can explain the present-day genetic variation patterns of these ancient populations.

Increased inbreeding and strong kinship structure in Taxus baccata estimated from both AFLP and SSR data

Habitat fragmentation can have severe genetic consequences for trees, such as increased inbreeding and decreased effective population size. In effect, local populations suffer from reduction of genetic variation, and thus loss of adaptive capacity, which consequently increases their risk of extinction. In Europe, Taxus baccata is among a number of tree species experiencing strong habitat fragmentation. However, there is little empirical data on the population genetic consequences of fragmentation for this species. This study aimed to characterize local genetic structure in two natural remnants of English yew in Poland based on both amplified fragment length polymorphism (AFLP) and microsatellite (SSR) markers. We introduced a Bayesian approach that estimates the average inbreeding coefficient using AFLP (dominant) markers. Results showed that, in spite of high dispersal potential (bird-mediated seed dispersal and wind-mediated pollen dispersal), English yew populations show strong kinship structure, with a spatial extent of 50-100 m, depending on the population. The estimated inbreeding levels ranged from 0.016 to 0.063, depending on the population and marker used. Several patterns were evident: (1) AFLP markers showed stronger kinship structure than SSRs; (2) AFLP markers provided higher inbreeding estimates than SSRs; and (3) kinship structure and inbreeding were more pronounced in denser populations regardless of the marker used. Our results suggest that, because both kinship structure and (bi-parental) inbreeding exist in populations of English yew, gene dispersal can be fairly limited in this species. Furthermore, at a local scale, gene dispersal intensity can be more limited in a dense population.

The Strait of Gibraltar as a major biogeographic barrier in Mediterranean conifers: a comparative phylogeographic survey

The Strait of Gibraltar (SG) is reputed for being both a bridge and a geographic barrier to biological exchanges between Europe and Africa. Major genetic breaks associated with this strait have been identified in various taxa, but it is unknown whether these disjunctions have been produced simultaneously or by independent biogeographic processes. Here, the genetic structure of five conifers distributed on both sides of the SG was investigated using mitochondrial (nad1 b/c, nad5-1, nad5-4 and nad7-1) and chloroplast (Pt1254, Pt15169, Pt30204, Pt36480, Pt71936 and Pt87268) DNA markers. The distribution of genetic variation was partially congruent between types of markers within the same species. Across taxa, there was a significant overlapping between the SG and the genetic breaks detected, especially for the four Tertiary species surveyed (Abies pinsapo complex, Pinus nigra, Pinus pinaster and Taxus baccata). For most of these taxa, the divergence of populations across the SG could date back to long before the Pleistocene glaciations. However, their strongly different cpDNA G(ST) and R(ST) values point out that they have had dissimilar population histories, which might include contrasting amounts of pollen-driven gene flow since their initial establishment in the region. The fifth species, Pinus halepensis, was genetically depauperated and homogenous on both sides of the SG. A further analysis of nuclear DNA sequences with coalescent-based isolation with migration models suggests a Pleistocene divergence of P. halepensis populations across the SG, which is in sharp contrast with the pre-Pleistocene divergence dates obtained for P. pinaster. Altogether, these results indicate that the genetic breaks observed across this putative biogeographical barrier have been produced by independent evolutionary processes related to the biological history of each individual species instead of a common vicariant phenomenon.