Plant Metabolic Engineering by Multigene Stacking: Synthesis of Diverse Mogrosides

Abstract

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.

1. Introduction

Recently, the low-sugar diet gradually became an attractive healthy lifestyle, and plant-derived sweeteners have drawn a tremendous amount of attention around the world. Up to now, many natural sweeteners with high potency have been approved by the U.S. Food and Drug Administration (FDA), including D-psicose (isolated from wheat, fig, and raisin) [1], Siraitia grosvenorii fruits extracts (mogrosides) [2], and stevia leaf extracts [3]. Notably, mogrosides, a group of sweet cucurbitane-type triterpene glycosides in Luo-han-guo, have an increasing market demand due to their superior sweetness compared to rebaudiana A [4]. Furthermore, at least so far, D-psicose is still not approved to be used as a food additive in China yet. Consequently, mogrosides are the representative of an emerging trend in the global natural sweetener market [5]. Mogrosides are commonly considered to be a mixture of cucurbitane-type triterpene glycosides with different numbers of glucoses, which have all kinds of pharmacological bioactivities including antitumor [6,7], anti-inflammatory [8], anti-diabetic [9,10], antiglycation and antioxidant [11]. For instance, mogroside V (five glucoses attached to mogrol) is an essential component at levels of 5.77–11.75 mg/g DW in the fruits of Siraitia grosvenorii. Siamenoside I (four glucoses attached to mogrol) is the sweetest component and has a preferable flavor than other mogrosides. However, the yield is merely at a lower level of 1 mg/g DW in the fruits of S. grosvenorii [12,13]. Additionally, in the mature fruits of Siraitia grosvenorii, most of which are seeds, account for approximately 70% of dry weight [14]. Thus, the plant extraction method is almost impossible to yield high concentrations of natural mogrosides, which has limited the extent of its marketability. Unfortunately, the complexity of mogrosides in special glycosylation sites makes the application of chemical synthesis extremely challenging. Furthermore, Siraitia grosvenorii is distributed only in China, which is limited by a planting environment of high humidity, high temperature and higher diurnal temperature variation. Additionally, to protect designation of origin, its seeds and genetic genes have already been legislated and are prohibited to be taken out of the country. Furthermore, mogrosides are only present in the fruits of Luo-han-guo, which are not detected in the leave, stem, root and seed [15]. This reduces the potential of market competitiveness for mogrosides which are prohibitively expensive and less economical. Hence, a traditional cultivation mode is apparently struggling to meet a surging global market demand of mogrosides, which urgently requires an environmental and efficient toolkit to be developed to produce mogrosides.

Alternatively, major extraordinary progress in synthetic biology has provided tremendous possibility for the mogrosides production in heterologous hosts. Heterologous biosynthesis of natural products in microbial or plant chassis is one of the research hotspots in the field of synthetic biology, which greatly reduces the natural resource damage, environmental pollution and economic benefits [16,17,18,19]. Researchers have published Siraitia grosvenorii genome data and mogrosides biosynthesis pathway has been investigated, which offers a solid theoretical basis for the heterologous production of mogrosides [15,20]. Microbial chassis were preferable in large-scale industrial production, which were comparatively well researched [17,21,22]. Several attempts have been made to explore the mogrosides in Saccharomyces cerevisiae. For example, mogroside V has been transformed into mogroside III E by β-glucosidases in the yeast [23]. Dekkera bruxellensis Exg1p enzyme has been used to hydrolyze mogroside V to siamenoside I, which contained 62.67 ± 3.71% in the fermentation products [24]. However, it has been not reported for the glycosylation resulting in mogroside V, siamenoside I, and mogroside III yet. This probably results from the lack of subcellular organelles required for functional expression of eukaryotic enzymes (e.g., cytochrome P450s), the absence of substrates or poor endogenous metabolic flux towards mogrosides. An alternative approach has potential to produce health-promoting natural products in plant chassis.

In the past decade, Nicotiana benthamiana and Arabidopsis thaliana are fast-growing and high-biomass model plants that can be used as a vehicle for metabolic engineering and have been a popular approach to improve the production of natural products, such as taxadiene, taxadiene-5α-ol [25], artemisinin [26], artemisinic acid [27], unusual fatty acid [28], triterpene [29], and DHA [30], etc. For example, artemisinin is an effective anti-malarial drug isolated from Artemisia annua with the low yields, which cannot meet the high market demand commonly. In 2011, researchers have reconstructed the artemisinin in tobacco, which paves the way for the development of a sustainable plant-based platform for the production of all kinds of secondary metabolites [26]. Additionally, they have a well- defined genetic background and efficient transformation system, which provides a valuable insight into the establishment of a bioreactor for mogrosides in plants. Taking advantage of the elucidated mogrosides biosynthesis pathway, all genes responsible for the biosynthesis have been characterized, and the precursor, 2,3-oxidosqualene commonly exists in plants. As shown in Figure 1, to de novo synthesize the mogrosides, at least 6 encoding genes are required to be transformed into the heterologous plants, including SgSQE1 (squalene epoxidase), SgCS (cucurbitadienol synthase), SgEPH2 (epoxide hydrolases), SgP450 (cytochrome P450 monooxygenase), SgUGT269-1 and SgUGT289-3 (UDP-glucosyltransferases). In fact, an appropriate multigene stacking strategy and transformation is used to produce mogrosides in heterologous plants.

Over the past two decades, there have been major advances in multigene stacking, including GATEWAY [31], Gibson assembly [32], In-Fusion [33], Cre recombinase/loxP-mediated recombination (TransGene Stacking II, TGS II) [34,35], and suppression thermo-interlaced PCR (STI PCR) [36]. Furthermore, the construction of fusion protein with Internal ribosomal entry sites (IRESs) and 2A peptides can be used to assemble multiple genes [37,38,39]. So far, there are abundant researches in the field of heterologous production with diverse plant chassis. β-carotene-biofortified rice endosperm [40], maize [41], canola [42] and soybean [43], anthocyanin-enriched rice endosperm [34], tomato [44], and maize [39]; ginsenoside aglycone rice [45]; astaxanthin-enriched rice [35] and tomato [46], and antimalarial drug, artemisinin, have been synthesized in Nicotiana benthamiana [26] and Physcomitrella patens [47]. Recently, multigene assembly strategy mainly focuses on 2-10 genes transformation in plants, and with the increasing number of transgenes, the vector gradually become unstable, thus leading to gene silencing. Therefore, an efficient multigene stacking strategy is important to improve the efficiency and applicability of multigene vector in all kinds of plant species for future research.

Here, in-fusion based gene stacking strategy (IGS) combined with 2A peptides was developed to assemble 6 mogrosides synthase genes and a mogrosides biosynthesis pathway has been engineered in tobacco and Arabidopsis. Our study also yielded a green toolkit for the development of multiple plant materials that are sources of mogrosides and have great market potential. This offered an alternative for production of mogrosides in perspective heterologous plants.

2. Results

2.1. Multigene Expression Vector Construction

