Evidence for hybridization and introgression within a species-rich oak (Quercus spp.) community

Hybridization has localized effect on genetic variation in closely related pine species

BACKGROUND

Hybridization is a known phenomenon in nature but its genetic impact on populations of parental species remains less understood. We investigated the evolutionary consequences of the interspecific gene flow in several contact zones of closely related pine species. Using a set of genetic markers from both nuclear and organellar genomes, we analyzed four hybrid zones (384 individuals) and a large panel of reference allopatric populations of parental taxa (2104 individuals from 96 stands).

RESULTS

We observed reduced genetic diversity in maternally transmitted mitochondrial genomes of pure pine species and hybrids from contact zones compared to reference allopatric populations. The distribution of mtDNA haplotypes followed geographic rather than species boundaries. Additionally, no new haplotypes emerged in the contact zones, instead these zones contained the most common local variants. However, species diverged significantly at nuclear genomes and populations in contact zones exhibited similar or higher genetic diversity compared to the reference stands. There were no signs of admixture in any allopatric population, while clear admixture was evident in the contact zones, indicating that hybridization has a geographically localized effect on the genetic variation of the analyzed pine species.

CONCLUSIONS

Our results suggest that hybrid zones act as sinks rather than melting pots of genetic diversity. Hybridization influences sympatric populations but is confined to contact zones. The spectrum of parental species ancestry in hybrids reflects the old evolutionary history of the sympatric populations. These findings also imply that introgression may play a crucial role in the adaptation of hybrids to specific environments.

The adaptive evolution of Quercus section Ilex using the chloroplast genomes of two threatened species

Chloroplast (cp) genome sequences have been extensively used for phylogenetic and evolutionary analyses, as many have been sequenced in recent years. Identification of Quercus is challenging because many species overlap phenotypically owing to interspecific hybridization, introgression, and incomplete lineage sorting. Therefore, we wanted to gain a better understanding of this genus at the level of the maternally inherited chloroplast genome. Here, we sequenced, assembled, and annotated the cp genomes of the threatened Quercus marlipoensis (160,995 bp) and Q. kingiana (161,167 bp), and mined these genomes for repeat sequences and codon usage bias. Comparative genomic analyses, phylogenomics, and selection pressure analysis were also performed in these two threatened species along with other species of Quercus. We found that the guanine and cytosine content of the two cp genomes were similar. All 131 annotated genes, including 86 protein-coding genes, 37 transfer RNA genes, and 8 ribosomal RNA genes, had the same order in the two species. A strong A/T bias was detected in the base composition of simple sequence repeats. Among the 59 synonymous codons, the codon usage pattern of the cp genomes in these two species was more inclined toward the A/U ending. Comparative genomic analyses indicated that the cp genomes of Quercus section Ilex are highly conserved. We detected eight highly variable regions that could be used as molecular markers for species identification. The cp genome structure was consistent and different within and among the sections of Quercus. The phylogenetic analysis showed that section Ilex was not monophyletic and was divided into two groups, which were respectively nested with section Cerris and section Cyclobalanopsis. The two threatened species sequenced in this study were grouped into the section Cyclobalanopsis. In conclusion, the analyses of cp genomes of Q. marlipoensis and Q. kingiana promote further study of the taxonomy, phylogeny and evolution of these two threatened species and Quercus.

Characterization of pollen tube development in distant hybridization of Chinese cork oak (Quercus variabilis L.)

MAIN CONCLUSION

This work mainly found that the stigma and style of Q. variabilis did not completely lose the specific recognition towards heterologous pollen, a fact which is different from previous studies. Quercus is the foundation species in the Northern Hemisphere, with extreme prevalence for interspecific hybridization. It is not yet entirely understood whether or how the pollen tube-female tissue interaction contributes to the "extensive hybridization" in oaks. Pollen storage conditions correlate with distant hybridization. We conducted hybridization experiments with Q. variabilis as female and Q. variabilis and Q. mongolica as male parents. And the differences in pollen tube (PT) development between intra- and distant interspecific hybridization were studied by fluorescence microscopy and scanning electron microscopy (SEM). Our results showed that -20 °C allowed pollen of both species to maintain some viability. Both Q. variabilis and Q. mongolica pollen germinated profusely on the stigmas. SEM results indicated that in the intraspecific hybridization, Q. variabilis pollen started to germinate at 6 h after pollination (hap), PTs elongated significantly at 12 hap, and entered the stigma at 24 hap. By contrast, Q. mongolica pollen germinated at 15 hap, and the PTs entered the stigma at 27 hap. By fluorescence microscopical studies it was observed that some PTs of Q. variabilis gathered at the style-joining at 96 hap, unlike the Q. mongolica which reached the style junction at 144 hap. The above results indicate that the abundant germination of heterologous pollen (HP) on the stigma and the "Feeble specificity recognition" of the stigma and transmitting tract to HP may create opportunities for the "extensive hybridization" of oaks. This work provides a sexual developmental reference for clarifying the causes of Quercus "extensive hybridization".

Comparison of chloroplast genomes and phylogenetic analysis of four species in Quercus section Cyclobalanopsis

The identification in Quercus L. species was considered to be difficult all the time. The fundamental phylogenies of Quercus have already been discussed by morphological and molecular means. However, the morphological characteristics of some Quercus groups may not be consistent with the molecular results (such as the group Helferiana), which may lead to blurring of species relationships and prevent further evolutionary researches. To understand the interspecific relationships and phylogenetic positions, we sequenced and assembled the CPGs (160,715 bp-160842 bp) of four Quercus section Cyclobalanopsis species by Illumina pair-end sequencing. The genomic structure, GC content, and IR/SC boundaries exhibited significant conservatism. Six highly variable hotspots were detected in comparison analysis, among which rpoC1, clpP and ycf1 could be used as molecular markers. Besides, two genes (petA, ycf2) were detected to be under positive selection pressure. The phylogenetic analysis showed: Trigonobalanus genus and Fagus genus located at the base of the phylogeny tree; The Quercus genus species were distincted to two clades, including five sections. All Compound Trichome Base species clustered into a single branch, which was in accordance with the results of the morphological studies. But neither of group Gilva nor group Helferiana had formed a monophyly. Six Compound Trichome Base species gathered together in pairs to form three branch respectively (Quercus kerrii and Quercus chungii; Quercus austrocochinchinensis with Quercus gilva; Quercus helferiana and Quercus rex). Due to a low support rate (0.338) in the phylogeny tree, the interspecies relationship between the two branches differentiated by this node remained unclear. We believe that Q. helferiana and Q. kerrii can exist as independent species due to their distance in the phylogeny tree. Our study provided genetic information in Quercus genus, which could be applied to further studies in taxonomy and phylogenetics.

Sequential hybridization may have facilitated ecological transitions in the Southwestern pinyon pine syngameon

Multispecies interbreeding networks, or syngameons, have been increasingly reported in natural systems. However, the formation, structure, and maintenance of syngameons have received little attention. Through gene flow, syngameons can increase genetic diversity, facilitate the colonization of new environments, and contribute to hybrid speciation. In this study, we evaluated the history, patterns, and consequences of hybridization in a pinyon pine syngameon using morphological and genomic data to assess genetic structure, demographic history, and geographic and climatic data to determine niche differentiation. We demonstrated that Pinus edulis, a dominant species in the Southwestern US and a barometer of climate change, is a core participant in the syngameon, involved in the formation of two drought-adapted hybrid lineages including the parapatric and taxonomically controversial fallax-type. We found that species remain morphologically and genetically distinct at range cores, maintaining species boundaries while undergoing extensive gene flow in areas of sympatry at range peripheries. Our study shows that sequential hybridization may have caused relatively rapid speciation and facilitated the colonization of different niches, resulting in the rapid formation of two new lineages. Participation in the syngameon may allow adaptive traits to be introgressed across species barriers and provide the changes needed to survive future climate scenarios.

