Synthesis of Crocin I and Crocin II by Multigene Stacking in Nicotiana benthamiana.

Crocins are a group of highly valuable water-soluble carotenoids that are reported to have many pharmacological activities, such as anticancer properties, and the potential for treating neurodegenerative diseases including Alzheimer's disease. Crocins are mainly biosynthesized in the stigmas of food-medicine herbs Crocus sativus L. and Gardenia jasminoides fruits. The distribution is narrow in nature and deficient in resources, which are scarce and expensive. Recently, the synthesis of metabolites in the heterologous host has opened up the potential for large-scale and sustainable production of crocins, especially for the main active compounds crocin I and crocin II. In this study, GjCCD4a, GjALDH2C3, GjUGT74F8, and GjUGT94E13 from G. jasminoides fruits were expressed in Nicotiana benthamiana. The highest total content of crocins in T1 generation tobacco can reach 78,362 ng/g FW (fresh weight) and the dry weight is expected to reach 1,058,945 ng/g DW (dry weight). Surprisingly, the primary effective constituents crocin I and crocin II can account for 99% of the total crocins in transgenic plants. The strategy mentioned here provides an alternative platform for the scale-up production of crocin I and crocin II in tobacco.

Post-Ripening and Key Glycosyltransferase Catalysis to Promote Sweet Mogrosides Accumulation of Siraitia grosvenorii Fruits.

Sweet mogrosides are not only the primary bioactive ingredient in Siraitia grosvenorii fruits that exhibit anti-tussive properties and expectorate phlegm, but they are also responsible for the fruit's sweetness. Increasing the content or proportion of sweet mogrosides in Siraitia grosvenorii fruits is significant for improving their quality and industrial production. Post-ripening is an essential step in the post-harvest processing of Siraitia grosvenorii fruits, but the underlying mechanism and condition of post-ripening on Siraitia grosvenorii quality improvement need to be studied systematically. Therefore, this study analyzed the mogroside metabolism in Siraitia grosvenorii fruits under different post-ripening conditions. We further examined the catalytic activity of glycosyltransferase UGT94-289-3 in vitro. The results showed that the post-ripening process of fruits could catalyze the glycosylation of bitter-tasting mogroside IIE and III to form sweet mogrosides containing four to six glucose units. After ripening at 35 °C for two weeks, the content of mogroside V changed significantly, with a maximum increase of 80%, while the increase in mogroside VI was over twice its initial amount. Furthermore, under the suitable catalytic condition, UGT94-289-3 could efficiently convert the mogrosides with less than three glucose units into structurally diverse sweet mogrosides, i.e., with mogroside III as the substrate, 95% of it can converted into sweet mogrosides. These findings suggest that controlling the temperature and related catalytic conditions may activate UGT94-289-3 and promote the accumulation of sweet mogrosides. This study provides an effective method for improving the quality of Siraitia grosvenorii fruits and the accumulation of sweet mogrosides, as well as a new economical, green, and efficient method for producing sweet mogrosides.

Root Metabolism and Effects of Root Exudates on the Growth of Ralstonia solanacearum and Fusarium moniliforme Were Significantly Different between the Two Genotypes of Peanuts.

Wild peanut species Arachis correntina (A. correntina) had a higher continuous cropping tolerance than peanut cultivars, closely correlating with the regulatory effects of its root exudates on soil microorganisms. To reveal the resistance mechanism of A. correntina to pathogens, we adopted transcriptomic and metabolomics approaches to analyze differentially expressed genes (DEGs) and differentially expressed metabolites (DEMs) between A. correntina and peanut cultivar Guihua85 (GH85) under hydroponic conditions. Interaction experiments of peanut root exudates with Ralstonia solanacearum (R. solanacearum) and Fusarium moniliforme (F. moniliforme) were carried out in this study. The result of transcriptome and metabolomics association analysis showed that there were fewer up-regulated DEGs and DEMs in A. correntina compared with GH85, which were closely associated with the metabolism of amino acids and phenolic acids. Root exudates of GH85 had stronger effects on promoting the growth of R. solanacearum and F. moniliforme than those of A. correntina under 1 and 5 percent volume (1% and 5%) of root exudates treatments. Thirty percent volume (30%) of A. correntina and GH85 root exudates significantly inhibited the growth of two pathogens. The exogenous amino acids and phenolic acids influenced R. solanacearum and F. moniliforme showing concentration effects from growth promotion to inhibition as with the root exudates. In conclusion, the greater resilience of A. correntina) to changes in metabolic pathways for amino acids and phenolic acids might aid in the repression of pathogenic bacteria and fungi.

