Long read sequencing of Toona sinensis (A. Juss) Roem: A chromosome-level reference genome for the family Meliaceae

Population structure and genetic diversity of Toona sinensis revealed by whole-genome resequencing

OBJECTIVES

Toona sinensis, commonly known as Chinese toon, is a perennial woody plant with significant economic and ecological importance. This study employed whole-genome resequencing of 180 T. sinensis samples collected from Shandong to analyze genetic variation and diversity, ultimately identifying 18,231 high-quality SNPs after rigorous quality control and linkage disequilibrium pruning. This comprehensive genomic resource provides novel insights into the genetic architecture of T. sinensis, facilitating the elucidation of population structure and supporting future breeding programs.

DATA DESCRIPTION

We performed whole-genome resequencing on 180 Toona sinensis samples, generating 1170.26 Gbp of clean data with a Q30 percentage of 93.69%. The average alignment rate to the reference genome was 96.72%, with an average coverage depth of 8 × and a genome coverage of 88.71%. Following data quality control and alignment, we performed SNP calling and filtering to identify high-quality SNPs across all samples. Population structure analyses were then conducted using the identified SNPs, including principal component analysis (PCA), structure analysis, and phylogenetic tree construction. These comprehensive analyses provide a foundation for understanding the genetic diversity and evolutionary dynamics of T. sinensis.

Integrated morphological, physiological and transcriptomic analyses reveal the responses of Toona sinensis seedlings to low-nitrogen stress

Nitrogen is one of the most essential elements for plant growth and development. In this study, the growth, physiology, and transcriptome of Toona sinensis (A. Juss) Roem seedlings were compared between low-nitrogen (LN) and normal-nitrogen (NN) conditions. These results indicate that LN stress adversely influences T. sinensis seedling growth. The activities of key enzymes related to nitrogen assimilation and phytohormone contents were altered by LN stress. A total of 2828 differentially expressed genes (DEGs) in roots and 1547 in leaves were identified between the LN and NN treatments. A differential enrichment analysis of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways indicated that nitrogen and sugar metabolism, flavonoid biosynthesis, plant hormone signal transduction, and ABC transporters, were strongly affected by LN stress. In summary, this research provides information for further understanding the response of T. sinensis to LN stress.

Integrating GC-MS and comparative transcriptome analysis reveals that TsERF66 promotes the biosynthesis of caryophyllene in Toona sinensis tender leaves

Introduction

The strong aromatic characteristics of the tender leaves of Toona sinensis determine their quality and economic value.

Methods and results

Here, GC-MS analysis revealed that caryophyllene is a key volatile compound in the tender leaves of two different T. sinensis varieties, however, the transcriptional mechanisms controlling its gene expression are unknown. Comparative transcriptome analysis revealed significant enrichment of terpenoid synthesis pathway genes, suggesting that the regulation of terpenoid synthesis-related gene expression is an important factor leading to differences in aroma between the two varieties. Further analysis of expression levels and genetic evolution revealed that TsTPS18 is a caryophyllene synthase, which was confirmed by transient overexpression in T. sinensis and Nicotiana benthamiana leaves. Furthermore, we screened an AP2/ERF transcriptional factor ERF-IX member, TsERF66, for the potential regulation of caryophyllene synthesis. The TsERF66 had a similar expression trend to that of TsTPS18 and was highly expressed in high-aroma varieties and tender leaves. Exogenous spraying of MeJA also induced the expression of TsERF66 and TsTPS18 and promoted the biosynthesis of caryophyllene. Transient overexpression of TsERF66 in T. sinensis significantly promoted TsTPS18 expression and caryophyllene biosynthesis.

Discussion

Our results showed that TsERF66 promoted the expression of TsTPS18 and the biosynthesis of caryophyllene in T. sinensis leaves, providing a strategy for improving the aroma of tender leaves.

Approaches to increase the validity of gene family identification using manual homology search tools

Identifying homologs is an important process in the analysis of genetic patterns underlying traits and evolutionary relationships among species. Analysis of gene families is often used to form and support hypotheses on genetic patterns such as gene presence, absence, or functional divergence which underlie traits examined in functional studies. These analyses often require precise identification of all members in a targeted gene family. Manual pipelines where homology search and orthology assignment tools are used separately are the most common approach for identifying small gene families where accurate identification of all members is important. The ability to curate sequences between steps in manual pipelines allows for simple and precise identification of all possible gene family members. However, the validity of such manual pipeline analyses is often decreased by inappropriate approaches to homology searches including too relaxed or stringent statistical thresholds, inappropriate query sequences, homology classification based on sequence similarity alone, and low-quality proteome or genome sequences. In this article, we propose several approaches to mitigate these issues and allow for precise identification of gene family members and support for hypotheses linking genetic patterns to functional traits.