A Small Set of Nuclear Markers for Reliable Differentiation of the Two Closely Related Oak Species Quercus Robur and Q. Petraea

Quercus robur and Q. petraea are, in addition to Fagus sylvatica, the main economically used deciduous tree species in Europe. Identification of these two species is crucial because they differ in their ecological demands. Because of a changing climate, foresters must know more than ever which species will perform better under given environmental conditions. The search for differentiating molecular markers between these two species has already lasted for decades. Until now, differentiation has only been possible in approaches with a combination of several molecular markers and a subsequent statistical analysis to calculate the probability of being one or the other species. Here, we used MiSeq Illumina data from pools of Q. robur and Q. petraea specimens and identified nuclear SNPs and small InDels versus the Q. robur reference genome. Selected sequence variants with 100% allele frequency difference between the two pools were further validated in an extended set of Q. robur and Q. petraea specimens, and then the number of markers was deliberately reduced to the smallest possible set for species differentiation. A combination of six markers from four nuclear regions is enough to identify Q. robur, Q. petraea or hybrids between these two species quite well and represents a marker set that is cost-efficient and useable in every laboratory.

Effect of Irrigation Dose on Powdery Mildew Incidence and Root Biomass of Sessile Oaks (Quercus petraea (Matt.) Liebl.)

The sessile oak is one of the most significant forest tree species in Europe. This species is vulnerable to various stresses, among which drought and powdery mildew have been the most serious threats. The aim of this study was to determine the influence of irrigation levels (overhead sprinklers) on the damage caused by powdery mildew to Quercus petraea growing in a nursery setting. Four irrigation rates were used: 100%, 75%, 50% and 25% of the full rate. The area of the leaves was measured and the ratio between the dry mass of the roots and the dry mass of the entire plant was calculated after the growing season in years’ 2015 and 2016. Limiting the total amount of water provided to a level between 53.6 mm × m⁻² and 83.6 mm × m⁻², particularly in the months when total precipitation was low (VII and VIII 2015), a supplemental irrigation rate between 3 and 9 mm × m⁻² resulted in a lower severity of oak powdery mildew on leaves and lead to a favorable allocation of the biomass of the sessile oak seedlings to the root system. The severity of infection on oak leaf blades was lower when irrigation rates were reduced. The greatest mean degree of infestation in 2015 was noted in the 100% irrigation rate (14.6%), 75% (6.25%), 50% (4.35%) and 25% (5.47%). In 2016, there was no significant difference between the mean area of leaves infected by powdery mildew depending on the applied irrigation rate. The shoot-root biomass rate showed greater variation under limited irrigation rates. Controlling the irrigation rate can become an effective component of integrated protection strategies against this pathogen.

Morphometric Analyses of Leaf Shapes in Four Sympatric Mediterranean Oaks and Hybrids in the Algerian Kabylie Forest

(1) Background: local morphological variation can provide useful information to clarify the role of hybridization in Mediterranean oaks. Accordingly, we have characterized putative hybrids and oak species coexisting in a highly diversified oak forest in Algeria with four native oak species (Quercus suber L., Q. ilex L., Q. canariensis Willd., and Q. afares Pomel). (2) Methods: sixteen plots, including 89 trees from the four native species and their putative hybrids were sampled. Leaves were scanned and their geometric morphometry analyzed by using 11 landmarks on the right side on their abaxial surface. Variation within and among species, and the relationship among oak species and their hybrids were analyzed, utilizing an ANOVA and a canonical analysis using morphoJ software. (3) Results: using the geometric morphometry analysis, we observed that Q. afares shape is intermediate between Q. suber and Q. canariensis, being Q. ilex very different from the others, and that there is no overlap among them. Putative hybrids are morphologically close to Q. afares and, to a lesser extent, to Q. suber and, finally, to Q. ilex and Q. canariensis. (4) Conclusions: the study opens the field for future molecular characterization of hybrids, and for determining their role in terms of adaptation to actual and predicted future climatic conditions. The morphological proximity of hybrids to Q. afares demonstrate this species advantage in the area of study and its importance for the future evolution of the species in the Mediterranean.

Assessment of Drought-Tolerant Provenances of Austria’s Indigenous Tree Species

(1) Background: Forestry will have to react to climate change because many tree species suffer. Mitigation can be realized either by planting non-native trees from regions with high climatic stress or by utilizing native tree provenances already adapted to stressful environments. Non-native trees have often generated problems in the past due to uncontrolled invasiveness. The use of native trees pre-adapted to the prospective climatic conditions is far less risky for the respective ecosystems. We offer a tool for selecting ecotypes of native trees as provenances for future forestry. (2) Methods: We propose the selection of tree species native to Middle Europe from a database of vegetation relevés of ± natural forest stands. By calculating the mean ecological indicator values of stands from their vegetation, cover sites can be elected that can provide seeds of provenances well adapted to future climatic conditions. (3) Results: By selecting the 10% partition of the most extreme stands of European tree species, seeds can be sampled and propagated for re-cultivating forests fit for future climate. (4) Conclusions: One can expect ecotypes of tree species that grow well on dry sites, since generations have faced evolutionary selection, for survival under stressful environments. This approach helps to avoid ecological risks of non-native trees.

Evidence for extensive hybridisation and past introgression events in feather grasses using genome-wide SNP genotyping

BACKGROUND

The proper identification of feather grasses in nature is often limited due to phenotypic variability and high morphological similarity between many species. Among plausible factors influencing this issue are hybridisation and introgression recently detected in the genus. Nonetheless, to date, only a bounded set of taxa have been investigated using integrative taxonomy combining morphological and molecular data. Here, we report the first large-scale study on five feather grass species across several hybrid zones in Russia and Central Asia. In total, 302 specimens were sampled in the field and classified based on the current descriptions of these taxa. They were then genotyped with high density genome-wide markers and measured based on a set of morphological characters to delimitate species and assess levels of hybridisation and introgression. Moreover, we tested species for past introgression and estimated divergence times between them.

RESULTS

Our findings demonstrated that 250 specimens represent five distinct species: S. baicalensis, S. capillata, S. glareosa, S. grandis and S. krylovii. The remaining 52 individuals provided evidence for extensive hybridisation between S. capillata and S. baicalensis, S. capillata and S. krylovii, S. baicalensis and S. krylovii, as well as to a lesser extent between S. grandis and S. krylovii, S. grandis and S. baicalensis. We detected past reticulation events between S. baicalensis, S. krylovii, S. grandis and inferred that diversification within species S. capillata, S. baicalensis, S. krylovii and S. grandis started ca. 130-96 kya. In addition, the assessment of genetic population structure revealed signs of contemporary gene flow between populations across species from the section Leiostipa, despite significant geographical distances between some of them. Lastly, we concluded that only 5 out of 52 hybrid taxa were properly identified solely based on morphology.

CONCLUSIONS

Our results support the hypothesis that hybridisation is an important mechanism driving evolution in Stipa. As an outcome, this phenomenon complicates identification of hybrid taxa in the field using morphological characters alone. Thus, integrative taxonomy seems to be the only reliable way to properly resolve the phylogenetic issue of Stipa. Moreover, we believe that feather grasses may be a suitable genus to study hybridisation and introgression events in nature.

Chromosome-scale genome assembly of Castanopsis tibetana provides a powerful comparative framework to study the evolution and adaptation of Fagaceae trees

