Outlier loci highlight the direction of introgression in oaks.

Loci considered to be under selection are generally avoided in attempts to infer past demographic processes as they do not fit neutral model assumptions. However, opportunities to better reconstruct some aspects of past demography might thus be missed. Here we examined genetic differentiation between two sympatric European oak species with contrasting ecological dynamics (Quercus robur and Quercus petraea) with both outlier (i.e. loci possibly affected by divergent selection between species or by hitchhiking effects with genomic regions under selection) and nonoutlier loci. We sampled 855 individuals in six mixed forests in France and genotyped them with a set of 262 SNPs enriched with markers showing high interspecific differentiation, resulting in accurate species delimitation. We identified between 13 and 74 interspecific outlier loci, depending on the coalescent simulation models and parameters used. Greater genetic diversity was predicted in Q. petraea (a late-successional species) than in Q. robur (an early successional species) as introgression should theoretically occur predominantly from the resident species to the invading species. Remarkably, this prediction was verified with outlier loci but not with nonoutlier loci. We suggest that the lower effective interspecific gene flow at loci showing high interspecific divergence has better preserved the signal of past asymmetric introgression towards Q. petraea caused by the species' contrasting dynamics. Using markers under selection to reconstruct past demographic processes could therefore have broader potential than generally recognized.

Detecting the footprints of divergent selection in oaks with linked markers.

Genome scans are increasingly used to study ecological speciation, providing a useful genome-wide perspective on divergent selection in the presence of gene flow. Here, we compare current approaches to detect footprints of divergent selection in closely related species. We analyzed 192 individuals from two interfertile European temperate oak species using 30 nuclear microsatellites from eight linkage groups. These markers present little intraspecific differentiation and can be used in combination to assign individual genotypes to species. We first show that different outlier detection tests give somewhat different results, possibly due to model constraints. Second, using linkage information for these markers, we further characterize the signature of divergent selection in the presence of gene flow. In particular, we show that recombination estimates for regions with outlier markers are lower than those for a control region, in line with a prediction from ecological speciation theory. Most importantly, we show that analyses at the haplotype level can distinguish between truly divergent (bi-directional) selection and positive selection in one of the two species, offering a new and improved method for characterizing the speciation process.

Multiplexed microsatellite markers for genetic studies of beech.

European beech (Fagus sylvatica L.) is one of the economically most important broadleaved tree species in Europe and has become a model for studying climate change effects on forests. Multiplex PCR of microsatellites is a fast and cost-effective technique allowing high-throughput genotyping. Here we present the procedure used to develop two multiplex kits (8-plexes) for European beech. We paid particular attention to quality control throughout all steps of the multiplex kits development (null allele detection, error rate measurements, linkage disequilibrium). Preliminary assays suggest that the 16 amplified loci are largely devoid of null alleles and allow rapid and cost-effective genotyping of beech with low error rates. The two kits, which differ in their levels of polymorphism, most likely due to marker origin, were also informative in seven other beech species tested.

Current trends in microsatellite genotyping.

Microsatellites have been popular molecular markers ever since their advent in the late eighties. Despite growing competition from new genotyping and sequencing techniques, the use of these versatile and cost-effective markers continues to increase, boosted by successive technical advances. First, methods for multiplexing PCR have considerably improved over the last years, thereby decreasing genotyping costs and increasing throughput. Second, next-generation sequencing technologies allow the identification of large numbers of microsatellite loci at reduced cost in non-model species. As a consequence, more stringent selection of loci is possible, thereby further enhancing multiplex quality and efficiency. However, current practices are lagging behind. By surveying recently published population genetic studies relying on simple sequence repeats, we show that more than half of the studies lack appropriate quality controls and do not make use of multiplex PCR. To make the most of the latest technical developments, we outline the need for a well-established strategy including standardized high-throughput bench protocols and specific bioinformatic tools, from primer design to allele calling.

Two highly validated multiplexes (12-plex and 8-plex) for species delimitation and parentage analysis in oaks (Quercus spp.).

