Distinct male reproductive strategies in two closely related oak species

Reproductive strategies of closely related species distributed along successional gradients should differ as a consequence of the trade‐off between competition and colonization abilities. We compared male reproductive strategies of Quercus robur and Q. petraea, two partly interfertile European oak species with different successional status. In the studied even‐aged stand, trees of the late‐successional species (Q. petraea) grew faster and suffered less from intertree competition than trees of the early‐successional species (Q. robur). A large‐scale paternity study and a spatially explicit individual‐based mating model were used to estimate parameters of pollen production and dispersal as well as sexual barriers between species. Male fecundity was found to be dependent both on a tree's circumference and on its environment, particularly so for Q. petraea. Pollen dispersal was greater and more isotropic in Q. robur than in Q. petraea. Premating barriers to hybridization were strong in both species, but more so in Q. petraea than in Q. robur. Hence, predictions based on the competition–colonization trade‐off are well supported, whereas the sexual barriers themselves seem to be shaped by colonization dynamics.     

Electrophoretic Analysis of Genetic Variation within and between Populations of Quercus cerris,Q. pubescens,Q. petraea and Q. robur (Fagaceae) from Eastern Austria

Allozyme variation at 14 loci was studied electrophoretically in 19 population samples of Quercus cerris, Q. pubescens, Q. petraea, and Q. robur. Genetic variation, measured as mean heterozygosity, was found to be high within populations but low among populations. A significant deficit of heterozygotes was observed in most of the populations of Q. pubescens, Q. petraea, and Q. robur. At the species level, Q. cerris (sect. Cerris) appears to be clearly separated from the rest (sect. Quercus). Q. pubescens, Q. petraea, and Q. robur form a cluster of closely-related species. The degree of genetic differentiation among these species is low, i.e. in the range normally found among con-specific populations. Nevertheless, allozyme data allow the elaboration of a dendrogram which is in reasonable accordance with the taxonomic classification, but also supports the importance of hybridization and introgression. This is also underlined by the fact that the allelic differentiation at the Got-2 locus corresponds with the phenotypic classification (petraea-like, hybrids, and robur-like) in samples from mixed populations.     

Hybridization and divergence in multi‐species oak (Quercus) communities

Oaks (Quercus: Fagaceae) commonly interbreed yet retain their morphological, genetic and ecological distinctiveness. Post‐zygotic isolation mechanisms, such as ecologically dependent selection on adaptive loci, may therefore limit introgression. To test this hypothesis, we quantified hybridization and genetic divergence across the contact zone of four red oaks (Quercus section Lobatae) in the Great Lakes region of North America using a suite of 259 amplified fragment length polymorphisms and 27 genic and genomic microsatellite markers. First, we identified hybrids using genetic structure analysis and confirmed the reliability of our assignments via simulations. Then, we identified candidate loci for species maintenance with three complementary tests for selection and obtained partial gene sequences linked to an outlier locus and three other loci. We detected evidence of recent hybridization among all species and considerable gene flow between Q. ellipsoidalis and Q. velutina. Overall, c. 20% of Q. velutina had recent ancestry from Q. ellipsoidalis, whereas nearly 30% of Q. ellipsoidalis had a Q. velutina ancestor. Most loci were negligibly to weakly differentiated among species, but two gene‐linked microsatellites deviated significantly from neutral expectations in multiple, complementary outlier tests. Both outlier loci were located in the same 15‐cM bin on an existing Q. robur linkage map, a region under divergent selection in other oak species. Adaptive loci in this highly differentiated genomic region may contribute to ecological divergence among species and limit introgression.     

Flooding effects on starch partitioning during early growth of two oak species

Pedunculate (Quercus robur L.) and sessile (Q. petraea [Matt.] Liebl.) oaks are the most common oak species in Western Europe. They are known to display different ecological requirements, particularly relative to root hypoxia induced by flooding: In a glasshouse study of seedlings, we quantified the effects of flooding on starch mobilization from cotyledons and starch partitioning. Growth and distribution of lateral roots were also measured. The above-ground growth of Q. robur was less affected by flooding than that of Q. petraea which failed to develop a second flush. Root growth was also severely inhibited, particularly in Q. petraea. In Q. robur, lateral root initiation as well as elongation was restricted to the soil surface layer. Flooding markedly reduced total growth and concentrations of in all components except stems. Starch mobilization from cotyledons was delayed by flooding, especially in Q. robur seedlings. Under flooding, the decrease of cotyledons dry mass and starch content in Q. robur was lower than in Q. petraea, whereas Q. robur displayed larger growth than Q. petraea. The features of carbohydrate management may be crucial in the observed differences in flooding tolerance of these species.     

