Dear Editor, Beneficial symbiotic fungi colonize plant tissues,delivering crucial ecosystem services such as carbon sequestration and plant fertilization[1].Specifically,trees that form a nutrient-acquiring symbiosis with mutualistic ectomycorrhizal(ECM)fungi gain advantages from these associations by experiencing enhanced growth rates and increased resilience to both biotic and abiotic stresses[2,3].Despite the vital role ECM fungi play in the nutrition and well-being of trees,identifying key regulators participating in the molecular communication between plant and fungal cells is still in its early stages[4].The mutualistic relationship between Laccaria bicolor and Populus spp.

Abiotic stresses such as drought change plant-pathogen interactions by affecting both hosts and pathogens. Here, we aimed to unravel the molecular mechanisms underlying forest tree-pathogen interactions under drought stress and subsequent rewatering. We conducted glasshouse experiments involving infection by the stem canker-causing fungal pathogen Chrysoporthe austroafricana under drought stress and rewatering in Eucalyptus grandis and investigated host and pathogen transcriptomic changes using RNA-seq data from our current combined stress experiment as well as previous single stress studies. We found that mild drought stress enhances disease progression while, upon rewatering, pathogen infection delays recovery of leaf stomatal conductance. Transcriptomic changes in the host support increased susceptibility to the pathogen while the in planta fungal transcriptome suggests prioritization of survival in the drought-stressed host. Upon rewatering, changes in the host transcriptome suggest allocation of resources to stress responses at the expense of growth and carbohydrate storage while that of the pathogen indicate downregulation of some fungal metabolic pathways potentially because the pathogen takes advantage of changes in the host. Our study identified key molecular processes and genes that provide mechanistic insights into tree-pathogen interactions under abiotic stresses. This enables prediction of tree resilience under a changing climate and contributes towards future tree health improvement endeavours.

The ionome represents elemental composition in plant tissues and can be an indicator of nutrient status as well as overall plant performance. Thus, identifying genetic determinants governing elemental uptake and storage is an important goal for breeding and engineering biomass feedstocks with improved performance. In this study, we coupled high-throughput ionome characterization of leaf tissues with high-resolution genome-wide association studies (GWAS) to uncover genetic loci that modulate ionomic composition in leaves of poplar (Populus trichocarpa). Significant agreement was observed across the three ionomic profiling platforms tested: inductively coupled plasma-mass spectrometry (ICP-MS), neutron activation analysis (NAA) and laser-induced breakdown spectroscopy (LIBS). Relative quantification of 20 elements using ICP-MS across a population of 584 genotypes, revealed larger variation in micro-nutrients and trace elements content than for macro-nutrients across genotypes. The GWAS performed using a set of high-density (>8.2 million) single nucleotide polymorphisms, identified over 600 loci significantly associated with variations in these mineral elements, pointing to numerous uncharacterized candidate genes. A significant enrichment for genes related to ion homeostasis and transport was observed, including several members of the cation-proton antiporters (CPA) family and MATE efflux transporters, previously reported to be critical for plant growth and fitness in other species. Our results also included a polymorphic copy of the high-affinity molybdenum transporter MOT1 found directly associated to molybdenum content. For the first time in a perennial plant, our results provide evidence of genetic control of mineral content in a model tree species.

While being the most planted tree worldwide, most Eucalyptus species are sensitive to frost. Long-term exposure to cold temperatures, among other abiotic cues, triggers modification of secondary xylem differentiation in Eucalyptus, but the molecular mechanisms behind it remain unknown. Overexpression of key players of cold signalling pathway, the CRT-repeat binding factors (CBF), not only causes the expected increase in freezing tolerance but also remodels wood in a similar trend to cold acclimation, making it a good candidate for wood functional adaptation to cold stress. To gain insight in CBF role in cold-induced secondary cell wall (SCW) remodelling, we used both targeted and untargeted methodologies to show that chilling and freezing temperatures induced the deposition of a thick SCW with alterations in lignin and polysaccharides composition as well as modifications of wood anatomy in a Eucalyptus cold-tolerant hybrid. Using co-expression network approaches, we identified CBF transcription factors (TFs) as a regulatory hub in xylem cold response. Direct targets of the CBF TFs were identified by DAPseq and unravelled promising candidates involved in SCW deposition and hormonal signalling pathways. Our results shed new light on the interplay between cold response and wood formation, bringing new evidence for the role of the cell wall in trees tolerance to abiotic stresses. ### Competing Interest Statement The authors have declared no competing interest.

