Harmful self-pollination drives gynodioecy in European chestnut, a self-incompatible tree

PREMISE

Gynodioecy is a rare sexual system in which two genders (sensu Lloyd, 1980), cosexuals and females, coexist. To survive, female plants must compensate for their lack of siring capacity and male attractiveness. In European chestnut (Castanea sativa), an outcrossing tree, self-pollination reduces fruit set in cosexual individuals because of late-acting self-incompatibility and early inbreeding depression. Could this negative sexual interaction explain the presence of females in this species?

METHODS

We studied gender variation in wild populations of European chestnut. In addition, we compared fruit set (the proportion of flowers giving fruits) and other key female fitness components as well as reproductive allocation between genders. We then performed emasculation experiments in cosexual trees, by removing nectar-producing fertile male inflorescences. We also removed sterile but nectar-producing male inflorescences from female trees, as a control.

RESULTS

We found a highly variable proportion of male-sterile individuals in the wild in European chestnut. In the experimental plot, trees from each gender had similar size, flower density, and burr set, but different fruit set. Removing nectar-producing male inflorescences from branches or entire trees increased fruit set in cosexual but not in female trees.

CONCLUSIONS

These results show that self-pollination impairs fruit set in cosexual trees. Female trees avoid these problems as they do not produce pollen but continue to attract pollinators thanks to their rewarding male-sterile inflorescences, resulting in a much higher fruit set than in cosexuals. This demonstrates that even outcrossed plants can benefit from the cessation of self-pollination, to the point that unisexuality can evolve.

Adaptive function of duodichogamy: Why do chestnut trees have two pollen emission phases?

PREMISE

Intersexual mating facilitation in flowering plants has been largely underexplored. Duodichogamy is a rare flowering system in which individual plants flower in the sequence male-female-male. We studied the adaptive advantages of this flowering system using chestnuts (Castanea spp., Fagaceae) as models. These insect-pollinated trees produce many unisexual male catkins responsible for a first staminate phase and a few bisexual catkins responsible for a second staminate phase. We hypothesized that duodichogamy increases female mating success by facilitating pollen deposition on stigmas of the rewardless female flowers through their proximity with attractive male flowers responsible for the minor staminate phase.

METHODS

We monitored insect visits to 11 chestnut trees during the entire flowering period and explored reproductive traits of all known duodichogamous species using published evidence.

RESULTS

In chestnuts, insects visited trees more frequently during the first staminate phase but visited female flowers more frequently during the second staminate phase. All 21 animal-pollinated duodichogamous species identified are mass-flowering woody plants at high risk of self-pollination. In 20 of 21 cases, gynoecia (female flower parts) are located close to androecia (male flower parts), typically those responsible for the second minor staminate phase, whereas androecia are often distant from gynoecia.

CONCLUSIONS

Our results suggest that duodichogamy increases female mating success by facilitating pollen deposition on stigmas by means of the attractiveness of the associated male flowers while effectively limiting self-pollination.

Sexual interference revealed by joint study of male and female pollination success in chestnut

Most seed plants produce both pollen and ovules. In principle, pollen export could interfere with pollen import through self-pollination, resulting in ovule usurpation and reduced fruit set. Evidence for such interference exists under experimental settings but its importance under natural conditions is unknown. To test for sexual interference in nature, it is necessary to study together mating system, through paternity analyses, and fruit set, the proportion of flowers giving seeds or fruits. We developed a new model combining both processes, using chestnut (Castanea) as case study. We carried out a paternity analysis in an intensively studied plot of 273 trees belonging to three interfertile chestnut species and including a range of individuals with more or less functional stamens, resulting in a large data set of 1924 mating events. We then measured fruit set on 216 of these trees. Fruit set of male-fertile trees was much lower than that of male-sterile trees. Our process-based model shows that pollen is not limiting in the study site and hence cannot account for reduced fruit set. It also indicates that self-pollination is high (74%) but selfing rate is low (4%). Self-pollen is less competitive than cross-pollen, reducing sexual interference, but not sufficiently, as many ovules end up being self-fertilized, 95% of which abort before fruit formation, resulting in the loss of 46% of the fruit crop. These results suggest that the main cause of reduced reproductive potential in chestnut is sexual interference by self-pollen, raising questions on its evolutionary origins.

Asymmetric character displacement in mixed oak stands

