WOODIV, a database of occurrences, functional traits, and phylogenetic data for all Euro-Mediterranean trees.

Trees play a key role in the structure and function of many ecosystems worldwide. In the Mediterranean Basin, forests cover approximately 22% of the total land area hosting a large number of endemics (46 species). Despite its particularities and vulnerability, the biodiversity of Mediterranean trees is not well known at the taxonomic, spatial, functional, and genetic levels required for conservation applications. The WOODIV database fills this gap by providing reliable occurrences, four functional traits (plant height, seed mass, wood density, and specific leaf area), and sequences from three DNA-regions (rbcL, matK, and trnH-psbA), together with modelled occurrences and a phylogeny for all 210 Euro-Mediterranean tree species. We compiled, homogenized, and verified occurrence data from sparse datasets and collated them on an INSPIRE-compliant 10 × 10 km grid. We also gathered functional trait and genetic data, filling existing gaps where possible. The WOODIV database can benefit macroecological studies in the fields of conservation, biogeography, and community ecology.

Functional Relationships of Wood Anatomical Traits in Norway Spruce.

The quantitative assessment of wood anatomical traits offers important insights into those factors that shape tree growth. While it is known that conduit diameter, cell wall thickness, and wood density vary substantially between and within species, the interconnection between wood anatomical traits, tree-ring width, tree height and age, as well as environment effects on wood anatomy remain unclear. Here, we measure and derived 65 wood anatomical traits in cross-sections of the five outermost tree rings (2008-2012) of 30 Norway spruce [Picea abies (L.) H. Karst.] trees growing along an altitudinal gradient (1,400-1,750 m a.s.l.) in the northern Apennines (Italy). We assess the relationship among each anatomical trait and between anatomical trait groups according to their function for (i) tree-ring growth, (ii) cell growth, (iii) hydraulic traits, and (iv) mechanical traits. The results show that tree height significantly affects wood hydraulic traits, as well as number and tangential diameter of tracheids, and ultimately the total ring width. Moreover, the amount of earlywood and latewood percentage influence wood hydraulic safety and efficiency, as well as mechanical traits. Mechanically relevant wood anatomical traits are mainly influenced by tree age, not necessarily correlated with tree height. An additional level of complexity is also indicated by some anatomical traits, such as latewood lumen diameter and the cell wall reinforcement index, showing large inter-annual variation as a proxy of phenotypic plasticity. This study unravels the complex interconnection of tree-ring tracheid structure and identifies anatomical traits showing a large inter-individual variation and a strong interannual coherency. Knowing and quantifying anatomical variation in cells of plant stem is crucial in ecological and biological studies for an appropriate interpretation of abiotic drivers of wood formation often related to tree height and/or tree age.

Individual reproductive success in Norway spruce natural populations depends on growth rate, age and sensitivity to temperature.

Quantifying the individual reproductive success and understanding its determinants is a central issue in evolutionary research for the major consequences that the transmission of genetic variation from parents to offspring has on the adaptive potential of populations. Here, we propose to distil the myriad of information embedded in tree-ring time series into a set of tree-ring-based phenotypic traits to be investigated as potential drivers of reproductive success in forest trees. By using a cross-disciplinary approach that combines parentage analysis and a thorough dendrophenotypic characterisation of putative parents, we assessed sex-specific relationships between such dendrophenotypic traits (i.e., age, growth rate and parameters describing sensitivity to climate and to extreme climatic events) and reproductive success in Norway spruce. We applied a full probability method for reconstructing parent-offspring relationships between 604 seedlings and 518 adult trees sampled within five populations from southern and central Europe. We found that individual female and male reproductive success was positively associated with tree growth rate and age. Female reproductive success was also positively influenced by the correlation between growth and the mean temperature of the previous vegetative season. Overall, our results showed that Norway spruce individuals with the highest fitness are those who are able to keep high-growth rates despite potential growth limitations caused by reproductive costs and climatic limiting conditions. Identifying such functional links between the individual ecophysiological behaviour and its evolutionary gain would increase our understanding on how natural selection shapes the genetic composition of forest tree populations over time.

Disentangling the effects of spatial proximity and genetic similarity on individual growth performances in Norway spruce natural populations.