Development of ten microsatellite markers for <Emphasis Type="Italic">Quercus mongolica</Emphasis> var. <Emphasis Type="Italic">crispula</Emphasis> by database mining

Simple sequence repeats (SSRs) in the NCBI dbEST database were surveyed to identify potential SSR markers for Quercus mongolica. In total, 2,691 gene sequences, mainly from expressed sequence tags (ESTs) for Q. robur and Q. petraea had been registered. Twenty-two PCR primers were designed for SSRs in these sequences and screened for polymorphisms in 16 Q. mongolica trees. Ten loci were easily genotyped and showed polymorphism, with numbers of alleles and expected heterozygosity ranging from 3 to 15 and 0.28 to 0.94, respectively. These EST-SSR markers should be useful for studying the genetic diversity of Quercus species.     

Insights into drought adaptation of two European oak species revealed by nucleotide diversity of candidate genes

Nucleotide diversity, correlations to environmental conditions, and deviations from standard neutral models were examined in eight candidate genes for drought response of European oaks Quercus petraea and Quercus robur. Baseline genetic parameters were calculated and deviation from neutrality was tested using Tajima’s D, and Fu and Li’s D and F statistics. Population structure was investigated using differentiation indices and Bayesian clustering. Similar values of total nucleotide diversity were found in both species with significant higher diversity at non-synonymous sites in Q. robur while the number of haplotypes was significantly higher in Q. petraea. Significant population differentiation was found for three genes in Q. petraea and for one gene in Q. robur. Within four genes, strong correlations were found between the local temperature–precipitation regime and the allele frequencies of six alleles, of which three were private to Q. petraea. Using various population genetic and Bayesian tests for neutrality, four outlier single nucleotide polymorphisms (SNPs) under putative selection were detected in two of the analyzed genes. Significant differentiation and strong allelic correlation to environmental conditions support preceding gene expression profiling experiments, where functional impact of candidate genes in drought response has been revealed. Q. petraea populations were found to be more differentiated as compared to Q. robur—this could be linked to the higher adaptive potential of this species under arid conditions.     

Genetic evidence of reproductive isolation in a remote enclave of Quercus pubescens in the presence of cross-fertile species

Peripheral populations may be crucial for understanding processes underlying adaptive genetic variation. Their evolution and ecology are driven by various genetic and demographic processes, such as selection, gene flow and bottleneck. Peripheral populations often experience a reduction in density resulting in the Allee effect. The presence of interfertile species increases the opportunity for hybridisation, which allows for a rescue from the Allee effect, but at the risk of genetic extinction through introgression. In this article we investigated a peripheral population of Quercus pubescens, a European tree species. The study population is located in Poland, several hundred kilometres northwards from the main species range. Due to geographic separation, the study population exists under strong pressure of introgression from potentially inter-fertile Q. petraea and Q. robur, which are the only common oaks in Poland. The intermediate morphology between typical Q. pubescens and a common oak species found in the study population supports the introgression hypothesis, which could be in line with the earlier studies of this species complex conducted in the main range of Q. pubescens. Alternatively, the intermediate morphology could reflect the founder effect or selection at an ecological extreme. We attempted to verify these hypotheses using microsatellites and a reference of common oak species. The results showed that the study population is genetically distinct from both Q. petraea and Q. robur. Additionally, the population is characterised by a low effective population size and limited gene dispersal. This suggests that the study population reveals strong reproductive isolation from common species, implying alternative sources of atypical morphology.     

Differential introgression from a sister species explains high FST outlier loci within a mussel species

