The Litsea genome and the evolution of the laurel family

Probing of plant transcriptomes reveals the hidden genetic diversity of the family Secoviridae

Secoviruses are single-stranded RNA viruses that infect plants. In the present study, we identified 61 putative novel secoviral genomes in various plant species by mining publicly available plant transcriptome data. These viral sequences represent the genomes of 13 monopartite and 48 bipartite secovirids. The genome sequences of 52 secovirids were coding-complete, and nine were partial. Except for small open reading frames (ORFs) determined in waikaviral genomes and RNA2 of torradoviruses, all of the recovered genomes/genome segments contained a large ORF encoding a polyprotein. Based on genome organization and phylogeny, all but three of the novel secoviruses were assigned to different genera. The genome organization of two identified waika-like viruses resembled that of the recently identified waika-like virus Triticum aestivum secovirus. Phylogenetic analysis revealed a pattern of host-virus co-evolution in a few waika- and waika-like viruses and increased phylogenetic diversity of nepoviruses. The study provides a basis for further investigation of the biological properties of these novel secoviruses.

The chromosome-level genomes of the herbal magnoliids Warburgia ugandensis and Saururus chinensis

Warburgia ugandensis and Saururus chinensis are two of the most important medicinal plants in magnoliids and are widely utilized in traditional Kenya and Chinese medicine, respectively. The absence of higher-quality reference genomes has hindered research on the medicinal compound biosynthesis mechanisms of these plants. We report the chromosome-level genome assemblies of W. ugandensis and S. chinensis, and generated 1.13 Gb and 0.53 Gb genomes from 74 and 27 scaffolds, respectively, using BGI-DIPSEQ, Nanopore, and Hi-C sequencing. The scaffold N50 lengths were 82.97 Mb and 48.53 Mb, and the assemblies were anchored to 14 and 11 chromosomes of W. ugandensis and S. chinensis, respectively. In total, 24,739 and 20,561 genes were annotated, and 98.5% and 98% of the BUSCO genes were fully represented, respectively. The chromosome-level genomes of W. ugandensis and S. chinensis will be valuable resources for understanding the genetics of these medicinal plants, studying the evolution of magnoliids and angiosperms and conserving plant genetic resources.

Chromosome-level genome assembly of the threatened resource plant Cinnamomum chago

Cinnamomum chago is a tree species endemic to Yunnan province, China, with potential economic value, phylogenetic importance, and conservation priority. We assembled the genome of C. chago using multiple sequencing technologies, resulting in a high-quality, chromosomal-level genome with annotation information. The assembled genome size is approximately 1.06 Gb, with a contig N50 length of 92.10 Mb. About 99.92% of the assembled sequences could be anchored to 12 pseudo-chromosomes, with only one gap, and 63.73% of the assembled genome consists of repeat sequences. In total, 30,497 genes were recognized according to annotation, including 28,681 protein-coding genes. This high-quality chromosome-level assembly and annotation of C. chago will assist us in the conservation and utilization of this valuable resource, while also providing crucial data for studying the evolutionary relationships within the Cinnamomum genus, offering opportunities for further research and exploration of its diverse applications.

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.

Alcohol dehydrogenases regulated by a MYB44 transcription factor underlie Lauraceae citral biosynthesis

Lineage-specific terpenoids have arisen throughout the evolution of land plants and are believed to play a role in interactions between plants and the environment. Species-specific gene clusters in plants have provided insight on the evolution of secondary metabolism. Lauraceae is an ecologically important plant family whose members are also of considerable economic value given their monoterpene contents. However, the gene cluster responsible for the biosynthesis of monoterpenes remains yet to be elucidated. Here, a Lauraceae-specific citral biosynthetic gene cluster (CGC) was identified and investigated using a multifaceted approach that combined phylogenetic, collinearity, and biochemical analyses. The CGC comprises MYB44 as a regulator and 2 alcohol dehydrogenases (ADHs) as modifying enzymes, which derived from species-specific tandem and proximal duplication events. Activity and substrate divergence of the ADHs has resulted in the fruit of mountain pepper (Litsea cubeba), a core Lauraceae species, consisting of more than 80% citral. In addition, MYB44 negatively regulates citral biosynthesis by directly binding to the promoters of the ADH-encoding genes. The aggregation of citral biosynthetic pathways suggests that they may form the basis of important characteristics that enhance adaptability. The findings of this study provide insights into the evolution of and the regulatory mechanisms involved in plant terpene biosynthesis.

The American Cherimoya Genome Reveals Insights into the Intra-Specific Divergence, the Evolution of Magnoliales, and a Putative Gene Cluster for Acetogenin Biosynthesis

Annona cherimola (cherimoya) is a species renowned for its delectable fruit and medicinal properties. In this study, we developed a chromosome-level genome assembly for the cherimoya 'Booth' cultivar from the United States. The genome assembly has a size of 794 Mb with a N50 = 97.59 Mb. The seven longest scaffolds account for 87.6% of the total genome length, which corresponds to the seven pseudo-chromosomes. A total of 45,272 protein-coding genes (≥30 aa) were predicted with 92.9% gene content completeness. No recent whole genome duplications were identified by an intra-genome collinearity analysis. Phylogenetic analysis supports that eudicots and magnoliids are more closely related to each other than to monocots. Moreover, the Magnoliales was found to be more closely related to the Laurales than the Piperales. Genome comparison revealed that the 'Booth' cultivar has 200 Mb less repeats than the Spanish cultivar 'Fino de Jete', despite their highly similar (>99%) genome sequence identity and collinearity. These two cultivars were diverged during the early Pleistocene (1.93 Mya), which suggests a different origin and domestication of the cherimoya. Terpene/terpenoid metabolism functions were found to be enriched in Magnoliales, while TNL (Toll/Interleukin-1-NBS-LRR) disease resistance gene has been lost in Magnoliales during evolution. We have also identified a gene cluster that is potentially responsible for the biosynthesis of acetogenins, a class of natural products found exclusively in Annonaceae. The cherimoya genome provides an invaluable resource for supporting characterization, conservation, and utilization of Annona genetic resources.

Anthropogenic disturbance driving population decline of a dominant tree in East Asia evergreen broadleaved forests over the last 11,000 years

Current biodiversity loss is generally considered to have been caused by anthropogenic disturbance, but it is unclear when anthropogenic activities began to affect biodiversity loss. One hypothesis suggests it began with the Industrial Revolution, whereas others propose that anthropogenic disturbance has been associated with biodiversity decline since the early Holocene. To test these hypotheses, we examined the unique vegetation of evergreen broadleaved forests (EBLFs) in East Asia, where humans have affected landscapes since the early Holocene. We adopted a genomic approach to infer the demographic history of a dominant tree (Litsea elongata) of EBLFs. We used Holocene temperature and anthropogenic disturbance factors to calculate the correlation between these variables and the historical effective population size of L. elongata with Spearman statistics and integrated the maximum-entropy niche model to determine the impact of climate change and anthropogenic disturbance on fluctuation in its effective population size. We identified 9 well-defined geographic clades for the populations of L. elongata. Based on the estimated historical population sizes of these clades, all the populations contracted, indicating persistent population decline over the last 11,000 years. Demographic history of L. elongata and human population change, change in cropland use, and change in irrigated rice area were significantly negatively correlated, whereas climate change in the Holocene was not correlated with demographic history. Our results support the early human impact hypothesis and provide comprehensive evidence that early anthropogenic disturbance may contribute to the current biodiversity crisis in East Asia.

