Antoine Kremer

Identification of refugia and post-glacial colonisation routes of European white oaks based on chloroplast DNA and fossil pollen evidence

Abstract The geographic distribution throughout Europe of each of 32 chloroplast DNA variants belonging to eight white oak species sampled from 2613 populations is presented. Clear-cut geographic patterns were revealed by the survey. These distributions, together with the available palynological information, were used to infer colonisation routes out of the glacial period refugia. In western Europe in particular, movements out of the Iberian and the Italian Peninsulas can be clearly identified. Separate refugia are also present in eastern Balkans, whereas further west in this peninsula similarities with Italy were evident. Movements resulting in the exchange of haplotypes between refugia both during the present interglacial and probably also during earlier glacial cycles were therefore inferred. The consequences of these past exchanges is that phylogenetically divergent haplotypes have sometimes followed very similar colonisation routes, limiting somewhat the phylogeographic structure. Cases of geographic disjunction in the present-day distribution of haplotypes are also apparent and could have been induced by the existence of rapid climatic changes at the end of the glacial period (specifically the Younger Dryas cold period), which resulted in range restriction following an early warm period during which oak first expanded from its primary refugia. This cold phase was followed by a new period of expansion at the outset of the Holocene, involving in some cases ‘secondary’ refugia. It is expected that these short climate oscillations would have led to a partial reshuffling of haplotype distribution. Early association between haplotypes and oak species are also suggested by the data, although extensive introgression among species has ultimately largely blurred the pattern. This implies that colonisation routes may have been initially constrained by the ecological characteristics of the species hosting each chloroplast variant. We suggest for instance that two oak species distributed in the north of the Iberian Peninsula ( Quercus petraea and Q. pubescens ) are recent post-glacial immigrants there. When considered together, conclusions on the location of glacial period refugia and the colonisation routes derived from molecular information and fossil pollen data appear to be both largely compatible and complementary.

Chloroplast DNA variation in European white oaks: Phylogeography and patterns of diversity based on data from over 2600 populations

A consortium of 16 laboratories have studied chloroplast DNA (cpDNA) variation in European white oaks. A common strategy for molecular screening, based on restriction analysis of four PCR-amplified cpDNA fragments, was used to allow comparison among the different laboratories. A total of 2613 oak populations (12,214 individual trees from eight species) were sampled from 37 countries, and analysed with the four fragments. They belong to eight related oak species: Quercus robur, Q. petraea, Q. pubescens, Q. frainetto, Q. faginea, Q. pyrenaica, Q. canariensis and Q. macranthera. During this survey, 45 chloroplast variants were detected and are described together with their phylogenetic relationships, but several of these haplotypes were pooled when there were some risks of confusion across laboratories during the survey, and finally 32 remained that were mapped and used in diversity analyses. A strong phylogeographic structure is apparent from the data, where related haplotypes have broadly similar geographic distributions. In total, six cpDNA lineages are identified, which have distinct geographic distributions, mainly along a longitudinal gradient. Most haplotypes found in northern Europe are also present in the south, whereas the converse is not true, suggesting that the majority of mutations observed were generated prior to postglacial recolonisation, corroborating the conclusions of earlier studies. The description of a new western European lineage constitutes a major finding, compared to earlier phylogenetic treatments. Although the eight oak species studied systematically share cpDNA variants when in sympatry, they partition cpDNA diversity differently, as a consequence of their different ecology and life history attributes. Regional differences in levels of differentiation also exist (either species-specific or general); these seem to be related to the intensity of past and present management of the forests across Europe but also to the level of fragmentation of the range within these regions.

Forest tree genomics: growing resources and applications

Over the past two decades, research in forest tree genomics has lagged behind that of model and agricultural systems. However, genomic research in forest trees is poised to enter into an important and productive phase owing to the advent of next-generation sequencing technologies, the enormous genetic diversity in forest trees and the need to mitigate the effects of climate change. Research on long-lived woody perennials is extending our molecular knowledge of complex life histories and adaptations to the environment — enriching a field that has traditionally drawn biological inference from a few short-lived herbaceous species.

