Genome-wide analysis of the C2H2-ZFP gene family in Stevia rebaudiana reveals involvement in abiotic stress response.

Stevia (Stevia rebaudiana Bertoni) is a natural sweetener plant that accumulates highly sweet steviol glycosides (SGs) especially in leaves. Stevia is native to humid areas and does not have a high tolerance to drought which is the most serious abiotic stress restricting its production worldwide. C2H2 zinc finger proteins (C2H2-ZFPs) are a group of well-known transcription factors that involves in various developmental, physiological and biochemical activities as well as in response to abiotic stresses. Here we analyzed C2H2-ZFP gene family in stevia and identified a total of 185 putative SrC2H2-ZF proteins from the genome sequence of S. rebaudiana. We further characterized the identified C2H2-ZF domains and their organization, additional domains and motifs and analyzed their physicochemical properties, localization and gene expression patterns. The cis-element analysis suggested multiple roles of SrC2H2-ZFPs in response to light, phytohormone, and abiotic stresses. In silico analysis revealed that the stevia C2H2-ZFP genes are interactively expressed in different tissues and developmental stages and some C2H2-ZFP genes are involved in response to drought stress. This study provides a background for future exploration of the functional, and regulatory aspects of the C2H2-ZFP gene family in S. rebaudiana.

Transcriptomic Analyses Reveal Insights into the Shared Regulatory Network of Phenolic Compounds and Steviol Glycosides in Stevia rebaudiana.

Stevia rebaudiana (Bertoni) is a highly valuable crop for the steviol glycoside content in its leaves, which are no-calorie sweeteners hundreds of times more potent than sucrose. The presence of health-promoting phenolic compounds, particularly flavonoids, in the leaf of S. rebaudiana adds further nutritional value to this crop. Although all these secondary metabolites are highly desirable in S. rebaudiana leaves, the genes regulating the biosynthesis of phenolic compounds and the shared gene network between the regulation of biosynthesis of steviol glycosides and phenolic compounds still need to be investigated in this species. To identify putative candidate genes involved in the synergistic regulation of steviol glycosides and phenolic compounds, four genotypes with different contents of these compounds were selected for a pairwise comparison RNA-seq analysis, yielding 1136 differentially expressed genes. Genes that highly correlate with both steviol glycosides and phenolic compound accumulation in the four genotypes of S. rebaudiana were identified using the weighted gene co-expression network analysis. The presence of UDP-glycosyltransferases 76G1, 76H1, 85C1, and 91A1, and several genes associated with the phenylpropanoid pathway, including peroxidase, caffeoyl-CoA O-methyltransferase, and malonyl-coenzyme A:anthocyanin 3-O-glucoside-6″-O-malonyltransferase, along with 21 transcription factors like SCL3, WRK11, and MYB111, implied an extensive and synergistic regulatory network involved in enhancing the production of such compounds in S. rebaudiana leaves. In conclusion, this work identified a variety of putative candidate genes involved in the biosynthesis and regulation of particular steviol glycosides and phenolic compounds that will be useful in gene editing strategies for increasing and steering the production of such compounds in S. rebaudiana as well as in other species.

In vitro propagation of Indonesian stevia (Stevia rebaudiana) genotype using axenic nodal segments.

OBJECTIVE: The high industrial demand for Stevia cultivation (Stevia rebaudiana) has increased due to its high stevioside content derived from the leaves. However, the low germination rate makes the cultivation of the plant become the main obstacle. Therefore, an efficient cultivation technique is required. This present work aims to analyze the effect of five combinations of Kinetin (Kin) and benzyladenine (BA) on stevia micropropagation using nodal segment explants. RESULTS: The micropropagation of stevia was performed using Murashige and Skoog (MS) medium supplemented with BA and Kin. We analyzed different organogenesis and callogenesis responses. In addition, the number of shoots and root formed during in vitro culture were also observed. Our results demonstrated that all treatments with Kin, both alone and in combination with BA, resulted in the development of callus on all nodal segment explants. Explants treated in MS with 1 mg L-1 BA exhibited the best average of shoot number (36.27). In contrast, the treatment without PGR resulted in the best root formation (2.6). The overall results suggested that different combination of BA and Kin resulted in distinct organogenesis responses, where 1 mg L-1 of BA was potentially used for boosting the number of shoots in micropropagation of stevia accession Mini.

