The chromosome-scale genome of Phoebe bournei reveals contrasting fates of terpene synthase (TPS)-a and TPS-b subfamilies

Identification and variation analysis of the composition and content of essential oil and fragrance compounds in Phoebe zhennan wood at different tree ages

Wood essential oil and wood products with special fragrances are high value-added forest products. Despite the availability of essential oil and volatile organic compounds (VOCs) from Phoebe zhennan wood, their variation and dependence on tree age have not been examined. After essential oil extraction and wood processing, the yields and compositions of essential oils and VOCs in wood from P. zhennan trees of different ages (10a, 30a, and 80a) were determined. The yield of essential oil from 30a wood was significantly greater than that from 10a and 80a wood. Liquid chromatography-mass spectrometry (LC-MS) and gas chromatography-mass spectrometry (GC-MS) revealed 672 and 41 volatile compounds, respectively, in the essential oil and wood, the majority of which exhibited large fluctuations in relative content and composition depending on tree age. Sesquiterpenoids, fatty acids and conjugates may greatly contribute to the main components of essential oil from wood. Almost all major sesquiterpenoid compounds, such as caryophyllene α-oxide, eudesmo, and cubebene, were identified in the essential oils from the 30a and 80a wood, and their relative contents were much greater than those in the 10a wood. The main components of the wood fragrance were sesquiterpenoids. The types and relative contents of sesquiterpenoids from wood increased with tree age. These results suggest that choosing wood from trees of a suitable age will significantly improve the efficiency of wood utilization.

Genome-Wide Identification, Characterization, and Expression Analysis of the BES1 Family Genes under Abiotic Stresses in Phoebe bournei

The BRI1 EMS suppressor 1(BES1) transcription factor is a crucial regulator in the signaling pathway of Brassinosteroid (BR) and plays an important role in plant growth and response to abiotic stress. Although the identification and functional validation of BES1 genes have been extensively explored in various plant species, the understanding of their role in woody plants-particularly the endangered species Phoebe bournei (Hemsl.) Yang-remains limited. In this study, we identified nine members of the BES1 gene family in the genome of P. bournei; these nine members were unevenly distributed across four chromosomes. In our further evolutionary analysis of PbBES1, we discovered that PbBES1 can be divided into three subfamilies (Class I, Class II, and Class IV) based on the evolutionary tree constructed with Arabidopsis thaliana, Oryza sativa, and Solanum lycopersicum. Each subfamily contains 2-5 PbBES1 genes. There were nine pairs of homologous BES1 genes in the synteny analysis of PbBES1 and AtBES1. Three segmental replication events and one pair of tandem duplication events were present among the PbBES1 family members. Additionally, we conducted promoter cis-acting element analysis and discovered that PbBES1 contains binding sites for plant growth and development, cell cycle regulation, and response to abiotic stress. PbBES1.2 is highly expressed in root bark, stem bark, root xylem, and stem xylem. PbBES1.3 was expressed in five tissues. Moreover, we examined the expression profiles of five representative PbBES1 genes under heat and drought stress. These experiments preliminarily verified their responsiveness and functional roles in mediating responses to abiotic stress. This study provides important clues to elucidate the functional characteristics of the BES1 gene family, and at the same time provides new insights and valuable information for the regulation of resistance in P. bournei.

Analysis of the UDP-Glucosyltransferase (UGT) Gene Family and Its Functional Involvement in Drought and Salt Stress Tolerance in Phoebe bournei

