
==== Front
PLoS One
PLoS One
plos
PLOS ONE
1932-6203
Public Library of Science San Francisco, CA USA

10.1371/journal.pone.0309978
PONE-D-24-13936
Research Article
Physical Sciences
Chemistry
Chemical Compounds
Alkaloids
Biology and Life Sciences
Biochemistry
Biosynthesis
Biology and Life Sciences
Genetics
Gene Expression
Biology and Life Sciences
Agriculture
Agrochemicals
Fertilizers
Physical Sciences
Chemistry
Chemical Compounds
Organic Compounds
Steroids
Physical Sciences
Chemistry
Organic Chemistry
Organic Compounds
Steroids
Biology and Life Sciences
Biochemistry
Proteins
Post-Translational Modification
Hydroxylation
Biology and Life Sciences
Computational Biology
Genome Analysis
Transcriptome Analysis
Biology and Life Sciences
Genetics
Genomics
Genome Analysis
Transcriptome Analysis
Physical Sciences
Chemistry
Chemical Compounds
Organic Compounds
Terpenes
Physical Sciences
Chemistry
Organic Chemistry
Organic Compounds
Terpenes
Transcriptomic analysis reveals effects of fertilization towards growth and quality of Fritillariae thunbergii bulbus
Transcriptomic analysis of Fritillariae thunbergii bulbus under different fertilization treatments
https://orcid.org/0009-0006-3455-787X
Huang Luman Writing – original draft 1
Liang Shuang Data curation 1
Luo Lei Formal analysis 2
Wu Mengmin Software 1
https://orcid.org/0009-0009-7194-124X
Fu Hongwei Data curation Supervision Validation 1 *
https://orcid.org/0000-0002-6952-1096
Zhong Zhuoheng Writing – review & editing 1 *
1 College of Life Sciences and Medicine, Zhejiang Sci-Tech University, Hangzhou, P.R. China
2 Zhejiang Institute for Food and Drug Control, Hangzhou, P.R. China
Li Minhui Editor
Baotou Medical College, MONGOLIA
Competing Interests: The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence this work reported in this paper.

* E-mail: fhw668@zju.edu.cn (HF); zhongzhh@zstu.edu.cn (ZZ)
20 9 2024
2024
19 9 e03099787 4 2024
17 8 2024
© 2024 Huang et al
2024
Huang et al
https://creativecommons.org/licenses/by/4.0/ This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Fritillariae thunbergii Bulbus (FTB) is a traditional Chinese medicine that has been widely cultivated for its expectorant, antitussive, antiasthmatic, antiviral, and anticancer properties. The yield and quality of F. thunbergii are influenced by cultivation conditions, such as the use of fertilizers. However, the optimal type of fertilizers for maximum quality and yield and underlying mechanisms are not clear. We collected F. thunbergii using raw chicken manure (RC), organic fertilizer (OF), and plant ash (PA) as the base fertilizer in Pan’an County, Jinhua City, Zhejiang Province as experimental materials. The combined results of HPLC-ELSD detection and yield statistics showed that the F. thunbergii with OF application was the best, with the content of peimine and peiminine reaching 0.0603% and 0.0502%, respectively. In addition, the yield was 2.70 kg/m2. Transcriptome analysis indicated that up-regulation of the ABA signaling pathway might promote bulb yield. Furthermore, putative key genes responsible for steroidal alkaloid accumulation were identified. These results provided guiding significance for the rational fertilization conditions of F. thunbergii as well as the basis for the exploration of functional genes related to the alkaloid biosynthesis pathway.

the Research Initiation Funding of Zhejiang Sci-Tech University No. 19042112-Y National Support Funding for Reform and Development of Local Universities No. 303013-2021-0007 Present work was supported by the Research Initiation Funding of Zhejiang Sci-Tech University (Grant number: 19042112-Y; 22042024-Y). Data AvailabilityThe raw data of RNA-Seq have been submitted to the the NCBI Sequence Read Archive (SRA) database (Accession: PRJNA950562).
Data Availability

The raw data of RNA-Seq have been submitted to the the NCBI Sequence Read Archive (SRA) database (Accession: PRJNA950562).
==== Body
pmcIntroduction

Fritillaria (“Beimu” in Chinese) genus, belonging to the Liliaceae family, was ranked as one kind of the most traditional medicinal plants in China [1]. So far, about 130 species of this genus are found to be distributed in temperate regions of the northern hemisphere [2]. In the 2020 edition of the Chinese Pharmacopoeia, the bulbus of 11 cultivars from Fritillaria were officially recorded and classified into five species, namely Fritillariae Cirrhosae Bulbus (FCB), Fritillariae Hupehensis Bulbus (FHB), Fritillariae Ussuriensis Bulbus (FUB), Fritillariae Pallidiflorae Bulbus (FPB), and Fritillariae Thunbergii Bulbus (FTB) [3]. F. thunbergii, also known as Zhe Bei, is one of the famous "Zhejiang Eight Flavors". Modern medical research reported that it had functions such as cough relieving, phlegm relieving, asthma relieving, analgesic and anti-inflammatory, anti-ulcer, antioxidant, anti-tumor, and reversing tumor cell resistance [4–8]. Its main chemical components include alkaloids, polysaccharides, total saponins, flavonoids, and volatile components [7]. According to research, steroidal alkaloids have important medicinal value in eliminating phlegm and relieving cough, reducing blood pressure and promoting blood circulation, relieving pain and ulcers, anti-inflammatory and antioxidant effects, and anti-tumor effects [9]. At present, the synthesis pathway of steroidal alkaloids in F. thunbergii is not very clear.

With the widespread application of FTB, wild resources are scarce. To meet the increasing market demand, it has been introduced as a domestic variety, and the planting areas are constantly expanding. Pan’an County is one of the main production areas of F. thunbergii [10]. However, differences in fertilization conditions, soil, and other factors have led to significant differences in the quality and yield of F. thunbergii. The "Technical Regulations for the Production of F. thunbergii" stipulated that raw chicken manure, commercial organic fertilizer, and plant ash could be used as fertilizer for F. thunbergii. However, there is still a lack of research on these three fertilizer applications on the yield and quality of F. thunbergii.

Transcriptomic information from many medicinal herbs has been acquired for investigation of differential alterations at the gene expression level as well as potential biosynthetic pathways of secondary metabolites [11]. According to previous analysis of transcriptome data, the steroid alkaloid biosynthetic pathway was divided into three stages: (1) Terpenoid backbone biosynthesis, where isopentenyl diphosphate (IPP) and dimethylallyl diphosphate (DMAPP) produced from mevalonate (MVA) and 2-C-methyl-D-erythritol-4-phosphate (MEP) pathways, respectively, were converted into farnesyl pyrophosphate (FPP) under the action of farnesyl pyrophaophate synthase (FPS) [12]; (2) Steroid circular skeleton precursor biosynthesis, two FPPs were catalyzed by squalene synthase (SQS) and squalene epoxidase (SQE) to form 2,3-oxidosqualene (OS) [13]; (3) Steroid compounds biosynthesis, several members of plant cytochrome P450 (CYP450) had also been confirmed to participate in the biosynthesis of steroid compounds [14], but the specific functions have yet to be testified. The varied enzymes involved in the biosynthesis of steroid alkaloids in FTB still need to be further studied.

In addition, many researchs have also proved that altitude have an significant impact on secondary metabolism in medicinal plants [15, 16]. Indeed, omics have been applied to the in vitro reproduction techniques of high-value Himalayan Fritillaria species [17]. Similarly, appropriate organic fertilizer was beneficial for the accumulation of secondary metabolites of Salvia miltiorrhiza Bunge and the improvement of yield [18], different fertilization conditions have an impact on the yield and quality of Ligusticum chuanxiong Hort [19], Herba Artimisiae Sieversianae [20]. However, there is little research on the effects of different fertilization conditions on the yield and alkaloid content of Bulbs of F. thunbergii. In this study, we collected FTB obtained with three kinds of fertilizer from Pan’an County in Zhejiang province in China. High-performance liquid chromatography-evaporative light scattering detector (HPLC-ELSD) and high-throughput RNA-Seq technology were performed to jointly analyze variations among different groups at both metabolite and gene expression levels, with the hypothesis of potential genes related to steroidal alkaloid biosynthesis pathway in F. thunbergii and providing a genetic basis for subsequent research on steroid alkaloid biosynthesis. It has a guiding significance for rational fertilization of F. thunbergii to improve yield and quality.

Materials and methods

Plant materials

Triennial Fritillariae thunbergii bulbs (cultivar ‘Zhebei 1’, identified by Prof. Bizeng Mao from Zhejiang University) were cultivated and collected on the same day in a natural habitat with fertile soil in Pan’an, Zhejiang, China (located at 28°99’N, 120°62’E, altitude 600 m). Raw chicken manure (RC), organic fertilizer (OF), and plant ash (PA) were separately applied as the base fertilizer and the applications were maintained constant during cultivation. Fritillaria thunbergii bulb yield was determined by harvesting the plants in an area of 2 m2 in each plot. For collecting samples for analysis, all the bulbs were collected biological replicates to randomly select different regions in the north, south, east, and west of the same production area, and mix well on the same day and washed with pure water, after which one portion was quickly frozen in liquid nitrogen and stored at -80°C until processing, and these samples were used for RNA extraction and library construction. For each fertilization conditions, three biological replicas were prepared for each sample.

Soil and fertilizer sources

The soil is sandy loam soil, with a turtle-shaped border. Organic fertilizer (SITEWO, Anhui, China). Chicken manure and plant ash were made by local farmers. The fertilization amount for different fertilizers was 5 kg/m2.

Chemical reagents

Strong ammonia solution (28%, Nanjing Chemical Reagent Co., Ltd), Chloroform (≥99.5%, Nanjing Chemical Reagent Co., Ltd), Methanol (≥99.5%, Hangzhou Gaojing Fine Chemical Co., Ltd).

Determination of alkaloid contents

A method based on the Chinese Pharmacopoeia was used to extract peimine and peiminine from FTB. Briefly, dried bulbs were ground into powder and screened through a 65-mesh sieve, 2 g of the dried powder was accurately weighed into a flask (with three biological replicates per sample), soaked in 4 mL strong ammonia solution for 1 h, and accurately added with 40 mL of a mixed solution of chloroform-methanol (4:1, v/v). The mixture was heated and refluxed in a water bath at 80°C for 2 h. Precisely measured supernatant (10 mL) was put in an evaporating dish and evaporated to dryness, 2 mL of methanol was added to dissolve the residue and accurately transferred to a 2 mL-volumetric-flask. Then the solution was centrifuged for 10 min at 13 400 g, and the supernatant was collected and filtrated through 0.22 μm membrane filters (Jinteng, Tianjin, China) before HPLC-ELSD analysis.

