Genomic landscape of the global oak phylogeny

The influence of post-glacial migration and hybridization on the gene pool of marginal Quercus pubescens populations in Central Europe

BACKGROUND AND AIMS

In Central Europe, the drought-tolerant downy oak (Quercus pubescens) is at the northern edge of its natural distribution range, often growing in small and spatially isolated populations. Here, we elucidate how the population genetic structure of Central European Q. pubescens was shaped by geographical barriers, genetic drift and introgression with the closely related sessile oak (Q. petraea).

METHODS

In total, 27 Q. pubescens populations from the northern margin of its natural distribution range were sampled. Based on 16 nuclear microsatellite markers (nSSRs), Bayesian clustering and distance-based analyses were performed to determine the intraspecific genetic structure and to identify genetic barriers. To identify drivers of introgression with Q. petraea, generalized linear models were applied to link levels of introgression with environmental conditions. To track post-glacial migration routes, the spatial distribution of haplotypes based on eight chloroplast microsatellite markers (cpSSRs) was investigated.

KEY RESULTS

Based on nSSRs, the study populations of Q. pubescens were divided into a western and an eastern genetic cluster. Within these clusters, more pronounced genetic substructure was observed in the west, probably due to a rugged topography and limited gene flow. Introgression from Q. petraea was more prevalent at wetter and north-exposed sites and in the west. The identified cpSSR haplotypes followed known migration pathways.

CONCLUSIONS

Our results suggest two late-glacial refugia in or near the southwestern Alps and the southeastern Alps as potential sources for post-glacial migration. Although some genetic exchange is evident in northern Italy, south of the Alps, the two clusters remain distinct at a large scale. Landscape features and introgression with Q. petraea shaped the genetic substructure at a smaller scale. Our study provides a comprehensive overview of the genetic structure of Q. pubescens in Central Europe, relevant for conservation.

A seasonal strategy for pollen tube growth and ovule development to overcome winter in Japanese stone oak (Lithocarpus edulis)

Delayed fertilization is commonly observed in many acorn-producing Fagaceae trees, yet its underlying mechanisms and adaptive significance remain poorly understood. In recent years, a new hypothesis has been proposed suggesting that the nearly year-long delay in fertilization is driven by an overwintering strategy, wherein female gametophyte development is delayed, and pollen tube growth is arrested before winter. This mechanism allows ovules to be fertilized and seeds to develop during more favorable seasons while avoiding adverse winter conditions. However, empirical evidence for this overwintering strategy has been limited. To address this, we observed the seasonal progression of pollen tube growth and ovule development in Lithocarpus edulis, a species with spring and autumn flowering seasons. Monthly observations of pistillate flowers from both seasons were conducted using microtome techniques and scanning confocal microscopy. Our findings revealed that pollen tubes were arrested at the style joining site, and ovules remained immature in both spring and autumn flowers prior to winter. Following winter, pollen tube regrowth and ovule maturation were synchronized in the subsequent spring, regardless of the flowering season. These results support the hypothesis that ovule development is delayed, leading to delayed fertilization, until after winter. This study highlights the importance of temporally coordinating fertilization phenology with flowering and fruiting phenology in seasonal environments to avoid unfavorable winter conditions.

A Vanishing Imprint? Modeling the Present and Future Distribution of the Enigmatic Quercus crenata Lam., a Mediterranean Sporadic Tree Species

Poorly known, rare species are important biodiversity elements; understanding the relationships between their effective numbers, geographical distribution, ecology, and adaptive potential is an unquestionable critical aspect to reverse biodiversity decline under climate change. Quercus crenata Lam. is a sporadic Mediterranean tree species with debated taxonomy and evolutionary history. Confusing identifications and a scattered distribution combine to create an incautious lack of comprehensive and reliable information on its spatial distribution, ecology, and genetic resources, thereby hindering correct management and conservation efforts. This work undertook a first decisive step to address these knowledge gaps, integrating all previous dispersive information and presenting a comprehensive map of Q. crenata occurrences, with 923 established records: 495 new field observations and 428 verified from all available literature and online databases. The taxon occurs with extremely low individual numbers across central and northern Italy, southern France, western Slovenia, and Croatia, mainly at altitudes between sea level and 1100 m a.s.l. A large part of the species records are outside current networks of protected areas. The MaxEnt-based distribution model highlights Q. crenata adaptation to mild Mediterranean climates with moderate temperature fluctuations, moderate-to-high water requirements, and diverse soil types. A broader present potential distribution than currently assessed is suggested, underlining the possibility of identifying new occurrences with accurate searches in targeted sites, thus refining the understanding of the species' actual present distribution. Future range projections under three carbon emission scenarios with increasing severity (SSP1, SSP3, and SSP5) predict substantial range losses by 2100, ranging from 32% to a drastic 99% reduction under the most severe scenario (SSP5). On these bases, our findings underline the urgent need to improve current conservation practices, which should be conveniently implemented by exhaustive genetic investigations.

Metabolomics and Antioxidant Activity of Valonea from Quercus variabilis Produced in Different Geographical Regions in China

The genus Quercus is widely distributed globally and serves as a potential source of phenolic compounds, which are renowned for their potent biological activities. The primary objective of this study was to determine the concentrations of metabolite components and evaluate the relative antioxidant activities of valonea (acorn cups) from Quercus variabilis (Q. variabilis) of different geographic origins using a UPLC-ESI-MS/MS-based metabolomics approach. A total of 791 metabolite components were identified, with significant variations in their concentrations observed among samples from different geographic locations. Among these, 1-O-galloyl-β-D-glucose was identified as a key active compound. The biosynthesis of galloyl sugars, galactose metabolism, and pathways for starch and sucrose metabolism represent the three pathways that correspond to the differential metabolites, encompassing 23, 11, and 7 metabolites, respectively. The variations in the antioxidant effectiveness of valonea could mainly be linked to the synthesis of galloyl sugars. These findings improve our knowledge of the composition of valonea and offer valuable resources for its extensive utilization and focused development.

An integrative framework reveals widespread gene flow during the early radiation of oaks and relatives in Quercoideae (Fagaceae)

Although the frequency of ancient hybridization across the Tree of Life is greater than previously thought, little work has been devoted to uncovering the extent, timeline, and geographic and ecological context of ancient hybridization. Using an expansive new dataset of nuclear and chloroplast DNA sequences, we conducted a multifaceted phylogenomic investigation to identify ancient reticulation in the early evolution of oaks (Quercus). We document extensive nuclear gene tree and cytonuclear discordance among major lineages of Quercus and relatives in Quercoideae. Our analyses recovered clear signatures of gene flow against a backdrop of rampant incomplete lineage sorting, with gene flow most prevalent among major lineages of Quercus and relatives in Quercoideae during their initial radiation, dated to the Early-Middle Eocene. Ancestral reconstructions including fossils suggest ancestors of Castanea + Castanopsis, Lithocarpus, and the Old World oak clade probably co-occurred in North America and Eurasia, while the ancestors of Chrysolepis, Notholithocarpus, and the New World oak clade co-occurred in North America, offering ample opportunity for hybridization in each region. Our study shows that hybridization-perhaps in the form of ancient syngameons like those seen today-has been a common and important process throughout the evolutionary history of oaks and their relatives. Concomitantly, this study provides a methodological framework for detecting ancient hybridization in other groups.

A haplotype-resolved reference genome of Quercus alba sheds light on the evolutionary history of oaks

