Fine-scale spatial genetic structure in mixed oak stands with different levels of hybridization

Genomic basis of the distinct biosynthesis of β-glucogallin, a biochemical marker for hydrolyzable tannin production, in three oak species

Hydrolyzable tannins (HTs), predominant polyphenols in oaks, are widely used in grape wine aging, feed additives, and human healthcare. However, the limited availability of a high-quality reference genome of oaks greatly hampered the recognition of the mechanism of HT biosynthesis. Here, high-quality reference genomes of three Asian oak species (Quercus variabilis, Quercus aliena, and Quercus dentata) that have different HT contents were generated. Multi-omics studies were carried out to identify key genes regulating HT biosynthesis. In vitro enzyme activity assay was also conducted. Dual-luciferase and yeast one-hybrid assays were used to reveal the transcriptional regulation. Our results revealed that β-glucogallin was a biochemical marker for HT production in the cupules of the three Asian oaks. UGT84A13 was confirmed as the key enzyme for β-glucogallin biosynthesis. The differential expression of UGT84A13, rather than enzyme activity, was the main reason for different β-glucogallin and HT accumulation. Notably, sequence variations in UGT84A13 promoters led to different trans-activating activities of WRKY32/59, explaining the different expression patterns of UGT84A13 among the three species. Our findings provide three high-quality new reference genomes for oak trees and give new insights into different transcriptional regulation for understanding β-glucogallin and HT biosynthesis in closely related oak species.

Genetic Diversity and Spatial Genetic Structure in Isolated Scots Pine (Pinus sylvestris L.) Populations Native to Eastern and Southern Carpathians

Small, isolated populations are more vulnerable to natural disturbances and loss of genetic diversity. Scots pine, an abundant tree species in the boreal forest of Eurasia, has a scattered natural distribution across Eastern and Southern Carpathian Mountains, where only a few relict populations still exist. We estimated genetic diversity and spatial genetic structure in Scots pine on the basis of microsatellite nuclear markers (nSSR) data. We found a relatively high level of genetic diversity (He = 0.697) within populations and no evidence of recent bottlenecks. Genetic diversity was lower in peat bog populations, as compared to populations that grow on rocky slopes or acidic soils and nutrient-poor sites. Population genetic structure was weak, and genetic discontinuities among populations were detected. Spatial genetic structure (SGS) was observed in nearly all Scots pine populations. The strength of SGS, quantified by Sp statistics, varied greatly among populations, ranging from 0.0011 to 0.0207, with an average of 0.01. Our study highlights that Eastern and Southern Carpathian populations still possess high within-population diversity in spite of the recent fragmentation and reduction of the Scots pine natural distribution range. We discuss the importance of spatial patterns of genetic diversity for developing strategies of conservation and sustainable use of Scots pine genetic resources in the Carpathian region.

Climate Change Synchronizes Growth and iWUE Across Species in a Temperate-Submediterranean Mixed Oak Forest

Tree species have good tolerance to a range of environmental conditions, though their ability to respond and persist to environmental changes is dramatically reduced at the rear-edge distribution limits. At those edges, gene flow conferring adaptation is impaired due to lack of populations at lower latitudes. Thus, trees mainly rely on phenotypic changes to buffer against long-term environmental changes. Interspecific hybridization may offer an alternative mechanism in the generation of novel genetic recombinants that could be particularly valuable to ensure persistence in geographically isolated forests. In this paper, we take advantage of the longevity of a temperate-submediterranean mixed-oak forest to explore the long-term impact of environmental changes on two different oak species and their hybrid. Individual trees were genetically characterized and classified into three groups: pure Quercus petraea (Matt.), Liebl, pure Q. pyrenaica Willd, and hybrids. We calculated basal area increment and intrinsic water-use efficiency (iWUE) from tree-ring width and δ13C per genetic group, respectively. Tree-growth drivers were assessed using correlation analyses and generalized linear mixed models for two contrasting climatic periods: (1880-1915, colder with [CO2] < 303 ppm; and 1980-2015, warmer with [CO2] > 338 ppm). The three genetic groups have increased radial growth and iWUE during the last decades, being the least drought-tolerant QuPe the most sensitive species to water stress. However, no significant differences were found among genetic groups neither in mean growth rate nor in mean iWUE. Furthermore, little differences were found in the response to climate among groups. Genetic groups only differed in the relationship between δ13C and temperature and precipitation during the earlier period, but such a difference disappeared during the recent decades. Climate change may have promoted species-level convergence as a response to environment-induced growth limitations, which translated in synchronized growth and response to climate as well as a tighter stomatal control and increased iWUE across coexisting oak species.