Uridine diphosphate glycosyltransferases (UDP-GTs, UGTs), which are regulated by UGT genes, play a crucial role in glycosylation. In vivo, the activity of UGT genes can affect the availability of metabolites and the rate at which they can be eliminated from the body. UGT genes can exert their regulatory effects through mechanisms such as post-transcriptional modification, substrate subtype specificity, and drug interactions. Phoebe bournei is an economically significant tree species that is endemic to southern China. Despite extensive studies on the UGT gene family in various species, a comprehensive investigation of the UGT family in P. bournei has not been reported. Therefore, we conducted a systematic analysis to identify 156 UGT genes within the entire P. bournei genome, all of which contained the PSPG box. The PbUGT family consists of 14 subfamilies, consistent with Arabidopsis thaliana. We observed varying expression levels of PbUGT genes across different tissues in P. bournei, with the following average expression hierarchy: leaf > stem xylem > stem bark > root xylem > root bark. Covariance analysis revealed stronger covariance between P. bournei and closely related species. In addition, we stressed the seedlings with 10% NaCl and 10% PEG-6000. The PbUGT genes exhibited differential expression under drought and salt stresses, with specific expression patterns observed under each stress condition. Our findings shed light on the transcriptional response of PbUGT factors to drought and salt stresses, thereby establishing a foundation for future investigations into the role of PbUGT transcription factors.

Genome Identification and Expression Profiling of the PIN-Formed Gene Family in Phoebe bournei under Abiotic Stresses

PIN-formed (PIN) proteins-specific transcription factors that are widely distributed in plants-play a pivotal role in regulating polar auxin transport, thus influencing plant growth, development, and abiotic stress responses. Although the identification and functional validation of PIN genes have been extensively explored in various plant species, their understanding in woody plants-particularly the endangered species Phoebe bournei (Hemsl.) Yang-remains limited. P. bournei is an economically significant tree species that is endemic to southern China. For this study, we employed bioinformatics approaches to screen and identify 13 members of the PIN gene family in P. bournei. Through a phylogenetic analysis, we classified these genes into five sub-families: A, B, C, D, and E. Furthermore, we conducted a comprehensive analysis of the physicochemical properties, three-dimensional structures, conserved motifs, and gene structures of the PbPIN proteins. Our results demonstrate that all PbPIN genes consist of exons and introns, albeit with variations in their number and length, highlighting the conservation and evolutionary changes in PbPIN genes. The results of our collinearity analysis indicate that the expansion of the PbPIN gene family primarily occurred through segmental duplication. Additionally, by predicting cis-acting elements in their promoters, we inferred the potential involvement of PbPIN genes in plant hormone and abiotic stress responses. To investigate their expression patterns, we conducted a comprehensive expression profiling of PbPIN genes in different tissues. Notably, we observed differential expression levels of PbPINs across the various tissues. Moreover, we examined the expression profiles of five representative PbPIN genes under abiotic stress conditions, including heat, cold, salt, and drought stress. These experiments preliminarily verified their responsiveness and functional roles in mediating responses to abiotic stress. In summary, this study systematically analyzes the expression patterns of PIN genes and their response to abiotic stresses in P. bournei using whole-genome data. Our findings provide novel insights and valuable information for stress tolerance regulation in P. bournei. Moreover, the study offers significant contributions towards unraveling the functional characteristics of the PIN gene family.

Co-expression network analysis reveals PbTGA4 and PbAPRR2 as core transcription factors of drought response in an important timber species Phoebe bournei

Phoebe bournei is one of the main afforestation tree species in subtropical regions of China and is famous for its timber. Its distribution and growth are significantly impaired by water conditions. Thus, it is essential to understand the mechanism of the stress response in P. bournei. Here, we analyzed the phenotypic changes and transcriptomic rearrangement in the leaves and roots of P. bournei seedlings grown for 0 h, 1 h, 24 h, and 72 h under simulated drought conditions (10% PEG 6000). The results showed that drought stress inhibited plant photosynthesis and increased oxidoreductase activity and abscisic acid (ABA) accumulation. Spatio-temporal transcriptomic analysis identified 2836 and 3704 differentially expressed genes (DEGs) in leaves and roots, respectively. The responsive genes in different organs presented various expression profiles at different times. Gene co-expression network analysis identified two core transcription factors, TGA4 and APRR2, from two modules that showed a strong positive correlation with ABA accumulation. Our study investigated the different responses of aboveground and belowground organs of P. bournei to drought stress and provides critical information for improving the drought resistance of this timber species.

Genome-Wide Identification and Expression Analysis of Calmodulin and Calmodulin-like Genes, Revealing CaM3 and CML13 Participating in Drought Stress in Phoebe bournei