White oak (Quercus alba) is an abundant forest tree species across eastern North America that is ecologically, culturally, and economically important. We report the first haplotype-resolved chromosome-scale genome assembly of Q. alba and conduct comparative analyses of genome structure and gene content against other published Fagaceae genomes. We investigate the genetic diversity of this widespread species and the phylogenetic relationships among oaks using whole genome data. Despite strongly conserved chromosome synteny and genome size across Quercus, certain gene families have undergone rapid changes in size, including defense genes. Unbiased annotation of resistance (R) genes across oaks revealed that the overall number of R genes is similar across species - as are the chromosomal locations of R gene clusters - but, gene number within clusters is more labile. We found that Q. alba has high genetic diversity, much of which predates its divergence from other oaks and likely impacts divergence time estimations. Our phylogenetic results highlight widespread phylogenetic discordance across the genus. The white oak genome represents a major new resource for studying genome diversity and evolution in Quercus. Additionally, we show that unbiased gene annotation is key to accurately assessing R gene evolution in Quercus.

Widespread slow growth of acquisitive tree species

Trees are an important carbon sink as they accumulate biomass through photosynthesis1. Identifying tree species that grow fast is therefore commonly considered to be essential for effective climate change mitigation through forest planting. Although species characteristics are key information for plantation design and forest management, field studies often fail to detect clear relationships between species functional traits and tree growth2. Here, by consolidating four independent datasets and classifying the acquisitive and conservative species based on their functional trait values, we show that acquisitive tree species, which are supposedly fast-growing species, generally grow slowly in field conditions. This discrepancy between the current paradigm and field observations is explained by the interactions with environmental conditions that influence growth. Acquisitive species require moist mild climates and fertile soils, conditions that are generally not met in the field. By contrast, conservative species, which are supposedly slow-growing species, show generally higher realized growth due to their ability to tolerate unfavourable environmental conditions. In general, conservative tree species grow more steadily than acquisitive tree species in non-tropical forests. We recommend planting acquisitive tree species in areas where they can realize their fast-growing potential. In other regions, where environmental stress is higher, conservative tree species have a larger potential to fix carbon in their biomass.

Intraspecific character displacement in oaks

Character displacement refers to the process by which species diverge more in sympatry due to competition for resources. This competition-driven speciation can also occur within populations, known as intraspecific character displacement (ICD). ICD can promote divergence within species by influencing intraspecific competition or encouraging the evolution of alternative phenotypes. Despite its significance, ICD remains understudied and requires further exploration. In this study, we investigate how competition influences genetic and morphological differentiation within species in sympatric and allopatric populations. We focused on Quercus serrata (in China and Japan) and Q. serrata var. brevipetiolata (found only in China), which belong to a small monophyletic group of oak species nested within Section Quercus (white oaks). Using genetic markers, we detected divergence between Chinese and Japanese populations and further diversification within China, with asymmetric historical gene flow primarily from Q. serrata (the earlier diverged species) to Q. serrata var. brevipetiolata (the later variety). Although genetic differentiation did not differ between sympatric and allopatric populations, leaf morphological variation, analyzed through the geometric morphometric method (GMM) and traditional morphological method, revealed greater trait variation in sympatry. In addition, we found an allometric growth relationship between leaf size and leaf mass of Q. serrata and Q. serrata var. brevipetiolata, with the leaf area of Q. serrata var. brevipetiolata decreasing more disproportionately to leaf mass. This suggests a resource trade-off, where Q. serrata var. brevipetiolata, the later diverged variety, adopts more resource-conservative traits in sympatry. Further analysis of trait variation with environmental factors supports these findings, while genetic variation along climate gradients showed significant responses primarily in Q. serrata, regardless of sympatric or allopatric conditions. Although neutral genetic markers are insufficient to capture selection-driven adaptive differentiation, we inferred that Q. serrata var. brevipetiolata is progressing towards ecological divergence from Q. serrata. Overall, our results highlight the role of ICD in driving morphological diversification and resource-use strategies within species in response to competitive pressures.

Adaptive Evolution of Freezing Tolerance in Oaks Is Key to Their Dominance in North America

Freezing tolerance plays a pivotal role in shaping the distribution and diversification of organisms. We investigated the dynamics of adaptation to climate and potential trade-offs between stem freezing tolerance and growth rate in 48 Quercus species. Species from colder regions exhibited higher freezing tolerance, lower growth rates and higher winter-acclimation potential than species from warmer climates. Despite an evolutionary lag, freezing tolerance in oaks is closely aligned with its optimal state. Deciduous species showed marked variability in freezing tolerance across their broad climatic range, while evergreen species, confined to warm climates, displayed low freezing tolerance. Annual growth rates were constrained in all deciduous species, but those that evolved in warm latitudes lost freezing tolerance, precluding a trade-off between freezing tolerance and growth. We provide evidence that the capacity to adapt to a wide range of thermal environments was critical to adaptive radiation and the current dominance of the North American oaks.

Response of Pedunculate Oak (Quercus robur L.) to Adverse Environmental Conditions in Genetic and Dendrochronological Studies

Pedunculate oak (Quercus robur L.) is widely distributed across Europe and serves critical ecological, economic, and recreational functions. Investigating its responses to stressors such as drought, extreme temperatures, pests, and pathogens provides valuable insights into its capacity to adapt to climate change. Genetic and dendrochronological studies offer complementary perspectives on this adaptability. Tree-ring analysis (dendrochronology) reveals how Q. robur has historically responded to environmental stressors, linking growth patterns to specific conditions such as drought or temperature extremes. By examining tree-ring width, density, and dynamics, researchers can identify periods of growth suppression or enhancement and predict forest responses to future climatic events. Genetic studies further complement this by uncovering adaptive genetic diversity and inheritance patterns. Identifying genetic markers associated with stress tolerance enables forest managers to prioritize the conservation of populations with higher adaptive potential. These insights can guide reforestation efforts and support the development of climate-resilient oak populations. By integrating genetic and dendrochronological data, researchers gain a holistic understanding of Q. robur's mechanisms of resilience. This knowledge is vital for adaptive forest management and sustainable planning in the face of environmental challenges, ultimately helping to ensure the long-term viability of oak populations and their ecosystems. The topics covered in this review are very broad. We tried to include the most relevant, important, and significant studies, but focused mainly on the relatively recent Eastern European studies because they include the most of the species' area. However, although more than 270 published works have been cited in this review, we have, of course, missed some published studies. We apologize in advance to authors of those relevant works that have not been cited.

Genomic insights into ecological adaptation of oaks revealed by phylogenomic analysis of multiple species

Understanding the ecological adaptation of tree species can not only reveal the evolutionary potential but also benefit biodiversity conservation under global climate change. Quercus is a keystone genus in Northern Hemisphere forests, and its wide distribution in diverse ecosystems and long evolutionary history make it an ideal model for studying the genomic basis of ecological adaptations. Here we used a newly sequenced genome of Quercus gilva, an evergreen oak species from East Asia, with 18 published Fagales genomes to determine how Fagaceae genomes have evolved, identify genomic footprints of ecological adaptability in oaks in general, as well as between evergreen and deciduous oaks. We found that oak species exhibited a higher degree of genomic conservation and stability, as indicated by the absence of large-scale chromosomal structural variations or additional whole-genome duplication events. In addition, we identified expansion and tandem repetitions within gene families that contribute to plant physical and chemical defense (e.g., cuticle biosynthesis and oxidosqualene cyclase genes), which may represent the foundation for the ecological adaptation of oak species. Circadian rhythm and hormone-related genes may regulate the habits of evergreen and deciduous oaks. This study provides a comprehensive perspective on the ecological adaptations of tree species based on phylogenetic, genome evolutionary, and functional genomic analyses.

Heterogeneous occurrence of evergreen broad-leaved forests in East Asia: Evidence from plant fossils