Fagaceae species are increasingly used as models to elucidate the process and mechanism of adaptation and speciation by integrating ecology, evolution and genomics. The genus Castanopsis belongs to the family Fagaceae and is mainly distributed across subtropical and tropical Asia. In the present study, we reported the first chromosome-scale genome assembly of Castanopsis tibetana, a common species of evergreen broadleaved forests in subtropical China. The combination of Nanopore sequencing and Hi-C technologies enabled a high-quality genome assembly. The final assembled genome size of C. tibetana was 878.6 Mb (97.6% of the estimated genome size), consisting of 477 contigs with an N50 length of 3.3 Mb. The benchmarking universal single-copy orthologue (BUSCO) assessment indicated a completeness of 93.0%. Hi-C scaffolding generated 12 pseudochromosomes, representing 98.7% of the assembled genome. Subsequently, 40,937 protein-coding genes were predicted and 90.04% of them were functionally annotated. More than 476.9 Mb of repetitive sequences (54.3% of the genome) were identified, and the percentage of the genome covered by TE elements was 39.98%. Comparative genomics analysis revealed that C. tibetana was most closely related to Castanea mollissima and diverged at 18.48 Ma, and that C. tibetana has undergone considerable gene family expansion and contraction. Evidence of positive selection was detected in 53 genes, which showed different arrangement pattern compared to Quercus robur. The chromosome-scale genome assembly of C. tibetana will expand Fagaceae genome resources across the family and provide a powerful comparative framework to study the adaptation and evolution of Fagaceae trees.

Impact of Gene Flow and Introgression on the Range Wide Genetic Structure of Quercus robur (L.) in Europe

As for most other temperate broadleaved tree species, large-scale genetic inventories of pedunculate oak (Quercus robur L.) have focused on the plastidial genome, which showed the impact of post-glacial recolonization and manmade seed transfer. However, how have pollen mediated gene flow and introgression impacted the large-scale genetic structure? To answer these questions, we did a genetic inventory on 1970 pedunculate oak trees from 197 locations in 13 European countries. All samples were screened with a targeted sequencing approach on a set of 381 polymorphic loci (356 nuclear SNPs, 3 nuclear InDels, 17 chloroplast SNPs, and 5 mitochondrial SNPs). In a former analysis with additional 1763 putative Quercus petraea trees screened for the same gene markers we obtained estimates on the species admixture of all pedunculate oak trees. We identified 13 plastidial haplotypes, which showed a strong spatial pattern with a highly significant autocorrelation up to a range of 1250 km. Significant spatial genetic structure up to 1250 km was also observed at the nuclear loci. However, the differentiation at the nuclear gene markers was much lower compared to the organelle gene markers. The matrix of genetic distances among locations was partially correlated between nuclear and organelle genomes. Bayesian clustering analysis revealed the best fit to the data for a sub-division into two gene pools. One gene pool is dominating the west and the other is the most abundant in the east. The western gene pool was significantly influenced by introgression from Quercus petraea in the past. In Germany, we identified a contact zone of pedunculate oaks with different introgression intensity, likely resulting from different historical levels of introgression in glacial refugia or during postglacial recolonization. The main directions of postglacial recolonization were south to north and south to northwest in West and Central Europe, and for the eastern haplotypes also east to west in Central Europe. By contrast, the pollen mediated gene flow and introgression from Q. petraea modified the large-scale structure at the nuclear gene markers with significant west–east direction.

Evolutionary patterns of nucleotide substitution rates in plastid genomes of Quercus

Molecular evolution, including nucleotide substitutions, plays an important role in understanding the dynamics and mechanisms of species evolution. Here, we sequenced whole plastid genomes (plastomes) of Quercus fabri, Quercus semecarpifolia, Quercus engleriana, and Quercus phellos and compared them with 14 other Quercus plastomes to explore their evolutionary relationships using 67 shared protein-coding sequences. While many previously identified evolutionary relationships were found, our findings do not support previous research which retrieve Quercus subg. Cerris sect. Ilex as a monophyletic group, with sect. Ilex found to be polyphyletic and composed of three strongly supported lineages inserted between sections Cerris and Cyclobalanposis. Compared with gymnosperms, Quercus plastomes showed higher evolutionary rates (Dn/Ds = 0.3793). Most protein-coding genes experienced relaxed purifying selection, and the high Dn value (0.1927) indicated that gene functions adjusted to environmental changes effectively. Our findings suggest that gene interval regions play an important role in Quercus evolution. We detected greater variation in the intergenic regions (trnH-psbA, trnK_UUU-rps16, trnfM_CAU-rps14, trnS_GCU-trnG_GCC, and atpF-atpH), intron losses (petB and petD), and pseudogene loss and degradation (ycf15). Additionally, the loss of some genes suggested the existence of gene exchanges between plastid and nuclear genomes, which affects the evolutionary rate of the former. However, the connective mechanism between these two genomes is still unclear.

Variations in morphological and epidermal features of shade and sun leaves of two species: Quercus bambusifolia and Q. myrsinifolia

PREMISE

Microclimatic differences between the periphery and the interior of tree crowns result in a variety of adaptive leaf macromorphological and anatomical features. Our research was designed to reveal criteria for sun/shade leaf identification in two species of evergreen oaks, applicable to both modern and fossil leaves. We compared our results with those in other species similarly studied.

METHODS

For both Quercus bambusifolia and Q. myrsinifolia (section Cyclobalanopsis), leaves from single mature trees with well-developed crowns were collected in the South China Botanical Garden, Guangzhou, China. We focus on leaf characters often preserved in fossil material. SVGm software was used for macromorphological measurement. Quantitative analyses were performed and box plots generated using R software with IDE Rstudio. Leaf cuticles were prepared using traditional botanical techniques.

RESULTS

Principal characters for distinguishing shade and sun leaves in the studied oaks were identified as leaf lamina length to width ratio (L/W), and the degree of development of venation networks. For Q. myrsinifolia, shade and sun leaves differ in tooth morphology and the ratio of toothed lamina length to overall lamina length. The main epidermal characters are ordinary cell size and anticlinal wall outlines. For both species, plasticity within shade leaves exceeds that of sun leaves.

CONCLUSIONS

Morphological responses to sun and shade in the examined oaks are similar to those in other plant genera, pointing to useful generalizations for recognizing common foliar polymorphisms that must be taken into account when determining the taxonomic position of both modern and fossil plants.

Exaptation Traits for Megafaunal Mutualisms as a Factor in Plant Domestication

Megafaunal extinctions are recurring events that cause evolutionary ripples, as cascades of secondary extinctions and shifting selective pressures reshape ecosystems. Megafaunal browsers and grazers are major ecosystem engineers, they: keep woody vegetation suppressed; are nitrogen cyclers; and serve as seed dispersers. Most angiosperms possess sets of physiological traits that allow for the fixation of mutualisms with megafauna; some of these traits appear to serve as exaptation (preadaptation) features for farming. As an easily recognized example, fleshy fruits are, an exaptation to agriculture, as they evolved to recruit a non-human disperser. We hypothesize that the traits of rapid annual growth, self-compatibility, heavy investment in reproduction, high plasticity (wide reaction norms), and rapid evolvability were part of an adaptive syndrome for megafaunal seed dispersal. We review the evolutionary importance that megafauna had for crop and weed progenitors and discuss possible ramifications of their extinction on: (1) seed dispersal; (2) population dynamics; and (3) habitat loss. Humans replaced some of the ecological services that had been lost as a result of late Quaternary extinctions and drove rapid evolutionary change resulting in domestication.

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.

Oaks: an evolutionary success story