Ecological character displacement (ECD) refers to a pattern of increased divergence at sites where species ranges overlap caused by competition for resources. Although ECD is believed to be common, there are few in-depth studies that clearly establish its existence, especially in plants. Thus, we have compared leaf traits in allopatric and sympatric populations of two East Asian deciduous oaks: Quercus dentata and Quercus aliena. In contrast to previous studies, we define sympatry and allopatry at a local scale, thereby comparing populations that can or cannot directly interact. Using genetic markers, we found greater genetic divergence between the two oak species growing in mixed stands and inferred that long-term gene flow has predominantly occurred asymmetrically from the cold-tolerant species (Q. dentata) to the warm-demanding later colonizing species (Q. aliena). Analysis of leaf traits revealed greater divergence in mixed than in pure oak stands. This was mostly due to the later colonizing species being characterized by more resource-conservative traits in the presence of the other species. Controlling for relevant environmental differences did not alter these conclusions. These results suggest that asymmetric trait divergence can take place where species coexist, possibly due to the imbalance in demographic history of species resulting in asymmetric inter-specific selection pressures.

Sixty Years from the First Disease Description, a Novel Badnavirus Associated with Chestnut Mosaic Disease

Although chestnut mosaic disease (ChMD) was described several decades ago, its etiology is still not clear. Using classical approaches and high-throughput sequencing (HTS) techniques, we identified a novel Badnavirus that is a strong etiological candidate for ChMD. Two disease sources from Italy and France were submitted to HTS-based viral indexing. Total RNAs were extracted, ribodepleted, and sequenced on an Illumina NextSeq500 (2 × 150 nt or 2 × 75 nt). In each source, we identified a single contig of ≈7.2 kb that corresponds to a complete circular viral genome and shares homologies with various badnaviruses. The genomes of the two isolates have an average nucleotide identity of 90.5%, with a typical badnaviral genome organization comprising three open reading frames. Phylogenetic analyses and sequence comparisons showed that this virus is a novel species; we propose the name Chestnut mosaic virus (ChMV). Using a newly developed molecular detection test, we systematically detected the virus in symptomatic graft-inoculated indicator plants (chestnut and American oak) as well in chestnut trees presenting typical ChMD symptoms in the field (100 and 87% in France and Italy surveys, respectively). Datamining of publicly available chestnut sequence read archive transcriptomic data allowed the reconstruction of two additional complete ChMV genomes from two Castanea mollissima sources from the United States as well as ChMV detection in C. dentata from the United States. Preliminary epidemiological studies performed in France and central eastern Italy showed that ChMV has a high incidence in some commercial orchards and low within-orchard genetic diversity.

The oak syngameon: more than the sum of its parts

One of Anthropocene's most daunting challenges for conservation biology is habitat extinction, caused by rapid global change. Tree diversity has persisted through previous episodes of rapid change, even global extinctions. Given the pace of current change, our management of extant diversity needs to facilitate and even enhance the natural ability of trees to adapt and diversify. Numerous processes contribute to this evolutionary flexibility, including introgression, a widespread yet under-studied process. Reproductive networks, in which species remain distinct despite interspecific gene flow, are called syngameons, a concept largely inspired from work focusing on Quercus. Delineating and analyzing such species groups, empirically and theoretically, will provide insights into the nonadditive effects on evolution of numerous partially interfertile species exchanging genetic material episodically under changing environmental conditions. To conserve tree diversity, crossing experiments designed with an empirical and theoretical understanding of the constituent syngameon should be set up to assist diversification and adaptation in the Anthropocene. Our increasingly detailed knowledge of the oak genome and of oak interspecific and intraspecific phenotypic variation will improve our ability to sustain the diversity of this tree through an unpredictable and unprecedented future.

Chloroplast DNA variation in a hyperdiverse tropical tree community

We investigate chloroplast DNA variation in a hyperdiverse community of tropical rainforest trees in French Guiana, focusing on patterns of intraspecific and interspecific variation. We test whether a species genetic diversity is higher when it has congeners in the community with which it can exchange genes and if shared haplotypes are more frequent in genetically diverse species, as expected in the presence of introgression.We sampled a total of 1,681 individual trees from 472 species corresponding to 198 genera and sequenced them at a noncoding chloroplast DNA fragment.Polymorphism was more frequent in species that have congeneric species in the study site than in those without congeners (30% vs. 12%). Moreover, more chloroplast haplotypes were shared with congeners in polymorphic species than in monomorphic ones (44% vs. 28%).Despite large heterogeneities caused by genus-specific behaviors in patterns of hybridization, these results suggest that the higher polymorphism in the presence of congeners is caused by local introgression rather than by incomplete lineage sorting. Our findings suggest that introgression has the potential to drive intraspecific genetic diversity in species-rich tropical forests.

Invoking adaptation to decipher the genetic legacy of past climate change