Potential for evolutionary responses to climate change – evidence from tree populations

Evolutionary responses are required for tree populations to be able to track climate change. Results of 250 years of common garden experiments show that most forest trees have evolved local adaptation, as evidenced by the adaptive differentiation of populations in quantitative traits, reflecting environmental conditions of population origins. On the basis of the patterns of quantitative variation for 19 adaptation-related traits studied in 59 tree species (mostly temperate and boreal species from the Northern hemisphere), we found that genetic differentiation between populations and clinal variation along environmental gradients were very common (respectively, 90% and 78% of cases). Thus, responding to climate change will likely require that the quantitative traits of populations again match their environments. We examine what kind of information is needed for evaluating the potential to respond, and what information is already available. We review the genetic models related to selection responses, and what is known currently about the genetic basis of the traits. We address special problems to be found at the range margins, and highlight the need for more modeling to understand specific issues at southern and northern margins. We need new common garden experiments for less known species. For extensively studied species, new experiments are needed outside the current ranges. Improving genomic information will allow better prediction of responses. Competitive and other interactions within species and interactions between species deserve more consideration. Despite the long generation times, the strong background in quantitative genetics and growing genomic resources make forest trees useful species for climate change research. The greatest adaptive response is expected when populations are large, have high genetic variability, selection is strong, and there is ecological opportunity for establishment of better adapted genotypes.

Within-population genetic structure in Quercus robur L. and Quercus petraea (Matt.) Liebl. assessed with isozymes and microsatellites

The spatial distribution of alleles is described in a forest stand of natural origin of 5 ha comprising 355 mature Quercus petraea and Q. robur trees. Each tree was genotyped for six microsatellite loci. Previous studies on the same population based on isozymes allowed a comparison of different markers for the detection of spatial genetic structure. Different statistics were used: differentiation measures at different spatial scales, and spatial autocorrelation analysis based on Moran’s index I. For microsatellites, differentiation and autocorrelation were calculated with unordered alleles (identity in state) and with alleles ordered according to their size. Results showed the same tendency of a significant, but low, spatial genetic structure for markers and different statistics. Some differences could, however, be detected. First, microsatellites interpreted as unordered alleles exhibit stronger spatial structure than isozymes or microsatellite interpreted as ordered alleles. Second, differentiation and autocorrelation values were higher in Q. petraea than in Q. robur. These differences were attributed to species differences in gene flow via pollen or seed.

Long‐distance gene flow and adaptation of forest trees to rapid climate change

Forest trees are the dominant species in many parts of the world and predicting how they might respond to climate change is a vital global concern. Trees are capable of long-distance gene flow, which can promote adaptive evolution in novel environments by increasing genetic variation for fitness. It is unclear, however, if this can compensate for maladaptive effects of gene flow and for the long-generation times of trees. We critically review data on the extent of long-distance gene flow and summarise theory that allows us to predict evolutionary responses of trees to climate change. Estimates of long-distance gene flow based both on direct observations and on genetic methods provide evidence that genes can move over spatial scales larger than habitat shifts predicted under climate change within one generation. Both theoretical and empirical data suggest that the positive effects of gene flow on adaptation may dominate in many instances. The balance of positive to negative consequences of gene flow may, however, differ for leading edge, core and rear sections of forest distributions. We propose future experimental and theoretical research that would better integrate dispersal biology with evolutionary quantitative genetics and improve predictions of tree responses to climate change.

Fine‐scale genetic structure and gene dispersal inferences in 10 Neotropical tree species

