The honeysuckle genome provides insight into the molecular mechanism of carotenoid metabolism underlying dynamic flower coloration

Construction of a high-efficiency GjCCD4a mutant and its application for de novo biosynthesis of five crocins in Escherichia coli

Crocins are bioactive natural products that rarely exist in plants. High costs and resource shortage severely limit its development and application. Synthetic biology studies on crocins are of considerable global interest. However, the lack of high-efficiency genetic tools and complex cascade biocatalytic systems have substantially hindered progress in crocin biosynthesis-related research. Based on mutagenesis, a high-efficiency GjCCD4a mutant (N212m) was constructed with a catalytic efficiency that was 25.08-fold higher than that of the wild-type. Solubilized GjCCD4a was expressed via fusion with an MBP tag. Moreover, N212m and ten other genes were introduced into Escherichia coli for the de novo biosynthesis of five crocins. The engineered E57 strain produced crocins III and V with a total yield of 11.50 mg/L, and the E579 strain produced crocins I-V with a total output of 8.43 mg/L at shake-flask level. This study identified a marvelous genetic element (N212m) for crocin biosynthesis and achieved its de novo biosynthesis in E. coli using glucose. This study provides a reference for the large-scale production of five crocins using E. coli cell factories.

Haplotype-resolved chromosome-level genome assembly of Ehretia macrophylla

Ehretia macrophylla Wall, known as wild loquat, is an ecologically, economically, and medicinally significant tree species widely grown in China, Japan, Vietnam, and Nepal. In this study, we have successfully generated a haplotype-resolved chromosome-scale genome assembly of E. macrophylla by integrating PacBio HiFi long-reads, Illumina short-reads, and Hi-C data. The genome assembly consists of two haplotypes, with sizes of 1.82 Gb and 1.58 Gb respectively, and contig N50 lengths of 28.11 Mb and 21.57 Mb correspondingly. Additionally, 99.41% of the assembly was successfully anchored into 40 pseudo-chromosomes. We predicted 58,886 protein-coding genes, of which 99.60% were functionally annotated from databases. We furthermore detected 2.65 Gb repeat sequences, 659,290 rRNAs, 4,931 tRNAs and 4,688 other ncRNAs. The high-quality assembly of the genome offers a solid basis for furthering the fields of molecular breeding and functional genomics of E. macrophylla.

The chromosome-level genome and functional database accelerate research about biosynthesis of secondary metabolites in Rosa roxburghii

Rosa roxburghii Tratt, a valuable plant in China with long history, is famous for its fruit. It possesses various secondary metabolites, such as L-ascorbic acid (vitamin C), alkaloids and poly saccharides, which make it a high nutritional and medicinal value. Here we characterized the chromosome-level genome sequence of R. roxburghii, comprising seven pseudo-chromosomes with a total size of 531 Mb and a heterozygosity of 0.25%. We also annotated 45,226 coding gene loci after masking repeat elements. Orthologs for 90.1% of the Complete Single-Copy BUSCOs were found in the R. roxburghii annotation. By aligning with protein sequences from public platform, we annotated 85.89% genes from R. roxburghii. Comparative genomic analysis revealed that R. roxburghii diverged from Rosa chinensis approximately 5.58 to 13.17 million years ago, and no whole-genome duplication event occurred after the divergence from eudicots. To fully utilize this genomic resource, we constructed a genomic database RroFGD with various analysis tools. Otherwise, 69 enzyme genes involved in L-ascorbate biosynthesis were identified and a key enzyme in the biosynthesis of vitamin C, GDH (L-Gal-1-dehydrogenase), is used as an example to introduce the functions of the database. This genome and database will facilitate the future investigations into gene function and molecular breeding in R. roxburghii.

Integrative analysis of transcriptome and target metabolites uncovering flavonoid biosynthesis regulation of changing petal colors in Nymphaea 'Feitian 2'

BACKGROUND

Nymphaea (waterlily) is known for its rich colors and role as an important aquatic ornamental plant globally. Nymphaea atrans and some hybrids, including N. 'Feitian 2,' are more appealing due to the gradual color change of their petals at different flower developmental stages. The petals of N. 'Feitian 2' gradually change color from light blue-purple to deep rose-red throughout flowering. The mechanism of the phenomenon remains unclear.

RESULTS

In this work, flavonoids in the petals of N. 'Feitian 2' at six flowering stages were examined to identify the influence of flavonoid components on flower color changes. Additionally, six cDNA libraries of N. 'Feitian 2' over two blooming stages were developed, and the transcriptome was sequenced to identify the molecular mechanism governing petal color changes. As a result, 18 flavonoid metabolites were identified, including five anthocyanins and 13 flavonols. Anthocyanin accumulation during flower development is the primary driver of petal color change. A total of 12 differentially expressed genes (DEGs) in the flavonoid biosynthesis pathway were uncovered, and these DEGs were significantly positively correlated with anthocyanin accumulation. Six structural genes were ultimately focused on, as their expression levels varied significantly across different flowering stages. Moreover, 104 differentially expressed transcription factors (TFs) were uncovered, and three MYBs associated with flavonoid biosynthesis were screened. The RT-qPCR results were generally aligned with high-throughput sequencing results.

CONCLUSIONS

This research offers a foundation to clarify the mechanisms underlying changes in the petal color of waterlilies.

Coping with alpine habitats: genomic insights into the adaptation strategies of Triplostegia glandulifera (Caprifoliaceae)

How plants find a way to thrive in alpine habitats remains largely unknown. Here we present a chromosome-level genome assembly for an alpine medicinal herb, Triplostegia glandulifera (Caprifoliaceae), and 13 transcriptomes from other species of Dipsacales. We detected a whole-genome duplication event in T. glandulifera that occurred prior to the diversification of Dipsacales. Preferential gene retention after whole-genome duplication was found to contribute to increasing cold-related genes in T. glandulifera. A series of genes putatively associated with alpine adaptation (e.g. CBFs, ERF-VIIs, and RAD51C) exhibited higher expression levels in T. glandulifera than in its low-elevation relative, Lonicera japonica. Comparative genomic analysis among five pairs of high- vs low-elevation species, including a comparison of T. glandulifera and L. japonica, indicated that the gene families related to disease resistance experienced a significantly convergent contraction in alpine plants compared with their lowland relatives. The reduction in gene repertory size was largely concentrated in clades of genes for pathogen recognition (e.g. CNLs, prRLPs, and XII RLKs), while the clades for signal transduction and development remained nearly unchanged. This finding reflects an energy-saving strategy for survival in hostile alpine areas, where there is a tradeoff with less challenge from pathogens and limited resources for growth. We also identified candidate genes for alpine adaptation (e.g. RAD1, DMC1, and MSH3) that were under convergent positive selection or that exhibited a convergent acceleration in evolutionary rate in the investigated alpine plants. Overall, our study provides novel insights into the high-elevation adaptation strategies of this and other alpine plants.

Reduction of flavonoid content in honeysuckle via Erysiphe lonicerae-mediated inhibition of three essential genes in flavonoid biosynthesis pathways