An aliquot of the distinct filtrates (10 μL) was injected into an LC-20AT HPLC system equipped with an Alltech 2000ES ELSD Evaporative Light Scattering Detector and an automatic sampler (Shimadzu, Japan). The chromatographic separation was applied on a kromasil C18 column (250 mm × 4.6 mm, 5 μm) at the column temperature of 35°C and the flow rate is 1.0 mL/min. The mobile phase consisted of methanol (A) and 0.03% diethylamine solution (B)for chromatographic separation using isocratic elution parameters.

Total RNA preparation

For RNA-seq of FTB from different fertilization conditions, the total RNA was extracted using the SteadyPure Plant RNA Extraction Kit (Accurate Biotechnology, Hunan, China). The integrity and purity of the total RNA quality were determined by the PowerPac™ Basic electrophoresis apparatus (Bio-Rad, California, USA), and the content was quantified using the NanoDrop spectrophotometer 1000 (Thermo Fisher, MA, USA). Samples from three fertilization conditions containing three replicates were applied for the construction sequencing library.

Library construction and sequencing

Library construction and sequencing were performed by Novogene Bioinformatics Technology Co. Ltd. (Beijing, China). The extracted total RNA as input material and reverse-transcribed to cDNA for building library by NEBNext® Ultra™ RNA Library Prep Kit for Illumina® (New England Biolabs, Beijing, China) following the manufacturer’s protocol. The modified cDNA fragments with a length of 370~420 bp were preferentially selected, and amplified by PCR, then the PCR products were purified using the AMPure XP beads (Beckman Coulter, Beverly, USA) to obtain cDNA libraries from different conditions of F. Thunbergii plants. Finally, the Illumina novaseq 6000 platform was used for transcriptome sequencing, and 150 bp paired-end reads were generated by paired-end sequencing technology.

Transcriptome assembly and annotation

For assessment of the quality and reliability of the transcriptome sequencing data, the raw data was filtered to remove polynucleotide sequences, reads with adapters, and reads with nebulous bases. Then the high-quality clean reads were de novo assembled using Trinity (v2.6.6). Based on Trinity splicing, Corset software was used to aggregate redundant transcripts for a better detection rate of differentially expressed genes [21]. BUSCO software was used to evaluate the splicing quality of the spliced files [22]. Finally, all assembled transcripts were mapped to public databases by RSEM (v12.15) [23], and the database contains NCBI non-redundant protein sequences (Nr), NCBI non-redundant nucleotide sequences (Nt), Protein family (Pfam), euKaryotic Orthologous Groups (KOG), Swiss-Prot, Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) databases.

Differentially expressed genes analysis

The relative expression levels of each transcript were calculated according to fragments per kilobase of transcript per million mapped reads (FPKM) [24]. Before screening the differentially expressed genes (DEGs), unigenes with unstable expression in sample transcriptomes were removed. The FPKM value of three samples with p-value < 0.05 was set as the threshold for differentially expressed genes (DEGs). Furthermore, GOseq (1.10.0) and KOBAS (v2.0.12) software were used for GO function enrichment analysis and KEGG metabolic pathway enrichment analysis of DEGs, respectively [25, 26]. Furthermore, the clustered expression pattern was drawn through R (Version 3.0.3), Package ggplot2 and Package pheatmap.

Quantitative real time-polymerase chain reaction (qRT-PCR) analysis

To confirm the accuracy of the gene expression levels of DEGs obtained by transcript sequencing, qRT-PCR analysis was used. The primers used for qRT-PCR are shown in S1 Table. Total RNA was reverse-transcribed into cDNA using the Evo M-MLV RT Mix kit (Accurate Biotechnology, Hunan, China), then cDNA was used as the template for qRT-PCR using the SYBR® Green Premix Pro Taq HS q-PCR kit (Accurate Biotechnology, Hunan, China) on Applied BiosystemsTM 7500 Fast Dx Real-Time PCR Instrument (ABI, Carlsbad, USA). Actin was used as the reference gene, and the relative expression levels of each group (containing three biological replicates) were calculated based on the 2-ΔΔCt method [27].

Statistical analysis

All statistical analyses including steroid alkaloid contents and relative gene expressions were performed One-way ANOVA using SPSS (25.0.0).

Results

The contents of peimine and peiminine varied in FTB from different fertilization conditions

As the steroidal alkaloid peimine and peiminine are important medicinal ingredients in FTB, the amount of the two alkaloids in different fertilization conditions was quantified based on HPLC-ELSD analysis. As shown in Table 1, the contents of the two alkaloids in PA were the highest, with the content of peimine and peiminine were 0.0713% and 0.0615%, respectively. However, compared with the other two groups, the yield was the lowest (2.13 Kg/m2). The content of the two alkaloids in RC was the lowest, which just reached the pharmacopoeia standard (The total quantity of peimine and peiminine based on dried material content serves as a crucial quality control standard, with a minimum requirement of 0.080%.), but the yield was the highest (2.72 Kg/m2). The results indicated that, apart from RC which had a low concentration of pharmacologically active compounds, other FTBs produced in two conditions of Zhejiang province were excellent. Interestingly, obtained through comprehensive evaluations, the optimal fertilization condition of F. thunbergii was OF.

10.1371/journal.pone.0309978.t001 Table 1 Contents of two steroidal alkaloids in FTB from different fertilization conditions.

Sample	Peimine	Peiminine	Total content%	Production(Kg/m2)	
Content percentage%	Content percentage%	
RC	0.0452 ± 0.001c	0.0362 ± 0.0012c	0.0814 ± 0.0022c	2.72	
OF	0.0603 ± 0.0008b	0.0502 ± 0.0012b	0.1105 ± 0.002b	2.70	
PA	0.0713 ± 0.0011a	0.0615 ± 0.0018a	0.1328 ± 0.0031a	2.13	
Note.

a, b and c indicate significant differences in alkaloid content between different fertilization conditions (P < 0.05).

Transcriptome de novo assembly and functional annotations of unigenes

To investigate the difference of FTB from three fertilization conditions at the transcriptomic level, 9 RNA-seq libraries of FTB from three groups each with three biological replicates were constructed. Each library obtained an average of approximately 21 to 22 million raw reads. After the removal of adapters, uncertain bases, and low-quality raw sequences, a total of 328.14 million clean reads of FTB were obtained, ranging from 20.13 million to 23.27 million for each library. The Q20 and Q30 values (the percentage of bases with quality values ≥ 20 or 30, which refers to the error probability of the identified bases in the base calling process and used to measure the sequencing accuracy) of the nucleotide sequences reached above 90%, and the GC contents of each sample clean reads were around 50% (S2 Table). PCA analysis showed that the transcriptome data within the sample group had good repeatability (S1 Fig).

After de novo assembly of clean reads from 9 samples of FTB using Trinity software, a total of 139, 639 Unigenes were obtained through Corset clustering for subsequent comparative analysis. Unigenes with sequence lengths ranging from 300 bp to 500 bp, from 500 bp to 1000 bp, from 1000 bp to 2000 bp, and greater than 2000 bp exhibit a gradual decrease in the number of unigenes within the corresponding interval as the sequence length interval increases (S3 Table).

All the assembled unigenes were annotated according to NR, NT, Swiss-Prot, Pfam, GO, KOG and KEGG databases. The total number of unigenes with at least one annotation was 59, 428, making up approximately 50% of all unigenes. Among them, 7,400 unigenes were functionally annotated in the seven databases (Table 2). 47, 853 unigenes were annotated at most in Nr database, and Nr homologous species distribution analysis indicated that the highly homologous species related to F. thunbergii, such as Elaeis guineensis (6, 272; 13.11%), Phoenix dactylifera (5, 926; 12.38%), Quercus suber (5, 743; 11.09%), and Asparagus officinalis (3, 303; 6.90%) (S2A Fig). According to GO clustering analysis, 33, 948 unigenes were annotated (S2B Fig).

10.1371/journal.pone.0309978.t002 Table 2 Unigenes functional annotation of the FTB transcriptome.

Database	Unigenes	Percentage%	
NR	47853	41.53%	
NT	24581	21.33%	
Swiss-Prot	35293	30.63%	
PFAM	33951	29.47%	
GO	33948	29.46%	
KOG	11079	9.62%	
KEGG	17286	15.00%	
All Database-Annotated	7400	6.42%	
At least one Database-Annotated	59428	51.58%	

To identify involved biological pathways of the transcripts, 17, 286 unigenes were identified in 351 KEGG pathways with the most abundant number in metabolic pathways (S2C Fig). Among the metabolic pathways, several unigenes of secondary metabolic biosynthetic pathways might be involved in regulating the biosynthesis of steroid alkaloids, including 134 unigenes of terpenoid backbone biosynthesis (ko00900), 16 unigenes of sesquiterpenoid and triterpenoid biosynthesis (ko00909), 68 unigenes of steroid biosynthesis (ko00100), 47 unigenes of diterpenoid biosynthesis (ko00904) and 81 unigenes of ubiquinone and other terpenoid-quinone biosynthesis (ko00130).

Analysis of differentially expressed genes (DEGs)

A total of 16, 810 DEGS were screened according to the FPKM of three groups with p-value < 0.05. Firstly, we analyzed these DEGs based on GO annotation. The statistical analysis results are shown (Fig 1), where differentially expressed genes are distributed in three major categories: "biological process", "cellular component", and" molecular function". According to the P-value enriched by GO, it was found that in the category of "biological processes", the enrichment degree of "protein metabolic process" (1637) and "ribonucleoprotein complex biogenesis" (375) was the most significant. In the category of "Cellular component", the enrichment degree of "ribonucleoprotein complex" (434) and "macromolecular complex" (1923) is the most significant and the main taxonomic unit. In the category of "Molecular function", the enrichment level of "protein binding function" (1157) and "RNA binding function" (316) is the most significant. The results indicate that the GO subclass with the most significant enrichment of differential genes is mainly related to protein function and the binding function of various genetic materials.

10.1371/journal.pone.0309978.g001 Fig 1 GO functional classification of total DEGs.

The y-axis indicates the different categories. The x-axis indicates the gene ratio of enrichment analysis, The size of circle indicates the number of differential genes. The color of circle indicates the P-value.

KEGG annotation analysis illustrated that DEGs were mainly related to pathways such as terpenoid backbone biosynthesis (ko00900), sulfur metabolism (ko00920), isoquinoline alkaloid biosynthesis (ko00950), steroid biosynthesis (ko00100), stilbenoid, diarylheptanoid and gingerol biosynthesis (ko00945), tropane, piperidine and pyridine alkaloid biosynthesis (ko00960), carotenoid biosynthesis (ko00906), and zeatin biosynthesis (ko00908), which indicated that expression levels of genes related to secondary metabolic pathways in different FTBs were discrepant (Fig 2).

