Transcriptome and Metabolite analysis reveal candidate genes of the cardiac glycoside biosynthetic pathway from Calotropis procera

Four enzymes control natural variation in the steroid core of Erysimum cardenolides

Plants commonly produce families of structurally related metabolites with similar defensive functions. This apparent redundancy raises the question of underlying molecular mechanisms and adaptive benefits of such chemical variation. Cardenolides, a class defensive compounds found in the wallflower genus Erysimum (L., Brassicaceae) and scattered across other plant families, show substantial structural variation, with glycosylation and hydroxylation being common modifications of a steroid core, which itself may vary in terms of stereochemistry and saturation. Through a combination of chemical mutagenesis and analysis of gene coexpression networks, we identified four enzymes involved in cardenolide biosynthesis in Erysimum that work together to determine stereochemistry at carbon 5 of the steroid core: Ec3βHSD, a 3β-hydroxysteroid dehydrogenase, Ec3KSI, a ketosteroid isomerase, EcP5βR2, a progesterone 5β-reductase, and EcDET2, a steroid 5α-reductase. We biochemically characterized the activity of these enzymes in vitro and generated CRISPR/Cas9 knockout lines to confirm activity in vivo. Cardenolide biosynthesis was not eliminated in any of the knockouts. Instead, mutant plants accumulated cardenolides with altered saturation and stereochemistry of the steroid core. Furthermore, we found variation in carbon 5 configuration among the cardenolides of 44 species of Erysimum, where the occurrence of some 5β-cardenolides is associated with the expression and sequence of P5βR2. This may have allowed Erysimum species to fine-tune their defensive profiles to target specific herbivore populations over the course of evolution.

SIGNIFICANCE STATEMENT

Plants use an array of toxic compounds to defend themselves from attack against insects and other herbivores. One mechanism through which plants may evolve more toxic compounds is through modifications to the structure of compounds they already produce. In this study, we show how plants in the wallflower genus Erysimum use four enzymes to fine-tune the structure of toxic metabolites called cardenolides. Natural variation in the sequence and expression of a single enzyme called progesterone 5β-reductase 2 partly explains the variation in cardenolides observed across the Erysimum genus. These alterations to cardenolide structure over the course of evolution suggests that there may be context-dependent benefits to Erysimum to invest in one cardenolide variant over another.

The Madagascar palm genome provides new insights on the evolution of Apocynaceae specialized metabolism

Specialized metabolites possess diverse interesting biological activities and some cardenolides- and monoterpene indole alkaloids- (MIAs) derived pharmaceuticals are currently used to treat human diseases such as cancers or hypertension. While these two families of biocompounds are produced by specific subfamilies of Apocynaceae, one member of this medicinal plant family, the succulent tree Pachypodium lamerei Drake (also known as Madagascar palm), does not produce such specialized metabolites. To explore the evolutionary paths that have led to the emergence and loss of cardenolide and MIA biosynthesis in Apocynaceae, we sequenced and assembled the P. lamerei genome by combining Oxford Nanopore Technologies long-reads and Illumina short-reads. Phylogenomics revealed that, among the Apocynaceae whose genomes have been sequenced, the Madagascar palm is so far the species closest to the common ancestor between MIA producers/non-MIA producers. Transposable elements, constituting 72.48% of the genome, emerge as potential key players in shaping genomic architecture and influencing specialized metabolic pathways. The absence of crucial MIA biosynthetic genes such as strictosidine synthase in P. lamerei and non-Rauvolfioideae species hints at a transposon-mediated mechanism behind gene loss. Phylogenetic analysis not only showcases the evolutionary divergence of specialized metabolite biosynthesis within Apocynaceae but also underscores the role of transposable elements in this intricate process. Moreover, we shed light on the low conservation of enzymes involved in the final stages of MIA biosynthesis in the distinct MIA-producing plant families, inferring independent gains of these specialized enzymes along the evolution of these medicinal plant clades. Overall, this study marks a leap forward in understanding the genomic dynamics underpinning the evolution of specialized metabolites biosynthesis in the Apocynaceae family, with transposons emerging as potential architects of genomics restructuring and gene loss.

Anti-Leishmania major activity of Calotropis procera extract by increasing ROS production and upregulating TNF-α, IFN-γ and iNOS mRNA expression under in vitro conditions

BACKGROUND

Leishmaniasis, caused by protozoan parasites of the genus Leishmania, is a neglected tropical disease with 700,000 to 1,000,000 global new cases annually. Adverse effects associated with expense, long-term treatment and drug resistance have made conventional therapies unfavorable, encouraging the search for alternative drugs based on plant products. In this study, the effect of Calotropis procera (Asclepiadaceae) extract against viability of promastigotes and amastigotes of Leishmania major was evaluated in vitro.

METHODS

The extract from the leaves of C. procera seedlings was prepared using a methanol maceration method. The colorimetric cell viability 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assay was used to determine the growth-inhibitory effect of the extract on promastigotes. The level of reactive oxygen species (ROS) in promastigote cultures was determined after treatment with the extract using the 2',7'-dichlorofluorescein diacetate (DCFH-DA) method and compared with untreated cultures (control). After exposure to the extract the expression levels of tumor necrosis factor-α (TNF-α), interferon gamma (IFN-γ) and inducible nitric oxide synthase (iNOS) genes were determined and compared to control in peripheral blood mononuclear cells (PBMCs) infected with L. major.

RESULTS

Based on the MTT assay, the C. procera extract significantly reduced the proliferation of L. major promastigotes with IC50 values of 377.28 and 222.44 μg/mL for 24 and 72 h, respectively (p < 0.01). After treatment with 222.44 and 377.28 μg/mL of C. procera extract, ROS production in L. major promastigote cultures increased 1.2- to 1.65-fold and 2- to 4-fold compared to the control, respectively (p < 0.05). C. procera extract induced significant increases in gene expression of TNF-α (2.76-14.83 fold), IFN-γ (25.63-threefold) and iNOS (16.32-3.97 fold) in infected PBMCs compared to control (p < 0.01).

CONCLUSIONS

On the basis of its anti-leishmanial activity, C. procera can be considered as a promising new plant source for the potential treatment of leishmaniasis.

Promiscuous CYP87A enzyme activity initiates cardenolide biosynthesis in plants

Cardenolides are specialized, steroidal metabolites produced in a wide array of plant families1,2. Cardenolides play protective roles in plants, but these molecules, including digoxin from foxglove (Digitalis spp.), are better known for treatment of congenital heart failure, atrial arrhythmia, various cancers and other chronic diseases3-9. However, it is still unknown how plants synthesize 'high-value', complex cardenolide structures from, presumably, a sterol precursor. Here we identify two cytochrome P450, family 87, subfamily A (CYP87A) enzymes that act on both cholesterol and phytosterols (campesterol and β-sitosterol) to form pregnenolone, the first committed step in cardenolide biosynthesis in the two phylogenetically distant plants Digitalis purpurea and Calotropis procera. Arabidopsis plants overexpressing these CYP87A enzymes ectopically accumulated pregnenolone, whereas silencing of CYP87A in D. purpurea leaves by RNA interference resulted in substantial reduction of pregnenolone and cardenolides. Our work uncovers the key entry point to the cardenolide pathway, and expands the toolbox for sustainable production of high-value plant steroids via synthetic biology.