Evergreen broad-leaved forests (EBLFs) are widely distributed in East Asia and play a vital role in ecosystem stability. The occurrence of these forests in East Asia has been a subject of debate across various disciplines. In this study, we explored the occurrence of East Asian EBLFs from a paleobotanical perspective. By collecting plant fossils from four regions in East Asia, we have established the evolutionary history of EBLFs. Through floral similarity analysis and paleoclimatic reconstruction, we have revealed a diverse spatio-temporal pattern for the occurrence of EBLFs in East Asia. The earliest occurrence of EBLFs in southern China can be traced back to the middle Eocene, followed by southwestern China during the late Eocene-early Oligocene. Subsequently, EBLFs emerged in Japan during the early Oligocene and eventually appeared in central-eastern China around the Miocene. Paleoclimate simulation results suggest that the precipitation of wettest quarter (PWetQ, mm) exceeding 600 mm is crucial for the occurrence of EBLFs. Furthermore, the heterogeneous occurrence of EBLFs in East Asia is closely associated with the evolution of the Asian Monsoon. This study provides new insights into the occurrence of EBLFs in East Asia.

Plastome characterization and its phylogenetic implications on Lithocarpus (Fagaceae)

BACKGROUND

The genus Lithocarpus is a species-rich dominant woody lineage in East Asian evergreen broad-leaved forests. Despite its ecological and economic significance, the plastome structure and evolutionary history of the genus remain poorly understood. In this study, we comprehensively analyzed the 34 plastomes representing 33 Lithocarpus species. Of which, 21 were newly assembled. The plastome-based phylogenomic tree was reconstructed to reveal the maternal evolutionary patterns of the genus.

RESULTS

The Lithocarpus plastomes exhibit a typical quadripartite structure, ranging in length from 161,010 to 161,476 bp, and containing 131 genes, including 86 protein-coding genes, 8 rRNA genes, and 37 tRNA genes. Remarkably, the infA gene was identified as a pseudogene in 17 species. Significant variability was observed in simple sequence repeats (SSRs) as well as in the boundary regions between the two single-copy regions and the inverted repeat region (SC/IR) across the plastomes. Additionally, four genes (accD, atpF, rpl32, and rps8) were found to be under positive selection. The monophyletic status of Lithocarpus was strongly supported by plastome-based phylogeny; however, the phylogenetic tree topology showed a significant difference from that obtained by the nuclear genome-based phylogeny.

CONCLUSIONS

The plastome of Fagaceae is generally conserved. Nevertheless, genes related to metabolism, photosynthesis, and energy were under strong positive selection in Lithocarpus, likely driven by environmental pressures and local adaptation. The plastome-based phylogeny confirmed the monophyletic status of Lithocarpus and revealed a phylogeographic pattern indicating limited seed-mediated gene flow in the ancestral lineage. The prevalence of cytonuclear discordance in Lithocarpus and other Fagaceae genera suggests that ancient introgression, incomplete lineage sorting, and asymmetrical seed- and pollen-mediated geneflow might contribute to this discordance. Future studies are essential to test these hypotheses and further elucidate the divergence patterns of this unique Asian evergreen lineage.

Stomatal closure as a driver of minimum leaf conductance declines at high temperature and vapor pressure deficit in Quercus

Rising global temperatures and vapor pressure deficits (VPDs) are increasing plant water demand and becoming major drivers of large-scale plant mortality. Controlling transient leaf water loss after stomatal closure (minimum stomatal conductance [gmin]) is recognized as a key trait determining how long plants survive during soil drought. Yet, substantial uncertainty remains regarding how gmin responds to elevated temperatures and VPD and the underlying mechanisms. We measured gmin in 24 Quercus species from temperate and Mediterranean climates to determine whether gmin was sensitive to a coupled temperature and VPD increase. We also explored mechanistic links to phenology, climate, evolutionary history, and leaf anatomy. We found that gmin in all species exhibited a nonlinear negative temperature and VPD dependence. At 25 °C (VPD = 2.2 kPa), gmin varied from 1.19 to 8.09 mmol m-2 s-1 across species but converged to 0.57 ± 0.06 mmol m-2 s-1 at 45 °C (VPD = 6.6 kPa). In a subset of species, the effect of temperature and VPD on gmin was reversible and linked to the degree of stomatal closure, which was greater at 45 °C than at 25 °C. Our results show that gmin is dependent on temperature and VPD, is highly conserved in Quercus species, and is linked to leaf anatomy and stomatal behavior.

Applications of CRISPR Technologies in Forestry and Molecular Wood Biotechnology

Forests worldwide are under increasing pressure from climate change and emerging diseases, threatening their vital ecological and economic roles. Traditional breeding approaches, while valuable, are inherently slow and limited by the long generation times and existing genetic variation of trees. CRISPR technologies offer a transformative solution, enabling precise and efficient genome editing to accelerate the development of climate-resilient and productive forests. This review provides a comprehensive overview of CRISPR applications in forestry, exploring its potential for enhancing disease resistance, improving abiotic stress tolerance, modifying wood properties, and accelerating growth. We discuss the mechanisms and applications of various CRISPR systems, including base editing, prime editing, and multiplexing strategies. Additionally, we highlight recent advances in overcoming key challenges such as reagent delivery and plant regeneration, which are crucial for successful implementation of CRISPR in trees. We also delve into the potential and ethical considerations of using CRISPR gene drive for population-level genetic alterations, as well as the importance of genetic containment strategies for mitigating risks. This review emphasizes the need for continued research, technological advancements, extensive long-term field trials, public engagement, and responsible innovation to fully harness the power of CRISPR for shaping a sustainable future for forests.

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.

New insights on the phylogeny, evolutionary history, and ecological adaptation mechanism in cycle-cup oaks based on chloroplast genomes

Cycle-cup oaks (Quercus section Cyclobalanopsis) are one of the principal components of forests in the tropical and subtropical climates of East and Southeast Asia. They have experienced relatively recent increases in the diversification rate, driven by changing climates and the Himalayan orogeny. However, the evolutionary history and adaptive mechanisms at the chloroplast genome level in cycle-cup oaks remain largely unknown. Therefore, we studied this problem by conducting chloroplast genomics on 50 of the ca. 90 species. Comparative genomics and other analyses showed that Quercus section Cyclobalanopsis had a highly conserved chloroplast genome structure. Highly divergent regions, such as the ndhF and ycf1 gene regions and the petN-psbM and rpoB-trnC-GCA intergenic spacer regions, provided potential molecular markers for subsequent analysis. The chloroplast phylogenomic tree indicated that Quercus section Cyclobalanopsis was not monophyletic, which mixed with the other two sections of subgenus Cerris. The reconstruction of ancestral aera inferred that Palaeotropics was the most likely ancestral range of Quercus section Cyclobalanopsis, and then dispersed to Sino-Japan and Sino-Himalaya. Positive selection analysis showed that the photosystem genes had the lowest ω values among the seven functional gene groups. And nine protein-coding genes containing sites for positive selection: ndhA, ndhD, ndhF, ndhH, rbcL, rpl32, accD, ycf1, and ycf2. This series of analyses together revealed the phylogeny, evolutionary history, and ecological adaptation mechanism of the chloroplast genome of Quercus section Cyclobalanopsis in the long river of earth history. These chloroplast genome data provide valuable information for deep insights into phylogenetic relationships and intraspecific diversity in Quercus.

Comparative and phylogenetic analysis of chloroplast genomes from ten species in Quercus section Cyclobalanopsis