Mogrosides biosynthesis pathway has been studied extensively in Siraitia grosvenorii, 2,3-oxidosqualene was used as the substrate, which is existed widely in all kinds of plants [15]. To achieve heterologous biosynthesis of mogrosides, at least 6 mogrosides synthase genes were required to be introducing into the heterologous plants. For this reason, we developed an in-fusion based multigene stacking strategy with P2A peptides linker. Firstly, SgSQE1, SgCS, SgEPH2, SgP450, SgUGT269-1 and SgUGT289-3 were selected and cloned into the PBI121 (Figure S2A) driven by AtPD7, AtUBQ10, and CaMV 35S promoters, respectively. Secondly, the sequence of promoter + target gene + terminator was cloned and inserted into the pCAMBIA1300 (Figure S2B). After that, single-gene expression cassette and the double-gene expression cassette were combined into the pCAMBIA1300, where the triple-gene expression cassette were generated. Finally, the region of UBQ10:SgUGT269-1:Tmas::35S:SgUGT289-3:Tnos::35S:SgP450 was amplified and ligated with the sequence of SgSQE1:Thsp::35S:SgCS:Tnos::35S:SgEPH2:Tnos by P2A peptides. The multigene vector U22p-SCE was approximately 21.5 kb (from the left border to the right border) (Figure S2C). This construct was used to produce mogrosides in the tobacco and Arabidopsis thaliana, which provided an alternative toolkit for mogrosides production in Nicotiana benthamiana and Arabidopsis thaliana. In our study, the amino acid sequence of Gly-Ser-Gly has been introduced in N-terminal of P2A peptides, which improved the splicing efficiency [48]. This strategy highlights the potential of the natural products’ production in many heterologous plants.

2.2. Transient Expression Assays

To investigate whether the transformation of a multigene vector can produce the mogrosides, transient expression was developed in Nicotiana benthamiana. Six mogrosides synthase genes were introduced and expressed in tobacco using Agrobacterium infiltration. HPLC-MS/MS analysis indicated that relatively less mogrol has accumulated in the leaves of tobacco infiltrated with U22p-SCE at 36 h. There was no mogrol found in the WT tobacco leaves. However, no mogrosides have been synthesized in the leaves of tobacco (Figure S3). The reasons for the lack of mogrosides production in transient expression assay are mainly as follows: the transient expression assay has a short reaction time, which may cause the shortage of substrate accumulation. There is another reason to this result which is the stability of multigene expression vector, which has an adverse effect on enzyme activities. Furthermore, the lack of 2,3-oxidosqualene has contributed to limited accumulation of mogrosides in tobacco. In this case, a stable transformation of multigene expression vectors is necessary for confirming the mogrosides production in Nicotiana tabacum and Arabidopsis thaliana further.

2.3. Generation of Transgenic Tobacco and Arabidopsis thaliana

To assess the availability of Nicotiana benthamiana and Arabidopsis thaliana as plant chassis for mogrosides production, the U22p-SCE plasmid harboring mogrosides synthase genes were transformed into plants via Agrobacterium tumefaciens (Supplementary Materials, Figure S4). Sixteen Hyg-resistant tobacco transgenic lines were obtained. The sequences of the SgSQE1, SgCS, SgEPH2, SgP450, SgUGT269-1, SgUGT289-3, and Hyg genes were isolated from transgenic plants (Figure 2A). All the mogrosides biosynthesis-related genes were detected in six (N16, N22, N30, N32, N45, and N47) out of sixteen transgenic tobacco lines via PCR, which demonstrated that all target genes in the multigene vector were integrated into the tobacco genome. Then, the expression levels of target genes were measured by qPCR. All the target genes were overexpressed in the transgenic tobacco lines; however, no target genes were detected in the WT plants (Figure 2B). These data showed that the mogrosides biosynthesis-related genes were successfully expressed in Nicotiana benthamiana. Likewise, to select and generate transgenic Arabidopsis thaliana lines containing mogrosides biosynthesis-related genes, Arabidopsis thaliana (Col-0) plants were subjected to the Agrobacterium-mediated floral dip method (Figure 3A). T0 seeds were selected on MS media including Hyg, and 11 independent transgenic lines were identified by using PCR-based detection (Figure 3B). Subsequently, 7 transgenic lines were selected for RT-PCR analysis (Figure 3C).

The PCR analysis showed that the ratio of positive transgenic plants was 37.5% and 63.6% for tobacco and Arabidopsis, respectively, which were relatively lower. Sometimes, gene silencing and gene loss occurred in the plant transformation. This was probably due to the existence of redundant repetitive sequences and that the T-DNA fragment was large [49]. To avoid the adverse effect on the gene expression level and enzyme catalytic efficiency in heterologous plants, a different promoter and terminator may attempt to be introduced in the multigene vector, and codon optimization is a good choice to improve gene expression. Additionally, antibiotic resistance gene in transgenic plants should be evaluated. Hygromycin B (Hyg) was used as a selectable marker protein in our study, and concentration of hygromycin was highly negatively correlated with the growth of shoots [50]. Explants were severely yellowed and withered due to a high level of hygromycin in the plant transformation. Conversely, lower concentration may cause an extremely difficult screening process for positive plants. Furthermore, to improve the transformation efficiency of a multigene vector, an antibiotic resistance gene should be replaced by kanamycin or glyphosate, which was commonly used in the plant’s genetic engineering.

2.4. Heterologous Production of Mogrosides in Transgenic Tobacco and Arabidopsis thaliana

Our goal was to produce and obtain transgenic lines containing mogrosides by multigene transformation, that is, to develop a promising bio-factory for mogrosides production. It is apparent from Figure 4A,B that, compared with mogrosides, there was a majority of MII-E and a few of MIII accumulated in the transgenic tobacco leaves, respectively (Figure 4A,D). The average content of MII-E was 399.27, 3745.54, 5283.57, and 5663.55 ng/g FW, and the average content of MIII was 148.30, 212.92, and 252.73 ng/g FW (Figure 4B,C). It seemed that MI-A was completely catalyzed into the MII-E, but the MII-E was redundant. Surprisingly, HPLC-MS/MS-based analysis indicated that MIII and SI were present in transgenic Arabidopsis lines (Figure 5A–C). The average yields of SI were 29.65, 224.57, and 1036.96 ng/g FW in the AA3, AA6, and AU7 transgenic lines, respectively (Figure 6). Moreover, MIII was found only in the AA3 transgenic line, the level of which was 202.75 ng/g FW (Figure 6). MII-E, MIII, and SI production in transgenic tobacco and Arabidopsis lines suggests that these lines could represent novel platforms for mogrosides production.

HPLC-MS/MS analysis indicated that a plethora of MIII and a trace of MII-E were detected in the transgenic tobacco lines, and slight MIII and SI accumulated in transgenic Arabidopsis lines. Generally, the sweetness of mogrosides mainly depends on the number of glucose moieties and glycosylation sites. Mogroside V with five glucose moieties, the sweetness is approximately 300 times higher than that of sucrose [12,51]. Interestingly, siamenoside I with four glucose moieties was the sweetest component in the mogrosides and is 563 times sweeter than 5% sucrose because of the different glycosylation sites. Although the sweetness of mogroside III (three glucoses attached to mogrol) is tasteless, its intestinal maltase inhibitory effect has been investigated previously [52]. Previous studies commonly obtained siamenoside I, mogroside III or other mogrosides by hydrolyzing glycosidic bonds in the yeast via the biotransformation of mogroside V [23,24,53]. As a result, the strategy in the present study has significant potential for the production of SI in transgenic Arabidopsis plants via biotechnology.

3. Discussion