Calmodulin (CaM) and calmodulin-like (CML) proteins are major Ca2+ sensors involved in the regulation of plant development and stress responses by converting Ca2+ signals into appropriate cellular responses. However, characterization and expression analyses of CaM/CML genes in the precious species, Phoebe bournei, remain limited. In this study, five PbCaM and sixty PbCML genes were identified that only had EF-hand motifs with no other functional domains. The phylogenetic tree was clustered into 11 subgroups, including a unique clade of PbCaMs. The PbCaMs were intron-rich with four EF-hand motifs, whereas PbCMLs had two to four EF-hands and were mostly intronless. PbCaMs/CMLs were unevenly distributed across the 12 chromosomes of P. bournei and underwent purifying selection. Fragment duplication was the main driving force for the evolution of the PbCaM/CML gene family. Cis-acting element analysis indicated that PbCaMs/CMLs might be related to hormones, growth and development, and stress response. Expression analysis showed that PbCaMs were generally highly expressed in five different tissues and under drought stress, whereas PbCMLs showed specific expression patterns. The expression levels of 11 candidate PbCaMs/CMLs were responsive to ABA and MeJA, suggesting that these genes might act through multiple signaling pathways. The overexpression of PbCaM3/CML13 genes significantly increased the tolerance of yeast cells to drought stress. The identification and characterization of the CaM/CML gene family in P. bournei laid the foundation for future functional studies of these genes.

LcMYB106 suppresses monoterpene biosynthesis by negatively regulating LcTPS32 expression in Litsea cubeba

Litsea cubeba, the core species of the Lauraceae family, is valuable for the production of essential oils due to its high concentration of monoterpenes (90%). The key monoterpene synthase and metabolic regulatory network of monoterpene biosynthesis have provided new insights for improving essential oil content. However, there are few studies on the regulation mechanism of monoterpenes in L. cubeba. In this study, we investigated LcTPS32, a member of the TPS-b subfamily, and identified its function as an enzyme for the synthesis of monoterpenes, including geraniol, α-pinene, β-pinene, β-myrcene, linalool and eucalyptol. The quantitative real-time PCR analysis showed that LcTPS32 was highly expressed in the fruits of L. cubeba and contributed to the characteristic flavor of its essential oil. Overexpression of LcTPS32 resulted in a significant increase in the production of monoterpenes in L. cubeba by activating both the MVA and MEP pathways. Additionally, the study revealed that LcMYB106 played a negative regulatory role in monoterpenes biosynthesis by directly binding to the promoter of LcTPS32. Our study indicates that LcMYB106 could serve as a crucial target for metabolic engineering endeavors, aiming at enhancing the monoterpene biosynthesis in L. cubeba.

Comparative Metabolomics and Transcriptome Analysis Reveal the Fragrance-Related Metabolite Formation in Phoebe zhennan Wood

Nanmu (Phoebe zhennan) has a unique fragrance and is a high-quality tree species for forest conservation. The types and contents of volatile compounds in different tissues of nanmu wood are different, and the study of its volatile metabolites can help us to understand the source of its fragrance and functions. In order to explore the metabolites related to the wood fragrance of nanmu and to find out the unique volatile substances in the heartwood, gas chromatography-mass spectrometry (GC-MS) was performed to analyze the non-targeted metabolomics in five radial tissues from the sapwood to the heartwood of nanmu. A total of 53 volatile metabolites belonging to 11 classes were detected in all tissues, including terpenes, aromatic hydrocarbons, organoheterocyclics, phenols, esters, organic acids, alcohols, alkaloids, alkane, indoles derivatives, and others. And most of the volatile metabolites were identified for the first time in nanmu wood. Among them, terpenes and aromatic hydrocarbons were the main volatile components. In addition, 22 differential metabolites were screened from HW and SW, HW, and TZ via metabolomic analysis. Among these DAMs, three volatile metabolites (cadinene, a sesquiterpenoid; p-cymene, a monoterpenoid; 1,3,5-triisopropylbenzene, an aromatic hydrocarbon) contributed heavily to the characteristic fragrance of the heartwood. Additionally, the expression of transcripts showed that the unigenes in the terpenoid biosynthesis pathway were especially up-regulated in the SW. Therefore, we speculated that fragrance-related metabolites were synthesized in SW and then deposited in heartwood during sapwood transformed to heartwood. The expression levels of transcription factors (e.g., WRKY, C2H2, NAC) acted as the major regulatory factors in the synthesis of terpenoid. The results lay the foundations for further studies on the formation mechanism of fragrance components in nanmu wood and also provide a reference for the further development and utilization of nanmu wood.