Persistence of natural populations during periods of climate change is likely to depend on migration (range shifts) or adaptation. These responses were traditionally considered discrete processes and conceptually divided into the realms of ecology and evolution. In a milestone paper, Davis and Shaw (2001) Science 292:673 argued that the interplay of adaptation and migration was central to biotic responses to Quaternary climate, but since then there has been no synthesis of efforts made to set up this research program. Here we review some of the salient findings from molecular genetic studies assessing ecological and evolutionary responses to Quaternary climate change. These studies have revolutionized our understanding of population processes associated with past species migration. However, knowledge remains limited about the role of natural selection for local adaptation of populations to Quaternary environmental fluctuations and associated range shifts, and for the footprints this might have left on extant populations. Next-generation sequencing technologies, high-resolution paleoclimate analyses, and advances in population genetic theory offer an unprecedented opportunity to test hypotheses about adaptation through time. Recent population genomics studies have greatly improved our understanding of the role of contemporary adaptation to local environments in shaping spatial patterns of genetic diversity across modern-day landscapes. Advances in this burgeoning field provide important conceptual and methodological bases to decipher the historical role of natural selection and assess adaptation to past environmental variation. We suggest that a process called "temporal conditional neutrality" has taken place: some alleles favored in glacial environments become selectively neutral in modern-day conditions, whereas some alleles that had been neutral during glacial periods become under selection in modern environments. Building on this view, we present a new integrative framework for addressing the interplay of demographic and adaptive evolutionary responses to Quaternary climate dynamics, the research agenda initially envisioned by Davis and Shaw (2001) Science 292:673.

ARE CHLOROPLAST AND MITOCHONDRIAL DNA VARIATION SPECIES INDEPENDENT IN OAKS?

Extensive introgression of cytoplasmic genomes across oak species is now a well-established fact. To distinguish between ancient hybridization events and ongoing introgression, a direct test for the existence of local exchanges is proposed. Such local exchanges must be comparatively recent, that is, contemporaneous with or later than the last postglacial recolonization. The test is applied to an extensive set of data comprising 377 pure or mixed populations (1744 individuals) of four white oak species in southern France. After demonstrating that local exchanges have occurred frequently between all species pairs, another test is performed to check if species status does nevertheless play some role in restricting cytoplasmic gene flow. The results vary according to the species pairs considered, and the observed pattern may be related to the ecology and/or compatibility of interspecific crosses. It is also shown that, for some of these oak species, the presence of related species in a population significantly influences the intraspecific diversity. Altogether, the results demonstrate that (1) intraspecific cytoplasmic gene flow varies according to the species, (2) interspecific cytoplasmic gene flow varies according to the species pair, and (3) both components of gene flow are at least partly related.

Impacts of local adaptation of forest trees on associations with herbivorous insects: implications for adaptive forest management

Disruption of species interactions is a key issue in climate change biology. Interactions involving forest trees may be particularly vulnerable due to evolutionary rate limitations imposed by long generation times. One mitigation strategy for such impacts is Climate matching - the augmentation of local native tree populations by input from nonlocal populations currently experiencing predicted future climates. This strategy is controversial because of potential cascading impacts on locally adapted animal communities. We explored these impacts using abundance data for local native gallwasp herbivores sampled from 20 provenances of sessile oak (Quercus petraea) planted in a common garden trial. We hypothesized that non-native provenances would show (i) declining growth performance with increasing distance between provenance origin and trial site, and (ii) phenological differences to local oaks that increased with latitudinal differences between origin and trial site. Under a local adaptation hypothesis, we predicted declining gallwasp abundance with increasing phenological mismatch between native and climate-matched trees. Both hypotheses for oaks were supported. Provenance explained significant variation in gallwasp abundance, but no gall type showed the relationship between abundance and phenological mismatch predicted by a local adaptation hypothesis. Our results show that climate matching would have complex and variable impacts on oak gall communities.

Within-range translocations and their consequences in European larch

In contrast to biological invasions, translocations of individuals within a species range are understudied, due to difficulties in systematically detecting them. This results in limited knowledge about the corresponding processes and uncertainties regarding the status of extant populations. European larch, a forest tree whose fragmented native distribution is restricted to the Alps and to other Central European mountains, has been massively planted for at least 300 years. Here we focus on the genetic characterization of translocations having taken place within its native range. Microsatellite variation at 13 nuclear loci and sequence data of two mitochondrial DNA fragments were analyzed on the basis of a comprehensive range-wide population sample. Two complementary methods (Geneclass and Structure) were used to infer translocation events based on nuclear data whereas mitochondrial data were used for validation of these inferences. Using Geneclass, we found translocation events in a majority of populations. Additional cases of translocation and many instances of admixture were identified using Structure, thanks to the clear-cut ancestral genetic structure detected in this species. In particular, a strong divide between Alpine and Central European populations, also apparent at mitochondrial markers, helped uncover details on translocation events and related processes. Translocations and associated admixture events were found to be heterogeneously distributed across the species range, with a particularly high frequency in Central Europe. Furthermore, translocations frequently involved multiple geographic sources, some of which were over-represented. Our study illustrates the importance of range-wide investigations for tracing translocations back to their origins and for revealing some of their consequences. It provides some first clues for developing suitable conservation and management strategies.