The genus Quercus L. is widely acknowledged as a significant assemblage within East Asia tropical and subtropical broadleaf evergreen forests, possessing considerable economic importance. Nevertheless, the differentiation of Quercus species is deemed arduous, and the interrelations among these species remain enigmatic. Leveraging Illumina sequencing, we undertook the sequencing and assembly of the chloroplast (cp) genomes of seven species belonging to Quercus section Cyclobalanopsis (Quercus argyrotricha, Q. augustinii, Q. bambusifolia, Q. bella, Q. edithiae, Q. jenseniana, and Q. poilanei). Furthermore, we collated three previously published cp genome sequences of Cyclobalanopsis species (Q. litseoides, Q. obovatifolia, and Q. saravanensis). Our primary objective was to conduct comparative genomics and phylogenetic analyses of the complete cp genomes of ten species from Quercus section Cyclobalanopsis. This investigation unveiled that Quercus species feature a characteristic circular tetrad structure, with genome sizes ranging from 160,707 to 160,999 base pairs. The genomic configuration, GC content, and boundaries of inverted repeats/single copy regions exhibited marked conservation. Notably, four highly variable hotspots were identified in the comparative analysis, namely trnK-rps16, psbC-trnS, rbcL-accD, and ycf1. Furthermore, three genes (atpF, rpoC1, and ycf2) displayed signals of positive selection pressure. Phylogenetic scrutiny revealed that the four sections of Cyclobalanopsis clustered together as sister taxa. The branch support values ranged from moderate to high, with most nodes garnering 100% support, underscoring the utility of cp genomic data in elucidating the relationships within the genus. Divergence time analysis revealed that Section Cyclobalanopsis represents the earliest type of Quercus genus. The outcomes of this investigation establish a foundation for forthcoming research endeavors in taxonomy and phylogenetics.

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.

Relationships within Bolbitis sinensis Species Complex Using RAD Sequencing

Species identification and phylogenetic relationship clarification are fundamental goals in species delimitation. However, these tasks pose challenges when based on morphologies, geographic distribution, and genomic data. Previously, two species of the fern genus Bolbitis, B. × multipinna and B. longiaurita were described based on morphological traits; they are phylogenetically intertwined with B. sinensis and fail to form monophyletic groups. To address the unclear phylogenetic relationships within the B. sinensis species complex, RAD sequencing was performed on 65 individuals from five populations. Our integrated analysis of phylogenetic trees, neighbor nets, and genetic structures indicate that the B. sinensis species complex should not be considered as separate species. Moreover, our findings reveal differences in the degree of genetic differentiation among the five populations, ranging from low to moderate, which might be influenced by geographical distance and gene flow. The Fst values also confirmed that genetic differentiation intensifies with increasing geographic distance. Collectively, this study clarifies the complex phylogenetic relationships within the B. sinensis species complex, elucidates the genetic diversity and differentiation across the studied populations, and offers valuable genetic insights that contribute to the broader study of evolutionary relationships and population genetics within the Bolbitis species.

Challenge accepted: Evolutionary lineages versus taxonomic classification of North American shrub willows (Salix)

PREMISE

The huge diversity of Salix subgenus Chamaetia/Vetrix clade in North America and the lack of phylogenetic resolution within this clade has presented a difficult but fascinating challenge for taxonomists to resolve. Here we tested the existing taxonomic classification with molecular tools.

METHODS

In this study, 132 samples representing 46 species from 22 described sections of shrub willows from the United States and Canada were analyzed and combined with 67 samples from Eurasia. The ploidy levels of the samples were determined using flow cytometry and nQuire. Sequences were produced using a RAD sequencing approach and subsequently analyzed with ipyrad, then used for phylogenetic reconstructions (RAxML, SplitsTree), dating analyses (BEAST, SNAPPER), and character evolution analyses of 14 selected morphological traits (Mesquite).

RESULTS

The RAD sequencing approach allowed the production of a well-resolved phylogeny of shrub willows. The resulting tree showed an exclusively North American (NA) clade in sister position to a Eurasian clade, which included some North American endemics. The NA clade began to diversify in the Miocene. Polyploid species appeared in each observed clade. Character evolution analyses revealed that adaptive traits such as habit and adaxial nectaries evolved multiple times independently.

CONCLUSIONS

The diversity in shrub willows was shaped by an evolutionary radiation in North America. Most species were monophyletic, but the existing sectional classification could not be supported by molecular data. Nevertheless, monophyletic lineages share several morphological characters, which might be useful in the revision of the taxonomic classification of shrub willows.

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.

Assembly and analysis of the genome of Notholithocarpus densiflorus

Tanoak (Notholithocarpus densiflorus) is an evergreen tree in the Fagaceae family found in California and southern Oregon. Historically, tanoak acorns were an important food source for Native American tribes, and the bark was used extensively in the leather tanning process. Long considered a disjunct relictual element of the Asian stone oaks (Lithocarpus spp.), phylogenetic analysis has determined that the tanoak is an example of convergent evolution. Tanoaks are deeply divergent from oaks (Quercus) of the Pacific Northwest and comprise a new genus with a single species. These trees are highly susceptible to "sudden oak death" (SOD), a plant pathogen (Phytophthora ramorum) that has caused widespread deaths of tanoaks. In this study, we set out to assemble the genome and perform comparative studies among a number of individuals that demonstrated varying levels of susceptibility to SOD. First, we sequenced and de novo assembled a draft reference genome of N. densiflorus using cobarcoded library processing methods and an MGI DNBSEQ-G400 sequencer. To increase the contiguity of the final assembly, we also sequenced Oxford Nanopore long reads to 30× coverage. To our knowledge, the draft genome reported here is one of the more contiguous and complete genomes of a tree species published to date, with a contig N50 of ∼1.2 Mb, a scaffold N50 of ∼2.1 Mb, and a complete gene score of 95.5% through BUSCO analysis. In addition, we sequenced 11 genetically distinct individuals and mapped these onto the draft reference genome, enabling the discovery of almost 25 million single nucleotide polymorphisms and ∼4.4 million small insertions and deletions. Finally, using cobarcoded data, we were able to generate a complete haplotype coverage of all 11 genomes.

Rising Temperatures, Falling Leaves: Predicting the Fate of Cyprus's Endemic Oak under Climate and Land Use Change

Endemic island species face heightened extinction risk from climate-driven shifts, yet standard models often underestimate threat levels for those like Quercus alnifolia, an iconic Cypriot oak with pre-adaptations to aridity. Through species distribution modelling, we investigated the potential shifts in its distribution under future climate and land-use change scenarios. Our approach uniquely combines dispersal constraints, detailed soil characteristics, hydrological factors, and anticipated soil erosion data, offering a comprehensive assessment of environmental suitability. We quantified the species' sensitivity, exposure, and vulnerability to projected changes, conducting a preliminary IUCN extinction risk assessment according to Criteria A and B. Our projections uniformly predict range reductions, with a median decrease of 67.8% by the 2070s under the most extreme scenarios. Additionally, our research indicates Quercus alnifolia's resilience to diverse erosion conditions and preference for relatively dry climates within a specific annual temperature range. The preliminary IUCN risk assessment designates Quercus alnifolia as Critically Endangered in the future, highlighting the need for focused conservation efforts. Climate and land-use changes are critical threats to the species' survival, emphasising the importance of comprehensive modelling techniques and the urgent requirement for dedicated conservation measures to safeguard this iconic species.

Karyotype and LTR-RTs analysis provide insights into oak genomic evolution

BACKGROUND

Whole-genome duplication and long terminal repeat retrotransposons (LTR-RTs) amplification in organisms are essential factors that affect speciation, local adaptation, and diversification of organisms. Understanding the karyotype projection and LTR-RTs amplification could contribute to untangling evolutionary history. This study compared the karyotype and LTR-RTs evolution in the genomes of eight oaks, a dominant lineage in Northern Hemisphere forests.

RESULTS