Genome-wide transcriptional profiling and physiological investigation elucidating the molecular mechanism of multiple abiotic stress response in Stevia rebaudiana Bertoni.

Considering the major source of plant-derived low/non-calorie steviol glycosides (SGs), comprehensive physiological, biochemical, and deep transcriptional investigations were conducted to explicit deeper insight into multiple abiotic stress responses in Stevia rebaudiana. The physiological indicators including photosynthesis, chlorophyll, relative water content, shoot growth, electrolyte leakage, and SG biosynthesis were negatively impacted under drought (DS), followed by salinity (SS) and waterlogging (WS). Global transcriptional analysis revealed significant upregulated expression of the genes encoding for ROS detoxification (GST, SOD, APX, glutathione peroxidase), osmotic adjustment (alpha-trehalose-phosphate and S-adenosylmethionine decarboxylase), ion transporters (CAX, NHX, CNGS, VPPase, VATPase), water channel (PIP1, TIP) and abiotic stress-responsive candidate genes (LEA, HSPs, and Dehydrins) regulating abiotic stress response in S. rebaudiana. These inferences were complemented with predicted interactome network that revealed regulation of energy metabolism by key stress-responsive genes (GST, HKT1, MAPKs, P5CSs, PIP), transcription factors (HSFA2, DREB1A, DREB2A), and abiotic stress responsive pathways (ABA, ethylene, ion stress). This is the first detailed study to comprehend the molecular regulation of stress response and their interplay under DS, SS, and WS. The key genes and regulators can be functionally validated, and will facilitate targeted gene editing for genetic improvement of crop sustainability under changing environmental conditions in S. rebaudiana.

Genome-Wide Identification and Expression Analysis of bZIP Family Genes in Stevia rebaudiana.

The basic (region) leucine zippers (bZIPs) are evolutionarily conserved transcription factors widely distributed in eukaryotic organisms. In plants, they are not only involved in growth and development, defense and stress responses and regulation of physiological processes but also play a pivotal role in regulating secondary metabolism. To explore the function related to the bZIP gene family in Stevia rebaudiana Bertoni, we identified 105 SrbZIP genes at the genome-wide level and classified them into 12 subfamilies using bioinformation methods. Three main classes of cis-acting elements were found in the SrbZIP promoter regions, including development-related elements, defense and stress-responsive elements and phytohormone-responsive elements. Through protein-protein interaction network of 105 SrbZIP proteins, SrbZIP proteins were mainly classified into four major categories: ABF2/ABF4/ABI5 (SrbZIP51/SrbZIP38/SrbZIP7), involved in phytohormone signaling, GBF1/GBF3/GBF4 (SrbZIP29/SrbZIP63/SrbZIP60) involved in environmental signaling, AREB3 (SrbZIP88), PAN (SrbZIP12), TGA1 (SrbZIP69), TGA4 (SrbZIP82), TGA7 (SrbZIP31), TGA9 (SrbZIP95), TGA10 (SrbZIP79) and HY5 (SrbZIP96) involved in cryptochrome signaling, and FD (SrbZIP72) promoted flowering. The transcriptomic data showed that SrbZIP genes were differentially expressed in six S. rebaudiana cultivars ('023', '110', 'B1188', '11-14', 'GP' and 'GX'). Moreover, the expression levels of selected 15 SrbZIP genes in response to light, abiotic stress (low temperature, salt and drought), phytohormones (methyl jasmonate, gibberellic acid and salicylic acid) treatment and in different tissues were analyzed utilizing qRT-PCR. Some SrbZIP genes were further identified to be highly induced by factors affecting glycoside synthesis. Among them, three SrbZIP genes (SrbZIP54, SrbZIP63 and SrbZIP32) were predicted to be related to stress-responsive terpenoid synthesis in S. rebaudiana. The protein-protein interaction network expanded the potential functions of SrbZIP genes. This study firstly provided the comprehensive genome-wide report of the SrbZIP gene family, laying a foundation for further research on the evolution, function and regulatory role of the bZIP gene family in terpenoid synthesis in S. rebaudiana.