Chromosome-scale genomes of commercially important mahoganies, Swietenia macrophylla and Khaya senegalensis

Mahogany species (family Meliaceae) are highly valued for their aesthetic and durable wood. Despite their economic and ecological importance, genomic resources for mahogany species are limited, hindering genetic improvement and conservation efforts. Here we perform chromosome-scale genome assemblies of two commercially important mahogany species: Swietenia macrophylla and Khaya senegalensis. By combining 10X sequencing and Hi-C data, we assemble high-quality genomes of 274.49 Mb (S. macrophylla) and 406.50 Mb (K. senegalensis), with scaffold N50 lengths of 8.51 Mb and 7.85 Mb, respectively. A total of 99.38% and 98.05% of the assembled sequences are anchored to 28 pseudo-chromosomes in S. macrophylla and K. senegalensis, respectively. We predict 34,129 and 31,908 protein-coding genes in S. macrophylla and K. senegalensis, respectively, of which 97.44% and 98.49% are functionally annotated. The chromosome-scale genome assemblies of these mahogany species could serve as a vital genetic resource, especially in understanding the properties of non-model woody plants. These high-quality genomes could support the development of molecular markers for breeding programs, conservation efforts, and the sustainable management of these valuable forest resources.

Integrated Transcriptomic and Metabolomic Analysis Reveal the Underlying Mechanism of Anthocyanin Biosynthesis in Toona sinensis Leaves

Toona sinensis, commonly known as Chinese Toon, is a plant species that possesses noteworthy value as a tree and vegetable. Its tender young buds exhibit a diverse range of colors, primarily determined by the presence and composition of anthocyanins and flavonoids. However, the underlying mechanisms of anthocyanin biosynthesis in Toona sinensis have been rarely reported. To explore the related genes and metabolites associated with composition of leaf color, we conducted an analysis of the transcriptome and metabolome of five distinct Toona clones. The results showed that differentially expressed genes and metabolites involved in anthocyanin biosynthesis pathway were mainly enriched. A conjoint analysis of transcripts and metabolites was carried out in JFC (red) and LFC (green), resulting in the identification of 510 genes and 23 anthocyanin-related metabolites with a positive correlation coefficient greater than 0.8. Among these genes and metabolites, 23 transcription factors and phytohormone-related genes showed strong coefficients with 13 anthocyanin derivates, which mainly belonged to the stable types of delphinidin, cyanidin, peonidin. The core derivative was found to be Cyanidin-3-O-arabinoside, which was present in JFC at 520.93 times the abundance compared to LFC. Additionally, the regulatory network and relative expression levels of genes revealed that the structural genes DFR, ANS, and UFGT1 might be directly or indirectly regulated by the transcription factors SOC1 (MADS-box), CPC (MYB), and bHLH162 (bHLH) to control the accumulation of anthocyanin. The expression of these genes was significantly higher in red clones compared to green clones. Furthermore, RNA-seq results accurately reflected the true expression levels of genes. Overall, this study provides a foundation for future research aimed at manipulating anthocyanin biosynthesis to improve plant coloration or to derive human health benefits.

Post-Cyclization Skeletal Rearrangements in Plant Triterpenoid Biosynthesis by a Pair of Branchpoint Isomerases

Triterpenoids possess potent biological activities, but their polycyclic skeletons are challenging to synthesize. The skeletal diversity of triterpenoids in plants is generated by oxidosqualene cyclases based on epoxide-triggered cationic rearrangement cascades. Normally, triterpenoid skeletons then remain unaltered during subsequent tailoring steps. In contrast, the highly modified triterpenoids found in Sapindales plants imply the existence of post-cyclization skeletal rearrangement enzymes that have not yet been found. We report here a biosynthetic pathway in Sapindales plants for the modification of already cyclized tirucallane triterpenoids, controlling the pathway bifurcation between different plant triterpenoid classes. Using a combination of bioinformatics, heterologous expression in plants and chemical analyses, we identified a cytochrome P450 monooxygenase and two isomerases which harness the epoxidation-rearrangement biosynthetic logic of triterpene cyclizations for modifying the tirucallane scaffold. The two isomerases share the same epoxide substrate made by the cytochrome P450 monooxygenase CYP88A154, but generate two different rearrangement products, one containing a cyclopropane ring. Our findings reveal a process for skeletal rearrangements of triterpenoids in nature that expands their scaffold diversity after the initial cyclization. In addition, the enzymes described here are crucial for the biotechnological production of limonoid, quassinoid, apoprotolimonoid, and glabretane triterpenoids.

