Uncovering the genomic signature of ancient introgression between white oak lineages (Quercus)

Introgression and incomplete lineage sorting blurred phylogenetic relationships across the genomes of sclerophyllous oaks from southwest China

Resolving evolutionary relationships among closely related species with interspecific gene flow is challenging. Genome-scale data provide opportunities to clarify complex evolutionary relationships in closely related species and to observe variations in species relationships across the genomes of such species. The Himalayan-Hengduan subalpine oaks have a nearly completely sympatric distribution in southwest China and probably constitute a syngameon. In this study, we mapped resequencing data from different species in this group to the Quercus aquifolioides reference genome to obtain a high-quality filtered single nucleotide polymorphism (SNP) dataset. We also assembled their plastomes. We reconstructed their phylogenetic relationships, explored the level and pattern of introgression among these species and investigated gene tree variation in the genomes of these species using sliding windows. The same or closely related plastomes were found to be shared extensively among different species within a specific geographical area. Phylogenomic analyses of genome-wide SNP data found that most oaks in the Himalayan-Hengduan subalpine clade showed genetic coherence, but several species were found to be connected by introgression. The gene trees obtained using sliding windows showed that the phylogenetic relationships in the genomes of oaks are highly heterogeneous and therefore highly obscured. Our study found that all the oaks of the Himalayan-Hengduan subalpine clade from southwest China form a syngameon. The obscured phylogenetic relationships observed empirically across the genome are best explained by interspecific gene flow in conjunction with incomplete lineage sorting.

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

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

Divergence and reticulation in the Mexican white oaks: ecological and phylogenomic evidence on species limits and phylogenetic networks in the Quercus laeta complex (Fagaceae)

BACKGROUND AND AIMS

Introgressive hybridization poses a challenge to taxonomic and phylogenetic understanding of taxa, particularly when there are high numbers of co-occurring, intercrossable species. The genus Quercus exemplifies this situation. Oaks are highly diverse in sympatry and cross freely, creating syngameons of interfertile species. Although a well-resolved, dated phylogeny is available for the American oak clade, evolutionary relationships within many of the more recently derived clades remain to be defined, particularly for the young and exceptionally diverse Mexican white oak clade. Here, we adopted an approach bridging micro- and macroevolutionary scales to resolve evolutionary relationships in a rapidly diversifying clade endemic to Mexico.

METHODS

Ecological data and sequences of 155 low-copy nuclear genes were used to identify distinct lineages within the Quercus laeta complex. Concatenated and coalescent approaches were used to assess the phylogenetic placement of these lineages relative to the Mexican white oak clade. Phylogenetic network methods were applied to evaluate the timing and genomic significance of recent or historical introgression among lineages.

KEY RESULTS

The Q. laeta complex comprises six well-supported lineages, each restricted geographically and with mostly divergent climatic niches. Species trees corroborated that the different lineages are more closely related to other species of Mexican white oaks than to each other, suggesting that this complex is polyphyletic. Phylogenetic networks estimated events of ancient introgression that involved the ancestors of three present-day Q. laeta lineages.

CONCLUSIONS

The Q. laeta complex is a morphologically and ecologically related group of species rather than a clade. Currently, oak phylogenetics is at a turning point, at which it is necessary to integrate phylogenetics and ecology in broad regional samples to figure out species boundaries. Our study illuminates one of the more complicated of the Mexican white oak groups and lays groundwork for further taxonomic study.

Chromosome-scale genomes of Quercus sichourensis and Quercus rex provide insights into the evolution and adaptation of Fagaceae

The Fagaceae, a plant family with a wide distribution and diverse adaptability, has garnered significant interest as a subject of study in plant speciation and adaptation. Meanwhile, certain Fagaceae species are regarded as highly valuable wood resources due to the exceptional quality of their wood. In this study, we present two high-quality, chromosome-scale genome sequences for Quercus sichourensis (848.75 Mb) and Quercus rex (883.46 Mb). Comparative genomics analysis reveals that the difference in the number of plant disease resistance genes and the nonsynonymous and synonymous substitution ratio (Ka/Ks) of protein-coding genes among Fagaceae species are related to different environmental adaptations. Interestingly, most genes related to starch synthesis in the investigated Quercoideae species are located on a single chromosome, as compared to the outgroup species, Fagus sylvatica. Furthermore, resequencing and population analysis of Q. sichourensis and Q. rex reveal that Q. sichourensis has lower genetic diversity and higher deleterious mutations compared to Q. rex. The high-quality, chromosome-level genomes and the population genomic analysis of the critically endangered Q. sichourensis and Q. rex will provide an invaluable resource as well as insights for future study in these two species, even the genus Quercus, to facilitate their conservation.

Ensemble species distribution modeling and multilocus phylogeography provide insight into the spatial genetic patterns and distribution dynamics of a keystone forest species, Quercus glauca

BACKGROUND

Forests are essential for maintaining species diversity, stabilizing local and global climate, and providing ecosystem services. Exploring the impact of paleogeographic events and climate change on the genetic structure and distribution dynamics of forest keystone species could help predict responses to future climate change. In this study, we combined an ensemble species distribution model (eSDM) and multilocus phylogeography to investigate the spatial genetic patterns and distribution change of Quercus glauca Thunb, a keystone of East Asian subtropical evergreen broad-leaved forest.

RESULTS

A total of 781 samples were collected from 77 populations, largely covering the natural distribution of Q. glauca. The eSDM showed that the suitable habitat experienced a significant expansion after the last glacial maximum (LGM) but will recede in the future under a general climate warming scenario. The distribution centroid will migrate toward the northeast as the climate warms. Using nuclear SSR data, two distinct lineages split between east and west were detected. Within-group genetic differentiation was higher in the West than in the East. Based on the identified 58 haplotypes, no clear phylogeographic structure was found. Populations in the Nanling Mountains, Wuyi Mountains, and the southwest region were found to have high genetic diversity.

CONCLUSIONS

A significant negative correlation between habitat stability and heterozygosity might be explained by the mixing of different lineages in the expansion region after LGM and/or hybridization between Q. glauca and closely related species. The Nanling Mountains may be important for organisms as a dispersal corridor in the west-east direction and as a refugium during the glacial period. This study provided new insights into spatial genetic patterns and distribution dynamics of Q. glauca.