The genus Quercus is among the most widespread and species-rich tree genera in the northern hemisphere. The extraordinary species diversity in America and Asia together with the continuous continental distribution of a limited number of European species raise questions about how macro- and microevolutionary processes made the genus Quercus an evolutionary success. Synthesizing conclusions reached during the past three decades by complementary approaches in phylogenetics, phylogeography, genomics, ecology, paleobotany, population biology and quantitative genetics, this review aims to illuminate evolutionary processes leading to the radiation and expansion of oaks. From opposing scales of time and geography, we converge on four overarching explanations of evolutionary success in oaks: accumulation of large reservoirs of diversity within populations and species; ability for rapid migration contributing to ecological priority effects on lineage diversification; high rates of evolutionary divergence within clades combined with convergent solutions to ecological problems across clades; and propensity for hybridization, contributing to adaptive introgression and facilitating migration. Finally, we explore potential future research avenues, emphasizing the integration of microevolutionary and macroevolutionary perspectives.

Oak symbolism in the light of genomics

Throughout the Northern Hemisphere, human societies, political systems, and religions have appropriated oaks in symbolic representations. In this review, we explore the possible associations between recent genetic and genomic findings and the symbolic representations of oaks. We first consider the ways in which evolutionary history during the Holocene has tightened links between humans and oaks in Europe, and how this may have led to symbolic representations. We then show how recent findings concerning the structure and evolution of the oak genome have provided additional knowledge about symbolic representations, such as longevity, cohesiveness, and robustness.

Adaptive introgression as a driver of local adaptation to climate in European white oaks

Latitudinal and elevational gradients provide valuable experimental settings for studies of the potential impact of global warming on forest tree species. The availability of long-term phenological surveys in common garden experiments for traits associated with climate, such as bud flushing for sessile oaks (Quercus petraea), provide an ideal opportunity to investigate this impact. We sequenced 18 sessile oak populations and used available sequencing data for three other closely related European white oak species (Quercus pyrenaica, Quercus pubescens, and Quercus robur) to explore the evolutionary processes responsible for shaping the genetic variation across latitudinal and elevational gradients in extant sessile oaks. We used phenotypic surveys in common garden experiments and climatic data for the population of origin to perform genome-wide scans for population differentiation and genotype-environment and genotype-phenotype associations. The inferred historical relationships between Q. petraea populations suggest that interspecific gene flow occurred between Q. robur and Q. petraea populations from cooler or wetter areas. A genome-wide scan of differentiation between Q. petraea populations identified single nucleotide polymorphisms (SNPs) displaying strong interspecific relative divergence between these two species. These SNPs followed genetic clines along climatic or phenotypic gradients, providing further support for the likely contribution of introgression to the adaptive divergence of Q. petraea populations. Overall, the results indicate that outliers and associated SNPs are Q. robur ancestry-informative. We discuss the results of this study in the framework of the postglacial colonization scenario, in which introgression and diversifying selection have been proposed as essential drivers of Q. petraea microevolution.

A species-discriminatory single-nucleotide polymorphism set reveals maintenance of species integrity in hybridizing European white oaks (Quercus spp.) despite high levels of admixture

BACKGROUND AND AIMS

Hybridization and introgression play an important role in the evolution and diversification of plants. To assess the degree of past and current hybridization, the level of genetic admixture in populations needs to be investigated. Ongoing hybridization and blurred species separation have made it challenging to assign European white oak taxa based on leaf morphology and/or genetic markers and to assess the level of admixture. Therefore, there is a need for powerful markers that differentiate between taxa. Here, we established a condensed set of single-nucleotide polymorphism (SNP) markers to reliably differentiate between the three most common oak species in temperate European forests (Quercus robur, Q. petraea, Q. pubescens) and to assess the degree of admixture in a large set of selected Swiss populations.

METHODS

A training set of 194 presumably pure reference samples from Switzerland and Europe was used to assign 633 test individuals with two different approaches (population genetic-based/Bayesian vs. assumption-free/discriminative classifier) using 58 selected SNPs from coding regions. Admixture was calculated at the individual and population level with the Shannon diversity index based on individual assignment probabilities.

KEY RESULTS

Depending on the approach, 97.5-100 % of training individuals were assigned correctly, and additional analyses showed that the established SNP set could be further reduced while maintaining its discriminatory power. The two assignment approaches showed high overlap (99 %) in assigning training individuals and slightly less overlap in test individuals (84 %). Levels of admixture varied widely among populations. Mixed stands of Q. petraea and Q. pubescens revealed much higher degrees of admixture than mixed stands of the other two taxon pairs, accentuating high levels of gene flow between these two taxa in Switzerland.

CONCLUSIONS

Our set of SNPs warrants reliable taxon discrimination with great potential for further applications. We show that the three European white oak taxa have largely retained their species integrity in Switzerland despite high levels of admixture.

Demographic Inference of Divergence and Gene Exchange Between Castanopsis fabri and Castanopsis lamontii

The cytoplasmic genome of one species may be replaced by that of another species without leaving any trace of past hybridization in its nuclear genome, which can thus confuse the inference of genealogical relationship and evolutionary history of many congeneric species. In this study, we used sequence variations of chloroplast DNA and restriction site-associated DNA to investigate gene exchange between Castanopsis fabri and Castanopsis lamontii, and to infer the divergence history of the two species by comparing different divergence models based on the joint allele frequency spectrum. We evaluated climatic niche similarity of the two species using climatic variables across their entire distribution range in subtropical China. Clear genetic differentiation was revealed between C. fabri and C. lamontii, and gene exchange between the two species was discovered as a consequence of secondary contact. The gene exchange rates were variable across the genome. Gene exchange could allow C. fabri to widen its habitat through pollen swamping and broaden its climatic niche, and the chloroplast genome of C. lamontii is captured by C. fabri during this process. These results further our understanding of the timing and contribution of gene exchange to species divergence in forests.

What nature separated, and human joined together: About a spontaneous hybridization between two allopatric dogwood species (Cornus controversa and C. alternifolia)

In this study, possible hybridization between two allopatric species, Cornus controversa and Cornus alternifolia, was explored using molecular and morphological approaches. Scanning electron microscope analyses of the adaxial and the abaxial leaf surfaces yielded a few new not yet described characters typical for the particular species and intermediate for hybrids. With the use of 14 Random Amplified Polymorphic DNA and 5 Amplified Fragment Length Polymorphism primer combinations, 44 fragments species specific to C. controversa and 51 species specific to C. alternifolia were obtained. Most of these bands were also found in putative hybrids. All clustering analyses based on binary data combined from both methods confirmed a separate and intermediate status of the hybrids. Hybrid index estimates for hybrids C1-C5 indicated that all were the first generation of offspring (F₁). Chloroplast intergenic spacers (trnF-trnL and psbC-trnS) were used to infer the hybridization direction. Based on the assumption of maternal inheritance of chloroplast DNA, C. controversa seems to be the maternal parent of the hybrid. Internal transcribed spacer sequences of the five hybrids analyzed here indicated higher similarity with the sequences of C. controversa (all shared the majority of its single nucleotide polymorphisms). Sequence analysis of PI-like genes fully confirmed the hybrid origin of C1-C5 hybrids. Our results also showed that two specimens in the C. alternifolia group, A1 and A3, are not free of introgression. They are probably repeated backcrosses toward C. alternifolia. Furthermore, molecular data seem to point not only to unidirectional introgression toward C. controversa (the presence of hybrids) but to bidirectional introgression as well, since the presence of markers specific for C. controversa in the profiles of C. alternifolia specimen A3 was observed.

Abiotic genetic adaptation in the Fagaceae

