Polyploidy underlies co-option and diversification of biosynthetic triterpene pathways in the apple tribe

Occurrence of White Flesh Color and Refreshing Flavor Following Phytoene Synthase 2A Gene Variation in Loquat Fruit

Loquat fruits from the apple tribe of Rosaceae ripen from late spring to early summer, when most fresh fruits are out of season. Compared with the orange-fleshed varieties, the white-fleshed varieties are usually preferentially chosen by consumers for a more favorable flavor. Though the breeding of white-fleshed cultivars with large fruit size greatly promoted the development of loquat production, the mechanisms of how fruits with lighter pigments are generated and how fruit flavor changes following flesh color shifting remain to be elucidated. Pigment measurements indicated carotenoids as the dominant pigment underlying the change in flesh color changing. Genotyping and haplotyping of 807 loquat accessions revealed that the rise of PSY2Ad-PSY2Ad genotype blocks carotenoids accumulation and confers to white flesh color of loquat fruit. Analysis of widely targeted metabolomes on 18 representative cultivars identified 1420 metabolites, with 223 differentially accumulated metabolites between the two groups. Further metabolite comparison demonstrated that low levels of bitter or astringent compounds, such as flavonoids, lignans and coumarins, phenolic acids, nucleotides, alkaloids, and terpenoids, may confer to a refreshing flavor of white-fleshed fruits. Furthermore, 18 metabolic biomarkers were identified by machine learning to distinguish fruits with diverse flesh colors. This work gives insights into the understanding of how fruit color variation associated with flavor and also promotes white-fleshed loquat breeding by genus-wide genotyping.

Divergent fatty acid desaturase 2 is essential for falcarindiol biosynthesis in carrot

Carrots produce diverse falcarin-type polyacetylenes (PAs), which function as important antipathogenic phytoalexins and have potential as anticancer agents. Despite their abundance in carrot root tissues, the biosynthesis and evolutionary origins of falcarindiol, a C17-PA, remain unknown. Fatty acid desaturase 2 (FAD2) enzymes diversify PAs by introducing various double and/or triple carbon-carbon bonds into fatty acid chains. Here, we use association analysis to identify candidate FAD2 genes involved in falcarindiol biosynthesis. Using a rapid tobacco transient expression system, we found that DcFAD2 enzymes are highly functionally redundant and promiscuous. Combinatorial assays also revealed unexpected synergistic and redirective effects among FAD2 enzymes, further complicating the biosynthetic pathway. CRISPR-Cas9-mediated mutagenesis and overexpression studies identified previously overlooked DcFAD2 hub genes as essential for falcarindiol production. Evolutionary analysis suggests that the expansion of DcFAD2 genes underpins the abundance of falcarindiol in carrot, independent of the biosynthetic gene cluster previously identified in tomato. This work underscores the complexity of the falcarin biosynthetic network and identifies hub genes essential for falcarindiol biosynthesis in carrot.

Metabolic Engineering of Corynebacterium glutamicum for Producing Different Types of Triterpenoids

Triterpenoids widely exist in nature with diverse structures and possess various functional properties and biological effects. However, research on triterpenoids biosynthesis in Corynebacterium glutamicum is still limited to squalene, which restricts the development of C. glutamicum to produce high-value triterpenoids. In this study, C. glutamicum was developed as an efficient and flexible platform for the biosynthesis of different types of triterpenoids. Squalene was synthesized and the titer was improved to 400.1 mg/L in flask combining strategies of metabolic engineering and fermentation optimization. Particularly, intracellular squalene accounted for more than 97%, addressing the problem of leaking squalene in C. glutamicum, which may restrict the subsequent synthesis of other triterpenoids derived from squalene. Furthermore, 201.9 mg/L (3S)-2,3-oxidosqualene (SQO) and 264.9 mg/L (3S,22S)-2,3,22,23-dioxidosqualene (SDO) were successfully synthesized in strains harboring heterogeneous squalene epoxidase from Arabidopsis thaliana with different expression strengths. Therefore, a platform for de novo triterpenoids synthesis based on SQO or SDO was constructed in C. glutamicum. For instance, biosynthesis of α-amyrin and α-onocerin was achieved for the first time by introducing oxidosqualene cyclases in SQO- and SDO-producing C. glutamicum strains, respectively. After optimization, the titer of α-amyrin and α-onocerin was improved to 65.3 and 136.85 mg/L, respectively. Furthermore, ursolic acid, derived from α-amyrin, was synthesized after expressing cytochrome P450 enzyme and its compatible cytochrome P450 reductases with a titer of 486 μg/L. For the first time, reactions of epoxidation, cyclization, and oxidation from squalene were achieved in C. glutamicum, leading to the production of different types of triterpenoids. Our study provides a new platform for the production of triterpenoids, which will be helpful for the large-scale production of triterpenoids employing C. glutamicum as a chassis strain.

Genome-wide identification of UDP-glycosyltransferases involved in flavonol glycosylation induced by UV-B irradiation in Eriobotrya japonica

Flavonol glycosides are secondary metabolites important for plant development and stress defense such as UV-B irradiation. UDP-glycosyltransferase (UGT) catalyzes the last step in the biosynthesis of flavonol glycosides. Eriobotrya japonica is abundant in flavonol glycosides, but UGTs responsible for accumulation of flavonol glycosides remain unknown. Here, 13 flavonol glycosides including monoglycosides and diglycosides were characterized in different tissues of loquat by LC-MS/MS. UV-B irradiation significantly increased the accumulation of four quercetin glycosides and two kaempferol glycosides in loquat fruit. Based on UGT gene family analysis, transcriptome analysis, enzyme assays of recombinant proteins as well as transient overexpression assays in Nicotiana benthamiana, three UGTs were identified, i.e. EjUGT78T4 as flavonol 3-O-galactosyltransferase, EjUGT78S3 as flavonol 3-O-glucosyltransferase, and EjUGT91AK7 as flavonol 1 → 6 rhamnosyltransferase. This work elucidates the formation of flavonol glycosides in loquat through UGT-mediated glycosylation.

Multi-omics analysis reveals the role of UGT72 family genes in arbutin biosynthesis in Pyrus and evolution driven by whole genome duplication

The UGT72 gene family encodes proteins that glycosylate phenylpropanoids, and thus contribute to the synthesis of various phenolic substances. However, their functional role and evolutionary history in Pyrus spp. remains poorly understood. Here we explored the evolution, amplification, coding region structural variation, and functional divergence of the UGT72 gene family and its subfamilies. Further, we identified functional genes involved in arbutin synthesis and functionally validated the key genes. 15 UGT72 genes were identified in the complete genome sequence and classified into two subfamilies of Pyrus betulifolia. Significant expansion of the UGT72 gene family occurred after genome duplication in P. betulifolia. 53.33 % of all UGT72 family genes were found to have undergone expansion via WGD/segmental duplication. A noteworthy discovery was that the amplification of functional genes such as PbUGT72B1714 during polyploidization, combined with the loss of vital motifs and variations at important sites within these genes, significantly impacted the diversification of arbutin metabolism. These findings offer novel insights into how gene gains and losses caused by WGDs have contributed to metabolic diversification and evolutionary adaptation in Pyrus, as well as a groundwork for more detailed investigations into the mechanisms of arbutin metabolism.

High-quality genome of Firmiana hainanensis provides insights into the evolution of Malvaceae subfamilies and the mechanism of their wood density formation

The Malvaceae family, the most diverse family in the order Malvales, consists of nine subfamilies. Within the Firmiana genus of the Sterculioideae subfamily, most species are considered globally vulnerable, yet their genomes remain unexplored. Here, we present a chromosome-level genome assembly for a representative Firmiana species, F. hainanensis, 2n = 40, totaling 1536 Mb. Phylogenomic analysis shows that F. hainanensis and Durio zibethinus have the closest evolutionary relationship, with an estimated divergence time of approximately 21 MYA and distinct polyploidization events in their histories. Evolutionary trajectory analyses indicate that fissions and fusions may play a crucial role in chromosome number variation (2n = 14 to 2n = 96). Analysis of repetitive elements among Malvaceae reveals that the Tekay subfamily (belonging to the Gypsy group) contributes to variation in genome size (ranging from 324 Mb to 1620 Mb). Additionally, genes associated with P450, peroxidase, and microtubules, and thereby related to cell wall biosynthesis, are significantly contracted in F. hainanensis, potentially leading to its lower wood density relative to Hopea hainanensis. Overall, our study provides insights into the evolution of chromosome number, genome size, and the genetic basis of cell wall biosynthesis in Malvaceae species.

