Transcriptome Analysis of Secondary Metabolism Pathway, Transcription Factors, and Transporters in Response to Methyl Jasmonate in Lycoris aurea

Distinct fungal communities affecting opposite galanthamine accumulation patterns in two Lycoris species

Lycoris radiata is the main source of galanthamine, a clinical drug used in Alzheimer's disease; however, the galanthamine content in L. radiata is low. Lycoris aurea is another Lycoris species with high galanthamine content. Fungal endophytes can enhance plant secondary metabolite accumulation; thus, we compared the fungal communities in these two Lycoris species to identify certain fungal taxa in L. aurea capable of enhancing galanthamine accumulation. Several fungal endophytes, which were enriched in, exclusively isolated from L. aurea, or showed significant correlations with galanthamine, were demonstrated to enhance the accumulation of only galanthamine but no other Amaryllidaceae alkaloids (AAs) in L. radiata. These fungal endophytes mainly upregulated the downstream genes in the biosynthesis pathways of AAs in L. radiata, suggesting that they may allocate more precursors for galanthamine biosynthesis. This study demonstrated that fungal endophytes from L. aurea with higher galanthamine content can specifically enhance the accumulation of this medicinal alkaloid in other Lycoris species, thereby increasing the galanthamine source and reducing galanthamine separation and purification costs. This study broadens our understanding of the complex interactions between plant secondary metabolites and fungal endophytes.

Identification of transcription factor genes responsive to MeJA and characterization of a LaMYC2 transcription factor positively regulates lycorine biosynthesis in Lycoris aurea

Jasmonates (JAs) are among the main phytohormones, regulating plant growth and development, stress responses, and secondary metabolism. As the major regulator of the JA signaling pathway, MYC2 also plays an important role in plant secondary metabolite synthesis and accumulation. In this study, we performed a comparative transcriptome analysis of Lycoris aurea seedlings subjected to methyl jasmonate (MeJA) at different treatment times. A total of 31,193 differentially expressed genes (DEGs) were identified by RNA sequencing. Among them, 732 differentially expressed transcription factors (TFs) comprising 51 TF families were characterized. The most abundant TF family was WRKY proteins (80), followed by AP2/ERF-EFR (67), MYB (59), bHLH (52), and NAC protein (49) families. Subsequently, by calculating the Pearson's correlation coefficient (PCC) between the expression level of TF DEGs and the lycorine contents, 41 potential TF genes (|PCC| >0.8) involved in lycorine accumulation were identified, including 36 positive regulators and 5 negative regulators. Moreover, a MeJA-inducible MYC2 gene (namely LaMYC2) was cloned on the basis of transcriptome sequencing. Bioinformatic analyses revealed that LaMYC2 proteins contain the bHLH-MYC_N domain and bHLH-AtAIB_like motif. LaMYC2 protein is localized in the cell nucleus, and can partly rescue the MYC2 mutant in Arabidopsis thaliana. LaMYC2 protein could interact with most LaJAZs (especially LaJAZ3 and LaJAZ4) identified previously. Transient overexpression of LaMYC2 increased lycorine contents in L. aurea petals, which might be associated with the activation of the transcript levels of tyrosine decarboxylase (TYDC) and phenylalanine ammonia lyase (PAL) genes. By isolating the 887-bp-length promoter fragment upstream of the start codon (ATG) of LaTYDC, we found several different types of E-box motifs (CANNTG) in the promoter of LaTYDC. Further study demonstrated that LaMYC2 was indeed able to bind the E-box (CACATG) present in the LaTYDC promoter, verifying that the pathway genes involved in lycorine biosynthesis could be regulated by LaMYC2, and that LaMYC2 has positive roles in the regulation of lycorine biosynthesis. These findings demonstrate that LaMYC2 is a positive regulator of lycorine biosynthesis and may facilitate further functional research of the LaMYC2 gene, especially its potential regulatory roles in Amaryllidaceae alkaloid accumulation in L. aurea.

Integrated analysis of transcriptome and miRNAome reveals the heat stress response of Pinellia ternata seedlings

Pinellia ternata (Thunb.) Briet., a valuable herb native to China, is susceptible to the "sprout tumble" phenomenon because of high temperatures, resulting in a significant yield reduction. However, the molecular regulatory mechanisms underlying the response of P. ternata to heat stress are not well understood. In this study, we integrated transcriptome and miRNAome sequencing to identify heat-response genes, microRNAs (miRNAs), and key miRNA-target pairs in P. ternata that differed between heat-stress and room-temperature conditions. Transcriptome analysis revealed extensive reprogramming of 4,960 genes across various categories, predominantly associated with cellular and metabolic processes, responses to stimuli, biological regulation, cell parts, organelles, membranes, and catalytic and binding activities. miRNAome sequencing identified 1,597 known/conserved miRNAs that were differentially expressed between the two test conditions. According to the analysis, genes and miRNAs associated with the regulation of transcription, DNA template, transcription factor activity, and sequence-specific DNA binding pathways may play a major role in the resistance to heat stress in P. ternata. Integrated analysis of the transcriptome and miRNAome expression data revealed 41 high-confidence miRNA-mRNA pairs, forming 25 modules. MYB-like proteins and calcium-responsive transcription coactivators may play an integral role in heat-stress resistance in P. ternata. Additionally, the candidate genes and miRNAs were subjected to quantitative real-time polymerase chain reaction to validate their expression patterns. These results offer a foundation for future studies exploring the mechanisms and critical genes involved in heat-stress resistance in P. ternata.

Characterization of the Heat Shock Transcription Factor Family in Lycoris radiata and Its Potential Roles in Response to Abiotic Stresses

Heat shock transcription factors (HSFs) are an essential plant-specific transcription factor family that regulates the developmental and growth stages of plants, their signal transduction, and their response to different abiotic and biotic stresses. The HSF gene family has been characterized and systematically observed in various species; however, research on its association with Lycoris radiata is limited. This study identified 22 HSF genes (LrHSFs) in the transcriptome-sequencing data of L. radiata and categorized them into three classes including HSFA, HSFB, and HSFC, comprising 10, 8, and 4 genes, respectively. This research comprises basic bioinformatics analyses, such as protein sequence length, molecular weight, and the identification of its conserved motifs. According to the subcellular localization assessment, most LrHSFs were present in the nucleus. Furthermore, the LrHSF gene expression in various tissues, flower developmental stages, two hormones stress, and under four different abiotic stresses were characterized. The data indicated that LrHSF genes, especially LrHSF5, were essentially involved in L. radiata development and its response to different abiotic and hormone stresses. The gene-gene interaction network analysis revealed the presence of synergistic effects between various LrHSF genes' responses against abiotic stresses. In conclusion, these results provided crucial data for further functional analyses of LrHSF genes, which could help successful molecular breeding in L. radiata.

Functional Characterization of AP2/ERF Transcription Factors during Flower Development and Anthocyanin Biosynthesis Related Candidate Genes in Lycoris

The APETALA2/ethylene-responsive transcription factor (AP2/ERF) family has been extensively investigated because of its significant involvement in plant development, growth, fruit ripening, metabolism, and plant stress responses. To date, there has been little investigation into how the AP2/ERF genes influence flower formation and anthocyanin biosynthesis in Lycoris. Herein, 80 putative LrAP2/ERF transcription factors (TFs) with complete open reading frames (ORFs) were retrieved from the Lycoris transcriptome sequence data, which could be divided into five subfamilies dependent on their complete protein sequences. Furthermore, our findings demonstrated that genes belonging to the same subfamily had structural similarities and conserved motifs. LrAP2/ERF genes were analyzed for playing an important role in plant growth, water deprivation, and flower formation by means of gene ontology (GO) enrichment analysis. The expression pattern of the LrAP2/ERF genes differed across tissues and might be important for Lycoris growth and flower development. In response to methyl jasmonate (MeJA) exposure and drought stress, the expression of each LrAP2/ERF gene varied across tissues and time. Moreover, a total of 20 anthocyanin components were characterized using ultra-performance liquid chromatography-electrospray ionization tandem mass spectrometry (UPLC-ESI-MS/MS) analysis, and pelargonidin-3-O-glucoside-5-O-arabinoside was identified as the major anthocyanin aglycone responsible for the coloration of the red petals in Lycoris. In addition, we mapped the relationships between genes and metabolites and found that LrAP2/ERF16 is strongly linked to pelargonidin accumulation in Lycoris petals. These findings provide the basic conceptual groundwork for future research into the molecular underpinnings and regulation mechanisms of AP2/ERF TFs in anthocyanin accumulation and Lycoris floral development.