The chromosome-scale genome of Magnolia sinica (Magnoliaceae) provides insights into the conservation of plant species with extremely small populations (PSESP)

Magnolia sinica (Magnoliaceae) is a highly threatened tree endemic to southeast Yunnan, China. In this study, we generated for the first time a high-quality chromosome-scale genome sequence from M. sinica, by combining Illumina and ONT data with Hi-C scaffolding methods. The final assembled genome size of M. sinica was 1.84 Gb, with a contig N50 of ca. 45 Mb and scaffold N50 of 92 Mb. Identified repeats constituted approximately 57% of the genome, and 43,473 protein-coding genes were predicted. Phylogenetic analysis shows that the magnolias form a sister clade with the eudicots and the order Ceratophyllales, while the monocots are sister to the other core angiosperms. In our study, a total of 21 individuals from the 5 remnant populations of M. sinica, as well as 22 specimens belonging to 8 related Magnoliaceae species, were resequenced. The results showed that M. sinica had higher genetic diversity (θw = 0.01126 and θπ = 0.01158) than other related species in the Magnoliaceae. However, population structure analysis suggested that the genetic differentiation among the 5 M. sinica populations was very low. Analyses of the demographic history of the species using different models consistently revealed that 2 bottleneck events occurred. The contemporary effective population size of M. sinica was estimated to be 10.9. The different patterns of genetic loads (inbreeding and numbers of deleterious mutations) suggested constructive strategies for the conservation of these 5 different populations of M. sinica. Overall, this high-quality genome will be a valuable genomic resource for conservation of M. sinica.

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.

The Association between BZIP Transcription Factors and Flower Development in Litsea cubeba

The basic leucine zipper (bZIP) family is one of the largest families of transcription factors among eukaryotic organisms. Members of the bZIP family play various roles in regulating the intricate process of flower development in plants. Litsea cubeba (Lour.) (family: Lauraceae) is an aromatic, dioecious plant used in China for a wide range of applications. However, no study to date has undertaken a comprehensive analysis of the bZIP gene family in L. cubeba. In this work, we identified 68 members of the bZIP gene family in L. cubeba and classified them into 12 subfamilies based on previous studies on Arabidopsis thaliana. Transcriptome data analysis revealed that multiple LcbZIP genes exhibit significantly high expression levels in the flowers of L. cubeba, while some also demonstrate distinct temporal specificity during L. cubeba flower development. In particular, some LcbZIP genes displayed specific and high expression levels during the stamen and pistil degradation process. Using differential gene expression analysis, weighted gene co-expression network analysis, and Gene Ontology enrichment analysis, we identified six candidate LcbZIP genes that potentially regulate stamen or pistil degradation during flower development. In summary, our findings provide a framework for future functional analysis of the LcbZIP gene family in L. cubeba and offer novel insights for investigating the mechanism underlying pistil and stamen degeneration in this plant.

Expression divergence of expansin genes drive the heteroblasty in Ceratopteris chingii

BACKGROUND

Sterile-fertile heteroblasty is a common phenomenon observed in ferns, where the leaf shape of a fern sporophyll, responsible for sporangium production, differs from that of a regular trophophyll. However, due to the large size and complexity of most fern genomes, the molecular mechanisms that regulate the formation of these functionally different heteroblasty have remained elusive. To shed light on these mechanisms, we generated a full-length transcriptome of Ceratopteris chingii with PacBio Iso-Seq from five tissue samples. By integrating Illumina-based sequencing short reads, we identified the genes exhibiting the most significant differential expression between sporophylls and trophophylls.

RESULTS

The long reads were assembled, resulting in a total of 24,024 gene models. The differential expressed genes between heteroblasty primarily involved reproduction and cell wall composition, with a particular focus on expansin genes. Reconstructing the phylogeny of expansin genes across 19 plant species, ranging from green algae to seed plants, we identified four ortholog groups for expansins. The observed high expression of expansin genes in the young sporophylls of C. chingii emphasizes their role in the development of heteroblastic leaves. Through gene coexpression analysis, we identified highly divergent expressions of expansin genes both within and between species.

CONCLUSIONS

The specific regulatory interactions and accompanying expression patterns of expansin genes are associated with variations in leaf shapes between sporophylls and trophophylls.

Main Habitat Factors Driving the Phenotypic Diversity of Litsea cubeba in China

Litsea cubeba (Lour.) Pers. is an important woody spice tree in southern China, and its fruit is a rich source of valuable essential oil. We surveyed and sampled L. cubeba germplasm resources from 36 provenances in nine Chinese provinces, and detected rich phenotypic diversity. The survey results showed that plants of SC-KJ, SC-HJ, and SC-LS provenance presented higher leaf area (LA); YN-SM and YN-XC plants had larger thousand-grain fresh weight (TFW); and HN-DX plants had the highest essential oil content (EOC). To explain the large differences in the phenotypes of L. cubeba among different habitats, we used Pearson's correlation analysis, multiple stepwise regression path analysis, and redundancy analysis to evaluate the phenotypic diversity of L. cubeba. It was found that compared to other traits, leaf and fruit traits had more significant geographical distributions, and that leaf phenotypes were correlated to fruit phenotypes. The results showed that elevation, latitude, longitude, total soil porosity (SP), soil bulk density (SBD), and average annual rainfall (AAR, mm) contributed significantly to the phenotypic diversity of L. cubeba. Geographical factors explained a higher percentage of variation in phenotypic diversity than did soil factors and climate factors. Plants of SC-KJ and HN-DX provenances could be important resources for domestication and breeding to develop new high-yielding varieties of this woody aromatic plant. This study describes significant phenotypic differences in L. cubeba related to adaptation to different environments, and provides a theoretical basis for the development of a breeding strategy and for optimizing L. cubeba cultivation.

A chromosome-level genome assembly provides insights into Cornus wilsoniana evolution, oil biosynthesis, and floral bud development

Cornus wilsoniana W. is a woody oil plant with high oil content and strong hypolipidemic effects, making it a valuable species for medicinal, landscaping, and ecological purposes in China. To advance genetic research on this species, we employed PacBio together with Hi-C data to create a draft genome assembly for C. wilsoniana. Based on an 11-chromosome anchored chromosome-level assembly, the estimated genome size was determined to be 843.51 Mb. The N50 contig size and N50 scaffold size were calculated to be 4.49 and 78.00 Mb, respectively. Furthermore, 30 474 protein-coding genes were annotated. Comparative genomics analysis revealed that C. wilsoniana diverged from its closest species ~12.46 million years ago (Mya). Furthermore, the divergence between Cornaceae and Nyssaceae occurred >62.22 Mya. We also found evidence of whole-genome duplication events and whole-genome triplication γ, occurring at ~44.90 and 115.86 Mya. We further inferred the origins of chromosomes, which sheds light on the complex evolutionary history of the karyotype of C. wilsoniana. Through transcriptional and metabolic analysis, we identified two FAD2 homologous genes that may play a crucial role in controlling the oleic to linoleic acid ratio. We further investigated the correlation between metabolites and genes and identified 33 MADS-TF homologous genes that may affect flower morphology in C. wilsoniana. Overall, this study lays the groundwork for future research aimed at identifying the genetic basis of crucial traits in C. wilsoniana.

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