Multiplex PCR is a fast and cost-effective technique allowing increased genotyping throughput of microsatellites. We developed two multiplexes for Quercus petraea and Q. robur, a 12-plex of EST-SSRs (eSSRs) and an 8-plex of genomic SSRs (gSSRs). We studied the origin of allele calling errors at the human reader and software levels. We showed that the robustness of allele identification can be improved by binning on raw peak sizes prior to genetic data analysis. We checked through simulation the power of these markers for species delimitation and hybrid detection. The resolution achieved with all 20 markers was greatly improved compared to that of previous studies based on a subset of the markers. Preliminary PCR tests suggest that these multiplexes might be useful to study other oak species as well. The strategy used for multiplex microsatellite development (from PCR conditions to the definition of allele calling rules) should be broadly applicable.

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.

Detection of hybrids in nature: application to oaks (Quercus suber and Q. ilex).

Powerful and accurate detection of first-generation (F1) hybrids and backcrosses in nature is needed to achieve a better understanding of the function and dynamics of introgression. To document the frequency of ongoing interspecific gene exchange between two Mediterranean evergreen oaks, the cork oak (Quercus suber) and the holm oak (Q. ilex), we analyzed 1487 individuals originating from across the range of the two species using eight microsatellite loci and two Bayesian clustering approaches (implemented in the programs STRUCTURE and NEWHYBRIDS). Simulated data were used to assess the differences between the two clustering methods and to back up the choice of the threshold value for the posterior probability to discriminate admixed from pure individuals. We found that the use of STRUCTURE resulted in the highest power to detect hybrids, whereas NEWHYBRIDS provided the highest accuracy. Irrespective of the approach, the two species were clearly distinguished as independent genetic entities without any prior information. In contrast with previous reports, we found no evidence for unidirectional introgression. The overall hybridization rate was very low (<2% of introgressed individuals). Only two individuals were identified as F1 hybrids and five as early backcrosses. This work shows that the combined application of the two complementary Bayesian approaches and their systematic validation with simulations, fit for the case at hand, helps gain resolution in the identification of admixed individuals.

Contrasting effects of long distance seed dispersal on genetic diversity during range expansion.

Currently many attempts are made to reconstruct the colonization history of plant species after the last ice age. A surprising finding is that during the colonization phase genetic diversity did not decrease as much as expected. In this paper we examine whether long distance seed dispersal events could play a role in the unexpected maintenance of genetic diversity during range expansion. This study is based on simulations carried out with a maternally inherited haploid locus using a cellular automaton. The simulations reveal a close relationship between the frequency of long distance seed dispersal events and the amount of genetic diversity preserved during colonization. In particular, when the colonized region is narrow, a complete loss of genetic diversity results from the occurrence of very rare long distance dispersal (LDD) events. We call this phenomenon the 'embolism effect'. However, slightly higher rates of LDD events reverse this effect, up to the point that diversity is better preserved than in a pure diffusion model. This phenomenon is linked to the reorganization of the genetic structure during colonization and is called the 'reshuffling effect'.

Chloroplast DNA variation of Quercus rubra L. in North America and comparison with other Fagaceae.

Quercus rubra is one of the most important timber and ornamental tree species from eastern North America. It is a widespread species growing under variable ecological conditions. Chloroplast DNA variation was studied by PCR-RFLP (polymerase chain reaction-restriction fragment length polymorphism) in 290 individuals from 66 populations sampled throughout the natural range. A total of 12 haplotypes were detected, with one found in 75% of the trees. Population differentiation is relatively low (G(ST) = 0.46), even when similarities between haplotypes are taken into account (N(ST) = 0.50), pointing to a weak phylogeographical structure. Furthermore, no spatial structure of genetic diversity could be detected. The genetic differentiation increased northwards, reflecting the postglacial history of Q. rubra. The unusual aspect of this study was the low level of chloroplast DNA genetic differentiation in Q. rubra compared to that typically observed in other oak species. Palynological evidence indicates that during the last glacial maximum, Q. rubra had one major distribution range with populations located relatively far to the north, resulting in only modest movement northwards when climate improved, whereas European white oaks were largely restricted to the southern European peninsulas and experienced extensive movements during the postglacial period. The contrasted geographical features and levels of tree species richness of both continents might further explain why congeneric species sharing similar life history traits have genetic structures that are so different.