Methyl jasmonate and salicylic acid as powerful elicitors for enhancing the production of secondary metabolites in medicinal plants: an updated review

Plant secondary metabolites are bioactive scaffolds that are crucial for plant survival in the environment and to maintain a defense mechanism from predators. These compounds are generally present in plants at a minimal level and interestingly, they are found to have a wide variety of therapeutic values for humans. Several medicinal plants are used for pharmaceutical purposes due to their affordability, fewer adverse effects, and vital role in traditional remedies. Owing to this reason, these plants are exploited at a high range worldwide and therefore many medicinal plants are on the threatened list. There is a need of the hour to tackle this major problem, one effective approach called elicitation can be used to enhance the level of existing and novel plant bioactive compounds using different types of elicitors namely biotic and abiotic. This process can be generally achieved by in vitro and in vivo experiments. The current comprehensive review provides an overview of biotic and abiotic elicitation strategies used in medicinal plants, as well as their effects on secondary metabolites enhancement. Further, this review mainly deals with the enhancement of biomass and biosynthesis of different bioactive compounds by methyl jasmonate (MeJA) and salicylic acid (SA) as elicitors of wide medicinal plants in in vitro by using different cultures. The present review was suggested as a significant groundwork for peers working with medicinal plants by applying elicitation strategies along with advanced biotechnological approaches.

System network analysis of Rosmarinus officinalis transcriptome and metabolome-Key genes in biosynthesis of secondary metabolites

Medicinal plants contain valuable compounds that have attracted worldwide interest for their use in the production of natural drugs. The presence of compounds such as rosmarinic acid, carnosic acid, and carnosol in Rosmarinus officinalis has made it a plant with unique therapeutic effects. The identification and regulation of the biosynthetic pathways and genes will enable the large-scale production of these compounds. Hence, we studied the correlation between the genes involved in biosynthesis of the secondary metabolites in R. officinalis using proteomics and metabolomics data by WGCNA. We identified three modules as having the highest potential for the metabolite engineering. Moreover, the hub genes highly connected to particular modules, TFs, PKs, and transporters were identified. The TFs of MYB, C3H, HB, and C2H2 were the most likely candidates associated with the target metabolic pathways. The results indicated that the hub genes including Copalyl diphosphate synthase (CDS), Phenylalanine ammonia lyase (PAL), Cineole synthase (CIN), Rosmarinic acid synthase (RAS), Tyrosine aminotransferase (TAT), Cinnamate 4-hydroxylase (C4H), and MYB58 are responsible for biosynthesis of important secondary metabolites. Thus, we confirmed these results using qRT-PCR after treating R. officinalis seedlings with methyl jasmonate. These candidate genes may be employed for genetic and metabolic engineering research to increase R. officinalis metabolite production.

Characterization of the WRKY Gene Family Related to Anthocyanin Biosynthesis and the Regulation Mechanism under Drought Stress and Methyl Jasmonate Treatment in Lycoris radiata

Lycoris radiata, belonging to the Amaryllidaceae family, is a well-known Chinese traditional medicinal plant and susceptible to many stresses. WRKY proteins are one of the largest families of transcription factors (TFs) in plants and play significant functions in regulating physiological metabolisms and abiotic stress responses. The WRKY TF family has been identified and investigated in many medicinal plants, but its members and functions are not identified in L. radiata. In this study, a total of 31 L. radiata WRKY (LrWRKY) genes were identified based on the transcriptome-sequencing data. Next, the LrWRKYs were divided into three major clades (Group I-III) based on the WRKY domains. A motif analysis showed the members within same group shared a similar motif component, indicating a conservational function. Furthermore, subcellular localization analysis exhibited that most LrWRKYs were localized in the nucleus. The expression pattern of the LrWRKY genes differed across tissues and might be important for Lycoris growth and flower development. There were large differences among the LrWRKYs based on the transcriptional levels under drought stress and MeJA treatments. Moreover, a total of 18 anthocyanin components were characterized using an ultra-performance liquid chromatography-electrospray ionization tandem mass spectrometry (UPLC-ESI-MS/MS) analysis and pelargonidin-3-O-glucoside-5-O-arabinoside as well as cyanidin-3-O-sambubioside were identified as the major anthocyanin aglycones responsible for the coloration of the red petals in L. radiata. We further established a gene-to-metabolite correlation network and identified LrWRKY3 and LrWRKY27 significant association with the accumulation of pelargonidin-3-O-glucoside-5-O-arabinoside in the Lycoris red petals. These results provide an important theoretical basis for further exploring the molecular basis and regulatory mechanism of WRKY TFs in anthocyanin biosynthesis and in response to drought stress and MeJA treatment.

Transcriptome-wide characterization of bHLH transcription factor genes in Lycoris radiata and functional analysis of their response to MeJA