Fagaceae can be found in tropical and temperate regions and contain species of major ecological and economic importance. In times of global climate change, tree populations need to adapt to rapidly changing environmental conditions. The predicted warmer and drier conditions will potentially result in locally maladapted populations. There is evidence that major genera of the Fagaceae are already negatively affected by climate change-related factors such as drought and associated biotic stressors. Therefore, knowledge of the mechanisms underlying adaptation is of great interest. In this review, we summarise current literature related to genetic adaptation to abiotic environmental conditions. We begin with an overview of genetic diversity in Fagaceae species and then summarise current knowledge related to drought stress tolerance, bud burst timing and frost tolerance in the Fagaceae. Finally, we discuss the role of hybridisation, epigenetics and phenotypic plasticity in adaptation.

Natural hybridisation among Quercus glabrescens, Q. rugosa and Q. obtusata (Fagaceae): Microsatellites and secondary metabolites markers

Natural hybridisation has significant ecological, genetic and evolutionary consequences altering morphological and chemical characters of individuals. Quercus glabrescens, Q. rugosa and Q. obtusata are white oak species well separated by their morphological characters when they occur in allopatry in Mexican temperate forests. However, in sympatry, individuals with atypical morphology have been observed, suggesting hybridisation events. In this study, we determined, with microsatellites and secondary metabolites, if interspecific gene flow occurs when these three oak species coexist in sympatry. In total, 180 individuals belonging to seven populations [three allopatric (one for each parental species) and four sympatric sites] were analysed. Allopatric populations represent well-defined genetic groups and the sympatric populations showed genetic evidence of hybridisation between Q. glabrescens × Q. rugosa and Q. glabrescens × Q. obtusata. The hybridisation percentage varied between sites and combination of involved species. We registered the presence of unique flavonoid compounds for Q. glabrescens (caffeic acid and flavonol 2), Q. rugosa (flavonol 5) and Q. obtusata (flavonol 1). Three compounds (quercetin rhamnoside, flavonol 3 and alkyl coumarate) were expressed in all taxa. Finally, the hybrid genotypes identified in this study (Q. glabrescens × Q. rugosa and Q. glabrescens × Q. obtusata) showed specific chemical profiles, resulting from a combination of those of their parental species. These results show that hybridisation events between these oak species alter chemical expression of secondary metabolites, creating a mosaic of resources and conditions that provide the substrate for different combinations of foliar-associated species such as herbivores, endophytic fungi or epiphyte plants.

Diversification, adaptation, and community assembly of the American oaks (Quercus), a model clade for integrating ecology and evolution

Contents Summary 669 I. Model clades for the study and integration of ecology and evolution 670 II. Oaks: an important model clade 671 III. Insights from the history of the American oaks for understanding community assembly and ecosystem dominance 673 IV. Bridging the gap between micro- and macroevolutionary processes relevant to ecology 679 V. How do we reconcile evidence for adaptive evolution with niche conservatism and long-term stasis? 682 VI. High plasticity and within-population genetic variation contribute to population persistence 683 VII. Emerging technologies for tracking functional change 685 VIII. Conclusions 685 Acknowledgements 686 References 686 SUMMARY: Ecologists and evolutionary biologists are concerned with explaining the diversity and composition of the natural world and are aware of the inextricable linkages between ecological and evolutionary processes that maintain the Earth's life support systems. Yet examination of these linkages remains challenging due to the contrasting nature of focal systems and research approaches. Model clades provide a critical means to integrate ecology and evolution, as illustrated by the oaks (genus Quercus), an important model clade, given their ecological dominance, remarkable diversity, and growing phylogenetic, genomic, and ecological data resources. Studies of the clade reveal that their history of sympatric parallel adaptive radiation continues to influence community assembly today, highlighting questions on the nature and extent of coexistence mechanisms. Flexible phenology and hydraulic traits, despite evolutionary stasis, may have enabled adaptation to a wide range of environments within and across species, contributing to their high abundance and diversity. The oaks offer fundamental insights at the intersection of ecology and evolution on the role of diversification in community assembly processes, on the importance of flexibility in key functional traits in adapting to new environments, on factors contributing to persistence of long-lived organisms, and on evolutionary legacies that influence ecosystem function.

Similarity analysis between species of the genus Quercus L. (Fagaceae) in southern Italy based on the fractal dimension

The fractal dimension (FD) is calculated for seven species of the genus Quercus L. in Calabria region (southern Italy), five of which have a marcescent-deciduous and two a sclerophyllous character. The fractal analysis applied to the leaves reveals different FD values for the two groups. The difference between the means and medians is very small in the case of the marcescent-deciduous group and very large when these differences are established between both groups: all this highlights the distance between the two groups in terms of similarity. Specifically, Q.crenata, which is hybridogenic in origin and whose parental species are Q.cerris and Q.suber, is more closely related to Q.cerris than to Q.suber, as also expressed in the molecular analysis. We consider that, in combination with other morphological, physiological and genetic parameters, the fractal dimension is a useful tool for studying similarities amongst species.

A nuclear DNA barcode for eastern North American oaks and application to a study of hybridization in an Arboretum setting

DNA barcoding has proved difficult in a number of woody plant genera, including the ecologically important oak genus Quercus. In this study, we utilized restrictionsite-associated DNA sequencing (RAD-seq) to develop an economical single nucleotide polymorphism (SNP) DNA barcoding system that suffices to distinguish eight common, sympatric eastern North American white oak species. Two de novo clustering pipelines, PyRAD and Stacks, were used in combination with postclustering bioinformatic tools to generate a list of 291 potential SNPs, 80 of which were included in a barcoding toolkit that is easily implemented using MassARRAY mass spectrometry technology. As a proof-of-concept, we used the genotyping toolkit to infer potential hybridization between North American white oaks transplanted outside of their native range (Q. michauxii, Q. montana, Q muehlenbergii/Q. prinoides, and Q. stellata) into a horticultural collection surrounded by natural forests of locally native trees (Q. alba and Q. macrocarpa) in the living collection at The Morton Arboretum (Lisle, IL, USA). Phylogenetic and clustering analyses suggested low rates of hybridization between cultivated and native species, with the exception of one Q. michauxii mother tree, the acorns of which exhibited high admixture from either Q. alba or Q. stellata and Q. macrocarpa, and a hybrid between Q. stellata that appears to have backcrossed almost exclusively to Q. alba. Together, RAD-seq and MassARRAY technologies allow for efficient development and implementation of a multispecies barcode for one of the more challenging forest tree genera.

The Application and Limitation of Universal Chloroplast Markers in Discriminating East Asian Evergreen Oaks

The East Asian subtropics mostly occupied by evergreen broad-leaved forests (EBLFs), is one of the global diversity centers for evergreen oaks. Evergreen oaks are keystone canopy trees in EBLFs with important ecosystem function and crucial significance for regional biodiversity conservation. However, the species composition and diversity of Asian evergreen oaks are poorly understood. Here, we test whether the four chloroplast markers atpI-atpH, matK, psbA-trnH, and ycf1, can discriminate the two evergreen oak sections in Asia - Cyclobalanopsis and Ilex. Two hundred and seventy-two individuals representing 57 species were scanned and 17 species from other oaks sections were included for phylogenetic reconstruction. The genetic diversity of the Quercus sections was also compared. Overall, we found that universal chloroplast DNA (cpDNA) barcoding markers could resolve two clades in Quercus, i.e., subgenus Cerris (Old World Clade) and subgenus Quercus (New World Clade). The chloroplast markers distinguished the main sections, with few exceptions. Each cpDNA region showed no barcoding gap and none of them provided good resolution at the species level. The best species resolution (27.78%) was obtained when three or four markers were combined and analyzed using BLAST. The high conservation of the cpDNA and complicated evolutionary patterns, due to incomplete lineage sorting, interspecific hybridization and introgressions may hinder the ability of cpDNA markers to discriminate different species. When comparing diversification pattern across Quercus sections (Cyclobalanopsis, Ilex, Cerris, Quercus, and Protobalanus), we found that section Ilex was the most genetically diverse, and section Cyclobalanopsis was lower genetically diverse. This diversification pattern may have resulted from the interplay of the Eurasia Cenozoic tectonic movements, climate changes and different niches of their ancestral lineages.