10.1371/journal.pone.0309978.g002 Fig 2 KEGG functional classification of total DEGs.

The y-axis indicates the different categories. The x-axis indicates the number of DEGs in a category.

Next, five enrichment pathways with the highest number of DEGs were selected and involved in regulating plant life activities. The plant hormone signal transduction pathway (ko04075) in the "Environmental Information Processing" branch is an important pathway for regulating plant growth [28] (Fig 3A, S4 Table). The "plant pathogen interaction pathway (ko04626)" in the "organic system" branch mainly includes immune responses activated by pathogen-related molecular features and immune responses activated by effector proteins [29] (Fig 3B, S5 Table). The plant immune system can resist the invasion of most natural pathogens through these two defense mechanisms [30]. In the "genetic information processing" branch, the "mRNA monitoring pathway (ko03015)" is also known as uninterrupted decay (Fig 3C, S6 Table), which can effectively prevent the synthesis of other abnormal proteins while eliminating some harmful effects [31]; the "Ubiquitin mediated protein hydrolysis (ko04120)" is the main pathway for intracellular protein degradation (Fig 3D, S7 Table) and participates in multiple processes throughout the cell life cycle to regulate protein stability, function, activity, and subcellular localization [32, 33]; The "terpenoid skeleton biosynthesis (ko00900)" is involved in the generation of steroid alkaloid precursors (Fig 3E, S8 Table) and is highly important for steroid alkaloid synthesis [34].

10.1371/journal.pone.0309978.g003 Fig 3 The clustered expression profiles of five enriched pathways with the highest number of DEGs in the FTB of three fertilization conditions.

FPKM of each gene were normalized according to Z-score. Red indicates high gene expression, blue indicates low gene expression.

Pathway mapping analysis of DEGs involved in plant hormone signal transduction

Plant hormones are small plant organic molecules, which mainly contain auxin, cytokinin (CK), gibberellin (GA), abscisic acid (ABA), ethylene (ET), brassinosteroid (BR), jasmonic acid (JA), salicylic acid (SA) and strigolactones. They play an indispensable role in regulating plant growth and abiotic stress responses at very low concentrations in plants [35]. Based on KEGG annotation and gene expression files, a total of 123 DEGs were mapped to the plant hormone signal transduction pathway (ko04075) (Fig 4, S4 Table).

10.1371/journal.pone.0309978.g004 Fig 4 Expression heat maps of DEGs for the structural genes involved in plant hormone signal transduction (ko04075) from three fertilization conditions.

FPKM of each gene were normalized according to Z-score. Red indicates high gene expression, blue indicates low gene expression. Detailed gene expression profiles were shown in S4 Table.

ABA is the prominent medium of these response processes, the gene expression of SnRK2 and JAZ in PA was the highest. In addition, ABF, as a master regulatory transcription factor, belongs to a subfamily of bZIP proteins [36] and is also up-regulated in PA. In the JA signaling pathway, the expression profile of JAZ in PA and OF was relatively higher, Similarly, the majority of transcripts enriched in the JA signaling pathway were increased in PA, which was Similar to the expression patterns in the ABA pathway (Fig 4).

Identification and expression analysis of DEGs associated with steroid biosynthesis

To distinguish the presumable genes involved in the steroidal alkaloids biosynthesis in F. thunbergii, a possible biosynthetic pathway containing DEGs expression level has been made according to all DEGs enrichment results of KEGG pathways, it was initially synthesized via terpenoid backbone biosynthesis (ko00900) followed by steroid biosynthesis (ko00100). The results indicated that 62 DEGs were involved in terpenoid backbone biosynthesis (ko00900). Interestingly, these DEGs mainly show the same expression trend: low expression level in RC, and the expression level increases successively in PA and OF (Fig 5, S8 Table). which is consistent with the trend of differences in FTB content. DXS can catalyze the production of 1-deoxyglucose 5-phosphate (DXP) from pyruvate and 3-phosphate glyceraldehyde, which is the first rate limiting step in the catalytic MEP pathway. DXR can catalyze the synthesis of various important intermediate metabolite enzymes such as IPP, pyridoxine, and thiamine by DXP. Which is the second-rate limiting enzyme in the MEP pathway. Compared with RC, DXS and DXR in PA and OF were significantly upregulated.

10.1371/journal.pone.0309978.g005 Fig 5 Expression heat maps of DEGs in terpenoid backbone biosynthesis (ko00900) in the FTB from three fertilization conditions.

FPKM of each gene were normalized according to Z-score. The dashed arrow represents presumed terminal biosynthesis steps. AACT, acetyl CoA acetyltransferase; HMGS, 3-hydroxy-3-methylglutaryl-CoA synthase; HMGR,3-hydroxy-3-methylglutaryl-CoA reductase; MVK, mevalonate kinase; PMK, 5-phosphomevalonate kinase; MVD, mevalonate 5-diphosphate decarboxylas;DXS,1-deoxy-D-xylulose-5-phosphate synthase; DXR, 1-deoxy-D-xylulose-5-phosphate reductoisomerase; MCT, 2-C-methyl-D-erythritol-4-phosphate cytidylyltransferase; CMK, 4-diphosphocytidyl-2-C-methyl-D-erythritol kinase; MCS, 2-C-methyl-D-erythritol 2,4-cyclodiphosphate synthase; HDS,1-hydroxy-2-methyl-2-(E)-butenyl-4-diphosphate synthase; HDR, 4-hydroxy-3-methylbut-2-enyl diphosphate reductase; IDI, isopentenyl-diphosphate Delta-isomerase; GPS, geranyl diphosphate synthase; FPS, farnesyl pyrophosphate synthase.

According to a previous study, cholesterol was the precursor for the biosynthesis of steroidal alkaloids, which is important for the downstream accumulation of steroid alkaloids [34, 37]. As shown in Fig 6 (S9 Table), the pathway from 2,3-Oxidosqualene toward cholesterol in the KEGG database was already clear, we can directly obtain relevant information on DEGs in this pathway through enrichment analysis. the expression level of DEGs mainly is up-regulated in OF and PA, and a downregulation trend in RC. To some extent, it provides reliable evidence for the upregulation of genes related to steroid biosynthesis pathways.

10.1371/journal.pone.0309978.g006 Fig 6 Putative biosynthesis pathway from cholesterol toward steroidal alkaloids and schematic diagram of fold change in expression level of genes in the FTB of three fertilization conditions.

The dashed arrow represents the presumed terminal biosynthesis steps. FPKM of each gene were normalized according to Z-score. Red indicates high gene expression, blue indicates low gene expression. SQS, squalene synthase; SQE, squalene monooxygenase; LSS, lanosterol synthase; NSDHL, sterol-4alpha-carboxylate 3-dehydrogenase; 3-KSR, 3-keto steroid reductase; DHCR24, delta24-sterol reductase; EBP, cholestenol Delta-isomerase; SC5DL, elta7-sterol 5-desaturase; DHCR7, 7-dehydrocholesterol reductase; CYPs, cytochrome P450s.

Peimine and peiminine are considered to be the principal steroidal alkaloids of FTB. Yet, there is limited information about the genes involved in the downstream biosynthesis from cholesterol to steroidal alkaloids in plants at the biochemical level. The biosynthetic pathway from cholesterol into verazine was confirmed to be compromised by a series of C-22 and C-26 hydroxylation/oxidation reactions and C-26 transamination reactions in Veratrum Californicum and Veratrum grandiflorum [38, 39]. Verazine has five circular structures and has formed a basic circular skeleton. It was speculated that it is a precursor of the C-nor D-homo structure formed by the reaction of Peimine and peiminine, plant cytochrome P450 (CYP450) might play prominent roles in the secondary metabolite biosynthesis processes [18]; Based on this knowledge, potential transcripts involved in steroid alkaloid biosynthesis within F. Thunbergii were analyzed (Fig 6, Table 3). Present speculation with different from Shan et al.’s findings [40], which speculated that peimine and peiminine in Fritillaria anhuiensis were biosynthesized by two precursor substances respectively. Considering the only difference in structure is C-5, we inferred that they might have common precursor cholesterol and shared some subsequent conversion reactions.

10.1371/journal.pone.0309978.t003 Table 3 Candidate CYPs associated with steroidal alkaloid biosynthesis.

Function	Unigenes_ID	Reference gene	Species	
C-6 oxidation	Cluster-73431.46232	CYP85A3 [68]	Solanum lycopersicum	
C-6 hydroxylation	Cluster-63258.0	CYP716A53v2 [68]	Panax ginseng	
C-6 hydroxylation	Cluster-61668.0	CYP716E26 [68]	Solanum lycopersicum	
C-11 hydroxylation	Cluster-73431.43439	CYP88D6 [68]	Glycyrrhiza glabra	
C-16 hydroxylation	Cluster-73431.41324	CYP716A141 [64]	Platycodon grandiflorus	
C-16 hydroxylation	Cluster-73431.51706	CYP86A2 [66]	Arabidopsis thaliana	
C-20 hydroxylation	Cluster-73431.28317	CYP71D3 [62]	Lotus japonicus	
C-22 oxidation	Cluster-73431.18721	CYP90G1 [38]	Veratrum californicum	
C-22 hydroxylation	Cluster-73431.13425	CYP90B1 [54, 55]	Arabidopsis thaliana	
C-22 hydroxylation	Cluster-73431.15630	CYP90B2 [56]	Oryza sativa	
C-22 hydroxylation	Cluster-39298.0	CYP724B1 [56]	Oryza sativa	
C-22 hydroxylation	Cluster-73431.27372	CYP72A188 [58]	Solanum tuberosum	
C-26 oxidation	Cluster-73431.15020	GAME4 [62]	Solanum lycopersicum	
C-26 hydroxylation	Cluster-73431.27371	CYP72A208 [58]	Solanum tuberosum	
C-26 hydroxylation	Cluster-73431.27483	CYP734A6 [38, 61]	Arabidopsis thaliana	
C-26 oxidation/hydroxylation	Cluster-73431.14169	CYP94N1 [38, 61]	Veratrum californicum	
C-26 transamination	Cluster-73431.19556	GABAT1 [38]	Veratrum californicum	
C-26 transamination	Cluster-73431.26754	GABAT2 [38]	Solanum lycopersicum	

Quantitative real-time PCR (qRT-PCR) validation