BACKGROUND

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

RESULTS

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

CONCLUSIONS

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

Genome-wide characterization of long terminal repeat retrotransposons provides insights into trait evolution of four cucurbit species

Cucurbits are a diverse plant family that includes economically important crops, such as cucumber, watermelon, melon, and pumpkin. Knowledge of the roles that long terminal repeat retrotransposons (LTR-RTs) have played in diversification of cucurbit species is limited; to add to understanding of the roles of LTR-RTs, we assessed their distributions in four cucurbit species. We identified 381, 578, 1086, and 623 intact LTR-RTs in cucumber (Cucumis sativus L. var. sativus cv. Chinese Long), watermelon (Citrullus lanatus subsp. vulgaris cv. 97103), melon (Cucumis melo cv. DHL92), and Cucurbita (Cucurbita moschata var. Rifu), respectively. Among these LTR-RTs, the Ale clade of the Copia superfamily was the most abundant in all the four cucurbit species. Insertion time and copy number analysis revealed that an LTR-RT burst occurred approximately 2 million years ago in cucumber, watermelon, melon, and Cucurbita, and may have contributed to their genome size variation. Phylogenetic and nucleotide polymorphism analyses suggested that most LTR-RTs were formed after species diversification. Analysis of gene insertions by LTR-RTs revealed that the most frequent insertions were of Ale and Tekay and that genes related to dietary fiber synthesis were the most commonly affected by LTR-RTs in Cucurbita. These results increase our understanding of LTR-RTs and their roles in genome evolution and trait characterization in cucurbits.

Independent rediploidization masks shared whole genome duplication in the sturgeon-paddlefish ancestor

Whole genome duplication (WGD) is a dramatic evolutionary event generating many new genes and which may play a role in survival through mass extinctions. Paddlefish and sturgeon are sister lineages that both show genomic evidence for ancient WGD. Until now this has been interpreted as two independent WGD events due to a preponderance of duplicate genes with independent histories. Here we show that although there is indeed a plurality of apparently independent gene duplications, these derive from a shared genome duplication event occurring well over 200 million years ago, likely close to the Permian-Triassic mass extinction period. This was followed by a prolonged process of reversion to stable diploid inheritance (rediploidization), that may have promoted survival during the Triassic-Jurassic mass extinction. We show that the sharing of this WGD is masked by the fact that paddlefish and sturgeon lineage divergence occurred before rediploidization had proceeded even half-way. Thus, for most genes the resolution to diploidy was lineage-specific. Because genes are only truly duplicated once diploid inheritance is established, the paddlefish and sturgeon genomes are thus a mosaic of shared and non-shared gene duplications resulting from a shared genome duplication event.

Variations in Essential Oils from the Leaves of Cinnamomum bodinieri in China

Cinnamomum plants are rich in natural essential oils, which are widely used as materials in the fragrance, insecticidal, antibacterial agent, pharmaceutical, and food industries; however, few studies have investigated the essential oil components of Cinnamomum bodinieri. Therefore, this study investigated the diversity of essential oils from the leaves of 885 individual C. bodinieri plants across 32 populations in five provinces. Essential oils were extracted by hydrodistillation, and then qualitative and quantitative analyses of the compounds were performed by GC-MS and GC-FID. A total of 87 chemical constituents were identified in the essential oils, including 33 monoterpenes, 48 sesquiterpenes, and six other compounds. The average oil yield was 0.75%, and individual oil yields ranged from 0.01% to 4.28%. A total of 16 chemotypes were classified according to variations in the essential oil chemical constituents of C. bodinieri, among which the camphor-type, citral-type, and eucalyptol-type were dominant. Moreover, the borneol-type, cymol-type, elemol-type, methylisoeugenol-type, and selina-6-en-4-ol-type were reported in C. bodinieri for the first time. The yield and principal components of the essential oils were mainly affected by altitude, temperature, and sunshine duration, among which altitude had the most significant effect; thus, low-altitude areas are more suitable for the synthesis and accumulation of essential oils. Based on the different characteristics of the essential oils in the leaves of C. bodinieri, several excellent populations and individuals were identified in this study. Moreover, the findings provide a foundation for breeding superior varieties and studying essential oil biosynthesis mechanisms in the future.

LcWRKY17, a WRKY Transcription Factor from Litsea cubeba, Effectively Promotes Monoterpene Synthesis

The WRKY gene family is one of the most significant transcription factor (TF) families in higher plants and participates in many secondary metabolic processes in plants. Litsea cubeba (Lour.) Person is an important woody oil plant that is high in terpenoids. However, no studies have been conducted to investigate the WRKY TFs that regulate the synthesis of terpene in L. cubeba. This paper provides a comprehensive genomic analysis of the LcWRKYs. In the L. cubeba genome, 64 LcWRKY genes were discovered. According to a comparative phylogenetic study with Arabidopsis thaliana, these L. cubeba WRKYs were divided into three groups. Some LcWRKY genes may have arisen from gene duplication, but the majority of LcWRKY evolution has been driven by segmental duplication events. Based on transcriptome data, a consistent expression pattern of LcWRKY17 and terpene synthase LcTPS42 was found at different stages of L. cubeba fruit development. Furthermore, the function of LcWRKY17 was verified by subcellular localization and transient overexpression, and overexpression of LcWRKY17 promotes monoterpene synthesis. Meanwhile, dual-Luciferase and yeast one-hybrid (Y1H) experiments showed that the LcWRKY17 transcription factor binds to W-box motifs of LcTPS42 and enhances its transcription. In conclusion, this research provided a fundamental framework for future functional analysis of the WRKY gene families, as well as breeding improvement and the regulation of secondary metabolism in L. cubeba.

LcERF134 increases the production of monoterpenes by activating the terpene biosynthesis pathway in Litsea cubeba

Litsea cubeba, an aromatic species of the Lauraceae family, produces a diverse array of monoterpenes. The biosynthesis of monoterpenes is regulated by transcriptional factors (TFs), such as APETALA2/ethylene response factor (AP2/ERF). However, the regulatory mechanisms that control the AP2/ERF gene responsible for the biosynthesis of monoterpenes in L. cubeba have yet to be elucidated. Here, we identified an AP2/ERF gene, LcERF134, as an activator for the accumulation of citral and other monoterpenes. The expression level of LcERF134 was consistent with terpene synthase LcTPS42 in the pericarp. The transient overexpression of LcERF134 significantly increased monoterpene production in L. cubeba as well as the expression of rate-limiting genes involved in the monoterpene biosynthesis pathway. Furthermore, yeast one-hybrid, dual-luciferase and electrophoretic mobility shift assays demonstrated that LcERF134 activated the monoterpene biosynthesis pathway by directly binding to the GCC-box elements of the LcTPS42 and LcGPPS.SSU1 promoters. However, the overexpression of LcERF134 in tomatoes had no impact on the synthesis of monoterpenes, thus indicating that LcERF134 is a species-specific TF. Our research demonstrated that LcERF134 significantly increased the biosynthesis of monoterpenes by inducing the expression of LcTPS42 and LcGPPS.SSU1, thus offering insight into how to enhance the flavor of L. cubeba essential oil.