Cambial growth is a phenotypic trait influenced by various physiological processes, numerous biotic and abiotic drivers, as well as by the genetic background. By archiving the outcome of such complex interplay, tree-rings are an exceptional resource for addressing individual long-term growth responses to changing environments and climate. Disentangling the effects of the different drivers of tree growth, however, remains challenging because of the lack of multidisciplinary data. Here, we combine individual dendrochronological, genetic and spatial data to assess the relative importance of genetic similarity and spatial proximity on Norway spruce (Picea abies (L.) Karst.) growth performances. We intensively sampled five plots from two populations in southern and central Europe, characterizing a total of 482 trees. A two-step analytical framework was developed. First, the effects of climate and tree age on tree-ring width (TRW) were estimated for each individual using a random slope linear mixed-effects model. Individual parameters were then tested against genetic and spatial variables by Mantel tests, partial redundancy analyses and variance partitioning. Our modelling approach successfully captured a large fraction of variance in TRW (conditional R2 values up to 0.94) which was largely embedded in inter-individual differences. All statistical approaches consistently showed that genetic similarity was not related to variation in the individual parameters describing growth responses. In contrast, up to 29% of the variance of individual parameters was accounted by spatial variables, revealing that microenvironmental features are more relevant than genetic similarity in determining similar growth patterns. Our study highlights both the advantages of modelling dendrochronological data at the individual level and the relevance of microenvironmental variation on individual growth patterns. These two aspects should be carefully considered in future multidisciplinary studies on growth dynamics in natural populations.

Genetic diversity and divergence at the Arbutus unedo L. (Ericaceae) westernmost distribution limit.

Mediterranean forests are fragile ecosystems vulnerable to recent global warming and reduction of precipitation, and a long-term negative effect is expected on vegetation with increasing drought and in areas burnt by fires. We investigated the spatial distribution of genetic variation of Arbutus unedo in the western Iberia Peninsula, using plastid markers with conservation and provenance regions design purposes. This species is currently undergoing an intense domestication process in the region, and, like other species, is increasingly under the threat from climate change, habitat fragmentation and wildfires. We sampled 451 trees from 15 natural populations from different ecological conditions spanning the whole species' distribution range in the region. We applied Bayesian analysis and identified four clusters (north, centre, south, and a single-population cluster). Hierarchical AMOVA showed higher differentiation among clusters than among populations within clusters. The relatively low within-clusters differentiation can be explained by a common postglacial history of nearby populations. The genetic structure found, supported by the few available palaeobotanical records, cannot exclude the hypothesis of two independent A. unedo refugia in western Iberia Peninsula during the Last Glacial Maximum. Based on the results we recommend a conservation strategy by selecting populations for conservation based on their allelic richness and diversity and careful seed transfer consistent with current species' genetic structure.

The eastern part of the Fertile Crescent concealed an unexpected route of olive (Olea europaea L.) differentiation.

Background and Aims: Olive is considered a native plant of the eastern side of the Mediterranean basin, from where it should have spread westward along the Mediterranean shores, while little is known about its diffusion in the eastern direction. Methods: Genetic diversity levels and population genetic structure of a wide set of olive ecotypes and varieties collected from several provinces of Iran, representing a high percentage of the entire olive resources present in the area, was screened with 49 chloroplast and ten nuclear simple sequence repeat markers, and coupled with archaeo-botanical and historical data on Mediterranean olive varieties. Approximate Bayesian Computation was applied to define the demographic history of olives including Iranian germplasm, and species distribution modelling was performed to understand the impact of the Late Quaternary on olive distribution. Key Results: The results of the present study demonstrated that: (1) the climatic conditions of the last glacial maximum had an important role on the actual olive distribution, (2) all Iranian olive samples had the same maternal inheritance as Mediterranean cultivars, and (3) the nuclear gene flow from the Mediterranean basin to the Iranian plateau was almost absent, as well as the contribution of subspecies cuspidata to the diversity of Iranian olives. Conclusions: Based on this evidence, a new scenario for the origin and distribution of this important fruit crop has been traced. The evaluation of olive trees growing in the eastern part of the Levant highlighted a new perspective on the spread and distribution of olive, suggesting two routes of olive differentiation, one westward, spreading along the Mediterranean basin, and another moving towards the east and reaching the Iranian plateau before its domestication.

Reconstructing the plant mitochondrial genome for marker discovery: a case study using Pinus.