To further confirm the accuracy of the RNA-seq data, qRT-PCR analysis of 12 genes related to the steroidal alkaloid biosynthesis pathway was performed. The relative gene expression levels of PP2C, SnRK2, ABF, JAZ, ACAT, DXS, MVD, FPS, DSDHL, CYP94N1, CYP90G1, and GABAT1 from three FTBs of different fertilization conditions were calculated, which were consistent with RNA-Seq FPKM values (Fig 7, S1 Table).

10.1371/journal.pone.0309978.g007 Fig 7 Comparative expression pattern validation related to the steroidal alkaloid biosynthesis pathway structural genes as obtained from RNA-seq data and qRT-PCR.

Discussion

Discrepancy of alkaloid contents in F. Thunbergii from different fertilization conditions

Research has shown that potassium is essential for plant growth. With sufficient nitrogen and phosphorus, more alkaloids and potassium are allocated to the underground as potassium application increases. The yield and quality of bulbs are positively correlated with the accumulation of underground potassium, so potassium application is an effective strategy to improve the yield and quality of F. Thunbergii [41, 42]. Chicken manure, organic fertilizer, and plant ash fertilizers all contain various elements necessary for plant growth. after HPLC-ELSD detection, The alkaloid content of F. Thunbergii under three cultivation methods in Pan’an County was detected by HPLC-ELSD. The total amount of steroidal alkaloids in FTB in RC was lower than that of other samples, but the yield under RC fertilization conditions was the highest (Table 1). The results in the plant ash group were exactly the opposite. Therefore, it is recommended to choose a cultivation method that combines organic fertilizer and plant ash, providing an experimental basis and theoretical guidance for improving yield and steroidal alkaloid content.

Revelation of plant hormone signal transduction pathway in F. Thunbergii

Plant hormones are rapid response signal compounds related to regulating growth, development in key stages, and resistance to environmental stress responses [29]. According to Fig 4, 123 DEGs are distributed in eight plant hormone signal transduction pathways, among which the abscisic acid (ABA) pathway is positively correlated with the alkaloid metabolism profile of FTB, thus ABA may be the uppermost hormone for regulating the growth of FTB. ABA is the prominent medium of these response processes owing to regulating stomatal closure [43] and seed dormancy [29], resulting in the maintenance of growth and survival of plants. The gene expression of SnRK2 and JAZ in PA was the highest, that SnRK2 (ABA-dependent phosphorylase) interacts with ABA that determines the variation of symplastic calcium ion concentrations, the increase of ion efflux, the reduction of cell expansion caused by water outflow from the guard cells, and finally leading to stomatal closure [43]. ABF induced by ABA and collaboratively regulate ABRE-dependent ABA signaling involved in various stress treatments [44]. ABA was also called “stress hormone”, which enhanced the stress resistance of plant tissue under environmental or abiotic stress such as cold [45] or drought [46].

Additionally, ABA had an antagonistic effect on ET, BR, and JA [47]. In the JA signaling pathway, almost all DEGs showed an upward trend in PA (Fig 4), majority of genes in the JA signaling pathway also demonstrated a similar expression spectrum. JAZ protein is a repressor of jasmonic acid (JA) signal transduction, JAZ has the highest concentration of DEGs in the JA signal transduction pathway. The gene expression level in PA samples is significantly higher than that in samples from other conditions, and the trend of gene expression in the ABA pathway is similar. Furthermore, In the BR signaling pathway, almost all DEGs showed a downward trend in PA, BR-activated transcription factors repressed the expression of desiccation transcription factor 26 (RD26), suggesting that antagonistic crosstalk between ABA signaling and BR signaling contributed to drought stress responses [48]. Crosstalk also existed between cytokinin and ABA, type A Arabidopsis response regulator 5 (ARR5) mediated by SnRK2 is a negative regulator of cytokinin signal transduction, which heightened the stability of its protein, thus lessening cytokinin reaction during drought stress [49]. Further attempts and innovations at stress and hormone action are required to improve the principal active ingredients of F. thunbergii. The effect of hormones on the growth and quality of F. thunbergii is currently not fully studied, and further experimental verification will be designed in the future.

Speculation of plant steroid alkaloids biosynthesis associated genes

Based on the present data, we screened out candidate genes that affected the accumulation of steroidal alkaloids in FTB, which could help improve the quality of F. thunbergii. Firstly, in the biosynthesis of terpenoid backbone, as though both MVA and MEP pathways existed in plants, only the genes in the MEP pathway of F. thunbergii are higher than F. cirrhosae, which indicating the MEP pathway may be the main way to control the biosynthesis of DMAPP/IPP [50]. We found that the DEGs were mainly distributed in MEP pathways, which was consistent with previous conclusions (Fig 5, S8 Table). DXS is the first rate-limiting step in the catalytic MEP pathway, which catalyzes the production of 1-deoxyglucose 5-phosphate (DXP) from pyruvate and 3-phosphate glyceraldehyde [51]. DXR is the second-rate limiting enzyme in the MEP pathway, which can catalyze the synthesis of various important intermediate metabolite enzymes such as IPP, pyridoxine, and thiamine by DXP [52]. Compared with RC samples, the expression level of DXS and DXR in PA and OF exceeded RC, which provided strong evidence for the enzyme to be the key rate-limiting enzyme for the biosynthesis of the terpenoid backbone.

Understanding the medicinal ingredients biosynthesis of FTB is essential to bioengineering and metabolic engineering for sustainable production [38, 53]. So far, we do not know the specific mechanism of the downstream pathway from cholesterol toward steroidal alkaloids in FTB, and the pathway prediction diagram was plotted based on previous literature (Fig 6). A large number of enzymes from other plants took part in the switch from cholesterol to steroidal alkaloids or glycoalkaloids have been reported. The steroid hydroxylation at C-22 position could be catalysed by CYP90B27 in Veratrum californicum [38], CYP90B1 (DWF4) in Arabidopsis thaliana [54, 55], CYP90B2 and CYP724B1 in Oryza sativa (rice) [56], CYP90B3 and CYP724B2 in Lycopersicon esculentum (tomato) [57], and CYP72A188 (PGA2) in Solanum tuberosum (potato) [58]. Next, CYP734A1 in A. thaliana [59], CYP734A7 in L. esculentum [60], and CYP72A208 (PGA1) in S. tuberosum [58] participated in C-26 hydroxylation reactions. In addition, CYP94N1 in V. californicum and CYP734A6 in O. sativa also participated in C-26 oxidation [38, 61]. Afterward, the products of the previous reactions were subject to transamination reaction at C-26 by GABAT1 [38] or GAME12 [62] and oxidation at C-22 by CYP94N1 or CYP90G1 [38], which is prepared for the formation of F-ring and the generation of verazine. The previous research revealed that cevanine, veratramine, and jervin were probably generated by the same C-nor-D-homo-steroidal intermediate [63], but some specific enzymes catalyzed the subsequent reaction of conversion from verazine to C-nor-D-homo-steroidal skeleton is still unknown. The synthesis mechanism may mainly include C-16 hydroxylation, C-22 reduction, wagner-meerwein rearrangement, and C-12 hydroxylation reactions. The pyrrolidine ring-opening reaction and cyclization reaction might be carried out again later [39] (Fig 6). Additionally, the ability to catalyze C-16 hydroxylation was also enlightened for CYP716s, such as PgCYP716A141 [64], BfCYP716AY1 [65], AtCYP86A2 [66] and St16DOX [67]. CYP716A47 in Panax ginseng partake in C-12 hydroxylation [68]. According to the structure of peimine and peiminine, it is possible to proceed with hydroxylation reactions at C-6 and C-20. The expression level of genes is related to the accumulation of alkaloids in F. thunbergii, but the functional function of related genes still needs further verification.

Conclusion

Different fertilization conditions have different effects on the yield and quality of F. thunbergii. This study revealed the potential mechanism of increased steroidal alkaloid content in the bulb yield of Fritillaria after applying different fertilizers. According to the analysis of differences in steroidal alkaloid content and gene expression levels, The PA and OF samples had multiple common genes upregulated in the plant hormone signal transduction pathway, terpenoid skeleton biosynthesis pathway, and cholesterol synthesis pathway, while RC samples showed the opposite trend. Within plant hormone signal transduction pathways, the ABA signal transduction pathway is positively correlated with the steroid alkaloid metabolism spectrum of F. thunbergii, suggesting that ABA might be the most important hormone regulating the synthesis of steroid alkaloids in F. thunbergii. Furthermore, compared with the MVA pathway, the MEP pathway contained a higher number of DEGs from different conditions, which might be a key synthesis pathway in F. thunbergii that affects the accumulation of downstream steroidal alkaloids by controlling DMAPP/IPP biosynthesis. Cholesterol has been characterized as a precursor substance of steroid alkaloids, however, the downstream biosynthesis of cholesterol to steroidal alkaloids in plants remained unclear. Conduct in-depth exploration of key genes in the biosynthesis pathway of steroidal alkaloids, we identified a total of 62 differentially expressed genes, including 9 of MVA pathways and 29 of MEP pathways, and the expression trend of 19 differentially expressed genes in the steroid synthesis pathway are consistent with the accumulation of steroid alkaloids, which may be involved in steroid alkaloid synthesis, laying the foundation for elucidating the biosynthetic pathway of steroid alkaloids in Streptomyces.

Supporting information

S1 Fig PCA principal component analysis diagram.

(TIF)

S2 Fig Functional annotations and classification of the assembled unigenes.

(A) Homologous species distribution unigenes about FTB compared against the Nr database. (B) GO function annotation and classification statistics of the assembled unigenes. The results are summarized in three main categories: BP, biological process; CC, cellular component; MF, molecular function. (C) KEGG function annotation and classification statistics of the assembled unigenes.

(ZIP)

S1 Table Primers for qRT-PCR analysis.

(DOCX)

S2 Table Numbers and quality of RNA-seq reads produced in each sample.

(DOCX)

S3 Table Length distribution of unigenes from the transcriptome data.

(DOCX)

S4 Table Transcripts and FPKM of genes involved in plant hormone signal transduction (ko04075).

(DOCX)

S5 Table Transcripts and FPKM of genes involved in plant-pathogen interaction (ko04626).

(DOCX)

S6 Table Transcripts and FPKM of genes involved in mRNA surveillance pathway (ko03015).

(DOCX)

S7 Table Transcripts and FPKM of genes involved in ubiquitin mediated proteolysis (ko04120).

(DOCX)

S8 Table Transcripts and FPKM of genes involved in terpenoid backbone biosynthesis (ko00900).

(DOCX)

S9 Table Transcripts and FPKM of genes involved in steroid biosynthesis (ko00100).

(DOCX)