Cryptic no more: soil macrofossils uncover Pleistocene forest microrefugia within a periglacial desert

Despite their critical importance for understanding the local effects of global climate change on biodiversity, glacial microrefugia are not well studied because they are difficult to detect by using classical palaeoecological or population genetics approaches. We used soil macrofossil charcoal analysis to uncover the presence of cryptic glacial refugia for European beech (Fagus sylvatica) and other tree species in the Landes de Gascogne (southwestern France). Using botanical identification and direct radiocarbon dating (140 (14) C-dates) of macrofossil charcoal extracted from mineral soils, we reconstructed the glacial and postglacial history of all extant beech stands in the region (n = 11). Soil charcoal macrofossils were found in all sites, allowing the identification of up to at least 14 distinct fire events per site. There was direct evidence of the presence of beech during the last glacial period at three sites. Beech was detected during Heinrich stadial-1, one of the coldest and driest intervals of the last glacial period in Western Europe. Together with previous results on the genetic structure of the species in the region, these findings suggest that beech persisted in situ in several microrefugia through full glacial and interglacial periods up to the present day.

Stronger spatial genetic structure in recolonized areas than in refugia in the European beech

Extant rear-edge populations located in former glacial refugia remain understudied despite their high conservation value. These populations should have experienced strong genetic drift due to their small size and long isolation. Moreover, the prolonged action of isolation by distance in refugial areas should result in stronger regional spatial genetic structure (SGS) than in recolonized areas, but empirical tests of this prediction are scarce. To fill this gap, we first used a set of 16 microsatellite markers to investigate the genetic structure of European beech in France in 65 populations from three refugial areas and one control recolonized (nonrefugial) area. Then, using the same approach, we reanalysed published isozyme data from 375 populations distributed across the entire species range. We found stronger genetic differentiation among populations in refugia than in recolonized areas. However, contrary to expectations, regional SGS was lower within refugia than within recolonized areas. Published studies presenting similar analyses suggest that our results could have generality across different biogeographical settings and types of organisms. Strong and prolonged genetic drift in refugial areas could have erased the signature of range expansions that is still visible in recolonized areas. Our results therefore suggest that Pleistocene population isolation has played a key role in increasing the genetic complexity of extant rear-edge populations.

Putting the biological species concept to the test: using mating networks to delimit species

Although interfertility is the key criterion upon which Mayr's biological species concept is based, it has never been applied directly to delimit species under natural conditions. Our study fills this gap. We used the interfertility criterion to delimit two closely related oak species in a forest stand by analyzing the network of natural mating events between individuals. The results reveal two groups of interfertile individuals connected by only few mating events. These two groups were largely congruent with those determined using other criteria (morphological similarity, genotypic similarity and individual relatedness). Our study, therefore, shows that the analysis of mating networks is an effective method to delimit species based on the interfertility criterion, provided that adequate network data can be assembled. Our study also shows that although species boundaries are highly congruent across methods of species delimitation, they are not exactly the same. Most of the differences stem from assignment of individuals to an intermediate category. The discrepancies between methods may reflect a biological reality. Indeed, the interfertility criterion is an environment-dependant criterion as species abundances typically affect rates of hybridization under natural conditions. Thus, the methods of species delimitation based on the interfertility criterion are expected to give results slightly different from those based on environment-independent criteria (such as the genotypic similarity criteria). However, whatever the criterion chosen, the challenge we face when delimiting species is to summarize continuous but non-uniform variations in biological diversity. The grade of membership model that we use in this study appears as an appropriate tool.

Efficient mitigation of founder effects during the establishment of a leading-edge oak population