Understanding metabolic diversification in plants: branchpoints in the evolution of specialized metabolism

Plants are chemical engineers par excellence. Collectively they make a vast array of structurally diverse specialized metabolites. The raw materials for building new pathways (genes encoding biosynthetic enzymes) are commonly recruited directly or indirectly from primary metabolism. Little is known about how new metabolic pathways and networks evolve in plants, or what key nodes contribute to branches that lead to the biosynthesis of diverse chemicals. Here we review the molecular mechanisms underlying the generation of biosynthetic branchpoints. We also consider examples in which new metabolites are formed through the joining of precursor molecules arising from different biosynthetic routes, a scenario that greatly increases both the diversity and complexity of specialized metabolism. Given the emerging importance of metabolic gene clustering in helping to identify new enzymes and pathways, we further cover the significance of biosynthetic gene clusters in relation to metabolic networks and dedicated biosynthetic pathways. In conclusion, an improved understanding of the branchpoints between metabolic pathways will be key in order to be able to predict and illustrate the complex structure of metabolic networks and to better understand the plasticity of plant metabolism. This article is part of the theme issue 'The evolution of plant metabolism'.

Insights into ACO genes across Rosaceae: evolution, expression, and regulatory networks in fruit development

BACKGROUND

ACO (1-aminocyclopropane-1-carboxylic acid) serves as a pivotal enzyme within the plant ethylene synthesis pathway, exerting influence over critical facets of plant biology such as flowering, fruit ripening, and seed development.

OBJECTIVE

This study aims to identify ACO genes from representative Rosaceae genomes, reconstruct their phylogenetic relationships by integrating synteny information, and investigate their expression patterns and networks during fruit development.

METHODS

we utilize a specialized Hidden Markov Model (HMM), crafted on the sequence attributes of ACO gene-encoded proteins, to systematically identify and analyze ACO gene family members across 12 representative species within the Rosaceae botanical family. Through transcriptome analysis, we delineate the expression patterns of ACO genes in six distinct Rosaceae fruits.

RESULTS

Our investigation reveals the presence of 62 ACO genes distributed among the surveyed Rosaceae species, characterized by hydrophilic proteins predominantly expressed within the cytoplasm. Phylogenetic analysis categorizes these ACO genes into three discernible classes, namely Class I, Class II, and Class III. Further scrutiny via collinearity assessment indicates a lack of collinearity relationships among these classes, highlighting variations in conserved motifs and promoter types within each class. Transcriptome analysis unveils significant disparities in both expression levels and trends of ACO genes in fruits exhibiting respiratory bursts compared to those that do not. Employing Weighted Gene Co-Expression Network Analysis (WGCNA), we discern that the co-expression correlation of ACO genes within loquat fruit notably differs from that observed in apples. Our findings, derived from Gene Ontology (GO) enrichment results, signify the involvement of ACO genes and their co-expressed counterparts in biological processes linked to terpenoid metabolism and carbohydrate synthesis in loquat. Moreover, our exploration of gene regulatory networks (GRN) highlights the potential pivotal role of the GNAT transcription factor (Ejapchr1G00010380) in governing the overexpression of the ACO gene (Ejapchr10G00001110) within loquat fruits.

CONCLUSION

The constructed HMM of ACO proteins offers a precise and systematic method for identifying plant ACO proteins, facilitating phylogenetic reconstruction. ACO genes from representative Rosaceae fruits exhibit diverse expression and regulative patterns, warranting further function characterizations.

Genome-wide identification of R2R3-MYB family in Eriobotrya japonica and functional analysis of EjMYB5 involved in proanthocyanidin biosynthesis

Loquat (Eriobotrya japonica Lindl.) is a popular fruit and medicinal plant. Proanthocyanidins (PAs), as one of the main types of flavonoids, are the key components of loquat fruit quality and medicinal properties. However, the identification of transcription factors (TFs) involved in PA accumulation in loquat remains limited. R2R3-MYB TFs play key regulatory role in PA accumulation in plants. In this study, 190 R2R3-MYB TFs were identified in loquat genome. Combined with transcriptome data, R2R3-MYB TF EjMYB5 involved in PA accumulation in loquat was isolated. EjMYB5 was transcriptional activator localized to nucleus. Expression of EjMYB5 was closely related to PA accumulation in loquat fruits. Heterogenous overexpression of EjMYB5 in tomato (Solanum lycopersicum) inhibited anthocyanin accumulation and promoted PA accumulation. Additionally, transient overexpression of EjMYB5 in tobacco (Nicotiana benthamiana) leaves promoted PA accumulation by upregulating flavonoid biosynthesis genes (NtDFR, NtANS, and NtLAR). Transcriptome analysis of EjMYB5-overexpressing tomato fruits suggested that EjMYB5 was involved in several biological pathways, including lipid metabolism, MAPK signaling, phenylpropanoid biosynthesis, and flavonoid biosynthesis. Collectively, our findings provided basic data for further analysis the function of R2R3-MYB TFs in loquat, and revealed that EjMYB5 functioned as PA accumulation in loquat.

Unraveling the role of cytochrome P450 enzymes in oleanane triterpenoid biosynthesis in arjuna tree

Triterpenoids (C30-isoprenoids) represent a major group of natural products with various physiological functions in plants. Triterpenoids and their derivatives have medicinal uses owing to diverse bioactivities. Arjuna (Terminalia arjuna) tree bark accumulates highly oxygenated β-amyrin-derived oleanane triterpenoids (e.g., arjunic acid, arjungenin, and arjunolic acid) with cardioprotective roles. However, biosynthetic routes and enzymes remain poorly understood. We mined the arjuna transcriptome and conducted cytochrome P450 monooxygenase (P450) assays using Saccharomyces cerevisiae and Nicotiana benthamiana to identify six P450s and two P450 reductases for oxidative modifications of oleanane triterpenoids. P450 assays using oleananes revealed a greater substrate promiscuity of C-2α and C-23 hydroxylases/oxidases than C-28 oxidases. CYP716A233 and CYP716A432 catalyzed β-amyrin/erythrodiol C-28 oxidation to produce oleanolic acid. C-2α hydroxylases (CYP716C88 and CYP716C89) converted oleanolic acid and hederagenin to maslinic acid and arjunolic acid. CYP716C89 also hydroxylated erythrodiol and oleanolic aldehyde. However, CYP714E107a and CYP714E107b catalyzed oleanolic acid/maslinic acid/arjunic acid, C-23 hydroxylation to form hederagenin, arjunolic acid and arjungenin, and hederagenin C-23 oxidation to produce gypsogenic acid, but at a lower rate than oleanolic acid C-23 hydroxylation. Overall, P450 substrate selectivity suggested that C-28 oxidation is the first P450-catalyzed oxidative modification in the arjuna triterpenoid pathway. However, the pathway might branch thereafter through C-2α/C-23 hydroxylation of oleanolic acid. Taken together, these results provided new insights into substrate range of P450s and unraveled biosynthetic routes of triterpenoids in arjuna. Moreover, complete elucidation and reconstruction of arjunolic acid pathway in S. cerevisiae and N. benthamiana suggested the utility of arjuna P450s in heterologous production of cardioprotective compounds.

Gene Duplication and Functional Diversification of MADS-Box Genes in Malus × domestica following WGD: Implications for Fruit Type and Floral Organ Evolution

The evolution of the MADS-box gene family is essential for the rapid differentiation of floral organs and fruit types in angiosperms. Two key processes drive the evolution of gene families: gene duplication and functional differentiation. Duplicated copies provide the material for variation, while advantageous mutations can confer new functions on gene copies. In this study, we selected the Rosaceae family, which includes a variety of fruit types and flower organs, as well as species that existed before and after whole-genome duplication (WGD). The results indicate that different fruit types are associated with different copies of MADS-box gene family duplications and WGD events. While most gene copies derived from WGD have been lost, MADS-box genes not only retain copies derived from WGD but also undergo further gene duplication. The sequences, protein structures, and expression patterns of these gene copies have undergone significant differentiation. This work provides a clear example of MADS-box genes in the context of gene duplication and functional differentiation, offering new insights into the evolution of fruit types and floral organs.

