RNA-Seq for gene identification and transcript profiling of three Stevia rebaudiana genotypes

Rational Design for the Complete Synthesis of Stevioside in Saccharomyces cerevisiae

Stevioside is a secondary metabolite of diterpenoid glycoside production in plants. It has been used as a natural sweetener in various foods because of its high sweetness and low-calorie content. In this study, we constructed a Saccharomyces cerevisiae strain for the complete synthesis of stevioside using a metabolic engineering strategy. Firstly, the synthesis pathway of steviol was modularly constructed in S. cerevisiae BY4742, and the precursor pathway was strengthened. The yield of steviol was used as an indicator to investigate the expression effect of different sources of diterpene synthases under different combinations, and the strains with further improved steviol yield were screened. Secondly, glycosyltransferases were heterologously expressed in this strain to produce stevioside, the sequence of glycosyltransferase expression was optimized, and the uridine diphosphate-glucose (UDP-Glc) supply was enhanced. Finally, the results showed that the strain SST-302III-ST2 produced 164.89 mg/L of stevioside in a shake flask experiment, and the yield of stevioside reached 1104.49 mg/L in an experiment employing a 10 L bioreactor with batch feeding, which was the highest yield reported. We constructed strains with a high production of stevioside, thus laying the foundation for the production of other classes of steviol glycosides and holding good prospects for application and promotion.

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.

Novel insights into maize (Zea mays) development and organogenesis for agricultural optimization

In maize, intrinsic hormone activities and sap fluxes facilitate organogenesis patterning and plant holistic development; these hormone movements should be a primary focus of developmental biology and agricultural optimization strategies. Maize (Zea mays) is an important crop plant with distinctive life history characteristics and structural features. Genetic studies have extended our knowledge of maize developmental processes, genetics, and molecular ecophysiology. In this review, the classical life cycle and life history strategies of maize are analyzed to identify spatiotemporal organogenesis properties and develop a definitive understanding of maize development. The actions of genes and hormones involved in maize organogenesis and sex determination, along with potential molecular mechanisms, are investigated, with findings suggesting central roles of auxin and cytokinins in regulating maize holistic development. Furthermore, investigation of morphological and structural characteristics of maize, particularly node ubiquity and the alternate attachment pattern of lateral organs, yields a novel regulatory model suggesting that maize organ initiation and subsequent development are derived from the stimulation and interaction of auxin and cytokinin fluxes. Propositions that hormone activities and sap flow pathways control organogenesis are thoroughly explored, and initiation and development processes of distinctive maize organs are discussed. Analysis of physiological factors driving hormone and sap movement implicates cues of whole-plant activity for hormone and sap fluxes to stimulate maize inflorescence initiation and organ identity determination. The physiological origins and biogenetic mechanisms underlying maize floral sex determination occurring at the tassel and ear spikelet are thoroughly investigated. The comprehensive outline of maize development and morphogenetic physiology developed in this review will enable farmers to optimize field management and will provide a reference for de novo crop domestication and germplasm improvement using genome editing biotechnologies, promoting agricultural optimization.

Physiological Analysis and Transcriptome Sequencing Reveal the Effects of Salt Stress on Banana (Musa acuminata cv. BD) Leaf

The salinization of soil is a widespread environmental problem. Banana (Musa acuminata L.) is a salt-sensitive plant whose growth, development, and production are constrained by salt stresses. However, the tolerance mechanism of this salt-sensitive banana to salt stress is still unclear. This study aimed to investigate the influence of NaCl treatment on phenotypic, physiological, and transcriptome changes in bananas. We found that the content of root activity, MDA, Pro, soluble sugar, soluble protein, and antioxidant enzymes activity in salt-stress treatment were significantly higher than the control in bananas. Transcriptome sequencing result identified an overall of 3,378 differentially expressed genes (DEGs) in banana leaves, and the Kyoto Encyclopedia of Genes and Genomes analysis indicated that these DEGs were involved in phenylpropanoid biosynthesis process, ribosome process, starch and sucrose metabolism, amino sugar process, and plant hormone signal transduction process that had simultaneously changed their expression under salt stress, which indicated these DEGs may play a role in promoting BD banana growth under salt treatments. The genes which were enriched in the phenylpropanoid biosynthesis process, starch and sucrose metabolism process, amino sugar process, and plant hormone signal transduction process were specifically regulated to respond to the salt stress treatments. Here, totally 48 differentially expressed transcription factors (TFs), including WRKY, MYB, NAC, and bHLH, were annotated in BD banana under salt stress. In the phenylpropane biosynthesis pathway, all transcripts encoding key enzymes were found to be significantly up-regulated, indicating that the genes in these pathways may play a significant function in the response of BD banana to salt stress. In conclusion, this study provides new insights into the mechanism of banana tolerance to salt stress, which provides a potential application for the genetic improvement of banana with salt tolerance.

Glycosyltransferases: Mining, engineering and applications in biosynthesis of glycosylated plant natural products

UDP-Glycosyltransferases (UGTs) catalyze the transfer of nucleotide-activated sugars to specific acceptors, among which the GT1 family enzymes are well-known for their function in biosynthesis of natural product glycosides. Elucidating GT function represents necessary step in metabolic engineering of aglycone glycosylation to produce drug leads, cosmetics, nutrients and sweeteners. In this review, we systematically summarize the phylogenetic distribution and catalytic diversity of plant GTs. We also discuss recent progress in the identification of novel GT candidates for synthesis of plant natural products (PNPs) using multi-omics technology and deep learning predicted models. We also highlight recent advances in rational design and directed evolution engineering strategies for new or improved GT functions. Finally, we cover recent breakthroughs in the application of GTs for microbial biosynthesis of some representative glycosylated PNPs, including flavonoid glycosides (fisetin 3-O-glycosides, astragalin, scutellarein 7-O-glucoside), terpenoid glycosides (rebaudioside A, ginsenosides) and polyketide glycosides (salidroside, polydatin).

Melatonin-induced physiology and transcriptome changes in banana seedlings under salt stress conditions

