Oak genome reveals facets of long lifespan

Comparative transcriptome reprogramming in oak galls containing asexual or sexual generations of gall wasps

Oak gall wasps have evolved strategies to manipulate the developmental pathways of their host to induce gall formation. This provides shelter and nutrients for the developing larva. Galls are entirely host tissue; however, the initiation, development, and physical appearance are controlled by the inducer. The underlying molecular mechanisms of gall formation, by which one or a small number of cells are reprogrammed and commit to a novel developmental path, are poorly understood. In this study, we sought a deeper insight into the molecular underpinnings of this process. Oak gall wasps have two generations each year, one sexual, and one asexual. Galls formed by these two generations exhibit a markedly different appearance. We sequenced transcriptomes of both the asexual and sexual generations of Neuroterus quercusbaccarum and Neuroterus numismalis. We then deployed Nanopore sequencing to generate long-read sequences to test the hypothesis that gall wasps introduce DNA insertions to determine gall development. We detected potential genome rearrangements but did not uncover any non-host DNA insertions. Transcriptome analysis revealed that transcriptomes of the sexual generations of distinct species of wasp are more similar than inter-generational comparisons from the same species of wasp. Our results highlight the intricate interplay between the host leaves and gall development, suggesting that season and requirements of the gall structure play a larger role than species in controlling gall development and structure.

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.

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.

Convergent and/or parallel evolution of RNA-binding proteins in angiosperms after polyploidization

Increasing studies suggest that the biased retention of stress-related transcription factors (TFs) after whole-genome duplications (WGDs) could rewire gene transcriptional networks, facilitating plant adaptation to challenging environments. However, the role of posttranscriptional factors (e.g. RNA-binding proteins, RBPs) following WGDs has been largely ignored. Uncovering thousands of RBPs in 21 representative angiosperm species, we integrate genomic, transcriptomic, regulatomic, and paleotemperature datasets to unravel their evolutionary trajectories and roles in adapting to challenging environments. We reveal functional enrichments of RBP genes in stress responses and identify their convergent retention across diverse angiosperms from independent WGDs, coinciding with global cooling periods. Numerous RBP duplicates derived from WGDs are then identified as cold-induced. A significant overlap of 29 orthogroups between WGD-derived and cold-induced RBP genes across diverse angiosperms highlights a correlation between WGD and cold stress. Notably, we unveil an orthogroup (Glycine-rich RNA-binding Proteins 7/8, GRP7/8) and relevant TF duplicates (CCA1/LHY, RVE4/8, CBF2/4, etc.), co-retained in different angiosperms post-WGDs. Finally, we illustrate their roles in rewiring circadian and cold-regulatory networks at both transcriptional and posttranscriptional levels during global cooling. Altogether, we underline the adaptive evolution of RBPs in angiosperms after WGDs during global cooling, improving our understanding of plants surviving periods of environmental turmoil.

Phylogenomic profiles of whole-genome duplications in Poaceae and landscape of differential duplicate retention and losses among major Poaceae lineages

Poaceae members shared a whole-genome duplication called rho. However, little is known about the evolutionary pattern of the rho-derived duplicates among Poaceae lineages and implications in adaptive evolution. Here we present phylogenomic/phylotranscriptomic analyses of 363 grasses covering all 12 subfamilies and report nine previously unknown whole-genome duplications. Furthermore, duplications from a single whole-genome duplication were mapped to multiple nodes on the species phylogeny; a whole-genome duplication was likely shared by woody bamboos with possible gene flow from herbaceous bamboos; and recent paralogues of a tetraploid Oryza are implicated in tolerance of seawater submergence. Moreover, rho duplicates showing differential retention among subfamilies include those with functions in environmental adaptations or morphogenesis, including ACOT for aquatic environments (Oryzoideae), CK2β for cold responses (Pooideae), SPIRAL1 for rapid cell elongation (Bambusoideae), and PAI1 for drought/cold responses (Panicoideae). This study presents a Poaceae whole-genome duplication profile with evidence for multiple evolutionary mechanisms that contribute to gene retention and losses.

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.

The genome assembly of Island Oak (Quercus tomentella), a relictual island tree species

Island oak (Quercus tomentella) is a rare relictual island tree species that exists only on six islands off the coast of California and Mexico, but was once widespread throughout mainland California. Currently, this species is endangered by threats such as non-native plants, grazing animals, and human removal. Efforts for conservation and restoration of island oak currently underway could benefit from information about its range-wide genetic structure and evolutionary history. Here we present a high-quality genome assembly for Q. tomentella, assembled using PacBio HiFi and Omni-C sequencing, developed as part of the California Conservation Genomics Project (CCGP). The resulting assembly has a length of 781 Mb, with a contig N50 of 22.0 Mb and a scaffold N50 of 63.4 Mb. This genome assembly will provide a resource for genomics-informed conservation of this rare oak species. Additionally, this reference genome will be the first one available for a species in Quercus section Protobalanus, a unique oak clade present only in western North America.

Low-frequency somatic mutations are heritable in tropical trees Dicorynia guianensis and Sextonia rubra

Somatic mutations potentially play a role in plant evolution, but common expectations pertaining to plant somatic mutations remain insufficiently tested. Unlike in most animals, the plant germline is assumed to be set aside late in development, leading to the expectation that plants accumulate somatic mutations along growth. Therefore, several predictions were made on the fate of somatic mutations: mutations have generally low frequency in plant tissues; mutations at high frequency have a higher chance of intergenerational transmission; branching topology of the tree dictates mutation distribution; and exposure to UV (ultraviolet) radiation increases mutagenesis. To provide insights into mutation accumulation and transmission in plants, we produced two high-quality reference genomes and a unique dataset of 60 high-coverage whole-genome sequences of two tropical tree species, Dicorynia guianensis (Fabaceae) and Sextonia rubra (Lauraceae). We identified 15,066 de novo somatic mutations in D. guianensis and 3,208 in S. rubra, surprisingly almost all found at low frequency. We demonstrate that 1) low-frequency mutations can be transmitted to the next generation; 2) mutation phylogenies deviate from the branching topology of the tree; and 3) mutation rates and mutation spectra are not demonstrably affected by differences in UV exposure. Altogether, our results suggest far more complex links between plant growth, aging, UV exposure, and mutation rates than commonly thought.

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.

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.