Phylogeny and Systematics of Sassafras (Lauraceae), an Interesting Genus with Disjunct Distributions in Eastern North America and East Asia

The Lauraceae is a family of the order Laurales, with 2500-3000 species comprising 50 genera, mainly distributed in tropical and subtropical evergreen broad-leaved forests. Up to two decades ago, the systematic classification of the Lauraceae was based on floral morphology, but molecular phylogenetic approaches have made tremendous advances in elucidating tribe- and genus-level relationships within the family in recent decades. Our review focused on the phylogeny and systematics of Sassafras, a genus of three species with highly disjunct distributions in eastern North America and East Asia, whose tribe affiliation within the Lauraceae has long been controversial. By synthesizing information on the floral biology and molecular phylogeny of Sassafras, this review aimed to explore the position of Sassafras within the Lauraceae, and to provide suggestions and implications for future phylogenetic studies. Our synthesis highlighted Sassafras as a transitional type between Cinnamomeae and Laureae with a closer genetic relationship with Cinnamomeae, as revealed by molecular phylogenetic evidence, while it shares many similar characteristics with Laureae in morphology. We therefore discovered that several molecular and morphological methods should be concurrently considered to illuminate the phylogeny and systematics of Sassafras in Lauraceae.

Unique gene duplications and conserved microsynteny potentially associated with resistance to wood decay in the Lauraceae

Wood decay resistance (WDR) is marking the value of wood utilization. Many trees of the Lauraceae have exceptional WDR, as evidenced by their use in ancient royal palace buildings in China. However, the genetics of WDR remain elusive. Here, through comparative genomics, we revealed the unique characteristics related to the high WDR in Lauraceae trees. We present a 1.27-Gb chromosome-level assembly for Lindera megaphylla (Lauraceae). Comparative genomics integrating major groups of angiosperm revealed Lauraceae species have extensively shared gene microsynteny associated with the biosynthesis of specialized metabolites such as isoquinoline alkaloids, flavonoid, lignins and terpenoid, which play significant roles in WDR. In Lauraceae genomes, tandem and proximal duplications (TD/PD) significantly expanded the coding space of key enzymes of biosynthesis pathways related to WDR, which may enhance the decay resistance of wood by increasing the accumulation of these compounds. Among Lauraceae species, genes of WDR-related biosynthesis pathways showed remarkable expansion by TD/PD and conveyed unique and conserved motifs in their promoter and protein sequences, suggesting conserved gene collinearity, gene expansion and gene regulation supporting the high WDR. Our study thus reveals genomic profiles related to biochemical transitions among major plant groups and the genomic basis of WDR in the Lauraceae.

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

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

The genome of Magnolia hypoleuca provides a new insight into cold tolerance and the evolutionary position of magnoliids

Magnolia hypoleuca Sieb. & Zucc, a member of the Magnoliaceae of magnoliids, is one of the most economically valuable, phylogenetic and ornamental tree species in Eastern China. Here, the 1.64 Gb chromosome-level assembly covers 96.64% of the genome which is anchored to 19 chromosomes, with a contig N50 value of 1.71 Mb and 33,873 protein-coding genes was predicted. Phylogenetic analyses between M. hypoleuca and other 10 representative angiosperms suggested that magnoliids were placed as a sister group to the eudicots, rather than sister to monocots or both monocots and eudicots. In addition, the relative timing of the whole-genome duplication (WGD) events about 115.32 Mya for magnoliid plants. M. hypoleuca was found to have a common ancestor with M. officinalis approximately 23.4 MYA, and the climate change of OMT (Oligocene-Miocene transition) is the main reason for the divergence of M. hypoleuca and M. officinalis, which was along with the division of Japanese islands. Moreover, the TPS gene expansion observed in M. hypoleuca might contribute to the enhancement of flower fragrance. Tandem and proximal duplicates of younger age that have been preserved have experienced more rapid sequence divergence and a more clustered distribution on chromosomes contributing to fragrance accumulation, especially phenylpropanoid, monoterpenes and sesquiterpenes and cold tolerance. The stronger selective pressure drived the evolution of tandem and proximal duplicates toward plant self-defense and adaptation. The reference M. hypoleuca genome will provide insights into the evolutionary process of M. hypoleuca and the relationships between the magnoliids with monocots and eudicots, and enable us to delve into the fragrance and cold tolerance produced by M. hypoleuca and provide more robust and deep insight of how the Magnoliales evolved and diversified.

Revisiting ancient polyploidy in leptosporangiate ferns

Ferns, and particularly homosporous ferns, have long been assumed to have experienced recurrent whole-genome duplication (WGD) events because of their substantially large genome sizes, surprisingly high chromosome numbers, and high degrees of polyploidy among many extant members. As the number of sequenced fern genomes is limited, recent studies have employed transcriptome data to find evidence for WGDs in ferns. However, they have reached conflicting results concerning the occurrence of ancient polyploidy, for instance, in the lineage of leptosporangiate ferns. Because identifying WGDs in a phylogenetic context is the foremost step in studying the contribution of ancient polyploidy to evolution, we here revisited earlier identified WGDs in leptosporangiate ferns, mainly the core leptosporangiate ferns, by building K_S -age distributions and applying substitution rate corrections and by conducting statistical gene tree-species tree reconciliation analyses. Our integrative analyses not only identified four ancient WGDs in the sampled core leptosporangiate ferns but also identified false positives and false negatives for WGDs that recent studies have reported earlier. In conclusion, we underscore the significance of substitution rate corrections and uncertainties in gene tree-species tree reconciliations in calling WGD events and advance an exemplar workflow to overcome such often-overlooked issues.

Genetic analysis and elite tree selection of the main resin components of slash pine

Pine resin, as a natural material, has been widely used in food, pharmaceutical, and chemical industries. Slash pine (Pinus elliottii Engelm var. elliottii) is the primary tree species for resin tapping due to its high resin yield, low resin crystallization rate, and high turpentine content. Current researches focuse on the targeted improvement of several significant components to meet industrial needs rather than just resin yield. The objective of this study was to examine the genetic variation and correlation of genetic and phenotype for four main resin components (α pinene, β pinene, abietic acid, and levoprimaric acid) of 219 half-sib progenies from 59 families. The results showed that the levopimaric acid had the largest content (mean value = 21.63%), while the β pinene content had the largest variation coefficient (CV = 0.42). The α pinene content has the highest heritability (h² = 0.67), while levopimaric acid has the lowest heritability (h² = 0.51). There was a significant negative correlation between α pinene and the other three components and a significant positive correlation between β pinene and the two diterpenes. The family ranking and genetic gain suggested that it is possible to improve the contents of main resin components of slash pine through genetic breeding selection.

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.