Soil salinization poses a serious threat to the ecological environment and agricultural production and is one of the most common abiotic stresses in global agricultural production. As a salt-sensitive plant, the growth, development, and production of bananas (Musa acuminata L.) are restricted by salt stress. Melatonin is known to improve the resistance of plants to stress. The study analyzed the effects of 100 μM melatonin on physiological and transcriptome changes in banana varieties (AAA group cv. Cavendish) under 60 mmol/l of NaCl salt stress situation. The phenotypic results showed that the application of exogenous melatonin could maintain banana plants' health growth and alleviate the damage caused by salt stress. The physiological data show that the application of exogenous melatonin can enhance salt tolerance of banana seedlings by increasing the content of proline content and soluble protein, slowing down the degradation of chlorophyll, reducing membrane permeability and recovery of relative water content, increasing the accumulation of MDA, and enhancing antioxidant defense activity. Transcriptome sequencing showed that melatonin-induced salt tolerance of banana seedlings involved biological processes, molecular functions, and cellular components. We also found that differentially expressed genes (DEGs) are involved in a variety of metabolic pathways, including amino sugar and nucleotide sugar metabolism, phenylalanine metabolism, cyanoamino acid metabolism, starch and sucrose metabolism, and linoleic acid metabolism. These major metabolism and biosynthesis may be involved in the potential mechanism of melatonin under salt stress. Furthermore, some members of the transcription factor family, such as MYB, NAC, bHLH, and WRKY, might contribute to melatonin alleviating salt stress tolerance of the banana plant. The result laid a basis for further clarifying the salt stress resistance mechanism of bananas mediated by exogenous melatonin and provides theoretical bases to utilize melatonin to improve banana salt tolerance in the future.

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 (HDM₁, HDM₂, HDM₃), 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.

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

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.

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.

Synthesis and production of steviol glycosides: recent research trends and perspectives

Abstract

Steviol glycosides (SvGls) are plant secondary metabolites belonging to a class of chemical compounds known as diterpenes. SvGls have been discovered only in a few plant species, including in the leaves of Stevia rebaudiana Bertoni. Over the last few decades, SvGls have been extensively researched for their extraordinary sweetness. As a result, the nutritional and pharmacological benefits of these secondary metabolites have grown increasingly apparent. In the near future, SvGls may become a basic, low-calorie, and potent sweetener in the growing natural foods market, and a natural anti-diabetic remedy, a highly competitive alternative to commercially available synthetic drugs. Commercial cultivation of stevia plants and the technologies of SvGls extraction and purification from plant material have already been introduced in many countries. However, new conventional and biotechnological solutions are still being sought to increase the level of SvGls in plants. Since many aspects related to the biochemistry and metabolism of SvGls in vivo, as well as their relationship to the overall physiology of S. rebaudiana are not yet understood, there is also a great need for in-depth scientific research on this topic. Such research may have positive impact on optimization of the profile and SvGls concentration in plants and thus lead to obtaining desired yield. This research summarizes the latest approaches and developments in SvGls production.

Key points

• Steviol glycosides (SvGls) are found in nature in S. rebaudiana plants.

• They exhibit nutraceutical properties.

• This review provides an insight on different approaches to produce SvGls.

• The areas of research that still need to be explored have been identified.

Characterization and analysis of the transcriptome in Opisina arenosella from different developmental stages using single-molecule real-time transcript sequencing and RNA-seq

Opisina arenosella is one of the main pests harming coconut trees. To date, there have been few studies on the molecular genetics, biochemistry and physiology of O. arenosella at the transcriptional level, and there are no available reference genomes. Here, Illumina RNA sequencing combined with PacBio single-molecule real-time analysis was applied to study the transcriptome of this pest at different developmental stages, providing reference data for transcript expression analysis. Twelve samples of O. arenosella from different stages of development were sequenced using Illumina RNA sequencing, and the pooled RNA samples were sequenced with PacBio technology (Iso-Seq). A full-length transcriptome with 41,938 transcripts was captured, and the N50 and N90 lengths were 3543 bp and 1646 bp, respectively. A total of 36,925 transcripts were annotated in public databases, 6493 of which were long noncoding RNAs, while 2510 represented alternative splicing events. There were significant differences in the gene expression profiles at different developmental stages, with high levels of differential gene expression associated with growth, development, carbohydrate metabolism and immunity. This work provides resources and information for the study of the transcriptome and gene function of O. arenosella and provides a valuable foundation for understanding the changes in gene expression during development.

Characterizing glycosyltransferases by a combination of sequencing platforms applied to the leaf tissues of Stevia rebaudiana

Background

Stevia rebaudiana (Bertoni) is considered one of the most valuable plants because of the steviol glycosides (SGs) that can be extracted from its leaves. Glycosyltransferases (GTs), which can transfer sugar moieties from activated sugar donors onto saccharide and nonsaccharide acceptors, are widely distributed in the genome of S. rebaudiana and play important roles in the synthesis of steviol glycosides.

Results

Six stevia genotypes with significantly different concentrations of SGs were obtained by induction through various mutagenic methods, and the contents of seven glycosides (stevioboside, Reb B, ST, Reb A, Reb F, Reb D and Reb M) in their leaves were considerably different. Then, NGS and single-molecule real-time (SMRT) sequencing were combined to analyse leaf tissue from these six different genotypes to generate a full-length transcriptome of S. rebaudiana. Two phylogenetic trees of glycosyltransferases (SrUGTs) were constructed by the neighbour-joining method and successfully predicted the functions of SrUGTs involved in SG biosynthesis. With further insight into glycosyltransferases (SrUGTs) involved in SG biosynthesis, the weighted gene co-expression network analysis (WGCNA) method was used to characterize the relationships between SrUGTs and SGs, and forty-four potential SrUGTs were finally obtained, including SrUGT85C2, SrUGT74G1, SrUGT76G1 and SrUGT91D2, which have already been reported to be involved in the glucosylation of steviol glycosides, illustrating the reliability of our results.

Conclusion

Combined with the results obtained by previous studies and those of this work, we systematically characterized glycosyltransferases in S. rebaudiana and forty-four candidate SrUGTs involved in the glycosylation of steviol glucosides were obtained. Moreover, the full-length transcriptome obtained in this study will provide valuable support for further research investigating S. rebaudiana.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07195-5.

Transcriptional and physiological analyses of reduced density in apple provide insight into the regulation involved in photosynthesis