Unlocking secrets of nature's chemists: Potential of CRISPR/Cas-based tools in plant metabolic engineering for customized nutraceutical and medicinal profiles

Plant species have evolved diverse metabolic pathways to effectively respond to internal and external signals throughout their life cycle, allowing adaptation to their sessile and phototropic nature. These pathways selectively activate specific metabolic processes, producing plant secondary metabolites (PSMs) governed by genetic and environmental factors. Humans have utilized PSM-enriched plant sources for millennia in medicine and nutraceuticals. Recent technological advances have significantly contributed to discovering metabolic pathways and related genes involved in the biosynthesis of specific PSM in different food crops and medicinal plants. Consequently, there is a growing demand for plant materials rich in nutrients and bioactive compounds, marketed as "superfoods". To meet the industrial demand for superfoods and therapeutic PSMs, modern methods such as system biology, omics, synthetic biology, and genome editing (GE) play a crucial role in identifying the molecular players, limiting steps, and regulatory circuitry involved in PSM production. Among these methods, clustered regularly interspaced short palindromic repeats-CRISPR associated protein (CRISPR/Cas) is the most widely used system for plant GE due to its simple design, flexibility, precision, and multiplexing capabilities. Utilizing the CRISPR-based toolbox for metabolic engineering (ME) offers an ideal solution for developing plants with tailored preventive (nutraceuticals) and curative (therapeutic) metabolic profiles in an ecofriendly way. This review discusses recent advances in understanding the multifactorial regulation of metabolic pathways, the application of CRISPR-based tools for plant ME, and the potential research areas for enhancing plant metabolic profiles.

Temporospatial pattern of flavonoid metabolites and potential regulatory pathway of PbMYB211-coordinated kaempferol-3-O-rhamnoside biosynthesis in Phoebe bournei

Phoebe is a well-known timber tree species that contains abundant metabolites characterized by flavonoids that are widely used in the pharmaceutical industry. Nevertheless, temporospatial flavonoid metabolism variations substantially impact the Phoebe industry. Thus, a metabolomics analysis was carried out and identified 465 metabolites (102 flavonoids) in P. bournei, revealing distinct distribution patterns among five studied organs, and most of the flavonoids were dominant in the leaves. Furthermore, three kaempferol glycoside derivatives were significantly accumulated in the leaves and showed higher contents in young leaves than in mature leaves and differences between spring and autumn. For instance, greater accumulation of kaempferol-7-O-rhamnoside was detected in spring, whereas higher contents of kaempferol-3-O-arabinofuranoside and kaempferol-3-O-rhamnoside were found in autumn. Integrated metabolomics and transcriptomics identified 20 transcription factors (TFs) and 12 structural genes that participate in kaempferol derivative synthesis and elucidated a potential regulatory mechanism in P. bournei. Of the identified genes, PbMYB211 might contribute significantly to the kaempferol-3-O-rhamnoside content by regulating the target structural gene PbUGT139, as revealed by transient overexpression analysis. Overall, this study illuminated the temporospatial accumulation of flavonoids among different organs, seasons, and developmental stages in P. bournei and elucidated a potential regulatory pathway of kaempferol-3-O-rhamnoside. The results provide important insights into harvest techniques and a theoretical basis for the comprehensive utilization of P. bournei.

Chromosome-scale genomes of commercial timber trees (Ochroma pyramidale, Mesua ferrea, and Tectona grandis)