Chloroplast genomes of Eriobotrya elliptica and an unknown wild loquat "YN-1"

The chloroplast genomes of wild loquat can help to determine their place in the history of evolution. Here, we sequenced and assembled two novel wild loquat's chloroplast genomes, one is Eriobotrya elliptica, and the other is an unidentified wild loquat, which we named "YN-1". Their sizes are 159,471 bp and 159,399 bp, respectively. We also assembled a cultivated loquat named 'JFZ', its chloroplast genome size is 159,156 bp. A comparative study was conducted with six distinct species of loquats, including five wild loquats and one cultivated loquat. The results showed that both E. elliptica and "YN-1" have 127 genes, one gene more than E. fragrans, which is psbK. Regions trnF-GAA-ndhJ, petG-trnP-UGG, and rpl32-trnL-UAG were found to exhibit high variability. It was discovered that there was a positive selection on rpl22 and rps12. RNA editing analysis found several chilling stress-specific RNA editing sites, especially in rpl2 gene. Phylogenetic analysis results showed that "YN-1" is closely related to E. elliptica, E. obovata and E. henryi.

Three Alternaria Species, Including a New Species, Causing Leaf Spot Disease of Loquat (Eriobotrya japonica) in China

Loquat (Eriobotrya japonica) is an economically important subtropical fruit crop in China. Field surveys conducted in different loquat orchards located in Chongqing, Sichuan, and Fujian provinces between 2017 and 2020 resulted in a collection of 56 Alternaria-like isolates from trees exhibiting symptoms of loquat leaf spot. Multigene phylogenetic analyses using seven gene regions, namely, ITS, gapdh, RPB2, tef1, Alt a 1, endoPG, and OPA10-2, showed that all the isolates belonged to the genus Alternaria, and supporting morphological analysis identified them as members of species A. alternata, A. gaisen, and A. chongqingensis sp. nov. In vitro and in vivo pathogenicity tests showed all the identified species to be pathogenic and able to cause leaf spot disease on loquat. Moreover, comprehensive phylogenetic analyses employing all combinations of the above seven gene sequences revealed the capability of Alt a 1-tef1-endoPG to provide a well-resolved gene tree for Alternaria spp. at the species level. This study adds to the current knowledge on an unknown species (A. chongqingensis sp. nov.) and is the first report of A. gaisen in loquat worldwide.

Integrated Metabolomics and Transcriptomics Analysis Reveals New Insights into Triterpene Biosynthesis in Rosa rugosa

Rosa rugosa is highly regarded for its aesthetic and therapeutic qualities. In particular, R. rugosa's flowers are known to produce essential oils containing a mixture of volatile terpenes, phenylpropanoids, and other compounds. Despite this, extensive research exists on volatile terpenes in flowers, while the knowledge of non-volatile terpenes in distinct tissues is still limited. Using UPLC-ESI-MS/MS, a comprehensive analysis of the terpene metabolites in five different tissues of R. rugosa was conducted. These metabolites accumulated in distinct tissues, and the majority of them were triterpenoids. Transcriptome data were collected from five tissues using RNA-seq. Transcriptomics and metabolomics were utilized to evaluate the triterpene biosynthesis pathway, resulting in new insights into its regulation and biosynthesis. The RrOSC10 was identified as a key enzyme in converting 2,3-oxidosqualene into α-amyrin, potentially contributing to the triterpene biosynthesis pathway. Furthermore, the expression of the RrOSC10 gene was upregulated by salinity for 0.5 h and 1 h, with subsequent downregulation at 2 h. This study lays a foundation for future research on the biosynthesis and accumulation of triterpenes in R. rugosa.

Signatures of Adaptation and Purifying Selection in Highland Populations of Dasiphora fruticosa

High mountains harbor a considerable proportion of biodiversity, but we know little about how diverse plants adapt to the harsh environment. Here we finished a high-quality genome assembly for Dasiphora fruticosa, an ecologically important plant distributed in the Qinghai-Tibetan Plateau and lowland of the Northern Hemisphere, and resequenced 592 natural individuals to address how this horticulture plant adapts to highland. Demographic analysis revealed D. fruticosa underwent a bottleneck after Naynayxungla Glaciation. Selective sweep analysis of two pairs of lowland and highland populations identified 63 shared genes related to cell wall organization or biogenesis, cellular component organization, and dwarfism, suggesting parallel adaptation to highland habitats. Most importantly, we found that stronger purging of estimated genetic load due to inbreeding in highland populations apparently contributed to their adaptation to the highest mountain. Our results revealed how plants could tolerate the extreme plateau, which could provide potential insights for species conservation and crop breeding.

Genomic basis of the distinct biosynthesis of β-glucogallin, a biochemical marker for hydrolyzable tannin production, in three oak species

Hydrolyzable tannins (HTs), predominant polyphenols in oaks, are widely used in grape wine aging, feed additives, and human healthcare. However, the limited availability of a high-quality reference genome of oaks greatly hampered the recognition of the mechanism of HT biosynthesis. Here, high-quality reference genomes of three Asian oak species (Quercus variabilis, Quercus aliena, and Quercus dentata) that have different HT contents were generated. Multi-omics studies were carried out to identify key genes regulating HT biosynthesis. In vitro enzyme activity assay was also conducted. Dual-luciferase and yeast one-hybrid assays were used to reveal the transcriptional regulation. Our results revealed that β-glucogallin was a biochemical marker for HT production in the cupules of the three Asian oaks. UGT84A13 was confirmed as the key enzyme for β-glucogallin biosynthesis. The differential expression of UGT84A13, rather than enzyme activity, was the main reason for different β-glucogallin and HT accumulation. Notably, sequence variations in UGT84A13 promoters led to different trans-activating activities of WRKY32/59, explaining the different expression patterns of UGT84A13 among the three species. Our findings provide three high-quality new reference genomes for oak trees and give new insights into different transcriptional regulation for understanding β-glucogallin and HT biosynthesis in closely related oak species.

Transposable elements in Rosaceae: insights into genome evolution, expression dynamics, and syntenic gene regulation

Transposable elements (TEs) exert significant influence on plant genomic structure and gene expression. Here, we explored TE-related aspects across 14 Rosaceae genomes, investigating genomic distribution, transposition activity, expression patterns, and nearby differentially expressed genes (DEGs). Analyses unveiled distinct long terminal repeat retrotransposon (LTR-RT) evolutionary patterns, reflecting varied genome size changes among nine species over the past million years. In the past 2.5 million years, Rubus idaeus showed a transposition rate twice as fast as Fragaria vesca, while Pyrus bretschneideri displayed significantly faster transposition compared with Crataegus pinnatifida. Genes adjacent to recent TE insertions were linked to adversity resistance, while those near previous insertions were functionally enriched in morphogenesis, enzyme activity, and metabolic processes. Expression analysis revealed diverse responses of LTR-RTs to internal or external conditions. Furthermore, we identified 3695 pairs of syntenic DEGs proximal to TEs in Malus domestica cv. 'Gala' and M. domestica (GDDH13), suggesting TE insertions may contribute to varietal trait differences in these apple varieties. Our study across representative Rosaceae species underscores the pivotal role of TEs in plant genome evolution within this diverse family. It elucidates how these elements regulate syntenic DEGs on a genome-wide scale, offering insights into Rosaceae-specific genomic evolution.

Genome-wide identification and expression analysis of the Eriobotrya japonica TIFY gene family reveals its functional diversity under abiotic stress conditions

BACKGROUND

Plant-specific TIFY proteins are widely found in terrestrial plants and play important roles in plant adversity responses. Although the genome of loquat at the chromosome level has been published, studies on the TIFY family in loquat are lacking. Therefore, the EjTIFY gene family was bioinformatically analyzed by constructing a phylogenetic tree, chromosomal localization, gene structure, and adversity expression profiling in this study.

RESULTS