Different densities have a great influence on the physiological process and growth of orchard plants. Exploring the molecular basis and revealing key candidate genes for different densities management of orchard has great significance for production capacity improvement. In this study, transcriptome sequencing of apple trees was carried out at three different sampling heights to determine gene expression patterns under high density(HD) and low density(LD) and the physiological indices were measured to determine the effect of density change on plants. As a result, physiological indexes showed that the content of Chlorophyll, ACC, RUBP and PEP in the LD was apparently higher than that in control group(high density, HD). While the content of PPO and AO in the LD was noticeably lower than that in the HD. There were 3808 differentially expressed genes (DEGs) were detected between HD and LD, of which 1935, 2390 and 1108 DEGs were found in the three comparisons(middle-upper, lower-outer and lower-inner), respectively. 274 common differentially expressed genes (co-DEGs) were contained in all three comparisons. Functional enrichment and KEGG pathway analysis found these genes were involved in Carbon fixation in photosynthetic organisms, Circadian rhythm, Photosynthesis - antenna proteins, Photosynthesis, chlorophyll metabolism, Porphyrin, sugar metabolism and so on. Among these genes, LHCB family participated in photosynthesis as parts of photosystem II. In addition, SPA1, rbcL, SNRK2, MYC2, BSK, SAUR and PP2C are involved in Circadian rhythm, the expression of genes related to glycometabolism and hormone signaling pathway is also changed. The results revealed that the decrease of plant density changed the photosynthetic efficiency of leaves and the expression of photosynthesis-related genes, which provide a theoretical basis for the actual production regulation of apples.

Terpenoids as Potential Geroprotectors

Terpenes and terpenoids are the largest groups of plant secondary metabolites. However, unlike polyphenols, they are rarely associated with geroprotective properties. Here we evaluated the conformity of the biological effects of terpenoids with the criteria of geroprotectors, including primary criteria (lifespan-extending effects in model organisms, improvement of aging biomarkers, low toxicity, minimal adverse effects, improvement of the quality of life) and secondary criteria (evolutionarily conserved mechanisms of action, reproducibility of the effects on different models, prevention of age-associated diseases, increasing of stress-resistance). The number of substances that demonstrate the greatest compliance with both primary and secondary criteria of geroprotectors were found among different classes of terpenoids. Thus, terpenoids are an underestimated source of potential geroprotectors that can effectively influence the mechanisms of aging and age-related diseases.

Development of EST-based SSR and SNP markers in Gastrodia elata (herbal medicine) by sequencing, de novo assembly and annotation of the transcriptome

Tianma (Gastrodia elata Blume) has unique biological characteristics and high medicinal value. The wild resource of G. elata is being overutilized and should be conserved as it is already included in the list of endangered species in China. The population size of cultivated G. elata is small because of domestication bottleneck. Therefore, it is of utmost importance to evolve high-quality varieties and conserve wild resources of G. elata. In this study, we sequenced tuber transcriptomes of three major cultivated sub-species of Gastrodia elata, namely G. elata BI. f. elata, G. elata Bl. f. glauca S. Chow, and G. elata Bl. f. Viridis, and obtained about 7.8G clean data. The assembled high-quality reads of three sub-species were clustered into 56,884 unigenes. Of these, 31,224 (54.89%), 25,733 (45.24%), 22,629 (39.78%), and 11,856 (20.84%) unigenes were annotated by Nr, Swiss-Port, Eukaryotic Ortholog Groups (KOG), and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases, respectively. Here, a total of 3766 EST-SSRs and 128,921 SNPs were identified from the unigenes. The results not only offer huge number of genes that were responsible for the growth, development, and metabolism of bioactive components, but also a large number of molecular markers were detected for future studies on the conservation genetics and molecular breeding of G. elata.

Efficient and High-Quality RNA Isolation from Metabolite-Rich Tissues of Stevia rebaudiana, an Important Commercial Crop

Stevia rebaudiana, a perennial herb native to northeastern Paraguay, has gained immense attention globally over the recent decades due to the natural sweetness of its leaves. Like in most plants, this particular species contains high amount of secondary metabolites, thus rendering the isolation of high quality and quantity RNA extract for molecular applications rather challenging. An effective, high-yield and high-quality RNA isolation protocol for this economically important plant species was devised here based on the cetyltrimethylammonium bromide (CTAB) extraction method, with an additional genomic DNA (gDNA) removal step. DNA and other contaminants that may affect downstream applications were effectively removed. Our results exhibited that RNA samples isolated from the leaves and stems of Stevia rebaudiana using this improvised method are high in integrity and quality with RNA integrity number (RIN) of more than 8 and low in contaminants.

Overexpression of SrDXS1 and SrKAH enhances steviol glycosides content in transgenic Stevia plants

BACKGROUND

Stevia rebaudiana produces sweet-tasting steviol glycosides (SGs) in its leaves which can be used as natural sweeteners. Metabolic engineering of Stevia would offer an alternative approach to conventional breeding for enhanced production of SGs. However, an effective protocol for Stevia transformation is lacking.

RESULTS

Here, we present an efficient and reproducible method for Agrobacterium-mediated transformation of Stevia. In our attempts to produce transgenic Stevia plants, we found that prolonged dark incubation is critical for increasing shoot regeneration. Etiolated shoots regenerated in the dark also facilitated subsequent visual selection of transformants by green fluorescent protein during Stevia transformation. Using this newly established transformation method, we overexpressed the Stevia 1-deoxy-d-xylulose-5-phosphate synthase 1 (SrDXS1) and kaurenoic acid hydroxylase (SrKAH), both of which are required for SGs biosynthesis. Compared to control plants, the total SGs content in SrDXS1- and SrKAH-overexpressing transgenic lines were enhanced by up to 42-54% and 67-88%, respectively, showing a positive correlation with the expression levels of SrDXS1 and SrKAH. Furthermore, their overexpression did not stunt the growth and development of the transgenic Stevia plants.

CONCLUSION

This study represents a successful case of genetic manipulation of SGs biosynthetic pathway in Stevia and also demonstrates the potential of metabolic engineering towards producing Stevia with improved SGs yield.

Physiological Analysis and Transcriptome Profiling of Inverted Cuttings of Populus yunnanensis Reveal That Cell Wall Metabolism Plays a Crucial Role in Responding to Inversion

Inverted cuttings of Populus yunnanensis remain alive by rooting from the original morphological apex and sprouting from the base, but the lateral branches exhibit less vigorous growth than those of the upright plant. In this study, we examined the changes in hormone contents, oxidase activities, and transcriptome profiles between upright and inverted cuttings of P. yunnanensis. The results showed that the indole-3-acetic acid (IAA) and gibberellic acid (GA₃) contents were significantly lower in inverted cuttings than in upright cuttings only in the late growth period (September and October), while the abscisic acid (ABA) level was always similar between the two direction types. The biosynthesis of these hormones was surprisingly unrelated to the inversion of P. yunnanensis during the vegetative growth stage (July and August). Increased levels of peroxidases (PODs) encoded by 13 differentially expressed genes (DEGs) served as lignification promoters that protected plants against oxidative stress. Kyoto encyclopedia of genes and genomes (KEGG) enrichment analysis showed that most DEGs (107) were related to carbohydrate metabolism. Furthermore, altered activities of uridine diphosphate (UDP)-sugar pyrophosphorylase (USP, 15 DEGs) for nucleotide sugars, pectin methylesterase (PME, 7 DEGs) for pectin, and POD (13 DEGs) for lignin were important factors in the response of the trees to inversion, and these enzymes are all involved cell wall metabolism.