Wood is the most important natural and endlessly renewable source of energy. Despite the ecological and economic importance of wood, many aspects of its formation have not yet been investigated. We performed chromosome-scale genome assemblies of three timber trees (Ochroma pyramidale, Mesua ferrea, and Tectona grandis) which exhibit different wood properties such as wood density, hardness, growth rate, and fiber cell wall thickness. The combination of 10X, stLFR, Hi-Fi sequencing and HiC data led us to assemble high-quality genomes evident by scaffold N50 length of 55.97 Mb (O. pyramidale), 22.37 Mb (M. ferrea) and 14.55 Mb (T. grandis) with >97% BUSCO completeness of the assemblies. A total of 35774, 24027, and 44813 protein-coding genes were identified in M. ferrea, T. grandis and O. pyramidale, respectively. The data generated in this study is anticipated to serve as a valuable genetic resource and will promote comparative genomic analyses, and it is of practical importance in gaining a further understanding of the wood properties in non-model woody species.

Genome-Wide Identification of GATA Family Genes in Phoebe bournei and Their Transcriptional Analysis under Abiotic Stresses

GATA transcription factors are crucial proteins in regulating transcription and are characterized by a type-IV zinc finger DNA-binding domain. They play a significant role in the growth and development of plants. While the GATA family gene has been identified in several plant species, it has not yet been reported in Phoebe bournei. In this study, 22 GATA family genes were identified from the P. bournei genome, and their physicochemical properties, chromosomal distribution, subcellular localization, phylogenetic tree, conserved motif, gene structure, cis-regulatory elements in promoters, and expression in plant tissues were analyzed. Phylogenetic analysis showed that the PbGATAs were clearly divided into four subfamilies. They are unequally distributed across 11 out of 12 chromosomes, except chromosome 9. Promoter cis-elements are mostly involved in environmental stress and hormonal regulation. Further studies showed that PbGATA11 was localized to chloroplasts and expressed in five tissues, including the root bark, root xylem, stem bark, stem xylem, and leaf, which means that PbGATA11 may have a potential role in the regulation of chlorophyll synthesis. Finally, the expression profiles of four representative genes, PbGATA5, PbGATA12, PbGATA16, and PbGATA22, under drought, salinity, and temperature stress, were detected by qRT-PCR. The results showed that PbGATA5, PbGATA22, and PbGATA16 were significantly expressed under drought stress. PbGATA12 and PbGATA22 were significantly expressed after 8 h of low-temperature stress at 10 °C. This study concludes that the growth and development of the PbGATA family gene in P. bournei in coping with adversity stress are crucial. This study provides new ideas for studying the evolution of GATAs, provides useful information for future functional analysis of PbGATA genes, and helps better understand the abiotic stress response of P. bournei.

Genome-Wide Identification, Expression and Stress Analysis of the GRAS Gene Family in Phoebe bournei

GRAS genes are important transcriptional regulators in plants that govern plant growth and development through enhancing plant hormones, biosynthesis, and signaling pathways. Drought and other abiotic factors may influence the defenses and growth of Phoebe bournei, which is a superb timber source for the construction industry and building exquisite furniture. Although genome-wide identification of the GRAS gene family has been completed in many species, that of most woody plants, particularly P. bournei, has not yet begun. We performed a genome-wide investigation of 56 PbGRAS genes, which are unequally distributed across 12 chromosomes. They are divided into nine subclades. Furthermore, these 56 PbGRAS genes have a substantial number of components related to abiotic stress responses or phytohormone transmission. Analysis using qRT-PCR showed that the expression of four PbGRAS genes, namely PbGRAS7, PbGRAS10, PbGRAS14 and PbGRAS16, was differentially increased in response to drought, salt and temperature stresses, respectively. We hypothesize that they may help P. bournei to successfully resist harsh environmental disturbances. In this work, we conducted a comprehensive survey of the GRAS gene family in P. bournei plants, and the results provide an extensive and preliminary resource for further clarification of the molecular mechanisms of the GRAS gene family in P. bournei in response to abiotic stresses and forestry improvement.