Nightmare or delight: Taxonomic circumscription meets reticulate evolution in the phylogenomic era

Phylogenetic studies in the phylogenomics era have demonstrated that reticulate evolution greatly impedes the accuracy of phylogenetic inference, and consequently can obscure taxonomic treatments. However, the systematics community lacks a broadly applicable strategy for taxonomic delimitation in groups characterized by pervasive reticulate evolution. The red-fruit genus, Stranvaesia, provides an ideal model to examine the influence of reticulation on generic circumscription, particularly where hybridization and allopolyploidy dominate the evolutionary history. In this study, we conducted phylogenomic analyses integrating data from hundreds of single-copy nuclear (SCN) genes and plastomes, and interrogated nuclear paralogs to clarify the inter/intra-generic relationship of Stranvaesia and its allies in the framework of Maleae. Analyses of phylogenomic discord and phylogenetic networks showed that allopolyploidization and introgression promoted the origin and diversification of the Stranvaesia clade, a conclusion further bolstered by cytonuclear and gene tree discordance. With a well-inferred phylogenetic backbone, we propose an updated generic delimitation of Stranvaesia and introduce a new genus, Weniomeles. This new genus is distinguished by its purple-black fruits, thorns trunk and/or branches, and a distinctive fruit core anatomy characterized by multilocular separated by a layer of sclereids and a cluster of sclereids at the top of the locules. Through this study, we highlight a broadly-applicable workflow that underscores the significance of reticulate evolution analyses in shaping taxonomic revisions from phylogenomic data.

Phylogenomics reveals extreme gene tree discordance in a lineage of dominant trees: hybridization, introgression, and incomplete lineage sorting blur deep evolutionary relationships despite clear species groupings in Eucalyptus subgenus Eudesmia

Eucalypts are a large and ecologically important group of plants on the Australian continent, and understanding their evolution is important in understanding evolution of the unique Australian flora. Previous phylogenies using plastome DNA, nuclear-ribosomal DNA, or random genome-wide SNPs, have been confounded by limited genetic sampling or by idiosyncratic biological features of the eucalypts, including widespread plastome introgression. Here we present phylogenetic analyses of Eucalyptus subgenus Eudesmia (22 species from western, northern, central and eastern Australia), in the first study to apply a target-capture sequencing approach using custom, eucalypt-specific baits (of 568 genes) to a lineage of Eucalyptus. Multiple accessions of all species were included, and target-capture data were supplemented by separate analyses of plastome genes (average of 63 genes per sample). Analyses revealed a complex evolutionary history likely shaped by incomplete lineage sorting and hybridization. Gene tree discordance generally increased with phylogenetic depth. Species, or groups of species, toward the tips of the tree are mostly supported, and three major clades are identified, but the branching order of these clades cannot be confirmed with confidence. Multiple approaches to filtering the nuclear dataset, by removing genes or samples, could not reduce gene tree conflict or resolve these relationships. Despite inherent complexities in eucalypt evolution, the custom bait kit devised for this research will be a powerful tool for investigating the evolutionary history of eucalypts more broadly.

Introgression between highly divergent fungal sister species

To understand how species evolve and adapt to changing environments, it is important to study gene flow and introgression due to their influence on speciation and radiation events. Here, we apply a novel experimental system for investigating these mechanisms using natural populations. The system is based on two fungal sister species with morphological and ecological similarities occurring in overlapping habitats. We examined introgression between these species by conducting whole genome sequencing of individuals from populations in North America and Europe. We assessed genome-wide nucleotide divergence and performed crossing experiments to study reproductive barriers. We further used ABBA-BABA statistics together with a network analysis to investigate introgression, and conducted demographic modelling to gain insight into divergence times and introgression events. The results revealed that the species are highly divergent and incompatible in vitro. Despite this, small regions of introgression were scattered throughout the genomes and one introgression event likely involves a ghost population (extant or extinct). This study demonstrates that introgression can be found among divergent species and that population histories can be studied without collections of all the populations involved. Moreover, the experimental system is shown to be a useful tool for research on reproductive isolation in natural populations.

Niche evolution in a northern temperate tree lineage: biogeographical legacies in cork oaks (Quercus section Cerris)

BACKGROUND AND AIMS

Cork oaks (Quercus section Cerris) comprise 15 extant species in Eurasia. Despite being a small clade, they display a range of leaf morphologies comparable to the largest sections (>100 spp.) in Quercus. Their fossil record extends back to the Eocene. Here, we explore how cork oaks achieved their modern ranges and how legacy effects might explain niche evolution in modern species of section Cerris and its sister section Ilex, the holly oaks.

METHODS

We inferred a dated phylogeny for cork and holly oaks using a reduced-representation next-generation sequencing method, restriction site-associated DNA sequencing (RAD-seq), and used D-statistics to investigate gene flow hypotheses. We estimated divergence times using a fossilized birth-death model calibrated with 47 fossils. We used Köppen profiles, selected bioclimatic parameters and forest biomes occupied by modern species to infer ancestral climatic and biotic niches.

KEY RESULTS

East Asian and Western Eurasian cork oaks diverged initially in the Eocene. Subsequently, four Western Eurasian lineages (subsections) differentiated during the Oligocene and Miocene. Evolution of leaf size, form and texture was correlated, in part, with multiple transitions from ancestral humid temperate climates to mediterranean, arid and continental climates. Distantly related but ecologically similar species converged on similar leaf traits in the process.

CONCLUSIONS

Originating in temperate (frost-free) biomes, Eocene to Oligocene ranges of the primarily deciduous cork oaks were restricted to higher latitudes (Siberia to north of Paratethys). Members of the evergreen holly oaks (section Ilex) also originated in temperate biomes but migrated southwards and south-westwards into then-(sub)tropical southern China and south-eastern Tibet during the Eocene, then westwards along existing pre-Himalayan mountain ranges. Divergent biogeographical histories and deep-time phylogenetic legacies (in cold and drought tolerance, nutrient storage and fire resistance) thus account for the modern species mosaic of Western Eurasian oak communities, which are composed of oaks belonging to four sections.

Deep reticulation: the long legacy of hybridization in vascular plant evolution