Scanning genomes for loci with high levels of population differentiation has become a standard of population genetics. F_ST outlier loci are most often interpreted as signatures of local selection, but outliers might arise for many other reasons too often left unexplored. Here, we tried to identify further the history and genetic basis underlying strong differentiation at F_ST outlier loci in a marine mussel. A genome scan of genetic differentiation has been conducted between Atlantic and Mediterranean populations of Mytilus galloprovincialis. The differentiation was low overall (F_ST = 0.03), but seven loci (2%) were strong F_ST outliers. We then analysed DNA sequence polymorphism at two outlier loci. The genetic structure proved to be the consequence of differential introgression of alleles from the sister‐hybridizing species Mytilus edulis. Surprisingly, the Mediterranean population was the most introgressed at these two loci, although the contact zone between the two species is nowadays localized along the Atlantic coasts of France and the British Isles. A historical contact between M. edulis and Mediterranean M. galloprovincialis should have happened during glacial periods. It proved difficult to disentangle two hypotheses: (i) introgression was adaptive, implying edulis alleles have been favoured in Mediterranean populations, or (ii) the genetic architecture of the barrier to edulis gene flow is different between the two M. galloprovincialis backgrounds. Five of the seven outliers between M. galloprovincialis populations were also outliers between M. edulis and Atlantic M. galloprovincialis, which would support the latter hypothesis. Differential introgression across semi‐permeable barriers to gene flow is a neglected scenario to interpret outlying loci that may prove more widespread than anticipated.     

Conservation of Nuclear SSR Loci Reveals High Affinity of <Emphasis Type="Italic">Quercus infectoria</Emphasis> ssp. <Emphasis Type="Italic">veneris</Emphasis> A. Kern (Fagaceae) to Section <Emphasis Type="Italic">Robur</Emphasis>

Conservation of 16 nuclear microsatellite loci, originally developed for Quercus macrocarpa (section Albae), Q. petraea, Q. robur (section Robur), and Q. myrsinifolia, (subgenus Cyclobalanopsis) was tested in a Q. infectoria ssp. veneris population from Cyprus. All loci could be amplified successfully and displayed allele size and diversity patterns that match those of oak species belonging to the section Robur. At least in one case, limited amplification and high levels of homozygosity support the occurrence of “null alleles” caused by a possible mutation in the highly conserved primer areas, thus hindering PCR. The sampled population exhibited high levels of diversity despite the very limited distribution of this species in Cyprus and extended population fragmentation. Allele sizes of Q. infectoria at locus QpZAG9 partially match those of Q. alnifolia and Q. coccifera from neighboring populations. However, sequencing showed homoplasy, excluding a case of interspecific introgression with the latter, phylogenetically remote species. Q. infectoria ssp. veneris sequences at this locus were concordant to those of other species of section Robur, while sequences of Quercus alnifolia and Quercus coccifera were almost identical to Q. cerris.     

Evidence for hybridization and introgression within a species-rich oak (<Emphasis Type="Italic">Quercus</Emphasis> spp.) community

Background  

Analysis of interspecific gene flow is crucial for the understanding of speciation processes and maintenance of species integrity. Oaks (genus Quercus, Fagaceae) are among the model species for the study of hybridization. Natural co-occurrence of four closely related oak species is a very rare case in the temperate forests of Europe. We used both morphological characters and genetic markers to characterize hybridization in a natural community situated in west-central Romania and which consists of Quercus robur, Q. petraea, Q. pubescen s, and Q. frainetto, respectively.     

The role of local ecology during hybridization at the initial stages of ecological speciation in a marine snail

Hybrid zones of ecologically divergent populations are ideal systems to study the interaction between natural selection and gene flow during the initial stages of speciation. Here, we perform an amplified fragment length polymorphism (AFLP) genome scan in parallel hybrid zones between divergent ecotypes of the marine snail Littorina saxatilis, which is considered a model case for the study of ecological speciation. Ridged‐Banded (RB) and Smooth‐Unbanded (SU) ecotypes are adapted to different shore levels and microhabitats, although they present a sympatric distribution at the mid‐shore where they meet and mate (partially assortatively). We used shell morphology, outlier and nonoutlier AFLP loci from RB, SU and hybrid specimens captured in sympatry to determine the level of phenotypic and genetic introgression. We found different levels of introgression at parallel hybrid zones and nonoutlier loci showed more gene flow with greater phenotypic introgression. These results were independent from the phylogeography of the studied populations, but not from the local ecological conditions. Genetic variation at outlier loci was highly correlated with phenotypic variation. In addition, we used the relationship between genetic and phenotypic variation to estimate the heritability of morphological traits and to identify potential Quantitative Trait Loci to be confirmed in future crosses. These results suggest that ecology (exogenous selection) plays an important role in this hybrid zone. Thus, ecologically based divergent natural selection is responsible, simultaneously, for both ecotype divergence and hybridization. On the other hand, genetic introgression occurs only at neutral loci (nonoutliers). In the future, genome‐wide studies and controlled crosses would give more valuable information about this process of speciation in the face of gene flow.     