Calotropis procera accumulates Uzarigenin and Calotropagenin in response to environmental lighting and drought

Calotropis procera (C. procera) was evaluated as a pharmaceutically useful plant and for its therapeutic effects in the most significant studies. Uzarigenin and Calotropagenin are significant components of this plant that have pharmacological effects on certain systems, including the digestive, immunological, and focal, and peripheral sensory systems. In this study, pathway genes are extracted from high throughput data acc.no. SRR1554320. Seven critical enzymes are involved in studying the effects of sunlight on the formation of Uzaragenin and Calotropagenin in C. procera before and after irrigation. Molecular identification and NCBI submission of six enzyme genes were successful; HSD (acc.no. OQ091761) for 3β-hydroxystroid dehydrogenase, OR (acc.no. OQ091762) for 5beta-pregnan oxidoreductase, MO (acc.no. OQ091763) for Pregnan monooxygenase, HOX (acc.no. OQ091764) for Steroid hydroxylase, MAT (acc.no. OQ091765) for Melonyletransferase, UHOX (acc.no. OQ091766) for Uzarigenin hydroxylase. During dawn after irrigation, the Uzargenin pathway showed the highest activity, however midday after irrigation was the lowest. The most period that showed high activity for the Uzargenin pathway was dawn after irrigation, however, midday after irrigation was the lowest. This data is confirmed by chromatography analysis (UPLC) to calculate the accumulation of Uzarigenin and Calotropagenin in different periods.

Cardenolide Increase in Foxglove after 2,1,3-Benzothiadiazole Treatment Reveals a Potential Link between Cardenolide and Phytosterol Biosynthesis

Cardenolides are steroidal metabolites in Digitalis lanata with potent cardioactive effects on animals. In plants, cardenolides are likely involved in various stress responses. However, the molecular mechanism of cardenolide increase during stresses is mostly unknown. Additionally, cardenolides are proposed to arise from cholesterol, but indirect results show that phytosterols may also be substrates for cardenolide biosynthesis. Here, we show that cardenolides increased after methyl jasmonate (MJ), sorbitol, potassium chloride (KCl) and salicylic acid analog [2,1,3-benzothiadiazole (BTH)] treatments. However, the expression of three known genes for cardenolide biosynthesis did not correlate well with these increases. Specifically, the expression of progesterone-5β-reductases (P5βR and P5βR2) did not correlate with the cardenolide increase. The expression of 3β-hydroxysteroid dehydrogenase (3βHSD) correlated with changes in cardenolide levels only during the BTH treatment. Mining the D. lanata transcriptome identified genes involved in cholesterol and phytosterol biosynthesis: C24 sterol sidechain reductase 1 (SSR1), C4 sterol methyl oxidase 1, and 3 (SMO1 and SMO3). Surprisingly, the expression of all three genes correlated well with the cardenolide increase after the BTH treatment. Phylogenetic analysis showed that SSR1 is likely involved in both cholesterol and phytosterol biosynthesis. In addition, SMO1 is likely specific to phytosterol biosynthesis, and SMO3 is specific to cholesterol biosynthesis. These results suggest that stress-induced increase of cardenolides in foxglove may correlate with cholesterol and phytosterol biosynthesis. In summary, this work shows that cardenolides are important for stress responses in D. lanata and reveals a potential link between phytosterol and cardenolide biosynthesis.

Comparative transcriptome and tissue-specific expression analysis of genes reveal tissue-cultured plants as an alternative source for phenylethanoids and phenylpropanoids in Rhodiola imbricata (Edgew.)

Rhodiola imbricata (Crassulaceae) is a traditional trans-Himalayan endangered medicinal herb with immense therapeutic applications. Over the years, over-exploitation, un-managed harvesting, and lack of captive cultivation procedures persuaded threat to its wild habitat. Plant tissue culture and RNA-Seq-based molecular bioprospection of key regulatory genes aid the understanding of molecular dynamics involved in specialized metabolites (phenylethanoids and phenylpropanoids) biosynthesis and its sustainable production. Hence, comparative transcriptomic analysis was performed using leaf and root tissues from the wild and tissue-cultured plants, revealing tissue-specific production of salidroside and rosavin. The transcriptome profiling resulted in 345 million high-quality reads yielding 92,380 unique transcripts with an N₅₀ of 1260 bp. Tissue-specific gene expression analysis revealed that both phenylethanoids and phenylpropanoids biosynthesis are predominantly associated with the shikimate pathway. In addition to RNA-Seq data, the downstream biosynthesis pathways genes viz., phospho-2-dehydro-3-deoxyheptonate aldolase (DAHPS), 3-dehydroquinate synthase (DHQS), shikimate kinase (SK), chorismate mutase (CM), arogenate dehydrogenase (TYRAAT), aromatic-L-amino-acid decarboxylase (TDC), phenylalanine ammonia-lyase (PAL), 4-coumarate-CoA ligase (4-CL), cinnamoyl-CoA reductase (CCR), and cinnamyl alcohol dehydrogenase (CAD) showed higher expression pattern in wild plant tissues compared to tissue-cultured plants. The transcript fold expression determined by RT-qPCR results followed similar patterns as those observed in RNA-seq and targeted metabolite profiling data. Salidroside and rosavin content in wild plants exhibited 2.40 fold and 1.77 fold increase accumulation compared to the tissue-cultured plant. The present investigation explained the tissue and condition-specific significant differences between the expression of proposed biosynthetic pathway genes and salidroside and rosavin content. Additionally, NAC, bHLH, and ARF were the most abundant transcription factor families found in the transcriptomic analysis of R. imbricata. The generated transcriptome dataset provides a valuable gene(s)/transcription factors hub that can be used for the sustainable production of salidroside and rosavin in R. imbricata under tissue culture conditions.

Transcriptome Analysis and Identification of the Cholesterol Side Chain Cleavage Enzyme BbgCYP11A1 From Bufo bufo gargarizans

Bufo bufo gargarizans Cantor are precious medicinal animals in traditional Chinese medicine (TCM). Bufadienolides as the major pharmacological components are generated from the venomous glands of B. bufo gargarizans. Bufadienolides are one type of cardiac aglycone with a six-member lactone ring and have properties of antitumor, cardiotonic, tonsillitis, and anti-inflammatory. The biosynthesis of bufadienolides is complex and unclear. This study explored the transcriptome of three different tissues (skin glands, venom glands, and muscles) of B. bufo gargarizans by high-throughput sequencing. According to the gene tissue–specific expression profile, 389 candidate genes were predicted possibly participating in the bufadienolides biosynthesis pathway. Then, BbgCYP11A1 was identified as a cholesterol side chain cleaving the enzyme in engineering yeast producing cholesterol. Furthermore, the catalytic activity of BbgCYP11A1 was studied with various redox partners. Interestingly, a plant NADPH-cytochrome P450 reductase (CPR) from Anemarrhena asphodeloides showed notably higher production than BbgAdx-2A-BbgAdR from B. bufo gargarizans. These results will provide certainly molecular research to reveal the bufadienolides biosynthesis pathway in B. bufo gargarizans.