Heterologous mogrosides biosynthesis in cucumber and tomato by genetic manipulation.

Mogrosides are widely used as high-value natural zero-calorie sweeteners that exhibit an array of biological activities and allow for vegetable flavour breeding by modern molecular biotechnology. In this study, we developed an In-fusion based gene stacking strategy for transgene stacking and a multi-gene vector harbouring 6 mogrosides biosynthesis genes and transformed it into Cucumis sativus and Lycopersicon esculentum. Here we show that transgenic cucumber can produce mogroside V and siamenoside I at 587 ng/g FW and 113 ng/g FW, respectively, and cultivated transgenic tomato with mogroside III. This study provides a strategy for vegetable flavour improvement, paving the way for heterologous biosynthesis of mogrosides.

Genome-Wide Identification and Expression Analysis of WRKY Transcription Factors in Siraitia siamensis.

WRKY transcription factors, as the largest gene family in higher plants, play an important role in various biological processes including growth and development, regulation of secondary metabolites, and stress response. In this study, we performed genome-wide identification and analysis of WRKY transcription factors in S. siamensis. A total of 59 SsWRKY genes were identified that were distributed on all 14 chromosomes, and these were classified into three major groups based on phylogenetic relationships. Each of these groups had similar conserved motifs and gene structures. We compared all the S. siamensis SsWRKY genes with WRKY genes identified from three diverse plant species, and the results implied that segmental duplication and tandem duplication play an important roles in the evolution processes of the WRKY gene family. Promoter region analysis revealed that SsWRKY genes included many cis-acting elements related to plant growth and development, phytohormone response, and both abiotic and biotic stress. Expression profiles originating from the transcriptome database showed expression patterns of these SsWRKY genes in four different tissues and revealed that most genes are expressed in plant roots. Fifteen SsWRKY genes with low-temperature response motifs were surveyed for their gene expression under cold stress, showing that most genes displayed continuous up-regulation during cold treatment. Our study provides a foundation for further study on the function and regulatory mechanism of the SsWRKY gene family.

Plant Metabolic Engineering by Multigene Stacking: Synthesis of Diverse Mogrosides.

Mogrosides are a group of health-promoting natural products that extracted from Siraitia grosvenorii fruit (Luo-han-guo or monk fruit), which exhibited a promising practical application in natural sweeteners and pharmaceutical development. However, the production of mogrosides is inadequate to meet the need worldwide, and uneconomical synthetic chemistry methods are not generally recommended for structural complexity. To address this issue, an in-fusion based gene stacking strategy (IGS) for multigene stacking has been developed to assemble 6 mogrosides synthase genes in pCAMBIA1300. Metabolic engineering of Nicotiana benthamiana and Arabidopsis thaliana to produce mogrosides from 2,3-oxidosqualene was carried out. Moreover, a validated HPLC-MS/MS method was used for the quantitative analysis of mogrosides in transgenic plants. Herein, engineered Arabidopsis thaliana produced siamenoside I ranging from 29.65 to 1036.96 ng/g FW, and the content of mogroside III at 202.75 ng/g FW, respectively. The production of mogroside III was from 148.30 to 252.73 ng/g FW, and mogroside II-E with concentration between 339.27 and 5663.55 ng/g FW in the engineered tobacco, respectively. This study provides information potentially applicable to develop a powerful and green toolkit for the production of mogrosides.

Heterologous Biosynthesis of Health-Promoting Baicalein in Lycopersicon esculentum.

Baicalein is a valuable flavonoid isolated from the medicinal plant Scutellaria baicalensis Georgi, which exhibits intensive biological activities, such as anticancer and antiviral activities. However, its production is limited in the root with low yield. In this study, In-Fusion and 2A peptide linker were developed to assemble SbCLL-7, SbCHI, SbCHS-2, SbFNSII-2 and SbCYP82D1.1 genes driven by the AtPD7, CaMV 35S and AtUBQ10 promoters with HSP, E9 and NOS terminators, and were used to engineer baicalein biosynthesis in transgenic tomato plants. The genetically modified tomato plants with this construct synthesized baicalein, ranging from 150 ng/g to 558 ng/g FW (fresh weight). Baicalein-fortified tomatoes have the potential to be health-promoting fresh vegetables and provide an alternative source of baicalein production, with great prospects for market application.