Hybridisation and colonisation dynamics in European oaks

Summary

Extensive hybridisation between the two sympatric species Quercus petraea and Q.robur is suggested by the near lack of genetic differentiation between the two species and supported by controlled crosses and mating system analysis in mixed stands. Further ecological and genetic evidence suggest that hybridisation does not impede the ecological specialisation of the two species, raising the issue of its evolutionary significance in oaks. Preferential unidirectional hybridization (pollen Q. petraea to ovule Q. robur) has been shown in various mixed stands and facilitates the introduction of sessile oak in existing pedunculate stands. If this unidirectional trend is reinforced in later backcrosses, then hybridisation leads to the dispersal of Q. petraea in existing stands of Q. robur. Hybridisation can therefore be seen as a ‘pollen-mediated’ dispersal mechanism, and has most likely contributed to the rapid migration of Q. petraea in Europe. Given the extant distribution of the species in Europe, migration through pollen swamping should be seen at the edges of the natural distribution of Q. petraea where the demographic imbalance of the two species will reinforce backcrosses.     

Bayesian clustering analyses for genetic assignment and study of hybridization in oaks: effects of asymmetric phylogenies and asymmetric sampling schemes

Bayesian clustering methods have been widely used for studying species delimitation and genetic introgression. In order to test the effect of phylogenetic relationships and sampling scheme on the inferred clustering solution and on the performance of Bayesian clustering analysis, I simulated genotypes of the interfertile oak species Quercus robur, Quercus petraea, and Quercus pubescens and I run analyses using two popular software programs, STRUCTURE and BAPS. First, based on purebred simulations, I compared clustering solutions resulting from different sample size configurations. While clustering solution generally reflected the taxonomic relationships when equal samples of each species were included, spurious partition was inferred by STRUCTURE when some species were represented by larger and others by smaller samples. In very unbalanced configurations, STRUCTURE failed to identify the three species, even if three subpopulations were assumed. By contrast, BAPS could properly identify the three species under any sampling scheme. Second, based on simulations of purebreds and hybrids, I tested the performance of individual assignments with variable number of loci. This analysis showed that STRUCTURE can detect introgressed individuals more efficiently than BAPS. However, BAPS could assign purebreds more efficiently with a lower number of loci. Method performance also depended on phylogenetic relationships. In the case of Q. petraea, Q. pubescens, and their hybrids, method performance was lower due to their phylogenetic affinity. Inclusion of three instead of two species into the analysis led to reduction of performance, and to misclassification of hybrids, which often reflected the phylogenetic affinity between Q. petraea and Q. pubescens.     

Population differentiation of sessile oak at the altitudinal front of migration in the French Pyrenees

To assess the effects of altitude on the level and structure of genetic diversity, a genetic survey was conducted in 12 populations of sessile oak (Quercus petraea) located between 130 and 1660 m in two parallel valleys on the northern side of the Pyrenees Mountains. Genetic diversity was monitored at 16 nuclear microsatellite loci and 5 chloroplast DNA (cpDNA) markers. The cpDNA survey suggested that extant populations in both valleys shared the same source populations from the plain. There was no visible trend of nuclear genetic diversity along altitude, even if indirect estimates of effective population sizes revealed a consistent reduction at higher altitudes. Population differentiation, although low, was mostly present among populations of the same valleys and reached similar levels than differentiation across the range of distribution of sessile oak. Contribution to the overall differentiation in the valleys was mostly due to the genetic divergence of the highest populations and the altitudinal variation of allelic frequencies at a few loci. Bayesian inference of migration between groups of populations showed that gene flow is preferentially unidirectional from lower altitudes in one valley to other groups of populations. Finally, we found evidence of clonal reproduction in high altitude populations. The introgression of Quercus robur and Quercus pubescens was also more frequent at the altitudinal margin suggesting that this mechanism may have contributed to the present migration and adaptation of Q. petraea and may also facilitate its future upslope shift in the context of climate change.     