Whole-genome-shotgun (WGS) sequencing of total genomic DNA was used to recover ~1 Mbp of novel mitochondrial (mtDNA) sequence from Pinus sylvestris (L.) and three members of the closely related Pinus mugo species complex. DNA was extracted from megagametophyte tissue from six mother trees from locations across Europe, and 100-bp paired-end sequencing was performed on the Illumina HiSeq platform. Candidate mtDNA sequences were identified by their size and coverage characteristics, and by comparison with published plant mitochondrial genomes. Novel variants were identified, and primers targeting these loci were trialled on a set of 28 individuals from across Europe. In total, 31 SNP loci were successfully resequenced, characterizing 15 unique haplotypes. This approach offers a cost-effective means of developing marker resources for mitochondrial genomes in other plant species where reference sequences are unavailable.

Multispecies genetic structure and hybridization in the Betula genus across Eurasia.

Boreal and cool temperate forests are the major land cover of northern Eurasia, and information about continental-scale genetic structure and past demographic history of forest species is important from an evolutionary perspective and has conservation implications. However, although many population genetic studies of forest tree species have been conducted in Europe or Eastern Asia, continental-scale genetic structure and past demographic history remain poorly known. Here, we focus on the birch genus Betula, which is commonly distributed in boreal and cool temperate forests, and examine 129 populations of two tetraploid and four diploid species collected from Iceland to Japan. All individuals were genotyped at seven to 18 nuclear simple sequence repeats (nSSRs). Pairwise FST' among the six species ranged from 0.285 to 0.903, and genetic differentiation among them was clear. structure analysis suggested that Betula pubescens is an allotetraploid and one of the parental species was Betula pendula. In both species pairs of B. pendula and B. plathyphylla, and B. pubescens and B. ermanii, genetic diversity was highest in central Siberia. A hybrid zone was detected around Lake Baikal for eastern and western species pairs regardless of ploidy level. Approximate Bayesian computation suggested that the divergence of B. pendula and B. platyphylla occurred around the beginning of the last ice age (36 300 years BP, 95% CI: 15 330-92 700) and hybridization between them was inferred to have occurred after the last glacial maximum (1614 years BP, 95% CI: 561-4710), with B. pendula providing a higher contribution to hybrids.

Local Adaptation in European Firs Assessed through Extensive Sampling across Altitudinal Gradients in Southern Europe.

BACKGROUND: Local adaptation is a key driver of phenotypic and genetic divergence at loci responsible for adaptive traits variations in forest tree populations. Its experimental assessment requires rigorous sampling strategies such as those involving population pairs replicated across broad spatial scales. METHODS: A hierarchical Bayesian model of selection (HBM) that explicitly considers both the replication of the environmental contrast and the hierarchical genetic structure among replicated study sites is introduced. Its power was assessed through simulations and compared to classical 'within-site' approaches (FDIST, BAYESCAN) and a simplified, within-site, version of the model introduced here (SBM). RESULTS: HBM demonstrates that hierarchical approaches are very powerful to detect replicated patterns of adaptive divergence with low false-discovery (FDR) and false-non-discovery (FNR) rates compared to the analysis of different sites separately through within-site approaches. The hypothesis of local adaptation to altitude was further addressed by analyzing replicated Abies alba population pairs (low and high elevations) across the species' southern distribution range, where the effects of climatic selection are expected to be the strongest. For comparison, a single population pair from the closely related species A. cephalonica was also analyzed. The hierarchical model did not detect any pattern of adaptive divergence to altitude replicated in the different study sites. Instead, idiosyncratic patterns of local adaptation among sites were detected by within-site approaches. CONCLUSION: Hierarchical approaches may miss idiosyncratic patterns of adaptation among sites, and we strongly recommend the use of both hierarchical (multi-site) and classical (within-site) approaches when addressing the question of adaptation across broad spatial scales.

The extent and meaning of hybridization and introgression between Siberian spruce (Picea obovata) and Norway spruce (Picea abies): cryptic refugia as stepping stones to the west?