Twenty-six EjTIFY genes were identified and categorized into four subfamilies (ZML, JAZ, PPD, and TIFY) based on their structural domains. Twenty-four EjTIFY genes were irregularly distributed on 11 of the 17 chromosomes, and the remaining two genes were distributed in fragments. We identified 15 covariate TIFY gene pairs in the loquat genome, 13 of which were involved in large-scale interchromosomal segmental duplication events, and two of which were involved in tandem duplication events. Many abiotic stress cis-elements were widely present in the promoter region. Analysis of the Ka/Ks ratio showed that the paralogous homologs of the EjTIFY family were mainly subjected to purifying selection. Analysis of the RNA-seq data revealed that a total of five differentially expressed genes (DEGs) were expressed in the shoots under gibberellin treatment, whereas only one gene was significantly differentially expressed in the leaves; under both low-temperature and high-temperature stresses, there were significantly differentially expressed genes, and the EjJAZ15 gene was significantly upregulated under both low- and high-temperature stress. RNA-seq and qRT-PCR expression analysis under salt stress conditions revealed that EjJAZ2, EjJAZ4, and EjJAZ9 responded to salt stress in loquat plants, which promoted resistance to salt stress through the JA pathway. The response model of the TIFY genes in the jasmonic acid pathway under salt stress in loquat was systematically summarized.

CONCLUSIONS

These results provide a theoretical basis for exploring the characteristics and functions of additional EjTIFY genes in the future. This study also provides a theoretical basis for further research on breeding for salt stress resistance in loquat. RT-qPCR analysis revealed that the expression of one of the three EjTIFY genes increased and the expression of two decreased under salt stress conditions, suggesting that EjTIFY exhibited different expression patterns under salt stress conditions.

High-quality genomes of Bombax ceiba and Ceiba pentandra provide insights into the evolution of Malvaceae species and differences in their natural fiber development

Members of the Malvaceae family, including Corchorus spp., Gossypium spp., Bombax spp., and Ceiba spp., are important sources of natural fibers. In the past decade, the genomes of several Malvaceae species have been assembled; however, the evolutionary history of Malvaceae species and the differences in their fiber development remain to be clarified. Here, we report the genome assembly and annotation of two natural fiber plants from the Malvaceae, Bombax ceiba and Ceiba pentandra, whose assembled genome sizes are 783.56 Mb and 1575.47 Mb, respectively. Comparative analysis revealed that whole-genome duplication and Gypsy long terminal repeat retroelements have been the major causes of differences in chromosome number (2n = 14 to 2n = 96) and genome size (234 Mb to 2676 Mb) among Malvaceae species. We also used comparative genomic analyses to reconstruct the ancestral Malvaceae karyotype with 11 proto-chromosomes, providing new insights into the evolutionary trajectories of Malvaceae species. MYB-MIXTA-like 3 is relatively conserved among the Malvaceae and functions in fiber cell-fate determination in the epidermis. It appears to perform this function in any tissue where it is expressed, i.e. in fibers on the endocarp of B. ceiba and in ovule fibers of cotton. We identified a structural variation in a cellulose synthase gene and a higher copy number of cellulose synthase-like genes as possible causes of the finer, less spinnable, weaker fibers of B. ceiba. Our study provides two high-quality genomes of natural fiber plants and offers insights into the evolution of Malvaceae species and differences in their natural fiber formation and development through multi-omics analysis.

Comparative transcriptome and metabolite profiling reveal diverse pattern of CYP-TS gene expression during corosolic acid biosynthesis in Lagerstroemia speciosa (L.) Pers

KEY MESSAGE

Extensive leaf transcriptome profiling and differential gene expression analysis of field grown and elicited shoot cultures of L. speciosa suggest that differential synthesis of CRA is mediated primarily by CYP and TS genes, showing functional diversity. Lagerstroemia speciosa L. is a tree species with medicinal and horticultural attributes. The pentacyclic triterpene, Corosolic acid (CRA) obtained from this species is widely used for the management of diabetes mellitus in traditional medicine. The high mercantile value of the compound and limited availability of innate resources entail exploration of alternative sources for CRA production. Metabolic pathway engineering for enhanced bioproduction of plant secondary metabolites is an attractive proposition for which, candidate genes in the pathway need to be identified and characterized. Therefore, in the present investigation, we focused on the identification of cytochrome P450 (CYP450) and oxidosqualene cyclases (OSC) genes and their differential expression during biosynthesis of CRA. The pattern of differential expression of these genes in the shoot cultures of L. speciosa, elicited with different epigenetic modifiers (azacytidine (AzaC), sodium butyrate (NaBu) and anacardic acid (AA)), was studied in comparison with field grown plant. Further, in vitro cultures with varying (low to high) concentrations of CRA were systematically assessed for the expression of CYP-TS and associated genes involved in CRA biosynthesis by transcriptome sequencing. The sequenced samples were de novo assembled into 180,290 transcripts of which, 92,983 transcripts were further annotated by UniProt. The results are collectively given in co-occurrence heat maps to identify the differentially expressed genes. The combined transcript and metabolite profiles along with RT-qPCR analysis resulted in the identification of CYP-TS genes with high sequence variation. Further, instances of concordant/discordant relation between CRA biosynthesis and CYP-TS gene expression were observed, indicating functional diversity in genes.

Multi-omics analysis provides new insights into the changes of important nutrients and fructose metabolism in loquat bud sport mutant

Bud sport is a common and stable somatic variation in perennial fruit trees, and often leads to significant modification of fruit traits and affects the breeding value. To investigate the impact of bud sport on the main metabolites in the fruit of white-fleshed loquat, we conducted a multi-omics analysis of loquat fruits at different developmental stages of a white-fleshed bud sport mutant of Dongting loquat (TBW) and its wild type (TBY). The findings from the detection of main fruit quality indices and metabolites suggested that bud sport resulted in a reduction in the accumulation of carotenoids, fructose, titratable acid and terpenoids at the mature stage of TBW, while leading to the accumulation of flavonoids, phenolic acids, amino acids and lipids. The comparably low content of titratable acid further enhances the balanced and pleasent taste profile of TBW. Expression patterns of differentially expressed genes involved in fructose metabolism exhibited a significant increase in the expression level of S6PDH (EVM0006243, EVM0044405) prior to fruit maturation. The comparison of protein sequences and promoter region of S6PDH between TBY and TBW revealed no structural variations that would impact gene function or expression, indicating that transcription factors may be responsible for the rapid up-regulation of S6PDH before maturation. Furthermore, correlation analysis helped to construct a comprehensive regulatory network of fructose metabolism in loquat, including 23 transcription factors, six structural genes, and nine saccharides. Based on the regulatory network and existing studies, it could be inferred that transcription factors such as ERF, NAC, MYB, GRAS, and bZIP may promote fructose accumulation in loquat flesh by positively regulating S6PDH. These findings improve our understanding of the nutritional value and breeding potential of white-fleshed loquat bud sport mutant, as well as serve as a foundation for exploring the genes and transcription factors that regulate fructose metabolism in loquat.

Comprehensive Evolutionary Analysis of the SMXL Gene Family in Rosaceae: Further Insights into Its Origin, Expansion, Diversification, and Role in Regulating Pear Branching

SMXL genes constitute a conserved gene family that is ubiquitous in angiosperms and involved in regulating various plant processes, including branching, leaf elongation, and anthocyanin biosynthesis, but little is known about their molecular functions in pear branching. Here, we performed genome-wide identification and investigation of the SMXL genes in 16 angiosperms and analyzed their phylogenetics, structural features, conserved motifs, and expression patterns. In total, 121 SMXLs genes were identified and were classified into four groups. The number of non-redundant SMXL genes in each species varied from 3 (Amborella trichopoda Baill.) to 18 (Glycine max Merr.) and revealed clear gene expansion events over evolutionary history. All the SMXL genes showed conserved structures, containing no more than two introns. Three-dimensional protein structure prediction revealed distinct structures between but similar structures within groups. A quantitative real-time PCR analysis revealed different expressions of 10 SMXL genes from pear branching induced by fruit-thinning treatment. Overall, our study provides a comprehensive investigation of SMXL genes in the Rosaceae family, especially pear. The results offer a reference for understanding the evolutionary history of SMXL genes and provide excellent candidates for studying fruit tree branching regulation, and in facilitating pear pruning and planting strategies.

Genome-wide identification, evolution, and expression analysis of the NAC gene family in chestnut (Castanea mollissima)