The Current Developments in Medicinal Plant Genomics Enabled the Diversification of Secondary Metabolites' Biosynthesis

Medicinal plants produce important substrates for their adaptation and defenses against environmental factors and, at the same time, are used for traditional medicine and industrial additives. Plants have relatively little in the way of secondary metabolites via biosynthesis. Recently, the whole-genome sequencing of medicinal plants and the identification of secondary metabolite production were revolutionized by the rapid development and cheap cost of sequencing technology. Advances in functional genomics, such as transcriptomics, proteomics, and metabolomics, pave the way for discoveries in secondary metabolites and related key genes. The multi-omics approaches can offer tremendous insight into the variety, distribution, and development of biosynthetic gene clusters (BGCs). Although many reviews have reported on the plant and medicinal plant genome, chemistry, and pharmacology, there is no review giving a comprehensive report about the medicinal plant genome and multi-omics approaches to study the biosynthesis pathway of secondary metabolites. Here, we introduce the medicinal plant genome and the application of multi-omics tools for identifying genes related to the biosynthesis pathway of secondary metabolites. Moreover, we explore comparative genomics and polyploidy for gene family analysis in medicinal plants. This study promotes medicinal plant genomics, which contributes to the biosynthesis and screening of plant substrates and plant-based drugs and prompts the research efficiency of traditional medicine.

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

Terpenoids, including aromatic volatile monoterpenoids and sesquiterpenoids, function in defense against pathogens and herbivores. Phoebe trees are remarkable for their scented wood and decay resistance. Unlike other Lauraceae species investigated to date, Phoebe species predominantly accumulate sesquiterpenoids instead of monoterpenoids. Limited genomic data restrict the elucidation of terpenoid variation and functions. Here, we present a chromosome-scale genome assembly of a Lauraceae tree, Phoebe bournei, and identify 72 full-length terpene synthase (TPS) genes. Genome-level comparison shows pervasive lineage-specific duplication and contraction of TPS subfamilies, which have contributed to the extreme terpenoid variation within Lauraceae species. Although the TPS-a and TPS-b subfamilies were both expanded via tandem duplication in P. bournei, more TPS-a copies were retained and constitutively expressed, whereas more TPS-b copies were lost. The TPS-a genes on chromosome 8 functionally diverged to synthesize eight highly accumulated sesquiterpenes in P. bournei. The essential oil of P. bournei and its main component, β-caryophyllene, exhibited antifungal activities against the three most widespread canker pathogens of trees. The TPS-a and TPS-b subfamilies have experienced contrasting fates over the evolution of P. bournei. The abundant sesquiterpenoids produced by TPS-a proteins contribute to the excellent pathogen resistance of P. bournei trees. Overall, this study sheds light on the evolution and adaptation of terpenoids in Lauraceae and provides valuable resources for boosting plant immunity against pathogens in various trees and crops.

Chromosome-level genome assembly and functional characterization of terpene synthases provide insights into the volatile terpenoid biosynthesis of Wurfbainia villosa

Wurfbainia villosa is a well-known medicinal and edible plant that is widely cultivated in the Lingnan region of China. Its dried fruits (called Fructus Amomi) are broadly used in traditional Chinese medicine for curing gastrointestinal diseases and are rich in volatile terpenoids. Here, we report a high-quality chromosome-level genome assembly of W. villosa with a total size of approximately 2.80 Gb, 42 588 protein-coding genes, and a very high percentage of repetitive sequences (87.23%). Genome analysis showed that W. villosa likely experienced a recent whole-genome duplication event prior to the W. villosa-Zingiber officinale divergence (approximately 11 million years ago), and a recent burst of long terminal repeat insertions afterward. The W. villosa genome enabled the identification of 17 genes involved in the terpenoid skeleton biosynthesis pathway and 66 terpene synthase (TPS) genes. We found that tandem duplication events have an important contribution to the expansion of WvTPSs, which likely drove the production of volatile terpenoids. In addition, functional characterization of 18 WvTPSs, focusing on the TPS-a and TPS-b subfamilies, showed that most of these WvTPSs are multi-product TPS and are predominantly expressed in seeds. The present study provides insights into the genome evolution and the molecular basis of the volatile terpenoids diversity in W. villosa. The genome sequence also represents valuable resources for the functional gene research and molecular breeding of W. villosa.

Chromosome-level genome assembly and resequencing of camphor tree (Cinnamomum camphora) provides insight into phylogeny and diversification of terpenoid and triglyceride biosynthesis of Cinnamomum

Cinnamomum species attract attentions owing to their scents, medicinal properties, and ambiguous relationship in the phylogenetic tree. Here, we report a high-quality genome assembly of Cinnamomum camphora, based on which two whole-genome duplication (WGD) events were detected in the C. camphora genome: one was shared with Magnoliales, and the other was unique to Lauraceae. Phylogenetic analyses illustrated that Lauraceae species formed a compact sister clade to the eudicots. We then performed whole-genome resequencing on 24 Cinnamomum species native to China, and the results showed that the topology of Cinnamomum species was not entirely consistent with morphological classification. The rise and molecular basis of chemodiversity in Cinnamomum were also fascinating issues. In this study, six chemotypes were classified and six main terpenoids were identified as major contributors of chemodiversity in C. camphora by the principal component analysis. Through in vitro assays and subcellular localization analyses, we identified two key terpene synthase (TPS) genes (CcTPS16 and CcTPS54), the products of which were characterized to catalyze the biosynthesis of two uppermost volatiles (i.e. 1,8-cineole and (iso)nerolidol), respectively, and meditate the generation of two chemotypes by transcriptional regulation and compartmentalization. Additionally, the pathway of medium-chain triglyceride (MCT) biosynthesis in Lauraceae was investigated for the first time. Synteny analysis suggested that the divergent synthesis of MCT and long-chain triglyceride (LCT) in Lauraceae kernels was probably controlled by specific medium-chain fatty acyl-ACP thioesterase (FatB), type-B lysophosphatidic acid acyltransferase (type-B LPAAT), and diacylglycerol acyltransferase 2b (DGAT 2b) isoforms during co-evolution with retentions or deletions in the genome.

Comparative Analysis of Complete Chloroplast Genomes of Nine Species of Litsea (Lauraceae): Hypervariable Regions, Positive Selection, and Phylogenetic Relationships