Karyotype projections showed that chromosomal evolution was relatively conservative in oaks, especially on chromosomes 1 and 7. Modern oak chromosomes formed through multiple fusions, fissions, and rearrangements after an ancestral triplication event. Species-specific chromosomal rearrangements revealed fragments preserved through natural selection and adaptive evolution. A total of 441,449 full-length LTR-RTs were identified from eight oak genomes, and the number of LTR-RTs for oaks from section Cyclobalanopsis was larger than in other sections. Recent amplification of the species-specific LTR-RTs lineages resulted in significant variation in the abundance and composition of LTR-RTs among oaks. The LTR-RTs insertion suppresses gene expression, and the suppressed intensity in gene regions was larger than in promoter regions. Some centromere and rearrangement regions indicated high-density peaks of LTR/Copia and LTR/Gypsy. Different centromeric regional repeat units (32, 78, 79 bp) were detected on different Q. glauca chromosomes.

CONCLUSION

Chromosome fusions and arm exchanges contribute to the formation of oak karyotypes. The composition and abundance of LTR-RTs are affected by its recent amplification. LTR-RTs random retrotransposition suppresses gene expression and is enriched in centromere and chromosomal rearrangement regions. This study provides novel insights into the evolutionary history of oak karyotypes and the organization, amplification, and function of LTR-RTs.

High-quality haplotype-resolved genome assembly for ring-cup oak (Quercus glauca) provides insight into oaks demographic dynamics

Quercus section Cyclobalanopsis represents a dominant woody lineage in East Asian evergreen broadleaved forests. Regardless of its ecological and economic importance, little is known about the genomes of species in this unique oak lineage. Quercus glauca is one of the most widespread tree species in the section Cyclobalanopsis. In this study, a high-quality haplotype-resolved reference genome was assembled for Q. glauca from PacBio HiFi and Hi-C reads. The genome size, contig N50, and scaffold N50 measured 902.88, 7.60, and 69.28 Mb, respectively, for haplotype1, and 913.28, 7.20, and 71.53 Mb, respectively, for haplotype2. A total of 37,457 and 38,311 protein-coding genes were predicted in haplotype1 and haplotype2, respectively. Homologous chromosomes in the Q. glauca genome had excellent gene pair collinearity. The number of R-genes in Q. glauca was similar to most East Asian oaks but less than oak species from Europe and America. Abundant structural variation in the Q. glauca genome could contribute to environmental stress tolerance in Q. glauca. Sections Cyclobalanopsis and Cerris diverged in the Oligocene, in agreement with fossil records for section Cyclobalanopsis, which document its presence in East Asia since the early Miocene. The demographic dynamics of closely related oak species were largely similar. The high-quality reference genome provided here for the most widespread species in section Cyclobalanopsis will serve as an essential genomic resource for evolutionary studies of key oak lineages while also supporting studies of interspecific introgression, local adaptation, and speciation in oaks.

Effect of geographical origin on the chemical characteristics of Q. canariensis acorns: profiling fatty acids, tocopherols, and phenolic composition

The impact of geographical origin on the chemical composition of acorns from Quercus species has significant attention. This study aimed to explore the phenolic composition of methanolic extracts, tocopherol content, and fatty acid composition of acorn oils from six different populations of Q. canariensis acorns. The obtained results revealed that acorn oil from BniMtir exhibited high levels of α-tocopherol (58 mg/kg). The fatty acids identified across all samples were Z-vaccenic + oleic acids (38.44-58.58%). In addition, the data highlighted the presence of quinic (32.514-60.216 µg/g DW) and gallic acids (1.674-10.849 µg/g DW), as well as catechin (+) (0.096-12.647 µg/g DW) in all populations. These variations in chemical composition from different regions are likely linked to their geographical origin. In conclusion, this study's finding should significance for the industry, offering valuable insights into the potential production of bioactive compounds from Q. canariensis acorns, which could have various applications.

Ex situ conservation of two rare oak species using microsatellite and SNP markers

Plant collections held by botanic gardens and arboreta are key components of ex situ conservation. Maintaining genetic diversity in such collections allows them to be used as resources for supplementing wild populations. However, most recommended minimum sample sizes for sufficient ex situ genetic diversity are based on microsatellite markers, and it remains unknown whether these sample sizes remain valid in light of more recently developed next-generation sequencing (NGS) approaches. To address this knowledge gap, we examine how ex situ conservation status and sampling recommendations differ when derived from microsatellites and single nucleotide polymorphisms (SNPs) in garden and wild samples of two threatened oak species. For Quercus acerifolia, SNPs show lower ex situ representation of wild allelic diversity and slightly lower minimum sample size estimates than microsatellites, while results for each marker are largely similar for Q. boyntonii. The application of missing data filters tends to lead to higher ex situ representation, while the impact of different SNP calling approaches is dependent on the species being analyzed. Measures of population differentiation within species are broadly similar between markers, but larger numbers of SNP loci allow for greater resolution of population structure and clearer assignment of ex situ individuals to wild source populations. Our results offer guidance for future ex situ conservation assessments utilizing SNP data, such as the application of missing data filters and the usage of a reference genome, and illustrate that both microsatellites and SNPs remain viable options for botanic gardens and arboreta seeking to ensure the genetic diversity of their collections.

Complete Chloroplast Genomes of Four Oaks from the Section Cyclobalanopsis Improve the Phylogenetic Analysis and Understanding of Evolutionary Processes in the Genus Quercus

Quercus is a valuable genus ecologically, economically, and culturally. They are keystone species in many ecosystems. Species delimitation and phylogenetic studies of this genus are difficult owing to frequent hybridization. With an increasing number of genetic resources, we will gain a deeper understanding of this genus. In the present study, we collected four Quercus section Cyclobalanopsis species (Q. poilanei, Q. helferiana, Q. camusiae, and Q. semiserrata) distributed in Southeast Asia and sequenced their complete genomes. Following analysis, we compared the results with those of other species in the genus Quercus. These four chloroplast genomes ranged from 160,784 bp (Q. poilanei) to 161,632 bp (Q. camusiae) in length, with an overall guanine and cytosine (GC) content of 36.9%. Their chloroplast genomic organization and order, as well as their GC content, were similar to those of other Quercus species. We identified seven regions with relatively high variability (rps16, ndhk, accD, ycf1, psbZ-trnG-GCC, rbcL-accD, and rpl32-trnL-UAG) which could potentially serve as plastid markers for further taxonomic and phylogenetic studies within Quercus. Our phylogenetic tree supported the idea that the genus Quercus forms two well-differentiated lineages (corresponding to the subgenera Quercus and Cerris). Of the three sections in the subgenus Cerris, the section Ilex was split into two clusters, each nested in the other two sections. Moreover, Q. camusiae and Q. semiserrata detected in this study diverged first in the section Cyclobalanopsis and mixed with Q. engleriana in the section Ilex. In particular, 11 protein coding genes (atpF, ndhA, ndhD, ndhF, ndhK, petB, petD, rbcL, rpl22, ycf1, and ycf3) were subjected to positive selection pressure. Overall, this study enriches the chloroplast genome resources of Quercus, which will facilitate further analyses of phylogenetic relationships in this ecologically important tree genus.

Molecular identification and functional characterization of an environmental stress responsive glutaredoxin gene ROXY1 in Quercus glauca

Quercus glauca is a valuable natural resource with both economic and ecological values. It is one of the dominant forest tree species widely distributed in Southern China. As a perennial broadleaf plant, Q. glauca inevitably encounters numerous stresses from environment. Glutaredoxins (GRXs) are a kind of small oxidoreductases that play an important role in response to oxidative stress. CC-type GRXs also known as ROXYs are specific to land plants. In this study, we isolated a CC-type GRX gene, QgROXY1, from Q. glauca. Expression of QgROXY1 is induced by a variety of environmental stimuli. QgROXY1 protein localizes to both cytoplasm and nucleus; whereas the nucleus localized QgROXY1 could physically interact with the basic region/leucine zipper motif (bZIP) transcription factor AtTGA2 from Arabidopsis thaliana. Transgenic A. thaliana ectopically expressing QgROXY1 is hypersensitive to exogenously applied salicylic acid. Induction of plant defense gene is significantly impaired in QgROXY1 transgenic plants that results in enhanced susceptibility to infection of Botrytis cinerea pathogen, indicating the evolutionary conserved function among ROXY homologs in weedy and woody plants. This is the first described function for the ROXYs in tree plants. Through this case study, we demonstrated the feasibility and efficacy of molecular technology applied to characterization of gene function in tree species.