Systematic characterization of flavonoids from Siraitia grosvenorii leaf extract using an integrated strategy of high-speed counter-current chromatography combined with ultra high performance liquid chromatography and electrospray ionization quadrupole time-of-flight mass spectrometry.

The chemical constituents of the Siraitia grosvenorii leaf extract were studied. Firstly, high-speed counter-current chromatography was applied to the one-step separation of four compounds from S. grosvenorii leaf extract with the solvent system composed of 0.01% acetic acid water/n-butanol/n-hexane/methanol (5:3:1:1, v/v/v/v). In this work, 270 mg of crude sample yielded four compounds, a new kaempferol O-glycoside derivative, kaempferol 3-O-α-L-[4-O-(4-carboxy-3-hydroxy-3-methylbutanoyl)]-rhamnopyranoside-7-O-α-L-rhamnopyranoside, named kaempferitrin A (2.1 mg, 90%), and three known compounds, grosvenorine (3.4 mg, 93%), kaempferitrin (14.4 mg, 99%) and afzelin (4 mg, 98%), and the structures of these compounds were identified by NMR spectroscopy and mass spectrometry. Then, ultra high performance liquid chromatography with electrospray ionization quadrupole time-of-flight mass spectrometry was used to illustrate the dominant flavonoids in S. grosvenorii leaf extract. 34 flavonoids including 19 kaempferol O-glycosides, 4 quercetin O-glycosides, 6 flavanone derivatives, and 5 polymethoxyflavones, were accurately or tentatively identified by carefully comparing their retention times, UV data, precise masses, the typical fragments of the standards and literature data. Most of these compounds were reported for the first time. This study establishes a foundation for the further development and utilization of S. grosvenorii leaves in future.

Complete chloroplast genomes of two Siraitia Merrill species: Comparative analysis, positive selection and novel molecular marker development.

Siraitia grosvenorii fruit, known as Luo-Han-Guo, has been used as a traditional Chinese medicine for many years, and mogrosides are its primary active ingredients. Unfortunately, Siraitia siamensis, its wild relative, might be misused due to its indistinguishable appearance, not only threatening the reliability of the medication but also partly exacerbating wild resource scarcity. Therefore, high-resolution genetic markers must be developed to discriminate between these species. Here, the complete chloroplast genomes of S. grosvenorii and S. siamensis were assembled and analyzed for the first time; they were 158,757 and 159,190 bp in length, respectively, and possessed conserved quadripartite circular structures. Both contained 134 annotated genes, including 8 rRNA, 37 tRNA and 89 protein-coding genes. Twenty divergences (Pi > 0.03) were found in the intergenic regions. Nine protein-coding genes, accD, atpA, atpE, atpF, clpP, ndhF, psbH, rbcL, and rpoC2, underwent selection within Cucurbitaceae. Phylogenetic relationship analysis indicated that these two species originated from the same ancestor. Finally, four pairs of molecular markers were developed to distinguish the two species. The results of this study will be beneficial for taxonomic research, identification and conservation of Siraitia Merrill wild resources in the future.

Modification of isoprene synthesis to enable production of curcurbitadienol synthesis in Saccharomyces cerevisiae.

Cucurbitane-type triterpenoids such as mogrosides and cucurbitacins that are present in the plants of Cucurbitaceae are widely used in Asian traditional medicine. Cucurbitadienol is the skeleton of cucurbitane-type triterpenoids. As an alternative production strategy, we developed baker's yeast Saccharomyces cerevisiae as a microbial host for the eventual transformation of cucurbitadienol. The synthetic pathway of cucurbitadienol was constructed in Saccharomyces cerevisiae by introducing the cucurbitadienol synthase gene from different plants, resulting in 7.80 mg cucurbitadienol from 1 L of fermentation broth. Improving supplies of isoprenoid precursors was then investigated for increasing cucurbitadienol production. Cucurbitadienol production increased to 21.47 mg/L through the overexpression of a global regulatory factor (UPC2) gene of triterpenoid synthase. In addition, knockout of the ERG7 gene increased cucurbitadienol production from 21.47 to 61.80 mg/L. Finally, fed-batch fermentation was performed, and 63.00 mg/L cucurbitadienol was produced. This work is an important step towards the total biosynthesis of valuable cucurbitane-type triterpenoids and demonstrates the potential for developing a sustainable and secure yeast biomanufacturing platform for triterpenoids.