Numerous plant species are shifting their range polewards in response to ongoing climate change. Range shifts typically involve the repeated establishment and growth of leading-edge populations well ahead of the main species range. How these populations recover from founder events and associated diversity loss remains poorly understood. To help fill this gap, we exhaustively investigated a newly established population of holm oak (Quercus ilex) growing more than 30 km ahead of the nearest larger stands. Pedigree reconstructions showed that plants belong to two non-overlapping generations and that the whole population originates from only two founder trees. The four first-generation trees that have reached maturity showed disparate mating patterns despite being full-sibs. Long-distance pollen immigration was notable despite the strong isolation of the stand: 6 per cent gene flow events in acorns collected on the trees (n = 255), and as much as 27 per cent among their established offspring (n = 33). Our results show that isolated leading-edge populations of wind-pollinated forest trees can rapidly restore their genetic diversity through the interacting effects of efficient long-distance pollen flow and purging of inbred individuals during recruitment. They imply that range expansions of these species are primarily constrained by initial propagule arrival rather than by subsequent gene flow.

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.

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.

Fine-scale environmental control of hybridization in oaks

Natural hybridization is attracting much interest in modern speciation and conservation biology studies, but the underlying mechanisms remain poorly understood. In particular, it is unclear why environmental changes often increase hybridization rates. To study this question, we surveyed mating events in a mixed oak stand and developed a spatially explicit individual-based hybridization model. This model, where hybridization is frequency-dependent, pollen is nonlimiting and which allows immigrant pollen to compete with local pollen, takes into account species-specific pollen dispersal and sexual barriers to hybridization. The consequences of pollen limitation on hybridization were studied using another simple model. The results indicate that environmental changes could increase hybridization rates through two distinct mechanisms. First, by disrupting the spatial organization of communities, they should decrease the proportion of conspecific pollen available for mating, thus increasing hybridization rates. Second, by decreasing the density of conspecifics, they should increase pollen limitation and thus hybridization rates, as a consequence of chance pollination predominating over deterministic pollen competition. Altogether, our results point to a need for considering hybridization events at the appropriate level of organization and provide new insights into why hybridization rates generally increase in disturbed environments.

Conserving biodiversity under climate change: the rear edge matters

Modern climate change is producing poleward range shifts of numerous taxa, communities and ecosystems worldwide. The response of species to changing environments is likely to be determined largely by population responses at range margins. In contrast to the expanding edge, the low-latitude limit (rear edge) of species ranges remains understudied, and the critical importance of rear edge populations as long-term stores of species' genetic diversity and foci of speciation has been little acknowledged. We review recent findings from the fossil record, phylogeography and ecology to illustrate that rear edge populations are often disproportionately important for the survival and evolution of biota. Their ecological features, dynamics and conservation requirements differ from those of populations in other parts of the range, and some commonly recommended conservation practices might therefore be of little use or even counterproductive for rear edge populations.

Direction and extent of organelle DNA introgression between two spruce species in the Qinghai-Tibetan Plateau

A recent model has shown that, during range expansion of one species in a territory already occupied by a related species, introgression should take place preferentially from the resident species towards the invading species and genome components experiencing low rates of gene flow should introgress more readily than those experiencing high rates of gene flow. Here, we use molecular markers from two organelle genomes with contrasted rates of gene flow to test these predictions by examining genetic exchanges between two morphologically distinct spruce Picea species growing in the Qinghai-Tibetan Plateau. The haplotypes from both mitochondrial (mt) DNA and chloroplast (cp) DNA cluster into two distinct lineages that differentiate allopatric populations of the two species. By contrast, in sympatry, the species share the same haplotypes, suggesting interspecific genetic exchanges. As predicted by the neutral model, all sympatric populations of the expanding species had received their maternally inherited mtDNA from the resident species, whereas for paternally inherited cpDNA introgression is more limited and not strictly unidirectional. Our results underscore cryptic introgressions of organelle DNAs in plants and the importance of considering rates of gene flow and range shifts to predict direction and extent of interspecific genetic exchanges.

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.

Exploring species limits in two closely related Chinese oaks

BACKGROUND

The species status of two closely related Chinese oaks, Quercus liaotungensis and Q. mongolica, has been called into question. The objective of this study was to investigate the species status and to estimate the degree of introgression between the two taxa using different approaches.

METHODOLOGY/PRINCIPAL FINDINGS

Using SSR (simple sequence repeat) and AFLP (amplified fragment length polymorphism) markers, we found that interspecific genetic differentiation is significant and higher than the differentiation among populations within taxa. Bayesian clusters, principal coordinate analysis and population genetic distance trees all classified the oaks into two main groups consistent with the morphological differentiation of the two taxa rather than with geographic locations using both types of markers. Nevertheless, a few individuals in Northeast China and many individuals in North China have hybrid ancestry according to Bayesian assignment. One SSR locus and five AFLPs are significant outliers against neutral expectations in the interspecific FST simulation analysis, suggesting a role for divergent selection in differentiating species.

MAIN CONCLUSIONS/SIGNIFICANCE