Species identification through deep learning and geometrical morphology in oaks (Quercus spp.): Pros and cons

Plant phenotypic characteristics, especially leaf morphology of leaves, are an important indicator for species identification. However, leaf shape can be extraordinarily complex in some species, such as oaks. The great variation in leaf morphology and difficulty of species identification in oaks have attracted the attention of scientists since Charles Darwin. Recent advances in discrimination technology have provided opportunities to understand leaf morphology variation in oaks. Here, we aimed to compare the accuracy and efficiency of species identification in two closely related deciduous oaks by geometric morphometric method (GMM) and deep learning using preliminary identification of simple sequence repeats (nSSRs) as a prior. A total of 538 Asian deciduous oak trees, 16 Q. aliena and 23 Q. dentata populations, were firstly assigned by nSSRs Bayesian clustering analysis to one of the two species or admixture and this grouping served as a priori identification of these trees. Then we analyzed the shapes of 2328 leaves from the 538 trees in terms of 13 characters (landmarks) by GMM. Finally, we trained and classified 2221 leaf-scanned images with Xception architecture using deep learning. The two species can be identified by GMM and deep learning using genetic analysis as a priori. Deep learning is the most cost-efficient method in terms of time-consuming, while GMM can confirm the admixture individuals' leaf shape. These various methods provide high classification accuracy, highlight the application in plant classification research, and are ready to be applied to other morphology analysis.

Tracking population structure and phenology through time using ancient genomes from waterlogged white oak wood

Whole genome characterizations of crop plants based on ancient DNA have provided unique keys for a better understanding of the evolutionary origins of modern cultivars, the pace and mode of selection underlying their adaptation to new environments and the production of phenotypes of interest. Although forests are among the most biologically rich ecosystems on earth and represent a fundamental resource for human societies, no ancient genome sequences have been generated for trees. This contrasts with the generation of multiple ancient reference genomes for important crops. Here, we sequenced the first ancient tree genomes using two white oak wood remains from Germany dating to the Last Little Ice Age (15th century CE, 7.3× and 4.0×) and one from France dating to the Bronze Age (1700 BCE, 3.4×). We assessed the underlying species and identified one medieval remains as a hybrid between two common oak species (Quercus robur and Q. petraea) and the other two remains as Q. robur. We found that diversity at the global genome level had not changed over time. However, exploratory analyses suggested that a reduction of diversity took place at different time periods. Finally, we determined the timing of leaf unfolding for ancient trees for the first time. The study extends the application of ancient wood beyond the classical proxies of dendroclimatology, dendrochronology, dendroarchaeology and dendroecology, thereby enhancing resolution of inferences on the responses of forest ecosystems to past environmental changes, epidemics and silvicultural practices.

Leaf morphometric analysis and potential distribution modelling contribute to taxonomic differentiation in the Quercus microphylla complex

Mexico is a major center of evolutionary radiation for the genus Quercus, with oak species occurring across different habitat types and showing a wide variation in morphology and growth form. Despite representing about 20% of Mexican species, scrub oaks have received little attention and even basic aspects of their taxonomy and geographic distribution remain unresolved. In this study, we analyzed the morphological and climatic niche differentiation of scrub oak populations forming a complex constituted by six named species, Quercus cordifolia, Quercus frutex, Quercus intricata, Quercus microphylla, Quercus repanda, Quercus striatula and a distinct morphotype of Q. striatula identified during field and herbarium work (hereafter named Q. striatula II). Samples were obtained from 35 sites covering the geographic distribution of the complex in northern and central Mexico. Morphological differentiation was analyzed through geometric morphometrics of leaf shape and quantification of trichome traits. Our results indicated the presence of two main morphological groups with geographic concordance. The first was formed by Q. frutex, Q. microphylla, Q. repanda and Q. striatula, distributed in the Trans-Mexican Volcanic Belt, the Sierra Madre Occidental and a little portion of the south of the Mexican Altiplano (MA). The second group consists of Q. cordifola, Q. intricata and Q. striatula II, found in the Sierra Madre Oriental and the MA. Therefore, our evidence supports the distinctness of the Q. striatula II morphotype, indicating the need for a taxonomic revision. Within the two groups, morphological differentiation among taxa varied from very clear to low or inexistent (i.e. Q. microphylla-Q. striatula and Q. cordifolia-Q. striatula II) but niche comparisons revealed significant niche differentiation in all pairwise comparisons, highlighting the relevance of integrative approaches for the taxonomic resolution of complicated groups such as the one studied here.

Nuclear and chloroplast DNA phylogeography reveals high genetic diversity and postglacial range expansion in Quercus mexicana

PREMISE

Phylogeographical studies are fundamental for understanding factors that influence the spatial distribution of genetic lineages within species. Population expansions and contractions, distribution shifts, and climate changes are among the most important factors shaping the genetic compositions of populations.

METHODS

We investigated the phylogeography of an endemic oak, Quercus mexicana (Fagaceae), which has a restricted distribution in northeastern Mexico along the Sierra Madre Oriental and adjacent areas. Nuclear and chloroplast DNA microsatellite markers were used to describe the genetic diversity and structure of 39 populations of Q. mexicana along its entire distribution area. We tested whether population expansion or contraction events influenced the genetic diversity and structure of the species. We also modeled the historical distributional range of Q. mexicana (for the Mid Holocene, the Last Glacial Maximum, and the Last Interglacial) to estimate the extent to which climate fluctuations have impacted the distribution of this oak species.

RESULTS

Our results revealed high genetic diversity and low genetic structure in Q. mexicana populations. Ecological niche models suggested historical fluctuations in the distributional range of Q. mexicana. Historical range changes, gene flow, and physical barriers seem to have played an important role in shaping the phylogeographic structure of Q. mexicana.

CONCLUSIONS

Our study indicates that the genetic structure of Q. mexicana may have been the result of responses of oak trees not only to heterogeneous environments present in the Sierra Madre Oriental and adjacent areas, but also to elevational and latitudinal shifts in response to climate changes in the past.

Angiosperm phylogenetic diversity is lower in Africa than South America

Although originating from a common Gondwanan flora, the diversity and composition of the floras of Africa and South America have greatly diverged since continental breakup of Africa from South America now having much higher plant species richness. However, the phylogenetic diversity of the floras and what this tells us about their evolution remained unexplored. We show that for a given species richness and considering land surface area, topography, and present-day climate, angiosperm phylogenetic diversity in South America is higher than in Africa. This relationship holds regardless of whether all climatically matched areas or only matched areas in tropical climates are considered. Phylogenetic diversity is high relative to species richness in refugial areas in Africa and in northwestern South America, once the gateway for immigration from the north. While species richness is strongly influenced by massive plant radiations in South America, we detect a pervasive influence of historical processes on the phylogenetic diversity of both the South American and African floras.

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.

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.

Cenozoic migration of a desert plant lineage across the North Atlantic