Hybridization has long been recognized as a fundamental evolutionary process in plants but, until recently, our understanding of its phylogenetic distribution and biological significance across deep evolutionary scales has been largely obscure. Over the past decade, genomic and phylogenomic datasets have revealed, perhaps not surprisingly, that hybridization, often associated with polyploidy, has been common throughout the evolutionary history of plants, particularly in various lineages of flowering plants. However, phylogenomic studies have also highlighted the challenges of disentangling signals of ancient hybridization from other sources of genomic conflict (in particular, incomplete lineage sorting). Here, we provide a critical review of ancient hybridization in vascular plants, outlining well-documented cases of ancient hybridization across plant phylogeny, as well as the challenges unique to documenting ancient versus recent hybridization. We provide a definition for ancient hybridization, which, to our knowledge, has not been explicitly attempted before. Further documenting the extent of deep reticulation in plants should remain an important research focus, especially because published examples likely represent the tip of the iceberg in terms of the total extent of ancient hybridization. However, future research should increasingly explore the macroevolutionary significance of this process, in terms of its impact on evolutionary trajectories (e.g. how does hybridization influence trait evolution or the generation of biodiversity over long time scales?), as well as how life history and ecological factors shape, or have shaped, the frequency of hybridization across geologic time and plant phylogeny. Finally, we consider the implications of ubiquitous ancient hybridization for how we conceptualize, analyze, and classify plant phylogeny. Networks, as opposed to bifurcating trees, represent more accurate representations of evolutionary history in many cases, although our ability to infer, visualize, and use networks for comparative analyses is highly limited. Developing improved methods for the generation, visualization, and use of networks represents a critical future direction for plant biology. Current classification systems also do not generally allow for the recognition of reticulate lineages, and our classifications themselves are largely based on evidence from the chloroplast genome. Updating plant classification to better reflect nuclear phylogenies, as well as considering whether and how to recognize hybridization in classification systems, will represent an important challenge for the plant systematics community.

The use of Anchored Hybrid Enrichment data to resolve higher-level phylogenetic relationships: A proof-of-concept applied to Asterales (Eudicotyledoneae; Angiosperms)

Anchored Hybrid Enrichment (AHE) is a tool for capturing orthologous regions of the nuclear genome shared in low or single copy across lineages. Despite the increasing number of studies using this method, its usefulness to estimate relationships at deeper taxonomic levels in plants has not been fully explored. Here we present a proof of concept about the performance of nuclear loci obtained with AHE to infer phylogenetic relationships and explore the use of gene sampling schemes to estimate divergence times in Asterales. We recovered low-copy nuclear loci using the AHE method from herbarium material and silica-preserved samples. Maximum likelihood, Bayesian inference, and coalescence approaches were used to reconstruct phylogenomic relationships. Dating analyses were conducted under a multispecies coalescent approach by jointly inferring species tree and divergence times with random gene sampling schemes and multiple calibrations. We recovered 403 low-copy nuclear loci for 63 species representing nine out of eleven families of Asterales. Phylogenetic hypotheses were congruent among the applied methods and previously published results. Analyses with concatenated datasets were strongly supported, but coalescence-based analyses showed low support for the phylogenetic position of families Argophyllaceae and Alseuosmiaceae. Estimated family ages were congruent among gene sampling schemes, with the mean age for Asterales around 130 Myr. Our study documents the usefulness of AHE for resolving phylogenetic relationships at deep phylogenetic levels in Asterales. Observed phylogenetic inconsistencies were possibly due to the non-inclusion of families Phellinceae and Pentaphragmataceae. Random gene sampling schemes produced consistent age estimates with coalescence and species tree relaxed clock approaches.

Insights from the genomes of 4 diploid Camelina spp

Plant evolution has been a complex process involving hybridization and polyploidization making understanding the origin and evolution of a plant's genome challenging even once a published genome is available. The oilseed crop, Camelina sativa (Brassicaceae), has a fully sequenced allohexaploid genome with 3 unknown ancestors. To better understand which extant species best represent the ancestral genomes that contributed to C. sativa's formation, we sequenced and assembled chromosome level draft genomes for 4 diploid members of Camelina: C. neglecta C. hispida var. hispida, C. hispida var. grandiflora, and C. laxa using long and short read data scaffolded with proximity data. We then conducted phylogenetic analyses on regions of synteny and on genes described for Arabidopsis thaliana, from across each nuclear genome and the chloroplasts to examine evolutionary relationships within Camelina and Camelineae. We conclude that C. neglecta is closely related to C. sativa's sub-genome 1 and that C. hispida var. hispida and C. hispida var. grandiflora are most closely related to C. sativa's sub-genome 3. Further, the abundance and density of transposable elements, specifically Helitrons, suggest that the progenitor genome that contributed C. sativa's sub-genome 3 maybe more similar to the genome of C. hispida var. hispida than that of C. hispida var. grandiflora. These diploid genomes show few structural differences when compared to C. sativa's genome indicating little change to chromosome structure following allopolyploidization. This work also indicates that C. neglecta and C. hispida are important resources for understanding the genetics of C. sativa and potential resources for crop improvement.

Hyb-Seq provides new insights into the phylogeny and evolution of the Chrysanthemum zawadskii species complex in China

A species complex is an assemblage of closely related species with blurred boundaries, and from which species could arise from different speciation processes and/or a speciation continuum. Such a complex can provide an opportunity to investigate evolutionary mechanisms acting on speciation. The Chrysanthemum zawadskii species complex in China, a monophyletic group of Chrysanthemum, consists of seven species with considerable morphological variation, diverse habitats and different distribution patterns. Here, we used Hyb-Seq data to construct a well-resolved phylogeny of the C. zawadskii complex. Then, we performed comparative analyses of variation patterns in morphology, ecology and distribution to investigate the roles of geography and ecology in this complex's diversification. Lastly, we implemented divergence time estimation, species distribution modelling and ancestral area reconstruction to trace the evolutionary history of this complex. We concluded that the C. zawadskii complex originated in the Qinling-Daba mountains during the early Pliocene and then spread west and northward along the mountain ranges to northern China. During this process, geographical and ecological factors imposing different influences resulted in the current diversification and distribution patterns of this species complex, which is composed of both well-diverged species and diverging lineages on the path of speciation.

Systematics of Thraupis (Aves, Passeriformes) reveals an extensive hybrid zone between T. episcopus (Blue-gray Tanager) and T. sayaca (Sayaca Tanager)