Hybridization and geographic distribution shapes the spatial genetic structure of two co-occurring orchid species

Multiple ecological and life-history traits shape the fine-scale spatial genetic structure (FSGS) of a given population. The occurrence in core versus peripheral populations, levels of outcrossing, pollen and seed dispersal, and hybridization are important biological properties that influence the kinship of individuals within populations. We examined spatial genetic structure within 15 populations of Epidendrum fulgens and E. puniceoluteum distributed along a linear gradient of Brazilian coastal vegetation, including both allopatric and sympatric populations where the two orchid species hybridize. We analyzed 581 mapped specimens using nine simple sequence repeat loci, aiming to investigate how geographic distribution and hybridization shape within-population FSGS. A significant increase in FSGS was found towards peripheral populations, compared to core populations. Analysis of short-distance and long-distance components of FSGS identified biparental inbreeding and higher levels of FSGS at peripheral populations, when compared to core populations. In contrast, the relatively high density of reproductive adults in core populations potentially leads to highly overlapping seed and pollen movement, decreasing FSGS. Hybridization was an important factor shaping within-population spatial genetic structure at sympatric sites, decreasing the FSGS observed in parental species. Our results indicate that different ecological forces act in concert to create a gradient of FSGS along species distribution ranges, shaped by extensive levels of intraspecific and interspecific gene exchange.

Quantitative evaluation of hybridization and the impact on biodiversity conservation

Anthropogenic hybridization is an increasing conservation threat worldwide. In South Africa, recent hybridization is threatening numerous ungulate taxa. For example, the genetic integrity of the near‐threatened bontebok (Damaliscus pygargus pygargus) is threatened by hybridization with the more common blesbok (D. p. phillipsi). Identifying nonadmixed parental and admixed individuals is challenging based on the morphological traits alone; however, molecular analyses may allow for accurate detection. Once hybrids are identified, population simulation software may assist in determining the optimal conservation management strategy, although quantitative evaluation of hybrid management is rarely performed. In this study, our objectives were to describe species‐wide and localized rates of hybridization in nearly 3,000 individuals based on 12 microsatellite loci, quantify the accuracy of hybrid assignment software (STRUCTURE and NEWHYBRIDS), and determine an optimal threshold of bontebok ancestry for management purposes. According to multiple methods, we identified 2,051 bontebok, 657 hybrids, and 29 blesbok. More than two‐thirds of locations contained at least some hybrid individuals, with populations varying in the degree of introgression. HYBRIDLAB was used to simulate four generations of coexistence between bontebok and blesbok, and to optimize a threshold of ancestry, where most hybrids will be detected and removed, and the fewest nonadmixed bontebok individuals misclassified as hybrids. Overall, a threshold Q‐value (admixture coefficient) of 0.90 would remove 94% of hybrid animals, while a threshold of 0.95 would remove 98% of hybrid animals but also 8% of nonadmixed bontebok. To this end, a threshold of 0.90 was identified as optimal and has since been implemented in formal policy by a provincial nature conservation agency. Due to widespread hybridization, effective conservation plans should be established and enforced to conserve native populations that are genetically unique.

Plastome data reveal multiple geographic origins of Quercus Group Ilex