Transcriptomic analysis reveals potential candidate pathways and genes involved in toxin biosynthesis in true toads

Synthesized chemical defenses have broadly evolved across countless taxa and are important in shaping evolutionary and ecological interactions within ecosystems. However, the underlying genomic mechanisms by which these organisms synthesize and utilize their toxins are relatively unknown. Herein, we use comparative transcriptomics to uncover potential toxin synthesizing genes and pathways, as well as interspecific patterns of toxin synthesizing genes across ten species of North American true toads (Bufonidae). Upon assembly and annotation of the ten transcriptomes, we explored patterns of relative gene expression and possible protein-protein interactions across the species to determine what genes and/or pathways may be responsible for toxin synthesis. We also tested our transcriptome dataset for signatures of positive selection to reveal how selection may be acting upon potential toxin producing genes. We assembled high quality transcriptomes of the bufonid parotoid gland, a tissue not often investigated in other bufonid related RNAseq studies. We found several genes involved in metabolic and biosynthetic pathways (e.g. steroid biosynthesis, terpenoid backbone biosynthesis, isoquinoline biosynthesis, glucosinolate biosynthesis) that were functionally enriched and/or relatively expressed across the ten focal species that may be involved in the synthesis of alkaloid and steroid toxins, as well as other small metabolic compounds that cause distastefulness in bufonids. We hope that our study lays a foundation for future studies to explore the genomic underpinnings and specific pathways of toxin synthesis in toads, as well as at the macroevolutionary scale across numerous taxa that produce their own defensive toxins.

Metabolomics integrated with transcriptomics reveals the distribution of iridoid and crocin metabolic flux in Gardenia jasminoides Ellis

Gardenia jasminoides Ellis (G. jasminoides) fruits are used as a resource for obtaining natural colorants and in traditional Chinese herbal medicine. However, G. jasminoides presents a relatively long flowering period and different ripening periods, so there are significant differences in the accumulation of metabolites in fruits of different colors. In addition, the complete metabolic pathways of iridoidsand crocins, which are used as medicinal composition of G. jasminoides, are poorly understood at present. In this research, we comprehensively compared the transcriptome and metabolites profiles of the developmental stages and locations of iridoid and crocin biosynthesis. A large number of differentially expressed genes (DEGs) and differentially accumulated metabolites (DAMs) were detected in four groups of samples, and clear variation in the pattern of metabolite abundance and gene expression were observed among different fruit colors and parts. Geniposide and gardenoside mainly accumulated in the sarcocarp of green fruit (GFS) and the sarcocarp of red fruit (FS), respectively. Crocin mainly accumulated in the peel and sarcocarp of red fruits. In the iridoid pathway, we hypothesized that there was a transport mechanism from the sarcocarp to the peel of G. jasminoides because of the inconsistent expression of G8O, 10-HGO and IS associated with differences in fruit ripening. UGTs play an important role in the biosynthesis of the active components of G. jasminoides. Combined transcriptome and metabonomics analysis showed a negative correlation between the biosynthesis of geniposide and crocin. The redirection of the metabolic flux and the regulation of key enzymes may be the main reasons for the changes in the biosynthesis of iridoid and crocin in G. jasminoides fruit. Our study expended valuable information for functional genomic library and provided new insights for metabolic engineering of secondary metabolite in G. Jasminoides.

Anti-proliferative, antioxidant effects of methanol extract of Calotropis procera leaf on lung cancer cells (H1299) and its ameliorative effect on expression of CD146 on blood cells

Background

Calotropis procera leaf is one of the plants commonly utilized in phytomedicine in Nigeria. The present investigation explored the use of the extracts on cell viability and apoptosis respectively. In this study, the expression of the Cluster of differentiation 146 (CD146) in the blood of lung cancer patients on regulatory T cells (Tregs) was determined. The antioxidant and anti-proliferative effects of methanol extracts of Caloropis procera leaf on lung cancer cell H1299 were investigated.

Methods

From the flow cytometry, the expression of the CD146+ in the T cells were evaluated using the healthy patient, adenocarcinoma, squamous, and small cell lung cancer respectively. The apoptosis of granulocytes, monocytes, lymphocytes, CD4+, and Treg were determined by 7-amino-actinomycin D/Annexin V-Allophycocyanin (APC) staining during the resting stage and after 24 h respectively. Immunofluorescence was conducted. Cell viability assay, hydroxyl (OH), hydrogen peroxide (H2O2) and nitric oxide (NO) scavenging radicals were conducted. Reducing power and flavonoid content of Calotropis procera were investigated. The effect of the Calotropis procera at different concentrations at 24hrs was determined.

Results

From the flow cytometry, the expression of the CD146+ on the T cells includes 4.60 % in healthy patients, 10.10, 12.20, 9.80 % in adenocarcinoma, squamous and small cell lung cancer. The apoptosis of granulocytes, monocytes, lymphocytes, CD4+ and Treg were determined by 7-amino-actinomycin D/Annexin V-APC staining during the resting stage and after 24 h which indicate that apoptosis also occurred on Treg. Immunofluorescence shows the presence of CD146 in lung cancer patient’s tissues. The methanol extracts of Calotropis procera leaf have antioxidant and anti-proliferative effects. Methanol extract of Calotropis procera leaf reduced CD146 expression on blood cells at 24 h.

Conclusion

Increased CD146 expression in the Treg of lung cancer patients indicates that it may be a possible target for the treatment of lung cancer by utilizing potent immunotherapy or natural products such as methanol extract of Calotropis procera leaf which may ameliorate the expression of CD146. Calotropis procera has antioxidant, inhibitory capacity on H1299 lung cancer cells, and the ability to scavenge OH, H2O2, and NO radicals. Hence, this investigation strengthens the phyto-medicinal properties of Calotropis procera.

Calotropis procera (Aiton) Dryand (Apocynaceae) as an anti-cancer agent against canine mammary tumor and osteosarcoma cells