All results based on SSRS and AFLPS reached the same conclusion: Q. liaotungensis and Q. mongolica maintain distinct gene pools in most areas of sympatry. They should therefore be considered as discrete taxonomic units. Yet, the degree of introgression varies between the two species in different contact zones, which might be caused by different population history or by local environmental factors.

Sex-biased dispersal promotes adaptive parental effects

Background

In heterogeneous environments, sex-biased dispersal could lead to environmental adaptive parental effects, with offspring selected to perform in the same way as the parent dispersing least, because this parent is more likely to be locally adapted. We investigate this hypothesis by simulating varying levels of sex-biased dispersal in a patchy environment. The relative advantage of a strategy involving pure maternal (or paternal) inheritance is then compared with a strategy involving classical biparental inheritance in plants and in animals.

Results

We find that the advantage of the uniparental strategy over the biparental strategy is maximal when dispersal is more strongly sex-biased and when dispersal distances of the least mobile sex are much lower than the size of the environmental patches. In plants, only maternal effects can be selected for, in contrast to animals where the evolution of either paternal or maternal effects can be favoured. Moreover, the conditions for environmental adaptive maternal effects to be selected for are more easily fulfilled in plants than in animals.

Conclusions

The study suggests that sex-biased dispersal can help predict the direction and magnitude of environmental adaptive parental effects. However, this depends on the scale of dispersal relative to that of the environment and on the existence of appropriate mechanisms of transmission of environmentally induced traits.

Relevance of genetics for conservation policies: the case of Minorcan cork oaks

BACKGROUND AND AIMS

Marginal populations of widely distributed species can be of high conservation interest when they hold a significant or unique portion of the genetic diversity of the species. However, such genetic information is frequently lacking. Here the relevance of genetic surveys to develop efficient conservation strategies for such populations is illustrated using cork oak (Quercus suber) from Minorca (Balearic Islands, Spain) as a case study. Cork oak is highly endangered on the island, where no more than 67 individuals live in small, isolated stands in siliceous sites. As a consequence, it was recently granted protected status.

METHODS

Two Bayesian clustering approaches were used to analyse the genetic structure of the Minorcan population, on the basis of nuclear microsatellite data. The different groups within the island were also compared with additional island and continental populations surrounding Minorca.

KEY RESULTS

Very high genetic diversity was found, with values comparable with those observed in continental parts of the species' range. Furthermore, the Minorcan oak stands were highly differentiated from one another and were genetically related to different continental populations of France and Spain.

CONCLUSIONS

The high levels of genetic diversity and inter-stands differentiation make Minorcan cork oak eligible for specific conservation efforts. The relationship of Minorcan stands to different continental populations in France and Spain probably reflects multiple colonization events. However, discrepancy between chloroplast DNA- and nuclear DNA-based groups does not support a simple scenario of recent introduction. Gene exchanges between neighbouring cork oak stands and with holm oak have created specific and exceptional genetic combinations. They also constitute a wide range of potential genetic resources for research on adaptation to new environmental conditions. Conservation guidelines that take into account these findings are provided.

Plant traits correlated with generation time directly affect inbreeding depression and mating system and indirectly genetic structure

Background

Understanding the mechanisms that control species genetic structure has always been a major objective in evolutionary studies. The association between genetic structure and species attributes has received special attention. As species attributes are highly taxonomically constrained, phylogenetically controlled methods are necessary to infer causal relationships. In plants, a previous study controlling for phylogenetic signal has demonstrated that Wright's F_ST, a measure of genetic differentiation among populations, is best predicted by the mating system (outcrossing, mixed-mating or selfing) and that plant traits such as perenniality and growth form have only an indirect influence on F_ST via their association with the mating system. The objective of this study is to further outline the determinants of plant genetic structure by distinguishing the effects of mating system on gene flow and on genetic drift. The association of biparental inbreeding and inbreeding depression with population genetic structure, mating system and plant traits are also investigated.

Results

Based on data from 263 plant species for which estimates of F_ST, inbreeding (F_IS) and outcrossing rate (t_m) are available, we confirm that mating system is the main influencing factor of F_ST. Moreover, using an alternative measure of F_ST unaffected by the impact of inbreeding on effective population size, we show that the influence of t_m on F_ST is due to its impact on gene flow (reduced pollen flow under selfing) and on genetic drift (higher drift under selfing due to inbreeding). Plant traits, in particular perenniality, influence F_ST mostly via their effect on the mating system but also via their association with the magnitude of selection against inbred individuals: the mean inbreeding depression increases from short-lived herbaceous to long-lived herbaceous and then to woody species. The influence of perenniality on mating system does not seem to be related to differences in stature, as proposed earlier, but rather to differences in generation time.