The Neotropical avian genus Thraupis (Passeriformes, Thraupidae) currently comprises seven species that are widespread and abundant throughout their ranges. However, no phylogenetic hypothesis with comprehensive intraspecific sampling is available for the group and, therefore, currently accepted species limits remain untested. We obtained sequence data for two mitochondrial (ND2, cyt-b) and three non-coding nuclear (TGFB2, MUSK, and βF5) markers from 118 vouchered museum specimens. We conducted population structure and coalescent-based species-tree analyses using a molecular clock calibration. We integrated these results with morphometric and coloration analyses of 1,003 museum specimens to assess species limits within Thraupis. Our results confirm that Thraupis is a monophyletic group and support its origin in the late Miocene and subsequent diversification during the Pleistocene. However, we found conflicts with previous phylogenies. We recovered Thraupis glaucocolpa to be sister to all other species in the genus, and T. cyanoptera to the remaining five species. Our phylogenetic trees and population structure analyses uncovered phylogeographic structure within Thraupis episcopus that is congruent with geographic patterns of phenotypic variation and distributions of some named taxa. The first genetic and phenotypic cluster in T. episcopus occurs east of the Andes and is diagnosed by the white patch on the lesser and median wing coverts, whereas the second group has a blue patch on the wing and distributes to the west of Colombia's eastern Andes. Finally, we present evidence of hybridization and ongoing gene flow between several taxa at different taxonomic levels and discuss its taxonomic implications.

A first complete phylogenomic hypothesis for diploid blueberries (Vaccinium section Cyanococcus)

PREMISE

The true blueberries (Vaccinium sect. Cyanococcus; Ericaceae), endemic to North America, have been intensively studied for over a century. However, with species estimates ranging from nine to 24 and much confusion regarding species boundaries, this ecologically and economically valuable group remains inadequately understood at a basic evolutionary and taxonomic level. As a first step toward understanding the evolutionary history and taxonomy of this species complex, we present the first phylogenomic hypothesis of the known diploid blueberries.

METHODS

We used flow cytometry to verify the ploidy of putative diploid taxa and a target-enrichment approach to obtain a genomic data set for phylogenetic analyses.

RESULTS

Despite evidence of gene flow, we found that a primary phylogenetic signal is present. Monophyly for all morphospecies was recovered, with two notable exceptions: one sample of V. boreale was consistently nested in the V. myrtilloides clade and V. caesariense was nested in the V. fuscatum clade. One diploid taxon, Vaccinium pallidum, is implicated as having a homoploid hybrid origin.

CONCLUSIONS

This foundational study represents the first attempt to elucidate evolutionary relationships of the true blueberries of North America with a phylogenomic approach and sets the stage for multiple avenues of future study such as a taxonomic revision of the group, the verification of a homoploid hybrid taxon, and the study of polyploid lineages within the context of a diploid phylogeny.

A snapshot of progenitor-derivative speciation in Iberodes (Boraginaceae)

Traditional classification of speciation modes has focused on physical barriers to gene flow. Allopatric speciation with complete reproductive isolation is viewed as the most common mechanism of speciation. Parapatry and sympatry, by contrast, entail speciation in the face of ongoing gene flow, making them more difficult to detect. The genus Iberodes (Boraginaceae, NW Europe) comprises five species with contrasting morphological traits, habitats, and species distributions. Based on the predominance of narrow and geographically distant endemic species, we hypothesized that geographic barriers were responsible for most speciation events in Iberodes. We undertook an integrative study including: (i) phylogenomics through restriction-site associated DNA sequencing, (ii) genetic structure analyses, (iii) demographic modeling, (iv) morphometrics, and (v) climatic niche modeling and niche overlap analysis. Results revealed a history of recurrent progenitor-derivative speciation manifested by a paraphyletic pattern of nested species differentiation. Budding speciation mediated by ecological differentiation is suggested for the coastal lineage, deriving from the inland widespread I. linifolia during Late Pliocene. Meanwhile, geographic isolation followed by niche shifts are suggested for the more recent differentiation of the coastland taxa. Our work provides a model for distinguishing speciation via ecological differentiation of peripheral, narrowly endemic I. kuzinskyanae and I. littoralis from a widespread extant ancestor, I. linifolia. Ultimately, our results illustrate a case of Pliocene speciation in the probable absence of geographic barriers and get away from the traditional cladistic perspective of speciation as producing two species from an extinct ancestor, thus reminding us that phylogenetic trees tell only part of the story.

A chromosome-scale genome assembly of the Mongolian oak (Quercus mongolica)

Mongolian oak (Quercus mongolica Fisch.) is an ecologically and economically important white oak species native to and widespread in the temperate zone of East Asia. Here, we present a chromosome-scale reference genome assembly of Q. mongolica, a representative white oak species, by combining Illumina and PacBio data with Hi-C mapping technologies that is the first reference genome created for an Asian oak. Our results showed that the PacBio draft genome size was 809.84 Mb, with a BUSCO complete gene percentage of 92.71%. Hi-C scaffolding anchored 774.59 Mb contigs (95.65% of draft assembly) onto 12 pseudochromosomes. The contig N50 and scaffold N50 were 2.64 Mb and 66.74 Mb, respectively. Of the 36,553 protein-coding genes predicted in the study, approximately 95% had functional annotations in public databases. A total of 435.34 Mb (53.75% of the genome) of repetitive sequences were predicted in the assembled genome. Genome evolution analysis showed that Q. mongolica is closely related to Q. robur from Europe, and they shared a common ancestor ~11.8 million years ago. Gene family evolution analysis of Q. mongolica revealed that the nucleotide-binding site (NBS)-encoding gene family related to disease resistance was significantly contracted, whereas the ECERIFERUM 1 (CER1) homologous genes related to cuticular wax biosynthesis was significantly expanded. This pioneering Asian oak genome resource represents an important supplement to the oak genomics community and will improve our understanding of Asian white oak biology and evolution.

Phylogenomic analyses highlight innovation and introgression in the continental radiations of Fagaceae across the Northern Hemisphere