Cyanidin-3-O-glucoside Contributes to Leaf Color Change by Regulating Two bHLH Transcription Factors in Phoebe bournei

Anthocyanins produce different-colored pigments in plant organs, which provide ornamental value. Thus, this study was conducted to understand the mechanism of anthocyanin synthesis in ornamental plants. Phoebe bournei, a Chinese specialty tree, has high ornamental and economic value due to its rich leaf color and diverse metabolic products. Here, the metabolic data and gene expression of red P. bournei leaves at the three developmental stages were evaluated to elucidate the color-production mechanism in the red-leaved P. bournei species. First, metabolomic analysis identified 34 anthocyanin metabolites showing high levels of cyanidin-3-O-glucoside (cya-3-O-glu) in the S1 stage, which may suggest that it is a characteristic metabolite associated with the red coloration of the leaves. Second, transcriptome analysis showed that 94 structural genes were involved in anthocyanin biosynthesis, especially flavanone 3'-hydroxy-lase (PbF3'H), and were significantly correlated with the cya-3-O-glu level. Third, K-means clustering analysis and phylogenetic analyses identified PbbHLH1 and PbbHLH2, which shared the same expression pattern as most structural genes, indicating that these two PbbHLH genes may be regulators of anthocyanin biosynthesis in P. bournei. Finally, overexpression of PbbHLH1 and PbbHLH2 in Nicotiana tabacum leaves triggered anthocyanin accumulation. These findings provide a basis for cultivating P. bournei varieties that have high ornamental value.

Phylogenetic studies of magnoliids: Advances and perspectives

Magnoliids are the largest flowering plant clades outside of the eudicots and monocots, which are distributed worldwide and have high economic, ornamental and ecological values. Eudicots, monocots and magnoliids are the three major clades of Mesangiospermae, and their phylogenetic relationship is one of the most interesting issues. In recent years, with the continuous accumulation of genomic information, the evolutionary status of magnoliids has become a hot spot in plant phylogenetic research. Although great efforts have been made to study the evolution of magnoliids using molecular data from several representative species such as nuclear genome, plastid genome, mitochondrial genome, and transcriptome, the results of current studies on the phylogenetic status of magnoliids are inconsistent. Here, we systematically describe the current understanding of the molecular research on magnoliid phylogeny and review the differences in the evolutionary state of magnoliids. Understanding the research approaches and limitations of magnoliid phylogeny can guide research strategies to further improve the study of the phylogenetic evolution of magnoliids.

Integrated Transcriptomic, Metabolomic, and Physiological Analyses Reveal New Insights into Fragrance Formation in the Heartwood of Phoebe hui

Phoebe hui is an extremely valuable tree that is the main source of the fragrant golden-thread nanmu wood. Although the fragrance of wood has been investigated in several trees, the potential substances and gene regulation mechanisms that are involved in fragrance formation are poorly understood. Here, three radial tissues, sapwood (SW), heartwood (HW), and the transition zone (TZ) in between them, were compared via integrative physiological, volatile-metabolomic, and transcriptomic analyses to identify the key metabolites and regulatory mechanisms involved in fragrance formation. During heartwood formation, gradual starch grain loss was accompanied by the deposition of lipids and extractives in the cell lumen. Extracts of terpenoids were synthesized and accumulated in the heartwood, including monoterpenoids (limonene and p-cymene) and sesquiterpenes (cubebene and guaiadiene); these were identified as being closely related to the special fragrance of the wood. Additionally, the expression of transcripts showed that the genes related to primary metabolism were specifically upregulated in the SW, whereas genes annotated in terpenoid biosynthesis were specifically upregulated in the HW. Therefore, we speculated that terpenoid biosynthesis occurs in situ in the HW via the HW formation model of Type-III (Santalum) using the precursors that were produced by primary metabolism in the SW. The expression levels of transcription factors (e.g., MYB, WRKY, and C2H2) acted as the major regulatory factors in the synthesis of terpenoids. Our results explain the special fragrance in P. hui and broaden the current knowledge of the regulatory mechanisms of fragrance formation. This work provides a framework for future research that is focused on improving wood quality and value.