Nucleotide sequences from the plastome are currently the main source for assessing taxonomic and phylogenetic relationships in flowering plants and their historical biogeography at all hierarchical levels. One major exception is the large and economically important genus Quercus (oaks). Whereas differentiation patterns of the nuclear genome are in agreement with morphology and the fossil record, diversity patterns in the plastome are at odds with established taxonomic and phylogenetic relationships. However, the extent and evolutionary implications of this incongruence has yet to be fully uncovered. The DNA sequence divergence of four Euro-Mediterranean Group Ilex oak species (Quercus ilex L., Q. coccifera L., Q. aucheri Jaub. & Spach., Q. alnifolia Poech.) was explored at three chloroplast markers (rbcL, trnK/matK, trnH-psbA). Phylogenetic relationships were reconstructed including worldwide members of additional 55 species representing all Quercus subgeneric groups. Family and order sequence data were harvested from gene banks to better frame the observed divergence in larger taxonomic contexts. We found a strong geographic sorting in the focal group and the genus in general that is entirely decoupled from species boundaries. High plastid divergence in members of Quercus Group Ilex, including haplotypes shared with related, but long isolated oak lineages, point towards multiple geographic origins of this group of oaks. The results suggest that incomplete lineage sorting and repeated phases of asymmetrical introgression among ancestral lineages of Group Ilex and two other main Groups of Eurasian oaks (Cyclobalanopsis and Cerris) caused this complex pattern. Comparison with the current phylogenetic synthesis also suggests an initial high- versus mid-latitude biogeographic split within Quercus. High plastome plasticity of Group Ilex reflects geographic area disruptions, possibly linked with high tectonic activity of past and modern distribution ranges, that did not leave imprints in the nuclear genome of modern species and infrageneric lineages.

Plastome data reveal multiple geographic origins of Quercus Group Ilex

Nucleotide sequences from the plastome are currently the main source for assessing taxonomic and phylogenetic relationships in flowering plants and their historical biogeography at all hierarchical levels. One major exception is the large and economically important genus Quercus (oaks). Whereas differentiation patterns of the nuclear genome are in agreement with morphology and the fossil record, diversity patterns in the plastome are at odds with established taxonomic and phylogenetic relationships. However, the extent and evolutionary implications of this incongruence has yet to be fully uncovered. The DNA sequence divergence of four Euro-Mediterranean Group Ilex oak species (Quercus ilex L., Q. coccifera L., Q. aucheri Jaub. & Spach., Q. alnifolia Poech.) was explored at three chloroplast markers (rbcL, trnK/matK, trnH-psbA). Phylogenetic relationships were reconstructed including worldwide members of additional 55 species representing all Quercus subgeneric groups. Family and order sequence data were harvested from gene banks to better frame the observed divergence in larger taxonomic contexts. We found a strong geographic sorting in the focal group and the genus in general that is entirely decoupled from species boundaries. High plastid divergence in members of Quercus Group Ilex, including haplotypes shared with related, but long isolated oak lineages, point towards multiple geographic origins of this group of oaks. The results suggest that incomplete lineage sorting and repeated phases of asymmetrical introgression among ancestral lineages of Group Ilex and two other main Groups of Eurasian oaks (Cyclobalanopsis and Cerris) caused this complex pattern. Comparison with the current phylogenetic synthesis also suggests an initial high- versus mid-latitude biogeographic split within Quercus. High plastome plasticity of Group Ilex reflects geographic area disruptions, possibly linked with high tectonic activity of past and modern distribution ranges, that did not leave imprints in the nuclear genome of modern species and infrageneric lineages.

Individual spatial aggregation correlates with between-population variation in fine-scale genetic structure of Silene ciliata (Caryophyllaceae)

Fine-scale genetic structure (FSGS) can vary among populations within species depending on multiple demographic and environmental factors. Theoretical models predict that FSGS should decrease in high-density populations and increase in populations where individuals are spatially aggregated. However, few empirical studies have compared FSGS between populations with different degrees of individual spatial aggregation and microhabitat heterogeneity. In this work, we studied the relationship between spatial and genetic structure in five populations of alpine specialist Silene ciliata Poiret (Caryophyllaceae). We mapped all individuals in each population and genotyped 96 of them using 10 microsatellite markers. We found significant FSGS consistent with an isolation-by-distance process in three of the five populations. The intensity of FSGS was positively associated with individual spatial aggregation. However, no association was found between FSGS and global population density or microhabitat heterogeneity. Overall, our results support theoretical studies indicating that stronger spatial aggregation tends to increase the magnitude of FSGS. They also highlight the relevance of characterizing local plant distribution and microhabitat to better understand the mechanisms that generate intraspecific variation in FSGS across landscapes.

High rates of gene flow by pollen and seed in oak populations across Europe

Gene flow is a key factor in the evolution of species, influencing effective population size, hybridisation and local adaptation. We analysed local gene flow in eight stands of white oak (mostly Quercus petraea and Q. robur, but also Q. pubescens and Q. faginea) distributed across Europe.