Northern Hemisphere forests changed drastically in the early Eocene with the diversification of the oak family (Fagaceae). Cooling climates over the next 20 million years fostered the spread of temperate biomes that became increasingly dominated by oaks and their chestnut relatives. Here we use phylogenomic analyses of nuclear and plastid genomes to investigate the timing and pattern of major macroevolutionary events and ancient genome-wide signatures of hybridization across Fagaceae. Innovation related to seed dispersal is implicated in triggering waves of continental radiations beginning with the rapid diversification of major lineages and resulting in unparalleled transformation of forest dynamics within 15 million years following the K-Pg extinction. We detect introgression at multiple time scales, including ancient events predating the origination of genus-level diversity. As oak lineages moved into newly available temperate habitats in the early Miocene, secondary contact between previously isolated species occurred. This resulted in adaptive introgression, which may have further amplified the diversification of white oaks across Eurasia.

Fagaceae are diverse family including trees of ecological and economic importance. This phylogenomic analysis of nuclear and plastid genomes reconstructs evolutionary history and finds evidence of multiple adaptive introgression events in this important plant family.

Influence of Pliocene and Pleistocene climates on hybridization patterns between two closely related oak species in China

BACKGROUND AND AIMS

Contemporary patterns of genetic admixture reflect imprints of both ancient and recent gene flow, which can provide us with valuable information on hybridization history in response to palaeoclimate change. Here, we examine the relationships between present admixture patterns and past climatic niche suitability of two East Asian Cerris oaks (Quercus acutissima and Q. chenii) to test the hypothesis that the mid-Pliocene warm climate promoted while the Pleistocene cool climate limited hybridization among local closely related taxa.

METHODS

We analyse genetic variation at seven nuclear microsatellites (1111 individuals) and three chloroplast intergenic spacers (576 individuals) to determine the present admixture pattern and ancient hybridization history. We apply an information-theoretic model selection approach to explore the associations of genetic admixture degree with past climatic niche suitability at multiple spatial scales.

KEY RESULTS

More than 70 % of the hybrids determined by Bayesian clustering analysis and more than 90 % of the individuals with locally shared chloroplast haplotypes are concentrated within a mid-Pliocene contact zone between ~30°N and 35°N. Climatic niche suitabilities for Q. chenii during the mid-Pliocene Warm Period [mPWP, ~3.264-3.025 million years ago (mya)] and during the Last Glacial Maximum (LGM, ~0.022 mya) best explain the admixture patterns across all Q. acutissima populations and across those within the ancient contact zone, respectively.

CONCLUSIONS

Our results highlight that palaeoclimate change shapes present admixture patterns by influencing the extent of historical range overlap. Specifically, the mid-Pliocene warm climate promoted ancient contact, allowing widespread hybridization throughout central China. In contrast, the Pleistocene cool climate caused the local extinction of Q. chenii, reducing the probability of interspecific gene flow in most areas except those sites having a high level of ecological stability.

Analysis of Paralogs in Target Enrichment Data Pinpoints Multiple Ancient Polyploidy Events in Alchemilla s.l. (Rosaceae)

Target enrichment is becoming increasingly popular for phylogenomic studies. Although baits for enrichment are typically designed to target single-copy genes, paralogs are often recovered with increased sequencing depth, sometimes from a significant proportion of loci, especially in groups experiencing whole-genome duplication (WGD) events. Common approaches for processing paralogs in target enrichment data sets include random selection, manual pruning, and mainly, the removal of entire genes that show any evidence of paralogy. These approaches are prone to errors in orthology inference or removing large numbers of genes. By removing entire genes, valuable information that could be used to detect and place WGD events is discarded. Here, we used an automated approach for orthology inference in a target enrichment data set of 68 species of Alchemilla s.l. (Rosaceae), a widely distributed clade of plants primarily from temperate climate regions. Previous molecular phylogenetic studies and chromosome numbers both suggested ancient WGDs in the group. However, both the phylogenetic location and putative parental lineages of these WGD events remain unknown. By taking paralogs into consideration and inferring orthologs from target enrichment data, we identified four nodes in the backbone of Alchemilla s.l. with an elevated proportion of gene duplication. Furthermore, using a gene-tree reconciliation approach, we established the autopolyploid origin of the entire Alchemilla s.l. and the nested allopolyploid origin of four major clades within the group. Here, we showed the utility of automated tree-based orthology inference methods, previously designed for genomic or transcriptomic data sets, to study complex scenarios of polyploidy and reticulate evolution from target enrichment data sets.[Alchemilla; allopolyploidy; autopolyploidy; gene tree discordance; orthology inference; paralogs; Rosaceae; target enrichment; whole genome duplication.].

Quercus Conservation Genetics and Genomics: Past, Present, and Future

Quercus species (oaks) have been an integral part of the landscape in the northern hemisphere for millions of years. Their ability to adapt and spread across different environments and their contributions to many ecosystem services is well documented. Human activity has placed many oak species in peril by eliminating or adversely modifying habitats through exploitative land usage and by practices that have exacerbated climate change. The goal of this review is to compile a list of oak species of conservation concern, evaluate the genetic data that is available for these species, and to highlight the gaps that exist. We compiled a list of 124 Oaks of Concern based on the Red List of Oaks 2020 and the Conservation Gap Analysis for Native U.S. Oaks and their evaluations of each species. Of these, 57% have been the subject of some genetic analysis, but for most threatened species (72%), the only genetic analysis was done as part of a phylogenetic study. While nearly half (49%) of published genetic studies involved population genetic analysis, only 16 species of concern (13%) have been the subject of these studies. This is a critical gap considering that analysis of intraspecific genetic variability and genetic structure are essential for designing conservation management strategies. We review the published population genetic studies to highlight their application to conservation. Finally, we discuss future directions in Quercus conservation genetics and genomics.

An Updated Infrageneric Classification of the North American Oaks (Quercus Subgenus Quercus): Review of the Contribution of Phylogenomic Data to Biogeography and Species Diversity

The oak flora of North America north of Mexico is both phylogenetically diverse and species-rich, including 92 species placed in five sections of subgenus Quercus, the oak clade centered on the Americas. Despite phylogenetic and taxonomic progress on the genus over the past 45 years, classification of species at the subsectional level remains unchanged since the early treatments by WL Trelease, AA Camus, and CH Muller. In recent work, we used a RAD-seq based phylogeny including 250 species sampled from throughout the Americas and Eurasia to reconstruct the timing and biogeography of the North American oak radiation. This work demonstrates that the North American oak flora comprises mostly regional species radiations with limited phylogenetic affinities to Mexican clades, and two sister group connections to Eurasia. Using this framework, we describe the regional patterns of oak diversity within North America and formally classify 62 species into nine major North American subsections within sections Lobatae (the red oaks) and Quercus (the white oaks), the two largest sections of subgenus Quercus. We also distill emerging evolutionary and biogeographic patterns based on the impact of phylogenomic data on the systematics of multiple species complexes and instances of hybridization.