Our goal was to evaluate phytochemical characterization and the antitumor potential of Calotropis procera. The phytochemical constitution of the crude extract (CE) revealed the presence of flavonoids, glycosides and cardenolide. The MTT assay was used to evaluate the cytotoxicity of CE, methanolic (MF) and ethyl acetate fractions (EAF) of C. procera in canine osteosarcoma cells (OST), canine mammary tumor (CMT), and canine skin fibroblasts (non-tumor cell). Doxorubicin was also used as a positive control. Results showed that CE, MF and EAF promoted a decrease in the viability of OST and CMT cells and did not alter the fibroblasts viability. C. procera also decreased the number of cells, corroborating to the decrease in proliferation and the cell cycle arrest in the G0/G1 phase. It was also evaluated the cell morphology by light and fluorescence microscopy, being demonstrated a reduction in cytoplasmic and cell rounding characteristic of programmed cell death. Moreover, flow cytometry data demonstrated that CE treatment promoted increase of caspase-3 and p53, showing that the cell death was activated in OST cells. In addition, there was a decrease in CD31, VEGF, osteopontin and TGF-β after CE treatment, suggesting that CE exerts its antitumor effect by reducing angiogenesis and tumor progression in OST cells. Moreover, CMT cells showed a reduction in PCNA after treatment with MF and CE. Analyzing the data together, C. procera, especially CE, showed an antitumor potential in both OST and CMT cells, encouraging us to continue investigating its use in cancer therapy.

Bioactivities and in silico study of Pergularia tomentosa L. phytochemicals as potent antimicrobial agents targeting type IIA topoisomerase, TyrRS, and Sap1 virulence proteins

Pergularia tomentosa L. (P. tomentosa) has been largely used in Tunisian folk medicine as remedies against skin diseases, asthma, and bronchitis. The main objectives of this study were to identify phytochemical compounds that have antioxidant and antimicrobial properties from the stem, leaves, and fruit crude methanolic extracts of P. tomentosa, and to search for tyrosyl-tRNA synthetase (TyrRS), topoisomerase type IIA, and Candidapepsin-1 (SAP1) enzyme inhibitors through molecular docking study. Phytochemical quantification revealed that fruit and leaves extracts displayed the highest total flavonoids (582 mg QE/g Ex; 219 mg QE/g Ex) and tannins content (375 mg TAE/g Ex; 216 mg TAE/g Ex), also exhibiting significant scavenging activity to decrease free radicals for ABTS, DPPH, β-carotene, and FRAP assay with IC₅₀ values (> 1 mg/mL). Additionally, promising antimicrobial activities towards different organs have been observed against several bacteria and Candida strains. From the liquid chromatography-mass spectrometry (LC-MS) analysis, five polyphenolic compounds, namely digitoxigenin, digitonin glycoside and calactina in the leaves, kaempferol in the fruit, and calotropagenin in the stems, were identified. They were also analyzed for their drug likeliness, based on computational methods. Molecular docking study affirmed that the binding affinity of calactin and actodigin to the active site of TyrRS, topoisomerase type IIA, and SAP1 target virulence proteins was the highest among the examined dominant compounds. Therefore, this study indicated that P. tomentosa methanolic extracts displayed great potential to become a potent antimicrobial agent and might be a promising source for therapeutic and nutritional functions. These phytocompounds could be further promoted as a candidate for drug discovery and development.

3D-surface MALDI mass spectrometry imaging for visualising plant defensive cardiac glycosides in Asclepias curassavica

Mass spectrometry-based imaging (MSI) has emerged as a promising method for spatial metabolomics in plant science. Several ionisation techniques have shown great potential for the spatially resolved analysis of metabolites in plant tissue. However, limitations in technology and methodology limited the molecular information for irregular 3D surfaces with resolutions on the micrometre scale. Here, we used atmospheric-pressure 3D-surface matrix-assisted laser desorption/ionisation mass spectrometry imaging (3D-surface MALDI MSI) to investigate plant chemical defence at the topographic molecular level for the model system Asclepias curassavica. Upon mechanical damage (simulating herbivore attacks) of native A. curassavica leaves, the surface of the leaves varies up to 700 μm, and cardiac glycosides (cardenolides) and other defence metabolites were exclusively detected in damaged leaf tissue but not in different regions of the same leaf. Our results indicated an increased latex flow rate towards the point of damage leading to an accumulation of defence substances in the affected area. While the concentration of cardiac glycosides showed no differences between 10 and 300 min after wounding, cardiac glycosides decreased after 24 h. The employed autofocusing AP-SMALDI MSI system provides a significant technological advancement for the visualisation of individual molecule species on irregular 3D surfaces such as native plant leaves. Our study demonstrates the enormous potential of this method in the field of plant science including primary metabolism and molecular mechanisms of plant responses to abiotic and biotic stress and symbiotic relationships.

Systematic Phytochemical Screening of Different Organs of Calotropis procera and the Ovicidal Effect of Their Extracts to the Foodstuff Pest Cadra cautella

In developing countries, crop deterioration is mainly caused by inappropriate storage conditions that promote insect infestation. Synthetic pesticides are associated with serious adverse effects on humans and the environment. Thus, finding alternative “green” insecticides is a very pressing need. Calotropis procera (Aiton) Dryand (Apocynaceae) growing in Saudi Arabia was selected for this purpose. LC-MS/MS analysis was applied to investigate the metabolic composition of different C. procera extracts. Particularly, C. procera latex and leaves showed a high presence of cardenolides including calactin, uscharidin, 15β-hydroxy-calactin, 16β-hydroxy-calactin, and 12β-hydroxy-calactin. The ovicidal activity of the extracts from different plant organs (flowers, leaves, branches, roots), and of the latex, against Cadra cautella (Walker) (Lepidoptera, Pyralidae) was assessed. Extracts of C. procera roots displayed the most potent activity with 50% of C. cautella eggs not hatching at 10.000 ppm (1%).

Sequencing, de novo assembly and annotation of Digitalis ferruginea subsp. schischkinii transcriptome

There is an increasing demand for elucidating the biosynthetic pathway of medicinal plants, which are capable of producing several metabolites with great potentials for industrial drug production. Digitalis species are important medicinal plants for the production of cardenolide compounds. Advancement on culture techniques is strictly related to our understanding of the genomic background of species. There are a limited number of genomic studies on Digitalis species. The goal of this study is to contribute to the genomic data of Digitalis ferruginea subsp. schischkinii by presenting transcriptome annotation. Digitalis ferruginea subsp. schischkinii has a limited distribution in Turkey and Transcaucasia, and has a high level of lanatoside C, an important cardenolide. In the study, we sequenced the cDNA library prepared from RNA pools of D. ferruginea subsp. schischkinii tissues treated with various stress conditions. Comprehensive bioinformatics approaches were used for de novo assembly and functional annotation of D. ferruginea subsp. schischkinii transcriptome sequence data along with TF families predictions and phylogenetic analysis. In the study, 58,369 unigenes were predicted and unigenes were annotated by analyzing the sequence data in the non-redundant (NR) protein database, the non-redundant nucleotide (NT) database, Gene Orthology (GO), EuKaryotic Orthologous Groups (KOG), Kyoto Encyclopedia of Genes and Genomes (KEGG), SwissProt, and InterPro databases. This study is the first transcriptome data for D. ferruginea subsp. schischkinii.

Localization of constituents for determining the age and parts of ginseng through ultraperfomance liquid chromatography quadrupole/time of flight-mass spectrometry combined with desorption electrospray ionization mass spectrometry imaging