Transcriptome analysis of Nicotiana benthamiana infected by Tobacco curly shoot virus

BACKGROUND

Tobacco curly shoot virus (TbCSV) is a monopartite begomovirus associated with betasatellite (Tobacco curly shoot betasatellite, TbCSB), which causes serious leaf curl disease on tomato and tobacco in China. It is interesting that TbCSV induced severe upward leaf curling in Nicotiana benthamiana, but in the presence of TbCSB, symptoms changed to be downward leaf curling. However, the mechanism of interactions between viral pathogenicity, host defense, viral-betasatellite interactions and virus-host interactions remains unclear.

METHODS

In this study, RNA-seq was used to analyze differentially expressed genes (DEGs) in N. benthamiana plants infected by TbCSV (Y35A) and TbCSV together with TbCSB (Y35AB) respectively.

RESULTS

Through mapping to N. benthamiana reference genome, 59,814 unigenes were identified. Transcriptome analysis revealed that a total of 4081 and 3196 DEGs were identified in Y35AB vs CK (control check) and Y35A vs CK, respectively. Both GO and KEGG analyses were conducted to classify the DEGs. Ten of the top 15 GO terms were enriched in both DEGs of Y35AB vs CK and Y35A vs CK, and these enriched GO terms mainly classified into three categories including biological process, cellular component and molecular function. KEGG pathway analysis indicated that 118 and 111 pathways were identified in Y35AB vs CK and Y35A vs CK, respectively, of which nine and six pathways were significantly enriched. Three major pathways in Y35AB vs CK involved in metabolic pathways, carbon metabolism and photosynthesis, while those in Y35A vs CK were related to Ribosome, Glyoxylate and dicarboxylate metabolism and DNA replication. We observed that 8 PR genes were significantly up-regulated and 44 LRR-RLK genes were significantly differentially expressed in Y35A treatment or in Y35AB treatment. In addition, 7 and 13 genes were identified to be significantly changed in biosynthesis and signal transduction pathway of brassinosteroid (BR) and jasmonic acid (JA) respectively.

CONCLUSIONS

These results presented here would be particularly useful to further elucidate the response of the host plant against virus infection.

Validation of reference genes for RT-qPCR studies in Stevia rebaudiana in response to elicitor agents

Stevia rebaudiana is an important source of natural steviol glycosides and is of increasing interest in various fields of study. Therefore, understanding the molecular processes regulating its metabolism is of great importance. In this study, the stability of seven reference genes (18S ribosomal RNA, Actin, Aquaporin, Calmodulin, Eukaryote elongation factor 1-α, Malate dehydrogenase, and Ubiquitin) under the effect of three stress-related elicitors (methyl jasmonate, salicylic acid, and spermidine) was evaluated in stevia plants. We used RefFinder software, which makes use of the four main currently available algorithms for reference gene selection: geNorm, NormFinder, BestKeeper, and the Comparative ∆Ct method. The results indicated that Ubiquitin and Actin can be used as reference genes under all tested experimental conditions. The genes, 18S ribosomal RNA, traditionally used as a reference gene, along with Calmodulin, showed the lowest stability. The expression of Deoxyxylulose-5-phosphate synthase and Kaurenoic acid hydroxylase genes was used to confirm the validated reference genes, showing that inadequacy of the reference gene may lead to erroneous results. This is the first study on the stability of reference genes in Stevia rebaudiana plants, and is of great relevance for further analysis of the gene expression of the steviol glycoside biosynthetic pathway.

Comparative analysis of whole flower transcriptomes in the Zingiberales

The advancement of next generation sequencing technologies (NGS) has revolutionized our ability to generate large quantities of data at a genomic scale. Despite great challenges, these new sequencing technologies have empowered scientists to explore various relevant biological questions on non-model organisms, even in the absence of a complete sequenced reference genome. Here, we analyzed whole flower transcriptome libraries from exemplar species across the monocot order Zingiberales, using a comparative approach in order to gain insight into the evolution of the molecular mechanisms underlying flower development in the group. We identified 4,153 coding genes shared by all floral transcriptomes analyzed, and 1,748 genes that are only retrieved in the Zingiberales. We also identified 666 genes that are unique to the ginger lineage, and 2,001 that are only found in the banana group, while in the outgroup species Dichorisandra thyrsiflora J.C. Mikan (Commelinaceae) we retrieved 2,686 unique genes. It is possible that some of these genes underlie lineage-specific molecular mechanisms of floral diversification. We further discuss the nature of these lineage-specific datasets, emphasizing conserved and unique molecular processes with special emphasis in the Zingiberales. We also briefly discuss the strengths and shortcomings of de novo assembly for the study of developmental processes across divergent taxa from a particular order. Although this comparison is based exclusively on coding genes, with particular emphasis in transcription factors, we believe that the careful study of other regulatory mechanisms, such as non-coding RNAs, might reveal new levels of complexity, which were not explored in this work.

Development of Photoaffinity Probe for the Discovery of Steviol Glycosides Biosynthesis Pathway in Stevia rebuadiana and Rapid Substrate Screening

Functional discovery and characterization of the target enzymes responsible for the biosynthesis pathway coded for the genes is ongoing, and the unknown functional diversity of this class of enzymes has been revealed by genome sequencing. Commonly, it is feasible in annotating of biosynthetic genes of prokaryotes due to the existence of gene clusters of secondary metabolites. However, in eukaryotes, the biosynthetic genes are not compactly clustered in the way of prokaryotes. Hence, it remains challenging to identify the biosynthetic pathways of newly discovered natural products in plants. Steviol glycosides are one class of natural sweeteners found in high abundance in the herb Stevia rebaudiana. Here, we applied the chemoproteomic strategy for the proteomic profiling of the biosynthetic enzymes of steviol glycosides in Stevia rebaudiana. We not only identified a steviol-catalyzing UDP-glycosyltransferase (UGT) UGT73E1 involved in steviol glycoside biosynthesis but also built up a probe-based platform for the screening of potential substrates of functional uncharacterized UGT rapidly. This approach would be a complementary tool in mining novel synthetic parts for assembling of synthetic biological systems for the biosynthesis of other complex natural products.