10.1371/journal.pone.0309978.r001
Decision Letter 0
Li Minhui Academic Editor
© 2024 Minhui Li
2024
Minhui Li
https://creativecommons.org/licenses/by/4.0/ This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Submission Version0
10 Jun 2024

PONE-D-24-13936Transcriptomic analysis reveals effects of fertilization towards growth and quality of Fritillariae thunbergii  bulbusPLOS ONE

Dear Dr. Zhong,

Thank you for submitting your manuscript to PLOS ONE. After careful consideration, we feel that it has merit but does not fully meet PLOS ONE’s publication criteria as it currently stands. Therefore, we invite you to submit a revised version of the manuscript that addresses the points raised during the review process.

Thank you for your contributions to the field. Suggest making revisions based on the reviewer's comments. All my decisions are made in accordance with PLOS ONE publication standards.

Please submit your revised manuscript by Jul 25 2024 11:59PM. If you will need more time than this to complete your revisions, please reply to this message or contact the journal office at plosone@plos.org. When you're ready to submit your revision, log on to https://www.editorialmanager.com/pone/ and select the 'Submissions Needing Revision' folder to locate your manuscript file.

Please include the following items when submitting your revised manuscript:A rebuttal letter that responds to each point raised by the academic editor and reviewer(s). You should upload this letter as a separate file labeled 'Response to Reviewers'.

A marked-up copy of your manuscript that highlights changes made to the original version. You should upload this as a separate file labeled 'Revised Manuscript with Track Changes'.

An unmarked version of your revised paper without tracked changes. You should upload this as a separate file labeled 'Manuscript'.

If you would like to make changes to your financial disclosure, please include your updated statement in your cover letter. Guidelines for resubmitting your figure files are available below the reviewer comments at the end of this letter.

If applicable, we recommend that you deposit your laboratory protocols in protocols.io to enhance the reproducibility of your results. Protocols.io assigns your protocol its own identifier (DOI) so that it can be cited independently in the future. For instructions see: https://journals.plos.org/plosone/s/submission-guidelines#loc-laboratory-protocols. Additionally, PLOS ONE offers an option for publishing peer-reviewed Lab Protocol articles, which describe protocols hosted on protocols.io. Read more information on sharing protocols at https://plos.org/protocols?utm_medium=editorial-email&utm_source=authorletters&utm_campaign=protocols.

We look forward to receiving your revised manuscript.

Kind regards,

Minhui Li, PhD

Academic Editor

PLOS ONE

Journal requirements:

When submitting your revision, we need you to address these additional requirements.

1. Please ensure that your manuscript meets PLOS ONE's style requirements, including those for file naming. The PLOS ONE style templates can be found at

https://journals.plos.org/plosone/s/file?id=wjVg/PLOSOne_formatting_sample_main_body.pdf and

https://journals.plos.org/plosone/s/file?id=ba62/PLOSOne_formatting_sample_title_authors_affiliations.pdf

2. We note that the grant information you provided in the ‘Funding Information’ and ‘Financial Disclosure’ sections do not match.

When you resubmit, please ensure that you provide the correct grant numbers for the awards you received for your study in the ‘Funding Information’ section.

3. PLOS requires an ORCID iD for the corresponding author in Editorial Manager on papers submitted after December 6th, 2016. Please ensure that you have an ORCID iD and that it is validated in Editorial Manager. To do this, go to ‘Update my Information’ (in the upper left-hand corner of the main menu), and click on the Fetch/Validate link next to the ORCID field. This will take you to the ORCID site and allow you to create a new iD or authenticate a pre-existing iD in Editorial Manager. Please see the following video for instructions on linking an ORCID iD to your Editorial Manager account: https://www.youtube.com/watch?v=_xcclfuvtxQ.

4. Please include captions for your Supporting Information files at the end of your manuscript, and update any in-text citations to match accordingly. Please see our Supporting Information guidelines for more information: http://journals.plos.org/plosone/s/supporting-information.

Additional Editor Comments:

Two reviewers have provided valuable feedback. Please respond to their comments and make appropriate additions to the article.

Additional Comments from PLOS Editorial Office:

We note that one or more reviewers has recommended that you cite specific previously published works. As always, we recommend that you please review and evaluate the requested works to determine whether they are relevant and should be cited. It is not a requirement to cite these works. We appreciate your attention to this request.

[Note: HTML markup is below. Please do not edit.]

Reviewers' comments:

Reviewer's Responses to Questions

Comments to the Author

1. Is the manuscript technically sound, and do the data support the conclusions?

The manuscript must describe a technically sound piece of scientific research with data that supports the conclusions. Experiments must have been conducted rigorously, with appropriate controls, replication, and sample sizes. The conclusions must be drawn appropriately based on the data presented.

Reviewer #1: Partly

Reviewer #2: Yes

**********

2. Has the statistical analysis been performed appropriately and rigorously?

Reviewer #1: I Don't Know

Reviewer #2: Yes

**********

3. Have the authors made all data underlying the findings in their manuscript fully available?

The PLOS Data policy requires authors to make all data underlying the findings described in their manuscript fully available without restriction, with rare exception (please refer to the Data Availability Statement in the manuscript PDF file). The data should be provided as part of the manuscript or its supporting information, or deposited to a public repository. For example, in addition to summary statistics, the data points behind means, medians and variance measures should be available. If there are restrictions on publicly sharing data—e.g. participant privacy or use of data from a third party—those must be specified.

Reviewer #1: Yes

Reviewer #2: Yes

**********

4. Is the manuscript presented in an intelligible fashion and written in standard English?

PLOS ONE does not copyedit accepted manuscripts, so the language in submitted articles must be clear, correct, and unambiguous. Any typographical or grammatical errors should be corrected at revision, so please note any specific errors here.

Reviewer #1: Yes

Reviewer #2: Yes

**********

5. Review Comments to the Author

Please use the space provided to explain your answers to the questions above. You may also include additional comments for the author, including concerns about dual publication, research ethics, or publication ethics. (Please upload your review as an attachment if it exceeds 20,000 characters)

Reviewer #1: The Abstract is lacking key results like the yield and alkaloids content obtained from the best fertilizer.

In the introduction, is is preferable to state the hypothesis tested rather than the objectives.

In the methods section, the authors should justify why they use only one cultiver rather than using many. the quantity and justification of the quantity of eah of the fertiliser used for the study. what was the growing conditions of the plants. Alkaloids content was evaluated of plant tissue of how old? You should use AOVA test with mean separation as a statistical means. Also provide The software used for clustering transcriptome data

Reviewer #2: The manuscript entitled "Transcriptomic analysis reveals effects of fertilization towards growth and quality of Fritillariae thunbergii bulbus" provides valuable insights into optimizing the cultivation conditions of Fritillariae thunbergii (FTB), a traditional Chinese medicine renowned for its medicinal properties. The study investigates the impact of different fertilizers on the yield and quality of FTB, focusing on raw chicken manure (RC), organic fertilizer (OF), and plant ash (PA). Through a combination of HPLC-ELSD detection and yield statistics, OF application emerged as the most effective fertilizer. Transcriptome analysis uncovered the up-regulation of the ABA signaling pathway, potentially enhancing bulb yield, along with the identification of key genes associated with steroidal alkaloid accumulation. However, the manuscript requires minor revisions. Specifically, replication information regarding biological and technical replicates should be provided, and transcription factor/gene names should be italicized throughout the manuscript. Additionally, the methodology section should include NCBI SRA submission information. Furthermore, the introduction and discussion sections would benefit from the inclusion of references such as (Kumar, V., Sharma, S., Kumar, P. (2024). Insight into the Genetics and Genomics Studies of the Fritillaria Species. In: Gahlaut, V., Jaiswal, V. (eds) Genetics and Genomics of High-Altitude Crops. Springer, Singapore. https://doi.org/10.1007/978-981-99-9175-4_4; https://doi.org/10.1016/j.indcrop.2023.117541; Kumar, V., Kumar, P., Bhargava, B. et al. Transcriptomic and Metabolomic Reprogramming to Explore the High-Altitude Adaptation of Medicinal Plants: A Review. J Plant Growth Regul 42, 7315–7329 (2023). https://doi.org/10.1007/s00344-023-11018-8; Kapoor, B., Kumar, A. & Kumar, P. Transcriptome repository of North-Western Himalayan endangered medicinal herbs: a paramount approach illuminating molecular perspective of phytoactive molecules and secondary metabolism. Mol Genet Genomics 296, 1177–1202 (2021). https://doi.org/10.1007/s00438-021-01821-x) to enrich the context and support the findings. Overall, these adjustments will enhance the clarity and completeness of the manuscript, contributing to its scientific significance.

**********

6. PLOS authors have the option to publish the peer review history of their article (what does this mean?). If published, this will include your full peer review and any attached files.

If you choose “no”, your identity will remain anonymous but your review may still be made public.

Do you want your identity to be public for this peer review? For information about this choice, including consent withdrawal, please see our Privacy Policy.

Reviewer #1: Yes: Eric Bertrand Kouam

Reviewer #2: No

**********

[NOTE: If reviewer comments were submitted as an attachment file, they will be attached to this email and accessible via the submission site. Please log into your account, locate the manuscript record, and check for the action link "View Attachments". If this link does not appear, there are no attachment files.]

While revising your submission, please upload your figure files to the Preflight Analysis and Conversion Engine (PACE) digital diagnostic tool, https://pacev2.apexcovantage.com/. PACE helps ensure that figures meet PLOS requirements. To use PACE, you must first register as a user. Registration is free. Then, login and navigate to the UPLOAD tab, where you will find detailed instructions on how to use the tool. If you encounter any issues or have any questions when using PACE, please email PLOS at figures@plos.org. Please note that Supporting Information files do not need this step.

Attachment Submitted filename: PONE-D-24-13936. review Kouam.pdf

10.1371/journal.pone.0309978.r002
Author response to Decision Letter 0
Submission Version1
11 Jul 2024

Dear Editor and Reviewers:

Thanks for your letter and the reviewers' comments concerning our manuscript, "Transcriptomic analysis reveals effects of fertilization towards growth and quality of Fritillariae thunbergii bulbus" (Manuscript ID: PONE-D-24-13936). Those comments are all valuable and helpful for revising and improving our paper. We have studied all comments carefully and have made conscientious corrections. Our response is given in standard font, and changes to the manuscript are given in red text. The responses to the reviewers' comments are as follows:

Response to the reviewers:

Reviewer 1

We feel great thanks for your professional review work on our article. As you are concerned, there are several problems that need to be addressed. According to your nice suggestions, we have made extensive corrections to our previous draft, the detailed corrections are listed below.

1. The Abstract is lacking key results like the yield and alkaloids content obtained from the best fertilizer.