Ginseng has been used for prevention and treatment of disease for thousands of years in China and many other Asian countries. Phytochemical studies have indicated that ginsenosides, polysaccharides, alkaloids, and phenolic acids are the active constituents of ginseng. Main and branch roots of ginseng exhibit distinct bioactive behavior. Furthermore, the bioactive behavior of ginseng depends on its age. Traditional analysis is complex preparation and provides inadequate of chemical information of the original distribution of analytes. Therefore, in this study, ultraperformance liquid chromatography quadrupole/time of flight-mass spectrometry (UPLC-QTOF MS) and desorption electrospray ionization mass spectrometry imaging (DESI-MSI) combined with orthogonal partial least squares discriminant analysis were used to discriminate ginseng in different age and parts of ginseng, and profiled distribution of selected markers. The results indicated that UPLC-QTOF-MS and DESI-MSI could be used to determine the parts and age of ginseng. Fifteen variables including five of protopanaxatriol (PPT), four of protopanaxadiol (PPD), and six of other types were assumed as markers for different parts of ginseng. Moreover, four variables of PPT, four of PPD, and ten of other types were used to determine the age of ginseng samples. An analysis of localization of markers indicated that malonyl ginsenoside, including malonyl-ginsenoside Rb₁, Rb₂, Rc, and Rd was mainly distributed in the corks. Neutral ginsenoside Rg₁, yesanchinoisde D, and chikusetsusaponin Iva were mainly distributed in the cork and phloem. Non-ginsenoside castanoside H, 20(S)-protopanaxatriol, unknown 2, saponin III and cistanoside C were distributed in all tissues. Ethyloleate, unknown 1 and monolinolein were distributed in the cork.

Comparative transcriptome analysis of root, stem, and leaf tissues of Entada phaseoloides reveals potential genes involved in triterpenoid saponin biosynthesis

Background

Entada phaseoloides (L.) Merr. is an important traditional medicinal plant. The stem of Entada phaseoloides is popularly used as traditional medicine because of its significance in dispelling wind and dampness and remarkable anti-inflammatory activities. Triterpenoid saponins are the major bioactive compounds of Entada phaseoloides. However, genomic or transcriptomic technologies have not been used to study the triterpenoid saponin biosynthetic pathway in this plant.

Results

We performed comparative transcriptome analysis of the root, stem, and leaf tissues of Entada phaseoloides with three independent biological replicates and obtained a total of 53.26 Gb clean data and 116,910 unigenes, with an average N50 length of 1218 bp. Putative functions could be annotated to 42,191 unigenes (36.1%) based on BLASTx searches against the Non-redundant, Uniprot, KEGG, Pfam, GO, KEGG and COG databases. Most of the unigenes related to triterpenoid saponin backbone biosynthesis were specifically upregulated in the stem. A total of 26 cytochrome P450 and 17 uridine diphosphate glycosyltransferase candidate genes related to triterpenoid saponin biosynthesis were identified. The differential expressions of selected genes were further verified by qPT-PCR.

Conclusions

The dataset reported here will facilitate the research about the functional genomics of triterpenoid saponin biosynthesis and genetic engineering of Entada phaseoloides.

Transcriptomic response to cold of thermophilous medicinal plant Marsdenia tenacissima

Extracts from Marsdenia tenacissima, involving tenacissoside H, I and G, have been used as remedies of cancer, inflammation and asthma. Low temperature serves as one of the main factors constrain the planting expansion and quality of M. tenacissima, but its functional mechanism has been known scarcely for the lack of genomic information and transcriptional profile. Here we investigated the transcriptomic responses of M. tenacissima under cold stress to gain insight into the molecular mechanism of low temperature sensitivity. Total RNAs were collected from samples obtained at 4-time points (after 0, 3, 6 and 48 h cold treatments with 4 °C, respectively), then used for library construction and sequenced on the Illumina Hiseq™ 4000 platform. Passing quality assessments, 500794 transcripts, and 206137 unigenes were de novo assembly out in Trinity v2.4.0, holding contig N50 of 2566 bp and unigene mean length of 754 bp. 44.20% of assembled unigenes were annotated to the well-known public protein database on a basis of sequence similarity. Using statistical comparison of the fragments per kilo base of transcript per million reads mapped (FPKM) values between conditions, 6082 group-specific differentially expressed genes (DEGs) were identified and considered as cold-responsive genes, which contained copious transcription factors and active secondary metabolism. Among them, 43 unigenes were constantly up-regulated expression along with cold time, which mainly implicated in the biosynthesis of secondary metabolites, carbon metabolism, RNA and DNA metabolism. Conversely, 21 unigenes involved in photosynthesis, cell wall, protein degradation, and transporters were downregulated continually with cold timescale. Experimentally, MtEF1α was chosen as the best housekeeping gene. Functional enrichments found that damaging of cold stress on M. tenacissima may be ascribed to inability of photosynthesis, ribsome processing, flavonoid biosynthesis and terpenoids degradation. Correlation analysis between cold induced transcription factors and tenacissoside biosynthesis-related genes indicated that 3β-HSD significant positively correlated with bHLH51, and 4-MSO with NF-YB, GRAS3, Trihelix, FAR1, MYB60, MYBS1, bZIP43. Further promoter clone found MYB-binding site in the promoter of 4-MSO. In view of the reported cold tolerance of MYB60, it is recommended as a potential candidate suitable for future molecular design of exaptation cultivation with high bioactive constituents.

High-throughput metabolomic and transcriptomic analyses vet the potential route of cerpegin biosynthesis in two varieties of Ceropegia bulbosa Roxb

Exploration with high-throughput transcriptomics and metabolomics of two varieties of Ceropegia bulbosa identifies candidate genes, crucial metabolites and a potential cerpegin biosynthetic pathway. Ceropegia bulbosa is an important medicinal plant, used in the treatment of various ailments including diarrhea, dysentery, and syphilis. This is primarily attributed to the presence of pharmaceutically active secondary metabolites, especially cerpegin. As this plant belongs to an endemic threatened category, genomic resources are not available hampering exploration on the molecular basis of cerpegin accumulation till now. Therefore, we undertook high-throughput metabolomic and transcriptomic analyses using different tissues from two varieties namely, C. bulbosa var. bulbosa and C. bulbosa var. lushii. Metabolomic analysis revealed spatial and differential accumulation of various metabolites. We chemically synthesized and characterized the cerpegin and its derivatives by liquid chromatography tandem-mass spectrometry (LC-MS/MS). Importantly, these comparisons suggested the presence of cerpegin and 5-allyl cerpegin in all C. bulbosa tissues. Further, de novo transcriptome analysis indicated the presence of significant transcripts for secondary metabolic pathways through the Kyoto encyclopedia of genes and genomes database. Tissue-specific profiling of transcripts and metabolites showed a significant correlation, suggesting the intricate mechanism of cerpegin biosynthesis. The expression of potential candidate genes from the proposed cerpegin biosynthetic pathway was further validated by qRT-PCR and NanoString nCounter. Overall, our findings propose a potential route of cerpegin biosynthesis. Identified transcripts and metabolites have built a foundation as new molecular resources that could facilitate future research on biosynthesis, regulation, and engineering of cerpegin or other important metabolites in such non-model plants.