As one of the biggest plant specific transcription factor (TF) families, basic helix-loop-helix (bHLH) protein, plays significant roles in plant growth, development, and abiotic stress responses. However, there has been minimal research about the effects of methyl jasmonate (MeJA) treatment on the bHLH gene family in Lycoris radiata (L'Her.) Herb. In this study, based on transcriptome sequencing data, 50 putative L. radiata bHLH (LrbHLH) genes with complete open reading frames (ORFs), which were divided into 20 bHLH subfamilies, were identified. The protein motif analyses showed that a total of 10 conserved motifs were found in LrbHLH proteins and motif 1 and motif 2 were the most highly conserved motifs. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of LrbHLH genes revealed their involvement in regulation of plant growth, jasmonic acid (JA) mediated signaling pathway, photoperiodism, and flowering. Furthermore, subcellular localization revealed that most LrbHLHs were located in the nucleus. Expression pattern analysis of LrbHLH genes in different tissues and at flower developmental stages suggested that their expression differed across lineages and might be important for plant growth and organ development in Lycoris. In addition, all LrbHLH genes exhibited specific spatial and temporal expression patterns under MeJA treatment. Moreover, protein-protein interaction (PPI) network analysis and yeast two-hybrid assay showed that numerous LrbHLHs could interact with jasmonate ZIM (zinc-finger inflorescence meristem) domain (JAZ) proteins. This research provides a theoretical basis for further investigation of LrbHLHs to find their functions and insights for their regulatory mechanisms involved in JA signaling pathway.

Interplay of transcription factors orchestrating the biosynthesis of plant alkaloids

Plants produce a range of secondary metabolites primarily as defence molecules. A plant has to invest considerable energy to synthesise alkaloids, and sometimes they are even toxic to themselves. Hence, the biosynthesis of alkaloids is a spatiotemporally regulated process under quantitative feedback regulation which is accomplished by the signal reception, transcriptional/translational regulation, transport, storage and accumulation. The transcription factors (TFs) initiate the biosynthesis of alkaloids after appropriate cues. The present study recapitulates last decade understanding of the role of TFs in alkaloid biosynthesis. The present review discusses TF families, viz. AP2/ERF, bHLH, WRKY, MYB involved in the biosynthesis of various types of alkaloids. It also highlights the role of the jasmonic acid cascade and post-translational modifications of TF proteins. A thorough understanding of TFs will help us to decide a strategy to facilitate successful pathway manipulation and in vitro production.

Comparative transcriptome and metabolome analyses identified the mode of sucrose degradation as a metabolic marker for early vegetative propagation in bulbs of Lycoris

Vegetative propagation (VP) is an important practice for production in many horticultural plants. Sugar supply constitutes the basis of VP in bulb flowers, but the underlying molecular basis remains elusive. By performing a combined sequencing technologies coupled with ultra-high performance liquid chromatography-quadrupole time-of-flight mass spectrometry approach for metabolic analyses, we compared two Lycoris species with contrasting regeneration rates: high-regeneration Lycoris sprengeri and low-regeneration Lycoris aurea. A comprehensive multi-omics analyses identified both expected processes involving carbohydrate metabolism and transcription factor networks, as well as the metabolic characteristics for each developmental stage. A higher abundance of the differentially expressed genes including those encoding ethylene responsive factors was detected at bulblet initiation stage compared to the late stage of bulblet development. High hexose-to-sucrose ratio correlated to bulblet formation across all the species examined, indicating its role in the VP process in Lycoris bulb. Importantly, a clear difference between cell wall invertase (CWIN)-catalyzed sucrose unloading in high-regeneration species and the sucrose synthase-catalyzed pathway in low-regeneration species was observed at the bulblet initiation stage, which was supported by findings from carboxyfluorescein tracing and quantitative real-time PCR analyses. Collectively, the findings indicate a sugar-mediated model of the regulation of VP in which high CWIN expression or activity may promote bulblet initiation via enhancing apoplasmic unloading of sucrose or sugar signals, whereas the subsequent high ratio of hexose-to-sucrose likely supports cell division characterized in the next phase of bulblet formation.

Aneuploidy promotes intraspecific diversification of the endemic East Asian herb Lycoris aurea complex

Polyploidy has received considerable interest in the past, but aneuploidy and partial rearrangements may also influence genomic divergence. In this study, we reported a comprehensive cytogeographic, morphological and genetic analysis of Lycoris aurea complex throughout its range and attempted to explore the association between aneuploidy and species diversification. The karyotypes of this complex presented aneuploidy variations mainly divided into four cytotypes: I (2n = 10m + 2T), II (2n = 8m + 6T), III (2n = 7m + 8T), and IV (2n = 6m + 10T). Cytotype distributions were highly structured geographically. Two main cytotypes, II and IV, are geographically allopatric. The populations with cytotype II are mainly distributed in central China and the southern islands of Japan. Cytotypes IV is disjunctly distributed in southwestern and southeastern China. The cytotypes with fewer chromosome numbers tend to occur at high latitudes. For analyzing the phylogeographic pattern and genetic structure of this complex, we sequenced four chloroplast DNA fragments (4,748 bp in total) of 241 individuals from 42 populations. Extremely high diversity of cpDNA haplotypes was found, with genetic diversity index (H _d) being 0.932 and 98.61% of the genetic variation occurring among populations, indicating that this complex has undergone strong intraspecific differentiation. The cytotype II had the highest haplotype diversity (H _d = 0.885), while cytotype IV harbored the highest nucleotide diversity (π = 4.09 × 10^-3). We detected significant leaf morphological differences not only between cytotype II and IV but also between west lineage and east lineage within cytotype IV. These results illustrated that aneuploidy contributed to extensive morphological and genetic differentiation in L. aurea complex. It was suggested that L. aurea complex should comprise multiple independent evolutionary lineages, and accurate species delimitation needs to be established further in an integrative taxonomic approach.

Biotechnological Approaches to Optimize the Production of Amaryllidaceae Alkaloids

Amaryllidaceae alkaloids (AAs) are plant specialized metabolites with therapeutic properties exclusively produced by the Amaryllidaceae plant family. The two most studied representatives of the family are galanthamine, an acetylcholinesterase inhibitor used as a treatment of Alzheimer’s disease, and lycorine, displaying potent in vitro and in vivo cytotoxic and antiviral properties. Unfortunately, the variable level of AAs’ production in planta restricts most of the pharmaceutical applications. Several biotechnological alternatives, such as in vitro culture or synthetic biology, are being developed to enhance the production and fulfil the increasing demand for these AAs plant-derived drugs. In this review, current biotechnological approaches to produce different types of bioactive AAs are discussed.

BrCNGC gene family in field mustard: genome-wide identification, characterization, comparative synteny, evolution and expression profiling

CNGCs are ligand-gated calcium signaling channels, which participate in important biological processes in eukaryotes. However, the CNGC gene family is not well-investigated in Brassica rapa L. (i.e., field mustard) that is economically important and evolutionary model crop. In this study, we systematically identified 29 member genes in BrCNGC gene family, and studied their physico-chemical properties. The BrCNGC family was classified into four major and two sub phylogenetic groups. These genes were randomly localized on nine chromosomes, and dispersed into three sub-genomes of B. rapa L. Both whole-genome triplication and gene duplication (i.e., segmental/tandem) events participated in the expansion of the BrCNGC family. Using in-silico bioinformatics approaches, we determined the gene structures, conserved motif compositions, protein interaction networks, and revealed that most BrCNGCs can be regulated by phosphorylation and microRNAs of diverse functionality. The differential expression patterns of BrCNGC genes in different plant tissues, and in response to different biotic, abiotic and hormonal stress types, suggest their strong role in plant growth, development and stress tolerance. Notably, BrCNGC-9, 27, 18 and 11 exhibited highest responses in terms of fold-changes against club-root pathogen Plasmodiophora brassicae, Pseudomonas syringae pv. maculicola, methyl-jasmonate, and trace elements. These results provide foundation for the selection of candidate BrCNGC genes for future breeding of field mustard.

Early Sucrose Degradation and the Dominant Sucrose Cleavage Pattern Influence Lycoris sprengeri Bulblet Regeneration In Vitro

Bulblet formation and development determine the quantitative and qualitative traits, respectively, of bulb yield for most flowering bulbs. For Lycoris species, however, the underlying molecular mechanism remains elusive. Here, clonal bulblets of Lycoris sprengeri (Ls) derived from the same probulb were used as explants to establish efficient and inefficient in vitro regeneration systems by adjusting the 6-benzyladenine (BA) concentrations in media. BA application did not change the biological processes among groups but led to earlier decreases in sucrose and total soluble sugar (TSS) contents. Correlation analyses showed that the BA treatments changed the interaction between carbohydrate and endogenous hormone contents during bulblet regeneration. We found that two sucrose degradation enzyme-related genes, cell wall invertase (CWIN) and sucrose synthase, exhibited exactly opposite expression patterns during the competence stage. In addition, the regeneration system that obtained more bulblets showed significantly higher expression of LsCWIN2 than those that obtained fewer bulblets. Our data demonstrate the essential role of BA in accelerating sucrose degradation and the selection of a dominant sucrose cleavage pattern at the competence stage of in vitro bulblet regeneration. We propose that a relatively active CWIN-catalyzed pathway at the competence stage might promote bulblet regeneration, thus influencing bulb yield.

An ATP-Binding Cassette Transporter, LaABCB11, Contributes to Alkaloid Transport in Lycoris aurea

As a kind of Amaryllidaceae alkaloid which is accumulated in the species of Lycoris plants, lycorine has a range of physiological effects. The biosynthesis pathway of lycorine has been partly revealed, but the transport and accumulation mechanisms of lycorine have rarely been studied. In this study, an ATP-binding cassette (ABC) transporter from Lycoris aurea (L'Hér) Herb., namely LaABCB11, was cloned and functionally characterized. Heterologous expression showed that LaABCB11 transported lycorine in an outward direction, increased the tolerance of yeast cells to lycorine, and caused a lower lycorine accumulation in transformants than control or mutant in yeast. LaABCB11 is associated with the plasma membrane, and in situ hybridization indicated that LaABCB11 was mainly expressed in the phloem of leaves and bulbs, as well as in the cortical cells of roots. These findings suggest that LaABCB11 functions as a lycorine transport and it might be related to the translocation and accumulation of lycorine from the leaves and bulbs to the roots.

Systemic Expression of Genes Involved in the Plant Defense Response Induced by Wounding in Senna tora

Wounds in tissues provide a pathway of entry for pathogenic fungi and bacteria in plants. Plants respond to wounding by regulating the expression of genes involved in their defense mechanisms. To analyze this response, we investigated the defense-related genes induced by wounding in the leaves of Senna tora using RNA sequencing. The genes involved in jasmonate and ethylene biosynthesis were strongly induced by wounding, as were a large number of genes encoding transcription factors such as ERFs, WRKYs, MYBs, bHLHs, and NACs. Wounding induced the expression of genes encoding pathogenesis-related (PR) proteins, such as PR-1, chitinase, thaumatin-like protein, cysteine proteinase inhibitor, PR-10, and plant defensin. Furthermore, wounding led to the induction of genes involved in flavonoid biosynthesis and the accumulation of kaempferol and quercetin in S. tora leaves. All these genes were expressed systemically in leaves distant from the wound site. These results demonstrate that mechanical wounding can lead to a systemic defense response in the Caesalpinioideae, a subfamily of the Leguminosae. In addition, a co-expression analysis of genes induced by wounding provides important information about the interactions between genes involved in plant defense responses.

Investigation of Chemical Profiles of Different Parts of Morus alba Using a Combination of Molecular Networking Methods with Mass Spectral Data from Two Ionization Modes of LC/MS

Plants produce numerous secondary metabolites with diverse physicochemical properties. Because different parts of a single plant produce various components, several spectroscopic methods are necessary to inspect their chemical profiles. Mass spectral data are recognized as one of the most useful tools for analyzing components with a wide range of polarities. However, interpreting mass spectral data generated from positive and negative ionization modes is a challenging task because of the diverse chemical profiles of secondary metabolites. Herein, we combine and analyze mass spectral data generated in two ionization modes to detect as many metabolites as possible using the molecular networking approach. We selected different parts of a single plant, Morus alba (Moraceae), which are used in the functional food and medicinal herb industries. The mass spectral data generated from two ionization modes were combined and analyzed using various molecular networking workflows. We confirmed that our approach could be applied to simultaneously analyze the different types of secondary metabolites with different physicochemical properties.

Transcriptome and Small RNA Combined Sequencing Analysis of Cold Tolerance in Non-heading Chinese Cabbage

Background

Non-heading Chinese cabbage (Brassica rapa ssp. chinensis) is an important leaf vegetable grown worldwide. However, there has currently been not enough transcriptome and small RNA combined sequencing analysis of cold tolerance, which hinders further functional genomics research.

Results

In this study, 63.43 Gb of clean data was obtained from the transcriptome analysis. The clean data of each sample reached 6.99 Gb, and the basic percentage of Q30 was 93.68% and above. The clean reads of each sample were sequence aligned with the designated reference genome (Brassica rapa, IVFCAASv1), and the efficiency of the alignment varied from 81.54 to 87.24%. According to the comparison results, 1,860 new genes were discovered in Pak-choi, of which 1,613 were functionally annotated. Among them, 13 common differentially expressed genes were detected in all materials, including seven upregulated and six downregulated. At the same time, we used quantitative real-time PCR to confirm the changes of these gene expression levels. In addition, we sequenced miRNA of the same material. Our findings revealed a total of 34,182,333 small RNA reads, 88,604,604 kinds of small RNAs, among which the most common size was 24 nt. In all materials, the number of common differential miRNAs is eight. According to the corresponding relationship between miRNA and its target genes, we carried out Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis on the set of target genes on each group of differentially expressed miRNAs. Through the analysis, it is found that the distributions of candidate target genes in different materials are different. We not only used transcriptome sequencing and small RNA sequencing but also used experiments to prove the expression levels of differentially expressed genes that were obtained by sequencing. Sequencing combined with experiments proved the mechanism of some differential gene expression levels after low-temperature treatment.

Conclusion

In all, this study provides a resource for genetic and genomic research under abiotic stress in Pak-choi.

De Novo Transcriptome Assembly of Two Microsorum Fern Species Identifies Enzymes Required for Two Upstream Pathways of Phytoecdysteroids

Microsorum species produce a high amount of phytoecdysteroids (PEs), which are widely used in traditional medicine in the Pacific islands. The PEs in two different Microsorum species, M. punctatum (MP) and M. scolopendria (MS), were examined using high-performance liquid chromatography (HPLC). In particular, MS produces a high amount of 20-hydroxyecdysone, which is the main active compound in PEs. To identify genes for PE biosynthesis, we generated reference transcriptomes from sterile frond tissues using the NovaSeq 6000 system. De novo transcriptome assembly after deleting contaminants resulted in 57,252 and 54,618 clean transcripts for MP and MS, respectively. The clean Microsorum transcripts for each species were annotated according to gene ontology terms, UniProt pathways, and the clusters of the orthologous group protein database using the MEGAN6 and Sma3s programs. In total, 1852 and 1980 transcription factors were identified for MP and MS, respectively. We obtained transcripts encoding for 38 and 32 enzymes for MP and MS, respectively, potentially involved in mevalonate and sterol biosynthetic pathways, which produce precursors for PE biosynthesis. Phylogenetic analyses revealed many redundant and unique enzymes between the two species. Overall, this study provides two Microsorum reference transcriptomes that might be useful for further studies regarding PE biosynthesis in Microsorum species.

Sweet Basil Has Distinct Synthases for Eugenol Biosynthesis in Glandular Trichomes and Roots with Different Regulatory Mechanisms

Production of a volatile phenylpropene; eugenol in sweet basil is mostly associated with peltate glandular trichomes (PGTs) found aerially. Currently only one eugenol synthase (EGS), ObEGS1 which belongs to PIP family is identified from sweet basil PGTs. Reports of the presence of eugenol in roots led us to analyse other EGSs in roots. We screened for all the PIP family reductase transcripts from the RNA-Seq data. In vivo functional characterization of all the genes in E. coli showed their ability to produce eugenol and were termed as ObEGS2-8. Among all, ObEGS1 displayed highest expression in PGTs and ObEGS4 in roots. Further, eugenol was produced only in the roots of soil-grown plants, but not in roots of aseptically-grown plants. Interestingly, eugenol production could be induced in roots of aseptically-grown plants under elicitation suggesting that eugenol production might occur as a result of environmental cues in roots. The presence of ObEGS4 transcript and protein in aseptically-grown plants indicated towards post-translational modifications (PTMs) of ObEGS4. Bioinformatics analysis showed possibility of phosphorylation in ObEGS4 which was further confirmed by in vitro experiment. Our study reveals the presence of multiple eugenol synthases in sweet basil and provides new insights into their diversity and tissue specific regulation.

Seasonal Variation in Transcriptomic Profiling of Tetrastigma hemsleyanum Fully Developed Tuberous Roots Enriches Candidate Genes in Essential Metabolic Pathways and Phytohormone Signaling

Tetrastigma hemsleyanum Diels et Gilg (Sanyeqing, SYQ) is a perennial climbing liana and an endemic plant to southern China. Its tuberous roots (TRs) are used in traditional Chinese medicine for treating some diseases such as high fever, pneumonia, asthma, hepatitis, and cancers. However, the mechanisms underlying the development of TR and the content of flavonoids and phenylpropanoids (FPs) are not well-understood. In this study, we performed a transcriptomic analysis of 12 fully developed TR (FD-TR) samples harvested in four seasons [spring (Sp), summer (Su), autumn (Au), and winter (Wi)] using the RNA-Sequencing (RNA-Seq). We obtained a total of 78.54 Gb raw data and 65,578 unigenes. Then, the unigenes were annotated by using six databases such as non-redundant protein database (NR), Pfam, eggNOG, SWISSProt, Kyoto Encyclopedia of Genes and Genomes (KEGG), and gene ontology (GO). The transcriptomic profiling showed closer relationships between the samples obtained in Su and Au than those obtained in Sp and Wi based on the results of both total unigenes and differentially expressed genes (DEGs). Three pathways, including the biosynthesis of FPs, metabolism of starch and sucrose, and signaling of phytohormones, were highly enriched, suggesting a gene-level seasonal variation. Based on the numbers of DEGs, brassinosteroid (BR) signal transduction factors appeared to play a key role in modulating the development of TRs while most of the auxin signaling genes were mainly activated in Wi and Sp FD-TRs. Most genes in the biosynthesis and biodegradation of starch and biodegradation of cellulose were activated in Wi FD-TRs. As determined by the high performance liquid chromatography (HPLC) and aluminum nitrate colorimetric method, the contents of total flavonoids and most detected FP components increased from Sp to Au but decreased in Wi. Enhanced expression levels of some genes in the biosynthetic pathways of FPs were detected in Su and Au samples, which corroborated well with metabolite content. Our findings provide the first transcriptomic and biochemical data on a seasonal variation in the composition of medically important metabolites in SYQ FD-TRs.

Transcriptomic and Metabolomic Analyses Reveals That Exogenous Methyl Jasmonate Regulates Galanthamine Biosynthesis in Lycoris longituba Seedlings

The Amaryllidaceae alkaloid galanthamine (Gal) in Lycoris longituba is a secondary metabolite that has been used to treat Alzheimer's disease. Plant secondary metabolism is affected by methyl jasmonate (MeJA) exposure, although the regulatory mechanisms of MeJA on L. longituba seedlings remains largely unknown. In the present study, 75, 150, and 300 μM MeJA were used as treatments on L. longituba seedlings for 7, 14, 21, and 28 days, while 0 μM MeJA was used as the control (MJ-0). The effect of exogenous MeJA on Gal synthesis in L. longituba was then investigated using transcriptomic sequencing and metabolite profiling via GC-MS and LC-MS analysis. Galanthamine (Gal), lycorine (Lyc), and lycoramine (Lycm) abundances were 2. 71-, 2. 01-, and 2.85-fold higher in 75 μM MeJA (MJ-75) treatment plants compared to MJ-0 treatment plants after 7 days of cultivation. Transcriptomic analysis further showed that MJ-75 treatment significantly induced the expression of norbelladine synthase (NBS) and norbelladine 4'-O-methyltransferase (OMT), which are involved in the Gal biosynthesis pathway. In addition, increased expression was observed in MJ-75 treatment plants for genes in the JA synthesis and JA signaling pathways including those of allene oxide cyclase (AOC), 12-oxo-phytodienoic acid reductase (OPR), jasmonic acid amino acid synthase (JAR), and transcription factor MYC. The L. longituba tyrosine decarboxylase (LlTYDC) enzyme was identified and proposed to be involved in the Gal biosynthetic pathway. Metabolomics results demonstrated that the accumulation of Amaryllidaceae alkaloids, and especially alkaloids in the Gal biosynthesis pathway, could be induced by MJ-75 treatment. Interestingly, metabolites in the JA synthesis pathway were also affected by MeJA treatment. Overall, this multi-omics study suggests that both the JA synthesis/JA signaling and Gal biosynthesis pathways were affected by exogenous MeJA treatment. This comprehensive study of gene expression and metabolite contents can help us better understand the molecular mechanisms underlying MeJA-mediated Gal biosynthesis in L. longituba.

Characterization, validation, and cross-species transferability of EST-SSR markers developed from Lycoris aurea and their application in genetic evaluation of Lycoris species

BACKGROUND

The Lycoris genus includes many ornamentally and medicinally important species. Polyploidization and hybridization are considered modes of speciation in this genus, implying great genetic diversity. However, the lack of effective molecular markers has limited the genetic analysis of this genus.

RESULTS

In this study, mining of EST-SSR markers was performed using transcriptome sequences of L. aurea, and 839 primer pairs for non-redundant EST-SSRs were successfully designed. A subset of 60 pairs was randomly selected for validation, of which 44 pairs could amplify products of the expected size. Cross-species transferability of the 60 primer pairs among Lycoris species were assessed in L. radiata Hreb, L. sprengeri Comes ex Baker, L. chinensis Traub and L. anhuiensis, of which between 38 to 77% of the primers were able to amplify products in these Lycoris species. Furthermore, 20 and 10 amplification products were selected for sequencing verification in L. aurea and L. radiata respectively. All products were validated as expected SSRs. In addition, 15 SSRs, including 10 sequence-verified and 5 unverified SSRs were selected and used to evaluate the genetic diversity of seven L. radiata lines. Among these, there were three sterile lines, three fertile lines and one line represented by the offspring of one fertile line. Unweighted pair group method with arithmetic mean analysis (UPGMA) demonstrated that the outgroup, L. aurea was separated from L. radiata lines and that the seven L. radiata lines were clustered into two groups, consistent with their fertility. Interestingly, even a dendrogram with 34 individuals representing the seven L. radiata lines was almost consistent with fertility.

CONCLUSIONS

This study supplies a pool of potential 839 non-redundant SSR markers for genetic analysis of Lycoris genus, that present high amplification rate, transferability and efficiency, which will facilitate genetic analysis and breeding program in Lycoris.

Complete Chloroplast Genomes and Comparative Analyses of L. chinensis, L. anhuiensis, and L. aurea (Amaryllidaceae)

The genus Lycoris (about 20 species) includes important medicinal and ornamental plants. Due to the similar morphological features and insufficient genomic resources, germplasm identification and molecular phylogeny analysis are very limited. Here, we sequenced the complete chloroplast genomes of L. chinensis, L. anhuiensis, and L. aurea; they have very similar morphological traits that make it difficult to identify. The full length of their cp genomes was nearly 158k bp with the same guanine-cytosine content of 37.8%. A total of 137 genes were annotated, including 87 protein-coding genes, 42 tRNAs, and eight rRNAs. A comparative analysis revealed the conservation in sequence size, GC content, and gene content. Some variations were observed in repeat structures, gene expansion on the IR-SC (Inverted Repeat-Single-Copy) boundary regions. Together with the cpSSR (chloroplast simple sequence repeats), these genetic variations are useful to develop molecular markers for germplasm identification. Phylogenetic analysis showed that seven Lycoris species were clustered into a monophyletic group, and closed to Narcissus in Amaryllidaceae. L. chinensis, L. anhuiensis, and L. longituba were clustered together, suggesting that they were very likely to be derived from one species, and had the same ancestor with L. squamigera. Our results provided information on the study of genetic diversity, origins or relatedness of native species, and the identification of cultivars.

Chemical and Biological Aspects of Montanine-Type Alkaloids Isolated from Plants of the Amaryllidaceae Family

Plants of the Amaryllidaceae family are promising therapeutic tools for human diseases and have been used as alternative medicines. The specific secondary metabolites of this plant family, called Amaryllidaceae alkaloids (AA), have attracted considerable attention due to their interesting pharmacological activities. One of them, galantamine, is already used in the therapy of Alzheimer’s disease as a long acting, selective, reversible inhibitor of acetylcholinesterase. One group of AA is the montanine-type, such as montanine, pancracine and others, which share a 5,11-methanomorphanthridine core. So far, only 14 montanine-type alkaloids have been isolated. Compared with other structural-types of AA, montanine-type alkaloids are predominantly present in plants in low concentrations, but some of them display promising biological properties, especially in vitro cytotoxic activity against different cancerous cell lines. The present review aims to summarize comprehensively the research that has been published on the Amaryllidaceae alkaloids of montanine-type.

Transcriptome and Metabolite Conjoint Analysis Reveals that Exogenous Methyl Jasmonate Regulates Monoterpene Synthesis in Grape Berry Skin

Monoterpene is one of the important sources of varietal aroma, which provides a strong floral and fruity aroma in wines. Methyl jasmonate (MeJA) affects plant secondary metabolism. However, the regulatory mechanisms of monoterpene biosynthesis after MeJA application on grapes are not illuminated. In the present study, 10 mM MeJA was used as treatments in Italian Riesling grape at the preveraison stage in different ways, including grape cluster soaking, foliar spraying, and whole vine spraying, designated as T1, T2, and T3, respectively, while a blank group was used as the control (CK). HS-SPME/GC-MS and transcriptome sequencing analysis were performed to investigate the effect of exogenous MeJA on monoterpene synthesis in grape berry skin. The results of GC-MS showed that the application of MeJA induced the accumulation of volatile monoterpenes in grape berry skin, especially linalool, α-terpineol, and oxides. In addition, transcriptome analysis showed that differentially expressed genes were increased from T2 to T3 to T1 compared with CK, and significantly enriched in JA and monoterpene synthesis pathways. T1 application significantly upregulated the mRNA expression levels of LOX2S, AOS, OPR, and JMT involved in the JA biosynthesis pathway, as well as DXS, HMGCR, TPS14, and α-terpineol synthesis genes involved in the monoterpene synthesis pathway compared with T2, T3, and CK. Thus, grape cluster soaking treatment with MeJA could greatly activate volatile monoterpene synthesis. The results will deeply increase our understanding of the monoterpene biosynthesis of grape berry skin in response to MeJA.

Identification, molecular characterization and expression of JAZ genes in Lycoris aurea

Jasmonates (JAs) are key phytohormones involved in regulation of plant growth and development, stress responses, and secondary metabolism. It has been reported that treatments with JAs could increase the contents of Amaryllidaceae alkaloids in Amaryllidaceae plants. Jasmonate ZIM (zinc-finger inflorescence meristem) domain (JAZ) proteins are key components in JA signal processes. However, JAZ proteins have not been characterized in genus Lycoris. In this study, we identified and cloned seven differentially expressed JAZ genes (namely LaJAZ1–LaJAZ7) from Lycoris aurea. Bioinformatic analyses revealed that these seven LaJAZ proteins contain the ZIM domain and JA-associated (Jas, also named CCT_2) motif. Quantitative reverse transcription polymerase chain reaction (qRT-PCR) analysis revealed that these LaJAZ genes display different expression patterns in L. aurea tissues, and most of them are inducible when treated with methyl jasmonate (MeJA) treatment. Subcellular localization assay demonstrated that LaJAZ proteins are localized in the cell nucleus or cytoplasm. In addition, LaJAZ proteins could interact with each other to form homodimer and/or heterodimer. The findings in this study may facilitate further functional research of the LaJAZ genes, especially the potential regulatory mechanism of plant secondary metabolites including Amaryllidaceae alkaloids in L. aurea.

Time-dependent behavior of phenylpropanoid pathway in response to methyl jasmonate in Scrophularia striata cell cultures

MeJA triggers a time-dependent behavior of the phenylpropanoid compounds. Plant cells produce a large number of metabolites in response to environmental factors. The cellular responses to environmental changes are orchestrated by signaling molecules, such as methyl jasmonate (MeJA). To understand how the MeJA changes the behavior of amino acids, carbohydrates, and phenylpropanoid compounds such as phenolic acids, phenylethanoid-glycosides, and flavonoids in Scrophularia striata cells; we monitored the metabolic responses for different times of exposure. In this study, we performed a time course analysis of metabolites and enzymes in S. striata cells exposed to MeJA (100 µM) and evaluated the metabolic flux towards carbon-rich secondary metabolites production. Moreover, we calculated the biosynthetic energy cost for free amino acids. Our results indicated that MeJA accelerates the sucrose degradation and directs the metabolic fluxes towards a pool of flavonoids and phenylethanoid glycosides through a change in enzyme behavior in the entry point and center of the phenylpropanoid pathway. MeJA also decreased and then raised the amino acid biosynthesis cost in S. striata cells in a time-dependent manner, indicating the cells evolve to utilize amino acids more economically by reducing cell growth. Finally, we classified the marked changes in the metabolites level and enzyme activities into three groups including early-, late-, and oscillatory-response groups to MeJA and summarized our findings as a model depicting pathway interactions during MeJA elicitation. Determination of metabolic levels in response to MeJA suggests that the changes in metabolic responses are time-dependent.

Transcriptome Analysis of Different Tissues Reveals Key Genes Associated With Galanthamine Biosynthesis in Lycoris longituba

L ycoris longituba is a traditional medicinal plant containing the bioactive compound galanthamine (Gal), a type of Amaryllidaceae alkaloid and can be used to treat Alzheimer's disease. However, research on its genome or transcriptome and associated genes in the biosynthetic pathway is incomplete. In this study, we estimated the nuclear genome size of this species to be 29.33 Gb by flow cytometry. Then, RNA sequencing of the leaves, roots, and bulbs of L. longituba was carried out. After de novo assembly, 474,589 all-transcripts and 333,440 all-unigenes were finally generated. In addition, the differentially expressed genes (DEGs) were identified, and genes involved in the galanthamine metabolic pathway encoding tyrosine decarboxylase (TYDC), phenylalanine ammonia-lyase (PAL), cinnamate 4-hydroxylase (C4H), p-coumarate 3-hydroxylase (C3H), norbelladine synthase (NBS), norbelladine 4'-O-methyltransferase (OMT), and noroxomaritidine synthase (CYP96T1) were detected and validated by real-time quantitative PCR analysis. One candidate gene, Lycoris longituba O-Methyltransferase (LlOMT), was identified in the proposed galanthamine biosynthetic pathway. Sequence analysis showed that LlOMT is a class I OMT. LlOMT is localized in the cytoplasm, and biochemical analysis indicated that the recombinant LlOMT catalyzes norbelladine to generate 4'-O-methylnorbelladine. The protoplast transformation result showed that the overexpression of LlOMT could increase the Gal content. Our results indicate that LlOMT may play a role in galanthamine biosynthesis in L. longituba. This work provides a useful resource for the metabolic engineering of Amaryllidaceae alkaloids.

De novo transcriptomic analysis of light-induced flavonoid pathway, transcription factors in the flower buds of Lonicera japonica

KEY MESSAGE

Transcriptomic analysis of the relationship between gene expression patterns and flavonoid contents in the flower buds of Lonicera japonica under light-induced conditions, especially the flavonoid pathway genes and transcription factors.

ABSTRACT

Flos Lonicerae Japonicae (FLJ), the flower buds of Lonicera japonica Thunb., has been used to treat some human diseases including severe respiratory syndromes and hand-foot-and-mouth diseases owing to its putative antibacterial, and antiviral effects. Luteoloside is a flavonoid that is used by the Chinese Pharmacopoeia to evaluate the quality of FLJ. Light is an important environmental factor that affects flavonoid biosynthesis in the flower buds of L. japonica. However, how light triggers increases in flavonoid production remains unclear. To enhance our understanding of the mechanism involved in light-regulated flavonoid biosynthesis, we sequenced the transcriptomes of L. japonica exposed to three different light conditions: 100% light intensity (CK), 50% light intensity (LI50), and 25% light intensity (LI25) using an Illumina HiSeq 4000 System. A total of 77,297 unigenes with an average length of 809 bp were obtained. Among them, 43,334 unigenes (56.06%) could be matched to at least one biomolecular database. Additionally, 4188, 1545 and 1023 differentially expressed genes (DEGs) were identified by comparative transcriptomics LI25-vs-CK, LI50-vs-CK, and LI25-vs-LI50, respectively. Of note, genes known to be involved in flavonoid biosynthesis, such as 4-coumarate coenzyme A ligase (4CL), and chalcone synthase (CHS) were up-regulated. In addition, a total of 1649 transcription factors (TFs) were identified and divided into 58 TF families; 98 TFs exhibited highly dynamic changes in response to light intensity. Quantitative real-time PCR (qRT-PCR) was used to test the expression profiles of the RNA sequencing (RNA-Seq) data. This study offers insight into how transcriptional expression pattern is influenced by light in the flower buds of L. japonica, and will enhance the understanding of molecular mechanisms of flavonoid biosynthesis in response to light in L. japonica.

Amaryllidaceae alkaloids: identification and partial characterization of montanine production in Rhodophiala bifida plant

Rhodophiala bifida (R. bifida) is a representative of the Amaryllidaceae plant family and is rich in montanine, an alkaloid with high pharmaceutical potential. Despite the interest in these compounds, many steps of the biosynthetic pathway have not been elucidated. In this study, we identified the alkaloids produced in different organs of R. bifida under different growth conditions, set up the conditions for in vitro R. bifida regeneration and initiated the molecular characterization of two R. bifida genes involved in alkaloids biosynthesis: the Norbelladine 4'-O-Methyltransferase (RbN4OMT) and the Cytochrome P450 (RbCYP96T). We show that montanine is the main alkaloid produced in the different R. bifida organs and developed a direct organogenesis regeneration protocol, using twin-scale explants cultivated on media enriched with naphthalene acetic acid and benzyladenine. Finally, we analyzed the RbN4OMT and RbCYP96T gene expressions in different organs and culture conditions and compared them to alkaloid production. In different organs of R. bifida young, adult and regenerated plants, as well as under various growing conditions, the transcripts accumulation was correlated with the production of alkaloids. This work provides new tools to improve the production of this important pharmaceutical compound and for future biotechnological studies.

Cloning and characterization of a tyrosine decarboxylase involved in the biosynthesis of galanthamine in Lycoris aurea

Background

Galanthamine, one kind of Amaryllidaceae alkaloid extracted from the Lycoris species, is used in the treatment of Alzheimer’s disease. In regards to medical and economic importance, the biosynthesis and regulatory mechanism of the secondary metabolites in Lycoris remain uninvestigated.

Methods

BLAST was used to identify the sequence of tyrosine decarboxylase in the transcriptome of Lycoris aurea (L’Hér) Herb. The enzyme activity of this TYDC was determined by using heterologous expressed protein in the Escherichia coli cells. The related productive contents of tyramine were detected using High Performance Liquid Chromatography (HPLC). According to the available micro RNA sequencing profiles and degradome database of L. aurea, microRNA396 were isolated, which targets to LaTYDC1 and RNA Ligase-Mediated-Rapid Amplification of cDNA Ends (RLM-RACE) were used to confirm the cleavage. The expression levels of miR396 and LaTYDC1 were measured using a quantitative real-time polymerase chain reaction (qRT-PCR).

Results

LaTYDC1 was mainly expressed in root, bulb, leaf and flower fitting the models for galanthamine accumulation. This decarboxylase efficiently catalyzes tyrosine to tyramine conversion. Under methyl jasmonate (MeJA) treatment, the expression of LaTYDC1 and the content of tyramine sharply increase. The use of RLM-RACE confirms that miR396 promotes the degradation of LaTYDC1 mRNA. Under MeJA treatment, the expression of miR396 was suppressed while the expression level of LaTYDC1 sharply increased. Following the increase of the miR396 transcriptional level, LaTYDC1 was significantly repressed.

Conclusion

LaTYDC1 participates in the biosynthesis of galanthamine, and is regulated by miR396. This finding also provides genetic strategy for improving the yield of galanthamine in the future.

Response of Bioactive Metabolite and Biosynthesis Related Genes to Methyl Jasmonate Elicitation in Codonopsis pilosula

Bioactive metabolites in Codonopsis pilosula are of particular interest as an immunostimulant. Methyl jasmonate (MeJA) plays an important role in the elicitation of metabolite biosynthesis. Here, we explored the response of metabolites to MeJA elicitation in C. pilosula adventitious roots and multiple shoots. The results showed that the biomass, polysaccharide, and lobetyolin content of adventitious roots exhibited the highest increases with 100 µmol·L⁻¹ MeJA at the 16th day of subculture, whereas the atractylenolide III (a terpenoid) content increased extremely with 50 µmol·L⁻¹ MeJA treatment at the 7th day of subculture. In addition, the biomass and lobetyolin content significantly increased at the 4th day after treatment. Similarly, the polysaccharide and lobetyolin content increased in multiple shoots. Further identification of different metabolites responding to MeJA by ¹H-NMR showed an extremely significant increase of the lobetyolinin level, which coincided with lobetyolin. Accordingly, the precursor, fatty acids, showed a highly significant decrease in their levels. Furthermore, a significant increase in β-d-fructose-butanol glycoside was detected, which was accompanied by a decrease in the sucrose level. Accordingly, the enzyme genes responsible for terpenoid and carbohydrate biosynthesis, CpUGPase, and CpPMK, were up regulated. In conclusion, MeJA promoted culture growth and accelerated bioactive metabolite accumulation by regulating the expression of the metabolite biosynthesis related genes, CpUGPase and CpPMK in C. pilosula.

Genome-Centered Metagenomics Analysis Reveals the Symbiotic Organisms Possessing Ability to Cross-Feed with Anammox Bacteria in Anammox Consortia

Although using anammox communities for efficient wastewater treatment has attracted much attention, the pure anammox bacteria are difficult to obtain, and the potential roles of symbiotic bacteria in anammox performance are still elusive. Here, we combined long-term reactor operation, genome-centered metagenomics, community functional structure, and metabolic pathway reconstruction to reveal multiple potential cross-feedings during anammox reactor start-up according to the 37 recovered metagenome-assembled genomes (MAGs). We found Armatimonadetes and Proteobacteria could contribute the secondary metabolites molybdopterin cofactor and folate for anammox bacteria to benefit their activity and growth. Chloroflexi-affiliated bacteria encoded the function of biosynthesizing exopolysaccharides for anammox consortium aggregation, based on the partial nucleotide sugars produced by anammox bacteria. Chlorobi-affiliated bacteria had the ability to degrade extracellular proteins produced by anammox bacteria to amino acids to affect consortium aggregation. Additionally, the Chloroflexi-affiliated bacteria harbored genes for a nitrite loop and could have a dual role in anammox performance during reactor start-up. Cross-feeding in anammox community adds a different dimension for understanding microbial interactions and emphasizes the importance of symbiotic bacteria in the anammox process for wastewater treatment.

De novo transcriptome sequencing and analysis of genes related to salt stress response in Glehnia littoralis

Soil salinity is one of the major environmental stresses affecting plant growth, development, and reproduction. Salt stress also affects the accumulation of some secondary metabolites in plants. Glehnia littoralis is an endangered medicinal halophyte that grows in coastal habitats. Peeled and dried Glehnia littoralis roots, named Radix Glehniae, have been used traditionally as a Chinese herbal medicine. Although Glehnia littoralis has great ecological and commercial value, salt-related mechanisms in Glehnia littoralis remain largely unknown. In this study, we analysed the transcriptome of Glehnia littoralis in response to salt stress by RNA-sequencing to identify potential salt tolerance gene networks. After de novo assembly, we obtained 105,875 unigenes, of which 75,559 were annotated in public databases. We identified 10,335 differentially expressed genes (DEGs; false discovery rate <0.05 and |log₂ fold-change| ≥ 1) between NaCl treatment (GL2) and control (GL1), with 5,018 upregulated and 5,317 downregulated DEGs. To further this investigation, we performed Gene Ontology (GO) analysis and the Kyoto Encyclopaedia of Genes and Genomes (KEGG) pathway analysis. DEGs involved in secondary metabolite biosynthetic pathways, plant signal transduction pathways, and transcription factors in response to salt stress were analysed. In addition, we tested the gene expression of 15 unigenes by quantitative real-time PCR (qRT-PCR) to confirm the RNA-sequencing results. Our findings represent a large-scale assessment of the Glehnia littoralis gene resource, and provide useful information for exploring its molecular mechanisms of salt tolerance. Moreover, genes enriched in metabolic pathways could be used to investigate potential biosynthetic pathways of active compounds by Glehnia littoralis.

Transcriptome Analysis of JA Signal Transduction, Transcription Factors, and Monoterpene Biosynthesis Pathway in Response to Methyl Jasmonate Elicitation in Mentha canadensis L

Mentha canadensis L. has important economic value for its abundance in essential oils. Menthol is the main component of M. canadensis essential oils, which is certainly the best-known monoterpene for its simple structure and wide applications. However, the regulation of menthol biosynthesis remains elusive in M. canadensis. In this study, transcriptome sequencing of M. canadensis with MeJA treatment was applied to illustrate the transcriptional regulation of plant secondary metabolites, especially menthol biosynthesis. Six sequencing libraries were constructed including three replicates for both control check (CK) and methyl jasmonate (MeJA) treatment and at least 8 Gb clean bases was produced for each library. After assembly, a total of 81,843 unigenes were obtained with an average length of 724 bp. Functional annotation indicated that 64.55% of unigenes could be annotated in at least one database. Additionally, 4430 differentially expressed genes (DEGs) with 2383 up-regulated and 2047 down-regulated transcripts were identified under MeJA treatment. Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment indicated that "Monoterpenoid biosynthesis" was one of the most significantly enriched pathways in metabolism. Subsequently, DEGs involved in JA signal transduction, transcription factors, and monoterpene biosynthesis were analyzed. 9 orthologous genes involved in menthol biosynthesis were also identified. This is the first report of a transcriptome study of M. canadensis and will facilitate the studies of monoterpene biosynthesis in the genus Mentha.

Agrobiotechnology Goes Wild: Ancient Local Varieties as Sources of Bioactives

The identification and use of species that have best adapted to their growth territory is of paramount importance to preserve biodiversity while promoting sustainable agricultural practices. Parameters including resistance to natural conditions (biotic and abiotic risk factors), biomass and fruit productivity, and phytochemical content with nutraceutical potential, could be used as quantitative markers of the adaptability of plants to wild environments characterized by minimal human impact. Ancient varieties, which are plant varieties growing in regional territories and not destined for market distribution, are a source of unique genetic characters derived from many years of adaptation to the original territory. These plants are often more resistant to biotic and abiotic stresses. In addition, these varieties have a high phytochemical (also known as bioactives) content considered health-beneficial. Notably, the content of these compounds is often lower in commercial cultivars. The use of selected territorial varieties according to the cultivation area represents an opportunity in the agricultural sector in terms of biodiversity preservation, environmental sustainability, and valorization of the final products. Our survey highlights the nutraceutical potential of ancient local varieties and stresses the importance of holistic studies (-omics) to investigate their physiology and secondary metabolism.

Molecular Cloning and Characterization of a meta/para-O-Methyltransferase from Lycoris aurea

O-methyltransferases (OMTs) have been demonstrated to play key roles in the biosynthesis of plant secondary metabolites, such as alkaloids, isoprenoids, and phenolic compounds. Here, we isolated and characterized an OMT gene from Lycoris aurea (namely LaOMT1), based on our previous transcriptome sequencing data. Sequence alignment and phylogenetic analysis showed that LaOMT1 belongs to the class I OMT, and shares high identity to other known plant OMTs. Also, LaOMT1 is highly identical in its amino acid sequence to NpN4OMT, a norbelladine 4′-OMT from Narcissus sp. aff. pseudonarcissus involved in the biosynthesis of Amaryllidaceae alkaloids. Biochemical analysis indicated that the recombinant LaOMT1 displayed both para and metaO-methylation activities with caffeic acid and 3,4-dihydroxybenzaldehyde, and showed a strong preference for the meta position. Besides, LaOMT1 also catalyzes the O-methylation of norbelladine to form 4′-O-methylnorbelladine, which has been demonstrated to be a universal precursor of all the primary Amaryllidaceae alkaloid skeletons. The results from quantitative real-time PCR assay indicated that LaOMT1 was ubiquitously expressed in different tissues of L. aurea, and its highest expression level was observed in the ovary. Meanwhile, the largest concentration of lycorine and galanthamine were found in the ovary, whereas the highest level of narciclasine was observed in the bulb. In addition, sodium chloride (NaCl), cold, polyethylene glycol (PEG), sodium nitroprusside (SNP), and methyl jasmonate (MeJA) treatments could significantly increase LaOMT1 transcripts, while abscisic acid (ABA) treatment dramatically decreased the expression level of LaOMT1. Subcellular localization showed that LaOMT1 is mainly localized in cytoplasm and endosome. Our results in this study indicate that LaOMT1 may play a multifunctional role, and lay the foundation for Amaryllidaceae alkaloid biosynthesis in L. aurea.

Transcriptome and miRNA analyses of the response to Corynespora cassiicola in cucumber

Cucumber (Cucumis sativus L.) target leaf spot (TLS), which is caused by the fungus Corynespora cassiicola (C. cassiicola), seriously endangers the production of cucumber. In this assay, we performed comprehensive sequencing of the transcriptome and microRNAs (miRNAs) of a resistant cucumber (Jinyou 38) during C. cassiicola inoculation using the Illumina NextSeq 500 platform. The possible genes related to the response to C. cassiicola were associated with plant hormones, transcription factors, primary metabolism, Ca²⁺ signaling pathways, secondary metabolism and defense genes. In total, 150 target genes of these differentially expressed miRNAs were predicted by the bioinformatic analysis. By analyzing the function of the target genes, several candidate miRNAs that may be related to the response to C. cassiicola stress were selected. We also predicted 7 novel miRNAs and predicted their target genes. Moreover, the expression patterns of the candidate genes and miRNAs were tested by quantitative real-time RT-PCR. According to the analysis, genes and miRNAs associated with secondary metabolism, particularly the phenylpropanoid biosynthesis pathway, may play a major role in the resistance to C. cassiicola stress in cucumber. These results offer a foundation for future studies exploring the mechanism and key genes of resistance to cucumber TLS.

Characterization of phenolics, amino acids, fatty acids and antioxidant activity in pulp and seeds of high altitude Himalayan crab apple fruits (Malus baccata)

Phytochemicals in fruits and vegetables have achieved immense significance owing to the increasing evidence which signifying their activity for antioxidant and prevention of chronic diseases. The amount of phloretin (88.39 µg mg^-1) and phloridzin (83.03 µg mg^-1) were found to be higher among other phenolics determined using UPLC. DPPH, ABTS⁺, metal chelating and ·OH radical assays were used to evaluate antioxidant activity. Malus baccata pulp portion showed higher antioxidant activity than seed portion. HPLC analysis for free amino acids showed that serine (9.06 µg mg^-1), alanine (8.03 µg mg^-1), tyrosine (10.33 µg mg^-1), and cysteine (76.86 µg mg^-1) were only detected in pulp portion while seed comprised of histidine (3.96 µg mg^-1) only. Seed portion was also determined for their fatty acid composition including palmitic acid (0.89%), ethyl palmitate (0.56%), methyl petroselinate (0.90%) and linolein (3.93%) using GC-MS analysis. HPAEC technique detected fructose and sucrose in a fair amount of 21 and 17.3 mg g^-1 in pulp, while 9.4 and 4.24 mg g^-1 in seed portion, respectively. The present study suggested that M. baccata fruit is a rich source of phenolic and other chemical components which can be used in food products and nutraceutical formulations.

A Transcriptome Approach Toward Understanding Fruit Softening in Persimmon

Persimmon (Diospyros kaki Thunb.), which is a climacteric fruit, softens in 3-5 weeks after harvest. However, little is known regarding the transcriptional changes that underlie persimmon ripening. In this study, high-throughput de novo RNA sequencing was performed to examine differential expression between freshly harvested (FH) and softened (ST) persimmon fruit peels. Using the Illumina HiSeq platform, we obtained 259,483,704 high quality reads and 94,856 transcripts. After the removal of redundant sequences, a total of 31,258 unigenes were predicted, 1,790 of which were differentially expressed between FH and ST persimmon (1,284 up-regulated and 506 down-regulated in ST compared with FH). The differentially expressed genes (DEGs) were further subjected to KEGG pathway analysis. Several pathways were found to be up-regulated in ST persimmon, including "amino sugar and nucleotide sugar metabolism." Pathways down-regulated in ST persimmon included "photosynthesis" and "carbon fixation in photosynthetic organisms." Expression patterns of genes in these pathways were further confirmed using quantitative real-time RT-PCR. Ethylene gas production during persimmon softening was monitored with gas chromatography and found to be correlated with the fruit softening. Transcription involved in ethylene biosynthesis, perception and signaling was up-regulated. On the whole, this study investigated the key genes involved in metabolic pathways of persimmon fruit softening, especially implicated in increased sugar metabolism, decreased photosynthetic capability, and increased ethylene production and other ethylene-related functions. This transcriptome analysis provides baseline information on the identity and modulation of genes involved in softening of persimmon fruits and can underpin the future development of technologies to delay softening in persimmon.