Rangewide phylogeography of a bird-dispersed Eurasian shrub: contrasting Mediterranean and temperate glacial refugia.

We studied the phylogeography of alder buckthorn (Frangula alnus), a bird-dispersed shrub or small tree distributed over most of Europe and West Asia and present in three of the four main refugia of West Palaearctic temperate woody plants: the Iberian Peninsula, the Balkans and Anatolia. A total of 78 populations from 21 countries were analysed for chloroplast DNA variation using polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP), and 21 different haplotypes were distinguished. We found a very strong overall population differentiation (GST = 0.81) and phylogeographical structure, and a sharp contrast between the haplotype-rich refugia and the almost completely uniform area of postglacial colonization. The haplotype network comprises three lineages made up of haplotypes from the Iberian Peninsula, Anatolia with the Caucasus, and temperate Europe. The Iberian and the Anatolian branches represent parts of a major lineage that spans over the whole northern Mediterranean Basin and some neighbouring areas and probably dates back to the Tertiary. Many haplotypes of this lineage are distributed locally and most populations are fixed for a single haplotype; these populations have apparently been very stable since their establishment, experiencing negligible gene flow and few mutations. The temperate European lineage consists of one very widespread and abundant plus six locally distributed haplotypes. Four of them are located in Southeast Europe, the putative refugium of all extant temperate European populations. Contrary to populations from Iberia and Anatolia, F. alnus populations from the southeastern European refugium have most genetic variation within populations. Bird-mediated seed dispersal has apparently allowed not only a very rapid postglacial expansion of F. alnus but also subsequent regular seed exchanges between populations of the largely continuous species range in temperate Europe. In contrast, the disjunct F. alnus populations persisting in Mediterranean mountain ranges seem to have experienced little gene flow and have therefore accumulated a high degree of differentiation, even at short distances. Populations from the southern parts of the glacial refugia have contributed little to the postglacial recolonization of Europe, but their long-term historical continuity has allowed them to maintain a unique store of genetic variation.

Phylogeography of maritime pine inferred with organelle markers having contrasted inheritance.

Range-wide variation of maritime pine was studied at maternally inherited and paternally inherited markers (mitochondrial DNA and chloroplast DNA). While chloroplast DNA exhibits the highest diversity, phylogeographic inferences from this marker are blurred by homoplasy and extensive pollen flow. In contrast, the only three mitochondrial haplotypes found provide a clear picture of nonoverlapping areas colonized from different refugia, with no single population having a mixed composition (GST = 1). Comparison of the genetic structure inferred from both organelle genomes allows the investigation of differential seed and pollen dispersal, pointing to pollen, but not seed, dispersal across the Strait of Gibraltar (from Morocco into Iberia). A comparison with already available genetic information, especially that of one of the maritime pine's most threatening insect pests, the bast scale Matsucoccus feytaudi, further completes the picture.

Checking the geographical origin of oak wood: molecular and statistical tools.

New methods for better identification of timber geographical origin would constitute an important technical element in the forest industry, for phytosanitary certification procedures or in the chain of custody developed for the certification of timber from sustainably managed forests. In the case of the European white oaks, a detailed reference map of chloroplast (cp) DNA variation across the range exists, and we propose here to use the strong geographical structure, characterized by a differentiation of western vs. eastern populations, for the purpose of oak wood traceability. We first developed cpDNA markers permitting the characterization of haplotype on degraded DNA obtained from wood samples. The techniques were subsequently validated by confirming the full correspondence between genotypes obtained from living tissues (buds) and from wood collected from the same individual oak. Finally, a statistical procedure was used to test if the haplotype composition of a lot of wood samples is consistent with its presumed geographical origin. Clearly, the technique cannot permit the unambiguous identification of wood products of unknown origin but can be used to check the conformity of genetic composition of wood samples with the region of alleged origin. This could lead to major applications not only in the forest industry but also in archaeology or in palaeobotany.