Boreal species were repeatedly exposed to ice ages and went through cycles of contraction and expansion while sister species alternated periods of contact and isolation. The resulting genetic structure is consequently complex, and demographic inferences are intrinsically challenging. The range of Norway spruce (Picea abies) and Siberian spruce (Picea obovata) covers most of northern Eurasia; yet their geographical limits and histories remain poorly understood. To delineate the hybrid zone between the two species and reconstruct their joint demographic history, we analysed variation at nuclear SSR and mitochondrial DNA in 102 and 88 populations, respectively. The dynamics of the hybrid zone was analysed with approximate Bayesian computation (ABC) followed by posterior predictive structure plot reconstruction and the presence of barriers across the range tested with estimated effective migration surfaces. To estimate the divergence time between the two species, nuclear sequences from two well-separated populations of each species were analysed with ABC. Two main barriers divide the range of the two species: one corresponds to the hybrid zone between them, and the other separates the southern and northern domains of Norway spruce. The hybrid zone is centred on the Urals, but the genetic impact of Siberian spruce extends further west. The joint distribution of mitochondrial and nuclear variation indicates an introgression of mitochondrial DNA from Norway spruce into Siberian spruce. Overall, our data reveal a demographic history where the two species interacted frequently and where migrants originating from the Urals and the West Siberian Plain recolonized northern Russia and Scandinavia using scattered refugial populations of Norway spruce as stepping stones towards the west.

Genomic resources notes accepted 1 October 2014-30 November 2014.

This article documents the public availability of (i) RAD sequencing data and validated SNPs for the American mink Neovison vison and (ii) Transcriptome resources for two nonmodel freshwater crustacean species, the copepod Eucyclops serrulatus and the amphipod Echinogammarus veneris.

Limited genetic variability and phenotypic plasticity detected for cavitation resistance in a Mediterranean pine.

Resistance to cavitation is a major determinant of plant survival under severe drought and can be used to quantify species adaptive potential. Interspecific variation in this key trait is well defined in woody species, but intraspecific variation (level and structure) resulting from standing genetic variation and phenotypic plasticity has never been determined. Combining for the first time in situ characterization of natural populations and two reciprocal common gardens in dry and wet sites, we estimated variance components (phenotypic, genetic, environmental, and genetic × environmental) of cavitation resistance based on 513 genotypes of a Mediterranean pine, Pinus pinaster. Despite the selected populations being climatically contrasted, phenotypic plasticity in resistance to cavitation remained low and was essentially attributed to family level. Between-population variation in cavitation resistance for both phenotypic and genetic variation was limited. These results strongly suggest that cavitation resistance is buffered against genetic and to a lesser extent environmental variation (canalization) in maritime pine. Consequently, in a drier world, the increasing drought tolerance of Pinus species might be severely constrained by the low level of cavitation resistance variation, resulting in a large-scale loss of productivity.

Micro- and macro-geographic scale effect on the molecular imprint of selection and adaptation in Norway spruce.

Forest tree species of temperate and boreal regions have undergone a long history of demographic changes and evolutionary adaptations. The main objective of this study was to detect signals of selection in Norway spruce (Picea abies [L.] Karst), at different sampling-scales and to investigate, accounting for population structure, the effect of environment on species genetic diversity. A total of 384 single nucleotide polymorphisms (SNPs) representing 290 genes were genotyped at two geographic scales: across 12 populations distributed along two altitudinal-transects in the Alps (micro-geographic scale), and across 27 populations belonging to the range of Norway spruce in central and south-east Europe (macro-geographic scale). At the macrogeographic scale, principal component analysis combined with Bayesian clustering revealed three major clusters, corresponding to the main areas of southern spruce occurrence, i.e. the Alps, Carpathians, and Hercynia. The populations along the altitudinal transects were not differentiated. To assess the role of selection in structuring genetic variation, we applied a Bayesian and coalescent-based F(ST)-outlier method and tested for correlations between allele frequencies and climatic variables using regression analyses. At the macro-geographic scale, the F(ST)-outlier methods detected together 11 F(ST)-outliers. Six outliers were detected when the same analyses were carried out taking into account the genetic structure. Regression analyses with population structure correction resulted in the identification of two (micro-geographic scale) and 38 SNPs (macro-geographic scale) significantly correlated with temperature and/or precipitation. Six of these loci overlapped with F(ST)-outliers, among them two loci encoding an enzyme involved in riboflavin biosynthesis and a sucrose synthase. The results of this study indicate a strong relationship between genetic and environmental variation at both geographic scales. It also suggests that an integrative approach combining different outlier detection methods and population sampling at different geographic scales is useful to identify loci potentially involved in adaptation.

Within-population genetic structure in beech (Fagus sylvatica L.) stands characterized by different disturbance histories: does forest management simplify population substructure?