Honeysuckle, valued for its wide-ranging uses in medicine, cuisine, and aesthetics, faces a significant challenge in cultivation due to powdery mildew, primarily caused by the Erysiphe lonicerae pathogen. The interaction between honeysuckle and E. lonicerae, especially concerning disease progression, remains insufficiently understood. Our study, conducted in three different locations, found that honeysuckle naturally infected with E. lonicerae showed notable decreases in total flavonoid content, with reductions of 34.7%, 53.5%, and 53.8% observed in each respective site. Controlled experiments supported these findings, indicating that artificial inoculation with E. lonicerae led to a 20.9% reduction in flavonoid levels over 21 days, worsening to a 54.8% decrease by day 42. Additionally, there was a significant drop in the plant's total antioxidant capacity, reaching an 81.7% reduction 56 days after inoculation. Metabolomic analysis also revealed substantial reductions in essential medicinal components such as chlorogenic acid, luteolin, quercetin, isoquercetin, and rutin. Investigating gene expression revealed a marked decrease in the relative expression of the LjPAL1 gene, starting as early as day 7 post-inoculation and falling to a minimal level (fold change = 0.29) by day 35. This trend was mirrored by a consistent reduction in phenylalanine ammonia-lyase activity in honeysuckle through the entire process, which decreased by 72.3% by day 56. Further analysis showed significant and sustained repression of downstream genes LjFNHO1 and LjFNGT1, closely linked to LjPAL1. We identified the mechanism by which E. lonicerae inhibits this pathway and suggest that E. lonicerae may strategically weaken the honeysuckle's disease resistance by targeting key biosynthetic pathways, thereby facilitating further pathogen invasion. Based on our findings, we recommend two primary strategies: first, monitoring medicinal constituent levels in honeysuckle from E. lonicerae-affected areas to ensure its therapeutic effectiveness; and second, emphasizing early prevention and control measures against honeysuckle powdery mildew due to the persistent decline in crucial active compounds.

Exploration of exogenous chlorogenic acid as a potential plant stimulant: enhancing physiochemical properties in Lonicera japonica

In this study, we applied exogenous chlorogenic acid (CGA) to Lonicera japonica (L. japonica) leaves via foliar sprays every Monday, Wednesday, and Friday for a period of 12 months. Our continuous monitoring over this period revealed a consistent increase in flavonoid levels from the second to the tenth month following the commencement of CGA treatment. This was accompanied by a notable upregulation in the expression of four secondary metabolite-related enzyme genes: LjPAL1, LjPAL2, LjPAL3, and LjISY1. Concurrently, there was a significant enhancement in the total activity of the enzyme phenylalanine ammonia-lyase. The total antioxidant capacity of the plants also showed a marked increase from the third to the seventh month post-treatment initiation, subsequently stabilizing. This increase was also reflected in the elevated activities of key antioxidant enzymes: peroxidase, polyphenol oxidase, and superoxide dismutase. Furthermore, the treatment notably enhanced various indicators of nutrient growth, such as total protein content, total sugar content, and leaf area. Notably, the relative expression of LjTF1, a kind of BZIP transcription factor gene known for its extensive regulatory effects, showed a significant and sustained increase after the start of exogenous CGA treatment. Subsequent metabolomic analysis revealed significant changes in L. japonica metabolites. Specifically, 172 differentially expressed metabolites (DEMs) showed a notable increase (Fold > 1), predominantly in pathways related to nutrient metabolism such as carbohydrate, amino acid, and energy metabolism. Notably, some of the highly expressed DEMs (Fold > 4) are key antioxidants and medicinal components in L. japonica. The experimental findings were in alignment with the metabolomics analysis, indicating that exogenous CGA can act as a stimulant for L. japonica. It promotes the significant accumulation of certain secondary metabolites, enhances nutritive growth, and boosts the plant's total antioxidant capacity.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12298-024-01435-8.

Integrated metabolic profiling and transcriptome analysis of Lonicera japonica flowers for chlorogenic acid, luteolin and endogenous hormone syntheses

The active ingredients of many medicinal plants are the secondary metabolites associated with the growth period. Lonicera japonica Thunb. is an important traditional Chinese medicine, and the flower development stage is an important factor that influences the quality of medicinal ingredients. In this study, transcriptomics and metabolomics were performed to reveal the regulatory mechanism of secondary metabolites during flowering of L. japonica. The results showed that the content of chlorogenic acid (CGA) and luteolin gradually decreased from green bud stage (Sa) to white flower stage (Sc), especially from white flower bud stage (Sb) to Sc. Most of the genes encoding the crucial rate-limiting enzymes, including PAL, C4H, HCT, C3'H, F3'H and FNSII, were down-regulated in three comparisons. Correlation analysis identified some members of the MYB, AP2/ERF, bHLH and NAC transcription factor families that are closely related to CGA and luteolin biosynthesis. Furthermore, differentially expressed genes (DEGs) involved in hormone biosynthesis, signalling pathways and flowering process were analysed in three flower developmental stage.

ERF5.1 modulates carotenoid accumulation by interacting with CCD4.1 in Lycium

Carotenoids are important natural pigments and have medical and health functions for humans. Carotenoid cleavage dioxygenase 4 (CCD4) and ethylene responsive factor (ERF) participate in carotenoid metabolism, but their roles in Lycium have not been discovered. Here, we annotated LbCCDs from the Lycium reference genome and found that LbCCD4.1 expression was significantly correlated with the carotenoid metabolites during Lycium five fruit developmental stages. Over-expression of LbCCD4.1 in NQ's leaves resulted in a series of significantly lower contents of carotenoid metabolites, including β-carotene and β-cryptoxanthin. Moreover, LbERF5.1, a transcription factor belonging to the ERF family that was located in the nucleus, was isolated. Significant reductions in the carotenoids, especially lutein, violaxanthin and their derivatives, were observed in over-expressing ERF5.1 transgenic NQ's leaves. Over-expression or virus-induced gene silencing of LbERF5.1 in NQ's leaves induced a consistent up- or down-expression, respectively, of LbCCD4.1. Furthermore, yeast one-hybrid and dual-luciferase reporter assays showed that ERF5.1 interacted with the promoter of CCD4.1 to increase its expression, and LbERF5.1 could bind to any one of the three predicted binding sites in the promoter of LbCCD4.1. A transcriptome analysis of LbERF5.1 and LbCCD4.1 over-expressed lines showed similar global transcript expression, and geranylgeranyl diphosphate synthase, phytoene synthase, lycopene δ-cyclase cytochrome, cytochrome P450-type monooxygenase 97A, cytochrome P450-type monooxygenase 97C, and zeaxanthin epoxidase in the carotenoid biosynthesis pathway were differentially expressed. In summary, we uncovered a novel molecular mechanism of carotenoid accumulation that involved an interaction between ERF5.1 and CCD4.1, which may be used to enhance carotenoid in Lycium.

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.

Magnesium accelerates changes in the fruit ripening and carotenoid accumulation in Satsuma Mandarin pulp

This study aims to further examine the effect of Magnesium (Mg) application on fruit quality and carotenoid metabolism in Satsuma mandarin pulp. For this, a field experiment was using 20-year-old Satsuma mandarin (C. unshiu Marc.) for two treatment; (1) CK treatment (without Mg), (2) Mg fertilizer treatment (200 g MgO plant-1). Compared with CK, Mg treatment substantially raised the Mg content in pulp at 90 to 150 DAF (the fruit expansion period), increasing by 15.69%-21.74%. Mg treatment also increased fruit TSS content by 15.84% and 9.88%, decreased fruit TA content in by 34.25% and 33.26% at 195 DAF and 210 DAF (the fruit ripening period). Moreover, at 120 to 195 DAF, Mg treatment significantly increased the levels of lutein, β-cryptoxanthin, zeaxanthin and violaxanthin in the pulp. This can be explained by the increased expression of important biosynthetic genes, including CitPSY, CitPDS, CitLCYb1, CitLCYb2, CitLCYe, CitHYb, and CitZEP, that played a role in altering the carotenoid composition. The findings of this research offer a novel approach for augmenting both the economic and nutritional worth of citrus fruits.

Chromosome-scale genome assembly provides insights into flower coloration mechanisms of Canna indica

Canna indica, which produce conspicuous and colorful flowers, are widely appreciated as ornamental plants. We used Pacific Biosciences sequencing (PacBio) and chromosome conformation capture (Hi-C) genome scaffolding to build a high-quality chromosome-scale genome assembly of C. indica and the genome assembly was 821Mb with a contig N50 of 48Mb assembled into nine chromosomes. The genome of C. indica was predicted to contain 31,130 genes and 30,816 genes were functionally annotated. Genome annotation identified 522 Mb (63.59 %) as repetitive sequences. Genome evolution analysis showed that whole-genome duplication occurred 53.4 million years ago. Transcriptome analysis revealed that petal coloration was linked with the expression of genes encoding enzymes involved in anthocyanin biosynthesis, carotenoid biosynthesis, and the methylerythritol phosphate (MEP) pathway. Furthermore, modules of co-expressed genes and hub genes were identified via weighted gene co-expression network analysis. These results suggested that, in Canna indica, deep red petal coloration was regulated by CHS2 and yellow petal coloration was associated with expression of ARF6 and NAC14. Considered together, the current study revealed a high-quality reference genome which may provide new insights into the molecular basis of flower coloration in Canna indica and help enhance the conservation and breeding of ornamental plants in general.