A case of chloroplast heteroplasmy in kiwifruit (Actinidia deliciosa) that is not transmitted during sexual reproduction.

We report the first case of plastid chimera within the Actinidia genus, where plastid inheritance was believed to be paternal. The heterogeneity of chloroplast DNA observed in the hexaploid Actinidia deliciosa cultivar D uno involves the presence or absence of a particular MspI restriction site in the region between the psbC gene and the tRNA-Ser(UGA) gene. The heterogeneity was first observed using restriction fragment length polymorphism and then confirmed through cloning and sequencing. The analysis of the cloned fragments revealed the presence of two haplotypes: the most frequent type was found in 123 (88.5%) out of a total of 139 colonies screened. Partial sequences of the psbC-trnS fragment from both haplotypes revealed that the polymorphism occurs within the coding region of the psbC gene and consists of a synonymous transition. A contamination-free cross involving D uno as the male parent produced only plants characterized by the most frequent haplotype, indicating either selection bias against the rare type or more likely fixation of the frequent type in tissues leading to the formation of the male gametes. The MspI restriction profiles performed on various tissues suggest that the rarer type is absent from the histogenic layer LII and that D uno is a periclinal plastid chimera.

Phylogeography of the common ivy (Hedera sp.) in Europe: genetic differentiation through space and time.

We studied the phylogeography of ivy (Hedera sp.), a liana widespread in Europe, throughout its natural range. The populations sampled belong to four closely related species differing by ploidy levels and morphological characters. Chloroplast (cp) markers were used and 13 haplotypes were detected, usually shared across species, contrary to ribosomal internal transcribed spacer (ITS) variants. We demonstrated the existence of a strong overall cpDNA phylogeographical structure. Several methods of data analysis were conducted to describe how this structure and the genetic diversity change through space and time. Southern populations, especially those from Spain, are the most divergent. Pairwise estimates of differentiation point to isolation by distance, and the existence of a latitudinal gradient of divergence was demonstrated using a regression procedure. Similarly, latitudinal differences in haplotype richness and diversity exist, as shown by population permutations ('differentiation through space'). Finally, we measured differentiation by taking into account successive levels of divergence between haplotypes ('differentiation through time'). Genetic differentiation turns out to be much greater when differences between closely related haplotypes are not considered. Further, these results indicate that the phylogeographical structure is essentially due to the relative distribution of the most similar haplotypes. Diversity decreases from south to north, whereas haplotype frequencies change longitudinally. It appears that Hedera survived in Spanish and Balkan refugia during the last ice age. A third refugium must have been present in the Alps or in Italy. During the northward expansion, the decrease in overall diversity was attenuated by some mixing of lineages at intermediate latitudes, resulting in comparatively higher levels of differentiation in the south.

Frequent cytoplasmic exchanges between oak species that are not closely related: Quercus suber and Q. ilex in Morocco.

Chloroplast (cp) and mitochondrial (mt) DNA variation were studied in 97 populations of cork oak (Quercus suber) in Morocco; in 31 of these populations, holm oak (Quercus ilex), a clearly distinct species, also occurred and was compared with Q. suber. Three cpDNA and one mtDNA primer pairs were used in the survey, each in combination with one restriction enzyme. Six haplotypes belonging to two very divergent lineages were detected; one lineage predominates in each species, and is probably ancestral, as inferred from comparisons with other oak species. In the mixed-species populations, cytoplasmic genomes were frequently shared across species, as indicated by an introgression ratio of 0.63. This index is a new measure of the propensity of species to share locally genetic markers, varying from zero (complete differentiation) to one (no differentiation). By contrast, more closely related deciduous oak species (Q. robur, Q. petraea and Q. pubescens) have introgression ratios varying from 0.82 to 0.97. The introgression events appear to have been more frequent in the direction Q. ilex (female) x Q. suber (male), a finding which seems attributable to the flowering phenology of these two species. This asymmetry may have favoured immigration of Q. suber beyond its main range, in regions already colonized by Q. ilex. There, rare hybridization and further introgression through long distance pollen flow have established populations that are morphologically indistinguishable from Q. suber but that have cytoplasmic genomes originating from the local Q. ilex populations.