Genome-wide identification of SrbHLH transcription factors highlights its potential role in rebaudioside A (RA) biosynthesis in Stevia rebaudiana.

Stevia rebaudiana Bertoni is a valuable medicinal plant and an essential source of natural sweetener, steviol glycosides (SGs), with rebaudioside A (RA) being one of the main components of SGs. bHLH transcription factors play a crucial role in plant development and secondary metabolism. In this study, 159 SrbHLH genes were identified from the S. rebaudiana genome, and each gene was named based on its chromosome location. The SrbHLH proteins were then clustered into 18 subfamilies through phylogenetic analysis. The analysis of conserved motifs and gene structure further supported the classification of the SrbHLH family. Chromosomal location and gene duplication events of SrbHLH genes were also studied. Moreover, based on the RNA-Seq data of different tissues of S. rebaudiana, 28 SrbHLHs were co-expressed with structural genes involved in RA biosynthesis. The expression pattern of candidate SrbHLH genes were confirmed by qPCR. Finally, dual luciferase reporter assays (DLAs) and subcellular localization analysis verified SrbHLH22, SrbHLH111, SrbHLH126, SrbHLH142, and SrbHLH152 are critical regulators of RA biosynthesis. This study provides new insights into the function of SrbHLHs in regulating SGs biosynthesis and lays the foundation for future applications of SrbHLH genes in molecular breeding of S. rebaudiana.

Enhancement of diterpenoid steviol glycosides by co-overexpressing SrKO and SrUGT76G1 genes in Stevia rebaudiana Bertoni.

Stevia rebaudiana (stevia) contains commercially important steviol glycosides, stevioside and rebaudioside A, these compounds have insulinotropic and anti-hyperglycemic effect. Steviol, stevioside and rebaudioside-A have taste modulation and insulin potentiation activity. Stevia leaves are composed of steviol (2-5%), stevioside (4-13%) and rebaudioside-A (1-6%). Stevioside has after-taste bitterness, rebaudioside-A is sweetest in taste among all the glycosides present. Therefore, lower ratio of rebaudioside-A to stevioside has bitter after-taste, which makes stevia plants unpalatable. By over-expressing the genes, SrUGT76G1 and SrKO, we propose to increase the ratio of RebA to stevioside in stevia. Various lines were generated and amongst them, seven lines had both the transgenes present. Co-overxpresion of SrUGT76G1 and SrKO led to the increased concentration of RebA in all the seven transgenic lines (KU1-KU7) than control plant and RebA to stevioside ratio also increased significantly. Steviol, stevioside and RebA showed a differential concentration in all the seven lines, but the pattern was the same in all of them and the ratio of RebA to stevioside increased dramatically. In transgenic line 2 (KU2), RebA showed a steep increase in concentration 52% the rebaudioside-A to stevioside ratio increased from 0.74 (control) to 2.83. In overall all the lines, RebA showed a positive correlation with steviol and stevioside. Overexpression of SrKO led to an increase in steviol which increased the stevioside, overexpression of SrUGT76G1 ultimately increased RebA concentration. In conclusion, concentration of RebA increased significantly with co- overexpression of SrUGT6G1 and SrKO genes. Lines with increased RebA are more palatable and commercially viable.