Genome-wide identification and in-silico expression analysis of carotenoid cleavage oxygenases gene family in Oryza sativa (rice) in response to abiotic stress

Rice constitutes a foundational cereal and plays a vital role in the culinary sector. However, the detriments of abiotic stress on rice quality and productivity are noteworthy. Carotenoid cleavage oxygenases (CCO) hold vital importance as they enable the particular breakdown of carotenoids and significantly contribute towards the growth and response to abiotic stress in rice. Due to the insufficient information regarding rice CCOs and their potential role in abiotic stress, their utilization in stress-resistant genetic breeding remains limited. The current research identified 16 CCO genes within the Oryza sativa japonica group. These OsCCO genes can be bifurcated into three categories based on their conserved sequences: NCEDs (9-Cis-epoxycarotenoid dioxygenases), CCDs (Carotenoid cleavage dioxygenases) and CCD-like (Carotenoid cleavage dioxygenases-like). Conserved motifs were found in the OsCCO gene sequence via MEME analysis and multiple sequence alignment. Stress-related cis-elements were detected in the promoter regions of OsCCOs genes, indicating their involvement in stress response. Additionally, the promoters of these genes had various components related to plant light, development, and hormone responsiveness, suggesting they may be responsive to plant hormones and involved in developmental processes. MicroRNAs play a pivotal role in the regulation of these 16 genes, underscoring their significance in rice gene regulation. Transcriptome data analysis suggests a tissue-specific expression pattern for rice CCOs. Only OsNCED6 and OsNCED10 significantly up-regulated during salt stress, as per RNA seq analyses. CCD7 and CCD8 levels were also higher in the CCD group during the inflorescence growth stage. This provides insight into the function of rice CCOs in abiotic stress response and identifies possible genes that could be beneficial for stress-resistant breeding.

The Functional Characteristics and Soluble Expression of Saffron CsCCD2

Crocins are important natural products predominantly obtained from the stigma of saffron, and that can be utilized as a medicinal compound, spice, and colorant with significant promise in the pharmaceutical, food, and cosmetic industries. Carotenoid cleavage dioxygenase 2 (CsCCD2) is a crucial limiting enzyme that has been reported to be responsible for the cleavage of zeaxanthin in the crocin biosynthetic pathway. However, the catalytic activity of CsCCD2 on β-carotene/lycopene remains elusive, and the soluble expression of CsCCD2 remains a big challenge. In this study, we reported the functional characteristics of CsCCD2, that can catalyze not only zeaxanthin cleavage but also β-carotene and lycopene cleavage. The molecular basis of the divergent functionality of CsCCD2 was elucidated using bioinformatic analysis and truncation studies. The protein expression optimization results demonstrated that the use of a maltose-binding protein (MBP) tag and the optimization of the induction conditions resulted in the production of more soluble protein. Correspondingly, the catalytic efficiency of soluble CsCCD2 was higher than that of the insoluble one, and the results further validated its functional verification. This study not only broadened the substrate profile of CsCCD2, but also achieved the soluble expression of CsCCD2. It provides a firm platform for CsCCD2 crystal structure resolution and facilitates the synthesis of crocetin and crocins.

Combinatorial transcriptomics and metabolomics analysis reveals the effects of the harvesting stages on the accumulation of phenylpropanoid metabolites in Lonicera japonica

The flower buds of Lonicera japonica are widely used for its high medicinal value. It is reported that the accumulation of phenylpropanoids in the buds of L. japonica is affected by the stage at which it is harvested. However, the changes of active components and the underlying mechanisms in flower buds at different harvesting stages have not been reported. Integrative analyses of transcriptomics and metabolomics was used to explore the underlying mechanism of harvesting stages (green bud, GB; and white bud, WB) on the phenylpropanoids metabolites accumulation in L. japonica . The result showed that 3735 differentially expressed genes were identified, and the genes related to glycolysis/gluconeogenesis and phenylalanine biosynthesis pathway were significantly upregulated in GB stage. A total of 510 differential metabolites were identified in GB stage. Among them, 14 phenylpropanoids were changed during the GB and WB, seven of which increased in GB, including caffeic acid, sauchinone, coniferin, secoisolariciresinol diglucoside, scopolin, methyl cinnamate, chlorogenic acid, 7-hydroxycoumarin, while others such as sibiricose A6, coumarin, eleutheroside E decreased. Further correlation analysis showed that the unigenes for CSE, CAD, bg1, ADH, ALDH, DLAT and ENO significantly correlated with the 10 phenylpropanoid. The above results would provide basic data for the selection of harvesting stages in the production of L. japonica .

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

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

A Comprehensive Quality Analysis of Different Colors of Medicinal and Edible Honeysuckle

Honeysuckle (the dried flower bud or opening flower of Lonicera japonica Thunb.), a medicinal and edible substance, has is greatly popular among consumers for its remarkable health effects, such as antioxidant, antibacterial, and anti-inflammatory effects. However, due to the influences of processing methods, storage conditions, and other factors, honeysuckles show different colors which can directly reflect the quality and the price on the market. In order to comprehensively compare the quality of different colors, 55 batches of honeysuckle samples were collected and analyzed. Their color parameters, chlorophyll content (chl), total phenol content (TPC), total flavonoid content (TFC), antioxidant activity (AA), main active compounds, and metabolites were measured. As a result, the initial green-white (GW) samples, a kind of highest-quality honeysuckle, had the smallest a* value, largest h*, chl, TPC, TFC, and AA values, and highest content of chlorogenic acid and cynaroside. There was a significant difference between GW samples and a series of discolored samples. As the color darkened or lightened, the quality gradually decreased. The yellow-brown (YB) samples were of the worst quality and were no longer available for clinical and health purposes. A series of differential metabolites, such as quercetin-7-O-glucoside and secologanoside, could be used as important references to evaluate the quality of differently colored samples. The metabolic profile of honeysuckle provided new insights into the process of color change and laid a foundation for further honeysuckle quality control. The correlation results showed that the a* and h* values significantly affect the abovementioned quality indicators and the 10 main active compounds. In other words, the color difference could directly reflect the quality and clinical efficacy. Multiple regression analysis was carried out using combined L*, a*, and b* values to predict the quality of honeysuckle. This is the first time the quality of different color honeysuckle samples on the post-harvest link has been systematically compared and a demonstration of medicinal and edible substances with different colors has been provided.

Genome-wide analysis of the MADS-box gene family in Lonicera japonica and a proposed floral organ identity model

BACKGROUND

Lonicera japonica Thunb. is widely used in traditional Chinese medicine. Medicinal L. japonica mainly consists of dried flower buds and partially opened flowers, thus flowers are an important quality indicator. MADS-box genes encode transcription factors that regulate flower development. However, little is known about these genes in L. japonica.

RESULTS

In this study, 48 MADS-box genes were identified in L. japonica, including 20 Type-I genes (8 Mα, 2 Mβ, and 10 Mγ) and 28 Type-II genes (26 MIKCc and 2 MIKC*). The Type-I and Type-II genes differed significantly in gene structure, conserved domains, protein structure, chromosomal distribution, phylogenesis, and expression pattern. Type-I genes had a simpler gene structure, lacked the K domain, had low protein structure conservation, were tandemly distributed on the chromosomes, had more frequent lineage-specific duplications, and were expressed at low levels. In contrast, Type-II genes had a more complex gene structure; contained conserved M, I, K, and C domains; had highly conserved protein structure; and were expressed at high levels throughout the flowering period. Eleven floral homeotic MADS-box genes that are orthologous to the proposed Arabidopsis ABCDE model of floral organ identity determination, were identified in L. japonica. By integrating expression pattern and protein interaction data for these genes, we developed a possible model for floral organ identity determination.