Populus is a promising lignocellulosic feedstock for biofuels and bioproducts. However, the cell wall biopolymer lignin is a major barrier in conversion of biomass to biofuels. To investigate the variability and underlying genetic basis of the complex structure of lignin, a population of 409 three-year-old, naturally varying Populus trichocarpa genotypes were characterized by heteronuclear single quantum coherence (HSQC) nuclear magnetic resonance (NMR). A subsequent genome-wide association study (GWAS) was conducted using approximately 8.3 million single nucleotide polymorphisms (SNPs), which identified 756 genes that were significantly associated (-log(10)(p-value)>6) with at least one lignin phenotype. Several promising candidate genes were identified, many of which have not previously been reported to be associated with lignin or cell wall biosynthesis. These results provide a resource for gaining insights into the molecular mechanisms of lignin biosynthesis and new targets for future genetic improvement in poplar.

Plant establishment requires the formation and development of an extensive root system with architecture modulated by complex genetic networks. Here, we report the identification of the PtrXB38 gene as an expression quantitative trait loci (eQTL) hotspot, mapped using 390 leaf and 444 xylem Populus trichocarpa transcriptomes. Among predicted targets of this trans-eQTL were genes involved in plant hormone responses and root development. Overexpression of PtrXB38 in Populus led to significant increases in callusing and formation of both stem-born roots and base-born adventitious roots. Omics studies revealed that genes and proteins controlling auxin transport and signaling were involved in PtrXB38-mediated adventitious root formation. Protein-protein interaction assays indicated that PtrXB38 interacts with components of endosomal sorting complexes required for transport machinery, implying that PtrXB38-regulated root development may be mediated by regulating endocytosis pathway. Taken together, this work identified a crucial root development regulator and sheds light on the discovery of other plant developmental regulators through combining eQTL mapping and omics approaches.

To improve our understanding of genetic mechanisms underlying complex traits in plants, a comprehensive analysis of gene variants is required. Eucalyptus is an important forest plantation genus that is highly outbred. Trait dissection and molecular breeding in eucalypts currently relies on biallelic single-nucleotide polymorphism (SNP) markers. These markers fail to capture the large amount of haplotype diversity in these species, and thus multi-allelic markers are required. We aimed to develop a gene-based haplotype mining panel for Eucalyptus species. We generated 17 999 oligonucleotide probe sets for targeted sequencing of selected regions of 6293 genes implicated in growth and wood properties, pest and disease resistance, and abiotic stress responses. We identified and phased 195 834 SNPs using a read-based phasing approach to reveal SNP-based haplotypes. A total of 8915 target regions (at 4637 gene loci) passed tests for Mendelian inheritance. We evaluated the haplotype panel in four Eucalyptus species (E. grandis, E. urophylla, E. dunnii and E. nitens) to determine its ability to capture diversity across eucalypt species. This revealed an average of 3.13-4.52 haplotypes per target region in each species, and 33.36% of the identified haplotypes were shared by at least two species. This haplotype mining panel will enable the analysis of haplotype diversity within and between species, and provide multi-allelic markers that can be used for genome-wide association studies and gene-based breeding approaches.