The fine-scale assessment of both spatially and non-spatially distributed genetic variation is crucial to preserve forest genetic resources through appropriate forest management. Cryptic within-population genetic structure may be more common than previously thought in forest tree populations, which has strong implications for the potential of forests to adapt to environmental change. The present study was aimed at comparing within-population genetic structure in European beech (Fagus sylvatica L.) plots experiencing different disturbance levels. Five plot pairs made up by disturbed and undisturbed plots having the same biogeographic history were sampled throughout Europe. Overall, 1298 individuals were analyzed using four highly polymorphic nuclear microsatellite markers (SSRs). Bayesian clustering within plots identified 3 to 11 genetic clusters (within-plot θ ST ranged from 0.025 to 0.124). The proportion of within-population genetic variation due to genetic substructuring (F CluPlot = 0.067) was higher than the differentiation among the 10 plots (F PlotTot = 0.045). Focusing on the comparison between managed and unmanaged plots, disturbance mostly explains differences in the complexity of within-population genetic structure, determining a reduction of the number of genetic clusters present in a standardized area. Our results show that: i) genetic substructuring needs to be investigated when studying the within-population genetic structure in forest tree populations, and ii) indices describing subtle characteristics of the within-population genetic structure are good candidates for providing early signals of the consequences of forest management, and of disturbance events in general.

Population structure of Cynara cardunculus complex and the origin of the conspecific crops artichoke and cardoon.

BACKGROUND AND AIMS: Globe artichoke and leafy cardoon, two crops within the same species Cynara cardunculus, are traditionally cultivated in the Mediterranean region and play a significant role in the agricultural economy of this area. The two cultigens have different reproductive systems: artichoke is generally vegetatively propagated, while leafy cardoon is seed propagated. The domestication events underlying the origin of both artichoke and cultivated cardoon from their wild relative and the area of occurrence are not yet fully understood. The aim of this study was to investigate population structure in wild cardoon, globe artichoke and leafy cardoon material and infer domestication events. METHODS: Thirty-five microsatellite (simple sequence repeat) markers, distributed in the C. cardunculus genome, and a large geographical and numerical sampling in southern Europe and North Africa were used to assess population structure and diversity. KEY RESULTS: The results suggest the presence of two distinct domestication events for artichoke and leafy cardoon, and also suggest a new possible scenario, with western wild cardoon having originated from cultivated cardoon escaped from cultivation. Evidence was found for a demographic bottleneck in the past history of globe artichoke. CONCLUSIONS: The results shed new light on the relationships between the three taxa of C. cardunculus and highlight relevant aspects on the evolution of domestication of two crops with a different reproductive system within the same species. It is proposed that the probable centre of origin of artichoke is located in southern Italy, probably Sicily.

Range-wide phylogeography and gene zones in Pinus pinaster Ait. revealed by chloroplast microsatellite markers.

Some 1339 trees from 48 Pinus pinaster stands were characterized by five chloroplast microsatellites, detecting a total of 103 distinct haplotypes. Frequencies for the 16 most abundant haplotypes (p(k) > 0.01) were spatially interpolated over a lattice made by 430 grid points. Fitting of spatially interpolated values on raw haplotype frequencies at the same geographical location was tested by regression analysis. A range-wide 'diversity map' based on interpolated haplotype frequencies allowed the identification of one 'hotspot' of diversity in central and southeastern Spain, and two areas of low haplotypic diversity located in the western Iberian peninsula and Morocco. Principal component analysis (PCA) carried out on haplotypes frequency surfaces allowed the construction of a colour-based 'synthetic' map of the first three PC components, enabling the detection of the main range-scale genetic trends and the identification of three main 'gene pools' for the species: (i) a 'southeastern' gene pool, including southeastern France, Italy, Corsica, Sardinia, Pantelleria and northern Africa; (ii) an 'Atlantic' gene pool, including all the western areas of the Iberian peninsula; and (iii) a 'central' gene pool, located in southeastern Spain. Multivariate and AMOVA analyses carried out on interpolated grid point frequency values revealed the existence of eight major clusters ('gene zones'), whose genetic relationships were related with the history of the species. In addition, demographic models showed more ancient expansions in the eastern and southern ranges of maritime pine probably associated to early postglacial recolonization. The delineation of the gene zones provides a baseline for designing conservation areas in this key Mediterranean pine.