Adult trees within a given area in each stand were exhaustively sampled (range [239, 754], mean 423), mapped, and acorns were collected ([17,147], 51) from several mother trees ([3], [47], 23). Seedlings ([65,387], 178) were harvested and geo-referenced in six of the eight stands. Genetic information was obtained from screening distinct molecular markers spread across the genome, genotyping each tree, acorn or seedling. All samples were thus genotyped at 5–8 nuclear microsatellite loci. Fathers/parents were assigned to acorns and seedlings using likelihood methods. Mating success of male and female parents, pollen and seed dispersal curves, and also hybridisation rates were estimated in each stand and compared on a continental scale.

On average, the percentage of the wind-borne pollen from outside the stand was 60%, with large variation among stands (21–88%). Mean seed immigration into the stand was 40%, a high value for oaks that are generally considered to have limited seed dispersal. However, this estimate varied greatly among stands (20–66%). Gene flow was mostly intraspecific, with large variation, as some trees and stands showed particularly high rates of hybridisation.

Our results show that mating success was unevenly distributed among trees. The high levels of gene flow suggest that geographically remote oak stands are unlikely to be genetically isolated, questioning the static definition of gene reserves and seed stands.

Genetic evidence for hybridization in red oaks (Quercus sect. Lobatae, Fagaceae)

PREMISE OF THE STUDY

Hybridization is pervasive in many plant taxa, with consequences for species taxonomy, local adaptation, and management. Oaks (Quercus spp.) are thought to hybridize readily yet retain distinct traits, drawing into question the biological species concept for such taxa, but the true extent of gene flow is controversial. Genetic data are beginning to shed new light on this issue, but red oaks (section Lobatae), an important component of North American forests, have largely been neglected. Moreover, gene flow estimates may be sensitive to the choice of life stage, marker type, or genetic structure statistic.

METHODS