The extent of gene dispersal is a fundamental factor of the population and evolutionary dynamics of tropical tree species, but directly monitoring seed and pollen movement is a difficult task. However, indirect estimates of historical gene dispersal can be obtained from the fine‐scale spatial genetic structure of populations at drift–dispersal equilibrium. Using an approach that is based on the slope of the regression of pairwise kinship coefficients on spatial distance and estimates of the effective population density, we compare indirect gene dispersal estimates of sympatric populations of 10 tropical tree species. We re‐analysed 26 data sets consisting of mapped allozyme, SSR (simple sequence repeat), RAPD (random amplified polymorphic DNA) or AFLP (amplified fragment length polymorphism) genotypes from two rainforest sites in French Guiana. Gene dispersal estimates were obtained for at least one marker in each species, although the estimation procedure failed under insufficient marker polymorphism, limited sample size, or inappropriate sampling area. Estimates generally suffered low precision and were affected by assumptions regarding the effective population density. Averaging estimates over data sets, the extent of gene dispersal ranged from 150 m to 1200 m according to species. Smaller gene dispersal estimates were obtained in species with heavy diaspores, which are presumably not well dispersed, and in populations with high local adult density. We suggest that limited seed dispersal could indirectly limit effective pollen dispersal by creating higher local tree densities, thereby increasing the positive correlation between pollen and seed dispersal distances. We discuss the potential and limitations of our indirect estimation procedure and suggest guidelines for future studies.

Comparison of microsatellites and amplified fragment length polymorphism markers for parentage analysis

This study compares the properties of dominant markers, such as amplified fragment length polymorphisms (AFLPs), with those of codominant multiallelic markers, such as microsatellites, in reconstructing parentage. These two types of markers were used to search for both parents of an individual without prior knowledge of their relationships, by calculating likelihood ratios based on genotypic data, including mistyping. Experimental data on 89 oak trees genotyped for six microsatellite markers and 159 polymorphic AFLP loci were used as a starting point for simulations and tests. Both sets of markers produced high exclusion probabilities, and among dominant markers those with dominant allele frequencies in the range 0.1–0.4 were more informative. Such codominant and dominant markers can be used to construct powerful statistical tests to decide whether a genotyped individual (or two individuals) can be considered as the true parent (or parent pair). Gene flow from outside the study stand (GFO), inferred from parentage analysis with microsatellites, overestimated the true GFO, whereas with AFLPs it was underestimated. As expected, dominant markers are less efficient than codominant markers for achieving this, but can still be used with good confidence, especially when loci are deliberately selected according to their allele frequencies.

IDENTIFICATION AND CHARACTERIZATION OF (GA/CT)N- MICROSATELLITE LOCI FROM QUERCUS PETRAEA

In this study a size selected genomic library from Quercus petraea was screened for (GA/CT)n-microsatellite sequences. The resulting loci were analysed by PCR for their usefulness as molecular markers in Q. petraea and Q. robur. 17 out of 52 tested primer pairs resulted in the amplification of a polymorphic single-locus pattern. The number of alleles found per locus varied from 6 to 16. Combining the genetic variation observed for the characterized loci provides a unique genotype for all the individuals tested. Using intraspecific controlled crosses of Q. robur trees Mendelian inheritance could be shown for five loci.

Species relative abundance and direction of introgression in oaks

Successful hybridisation and subsequent introgression lead to the transfer of genetic material across species boundaries. In this process, species relative abundance can play a significant role. If one species is less abundant than the other, its females will receive many heterospecific gametes, increasing mate‐recognition errors and thus hybridisation rate. Moreover, first‐generation hybrids will also more likely mate with the more abundant species, leading to asymmetric introgression. These predictions have important fundamental consequences, especially during biological invasions or when a rare species threatened by extinction is surrounded by individuals from a related species. However, experimental tests in nature of the importance of the relative abundance of each species on hybridisation dynamics remain scarce. We assess here the impact of species relative abundance on hybridisation dynamics among four species from the European white oak species complex. A total of 2107 oak trees were genotyped at 10 microsatellite markers and Bayesian clustering methods were used to identify reference trees of each species. We then used these reference trees to simulate purebred and hybrid genotypes to determine optimal threshold for genetic assignment. With this approach, we found widespread evidence of hybridisation between all studied oak species, with high occurrence of hybrids, varying from 11% to 31% according to stand and sampling strategies. This finding suggests that hybridisation is a common phenomenon that plays a significant role in evolution of this oak species complex. In addition, we demonstrate a strong impact of species abundance on both hybridisation rate and introgression directionality.