Plant secondary metabolites involved in the stress tolerance of long-lived trees

Ancient trees are natural wonders because of their longevity, having lived for hundreds or thousands of years, and their ability to withstand changing environments and a variety of stresses. These long-lived trees have sophisticated defense mechanisms, such as the production of specialized plant metabolites (SPMs). In this review, we provide an overview of the major biotic and abiotic stresses that long-lived trees often face, as well as an analysis of renowned ancient tree species and their unique protective SPMs against environmental stressors. We also discuss the synthesis and accumulation of defensive SPMs induced by environmental factors and endophytes in these trees. Furthermore, we conducted a comparative genomic analysis of 17 long-lived tree species and discovered significant expansions of SPM biosynthesis gene families in these species. Our comprehensive review reveals the crucial role of SPMs in high resistance in long-lived trees, providing a novel natural resource for plant defense, crop improvement and even the pharmaceutical industry.

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.

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.

Signatures of local adaptation to current and future climate in phenology-related genes in natural populations of Quercus robur

BACKGROUND

Local adaptation is a key evolutionary process that enhances the growth of plants in their native habitat compared to non-native habitats, resulting in patterns of adaptive genetic variation across the entire geographic range of the species. The study of population adaptation to local environments and predicting their response to future climate change is important because of climate change.

RESULTS

Here, we explored the genetic diversity of candidate genes associated with bud burst in pedunculate oak individuals sampled from 6 populations in Poland. Single nucleotide polymorphism (SNP) diversity was assessed in 720 candidate genes using the sequence capture technique, yielding 18,799 SNPs. Using landscape genomic approaches, we identified 8 FST outliers and 781 unique SNPs in 389 genes associated with geography, climate, and phenotypic variables (individual/family spring and autumn phenology, family diameter at breast height (DBH), height, and survival) that are potentially involved in local adaptation. Then, using a nonlinear multivariate model, Gradient Forests, we identified vulnerable areas of the pedunculate oak distribution in Poland that are at risk from climate change.

CONCLUSIONS

The model revealed that pedunculate oak populations in the eastern part of the analyzed geographical region are the most sensitive to climate change. Our results might offer an initial evaluation of a potential management strategy for preserving the genetic diversity of pedunculate oak.

PPGR: a comprehensive perennial plant genomes and regulation database