CD146 T cells in lung cancer: its function, detection, and clinical implications as a biomarker and therapeutic target

CD146 alternatively called melanoma cell adhesion molecule (MCAM), is a biomarker and therapeutic target of clinical significance. It is found on different cells including the endothelial cells and lymphocytes which participate in heterotypic and homotypic ligand-receptor. This review concentrated on the CD146 expression T cells (or lymphocytes) centering on Treg in lung cancer. Here, we have also considered the vigorous investigation of CD146 mainly acknowledged new roles, essential mechanisms and clinical implications of CD146 in cancer. CD146 has progressively become a significant molecule, particularly recognized as a novel biomarker, prognosis and therapy for cancer. Hence, targeting CD146 expression by utilization of methanol extracts of Calotropis procera leaf may be useful for the treatment of carcinogenesis.

A draft genome and transcriptome of common milkweed (Asclepias syriaca) as resources for evolutionary, ecological, and molecular studies in milkweeds and Apocynaceae

Milkweeds (Asclepias) are used in wide-ranging studies including floral development, pollination biology, plant-insect interactions and co-evolution, secondary metabolite chemistry, and rapid diversification. We present a transcriptome and draft nuclear genome assembly of the common milkweed, Asclepias syriaca. This reconstruction of the nuclear genome is augmented by linkage group information, adding to existing chloroplast and mitochondrial genomic resources for this member of the Apocynaceae subfamily Asclepiadoideae. The genome was sequenced to 80.4× depth and the draft assembly contains 54,266 scaffolds ≥1 kbp, with N50 = 3,415 bp, representing 37% (156.6 Mbp) of the estimated 420 Mbp genome. A total of 14,474 protein-coding genes were identified based on transcript evidence, closely related proteins, and ab initio models, and 95% of genes were annotated. A large proportion of gene space is represented in the assembly, with 96.7% of Asclepias transcripts, 88.4% of transcripts from the related genus Calotropis, and 90.6% of proteins from Coffea mapping to the assembly. Scaffolds covering 75 Mbp of the Asclepias assembly formed 11 linkage groups. Comparisons of these groups with pseudochromosomes in Coffea found that six chromosomes show consistent stability in gene content, while one may have a long history of fragmentation and rearrangement. The progesterone 5β-reductase gene family, a key component of cardenolide production, is likely reduced in Asclepias relative to other Apocynaceae. The genome and transcriptome of common milkweed provide a rich resource for future studies of the ecology and evolution of a charismatic plant family.

Reference genes selection for Calotropis procera under different salt stress conditions

Calotropis procera is a perennial Asian shrub with significant adaptation to adverse climate conditions and poor soils. Given its increased salt and drought stress tolerance, C. procera stands out as a powerful candidate to provide alternative genetic resources for biotechnological approaches. The qPCR (real-time quantitative polymerase chain reaction), widely recognized among the most accurate methods for quantifying gene expression, demands suitable reference genes (RGs) to avoid over- or underestimations of the relative expression and incorrect interpretation. This study aimed at evaluating the stability of ten RGs for normalization of gene expression of root and leaf of C. procera under different salt stress conditions and different collection times. The selected RGs were used on expression analysis of three target genes. Three independent experiments were carried out in greenhouse with young plants: i) Leaf₁₀₀ = leaf samples collected 30 min, 2 h, 8 h and 45 days after NaCl-stress (100 mM NaCl); ii) Root₅₀ and iii) Root₂₀₀ = root samples collected 30 min, 2 h, 8 h and 1day after NaCl-stress (50 and 200 mM NaCl, respectively). Stability rank among the three algorithms used showed high agreement for the four most stable RGs. The four most stable RGs showed high congruence among all combination of collection time, for each software studied, with minor disagreements. CYP23 was the best RG (rank of top four) for all experimental conditions (Leaf₁₀₀, Root_50, and Root₂₀₀). Using appropriated RGs, we validated the relative expression level of three differentially expressed target genes (NAC78, CNBL4, and ND1) in Leaf₁₀₀ and Root₂₀₀ samples. This study provides the first selection of stable reference genes for C. procera under salinity. Our results emphasize the need for caution when evaluating the stability RGs under different amplitude of variable factors.

UPLC and ESI-MS analysis of metabolites of Rauvolfia tetraphylla L. and their spatial localization using desorption electrospray ionization (DESI) mass spectrometric imaging

Rauvolfia tetraphylla L. (family Apocynaceae), often referred to as the wild snakeroot plant, is an important medicinal plant and produces a number of indole alkaloids in its seeds and roots. The plant is often used as a substitute for Ravuolfia serpentine (L.) Benth. ex Kurz known commonly as the Indian snakeroot plant or sarphagandha in the preparation of Ayurvedic formulations for a range of diseases including hypertension. In this study, we examine the spatial localization of the various indole alkaloids in developing fruits and plants of R. tetraphylla using desorption electrospray ionization mass spectrometry imaging (DESI-MSI). A semi-quantitative analysis of the various indole alkaloids was performed using UPLC-ESI/MS. DESI-MS images showed that the distribution of ajmalcine, yohimbine, demethyl serpentine and mitoridine are largely localized in the fruit coat while that for ajmaline is restricted to mesocarp of the fruit. At a whole plant level, the ESI-MS intensities of many of the ions were highest in the roots and lesser in the shoot region. Within the root tissue, except sarpagine and ajmalcine, all other indole alkaloids occurred in the epidermal and cortex tissues. In leaves, only serpentine, ajmalcine, reserpiline and yohimbine were present. Serpentine was restricted to the petiolar region of leaves. Principal component analysis based on the presence of the indole alkaloids, clearly separated the four tissues (stem, leaves, root and fruits) into distinct clusters. In summary, the DESI-MSI results indicated a clear tissue localization of the various indole alkaloids, in fruits, leaves and roots of R. tetraphylla. While it is not clear of how such localization is attained, we discuss the possible pathways of indole alkaloid biosynthesis and translocation during fruit and seedling development in R. tetraphylla. We also briefly discuss the functional significance of the spatial patterns in distribution of metabolites.

Identification, expression, and phylogenetic analyses of terpenoid biosynthesis-related genes in secondary xylem of loblolly pine (Pinus taeda L.) based on transcriptome analyses