We coupled genetic structure data with parentage analyses for two mixed-species stands in North Carolina. Genetic structure analyses of adults (including F(ST), R(ST), G'(ST), and structure) reflect long-term patterns of gene flow, while the percentage of seedlings with parents of two different species reflect current levels of gene flow.

KEY RESULTS

Genetic structure analyses revealed low differentiation in microsatellite allele frequencies between co-occurring species, suggesting past gene flow. However, methods differed in their sensitivity to differentiation, indicating a need for caution when drawing conclusions from a single method. Parentage analyses identified >20% of seedlings as potential hybrids. The species examined exhibit distinct morphologies, suggesting selection against intermediate phenotypes.

CONCLUSIONS

Our results suggest that hybridization between co-occurring red oaks occurs, but that selection may limit introgression, especially at functional loci. However, by providing a source of genetic variation, hybridization could influence the response of oaks and other hybridizing taxa to environmental change.

Patterns of hybridization and asymmetrical gene flow in hybrid zones of the rare Eucalyptus aggregata and common E. rubida

The patterns of hybridization and asymmetrical gene flow among species are important for understanding the processes that maintain distinct species. We examined the potential for asymmetrical gene flow in sympatric populations of Eucalyptus aggregata and Eucalyptus rubida, both long-lived trees of southern Australia. A total of 421 adults from three hybrid zones were genotyped with six microsatellite markers. We used genealogical assignments, admixture analysis and analyses of spatial genetic structure and spatial distribution of individuals, to assess patterns of interspecific gene flow within populations. A high number of admixed individuals were detected (13.9-40% of individuals), with hybrid populations consisting of F(1) and F(2) hybrids and backcrosses in both parental directions. Across the three sites, admixture proportions were skewed towards the E. aggregata genetic cluster (x=0.56-0.65), indicating that backcrossing towards E. aggregata is more frequent. Estimates of long-term migration rates also indicate asymmetric gene flow, with higher migration rates from E. aggregata to hybrids compared with E. rubida. Taken together, these results indicate a greater genetic input from E. aggregata into the hybrid populations. This asymmetry probably reflects differences in style lengths (E. rubida: ~7 mm, E. aggregata: ~4 mm), which can prevent pollen tubes of smaller-flowered species from fertilizing larger-flowered species. However, analyses of fine-scale genetic structure suggest that localized seed dispersal (<40 m) and greater clustering between hybrid and E. aggregata individuals may also contribute to directional gene flow. Our study highlights that floral traits and the spatial distributions of individuals can be useful predictors of the directionality of interspecific gene flow in plant populations.

Fine-scale spatial genetic structure and gene flow in a small, fragmented population of Sinojackia rehderiana (Styracaceae), an endangered tree species endemic to China

Populations of Sinojackia rehderiana are highly threatened and have small and scattered distribution due to habitat fragmentation and human activities. Understanding changes in genetic diversity, the fine-scale spatial genetic structure (SGS) at different life stages and gene flow of S. rehderiana is critical for developing successful conservation strategies for fragmented populations of this endangered species. In this study, 208 adults, 114 juveniles and 136 seedlings in a 50 × 100-m transect within an old-growth forest were mapped and genotyped using eight microsatellite makers to investigate the genetic diversity and SGS of this species. No significant differences in genetic diversity among different life-history stages were found. However, a significant heterozygote deficiency in adults and seedlings may result from substantial biparental inbreeding. Significant fine-scale spatial structure was found in different life-history stages within 19 m, suggesting that seed dispersal mainly occurred near a mother tree. Both historical and contemporary estimates of gene flow (13.06 and 16.77 m) indicated short-distance gene dispersal in isolated populations of S. rehderiana. The consistent spatial structure revealed in different life stages is most likely the result of limited gene flow. Our results have important implications for conservation of extant populations of S. rehderiana. Measures for promoting pollen flow should be taken for in situ conservation. The presence of a SGS in fragmented populations implies that seeds for ex situ conservation should be collected from trees at least 19-m apart to reduce genetic similarity between neighbouring individuals.

Reproductive patterns shape introgression dynamics and species succession within the European white oak species complex

The reproductive system of hybrids is an important factor shaping introgression dynamics within species complexes. We combined paternity and parentage analyses with previous species characterization by genetic assignment, to directly identify reproductive events that occurred within a stand comprising four European white oak species. Comparing species status of parent pairs provided a precise quantification of hybridization rate, backcrosses, and intraspecific matings in two life stages. The detailed mating system analysis revealed new findings on the dynamics of interspecific gene flow. First, hybrids acted successfully as both male and female during reproduction. They produced acorns and seedlings that were as viable as those sired by purebreds. Second, species maintenance could be due to a relatively low level of interspecific mating contrasting with a large proportion of intraspecific crosses and backcrosses. Despite a high proportion of hybrids and extensive interspecific gene flow, partial species integrity is maintained by genetically controlled pollen discrimination, ensuring preferential matings within purebreds and high parental species fidelity in hybrid reproduction, which impedes complete collapse into a continuous hybrid swarm. Finally, we showed that pollen from the different species had unequal contributions to reproduction suggesting that introgression processes could ultimately lead to extirpation or expansion of some species.

Spatial and temporal genetic structure in a hybrid cordgrass invasion

Invasive hybrids and their spread dynamics pose unique opportunities to study evolutionary processes. Invasive hybrids of native Spartina foliosa and introduced S. alterniflora have expanded throughout San Francisco Bay intertidal habitats within the past 35 years by deliberate plantation and seeds floating on the tide. Our goals were to assess spatial and temporal scales of genetic structure in Spartina hybrid populations within the context of colonization history. We genotyped adult and seedling Spartina using 17 microsatellite loci and mapped their locations in three populations. All sampled seedlings were hybrids. Bayesian ordination analysis distinguished hybrid populations from parent species, clearly separated the population that originated by plantation from populations that originated naturally by seed and aligned most seedlings within each population. Population genetic structure estimated by analysis of molecular variance was substantial (F(ST)=0.21). Temporal genetic structure among age classes varied highly between populations. At one population, the divergence between adults and 2004 seedlings was low (F(ST)=0.02) whereas at another population this divergence was high (F(ST)=0.26). This latter result was consistent with local recruitment of self-fertilized seed produced by only a few parental plants. We found fine-scale spatial genetic structure at distances less than ∼200 m, further supporting local seed and/or pollen dispersal. We posit a few self-fertile plants dominating local recruitment created substantial spatial genetic structure despite initial long-distance, human dispersal of hybrid Spartina through San Francisco Bay. Fine-scale genetic structure may more strongly develop when local recruits are dominated by the offspring of a few self-fertile plants.

Fine-scale spatial genetic structure and gene dispersal in Silene latifolia

Plants are sessile organisms, often characterized by limited dispersal. Seeds and pollen are the critical stages for gene flow. Here we investigate spatial genetic structure, gene dispersal and the relative contribution of pollen vs seed in the movement of genes in a stable metapopulation of the white campion Silene latifolia within its native range. This short-lived perennial plant is dioecious, has gravity-dispersed seeds and moth-mediated pollination. Direct measures of pollen dispersal suggested that large populations receive more pollen than small isolated populations and that most gene flow occurs within tens of meters. However, these studies were performed in the newly colonized range (North America) where the specialist pollinator is absent. In the native range (Europe), gene dispersal could fall on a different spatial scale. We genotyped 258 individuals from large and small (15) subpopulations along a 60 km, elongated metapopulation in Europe using six highly variable microsatellite markers, two X-linked and four autosomal. We found substantial genetic differentiation among subpopulations (global F(ST)=0.11) and a general pattern of isolation by distance over the whole sampled area. Spatial autocorrelation revealed high relatedness among neighboring individuals over hundreds of meters. Estimates of gene dispersal revealed gene flow at the scale of tens of meters (5-30 m), similar to the newly colonized range. Contrary to expectations, estimates of dispersal based on X and autosomal markers showed very similar ranges, suggesting similar levels of pollen and seed dispersal. This may be explained by stochastic events of extensive seed dispersal in this area and limited pollen dispersal.

Natural hybridisation between kermes (Quercus coccifera L.) and holm oaks (Q. ilex L.) revealed by microsatellite markers

Hybridisation between species of the genus Quercus is a common phenomenon as a result of weak reproductive isolation mechanisms between phylogenetically close species that frequently co-occur in mixed stands. In this study, we use microsatellite markers to analyse introgression between kermes (Quercus coccifera L.) and holm (Q. ilex L.) oak, two closely related taxa that frequently dominate the landscape in extensive areas in the Mediterranean region. All tested microsatellites amplified and were polymorphic in both kermes and holm oaks. Bayesian admixture analyses showed a good correspondence between each species and one of the two inferred genetic clusters. Five sampled individuals were a priori tentatively identified as hybrids on the basis of intermediate morphological characteristics, and it was confirmed that they also presented mixed genotypes. However, we also detected different levels of genetic introgression among morphologically pure individuals, suggesting that successful backcrossing and/or reduced phenotypic expression of genetic variance in certain individuals may have resulted in strong convergence towards a single species phenotype.

Patterns of contemporary hybridization inferred from paternity analysis in a four-oak-species forest

Background

Few studies address the issue of hybridization in a more than two-species context. The species-rich Quercus complex is one of the systems which can offer such an opportunity. To investigate the contemporary pattern of hybridization we sampled and genotyped 320 offspring from a natural mixed forest comprising four species of the European white oak complex: Quercus robur, Q. petraea, Q. pubescens, and Q. frainetto.

Results

A total of 165 offspring were assigned unambiguously to one of the pollen donors within the study plot. The minimum amount of effective pollen originating from outside the plot varied markedly among the seed parents, ranging from 0.18 to 0.87. The majority of the successful matings (64.1%) occurred between conspecific individuals indicating the existence of reproductive barriers between oak species. However, the isolation was not complete since we found strong evidence for both first-generation (8.4%) and later-generation hybrids (27.5%). Only two out of eight seed parents, belonging to Q. petraea and Q. robur, showed a high propensity to hybridize with Q. pubescens and Q. petraea, respectively. Significant structure of the effective pollen pools (Φ_pt = 0.069, P = 0.01) was detected in our sample. However, no support was found for the isolation by distance hypothesis. The proportion of hybrids was much higher (79%) in the seed generation when compared to the adult tree generation.

Conclusion

First-generation hybrids were observed only between three out of six possible species combinations. Hybrids between one pair of species preferred to mate with one of their parental species. The observation of first and later-generation hybrids in higher frequency in acorns than in adults might be explained by selection against hybrid genotypes, the history of this uneven-aged forest or past introgression between species.

Ancient and current gene flow between two distantly related Mediterranean oak species, Quercus suber and Q. ilex

BACKGROUND AND AIMS

Quercus suber and Q. ilex are distantly related and their distributions partially overlap. They hybridize occasionally, but the complete replacement of Q. suber chloroplast DNA (cpDNA) by that of Q. ilex was identified in two specific geographical areas. The objective of this study was to determine whether the contrasting situation reflected current or recent geographical interspecies gene flow variation or was the result of ancient introgression.

METHODS

cpDNA PCR-RFLPs (restriction fragment length polymorphisms) and variation at ten nuclear microsatellite loci were analysed in populations of each species, in 16 morphologically intermediate individuals and the progeny of several of them. Interspecies nuclear introgression was based on individual admixture rates using a Bayesian approach with no a priori species assignment, and on a maximum-likelihood (ML) method, using allele frequencies in the allopatric populations of each species as controls. Gene flow was compared specifically between populations located within and outside the specific areas.

KEY RESULTS

High interspecies nuclear genetic differentiation was observed, with twice the number of alleles in Q. ilex than in Q. suber. According to Bayesian assignment, approx. 1 % of individuals had a high probability of being F(1) hybrids, and bidirectional nuclear introgression affected approx. 4 % of individuals in each species. Hybrid and introgressed individuals were identified predominantly in mixed stands and may have a recent origin. Higher proportions including allospecific genes recovered from past hybridization were obtained using the ML method. Similar rates of hybridization and of nuclear introgression, partially independent of cpDNA interspecies transfer suggestive of gene filtering, were obtained in the populations located within and outside the areas of complete cpDNA replacement.

CONCLUSIONS

The results did not provide evidence for geographical variation in interspecies gene flow. In contrast, historical introgression is supported by palynological records and constitutes the more reliable origin of cpDNA replacement in specific regions.

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.

Contrasting nuclear and cytoplasmic exchanges between phylogenetically distant oak species (Quercus suber L. and Q. ilex L.) in Southern France: inferring crosses and dynamics

Gene flow is particularly frequent in the genus Quercus (oaks), especially between closely related species. We focus here on Quercus ilex and the cork-producing Quercus suber, which occasionally hybridize although they are phylogenetically markedly separated. Morphological observations were combined with both allozymic and chloroplastic diagnostic markers to characterize hybridization and introgression and to infer their dynamics in two French regions (French Catalonia and Provence), which are separated by several hundred kilometres. Some hybrids were found in both regions, indicating recent hybridization events. As expected from previous studies, most hybrids resulted from female symbol Q. ilex x male symbol Q. suber crosses, but our data showed that the reciprocal cross is also possible. Partial independence between nuclear and chloroplastic introgression was observed in the two species. Nuclear introgression was limited in both species and both regions, with no preferred direction. In Provence, chloroplastic introgression was very rare in both species. Conversely, all Q. suber individuals from French Catalonia were introgressed by Q. ilex chlorotypes. This might be explained by introgression in the Iberian Peninsula antedating the first occurrence of the two species in French Catalonia. We also observed a new chlorotype that was created locally, and was exchanged between the two species. However, the two species still remain genetically differentiated. The dynamics and complexity of exchanges and the factors determining them (including human management of Q. suber) are discussed.

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.

Clonality and spatial genetic structure in Populus x canescens and its sympatric backcross parent P. alba in a Central European hybrid zone

Spatial genetic structure (SGS) holds the key to understanding the role of clonality in hybrid persistence, but multilocus SGS in hybrid zones has rarely been quantified. Here, the aim was to fill this gap for natural hybrids between two diploid, ecologically divergent European tree species with mixed sexual/asexual reproduction, Populus alba and P. tremula. Nuclear microsatellites were used to quantify clonality, SGS, and historical gene dispersal distances in up to 407 trees from an extensive Central European hybrid zone including three subpopulation replicates. The focus was on P. x canescens and its backcross parent P. alba, as these two genotypic classes co-occur and interact directly. Sexual recombination in both taxa was more prominent than previously thought, but P. x canescens hybrids tended to build larger clones extending over larger areas than P. alba. The 3.4 times stronger SGS in the P. x canescens genet population was best explained by a combination of interspecific gene flow, assortative mating, and increased clonality in hybrids. Clonality potentially contributes to the maintenance of hybrid zones of P. alba and P. tremula in time and space. Both clonality and SGS need to be taken into account explicitly when designing population genomics studies of locus-specific effects in hybrid zones.