Seeds of Stevia rebaudiana Bertoni as a Source of Plant Growth-Promoting Endophytic Bacteria with the Potential to Synthesize Rebaudioside A.

In this study, a new strain of Pantoea vagans, SRS89, was isolated from surface-sterilized stevia seeds. The isolate was evaluated using morphological, molecular, and biochemical methods. The bacterium was 1.5 µm long, yellowish in color, and classified as Gram-negative. Whole genome sequencing of our strain revealed the presence of a 4,610,019 bp chromosome, and genome annotation resulted in the detection of 4283 genes encoding 4204 putative coding sequences. Phylogenic analysis classified the genome of our strain close to the MP7 and LMG 24199 strains of P. vagans. Functional analysis showed that the highest number of genes within the analyzed bacterium genome were involved in transcription, amino acid transport and metabolism, and carbohydrate transport and metabolism. We also identified genes for enzymes involved in the biosynthesis of carotenoids and terpenoids. Furthermore, we showed the presence of growth regulators, with the highest amount noted for gibberellic acid A3, indole-3-acetic acid, and benzoic acid. However, the most promising property of this strain is its ability to synthesize rebaudioside A; the estimated amount quantified using reversed-phase (RP)-HPLC was 4.39 mg/g of the dry weight of the bacteria culture. The isolated endophytic bacterium may be an interesting new approach to the production of this valuable metabolite.

Multi-walled carbon nanotubes interact with light intensity to affect morpho-biochemical, nutrient uptake, DNA damage, and secondary metabolism of Stevia rebaudiana.

In this study, the interaction between nanoparticles (0, 50, 100, and 150 mg L-1) and light intensity (100, 200, and 400 µmol·m-2·s-1) was evaluated for effectiveness in improving stevia shoot induction by measuring morphological traits, nutrient absorption, total carbohydrates, steviol glycosides (SVglys), and DNA damage in two DNA sequence regions (promoter and sequence of the UGT76G1 gene). MWCNTs at a concentration of 50 mg L-1 in interaction with the light intensity of 200 µmol·m-2·s-1 improved the morphological traits and absorption of nutrients such as nitrogen (N), phosphorous (P), potassium (K), calcium (Ca), iron (Fe), and Manganese (Mn), compared to other treatments. Also, under this interaction, the accumulation of total carbohydrates and SVglys was elevated. Moreover, DNA damage in both regions of the DNA sequence under light intensity at low concentrations of MWCNTs (0 and 50 mg L-1) did not show a significant change but increased with increasing MWCNT concentration at high light intensities (200 and 400 µmol·m-2·s-1). These results demonstrate that the advantages and phytotoxicity of MWCNTs in the in vitro culture of stevia are dose-dependent and are affected by light intensity. Based on this, the interaction of 50 mg L-1 of MWCNTs with the light intensity of 200 µmol·m-2·s-1 is recommended to improve stevia micropropagation and subsequent growth and metabolism.

The effect of alginate as an elicitor on transcription of steviol glycosides biosynthesis pathway related key genes and sweeteners content in in vitro cultured Stevia rebaudiana.

BACKGROUND: Stevia rebaudiana is a medicinal herb that accumulates non-caloric sweeteners called steviol glycosides (SGs) which are approximately 300 times sweeter than sucrose. This study used alginate (ALG) as an elicitor to increase steviol glycosides accumulation and elucidate gene transcription in the steviol glycosides biosynthesis pathway. METHODS AND RESULTS: To minimize the grassy taste associated with stevia sweeteners, plantlets were grown in complete darkness. ALG was applied to stevia plants grown in suspension culture with a Murashige and Skoog (MS) medium to determine its effect on SGs' content and the transcription profile of SG-related genes using the HPLC and RT-qPCR methods, respectively. Treatment with alginate did not significantly affect plantlet growth parameters such as shoot number, dry and fresh weight. Rebaudioside A (Reb A) content increased approximately sixfold in the presence of 1g L-1 alginate and KS, KAH, and UGT74G1 genes showed significant up-regulation. When the concentration was increased to 2g L-1, the transcription of KO and UGT76G1, responsible for the conversion of stevioside to Reb A, was increased about twofold. CONCLUSIONS: The current study proposes that adding alginate to the MS suspension medium can increase Reb A levels by altering the SG biosynthesize pathway's transcription profile. The present experiment provides new insights into the biochemical and transcriptional response mechanisms of suspension-cultured stevia plants to alginate.