Hybridization and introgression in sympatric and allopatric populations of four oak species

BACKGROUND

Hybridization and introgression are vital sources of novel genetic variation driving diversification during reticulated evolution. Quercus is an important model clade, having extraordinary diverse and abundant members in the Northern hemisphere, that are used to studying the introgression of species boundaries and adaptive processes. China is the second-largest distribution center of Quercus, but there are limited studies on introgressive hybridization.

RESULTS

Here, we screened 17 co-dominant nuclear microsatellite markers to investigate the hybridization and introgression of four oaks (Quercus acutissima, Quercus variabilis, Quercus fabri, and Quercus serrata) in 10 populations. We identified 361 alleles in the four-oak species across 17 loci, and all loci were characterized by high genetic variability (HE = 0.844-0.944) and moderate differentiation (FST = 0.037-0.156) levels. A population differentiation analysis revealed the following: allopatric homologous (FST = 0.064) < sympatric heterogeneous (FST = 0.071) < allopatric heterogeneous (FST = 0.084). A Bayesian admixture analysis determined four types of hybrids (Q. acutissima × Q. variabilis, Q. fabri × Q. serrata, Q. acutissima × Q. fabri, and Q. acutissima × Q. variabilis × Q. fabri) and their asymmetric introgression. Our results revealed that interspecific hybridization is commonly observed within the section Quercus, with members having tendency to hybridize.

CONCLUSIONS

Our study determined the basic hybridization and introgression states among the studied four oak species and extended our understanding of the evolutionary role of hybridization. The results provide useful theoretical data for formulating conservation strategies.

Phylogenetic relationships in Chinese oaks (Fagaceae, Quercus): Evidence from plastid genome using low-coverage whole genome sequencing

China is a second center of oak diversity but with less intensively systematic studies. Here, with 49 species representing all four sections in China, we firstly gave insight into the comprehensive phylogenetic relationships of Chinese oaks based on 54 complete plastid genomes. Our results recovered a robust phylogenetic framework and provided strong support for most nodes. The phylogenetic tree supported Quercus section Ilex as not monophyletic, in which Quercus section Cyclobalanopsis and Quercus section Cerris were nested. Most likely, incomplete lineage sorting and/or introgression among ancestral lineages in these three sections resulted in this complex pattern. The current distribution, diversification and molecular differentiation of Q. sect. Ilex in China are likely consequences of local adaptation to the geographic and paleoclimatic changes, which were driven by the uplift of Tibetan Plateau, the Hengduan Mountains and the Himalayas.

Modern Strategies to Assess and Breed Forest Tree Adaptation to Changing Climate

Studying the genetics of adaptation to new environments in ecologically and industrially important tree species is currently a major research line in the fields of plant science and genetic improvement for tolerance to abiotic stress. Specifically, exploring the genomic basis of local adaptation is imperative for assessing the conditions under which trees will successfully adapt in situ to global climate change. However, this knowledge has scarcely been used in conservation and forest tree improvement because woody perennials face major research limitations such as their outcrossing reproductive systems, long juvenile phase, and huge genome sizes. Therefore, in this review we discuss predictive genomic approaches that promise increasing adaptive selection accuracy and shortening generation intervals. They may also assist the detection of novel allelic variants from tree germplasm, and disclose the genomic potential of adaptation to different environments. For instance, natural populations of tree species invite using tools from the population genomics field to study the signatures of local adaptation. Conventional genetic markers and whole genome sequencing both help identifying genes and markers that diverge between local populations more than expected under neutrality, and that exhibit unique signatures of diversity indicative of "selective sweeps." Ultimately, these efforts inform the conservation and breeding status capable of pivoting forest health, ecosystem services, and sustainable production. Key long-term perspectives include understanding how trees' phylogeographic history may affect the adaptive relevant genetic variation available for adaptation to environmental change. Encouraging "big data" approaches (machine learning-ML) capable of comprehensively merging heterogeneous genomic and ecological datasets is becoming imperative, too.

Conflicting phylogenetic signals in plastomes of the tribe Laureae (Lauraceae)

Background

Gene tree discordance is common in phylogenetic analyses. Many phylogenetic studies have excluded non-coding regions of the plastome without evaluating their impact on tree topology. In general, plastid loci have often been treated as a single unit, and tree discordance among these loci has seldom been examined. Using samples of Laureae (Lauraceae) plastomes, we explored plastome variation among the tribe, examined the influence of non-coding regions on tree topology, and quantified intra-plastome conflict.

Results

We found that the plastomes of Laureae have low inter-specific variation and are highly similar in structure, size, and gene content. Laureae was divided into three groups, subclades I, II and III. The inclusion of non-coding regions changed the phylogenetic relationship among the three subclades. Topologies based on coding and non-coding regions were largely congruent except for the relationship among subclades I, II and III. By measuring the distribution of phylogenetic signal across loci that supported different topologies, we found that nine loci (two coding regions, two introns and five intergenic spacers) played a critical role at the contentious node.

Conclusions

Our results suggest that subclade III and subclade II are successively sister to subclade I. Conflicting phylogenetic signals exist between coding and non-coding regions of Laureae plastomes. Our study highlights the importance of evaluating the influence of non-coding regions on tree topology and emphasizes the necessity of examining discordance among different plastid loci in phylogenetic studies.

High-throughput sequencing of 5S-IGS in oaks: Exploring intragenomic variation and algorithms to recognize target species in pure and mixed samples