Transcriptome analysis highlights key differentially expressed genes involved in cellulose and lignin biosynthesis of sugarcane genotypes varying in fiber content

Sugarcane (Saccharum spp. hybrids) is a potential lignocellulosic feedstock for biofuel production due to its exceptional biomass accumulation ability, high convertible carbohydrate content and a favorable energy input/output ratio. Genetic modification of biofuel traits to improve biomass conversion requires an understanding of the regulation of carbohydrate and lignin biosynthesis. RNA-Seq was used to investigate the transcripts differentially expressed between the immature and mature tissues of the sugarcane genotypes varying in fiber content. Most of the differentially expressed transcripts were found to be down-regulated during stem maturation, highlighting their roles in active secondary cell-wall development in the younger tissues of both high and low fiber genotypes. Several cellulose synthase genes (including CesA2, CesA4, CesA7 and COBRA-like protein), lignin biosynthesis-related genes (ρ-coumarate 3-hydroxylase, ferulate 5-hydroxylase, cinnamyl alcohol dehydrogenase and gentiobiase) and transcription regulators for the secondary cell-wall synthesis (including LIM, MYB, PLATZ, IAA24, C2H2 and C2C2 DOF zinc finger gene families) were exclusively differentially expressed between immature and mature tissues of high fiber genotypes. These findings reveal target genes for subsequent research on the regulation of cellulose and lignin metabolism.

Molecular dissection of transcriptional reprogramming of steviol glycosides synthesis in leaf tissue during developmental phase transitions in Stevia rebaudiana Bert

Stevia is a natural source of commercially important steviol glycosides (SGs), which share biosynthesis route with gibberellic acids (GAs) through plastidal MEP and cytosolic MVA pathways. Ontogeny-dependent deviation in SGs biosynthesis is one of the key factor for global cultivation of Stevia, has not been studied at transcriptional level. To dissect underlying molecular mechanism, we followed a global transcriptome sequencing approach and generated more than 100 million reads. Annotation of 41,262 de novo assembled transcripts identified all the genes required for SGs and GAs biosynthesis. Differential gene expression and quantitative analysis of important pathway genes (DXS, HMGR, KA13H) and gene regulators (WRKY, MYB, NAC TFs) indicated developmental phase dependent utilization of metabolic flux between SGs and GAs synthesis. Further, identification of 124 CYPs and 45 UGTs enrich the genomic resources, and their PPI network analysis with SGs/GAs biosynthesis proteins identifies putative candidates involved in metabolic changes, as supported by their developmental phase-dependent expression. These putative targets can expedite molecular breeding and genetic engineering efforts to enhance SGs content, biomass and yield. Futuristically, the generated dataset will be a useful resource for development of functional molecular markers for diversity characterization, genome mapping and evolutionary studies in Stevia.

Analysis of the first Taraxacum kok-saghyz transcriptome reveals potential rubber yield related SNPs

Taraxacum kok-saghyz (TK) is a potential alternative crop for natural rubber (NR) production, due to its high molecular weight rubber, short breeding cycle, and diverse environmental adaptation. However, improvements in rubber yield and agronomically relevant traits are still required before it can become a commercially-viable crop. An RNA-Seq based transcriptome was developed from a pool of roots from genotypes with high and low rubber yield. A total of 55,532 transcripts with lengths over 200 bp were de novo assembled. As many as 472 transcripts were significantly homologous to 49 out of 50 known plant putative rubber biosynthesis related genes. 158 transcripts were significantly differentially expressed between high rubber and low rubber genotypes. 21,036 SNPs were different in high and low rubber TK genotypes. Among these, 50 SNPs were found within 39 transcripts highly homologous to 49 publically-searched rubber biosynthesis related genes. 117 SNPs were located within 36 of the differentially expressed gene sequences. This comprehensive TK transcriptomic reference, and large set of SNPs including putative exonic markers associated with rubber related gene homologues and differentially expressed genes, provides a solid foundation for further genetic dissection of rubber related traits, comparative genomics and marker-assisted selection for the breeding of TK.

Association of gene expression with biomass content and composition in sugarcane

About 64% of the total aboveground biomass in sugarcane production is from the culm, of which ~90% is present in fiber and sugars. Understanding the transcriptome in the sugarcane culm, and the transcripts that are associated with the accumulation of the sugar and fiber components would facilitate the modification of biomass composition for enhanced biofuel and biomaterial production. The Sugarcane Iso-Seq Transcriptome (SUGIT) database was used as a reference for RNA-Seq analysis of variation in gene expression between young and mature tissues, and between 10 genotypes with varying fiber content. Global expression analysis suggests that each genotype displayed a unique expression pattern, possibly due to different chromosome combinations and maturation amongst these genotypes. Apart from direct sugar- and fiber-related transcripts, the differentially expressed (DE) transcripts in this study belonged to various supporting pathways that are not obviously involved in the accumulation of these major biomass components. The analysis revealed 1,649 DE transcripts between the young and mature tissues, while 555 DE transcripts were found between the low and high fiber genotypes. Of these, 151 and 23 transcripts respectively, were directly involved in sugar and fiber accumulation. Most of the transcripts identified were up-regulated in the young tissues (2 to 22-fold, FDR adjusted p-value <0.05), which could be explained by the more active metabolism in the young tissues compared to the mature tissues in the sugarcane culm. The results of analysis of the contrasting genotypes suggests that due to the large number of genes contributing to these traits, some of the critical DE transcripts could display less than 2-fold differences in expression and might not be easily identified. However, this transcript profiling analysis identified full-length candidate transcripts and pathways that were likely to determine the differences in sugar and fiber accumulation between tissue types and contrasting genotypes.

Transcriptomic profiling of Melilotus albus near-isogenic lines contrasting for coumarin content