An interdisciplinary approach towards sustainable and higher steviol glycoside production from in vitro cultures of Stevia rebaudiana.

Stevia rebaudiana is one of the vastly acclaimed commercial plant in the world and belongs to Asteraceae family. The exclusive advantage of Stevia over artificial sweeteners is impeccable and targets its potentiality to the presence of diterpene glycosides. Moreover, the flaunting sweetness of steviol glycosides with associated medicinal benefits, turns the plant to be one of the most economic assets, globally. As compared to vegetative propagation through stem-cuttings, plant tissue culture is the most suitable approach in obtaining true-to-type plants of superior quality. During last few decades, significant in vitro propagation methods have been developed and still the research is ongoing. The present review discusses the tissue culture perspectives of S. rebaudiana, primarily focusing on the mineral nutrition, growth regulators and other accessory factors, motioning the optimum growth and development of the plant. Another crucial aspect is the generation of sweeter varieties in order to reduce the bitter-off taste, which is noticed after the consumption of the leaves. The in vitro cultures pose an efficient alternative system for production of steviol glycosides, with higher rebaudioside(s) content. Moreover, the review also covers the recent approaches pertaining to scale-up studies and genome editing perspectives.

Enhanced specialized metabolite, trichome density, and biosynthetic gene expression in Stevia rebaudiana (Bertoni) Bertoni plants inoculated with endophytic bacteria Enterobacter hormaechei.

Stevia rebaudiana (Bertoni) Bertoni is a plant of economic interest in the food and pharmaceutical industries due its steviol glycosides (SG), which are rich in metabolites that are 300 times sweeter than sucrose. In addition, S. rebaudiana plants contain phenolic compounds and flavonoids with antioxidant activity. Endophytic bacteria promote the growth and development and modulate the metabolism of the host plant. However, little is known regarding the role of endophytic bacteria in the growth; synthesis of SG, flavonoids and phenolic compounds; and the relationship between trichome development and specialized metabolites in S. rebaudiana, which was the subject of this study. The 12 bacteria tested did not increase the growth of S. rebaudiana plants; however, the content of SG increased with inoculation with the bacteria Enterobacter hormaechei H2A3 and E. hormaechei H5A2. The SG content in leaves paralleled an increase in the density of glandular, short, and large trichome. The image analysis of S. rebaudiana leaves showed the presence of SG, phenolic compounds, and flavonoids principally in glandular and short trichomes. The increase in the transcript levels of the KO, KAH, UGT74G1, and UGT76G1 genes was related to the SG concentration in plants of S. rebaudiana inoculated with E. hormaechei H2A3 and E. hormaechei H5A2. In conclusion, inoculation with the stimulating endophytes E. hormaechei H2A3 and E. hormaechei H5A2 increased SG synthesis, flavonoid content and flavonoid accumulation in the trichomes of S. rebaudiana plants.

Lelliottia steviae sp. nov. isolated from Stevia rebaudiana Bertoni.