In the last few decades, plant synthetic metabolic engineering is representative of an emerging trend in the production of natural compounds [54]. Firstly, it is based on the identification of key synthase genes, gene flux, and regulatory mechanism for biosynthetic pathway of target compounds; then, efficient multigene stacking and transformation are required to establish in the organisms, which is a pivotal step in this research [55]. Generally, one or more exogenous genes often need to be integrated into the heterologous host genome, which is required to develop all kinds of multigene stacking strategies, such as crossbreeding, re-transformation, and co-transformation. However, these strategies are not only time-consuming but also require randomness and blindness. Currently, multigene vector transformation, that is, multiple genes which are linked in a single T-DNA region, has a great advantage over the other traditional methods [56]. All genes are inserted into a chromosomal site in the plant host genomes by Agrobacterium-mediated transformation, and inherited as a single unit. Therefore, the loss of exogenous genes can be reduced and can improve the transgenic efficiency. Due to the significant advantages of the multigene vector transformation, this method is more widely used in current synthetic biology, and a variety of multigene stacking strategies have been developed recently [57,58,59,60,61,62,63,64]. Researchers developed a TGS II system for multigene assembly, and astaxanthin, β-carotene, and anthocyanins have been produced in rice endosperm by the transformation of 4-10 structural genes, respectively [34,35]. Protopanaxadiol was successfully synthesized in the transgenic rice by the integration of ginseng PgDDS and PgCYP716A47 genes using GATEWAY [45]. Gibson assembly has been used to engineer the taxadiene and taxadiene-5α-ol biosynthesis pathway in tobacco [25]. These multigene stacking strategies were rarely used for 6 or more genes transformation except for patented system TGS II. In this case, to achieve the de novo synthesis of mogrosides in heterologous plants, a simple and efficient multigene stacking needed to be developed. Therefore, we attempted to assemble 6 mogrosides synthase genes via in-fusion technology and self-cleaving 2A peptides, which provided an alternative toolkit for mogrosides production in Nicotiana benthamiana and Arabidopsis thaliana. In our study, the amino acid sequence of Gly-Ser-Gly has been introduced in N-terminal of P2A peptides, which improved the splicing efficiency [48]. This strategy highlights the potential of the natural products’ production in many heterologous plants.

The PCR analysis showed the ratio of positive transgenic plants were 37.5% and 63.6% for tobacco and Arabidopsis, respectively. The multigene assembly and transformation carrying a large multigene vector cause the unstable expression of multigene in the host plants which are in addition to the existence of redundant repetitive sequences which commonly are the initiator of partial or all foreign genes silencing [49]. These make the transgenic efficiency of a multigene vector relatively lower. To avoid the adverse effect on the gene expression level and enzyme catalytic efficiency in heterologous plants, a different promoter and terminator should be introduced in the multigene vector. Additionally, codon optimization is a good choice to improve gene expression and the antibiotic resistance gene in transgenic plants should be evaluated. Hygromycin B (Hyg) was used as a selectable marker protein in our study, and concentration of hygromycin was highly negatively correlated with the growth of shoots [50]. Furthermore, to improve the transformation efficiency of a multigene vector, an antibiotic resistance gene should be replaced by kanamycin or glyphosate, which was commonly used in the plant’s genetic engineering.

HPLC-MS/MS analysis indicated that a plethora of MII-E and a trace of MII-E were detected in the transgenic tobacco lines, and less MIII and SI accumulated in transgenic Arabidopsis lines. Generally, the sweetness of mogrosides differs depending on the number of glycosylation and glycosylation sites. Previous studies have indicated that in mogroside V with five glucose moieties, the sweetness is approximately 300 times higher than that of sucrose [12,51]. Nevertheless, siamenoside I with four glucose moieties was the sweetest component in the mogrosides and is 563 times sweeter than 5% sucrose because of the different glycosylation sites. In addition, the flavor of siamenoside I was better than that of mogroside V [24]. Additionally, mogroside III has three glucose moieties, which is somehow less sweet than mogroside V [65]. Studies found that the extremely low level of siamenoside I in natural extracts has limited the extent of its marketability [25]. Unfortunately, the complexity of siamenoside I in special glycosylation sites makes the application of chemical synthesis extremely challenging. As a result, the strategy in the present study has significant potential for the production of SI in transgenic Arabidopsis plants via biotechnology.

Besides, multigene vector U22p-SCE harbouring 6 mogrosides biosynthesis gene has been assembled via in-fusion technology and self-cleaving 2A peptides, and transformed into the tobacco and Arabidopsis. Although all candidate genes involved in mogroside biosynthesis pathway were highly expressed in the transgenic lines, mogroside V was not detected in the transgenic lines. There are several possible explanations for this result. Firstly, a possible explanation for these results may be the lack of an adequate precursor accumulation in the tobacco and Arabidopsis. To reduce the consumption of substrates, there was an introduction of various rate-limiting enzymes, such as 3-hydroxy-3-methyl glutaryl coenzyme A reductase (HMGR) and squalene synthase (SQS), which is beneficial to the accumulation of substrates. Previous research has found that overexpression of PnSS and PnHMGR could significantly improve the accumulation of total saponins in Panax notoginseng. Additionally, the yield of co-overexpression of PnSS and PnHMGR was 3-fold and 1.5-fold higher than the control and PnHMGR-overexpressed lines, respectively [64]. Therefore, the introduction of rate-limiting enzymes can facilitate the yield of mogrosides in the plant hosts. Another possible explanation for this might be the protein interaction and localization in transgenic plants, which have impacted greatly on the catalytic efficiency of enzyme. We need to attempt to predict the signal peptides and localization for target proteins, and compartmentalized engineering can be used to produce mogrosides in heterologous plants.

Together, our study is highly important because it marks the success of the transformation of mogroside biosynthesis genes in tobacco and Arabidopsis, which offers a perspective on producing plant materials that could serve as a source of multiple mogrosides. This strategy offers an excellent alternative to the traditional extraction methods, which is high-valued toolkit for the production of mogrosides. Furthermore, this study provided a necessary theoretical foundation for mogrosides production in heterologous plants.

4. Materials and Methods

4.1. Plant Material, Chemicals, and Strains

Nicotiana benthamiana and Arabidopsis thaliana (ecotype Columbia-0) were used for transformation assay. The binary plasmid PBI121 and pCAMBIA1300 were used for multigene vector construction. E. coli competent cells including DH5α and XL10-Gold (WeidiBio, Shanghai, China) were used in this study, and Agrobacterium tumefaciens competent cell GV3101 (WeidiBio, Shanghai, China) was used for multigene vector transformations.

HPLC-grade methanol, acetonitrile, and formic acid were obtained from Fisher (Emerson, IA, USA). ClonExpress II One Step Cloning Kit, ClonExpress MultiS One Step Cloning Kit, ClonExpress Ultra One Step Cloning Kit, ClonExpress Ultra One Step Cloning Kit, HiScript Ⅲ 1st Strand cDNA Synthesis Kit (+gDNA wiper), and Taq Pro Universal SYBR qPCR Master Mix were obtained from Vazyme Biotech Co., Ltd. (Nanjing, China). Ultrapure RNA Kit (DNase I) was obtained from CWBIO. Co., Ltd. (Beijing, China). Plant genomic DNA kits, restriction enzymes, and DNA Marker were purchased from TianGen Biotech Co., Ltd. (Beijing, China). KOD One PCR master Mix was obtained from TOYOBO Biotech Co., Ltd. (Shanghai, China). The PBI121 and pCAMBIA1300 plasmids were stored in the laboratory. Mogrosides standards, including mogrol, mogroside II-E (MII-E), mogroside III (MIII) and siamenoside I (SI) were purchased from Chengdu Must Bio-Technology Co., Ltd. (Sichuan, China). Other chemical reagents were obtained from Beijing Chemical Corporation (Beijing, China) unless otherwise specified.

4.2. Multigene Stacking

The coding genes of SgSQE1 (squalene epoxidase), SgCS (cucurbitadienol synthase), SgEPH2 (epoxide hydrolases), SgP450 (cytochrome P450 monooxygenase), SgUGT269-1, and SgUGT289-3 (UDP-glucosyltransferases) were isolated from Siraitia grosvenorii fruit. Ubiquitin 10 (AtUBQ10) and serine carboxypeptidase-like AtSCPL30 (PD7, 456 bp) promoters were cloned from Arabidopsis thaliana using KOD One PCR Master Mix. CaMV 35S promoter and NOS terminator were amplified from the PBI121 plasmid, and mannopine synthase (MAS) terminators and heat-shock protein (HSP) 18.2 terminators were chemically synthesized by GENEWIZ (Suzhou, China).

