Chromosome-scale assembly of the Kandelia obovata genome

In silico analysis of trehalose biosynthesis genes provides clues to reveal its roles in Avicennia marina adaptation to tidal submergence

Trehalose has an important function for alleviating various abiotic stress in plants. Nevertheless, the functional and evolutionary characteristics of trehalose biosynthesis genes in mangrove plants is not documented. Here, using typical mangrove Avicennia marina, we found the trehalose content decreased in the roots and leaves and T6P increased significantly in the leaves under tidal submergence. Then, the basic physicochemical properties and gene structure of trehalose biosynthesis genes (AmTPS and AmTPP), and the conserved domain and motifs of AmTPS and AmTPP proteins were analyzed. The collinearity analysis and Ka/Ks values indicated that AmTPS and AmTPP are evolutionarily conserved. Tissue-specific expression profiling showed that most AmTPS and AmTPP genes have tissue specificity. RNA-Seq analysis showed that five AmTPS genes were markedly up-regulated in A. marina treated with tidal submergence. Subcellular localization analysis revealed three genes including AmTPS10B, AmTPS11A and AmTPS11C out of these five up-regulated AmTPS genes work in plasma membrane, cytoplasm and nucleus. Finally, integrative analysis of bioinformatics and RNA-Seq analysis were performed to excavate transcription factors that may regulate AmTPS and AmTPP genes in A. marina response to submergence. Taken together, these findings provide new insights into the response to tidal submergence in A. marina at the aspect of trehalose.

MangroveDB: A Comprehensive Online Database for Mangroves Based on Multi-Omics Data