Functional expression of two NADPH-cytochrome P450 reductases from Siraitia grosvenorii.

Cytochrome P450 reductase (CPR) is the redox partner of various P450s involved in primary and secondary metabolism. Here, we identified and characterized two paralogs of cytochrome P450 reductase from Siraitia grosvenorii. There were two full-length CPR isoforms in the S. grosvenorii fruit transcriptome dataset. They had the same open reading frames of 2, 124 bp, encoding 707 amino acids. A phylogenetic analysis characterized both SgCPR1 and SgCPR2 as Class II dicotyledonous CPRs. The recombinant proteins SgCPR1 and SgCPR2 could reduce cytochrome c and ferricyanide in a NADPH-dependent manner. The SgCPR1 and SgCPR2 transcripts were detected in all examined tissues of S. grosvenorii, and in fresh fruit, they had expression patterns similar to several key enzymes that require CPR as a partner during their biosynthesis. The expression levels of the SgCPRs were induced after a methyl jasmonate treatment. The extracts from yeast co-expressing SgCPR1/SgCPR2 and the cytochrome P450 enzyme CYP76AH1 produced ferruginol, indicating the positive effects of SgCPR1/SgCPR2 on the CYP76AH1 activity. A docking analysis confirmed the experimentally deduced functional activities of SgCPR1 and SgCPR2 for NADPH, FAD and FMN. Thus, SgCRP1 and SgCPR2 are both likely to participate in secondary metabolism, especially mogroside biosynthesis in S. grosvenorii.

Rapid Evaluation of Chemical Consistency of Artificially Induced and Natural Resina Draconis Using Ultra-Performance Liquid Chromatography Quadrupole-Time-of-Flight Mass Spectrometry-Based Chemical Profiling.

Resina Draconis is a highly valued traditional medicine widely used in Arabia since ancient times, and it has been commonly used as an antidiarrheic, antimicrobial, antiulcer, blood circulation promoter as well as an anti-inflammatory agent. The tree source from which this medicine orignates grows extremely slowly, producing a very low yield of Resina Draconis. To meet the increasing market demand, artificial methods for stimulating Resina Draconis formation have been developed and applied. However, the chemical differences between artificially induced Resina Draconis (AIRD) and natural Resina Draconis (NRD) have been rarely studied. The aim of this research was to explore and identify the chemical constituents of AIRD and NRD using ultra-performance liquid chromatography-quadrupole-time-of-flight mass spectrometry (UHPLC-QTOF-MS/MS) based chemical profiling. A total of 56 chromatographic peaks were detected in AIRD, of these, 44 peaks have had their structures tentatively characterized based on high-resolution mass spectra (HRMS) data, fragmentation ions information, reference standards data and literature review. In total, 40 peaks were found both in AIRD and NRD. The potential chemical transformation mechanisms active in Resina Draconis during formation were explored. To the best of our knowledge, this is the first evaluation of the chemical profiles of both AIRD and NRD. Furthermore, these findings are expected to provide a rational basis for the quality assessment of AIRD and the use of AIRD as a substitute for NRD.

[Prospects of molecular breeding in medical plants].

The molecular-assisted breeding, transgenic breeding and molecular designing breeding are three development directions of plant molecular breeding. Base on these three development directions, this paper summarizes developing status and new tendency of research field of genetic linkage mapping, QTL mapping, association mapping, molecular-assisted selections, pollen-mediated transformations, agrobacterium-mediated transformations, particle gun-mediated transformations, genome editing technologies, whole-genome sequencing, transcriptome sequencing, proteome sequencing and varietal molecular designing. The objective and existing problem of medical plant molecular breeding were discussed the prospect of these three molecular breeding technologies application on medical plant molecular breeding was outlooked.

The effect of 8-OH-DPAT and dapoxetine on gene expression in the brain of male rats during ejaculation.