The strategy of the multigene vector construction in this study is shown in Figure S1. Firstly, SgCS, SgEPH2, SgP450, and SgUGT289-3 were subcloned into the BamHI and SacI sites of PBI121 vector with CaMV 35S promoter and the NOS terminator to generate 35S:SgCS:Tnos, 35S:SgEPH2:Tnos, 35S:SgP450:Tnos, and 35S:SgUGT289-3:Tnos using a ClonExpress II One Step Cloning Kit. Simultaneously, AtPD7:SgSQE1:Thsp and UBQ10:SgUGT269-1:Tmas were constructed using a ClonExpress MultiS One Step Cloning Kit. Secondly, PD7:SgSQE1:Thsp and 35S:SgCS:Tnos were fused into the pCAMBIA1300 at the EcoRI and HindIII sites to generate PD7:SgSQE1:Thsp::35S:SgCS:Tnos. Similarly, another double-gene expression cassette, UBQ:SgUGT269-1:Tmas::35S::SgUGT289-3:Tnos was constructed. Thirdly, PD7:SgSQE1:Thsp::35S:SgCS:Tnos and 35S:SgEPH2:Tnos, UBQ10:SgUGT269-1:Tmas::35S:SgUGT289-3:Tnos and 35S:SgP450:Tnos were ligated into the EcoRI/HindIII sites of the pCAMBIA1300 plasmid, respectively. Finally, UBQ10:SgUGT269-1:Tmas::35S:SgUGT289-3:Tnos::35S:SgP450 (6.1 kb) and SgSQE1:Thsp::35S:SgCS:Tnos::35S:SgEPH2:Tnos (6.5-kb) were ligated by a 2A peptide linker from porcine teschovirus (amino acid sequence is GSGATNFSLLKQAGDVEENPGP) (P2A) via a ClonExpress Ultra One Step Cloning Kit. The multigene expression vector U22p-SCE was assembly. All primers used for multi-gene vector construction are listed in the Table S1.

4.3. Transient Expression Assay

To analyze the availability of multigene expression vector, multigene vector U22p-SCE was introduced into Agrobacterium tumefaciens strain GV3101 for transient expression assay in Nicotiana benthamiana. The Agrobacterium strain with U22p-SCE was cultured in LB medium (25 mg/L rifampicin and 50 mg/L kanamycin) at 28 °C. Then, cells were centrifuged for 15 min at 5000× g and resuspended in buffer with 10 mM 2-(N-morpholino) ethanesulfonic acid (MES), 10 mM MgSO₄, and 200 μM acetosyringone (AS) to the OD₆₀₀ = 0.6. They were kept in the dark about 2-4 h. The strains were infiltrated into the leaves of Nicotiana benthamiana using a needleless syringe. Non-infiltrated plants were as the negative control. The assay was performed independently six times.

4.4. Agrobacterium-Transformation in Tobacco

Tobacco transformation was performed according to the leaf disc method, with slight modifications. In short, it was primarily necessary to aseptically cut the leaves of one-month-old Nicotiana benthamiana plants for the leaves to be used as explants. These leaf explants were pre-cultivated in MS media consisting of 3% sucrose, 0.7% agar, 1 mg/L 6-BA, and 0.2 mg/L naphthylacetic acid (NAA) for 2 d in the dark at 25 °C. Subsequently, leaf explants of Nicotiana benthamiana were inoculated with multigene expression vectors via Agrobacterium tumefaciens strain GV3101. Next, the explants were transferred to cocultivation media (MS media including 3% sucrose, 0.7% agar, 1 mg/L 6-BA, 0.2 mg/L NAA, and 100 μM AS) for 3 d in the dark at 25 °C. The explants were incubated on regeneration media (MS media including 3% sucrose, 0.7% agar, 1 mg/L 6-BA, 0.2 mg/L NAA, 400 mg/L cefotaxime, and 7 mg/L Hyg) at 25 °C under a 16 h/8 h (light/dark) photoperiod. When the shoots were approximately 2-3 cm in length, rooting media (MS media including 3% sucrose, 0.7% agar, 0.2 mg/L NAA, 400 mg/L cefotaxime, and 7 mg/L Hyg) were used to cultivate the regenerated shoots. The resistant plants were identified by PCR for the 6 candidate genes and Hyg-resistant gene using specific primers (Table S2). Wild-type (WT) plant was used as a negative control.

4.5. Arabidopsis thaliana Transformation

For Arabidopsis thaliana transformation, the Col-0 ecotype was used. The multigene vector U22p-SCE was introduced into the Agrobacterium tumefaciens strain GV3101 through the freeze-thaw method, and then the floral-dip method involving Agrobacterium harboring the multigene expression vector was used as previously described [66]. To screen for the transgenic lines, T0 seeds were harvested and grown on selective media including Hyg (30 mg/L), and then the T1 transgenic lines were grown in a greenhouse under nonstress conditions for further study. The resistant plants were detected by PCR using specific primers (Table S2). In addition, genomic DNA from WT plants was a negative control.

4.6. Gene Expression Level Analysis

Total RNA was extracted from 0.1 g transgenic plant leaves using a CWBIO RNA extraction kit. The first-strand cDNA was s reverse transcribed from 1 μg of total RNA via HiScript Ⅲ 1st Strand cDNA Synthesis Kit (+gDNA wiper). The relative gene expression level of SgSQE1, SgCS, SgEPH2, SgP450, SgUGT269-1, and SgUGT289-3 were measured by real-time quantitative PCR (qRT-PCR) using Taq Pro Universal SYBR qPCR Master Mix. An ABI CFX96^TM Real-Time System (Waltham, Massachusetts, USA) with the following reaction procedure: 95 °C for 30 s; 40 cycles of 95 °C for 3 s, 55 °C for 10 s. Nbactin (Accession number: XM_016619439) was internal control. This assay was conducted for multiple technical replicates. All the primers used for qRT-PCR are listed in Table S3.

For the gene expression level analysis of Arabidopsis, Reverse Transcription PCR (RT-PCR) was performed according to the protocol. PCR detection with equivalent cDNA as template was used to determine the gene expression level of SgSQE1, SgCS, SgEPH2, SgP450, SgUGT269-1, and SgUGT289-3 using the specific primers (Table S3) via KOD One PCR Master Mix. Atactin (Accession number: NM_001338359) was used as positive control. WT plants was the negative control. The gene expression of transgenes in Arabidopsis were analyzed with DNA gel electrophoresis.

4.7. HPLC-MS/MS Analysis of Mogrosides

The quantitative analysis of mogrosides was performed as previously described with the HPLC-MS/MS method with minor modification [67]. Mogrosides were extracted from 4 g transgenic tobacco leaves and 5 g Arabidopsis leaves using 80% methanol by ultrasound extraction. The extraction conditions were as follows: 40 kHz frequency for 1 h at room temperature, and then centrifuged at 5000× g for 20 min. After, the collected supernatant was concentrated with nitrogen and diluted in 1 mL methanol.