Coumarin and its derivatives are widely used as fragrances in industrial products and have medical value. The goal of the present study was to discover genes and pathways related to coumarin biosynthesis in Melilotus albus using transcriptome analysis. The genes of five M. albus near-isogenic lines (NILs) that had different coumarin content and β-glucosidase activity according to the investigation of pedigree were quantified and then analysed by RNA-Seq. Using transcriptome analysis, differentially expressed genes (DEGs) were identified in two pairwise comparisons that differed in coumarin content as well as in two pairwise comparisons that differed in β-glucosidase activity. Gene expression pattern analysis suggested similar transcriptional trends in the genotypes with the same coumarin levels. Furthermore, the Kyoto Encyclopedia of Genes and Genomes (KEGG) database of DEGs was used to identify functional pathways associated with coumarin biosynthesis. We identified 111 unigenes, with several DEGs among them possibly being related to coumarin synthesis pathways. Unigenes encoding a hexokinase, an abscisic acid receptor, a phenylalanine ammonia-lyase (PAL) and two peroxidases particularly showed correspondence with the coumarin content of different genotypes. These results will contribute to a better understanding of the coumarin biosynthesis in M. albus.

A survey of the complex transcriptome from the highly polyploid sugarcane genome using full-length isoform sequencing and de novo assembly from short read sequencing

BACKGROUND

Despite the economic importance of sugarcane in sugar and bioenergy production, there is not yet a reference genome available. Most of the sugarcane transcriptomic studies have been based on Saccharum officinarum gene indices (SoGI), expressed sequence tags (ESTs) and de novo assembled transcript contigs from short-reads; hence knowledge of the sugarcane transcriptome is limited in relation to transcript length and number of transcript isoforms.

RESULTS

The sugarcane transcriptome was sequenced using PacBio isoform sequencing (Iso-Seq) of a pooled RNA sample derived from leaf, internode and root tissues, of different developmental stages, from 22 varieties, to explore the potential for capturing full-length transcript isoforms. A total of 107,598 unique transcript isoforms were obtained, representing about 71% of the total number of predicted sugarcane genes. The majority of this dataset (92%) matched the plant protein database, while just over 2% was novel transcripts, and over 2% was putative long non-coding RNAs. About 56% and 23% of total sequences were annotated against the gene ontology and KEGG pathway databases, respectively. Comparison with de novo contigs from Illumina RNA-Sequencing (RNA-Seq) of the internode samples from the same experiment and public databases showed that the Iso-Seq method recovered more full-length transcript isoforms, had a higher N50 and average length of largest 1,000 proteins; whereas a greater representation of the gene content and RNA diversity was captured in RNA-Seq. Only 62% of PacBio transcript isoforms matched 67% of de novo contigs, while the non-matched proportions were attributed to the inclusion of leaf/root tissues and the normalization in PacBio, and the representation of more gene content and RNA classes in the de novo assembly, respectively. About 69% of PacBio transcript isoforms and 41% of de novo contigs aligned with the sorghum genome, indicating the high conservation of orthologs in the genic regions of the two genomes.

CONCLUSIONS

The transcriptome dataset should contribute to improved sugarcane gene models and sugarcane protein predictions; and will serve as a reference database for analysis of transcript expression in sugarcane.

Transcriptomic Comparison Reveals Candidate Genes for Triterpenoid Biosynthesis in Two Closely Related Ilex Species

Native to Southern China, Ilex pubescens and Ilex asprella are frequently used in traditional Chinese medicine. Both of them produce a large variety of ursane-type triterpenoid saponins, which have been demonstrated to have different pharmacological effects. However, little is known about their biosynthesis. In this study, transcriptomic analysis of I. pubescens and comparison with its closely related specie I. asprella were carried out to identify potential genes involved in triterpenoid saponin biosynthesis. Through RNA sequencing (RNA-seq) and de novo transcriptome assembly of I. pubescens, a total of 68,688 UniGene clusters are obtained, of which 32,184 (46.86%) are successfully annotated by comparison with the sequences in major public databases (NCBI, Swiss-Prot, and KEGG). It includes 128 UniGenes related to triterpenoid backbone biosynthesis, 11 OSCs (oxidosqualene cyclases), 233 CYPs (cytochrome P450), and 269 UGTs (UDP-glycosyltransferases). By homology-based blast and phylogenetic analysis with well-characterized genes involved in triterpenoid saponin biosynthesis, 5 OSCs, 14 CYPs, and 1 UGT are further proposed as the most promising candidate genes. Transcriptomic comparison between two Ilex species using blastp and OrthoMCL method reveals high sequence similarity. All OSCs and UGTs as well as most CYPs are classified as orthologous genes, while only 5 CYPs in I. pubescens and 3 CYPs in I. asprella are species-specific. One of OSC candidates, named as IpAS1, was successfully cloned and expressed in Saccharomyces cerevisiae INVSc1. Analysis of the yeast extract by gas chromatography (GC) and gas chromatography-mass spectrometry (GC-MS) shows IpAS1 is a mixed amyrin synthase, producing α-amyrin and β-amyrin at ratio of 5:1, which is similar to its ortholog IaAS1 from I. asprella. This study is the first exploration to profile the transcriptome of I. pubescens, the generated data and gene models will facilitate further molecular studies on the physiology and metabolism in this plant. By comparative transcriptomic analysis, a series of candidate genes involved in the biosynthetic pathway of triterpenoid saponins are identified, providing new insight into their biosynthesis at transcriptome level.

Transcriptome Analysis to Identify the Putative Biosynthesis and Transport Genes Associated with the Medicinal Components of Achyranthes bidentata Bl

Achyranthes bidentata is a popular perennial medicine herb used for 1000s of years in China to treat various diseases. Although this herb has multiple pharmaceutical purposes in China, no transcriptomic information has been reported for this species. In addition, the understanding of several key pathways and enzymes involved in the biosynthesis of oleanolic acid and ecdysterone, two pharmacologically active classes of metabolites and major chemical constituents of A. bidentata root extracts, is limited. The aim of the present study was to characterize the transcriptome profile of the roots and leaves of A. bidentata to uncover the biosynthetic and transport mechanisms of the active components. In this study, we identified 100,987 transcripts, with an average length of 1146.8 base pairs. A total of 31,634 (31.33%) unigenes were annotated, and 12,762 unigenes were mapped to 303 pathways according to the Kyoto Encyclopedia of Genes and Genomes pathway database. Moreover, we identified a total of 260 oleanolic acid and ecdysterone genes encoding biosynthetic enzymes. Furthermore, the key enzymes involved in the oleanolic acid and ecdysterone synthesis pathways were analyzed using quantitative real-time polymerase chain reaction, revealing that the roots expressed these enzymes to a greater extent than the leaves. In addition, we identified 85 ATP-binding cassette transporters, some of which might be involved in the translocation of secondary metabolites.

Stevia Glycosides: Chemical and Enzymatic Modifications of Their Carbohydrate Moieties to Improve the Sweet-Tasting Quality