Species Boundaries Between Three Sympatric Oak Species: Quercus aliena, Q. dentata, and Q. variabilis at the Northern Edge of Their Distribution in China

Oaks are important timber trees with wide distributions in China, but few genetic studies have been conducted on a fine scale. In this study, we seek to investigate the genetic diversity and differentiation of three sympatric oak species (Quercus aliena Blume, Quercus dentata Thunb. ex Murray, and Quercus variabilis Blume) in their northern distribution in China using 17 bi-parentally inherited nSSRs markers and five maternally inherited chloroplast DNA (cpDNA) fragments. Both the cpDNA and the nSSRs show a high level of genetic differentiation between different oak sections. The chloroplast haplotypes are clustered into two lineages. Clear species boundaries are detected between Q. variabilis and either Q. aliena or Q. dentata. The sharing of chloroplast haplotype H1 between Q. aliena and Q. dentata suggests very recent speciation and incomplete lineage sorting or introgression of H1 from one species to another. The nSSRs data indicate a complete fixation of variation within sites for all three oak species, and that extensive gene flow occurs within species whereas only limited gene flow is detected between Q. aliena and Q. dentata and nearly no gene flow can be detected between Q. aliena and Q. variabilis and between Q. dentata and Q. variabilis. Prezygotic isolation may have contributed to the species boundaries of these three sympatric oak species.

A genetic legacy of introgression confounds phylogeny and biogeography in oaks

Oaks (Quercus L.) have long been suspected to hybridize in nature, and widespread genetic exchange between morphologically defined species is well documented in two- to six-species systems, but the phylogenetic consequences of hybridization in oaks have never been demonstrated in a phylogenetically diverse sample. Here, we present phylogenomic analyses of a ca 30 Myr clade that strongly support morphologically defined species and the resolution of novel clades of white oaks; however, historical hybridization across clade boundaries is detectable and, undiagnosed, would obscure the imprint of biogeographic history in the phylogeny. Phylogenetic estimation from restriction-site-associated DNA sequencing data for 156 individuals representing 81 species supports two successive intercontinental disjunctions of white oaks: an early vicariance between the Eurasian and American white oaks, and a second, independent radiation represented by two relictual species. A suite of subsampled and partitioned analyses, however, supports a more recent divergence of the Eurasian white oaks from within the American white oaks and suggests that historic introgression between the Eurasian white oaks and a now-relictual lineage biases concatenated phylogenetic estimates. We demonstrate how divergence and reticulation both influence our understanding of the timing and nature of diversification and global colonization in these ecologically and economically important taxa.

Exploring the potential of gametic reconstruction of parental genotypes by F1 hybrids as a bridge for rapid introgression

Interspecific hybridization and genetic introgression are commonly observed in natural populations of many species, especially trees. Among oaks, gene flow between closely related species has been well documented. And yet, hybridization does not lead to a "melting pot", i.e., the homogenization of phenotypic traits. Here, we explore how the combination of several common reproductive and genomic traits could create an avenue for interspecific gene flow that partially explains this apparent paradox. During meiosis, F₁ hybrids will produce approximately (½)ⁿ "reconstructed" parental gametes, where n equals the number of chromosomes. Crossing over would introduce a small amount of introgressive material. The resulting parental-type gametophytes would probably possess a similar fertilization advantage as conspecific pollen. The resulting "backcross" would actually be the genetic equivalent of a conspecific out-cross, with a small amount of heterospecific DNA captured through crossing over. Even with detailed genomic analysis, the resulting offspring would not appear to be a backcross. This avenue for rapid introgression between species through the F₁ hybrid will be viable for organisms that meet certain conditions: low base chromosome number, conserved genomic structure and size, production of billions of gametes/gametophytes during each reproductive event, and conspecific fertilization advantage.

Extensive recent secondary contacts between four European white oak species

Historical trajectories of tree species during the late Quaternary have been well reconstructed through genetic and palaeobotanical studies. However, many congeneric tree species are interfertile, and the timing and contribution of introgression to species divergence during their evolutionary history remains largely unknown. We quantified past and current gene flow events between four morphologically divergent oak species (Quercus petraea, Q. robur, Q. pyrenaica, Q. pubescens), by two independent inference methods: diffusion approximation to the joint frequency spectrum (∂a∂i) and approximate Bayesian computation (ABC). For each pair of species, alternative scenarios of speciation allowing gene flow over different timescales were evaluated. Analyses of 3524 single nucleotide polymorphisms (SNPs) randomly distributed in the genome, showed that these species evolved in complete isolation for most of their history, but recently came into secondary contact, probably facilitated by the most recent period of postglacial warming. We demonstrated that: there was sufficient genetic differentiation before secondary contact for the accumulation of barriers to gene flow; and current European white oak genomes are a mosaic of genes that have crossed species boundaries and genes impermeable to gene flow.

Introgression Threatens the Genetic Diversity of Quercus austrocochinchinensis (Fagaceae), an Endangered Oak: A Case Inferred by Molecular Markers

Natural introgression can cause negative effects where rare species experience genetic assimilation and invade by their abundant congeners. Quercus austrocochinchinensis and Q. kerrii (subgenus Cyclobalanopsis) are a pair of closely related species in the Indo-China area. Morphological intermediates of the two species have been reported in this region. In this study, we used AFLP, SSR and two key leaf morphological diagnostic traits to study the two Q. austrocochinchinensis populations, two pure Q. kerrii and two putative hybrid populations in China. Rates of individual admixture were examined using the Bayesian clustering programs STRUCTURE and NewHybrids, with no a priori species assignment. In total, we obtained 151 SSR alleles and 781 polymorphic loci of AFLP markers. Population differentiation inferred by SSR and AFLP was incoherent with recognized species boundaries. Bayesian admixture analyses and principal coordinate analysis identified more hybrids and backcrossed individuals than morphological intermediates in the populations. SSR inferred a wide genetic assimilation in Q. austrocochinchinensis, except for subpopulation D2 in the core area of Xi-Shuang-Ban-Na Nature Reserve (XSBN). However, AFLP recognized more Q. austrocochinchinensis purebreds than SSR. Analysis using NewHybrids on AFLP data indicated that these hybridized individuals were few F2 and predominantly backcrosses with both parental species. All these evidences indicate the formation of a hybrid swarm at XSBN where the two species co-exist. Both AFLP and SSR recognized that the core protected area of XSBN (D2) has a high percentage of Q. austrocochinchinensis purebreds and a unique germplasm. The Hainan population and the other subpopulations of XSBN of the species might have lost their genetic integrity. Our results revealed a clear genetic differentiation in the populations and subpopulations of Q. austrocochinchinensis and ongoing introgression between Q. austrocochinchinensis and Q. kerrii at the disturbed contact areas. Combining the results from genetic and morphological analyses, the conservation of subpopulation D2 should be prioritized. Conservation and restoration of the integrity of tropical ravine rainforest is an important long-term goal for the successful conservation of Q. austrocochinchinensis. The fine-scale landscape might play an essential role in shaping the spatial patterns of hybridization. Further studies are needed to evaluate these patterns and dynamics.