Answer: Line 25-27, the yield and alkaloids content of “The combined results of HPLC-ELSD detection and yield statistics showed that the F. thunbergii with OF application was the best, with the content of peimine and peiminine reaching 0.0603% and 0.0502%, respectively. In addition, the yield was 2.70 kg/m². ” was added.

2. In the introduction, it is preferable to state the hypothesis tested rather than the objectives.

Answer: Line 71-77, and Line 82-84, based on existing literatures on the correlation between altitude and secondary metabolic yield of medicinal plants, and the application to in vitro cultivation of Himalayan Fritillaria species, it is speculated that different fertilization methods may be related to the alkaloid content and yield of Fritillaria.

3. In the methods section, the authors should justify why they use only one cultiver rather than using many. the quantity and justification of the quantity of each of the fertiliser used for the study. what was the growing conditions of the plants. Alkaloids content was evaluated of plant tissue of how old? You should use AOVA test with mean separation as a statistical means. Also provide The software used for clustering transcriptome data.

Answer: To control variables, we selected the same batch of Fritillariae thunbergii bulbs for different fertilization conditions to prevent the influence of different varieties on their alkaloid contents. Among them, ‘Zhebei 1’ has a wide range of applications, which is more representative and practical.

The growth environment of Fritillariae thunbergii bulbs in the three experimental fields was the same natural environment, with the same altitude 600 m.

Three year old bulbs of Fritillariae thunbergii were used in this experiment.

Line 173-174, and Line 194-196, we conducted One-way ANOVA of alkaloid content of different fertilization conditions using SPSS (25.0.0).

Line 161-162, “Furthermore, the clustered expression pattern was drawn through R (Version 3.0.3), Package ggplot2 and Package pheatmap.” were added.

Additionally, Line 184-193, we have provided a more detailed description of the experimental results.

4.There is no standard deviation?- yield

Answer: Thank you for your reminder. Line 93-94, Fritillaria thunbergii bulb yield was determined by harvesting the plants in an area of 2 m2 in each plot. We collected bulbs from a whole experimental field and calculated the yield.

Reviewer 2

The manuscript entitled "Transcriptomic analysis reveals effects of fertilization towards growth and quality of Fritillariae thunbergii bulbus" provides valuable insights into optimizing the cultivation conditions of Fritillariae thunbergii (FTB), a traditional Chinese medicine renowned for its medicinal properties. The study investigates the impact of different fertilizers on the yield and quality of FTB, focusing on raw chicken manure (RC), organic fertilizer (OF), and plant ash (PA). Through a combination of HPLC-ELSD detection and yield statistics, OF application emerged as the most effective fertilizer. Transcriptome analysis uncovered the up-regulation of the ABA signaling pathway, potentially enhancing bulb yield, along with the identification of key genes associated with steroidal alkaloid accumulation. However, the manuscript requires minor revisions.

Answer: Thank you for your valuable comments, revisions have been made to improve the readability of the present manuscript. Detailed changes and corresponding answers have been displayed below.

1. Replication information regarding biological and technical replicates should be provided.

Answer: Line 95-96, “ all the bulbs were collected biological replicates to randomly select different regions in the north, south, east, and west of the same production area” was added. Line 99-100, and Line 112-113, the number biological replicates has been explained.

2. Transcription factor/gene names should be italicized throughout the manuscript.

Answer: We were really sorry for our careless mistakes. Thank you for noticing. Transcription factor/gene names were corrected to Italic format.

3. The methodology section should include NCBI SRA submission information.

Answer: Line 175-177, “Data avaliability: The raw data of RNA-Seq have been submitted to the the NCBI Sequence Read Archive (SRA) database (Accession: PRJNA950562).” was added.