CONCLUSION

This study genome-widely identified and characterized the MADS-box gene family in L. japonica. Eleven floral homeotic MADS-box genes were identified and a possible model for floral organ identity determination was also developed. This study contributes to our understanding of the MADS-box gene family and its possible involvement in floral organ development in L. japonica.

The first chromosome-level Fallopia multiflora genome assembly provides insights into stilbene biosynthesis

Fallopia multiflora (Thunb.) Harald, a vine belonging to the Polygonaceae family, is used in traditional medicine. The stilbenes contained in it have significant pharmacological activities in anti-oxidation and anti-aging. This study describes the assembly of the F. multiflora genome and presents its chromosome-level genome sequence containing 1.46 gigabases of data (with a contig N50 of 1.97 megabases), 1.44 gigabases of which was assigned to 11 pseudochromosomes. Comparative genomics confirmed that F. multiflora shared a whole-genome duplication event with Tartary buckwheat and then underwent different transposon evolution after separation. Combining genomics, transcriptomics, and metabolomics data to map a network of associated genes and metabolites, we identified two FmRS genes responsible for the catalysis of one molecule of p-coumaroyl-CoA and three molecules of malonyl-CoA to resveratrol in F. multiflora. These findings not only serve as the basis for revealing the stilbene biosynthetic pathway but will also contribute to the development of tools for increasing the production of bioactive stilbenes through molecular breeding in plants or metabolic engineering in microbes. Moreover, the reference genome of F. multiflora is a useful addition to the genomes of the Polygonaceae family.

A phylogenomic analysis of Lonicera and its bearing on the evolution of organ fusion

PREMISE

The ~140 species of Lonicera are characterized by variously fused leaves, bracteoles, and ovaries, making it a model system for studying the evolution and development of organ fusion. However, previous phylogenetic analyses, based mainly on chloroplast DNA markers, have yielded uncertain and conflicting results. A well-supported phylogeny of Lonicera will allow us to trace the evolutionary history of organ fusion.

METHODS

We inferred the phylogeny of Lonicera using restriction site-associated DNA sequencing (RADSeq), sampling all major clades and 18 of the 23 subsections. This provided the basis for inferring the evolution of five fusion-related traits.

RESULTS

RADSeq data yielded a well-resolved and well-supported phylogeny. The two traditionally recognized subgenera (Periclymenum and Chamaecerasus), three of the four sections (Isoxylosteum, Coeloxylosteum, and Nintooa), and half of the subsections sampled were recovered as monophyletic. However, the large and heterogeneous section Isika was strongly supported as paraphyletic. Nintooa, a clade of ~22 mostly vine-forming species, including L. japonica, was recovered in a novel position, raising the possibility of cytonuclear discordance. We document the parallel evolution of fused leaves, bracteoles, and ovaries, with rare reversals. Most strikingly, complete cupules, in which four fused bracteoles completely enclose two unfused ovaries, arose at least three times. Surprisingly, these appear to have evolved directly from ancestors with free bracteoles instead of partial cupules.

CONCLUSIONS

We provide the most comprehensive and well-supported phylogeny of Lonicera to date. Our inference of multiple evolutionary shifts in organ fusion provides a solid foundation for in depth developmental and functional analyses.

A reference-grade genome assembly for Astragalus mongholicus and insights into the biosynthesis and high accumulation of triterpenoids and flavonoids in its roots

Astragalus membranaceus var. mongholicus (AMM), a member of the Leguminosae, is one of the most important medicinal plants worldwide. The dried roots of AMM have a wide range of pharmacological effects and are a traditional Chinese medicine. Here, we report the first chromosome-level reference genome of AMM, comprising nine pseudochromosomes with a total size of 1.47 Gb and 27 868 protein-encoding genes. Comparative genomic analysis reveals that AMM has not experienced an independent whole-genome duplication (WGD) event after the WGD event shared by the Papilionoideae species. Analysis of long terminal repeat retrotransposons suggests a recent burst of these elements at approximately 0.13 million years ago, which may explain the large size of the AMM genome. Multiple gene families involved in the biosynthesis of triterpenoids and flavonoids were expanded, and our data indicate that tandem duplication has been the main driver for expansion of these families. Among the expanded families, the phenylalanine ammonia-lyase gene family was primarily expressed in the roots of AMM, suggesting their roles in the biosynthesis of phenylpropanoid compounds. The functional versatility of 2,3-oxidosqualene cyclase genes in cluster III may play a critical role in the diversification of triterpenoids in AMM. Our findings provide novel insights into triterpenoid and flavonoid biosynthesis and can facilitate future research on the genetics and medical applications of AMM.

A lack of population structure characterizes the invasive Lonicera japonica in West Virginia and across eastern North America

Invasive plant species cause massive ecosystem damage globally, yet represent powerful case studies in population genetics and rapid adaptation to new habitats. The availability of digitized herbarium collections data, and the ubiquity of invasive species across the landscape make them highly accessible for studies of invasion history and population dynamics associated with their introduction, establishment, spread, and ecological interactions. Here we focus on Lonicera japonica, one of the most damaging invasive vine species in North America. We leveraged digitized collections data and contemporary field collections to reconstruct the invasion history and characterize patterns of genomic variation in the eastern USA, using a straightforward method for generating nucleotide polymorphism data and a recently published, chromosome-level genome for the species. We found an overall lack of population structure among sites in northern West Virginia, USA, as well as across sites in the central and eastern USA. Heterozygosity and population differentiation were both low based on Fst, analysis of molecular variance, principal components analysis, and cluster-based analyses. We also found evidence of high inbreeding coefficients and significant linkage disequilibrium, in line with the ability of this otherwise outcrossing, perennial species to propagate vegetatively. Our findings corroborate earlier studies based on allozyme data, and suggest that intentional, human-assisted spread explains the lack of population structure, as this species was planted for erosion control and as an ornamental, escaping cultivation repeatedly across the USA. Finally, we discuss how plant invasion genomics can be incorporated into experiential undergraduate education as a way to integrate teaching and research.

Light plays a critical role in the accumulation of chlorogenic acid in Lonicera macranthoides Hand.-Mazz

Light has important effects on plant metabolism. However, the relationship between the chlorogenic acid (CGA) content and light in plants remains unclear. Here, we investigated the effects of shading treatment on gene expression and CGA content in Lonicera macranthoides Hand.-Mazz. (LM), a widely used medicinal plant. A total of 1891 differentially expressed genes (DEGs) were obtained in flower buds and 819 in leaves in response to light in shading treatment compared to the control sample by RNA-Seq. After shading treatment, the content of CGA in LM leaves decreased significantly by 1.78-fold, the carotenoid content increased, and the soluble sugar and starch contents significantly decreased. WGCNA and the expression of related genes verified by qRT‒PCR revealed that CGA synthesis pathway enzyme genes form a co-expression network with genes for carbohydrate synthesis, photosynthesis, light signalling elements, and transcription factor genes (TFs) that affect the accumulation of CGA. Through a virus-induced gene silencing (VIGS) system and CGA assay in Nicotiana benthamiana (NB), we determined that downregulation of NbHY5 expression decreased the CGA content in NB leaves. In this study, we found that light provides energy and material for the accumulation of CGA in LM, and light affects the expression of CGA accumulation-related genes. Our results show that different light intensities have multiple effects on leaves and flower buds in LM and are able to coregulate LmHY5 expression and CGA synthesis.

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.

Whole-Genome Sequence, Assembly and Annotation of an Invasive Plant, Lonicera maackii (Amur Honeysuckle)