A multi-tiered transcriptional network regulates xylem differentiation and secondary cell wall (SCW) formation in plants, with evidence of both conserved and lineage-specific SCW network architecture. We aimed to elucidate the roles of selected R2R3-MYB transcription factors (TFs) linked to Eucalyptus wood formation by identifying genome-wide TF binding sites and direct target genes through an improved DAP-seq protocol combined with machine learning for target gene assignment (DAP-seq-ML). We applied this to five TFs including a well-studied SCW master regulator (EgrMYB2; homolog of AtMYB83), a repressor of lignification (EgrMYB1; homolog of AtMYB4), a TF affecting SCW thickness and vessel density (EgrMYB137; homolog of PtrMYB074) and two TFs with unclear roles in SCW regulation (EgrMYB135 and EgrMYB122). Each DAP-seq TF peak set (average 12,613 peaks) was enriched for canonical R2R3-MYB binding motifs. To improve the reliability of target gene assignment to peaks, a random forest classifier was developed from Arabidopsis DAP-seq, RNA-seq, chromatin, and conserved noncoding sequence data which demonstrated significantly higher precision and recall to the baseline method of assigning genes to proximal peaks. EgrMYB1, EgrMYB2 and EgrMYB137 predicted targets showed clear enrichment for SCW-related biological processes. As validation, EgrMYB137 overexpression in transgenic Eucalyptus hairy roots increased xylem lignification, while its dominant repression in transgenic Arabidopsis and Populus reduced xylem lignification, stunted growth, and caused downregulation of SCW genes. EgrMYB137 targets overlapped significantly with those of EgrMYB2, suggesting partial functional redundancy. Our results show that DAP-seq-ML identified biologically relevant R2R3-MYB targets supported by the finding that EgrMYB137 promotes SCW lignification in planta.

The importance of drought as a constraint to agriculture and forestry is increasing with climate change. Genetic improvement of plants' resilience is one of the mitigation strategies to curb this threat. Although recovery from drought stress is important to long-term drought adaptation and has been considered as an indicator of dehydration tolerance in annual crops, this has not been well explored in forest trees. Thus, we aimed to investigate the physiological and transcriptional changes during drought stress and rewatering in Eucalyptus grandis W. Hill ex Maiden. We set up a greenhouse experiment where we imposed drought stress on 2-year-old seedlings and rewatered the recovery group after 17 days of drought. Our measurement of leaf stomatal conductance (gs) showed that, while gs was reduced by drought stress, it fully recovered after 5 days of rewatering. The RNA-seq analysis from stem samples revealed that genes related to known stress responses such as phytohormone and reactive oxygen species signaling were upregulated, while genes involved in metabolism and growth were downregulated due to drought stress. We observed reprogramming of signal transduction pathways and metabolic processes at 1 day of rewatering, indicating a quick response to rewatering. Our results suggest that recovery from drought stress may entail alterations in the jasmonic acid, salicylic acid, ethylene and brassinosteroid signaling pathways. Using co-expression network analysis, we identified hub genes, including the putative orthologs of ABI1, ABF2, ABF3, HAI2, BAM1, GolS2 and SIP1 during drought and CAT2, G6PD1, ADG1 and FD-1 during recovery. Taken together, by highlighting the molecular processes and identifying key genes, this study gives an overview of the mechanisms underlying the response of E. grandis to drought stress and recovery that trees may face repeatedly throughout their long life cycle. This provides a useful reference to the identification and further investigation of signaling pathways and target genes for future tree improvement.

The appearance of lignin, 400 million-years ago, has been crucial for the successful land colonization by plants. This complex phenolic polymer which confers hydrophobicity and resistance to degradation to secondary cell walls, is also neosynthesized in response to both biotic and abiotic stresses. Lignin biosynthesis during development is regulated by a complex hierarchical transcriptional gene regulatory network. Notably, distinct abiotic stresses can co-opt different genes of this network, leading to different lignification patterns. Co-optation of this developmental network is a mean to facilitate adaptation to stress by promoting functional adaptation. Supporting this hypothesis, a growing number of transcription factors known to be involved in the regulation of secondary cell wall, were shown to integrate both developmental and environmental signals leading to an increased tolerance to abiotic stresses, such as drought or salinity. On the other hand, some transcription factors known to be involved in stress response were shown to impact lignin biosynthesis as part of the functional adaptation performed by plants to cope with natural constraints. In this review, we focus on the most recent findings on the cross-talk between stress-signaling pathways and regulation of lignin biosynthesis in response to abiotic stress and to a lesser extent to biotic stresses. We consider the transcriptional regulation level which has received much attention and the emerging roles of posttranscriptional and posttranslational regulations in the control of lignin biosynthesis in response to stresses. We end discussing future outlines and challenges to fill gaps in this emerging and exciting research field.