Conclusion

Plant traits correlated with generation time affect both inbreeding depression and mating system. These in turn modify genetic drift and gene flow and ultimately genetic structure.

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.

The hidden side of invasions: massive introgression by local genes

Despite hundreds of reports involving both plants and animals, the mechanisms underlying introgression remain obscure, even if some form of selection is frequently invoked. Introgression has repeatedly been reported in species that have recently colonized a new habitat, suggesting that demographic processes should be given more attention for understanding the mechanisms of introgression. Here we show by spatially explicit simulations that massive introgression of neutral genes takes place during the invasion of an occupied territory if interbreeding is not severely prevented between the invading and the local species. We also demonstrate that introgression occurs almost exclusively from the local to the invading species, especially for populations located far away from the source of the invasion, and this irrespective of the relative densities of the two species. This pattern is strongest at markers experiencing reduced gene flow, in keeping with the observation that organelle genes are often preferentially introgressed across species boundaries. A survey of the literature shows that a majority of published empirical studies of introgression during range expansions, in animals and in plants, follow the predictions of our model. Our results imply that speciation genes can be identified by comparing genomes of interfertile native and invading species pairs.

Ever deeper phylogeographies: trees retain the genetic imprint of Tertiary plate tectonics

Changes in species distributions after the last glacial maximum (c. 18 000 years bp) are beginning to be understood, but information diminishes quickly as one moves further back in time. In this issue of Molecular Ecology, Magri et al. (2007) present the fascinating case of a Mediterranean tree species whose populations preserve the genetic imprints of plate tectonic events that took place between 25 million years and 15 million years ago. The study provides a unique insight into the pace of evolution of trees, which, despite interspecific gene flow, can retain a cohesive species identity over timescales long enough to allow the diversification of entire plant and animal genera.

The coup de grâce for the nested clade phylogeographic analysis?

Nested clade phylogeographic analysis (NCPA) has become a popular method for reconstructing the history of populations across species ranges. Ever since its invention in 1995, criticisms have been formulated, but the method, which has been regularly updated, continues to attract investigators. Molecular Ecology has published a large fraction of the literature on the topic--both pro and con. A recent study by Panchal and Beaumont (2007) finally allows a precise evaluation of the method by developing software that automates the somewhat complicated NCPA procedure. Using simulations of random-mating populations, Panchal and Beaumont find a high frequency of false-positives with their automated NCPA procedure (over 75%). These findings, which echo and amplify earlier warnings, appear serious enough to suggest to researchers to await further evaluation of the method. Although no other all-encompassing method such as the NCAP currently exists to evaluate phylogeographic data sets, researchers have many alternative methods to test ever more refined hypotheses.

Can population genetic structure be predicted from life-history traits?

Population genetic structure is a key parameter in evolutionary biology. Earlier comparative studies have shown that genetic structure depends on species ecological attributes and life-history traits, but species phylogenetic relatedness had not been accounted for. Here we reevaluate the relationships between genetic structure and species traits in seed plants. Each species is characterized by a set of life-history and ecological features as well as by its geographic range size, its heterozygote deficit, and its genetic structure at nuclear and organelle markers to distinguish between pollen- and seed-mediated gene flow. We use both a conventional regression approach and a method that controls for phylogenetic relationships. Once phylogenetic conservatism and covariation among traits are taken into account, genetic structure is shown to be related with only a few synthetic traits, such as mating system for nuclear markers and seed dispersal mode or geographic range size for organelle markers. Along with other studies on invasiveness or rarity, our work illustrates the fact that predicting the fate of species across a broad taxonomic assemblage on the basis of simple traits is rarely possible, a testimony of the highly contingent nature of evolution.

Can Population Genetic Structure Be Predicted from Life‐History Traits?

Population genetic structure is a key parameter in evolutionary biology. Earlier comparative studies have shown that genetic structure depends on species ecological attributes and life-history traits, but species phylogenetic relatedness had not been accounted for. Here we reevaluate the relationships between genetic structure and species traits in seed plants. Each species is characterized by a set of life-history and ecological features as well as by its geographic range size, its heterozygote deficit, and its genetic structure at nuclear and organelle markers to distinguish between pollen- and seed-mediated gene flow. We use both a conventional regression approach and a method that controls for phylogenetic relationships. Once phylogenetic conservatism and covariation among traits are taken into account, genetic structure is shown to be related with only a few synthetic traits, such as mating system for nuclear markers and seed dispersal mode or geographic range size for organelle markers. Along with other studies on invasiveness or rarity, our work illustrates the fact that predicting the fate of species across a broad taxonomic assemblage on the basis of simple traits is rarely possible, a testimony of the highly contingent nature of evolution.

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'.