A Gram-negative, aerobic, chemoheterotrophic, rod-shaped, and motile bacterium, designated as LST-1T, was isolated from wild Stevia rebaudiana Bertoni and subjected to a polyphasic taxonomic analysis. The LST-1 strain grew optimally at 37 °C and pH 6.0-7.0 in the presence of 0.5% (w/v) NaCl. Phylogenetic analysis based on the 16S rDNA sequence indicated that LST-1 is closely related to Lelliottia jeotgali PFL01T (99.85%), Lelliottia nimipressuralis LMG10245T (98.82%), and Lelliottia amnigena LMG2784T (98.54%). Multi-locus sequence typing of concatenated partial atpD, infB, gyrB, and rpoB genes was performed to improve the resolution, and clear distinctions between the closest related type strains were observed. The results of average nucleotide identify analyses and DNA-DNA hybridization with four species (16S rDNA similarity > 98.65%) were less than 90 and 40%, respectively, verifying the distinct characteristics from other species of Lelliottia. The cellular fatty acid profile of the strain consisted of C16:0, Summed Feature3, and Summed Feature8 (possibly 16:1 w6c/16:1 w7c and 18:1 w6c) as major components. The major polar lipids included phosphatidylethanolamine, phosphatidylglycerol, an aminophospholipid, three non-characteristic phospholipids, and a non-characteristic lipid. The genome of LST-1T was 4,611,055 bp in size, with a G + C content of 55.02%. The unique combination of several phenotypic, chemotaxonomic, and genomic characteristics proved that strain LST-1T belongs to a novel species, for which the name Lelliottia steviae sp. nov. is proposed. The type strain is LST-1T (= CGMCC 1.19175T = JCM 34938T).Repositories: The genbank accession numbers for the 16S rRNA gene and genome sequences of strain LST-1T are MZ497264 and CP063663, respectively.

Characterisation of LTR-Retrotransposons of Stevia rebaudiana and Their Use for the Analysis of Genetic Variability.

Stevia rebaudiana is one of the most important crops belonging to the Asteraceae family. Stevia is cultivated all over the world as it represents a valid natural alternative to artificial sweeteners thanks to its leaves, which produce steviol glycosides that have high sweetening power and reduced caloric value. In this work, the stevia genome sequence was used to isolate and characterise full-length long-terminal repeat retrotransposons (LTR-REs), which account for more than half of the genome. The Gypsy retrotransposons were twice as abundant as the Copia ones. A disproportionate abundance of elements belonging to the Chromovirus/Tekay lineage was observed among the Gypsy elements. Only the SIRE and Angela lineages represented significant portions of the genome among the Copia elements. The dynamics with which LTR-REs colonised the stevia genome were also estimated; all isolated full-length elements turned out to be relatively young, with a proliferation peak around 1-2 million years ago. However, a different analysis conducted by comparing sequences encoding retrotranscriptase showed the occurrence of an older period in which there was a lot of LTR-RE proliferation. Finally, a group of isolated full-length elements belonging to the lineage Angela was used to analyse the genetic variability in 25 accessions of S. rebaudiana using the Inter-Retrotransposon Amplified Polymorphism (IRAP) protocol. The obtained fingerprints highlighted a high degree of genetic variability and were used to study the genomic structures of the different accessions. It was hypothesised that there are four ancestral subpopulations at the root of the analysed accessions, which all turned out to be admixed. Overall, these data may be useful for genome sequence annotations and for evaluating genetic variability in this species, which may be useful in stevia breeding.

Enhancement of Rebaudioside M Production by Structure-Guided Engineering of Glycosyltransferase UGT76G1.

Owing to zero-calorie and advanced organoleptic properties similar to sucrose, the plant-derived rebaudioside M (Reb M) has been considered as a next generation sweetener. However, a low content of Reb M in Stevia rebaudiana Bertoni and low enzymatic activity of UGT76G1, which is an uridine diphosphate glucose (UDPG)-dependent glycosyltransferase with the ability to glycosylate rebaudioside D (Reb D) to produce Reb M through the formation of ß-1,3 glycosidic bond, restrict its commercial usage. To improve the catalytic activity of UGT76G1, a variant UGT76G1-T284S/M88L/L200A was obtained by structure-guided evolution, whose catalytic activity toward Reb D increased by 2.38 times compared with UGT76G1-T284S. This allowed us to prepare Reb M on a large-scale with a great yield of 90.50%. Moreover, molecular dynamics simulation illustrated that UGT76G1-T284S/M88L/L200A reduced distances from Reb D to catalytic residues and UDPG. Hence, we report an efficient method for the potential scale production of Reb M in this study.