Litsea is a group of evergreen trees or shrubs in the laurel family, Lauraceae. Species of the genus are widely used for a wide range of medicinal and industrial aspects. At present, most studies related to the gene resources of Litsea are restricted to morphological analyses or features of individual genomes, and currently available studies of select molecular markers are insufficient. In this study, we assembled and annotated the complete chloroplast genomes of nine species in Litsea, carried out a series of comparative analyses, and reconstructed phylogenetic relationships within the genus. The genome length ranged from 152,051 to 152,747 bp and a total of 128 genes were identified. High consistency patterns of codon bias, repeats, divergent analysis, single nucleotide polymorphisms (SNP) and insertions and deletions (InDels) were discovered across the genus. Variations in gene length and the presence of the pseudogene ycf1^Ψ, resulting from IR contraction and expansion, are reported. The hyper-variable gene rpl16 was identified for its exceptionally high Ka/Ks and Pi values, implying that those frequent mutations occurred as a result of positive selection. Phylogenetic relationships were recovered for the genus based on analyses of full chloroplast genomes and protein-coding genes. Overall, both genome sequences and potential molecular markers provided in this study enrich the available genomic resources for species of Litsea. Valuable genomic resources and divergent analysis are also provided for further research of the evolutionary patterns, molecular markers, and deeper phylogenetic relationships of Litsea.

Genomic insights into the evolution of plant chemical defense

Plant trait evolution can be impacted by common mechanisms of genome evolution, including whole-genome and small-scale duplication, rearrangement, and selective pressures. With the increasing accessibility of genome sequencing for non-model species, comparative studies of trait evolution among closely related or divergent lineages have supported investigations into plant chemical defense. Plant defensive compounds include major chemical classes, such as terpenoids, alkaloids, and phenolics, and are used in primary and secondary plant functions. These include the promotion of plant health, facilitation of pollination, defense against pathogens, and responses to a rapidly changing climate. We discuss mechanisms of genome evolution and use examples from recent studies to impress a stronger understanding of the link between genotype and phenotype as it relates to the evolution of plant chemical defense. We conclude with considerations for how to leverage genomics, transcriptomics, metabolomics, and functional assays for studying the emergence and evolution of chemical defense systems.

Genome of Lindera glauca provides insights into the evolution of biosynthesis genes for aromatic compounds

Lindera glauca is a crucial source of diverse industrial oil and medicines. The spicy aroma of tender leaves is caused by the presence of abundant aromatic compounds. Here, we present its chromosome-level genome assembly comprising 12 pseudochromosomes (2,092.2 Mb; scaffold N50: 186.5 Mb), which was predicted to have 65,145 protein-coding genes. Comparative genomic analyses indicated two whole-genome duplication (WGD) events in the Lauraceae family, contributing to the production of numerous terpene synthase (TPS) genes. We identified 138 TPS genes in L. glauca. Comparative transcriptomic analyses revealed high expression of genes Lg03G2346 and Lg08G140 in TPS-a and Lg07G2961 and Lg12G971 in TPS-b subfamilies, which regulated the biosynthesis of the monoterpenoid β-ocimene and sesquiterpenoid D-germacrene in L. glauca. The results suggested a molecular basis for species-specific terpenoid biosynthesis and provided a foundation for molecular breeding to produce desired characteristics and a valuable reference genome.

•

We provide the first chromosome-level genome for Lindera glauca

•

We explore the phylogenetic position and identify three WGD events of L. glauca

•

We identify genes involved in the main aromatic compounds

•

Analysis of metabolites in fruits and leaves by GC-MS is reported

Pervasive genome duplications across the plant tree of life and their links to major evolutionary innovations and transitions

Whole-genome duplication (WGD) has occurred repeatedly during plant evolution and diversification, providing genetic layers for evolving new functions and phenotypes. Advances in long-read sequencing technologies have enabled sequencing and assembly of over 1000 plant genomes spanning nearly 800 species, in which a large set of ancient WGDs has been uncovered. Here, we review the recently reported WGDs that occurred in major plant lineages and key evolutionary positions, and highlight their contributions to morphological innovation and adaptive evolution. Current gaps and challenges in integrating enormous volumes of sequenced plant genomes, accurately inferring WGDs, and developing web-based analysis tools are emphasized. Looking to the future, ambitious genome sequencing projects and global efforts may substantially recapitulate the plant tree of life based on broader sampling of phylogenetic diversity, reveal much of the timetable of ancient WGDs, and address the biological significance of WGDs in plant adaptation and radiation.

High-quality chromosome-level genome assembly of Litsea coreana L. provides insights into Magnoliids evolution and flavonoid biosynthesis

The magnoliid Litsea coreana has been the subject of a substantial amount of research owing to its production of many flavonoid metabolites, high food processing value, and a controversial phylogenetic position. For this study, we assembled a high-grade genome at the chromosome scale and annotation of L. coreana that was anchored to 12 chromosomes. The total genome was 1139.45 Mb, while the N50 scaffold was 97.18 Mb long. The analysis of phylogenetic trees constructed by different methods show that the phylogeny of Magnoliids is inconsistent, indicating that the differentiation process of monocots, eudicots, and Magnoliids still remains in dispute. An ancient whole-genome duplication (WGD) event was shown to have occurred before the Magnoliales and Laurels had differentiated. Subsequently, an independent WGD appeared in the Lauralean lineage. A total of 27 types of flavonoids were detected in all five tissues of L. coreana. Chalcone synthases (CHSs) that are responsible for production of flavonoids have been validated at the bioinformatics level. The retention of comparative genomic analyses of the CHS gene family showed that this family had contracted significantly in L. coreana. Our research further elaborated the evolution of Lauraceae and perfected the genetic basis of flavonoid biosynthesis in L. coreana. SIGNIFICANCE STATEMENT: Provides evidence that determines the evolutionary status of Magnoliids. The chalcone synthase gene family was significantly contracted in Litsea coreana.

LcERF19, an AP2/ERF transcription factor from Litsea cubeba, positively regulates geranial and neral biosynthesis

The APETALA2/ETHYLENE RESPONSE FACTOR (AP2/ERF) transcription factors (TFs) are involved in the regulation of specialized terpenoid biosynthesis. However, the AP2/ERF TFs in Litsea cubeba have not been characterized and their role in the biosynthesis of terpenoids is unknown. Here, 174 LcAP2/ERF TFs were identified in L. cubeba and categorized into four subfamilies: 27 AP2, 7 RAV, 1 Soloist, and 139 ERF. Transcriptomic and qRT-PCR assays both showed that the expression levels of LcERF19 were similar to that of terpene synthase LcTPS42 in the pericarp, which is related to the synthesis of geranial and neral in L. cubeba. LcERF19 was further shown to encode a nuclear-localized protein and its expression was strongly induced by jasmonate. Yeast one-hybrid and dual-luciferase assays showed that LcERF19 associated with GCC box elements of the LcTPS42 promoter and promoted its activity. Transient overexpression of LcERF19 in L. cubeba and overexpression of LcERF19 in tomato resulted in a significant increase in geranial and neral. Our findings show that LcERF19 enhances geranial and neral biosynthesis through activation of LcTPS42 expression, which provides a strategy to improve the flavor of tomato and other fruits.

A Chromosome-Level Genome of the Camphor Tree and the Underlying Genetic and Climatic Factors for Its Top-Geoherbalism