Previous paleobotanical work concluded that Paleogene elements of the sclerophyllous subhumid vegetation of western Eurasia and western North America were endemic to these disjunct regions, suggesting that the southern areas of the Holarctic flora were isolated at that time. Consequently, molecular studies invoked either parallel adaptation to dry climates from related ancestors, or long-distance dispersal in explaining disjunctions between the two regions, dismissing the contemporaneous migration of dry-adapted lineages via land bridges as unlikely. We report Vauquelinia (Rosaceae), currently endemic to western North America, in Cenozoic strata of western Eurasia. Revision of North American fossils previously assigned to Vauquelinia confirmed a single fossil-species of Vauquelinia and one of its close relative Kageneckia. We established taxonomic relationships of fossil-taxa using diagnostic character combinations shared with modern species and constructed a time-calibrated phylogeny. The fossil record suggests that Vauquelinia, currently endemic to arid and subdesert environments, originated under seasonally arid climates in the Eocene of western North America and subsequently crossed the Paleogene North Atlantic land bridge (NALB) to Europe. This pattern is replicated by other sclerophyllous, dry-adapted and warmth-loving plants, suggesting that several of these taxa potentially crossed the North Atlantic via the NALB during Eocene times.

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.

Disentangling leaf structural and material properties in relationship to their anatomical and chemical compositional traits in oaks (Quercus L.)

BACKGROUND AND AIMS

The existence of sclerophyllous plants has been considered an adaptive strategy against different environmental stresses. Given that it literally means 'hard-leaved', it is essential to quantify the leaf mechanical properties to understand sclerophylly. However, the relative importance of each leaf trait for mechanical properties is not yet well established.

METHODS

Genus Quercus is an excellent system to shed light on this because it minimizes phylogenetic variation while having a wide variation in sclerophylly. We measured leaf anatomical traits and cell wall composition, analysing their relationship with leaf mass per area and leaf mechanical properties in a set of 25 oak species.

KEY RESULTS

The upper epidermis outer wall makes a strong and direct contribution to the leaf mechanical strength. Moreover, cellulose plays a crucial role in increasing leaf strength and toughness. The principal component analysis plot based on leaf trait values clearly separates Quercus species into two groups corresponding to evergreen and deciduous species.

CONCLUSIONS

Sclerophyllous Quercus species are tougher and stronger owing to their thicker epidermis outer wall and/or higher cellulose concentration. Furthermore, section Ilex species share common traits, although they occupy different climates. In addition, evergreen species living in mediterranean-type climates share common leaf traits irrespective of their different phylogenetic origin.

Climate-influenced boreotropical survival and rampant introgressions explain the thriving of New World grapes in the north temperate zone

The north temperate region was characterized by a warm climate and a rich thermophilic flora before the Eocene, but early diversifications of the temperate biome under global climate change and biome shift remain uncertain. Moreover, it is becoming clear that hybridization/introgression is an important driving force of speciation in plant diversity. Here, we applied analyses from biogeography and phylogenetic networks to account for both introgression and incomplete lineage sorting based on genomic data from the New World Vitis, a charismatic component of the temperate North American flora with known and suspected gene flow among species. Biogeographic inference and fossil evidence suggest that the grapes were widely distributed from North America to Europe during the Paleocene to the Eocene, followed by widespread extinction and survival of relicts in the tropical New World. During the climate warming in the early Miocene, a Vitis ancestor migrated northward from the refugia with subsequent diversification in the North American region. We found strong evidence for widespread incongruence and reticulate evolution among nuclear genes within both recent and ancient lineages of the New World Vitis. Furthermore, the organellar genomes showed strong conflicts with the inferred species tree from the nuclear genomes. Our phylogenomic analyses provided an important assessment of the wide occurrence of reticulate introgression in the New World Vitis, which potentially represents one of the most important mechanisms for the diversification of Vitis species in temperate North America and even the entire temperate Northern Hemisphere. The scenario we report here may be a common model of temperate diversification of flowering plants adapted to the global climate cooling and fluctuation in the Neogene.

Solutions to fire and shade: resprouting, growing tall and the origin of Eurasian temperate broadleaved forest

Tree species of Eurasian broadleaved forest possess two divergent trait syndromes with contrasting patterns of resource allocation adapted to different selection environments: short-stature basal resprouters that divert resources to a bud bank adapted to frequent and severe disturbances such as fire and herbivory, and tall trees that delay reproduction by investing in rapid height growth to escape shading. Drawing on theory developed in savanna ecosystems, we propose a conceptual framework showing that the possession of contrasting trait syndromes is essential for the persistence of broadleaved trees in an open ecosystem that burns. Consistent with this hypothesis, trees of modern Eurasian broadleaved forest bear a suite of traits that are adaptive to surface and crown-fire regimes. We contend that limited opportunities in grassland restricts recruitment to disturbance-free refugia, and en masse establishment creates a wooded environment where shade limits the growth of light-demanding savanna plants. Rapid height growth, which involves investment in structural support and the switch from a multi-stemmed to a monopodial growth form, is adaptive in this shaded environment. Although clustering reduces surface fuel loads, these establishment nuclei are vulnerable to high-intensity crown fires. The lethal effects of canopy fire are avoided by seasonal leaf shedding, and aerial resprouting enhances rapid post-fire recovery of photosynthetic capacity. While these woody formations satisfy the structural definition of forest, their constituents are clearly derived from savanna. Contrasting trait syndromes thus represent the shift from consumer to resource regulation in savanna ecosystems. Consistent with global trends, the diversification of most contemporary broadleaved taxa coincided with the spread of grasslands, a surge in fire activity and a decline in wooded ecosystems in the late Miocene-Pliocene. Recognition that Eurasian broadleaved forest has savanna origins and persists as an alternative state with adjacent grassy ecosystems has far-reaching management implications in accordance with functional rather than structural criteria. Shade is a severe constraint to the regeneration and growth of both woody and herbaceous growth forms in consumer-regulated ecosystems. However, these ecosystems are highly resilient to disturbance, an essential process that maintains diversity especially among the species-rich herbaceous component that is vulnerable to shading when consumer behaviour is altered.

Applying molecular and genetic methods to trees and their fungal communities

Forests provide invaluable economic, ecological, and social services. At the same time, they are exposed to several threats, such as fragmentation, changing climatic conditions, or increasingly destructive pests and pathogens. Trees, the inherent species of forests, cannot be viewed as isolated organisms. Manifold (micro)organisms are associated with trees playing a pivotal role in forest ecosystems. Of these organisms, fungi may have the greatest impact on the life of trees. A multitude of molecular and genetic methods are now available to investigate tree species and their associated organisms. Due to their smaller genome sizes compared to tree species, whole genomes of different fungi are routinely compared. Such studies have only recently started in forest tree species. Here, we summarize the application of molecular and genetic methods in forest conservation genetics, tree breeding, and association genetics as well as for the investigation of fungal communities and their interrelated ecological functions. These techniques provide valuable insights into the molecular basis of adaptive traits, the impacts of forest management, and changing environmental conditions on tree species and fungal communities and can enhance tree-breeding cycles due to reduced time for field testing. It becomes clear that there are multifaceted interactions among microbial species as well as between these organisms and trees. We demonstrate the versatility of the different approaches based on case studies on trees and fungi. KEY POINTS: • Current knowledge of genetic methods applied to forest trees and associated fungi. • Genomic methods are essential in conservation, breeding, management, and research. • Important role of phytobiomes for trees and their ecosystems.

Genomic data and common garden experiments reveal climate-driven selection on ecophysiological traits in two Mediterranean oaks