The NAC gene family is one of the most important transcription factor families specific to plants, responsible for regulating many biological processes, including development, stress response, and signal transduction. However, it has not yet been characterized in chestnut, an important nut tree species. Here, we identified 115 CmNAC genes in the chestnut genome, which were divided into 16 subgroups based on the phylogenetic analysis. Numerous cis-acting elements related to auxin, gibberellin, and abscisic acid were identified in the promoter region of CmNACs, suggesting that they play an important role in the growth and development of chestnut. The results of the collinear analysis indicated that dispersed duplication and whole-genome-duplication were the main drivers of CmNAC gene expansion. RNA-seq data of developmental stages of chestnut nut, bud, and ovule revealed the expression patterns of CmNAC genes. Additionally, qRT-PCR experiments were used to verify the expression levels of some CmNAC genes. The comprehensive analysis of the above results revealed that some CmNAC members may be related to chestnut bud and nut development, as well as ovule fertility. The systematic analysis of this study will help to increase understanding of the potential functions of the CmNAC genes in chestnut growth and development.

Genome-Wide Identification and Expression Analysis of the SBP-Box Gene Family in Loquat Fruit Development

The loquat (Eriobotrya japonica L.) is a special evergreen tree, and its fruit is of high medical and health value as well as having stable market demand around the world. In recent years, research on the accumulation of nutrients in loquat fruit, such as carotenoids, flavonoids, and terpenoids, has become a hotspot. The SBP-box gene family encodes transcription factors involved in plant growth and development. However, there has been no report on the SBP-box gene family in the loquat genome and their functions in carotenoid biosynthesis and fruit ripening. In this study, we identified 28 EjSBP genes in the loquat genome, which were unevenly distributed on 12 chromosomes. We also systematically investigated the phylogenetic relationship, collinearity, gene structure, conserved motifs, and cis-elements of EjSBP proteins. Most EjSBP genes showed high expression in the root, stem, leaf, and inflorescence, while only five EjSBP genes were highly expressed in the fruit. Gene expression analysis revealed eight differentially expressed EjSBP genes between yellow- and white-fleshed fruits, suggesting that the EjSBP genes play important roles in loquat fruit development at the breaker stage. Notably, EjSBP01 and EjSBP19 exhibited completely opposite expression patterns between white- and yellow-fleshed fruits during fruit development, and showed a close relationship with SlCnr involved in carotenoid biosynthesis and fruit ripening, indicating that these two genes may participate in the synthesis and accumulation of carotenoids in loquat fruit. In summary, this study provides comprehensive information about the SBP-box gene family in the loquat, and identified two EjSBP genes as candidates involved in carotenoid synthesis and accumulation during loquat fruit development.

Phylogenomics insights into gene evolution, rapid species diversification, and morphological innovation of the apple tribe (Maleae, Rosaceae)

Maleae is one of the most widespread tribes of Rosaceae and includes several important fruit crops and ornamental plants. We used nuclear genes from 62 transcriptomes/genomes, including 26 newly generated transcriptomes, to reconstruct a well-supported phylogeny and study the evolution of fruit and leaf morphology and the possible effect of whole genome duplication (WGD). Our phylogeny recovered 11 well-supported clades and supported the monophyly of most genera (except Malus, Sorbus, and Pourthiaea) with at least two sampled species. A WGD was located to the most recent common ancestor (MRCA) of Maleae and dated to c. 54 million years ago (Ma) near the Early Eocene Climatic Optimum, supporting Gillenieae (x = 9) being a parental lineage of Maleae (x = 17) and including duplicate regulatory genes related to the origin of the fleshy pome fruit. Whole genome duplication-derived paralogs that are retained in specific lineages but lost in others are predicted to function in development, metabolism, and other processes. An upshift of diversification and innovations of fruit and leaf morphologies occurred at the MRCA of the Malinae subtribe, coinciding with the Eocene-Oligocene transition (c. 34 Ma), following a lag from the time of the WGD event. Our results provide new insights into the Maleae phylogeny, its rapid diversification, and morphological and molecular evolution.

Comparative genomics of the medicinal plants Lonicera macranthoides and L. japonica provides insight into genus genome evolution and hederagenin-based saponin biosynthesis

Lonicera macranthoides (LM) and L. japonica (LJ) are medicinal plants widely used in treating viral diseases, such as COVID-19. Although the two species are morphologically similar, their secondary metabolite profiles are significantly different. Here, metabolomics analysis showed that LM contained ~86.01 mg/g hederagenin-based saponins, 2000-fold higher than LJ. To gain molecular insights into its secondary metabolite production, a chromosome-level genome of LM was constructed, comprising 9 pseudo-chromosomes with 40 097 protein-encoding genes. Genome evolution analysis showed that LM and LJ were diverged 1.30-2.27 million years ago (MYA). The two plant species experienced a common whole-genome duplication event that occurred ∼53.9-55.2 MYA before speciation. Genes involved in hederagenin-based saponin biosynthesis were arranged in clusters on the chromosomes of LM and they were more highly expressed in LM than in LJ. Among them, oleanolic acid synthase (OAS) and UDP-glycosyltransferase 73 (UGT73) families were much more highly expressed in LM than in LJ. Specifically, LmOAS1 was identified to effectively catalyse the C-28 oxidation of β-Amyrin to form oleanolic acid, the precursor of hederagenin-based saponin. LmUGT73P1 was identified to catalyse cauloside A to produce α-hederin. We further identified the key amino acid residues of LmOAS1 and LmUGT73P1 for their enzymatic activities. Additionally, comparing with collinear genes in LJ, LmOAS1 and LmUGT73P1 had an interesting phenomenon of 'neighbourhood replication' in LM genome. Collectively, the genomic resource and candidate genes reported here set the foundation to fully reveal the genome evolution of the Lonicera genus and hederagenin-based saponin biosynthetic pathway.

Recent advances in triterpenoid pathway elucidation and engineering

Triterpenoids are among the most assorted class of specialized metabolites found in all the taxa of living organisms. Triterpenoids are the leading active ingredients sourced from plant species and are utilized in pharmaceutical and cosmetic industries. The triterpenoid precursor 2,3-oxidosqualene, which is biosynthesized via the mevalonate (MVA) pathway is structurally diversified by the oxidosqualene cyclases (OSCs) and other scaffold-decorating enzymes such as cytochrome P450 monooxygenases (P450s), UDP-glycosyltransferases (UGTs) and acyltransferases (ATs). A majority of the bioactive triterpenoids are harvested from the native hosts using the traditional methods of extraction and occasionally semi-synthesized. These methods of supply are time-consuming and do not often align with sustainability goals. Recent advancements in metabolic engineering and synthetic biology have shown prospects for the green routes of triterpenoid pathway reconstruction in heterologous hosts such as Escherichia coli, Saccharomyces cerevisiae and Nicotiana benthamiana, which appear to be quite promising and might lead to the development of alternative source of triterpenoids. The present review describes the biotechnological strategies used to elucidate complex biosynthetic pathways and to understand their regulation and also discusses how the advances in triterpenoid pathway engineering might aid in the scale-up of triterpenoid production in engineered hosts.

Characterization of the horse chestnut genome reveals the evolution of aescin and aesculin biosynthesis

Horse chestnut (Aesculus chinensis) is an important medicinal tree that contains various bioactive compounds, such as aescin, barrigenol-type triterpenoid saponins (BAT), and aesculin, a glycosylated coumarin. Herein, we report a 470.02 Mb genome assembly and characterize an Aesculus-specific whole-genome duplication event, which leads to the formation and duplication of two triterpenoid biosynthesis-related gene clusters (BGCs). We also show that AcOCS6, AcCYP716A278, AcCYP716A275, and AcCSL1 genes within these two BGCs along with a seed-specific expressed AcBAHD6 are responsible for the formation of aescin. Furthermore, we identify seven Aesculus-originated coumarin glycoside biosynthetic genes and achieve the de novo synthesis of aesculin in E. coli. Collinearity analysis shows that the collinear BGC segments can be traced back to early-diverging angiosperms, and the essential gene-encoding enzymes necessary for BAT biosynthesis are recruited before the splitting of Aesculus, Acer, and Xanthoceras. These findings provide insight on the evolution of gene clusters associated with medicinal tree metabolites.

Insights into the Evolution of Ohnologous Sequences and Their Epigenetic Marks Post-WGD in Malus Domestica