Camphor tree [Cinnamomum camphora (L.) J. Presl], a species in the magnoliid family Lauraceae, is known for its rich volatile oils and is used as a medical cardiotonic and as a scent in many perfumed hygiene products. Here, we present a high-quality chromosome-scale genome of C. camphora with a scaffold N50 of 64.34 Mb and an assembled genome size of 755.41 Mb. Phylogenetic inference revealed that the magnoliids are a sister group to the clade of eudicots and monocots. Comparative genomic analyses identified two rounds of ancient whole-genome duplication (WGD). Tandem duplicated genes exhibited a higher evolutionary rate, a more recent evolutionary history and a more clustered distribution on chromosomes, contributing to the production of secondary metabolites, especially monoterpenes and sesquiterpenes, which are the principal essential oil components. Three-dimensional analyses of the volatile metabolites, gene expression and climate data of samples with the same genotype grown in different locations showed that low temperature and low precipitation during the cold season modulate the expression of genes in the terpenoid biosynthesis pathways, especially TPS genes, which facilitates the accumulation of volatile compounds. Our study lays a theoretical foundation for policy-making regarding the agroforestry applications of camphor tree.

Genomes of leafy and leafless Platanthera orchids illuminate the evolution of mycoheterotrophy

To improve our understanding of the origin and evolution of mycoheterotrophic plants, we here present the chromosome-scale genome assemblies of two sibling orchid species: partially mycoheterotrophic Platanthera zijinensis and holomycoheterotrophic Platanthera guangdongensis. Comparative analysis shows that mycoheterotrophy is associated with increased substitution rates and gene loss, and the deletion of most photoreceptor genes and auxin transporter genes might be linked to the unique phenotypes of fully mycoheterotrophic orchids. Conversely, trehalase genes that catalyse the conversion of trehalose into glucose have expanded in most sequenced orchids, in line with the fact that the germination of orchid non-endosperm seeds needs carbohydrates from fungi during the protocorm stage. We further show that the mature plant of P. guangdongensis, different from photosynthetic orchids, keeps expressing trehalase genes to hijack trehalose from fungi. Therefore, we propose that mycoheterotrophy in mature orchids is a continuation of the protocorm stage by sustaining the expression of trehalase genes. Our results shed light on the molecular mechanism underlying initial, partial and full mycoheterotrophy.

A haplotype resolved chromosomal level avocado genome allows analysis of novel avocado genes

Avocado (Persea americana) is a member of the magnoliids, an early branching lineage of angiosperms that has high value globally with the fruit being highly nutritious. Here, we report a chromosome-level genome assembly for the commercial avocado cultivar Hass, which represents 80% of the world's avocado consumption. The DNA contigs produced from Pacific Biosciences HiFi reads were further assembled using a previously published version of the genome supported by a genetic map. The total assembly was 913 Mb with a contig N50 of 84 Mb. Contigs assigned to the 12 chromosomes represented 874 Mb and covered 98.8% of benchmarked single-copy genes from embryophytes. Annotation of protein coding sequences identified 48 915 avocado genes of which 39 207 could be ascribed functions. The genome contained 62.6% repeat elements. Specific biosynthetic pathways of interest in the genome were investigated. The analysis suggested that the predominant pathway of heptose biosynthesis in avocado may be through sedoheptulose 1,7 bisphosphate rather than via alternative routes. Endoglucanase genes were high in number, consistent with avocado using cellulase for fruit ripening. The avocado genome appeared to have a limited number of translocations between homeologous chromosomes, despite having undergone multiple genome duplication events. Proteome clustering with related species permitted identification of genes unique to avocado and other members of the Lauraceae family, as well as genes unique to species diverged near or prior to the divergence of monocots and eudicots. This genome provides a tool to support future advances in the development of elite avocado varieties with higher yields and fruit quality.

Antibacterial Activity of the Essential Oil From Litsea cubeba Against Cutibacterium acnes and the Investigations of Its Potential Mechanism by Gas Chromatography-Mass Spectrometry Metabolomics

Cutibacterium acnes (C. acnes) is an anaerobic Gram-positive bacterium generally considered as a human skin commensal, but is also involved in different infections, such as acne and surgical infections. Although there are a variety of treatments, the side effects and the problem of bacterial drug resistance still limit their clinical usage. In this study, we found that essential oil (EO) distilled from fresh mature Litsea cubeba possessed promising antibacterial activity against C. acnes. In order to elucidate its potential mechanism, bacteriostatic activity test, Live/Dead kit assay, scanning electron microscope (SEM), transmission electron microscope (TEM), and metabolomics were employed. In addition, the content of adenosine triphosphate (ATP) in bacterium and the activities of key enzymes involved in critical metabolic pathways were detected using a variety of biochemical assays. The results showed that EO exhibited significant antibacterial activity against C. acnes at a minimum inhibitory concentration (MIC) of 400 μg/mL and a minimum bactericidal concentration (MBC) of 800 μg/mL, and EO could destroy C. acnes morphology and inhibit its growth. Moreover, results from our study showed that EO had a significant effect on the C. acnes normal metabolism. In total, 86 metabolites were altered, and 34 metabolic pathways related to the carbohydrate metabolism, energy metabolism, amino acid metabolism, as well as cell wall and cell membrane synthesis were perturbed after EO administration. The synthesis of ATP in bacterial cells was also severely inhibited, and the activities of key enzymes of the glycolysis and Wood-Werkman cycle were significantly affected (Pyruvate Carboxylase, Malate Dehydrogenase and Pyruvate kinase activities were decreased, and Hexokinase was increased). Taken together, these results illustrated that the bacteriostatic effect of EO against C. acnes by breaking the bacterial cell morphology and perturbing cell metabolism, including inhibition of key enzyme activity and ATP synthesis. The results from our study may shed new light on the discovery of novel drugs with more robust efficacy.

Chromosome-level genome assembly of Aristolochia contorta provides insights into the biosynthesis of benzylisoquinoline alkaloids and aristolochic acids

Aristolochic acids (AAs) and their derivatives exist in multiple Aristolochiaceae species which had been or are being used as medicinal materials. During the past decades, AAs have received increasing attention due to their nephrotoxicity and carcinogenecity. Elimination of AAs in medicinal materials using biotechnological approaches is important to improve medication safety. However, it has not been achieved because of the limited information of AA biosynthesis available. Here, we report a high-quality reference-grade genome assembly of the AA-containing vine, Aristolochia contorta. Total size of the assembly is 209.27 Mb, which is assembled into 7 pseudochromosomes. Synteny analysis, Ks distribution and 4DTv suggest absences of whole-genome duplication events in A. contorta after the angiosperm-wide WGD. Based on genomic, transcriptomic and metabolic data, pathways and candidate genes of benzylisoquinoline alkaloid (BIA) and AA biosynthesis in A. contorta were proposed. Five O-methyltransferase genes, including AcOMT1-3, AcOMT5 and AcOMT7, were cloned and functionally characterized. The results provide a high-quality reference genome for AA-containing species of Aristolochiaceae. It lays a solid foundation for further elucidation of AA biosynthesis and regulation and molecular breeding of Aristolochiaceae medicinal materials.

The Diesel Tree Sindora glabra Genome Provides Insights Into the Evolution of Oleoresin Biosynthesis