Mangroves are dominant flora of intertidal zones along tropical and subtropical coastline around the world that offer important ecological and economic value. Recently, the genomes of mangroves have been decoded, and massive omics data were generated and deposited in the public databases. Reanalysis of multi-omics data can provide new biological insights excluded in the original studies. However, the requirements for computational resource and lack of bioinformatics skill for experimental researchers limit the effective use of the original data. To fill this gap, we uniformly processed 942 transcriptome data, 386 whole-genome sequencing data, and provided 13 reference genomes and 40 reference transcriptomes for 53 mangroves. Finally, we built an interactive web-based database platform MangroveDB (https://github.com/Jasonxu0109/MangroveDB), which was designed to provide comprehensive gene expression datasets to facilitate their exploration and equipped with several online analysis tools, including principal components analysis, differential gene expression analysis, tissue-specific gene expression analysis, GO and KEGG enrichment analysis. MangroveDB not only provides query functions about genes annotation, but also supports some useful visualization functions for analysis results, such as volcano plot, heatmap, dotplot, PCA plot, bubble plot, population structure, and so on. In conclusion, MangroveDB is a valuable resource for the mangroves research community to efficiently use the massive public omics datasets.

Unveiling the Role of β-Glucosidase Genes in Bletilla striata's Secondary Metabolism: A Genome-Wide Analysis

β-glucosidases (BGLUs) are abundant enzymes in plants that play pivotal roles in cell wall modification, hormone signal transduction, secondary metabolism, defense against herbivores, and volatile compound release. Bletilla striata, a perennial herb revered for its therapeutic properties, lacks a comprehensive analysis of its BGLU gene family despite the critical role these genes play in plant secondary metabolism. This study aims to perform a genome-wide analysis of the BGLU gene family in B. striata (BsBGLU) to elucidate their functions and regulatory mechanisms in secondary metabolite biosynthesis. We conducted a genome-wide screening to identify BsBGLU, followed by phylogenetic analysis to classify these genes into groups. Sequence characteristics were analyzed to predict functional roles. Simple sequence repeat (SSR) markers were examined to assess conservation and polymorphism among different landraces. Expression profiles of BsBGLUs were evaluated under sodium acetate and salicylic acid elicitor treatments and across different tissues. The accumulation of phylogenetic metabolites in different treatments and tissues was also analyzed by HPLC and LCMS detection to explore the correlation between gene expression and metabolite accumulation. A total of 23 BsBGLU genes were identified and classified into eight distinct groups. Sequence analysis suggested diverse functions related to hormone responses, secondary metabolism, and stress resistance. BsBGLUs with SSR sequences were conserved yet showed polymorphism among different B. striata landraces. Under elicitor treatments, expression profiling revealed that BsBGLUs significantly modulate the synthesis of secondary metabolites such as dactylorhin A and militarine. Tissue-specific expression analysis indicated that BsBGLU15 and BsBGLU28 were highly expressed in tubers compared to other tissues, suggesting their central role and a potential negative regulatory effect in metabolite accumulation. The elicitor NaAc can regulate metabolite synthesis by modulating the expression of BsBGLUs. The BsBGLU gene family in B. striata is integral to the modulation of secondary metabolite biosynthesis and accumulation and can respond to elicitors to promote the synthesis of militarine. These findings provide a theoretical foundation for the further exploration of BsBGLU gene functions and their regulatory mechanisms, advancing the production of medicinally active compounds in B. striata.

Genome-Wide Identification of Pentatricopeptide Repeat (PPR) Gene Family and Multi-Omics Analysis Provide New Insights Into the Albinism Mechanism of Kandelia obovata Propagule Leaves

Pentatricopeptide repeat (PPR) gene family constitutes one of the largest gene families in plants, which mainly participate in RNA editing and RNA splicing of organellar RNAs, thereby affecting the organellar development. Recently, some evidence elucidated the important roles of PPR proteins in the albino process of plant leaves. However, the functions of PPR genes in the woody mangrove species have not been investigated. In this study, using a typical true mangrove Kandelia obovata, we systematically identified 298 PPR genes and characterized their general features and physicochemical properties, including evolutionary relationships, the subcellular localization, PPR motif type, the number of introns and PPR motifs, and isoelectric point, and so forth. Furthermore, we combined genome-wide association studies (GWAS) and transcriptome analysis to identify the genetic architecture and potential PPR genes associated with propagule leaves colour variations of K. obovata. As a result, we prioritized 16 PPR genes related to the albino phenotype using different strategies, including differentially expressed genes analysis and genetic diversity analysis. Further analysis discovered two genes of interest, namely Maker00002998 (PLS-type) and Maker00003187 (P-type), which were differentially expressed genes and causal genes detected by GWAS analysis. Moreover, we successfully predicted downstream target chloroplast genes (rps14, rpoC1 and rpoC2) bound by Maker00002998 PPR proteins. The experimental verification of RNA editing sites of rps14, rpoC1, and rpoC2 in our previous study and the verification of interaction between Maker00002998 and rps14 transcript using in vitro RNA pull-down assays revealed that Maker00002998 PPR protein might be involved in the post-transcriptional process of chloroplast genes. Our result provides new insights into the roles of PPR genes in the albinism mechanism of K. obovata propagule leaves.

SOS1 gene family in mangrove (Kandelia obovata): Genome-wide identification, characterization, and expression analyses under salt and copper stress

BACKGROUND

Salt Overly Sensitive 1 (SOS1), a plasma membrane Na+/H+ exchanger, is essential for plant salt tolerance. Salt damage is a significant abiotic stress that impacts plant species globally. All living organisms require copper (Cu), a necessary micronutrient and a protein cofactor for many biological and physiological processes. High Cu concentrations, however, may result in pollution that inhibits the growth and development of plants. The function and production of mangrove ecosystems are significantly impacted by rising salinity and copper contamination.

RESULTS

A genome-wide analysis and bioinformatics techniques were used in this study to identify 20 SOS1 genes in the genome of Kandelia obovata. Most of the SOS1 genes were found on the plasma membrane and dispersed over 11 of the 18 chromosomes. Based on phylogenetic analysis, KoSOS1s can be categorized into four groups, similar to Solanum tuberosum. Kandelia obovata's SOS1 gene family expanded due to tandem and segmental duplication. These SOS1 homologs shared similar protein structures, according to the results of the conserved motif analysis. The coding regions of 20 KoSOS1 genes consist of amino acids ranging from 466 to 1221, while the exons include amino acids ranging from 3 to 23. In addition, we found that the 2.0 kb upstream promoter region of the KoSOS1s gene contains several cis-elements associated with phytohormones and stress responses. According to the expression experiments, seven randomly chosen genes experienced up- and down-regulation of their expression levels in response to copper (CuCl2) and salt stressors.

CONCLUSIONS

For the first time, this work systematically identified SOS1 genes in Kandelia obovata. Our investigations also encompassed physicochemical properties, evolution, and expression patterns, thereby furnishing a theoretical framework for subsequent research endeavours aimed at functionally characterizing the Kandelia obovata SOS1 genes throughout the life cycle of plants.

miR397 regulates cadmium stress response by coordinating lignin polymerization in the root exodermis in Kandelia obovata

Secondary lignification of the root exodermis of Kandelia obovata is crucial for its response to adversity such as high salinity and anaerobic environment, and this lignification is also effective in blocking cadmium transport to the roots. However, how the differences in lignification of root exodermis at different developmental stages respond to Cd stress and its regulatory mechanisms have not been revealed. In this study, after analyzing the root structure and cell wall thickness using a Phenom scanning electron microscope as well as measuring cadmium content in the root cell wall, we found that the exodermis of young and mature roots of K. obovata responded to Cd stress through the polymerization of different lignin monomers, forming two different mechanisms: chelation and blocking. Through small RNA sequencing, RLM-5'-RACE and dual luciferase transient expression system, we found that miR397 targets and regulates KoLAC4/17/7 expression. The expression of KoLAC4/17 promoted the accumulation of guaiacyl lignin during lignification and enhanced the binding of cadmium to the cell wall. Meanwhile, KoLAC7 expression promotes the accumulation of syringyl lignin during lignification, which enhances the obstruction of cadmium and improves the tolerance to cadmium. These findings enhance our understanding of the molecular mechanisms underlying the differential lignification of the root exodermis of K. obovata in response to cadmium stress, and provide scientific guidance for the conservation of mangrove forests under heavy metal pollution.

Stomatal dynamics are regulated by leaf hydraulic traits and guard cell anatomy in nine true mangrove species

Stomatal regulation is critical for mangroves to survive in the hyper-saline intertidal zone where water stress is severe and water availability is highly fluctuant. However, very little is known about the stomatal sensitivity to vapour pressure deficit (VPD) in mangroves, and its co-ordination with stomatal morphology and leaf hydraulic traits. We measured the stomatal response to a step increase in VPD in situ, stomatal anatomy, leaf hydraulic vulnerability and pressure-volume traits in nine true mangrove species of five families and collected the data of genome size. We aimed to answer two questions: (1) Does stomatal morphology influence stomatal dynamics in response to a high VPD in mangroves? with a consideration of possible influence of genome size on stomatal morphology; and (2) do leaf hydraulic traits influence stomatal sensitivity to VPD in mangroves? We found that the stomata of mangrove plants were highly sensitive to a step rise in VPD and the stomatal responses were directly affected by stomatal anatomy and hydraulic traits. Smaller, denser stomata was correlated with faster stomatal closure at high VPD across the species of Rhizophoraceae, and stomata size negatively and vein density positively correlated with genome size. Less negative leaf osmotic pressure at the full turgor (πo) was related to higher operating steady-state stomatal conductance (gs); and a higher leaf capacitance (Cleaf) and more embolism resistant leaf xylem were associated with slower stomatal responses to an increase in VPD. In addition, stomatal responsiveness to VPD was indirectly affected by leaf morphological traits, which were affected by site salinity and consequently leaf water status. Our results demonstrate that mangroves display a unique relationship between genome size, stomatal size and vein packing, and that stomatal responsiveness to VPD is regulated by leaf hydraulic traits and stomatal morphology. Our work provides a quantitative framework to better understand of stomatal regulation in mangroves in an environment with high salinity and dynamic water availability.

Genomic analysis of Nypa fruticans elucidates its intertidal adaptations and early palm evolution

Nypa fruticans (Wurmb), a mangrove palm species with origins dating back to the Late Cretaceous period, is a unique species for investigating long-term adaptation strategies to intertidal environments and the early evolution of palms. Here, we present a chromosome-level genome sequence and assembly for N. fruticans. We integrated the genomes of N. fruticans and other palm family members for a comparative genomic analysis, which confirmed that the common ancestor of all palms experienced a whole-genome duplication event around 89 million years ago, shaping the distinctive characteristics observed in this clade. We also inferred a low mutation rate for the N. fruticans genome, which underwent strong purifying selection and evolved slowly, thus contributing to its stability over a long evolutionary period. Moreover, ancient duplicates were preferentially retained, with critical genes having experienced positive selection, enhancing waterlogging tolerance in N. fruticans. Furthermore, we discovered that the pseudogenization of Early Methionine-labelled 1 (EM1) and EM6 in N. fruticans underly its crypto-vivipary characteristics, reflecting its intertidal adaptation. Our study provides valuable genomic insights into the evolutionary history, genome stability, and adaptive evolution of the mangrove palm. Our results also shed light on the long-term adaptation of this species and contribute to our understanding of the evolutionary dynamics in the palm family.

NRAMP gene family in Kandelia obovata: genome-wide identification, expression analysis, and response to five different copper stress conditions

Natural resistance-associated macrophage proteins (NRAMPs) are a class of metal transporters found in plants that exhibit diverse functions across different species. Transporter proteins facilitate the absorption, distribution, and sequestration of metallic elements within various plant tissues. Despite the extensive identification of NRAMP family genes in various species, a full analysis of these genes in tree species is still necessary. Genome-wide identification and bioinformatics analysis were performed to understand the roles of NRAMP genes in copper (CuCl2) stress in Kandelia obovata (Ko). In Arachis hypogaea L., Populus trichocarpa, Vitis vinifera, Phaseolus vulgaris L., Camellia sinensis, Spirodela polyrhiza, Glycine max L. and Solanum lycopersicum, a genome-wide study of the NRAMP gene family was performed earlier. The domain and 3D structural variation, phylogenetic tree, chromosomal distributions, gene structure, motif analysis, subcellular localization, cis-regulatory elements, synteny and duplication analysis, and expression profiles in leaves and CuCl2 were all investigated in this research. In order to comprehend the notable functions of the NRAMP gene family in Kandelia obovata, a comprehensive investigation was conducted at the genomic level. This study successfully found five NRAMP genes, encompassing one gene pair resulting from whole-genome duplication and a gene that had undergone segmental duplication. The examination of chromosomal position revealed an unequal distribution of the KoNRAMP genes across chromosomes 1, 2, 5, 7, and 18. The KoNRAMPs can be classified into three subgroups (I, II, and SLC) based on phylogeny and synteny analyses, similar to Solanum lycopersicum. Examining cis-regulatory elements in the promoters revealed five hormone-correlated responsive elements and four stress-related responsive elements. The genomic architecture and properties of 10 highly conserved motifs are similar among members of the NRAMP gene family. The conducted investigations demonstrated that the expression levels of all five genes exhibited alterations in response to different levels of CuCl2 stress. The results of this study offer crucial insights into the roles of KoNRAMPs in the response of Kandelia obovata to CuCl2 stress.

Genome-Wide Identification, Characterization, and Expression Analysis of the Copper-Containing Amine Oxidase Gene Family in Mangrove Kandelia obovata

Copper-containing amine oxidases (CuAOs) are known to have significant involvement in the process of polyamine catabolism, as well as serving crucial functions in plant development and response to abiotic stress. A genome-wide investigation of the CuAO protein family was previously carried out in sweet orange (Citrus sinensis) and sweet cherry (Prunus avium L.). Six CuAO (KoCuAO1-KoCuAO6) genes were discovered for the first time in the Kandelia obovata (Ko) genome through a genome-wide analysis conducted to better understand the key roles of the CuAO gene family in Kandelia obovata. This study encompassed an investigation into various aspects of gene analysis, including gene characterization and identification, subcellular localization, chromosomal distributions, phylogenetic tree analysis, gene structure analysis, motif analysis, duplication analysis, cis-regulatory element identification, domain and 3D structural variation analysis, as well as expression profiling in leaves under five different treatments of copper (CuCl2). Phylogenetic analysis suggests that these KoCuAOs, like sweet cherry, may be subdivided into three subgroups. Examining the chromosomal location revealed an unequal distribution of the KoCuAO genes across four out of the 18 chromosomes in Kandelia obovata. Six KoCuAO genes have coding regions with 106 and 159 amino acids and exons with 4 and 12 amino acids. Additionally, we discovered that the 2.5 kb upstream promoter region of the KoCuAOs predicted many cis elements linked to phytohormones and stress responses. According to the expression investigations, CuCl2 treatments caused up- and downregulation of all six genes. In conclusion, our work provides a comprehensive overview of the expression pattern and functional variety of the Kandelia obovata CuAO gene family, which will facilitate future functional characterization of each KoCuAO gene.

Genome-Wide Analysis of bHLH Family Genes and Identification of Members Associated with Cold/Drought-Induced Photoinhibition in Kandelia obovata

Plant basic helix-loop-helix (bHLH) transcription factors play pivotal roles in responding to stress, including cold and drought. However, it remains unclear how bHLH family genes respond to these stresses in Kandelia obovata. In this study, we identified 75 bHLH members in K. obovata, classified into 11 subfamilies and unevenly distributed across its 18 chromosomes. Collineation analysis revealed that segmental duplication primarily drove the expansion of KobHLH genes. The KobHLH promoters were enriched with elements associated with light response. Through RNA-seq, we identified several cold/drought-associated KobHLH genes. This correlated with decreased net photosynthetic rates (Pn) in the leaves of cold/drought-treated plants. Weighted gene co-expression network analysis (WGCNA) confirmed that 11 KobHLH genes were closely linked to photoinhibition in photosystem II (PS II). Among them, four Phytochrome Interacting Factors (PIFs) involved in chlorophyll metabolism were significantly down-regulated. Subcellular localization showed that KobHLH52 and KobHLH30 were located in the nucleus. Overall, we have comprehensively analyzed the KobHLH family and identified several members associated with photoinhibition under cold or drought stress, which may be helpfulfor further cold/drought-tolerance enhancement and molecular breeding through genetic engineering in K. obovata.

Genome-Wide Identification and Expression Analysis of the Copper Transporter (COPT/Ctr) Gene Family in Kandelia obovata, a Typical Mangrove Plant

The copper transporter (COPT/Ctr) gene family plays a critical part in maintaining the balance of the metal, and many diverse species depend on COPT to move copper (Cu) across the cell membrane. In Arabidopsis thaliana, Oryza sativa, Medicago sativa, Zea mays, Populus trichocarpa, Vitis vinifera, and Solanum lycopersicum, a genome-wide study of the COPT protein family was performed. To understand the major roles of the COPT gene family in Kandelia obovata (Ko), a genome-wide study identified four COPT genes in the Kandelia obovata genome for the first time. The domain and 3D structural variation, phylogenetic tree, chromosomal distributions, gene structure, motif analysis, subcellular localization, cis-regulatory elements, synteny and duplication analysis, and expression profiles in leaves and Cu were all investigated in this research. Structural and sequence investigations show that most KoCOPTs have three transmembrane domains (TMDs). According to phylogenetic research, these KoCOPTs might be divided into two subgroups, just like Populus trichocarpa. KoCOPT gene segmental duplications and positive selection pressure were discovered by universal analysis. According to gene structure and motif analysis, most KoCOPT genes showed consistent exon-intron and motif organization within the same group. In addition, we found five hormones and four stress- and seven light-responsive cis-elements in the KoCOPTs promoters. The expression studies revealed that all four genes changed their expression levels in response to copper (CuCl2) treatments. In summary, our study offers a thorough overview of the Kandelia obovata COPT gene family's expression pattern and functional diversity, making it easier to characterize each KoCOPT gene's function in the future.

Description and genomic characterization of Gallaecimonas kandeliae sp. nov., isolated from the sediments of mangrove plant Kandelia obovate

The genus Gallaecimonas, proposed by Rodríguez-Blanco et al. (Int J Syst Evol Microbiol 60:504-509, 2010), is mainly isolated from marine environments. So far, only three species have been identified and characterized in this genus. In this study, a new Gallaecimonas strain named Q10^T was isolated from the sediments of mangrove plant Kandelia obovate taken from Dapeng district, Shenzhen, China. Strain Q10^T was a Gram-stain-negative, non-motile, strictly aerobic, rod-shaped bacterium, and grew with 0-8.0% (w/v) NaCl, at 10-45 °C and at pH 5.5-8.5. Phylogenetic analysis indicated that strain Q10^T and the three Gallaecimonas species formed a clade in the tree, with 16S rRNA gene sequence similarities ranging from 96.0 to 97.0%. The major respiratory quinone is Q8. The polar lipids comprised aminolipid, aminophospholipid, diphosphatidylglycerol, glycolipid, phosphatidylethanolamine, phosphatidylglycerol, glycophospholipid and phospholipid. The predominant fatty acids are C_16:0, C_17:1ω8c, summed feature 3 (C_16:1ω7c/C_16:1ω6c), and iso-C_16:0. The complete genome of strain Q10^T is 3,836,841 bp with a G+C content of 62.6 mol%. The orthologous proteins analysis revealed 55 unique proteins in strain Q10^T related to important biological processes, especially three frataxins related to iron-sulfur cluster assembly, which may play a pivotal role in environmental adaptability of this species. Based on polyphasic taxonomic data, strain Q10^T is considered to represent a novel species within the genus Gallaecimonas, for which the name Gallaecimonas kandelia sp. nov. is proposed. The type strain is Q10^T (=KCTC 92860^T=MCCC 1K08421^T). These results contribute to a better understanding of general features and taxonomy of the genus Gallaecimonas.

Expansion and adaptive evolution of the WRKY transcription factor family in Avicennia mangrove trees

Mangroves are adapted to intertidal zones, which present extreme environmental conditions. WRKYs are among the most prominent transcription factors (TFs) in higher plants and act through various interconnected networks to regulate responses to multiple abiotic stressors. Here, based on omic data, we investigated the landscape and evolutionary patterns of WRKYs in the main mangrove genus Avicennia. We found that both the number and the proportion of TFs and WRKYs in Avicennia species exceeded their inland relatives, indicating a significant expansion of WRKYs in Avicennia. We identified 109 WRKY genes in the representative species Avicennia marina. Comparative genomic analysis showed that two recent whole-genome duplication (WGD) events played a critical role in the expansion of WRKYs, and 88% of Avicennia marina WRKYs (AmWRKYs) have been retained following these WGDs. Applying comparative transcriptomics on roots under experimental salt gradients, we inferred that there is high divergence in the expression of WGD-retained AmWRKYs. Moreover, we found that the expression of 16 AmWRKYs was stable between freshwater and moderately saline water but increased when the trees were exposed to high salinity. In particular, 14 duplicates were retained following the two recent WGD events, indicating potential neo- and sub-functionalization. We also found that WRKYs could interact with other upregulated genes involved in signalling pathways and natural antioxidant biosynthesis to enhance salt tolerance, contributing to the adaptation to intertidal zones. Our omic data of the WRKY family in A. marina broadens the understanding of how a TF family relates to the adaptive evolution of mangroves.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-023-00177-y.

Genome-Wide Identification and Bioinformatics Analyses of Host Defense Peptides Snakin/GASA in Mangrove Plants

Host defense peptides (HDPs) are components of plant defensive barriers that resist microbial infection. Members of the Snakin/GASA protein family in plants have functions of regulating plant growth, defense, and bacteriostasis. Most mangrove plants grow in coastal zones. In order to survive in harsh environments, mangrove plants have evolved complex adaptations against microbes. In this study, Snakin/GASA family members were identified and analyzed in the genomes of three mangrove species. Twenty-seven, thirteen, and nine candidate Snakin/GASA family members were found in Avicennia marina, Kandelia obovata, and Aegiceras corniculatum, respectively. These Snakin/GASA family members were identified and categorized into three subfamilies via phylogenetic analysis. The genes coding for the Snakin/GASA family members were unevenly distributed on chromosomes. Collinearity and conservative motif analyses showed that the Snakin/GASA family members in K. obovata and A. corniculatum underwent multiple gene duplication events. Snakin/GASA family member expression in normal leaves and leaves infected with pathogenic microorganisms of the three mangrove species was verified using real-time quantitative polymerase chain reaction. The expression of KoGASA3 and 4, AcGASA5 and 10, and AmGASA1, 4, 5, 15, 18, and 23 increased after microbial infection. This study provides a research basis for the verification of HDPs from mangrove plants and suggests directions for the development and utilization of marine biological antimicrobial peptides.

Diverse mangroves deviate from other angiosperms in their genome size, leaf cell size and cell packing density relationships

BACKGROUND AND AIMS

While genome size limits the minimum sizes and maximum numbers of cells that can be packed into a given leaf volume, mature cell sizes can be substantially larger than their meristematic precursors and vary in response to abiotic conditions. Mangroves are iconic examples of how abiotic conditions can influence the evolution of plant phenotypes.

METHODS

Here, we examined the coordination between genome size, leaf cell sizes, cell packing densities and leaf size in 13 mangrove species across four sites in China. Four of these species occurred at more than one site, allowing us to test the effect of climate on leaf anatomy.

RESULTS

We found that genome sizes of mangroves were very small compared to other angiosperms, but, like other angiosperms, mangrove cells were always larger than the minimum size defined by genome size. Increasing mean annual temperature of a growth site led to higher packing densities of veins (Dv) and stomata (Ds) and smaller epidermal cells but had no effect on stomatal size. In contrast to other angiosperms, mangroves exhibited (1) a negative relationship between guard cell size and genome size; (2) epidermal cells that were smaller than stomata; and (3) coordination between Dv and Ds that was not mediated by epidermal cell size. Furthermore, mangrove epidermal cell sizes and packing densities covaried with leaf size.

CONCLUSIONS

While mangroves exhibited coordination between veins and stomata and attained a maximum theoretical stomatal conductance similar to that of other angiosperms, the tissue-level tradeoffs underlying these similar relationships across species and environments were markedly different, perhaps indicative of the unique structural and physiological adaptations of mangroves to their stressful environments.

In silico analysis of NAC gene family in the mangrove plant Avicennia marina provides clues for adaptation to intertidal habitats

NAC (NAM, ATAF1/2, CUC2) transcription factors (TFs) constitute a plant-specific gene family. It is reported that NAC TFs play important roles in plant growth and developmental processes and in response to biotic/abiotic stresses. Nevertheless, little information is known about the functional and evolutionary characteristics of NAC TFs in mangrove plants, a group of species adapting coastal intertidal habitats. Thus, we conducted a comprehensive investigation for NAC TFs in Avicennia marina, one pioneer species of mangrove plants. We totally identified 142 NAC TFs from the genome of A. marina. Combined with NAC proteins having been functionally characterized in other organisms, we built a phylogenetic tree to infer the function of NAC TFs in A. marina. Gene structure and motif sequence analyses suggest the sequence conservation and transcription regulatory regions-mediated functional diversity. Whole-genome duplication serves as the driver force to the evolution of NAC gene family. Moreover, two pairs of NAC genes were identified as positively selected genes of which AmNAC010/040 may be imposed on less constraint toward neofunctionalization. Quite a few stress/hormone-related responsive elements were found in promoter regions indicating potential response to various external factors. Transcriptome data revealed some NAC TFs were involved in pneumatophore and leaf salt gland development and response to salt, flooding and Cd stresses. Gene co-expression analysis found a few NAC TFs participates in the special biological processes concerned with adaptation to intertidal environment. In summary, this study provides detailed functional and evolutionary information about NAC gene family in mangrove plant A. marina and new perspective for adaptation to intertidal habitats.

Role of exogenous abscisic acid in freezing tolerance of mangrove Kandelia obovata under natural frost condition at near 32°N

BACKGROUND

Mangroves possess substantial ecological, social, and economic functions in tropical and subtropical coastal wetlands. Kandelia obovata is the most cold-resistance species among mangrove plants, with a widespread distribution in China that ranges from Sanya (18° 12' N) to Wenzhou (28° 20' N). Here, we explored the temporal variations in physiological status and transcriptome profiling of K. obovata under natural frost conditions at ~ 32^oN, as well as the positive role of exogenous abscisic acid (ABA) in cold resistance.

RESULTS

The soluble sugar (SS) and proline (Pro) functioned under freezing stress, of which SS was more important for K. obovata. Consistently, up-regulated DEGs responding to low temperature were significantly annotated to glycometabolism, such as starch and sucrose metabolism and amino sugar and nucleotide sugar metabolism. Notably, the top 2 pathways of KEGG enrichment were phenylpropanoid biosynthesis and flavonoid biosynthesis. For the antioxidant system, POD in conjunction with CAT removed hydrogen peroxide, and CAT appeared to be more important. The up-regulated DEGs responding to low temperature and ABA were also found to be enriched in arginine and proline metabolism, starch and sucrose metabolism, and peroxisome. Moreover, ABA triggered the expression of P5CS and P5CR, but inhibited the ProDH expression, which might contribute to Pro accumulation. Interestingly, there was no significant change in malondialdehyde (MDA) content during the cold event (P > 0.05), suggesting foliar application of ABA effectively alleviated the adverse effects of freezing stress on K. obovata by activating the antioxidant enzyme activity and increasing osmolytes accumulation, such as Pro, and the outcome was proportional to ABA concentration.

CONCLUSIONS

This study deepened our understanding of the physiological characters and molecular mechanisms underlying the response of K. obovata to natural frost conditions and exogenous ABA at the field level, which could provide a sound theoretical foundation for expanding mangroves plantations in higher latitudes, as well as the development coastal landscape.

Comprehensive analysis of the NAC transcription factor gene family in Kandelia obovata reveals potential members related to chilling tolerance

BACKGROUND

Kandelia obovata is an important mangrove species extensively distributed in Eastern Asia that is susceptible to low-temperature stress. NAC (NAM, ATAF1/2 and CUC2) domain proteins are transcription factors (TFs) that play various roles in plant growth and development and in the plant response to environmental stresses. Nevertheless, genome-wide analyses of K. obovata NAC genes (KoNACs) and their responses to chilling stress have rarely been studied.

METHODS

The KoNAC gene family was identified and characterized using bioinformatic analysis, the subcellular location of some NAC proteins was confirmed using confocal microscopy analysis, and the KoNACs that responded to chilling stress were screened using RNA-seq and qRT-PCR analysis.

RESULTS

A total of 79 KoNACs were identified, and they were unequally distributed across all 18 chromosomes of K. obovata. The KoNAC proteins could be divided into 16 subgroups according to the phylogenetic tree based on NAC family members of Arabidopsis thaliana. The KoNACs exhibited greater synteny with A. thaliana sequences than with Oryza sativa sequences, indicating that KoNACs underwent extensive evolution after the divergence of dicotyledons and monocotyledons. Segmental duplication was the main driving force of the expansions of KoNAC genes. Confocal microscopy analysis verified that the four randomly selected KoNACs localized to the nucleus, indicating the accuracy of the bioinformatic predictions. Tissue expression pattern analysis demonstrated that some KoNAC genes showed tissue-specific expression, suggesting that these KoNACs might be important for plant development and growth. Additionally, the expression levels of 19 KoNACs were significantly (15 positively and 4 negatively) induced by cold treatment, demonstrating that these KoNACs might play important roles during cold stress responses and might be candidate genes for the genetic engineering of K. obovata with enhanced chilling stress tolerance. Coexpression network analysis revealed that 381 coexpressed pairs (between 13 KoNACs and 284 other genes) were significantly correlated.

CONCLUSIONS

Seventy-nine KoNACs were identified in K. obovata, nineteen of which displayed chilling-induced expression patterns. These genes may serve as candidates for functional analyses of KoNACs engaged in chilling stress. Our results lay the foundation for evolutionary analyses of KoNACs and their molecular mechanisms in response to environmental stress.

Genome-Wide Identification and Expression Analysis of the NAC Gene Family in Kandelia obovata, a Typical Mangrove Plant

The NAC (NAM, ATAF1/2, and CUC2) gene family, one of the largest transcription factor families in plants, acts as positive or negative regulators in plant response and adaption to various environmental stresses, including cold stress. Multiple reports on the functional characterization of NAC genes in Arabidopsis thaliana and other plants are available. However, the function of the NAC genes in the typical woody mangrove (Kandelia obovata) remains poorly understood. Here, a comprehensive analysis of NAC genes in K. obovata was performed with a pluri-disciplinary approach including bioinformatic and molecular analyses. We retrieved a contracted NAC family with 68 genes from the K. obovata genome, which were unevenly distributed in the chromosomes and classified into ten classes. These KoNAC genes were differentially and preferentially expressed in different organs, among which, twelve up-regulated and one down-regulated KoNAC genes were identified. Several stress-related cis-regulatory elements, such as LTR (low-temperature response), STRE (stress response element), ABRE (abscisic acid response element), and WUN (wound-responsive element), were identified in the promoter regions of these 13 KoNAC genes. The expression patterns of five selected KoNAC genes (KoNAC6, KoNAC15, KoNAC20, KoNAC38, and KoNAC51) were confirmed by qRT-PCR under cold treatment. These results strongly implied the putative important roles of KoNAC genes in response to chilling and other stresses. Collectively, our findings provide valuable information for further investigations on the function of KoNAC genes.

Catalytic innovation underlies independent recruitment of polyketide synthases in cocaine and hyoscyamine biosynthesis

Tropane alkaloids such as hyoscyamine and cocaine are of importance in medicinal uses. Only recently has the hyoscyamine biosynthetic machinery become complete. However, the cocaine biosynthesis pathway remains only partially elucidated. Here we characterize polyketide synthases required for generating 3-oxo-glutaric acid from malonyl-CoA in cocaine biosynthetic route. Structural analysis shows that these two polyketide synthases adopt distinctly different active site architecture to catalyze the same reaction as pyrrolidine ketide synthase in hyoscyamine biosynthesis, revealing an unusual parallel/convergent evolution of biochemical function in homologous enzymes. Further phylogenetic analysis suggests lineage-specific acquisition of polyketide synthases required for tropane alkaloid biosynthesis in Erythroxylaceae and Solanaceae species, respectively. Overall, our work elucidates not only a key unknown step in cocaine biosynthesis pathway but also, more importantly, structural and biochemical basis for independent recruitment of polyketide synthases in tropane alkaloid biosynthesis, thus broadening the understanding of conservation and innovation of biosynthetic catalysts.

Evolution of coastal forests based on a full set of mangrove genomes

Genomic studies are now poised to explore whole communities of species. The ~70 species of woody plants that anchor the coastal ecosystems of the tropics, collectively referred to as mangroves, are particularly suited to this exploration. In this study, we de novo sequenced the whole genomes of 32 mangroves, which we combined with other sequences of 30 additional species, comprising almost all mangroves globally. These community-wide genomic data will be valuable for ecology, evolution and biodiversity research. While the data revealed 27 independent origins of mangroves, the total phylogeny shows only modest increases in species number, even in coastal areas of active speciation, suggesting that mangrove extinction is common. A possible explanation for common extinction is the frequent sea-level rises and falls (SLRs and SLFs) documented in the geological record. Indeed, near-extinctions of species with extremely small population size (N) often happened during periods of rapid SLR, as revealed by the genome-wide heterozygosity of almost all mangroves. Reduction in N has possibly been further compounded by population fragmentation and the subsequent accumulation of deleterious mutations, thus pushing mangroves even closer to extinction. Crucially, the impact of the next SLR will be exacerbated by human encroachment into these mangrove habitats, potentially altering the ecosystems of tropical coasts irreversibly.

A de novo reference assembly of the yellow mangrove Ceriops zippeliana genome

Mangroves are of great ecological and economical importance, providing shelters for a wide range of species and nursery habitats for commercially important marine species. Ceriops zippeliana (yellow mangrove) belongs to Rhizophoraceae family and is commonly distributed in the tropical and subtropical coastal communities. In this study, we present a high-quality assembly of the C. zippeliana genome. We constructed an initial draft assembly of 240,139,412 bases with an N50 contig length of 564,761 bases using the 10x Genomics linked-read technology. This assembly was further scaffolded with RagTag using a chromosome-scale assembly of a closely related Ceriops species as a reference. The final assembly contained 243,228,612 bases with an N50 scaffold length of 10,559,178 Mb. The size of the final assembly was close to those estimated using DNA flow cytometry (248 Mb) and the k-mer distribution analysis (246 Mb). We predicted a total of 23,474 gene models and 21,724 protein-coding genes in the C. zippeliana genome, of which 16,002 were assigned gene ontology terms. We recovered 97.1% of the highly conserved orthologs based on the Benchmarking Universal Single-Copy Orthologs analysis. The phylogenetic analysis based on single-copy orthologous genes illustrated that C. zippeliana and Ceriops tagal diverged approximately 10.2 million years ago (MYA), and their last common ancestor and Kandelia obovata diverged approximately 29.9 MYA. The high-quality assembly of C. zippeliana presented in this work provides a useful genomic resource for studying mangroves' unique adaptations to stressful intertidal habitats and for developing sustainable mangrove forest restoration and conservation programs.

Genome-Wide Identification of WRKY Genes and Their Responses to Chilling Stress in Kandelia obovata

Background:Kandelia obovata, a dominant mangrove species, is widely distributed in tropical and subtropical areas. Low temperature is the major abiotic stress that seriously limits the survival and growth of mangroves. WRKY transcription factors (TFs) play vital roles in responses to biotic and abiotic stresses. However, genome-wide analysis of WRKY genes in K. obovata and their responses to chilling stress have not been reported.

Methods: Bioinformatic analysis was used to identify and characterize the K. obovata WRKY (KoWRKY) gene family, RNA-seq and qRT–PCR analyses were employed to screen KoWRKYs that respond to chilling stress.

Results: Sixty-four KoWRKYs were identified and they were unevenly distributed across all 18 K. obovata chromosomes. Many orthologous WRKY gene pairs were identified between Arabidopsis thaliana and K. obovata, showing high synteny between the two genomes. Segmental duplication events were found to be the major force driving the expansion for the KoWRKY gene family. Most of the KoWRKY genes contained several kinds of hormone- and stress-responsive cis-elements in their promoter. KoWRKY proteins belonged to three groups (I, II, III) according to their conserved WRKY domains and zinc-finger structure. Expression patterns derived from the RNA-seq and qRT–PCR analyses revealed that 9 KoWRKYs were significantly upregulated during chilling acclimation in the leaves. KEGG pathway enrichment analysis showed that the target genes of KoWRKYs were significantly involved in 11 pathways, and coexpression network analysis showed that 315 coexpressed pairs (KoWRKYs and mRNAs) were positively correlated.

Conclusion: Sixty-four KoWRKYs from the K. obovata genome were identified, 9 of which exhibited chilling stress-induced expression patterns. These genes represent candidates for future functional analysis of KoWRKYs involved in chilling stress related signaling pathways in K. obovata. Our results provide a basis for further analysis of KoWRKY genes to determine their functions and molecular mechanisms in K. obovata in response to chilling stress.

De novo genome assembly and in natura epigenomics reveal salinity-induced DNA methylation in the mangrove tree Bruguiera gymnorhiza

Mangroves are adapted to harsh environments, such as high ultraviolet (UV) light, low nutrition, and fluctuating salinity in coastal zones. However, little is known about the transcriptomic and epigenomic basis of the resilience of mangroves due to limited available genome resources. We performed a de novo genome assembly and in natura epigenome analyses of the mangrove Bruguiera gymnorhiza, one of the dominant mangrove species. We also performed the first genome-guided transcriptome assembly for mangrove species. The 309 Mb of the genome is predicted to encode 34 403 genes and has a repeat content of 48%. Depending on its growing environment, the natural B. gymnorhiza population showed drastic morphological changes associated with expression changes in thousands of genes. Moreover, high-salinity environments induced genome-wide DNA hypermethylation of transposable elements (TEs) in the B. gymnorhiza. DNA hypermethylation was concurrent with the transcriptional regulation of chromatin modifier genes, suggesting robust epigenome regulation of TEs in the B. gymnorhiza genome under high-salinity environments. The genome and epigenome data in this study provide novel insights into the epigenome regulation of mangroves and a better understanding of the adaptation of plants to fluctuating, harsh natural environments.

Comparative Analysis and Phylogenetic Relationships of Ceriops Species (Rhizophoraceae) and Avicennia lanata (Acanthaceae): Insight into the Chloroplast Genome Evolution between Middle and Seaward Zones of Mangrove Forests

Ceriops and Avicennia are true mangroves in the middle and seaward zones of mangrove forests, respectively. The chloroplast genomes of Ceriops decandra, Ceriops zippeliana, and Ceriops tagal were assembled into lengths of 166,650, 166,083 and 164,432 bp, respectively, whereas Avicennia lanata was 148,264 bp in length. The gene content and gene order are highly conserved among these species. The chloroplast genome contains 125 genes in A. lanata and 129 genes in Ceriops species. Three duplicate genes (rpl2, rpl23, and trnM-CAU) were found in the IR regions of the three Ceriops species, resulting in expansion of the IR regions. The rpl32 gene was lost in C. zippeliana, whereas the infA gene was present in A. lanata. Short repeats (<40 bp) and a lower number of SSRs were found in A. lanata but not in Ceriops species. The phylogenetic analysis supports that all Ceriops species are clustered in Rhizophoraceae and A. lanata is in Acanthaceae. In a search for genes under selective pressures of coastal environments, the rps7 gene was under positive selection compared with non-mangrove species. Finally, two specific primer sets were developed for species identification of the three Ceriops species. Thus, this finding provides insightful genetic information for evolutionary relationships and molecular markers in Ceriops and Avicennia species.

Mangrove Forests: Natural Laboratories for Studying Epigenetic and Climate Changes

In the adaptation market, plants cash in the changes of their DNA (either genetic or epigenetic) to purchase fitness. Fitness is negatively affected by stressful conditions caused by climate change and well-designed studies are required to investigate the fine-tuning cooperation of epigenetic and genetic changes in response to those stresses. Mangrove trees are promising model systems for studying climate change because the effects of environmental changes are already evident in their natural habitats where they are exposed to different salinity levels ranging from saltwater to freshwater. In addition, as mangrove species are known to have very low genetic diversity caused by their stressful living conditions, epigenetic variation is likely to be a vital source for them to respond to environmental changes. This mini review aims to provide an overview of available studies on epigenetic regulation and adaptation of mangroves.

De novo reference assembly of the upriver orange mangrove (Bruguiera sexangula) genome

Upriver orange mangrove (Bruguiera sexangula) is a member of the most mangrove-rich taxon (Rhizophoraceae family) and is commonly distributed in the intertidal zones in tropical and subtropical latitudes. In this study, we employed the 10× Genomics linked-read technology to obtain a preliminary de novo assembly of the B. sexangula genome, which was further scaffolded to a pseudomolecule level using the Bruguiera parviflora genome as a reference. The final assembly of the B. sexangula genome contained 260 Mb with an N50 scaffold length of 11,020,310 bases. The assembly comprised 18 pseudomolecules (corresponding to the haploid chromosome number in B. sexangula), covering 204,645,832 bases or 78.6% of the 260-Mb assembly. We predicted a total of 23,978 protein-coding sequences, 17,598 of which were associated with gene ontology terms. Our gene prediction recovered 96.6% of the highly conserved orthologs based on the Benchmarking Universal Single-Copy Orthologs (BUSCO) analysis. The chromosome-level assembly presented in this work provides a valuable genetic resource to help strengthen our understanding of mangroves’ physiological and morphological adaptations to the intertidal zones.

Genetic and molecular mechanisms underlying mangrove adaptations to intertidal environments

Mangroves are halophytic plants belonging to diverse angiosperm families that are adapted to highly stressful intertidal zones between land and sea. They are special, unique, and one of the most productive ecosystems that play enormous ecological roles and provide a large number of benefits to the coastal communities. To thrive under highly stressful conditions, mangroves have innovated several key morphological, anatomical, and physio-biochemical adaptations. The evolution of the unique adaptive modifications might have resulted from a host of genetic and molecular changes and to date we know little about the nature of these genetic and molecular changes. Although slow, new information has accumulated over the last few decades on the genetic and molecular regulation of the mangrove adaptations, a comprehensive review on it is not yet available. This review provides up-to-date consolidated information on the genetic, epigenetic, and molecular regulation of mangrove adaptive traits.

The First De Novo Transcriptome Assembly and Transcriptomic Dynamics of the Mangrove Tree Rhizophora stylosa Griff. (Rhizophoraceae)

Mangroves are salt-tolerant plant species that grow in coastal saline water and are adapted to harsh environmental conditions. In this study, we de novo assembled and functionally annotated the transcriptome of Rhizophora stylosa, the widely distributed mangrove from the largest mangrove family (Rhizophoraceae). The final transcriptome consists of 200,491 unigenes with an average length, and N50 of 912.7 and 1334 base pair, respectively. We then compared the genome-wide expression profiles between the two morphologically distinct natural populations of this species growing under different levels of salinity depending on their distance from the ocean. Among the 200,491 unigenes, 40,253 were identified as differentially expressed between the two populations, while 15,741 and 24,512 were up- and down-regulated, respectively. Functional annotation assigned thousands of upregulated genes in saline environment to the categories related to abiotic stresses such as response to salt-, osmotic-, and oxidative-stress. Validation of those genes may contribute to a better understanding of adaptation in mangroves species. This study reported, for the first time, the transcriptome of R. stylosa, and the dynamic of it in response to salt stress and provided a valuable resource for elucidation of the molecular mechanism underlying the salt stress response in mangroves and other plants that live under stress.

Chromosome-scale assembly of the Dendrobium chrysotoxum genome enhances the understanding of orchid evolution

As one of the largest families of angiosperms, the Orchidaceae family is diverse. Dendrobium represents the second largest genus of the Orchidaceae. However, an assembled high-quality genome of species in this genus is lacking. Here, we report a chromosome-scale reference genome of Dendrobium chrysotoxum, an important ornamental and medicinal orchid species. The assembled genome size of D. chrysotoxum was 1.37 Gb, with a contig N50 value of 1.54 Mb. Of the sequences, 95.75% were anchored to 19 pseudochromosomes. There were 30,044 genes predicted in the D. chrysotoxum genome. Two whole-genome polyploidization events occurred in D. chrysotoxum. In terms of the second event, whole-genome duplication (WGD) was also found to have occurred in other Orchidaceae members, which diverged mainly via gene loss immediately after the WGD event occurred; the first duplication was found to have occurred in most monocots (tau event). We identified sugar transporter (SWEET) gene family expansion, which might be related to the abundant medicinal compounds and fleshy stems of D. chrysotoxum. MADS-box genes were identified in D. chrysotoxum, as well as members of TPS and Hsp90 gene families, which are associated with resistance, which may contribute to the adaptive evolution of orchids. We also investigated the interplay among carotenoid, ABA, and ethylene biosynthesis in D. chrysotoxum to elucidate the regulatory mechanisms of the short flowering period of orchids with yellow flowers. The reference D. chrysotoxum genome will provide important insights for further research on medicinal active ingredients and breeding and enhances the understanding of orchid evolution.

A reference-grade genome identifies salt-tolerance genes from the salt-secreting mangrove species Avicennia marina

Water scarcity and salinity are major challenges facing agriculture today, which can be addressed by engineering plants to grow in the boundless seawater. Understanding the mangrove plants at the molecular level will be necessary for developing such highly salt-tolerant agricultural crops. With this objective, we sequenced the genome of a salt-secreting and extraordinarily salt-tolerant mangrove species, Avicennia marina, that grows optimally in 75% seawater and tolerates >250% seawater. Our reference-grade ~457 Mb genome contains 31 scaffolds corresponding to its chromosomes. We identified 31,477 protein-coding genes and a salinome consisting of 3246 salinity-responsive genes and homologs of 614 experimentally validated salinity tolerance genes. The salinome provides a strong foundation to understand the molecular mechanisms of salinity tolerance in plants and breeding crops suitable for seawater farming.

Chromosome-level assembly of the mangrove plant Aegiceras corniculatum genome generated through Illumina, PacBio and Hi-C sequencing technologies

Aegiceras corniculatum is a major mangrove plant species adapted to waterlogging and saline conditions, grows in the coastal intertidal zone of tropical and subtropical regions. Here, we present a chromosome-level genome assembly of A. corniculatum by incorporating PacBio long-read sequencing and Hi-C technology. The results showed that the PacBio draft genome size is 906.63 Mb. Hi-C scaffolding anchored 885.06 Mb contigs (97.62% of draft assembly) onto 24 pseudochromosomes. The contig N50 and scaffold N50 were 7.1 Mb and 37.74 Mb, respectively. Out of 40,727 protein-coding genes predicted in the study, 89% have functional annotations in public databases. We also showed that of the 603.93 Mb repetitive sequences predicted in the assembled genome, long terminal repeat retrotransposons constitute 41.52%. The genome evolution analysis showed that the A. corniculatum genome experienced two whole-genome duplication events and shared the ancient γ whole-genome triplication event. A comparative genomic analysis revealed an incidence of expansion in 1,488 gene families associated with essential metabolism and biosynthetic pathways, including photosynthesis, oxidative phosphorylation, phenylalanine, glyoxylate, dicarboxylate metabolism, and DNA replication, which probably constitute adaptation traits that allow the A. corniculatum to survive in the intertidal zone. Also, the systematic characterization of genes associated with flavonoid biosynthesis pathway and the AcNHX gene family conducted in this study will provide insight into the adaptation mechanism of A. corniculatum to intertidal environments. The high-quality genome reported here can provide historical insights into genomic transformations that support the survival of A. corniculatum under harsh intertidal habitats.

The Phoebe genome sheds light on the evolution of magnoliids

Lauraceae includes the genus Phoebe, and the family is linked to the evolution of magnoliids. We sequenced the genome of Phoebe bournei Nanmu. The assembled genome size was 989.19 Mb, with a contig N50 value of 2.05 Mb. A total of 28,198 protein-coding genes were annotated in P. bournei. Whole-genome duplication (WGD) analysis showed that Lauraceae has experienced two WGD events; the older WGD event occurred just before the divergence of Lauraceae and Magnoliales, and the more recent WGD was shared by all lineages of Lauraceae. The phylogenetic tree showed that magnoliids form a sister clade to monocots and eudicots. We also identified 63 MADS-box genes, including AGL12-like genes that may be related to the regulation of P. bournei roots and FIN219-like genes encoding GH3 proteins, which are involved in photomorphogenesis. SAUR50-like genes involved in light signal-mediated pedicel or stem development were also identified. Four ATMYB46- and three PtrEPSP-homologous genes related to lignin biosynthesis were identified. These genes may be associated with the formation of straight trunks in P. bournei. Overall, the P. bournei reference genome provides insight into the origin, evolution, and diversification of Phoebe and other magnoliids.