Perennial woody plants hold vital ecological significance, distinguished by their unique traits. While significant progress has been made in their genomic and functional studies, a major challenge persists: the absence of a comprehensive reference platform for collection, integration and in-depth analysis of the vast amount of data. Here, we present PPGR (Resource for Perennial Plant Genomes and Regulation; https://ngdc.cncb.ac.cn/ppgr/) to address this critical gap, by collecting, integrating, analyzing and visualizing genomic, gene regulation and functional data of perennial plants. PPGR currently includes 60 species, 847 million protein-protein/TF (transcription factor)-target interactions, 9016 transcriptome samples under various environmental conditions and genetic backgrounds. Noteworthy is the focus on genes that regulate wood production, seasonal dormancy, terpene biosynthesis and leaf senescence representing a wealth of information derived from experimental data, literature mining, public databases and genomic predictions. Furthermore, PPGR incorporates a range of multi-omics search and analysis tools to facilitate browsing and application of these extensive datasets. PPGR represents a comprehensive and high-quality resource for perennial plants, substantiated by an illustrative case study that demonstrates its capacity in unraveling gene functions and shedding light on potential regulatory processes.

A Survey of Xylella fastidiosa in the U.S. State of Virginia Reveals Wide Distribution of Both Subspecies fastidiosa and multiplex in Grapevine

Global travel and trade in combination with climate change are expanding the geographic distribution of plant pathogens. The bacterium Xylella fastidiosa is a prime example. Native to the Americas, it has spread to Europe, Asia, and the Middle East. To assess the risk that pathogen introductions pose to crops in newly invaded areas, it is key to survey their diversity, host range, and disease incidence in relation to climatic conditions where they are already present. We performed a survey of X. fastidiosa in grapevine in Virginia using a combination of quantitative PCR, multilocus sequencing, and metagenomics. We also analyzed samples from deciduous trees with leaf scorch symptoms. X. fastidiosa subspecies fastidiosa was identified in grapevines in all regions of the state, even in Northern Virginia, where the temperature was below -9°C for 10 days per year on average in the years preceding sampling. Unexpectedly, we also found for the first time grapevine samples infected with X. fastidiosa subspecies multiplex (Xfm). The Xfm lineage found in grapevines had been previously isolated from blueberries in the Southeastern United States and was distinct from that found in deciduous trees in Virginia. The obtained results will be important for risk assessment of X. fastidiosa introductions in other parts of the world.

Genome-Wide Identification of the CBF Gene Family and ICE Transcription Factors in Walnuts and Expression Profiles under Cold Conditions

Cold stress impacts woody tree growth and perennial production, especially when the temperature rapidly changes in late spring. To address this issue, we conducted the genome-wide identification of two important transcription factors (TFs), CBF (C-repeat binding factors) and ICE (inducers of CBF expression), in three walnut (Juglans) genomes. Although the CBF and ICE gene families have been identified in many crops, very little systematic analysis of these genes has been carried out in J. regia and J. sigillata. In this study, we identified a total of 16 CBF and 12 ICE genes in three Juglans genomes using bioinformatics analysis. Both CBF and ICE had conserved domains, motifs, and gene structures, which suggests that these two TFs were evolutionarily conserved. Most ICE genes are located at both ends of the chromosomes. The promoter cis-regulatory elements of CBF and ICE genes are largely involved in light and phytohormone responses. Based on 36 RNA sequencing of leaves from four walnut cultivars ('Zijing', 'Lvling', 'Hongren', and 'Liao1') under three temperature conditions (8 °C, 22 °C, and 5 °C) conditions in late spring, we found that the ICE genes were expressed more highly than CBFs. Both CBF and ICE proteins interacted with cold-related proteins, and many putative miRNAs had interactions with these two TFs. These results determined that CBF1 and ICE1 play important roles in the tolerance of walnut leaves to rapid temperature changes. Our results provide a useful resource on the function of the CBF and ICE genes related to cold tolerance in walnuts.

Chromosome-level genome assembly of Quercus variabilis provides insights into the molecular mechanism of cork thickness

Quercus variabilis is a deciduous woody species with high ecological and economic value, and is a major source of cork in East Asia. Cork from thick softwood sheets have higher commercial value than those from thin sheets. It is extremely difficult to genetically improve Q. variabilis to produce high quality softwood due to the lack of genomic information. Here, we present a high-quality chromosomal genome assembly for Q. variabilis with length of 791,89 Mb and 54,606 predicted genes. Comparative analysis of protein sequences of Q. variabilis with 11 other species revealed that specific and expanded gene families were significantly enriched in the "fatty acid biosynthesis" pathway in Q. variabilis, which may contribute to the formation of its unique cork. Based on weighted correlation network analysis of time-course (i.e., five important developmental ages) gene expression data in thick-cork versus thin-cork genotypes of Q. variabilis, we identified one co-expression gene module associated with the thick-cork trait. Within this co-expression gene module, 10 hub genes were associated with suberin biosynthesis. Furthermore, we identified a total of 198 suberin biosynthesis-related new candidate genes that were up-regulated in trees with a thick cork layer relative to those with a thin cork layer. Also, we found that some genes related to cell expansion and cell division were highly expressed in trees with a thick cork layer. Collectively, our results revealed that two metabolic pathways (i.e., suberin biosynthesis, fatty acid biosynthesis), along with other genes involved in cell expansion, cell division, and transcriptional regulation, were associated with the thick-cork trait in Q. variabilis, providing insights into the molecular basis of cork development and knowledge for informing genetic improvement of cork thickness in Q. variabilis and closely related species.

Identification of Functional Brassinosteroid Receptor Genes in Oaks and Functional Analysis of QmBRI1

Brassinosteroids (BRs) play important regulatory roles in plant growth and development, with functional BR receptors being crucial for BR recognition or signaling. Although functional BR receptors have been extensively studied in herbaceous plants, they remain largely under-studied in forest tree species. In this study, nine BR receptors were identified in three representative oak species, of which BRI1s and BRL1s were functional BR receptors. Dispersed duplications were a driving force for oak BR receptor expansion, among which the Brassinosteroid-Insensitive-1 (BRI1)-type genes diverged evolutionarily from most rosids. In oak BRI1s, we identified that methionine in the conserved Asn-Gly-Ser-Met (NGSM) motif was replaced by isoleucine and that the amino acid mutation occurred after the divergence of Quercus and Fagus. Compared with QmBRL1, QmBRI1 was relatively highly expressed during BR-induced xylem differentiation and in young leaves, shoots, and the phloem and xylem of young stems of Quercus mongolica. Based on Arabidopsis complementation experiments, we proved the important role of QmBRI1 in oak growth and development, especially in vascular patterning and xylem differentiation. These findings serve as an important supplement to the findings of the structural, functional and evolutionary studies on functional BR receptors in woody plants and provide the first example of natural mutation occurring in the conserved BR-binding region (NGSM motif) of angiosperm BRI1s.

Telomere-to-telomere and haplotype-resolved genome assembly of the Chinese cork oak (Quercus variabilis)

The Quercus variabilis, a deciduous broadleaved tree species, holds significant ecological and economical value. While a chromosome-level genome for this species has been made available, it remains riddled with unanchored sequences and gaps. In this study, we present a nearly complete comprehensive telomere-to-telomere (T2T) and haplotype-resolved reference genome for Q. variabilis. This was achieved through the integration of ONT ultra-long reads, PacBio HiFi long reads, and Hi-C data. The resultant two haplotype genomes measure 789 Mb and 768 Mb in length, with a contig N50 of 65 Mb and 56 Mb, and were anchored to 12 allelic chromosomes. Within this T2T haplotype-resolved assembly, we predicted 36,830 and 36,370 protein-coding genes, with 95.9% and 96.0% functional annotation for each haplotype genome. The availability of the T2T and haplotype-resolved reference genome lays a solid foundation, not only for illustrating genome structure and functional genomics studies but also to inform and facilitate genetic breeding and improvement of cultivated Quercus species.

A haplotype-resolved chromosome-scale genome for Quercus rubra L. provides insights into the genetics of adaptive traits for red oak species

Northern red oak (Quercus rubra L.) is an ecologically and economically important forest tree native to North America. We present a chromosome-scale genome of Q. rubra generated by the combination of PacBio sequences and chromatin conformation capture (Hi-C) scaffolding. This is the first reference genome from the red oak clade (section Lobatae). The Q. rubra assembly spans 739 Mb with 95.27% of the genome in 12 chromosomes and 33,333 protein-coding genes. Comparisons to the genomes of Quercus lobata and Quercus mongolica revealed high collinearity, with intrachromosomal structural variants present. Orthologous gene family analysis with other tree species revealed that gene families associated with defense response were expanding and contracting simultaneously across the Q. rubra genome. Quercus rubra had the most CC-NBS-LRR and TIR-NBS-LRR resistance genes out of the 9 species analyzed. Terpene synthase gene family comparisons further reveal tandem gene duplications in TPS-b subfamily, similar to Quercus robur. Phylogenetic analysis also identified 4 subfamilies of the IGT/LAZY gene family in Q. rubra important for plant structure. Single major QTL regions were identified for vegetative bud break and marcescence, which contain candidate genes for further research, including a putative ortholog of the circadian clock constituent cryptochrome (CRY2) and 8 tandemly duplicated genes for serine protease inhibitors, respectively. Genome-environment associations across natural populations identified candidate abiotic stress tolerance genes and predicted performance in a common garden. This high-quality red oak genome represents an essential resource to the oak genomic community, which will expedite comparative genomics and biological studies in Quercus species.

A first draft genome of holm oak (Quercus ilex subsp. ballota), the most representative species of the Mediterranean forest and the Spanish agrosylvopastoral ecosystem "dehesa"

The holm oak (Quercus ilex subsp. ballota) is the most representative species of the Mediterranean Basin and the agrosylvopastoral Spanish "dehesa" ecosystem. Being part of our life, culture, and subsistence since ancient times, it has significant environmental and economic importance. More recently, there has been a renewed interest in using the Q. ilex acorn as a functional food due to its nutritional and nutraceutical properties. However, the holm oak and its related ecosystems are threatened by different factors, with oak decline syndrome and climate change being the most worrying in the short and medium term. Breeding programs informed by the selection of elite genotypes seem to be the most plausible biotechnological solution to rescue populations under threat. To achieve this and other downstream analyses, we need a high-quality and well-annotated Q. ilex reference genome. Here, we introduce the first draft genome assembly of Q. ilex using long-read sequencing (PacBio). The assembled nuclear haploid genome had 530 contigs totaling 842.2 Mbp (N50 = 3.3 Mbp), of which 448.7 Mb (53%) were repetitive sequences. We annotated 39,443 protein-coding genes of which 94.80% were complete and single-copy genes. Phylogenetic analyses showed no evidence of a recent whole-genome duplication, and high synteny of the 12 chromosomes between Q. ilex and Quercus lobata and between Q. ilex and Quercus robur. The chloroplast genome size was 142.3 Kbp with 149 protein-coding genes successfully annotated. This first draft should allow for the validation of omics data as well as the identification and functional annotation of genes related to phenotypes of interest such as those associated with resilience against oak decline syndrome and climate change and higher acorn productivity and nutraceutical value.

Re-evaluating the driving force behind mutations

Experiments on tropical trees suggest that new mutations in plants are driven by age rather than number of cell divisions during growth.

Genome and whole-genome resequencing of Cinnamomum camphora elucidate its dominance in subtropical urban landscapes

BACKGROUND

Lauraceae is well known for its significant phylogenetic position as well as important economic and ornamental value; however, most evergreen species in Lauraceae are restricted to tropical regions. In contrast, camphor tree (Cinnamomum camphora) is the most dominant evergreen broadleaved tree in subtropical urban landscapes.

RESULTS

Here, we present a high-quality reference genome of C. camphora and conduct comparative genomics between C. camphora and C. kanehirae. Our findings demonstrated the significance of key genes in circadian rhythms and phenylpropanoid metabolism in enhancing cold response, and terpene synthases (TPSs) improved defence response with tandem duplication and gene cluster formation in C. camphora. Additionally, the first comprehensive catalogue of C. camphora based on whole-genome resequencing of 75 accessions was constructed, which confirmed the crucial roles of the above pathways and revealed candidate genes under selection in more popular C. camphora, and indicated that enhancing environmental adaptation is the primary force driving C. camphora breeding and dominance.

CONCLUSIONS

These results decipher the dominance of C. camphora in subtropical urban landscapes and provide abundant genomic resources for enlarging the application scopes of evergreen broadleaved trees.

The Reemergence of Phycopathology: When Algal Biology Meets Ecology and Biosecurity

Viruses, bacteria, and eukaryotic symbionts interact with algae in a variety of ways to cause disease complexes, often shaping marine and freshwater ecosystems. The advent of phyconomy (a.k.a. seaweed agronomy) represents a need for a greater understanding of algal disease interactions, where underestimated cryptic diversity and lack of phycopathological basis are prospective constraints for algal domestication. Here, we highlight the limited yet increasing knowledge of algal pathogen biodiversity and the ecological interaction with their algal hosts. Finally, we discuss how ecology and cultivation experience contribute to and reinforce aquaculture practice, with the potential to reshape biosecurity policies of seaweed cultivation worldwide.

The genome of the Wollemi pine, a critically endangered "living fossil" unchanged since the Cretaceous, reveals extensive ancient transposon activity

We present the genome of the living fossil, Wollemia nobilis, a southern hemisphere conifer morphologically unchanged since the Cretaceous. Presumed extinct until rediscovery in 1994, the Wollemi pine is critically endangered with less than 60 wild adults threatened by intensifying bushfires in the Blue Mountains of Australia. The 12 Gb genome is among the most contiguous large plant genomes assembled, with extremely low heterozygosity and unusual abundance of DNA transposons. Reduced representation and genome re-sequencing of individuals confirms a relictual population since the last major glacial/drying period in Australia, 120 ky BP. Small RNA and methylome sequencing reveal conservation of ancient silencing mechanisms despite the presence of thousands of active and abundant transposons, including some transferred horizontally to conifers from arthropods in the Jurassic. A retrotransposon burst 8-6 my BP coincided with population decline, possibly as an adaptation enhancing epigenetic diversity. Wollemia, like other conifers, is susceptible to Phytophthora, and a suite of defense genes, similar to those in loblolly pine, are targeted for silencing by sRNAs in leaves. The genome provides insight into the earliest seed plants, while enabling conservation efforts.

Low baseline intraspecific variation in leaf pressure-volume traits: Biophysical basis and implications for spectroscopic sensing

Intra-specific trait variation (ITV) plays a role in processes at a wide range of scales from organs to ecosystems across climate gradients. Yet, ITV remains rarely quantified for many ecophysiological traits typically assessed for species means, such as pressure volume (PV) curve parameters including osmotic potential at full turgor and modulus of elasticity, which are important in plant water relations. We defined a baseline "reference ITV" (ITVref ) as the variation among fully exposed, mature sun leaves of replicate individuals of a given species grown in similar, well-watered conditions, representing the conservative sampling design commonly used for species-level ecophysiological traits. We hypothesized that PV parameters would show low ITVref relative to other leaf morphological traits, and that their intraspecific relationships would be similar to those previously established across species and proposed to arise from biophysical constraints. In a database of novel and published PV curves and additional leaf structural traits for 50 diverse species, we found low ITVref for PV parameters relative to other morphological traits, and strong intraspecific relationships among PV traits. Simulation modeling showed that conservative ITVref enables the use of species-mean PV parameters for scaling up from spectroscopic measurements of leaf water content to enable sensing of leaf water potential.

Genetic and Epigenetic Mechanisms of Longevity in Forest Trees

Trees are unique in terms of development, sustainability and longevity. Some species have a record lifespan in the living world, reaching several millennia. The aim of this review is to summarize the available data on the genetic and epigenetic mechanisms of longevity in forest trees. In this review, we have focused on the genetic aspects of longevity of a few well-studied forest tree species, such as Quercus robur, Ginkgo biloba, Ficus benghalensis and F. religiosa, Populus, Welwitschia and Dracaena, as well as on interspecific genetic traits associated with plant longevity. A key trait associated with plant longevity is the enhanced immune defense, with the increase in gene families such as RLK, RLP and NLR in Quercus robur, the expansion of the CC-NBS-LRR disease resistance families in Ficus species and the steady expression of R-genes in Ginkgo biloba. A high copy number ratio of the PARP1 family genes involved in DNA repair and defense response was found in Pseudotsuga menziesii, Pinus sylvestris and Malus domestica. An increase in the number of copies of the epigenetic regulators BRU1/TSK/MGO3 (maintenance of meristems and genome integrity) and SDE3 (antiviral protection) was also found in long-lived trees. CHG methylation gradually declines in the DAL 1 gene in Pinus tabuliformis, a conservative age biomarker in conifers, as the age increases. It was shown in Larix kaempferi that grafting, cutting and pruning change the expression of age-related genes and rejuvenate plants. Thus, the main genetic and epigenetic mechanisms of longevity in forest trees were considered, among which there are both general and individual processes.

Pervasive tandem duplications and convergent evolution shape coral genomes

BACKGROUND

Over the last decade, several coral genomes have been sequenced allowing a better understanding of these symbiotic organisms threatened by climate change. Scleractinian corals are reef builders and are central to coral reef ecosystems, providing habitat to a great diversity of species.

RESULTS

In the frame of the Tara Pacific expedition, we assemble two coral genomes, Porites lobata and Pocillopora cf. effusa, with vastly improved contiguity that allows us to study the functional organization of these genomes. We annotate their gene catalog and report a relatively higher gene number than that found in other public coral genome sequences, 43,000 and 32,000 genes, respectively. This finding is explained by a high number of tandemly duplicated genes, accounting for almost a third of the predicted genes. We show that these duplicated genes originate from multiple and distinct duplication events throughout the coral lineage. They contribute to the amplification of gene families, mostly related to the immune system and disease resistance, which we suggest to be functionally linked to coral host resilience.

CONCLUSIONS

At large, we show the importance of duplicated genes to inform the biology of reef-building corals and provide novel avenues to understand and screen for differences in stress resilience.

Chromosome-scale genome assembly and insights into the metabolome and gene regulation of leaf color transition in an important oak species, Quercus dentata

Quercus dentata Thunb., a dominant forest tree species in northern China, has significant ecological and ornamental value due to its adaptability and beautiful autumn coloration, with color changes from green to yellow into red resulting from the autumnal shifts in leaf pigmentation. However, the key genes and molecular regulatory mechanisms for leaf color transition remain to be investigated. First, we presented a high-quality chromosome-scale assembly for Q. dentata. This 893.54 Mb sized genome (contig N50 = 4.21 Mb, scaffold N50 = 75.55 Mb; 2n = 24) harbors 31 584 protein-coding genes. Second, our metabolome analyses uncovered pelargonidin-3-O-glucoside, cyanidin-3-O-arabinoside, and cyanidin-3-O-glucoside as the main pigments involved in leaf color transition. Third, gene co-expression further identified the MYB-bHLH-WD40 (MBW) transcription activation complex as central to anthocyanin biosynthesis regulation. Notably, transcription factor (TF) QdNAC (QD08G038820) was highly co-expressed with this MBW complex and may regulate anthocyanin accumulation and chlorophyll degradation during leaf senescence through direct interaction with another TF, QdMYB (QD01G020890), as revealed by our further protein-protein and DNA-protein interaction assays. Our high-quality genome assembly, metabolome, and transcriptome resources further enrich Quercus genomics and will facilitate upcoming exploration of ornamental values and environmental adaptability in this important genus.

Multifeature analysis of age-related microbiome structures reveals defense mechanisms of Populus tomentosa trees

Root microbiota composition shifts during the development of most annual plants. Although some perennial plants can live for centuries, the host-microbiome partnerships and interaction mechanisms underlying their longevity remain unclear. To address this gap, we investigated age-related changes in the root metabolites, transcriptomes, and microbiome compositions of 1- to 35-yr-old Populus tomentosa trees. Ten co-response clusters were obtained according to their accumulation patterns, and members of each cluster displayed a uniform and clear pattern of abundance. Multi-omics network analysis demonstrated that the increased abundance of Actinobacteria with tree age was strongly associated with the flavonoid biosynthesis. Using genetic approaches, we demonstrate that the flavonoid biosynthesis regulator gene Transparent Testa 8 is associated with the recruitment of flavonoid-associated Actinobacteria. Further inoculation experiments of Actinobacteria isolates indicated that their colonization could significantly improve the host's phenotype. Site-directed mutagenesis revealed that the hyBl gene cluster, involved in biosynthesis of an aminocyclitol hygromycin B analog in Streptomyces isolate bj1, is associated with disease suppression. We hypothesize that interactions between perennial plants and soil microorganisms lead to gradual enrichment of a subset of microorganisms that may harbor a wealth of currently unknown functional traits.

A chromosome-scale genome assembly of Castanopsis hystrix provides new insights into the evolution and adaptation of Fagaceae species

Fagaceae species dominate forests and shrublands throughout the Northern Hemisphere, and have been used as models to investigate the processes and mechanisms of adaptation and speciation. Compared with the well-studied genus Quercus, genomic data is limited for the tropical-subtropical genus Castanopsis. Castanopsis hystrix is an ecologically and economically valuable species with a wide distribution in the evergreen broad-leaved forests of tropical-subtropical Asia. Here, we present a high-quality chromosome-scale reference genome of C. hystrix, obtained using a combination of Illumina and PacBio HiFi reads with Hi-C technology. The assembled genome size is 882.6 Mb with a contig N50 of 40.9 Mb and a BUSCO estimate of 99.5%, which are higher than those of recently published Fagaceae species. Genome annotation identified 37,750 protein-coding genes, of which 97.91% were functionally annotated. Repeat sequences constituted 50.95% of the genome and LTRs were the most abundant repetitive elements. Comparative genomic analysis revealed high genome synteny between C. hystrix and other Fagaceae species, despite the long divergence time between them. Considerable gene family expansion and contraction were detected in Castanopsis species. These expanded genes were involved in multiple important biological processes and molecular functions, which may have contributed to the adaptation of the genus to a tropical-subtropical climate. In summary, the genome assembly of C. hystrix provides important genomic resources for Fagaceae genomic research communities, and improves understanding of the adaptation and evolution of forest trees.

Modelling somatic mutation accumulation and expansion in a long-lived tree with hierarchical modular architecture

In a long-lived organism with a modular architecture, such as trees, somatic mutations accumulate throughout the long lifespan and result in genetic mosaicism in each module within the same individual. In recent years, next-generation sequencing technology has provided a snapshot of such intra-organismal genetic variability. However, the dynamic processes underlying the accumulation and expansion of somatic mutations during the growth remain poorly understood. In this study, we constructed a model to describe these processes in a form that can be applied to a real tree. Given that the proliferation dynamics of meristematic cells vary across plant species, multiple possible processes for elongation and branching were comprehensively expressed in our model. Using published data from a poplar tree, we compared the prediction of the models with the observation and explained the cell lineage dynamics underlying somatic mutations accumulation that were not evident from the snapshot of the sequenced data. We showed that the somatic genetic drift during growth increases inter-meristem mosaicism, resulting in genetically distinct branches and less integrity within an individual tree. We also showed that the somatic genetic drift during branching leads to the mutation accumulation pattern that does not reflect the tree topology. Our modelling framework can help interpret and provide further insights into the empirical findings of genetic mosaicism in long-lived trees.

Dual domestications and origin of traits in grapevine evolution

We elucidate grapevine evolution and domestication histories with 3525 cultivated and wild accessions worldwide. In the Pleistocene, harsh climate drove the separation of wild grape ecotypes caused by continuous habitat fragmentation. Then, domestication occurred concurrently about 11,000 years ago in Western Asia and the Caucasus to yield table and wine grapevines. The Western Asia domesticates dispersed into Europe with early farmers, introgressed with ancient wild western ecotypes, and subsequently diversified along human migration trails into muscat and unique western wine grape ancestries by the late Neolithic. Analyses of domestication traits also reveal new insights into selection for berry palatability, hermaphroditism, muscat flavor, and berry skin color. These data demonstrate the role of the grapevines in the early inception of agriculture across Eurasia.

Traits variation of acorns and cupules during maturation process in Quercus variabilis and Quercus aliena

Quercus variabilis and Quercus aliena are two native tree species in China, which have similar habitats, and their regeneration mainly depends on acorn dispersal. This study analyzed the contents of water, soluble sugar, starch, soluble protein, and total phenolics in acorns and cupules during the whole development process to explore the difference between species. Thereinto, starch and total phenol occupied the dominant roles as their high contents. The acorn starch contents increased sharply during development in both species, but the contents in Q. variabilis were almost twice those of Q. aliena when mature. Similarly, high expression levels of starch synthase, soluble starch synthase 2 (SSS2) were also found in the acorns of Q. variabilis. The total phenol contents in Q. variabilis acorns were high at the early stages, and decreased sharply to similar contents in Q. aliena when mature. Additionally, the cupules in Q. variabilis had high contents of total phenols during the whole development period. Similar trends were also found in the expression patterns of UGT84A13 and SDH. The high total phenols in acorns and cupules of Q. variabilis probably protect the acorns from Mechoris ursulus, as only Q. aliena suffered a severe pest infestation in the early development stages. This study not only clarifies the interspecific difference between storage and defense substances during the development process in acorns and cupules, but also deepens understanding the specialized mechanisms of plant-pest/animal interactions in Quercus.

Oak stands along an elevation gradient have different molecular strategies for regulating bud phenology

BACKGROUND

Global warming raises serious concerns about the persistence of species and populations locally adapted to their environment, simply because of the shift it produces in their adaptive landscape. For instance, the phenological cycle of tree species may be strongly affected by higher winter temperatures and late frost in spring. Given the variety of ecosystem services they provide, the question of forest tree adaptation has received increasing attention in the scientific community and catalyzed research efforts in ecology, evolutionary biology and functional genomics to study their adaptive capacity to respond to such perturbations.

RESULTS

In the present study, we used an elevation gradient in the Pyrenees Mountains to explore the gene expression network underlying dormancy regulation in natural populations of sessile oak stands sampled along an elevation cline and potentially adapted to different climatic conditions mainly driven by temperature. By performing analyses of gene expression in terminal buds we identified genes displaying significant dormancy, elevation or dormancy-by-elevation interaction effects. Our Results highlighted that low- and high-altitude populations have evolved different molecular strategies for minimizing late frost damage and maximizing the growth period, thereby increasing potentially their respective fitness in these contrasting environmental conditions. More particularly, population from high elevation overexpressed genes involved in the inhibition of cell elongation and delaying flowering time while genes involved in cell division and flowering, enabling buds to flush earlier were identified in population from low elevation.

CONCLUSION

Our study made it possible to identify key dormancy-by-elevation responsive genes revealing that the stands analyzed in this study have evolved distinct molecular strategies to adapt their bud phenology in response to temperature.

Translational research from models to crops: comparative genomics for plant breeding

The concept of translational research, which originated in the medical field in the 1980s, consists in improving the efficient transfer of research results obtained in a species (which can be considered as a model or pivot) to all the species for which these results are of interest for its improvement in Agriculture. In this context, comparative genomics is an important tool for translational research, effectively identifying genes controlling common functions between species. Editing and phenotyping tools must thus allow the functional validation of the gene conserved within the species for which the knowledge has been extrapolated, that is to say transferred, and the identification of the best alleles and associated genotypes for exploitation in current breeding programs.

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

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

How ancient forest fragmentation and riparian connectivity generate high levels of genetic diversity in a microendemic Malagasy tree

Understanding landscape changes is central to predicting evolutionary trajectories and defining conservation practices. While human-driven deforestation is intense throughout Madagascar, exceptions in areas such as the Loky-Manambato region (north) raise questions regarding the causes and age of forest fragmentation. The Loky-Manambato region also harbours a rich and endemic flora, whose evolutionary origin remains poorly understood. We assessed the genetic diversity of an endangered microendemic Malagasy olive species (Noronhia spinifolia Hong-Wa) to better understand the vegetation dynamics in the Loky-Manambato region and its influence on past evolutionary processes. We characterized 72 individuals sampled across eight forests through nuclear and mitochondrial restriction-associated DNA sequencing data and chloroplast microsatellites. Combined population and landscape genetics analyses indicate that N. spinifolia diversity is largely explained by the current forest cover, highlighting a long-standing habitat mosaic in the region. This sustains a major and long-term role of riparian corridors in maintaining connectivity across these antique mosaic habitats, calling for the study of organismal interactions that promote gene flow.

High level of somatic mutations detected in a diploid banana wild relative Musa basjoo

Plants are thought to lack an early segregating germline and often retain both asexual and sexual reproduction, both of which may allow somatic mutations to enter the gametes or clonal progeny, and thereby impact plant evolution. It is yet unclear how often these somatic mutations occur during plant development and what proportion is transmitted to their sexual or cloned offspring. Asexual "seedless" propagation has contributed greatly to the breeding in many fruit crops, such as citrus, grapes and bananas. Whether plants in these lineages experience substantial somatic mutation accumulation is unknown. To estimate the somatic mutation accumulation and inheritance among a clonal population of plant, here we assess somatic mutation accumulation in Musa basjoo, a diploid banana wild relative, using 30 whole-genome resequenced samples collected from five structures, including leaves, sheaths, panicle, roots and underground rhizome connecting three clonal individuals. We observed 18.5 high proportion de novo somatic mutations on average between each two adjacent clonal suckers, equivalent to ~ 2.48 × 10^-8 per site per asexual generation, higher than the per site per sexual generation rates (< 1 × 10^-8) reported in Arabidopsis and peach. Interestingly, most of these inter-ramet somatic mutations were shared simultaneously in different tissues of the same individual with a high level of variant allele fractions, suggesting that these somatic mutations arise early in ramet development and that each individual may develop only from a few apical stem cells. These results thus suggest substantial mutation accumulation in a wild relative of banana. Our work reveals the significance of somatic mutation in Musa basjoo genetics variations and contribute to the trait improvement breeding of bananas and other asexual clonal crops.

Molecular plasticity to soil water deficit differs between sessile oak (Quercus Petraea (Matt.) Liebl.) high- and low-water use efficiency genotypes

Water use efficiency (WUE) is an important adaptive trait for soil water deficit. The molecular and physiological bases of WUE regulation in crops have been studied in detail in the context of plant breeding. Knowledge for most forest tree species lags behind, despite the need to identify populations or genotypes able to cope with the longer, more intense drought periods likely to result from climate warming. We aimed to bridge this gap in knowledge for sessile oak (Quercus petraea (Matt.) Liebl.), one of the most ecologically and economically important tree species in Europe, using a factorial design including trees with contrasted phenotypic values (low and high WUE) and two watering regimes (control and drought). By monitoring the ecophysiological response, we first qualified genotypes for their WUE (by using instantaneous and long-term measures). We then performed RNA-seq to quantify gene expression for the three most extreme genotypes exposed to the two watering regimes. By analyzing the interaction term, we were able to capture the molecular strategy of each group of plants for coping with drought. We identified putative candidate genes potentially involved in the regulation of transpiration rate in high-WUE phenotypes. Regardless of water availability, trees from the high-WUE phenotypic class overexpressed genes associated with drought responses, and in the control of stomatal density and distribution, and displayed a downregulation of genes associated with early stomatal closure and high transpiration rate. Fine physiological screening of sessile oaks with contrasting WUE, and their molecular characterization (i) highlighted subtle differences in transcription between low- and high-WUE genotypes, identifying key molecular players in the genetic control of this trait and (ii) revealed the genes underlying the molecular strategy that evolved in each group to potentially cope with water deficit, providing new insight into the within-species diversity in drought adaptation strategies.

Strategies to mitigate shifts in red oak (Quercus sect. Lobatae) distribution under a changing climate

Red oaks (Quercus sect. Lobatae) are a taxonomic group of hardwood trees, which occur in swamp forests, subtropical chaparral and savannahs from Columbia to Canada. They cover a wide range of ecological niches, and many species are thought to be able to cope with current trends in climate change. Genus Quercus encompasses ca. 500 species, of which ca. 80 make up sect. Lobatae. Species diversity is greatest within the southeastern USA and within the northern and eastern regions of Mexico. This review discusses the weak reproductive barriers between species of red oaks and the effects this has on speciation and niche range. Distribution and diversity have been shaped by drought adaptations common to the species of sect. Lobatae, which enable them to fill various xeric niches across the continent. Drought adaptive traits of this taxonomic group include deciduousness, deep tap roots, ring-porous xylem, regenerative stump sprouting, greater leaf thickness and smaller stomata. The complex interplay between these anatomical and morphological traits has given red oaks features of drought tolerance and avoidance. Here, we discuss physiological and genetic components of these adaptations to address how many species of sect. Lobatae reside within xeric sites and/or sustain normal metabolic function during drought. Although extensive drought adaptation appears to give sect. Lobatae a resilience to climate change, aging tree stands, oak life history traits and the current genetic structures place many red oak species at risk. Furthermore, oak decline, a complex interaction between abiotic and biotic agents, has severe effects on red oaks and is likely to accelerate species decline and fragmentation. We suggest that assisted migration can be used to avoid species fragmentation and increase climate change resilience of sect. Lobatae.

Gene expression and genetic divergence in oak species highlight adaptive genes to soil water constraints

Drought and waterlogging impede tree growth and may even lead to tree death. Oaks, an emblematic group of tree species, have evolved a range of adaptations to cope with these constraints. The two most widely distributed European species, pedunculate (PO; Quercus robur L.) and sessile oak (SO; Quercus petraea Matt. Lieb), have overlapping ranges, but their respective distribution are highly constrained by local soil conditions. These contrasting ecological preferences between two closely related and frequently hybridizing species constitute a powerful model to explore the functional bases of the adaptive responses in oak. We exposed oak seedlings to waterlogging and drought, conditions typically encountered by the two species in their respective habitats, and studied changes in gene expression in roots using RNA-seq. We identified genes that change in expression between treatments differentially depending on species. These "species × environment"-responsive genes revealed adaptive molecular strategies involving adventitious and lateral root formation, aerenchyma formation in PO, and osmoregulation and ABA regulation in SO. With this experimental design, we also identified genes with different expression between species independently of water conditions imposed. Surprisingly, this category included genes with functions consistent with a role in intrinsic reproductive barriers. Finally, we compared our findings with those for a genome scan of species divergence and found that the expressional candidate genes included numerous highly differentiated genetic markers between the two species. By combining transcriptomic analysis, gene annotation, pathway analyses, as well as genome scan for genetic differentiation among species, we were able to highlight loci likely involved in adaptation of the two species to their respective ecological niches.

Chromosome-level genome assembly provides new insights into Japanese chestnut (Castanea crenata) genomes

Japanese chestnut (Castanea crenata Sieb. et Zucc) is an economically and ecologically important chestnut species in East Asia. Here, we presented a high-quality chromosome-level reference genome of the Japanese chestnut cultivar 'Tsukuba' by combining Nanopore long reads and Hi-C sequencing. The final assembly has a size of 718.30 Mb and consists of 12 pseudochromosomes ranging from 41.03 to 92.03 Mb, with a BUSCO complete gene percentage of 97.6%. A total of 421.37 Mb repetitive sequences and 46,744 gene models encoding 46,463 proteins were predicted in the genome. Genome evolution analysis showed that Japanese chestnut is closely related to Chinese chestnut and these species shared a common ancestor ~6.5 million years ago. This high-quality Japanese chestnut genome represents an important resource for the chestnut genomics community and will improve our understanding of chestnut biology and evolution.

Genome sequencing and comparative analysis of Ficus benghalensis and Ficus religiosa species reveal evolutionary mechanisms of longevity

Ficus benghalensis and Ficus religiosa are large woody trees well known for their long lifespan, ecological and traditional significance, and medicinal properties. To understand the genomic and evolutionary aspects of these characteristics, the whole genomes of these Ficus species were sequenced using 10x Genomics linked reads and Oxford Nanopore long reads. The draft genomes of F. benghalensis and F. religiosa comprised of 392.89 Mbp and 332.97 Mbp, respectively. We established the genome-wide phylogenetic positions of the two Ficus species with respect to 50 other Angiosperm species. Comparative evolutionary analyses with other phylogenetically closer Eudicot species revealed adaptive evolution in genes involved in key cellular mechanisms associated with prolonged survival including phytohormones signaling, senescence, disease resistance, and abiotic stress tolerance, which provide genomic insights into the mechanisms conferring longevity and suggest that longevity is a multifaceted phenomenon. This study also provides clues on the existence of CAM pathway in these Ficus species.

•

First whole genome assemblies of Ficus benghalensis and Ficus religiosa trees

•

Genome-wide phylogeny with 50 other Angiosperm species

•

Evolution of genes in phytohormone signaling, senescence, and stress tolerance

•

Genomic insights into longevity-regulating mechanisms

Anthocyanin Biosynthesis Associated with Natural Variation in Autumn Leaf Coloration in Quercus aliena Accessions

Quercus aliena is an economically important tree species and one of the dominant native oak species in China. Although its leaves typically turn yellow in autumn, we observed natural variants with red leaves. It is important to understand the mechanisms involved in leaf color variation in this species. Therefore, we compared a Q. aliena tree with yellow leaves and three variants with red leaves at different stages of senescence in order to determine the causes of natural variation. We found that the accumulation of anthocyanins such as cyanidin 3-O-glucoside and cyanidin 3-O-sambubiglycoside had a significant effect on leaf coloration. Gene expression analysis showed upregulation of almost all genes encoding enzymes involved in anthocyanin synthesis in the red-leaved variants during the early and main discoloration stages of senescence. These findings are consistent with the accumulation of anthocyanin in red variants. Furthermore, the variants showed significantly higher expression of transcription factors associated with anthocyanin synthesis, such as those encoded by genes QaMYB1 and QaMYB3. Our findings provide new insights into the physiological and molecular mechanisms involved in autumn leaf coloration in Q. aliena, as well as provide genetic resources for further development and cultivation of valuable ornamental variants of this species.