Federico Sebastiani

Relaxed Molecular Clock Provides Evidence for Long-Distance Dispersal of Nothofagus (Southern Beech)

Nothofagus (southern beech), with an 80-million-year-old fossil record, has become iconic as a plant genus whose ancient Gondwanan relationships reach back into the Cretaceous era. Closely associated with Wegeners theory of “Kontinentaldrift”, Nothofagus has been regarded as the “key genus in plant biogeography”. This paradigm has the New Zealand species as passengers on a Moas Ark that rafted away from other landmasses following the breakup of Gondwana. An alternative explanation for the current transoceanic distribution of species seems almost inconceivable given that Nothofagus seeds are generally thought to be poorly suited for dispersal across large distances or oceans. Here we test the Moas Ark hypothesis using relaxed molecular clock methods in the analysis of a 7.2-kb fragment of the chloroplast genome. Our analyses provide the first unequivocal molecular clock evidence that, whilst some Nothofagus transoceanic distributions are consistent with vicariance, trans-Tasman Sea distributions can only be explained by long-distance dispersal. Thus, our analyses support the interpretation of an absence of Lophozonia and Fuscospora pollen types in the New Zealand Cretaceous fossil record as evidence for Tertiary dispersals of Nothofagus to New Zealand. Our findings contradict those from recent cladistic analyses of biogeographic data that have concluded transoceanic Nothofagus distributions can only be explained by vicariance events and subsequent extinction. They indicate that the biogeographic history of Nothofagus is more complex than envisaged under opposing polarised views expressed in the ongoing controversy over the relevance of dispersal and vicariance for explaining plant biodiversity. They provide motivation and justification for developing more complex hypotheses that seek to explain the origins of Southern Hemisphere biota.

The distribution of Quercus suber chloroplast haplotypes matches the palaeogeographical history of the western Mediterranean

Combining molecular analyses with geological and palaeontological data may reveal timing and modes for the divergence of lineages within species. The Mediterranean Basin is particularly appropriate for this kind of multidisciplinary studies, because of its complex geological history and biological diversity. Here, we investigated chloroplast DNA of Quercus suber populations in order to detect possible relationships between their geographical distribution and the palaeogeographical history of the western Mediterranean domain. We analysed 110 cork oak populations, covering the whole distribution range of the species, by 14 chloroplast microsatellite markers, among which eight displayed variation among populations. We identified five haplotypes whose distribution is clearly geographically structured. Results demonstrated that cork oak populations have undergone a genetic drift geographically consistent with the Oligocene and Miocene break‐up events of the European–Iberian continental margin and suggested that they have persisted in a number of separate microplates, currently found in Tunisia, Sardinia, Corsica, and Provence, without detectable chloroplast DNA modifications for a time span of over 15 million years. A similar distribution pattern of mitochondrial DNA of Pinus pinaster supports the hypothesis of such long‐term persistence, in spite of Quaternary climate oscillations and of isolation due to insularity, and suggests that part of the modern geographical structure of Mediterranean populations may be traced back to the Tertiary history of taxa.

GENETICALLY DEPAUPERATE BUT WIDESPREAD: THE CASE OF AN EMBLEMATIC MEDITERRANEAN PINE

Abstract Genetic variation is generally considered a prerequisite for adaptation to new environmental conditions. Thus the discovery of genetically depauperate but geographically widespread species is unexpected. We used 12 paternally inherited chloroplast microsatellites to estimate population genetic variation across the full range of an emblematic circum-Mediterranean conifer, stone pine (Pinus pinea L.). The same chloroplast DNA haplotype is fixed in nearly all of the 34 investigated populations. Such a low level of variation is consistent with a previous report of very low levels of diversity at nuclear loci in this species. Stone pine appears to have passed through a severe and prolonged demographic bottleneck, followed by subsequent natural- and human-mediated dispersal across the Mediterranean Basin. No other abundant and widespread plant species has as little genetic diversity as P. pinea at both chloroplast and nuclear markers. However, the species harbors a nonnegligible amount of variation at adaptive traits. Thus a causal relationship between genetic diversity, as measured by marker loci, and the evolutionary precariousness of a species, cannot be taken for granted.