The invasive species Lonicera maackii (Amur Honeysuckle) is an increasing problem sweeping from the eastern United States toward the west, impacting normal forest development and animal survival across multiple taxa. Little is known about the genomics of this species, although a related invasive, Lonicera japonica, has been sequenced. Understanding the genomic foundation of the Lonicera maackii species could help us understand the biochemistry and life history that are the underpinnings of invasive success, as well as potential vulnerabilities and strengths which could guide research and development to control its spread. Here we present a draft, but high-quality, short-read whole-genome sequence, assembly, and annotation of Lonicera maackii, demonstrating that inexpensive and rapid short-read technologies can be successfully used in invasive species research. Despite being a short-read assembly, the genome length (7.93 × 10⁸) and completeness (estimated as 90.2-92.1% by BUSCO and Merqury) are close to the previously published chromosome-level sequencing of L. japonica. No bias, by means of a Gene Ontology analysis, was identified among missing BUSCOs. A duplication of the 5-enolpyruvylshikimate-3-phosphate (EPSP) synthase gene in both Lonicera species is identified, and the potential impact on controlling these invasive species is discussed. Future prospects for a diversity analysis of invasive species is also discussed.

The genome of single-petal jasmine (Jasminum sambac) provides insights into heat stress tolerance and aroma compound biosynthesis

Jasmine [Jasminum sambac (L.) Aiton] is a commercially important cultivated plant species known for its fragrant flowers used in the perfume industry, medicine and cosmetics. In the present study, we obtained a draft genome for the J. sambac cultivar 'Danbanmoli' (JSDB, a single-petal phenotype). We showed that the final genome of J. sambac was 520.80 Mb in size (contig N50 = 145.43 kb; scaffold N50 = 145.53 kb) and comprised 35,363 genes. Our analyses revealed that the J. sambac genome has undergone only an ancient whole-genome duplication (WGD) event. We estimated that the lineage that has given rise to J. sambac diverged from the lineage leading to Osmanthus fragrans and Olea europaea approximately 31.1 million years ago (Mya). On the basis of a combination of genomic and transcriptomic analyses, we identified 92 transcription factors (TFs) and 206 genes related to heat stress response. Base on a combination of genomic, transcriptomic and metabolomic analyses, a range of aroma compounds and genes involved in the benzenoid/phenylpropanoid and terpenoid biosynthesis pathways were identified. In the newly assembled J. sambac genome, we identified a total of 122 MYB, 122 bHLH and 69 WRKY genes. Our assembled J. sambac JSDB genome provides fundamental knowledge to study the molecular mechanism of heat stress tolerance, and improve jasmine flowers and dissect its fragrance.

Integration of genome and transcriptome reveal molecular regulation mechanism of early flowering trait in Prunus genus (Prunus mume and Prunus persica)

Flowering time is crucial for the survival and reproduction. Prunus genus belongs to the Rosaceae family and includes several hundred species of flowering trees and shrubs with important ornamental and economic values. However, the molecular mechanism underlying early flowering in Prunus genus is unclear. Here, we utilized the genome and transcriptome of P. mume and P. persica to explore the transcriptional regulation mechanism of early flowering. Comparative genomics found that genes accounting for 92.4% of the total P. mume genome and 91.2% of the total P. persica genome belonged to orthogroups. A total of 19,169 orthogroups were found between P. mume and P. persica, including 20,431 corresponding orthologues and 20,080 collinearity gene pairs. A total of 305 differentially expressed genes (DEGs) associated with early flowering were found, among which FT, TLI65, and NAP57 were identified as hub genes in the early flowering regulation pathway. Moreover, we identified twenty-five transcription factors (TFs) from nine protein families, including MADS-box, AP2/ERF, and MYB. Our results provide insights into the underlying molecular model of flowering time regulation in Prunus genus and highlight the utility of multi-omics in deciphering the properties of the inter-genus plants.

Genetic mapping, transcriptomic sequencing and metabolic profiling indicated a glutathione S-transferase is responsible for the red-spot-petals in Gossypium arboreum

A GST for red-spot-petals in Gossypium arboreum was identified as the candidate under the scope of multi-omics approaches. Colored petal spots are correlated with insect pollination efficiency in Gossypium species. However, molecular mechanisms concerning the formation of red spots on Gossypium arboreum flowers remain elusive. In the current study, the Shixiya1-R (SxyR, with red spots) × Shixiya1-W (SxyW, without red spots) segregating population was utilized to determine that the red-spot-petal phenotype was levered by a single dominant locus. This phenotype was expectedly related to the anthocyanin metabolites, wherein the cyanidin and delphinidin derivatives constituted the major partition. Subsequently, this dominant locus was narrowed to a 3.27 Mb range on chromosome 7 by genomic resequencing from the two parents and the two segregated progeny bulks that have spotted petals or not. Furthermore, differential expressed genes generated from the two bulks at either of three sequential flower developmental stages that spanning the spot formation were intersected with the annotated ones that allocated to the 3.27 Mb interval, which returned eight genes. A glutathione S-transferase-coding gene (Gar07G08900) out of the eight was the only one that exhibited simultaneously differential expression among all three developmental stages, and it was therefore considered to be the probable candidate. Finally, functional validation upon this candidate was achieved by the appearance of scattered petal spots with inhibited expression of Gar07G08900. In conclusion, the current report identified a key gene for the red spotted petal in G. arboreum under the scope of multi-omics approaches, such efforts and embedded molecular resources would benefit future applications underlying the flower color trait in cotton.

Natural products of medicinal plants: biosynthesis and bioengineering in post-genomic era

Globally, medicinal plant natural products (PNPs) are a major source of substances used in traditional and modern medicine. As we human race face the tremendous public health challenge posed by emerging infectious diseases, antibiotic resistance and surging drug prices etc., harnessing the healing power of medicinal plants gifted from mother nature is more urgent than ever in helping us survive future challenge in a sustainable way. PNP research efforts in the pre-genomic era focus on discovering bioactive molecules with pharmaceutical activities, and identifying individual genes responsible for biosynthesis. Critically, systemic biological, multi- and inter-disciplinary approaches integrating and interrogating all accessible data from genomics, metabolomics, structural biology, and chemical informatics are necessary to accelerate the full characterization of biosynthetic and regulatory circuitry for producing PNPs in medicinal plants. In this review, we attempt to provide a brief update on the current research of PNPs in medicinal plants by focusing on how different state-of-the-art biotechnologies facilitate their discovery, the molecular basis of their biosynthesis, as well as synthetic biology. Finally, we humbly provide a foresight of the research trend for understanding the biology of medicinal plants in the coming decades.

Haskap Berry Leaves (Lonicera caerulea L.)—The Favorable Potential of Medical Use

The presented research evaluates the medical use potential of Lonicera caerulea leaves, which are waste plants in cultivating berries. The study’s screening activity included the leaves of five varieties of Lonicera caerulea: Atut, Duet, Wojtek, Zojka, and Jugana. The microbiological analysis confirmed the safety of using Lonicera caerulea leaves without significant stabilization. Lonicera caerulea leaves standardization was carried out based on the results of the chromatographic analysis, and it showed differences in the contents of active compounds (loganic, chlorogenic and caffeic acids, and rutin), which are attributed to biological activity. For the Lonicera caerulea leaves varieties tested, the differences in the content of total polyphenol content, chlorophylls, and carotenoids were also confirmed. The screening of biological activity of five Lonicera caerulea leaf varieties was carried out concerning the possibility of inhibiting the activity of α-glucosidase, lipase, and hyaluronidase as well, and the antioxidant potential was determined. The defined profile of the biological activity of Lonicera caerulea leaves makes it possible to indicate this raw material as an essential material supporting the prevention and treatment of type II diabetes. However, this research showed that tested enzymes were strongly inhibited by the variety Jugana. The health-promoting potential of Lonicera caerulea leaves was correlated with the highest chlorogenic acid and rutin content in the variety Jugana.

Integrated volatile metabolomic and transcriptomic analysis provides insights into the regulation of floral scents between two contrasting varieties of Lonicera japonica

