Leaf morphological analyses in four European oak species (Quercus) and their hybrids: A comparison of traditional and geometric morphometric methods

Abstract

In this study, leaf morphology was assessed in a mixed oak stand (western France) using two geometric morphometric (landmark and outline) datasets and one dataset of 19 leaf measures. Adult oaks (817 oaks), comprising four white oak species (Quercus petraea, Q. robur, Q. pubescens and Q. pyrenaica), were sampled for DNA extraction and genetic analysis (nuclear microsatellites). Leaf morphology was assessed on 336 oaks, comprising pure species and hybrids as determined by genetic assignment. This comparative study of oak leaf morphology, based on the use of two free size geometric morphometric methods and a set of leaf measurements, combined with the genetic assignment of individuals to pure species or hybrids, provided information about the differences among species and the intermediate leaf morphology of their hybrids.     

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.     

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.     

Changes in the Dynamics of Foliar N Metabolites in Oak Saplings by Drought and Air Warming Depend on Species and Soil Type

Climate change poses direct or indirect influences on physiological mechanisms in plants. In particular, long living plants like trees have to cope with the predicted climate changes (i.e. drought and air warming) during their life span. The present study aimed to quantify the consequences of simulated climate change for foliar N metabolites over a drought-rewetting-drought course. Saplings of three Central European oak species (i.e. Quercus robur, Q. petraea, Q. pubescens) were tested on two different soil types (i.e. acidic and calcareous). Consecutive drought periods increased foliar amino acid-N and soluble protein-N concentrations at the expense of structural N in all three oak species. In addition, transient effects on foliar metabolite dynamics were observed over the drought-rewetting-drought course. The lowest levels of foliar soluble protein-N, amino acid-N and potassium cation with a minor response to drought and air warming were found in the oak species originating from the driest/warmest habitat (Q. pubescens) compared to Q. robur and Q. petraea. Higher foliar osmolyte-N and potassium under drought and air warming were observed in all oak species when grown on calcareous versus acidic soil. These results indicate that species-specific differences in physiological mechanisms to compensate drought and elevated temperature are modified by soil acidity.     

Extent of nuclear genome sharing among white oak species (Quercus L. subgen. Lepidobalanus (Endl.) Oerst.) in Slovakia estimated by allozymes

European white oak species (Quercus subg. Lepidobalanus) are known to share a substantial part of the nuclear and cytoplasmic genome. In spite of this fact, they have preserved their integrity. We analyzed genetic variation at thirteen isozyme loci in nine Quercus robur, ten Q. petraea, two Q. pubescens and two mixed (two-species) populations in Slovakia to assess the extent of allele sharing. Moreover, a detailed morphometric study based on thirteen leaf traits was done in the mixed stands. A Bayesian analysis of population structure in the set of pure populations revealed the existence of three clusters, which coincided with the supposed taxonomical identity of populations. However, the distribution of posterior probabilities of data was bimodal with a second peak at six groups, what indicates a possible substructure in the gene pools of the investigated oak species. Moreover, gene admixture was close to 50% in all investigated populations. Morphometric analysis of the two mixed stands revealed a clear morphological separation of the studied oak species, sharply contrasting with a high degree of gene admixture. Differentiation is not uniformly distributed over the genome. High levels of differentiation (F _ST > 0.10) were observed at only three out of thirteen loci. The hypotheses explaining extensive allele sharing (shared ancestral polymorphisms vs. high levels of interspecific gene flow) are briefly discussed.     

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.     

Osmotic stress responses of individual white oak (Quercus section, Quercus subgenus) genotypes cultured in vitro