Meliaceae genomes provide insights into wood development and limonoids biosynthesis

Meliaceae is a useful plant family owing to its high-quality timber and its many limonoids that have pharmacological and biological activities. Although some genomes of Meliaceae species have been reported, many questions regarding their unique family features, namely wood quality and natural products, have not been answered. In this study, we provide the whole-genome sequence of Melia azedarach comprising 237.16 Mb with a contig N50 of 8.07 Mb, and an improved genome sequence of Azadirachta indica comprising 223.66 Mb with a contig N50 of 8.91 Mb. Moreover, genome skimming data, transcriptomes and other published genomes were comprehensively analysed to determine the genes and proteins that produce superior wood and valuable limonoids. Phylogenetic analysis of chloroplast genomes, single-copy gene families and single-nucleotide polymorphisms revealed that Meliaceae should be classified into two subfamilies: Cedreloideae and Melioideae. Although the Meliaceae species did not undergo additional whole-genome duplication events, the secondary wall biosynthetic genes of the woody Cedreloideae species, Toona sinensis, expanded significantly compared to those of A. indica and M. azedarach, especially in downstream transcription factors and cellulose/hemicellulose biosynthesis-related genes. Moreover, expanded special oxidosqualene cyclase catalogues can help diversify Sapindales skeletons, and the clustered genes that regulate terpene chain elongation, cyclization and modification would support their roles in limonoid biosynthesis. The expanded clans of terpene synthase, O-methyltransferase and cytochrome P450, which are mainly derived from tandem duplication, are responsible for the different limonoid classes among the species. These results are beneficial for further investigations of wood development and limonoid biosynthesis.

Genome-wide analysis of WRKY transcription factor genes in Toona sinensis: An insight into evolutionary characteristics and terpene synthesis

WRKY transcription factors (TFs), one of the largest TF families, serve critical roles in the regulation of secondary metabolite production. However, little is known about the expression pattern of WRKY genes during the germination and maturation processes of Toona sinensis buds. In the present study, the new assembly of the T. sinensis genome was used for the identification of 78 TsWRKY genes, including gene structures, phylogenetic features, chromosomal locations, conserved protein domains, cis-regulatory elements, synteny, and expression profiles. Gene duplication analysis revealed that gene tandem and segmental duplication events drove the expansion of the TsWRKYs family, with the latter playing a key role in the creation of new TsWRKY genes. The synteny and evolutionary constraint analyses of the WRKY proteins among T. sinensis and several distinct species provided more detailed evidence of gene evolution for TsWRKYs. Besides, the expression patterns and co-expression network analysis show TsWRKYs may multi-genes co-participate in regulating terpenoid biosynthesis. The findings revealed that TsWRKYs potentially play a regulatory role in secondary metabolite synthesis, forming the basis for further functional characterization of WRKY genes with the intention of improving T. sinensis.

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.

Transcriptome profiling of Toona ciliata young stems in response to Hypsipyla robusta Moore

Toona ciliata is a traditional woody plant that can be used as a medicinal material in China. The extracts of its roots, stems, leaves, and flowers all have a wide range of bioactive compounds. However, T. ciliata has been facing an unresolved pest problem caused by Hypsipyla robusta Moore (HRM), which seriously affects its growth and development. In this study, the expression level of TcMYB3 gene reached the maximum (28-fold) at 12 h and transcriptome sequencing of young stems eaten by HRM for 0, 3, 12, and 21 h were performed. A large number of differentially expressed genes (DEGs) were identified including jointly up-regulated genes (263) and down-regulated genes (378). JA synthesis and signaling transduction, terpene biosynthesis, and MAPKs signaling pathway were analyzed in depth and found that TcOPR3, TcJAR1, TcJAZs, and TcTPS9 genes possessed anti-insect potential. Moreover, MYB and ERF transcription factor (TF) families were significantly strengthened to the point that they may participate in induced defense mechanisms in T. ciliata. These data not only provide insights into the molecular mechanisms in resistance of T. ciliata to HRM but also helps to explore the new biocontrol strategies against insects in eco-friendly woody plants.

A Chromosome-Level Genome Assembly of Toona ciliata (Meliaceae)