Pollen- versus seed-mediated gene flow in a scattered forest tree species.

We examined the spatial distribution of maternally inherited chloroplast DNA markers over the French part of the range of Sorbus torminalis, a scattered temperate forest tree native to most of Europe. The survey by restriction analysis of polymerase-chain-reaction amplified fragments for 880 individuals distributed among 55 populations allowed the detection of 25 haplotypes. The coefficient of differentiation among populations computed on the basis of haplotype frequency (G(STc) = 0.34) was one of the lowest found in forest trees so far, and the mean within-population diversity was relatively high, indicating multiple-mother foundation events. A significant but slight geographical pattern was observed, up to distances of about 100 km. This pattern of differentiation was compared to the genetic structure of the same populations revealed by biparentally inherited markers (isoenzymes), and a new method to quantify the relative importance of seed and pollen dispersal was derived, based on isolation-by-distance models. Neither pollen- nor seed-mediated gene flow was predominant in S. torminalis, a finding that differs from those for the majority of tree species studied so far. This result was most likely due to an extinction-recolonization dynamics based on efficient seed dispersal strategies. The joint screening of 31 individuals of the related Sorbus aria and of 163 hybrid individuals shows that hybridization occurs predominantly in one direction and is rarely followed by cytoplasmic introgression. As a consequence, interspecific gene flow should not significantly affect the diversity dynamics within S. torminalis.

Diverging trends between heterozygosity and allelic richness during postglacial colonization in the European beech.

Variation at 12 polymorphic isozyme loci was studied in the European beech on the basis of an extensive sample of 389 populations distributed throughout the species range. Special emphasis was given to the analysis of the pattern of geographic variation on the basis of two contrasting measures of genetic diversity, gene diversity (H) and allelic richness, and to their relationship. Measures of allelic richness were corrected for variation in sample size by using the rarefaction method. As expected, maximum allelic richness was found in the southeastern part of the range (southern Italy and the Balkans), where beech was confined during the last ice age. Surprisingly, H was lower in refugia than in recently colonized regions, resulting in a negative correlation between the two diversity measures. The decrease of allelic richness and the simultaneous increase of H during postglacial recolonization was attributed to several processes that differentially affect the two diversity parameters, such as bottlenecks due to long-distance founding events, selection during population establishment, and increased gene flow at low population densities.

Spatial and temporal distribution of chloroplast DNA polymorphism in a tropical tree species.

The level and the spatial organization of chloroplast DNA polymorphism were investigated in Dicorynia guianensis Hamshoff (Caesalpiniaceae) at different spatial and temporal scales. D. guianensis is a canopy tree of the rain forest that is distributed throughout the Guiana plateau in small aggregates. Twelve different haplotypes were identified using restriction analysis of polymerase chain reaction (PCR) amplified fragments of the chloroplast genome. When populations from different areas of French Guiana were compared, a clear geographical pattern of haplotype frequencies was identified along the Atlantic coast. This pattern is most likely the result of the restriction-expansion dynamics of the tropical forest during the Quaternary. At the local level, D. guianensis was characterized by a high level of within population diversity. Maintenance of within population diversity results from the dynamics of the aggregates; stochastic demography associated with the turnover of aggregates generates genetic differentiation among them. At the stand level, a strong spatial aggregation of haplotypes persisted from the adult to the seedling cohort indicating limited seed flow. There was also a strong difference in levels of diversity between the cohorts which suggested that recruitment over several years is needed in order to maintain genetic diversity during regeneration.