Lonicera japonica Thunb., belonging to the Caprifoliaceae family, is an important traditional Chinese medicinal plant. The L. japonica flower (LJF) is widely used in medicine, cosmetics, drinks, and food due to its medicinal and sweet-smelling properties. Considerable efforts have been devoted to investigating the pharmacological activities of LJF; however, the regulatory mechanism of the floral scents remains unknown. We previously selected and bred an elite variety of L. japonica var. chinensis Thunb. called 'Yujin2', which has a strong aroma and is used in functional drinks and cosmetics. In order to reveal the regulatory mechanism of the floral scents of LJF, volatile metabolomic and transcriptomic analyses of the LJF at the silver flowering stage of 'Yujin2' (strong aroma) and 'Fengjin1' (bland odor) were performed. Our results revealed that a total of 153 metabolites and 9,523 genes were differentially regulated in LJF between 'Yujin2' and 'Fengjin1'. The integrated analysis of omics data indicated that the biosynthetic pathways of terpenoids (i.e., monoterpenoids, including geraniol and alpha-terpineol; sesquiterpenoids, including farnesol, farnesal, and alpha-farnesene; triterpenoid squalene), tryptophan and its derivatives (methyl anthranilate), and fatty acid derivatives, were major contributors to the stronger aroma of 'Yujin2' compared to 'Fengjin1'. Moreover, several genes involved in the terpenoid biosynthetic pathway were characterized using quantitative real-time PCR. These results provide insights into the metabolic mechanisms and molecular basis of floral scents in LJF, enabling future screening of genes related to the floral scent regulation, such as alpha-terpineol synthase, geranylgeranyl diphosphate synthase, farnesyl pyrophosphate synthase, anthranilate synthase, as well as transcription factors such as MYB, WRKY, and LFY. The knowledge from this study will facilitate the breeding of quality-improved and more fragrant variety of L. japonica for ornamental purpose and functional beverages and cosmetics.

GAGE is a method for identification of plant species based on whole genome analysis and genome editing

Whole genomes of plants should be ideal databases for their species identification, but unfortunately there was no such method before this exploration. Here we report a plant species identification method based on the whole Genome Analysis and Genome Editing (GAGE). GAGE searches for target sequences from the whole genome of the subject plant and specifically detects them by employing a CRISPR/Cas12a system. Similar to how Mendel chose Pisum sativum (pea), we selected Crocus sativus (saffron) to establish GAGE, in which we constructed a library containing all candidate target sequences. Taking a target sequence in the ITS2 region as an example, we confirmed the feasibility, specificity, and sensitivity of GAGE. Consequently, we succeeded in not only using GAGE to identify Cr. sativus and its adulterants, but also executing GAGE in the plants from different classes including angiosperms, gymnosperms, ferns, and lycophytes. This sensitive and rapid method is the first plant species identification method based on the whole genome and provides new insights into the application of the whole genome in species identification.

A plant species identification method, GAGE, is reported that searches for target sequences from the whole genome of the subject plant and specifically detects them by employing a CRISPR/Cas12a system.

Transcriptional regulation mechanism of flavonoids biosynthesis gene during fruit development in astragalus membranaceus

Astragalus membranaceus, as an important medicinal plant, are an excellent source of flavonoids. Flavonoid compounds in A. membranaceus have been widely used in medicine and supplement, but known of the molecular mechanism of flavonoid biosynthesis is still very few. Here, we analyzed the association between flavonoid content and gene expression pattern during six different fruit developmental stages. Sixteen gene expression trends were significantly identified, involving 8,218 genes. The gene expression trend in profile 0 was positively correlated with flavonoid content, while the gene expression trend in profile 79 was negatively correlated with flavonoid content at six developmental stages. The expression level of genes involved in the general phenylpropane pathway was higher than that of genes involved in the flavonoid biosynthesis pathway. A total of 37 genes involved in flavonoid synthesis were identified in A. membranaceus. The expression pattern of flavonoid-related genes was highly correlated with flavonoid content. Our study deepened the understanding of the flavonoid synthesis mechanism and provided useful resources for future studies on the high flavonoid molecular breeding of A. membranaceus.

Flower color mutation, pink to orange, through CmGATA4 - CCD4a-5 module regulates carotenoids degradation in chrysanthemum

The carotenoids biosynthesis pathway in plants has been studied extensively, yet little is known about the regulatory mechanisms underlying this process, especially for ornamental horticulture plants. In this study, a natural variation of chrysanthemum with orange coloration was identified and compared with the wild type with pink coloration; the content and component of carotenoids were largely enriched in the mutant with orange coloration. CmCCD4a-5, the DNA sequence in both 'Pink yan' and the mutant, was identified and shown to function as a carotenoid degradation enzyme. Compared with 'Pink yan', the mutant shows lower expression level of CmCCD4a-5. Furthermore, CmGATA4 was found to have an opposite expression trend to CmCCD4a-5, and it could directly bind with the CmCCD4a-5 promoter. Taken together, this study demonstrates that CmGATA4 acts as a negative regulator of CmCCD4a-5 and, furthermore, low expression of CmCCD4a-5 resulted in carotenoid accumulation in the mutant.

Transcriptome profiles of yellowish-white and fuchsia colored flowers in the Rheum palmatum complex reveal genes related to color polymorphism

Flower color variation is ubiquitous in many plant species, and several studies have been conducted to elucidate the underlying molecular mechanism. There are two flower color variants (yellowish-white and fuchsia) in the Rheum palmatum complex, however, few studies have investigated this phenomenon. Here, we used transcriptome sequencing of the two color variants to shed light on the molecular and biochemical basis for these color morphs. Comparison of the two transcriptomes identified 9641 differentially expressed unigenes (DEGs), including 6477 up-regulated and 3163 down-regulated genes. Functional analyses indicated that several DEGs were related to the anthocyanin biosynthesis pathway, and the expression profiles of these DEGs were coincident with the qRT-PCR validation results, indicating that expression levels of structural genes have a profound effect on the color variation in the R. palmatum complex. Our results suggested that the interaction of transcription factors (MYB, bHLH and WRKY) also regulated the anthocyanin biosynthesis in the R. palmatum complex. Estimation of selection pressures using the dN/dS ratio showed that 1106 pairs of orthologous genes have undergone positive selection. Of these positively selected genes, 21 were involved in the anthocyanin biosynthetic pathway, indicating that they may encode the proteins for structural alteration and affect flower color in the R. palmatum complex.

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

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

Lilac (Syringa oblata) genome provides insights into its evolution and molecular mechanism of petal color change

Color change during flower opening is common; however, little is understood on the biochemical and molecular basis related. Lilac (Syringa oblata), a well-known woody ornamental plant with obvious petal color changes, is an ideal model. Here, we presented chromosome-scale genome assembly for lilac, resolved the flavonoids metabolism, and identified key genes and potential regulatory networks related to petal color change. The genome assembly is 1.05 Gb anchored onto 23 chromosomes, with a BUSCO score of 96.6%. Whole-genome duplication (WGD) event shared within Oleaceae was revealed. Metabolome quantification identified delphinidin-3-O-rutinoside (Dp3Ru) and cyanidin-3-O-rutinoside (Cy3Ru) as the major pigments; gene co-expression networks indicated WRKY an essential regulation factor at the early flowering stage, ERF more important in the color transition period (from violet to light nearly white), while the MBW complex participated in the entire process. Our results provide a foundation for functional study and molecular breeding in lilac.

New molecules in plant defence against pathogens