Toona ciliata Roem is an important timber species in the Toona genus of the Meliaceae family and an endangered species due to over-cutting and a low rate of natural regeneration in China. Although molecular markers have been applied to studying population genetic diversity, the absence of a reliable reference genome limits in-depth genetic conservation and evolutionary studies of this species. Here, we reported a high-quality assembly of the whole genome sequence of T. ciliata. The total assembled genome has 520.64 Mb in length anchored on 28 chromosomes (contig N50 = 4.48 Mb). A total of 42,159 genes were predicted after the ab initio, homology-based, and transcriptome analyses. A total of 41,284 protein-encoding genes (97.92%) were functionally annotated and 1,246 non-coding RNAs were identified in the T. ciliata genome. Phylogenomic analysis showed that T. ciliata was divergent at 15.06 (6-25) Ma from T. sinensis of the same genus Toona. This whole genome sequence provides a valuable resource to study the genetic conservation and molecular evolution of T. ciliata in the future.

Research progress of meliaceous limonoids from 2011 to 2021

Covering: July 2010 to December 2021Limonoids, a kind of natural tetranortriterpenoids with diverse skeletons and valuable insecticidal and medicinal bioactivities, are the characteristic metabolites of most plants of the Meliaceae family. The chemistry and bioactivities of meliaceous limonoids are a continuing hot area of natural products research; to date, about 2700 meliaceous limonoids have been identified. In particular, more than 1600, including thirty kinds of novel rearranged skeletons, have been isolated and identified in the past decade due to their wide distribution and abundant content in Meliaceae plants and active biosynthetic pathways. In addition to the discovery of new structures, many positive medicinal bioactivities of meliaceous limonoids have been investigated, and extensive achievements regarding the chemical and biological synthesis have been made. This review summarizes the recent research progress in the discovery of new structures, medicinal and agricultural bioactivities, and chem/biosynthesis of limonoids from the plants of the Meliaceae family during the past decade, with an emphasis on the discovery of limonoids with novel skeletons, the medicinal bioactivities and mechanisms, and chemical synthesis. The structures, origins, and bioactivities of other new limonoids were provided as ESI. Studies published from July 2010 to December 2021 are reviewed, and 482 references are cited.

Genomic Analysis Based on Chromosome-Level Genome Assembly Reveals an Expansion of Terpene Biosynthesis of Azadirachta indica

Azadirachta indica (neem), an evergreen tree of the Meliaceae family, is a source of the potent biopesticide azadirachtin. The lack of a chromosome-level assembly impedes an in-depth understanding of its genome architecture and the comparative genomic analysis of A. indica. Here, a high-quality genome assembly of A. indica was constructed using a combination of data from Illumina, PacBio, and Hi-C technology, which is the first chromosome-scale genome assembly of A. indica. Based on the length of our assembly, the genome size of A. indica is estimated to be 281 Mb anchored to 14 chromosomes (contig N50 = 6 Mb and scaffold N50 = 19 Mb). The genome assembly contained 115 Mb repetitive elements and 25,767 protein-coding genes. Evolutional analysis revealed that A. indica didn't experience any whole-genome duplication (WGD) event after the core eudicot γ event, but some genes and genome segment might likely experienced recent duplications. The secondary metabolite clusters, TPS genes, and CYP genes were also identified. Comparative genomic analysis revealed that most of the A. indica-specific TPS genes and CYP genes were located on the terpene-related clusters on chromosome 13. It is suggested that chromosome 13 may play an important role in the specific terpene biosynthesis of A. indica. The gene duplication events may be responsible for the terpene biosynthesis expansion in A. indica. The genomic dataset and genomic analysis created for A. indica will shed light on terpene biosynthesis in A. indica and facilitate comparative genomic research of the family Meliaceae.

Review on the Development and Applications of Medicinal Plant Genomes

With the development of sequencing technology, the research on medicinal plants is no longer limited to the aspects of chemistry, pharmacology, and pharmacodynamics, but reveals them from the genetic level. As the price of next-generation sequencing technology becomes affordable, and the long-read sequencing technology is established, the medicinal plant genomes with large sizes have been sequenced and assembled more easily. Although the review of plant genomes has been reported several times, there is no review giving a systematic and comprehensive introduction about the development and application of medicinal plant genomes that have been reported until now. Here, we provide a historical perspective on the current situation of genomes in medicinal plant biology, highlight the use of the rapidly developing sequencing technologies, and conduct a comprehensive summary on how the genomes apply to solve the practical problems in medicinal plants, like genomics-assisted herb breeding, evolution history revelation, herbal synthetic biology study, and geoherbal research, which are important for effective utilization, rational use and sustainable protection of medicinal plants.