No heat-stimulated germination found in herbaceous species from burned subtropical grassland

The inverse relationship between numbers of stomata (stomatal frequency) on tree leaves and ambient CO₂ concentration is increasingly applied for reconstructing past atmospheric CO₂ levels. The abundance of leaf remains of Quercus robur in Holocene peat and lake deposits in Europe makes this species potentially suitable for high-resolution stomatal frequency analysis. In order to quantify the CO₂ responsiveness of the species, the behavior of the stomatal index for Q. robur during the current anthropogenic CO₂ increase is determined on the basis of buried, herbarium and modern leaf material from the Netherlands. The stomatal index (SI), expressing the ratio of the number of stomata in a given area divided by the total number of stomata and other epidermal cells in that same area, is used in order to minimize influences on stomatal frequency of environmental conditions other than CO₂. The sigmoid SI response pattern recorded for Q. robur resembles that of the closely related species Q. petraea, although there is a difference in the timing of the response limitation of the two species to increasing atmospheric CO₂. For calibration purposes only the linear phase of the sigmoidal response curve is taken into consideration in the presented CO₂ response model, which allows confident combination of Q. robur and Q. petraea over the interval from 290 to 325 ppmv CO₂. The model is conservative in reconstructing past CO₂ mixing ratios outside the range of monitored response. As a result of the observed SI response limit, the model predicts CO₂ levels below 325 ppmv with a mean error of 10.2 ppmv, whereas higher CO₂ levels are underestimated.     

Single‐nucleotide polymorphism discovery and validation in high‐density SNP array for genetic analysis in European white oaks

An Illumina Infinium SNP genotyping array was constructed for European white oaks. Six individuals of Quercus petraea and Q. robur were considered for SNP discovery using both previously obtained Sanger sequences across 676 gene regions (1371 in vitro SNPs) and Roche 454 technology sequences from 5112 contigs (6542 putative in silico SNPs). The 7913 SNPs were genotyped across the six parental individuals, full‐sib progenies (one within each species and two interspecific crosses between Q. petraea and Q. robur) and three natural populations from south‐western France that included two additional interfertile white oak species (Q. pubescens and Q. pyrenaica). The genotyping success rate in mapping populations was 80.4% overall and 72.4% for polymorphic SNPs. In natural populations, these figures were lower (54.8% and 51.9%, respectively). Illumina genotype clusters with compression (shift of clusters on the normalized x‐axis) were detected in ~25% of the successfully genotyped SNPs and may be due to the presence of paralogues. Compressed clusters were significantly more frequent for SNPs showing a priori incorrect Illumina genotypes, suggesting that they should be considered with caution or discarded. Altogether, these results show a high experimental error rate for the Infinium array (between 15% and 20% of SNPs potentially unreliable and 10% when excluding all compressed clusters), and recommendations are proposed when applying this type of high‐throughput technique. Finally, results on diversity levels and shared polymorphisms across targeted white oaks and more distant species of the Quercus genus are discussed, and perspectives for future comparative studies are proposed.     

Tracking the progression of speciation: variable patterns of introgression across the genome provide insights on the species delimitation between progenitor–derivative spruces (Picea mariana × P. rubens)