Chloroplast DNA phylogeography of European ashes, Fraxinus sp. (Oleaceae): roles of hybridization and life history traits.

We investigated range‐wide phylogeographic variation in three European ash species (Fraxinus sp., Oleaceae). Chloroplast DNA (cpDNA) microsatellites were typed in the thermophilous Fraxinus angustifolia and Fraxinus ornus and the observed haplotypes and the geographic distribution of diversity were compared to cpDNA data previously obtained in the more cold‐tolerant Fraxinus excelsior. We found wide‐ranging haplotype sharing between the phylogenetically close F. angustifolia and F. excelsior, suggesting hybridization (i) in common glacial refuges in the Iberian Peninsula, northern Italy, the eastern and/or Dinaric Alps and the Balkan Peninsula, and/or (ii) during postglacial recolonization. The data allowed us to propose additional glacial refuges for F. angustifolia in southern Italy and in Turkey, and populations from the latter region were particularly polymorphic. There was evidence for refuge areas in Italy, the Balkan Peninsula and Turkey for F. ornus, which did not share any single chloroplast haplotype with the other species. In both F. angustifolia and F. ornus, cpDNA diversity (hS = 0.027 and hS = 0.009, respectively) was lower and fixation levels (GST = 0.964 and GST = 0.983, respectively) higher than in sympatric F. excelsior (hS = 0.096, GST = 0.870). These diversity patterns could be due to temperature tolerance or the demographic history.

A fast and cost-effective approach to develop and map EST-SSR markers: oak as a case study

BackgroundExpressed Sequence Tags (ESTs) are a source of simple sequence repeats (SSRs) that can be used to develop molecular markers for genetic studies. The availability of ESTs for Quercus robur and Quercus petraea provided a unique opportunity to develop microsatellite markers to accelerate research aimed at studying adaptation of these long-lived species to their environment. As a first step toward the construction of a SSR-based linkage map of oak for quantitative trait locus (QTL) mapping, we describe the mining and survey of EST-SSRs as well as a fast and cost-effective approach (bin mapping) to assign these markers to an approximate map position. We also compared the level of polymorphism between genomic and EST-derived SSRs and address the transferability of EST-SSRs in Castanea sativa (chestnut).ResultsA catalogue of 103,000 Sanger ESTs was assembled into 28,024 unigenes from which 18.6% presented one or more SSR motifs. More than 42% of these SSRs corresponded to trinucleotides. Primer pairs were designed for 748 putative unigenes. Overall 37.7% (283) were found to amplify a single polymorphic locus in a reference full-sib pedigree of Quercus robur. The usefulness of these loci for establishing a genetic map was assessed using a bin mapping approach. Bin maps were constructed for the male and female parental tree for which framework linkage maps based on AFLP markers were available. The bin set consisting of 14 highly informative offspring selected based on the number and position of crossover sites. The female and male maps comprised 44 and 37 bins, with an average bin length of 16.5 cM and 20.99 cM, respectively. A total of 256 EST-SSRs were assigned to bins and their map position was further validated by linkage mapping. EST-SSRs were found to be less polymorphic than genomic SSRs, but their transferability rate to chestnut, a phylogenetically related species to oak, was higher.ConclusionWe have generated a bin map for oak comprising 256 EST-SSRs. This resource constitutes a first step toward the establishment of a gene-based map for this genus that will facilitate the dissection of QTLs affecting complex traits of ecological importance.

Molecular footprints of local adaptation in two Mediterranean conifers