Stevia glycosides, extracted from the leaves of the plant Stevia rebaudiana Bertoni, display an amazing high degree of sweetness. As processed plant products, they are considered as excellent bio-alternatives for sucrose and artificial sweeteners. Being noncaloric and having beneficial properties for human health, they are the subject of an increasing number of studies for applications in food and pharmacy. However, one of the main obstacles for the successful commercialization of Stevia sweeteners, especially in food, is their slight bitter aftertaste and astringency. These undesirable properties may be reduced or eliminated by modifying the carbohydrate moieties of the steviol glycosides. A promising procedure is to subject steviol glycosides to enzymatic glycosylation, thereby introducing additional monosaccharide residues into the molecules. Depending on the number and positions of the monosaccharide units, the taste quality and sweetness potency of the compounds will vary. Many studies have been performed already, and this review summarizes the structures of native steviol glycosides and the recent data of modifications of the carbohydrate moieties that have been published to provide an overview of the current progress.

Transcriptome Analysis of Differentially Expressed Genes Involved in Proanthocyanidin Accumulation in the Rhizomes of Fagopyrum dibotrys and an Irradiation-Induced Mutant

The rhizome of Fagopyrum dibotrys is a traditional Chinese medicine that has recently gained attention due to substantial findings regarding its bioactive proanthocyanidin (PA) compounds. However, the molecular mechanism underlying PA accumulation in F. dibotrys remains elusive. We previously obtained an irradiation-induced mutant (RM_R) of F. dibotrys that had a higher PA content compared to that of the wild-type (CK_R). The present study aimed to elucidate the molecular mechanism underlying PA accumulation in F. dibotrys by comparing the rhizome transcriptomes of the irradiation-induced mutant and wild-type using RNA-seq analysis. A total of 53,540 unigenes were obtained, of which 29,901 (55.84%) were annotated based on BLAST searches against public databases, and 501 unique sequences were differentially expressed between the two samples, which consisted of 204 up-regulated and 297 down-regulated unigenes. Further analysis showed that the expression patterns of some unigenes encoding enzymes involved in PAs biosynthesis in F. dibotrys rhizomes differed between RM_R and CK_R. In addition, we identified transcription factor families and several cytochrome P450s that may be involved in PA regulation in F. dibotrys. Finally, 12 unigenes that encode PA biosynthetic enzymes were confirmed by qRT-PCR analysis. This study sheds light on the molecular mechanism underlying radiation-mediated flavonoid accumulation and regulation in F. dibotrys rhizomes. These results will also provide a platform for further functional genomic research on this particular species.

De novo Assembly and Characterization of the Transcriptome of Broomcorn Millet (Panicum miliaceum L.) for Gene Discovery and Marker Development

Broomcorn millet (Panicum miliaceum L.) is one of the world's oldest cultivated cereals, which is well-adapted to extreme environments such as drought, heat, and salinity with an efficient C4 carbon fixation. Discovery and identification of genes involved in these processes will provide valuable information to improve the crop for meeting the challenge of global climate change. However, the lack of genetic resources and genomic information make gene discovery and molecular mechanism studies very difficult. Here, we sequenced and assembled the transcriptome of broomcorn millet using Illumina sequencing technology. After sequencing, a total of 45,406,730 and 51,160,820 clean paired-end reads were obtained for two genotypes Yumi No. 2 and Yumi No. 3. These reads were mixed and then assembled into 113,643 unigenes, with the length ranging from 351 to 15,691 bp, of which 62,543 contings could be assigned to 315 gene ontology (GO) categories. Cluster of orthologous groups and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses assigned could map 15,514 unigenes into 202 KEGG pathways and 51,020 unigenes to 25 COG categories, respectively. Furthermore, 35,216 simple sequence repeats (SSRs) were identified in 27,055 unigene sequences, of which trinucleotides were the most abundant repeat unit, accounting for 66.72% of SSRs. In addition, 292 differentially expressed genes were identified between the two genotypes, which were significantly enriched in 88 GO terms and 12 KEGG pathways. Finally, the expression patterns of four selected transcripts were validated through quantitative reverse transcription polymerase chain reaction analysis. Our study for the first time sequenced and assembled the transcriptome of broomcorn millet, which not only provided a rich sequence resource for gene discovery and marker development in this important crop, but will also facilitate the further investigation of the molecular mechanism of its favored agronomic traits and beyond.

Transcriptome Analysis of Pepper (Capsicum annuum) Revealed a Role of 24-Epibrassinolide in Response to Chilling

Brassinosteroids (BRs) have positive effects on many processes during plant growth, development, and various abiotic stress responses. However, little information is available regarding the global gene expression of BRs in response to chilling stress in pepper. In this study, we used RNA sequencing to determine the molecular roles of 24-epibrassinolide (EBR) during a chilling stress response. There were 39,829 transcripts, and, among them, 656 were differently-expressed genes (DEGs) following EBR treatment (Chill+EBR) compared with the control (Chill only), including 335 up-regulated and 321 down-regulated DEGs. We selected 20 genes out of the 656 DEGs for RT-qPCR analysis to confirm the RNA-Seq. Based on GO enrich and KEGG pathway analysis, we found that photosynthesis was significantly up-enriched in biological processes, accompanied by significant increases in the net photosynthetic rate (Pn), Fv/Fm, and chlorophyll content. Furthermore, the results indicate that EBR enhanced endogenous levels of salicylic acid (SA) and jasmonic acid (JA) while suppressing the ethylene (ETH) biosynthesis pathway, suggesting that BRs function via a synergistic cross-talk with SA, JA, and ETH signaling pathways in response to chilling stress. In addition, EBR induced cellulose synthase-like protein and UDP-glycosyltransferase, suggesting a contribution to the formation of cell wall and hormone metabolism. EBR also triggered the calcium signaling transduction in cytoplasm, and activated the expression of cellular redox homeostasis related genes, such as GSTX1, PER72, and CAT2. This work, therefor, identified the specific genes showed different expression patterns in EBR-treated pepper and associated with the processes of hormone metabolism, redox, signaling, transcription, and defense. Our study provides the first evidence of the potent roles of BRs, at the transcription level, to induce the tolerance to chilling stress in pepper as a function of the combination of the transcriptional activities, signaling transduction, and metabolic homeostasis.

Comparative Transcriptomics Unravel Biochemical Specialization of Leaf Tissues of Stevia for Diterpenoid Production