Improving our knowledge of how past climate-driven selection has acted on present-day trait population divergence is essential to understand local adaptation processes and improve our predictions of evolutionary trajectories in the face of altered selection pressures resulting from climate change. In this study, we investigated signals of selection on traits related to drought tolerance and growth rates in two Mediterranean oak species (Quercus faginea and Q. lusitanica) with contrasting distribution ranges and climatic niches. We genotyped 182 individuals from 24 natural populations of the two species using restriction-site-associated DNA sequencing and conducted a thorough functional characterization in 1602 seedlings from 21 populations cultivated in common garden experiments under contrasting watering treatments. Our genomic data revealed that both Q. faginea and Q. lusitanica have very weak population genetic structure, probably as a result of high rates of pollen-mediated gene flow among populations and large effective population sizes. In contrast, common garden experiments showed evidence of climate-driven divergent selection among populations on traits related to leaf morphology, physiology and growth in both species. Overall, our study suggests that climate is an important selective factor for Mediterranean oaks and that ecophysiological traits have evolved in drought-prone environments even in a context of very high rates of gene flow among populations.

Phylogenomic insights into the origin and evolutionary history of evergreen broadleaved forests in East Asia under Cenozoic climate change

The evergreen versus deciduous leaf habit is an important functional trait for adaptation of forest trees and has been hypothesized to be related to the evolutionary processes of the component species under paleoclimatic change, and potentially reflected in the dynamic history of evergreen broadleaved forests (EBLFs) in East Asia. However, knowledge about the shift of evergreen versus deciduous leaf with the impact of paleoclimatic change using genomic data remains rare. Here, we focus on the Litsea complex (Lauraceae), a key lineage with dominant species of EBLFs, to gain insights into how evergreen versus deciduous trait shifted, providing insights into the origin and historical dynamics of EBLFs in East Asia under Cenozoic climate change. We reconstructed a robust phylogeny of the Litsea complex using genome-wide single-nucleotide variants (SNVs) with eight clades resolved. Fossil-calibrated analyses, diversification rate shifts, ancestral habit, ecological niche modelling and climate niche reconstruction were employed to estimate its origin and diversification pattern. Taking into account studies on other plant lineages dominating EBLFs of East Asia, it was revealed that the prototype of EBLFs in East Asia probably emerged in the Early Eocene (55-50 million years ago [Ma]), facilitated by the greenhouse warming. As a response to the cooling and drying climate in the Middle to Late Eocene (48-38 Ma), deciduous habits were evolved in the dominant lineages of the EBLFs in East Asia. Up to the Early Miocene (23 Ma), the prevailing of East Asian monsoon increased the extreme seasonal precipitation and accelerated the emergence of evergreen habits of the dominant lineages, and ultimately shaped the vegetation resembling that of today.

Stem hydraulic conductivity and embolism resistance of Quercus species are associated with their climatic niche

The hydraulic traits of a plant species may reflect its climate adaptations. Southwest China is considered as a biodiversity hotpot of the genus Quercus (oak). However, the hydraulic adaptations of Asian oaks to their climate niches remain unclear. Ten common garden-grown oak species with distinct natural distributions in eastern Asia were used to determine their stem xylem embolism resistance (water potential at 50% loss of hydraulic conductivity, P50), stem hydraulic efficiency (vessel anatomy and sapwood specific hydraulic conductivity (Ks)) and leaf anatomical traits. We also compiled four key functional traits: wood density, hydraulic-weighted vessel diameter, Ks and P50 data for 31 oak species from previous literature. We analyzed the relationship between hydraulic traits and climatic factors over the native ranges of 41 oak species. Our results revealed that the 10 Asian oak species, which are mainly distributed in humid subtropical habitats, possessed a stem xylem with low embolism resistance and moderate hydraulic efficiency. The deciduous and evergreen species of the 10 Asian oaks differed in the stem and leaf traits related to hydraulic efficiency. Ks differed significantly between the two phenological groups (deciduous and evergreens) in the 41-oak dataset. No significant difference in P50 between the two groups was found for the 10 Asian oaks or the 41-oak dataset. The oak species that can distribute in arid habitats possessed a stem xylem with high embolism resistance. Ks negatively related to the humidity of the native range of the 10 Asian oaks, but showed no trend when assessing the entire global oak dataset. Our study suggests that stem hydraulic conductivity and embolism resistance in Quercus species are shaped by their climate niche. Our findings assist predictions of oak drought resistance with future climate changes for oak forest management.

Reference-Based RADseq Unravels the Evolutionary History of Polar Species in 'the Crux Lichenologorum' Genus Usnea (Parmeliaceae, Ascomycota)

Nearly 90% of fungal diversity, one of the most speciose branches in the tree of life, remains undescribed. Lichenized fungi as symbiotic associations are still a challenge for species delimitation, and current species diversity is vastly underestimated. The ongoing democratization of Next-Generation Sequencing is turning the tables. Particularly, reference-based RADseq allows for metagenomic filtering of the symbiont sequence and yields robust phylogenomic trees of closely related species. We implemented reference-based RADseq to disentangle the evolution of neuropogonoid lichens, which inhabit harsh environments and belong to Usnea (Parmeliaceae, Ascomycota), one of the most taxonomically intriguing genera within lichenized fungi. Full taxon coverage of neuropogonoid lichens was sampled for the first time, coupled with phenotype characterizations. More than 20,000 loci of 126 specimens were analyzed through concatenated and coalescent-based methods, including time calibrations. Our analysis addressed the major taxonomic discussions over recent decades. Subsequently, two species are newly described, namely U. aymondiana and U. fibriloides, and three species names are resurrected. The late Miocene and Pliocene-Pleistocene boundary is inferred as the timeframe for neuropogonoid lichen diversification. Ultimately, this study helped fill the gap of fungal diversity by setting a solid backbone phylogeny which raises new questions about which factors may trigger complex evolutionary scenarios.

Genomic vulnerability to climate change in Quercus acutissima, a dominant tree species in East Asian deciduous forests

Understanding the evolutionary processes that shape the landscape of genetic variation and influence the response of species to future climate change is critical for biodiversity conservation. Here, we sampled 27 populations across the distribution range of a dominant forest tree, Quercus acutissima, in East Asia, and applied genome-wide analyses to track the evolutionary history and predict the fate of populations under future climate. We found two genetic groups (East and West) in Q. acutissima that diverged during Pliocene. We also found a heterogeneous landscape of genomic variation in this species, which may have been shaped by population demography and linked selections. Using genotype-environment association analyses, we identified climate-associated SNPs in a diverse set of genes and functional categories, indicating a model of polygenic adaptation in Q. acutissima. We further estimated three genetic offset metrics to quantify genomic vulnerability of this species to climate change due to the complex interplay between local adaptation and migration. We found that marginal populations are under higher risk of local extinction because of future climate change, and may not be able to track suitable habitats to maintain the gene-environment relationships observed under the current climate. We also detected higher reverse genetic offsets in northern China, indicating that genetic variation currently present in the whole range of Q. acutissima may not adapt to future climate conditions in this area. Overall, this study illustrates how evolutionary processes have shaped the landscape of genomic variation, and provides a comprehensive genome-wide view of climate maladaptation in Q. acutissima.

The complete chloroplast genome sequence of Quercus kerrii (Fagaceae), and comparative analysis with related species

Quercus kerrii Craib 1911 (section Cyclobalanopsis) is a widespread tree species in the tropical seasonal forests of southwest China and Northern Indo-China areas. In this study, we sequenced, assembled and annotated the complete chloroplast genome of Q. kerrii. The circular genome was 160,743 bp in length and had a GC content of 36.89%. The Q. kerrii chloroplast genome has a typical quadripartite structure, including two inverted repeat regions (length, 25,825 bp; GC content, 42.76%), a large single-copy region (length, 90,196 bp; GC content, 34.74%), and a small single-copy region (length, 18,897 bp; GC content, 30.60%). Genome annotation has indicated that the Q. kerrii chloroplast genome contained 131 genes, including 86 protein-coding genes, 37 tRNA, and eight rRNA. The phylogenetic tree showed that Q. kerrii had a close relationship with Q. schottkyana Rehder & E.H.Wilson 1916.