Loblolly pine (Pinus taeda L.) is one of the most important species for oleoresin (a mixture of terpenoids) in South China. The high oleoresin content of loblolly pine is associated with resistance to bark beetles and other economic benefits. In this study, we conducted transcriptome analyses of loblolly pine secondary xylem to gain insight into the genes involved in terpenoid biosynthesis. A total of 372 unigenes were identified as being critical for oleoresin production, including genes for ATP-binding cassette (ABC) transporters, the cytochrome P450 (CYP) protein family, and terpenoid backbone biosynthesis enzymes. Six key genes involved in terpenoid biosynthetic pathways were selected for multiple sequence alignment, conserved motif prediction, and phylogenetic and expression profile analyses. The protein sequences of all six genes exhibited a higher degree of sequence conservation, and upstream genes were relatively more conserved than downstream genes in terpenoid biosynthetic pathways. The N-terminal regions of these sequences were less conserved than the C-terminal ends, as the N-terminals were quite diverse in both length and composition. The phylogenetic analyses revealed that most genes originated from gene duplication after species divergence, and partial genes exhibited incomplete lineage sorting. In addition, the expression profile analyses showed that all six genes exhibited high expression levels during the high-oleoresin-yielding phase.

Metabolic Profiling of Pleurotus tuoliensis During Mycelium Physiological Maturation and Exploration on a Potential Indicator of Mycelial Maturation

Pleurotus tuoliensis is a valuable and rare edible fungus with extremely high nutritional and medicinal value. However, the relative immaturity of P. tuoliensis cultivation technology leads to fluctuating yields and quality. The main difficulty in P. tuoliensis cultivation is estimate of mycelial maturity. There is currently no measurable indicator that clearly characterizes the physiological maturation of mycelia. The aim of this study was to identify potential indicators of physiological maturation for P. tuoliensis mycelia by using metabolomics approach. A metabolite profiling strategy involving gas chromatography-mass spectrometry (GC/MS) was used to analyze changes to extracellular metabolites in mycelia collected at mycelium physiological maturation period (MPMP) day 0, MPMP day 35 at 17°C and MPMP day 35 at 29°C. 72 differential metabolites (37.8% up-regulated and 62.2% down-regulated) were identified based on the selected criteria [variable important in projection (VIP) greater than 1.0 and p < 0.01]. Metabolic pathways enrichment analysis showed that these metabolites are involved in glycolysis, organic acid metabolism, amino acid metabolism, tricarboxylic acid cycle (TCA), sugar metabolism, nicotinate and nicotinamide metabolism, and oxidative phosphorylation. In addition, the pyrimidine synthesis pathway was significantly activated during mycelium physiological maturation of P. tuoliensis. The abundance of N-carbamoyl-L-aspartate (CA-asp), a component of this pathway, was significantly increased at MPMP day 35, which motivated us to explore its potential as an indicator for physiological maturation of mycelia. The content of CA-asp in mycelia changed in a consistent manner during physiological maturation. The feasibility of using CA-asp as an indicator for mycelial maturation requires further investigation.

Control of β-sitosterol biosynthesis under light and watering in desert plant Calotropis procera

Most scientific studies on Calotropis procera refer to the plant as an important source of pharmaceutical compounds and its valuable benefits in medicine. One of the most important substances in this plant is the potential immunostimulant β-sitosterol (BS) that acts in improving human health. This study focused on the effects of lighting before and after irrigation on the BS accumulation pathway namely steroid biosynthesis. Studying the enzymes in BS biosynthetic pathway indicated the upregulation at dawn and predusk of the SMT2 and SMO2 genes encoding sterol methyltransferase 2 and methylsterol monooxygenase, two key enzymes in BS accumulation in C. procera. The results almost indicated no regulation at the different time points of the CYP710A gene encoding sterol 22-desaturase, an enzyme that acts in depleting β-sitosterol towards the biosynthesis of stigmasterol. RNA-Seq data was validated via quantitative RT-PCR and results were positive. The data of ultra-performance liquid chromatography-tandem mass spectrometry analysis with regard to BS accumulation also aligned with those of RNA-Seq analysis. We focused on the effects of light before and after watering on BS accumulation in C. procera. Our results show that BS accumulation is high at dawn in both dehydrated and well-watered condition. While, the BS was dramatically decrease at midday in well-watered plants. This increase/decrease in BS content is correlated with rates of expression of SMT 2 gene. This gene is a key convertor between the different branches in the cardiac glycoside biosynthesis. Accordingly, it could be suggested that BS (or one of the descendent product) may play an important role in C. procera tolerance to drought/light intensity conditions.

Protection afforded by methanol extract of Calotropis procera latex in experimental model of colitis is mediated through inhibition of oxidative stress and pro-inflammatory signaling

Calotropis procera, a latex producing plant is known to possess medicinal properties including its beneficial effect in gastrointestinal disorders. The anti-inflammatory effect of its latex in various experimental models is noteworthy and in light of this the present study was carried out with an objective to evaluate its efficacy in ulcerative colitis, an inflammatory condition of the colon. Colitis was induced in rats by acetic acid and the rats were divided into four groups where one group served as experimental control and the other groups were treated with two doses of methanol extract of dried latex of C. procera (MeDL; 50 and 150 mg/kg) and mesalazine (MSZ; 300 mg/kg). The study also included normal control (NC) group for comparison of various parameters related to colon like macroscopic changes, ulcer score, adherent mucus content, weight/length ratio, small intestinal transit, oxidative stress and inflammatory markers, tissue histology and immunoreactivity of cyclooxygenase-2 (COX-2), inducible nitric oxide synthase (iNOS) and nuclear factor kappa beta (NFκB) subunit p65. Treatment of colitic rats with MeDL produced a significant reduction in colonic mucosal damage as revealed by macroscopic and microscopic evaluation and normalization of tissue levels of oxidative stress markers and pro-inflammatory mediators. The protection afforded by MeDL was also evident from its restorative effect on tissue histology and expression of COX-2, iNOS and NFκB(p65). This study shows that by targeting oxidative stress and NFκB(p65) mediated pro-inflammatory signaling, the latex of C. procera affords protection in colitis and its effect was comparable to that of mesalazine. This study suggests that latex of C. procera could serve as a promising therapeutic option for treating inflammatory conditions of the colon.

De novo transcriptome analysis deciphered polyoxypregnane glycoside biosynthesis pathway in Gymnema sylvestre

Gymnema sylvestre is an important medicinal plant containing antidiabetic activity. Through de novo transcriptomic study, the pathways of polyoxypregnane glycosides were explored and candidate genes of these pathways were identified in G. sylvestre. High-quality raw reads were assembled into transcripts which resulted in 193,615 unigenes. These unigenes further decoded 58,274 coding DNA sequences (CDSs). Functional annotation of predicted CDSs was carried out using the protein databases, i.e., NCBI's non-redundant, Uniprot and Pfam. Eukaryotic orthologous group (KOG) classification and transcription factor analysis has revealed most CDS-enriched categories as "Signal transduction mechanism" and "Basic Helix loop helix" (bHLH) transcription factor family, respectively. A total of 16,569 CDSs were assigned minimum one Gene Ontology (GO) term. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis disclosed 235 CDSs which represented total 27 genes of pregnane glycoside pathways and 19 CDSs represented 10 important enzymes of polyoxypregnane glycoside biosynthesis, i.e., sterol 24-C-methyltransferase, cycloeucalenol cycloisomerase, Δ¹⁴-sterol reductase, C-8,7 sterol isomerase, sterol methyltransferase 2, C-5 sterol desaturase, sterol Δ⁷ reductase, Δ²⁴ sterol reductase, 3β-hydroxysteroid dehydrogenase and progesterone 5β reductase (5βPOR). This transcriptome analysis provided an important resource for future functional genomic studies in G. sylvestre.