Plants host a multipart immune signalling network to ward off pathogens. Pathogen attack upon plant tissues can often lead to an amplified state of (induced) defence against subsequent infections in distal tissues; this is known as systemic acquired resistance (SAR). The interaction of plants with beneficial microbes of the rhizosphere microbiome can also lead to an induced resistance in above-ground plant tissues, known as induced systemic resistance. Second messengers such as calcium (Ca2+), reactive oxygen species (ROS), and nitric oxide (NO) are necessary for cell-to-cell signal propagation during SAR and show emergent roles in the mediation of other SAR metabolites. These include the lysine-derived signals pipecolic acid (Pip) and N-hydroxypipecolic acid (NHP), which are key signalling metabolites in SAR. Emerging evidence additionally pinpoints plant volatiles as modulators of defence signalling within and between plants. Plant volatile organic compounds (VOCs) such as monoterpenes can promote SAR by functioning through ROS. Furthermore, plant-derived and additionally also microbial VOCs can target both salicylic acid and jasmonic acid signalling pathways in plants and modulate defence against pathogens. In this review, an overview of recent findings in induced defence signalling, with a particular focus on newer signalling molecules and how they integrate into these networks is discussed.

Fingerprinting, structure, and genetic relationships among selected accessions of blue honeysuckle (Lonicera caerulea L.) from European collections

•

Scarcely on December 13, 2018, L. caerulea fruits were included in the list of novel foods in EU.

•

Hence, the growing interest in L. caerulea extends to its genome.

•

R-ISSR explores other than the RAPD, ISSR, AFLP, and RFLP techniques part of genome and analyze other range of genetic variability.

•

R-ISSR markers could be used for Lonicera core germplasm collection, development of SCARs, genetic map construction, barcoding, protection of variety rights, MAS, and genomic selection.

Due to its value and economic importance, the genome of Lonicera caerulea L. has been widely studied in various fields of science. In this study the genetic structure and relationships between 24 accessions of L. caerulea of different origins were assessed. A total of 692, 814, and 258 loci were amplified using 43 RAPD (random amplified polymorphic DNA), 40 ISSR (intersimple sequence repeat), and 20 R-ISSR (RAPD+ISSR) primers, respectively. Among the amplified loci, 66–78% were polymorphic and 12–20% were private. Selected R-ISSR sequences were detected in Lonicera japonica transcripts. Cluster and STRUCTURE analyses performed for each of the techniques revealed the existing differences and unknown similarities between the genotypes. The r-factor values calculated in the Mantel test indicated highly significant positive correlations between the Nei distance matrices, similar to the F_ST values (F_{ST_RAPD} = 0.223, F_{ST_ISSR} = 0.279, F_{ST_R-ISSR} = 0.363) determined in the analysis of molecular variance. It was found that 78%, 72%, and 64% of the genetic variations were related to the differences observed within the populations, which suggest that the variations are mainly reflected in the differences among the genotypes. The principal coordinate analysis showed greater differences between the mean distances of the Lonicera genotype pair and the actual distances of the same pairs on the Nei matrix compared to multidimensional scaling. These differences were 45%, 56%, and 42% higher for RAPD, ISSR and R-ISSR, respectively.

Chromosome-scale genome assembly of an important medicinal plant honeysuckle

Lonicera japonica (honeysuckle) is one of the most important medicinal plants and widely utilized in traditional Chinese medicine. At present, there are many varieties of honeysuckle used in cultivation, among which Sijihua variety are widely cultivated due to its wide adaptability, stress resistance, early flowering and high yield. In this study, we assembled the genome of Sijihua, which was approximately 886.04 Mb in size with a scaffold N50 of 79.5 Mb. 93.28% of the total assembled sequences were anchored to 9 pseudo-chromosomes by using PacBio long reads and Hi-C sequencing data. We predicted 39,320 protein-coding genes and 92.87% of them could be annotated in NR, GO, KOG, KEGG and other databases. In addition, we identified 644 tRNAs, 2,156 rRNAs, 109 miRNAs and 5,502 pseudogenes from the genome. The chromosome-scale genome of Sijihua will be a significant resource for understanding the genetic basis of high stress-resistance, which will facilitate further study of the genetic diversity and accelerate the genetic improvement and breeding of L. japonica.

Measurement(s)	Lonicera japonica • RNA sequencing • genome assembly • sequence annotation
Technology Type(s)	SMRT Sequencing • RNA sequencing • Hi-C • biomolecular annotation design
Factor Type(s)	Genotype
Sample Characteristic - Organism	Lonicera japonica
Sample Characteristic - Environment	occurrence
Sample Characteristic - Location	Shandong Province

Low-Input High-Molecular-Weight DNA Extraction for Long-Read Sequencing From Plants of Diverse Families

Long-read DNA sequencing technologies require high molecular weight (HMW) DNA of adequate purity and integrity, which can be difficult to isolate from plant material. Plant leaves usually contain high levels of carbohydrates and secondary metabolites that can impact DNA purity, affecting downstream applications. Several protocols and kits are available for HMW DNA extraction, but they usually require a high amount of input material and often lead to substantial DNA fragmentation, making sequencing suboptimal in terms of read length and data yield. We here describe a protocol for plant HMW DNA extraction from low input material (0.1 g) which is easy to follow and quick (2.5 h). This method successfully enabled us to extract HMW from four species from different families (Orchidaceae, Poaceae, Brassicaceae, Asteraceae). In the case of recalcitrant species, we show that an additional purification step is sufficient to deliver a clean DNA sample. We demonstrate the suitability of our protocol for long-read sequencing on the Oxford Nanopore Technologies PromethION® platform, with and without the use of a short fragment depletion kit.

Genomic Data Reveals Profound Genetic Structure and Multiple Glacial Refugia in Lonicera oblata (Caprifoliaceae), a Threatened Montane Shrub Endemic to North China

Characterizing genetic diversity and structure and identifying conservation units are both crucial for the conservation and management of threatened species. The development of high-throughput sequencing technology provides exciting opportunities for conservation genetics. Here, we employed the powerful SuperGBS method to identify 33, 758 high-quality single-nucleotide polymorphisms (SNP) from 134 individuals of a critically endangered montane shrub endemic to North China, Lonicera oblata. A low level of genetic diversity and a high degree of genetic differentiation among populations were observed based on the SNP data. Both principal component and phylogenetic analyses detected seven clusters, which correspond exactly to the seven geographic populations. Under the optimal K = 7, Admixture suggested the combination of the two small and geographically neighboring populations in the Taihang Mountains, Dongling Mountains, and Lijiazhuang, while the division of the big population of Jiankou Great Wall in the Yan Mountains into two clusters. High population genetic diversity and a large number of private alleles were detected in the four large populations, while low diversity and non-private alleles were observed for the remaining three small populations, implying the importance of these large populations as conservation units in priority. Demographic history inference suggested two drastic contractions of population size events that occurred after the Middle Pleistocene Transition and the Last Glacial Maximum, respectively. Combining our previous ecological niche modeling results with the present genomic data, there was a possible presence of glacial refugia in the Taihang and Yan Mountains, North China. This study provides valuable data for the conservation and management of L. oblata and broadens the understanding of the high biodiversity in the Taihang and Yan Mountains.

Genome-wide identification, molecular characterization, and gene expression analyses of honeysuckle NHX antiporters suggest their involvement in salt stress adaptation

Background

Ion homeostasis is an essential process for the survival of plants under salt stress. Na⁺/H⁺ antiporters (NHXs) are secondary ion transporters that regulate Na⁺ compartmentalization or efflux reduce Na⁺ toxicity and play a critical role during plant development and stress responses.

Methods and Results

To gain insight into the functional divergence of NHX genes in honeysuckle, a total of seven LjNHX genes were identified on the whole genome level and were renamed according to their chromosomal positions. All LjNHXs possessed the Na⁺/H⁺ exchanger domain and the amiloride-binding site was presented in all NHX proteins except LjNHX4. The phylogenetic analysis divided the seven NHX genes into Vac-clade (LjNHX1/2/3/4/5/7) and PM-clade (LjNHX6) based on their subcellular localization and validated by the distribution of conserved protein motifs and exon/intron organization analysis. The protein-protein interaction network showed that LjNHX4/5/6/7 shared the same putatively interactive proteins, including SOS2, SOS3, HKT1, and AVP1. Cis-acting elements and gene ontology (GO) analysis suggested that most LjNHXs involve in the response to salt stress through ion transmembrane transport. The expression profile analysis revealed that the expression levels of LjNHX3/7 were remarkably affected by salinity. These results suggested that LjNHXs play significant roles in honeysuckle development and response to salt stresses.