Sindora glabra is an economically important tree that produces abundant oleoresin in the trunk. Here, we present a high-quality chromosome-scale assembly of S. glabra genome by combining Illumina HiSeq, Pacific Biosciences sequencing, and Hi-C technologies. The size of S. glabra genome was 1.11 Gb, with a contig N50 of 1.27 Mb and 31,944 predicted genes. This is the first sequenced genome of the subfamily Caesalpinioideae. As a sister taxon to Papilionoideae, S. glabra underwent an ancient genome triplication shared by core eudicots and further whole-genome duplication shared by early-legume in the last 73.3 million years. S. glabra harbors specific genes and expanded genes largely involved in stress responses and biosynthesis of secondary metabolites. Moreover, 59 terpene backbone biosynthesis genes and 64 terpene synthase genes were identified, which together with co-expressed transcription factors could contribute to the diversity and specificity of terpene compounds and high terpene content in S. glabra stem. In addition, 63 disease resistance NBS-LRR genes were found to be unique in S. glabra genome and their expression levels were correlated with the accumulation of terpene profiles, suggesting potential defense function of terpenes in S. glabra. These together provide new resources for understanding genome evolution and oleoresin production.

Comprehensive identification of bHLH transcription factors in Litsea cubeba reveals candidate gene involved in the monoterpene biosynthesis pathway

Litsea cubeba (Lour.) Person, an economically important aromatic plant producing essential oils, has lemon-like fragrance and 96.44-98.44% monoterpene contents. bHLH transcription factor plays an important role in plant secondary metabolism and terpene biosynthesis. In this study, we used bioinformatics to identify bHLH transcription factors in L. cubeba, 173 bHLH genes were identified from L. cubeba and divided these into 26 subfamilies based on phylogenetic analysis. The majority of bHLHs in each subfamily shared comparable structures and motifs. While LcbHLHs were unevenly distributed across 12 chromosomes, 10 tandem repeats were discovered. Expression profiles of bHLH genes in different tissues demonstrated that LcbHLH78 is a potential candidate gene for regulating monoterpene biosynthesis. LcbHLH78 and the terpene synthase LcTPS42 showed comparable expression patterns in various tissues and fruit development stages of L. cubeba. Subcellular localization analysis revealed that LcbHLH78 protein localizes to the nucleus, consistent with a transcription factor function. Importantly, transient overexpression of LcbHLH78 increased geraniol and linalol contents. Our research demonstrates that LcbHLH78 enhances terpenoid biosynthesis. This finding will be beneficial for improving the quality of L. cubeba and provides helpful insights for further research into the control mechanism of LcbHLH genes over terpenoid biosynthesis.

Chromosome-level genome assembly of Amomum tsao-ko provides insights into the biosynthesis of flavor compounds

Amomum tsao-ko is an economically important spice plant in the ginger family (Zingiberaceae). The dried ripe fruit has been widely used as spice and medicine in Southeast Asia due to its distinct flavor metabolites. However, there is little genomic information available to understand the biosynthesis of its characteristic flavor compounds. Here, we present a high-quality chromosome-level genome of A. tsao-ko with a total length of 2.08 Gb assembled into 24 chromosomes. Potential relationships between genetic variation and chemical constituents were analyzed by a genome-wide association study of 119 representative A. tsao-ko specimens in China. Metabolome and transcriptome correlation analysis of different plant organs and fruit developmental stages revealed the proposed biosynthesis of the characteristic bicyclononane aldehydes and aromatic metabolites in A. tsao-ko fruit. Transcription factors of 20 families may be involved in the regulatory network of terpenoids. This study provides genomic and chemical insights into the biosynthesis of characteristic aroma and flavor constituents, which can be used to improve the quality of A. tsao-ko as food and medicine.

NmRF: identification of multispecies RNA 2'-O-methylation modification sites from RNA sequences

2'-O-methylation (Nm) is a post-transcriptional modification of RNA that is catalyzed by 2'-O-methyltransferase and involves replacing the H on the 2'-hydroxyl group with a methyl group. The 2'-O-methylation modification site is detected in a variety of RNA types (miRNA, tRNA, mRNA, etc.), plays an important role in biological processes and is associated with different diseases. There are few functional mechanisms developed at present, and traditional high-throughput experiments are time-consuming and expensive to explore functional mechanisms. For a deeper understanding of relevant biological mechanisms, it is necessary to develop efficient and accurate recognition tools based on machine learning. Based on this, we constructed a predictor called NmRF based on optimal mixed features and random forest classifier to identify 2'-O-methylation modification sites. The predictor can identify modification sites of multiple species at the same time. To obtain a better prediction model, a two-step strategy is adopted; that is, the optimal hybrid feature set is obtained by combining the light gradient boosting algorithm and incremental feature selection strategy. In 10-fold cross-validation, the accuracies of Homo sapiens and Saccharomyces cerevisiae were 89.069 and 93.885%, and the AUC were 0.9498 and 0.9832, respectively. The rigorous 10-fold cross-validation and independent tests confirm that the proposed method is significantly better than existing tools. A user-friendly web server is accessible at http://lab.malab.cn/∼acy/NmRF.

Chloranthus genome provides insights into the early diversification of angiosperms

Chloranthales remain the last major mesangiosperm lineage without a nuclear genome assembly. We therefore assemble a high-quality chromosome-level genome of Chloranthus spicatus to resolve enigmatic evolutionary relationships, as well as explore patterns of genome evolution among the major lineages of mesangiosperms (eudicots, monocots, magnoliids, Chloranthales, and Ceratophyllales). We find that synteny is highly conserved between genomic regions of Amborella, Vitis, and Chloranthus. We identify an ancient single whole-genome duplication (WGD) (κ) prior to the divergence of extant Chloranthales. Phylogenetic inference shows Chloranthales as sister to magnoliids. Furthermore, our analyses indicate that ancient hybridization may account for the incongruent phylogenetic placement of Chloranthales + magnoliids relative to monocots and eudicots in nuclear and chloroplast trees. Long genes and long introns are found to be prevalent in both Chloranthales and magnoliids compared to other angiosperms. Overall, our findings provide an improved context for understanding mesangiosperm relationships and evolution and contribute a valuable genomic resource for future investigations.

The Chloranthus sessilifolius genome provides insight into early diversification of angiosperms

Most extant angiosperms belong to Mesangiospermae, which comprises eudicots, monocots, magnoliids, Chloranthales and Ceratophyllales. However, phylogenetic relationships between these five lineages remain unclear. Here, we report the high-quality genome of a member of the Chloranthales lineage (Chloranthus sessilifolius). We detect only one whole genome duplication within this species and find that polyploidization events in different Mesangiospermae lineage are mutually independent. We also find that the members of all floral development-related gene lineages are present in C. sessilifolius despite its extremely simplified flower. The AP1 and PI genes, however, show a weak floral tissue-specialized expression. Our phylogenomic analyses suggest that Chloranthales and magnoliids are sister groups, and both are together sister to the clade comprising Ceratophyllales and eudicots, while the monocot lineage is sister to all other Mesangiospermae. Our findings suggest that in addition to hybridization, incomplete lineage sorting may largely account for phylogenetic inconsistencies between the observed gene trees.