Background: Affordable high-throughput DNA and RNA sequencing technologies are allowing genomic analysis of plant and animal populations and as a result empowering new systems genetics approaches to study complex traits. The availability of intuitive tools to browse and analyze the resulting large-scale genetic and genomic datasets remain a significant challenge. Furthermore, these integrative genomics approaches require innovative methods to dissect the flow and interconnectedness of biological information underlying complex trait variation. The Plant Genome Integrative Explorer (PlantGenIE.org ) is a multi-species database and domain that houses online tools for model and woody plant species including Eucalyptus. Since the Eucalyptus Genome Integrative Explorer (EucGenIE) is integrated within PlantGenIE, it shares genome and expression analysis tools previously implemented within the various subdomains (ConGenIE, PopGenIE and AtGenIE). Despite the success in setting up integrative genomics databases, online tools for systems genetics modelling and high-resolution dissection of complex trait variation in plant populations have been lacking. Results: We have developed qtlXplorer (https://eucgenie.org/QTLXplorer) for visualizing and exploring systems genetics data from genome-wide association studies including quantitative trait loci (QTLs) and expression-based QTL (eQTL) associations. This module allows users to, for example, find co-located QTLs and eQTLs using an interactive version of Circos, or explore underlying genes using JBrowse. It provides users with a means to build systems genetics models and generate hypotheses from large-scale population genomics data. We also substantially upgraded the EucGenIE resource and show how it enables users to combine genomics and systems genetics approaches to discover candidate genes involved in biotic stress responses and wood formation by focusing on two multigene families, laccases and peroxidases. Conclusions: qtlXplorer adds a new dimension, population genomics, to the EucGenIE and PlantGenIE environment. The resource will be of interest to researchers and molecular breeders working in Eucalyptus and other woody plant species. It provides an example of how systems genetics data can be integrated with functional genetics data to provide biological insight and formulate hypotheses. Importantly, integration within PlantGenIE enables novel comparative genomics analyses to be performed from population-scale data.

Eucalypts are the most planted trees worldwide, but most of them are frost sensitive. Overexpressing transcription factors for CRT-repeat binding factors (CBFs) in transgenic Eucalyptus confer cold resistance both in leaves and stems. While wood plays crucial roles in trees and is affected by environmental cues, its potential role in adaptation to cold stress has been neglected. Here, we addressed this question by investigating the changes occurring in wood in response to the overexpression of two CBFs, taking advantage of available transgenic Eucalyptus lines. We performed histological, biochemical, and transcriptomic analyses on xylem samples. CBF ectopic expression led to a reduction of both primary and secondary growth, and triggered changes in xylem architecture with smaller and more frequent vessels and fibers exhibiting reduced lumens. In addition, lignin content and syringyl/guaiacyl (S/G) ratio increased. Consistently, many genes of the phenylpropanoid and lignin branch pathway were upregulated. Most of the features of xylem remodeling induced by CBF overexpression are reminiscent of those observed after long exposure of Eucalyptus trees to chilling temperatures. Altogether, these results suggest that CBF plays a central role in the cross-talk between response to cold and wood formation and that the remodeling of wood is part of the adaptive strategies to face cold stress.

The Eucalyptus Genome Integrative Explorer : an online resource for systems genetics in forest tree species

Human life intimately depends on plants for food, biomaterials, health, energy, and a sustainable environment. Various plants have been genetically improved mostly through breeding, along with limited modification via genetic engineering, yet they are still not able to meet the ever-increasing needs, in terms of both quantity and quality, resulting from the rapid increase in world population and expected standards of living. A step change that may address these challenges would be to expand the potential of plants using biosystems design approaches. This represents a shift in plant science research from relatively simple trial-and-error approaches to innovative strategies based on predictive models of biological systems. Plant biosystems design seeks to accelerate plant genetic improvement using genome editing and genetic circuit engineering or create novel plant systems through de novo synthesis of plant genomes. From this perspective, we present a comprehensive roadmap of plant biosystems design covering theories, principles, and technical methods, along with potential applications in basic and applied plant biology research. We highlight current challenges, future opportunities, and research priorities, along with a framework for international collaboration, towards rapid advancement of this emerging interdisciplinary area of research. Finally, we discuss the importance of social responsibility in utilizing plant biosystems design and suggest strategies for improving public perception, trust, and acceptance.

While potassium fertilization increases growth yield in Brazilian eucalyptus plantations, it could also increase water requirements, making trees more vulnerable to drought. Sodium fertilization, which has been shown to promote eucalyptus growth compared to K-deficient trees, could partially mitigate this adverse effect of potassium. However, little is known about the influence of K and Na fertilization on the tree metabolic response to water deficit. The aim of the present study was thus to analyze the transcriptome of leaves sampled from Eucalyptus grandis trees subjected to 37% rainfall reduction, and fertilized with potassium (K), sodium (Na), compared to control trees (C). The multifactorial experiment was set up in a field with a throughfall exclusion system. Transcriptomic analysis was performed on leaves from two-year-old trees, and data analyzed using multifactorial statistical analysis and weighted gene co-expression network analysis (WGCNA). Significant sets of genes were seen to respond to rainfall reduction, in interaction with K or Na fertilization, or to fertilization only (regardless of the water supply regime). The genes were involved in stress signaling, primary and secondary metabolism, secondary cell wall formation and photosynthetic activity. Our focus on key genes related to cation transporters and aquaporins highlighted specific regulation of ion homeostasis, and plant adjustment to water deficit. While water availability significantly affects the transcriptomic response of eucalyptus species, this study points out that the transcriptomic response is highly dependent on the fertilization regime. Our study is based on the first large-scale field trial in a tropical region, specifically designed to study the interaction between water availability and nutrition in eucalyptus. To our knowledge, this is the first global transcriptomic analysis to compare the influence of K and Na fertilization on tree adaptive traits in water deficit conditions.

Wood production in fast-growing Eucalyptus grandis trees is highly dependent on both potassium (K) fertilization and water availability but the molecular processes underlying wood formation in response to the combined effects of these two limiting factors remain unknown. E. grandis trees were submitted to four combinations of K-fertilization and water supply. Weighted gene co-expression network analysis and MixOmics-based co-regulation networks were used to integrate xylem transcriptome, metabolome and complex wood traits. Functional characterization of a candidate gene was performed in transgenic E. grandis hairy roots. This integrated network-based approach enabled us to identify meaningful biological processes and regulators impacted by K-fertilization and/or water limitation. It revealed that modules of co-regulated genes and metabolites strongly correlated to wood complex traits are in the heart of a complex trade-off between biomass production and stress responses. Nested in these modules, potential new cell-wall regulators were identified, as further confirmed by the functional characterization of EgMYB137. These findings provide new insights into the regulatory mechanisms of wood formation under stressful conditions, pointing out both known and new regulators co-opted by K-fertilization and/or water limitation that may potentially promote adaptive wood traits.

Wood, which represents the most abundant lignocellulosic biomass on earth, fulfils key roles in trees, and is also a raw material for multiple end-uses by mankind. The differentiation of this complex vascular tissue starts with cell division in the vascular cambium and is characterized by a massive deposition of lignified secondary cell walls mainly in fibres and vessels. A transcriptional network underlying this differentiation process ensures a tight regulation of the expression of thousands of genes both at the spatial and temporal levels. Most of our current knowledge of this hierarchical network was extrapolated from studies performed in Arabidopsis. Here, we review recent findings on the regulation of wood formation in angiosperm trees species highlighting conserved and distinct mechanisms with Arabidopsis. We provide examples shedding light on the central role of auxin and its cross-talk with other hormones at different stages of secondary xylem differentiation. Functional studies of trees' wood-associated transcription factors revealed diversified functions as compared to their Arabidopsis orthologs. Sophisticated mechanisms of alternative splicing and cross-regulation between the two distinct groups of top level NAC-domain master regulators were uncovered. These findings underlie the high level of complexity of wood formation in trees and provide a framework for future lines of research in this exiting research field.

Although eucalypts are the most planted hardwood trees worldwide, the majority of them are frost sensitive. The recent creation of frost-tolerant hybrids such as Eucalyptus gundal plants (E. gunnii × E. dalrympleana hybrids), now enables the development of industrial plantations in northern countries. Our objective was to evaluate the impact of cold on the wood structure and composition of these hybrids, and on the biosynthetic and regulatory processes controlling their secondary cell-wall (SCW) formation. We used an integrated approach combining histology, biochemical characterization and transcriptomic profiling as well as gene co-expression analyses to investigate xylem tissues from Eucalyptus hybrids exposed to cold conditions. Chilling temperatures triggered the deposition of thicker and more lignified xylem cell walls as well as regulation at the transcriptional level of SCW genes. Most genes involved in lignin biosynthesis, except those specifically dedicated to syringyl unit biosynthesis, were up-regulated. The construction of a co-expression network enabled the identification of both known and potential new SCW transcription factors, induced by cold stress. These regulators at the crossroads between cold signalling and SCW formation are promising candidates for functional studies since they may contribute to the tolerance of E. gunnii × E. dalrympleana hybrids to cold.

Fast growing Eucalyptus tree is well adapted to various soils and climate environment, but its growth varies strongly according to these factors. Thanks to potassium fertilization, productivity in the south of Brazil is one of the highest in the world. The development of more sustainable cultural practices requires an improved understanding of mineral nutrition, especially in interaction with water stress, which is more and more threatening plant culture. We aimed to characterize the effect of water availability and nutrition supply on wood formation and quality. An experimental design was set up on field with a highly productive Eucalyptus grandis clone planted in a split-plot design, with 2 factors tested in interaction: water availability set up with rainfall exclusion system, and K+ fertilization. We analyzed wood properties and performed large scale analysis of transcriptome (RNAseq) and metabolome in developing xylem. These data were integrated using multivariate statistical analyses and co-regulation networks. We identified promising transcription factors potentially involved in the regulation of wood formation. The functional characterization of one candiate in E. grandis transgenic roots demonstrated its implication in secondary cell wall biosynthesis, confirming the potential of our system biology approach to indentify new key regulators of wood formation in woody plants.

Annotation of the Eucalyptus grandis genome showed a large amplification of the dehydration-responsive element binding 1/C-repeat binding factor (DREB1/CBF) group without recent DREB2 gene duplication compared with other plant species. The present annotation of the CBF and DREB2 genes from a draft of the Eucalyptus gunnii genome sequence reveals at least one additional CBF copy in the E. gunnii genome compared with E. grandis, suggesting that this group is still evolving, unlike the DREB2 group. This study aims to investigate the redundancy/neo- or sub-functionalization of the duplicates and the relative involvement of the two groups in abiotic stress responses in both E. grandis and E. gunnii (lower growth but higher cold resistance). A comprehensive transcriptional analysis using high-throughput quantitative real-time polymerase chain reaction (qRT-PCR) was performed on leaves, stems and roots from the two Eucalyptus species after cold, heat or drought treatment. A large CBF cluster accounted for most of the cold response in all the organs, whereas heat and drought responses mainly involved a small CBF cluster and the DREB2 genes. In addition, CBF putative target genes, known to be involved in plant tolerance and development, were found to be cold-regulated. The higher transcript amounts of both the CBF and target genes in the cold tolerant E. gunnii contrasted with the higher CBF induction rates in the fast growing E. grandis. Altogether, the present results, in agreement with previous data about Eucalyptus transgenic lines over-expressing CBF, suggest that these factors, which promote both stress protection and growth limitation, participate in the trade-off between growth and resistance in this woody species.