Measuring biological diversity is a crucial but difficult undertaking, as exemplified in oaks where complex patterns of morphological, ecological, biogeographical and genetic differentiation collide with traditional taxonomy, which measures biodiversity in number of species (or higher taxa). In this pilot study, we generated high-throughput sequencing amplicon data of the intergenic spacer of the 5S nuclear ribosomal DNA cistron (5S-IGS) in oaks, using six mock samples that differ in geographical origin, species composition and pool complexity. The potential of the marker for automated genotaxonomy applications was assessed using a reference data set of 1,770 5S-IGS cloned sequences, covering the entire taxonomic breadth and distribution range of western Eurasian Quercus, and applying similarity (blast) and evolutionary approaches (maximum-likelihood trees and Evolutionary Placement Algorithm). Both methods performed equally well, allowing correct identification of species in sections Ilex and Cerris in the pure and mixed samples, and main lineages shared by species of sect. Quercus. Application of different cut-off thresholds revealed that medium- to high-abundance (>10 or 25) sequences suffice for a net species identification of samples containing one or a few individuals. Lower thresholds identify phylogenetic correspondence with all target species in highly mixed samples (analogous to environmental bulk samples) and include rare variants pointing towards reticulation, incomplete lineage sorting, pseudogenic 5S units and in situ (natural) contamination. Our pipeline is highly promising for future assessments of intraspecific and interpopulation diversity, and of the genetic resources of natural ecosystems, which are fundamental to empower fast and solid biodiversity conservation programmes worldwide.

Framework Phylogeny, Evolution and Complex Diversification of Chinese Oaks

Oaks (Quercus L.) are ideal models to assess patterns of plant diversity. We integrated the sequence data of five chloroplast and two nuclear loci from 50 Chinese oaks to explore the phylogenetic framework, evolution and diversification patterns of the Chinese oak's lineage. The framework phylogeny strongly supports two subgenera Quercus and Cerris comprising four infrageneric sections Quercus, Cerris, Ilex and Cyclobalanopsis for the Chinese oaks. An evolutionary analysis suggests that the two subgenera probably split during the mid-Eocene, followed by intergroup divergence within the subgenus Cerris around the late Eocene. The initial diversification of sections in the subgenus Cerris was dated between the mid-Oligocene and the Oligocene-Miocene boundary, while a rapid species radiation in section Quercus started in the late Miocene. Diversification simulations indicate a potential evolutionary shift on section Quercus, while several phenotypic shifts likely occur among all sections. We found significant negative correlations between rates of the lineage diversification and phenotypic turnover, suggesting a complex interaction between the species evolution and morphological divergence in Chinese oaks. Our infrageneric phylogeny of Chinese oaks accords with the recently proposed classification of the genus Quercus. The results point to tectonic activity and climatic change during the Tertiary as possible drivers of evolution and diversification in the Chinese oak's lineage.

Correlated evolution of morphology, gas exchange, growth rates and hydraulics as a response to precipitation and temperature regimes in oaks (Quercus)

It is hypothesised that tree distributions in Europe are largely limited by their ability to cope with the summer drought imposed by the Mediterranean climate in the southern areas and by their competitive potential in central regions with more mesic conditions. We investigated the extent to which leaf and plant morphology, gas exchange, leaf and stem hydraulics and growth rates have evolved in a coordinated way in oaks (Quercus) as a result of adaptation to contrasting environmental conditions in this region. We implemented an experiment in which seedlings of 12 European/North African oaks were grown under two watering treatments, a well-watered treatment and a drought treatment in which plants were subjected to three cycles of drought. Consistent with our hypothesis, species from drier summers had traits conferring more tolerance to drought such as small sclerophyllous leaves and lower percent loss of hydraulic conductivity. However, these species did not have lower growth rates as expected by a trade-off with drought tolerance. Overall, our results revealed that climate is an important driver of functional strategies in oaks and that traits have evolved along two coordinated functional axes to adapt to different precipitation and temperature regimes.

Disentangling Sources of Gene Tree Discordance in Phylogenomic Data Sets: Testing Ancient Hybridizations in Amaranthaceae s.l

Gene tree discordance in large genomic data sets can be caused by evolutionary processes such as incomplete lineage sorting and hybridization, as well as model violation, and errors in data processing, orthology inference, and gene tree estimation. Species tree methods that identify and accommodate all sources of conflict are not available, but a combination of multiple approaches can help tease apart alternative sources of conflict. Here, using a phylotranscriptomic analysis in combination with reference genomes, we test a hypothesis of ancient hybridization events within the plant family Amaranthaceae s.l. that was previously supported by morphological, ecological, and Sanger-based molecular data. The data set included seven genomes and 88 transcriptomes, 17 generated for this study. We examined gene-tree discordance using coalescent-based species trees and network inference, gene tree discordance analyses, site pattern tests of introgression, topology tests, synteny analyses, and simulations. We found that a combination of processes might have generated the high levels of gene tree discordance in the backbone of Amaranthaceae s.l. Furthermore, we found evidence that three consecutive short internal branches produce anomalous trees contributing to the discordance. Overall, our results suggest that Amaranthaceae s.l. might be a product of an ancient and rapid lineage diversification, and remains, and probably will remain, unresolved. This work highlights the potential problems of identifiability associated with the sources of gene tree discordance including, in particular, phylogenetic network methods. Our results also demonstrate the importance of thoroughly testing for multiple sources of conflict in phylogenomic analyses, especially in the context of ancient, rapid radiations. We provide several recommendations for exploring conflicting signals in such situations. [Amaranthaceae; gene tree discordance; hybridization; incomplete lineage sorting; phylogenomics; species network; species tree; transcriptomics.]

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.

Genomic landscape of the global oak phylogeny

The tree of life is highly reticulate, with the history of population divergence emerging from populations of gene phylogenies that reflect histories of introgression, lineage sorting and divergence. In this study, we investigate global patterns of oak diversity and test the hypothesis that there are regions of the oak genome that are broadly informative about phylogeny. We utilize fossil data and restriction-site associated DNA sequencing (RAD-seq) for 632 individuals representing nearly 250 Quercus species to infer a time-calibrated phylogeny of the world's oaks. We use a reversible-jump Markov chain Monte Carlo method to reconstruct shifts in lineage diversification rates, accounting for among-clade sampling biases. We then map the > 20 000 RAD-seq loci back to an annotated oak genome and investigate genomic distribution of introgression and phylogenetic support across the phylogeny. Oak lineages have diversified among geographic regions, followed by ecological divergence within regions, in the Americas and Eurasia. Roughly 60% of oak diversity traces back to four clades that experienced increases in net diversification, probably in response to climatic transitions or ecological opportunity. The strong support for the phylogeny contrasts with high genomic heterogeneity in phylogenetic signal and introgression. Oaks are phylogenomic mosaics, and their diversity may in fact depend on the gene flow that shapes the oak genome.

Phylogenomics of the genus Populus reveals extensive interspecific gene flow and balancing selection