Biotechnological interventions of in vitro propagation and production of valuable secondary metabolites in Stevia rebaudiana.

Plant cell and tissue culture makes provision of a sustainable and nature-friendly strategy for the production of secondary metabolites, and modern progress in gene editing and genome engineering provides novel possibilities to improve both the qualitative and quantitative aspects of such phytochemicals. The ever-expanding quest for plant-based medicine to treat diabetes facilitates large-scale cultivation of Stevia rebaudiana to enhance the yield of its much-coveted low-calorie sweetener glycosides. The potential to process stevia as a "natural" product should enhance the acceptance of steviosides as a natural calorie-free sweetener especially suitable for use in diabetic and weight control drinks and foods. Besides sweetener agents, S. rebaudiana is a potent source of many antioxidant compounds and is used to cure immunodeficiencies, neurologic disorders, inflammation, diabetes mellitus, Parkinson's disease, and Alzheimer's disease. This comprehensive review presents the research outcomes of the many biotechnological interventions implicated to upscale the yield of steviol glycosides and its derivatives in in vitro cell, callus, tissue, and organ cultures with notes on the use of bioreactor and genetic engineering in relation to the production of these valuable compounds in S. rebaudiana. KEY POINTS: • Critical and updated assessment on sustainable production of steviol glycosides from Stevia rebaudiana. • In vitro propagation of S. rebaudiana and elicitation of steviol glycosides production. • Genetic fidelity and diversity assessment of S. rebaudiana using molecular markers.

Comparative transcriptomic of Stevia rebaudiana provides insight into rebaudioside D and rebaudioside M biosynthesis.

Rebaudioside D (Reb D) and rebaudioside M (Reb M) are commercially important low/no-calorie natural sweeteners. However, they are present in a minor proportion of all steviol glycosides (SGs) in Stevia rebaudiana Bertoni (S. rebaudiana). Strain-dependent deviation in Reb D and Reb M biosynthesis is one key breach for breeding of S. rebaudiana, which has not been studied at the transcriptional level. Herein, five different S. rebaudiana varieties with distinct SGs contents, one cultivar having high stevioside content (HST), one cultivar having high Reb A content (HRA) and three cultivars having high Reb D and Reb M content (HDM1, HDM2, HDM3), were selected for RNA-seq analysis. In total, 131,655 de novo assembled unigenes were found in the RNA-seq data. According to Reb D and Reb M content divergence of S. rebaudiana accessions, 2186 differentially expressed genes (DEGs) were selected as potential genes related to Reb D and Reb M biosynthesis. Weighted Gene Co-expression Network Analysis (WGCNA) was used to explore the genes associated with the Reb D and Reb M biosynthesis. The unigenes from the positively associated turquoise module formed a layered co-expression network. There are 7 UDP-dependent glycosyltransferases (UGT) and 76 transcription factors (TFs) distributing at different regions which represented varying coherence of Reb D and Reb M biosynthesis. Particularly, two TFs having a strong correlation with two UGTs in the network were also discovered. The present study provided a comprehensive insight into networks for regulation of Reb D and Reb M contents in S. rebaudiana.

Transcriptomic Characterization of Nitrate-Enhanced Stevioside Glycoside Synthesis in Stevia (Stevia rebaudiana) Bertoni.