De novo transcriptomic analysis of leaf and fruit tissue of Cornus officinalis using Illumina platform

Cornus officinalis is one of the most widely used medicinal plants in China and other East Asian countries to cure diseases such as liver, kidney, cardiovascular diseases and frequent urination for thousands of years. It is a Level 3 protected species, and is one of the 42 national key protected wild species of animals and plants in China. However, the genetics and molecular biology of C. officinalis are poorly understood, which has hindered research on the molecular mechanism of its metabolism and utilization. Hence, enriching its genomic data and information is very important. In recent years, the fast-growing technology of next generation sequencing has provided an effective path to gain genomic information from nonmodel species. This study is the first to explore the leaf and fruit tissue transcriptome of C. officinalis using the Illumina HiSeq 4000 platform. A total of 57,954,134 and 60,971,652 clean reads from leaf and fruit were acquired, respectively (GenBank number SRP115440). The pooled reads from all two libraries were assembled into 56,392 unigenes with an average length 856 bp. Among these, 41,146 unigenes matched with sequences in the NCBI nonredundant protein database. The Gene Ontology database assigned 24,336 unigenes with biological process (83.26%), cellular components (53.58%), and molecular function (83.93%). In addition, 10,808 unigenes were assigned a KOG functional classification by the KOG database. Searching against the KEGG pathway database indicated that 18,435 unigenes were mapped to 371 KEGG pathways. Moreover, the edgeR database identified 4,585 significant differentially expressed genes (DEGs), of which 1,392 were up-regulated and 3,193 were down-regulated in fruit tissue compared with leaf tissue. Finally, we explored 581 transcription factors with 50 transcription factor gene families. Most DEGs and transcription factors were related to terpene biosynthesis and secondary metabolic regulation. This study not only represented the first de novo transcriptomic analysis of C. officinalis but also provided fundamental information on its genes and biosynthetic pathway. These findings will help us explore the molecular metabolism mechanism of terpene biosynthesis in C. officinalis.

Genome Assembly and Annotation of the Medicinal Plant Calotropis gigantea, a Producer of Anticancer and Antimalarial Cardenolides

Calotropis gigantea produces specialized secondary metabolites known as cardenolides, which have anticancer and antimalarial properties. Although transcriptomic studies have been conducted in other cardenolide-producing species, no nuclear genome assembly for an Asterid cardenolide-producing species has been reported to date. A high-quality de novo assembly was generated for C. gigantea, representing 157,284,427 bp with an N50 scaffold size of 805,959 bp, for which quality assessments indicated a near complete representation of the genic space. Transcriptome data in the form of RNA-sequencing libraries from a developmental tissue series was generated to aid the annotation and construction of a gene expression atlas. Using an ab initio and evidence-driven gene annotation pipeline, 18,197 high-confidence genes were annotated. Homologous and syntenic relationships between C. gigantea and other species within the Apocynaceae family confirmed previously identified evolutionary relationships, and suggest the emergence or loss of the specialized cardenolide metabolites after the divergence of the Apocynaceae subfamilies. The C. gigantea genome assembly, annotation, and RNA-sequencing data provide a novel resource to study the cardenolide biosynthesis pathway, especially for understanding the evolutionary origin of cardenolides and the engineering of cardenolide production in heterologous organisms for existing and novel pharmaceutical applications.

Transcriptomic and metabolic responses of Calotropis procera to salt and drought stress

BACKGROUND

Calotropis procera is a wild plant species in the family Apocynaceae that is able to grow in harsh, arid and heat stressed conditions. Understanding how this highly adapted plant persists in harsh environments should inform future efforts to improve the hardiness of crop and forage plant species. To study the plant response to droμght and osmotic stress, we treated plants with polyethylene glycol and NaCl and carried out transcriptomic and metabolomics measurements across a time-course of five days.

RESULTS

We identified a highly dynamic transcriptional response across the time-course including dramatic changes in inositol signaling, stress response genes and cytokinins. The resulting metabolome changes also involved sharp increases of myo-inositol, a key signaling molecule and elevated amino acid metabolites at later times.

CONCLUSIONS

The data generated here provide a first glimpse at the expressed genome of C. procera, a plant that is exceptionally well adapted to arid environments. We demonstrate, through transcriptome and metabolome analysis that myo-inositol signaling is strongly induced in response to drought and salt stress and that there is elevation of amino acid concentrations after prolonged osmotic stress. This work should lay the foundations of future studies in adaptation to arid environments.

Comparative transcriptome analysis of shoot and root tissue of Bacopa monnieri identifies potential genes related to triterpenoid saponin biosynthesis

Background

Bacopa monnieri commonly known as Brahmi is utilized in Ayurveda to improve memory and many other human health benefits. Bacosides enriched standardized extract of Bacopa monnieri is being marketed as a memory enhancing agent. In spite of its well known pharmacological properties it is not much studied in terms of transcripts involved in biosynthetic pathway and its regulation that controls the secondary metabolic pathway in this plant. The aim of this study was to identify the potential transcripts and provide a framework of identified transcripts involved in bacosides production through transcriptome assembly.

Results

We performed comparative transcriptome analysis of shoot and root tissue of Bacopa monnieri in two independent biological replicate and obtained 22.48 million and 22.0 million high quality processed reads in shoot and root respectively. After de novo assembly and quantitative assessment total 26,412 genes got annotated in root and 18,500 genes annotated in shoot sample. Quality of raw reads was determined by using SeqQC-V2.2. Assembled sequences were annotated using BLASTX against public database such as NR or UniProt. Searching against the KEGG pathway database indicated that 37,918 unigenes from root and 35,130 unigenes from shoot were mapped to 133 KEGG pathways. Based on the DGE data we found that most of the transcript related to CYP450s and UDP-glucosyltransferases were specifically upregulated in shoot tissue as compared to root tissue. Finally, we have selected 43 transcripts related to secondary metabolism including transcription factor families which are differentially expressed in shoot and root tissues were validated by qRT-PCR and their expression level were monitored after MeJA treatment and wounding for 1, 3 and 5 h.

Conclusions

This study not only represents the first de novo transcriptome analysis of Bacopa monnieri but also provides information about the identification, expression and differential tissues specific distribution of transcripts related to triterpenoid sapogenin which is one of the most important pharmacologically active secondary metabolite present in Bacopa monnieri. The identified transcripts in this study will establish a foundation for future studies related to carrying out the metabolic engineering for increasing the bacosides biosynthesis and its regulation for human health benefits.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3865-5) contains supplementary material, which is available to authorized users.