A Whole Genome Duplication (WGD) event occurred several Ma in a Rosaceae ancestor, giving rise to the Maloideae subfamily which includes today many pome fruits such as pear (Pyrus communis) and apple (Malus domestica). This complete and well-conserved genome duplication makes the apple an organism of choice to study the early evolutionary events occurring to ohnologous chromosome fragments. In this study, we investigated gene sequence evolution and expression, transposable elements (TE) density, and DNA methylation level. Overall, we identified 16,779 ohnologous gene pairs in the apple genome, confirming the relatively recent WGD. We identified several imbalances in QTL localization among duplicated chromosomal fragments and characterized various biases in genome fractionation, gene transcription, TE densities, and DNA methylation. Our results suggest a particular chromosome dominance in this autopolyploid species, a phenomenon that displays similarities with subgenome dominance that has only been described so far in allopolyploids.

The Identification of SQS/SQE/OSC Gene Families in Regulating the Biosynthesis of Triterpenes in Potentilla anserina

The tuberous roots of Potentilla anserina (Pan) are an edible and medicinal resource in Qinghai-Tibetan Plateau, China. The triterpenoids from tuberous roots have shown promising anti-cancer, hepatoprotective, and anti-inflammatory properties. In this study, we carried out phylogenetic analysis of squalene synthases (SQSs), squalene epoxidases (SQEs), and oxidosqualene cyclases (OSCs) in the pathway of triterpenes. In total, 6, 26, and 20 genes of SQSs, SQEs, and OSCs were retrieved from the genome of Pan, respectively. Moreover, 6 SQSs and 25 SQEs genes expressed in two sub-genomes (A and B) of Pan. SQSs were not expanded after whole-genome duplication (WGD), and the duplicated genes were detected in SQEs. Twenty OSCs were divided into two clades of cycloartenol synthases (CASs) and β-amyrin synthases (β-ASs) by a phylogenetic tree, characterized with gene duplication and evolutionary divergence. We speculated that β-ASs and CASs may participate in triterpenes synthesis. The data presented act as valuable references for future studies on the triterpene synthetic pathway of Pan.

Comparative genomics reveals the diversification of triterpenoid biosynthesis and origin of ocotillol-type triterpenes in Panax

Gene duplication is assumed to be the major force driving the evolution of metabolite biosynthesis in plants. Freed from functional burdens, duplicated genes can mutate toward novelties until fixed due to selective fitness. However, the extent to which this mechanism has driven the diversification of metabolite biosynthesis remains to be tested. Here we performed comparative genomics analysis and functional characterization to evaluate the impact of gene duplication on the evolution of triterpenoid biosynthesis using Panax species as models. We found that whole-genome duplications (WGDs) occurred independently in Araliaceae and Apiaceae lineages. Comparative genomics revealed the evolutionary trajectories of triterpenoid biosynthesis in plants, which was mainly promoted by WGDs and tandem duplication. Lanosterol synthase (LAS) was likely derived from a tandem duplicate of cycloartenol synthase that predated the emergence of Nymphaeales. Under episodic diversifying selection, the LAS gene duplicates produced by γ whole-genome triplication have given rise to triterpene biosynthesis in core eudicots through neofunctionalization. Moreover, functional characterization revealed that oxidosqualene cyclases (OSCs) responsible for synthesizing dammarane-type triterpenes in Panax species were also capable of producing ocotillol-type triterpenes. Genomic and biochemical evidence suggested that Panax genes encoding the above OSCs originated from the specialization of one OSC gene duplicate produced from a recent WGD shared by Araliaceae (Pg-β). Our results reveal the crucial role of gene duplication in diversification of triterpenoid biosynthesis in plants and provide insight into the origin of ocotillol-type triterpenes in Panax species.

Organic Acid Accumulation and Associated Dynamic Changes in Enzyme Activity and Gene Expression during Fruit Development and Ripening of Common Loquat and Its Interspecific Hybrid

Loquats have gained increasing attention from consumers and growers for their essential nutrients and unusual phenology, which could help plug a gap period at market in early spring. Fruit acid is a critical contributor to fruit quality. The dynamic changes in organic acid (OA) during fruit development and ripening of common loquat (Dawuxing, DWX) and its interspecific hybrid (Chunhua, CH) were compared, as well as the corresponding enzyme activity and gene expression. At harvest, titratable acid was significantly lower (p ≤ 0.01) in CH (0.11%) than in DWX loquats (0.35%). As the predominant OA compound, malic acid accounted for 77.55% and 48.59% of the total acid of DWX and CH loquats at harvest, followed by succinic acid and tartaric acid, respectively. PEPC and NAD-MDH are key enzymes that participate in malic acid metabolism in loquat. The OA differences in DWX loquat and its interspecific hybrid could be attributed to the coordinated regulation of multiple genes and enzymes associated with OA biosynthesis, degradation, and transport. The data obtained in this work will serve as a fundamental and important basis for future loquat breeding programs and even for improvements in loquat cultural practices.

Genome assembly of wild loquat (Eriobotrya japonica) and resequencing provide new insights into the genomic evolution and fruit domestication in loquat

Wild loquats (Eriobotrya japonica Lindl.) provide remarkable genetic resources for studying domestication and breeding improved varieties. Herein, we generate the first high-quality chromosome-level genome assembly of wild loquat, with 45 791 predicted protein-coding genes. Analysis of comparative genomics indicated that loquat shares a common ancestor with apple and pear, and a recent whole-genome duplication event occurred in loquat prior to its divergence. Genome resequencing showed that the loquat germplasms can be distinctly classified into wild and cultivated groups, and the commercial cultivars have experienced allelic admixture. Compared with cultivated loquats, the wild loquat genome showed very few selected genomic regions and had higher levels of genetic diversity. However, whole-genome scans of selective sweeps were mainly related to fruit quality, size, and flesh color during the domestication process. Large-scale transcriptome and metabolome analyses were further performed to identify differentially expressed genes (DEGs) and differentially accumulated metabolites (DAMs) in wild and cultivated loquats at various fruit development stages. Unlike those in wild loquat, the key DEGs and DAMs involved in carbohydrate metabolism, plant hormone signal transduction, flavonoid biosynthesis, and carotenoid biosynthesis were significantly regulated in cultivated loquats during fruit development. These high-quality reference genome, resequencing, and large-scale transcriptome/metabolome data provide valuable resources for elucidating fruit domestication and molecular breeding in loquat.

Genome-Wide Analysis of SAUR Gene Family Identifies a Candidate Associated with Fruit Size in Loquat (Eriobotrya japonica Lindl.)

Fruit size is an important fruit quality trait that influences the production and commodity values of loquats (Eriobotrya japonica Lindl.). The Small Auxin Upregulated RNA (SAUR) gene family has proven to play a vital role in the fruit development of many plant species. However, it has not been comprehensively studied in a genome-wide manner in loquats, and its role in regulating fruit size remains unknown. In this study, we identified 95 EjSAUR genes in the loquat genome. Tandem duplication and segmental duplication contributed to the expansion of this gene family in loquats. Phylogenetic analysis grouped the SAURs from Arabidopsis, rice, and loquat into nine clusters. By analyzing the transcriptome profiles in different tissues and at different fruit developmental stages and comparing two sister lines with contrasting fruit sizes, as well as by functional predictions, a candidate gene (EjSAUR22) highly expressed in expanding fruits was selected for further functional investigation. A combination of Indoleacetic acid (IAA) treatment and virus-induced gene silencing revealed that EjSAUR22 was not only responsive to auxin, but also played a role in regulating cell size and fruit expansion. The findings from our study provide a solid foundation for understanding the molecular mechanisms controlling fruit size in loquats, and also provide potential targets for manipulation of fruit size to accelerate loquat breeding.

Plant metabolic gene clusters in the multi-omics era

Secondary metabolism in plants gives rise to a vast array of small-molecule natural products. The discovery of operon-like gene clusters in plants has provided a new perspective on the evolution of specialized metabolism and the opportunity to rapidly advance the metabolic engineering of natural product production. Here, we review historical aspects of the study of plant metabolic gene clusters as well as general strategies for identifying plant metabolic gene clusters in the multi-omics era. We also emphasize the exploration of their natural variation and evolution, as well as new strategies for the prospecting of plant metabolic gene clusters and a deeper understanding of how their structure influences their function.