White oaks (Quercus section, Quercus subgenus) are widely distributed in Europe. Quercus petraea (sessile oak), an economically important species is predicted to be affected by climate change. Q. pubescens (pubescent oak) and Q. virgiliana (Italian pubescent oak) are economically less important, drought tolerant species. Frequent hybridization of white oaks was observed and currently the introgression of Q. pubescens and Q. virgiliana in non-mediterranean regions of Europe has been reported. Our goal was to use tissue cultures established from individual trees of the above taxa and their putative hybrids, all present in the forest stand of Síkf?kút LTER Research Area (NE Hungary) as simple experimental model systems for studying drought/osmotic stress tolerance. Tissue cultures are more suitable models for such studies, than seedlings, because they are genetically identical to the parent plants. Polyethylene glycol (PEG6000) treatments were used for this purpose. The identification of taxa was based on leaf morphological traits and microsatellite analysis and showed that Q. petraea is genetically distinct to all other taxa examined. We established six callus lines of Quercus. As expected, in Q. petraea cultures PEG6000 induced severe loss of fresh weight and the ability to recover after removal of the osmoticum, which was not characteristic for Q. pubescens and Q. virgiliana. Putative hybrids exhibited an intermediate response to osmotic stress. Activity gels showed the increase of single-strand preferring (SSP) nuclease and no significant change of guaiacol-peroxidase activities in drought-sensitive genotypes/cultures and no significant increase of SSP nuclease activities accompanied with increases of guaiacol-peroxidase activities in drought-tolerant ones. This indicates that drought/osmotic stress tolerance is associated to increased capacity of scavenging reactive oxygen species and hence less susceptibility to DNA damage. Our results confirm that tissue cultures of oak are suitable model systems for studying drought/osmotic stress responses.     

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.     

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.     

Development and Characterization of Genomic and Gene-Based Microsatellite Markers in North American Red Oak Species

Oaks (Quercus: Fagaceae) are ecological and economic keystones of many forested ecosystems but effective genetic management strategies are hindered by high levels of phenotypic plasticity within species and frequent hybridization among them. These same features, however, make oak communities suited for the study of speciation, hybridization, and genetic adaptation. Efforts to develop new and to adapt existing genomic resources to less-studied members of this genus should not only improve oak conservation and management but also aid the study of fundamental evolutionary processes. Here, we present a suite of 27 highly polymorphic simple sequence repeat (SSR) markers tested in four North American red oak (Quercus section Lobatae) species: Q. rubra, Q. ellipsoidalis, Q. coccinea, and Q. velutina. Five markers are genomic SSRs (gSSRs) — four novel and one previously transferred from Q. petraea — and 22 are gene-based SSRs derived from Q. robur and Q. petraea expressed sequence tags (EST-SSRs). Overall, levels of polymorphism detected with these primer pairs were high, with gene diversity (H _e) averaging 0.66 across all loci in natural populations. In addition, we show that EST-SSR markers may have the potential to detect divergent selection at stress-resistance candidate genes among closely related oak species.     

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.     

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.     

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.     

Utility of multilocus genotypes for taxon assignment in stands of closely related European white oaks from Switzerland

Background and Aims

European white oaks (Quercus petraea, Q. pubescens, Q. robur) have long puzzled plant biologists owing to disputed species differentiation. Extensive hybridization or shared ancestry have been proposed as alternative hypotheses to explain why genetic differentiation between these oak species is low. Species delimitation is usually weak and often shows gradual transitions in leaf morphology. Hence, individual identification may be difficult, but remains a critical step for both scientific work and practical management.

Methods

Multilocus genotype data (five nuclear microsatellites) were used from ten Swiss oak stands for taxon identification without a priori grouping of individuals or populations, using model-based Bayesian assignment tests.

Key Results

Three groups best structured the data, indicating that the taxonomical signal was stronger than the spatial signal. Most individuals showed high posterior probabilities for either of three genetic groups that were best circumscribed as taxonomical units. The assignment of a subset of trees, whose taxonomic status had been previously characterized in detail, supported this classification scheme.

Conclusions

Molecular-genetic assignment tests are useful in the identification of species status in critical taxon complexes such as the European white oaks. Such an approach is of practical importance for forest management, e.g. for stand certification or in seed trade to trace the origin of forest products.Key words: Assignment test, Bayesian inference, multilocus genotype, nuclear microsatellites, Quercus sp., species complex     

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.     

<Emphasis Type="Italic">Sporothrix inflata</Emphasis>, a root-inhabiting fungus of <Emphasis Type="Italic">Quercus robur</Emphasis> and <Emphasis Type="Italic">Q. petraea</Emphasis>

The anamorphic fungus Sporothrix inflata, known as a soil-borne fungus with worldwide distribution, was isolated for the first time from the cortex and central cylinder of living and dead roots of healthy and diseased oak trees (Quercus robur and Q. petraea). Isolation frequencies of S. inflata from oak roots varied according to the health status of trees, oak species, study sites, soil depth and root diameter. Colony morphology and growth rate of isolates are influenced by colony age and type of culture medium.     

Role of micromorphological leaf traits and molecular data in taxonomy of three sympatric white oak species and their hybrids (Quercus L.)