The 5-HT1A receptor agonist 8-hydroxy-2-[di-n-propylamino] tetralin (8-OH-DPAT) promotes ejaculation of male rats, whereas dapoxetine delays this process. However, the gene expression profile of the brain at ejaculation following administrationof these two compounds has not been fully elucidated. In the present study, a transcriptomic BodyMap was generated by conducting mRNA-Seq on brain samples of male Sprague-Dawley rats. The study included four groups: pre-copulatory control (CK) group, ejaculation (EJ) group, 0.5 mg/kg 8-OH-DPAT-ejaculation group (DPAT), and 60 mg/kg dapoxetine-ejaculation (DAP) group. The resulting analysis generated an average of approximately 47 million sequence reads. Significant differences in the gene expression profiles of the aforementioned groups were observed in the EJ (257 genes), DPAT (349 genes) and the DAP (207 genes) compared with the control rats. The results indicate that the expression of Drd1 and Slc6a3 was significantly different after treatment with 8-OH-DPAT, whereas the expression of Drd4 was significantly different after treatment with dapoxetine. Other genes, such as Wnt9b, Cdkn1a and Fosb, exhibited significant differences in expression after the two treatments and are related to bladder cancer, renal cell carcinoma and sexual addiction. The present study reveals the basic pattern of gene expression that was activated at ejaculation in the presence of 8-OH-DPAT or dapoxetine, providing preliminary gene expression information during rat ejaculation.

Cloning and characterization of squalene synthase and cycloartenol synthase from Siraitia grosvenorii.

Mogrosides and steroid saponins are tetracyclic triterpenoids found in Siraitia grosvenorii. Squalene synthase (SQS) and cycloartenol synthase (CAS) are key enzymes in triterpenoid and steroid biosynthesis. In this study, full-length cDNAs of SgSQS and SgCAS were cloned by a rapid amplification of cDNA-ends with polymerase chain reaction (RACE-PCR) approach. The SgSQS cDNA has a 1254 bp open reading frame (ORF) encoding 417 amino acids, and the SgCAS cDNA contains a 2298 bp ORF encoding 765 amino acids. Bioinformatic analysis showed that the deduced SgSQS protein has two transmembrane regions in the C-terminal. Both SgSQS and SgCAS have significantly higher levels in fruits than in other tissues, suggesting that steroids and mogrosides are competitors for the same precursors in fruits. Combined in silico prediction and subcellular localization, experiments in tobacco indicated that SgSQS was probably in the cytoplasm or on the cytoskeleton, and SgCAS was likely located in the nucleus or cytosol. These results will provide a foundation for further study of SgSQS and SgCAS gene functions in S. grosvenorii, and may facilitate improvements in mogroside content in fruit by regulating gene expression.

Systematic identification of flavonols, flavonol glycosides, triterpene and siraitic acid glycosides from Siraitia grosvenorii using high-performance liquid chromatography/quadrupole-time-of-flight mass spectrometry combined with a screening strategy.

The fruits of Siraitia grosvenorii are considered to be health-promoting because of the diversity of their bioactive ingredients. In the present study, a screening method, using high-performance liquid chromatography/quadrupole-time-of-flight mass spectrometry (HPLC-Q-TOF-MS) combined with a screening strategy, has been established. The technology was used to systematically screening the targeted metabolites, primarily from the complex matrix of S. grosvenorii. The compounds were then identified by their exact masses and characteristic fragment ions, in comparison with the fragmentation behaviors of 19 references. Finally, 122 compounds, including 53 flavonols and flavonol glycosides, 59 triterpene glycosides and 10 siraitic acid glycosides, were screened and identified in 10-, 50- and 80-day fruits, roots, stems and leaves of S. grosvenorii. 98 of them were reported for the first time. Additionally, the distribution of all identified components in different parts of the plant was determined and metabolic networks for flavonol and triterpene glycosides were proposed.

Profiles of brain central nervous system gene expression associated with ejaculation behavior in male rats.

The male rat has been used extensively as a model for evaluating the neurophysiology of sexual behavior. However, gene expression in the brain throughout the process of sexual intercourse has yet to be elucidated. In the present study, we created a transcriptomic BodyMap by performing mRNA-Seq on brain samples from pre-copulatory control (CK), fourth intromission (CR4), ejaculation (EJ) and post-ejaculatory interval 1-min (PEI1) Sprague Dawley rats (n=40, all male, each 10). The resulting analysis generated an average of approximately 47 million sequence reads, indicating changes in roughly 21,255 genes for each sample. Among of them, significant differences in gene expression relative to control rats were observed in the CR4 (139 genes), EJ (257 genes), and PEI1 (130 genes) groups. KEGG (Kyoto Encyclopedia of Genes and Genomes) enrichment analysis identified 22 pathway-related genes. We further identified eight important genes related to neural pathways using RT-qPCR and Western blot, ruling out the possibility of false positives. The results of the present study not only revealed the basic pattern of gene expression during male rate sexual activity but also provide preliminary data and methodology for further research regarding animal sexual activity.