A 4500 QTRAP LC-MS/MS (AB SCIEX, Toronto, ON, Canada) with an Agilent Poroshell 120 SB C18 column (100 mm × 2.1 mm, 2.7 µm) were used for HPLC-MS/MS analysis. The mobile phase was 0.1% formic acid-water (A) and acetonitrile (B) with a flow rate of 0.2 mL/min. The HPLC condition of mogrosides analysis was as following: 20% B for 0 min; 23% B for 3–5 min; 40% B for 18 min; 20% B for 18.01–20.10 min. For mogrol analysis, 20% B for 0 min; 30% B for 0.5 min; 88% B for 2–4 min, and 20% B for 5.50–8.00 min. The column effluent was monitored by mass spectrometry with electrospray ionization in positive mode and the multiple reaction monitoring (MRM) scanning was employed for quantification. The optimized mass spectrometric conditions and parameters were shown in Table 1. In the current analysis conditions, good linearity was achieved within the investigated ranges (50–5000 ng/mL) for MII-E, MIII, and SI. Typical equations for the standard curves of MII-E, MIII, and SI were y = 173.79x + 4960 (r² = 0.9998), y = 193.16x + 5616.2 (r² = 0.9997), and y = 48.64x + 75.78 (r² = 0.9976), respectively. Therefore, the contents of mogrosides in transgenic plants were calculated by the corresponding equation.

4.8. Data Analysis

SPSS 16.0 statistics program (IBM Co., Armonk, NY, USA) was used to analyze the data in this study. All graphs were illustrated by Origin 2019b (OriginLab Co., Northampton, MA, USA). Three biological and technical replicates were conducted in all experiments, respectively.

Author Contributions

Conceptualization, X.M. and J.L.; methodology, J.L., L.X., T.L., C.M., X.J. and Z.L.; formal analysis, L.X., T.L., X.H., S.C. and Z.L.; resources, F.L.; writing—original draft preparation, J.L.; writing—review and editing, J.L., T.L., Z.L., and X.M.; supervision, X.M.; project administration, X.M.; funding acquisition, X.M., Z.L., and J.L. All authors have read and agreed to the published version of the manuscript.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The authors confirm that all data in this experiment are available in the main text and Supplementary Information.

Conflicts of Interest

The authors declare no conflict of interest.

Sample Availability

All samples are available from the authors.

Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Figures and Table

Figure 1. The Flow Chart of mogrosides biosynthetic pathway in Siraitia grosvenorii fruits. Mogrosides synthase genes which were transformed in this study are marked in blue, including SgSQE1, squalene epoxidases; SgCS, curbitadienol synthase; SgEPH2, epoxide hydrolases; SgP450, cytochrome P450 mono-oxygenase; SgUGT269-1 and SgUGT289-3, UDP-glucosyltransferases; The substrate, 2,3-oxidosqualene is shown with a yellow background. MI-A1, mogroside I-A1; MII-E, mogroside II-E; MIII, mogroside III; SI, siamenoside I; MV, mogroside V (in the dotted bordered rectangle).

Figure 2. Molecular analysis of transgenic tobacco lines. (A) PCR-based analysis of the transgenic tobacco lines. The lanes from left to right represent the WT, N16, N22, N30, N31, N32, N45, N46, N47, and N48. An image of the DNA marker (4.5 kb) is shown in the upper-left corner of the figure. (B) Relative expression level analysis of 6 mogrosides biosynthesis genes in transgenic tobacco lines (N16, N22, N30, N32, N45, N47). The Nbactin is used as an internal control. Expression of tobacco WT plants was set to 1. The data are presented as the mean values ± SDs, n = 3 biologically independent samples, ** represents significant difference at p < 0.01 (LSD test).

Figure 3. Molecular analysis of transgenic Arabidopsis lines. (A) Arabidopsis WT plants and transgenic lines. (B) PCR-based analysis of transgenic Arabidopsis lines. The lanes from left to right represent the WT, AA3, AA5, AA6, AA7, AU6, AU7, AU8, AU10, AU11, AU12, and AU13. DNA marker (4500 bp) is used in the figure. (C) RT-PCR detection of 6 mogrosides biosynthesis genes in the WT and transgenic Arabidopsis lines AA3, AA5, AA6, AU7, AU10, AU11, A12 (from left to right). The Atactin is used as an internal control.

Figure 4. Production of mogrosides in transgenic tobacco lines. (A) HPLC-MS/MS analysis of MIII and MII-E in transgenic tobacco lines. (B) Accumulation of MIII in transgenic tobacco lines. (C) Accumulation of MII-E in transgenic tobacco lines. n.d., not detected. The data are presented as the mean values ± SDs, n = 3 biologically independent samples. The black arrows indicate the peak of MIII and MIIE. (D) Full-scan product ion of MIII and MII-E.5.0 × 104.

Figure 4. Production of mogrosides in transgenic tobacco lines. (A) HPLC-MS/MS analysis of MIII and MII-E in transgenic tobacco lines. (B) Accumulation of MIII in transgenic tobacco lines. (C) Accumulation of MII-E in transgenic tobacco lines. n.d., not detected. The data are presented as the mean values ± SDs, n = 3 biologically independent samples. The black arrows indicate the peak of MIII and MIIE. (D) Full-scan product ion of MIII and MII-E.5.0 × 104.

Figure 5. Production of mogrosides in transgenic Arabidopsis lines. (A) HPLC-MS/MS analysis of MIII in transgenic Arabidopsis lines. (B) HPLC-MS/MS analysis of SI in transgenic Arabidopsis lines. The black arrows indicate the peak of MIII and SI. (C) Full-scan product ion of SI.

Figure 6. Accumulation of mogrosides in transgenic Arabidopsis lines AA3, AA6, and AU7. n.d., not detected. The data are presented as the mean values ± SDs, n = 3 biologically independent samples.

Table 1

HPLC-MS/MS parameters.

Analytes	Molecular Formula	Retention Time (min)	Production (m/z)	DP (V)	CE (eV)
SI	C54H92O24	9.62	1123.6/961.61123.6/799.2	−220	−75
MIII	C48H82O19	11.86	961.6/799.4961.6/637.3	−170	−70
MII-E	C42H72O14	13.82	799.5/637.5799.5/475.5	−170	−65
Mogrol	C30H52O4	4.01	459.3/441.2459.3/423.3	80	20
MS parameters		Mogrosides		Mogrol
Ion mode		Negative		Positive
Source temperature (℃)		550		550
Ionization voltage (V)		−4500		5500
GS1 (psi)		55		18
GS2 (psi)		55		20
CUR (psi)		20		20
CAD		Medium		Medium
Dwell time (ms)		100		200
EP (V)		−10		10
CXP (V)		−15		10

Supplementary Materials

Supplementary materials can be found at https://www.mdpi.com/article/10.3390/ijms231810422/s1.

References

1. Food and Drug Administration, U.S.A GRAS Notice 000693: Generally Recognized as Safe (GRAS) Status of D-Allulose (D-Psicose). Available online: https://www.fda.gov/Food/IngredientsPackagingLabeling/GRAS/NoticeInventory (accessed on 2 February 2017).

2. Food and Drug Administration, U.S.A GRAS Notice 000253: Rebaudioside A for Use as a General Purpose Sweetener. Available online: https://www.fda.gov/Food/IngredientsPackagingLabeling/GRAS/NoticeInventory (accessed on 14 March 2008).

3. Food and Drug Administration, U.S.A GRAS Notice 000301: Luo Han Fruit. Available online: https://www.fda.gov/Food/IngredientsPackagingLabeling/GRAS/NoticeInventory (accessed on 15 January 2010).