Impacts of human-induced environmental disturbances on hybridization between two ecologically differentiated Californian oak species

Natural hybridization, which can be involved in local adaptation and in speciation processes, has been linked to different sources of anthropogenic disturbance. Here, we use genotypic data to study range-wide patterns of genetic admixture between the serpentine-soil specialist leather oak (Quercus durata) and the widespread Californian scrub oak (Quercus berberidifolia). First, we estimated hybridization rates and the direction of gene flow. Second, we tested the hypothesis that genetic admixture increases with different sources of environmental disturbance, namely anthropogenic destruction of natural habitats and wildfire frequency estimated from long-term records of fire occurrence. Our analyses indicate considerable rates of hybridization (> 25%), asymmetric gene flow from Q. durata into Q. berberidifolia, and a higher occurrence of hybrids in areas where both species live in close parapatry. In accordance with the environmental disturbance hypothesis, we found that genetic admixture increases with wildfire frequency, but we did not find a significant effect of other sources of human-induced habitat alteration (urbanization, land clearing for agriculture) or a suite of ecological factors (climate, elevation, soil type). Our findings highlight that wildfires constitute an important source of environmental disturbance, promoting hybridization between two ecologically well-differentiated native species.

Comparative Pollen Morphological Analysis and Its Systematic Implications on Three European Oak (Quercus L., Fagaceae) Species and Their Spontaneous Hybrids

Pollen morphology of three parental Quercus species (Q. robur L., Q. petraea (Matt) Liebl, Q. pubescens Willd.) and two spontaneous hybrids of these species (Q. ×calvescens Vuk. = Q. petraea × Q. pubescens and Q. ×rosacea Bechst. = Q. robur × Q. petraea) was investigated in this study. The pollen originated from 18 natural oak sites and 67 individuals (oak trees). Each individual was represented by 30 pollen grains. In total, 2010 pollen grains were measured. They were analysed for nine quantitative and four qualitative features. Pollen size and shape were important features to diagnosing Quercus parental species and hybrids. On the basis of exine ornamentation, it was possible to identify only Q. pubescens, while the remaining species and hybrids did not differ significantly with respect to this feature. The determination of the diagnostic value of endoaperture features requires further palynological studies. On the basis of pollen size and shape Q. robur × Q. petraea was clearly separated. Grouping of 67 oak trees on the basis of pollen grain features has shown that individuals from different as well as same taxa occurred in the same groups. Likewise, with respect to natural sites, oak trees originating from the same places as well as from geographically distant ones, grouped together. Pollen morphological features allow to distinguish a part of the studied Quercus taxa. Therefore, it can be used as an auxiliary feature in the taxonomy.

Comparative Analysis of the Complete Chloroplast Genomes of Five Quercus Species

Quercus is considered economically and ecologically one of the most important genera in the Northern Hemisphere. Oaks are taxonomically perplexing because of shared interspecific morphological traits and intraspecific morphological variation, which are mainly attributed to hybridization. Universal plastid markers cannot provide a sufficient number of variable sites to explore the phylogeny of this genus, and chloroplast genome-scale data have proven to be useful in resolving intractable phylogenetic relationships. In this study, the complete chloroplast genomes of four Quercus species were sequenced, and one published chloroplast genome of Quercus baronii was retrieved for comparative analyses. The five chloroplast genomes ranged from 161,072 bp (Q. baronii) to 161,237 bp (Q. dolicholepis) in length, and their gene organization and order, and GC content, were similar to those of other Fagaceae species. We analyzed nucleotide substitutions, indels, and repeats in the chloroplast genomes, and found 19 relatively highly variable regions that will potentially provide plastid markers for further taxonomic and phylogenetic studies within Quercus. We observed that four genes (ndhA, ndhK, petA, and ycf1) were subject to positive selection. The phylogenetic relationships of the Quercus species inferred from the chloroplast genomes obtained moderate-to-high support, indicating that chloroplast genome data may be useful in resolving relationships in this genus.

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.

High-density linkage mapping and distribution of segregation distortion regions in the oak genome

We developed the densest single-nucleotide polymorphism (SNP)-based linkage genetic map to date for the genus Quercus An 8k gene-based SNP array was used to genotype more than 1,000 full-sibs from two intraspecific and two interspecific full-sib families of Quercus petraea and Quercus robur A high degree of collinearity was observed between the eight parental maps of the two species. A composite map was then established with 4,261 SNP markers spanning 742 cM over the 12 linkage groups (LGs) of the oak genome. Nine genomic regions from six LGs displayed highly significant distortions of segregation. Two main hypotheses concerning the mechanisms underlying segregation distortion are discussed: genetic load vs. reproductive barriers. Our findings suggest a predominance of pre-zygotic to post-zygotic barriers.

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.

Evolutionary and demographic history of the Californian scrub white oak species complex: an integrative approach

Understanding the factors promoting species formation is a major task in evolutionary research. Here, we employ an integrative approach to study the evolutionary history of the Californian scrub white oak species complex (genus Quercus). To infer the relative importance of geographical isolation and ecological divergence in driving the speciation process, we (i) analysed inter- and intraspecific patterns of genetic differentiation and employed an approximate Bayesian computation (ABC) framework to evaluate different plausible scenarios of species divergence. In a second step, we (ii) linked the inferred divergence pathways with current and past species distribution models (SDMs) and (iii) tested for niche differentiation and phylogenetic niche conservatism across taxa. ABC analyses showed that the most plausible scenario is the one considering the divergence of two main lineages followed by a more recent pulse of speciation. Genotypic data in conjunction with SDMs and niche differentiation analyses support that different factors (geography vs. environment) and modes of speciation (parapatry, allopatry and maybe sympatry) have played a role in the divergence process within this complex. We found no significant relationship between genetic differentiation and niche overlap, which probably reflects niche lability and/or that multiple factors, have contributed to speciation. Our study shows that different mechanisms can drive divergence even among closely related taxa representing early stages of species formation and exemplifies the importance of adopting integrative approaches to get a better understanding of the speciation process.

A linkage disequilibrium perspective on the genetic mosaic of speciation in two hybridizing Mediterranean white oaks

We analyzed the genetic mosaic of speciation in two hybridizing Mediterranean white oaks from the Iberian Peninsula (Quercus faginea Lamb. and Quercus pyrenaica Willd.). The two species show ecological divergence in flowering phenology, leaf morphology and composition, and in their basic or acidic soil preferences. Ninety expressed sequence tag-simple sequence repeats (EST-SSRs) and eight nuclear SSRs were genotyped in 96 trees from each species. Genotyping was designed in two steps. First, we used 69 markers evenly distributed over the 12 linkage groups (LGs) of the oak linkage map to confirm the species genetic identity of the sampled genotypes, and searched for differentiation outliers. Then, we genotyped 29 additional markers from the chromosome bins containing the outliers and repeated the multilocus scans. We found one or two additional outliers within four saturated bins, thus confirming that outliers are organized into clusters. Linkage disequilibrium (LD) was extensive; even for loosely linked and for independent markers. Consequently, score tests for association between two-marker haplotypes and the 'species trait' showed a broad genomic divergence, although substantial variation across the genome and within LGs was also observed. We discuss the influence of several confounding effects on neutrality tests and review the evolutionary processes leading to extensive LD. Finally, we examine how LD analyses within regions that contain outlier clusters and quantitative trait loci can help to identify regions of divergence and/or genomic hitchhiking in the light of predictions from ecological speciation theory.

Infrared spectroscopy of pollen identifies plant species and genus as well as environmental conditions

BACKGROUND