Conclusions

The theoretical foundation was established in the present study for the further functional characterization of the NHX gene family in honeysuckle.

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.

Comparative genomics reveal the convergent evolution of CYP82D and CYP706X members related to flavone biosynthesis in Lamiaceae and Asteraceae

Distant species producing the same secondary metabolites is an interesting and common phenomenon in nature. A classic example of this is scutellarein whose derivatives have been used clinically for more than 30 years. Scutellarein occurs in significant amounts in species of two different orders, Scutellaria baicalensis and Erigeron breviscapus, which diverged more than 100 million years ago. Here, according to the genome-wide selection and functional identification of 39 CYP450 genes from various angiosperms, we confirmed that only seven Scutellaria-specific CYP82D genes and one Erigeron CYP706X gene could perform the catalytic activity of flavone 6-hydroxylase (F6H), suggesting that the convergent evolution of scutellarein production in these two distant species was caused by two independently evolved CYP450 families. We also identified seven Scutellaria-specific CYP82D genes encoding flavone 8-hydroxylase (F8H). The evolutionary patterns of CYP82 and CYP706 families via kingdom-wide comparative genomics highlighted the evolutionary diversity of CYP82D and the specificity of CYP706X in angiosperms. Multi-collinearity and phylogenetic analysis of CYP82D in Scutellaria confirmed that the function of F6H evolved from F8H. Furthermore, the SbaiCYP82D1^A319D , EbreCYP706X^R130A , EbreCYP706X^F312D and EbreCYP706X^A318D mutants can significantly decrease the catalytic activity of F6H, revealing the contribution of crucial F6H amino acids to the scutellarein biosynthesis of distant species. This study provides important insights into the multi-origin evolution of the same secondary metabolite biosynthesis in the plant kingdom.

Comparative transcriptome analysis reveals variations of bioactive constituents in Lonicera japonica flowers under salt stress

Lonicera japonica flowers (LJF) is a traditional Chinese medicine packed with phenols constituents and widely used in the treatments of various diseases throughout the world. However, there is still very little known on how LJF identifies and resists salt stress. Here in, we systematically investigated the effect of salt on the phenotypic, metabolite, and transcriptomic in LJF. During long term stress (35 days), 1055 differential expression genes (DEGs) involved in the biosynthesis of secondary metabolites were screened through transcriptome analysis, among which the candidate genes and pathways involved in phenols biosynthesis were highlighted; and performed by phylogenetic tree analysis and multiple nucleotide sequence alignment. Ninety compounds were identified and their relative levels were compared between the control and stressed groups based on the LC-MS analysis, Putative biosynthesis networks of phenolic acid and flavonoid were con-structed with structural DEGs. Strikingly, the expression patterns of structural DEGs were mostly consistent with the variations of phenols under salt stress. Notably, the upregulation of UDP-glycosyl transferases under salt stress indicated post-modification of glycosyl transferases may participate in downstream flavonoids synthesis. This study reveals the relationships of the gene regulation and the phenols biosynthesis in LJF under salt stress, paving the way for the use of gene-specific expression to improve the yield of biocomponent.

A chromosome-scale genome assembly for the holly (Ilex polyneura) provides insights into genomic adaptations to elevation in Southwest China

Southwest China is a plant diversity hotspot. The near-cosmopolitan genus Ilex (c. 664 spp., Aquifoliaceae) reaches its maximum diversity in this region, with many narrow-range and a few widespread species. Divergent selection on widespread species leads to local adaptation, with consequences for both conservation and utilization, but is counteracted by geneflow. Many Ilex species are utilized as teas, medicines, ornamentals, honey plants, and timber, but variation below the species level is largely uninvestigated. We therefore studied the widespread Ilex polyneura, which occupies most of the elevational range available and is cultivated for its decorative leafless branches with persistent red fruits. We assembled a chromosome-scale genome using approximately 100x whole genome long-read and short-read sequencing combined with Hi-C sequencing. The genome is approximately 727.1 Mb, with a contig N50 size of 5 124 369 bp and a scaffold N50 size of 36 593 620 bp, for which the BUSCO score was 97.6%, and 98.9% of the assembly was anchored to 20 pseudochromosomes. Out of 32 838 genes predicted, 96.9% were assigned functions. Two whole genome duplication events were identified. Using this genome as a reference, we conducted a population genomics study of 112 individuals from 21 populations across the elevation range using restriction site-associated DNA sequencing (RADseq). Most populations clustered into four clades separated by distance and elevation. Selective sweep analyses identified 34 candidate genes potentially under selection at different elevations, with functions related to responses to abiotic and biotic stresses. This first high-quality genome in the Aquifoliales will facilitate the further domestication of the genus.

A Chromosome-Level Genome Assembly of the Pygmy Mole Cricket Xya riparia

The pygmy mole cricket Xya riparia (Orthoptera: Tridactyloidea) is rarely studied or widely known. Some species of pygmy mole crickets, however, not only have a potential ecological value but are also important in the study of the evolution of the orthopteran genome and its phylogenetic relationships. The genome resources of pygmy crickets are limited and there are currently no publications referencing this species' genome. In this study, we assembled a reference genome of X. riparia at the chromosomal level using nanopore sequencing and Hi-C technology. An X. riparia genome of 1.67 Gb was successfully assembled from 164.01 Gb of nanopore sequencing data. The genome assembly showed a completeness of 98.97% benchmarking universal single-copy orthologs with a contig N50 of 4.18 Mb and the longest contig being 18.84 Mb. The contigs were clustered, ordered, and correctly oriented on six pseuchromosomes, which covered 95.63% of the genome assembly through Hi-C data with a scaffold N50 of 319.1 Mb and the longest scaffold being 397.8 Mb. Repeat sequences accounted for 42.88% of the whole-genome assembly. A total of 60,847 noncoding RNAs were detected. Moreover, 16,468 (87.91%) of the genes were functionally annotated. As this is the first high-quality reference genome of X. riparia at the chromosomal level, it will undoubtedly serve as a valuable resource for ecological, biological, and genetic research on pygmy mole crickets as well as for general research on Orthoptera's genome evolution and phylogenetic relationships.

Transcriptomic analysis unravels the molecular response of Lonicera japonica leaves to chilling stress

Lonicera japonica is not only an important resource of traditional Chinese medicine, but also has very high horticultural value. Studies have been performed on the physiological responses of L. japonica leaves to chilling, however, the molecular mechanism underlying the low temperature-induced leaves morphological changes remains unclear. In this study, it has been demonstrated that the ratio of pigments content including anthocyanins, chlorophylls, and carotenoids was significantly altered in response to chilling condition, resulting in the color transformation of leaves from green to purple. Transcriptomic analysis showed there were 10,329 differentially expressed genes (DEGs) co-expressed during chilling stress. DEGs were mainly mapped to secondary metabolism, cell wall, and minor carbohydrate. The upregulated genes (UGs) were mainly enriched in protein metabolism, transport, and signaling, while UGs in secondary metabolism were mainly involved in phenylpropaoids-flavonoids pathway (PFP) and carotenoids pathway (CP). Protein-protein interaction analysis illustrated that 21 interacted genes including CAX3, NHX2, ACA8, and ACA9 were enriched in calcium transport/potassium ion transport. BR biosynthesis pathway related genes and BR insensitive (BRI) were collectively induced by chilling stress. Furthermore, the expression of genes involved in anthocyanins and CPs as well as the content of chlorogenic acid (CGA) and luteoloside were increased in leaves of L. japonica under stress. Taken together, these results indicate that the activation of PFP and CP in leaves of L. japonica under chilling stress, largely attributed to the elevation of calcium homeostasis and stimulation of BR signaling, which then regulated the PFP/CP related transcription factors.