4. Sun-Waterhouse, D.; Wadhwa, S.S. Industry-relevant approaches for minimising the bitterness of bioactive compounds in functional foods: A review. Food Bioprocess Technol.; 2013; 6, pp. 607-627. [DOI: https://dx.doi.org/10.1007/s11947-012-0829-2]

5. Pawar, R.S.; Krynitsky, A.J.; Rader, J.I. Sweeteners from Plants-with Emphasis on Stevia Rebaudiana (Bertoni) and Siraitia Grosvenorii (Swingle). Anal. Bioanal. Chem.; 2013; 405, pp. 4397-4407. [DOI: https://dx.doi.org/10.1007/s00216-012-6693-0] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23341001]

6. Matsumoto, S.; Jin, M.; Dewa, Y.; Nishimura, J.; Moto, M.; Murata, Y.; Shibutani, M.; Mitsumori, K. Suppressive Effect of Siraitia Grosvenorii extract on dicyclanil-promoted hepatocellular proliferative lesions in male mice. J. Toxicol. Sci.; 2009; 34, pp. 109-118. [DOI: https://dx.doi.org/10.2131/jts.34.109] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19182440]

7. Liu, C.; Cai, T.Y.; Zhao, X.M.; Ma, L.Q.; Dou, D.Q.; Sun, Y.X. Apoptosis effect of Siraitia Grosvenorii Extracts on lung cancer cells A549 and its mechanisms. Chin. Pharm. Bull.; 2015; 31, pp. 1310-1314. [DOI: https://dx.doi.org/10.3969/j.issn.1001-1978.2015.09.025]

8. Tao, L.; Cao, F.; Xu, G.; Xie, H.; Zhang, M.; Zhang, C. Mogroside IIIE attenuates LPS-induced acute lung injury in mice partly through regulation of the TLR4/MAPK/NF-ΚB axis via AMPK activation. Phytother. Res.; 2017; 31, pp. 1097-1106. [DOI: https://dx.doi.org/10.1002/ptr.5833]

9. Suzuki, Y.A.; Tomoda, M.; Muratal, Y.; Inui, H.; Sugiura, M.; Nakano, Y. Antidiabetic effect of long-term supplementation with Siraitia Grosvenori on the spontaneously diabetic goto-kakizaki rat. Brit. J. Nutr.; 2007; 97, pp. 770-775. [DOI: https://dx.doi.org/10.1017/S0007114507381300]

10. Zhou, G.; Zhang, Y.; Li, Y.; Wang, M.; Li, X. The metabolism of a natural product mogroside V, in healthy and type 2 diabetic rats. J. Chromatogr. B; 2018; 1079, pp. 25-33. [DOI: https://dx.doi.org/10.1016/j.jchromb.2018.02.002]

11. Liu, H.; Wang, C.; Qi, X.; Zou, J.; Sun, Z. Antiglycation and antioxidant activities of mogroside extract from Siraitia Grosvenorii (Swingle) fruits. J. Food Sci. Tech.; 2018; 55, pp. 1880-1888. [DOI: https://dx.doi.org/10.1007/s13197-018-3105-2]

12. Makapugay, H.C.; Dhammika Nanayakkara, N.P.; Soejarto, D.D.; Kinghorn, A.D. High-performace liquid chromatographic analysis of the major sweet principle of lo han kuo fruits. J. Agric. Food Chem.; 1985; 33, pp. 348-350. [DOI: https://dx.doi.org/10.1021/jf00063a007]

13. Li, D.; Ikeda, T.; Huang, Y.; Liu, J.; Nohara, T.; Sakamoto, T.; Nonaka, G.I. Seasonal variation of mogrosides in lo han kuo (Siraitia Grosvenori) fruits. J. Nat. Med.; 2007; 61, pp. 307-312. [DOI: https://dx.doi.org/10.1007/s11418-006-0130-7]

14. Wang, B.; Yang, Z.; Xin, Z.; Ma, G.; Qian, Y.; Xie, T.; Prakash, I. Analysis of mogrosides in Siraitia Grosvenorii fruits at different stages of maturity. Nat. Prod. Commun.; 2019; 14, pp. 1-5. [DOI: https://dx.doi.org/10.1177/1934578X19878621]

15. Itkin, M.; Davidovich-Rikanati, R.; Cohen, S.; Portnoy, V.; Doron-Faigenboim, A.; Oren, E.; Freilich, S.; Tzuri, G.; Baranes, N.; Shen, S. et al. The biosynthetic pathway of the nonsugar, high-intensity sweetener mogroside V from Siraitia Grosvenorii. Proc. Natl. Acad. Sci. USA; 2016; 113, pp. E7619-E7628. [DOI: https://dx.doi.org/10.1073/pnas.1604828113] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27821754]

16. Srinivasan, P.; Smolke, C.D. Engineering a microbial biosynthesis platform for de novo production of tropane alkaloids. Nat. Commun.; 2019; 10, 3634. [DOI: https://dx.doi.org/10.1038/s41467-019-11588-w] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31406117]

17. Cravens, A.; Payne, J.; Smolke, C.D. Synthetic biology strategies for microbial biosynthesis of plant natural products. Nat. Commun.; 2019; 10, 2142. [DOI: https://dx.doi.org/10.1038/s41467-019-09848-w] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31086174]

18. Zhu, X.; Liu, X.; Liu, T.; Wang, Y.; Ahmed, N.; Li, Z.; Jiang, H. Synthetic biology of plant natural products: From pathway elucidation to engineered biosynthesis in plant cells. Plant Commun.; 2021; 2, 100229. [DOI: https://dx.doi.org/10.1016/j.xplc.2021.100229]

19. Huo, L.; Hug, J.J.; Fu, C.; Bian, X.; Zhang, Y.; Müller, R. Heterologous expression of bacterial natural product biosynthetic pathways. Nat. Prod. Rep.; 2019; 36, pp. 1412-1436. [DOI: https://dx.doi.org/10.1039/C8NP00091C]

20. Xia, M.; Han, X.; He, H.; Yu, R.; Zhen, G.; Jia, X.; Cheng, B.; Deng, X.W. Improved de novo genome assembly and analysis of the Chinese cucurbit Siraitia Grosvenorii, also known as monk fruit or Luo-Han-Guo. Gigascience; 2018; 7, giy067. [DOI: https://dx.doi.org/10.1093/gigascience/giy067]

21. Yang, D.; Park, S.Y.; Park, Y.S.; Eun, H.; Lee, S.Y. Metabolic engineering of Escherichia Coli for natural product biosynthesis. Trends Biotechnol.; 2020; 38, pp. 745-765. [DOI: https://dx.doi.org/10.1016/j.tibtech.2019.11.007]

22. Lane, S.; Zhang, Y.; Yun, E.J.; Ziolkowski, L.; Zhang, G.; Jin, Y.S.; Avalos, J.L. Xylose assimilation enhances the production of isobutanol in engineered Saccharomyces Cerevisiae. Biotechnol. Bioeng.; 2020; 117, pp. 372-381. [DOI: https://dx.doi.org/10.1002/bit.27202]

23. Chiu, C.H.; Wang, R.; Lee, C.C.; Lo, Y.C.; Lu, T.J. Biotransformation of mogrosides from Siraitia Grosvenorii Swingle by Saccharomyces Cerevisiae. J. Agric. Food Chem.; 2013; 61, pp. 7127-7134. [DOI: https://dx.doi.org/10.1021/jf402058p]

24. Wang, R.; Chen, Y.C.; Lai, Y.J.; Lu, T.J.; Huang, S.T.; Lo, Y.C. Dekkera Bruxellensis, a beer yeast that specifically bioconverts mogroside extracts into the intense natural sweetener siamenoside I. Food Chem.; 2019; 276, pp. 43-49. [DOI: https://dx.doi.org/10.1016/j.foodchem.2018.09.163]

25. Li, J.; Mutanda, I.; Wang, K.; Yang, L.; Wang, J.; Wang, Y. Chloroplastic metabolic engineering coupled with isoprenoid pool enhancement for committed taxanes biosynthesis in Nicotiana Benthamiana. Nat. Commun.; 2019; 10, 4850. [DOI: https://dx.doi.org/10.1038/s41467-019-12879-y] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31649252]

26. Farhi, M.; Marhevka, E.; Ben-Ari, J.; Algamas-Dimantov, A.; Liang, Z.; Zeevi, V.; Edelbaum, O.; Spitzer-Rimon, B.; Abeliovich, H.; Schwartz, B. et al. Generation of the potent anti-malarial drug artemisinin in tobacco. Nat. Biotechnol.; 2011; 29, pp. 1072-1074. [DOI: https://dx.doi.org/10.1038/nbt.2054] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/22158354]

27. Zhang, Y.; Nowak, G.; Reed, D.W.; Covello, P.S. The production of artemisinin precursors in tobacco. Plant Biotechnol. J.; 2011; 9, pp. 445-454. [DOI: https://dx.doi.org/10.1111/j.1467-7652.2010.00556.x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/20723135]

28. Bates, P.D.; Browse, J. The pathway of triacylglycerol synthesis through phosphatidylcholine in Arabidopsis produces a bottleneck for the accumulation of unusual fatty acids in transgenic seeds. Plant J.; 2011; 68, pp. 387-399. [DOI: https://dx.doi.org/10.1111/j.1365-313X.2011.04693.x]

29. Kempinski, C.; Chappell, J. Engineering triterpene metabolism in the oilseed of Arabidopsis Thaliana. Plant Biotechnol. J.; 2019; 17, pp. 386-396. [DOI: https://dx.doi.org/10.1111/pbi.12984]

30. Robert, S.S.; Singh, S.P.; Zhou, X.R.; Petrie, J.R.; Blackburn, S.I.; Mansour, P.M.; Nichols, P.D.; Liu, Q.; Green, A.G. Metabolic engineering of Arabidopsis to produce nutritionally important DHA in seed oil. Funct. Plant Biol.; 2005; 32, pp. 473-479. [DOI: https://dx.doi.org/10.1071/FP05084]

31. Chen, Q.J.; Zhou, H.M.; Chen, J.; Wang, X.C. A Gateway-Based platform for multigene plant transformation. Plant Mol. Biol.; 2006; 62, pp. 927-936. [DOI: https://dx.doi.org/10.1007/s11103-006-9065-3]

32. Gibson, D.G.; Benders, G.A.; Axelrod, K.C.; Zaveri, J.; Algire, M.A.; Moodie, M.; Montague, M.G.; Venter, J.C.; Smith, H.O.; Hutchison, C.A. One-step assembly in yeast of 25 overlapping DNA fragments to form a complete synthetic Mycoplasma genitalium genome. Proc. Natl. Acad. Sci. USA; 2008; 105, pp. 20404-20409. [DOI: https://dx.doi.org/10.1073/pnas.0811011106]

33. Sleight, S.C.; Bartley, B.A.; Lieviant, J.A.; Sauro, H.M. In-Fusion biobrick assembly and re-engineering. Nucleic Acids. Res.; 2010; 38, pp. 2624-2636. [DOI: https://dx.doi.org/10.1093/nar/gkq179]

34. Zhu, Q.; Yu, S.; Zeng, D.; Liu, H.; Wang, H.; Yang, Z.; Xie, X.; Shen, R.; Tan, J.; Li, H. et al. Development of “Purple Endosperm Rice” by engineering anthocyanin biosynthesis in the endosperm with a high-efficiency transgene stacking system. Mol. Plant; 2017; 10, pp. 918-929. [DOI: https://dx.doi.org/10.1016/j.molp.2017.05.008] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/28666688]

35. Zhu, Q.; Zeng, D.; Yu, S.; Cui, C.; Li, J.; Li, H.; Chen, J.; Zhang, R.; Zhao, X.; Chen, L. et al. From golden rice to ASTARice: Bioengineering astaxanthin biosynthesis in rice endosperm. Mol. Plant; 2018; 11, pp. 1440-1448. [DOI: https://dx.doi.org/10.1016/j.molp.2018.09.007] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/30296601]

36. Zhao, Z.; Xie, X.; Liu, W.; Huang, J.; Tan, J.; Yu, H.; Zong, W.; Tang, J.; Zhao, Y.; Xue, Y. et al. STI PCR: An efficient method for amplification and de novo synthesis of long DNA sequences. Mol. Plant; 2022; 15, pp. 620-629. [DOI: https://dx.doi.org/10.1016/j.molp.2021.12.018]

37. Douin, V.; Bornes, S.; Creancier, L.; Rochaix, P.; Favre, G.; Prats, A.C.; Couderc, B. Use and comparison of different internal ribosomal entry sites (IRES) in tricistronic retroviral vectors. BMC Biotechnol.; 2004; 4, 16. [DOI: https://dx.doi.org/10.1186/1472-6750-4-16] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/15279677]

38. He, Y.; Zhang, T.; Sun, H.; Zhan, H.; Zhao, Y. A reporter for noninvasively monitoring gene expression and plant transformation. Hortic. Res.; 2020; 7, pp. 1-6. [DOI: https://dx.doi.org/10.1038/s41438-020-00390-1]

39. Liu, X.; Yang, W.; Mu, B.; Li, S.; Li, Y.; Zhou, X.; Zhang, C.; Fan, Y.; Chen, R. Engineering of ‘Purple Embryo Maize’ with a multigene expression system derived from a bidirectional promoter and self-cleaving 2A peptides. Plant Biotechnol. J.; 2018; 16, pp. 1107-1109. [DOI: https://dx.doi.org/10.1111/pbi.12883]

40. Paine, J.A.; Shipton, C.A.; Chaggar, S.; Howells, R.M.; Kennedy, M.J.; Vernon, G.; Wright, S.Y.; Hinchliffe, E.; Adams, J.L.; Silverstone, A.L. et al. Improving the nutritional value of golden rice through increased pro-vitamin A content. Nat. Biotechnol.; 2005; 23, pp. 482-487. [DOI: https://dx.doi.org/10.1038/nbt1082]

41. Zhu, C.; Naqvi, S.; Breitenbach, J.; Sandmann, G.; Christou, P.; Capell, T. Combinatorial genetic transformation generates a library of metabolic phenotypes for the carotenoid pathway in maize. Proc. Natl. Acad. Sci. USA; 2008; 105, pp. 18232-18237. [DOI: https://dx.doi.org/10.1073/pnas.0809737105] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/19011084]

42. Shewmaker, C.K.; Sheehy, J.A.; Daley, M.; Colburn, S.; Ke, D.Y. Seed-specific overexpression of phytoene synthase: Increase in carotenoids and other metabolic effects. Plant J.; 1999; 20, pp. 401-412. [DOI: https://dx.doi.org/10.1046/j.1365-313x.1999.00611.x]

43. Kim, M.J.; Kim, J.K.; Kim, H.J.; Pak, J.H.; Lee, J.H.; Kim, D.H.; Choi, H.K.; Jung, H.W.; Lee, J.D.; Chung, Y.S. et al. Genetic modification of the soybean to enhance the β-carotene content through seed-specific expression. PLoS ONE; 2012; 7, e48287. [DOI: https://dx.doi.org/10.1371/journal.pone.0048287]

44. Butelli, E.; Titta, L.; Giorgio, M.; Mock, H.P.; Matros, A.; Peterek, S.; Schijlen, E.G.W.M.; Hall, R.D.; Bovy, A.G.; Luo, J. et al. Enrichment of tomato fruit with health-promoting anthocyanins by expression of select transcription factors. Nat. Biotechnol.; 2008; 26, pp. 1301-1308. [DOI: https://dx.doi.org/10.1038/nbt.1506] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/18953354]

45. Han, J.Y.; Baek, S.H.; Jo, H.J.; Yun, D.W.; Choi, Y.E. Genetically modified rice produces ginsenoside aglycone (protopanaxadiol). Planta; 2019; 250, pp. 1103-1110. [DOI: https://dx.doi.org/10.1007/s00425-019-03204-4]

46. Huang, J.C.; Zhong, Y.J.; Liu, J.; Sandmann, G.; Chen, F. Metabolic engineering of tomato for high-yield production of astaxanthin. Meta Eng.; 2013; 17, pp. 59-67. [DOI: https://dx.doi.org/10.1016/j.ymben.2013.02.005] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/23511430]

47. Khairul Ikram, N.K.B.; Beyraghdar Kashkooli, A.; Peramuna, A.V.; van der Krol, A.R.; Bouwmeester, H.; Simonsen, H.T. Stable production of the antimalarial drug artemisinin in the moss Physcomitrella Patens. Front. Bioeng. Biotech.; 2017; 5, pp. 1-8. [DOI: https://dx.doi.org/10.3389/fbioe.2017.00047]

48. Holst, J.; Vignali, K.M.; Burton, A.R.; Vignali, D.A.A. Rapid analysis of T-cell selection in vivo using T cell-receptor retrogenic mice. Nat. Methods; 2006; 3, pp. 191-197. [DOI: https://dx.doi.org/10.1038/nmeth858] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/16489336]

49. Ikeda, Y.; Pélissier, T.; Bourguet, P.; Becker, C.; Pouch-Pélissier, M.N.; Pogorelcnik, R.; Weingartner, M.; Weigel, D.; Deragon, J.M.; Mathieu, O. Arabidopsis proteins with a transposon-related domain act in gene silencing. Nat. Commun.; 2017; 8, 15122. [DOI: https://dx.doi.org/10.1038/ncomms15122]

50. Ee, S.F.; Khairunnisa, M.B.; Zeti-Azura, M.H.; Noor Azmi, S.; Zamri, Z. Effective hygromycin concentration for selection of Agrobacterium-mediated transgenic Arabidopsis Thaliana. Malays. Appl. Biol.; 2014; 43, pp. 119-123.

51. Kim, N.C.; Kinghorn, A.D. Sweet-tasting and sweetness-modifying constituents of plants. Stud. Nat. Prod. Chem.; 2002; 27, pp. 3-57. [DOI: https://dx.doi.org/10.1016/S1572-5995(02)80033-3]

52. Suzuki, Y.A.; Murata, Y.; Inui, H.; Sugiura, M.; Nakano, Y. Triterpene glycosides of Siraitia Grosvenori inhibit rat intestinal maltase and suppress the rise in blood glucose level after a single oral administration of maltose in rats. J. Agric. Food Chem.; 2005; 53, pp. 2941-2946. [DOI: https://dx.doi.org/10.1021/jf0478105]

53. Wang, R.; Lin, P.Y.; Huang, S.T.; Chiu, C.H.; Lu, T.J.; Lo, Y.C. Hyperproduction of β-glucanase Exg1 promotes the bioconversion of mogrosides in Saccharomyces Cerevisiae mutants defective in mannoprotein deposition. J. Agric. Food Chem.; 2015; 63, pp. 10271-10279. [DOI: https://dx.doi.org/10.1021/acs.jafc.5b03909]

54. Zhu, Q.; Wang, B.; Tan, J.; Liu, T.; Li, L.; Liu, Y.G. Plant synthetic metabolic engineering for enhancing crop nutritional quality. Plant Commun.; 2020; 1, 100017. [DOI: https://dx.doi.org/10.1016/j.xplc.2019.100017] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/33404538]

55. Dafny-Yelin, M.; Chung, S.M.; Frankman, E.L.; Tzfira, T. PSAT RNA interference vectors: A modular series for multiple gene down-regulation in plants. Plant Physiol.; 2007; 145, pp. 1272-1281. [DOI: https://dx.doi.org/10.1104/pp.107.106062] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/17766396]

56. Halpin, C. Gene stacking in transgenic plants-The challenge for 21st century plant biotechnology. Plant Biotechnol. J.; 2005; 3, pp. 141-155. [DOI: https://dx.doi.org/10.1111/j.1467-7652.2004.00113.x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/17173615]

57. De La Garza, R.I.D.; Gregory, J.F.; Hanson, A.D. Folate biofortification of tomato fruit. Proc. Natl. Acad. Sci. USA; 2007; 104, pp. 4218-4222. [DOI: https://dx.doi.org/10.1073/pnas.0700409104] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/17360503]

58. Nazhand, A.; Durazzo, A.; Lucarini, M.; Mobilia, M.A.; Omri, B.; Santini, A. Rewiring cellular metabolism for heterologous biosynthesis of taxol. Nat. Prod. Res.; 2020; 34, pp. 110-121. [DOI: https://dx.doi.org/10.1080/14786419.2019.1630122] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/31298589]

59. Breitel, D.; Brett, P.; Alseekh, S.; Fernie, A.R.; Butelli, E.; Martin, C. Metabolic engineering of tomato fruit enriched in L-DOPA. Meta Eng.; 2021; 65, pp. 185-196. [DOI: https://dx.doi.org/10.1016/j.ymben.2020.11.011]

60. Grützner, R.; Schubert, R.; Horn, C.; Yang, C.; Vogt, T.; Marillonnet, S. Engineering betalain biosynthesis in tomato for high level betanin production in fruits. Front. Plant Sci.; 2021; 12, 682443. [DOI: https://dx.doi.org/10.3389/fpls.2021.682443]

61. Gutensohn, M.; Dudareva, N. Tomato Fruits—A Platform for Metabolic Engineering of Terpenes. Methods Enzymol.; 2016; 576, pp. 333-359. [DOI: https://dx.doi.org/10.1016/bs.mie.2016.03.012]

62. Firsov, A.; Mitiouchkina, T.; Shaloiko, L.; Pushin, A.; Vainstein, A.; Dolgov, S. Agrobacterium-mediated transformation of chrysanthemum with artemisinin biosynthesis pathway genes. Plants; 2020; 9, 537. [DOI: https://dx.doi.org/10.3390/plants9040537]

63. Gou, Y.; Zhang, F.; Tang, Y.; Jiang, C.; Bai, G.; Xie, H.; Chen, M.; Liao, Z. Engineering nootkatone biosynthesis in Artemisia Annua. ACS Synth. Biol.; 2021; 10, pp. 957-963. [DOI: https://dx.doi.org/10.1021/acssynbio.1c00016]

64. Deng, B.; Zhang, P.; Ge, F.; Liu, D.Q.; Chen, C.Y. Enhancement of triterpenoid saponins biosynthesis in Panax Notoginseng cells by co-overexpressions of 3-Hydroxy-3-Methylglutaryl CoA reductase and squalene synthase genes. Biochem. Eng. J.; 2017; 122, pp. 38-46. [DOI: https://dx.doi.org/10.1016/j.bej.2017.03.001]

65. Matsumoto, K.; Kasai, R.; Ohtani, K.; Tanaka, O. Minor cucurbitane-glycosides from fruits of Siraitia Grosvenori (Cucurbitaceae). Chem. Pharm. Bull.; 1990; 38, pp. 2030-2032. [DOI: https://dx.doi.org/10.1248/cpb.38.2030]

66. Clough, S.J.; Bent, A.F. Floral Dip: A simplified method for Agrobacterium-mediated transformation of Arabidopsis Thaliana. Plant J.; 1998; 16, pp. 735-743. [DOI: https://dx.doi.org/10.1046/j.1365-313x.1998.00343.x] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/10069079]

67. Luo, Z.; Shi, H.; Zhang, K.; Qin, X.; Guo, Y.; Ma, X. Liquid chromatography with tandem mass spectrometry method for the simultaneous determination of multiple sweet mogrosides in the fruits of Siraitia Grosvenorii and its marketed sweeteners. J. Sep. Sci.; 2016; 39, pp. 4124-4135. [DOI: https://dx.doi.org/10.1002/jssc.201600563] [PubMed: https://www.ncbi.nlm.nih.gov/pubmed/27593208]