It is imperative to have reliable and timely methodologies for analysis and monitoring of seed plants in order to determine climate-related plant processes. Moreover, impact of environment on plant fitness is predominantly based on studies of female functions, while the contribution of male gametophytes is mostly ignored due to missing data on pollen quality. We explored the use of infrared spectroscopy of pollen for an inexpensive and rapid characterization of plants.

METHODOLOGY

The study was based on measurement of pollen samples by two Fourier transform infrared techniques: single reflectance attenuated total reflectance and transmission measurement of sample pellets. The experimental set, with a total of 813 samples, included five pollination seasons and 300 different plant species belonging to all principal spermatophyte clades (conifers, monocotyledons, eudicots, and magnoliids).

RESULTS

The spectroscopic-based methodology enables detection of phylogenetic variations, including the separation of confamiliar and congeneric species. Furthermore, the methodology enables measurement of phenotypic plasticity by the detection of inter-annual variations within the populations. The spectral differences related to environment and taxonomy are interpreted biochemically, specifically variations of pollen lipids, proteins, carbohydrates, and sporopollenins. The study shows large variations of absolute content of nutrients for congenital species pollinating in the same environmental conditions. Moreover, clear correlation between carbohydrate-to-protein ratio and pollination strategy has been detected. Infrared spectral database with respect to biochemical variation among the range of species, climate and biogeography will significantly improve comprehension of plant-environment interactions, including impact of global climate change on plant communities.

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.

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.

Nuclear and plastid markers reveal the persistence of genetic identity: a new perspective on the evolutionary history of Petunia exserta

Recently divergent species that can hybridize are ideal models for investigating the genetic exchanges that can occur while preserving the species boundaries. Petunia exserta is an endemic species from a very limited and specific area that grows exclusively in rocky shelters. These shaded spots are an inhospitable habitat for all other Petunia species, including the closely related and widely distributed species P. axillaris. Individuals with intermediate morphologic characteristics have been found near the rocky shelters and were believed to be putative hybrids between P. exserta and P. axillaris, suggesting a situation where Petunia exserta is losing its genetic identity. In the current study, we analyzed the plastid intergenic spacers trnS/trnG and trnH/psbA and six nuclear CAPS markers in a large sampling design of both species to understand the evolutionary process occurring in this biological system. Bayesian clustering methods, cpDNA haplotype networks, genetic diversity statistics, and coalescence-based analyses support a scenario where hybridization occurs while two genetic clusters corresponding to two species are maintained. Our results reinforce the importance of coupling differentially inherited markers with an extensive geographic sample to assess the evolutionary dynamics of recently diverged species that can hybridize.

Seeing the forest through the trees: comprehensive inference on individual mating patterns in a mixed stand of Quercus robur and Q. petraea

BACKGROUND AND AIMS

Sexual reproduction is one of the most important moments in a life cycle, determining the genetic composition of individual offspring. Controlled pollination experiments often show high variation in the mating system at the individual level, suggesting a persistence of individual variation in natural populations. Individual variation in mating patterns may have significant adaptive implications for a population and for the entire species. Nevertheless, field data rarely address individual differences in mating patterns, focusing rather on averages. This study aimed to quantify individual variation in the different components of mating patterns.

METHODS

Microsatellite data were used from 421 adult trees and 1911 seeds, structured in 72 half-sib families collected in a single mixed stand of Quercus robur and Q. petraea in northern Poland. Using a Bayesian approach, mating patterns were investigated, taking into account pollen dispersal, male fecundity, possible hybridization and heterogeneity in immigrant pollen pools.

KEY RESULTS

Pollen dispersal followed a heavy-tailed distribution (283 m on average). In spite of high pollen mobility, immigrant pollen pools showed strong genetic structuring among mothers. At the individual level, immigrant pollen pools showed highly variable divergence rates, revealing that sources of immigrant pollen can vary greatly among particular trees. Within the stand, the distribution of male fecundity appeared highly skewed, with a small number of dominant males, resulting in a ratio of census to effective density of pollen donors of 5·3. Male fecundity was not correlated with tree diameter but showed strong cline-like spatial variation. This pattern can be attributed to environmental variation. Quercus petraea revealed a greater preference (74 %) towards intraspecific mating than Q. robur (36 %), although mating preferences varied among trees.

CONCLUSIONS

Mating patterns can reveal great variation among individuals, even within a single even-age stand. The results show that trees can mate assortatively, with little respect for spatial proximity. Such selective mating may be a result of variable combining compatibility among trees due to genetic and/or environmental factors.

Inter- and intra-specific variability in isoprene production and photosynthesis of Central European oak species

European deciduous oaks are closely related and are known for their strong emission of volatile isoprenoids. They are chemo-taxonomically diverse, but hybridise frequently. Four-year-old oak seedlings growing together in a model ecosystem facility under near-natural conditions were studied. The leaves were morphologically classified in the three oak species Quercus robur, Q. pubescens and Q. petraea (with four provenances each) and further investigated by a molecular-genetic approach. Q. robur was morphologically and genetically clearly different from Q. pubescens and Q. petraea, whereas Q. pubescens and Q. petraea individuals used in this study were morphologically and genetically more similar. There was a minor impact of among and within species variability on isoprene synthesis, isoprene emission and photosynthesis. Isoprene emission rates normalised to 25 °C leaf temperature ranged from 5.78 to 10.66 nmol m(-2)  s(-1) , whereas photosynthesis ranged from 12.8 to 17.6 μmol m(-2)  s(-1) . On cloudy days, among the provenances of each species, only net photosynthesis of the Q. robur provenance Hünenberg was reduced and isoprene synthase activity of the Q. pubescens provenance Promotogno increased. On sunny days, photosynthesis did not differ among the provenances. Over all provenances, gas exchange on cloudy days did not differ significantly from sunny days. In the combined data of cloudy and sunny days, no differences between the studied provenances and oak species were detected in isoprene emission and photosynthesis. Thus, isoprene emission and photosynthesis rates were remarkably stable among oak species and provenances. The results indicate that taxonomic differences in the studied oak species are not reflected in isoprene emission and photosynthesis, probably because of the high plasticity of gene expression resulting in high phenotypic flexibility.

Novel SNP development and analysis at a NADP+ -specific IDH enzyme gene in a four species mixed oak forest

Closely related Quercus species generally exhibit low levels of genetic differentiation despite their ecological and morphological differences. However, at a few so-called 'outlier' loci they seem to remain genetically distinct. Isocitrate dehydrogenases (IDH) are key enzymes involved in the metabolic pathway of the citrate cycle. IDH has also been characterised as an 'outlier' marker, significantly differentiating the closely related Q. robur and Q. petraea with the isozyme technique. This ability to differentiate the species was tested here at molecular level: 13 single nucleotide polymorphism (SNP) markers were identified and developed within a NADP(+) -specific IDH gene in Quercus spp. and applied as molecular markers in a four species mixed oak forest in eastern Europe, where Q. robur, Q. petraea, Q. pubescens and Q. frainetto naturally co-exist. From the 13 developed SNPs, three groups were formed: non-synonymous, synonymous and non-coding SNPs. The levels of total gene diversity were moderate for all species investigated. The non-synonymous SNPs showed lower levels of gene diversity. Overall, the four closely related Quercus spp. were significantly differentiated (except Q. petraea with Q. frainetto). Analysis of non-random association of alleles revealed no clear physical clustering of the SNP sites in significant linkage disequilibrium (LD). However, separate LD analysis for each species showed a lower number of sites in significant LD for Q. robur than for the other species, possibly reflecting the history of the species in this specific geographical site and less efficient recombination effect due to the larger effective population size of Q. robur. Eleven statistically significant associations were found between seven SNPs and morphological traits that are commonly used to differentiate oak species.