Hybridisation and introgression occur with high frequency in the genus Quercus and interspecific hybrid individuals show patterns of morphological traits that might be influenced in different ways. Micromorphological leaf traits appear to be positive and stable in Quercus species, and by combining genetic and micromorphological analyses, it is possible to compare the patterns of variation in micromorphological leaf traits of pure and hybrid individuals. Trichomes and stomatal traits were examined using scanning electron microscopy at 150–2000 × magnification in sympatric oak species collected in a natural deciduous wood. Q. frainetto, Q. petraea and Q. pubescens appear to have a relatively predictable complement of trichome types. Both the pattern and quantitative values of each micromorphological trait examined (stomata and trichomes) have an important role in identifying hybrids and pure species; putative hybrids show a pattern of trichomes that is a combination of the parental types. These results, combined with the fact that micromorphological traits generally exhibit higher consistency, indicate that this source of information can be an excellent clue to hybridisation and introgression and useful in taxonomical, systematic and evolutionary studies on the European white oaks.     

Strength and variability of postmating reproductive isolating barriers between four European white oak species

The identification and quantification of the relative importance of reproductive isolating barriers is of fundamental importance to understand species maintenance in the face of interspecific gene flow between hybridising species. Yet, such assessments require extensive experimental fertilisations that are particularly difficult when dealing with more than two hybridising and long-generation-time species such as oaks. Here, we quantify the relative contribution of four postmating reproductive isolating barriers consisting of two prezygotic barriers (gametic incompatibility, conspecific pollen precedence) and two postzygotic barriers (germination rate, early survival) from extensively controlled pollinations between four oak species (Quercus robur, Quercus petraea, Quercus pubescens and Quercus pyrenaica) that have been shown to frequently hybridise in natural populations. We found high variation in the strength of total reproductive isolation between species, ranging from total reproductive isolation to advantage toward hybrid formation. As previously found, Q. robur pollen was unable to fertilise Q. petraea due to a strong reproductive isolating mechanism. On the contrary, Q. pubescens pollen was more efficient at fertilising Q. petraea than conspecific pollen. Overall, prezygotic barriers contribute far more than postzygotic barriers to isolate species reproductively, suggesting a role for reinforcement in the development of prezygotic barriers. Conspecific pollen precedence reduced hybrid formation when pollen competition was allowed; however, presence of conspecific pollen did not totally prevent hybridization. Our results suggest that pollen competition depends on multiple ecological and environmental parameters, including species local abundance, and that it may be of uppermost importance to understand interspecific gene flow among natural multispecies populations.     

Ectomycorrhizal fungal communities of pedunculate and sessile oak seedlings from bare-root forest nurseries

In this study, we present the detailed molecular investigation of the ectomycorrhizal (ECM) community of Quercus petraea and Quercus robur seedlings grown in bare-root forest nurseries. In all tested oak samples, mycorrhizal colonization was nearly 100%. Morphological observation and molecular investigations (sequencing of fungal ITS rDNA) revealed a total of 23 mycorrhizal taxa. The most frequent and abundant fungal taxa were Hebeloma sacchariolens, Tuber sp., and Peziza sp.; from the detected fungal taxa, 20 were noted for Q. petraea and 23 for Q. robur. Depending on the nursery, the species richness of identified ECM fungal taxa for both oak species ranged from six to 11 taxa. The mean species richness for all nurseries was 5.36 and 5.82 taxa per Q. petraea and Q. robur sample, respectively. According to the analysis of similarity, ECM fungal communities were similar for Q. petraea and Q. robur (R = 0.019; p = 0.151). On the other hand, detected fungal communities were significantly different between nurseries (R = 0.927; p < 0.0001). Using the Spearman rank correlation, it was determined that the ectomycorrhizal diversity (in terms of richness, the Shannon diversity, evenness, and Simpson dominance indices) is significantly related to the soil parameters of each nursery. We conclude that individual nursery may be considered as separate ecological niches that strongly discriminate diversity of ECM fungi.     

A set of primers for the amplification of chloroplast microsatellites in Quercus

The increase in demand for the certification of oak seed lots, as well as control of the geographical origin of oak wood, has led us to develop powerful genetic markers permitting us to discriminate among provenance regions. With the aim of detecting new chloroplast variants, we have identified 17 potential cpSSRs motifs from available oak sequences and tested their variability among French oak populations. Six loci were polymorphic at the intraspecific level in Quercus petraea and Q. robur. Moreover, conservation of the primer pairs was checked on a set of 21 forest tree species and they were all shown to work well on several Quercus species, and even within Fagacaea.