Nitrogen forms (nitrate (NO3-) or ammonium (NH4+)) are vital to plant growth and metabolism. In stevia (Stevia rebaudiana), it is important to assess whether nitrogen forms can influence the synthesis of the high-value terpene metabolites-steviol glycosides (SGs), together with the underlying mechanisms. Field and pot experiments were performed where stevia plants were fertilized with either NO3- or NH4+ nutrition to the same level of nitrogen. Physiological measurements suggested that nitrogen forms had no significant impact on biomass and the total nitrogen content of stevia leaves, but NO3--enhanced leaf SGs contents. Transcriptomic analysis identified 397 genes that were differentially expressed (DEGs) between NO3- and NH4+ treatments. Assessment of the DEGs highlighted the responses in secondary metabolism, particularly in terpenoid metabolism, to nitrogen forms. Further examinations of the expression patterns of SGs synthesis-related genes and potential transcription factors suggested that GGPPS and CPS genes, as well as the WRKY and MYB transcription factors, could be driving N form-regulated SG synthesis. We concluded that NO3-, rather than NH4+, can promote leaf SG synthesis via the NO3--MYB/WRKY-GGPPS/CPS module. Our study suggests that insights into the molecular mechanism of how SG synthesis can be affected by nitrogen forms.

Metabolic engineering for the synthesis of steviol glycosides: current status and future prospects.

With the pursuit of natural non-calorie sweeteners, steviol glycosides (SGs) have become one of the most popular natural sweeteners in the market. The SGs in Stevia are a mixture of SGs synthesized from steviol (a terpenoid). SGs are diterpenoids. Different SGs depend on the number and position of sugar groups on the core steviol backbone. This diversity comes from the processing of glycoside steviol by various glycosyltransferases. Due to the differences in glycosylation, each SG has unique sensory properties. At present, it is more complicated to extract high-quality SGs from plants, so the excavation of the metabolic pathways of engineered microorganisms to synthesize SGs has been extensively studied. Specifically, the expression of different glycosyltransferases in microbes is key to the synthesis of various SGs by engineered microorganisms. To trigger more researches on the functional characterization of the enzymes encoded by these genes, this review describes the latest research progresses of the related enzymes involved in SG biosynthesis and metabolic engineering.Key points• Outlines the research progress of key enzymes in the biosynthetic pathway of SGs• Factors affecting the catalytic capacity of stevia glucosyltransferase• Provide guidance for the efficient synthesis of SGs in microbial cell factories.

Comparative transcriptome analysis provides insights into steviol glycoside synthesis in stevia (Stevia rebaudiana Bertoni) leaves under nitrogen deficiency.

KEY MESSAGE: Transcriptome analysis revealed the potential mechanism of nitrogen regulating steviol glycosides synthesis via shifting of leaf carbon metabolic flux or inducing certain transcription factors. Nitrogen (N) plays key regulatory roles in both stevia (Stevia rebaudiana) growth and the synthesis of its functional metabolite steviol glycosides (SGs), but the mechanism by which this nutrient regulates SGs synthesis remains to be elucidated. To address this question, a pot experiment was performed in a greenhouse where stevia plants fertilized with N (the control as CK plants) and compared with plants without the supply of N. Physiological and biochemical analyses were conducted to test the growth and metabolic responses of plants to N regimes. Our results showed that N deficiency significantly inhibited plant growth and leaf photosynthesis, while increased leaf SGs contents in stevia (49.97, 46.64 and 84.80% respectively for rebaudioside A, stevioside, and rebaudioside C), which may be partly due to "concentration effect". Then, transcriptome analysis was conducted to understand the underlying mechanisms. A total of 535 differentially expressed genes were identified, and carbon metabolism-related events were highlighted by Gene Ontology and Kyoto Encyclopedia of Genes and Genomes. Many of these genes were significantly upregulated by N-deficiency, including those involved in "phenylpropanoid biosynthesis", "flavonoid biosynthesis" and "starch and sucrose metabolism". Our study also analyzed the expression patterns of SGs synthesis-related genes under two N regimes and the potential transcription factors linking N nutrition and SG metabolism. N-deficiency may promote SGs synthesis by changing the carbon metabolism flux or inducing certain transcription factors. Our results provide deeper insight into the relationship between N nutrition and SGs synthesis in stevia plants.