Optimization of Ultrasonic Extraction of Triterpenes from Loquat Peel and Pulp and Determination of Antioxidant Activity and Triterpenoid Components

The aim of this paper was to study the optimal extraction process of total triterpenes from loquat peel and pulp assisted by ultrasound. The effects of solid–liquid ratio, ethanol concentration, ultrasonic time, ultrasonic power, and ultrasonic temperature on the yield of triterpenoid acid in loquat were investigated by single-factor and response surface methodology. FRAP (Ferric ion reducing antioxidant power) method, ABTS (2,2′-Azino-bis(3-ethylbenzthiazoline-6-sulfonic acid)) method, and DPPH (1,1-Diphenyl-2-picrylhydrazyl) method were used to determine the antioxidant capacity of peel and pulp at different stages. LC-MS (Liquid Chromatograph Mass Spectrometer) was used to qualitatively analyze different tissues of loquat. The optimal extraction conditions were as follows: ethanol concentration 71%, ultrasonic time 45 min, ultrasonic power 160 W, solid–liquid ratio 1:10, and ultrasonic temperature 30 °C. The total triterpenoid content of loquat peel was 13.92 ± 0.20 mg/g. The optimal extraction conditions were ethanol concentration 85%, ultrasonic time 51 min, ultrasonic power 160 W, solid–liquid ratio 1:8, and ultrasonic temperature 43 °C. The total triterpenoid content of loquat pulp was 11.69 ± 0.25 mg/g. The contents of triterpenes and antioxidant capacity in the peel and pulp of loquat at the three stages were the highest in the fruit ripening stage (S3). LC-MS analysis showed that most of the triterpenes belonged to ursolic acid derivatives and oleanolic acid derivatives, which laid the foundation for further utilization and development of loquat peel and pulp.

Plant biosynthetic gene clusters in the context of metabolic evolution

Covering: up to 2022Plants produce a wide range of structurally and biosynthetically diverse natural products to interact with their environment. These specialised metabolites typically evolve in limited taxonomic groups presumably in response to specific selective pressures. With the increasing availability of sequencing data, it has become apparent that in many cases the genes encoding biosynthetic enzymes for specialised metabolic pathways are not randomly distributed on the genome. Instead they are physically linked in structures such as arrays, pairs and clusters. The exact function of these clusters is debated. In this review we take a broad view of gene arrangement in plant specialised metabolism, examining types of structures and variation. We discuss the evolution of biosynthetic gene clusters in the wider context of metabolism, populations and epigenetics. Finally, we synthesise our observations to propose a new hypothesis for biosynthetic gene cluster formation in plants.

Using interdisciplinary, phylogeny-guided approaches to understand the evolution of plant metabolism

To cope with relentless environmental pressures, plants produce an arsenal of structurally diverse chemicals, often called specialized metabolites. These lineage-specific compounds are derived from the simple building blocks made by ubiquitous core metabolic pathways. Although the structures of many specialized metabolites are known, the underlying metabolic pathways and the evolutionary events that have shaped the plant chemical diversity landscape are only beginning to be understood. However, with the advent of multi-omics data sets and the relative ease of studying pathways in previously intractable non-model species, plant specialized metabolic pathways are now being systematically identified. These large datasets also provide a foundation for comparative, phylogeny-guided studies of plant metabolism. Comparisons of metabolic traits and features like chemical abundances, enzyme activities, or gene sequences from phylogenetically diverse plants provide insights into how metabolic pathways evolved. This review highlights the power of studying evolution through the lens of comparative biochemistry, particularly how placing metabolism into a phylogenetic context can help a researcher identify the metabolic innovations enabling the evolution of structurally diverse plant metabolites.

The origin and evolution of RNase T2 family and gametophytic self-incompatibility system in plants

Ribonuclease (RNase) T2 genes are found widely in both eukaryotes and prokaryotes, and genes from this family have been revealed to have various functions in plants. In particular, S-RNase is known to be the female determinant in the S-RNase-based gametophytic self-incompatibility (GSI) system. However, the origin and evolution of the RNase T2 gene family and GSI system are not well understood. In this study, 785 RNase T2 genes were identified in 81 sequenced plant genomes representing broad-scale diversity and divided into three subgroups (Class I, II, and III) based on phylogenetic and synteny network analysis. Class I was found to be of ancient origin and to emerge in green algae, Class II was shown to originate with the appearance of angiosperms, while Class III was discovered to be eudicot-specific. Each of the three major classes could be further classified into several subclasses of which some subclasses were found to be lineage-specific. Furthermore, duplication, deletion, or inactivation of the S/S-like-locus was revealed to be linked to repeated loss and gain of self-incompatibility in different species from distantly related plant families with GSI. Finally, the origin and evolutionary history of S-locus in Rosaceae species was unraveled with independent loss and gain of S-RNase occurred in different subfamilies of Rosaceae. Our findings provide insights into the origin and evolution of the RNase T2 family and the GSI system in plants.

Deletion and tandem duplications of biosynthetic genes drive the diversity of triterpenoids in Aralia elata

Araliaceae species produce various classes of triterpene and triterpenoid saponins, such as the oleanane-type triterpenoids in Aralia species and dammarane-type saponins in Panax, valued for their medicinal properties. The lack of genome sequences of Panax relatives has hindered mechanistic insight into the divergence of triterpene saponins in Araliaceae. Here, we report a chromosome-level genome of Aralia elata with a total length of 1.05 Gb. The loss of 12 exons in the dammarenediol synthase (DDS)-encoding gene in A. elata after divergence from Panax might have caused the lack of dammarane-type saponin production, and a complementation assay shows that overexpression of the PgDDS gene from Panax ginseng in callus of A. elata recovers the accumulation of dammarane-type saponins. Tandem duplication events of triterpene biosynthetic genes are common in the A. elata genome, especially for AeCYP72As, AeCSLMs, and AeUGT73s, which function as tailoring enzymes of oleanane-type saponins and aralosides. More than 13 aralosides are de novo synthesized in Saccharomyces cerevisiae by overexpression of these genes in combination. This study sheds light on the diversity of saponins biosynthetic pathway in Araliaceae and will facilitate heterologous bioproduction of aralosides.

GURFAP: A Platform for Gene Function Analysis in Glycyrrhiza Uralensis

Glycyrrhiza uralensis (Licorice), which belongs to Leguminosae, is famous for the function of pharmacologic action and natural sweetener with its dried roots and rhizomes. In recent years, the whole-genome sequence of G. uralensis has been completed, which will help to lay the foundation for the study of gene function. Here, we integrated the available genomic and transcriptomic data of G. uralensis and constructed the G. uralensis gene co-expression network. We then annotated gene functions of G. uralensis via aligning with public databases. Furthermore, gene families of G. uralensis were predicted by tools including iTAK (Plant Transcription factor and Protein kinase Identifier and Classifier), HMMER (hidden Markov models), InParanoid, and PfamScan. Finally, we constructed a platform for gene function analysis in G. uralensis (GURFAP, www.gzybioinfoormatics.cn/GURFAP). For analyzed and predicted gene function, we introduced various tools including BLAST (Basic local alignment search tool), GSEA (Gene set enrichment analysis), Motif, Heatmap, and JBrowse. Our analysis based on this platform indicated that the biosynthesis of glycyrrhizin might be regulated by MYB and bHLH. We also took CYP88D6, CYP72A154, and bAS gene in the synthesis pathway of glycyrrhizin as examples to demonstrate the reliability and availability of our platform. Our platform GURFAP will provide convenience for researchers to mine the gene function of G. uralensis and thus discover more key genes involved in the biosynthetic pathway of active ingredients.

Unlocking plant metabolic diversity: A (pan)-genomic view

Plants produce a remarkable diversity of structurally and functionally diverse natural chemicals that serve as adaptive compounds throughout their life cycles. However, unlocking this metabolic diversity is significantly impeded by the size, complexity, and abundant repetitive elements of typical plant genomes. As genome sequencing becomes routine, we anticipate that links between metabolic diversity and genetic variation will be strengthened. In addition, an ever-increasing number of plant genomes have revealed that biosynthetic gene clusters are not only a hallmark of microbes and fungi; gene clusters for various classes of compounds have also been found in plants, and many are associated with important agronomic traits. We present recent examples of plant metabolic diversification that have been discovered through the exploration and exploitation of various genomic and pan-genomic data. We also draw attention to the fundamental genomic and pan-genomic basis of plant chemodiversity and discuss challenges and future perspectives for investigating metabolic diversity in the coming pan-genomics era.