Phylogenetic analysis is complicated by interspecific gene flow and the presence of shared ancestral polymorphisms, particularly those maintained by balancing selection. In this study, we aimed to examine the prevalence of these factors during the diversification of Populus, a model tree genus in the Northern Hemisphere. We constructed phylogenetic trees of 29 Populus taxa using 80 individuals based on re-sequenced genomes. Our species tree analyses recovered four main clades in the genus based on consensus nuclear phylogenies, but in conflict with the plastome phylogeny. A few interspecific relationships remained unresolved within the multiple-species clade because of inconsistent gene trees. Our results indicated that gene flow has been widespread within each clade and also occurred among the four clades during their early divergence. We identified 45 candidate genes with ancient polymorphisms maintained by balancing selection. These genes were mainly associated with mating compatibility, growth and stress resistance. Both gene flow and selection-mediated ancient polymorphisms are prevalent in the genus Populus. These are potentially important contributors to adaptive variation. Our results provide a framework for the diversification of model tree genus that will facilitate future comparative studies.

Target Nuclear and Off-Target Plastid Hybrid Enrichment Data Inform a Range of Evolutionary Depths in the Orchid Genus Epidendrum

Universal angiosperm enrichment probe sets designed to enrich hundreds of putatively orthologous nuclear single-copy loci are increasingly being applied to infer phylogenetic relationships of different lineages of angiosperms at a range of evolutionary depths. Studies applying such probe sets have focused on testing the universality and performance of the target nuclear loci, but they have not taken advantage of off-target data from other genome compartments generated alongside the nuclear loci. Here we do so to infer phylogenetic relationships in the orchid genus Epidendrum and closely related genera of subtribe Laeliinae. Our aims are to: 1) test the technical viability of applying the plant anchored hybrid enrichment (AHE) method (Angiosperm v.1 probe kit) to our focal group, 2) mine plastid protein coding genes from off-target reads; and 3) evaluate the performance of the target nuclear and off-target plastid loci in resolving and supporting phylogenetic relationships along a range of taxonomical depths. Phylogenetic relationships were inferred from the nuclear data set through coalescent summary and site-based methods, whereas plastid loci were analyzed in a concatenated partitioned matrix under maximum likelihood. The usefulness of target and flanking non-target nuclear regions and plastid loci was assessed through the estimation of their phylogenetic informativeness. Our study successfully applied the plant AHE probe kit to Epidendrum, supporting the universality of this kit in angiosperms. Moreover, it demonstrated the feasibility of mining plastome loci from off-target reads generated with the Angiosperm v.1 probe kit to obtain additional, uniparentally inherited sequence data at no extra sequencing cost. Our analyses detected some strongly supported incongruences between nuclear and plastid data sets at shallow divergences, an indication of potential lineage sorting, hybridization, or introgression events in the group. Lastly, we found that the per site phylogenetic informativeness of the ycf1 plastid gene surpasses that of all other plastid genes and several nuclear loci, making it an excellent candidate for assessing phylogenetic relationships at medium to low taxonomic levels in orchids.

ddRAD Sequencing-Based Identification of Genomic Boundaries and Permeability in Quercus ilex and Q. suber Hybrids

Hybridization and its relevance is a hot topic in ecology and evolutionary biology. Interspecific gene flow may play a key role in species adaptation to environmental change, as well as in the survival of endangered populations. Despite the fact that hybridization is quite common in plants, many hybridizing species, such as Quercus spp., maintain their integrity, while precise determination of genomic boundaries between species remains elusive. Novel high throughput sequencing techniques have opened up new perspectives in the comparative analysis of genomes and in the study of historical and current interspecific gene flow. In this work, we applied ddRADseq technique and developed an ad hoc bioinformatics pipeline for the study of ongoing hybridization between two relevant Mediterranean oaks, Q. ilex and Q. suber. We adopted a local scale approach, analyzing adult hybrids (sensu lato) identified in a mixed stand and their open-pollinated progenies. We have identified up to 9,435 markers across the genome and have estimated individual introgression levels in adults and seedlings. Estimated contribution of Q. suber to the genome is higher, on average, in hybrid progenies than in hybrid adults, suggesting preferential backcrossing with this parental species, maybe followed by selection during juvenile stages against individuals with higher Q. suber genomic contribution. Most discriminating markers seem to be scattered throughout the genome, suggesting that a large number of small genomic regions underlie boundaries between these species. A noticeable proportion of the markers (26%) showed allelic frequencies in adult hybrids very similar to one of the parental species, and very different from the other; a finding that seems relevant for understanding the hybridization process and the occurrence of adaptive introgression. Candidate marker databases developed in this study constitute a valuable resource to design large scale re-sequencing experiments in Mediterranean sclerophyllous oak species and could provide insight in species boundaries and on adaptive introgression between Q. suber and Q. ilex.

Molecular evolution of the internal transcribed spacers in red oaks (Quercus sect. Lobatae)

Previous studies of the Internal Transcribed Spacers of the nuclear ribosomal DNA (ITS) in sections Quercus (white oaks), Protobalanus (intermediate or golden cup oaks), Cerris (Cerris oaks), and Ilex (Ilex oaks) suggest that ITS regions undergo full concerted evolution in oaks; however, ITS evolution patterns in red oaks (section Lobatae) are unknown due to scant representation in published work. To determine whether full concerted evolution occurs in red oaks, the purpose of this study was to examine ITS sequences from 40 red oak species. The results show incomplete concerted evolution and the presence of three ITS ribotypes of lengths 505, 609, 601 bp, hereafter referred to as ITS-S (small), I ITS-M (medium), and ITS-L (large), respectively. Thirty species had only one ribotype (ITS-M), nine species had two ribotypes (different combinations of ITS-L, ITS-M, and ITS-S), and only one species had all three ribotypes. Furthermore, examination of these three ribotypes showed that only ITS-M is putatively functional and ITS-L and ITS-S are pseudogenes. Bayesian analysis strongly supported (100%) two pseudogenes clades but provided weak support for the monophyly of a putative functional clade (ITS-M); moreover, within the "functional" clade, species relationships were uncertain and, in most cases, sequences from the same species failed to group together. The results of the current study suggest that ITS may not be appropriate for phylogeny reconstruction of red oaks due to low levels of interspecific variation and incomplete concerted evolution.