Transcriptome Analysis of Taxillusi chinensis (DC.) Danser Seeds in Response to Water Loss.

BACKGROUND: Taxillus chinensis (DC.) Danser, the official species of parasitic loranthus that grows by parasitizing other plants, is used in various traditional Chinese medicine prescriptions. ABA-dependent and ABA-independent pathways are two major pathways in response to drought stress for plants and some genes have been reported to play a key role during the dehydration including dehydration-responsive protein RD22, late embryogenesis abundant (LEA) proteins, and various transcription factors (TFs) like MYB and WRKY. However, genes responding to dehydration are still unknown in loranthus. METHODS AND RESULTS: Initially, loranthus seeds were characterized as recalcitrant seeds. Then, biological replicates of fresh loranthus seeds (CK), and seeds after being dehydrated for 16 hours (Tac-16) and 36 hours (Tac-36) were sequenced by RNA-Seq, generating 386,542,846 high quality reads. A total of 164,546 transcripts corresponding to 114,971 genes were assembled by Trinity and annotated by mapping them to NCBI non-redundant (NR), UniProt, GO, KEGG pathway and COG databases. Transcriptome profiling identified 60,695, 56,027 and 66,389 transcripts (>1 FPKM) in CK, Tac-16 and Tac-36, respectively. Compared to CK, we obtained 2,102 up-regulated and 1,344 down-regulated transcripts in Tac-16 and 1,649 up-regulated and 2,135 down-regulated transcripts in Tac-36 by using edgeR. Among them some have been reported to function in dehydration process, such as RD22, heat shock proteins (HSP) and various TFs (MYB, WRKY and ethylene-responsive transcription factors). Interestingly, transcripts encoding ribosomal proteins peaked in Tac-16. It is indicated that HSPs and ribosomal proteins may function in early response to drought stress. Raw sequencing data can be accessed in NCBI SRA platform under the accession number SRA309567. CONCLUSIONS: This is the first time to profile transcriptome globally in loranthus seeds. Our findings provide insights into the gene regulations of loranthus seeds in response to water loss and expand our current understanding of drought tolerance and germination of seeds.

Cloning and expression of SgCYP450-4 from Siraitia grosvenorii.

CYP450 plays an essential role in the development and growth of the fruits of Siraitia grosvenorii. However, little is known about the SgCYP450-4 gene in S. grosvenorii. Here, based on transcriptome data, a full-length cDNA sequence of SgCYP450-4 was cloned by reverse transcriptase-polymerase chain reaction (RT-PCR) and rapid-amplification of cDNA ends (RACE) strategies. SgCYP450-4 is 1677 bp in length (GenBank accession No. AEM42985.1) and contains a complete open reading frame (ORF) of 1422 bp. The deduced protein was composed of 473 amino acids, the molecular weight is 54.01 kDa, the theoretical isoelectric point (PI) is 8.8, and the protein was predicted to possess cytochrome P450 domains. SgCYP450-4 gene was highly expressed in root, diploid fruit and fruit treated with hormone and pollination. At 10 days after treatment with pollination and hormones, the expression of SgCYP450-4 had the highest level and then decreased over time, which was consistent with the development of fruits of S. Grosvenorii. Hormonal treatment could significantly induce the expression of SgCYP450-4. These results provide a reference for regulation of fruit development and the use of parthenocarpy to generate seedless fruit, and provide a scientific basis for the production of growth regulator application agents.

[Selection of reference genes of Siraitia grosvenorii by real-time PCR].

Siraitia grosvenorii is a traditional Chinese medicine also as edible food. This study selected six candidate reference genes by real-time quantitative PCR, the expression stability of the candidate reference genes in the different samples was analyzed by using the software and methods of geNorm, NormFinder, BestKeeper, Delta CT method and RefFinder, reference genes for S. grosvenorii were selected for the first time. The results showed that 18SrRNA expressed most stable in all samples, was the best reference gene in the genetic analysis. The study has a guiding role for the analysis of gene expression using qRT-PCR methods, providing a suitable reference genes to ensure the results in the study on differential expressed gene in synthesis and biological pathways, also other genes of S. grosvenorii.