Stevia (Stevia rebaudiana) produces not only a group of diterpenoid glycosides known as steviol glycosides (SGs), but also other labdane-type diterpenoids that may be spatially separated from SGs. However, their biosynthetic routes and spatial distribution in leaf tissues have not yet been elucidated. Here, we integrate metabolome and transcriptome analyses of Stevia to explore the biosynthetic capacity of leaf tissues for diterpenoid metabolism. Tissue-specific chemical analyses confirmed that SGs were accumulated in leaf cells but not in trichomes. On the other hand, Stevia leaf trichomes stored other labdane-type diterpenoids such as oxomanoyl oxide and agatholic acid. RNA sequencing analyses from two different tissues of Stevia provided a comprehensive overview of dynamic metabolic activities in trichomes and leaf without trichomes. These metabolite-guided transcriptomics and phylogenetic and gene expression analyses clearly identified specific gene members encoding enzymes involved in the 2-C-methyl-d-erythritol 4-phosphate pathway and the biosynthesis of steviol or other labdane-type diterpenoids. Additionally, our RNA sequencing analysis uncovered copalyl diphosphate synthase (SrCPS) and kaurene synthase1 (SrKS1) homologs, SrCPS2 and KS-like (SrKSL), which were specifically expressed in trichomes. In vitro and in planta assays showed that unlike SrCPS and SrKS1, SrCPS2 synthesized labda-13-en-8-ol diphosphate and successively catalyzed the formation of manoyl oxide and epi-manoyl oxide in combination with SrKSL. Our findings suggest that Stevia may have evolved to use distinct metabolic pathways to avoid metabolic interferences in leaf tissues for efficient production of diverse secondary metabolites.

Effects of growth stage and fulvic acid on the diversity and dynamics of endophytic bacterial community in Stevia rebaudiana Bertoni leaves

The aim of this study was to learn the interactions among the endophytic bacteria, the plant growth, the foliar spray of fulvic acid, and the accumulation of steviol glycosides in the leaves of Stevia rebaudiana. Metagenomic DNA was extracted from the Stevia leaves at different growth stages with or without the fulvic acid treatment; and the diversity of endophytic bacteria in Stevia leaves was estimated by pyrosequencing of 16S rRNA genes. As results, Proteobacteria, Actinobacteria, Bacteroidetes, and Firmicutes were found to be the dominant phyla despite the growth stages and fulvic acid application. Stevia growth stages strongly regulated composition of endophytic community. The genera Agrobacterium (12.3%) and Erwinia (7.2%) dominated in seedling stage were apparently declined in the vegetable and initial flowering stages, while Sphingomonas and Methylobacterium increased in mature leaves at harvest time, which showed that the mature leaves of Stevia preferred to accumulate some certain endophytic bacteria. Sphingomonas and Methylobacterium constituted an important part of the core endophytic community and were positively correlated with the stevioside content and UGT74G1 gene expression, respectively; while Erwinia, Agrobacterium, and Bacillus were negatively correlated with the stevioside accumulation. Fulvic acid treatment accelerated the variation of endophytes along the growth stages and increased the steviol glycosides content. This is the first study to reveal the community composition of endophytic bacteria in the Stevia leaves, to evidence the strong effects of growth stage and fulvic acid application on the endophytes of Stevia, and to demonstrate the correlation between the endophytic bacteria and the steviol glycosides accumulation.

De novo assembly of the Japanese lawngrass (Zoysia japonica Steud.) root transcriptome and identification of candidate unigenes related to early responses under salt stress

Japanese lawngrass (Zoysia japonica Steud.) is an important warm-season turfgrass that is able to survive in a range of soils, from infertile sands to clays, and to grow well under saline conditions. However, little is known about the molecular mechanisms involved in its resistance to salt stress. Here, we used high-throughput RNA sequencing (RNA-seq) to investigate the changes in gene expression of Zoysia grass at high NaCl concentrations. We first constructed two sequencing libraries, including control and NaCl-treated samples, and sequenced them using the Illumina HiSeq™ 2000 platform. Approximately 157.20 million paired-end reads with a total length of 68.68 Mb were obtained. Subsequently, 100,800 unigenes with an N50 length of 1104 bp were assembled using Trinity, among which 70,127 unigenes were functionally annotated (E ≤ 10(-5)) in the non-redundant protein (NR) database. Furthermore, three public databases, the Kyoto Encyclopedia of Genes and Genomes (KEGG), Swiss-prot, and Clusters of Orthologous Groups (COGs), were used for gene function analysis and enrichment. The annotated genes included 46 Gene Ontology (GO) terms, 120 KEGG pathways, and 25 COGs. Compared with the control, 6035 genes were significantly different (false discovery rate ≤0.01, |log2Ratio|≥1) in the NaCl-treated samples. These genes were enriched in 10 KEGG pathways and 58 GO terms, and subjected to 25 COG categories. Using high-throughput next-generation sequencing, we built a database as a global transcript resource for Z. japonica Steud. roots. The results of this study will advance our understanding of the early salt response in Japanese lawngrass roots.

Benchmarking Transcriptome Quantification Methods for Duplicated Genes in Xenopus laevis

Xenopus is an important model organism for the study of genome duplication in vertebrates. With the full genome sequence of diploid Xenopus tropicalis available, and that of allotetraploid X. laevis close to being finished, we will be able to expand our understanding of how duplicated genes have evolved. One of the key features in the study of the functional consequence of gene duplication is how their expression patterns vary across different conditions, and RNA-seq seems to have enough resolution to discriminate the expression of highly similar duplicated genes. However, most of the current RNA-seq analysis methods were not designed to study samples with duplicate genes such as in X. laevis. Here, various computational methods to quantify gene expression in RNA-seq data were evaluated, using 2 independent X. laevis egg RNA-seq datasets and 2 reference databases for duplicated genes. The fact that RNA-seq can measure expression levels of similar duplicated genes was confirmed, but long paired-end reads are more informative than short single-end reads to discriminate duplicated genes. Also, it was found that bowtie, one of the most popular mappers in RNA-seq analysis, reports significantly smaller numbers of unique hits according to a mapping quality score compared to other mappers tested (BWA, GSNAP, STAR). Calculated from unique hits based on a mapping quality score, both expression levels and the expression ratio of duplicated genes can be estimated consistently among biological replicates, demonstrating that this method can successfully discriminate the expression of each copy of a duplicated gene pair. This comprehensive evaluation will be a useful guideline for studying gene expression of organisms with genome duplication using RNA-seq in the future.