Integration of genomics, transcriptomics and metabolomics identifies candidate loci underlying fruit weight in loquat

Fruit weight is an integral part of fruit-quality traits and directly influences commodity values and economic returns of fruit crops. Despite its importance, the molecular mechanisms underlying fruit weight remain understudied, especially for perennial fruit tree crops such as cultivated loquat (Eriobotrya japonica Lindl.). Auxin is known to regulate fruit development, whereas its role and metabolism in fruit development remain obscure in loquat. In this study, we applied a multi-omics approach, integrating whole-genome resequencing-based quantitative trait locus (QTL) mapping with an F1 population, population genomics analysis using germplasm accessions, transcriptome analysis, and metabolic profiling to identify the genomic regions potentially associated with fruit weight in loquat. We identified three major loci associated with fruit weight, supported by both QTL mapping and comparative genomic analysis between small- and big-fruited loquat cultivars. Comparison between two genotypes with contrasting fruit weight performance through transcriptomic and metabolic profiling revealed an important role of auxin in regulating fruit development, especially at the fruit enlarging stage. The multi-omics approach identified two homologs of ETHYLENE INSENSITIVE 4 (EjEIN4) and TORNADO 1 (EjTRN1) as promising candidates controlling fruit weight. Moreover, three single nucleotide polymorphism (SNP) markers were closely associated with fruit weight. Results from this study provided insights from multiple perspectives into the genetic and metabolic controls of fruit weight in loquat. The candidate genomic regions, genes, and sequence variants will facilitate understanding the molecular basis of fruit weight and lay a foundation for future breeding and manipulation of fruit weight in loquat.

Multi-Omics-Based Discovery of Plant Signaling Molecules

Plants produce numerous structurally and functionally diverse signaling metabolites, yet only relatively small fractions of which have been discovered. Multi-omics has greatly expedited the discovery as evidenced by increasing recent works reporting new plant signaling molecules and relevant functions via integrated multi-omics techniques. The effective application of multi-omics tools is the key to uncovering unknown plant signaling molecules. This review covers the features of multi-omics in the context of plant signaling metabolite discovery, highlighting how multi-omics addresses relevant aspects of the challenges as follows: (a) unknown functions of known metabolites; (b) unknown metabolites with known functions; (c) unknown metabolites and unknown functions. Based on the problem-oriented overview of the theoretical and application aspects of multi-omics, current limitations and future development of multi-omics in discovering plant signaling metabolites are also discussed.

Gibberellin Induced Transcriptome Profiles Reveal Gene Regulation of Loquat Flowering

Flowering is an integral part of the life cycle of flowering plants, which is essential for plant survival and crop production. Most woody fruit trees such as apples and pears bloom in spring, but loquat blooms in autumn and winter. Gibberellin (GA) plays a key role in the regulation of plant flower formation. In this study, we sprayed loquat plants with exogenous GA₃, which resulted in vigorous vegetative growth rather than floral bud formation. We then performed a comprehensive RNA-seq analysis on GA₃-treated and control-treated leaves and buds over three time periods to observe the effects of exogenous GA₃ application on floral initiation and development. The results showed that 111 differentially expressed genes (DEGs) and 563 DEGs were down-regulated, and 151 DEGs and 506 DEGs were up-regulated in buds and leaves, respectively, upon treatment with GA₃. Among those that are homologs of the DELLA-mediated GA signal pathway genes, some may be involved in the positive regulation of flower development, including EjWRKY75, EjFT, EjSOC1, EjAGL24, EjSPL, EjLFY, EjFUL, and EjAP1; while some may be involved in the negative regulation of flower development, including EjDELLA, EjMYC3, EjWRKY12, and EjWRKY13. Finally, by analyzing the co-expression of DEGs and key floral genes EjSOC1s, EjLFYs, EjFULs, EjAP1s, 330 candidate genes that may be involved in the regulation of loquat flowering were screened. These genes belong to 74 gene families, including Cyclin_C, Histone, Kinesin, Lipase_GDSL, MYB, P450, Pkinase, Tubulin, and ZF-HD_dimer gene families. These findings provide new insights into the regulation mechanism of loquat flowering.

Genotyping of polyploid plants using quantitative PCR: application in the breeding of white-fleshed triploid loquats (Eriobotrya japonica)

BACKGROUND

Ploidy manipulation is effective in seedless loquat breeding, in which flesh color is a key agronomic and economic trait. Few techniques are currently available for detecting the genotypes of polyploids in plants, but this ability is essential for most genetic research and molecular breeding.

RESULTS

We developed a system for genotyping by quantitative PCR (qPCR) that allowed flesh color genotyping in multiple tetraploid and triploid loquat varieties (lines). The analysis of 13 different ratios of DNA mixtures between two homozygous diploids (AA and aa) showed that the proportion of allele A has a high correlation (R² = 0.9992) with parameter b [b = a₁/(a₁ + a₂)], which is derived from the two normalized allele signals (a₁ and a₂) provided by qPCR. Cluster analysis and variance analysis from simulating triploid and tetraploid hybrids provided completely correct allelic configurations. Four genotypes (AAA, AAa, Aaa, aaa) were found in triploid loquats, and four (AAAA, AAAa, AAaa, Aaaa; absence of aaaa homozygotes) were found in tetraploid loquats. DNA markers analysis showed that the segregation of flesh color in all F₁ hybrids conformed to Mendel's law. When tetraploid B431 was the female parent, more white-fleshed triploids occurred among the progeny.

CONCLUSIONS

qPCR can detect the flesh color genotypes of loquat polyploids and provides an alternative method for analyzing polyploid genotype and breeding, dose effects and allele-specific expression.

Evaluation of phenolic composition and antioxidant properties of different varieties of Chinese citrus

Citrus peels have health-promoting effects and are a rich source of antioxidant substances. This study evaluated the compositions of phenolic compounds and antioxidant activities in the peels of 52 citrus varieties with consistent planting time and management. The highest levels of total phenols (72.95 ± 37.60 mg/g DW) and total flavonoids (71.43 ± 37.64 mg/g DW) were found in mandarin. The highest phenolic acid content (18.78 ± 0.38 mg/g DW), dominated by protocatechuic acid, was found in kumquat. The antioxidant potency composite index was 6.23-94.56, suggesting mandarin varieties HJ, TWPG, TTPG, AY28, BZH and TCJC had the highest antioxidant activity. Statistics analysis indicated phenolic compounds and antioxidant activity were positively correlated. Principal component analysis and hierarchical cluster analysis suggested a strong relationship between phenolic compound composition and genetic background. This study indicated significant differences in the biological properties of various types of citrus peels; which are valuable for future utilization and research of citrus peels.

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.

EjRAV1/2 Delay Flowering Through Transcriptional Repression of EjFTs and EjSOC1s in Loquat

Most species in Rosaceae usually need to undergo several years of juvenile phase before the initiation of flowering. After 4-6 years' juvenile phase, cultivated loquat (Eriobotrya japonica), a species in Rosaceae, enters the reproductive phase, blooms in the autumn and sets fruits during the winter. However, the mechanisms of the transition from a seedling to an adult tree remain obscure in loquat. The regulation networks controlling seasonal flowering are also largely unknown. Here, we report two RELATED TO ABI3 AND VP1 (RAV) homologs controlling juvenility and seasonal flowering in loquat. The expressions of EjRAV1/2 were relatively high during the juvenile or vegetative phase and low at the adult or reproductive phase. Overexpression of the two EjRAVs in Arabidopsis prolonged (about threefold) the juvenile period by repressing the expressions of flowering activator genes. Additionally, the transformed plants produced more lateral branches than the wild type plants. Molecular assays revealed that the nucleus localized EjRAVs could bind to the CAACA motif of the promoters of flower signal integrators, EjFT1/2, to repress their expression levels. These findings suggest that EjRAVs play critical roles in maintaining juvenility and repressing flower initiation in the early life cycle of loquat as well as in regulating seasonal flowering. Results from this study not only shed light on the control and maintenance of the juvenile phase, but also provided potential targets for manipulation of flowering time and accelerated breeding in loquat.