4. The introduction and discussion sections would benefit from the inclusion of references such as (Kumar, V., Sharma, S., Kumar, P. (2024). Insight into the Genetics and Genomics Studies of the Fritillaria Species. In: Gahlaut, V., Jaiswal, V. (eds) Genetics and Genomics of High-Altitude Crops. Springer, Singapore. https://doi.org/10.1007/978-981-99-9175-4_4;https://doi.org/10.1016/j.indcrop.2023.117541; Kumar, V., Kumar, P., Bhargava, B. et al. Transcriptomic and Metabolomic Reprogramming to Explore the High-Altitude Adaptation of Medicinal Plants: A Review. J Plant Growth Regul 42, 7315–7329 (2023). https://doi.org/10.1007/s00344-023-11018-8; Kapoor, B., Kumar, A. & Kumar, P. Transcriptome repository of North-Western Himalayan endangered medicinal herbs: a paramount approach illuminating molecular perspective of phytoactive molecules and secondary metabolism. Mol Genet Genomics 296, 1177–1202 (2021). https://doi.org/10.1007/s00438-021-01821-x) to enrich the context and support the findings. Overall, these adjustments will enhance the clarity and completeness of the manuscript, contributing to its scientific significance.

Answer: We sincerely appreciate the valuable comments. We have checked the literature carefully and added the four references on 11 (Line 58-60), 15 (Line 71-72), 16 (Line 71-72) and 17 (Line 72-73) into the Introduction part in the revised manuscript.

Thank you very much for your attention and time. Look forward to hearing from you.

Yours sincerely,

Corresponding author:

Name: Zhuoheng Zhong

E-mail: zhongzhh@zstu.edu.cn

Attachment Submitted filename: Response.docx

10.1371/journal.pone.0309978.r003
Decision Letter 1
Li Minhui Academic Editor
© 2024 Minhui Li
2024
Minhui Li
https://creativecommons.org/licenses/by/4.0/ This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Submission Version1
22 Aug 2024

Transcriptomic analysis reveals effects of fertilization towards growth and quality of Fritillariae thunbergii  bulbus

PONE-D-24-13936R1

Dear Dr. Fu,

We’re pleased to inform you that your manuscript has been judged scientifically suitable for publication and will be formally accepted for publication once it meets all outstanding technical requirements.

Within one week, you’ll receive an e-mail detailing the required amendments. When these have been addressed, you’ll receive a formal acceptance letter and your manuscript will be scheduled for publication.

An invoice will be generated when your article is formally accepted. Please note, if your institution has a publishing partnership with PLOS and your article meets the relevant criteria, all or part of your publication costs will be covered. Please make sure your user information is up-to-date by logging into Editorial Manager at Editorial Manager® and clicking the ‘Update My Information' link at the top of the page. If you have any questions relating to publication charges, please contact our Author Billing department directly at authorbilling@plos.org.

If your institution or institutions have a press office, please notify them about your upcoming paper to help maximize its impact. If they’ll be preparing press materials, please inform our press team as soon as possible -- no later than 48 hours after receiving the formal acceptance. Your manuscript will remain under strict press embargo until 2 pm Eastern Time on the date of publication. For more information, please contact onepress@plos.org.

Kind regards,

Minhui Li, PhD

Academic Editor

PLOS ONE

Additional Editor Comments (optional):

Thank you for your contributions to the field.

Reviewers' comments:

Reviewer's Responses to Questions

Comments to the Author

1. If the authors have adequately addressed your comments raised in a previous round of review and you feel that this manuscript is now acceptable for publication, you may indicate that here to bypass the “Comments to the Author” section, enter your conflict of interest statement in the “Confidential to Editor” section, and submit your "Accept" recommendation.

Reviewer #1: All comments have been addressed

Reviewer #2: All comments have been addressed

Reviewer #3: All comments have been addressed

**********

2. Is the manuscript technically sound, and do the data support the conclusions?

The manuscript must describe a technically sound piece of scientific research with data that supports the conclusions. Experiments must have been conducted rigorously, with appropriate controls, replication, and sample sizes. The conclusions must be drawn appropriately based on the data presented.

Reviewer #1: Yes

Reviewer #2: Yes

Reviewer #3: Yes

**********

3. Has the statistical analysis been performed appropriately and rigorously?

Reviewer #1: Yes

Reviewer #2: Yes

Reviewer #3: Yes

**********

4. Have the authors made all data underlying the findings in their manuscript fully available?

The PLOS Data policy requires authors to make all data underlying the findings described in their manuscript fully available without restriction, with rare exception (please refer to the Data Availability Statement in the manuscript PDF file). The data should be provided as part of the manuscript or its supporting information, or deposited to a public repository. For example, in addition to summary statistics, the data points behind means, medians and variance measures should be available. If there are restrictions on publicly sharing data—e.g. participant privacy or use of data from a third party—those must be specified.

Reviewer #1: Yes

Reviewer #2: Yes

Reviewer #3: Yes

**********

5. Is the manuscript presented in an intelligible fashion and written in standard English?

PLOS ONE does not copyedit accepted manuscripts, so the language in submitted articles must be clear, correct, and unambiguous. Any typographical or grammatical errors should be corrected at revision, so please note any specific errors here.

Reviewer #1: Yes

Reviewer #2: Yes

Reviewer #3: Yes

**********

6. Review Comments to the Author

Please use the space provided to explain your answers to the questions above. You may also include additional comments for the author, including concerns about dual publication, research ethics, or publication ethics. (Please upload your review as an attachment if it exceeds 20,000 characters)

Reviewer #1: Susbtancial improvement has been carried out in the revised manuscript. I am satisfied with the answers provided to the different questions that was raised.

Reviewer #2: The authors have carefully addressed the comments and duly incorporated in the revised manuscript. Hence manuscript can be accepted in its current form.

Reviewer #3: The author has revised the paper as requested and responded to the questions. I agree to publish this manuscript in Plos one

**********

7. PLOS authors have the option to publish the peer review history of their article (what does this mean?). If published, this will include your full peer review and any attached files.

If you choose “no”, your identity will remain anonymous but your review may still be made public.

Do you want your identity to be public for this peer review? For information about this choice, including consent withdrawal, please see our Privacy Policy.

Reviewer #1: Yes: Eric Bertrand Kouam

Reviewer #2: No

Reviewer #3: Yes: Jingmao You

**********

10.1371/journal.pone.0309978.r004
Acceptance letter
Li Minhui Academic Editor
© 2024 Minhui Li
2024
Minhui Li
https://creativecommons.org/licenses/by/4.0/ This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
10 Sep 2024

PONE-D-24-13936R1

PLOS ONE

Dear Dr. Fu,

I'm pleased to inform you that your manuscript has been deemed suitable for publication in PLOS ONE. Congratulations! Your manuscript is now being handed over to our production team.

At this stage, our production department will prepare your paper for publication. This includes ensuring the following:

* All references, tables, and figures are properly cited

* All relevant supporting information is included in the manuscript submission,

* There are no issues that prevent the paper from being properly typeset

If revisions are needed, the production department will contact you directly to resolve them. If no revisions are needed, you will receive an email when the publication date has been set. At this time, we do not offer pre-publication proofs to authors during production of the accepted work. Please keep in mind that we are working through a large volume of accepted articles, so please give us a few weeks to review your paper and let you know the next and final steps.

Lastly, if your institution or institutions have a press office, please let them know about your upcoming paper now to help maximize its impact. If they'll be preparing press materials, please inform our press team within the next 48 hours. Your manuscript will remain under strict press embargo until 2 pm Eastern Time on the date of publication. For more information, please contact onepress@plos.org.

If we can help with anything else, please email us at customercare@plos.org.

Thank you for submitting your work to PLOS ONE and supporting open access.

Kind regards,

PLOS ONE Editorial Office Staff

on behalf of

Dr. Minhui Li

Academic Editor

PLOS ONE
==== Refs
References

1 Nile SH , Su J , Wu D , Wang L , Hu J , Sieniawska E , et al . Fritillaria thunbergii Miq. (Zhe Beimu): A review on its traditional uses, phytochemical profile and pharmacological properties. Food and Chemical Toxicology. 2021; 153 . doi: 10.1016/j.fct.2021.112289 34029669
2 Hao DC , Gu XJ , Xiao PG , Peng Y . Phytochemical and biological research of Fritillaria medicine resources. Chinese Journal of Natural Medicines. 2013; 11 (4 ): 330–344. doi: 10.1016/s1875-5364(13)60050-3 23845541
3 Chinese Pharmacopoeia Commission. Pharmacopoeia of People’s Republic of China: One Edition. Beijing: China Medical Science Press; 2020: 126–127.
4 Wu F , Tian M , Sun Y , Wu C , Liu X . Efficacy, chemical composition, and pharmacological effects of herbal drugs derived from Fritillaria cirrhosa D. Don and Fritillaria thunbergii Miq. Front Pharmacol. 2022; 13 : 985935. doi: 10.3389/fphar.2022.985935 36532788
5 Lin BQ , Ji H , Li P , Jiang Y , Fang W . Selective antagonism activity of alkaloids from bulbs Fritillariae at muscarinic receptors: functional studies. European journal of pharmacology. 2006; 551 (1–3 ): 125–130. doi: 10.1016/j.ejphar.2006.08.066 17027747
6 Li ZH , An C , Hu KW , Zhou KH , Duan HH , Tang MK . Multidrug resistance reversal activity of total alkaloid from Fritillaria thunbergii on cisplatin-resistant human lung adenocarcinoma A549/DDP cells. Chinese Journal of Pharmacology and Toxicology. 2013; 27 (03 ): 315–320.
7 Li H , Hung A , Li M , Yang A . Fritillariae Thunbergii Bulbus: Traditional Uses, Phytochemistry, Pharmacodynamics, Pharmacokinetics and Toxicity. International Journal of Molecular Sciences. 2019; 20 (7 ). doi: 10.3390/ijms20071667 30987173
8 Suh WS , Lee SY , Park JE , Kim DH , Kim SY , Lee KR . Two new steroidal alkaloids from the bulbs of Fritillaria thunbergii. Heterocycles. 2018; 96 (5 ): 925–934. doi: 10.3987/com-18-13891
9 Wu QL , Ji YY , Dong LL , Shen XX , Wang ZA , Wang ZH . Study on the alkaloid content and synthesis related gene expression of Fritillaria thunbergii from different regions. Guangxi Botanical Journal. 2020; 40 (12 ): 1755–1763.
10 Lu X , Bi YM , Huo YZ , Yi SF , Jiao XL , Zhang J , et al . The effects of different fertilization and coverage measures on the yield and effective components of Fritillaria thunbergii. Modern Chinese Medicine, 2018; 20 (05 ): 576–580.
11 Kapoor B , Kumar A , Kumar P . Transcriptome repository of North-Western Himalayan endangered medicinal herbs: a paramount approach illuminating molecular perspective of phytoactive molecules and secondary metabolism. Molecular Genetics And Genomics. 2021; 296 (6 ): 1177–1202. doi: 10.1007/s00438-021-01821-x 34557965
12 Schaller H. New aspects of sterol biosynthesis in growth and development of higher plants. Plant Physiology and Biochemistry. 2004; 42 (6 ): 465–476. doi: 10.1016/j.plaphy.2004.05.012 15246059
13 Liu M , Li LN , Pan YT , Kong JQ . cDNA isolation and functional characterization of squalene synthase gene from Ornithogalum caudatum. Protein Expression and Purification. 2017; 130 : 63–72. doi: 10.1016/j.pep.2016.10.002 27725246
14 Cheng Y , Liu H , Tong X , Liu Z , Zhang X , Li D , et al . Identification and analysis of CYP450 and UGT supergene family members from the transcriptome of Aralia elata (Miq.) seem reveal candidate genes for triterpenoid saponin biosynthesis. BMC Plant Biology. 2020; 20 (1 ): 214. doi: 10.1186/s12870-020-02411-6 32404131
15 Vijay G , Vandana J . Genetics and Genomics of High-Altitude Crops. 1st ed. Singapore: Springer Nature Singapore; 2024.
16 Kumar V , Kumar P , Bhargava B , Sharma R , Irfan M , Chandora R . Transcriptomic and metabolomic reprogramming to explore the High-Altitude adaptation of medicinal plants: A Review. Journal of Plant Growth Regulation. 2023; 42 (12 ): 7315–29. doi: 10.1007/s00344-023-11018-8
17 Sharma S , Kumar P , Sharma R , Warghat AR . In vitro propagation and omics breakthroughs for understanding specialized metabolite production in high-value Himalayan Fritillaria species. Industrial Crops and Products. 2023; 205 : 117541. doi: 10.1016/j.indcrop.2023.117541
18 Hu JD , Mao G , Zhang ZW , Ma CD , Liang ZS , Xia GD , et al . Effect of different fertilization treatments on yield and secondary metabolites of Codonopsis pilosula. Zhongguo Zhong Yao Za Zhi. 2017; 42 (15 ): 2946–2953.29139262
19 Liang Q , Chen XF , Li Y , Zhang J , Meng J , Peng SM . Effects of Chemical Fertilizers and Organic Fertilizer on Yield of Ligusticum chuanxiong Rhizome. Zhong Yao Cai. 2015; 38 (10 ): 2015–2020.27254909
20 Thepsilvisut O , Chutimanukul P , Sae-Tan S , Ehara H . Effect of chicken manure and chemical fertilizer on the yield and qualities of white mugwort at dissimilar harvesting times. PLoS One. 2022; 17 (4 ): e0266190. doi: 10.1371/journal.pone.0266190 35472063
21 Davidson NM , Oshlack A . Corset: enabling differential gene expression analysis for de novo assembled transcriptomes. Genome Biology. 2014; 15 (7 ): 410. doi: 10.1186/s13059-014-0410-6 25063469
22 Waterhouse RM , Seppey M , Simão FA , Manni M , Ioannidis P , Klioutchnikov G , et al . BUSCO Applications from Quality Assessments to Gene Prediction and Phylogenomics. Molecular Biology and Evolution. 2018; 35 (3 ): 543–548. doi: 10.1093/molbev/msx319 29220515
23 Li B , Dewey CN . RSEM: accurate transcript quantification from RNA-Seq data with or without a reference genome. BMC Bioinformatics. 2011; 12 : 323. doi: 10.1186/1471-2105-12-323 21816040
24 Mortazavi A , Williams BA , McCue K , Schaeffer L , Wold B . Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nature Methods. 2008; 5 (7 ): 621–628. doi: 10.1038/nmeth.1226 18516045
25 Kanehisa M , Goto S . KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Research. 2000; 28 (1 ): 27–30. doi: 10.1093/nar/28.1.27 10592173
26 Wang L , Feng Z , Wang X , Wang X , Zhang X . DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics. 2010; 26 (1 ): 136–138. doi: 10.1093/bioinformatics/btp612 19855105
27 Livak KJ , Schmittgen TD . Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) Method. Methods. 2001; 25 (4 ): 402–408. doi: 10.1006/meth.2001.1262 11846609
28 Sun Y , Fan XY , Cao DM , Tang W , He K , Zhu JY , et al . Integration of brassinosteroid signal transduction with the transcription network for plant growth regulation in Arabidopsis. Developmental Cell. 2010; 19 (5 ): 765–777. doi: 10.1016/j.devcel.2010.10.010 21074725
29 Verma V , Ravindran P , Kumar PP . Plant hormone-mediated regulation of stress responses. BMC Plant Biology. 2016; 16 : 86. doi: 10.1186/s12870-016-0771-y 27079791
30 Muthamilarasan M , Prasad M . Plant innate immunity: An updated insight into defense mechanism. Journal of Biosciences. 2013; 38 (2 ): 433–449. doi: 10.1007/s12038-013-9302-2 23660678
31 Wagner E , Lykke-Andersen J . mRNA surveillance: the perfect persist. Journal of cell science. 2002; 115 (Pt 15 ): 3033–3038. doi: 10.1242/jcs.115.15.3033 12118059
32 Liu X , Wang Q , Chen W , Wang C . Dynamic regulation of innate immunity by ubiquitin and ubiquitin-like proteins. Cytokine & Growth Factor Reviews. 2013; 24 (6 ): 559–570. doi: 10.1016/j.cytogfr.2013.07.002 23953672
33 Shaid S , Brandts CH , Serve H , Dikic I . Ubiquitination and selective autophagy. Cell Death and Differentiation. 2013; 20 (1 ): 21–30. doi: 10.1038/cdd.2012.72 22722335
34 Sonawane PD , Pollier J , Panda S , Szymanski J , Massalha H , Yona M , et al . Plant cholesterol biosynthetic pathway overlaps with phytosterol metabolism. Nature Plants. 2016; 3 : 16205. doi: 10.1038/nplants.2016.205 28005066
35 Houseley J , Tollervey D . The many pathways of RNA degradation. Cell. 2009; 136 (4 ):763–776. doi: 10.1016/j.cell.2009.01.019 19239894
36 Choi H , Hong J , Ha J , Kang J , Kim SY . ABFs, a family of ABA-responsive element binding factors. Journal of Biological Chemistry. 2000; 275 (3 ): 1723–1730. doi: 10.1074/jbc.275.3.1723 10636868
37 Sawai S , Ohyama K , Yasumoto S , Seki H , Sakuma T , Yamamoto T , et al . Sterol side chain reductase 2 is a key enzyme in the biosynthesis of cholesterol, the common precursor of toxic steroidal glycoalkaloids in potato. Plant Cell. 2014; 26 (9 ): 3763–3774. doi: 10.1105/tpc.114.130096 25217510
38 Augustin MM , Ruzicka DR , Shukla AK , Augustin JM , Starks CM , O’Neil-Johnson M , et al . Elucidating steroid alkaloid biosynthesis in Veratrum californicum: production of verazine in Sf9 cells. The Plant Journal. 2015; 82 (6 ): 991–1003. doi: 10.1111/tpj.12871 25939370
39 Heretsch P , Tzagkaroulaki L , Giannis A . Cyclopamine and hedgehog signaling: chemistry, biology, medical perspectives. Angewandte Chemie International Edition. 2010; 49 (20 ): 3418–3427. doi: 10.1002/anie.200906967 20429080
40 Shan T , Shi Y , Xu J , Zhao L , Tao Y , Wu J . Transcriptome analysis reveals candidate genes related to steroid alkaloid biosynthesis in Fritillaria anhuiensis. Physiologia Plantarum. 2022; 174 (4 ): e13755.
41 Sui N , Wang L , Sun J , Wu J , Chen M , Zhang S , et al . Estimating Fritillaria thunbergii Miq. yield, quality, and potassium use efficiency in response to potassium application rate. Industrial Crops and Products. 2021; 164 : 113409. doi: 10.1016/j.indcrop.2021.113409
42 Wang L , Sun J , Kai G , Sui N . Screening of the optimal potassium fertilization rate for maximizing bulb yield and potassium nutrient evaluation in Fritillaria thunbergii miq. Clinical Complementary Medicine and Pharmacology. 2021; 1 (1 ): 100003. doi: 10.1016/j.ccmp.2021.100003
43 de Ollas C , Dodd IC . Physiological impacts of ABA-JA interactions under water-limitation. Plant Molecular Biology. 2016; 91 (6 ): 641–650. doi: 10.1007/s11103-016-0503-6 27299601
44 Yoshida T , Fujita Y , Sayama H , Kidokoro S , Maruyama K , Mizoi J , et al . AREB1, AREB2, and ABF3 are master transcription factors that cooperatively regulate ABRE-dependent ABA signaling involved in drought stress tolerance and require ABA for full activation. The Plant Journal. 2010; 61 (4 ): 672–685. doi: 10.1111/j.1365-313X.2009.04092.x 19947981
45 Li XY , Yang Y , Zhang LF , Zuo SY , Li LJ , Jiao J , et al . Regulation on contents of endogenous hormones and asr1 gene expression of Maize seedling by exogenous ABA under low-temperature stress. Acta Agronomica Sinica. 2017; 43 (01 ): 141–148.
46 Urano K , Maruyama K , Jikumaru Y , Kamiya Y , Yamaguchi-Shinozaki K , Shinozaki K . Analysis of plant hormone profiles in response to moderate dehydration stress. The Plant Journal. 2017; 90 (1 ): 17–36. doi: 10.1111/tpj.13460 27995695
47 Anderson JP , Badruzsaufari E , Schenk PM , Manners JM , Desmond OJ , Ehlert C , et al . Antagonistic interaction between abscisic acid and jasmonate-ethylene signaling pathways modulates defense gene expression and disease resistance in Arabidopsis. Plant Cell. 2004; 16 (12 ): 3460–3479. doi: 10.1105/tpc.104.025833 15548743
48 Zhao H , Yin CC , Ma B , Chen SY , Zhang JS . Ethylene signaling in rice and Arabidopsis: New regulators and mechanisms. Journal of Integrative Plant Biology. 2021; 63 (1 ): 102–125. doi: 10.1111/jipb.13028 33095478
49 Huang X , Hou L , Meng J , You H , Li Z , Gong Z , et al . The antagonistic action of abscisic acid and cytokinin signaling mediates drought stress response in Arabidopsis. Molecular Plant. 2018; 11 (7 ): 970–982. doi: 10.1016/j.molp.2018.05.001 29753021
50 Lu Q , Li R , Liao J , Hu Y , Gao Y , Wang M , et al . Integrative analysis of the steroidal alkaloids distribution and biosynthesis of bulbs Fritillariae Cirrhosae through metabolome and transcriptome analyses. BMC Genomics. 2022; 23 (1 ): 511. doi: 10.1186/s12864-022-08724-0 35836113
51 Estévez JM , Cantero A , Romero C , Kawaide H , Jiménez LF , Kuzuyama T , et al . Analysis of the expression of CLA1, a gene that encodes the 1-deoxyxylulose 5-phosphate synthase of the 2-C-methyl-D-erythritol-4-phosphate pathway in Arabidopsis. Plant Physiology. 2000; 124 (1 ): 95–104. doi: 10.1104/pp.124.1.95 10982425
52 Carretero-Paulet L , Ahumada I , Cunillera N , Rodríguez-Concepción M , Ferrer A , Boronat A , et al . Expression and molecular analysis of the Arabidopsis DXR gene encoding 1-deoxy-D-xylulose 5-phosphate reductoisomerase, the first committed enzyme of the 2-C-methyl-D-erythritol 4-phosphate pathway. Plant Physiology. 2002; 129 (4 ): 1581–1591. doi: 10.1104/pp.003798 12177470
53 Kumar P , Ashrita , Acharya V , Warghat AR . Comparative transcriptome analysis infers bulb derived in vitro cultures as a promising source for sipeimine biosynthesis in Fritillaria cirrhosa D. Don (Liliaceae, syn. Fritillaria roylei Hook.)—High value Himalayan medicinal herb. Phytochemistry. 2021; 183 : 112631. doi: 10.1016/j.phytochem.2020.112631 33370713
54 Choe S , Dilkes BP , Fujioka S , Takatsuto S , Sakurai A , Feldmann KA . The DWF4 gene of Arabidopsis encodes a cytochromecpha-hydroxylation steps in brassinosteroid biosynthesis. Plant Cell. 1998; 10 (2 ): 231–243. doi: 10.1105/tpc.10.2.231 9490746
55 Fujita S , Ohnishi T , Watanabe B , Yokota T , Takatsuto S , Fujioka S , et al . Arabidopsis CYP90B1 catalyses the early C-22 hydroxylation of C27, C28 and C29 sterols. The Plant Journal. 2006; 45 (5 ): 765–774. doi: 10.1111/j.1365-313X.2005.02639.x 16460510
56 Sakamoto T , Morinaka Y , Ohnishi T , Sunohara H , Fujioka S , Ueguchi-Tanaka M , et al . Erect leaves caused by brassinosteroid deficiency increase biomass production and grain yield in rice. Nature Biotechnology. 2006; 24 (1 ): 105–109. doi: 10.1038/nbt1173 16369540
57 Ohnishi T , Watanabe B , Sakata K , Mizutani M . CYP724B2 and CYP90B3 function in the early C-22 hydroxylation steps of brassinosteroid biosynthetic pathway in tomato. Bioscience, Biotechnology and Biochemistry. 2006; 70 (9 ): 2071–2080. doi: 10.1271/bbb.60034 16960392
58 Umemoto N , Nakayasu M , Ohyama K , Yotsu-Yamashita M , Mizutani M , Seki H , et al . Two cytochrome P450 monooxygenases catalyze early hydroxylation steps in the potato steroid glycoalkaloid biosynthetic pathway. Plant Physiology. 2016; 171 (4 ): 2458–2467. doi: 10.1104/pp.16.00137 27307258
59 Thornton LE , Rupasinghe SG , Peng H , Schuler MA , Neff MM . Arabidopsis CYP72C1 is an atypical cytochrome P450 that inactivates brassinosteroids. Plant Molecular Biology. 2010; 74 (1–2 ): 167–181. doi: 10.1007/s11103-010-9663-y 20669042
60 Ohnishi T , Nomura T , Watanabe B , Ohta D , Yokota T , Miyagawa H , et al . Tomato cytochrome P450 CYP734A7 functions in brassinosteroid catabolism. Phytochemistry. 2006; 67 (17 ): 1895–1906. doi: 10.1016/j.phytochem.2006.05.042 16872648
61 Sakamoto T , Kawabe A , Tokida-Segawa A , Shimizu B , Takatsuto S , Shimada Y , et al . Rice CYP734As function as multisubstrate and multifunctional enzymes in brassinosteroid catabolism. The Plant Journal. 2011; 67 (1 ): 1–12. doi: 10.1111/j.1365-313X.2011.04567.x 21418356
62 Itkin M , Heinig U , Tzfadia O , Bhide AJ , Shinde B , Cardenas PD , et al . Biosynthesis of antinutritional alkaloids in solanaceous crops is mediated by clustered genes. Science. 2013; 341 (6142 ):175–9. doi: 10.1126/science.1240230 23788733
63 Kaneko K , Mitsuhashi H , Hirayama K , Ohmori S . 11-Deoxojervine as a precursor for jervine biosynthesis in Veratrum grandiflorum. Phytochemistry. 1970; 9 : 2497–501. doi: 10.1016/S0031-9422(00)85769-6
64 Tamura K , Teranishi Y , Ueda S , Suzuki H , Kawano N , Yoshimatsu K , et al . Cytochrome P450 monooxygenase CYP716A141 is a unique b-Amyrin C-16b oxidase involved in triterpenoid saponin biosynthesis in Platycodon grandiflorus (vol 57, pg 874, 2017). Plant and Cell Physiology. 2017; 58 (6 ): 1119–1119. doi: 10.1093/pcp/pcx067 28575525
65 Moses T , Pollier J , Almagro L , Buyst D , Van Montagu M , Pedreño MÁ , et al . Combinatorial biosynthesis of sapogenins and saponins in Saccharomyces cerevisiae using a C-16α hydroxylase from Bupleurum falcatum. Proceedings of the National Academy of Sciences. 2014; 111 : 1634–1639. doi: 10.1073/pnas.1323369111 24434554
66 Duan H , Schuler MA . Differential expression and evolution of the Arabidopsis CYP86A subfamily. Plant Physiology. 2005; 137 (3 ): 1067–1081. doi: 10.1104/pp.104.055715 15709153
67 Nakayasu M , Umemoto N , Ohyama K , Fujimoto Y , Lee HJ , Watanabe B , et al . A dioxygenase catalyzes steroid 16α-hydroxylation in steroidal glycoalkaloid biosynthesis. Plant Physiology. 2017; 175 (1 ): 120–133. doi: 10.1104/pp.17.00501 28754839
68 Han JY , Kim HJ , Kwon YS , Choi YE . The Cyt P450 enzyme CYP716A47 catalyzes the formation of protopanaxadiol from dammarenediol-II during ginsenoside biosynthesis in Panax ginseng. Plant and Cell Physiology. 2011; 52 (12 ): 2062–2073. doi: 10.1093/pcp/pcr150 22039120