This study combines neutrality tests and environmental correlations to identify nonneutral patterns of evolution in candidate genes related to drought stress in two closely related Mediterranean conifers, Pinus pinaster Ait. and P. halepensis Mill. Based on previous studies, we selected twelve amplicons covering six candidate genes that were sequenced in a large sample spanning the full range of these two species. Neutrality tests relatively robust to demography (DHEW compound test and maximum likelihood multilocus Hudson-Kreitman-Aguadé test) were used to detect selection events at different temporal scales. Environmental associations between variation at candidate genes and climatic variables were also examined. These combined approaches detected distinct genes that may be targeted by selection, most of them specific to only one of the two conifers, despite their recent divergence (<10 Ma). An exception was 4-coumarate: CoA ligase, a gene involved in the production of various important secondary products that appeared to play a role in local adaptation processes of both pines. Another remarkable result was that all significant environmental correlations involved temperature indices, highlighting the importance of this climatic factor as a selective driver on Mediterranean pines. The ability to detect natural selection at the DNA sequence level depends on the nature and the strength of the selection events, on the timescale at which they occurred, and on the sensitivity of the methods to other evolutionary forces that can mimic selection (e.g., demography and population structure). Using complementary approaches can help to capture different aspects of the evolutionary processes that govern molecular variation at both intra- and interspecific levels.

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 (pk > 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.

Genetic effects of chronic habitat fragmentation revisited: Strong genetic structure in a temperate tree, Taxus baccata (Taxaceae), with great dispersal capability

Tree species are thought to be relatively resistant to habitat fragmentation because of their longevity and their aptitude for extensive gene flow, although recent empirical studies have reported negative genetic consequences, in particular after long-term habitat fragmentation in European temperate regions. Yet the response of each species to habitat loss may differ greatly depending on their biological attributes, in particular seed dispersal ability. In this study, we used demographic and molecular data to investigate the genetic consequences of chronic habitat fragmentation in remnant populations of Taxus baccata in the Montseny Mountains, northeast Spain. The age structure of populations revealed demographic bottlenecks and recruitment events associated with exploitation and management practices. We found a strong genetic structure, both at the landscape and within-population levels. We also detected high levels of inbreeding for a strictly outcrossing species. Chronic forest fragmentation resulting from long-term exploitation in the Montseny Mountains seems the most plausible explanation for the strong genetic structure observed. Our results support the view that, contrary to some predictions, tree species are not buffered from the adverse effects of habitat fragmentation, even in the case of species with a high dispersal potential.

Patterns of polymorphism resulting from long-range colonization in the Mediterranean conifer Aleppo pine.

*The evolutionary outcomes of range expansion/contraction depend on the biological system considered and the interactions among the evolutionary forces in place. In this study, we examined the demographic history and the local polymorphism patterns of candidate genes linked to drought tolerance of a widespread Mediterranean conifer (Pinus halepensis). *To that end, we used cpSSRs and coalescence modelling of nuclear genes to infer the demographic history of natural populations covering the species range. Ten drought-response candidate genes were then examined for their patterns of polymorphism and tested for selection considering plausible demographic scenarios. *Our results revealed a marked loss of genetic diversity from the relictual Greek population towards the western range of the species, as well as molecular signatures of intense bottlenecks. Moreover, we found an excess of derived polymorphisms in several genes sampled in the western part of the range - a potential result of the action of natural selection on populations confronted with new environments following long-range colonization. *Wide-range expansions-contractions of forest trees are accompanied by strong selective pressures, resulting in distinct evolutionary units. This knowledge is of crucial importance for the conservation and management of forests in the face of climate change.

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

English yew (Taxus baccata L., Taxaceae), a Tertiary relict, provides a seminal example of a widespread albeit locally endangered (often close to extinction) tree species. In order to gain detailed insights into the evolutionary dynamics of the species on a broad geographical scale, over 1000 trees from 91 populations of English yew in the western Mediterranean were analyzed using seven nuclear microsatellite markers. Our results revealed contrasting patterns of genetic structure at different spatial scales: genetic variation was highly structured at the local scale, while only a low proportion of the observed variation was attributed to regional differences. We also found a geographic gradient of decreasing diversity and increasing population divergence from northwest (central Europe and northern Iberian Peninsula) to southeast (Mediterranean Iberia and North Africa). The patterns revealed in this study probably reflect the combined effects of Quaternary climatic changes and recent impact of human activities, and potentially also more ancient events dating back to the Tertiary. Both climatic and anthropogenic factors seem to have conducted to a long history of population isolation, which may have contributed significantly to enhance population divergence through restricted gene flow and genetic drift.