Blind population genetics survey of tropical rainforest trees

Rainforest tree species can be difficult to identify outside of their period of reproduction. Vascular tissues from Carapa spp. individuals were collected during a short field trip in French Guiana and analysed in the laboratory with nuclear and chloroplast markers. Using a Bayesian approach, > 90% of the samples could be assigned to one of two distinct clusters corresponding to previously described species, making it possible to estimate the genetic structure of each species and to identify cases of introgression. We argue that this blind procedure represents a first-choice rather than a fallback option whenever related taxa are investigated.

Ice-age endurance: DNA evidence of a white spruce refugium in Alaska

Paleorecords offer key information for evaluating model simulations of species migration in response to forecast climatic change. However, their utility can be greatly compromised by the existence of glacial refugia that are undetectable in fossil records (cryptic refugia). Despite several decades of investigation, it remains controversial whether Beringia, the largely unglaciated area extending from northeastern Siberia to the Yukon Territory, harbored small populations of certain boreal tree species during the last glaciation. Here, we present genetic evidence for the existence of a glacial refuge in Alaska that helps to resolve this long-standing controversy. We sequenced chloroplast DNA (cpDNA) of white spruce (Picea glauca), a dominant boreal tree species, in 24 forest stands across northwestern North America. The majority of cpDNA haplotypes are unique, and haplotype diversity is relatively high in Alaska, arguing against the possibility that this species migrated into the region from areas south of the Laurentide Ice Sheet after the end of the last glaciation. Thus, white spruce apparently survived long glacial episodes under climatic extremes in a heterogeneous landscape matrix. These results suggest that estimated rates of tree migration from fossil records may be too high and that the ability of trees to track anthropogenic warming may be more limited than previously thought.

A new scenario for the quaternary history of European beech populations: palaeobotanical evidence and genetic consequences

Here, palaeobotanical and genetic data for common beech (Fagus sylvatica) in Europe are used to evaluate the genetic consequences of long-term survival in refuge areas and postglacial spread. Four large datasets are presented, including over 400 fossil-pollen sites, 80 plant-macrofossil sites, and 450 and 600 modern beech populations for chloroplast and nuclear markers, respectively. The largely complementary palaeobotanical and genetic data indicate that: (i) beech survived the last glacial period in multiple refuge areas; (ii) the central European refugia were separated from the Mediterranean refugia; (iii) the Mediterranean refuges did not contribute to the colonization of central and northern Europe; (iv) some populations expanded considerably during the postglacial period, while others experienced only a limited expansion; (v) the mountain chains were not geographical barriers for beech but rather facilitated its diffusion; and (vi) the modern genetic diversity was shaped over multiple glacial-interglacial cycles. This scenario differs from many recent treatments of tree phylogeography in Europe that largely focus on the last ice age and the postglacial period to interpret genetic structure and argue that the southern peninsulas (Iberian, Italian and Balkan) were the main source areas for trees in central and northern Europe.

Ancient plant DNA: review and prospects

Ancient DNA has received much attention since the mid-1980s, when the first sequence of an extinct animal species was recovered from a museum specimen. Since then, the majority of ancient DNA studies have focused predominantly on animal species, while studies in plant palaeogenetics have been rather limited, with the notable exception of cultivated species found in archaeological sites. Here, we outline the recent developments in the analysis of plant ancient DNA. We emphasize the trend from species identification to population-level investigation and highlight the potential and the difficulties in this field, related to DNA preservation and to risks of contamination. Further efforts towards the analysis of ancient DNA from the abundant store of fossil plant remains should provide new research opportunities in palaeoecology and phylogeography. In particular, intraspecific variation should be considered not only in cultivated plants but also in wild taxa if palaeogenetics is to become a fully emancipated field of plant research.

Standardizing for microsatellite length in comparisons of genetic diversity

Mutation rates at microsatellites tend to increase with the number of repeats of the motif, leading to higher levels of polymorphism at long microsatellites. To standardize levels of diversity when microsatellites differ in size, we investigate the relationship between tract length and variation and provide a formula to adjust allelic richness to a fixed mean number of repeats in the specific case of chloroplast microsatellites. A comparison between 39 loci from eight species of conifers (where chloroplast DNA is paternally inherited) and 64 loci from 12 species of angiosperms (where chloroplast DNA is generally predominantly maternally inherited) indicates that the greater allelic richness found in conifers remains significant after controlling for number of repeats. The approach stresses the advantage of reporting variation in number of repeats instead of relative fragment sizes.

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.