The genic species concept implies that while most of the genome can be exchanged somewhat freely between species through introgression, some genomic regions remain impermeable to interspecific gene flow. Hence, interspecific differences can be maintained despite ongoing gene exchange within contact zones. This study assessed the heterogeneous patterns of introgression at gene loci across the hybrid zone of an incipient progenitor–derivative species pair, Picea mariana (black spruce) and Picea rubens (red spruce). The spruce taxa likely diverged in geographic isolation during the Pleistocene and came into secondary contact during late Holocene. A total of 300 SNPs distributed across the 12 linkage groups (LG) of black spruce were genotyped for 385 individual trees from 33 populations distributed across the allopatric zone of each species and within the zone of sympatry. An integrative framework combining three population genomic approaches was used to scan the genomes, revealing heterogeneous patterns of introgression. A total of 23 SNPs scattered over 10 LG were considered impermeable to introgression and putatively under diverging selection. These loci revealed the existence of impermeable genomic regions forming the species boundary and are thus indicative of ongoing speciation between these two genetic lineages. Another 238 SNPs reflected selectively neutral diffusion across the porous species barrier. Finally, 39 highly permeable SNPs suggested ancestral polymorphism along with balancing selection. The heterogeneous patterns of introgression across the genome indicated that the speciation process between black spruce and red spruce is young and incomplete, albeit some interspecific differences are maintained, allowing ongoing species divergence even in sympatry. The approach developed in this study can be used to track the progression of ongoing speciation processes.     

Pure species in a continuum of genetic and morphological variation: sympatric oaks at the edge of their range

Background and Aims Studies on oaks (Quercus spp.) have often been hampered by taxonomic confusion, a situation further compounded by the occurrence of extensive interspecific hybridization. In the present study, a combination of genetic and morphological analyses was used to examine sympatric populations of Q. petraea and Q. robur at the north-western edge of their ranges in Northern Ireland, since it had previously been suggested that hybridization could facilitate the apparent rapid, long-distance dispersal of oaks following the glaciations.Methods Samples were collected from 24 sites across Northern Ireland that had been previously designated as ancient or semi-natural woodland. Genotypes were obtained from a total of 950 trees using 12 nuclear microsatellite loci, and admixture coefficients were calculated based on a Bayesian clustering approach. Individuals were also classified as Q. petraea, Q. robur or hybrids based on two objective morphometric characters shown previously to delineate pure individuals effectively. Genetically ‘pure’ individuals of both species, as defined by the Bayesian clustering, were also genotyped for five chloroplast microsatellites.Key Results Genetic and morphological analyses both indicated the presence of pure individuals of both species, as well as a continuum of intermediates. There was a good agreement between the molecular and morphological classification, with a generally clear separation between pure individuals.Conclusions Despite millennia of hybridization and introgression, genetically and morphologically pure individuals of both Q. petraea and Q. robur can be found at the edge of their range, where both species occur sympatrically. The high proportion of individuals exhibiting introgression compared with previous studies may reflect the historical role of hybridization in facilitating dispersal following the glaciations. This is further supported by the significantly higher chloroplast diversity in Q. robur compared with Q. petraea.     

Gene flow,divergent selection and resistance to introgression in two species of morning glories (Ipomoea)

Gene flow is thought to impede genetic divergence and speciation by homogenizing genomes. Recent theory and research suggest that sufficiently strong divergent selection can overpower gene flow, leading to loci that are highly differentiated compared to others. However, there are also alternative explanations for this pattern. Independent evidence that loci in highly differentiated regions are under divergent selection would allow these explanations to be distinguished, but such evidence is scarce. Here, we present multiple lines of evidence that many of the highly divergent SNPs in a pair of sister morning glory species, Ipomoea cordatotriloba and I. lacunosa, are the result of divergent selection in the face of gene flow. We analysed a SNP data set across the genome to assess the amount of gene flow, resistance to introgression and patterns of selection on loci resistant to introgression. We show that differentiation between the two species is much lower in sympatry than in allopatry, consistent with interspecific gene flow in sympatry. Gene flow appears to be substantially greater from I. lacunosa to I. cordatotriloba than in the reverse direction, resulting in sympatric and allopatric I. cordatotriloba being substantially more different than sympatric and allopatric I. lacunosa. Many SNPs highly differentiated in allopatry have experienced divergent selection, and, despite gene flow in sympatry, resist homogenization in sympatry. Finally, five out of eight floral and inflorescence characteristics measured exhibit asymmetric convergence in sympatry. Consistent with the pattern of gene flow, I. cordatotriloba traits become much more like those of I. lacunosa than the reverse. Our investigation reveals the complex interplay between selection and gene flow that can occur during the early stages of speciation.