High

Scientific Reports volume 12, Article number: 9347 (2022) Cite this article

2130 Accesses

2 Altmetric

Metrics details

Eleusine indica, which is used in traditional medicine, exhibits antiproliferative activity against several cancer cell lines. However, metabolomic studies to evaluate the metabolite changes induced by E. indica in cancer cells are still lacking. The present study investigated the anticancer effects of a root fraction of E. indica (R-S5-C1-H1) on H1299, MCF-7, and SK-HEP-1 cell lines and analyzed metabolic changes in the treated cancer cells using ultra-high-performance liquid chromatography high-resolution mass spectrometry (UHPLC-HRMS). Cell metabolic activity assays demonstrated that the cell viability of the three cancer cell lines was significantly reduced following treatment with R-S5-C1-H1, with half-maximal inhibitory concentrations values of 12.95 µg/mL, 15.99 µg/mL, and 13.69 µg/mL at 72 h, respectively. Microscopy analysis using Hoechst 33342 and Annexin V fluorescent dyes revealed that cells treated with R-S5-C1-H1 underwent apoptotic cell death, while chemometric analysis suggested that apoptosis was triggered 48 h after treatment with R-S5-C1-H1. Deconvoluted cellular metabolomics revealed that hydrophobic metabolites were significantly altered, including triacylglycerols, phosphatidylcholine, phosphatidylethanolamine, sphingomyelin, and ceramide, suggesting that apoptosis induction by R-S5-C1-H1 potentially occurred through modulation of phospholipid synthesis and sphingolipid metabolism. These metabolomic profiling results provide new insights into the anticancer mechanisms of E. indica and facilitate the overall understanding of molecular events following therapeutic interventions.

In the continuing battle of humans against cancer, which dates as far back as 3000 BC, the earliest written account of cancer is the description of breast cancer found in the Edwin Smith Papyrus1. Currently, cancer remains one of the leading causes of premature death among people aged between 30 and 69 years worldwide2,3. It is estimated that there were 19.3 million new cases in 20202 and almost 10.0 million deaths from cancer. Furthermore, the incidence of all cancers in Malaysia has been projected to almost double by 2040, from 48,639 to 86,666 new cases (2020)." href="/articles/s41598-022-13575-6#ref-CR4" id="ref-link-section-d279179327e580">4. Despite recent therapeutic advancements5, this disease continues to affect the quality of life of patients due to the limitations of current cancer treatments.

According to the World Health Organization, there are still many people who rely on herbal medicine as their primary healthcare resource6. Therefore, plants have an important role to play in the healthcare revolution. In traditional Chinese medicine, herbs have been used to prevent and treat a wide range of diseases7. For instance, the authorities in China encouraged the combined use of traditional Chinese medicine and Western medicine to protect and treat patients infected with the highly infectious novel coronavirus (SARS-CoV2), an approach that has proven to be effective to date8. In plants, secondary metabolites are distinct from the components of primary metabolism as they are not involved in regulatory metabolism; rather, they function to protect the plant against predation and microorganism invasion9. Decades of phytochemical research on over 150,000 plant species have revealed that secondary metabolites including phenolics, flavonoids, terpenoids, alkaloids, and sulfur-containing compounds, exert a variety of interesting biological effects10. Since 1981, nearly three-quarters (64.9 %) of prescribed anticancer drugs have been plant-derived11. Thus, medicinal plants are regarded as important and reliable sources for the discovery of anticancer drugs12,13,14.

Metabolomics is the comprehensive study of a complete set of small metabolites (<1500 Da) within a biological sample15. This technique broadens the systems biology view of organisms by bridging the gap between genotype and phenotype16,17. Metabolomics is widely applied in studies on plants, animals, medicine, and food. Among the technologies applied in metabolomics, including nuclear magnetic resonance (NMR), gas chromatography (GC), and liquid chromatography (LC) coupled with mass spectrometry (MS), LC-MS is the technique of choice for detecting and determining the elemental composition and molecular formula of an analyte of interest, as it can be performed without a derivatization step. Furthermore, the soft ionization and non-extensive heat application in LC for chromatographic separation reduce the chance of compound degradation during profiling18. Endogenous and exogenous metabolic characterizations of therapeutic agents using high-throughput metabolomics allow researchers to evaluate the efficacy and safety of these agents19. In cancer research, metabolomic studies have provided comprehensive information and improved the understanding of the underlying mechanisms of cancer pathogenesis and drug effects through the assessment of metabolic changes in cancer cells20,21. Thus, identification and quantification of extracted metabolites enable the monitoring of responses to external stimuli in test samples22,23.

Eleusine indica, also known as wiregrass or goosegrass, is one of the six diploid species of the genus Eleusine belonging to the family Poaceae24,25. It is commonly found in tropical and sub-tropical regions26,27 and is known to possess depurative, febrifugal, diuretic, and laxative properties25,28,29. Traditionally, the plant has been frequently used to treat hypertension, influenza, oliguria, and urine retention25,30,31, while decoctions of the whole plant are commonly used as antihelminthics and febrifuges32. The seeds are sometimes used as a famine food, as well as for the treatment of liver complaints26. Extracts of E. indica have been reported to exert antiproliferative effects against several types of human cancer cells including Rhabdomyosarcoma cells33, lung cancer cells (A549), and cervical cancer cells (HeLa)34. Although previous research has revealed the anticancer potential of E. indica, there remains a lack of scientific evidence regarding the effectiveness of the plant extracts against other human cancer cell lines, as well as metabolomic studies evaluating the metabolite changes in cancer cells induced by E. indica.

Therefore, this study aimed to i) evaluate the effects of E. indica root fractions on cell viability and cellular morphology in human non-small cell lung carcinoma (H1299), breast adenocarcinoma (MCF-7), and liver adenocarcinoma (SK-HEP-1) cells; and ii) examine the influence of the most active fraction on the metabolomic profiles of these cells using an LC-MS-based approach.

The experimental flowchart is illustrated in Fig. 1. Briefly, E. Indica was harvested and prepared for the extraction of compounds to be tested against different cancer cell lines. Half-maximal inhibitory concentrations (IC50) of the fractions were determined for each cell line using cell viability assays conducted 24, 48 and 72 h after treatment. Subsequently, staining with Hoechst 33342 and Annexin V staining was performed to determine the mechanism of cell death. Concurrently, cell metabolomics was carried out in order to elucidate the cell death-triggering mechanisms; in brief, snap-frozen cells underwent a Bligh and Dyer extraction protocol with modifications to separate polar and non-polar metabolomes. Following this, chromatographic separation of the extracted metabolomes was carried out using pentafluorophenyl columns. Metabolome profiling was carried out at an m/z range of 50–1500 using heated electrospray ionization, in both positive and negative modes. Acquired data were pre-processed using MZmine 235, followed by chemometric analysis. Each statistically significant metabolome underwent further compound matching and analysis.

Overall experimental workflow.

Bioassay-guided screening of E. indica extract fractions against H1299, MCF-7, and SK-HEP-1 cells using colorimetric 3-(4,5-Dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assays showed that the root fraction R-S5-C1-H1 exhibited the highest level of activity against all cell lines (Figure S1); this fraction was therefore selected for further analysis. As shown in Fig. 2, treatment with increasing concentrations of R-S5-C1-H1 resulted in significantly reduced cell viability in all three cancer cell lines after 24 h, 48 h, and 72 h, with the most marked effects observed at a concentration of 25 µg/mL. For H1299 and MCF-7 cell lines, numbers of viable cells decreased gradually with treatment time, while SK-HEP-1 cells showed a sharp decrease in viability at 72 h. Table 1 shows the IC50 values for H1299, MCF-7, and SK-HEP-1 cells after treating for 24 h, 48 h, and 72 h with fraction R-S5-C1-H1; doxorubicin (Dox) was used as a positive control. The highest level of R-S5-C1-H1 activity after 24 h was observed in H1299 cells (IC50 = 20.73 µg/mL), followed by MCF-7 (IC50 = 29.89 µg/mL) and SK-HEP-1 (IC50 = 44.36 µg/mL) cells. Prolonging the treatment up to 72 h resulted in an approximately twofold decrease in IC50 values for H1299 and MCF-7 cell lines, and a more than threefold decrease for SK-HEP-1 cells.

Viability of H1299 (a), MCF-7 (b), and SK-HEP-1 (c) cells following treatment with varying concentrations (0–100 µg/mL) of fraction R-S5-C1-H1 for 24 h, 48 h, and 72 h. Values are expressed as mean ± SEM of three independent experiments; differences were considered significant at p ≤ 0.05 (*) and p ≤ 0.001 (**).

When observed under a light microscope, H1299, MCF-7, and SK-HEP-1 cells treated with R-S5-C1-H1 displayed notable cell shrinkage and loss of adherence, while the cells were also observed to be rounded, with irregular and rough surfaces (Fig. 3a–c ii–iv). In contrast, untreated controls for all three cell lines were adherent and exhibited smooth surfaces (Fig. 3a–c i). More noticeable morphological changes were observed when the cells were exposed to longer treatment, in agreement with the results of the cell viability assay.

Light microscopy observation of morphological changes, indicated by red arrows, in H1299 (a), MCF-7 (b), and SK-HEP-1 (c) cells treated with 12.5 µg/mL (iv) and 25 µg/mL (ii and iii) of R-S5-C1-H1 for 24 h (ii), 48 h (iii), and 72 h (iv). Untreated control cells appeared adherent with smooth surfaces (i).

To examine the ability of R-S5-C1-H1 to induce apoptotic cell death in H1299, MCF-7, and SK-HEP-1 cells, Hoechst 33342 and Annexin V fluorescent dyes were used. Hoechst 33342, a DNA-binding dye, was used to detect nuclear condensation and fragmentation, which are major characteristics of apoptosis. Meanwhile, Annexin V, a 36-kDa calcium-binding protein, binds to phosphatidylserine (PS), which becomes translocated from the inner plasma membrane to the external cell surface during apoptosis. As shown in Fig. 4, apoptosis was observed in all three cancer cell lines following treatment with the IC50 of R-S5-C1-H1 for 48 h. Morphological changes, namely nuclear condensation and PS externalization, were revealed by Hoechst 33342 staining (blue) and Annexin V staining (green), respectively. The treated cells exhibited bright blue apoptotic nuclei (Fig. 4a–c iii) and green fluorescence (Fig. 4a–c iv), whereas, in the untreated control group, most cell nuclei showed weak homogeneous blue staining (Fig. 4a–c i) and absence of green fluorescence (Fig. 4a–c ii).

Morphological changes in apoptotic H1299 (a), MCF-7 (b), and SK-HEP-1 (c) cells following treatment with IC50 of R-S5-C1-H1 for 48 h, observed by fluorescence microscopy. Hoechst 33342 staining (blue) and Annexin V staining (green) were used to detect apoptotic cells in untreated control (i and ii) and treatment (iii and iv) groups. Red arrows indicate nuclear condensation (bright blue), while yellow arrows indicate phosphatidylserine externalization.

Principal component analysis (PCA), an unsupervised pattern-recognition technique, was used to inspect data distribution and identify outliers throughout the analysis. For non-polar metabolomes profiled in positive ionization mode, PCA revealed distinct differences between treatment and control groups for H1299, MCF-7, and SK-HEP-1 cells after 48 h treatment with the IC50 of R-S5-C1-H1 (Fig. 5a–c i). However, in negative ionization mode, the metabolomes of the treatment and control groups of the three cell lines profiled were co-clustered, indicating that the profiled metabolomes showed no significant differences. Therefore, these data were excluded from the analysis (data not shown). Next, the cell metabolites that showed the greatest changes after R-S5-C1-H1 treatment were identified through the use of S-plots (Fig. 5a–c ii). The S-plots were constructed by comparing the treatment and control groups of each cell line via orthogonal projection to latent structure-discriminant analysis (OPLS-DA). The values of parameters R2Y and Q2 were used to describe the quality of the model, representing the model’s goodness-of-fit (R2Y) and predictability (Q2). The parameters of each model are listed in Supplementary Table S1. In each of the S-plots, the scatter plot p[1] represents the magnitude of each variable within a model, while p(corr)[1] represents the reliability of each variable (modeled correlation). The most significantly altered variables are plotted at either the lower left or upper right of the S-plot, and both are highly correlated with group separation. In this study, compounds at the lower left and upper right of the S-plots were selected as potentially relevant metabolites that were altered as a result of R-S5-C1-H1 treatment, based on a p(corr)[1] cutoff value of 0.8. These compounds are labeled as per Fig. 5a–c ii. The validity of the OPLS-DA model was demonstrated by a permutation test (Figure S2). S-plots (Fig. 5a–c ii) and heat map analysis (Fig. 6) revealed differing metabolomic profiles for control and treatment groups, suggesting that R-S5-C1-H1 altered the levels of non-polar (mainly lipid) metabolites in the cells. Even between the cancer cell lines, R-S5-C1-H1 affected the metabolites differently. In general, membrane lipids, such as phosphatidylcholines (PC), phosphorylethanolamines (PE), and sphingomyelins (SM) were dysregulated in cells treated with the IC50 of R-S5-C1-H1 for 48 h. In contrast, triacylglycerols (TGs), the main lipids used for energy storage in cells, and ceramide (Cer), an apoptosis inducer, were found to be upregulated in the treated H1299, MCF-7, and SK-HEP-1 cells. Details of the identified hydrophobic metabolites affected by treatment with R-S5-C1-H1 are summarized in Tables S4–S6. An overview of metabolic pathways in the R-S5-C1-H1-treated cell lines is shown in Fig. 7.

Principal component analysis (PCA) (i) and S-plots (ii) derived from the control and treated non-polar metabolomes of H1299 (a), MCF-7 (b), and SK-HEP-1 (c) cells analyzed in positive ionization mode. Eclipses show the mean ± SD of the sample replicates. For the PCA, red dots represent the control group’s metabolites, while green dots represent the treatment (IC50 R-S5-C1-H1, 48 h) group’s metabolites. The labeled metabolites are the significantly perturbed metabolomes in the S-plots.

Heat map analysis of the extracted non-polar metabolomes of H1299 (a), MCF-7 (b), and SK-HEP-1 (c) cells analyzed in positive ionization mode. Groups are indicated at the top of each map by red (control, n = 6) and green (treatment [IC50 R-S5-C1-H1, 48 h], n = 6). The maps show 50 of the most significantly perturbed non-polar metabolites based on Euclidean distances and Ward clustering. The rows represent metabolites, while the columns represent the samples. The metabolite concentrations are represented on a log scale. The scale bar represents the normalized intensity of features. Decreased and increased metabolites are displayed in blue and red, respectively.

Overview of the metabolic pathways in cell lines 48 h after treatment with R-S5-C1-H1 E.indica root fraction. SM sphingomyelin, PC phosphatidylcholine, LPC lysophosphatidylcholine, PE phosphatidylethanolamine, PS phosphatidylserine, TG triacylglyceride, DG diacylglyceride, Cer ceramide, GlcCer glucosylceramide. The figure was created with BioRender.com.

Secondary plant metabolites consist of numerous chemical compounds, including phenolics, alkaloids, saponins, terpenes, lipids, and complex carbohydrates, and are produced by plant cells via various metabolic pathways. These pathways are derived from primary metabolic pathways and involve secondary modifications such as deamination. Secondary metabolites have been reported to exert antibiotic, antifungal, and antiviral activities in order to protect plants from a wide range of pathogenic invasions. These compounds also possess other important biological properties, including anti-inflammatory, antibacterial, antioxidant, and selective anticancer30,36activities, which provide a scientific basis for the use of medicinal plants in traditional medicine.

In this study, the effect of the most active fraction of E. indica root extract, R-S5-C1-H1, on the viability of H1299, MCF-7, and SK-HEP-1 cells was evaluated using MTT assays. R-S5-C1-H1 was found to affect the viability of all the studied cancer cell lines in a dose- and time-dependent manner, in agreement with other studies investigating the anticancer activity of E. indica33,34 in different cell lines. At 72 h post-treatment, R-S5-C1-H1 IC50 values were below 20 µg/mL for all three cancer cell lines. This suggests that R-S5-C1-H1 was highly active against the tested cell lines according to the criteria established by the US National Cancer Institute and Geran protocol (where IC50 ≤ 20 µg/mL = highly active, IC50 21–200 µg/mL = moderately active, IC50 201–500 µg/mL = weakly active and IC50 ≥ 501 µg/mL = not active37). Cell viability was significantly reduced by higher concentrations of R-S5-C1-H1 or by longer treatment times. Prolonging treatment time while reducing drug dose can be a useful way of potentially minimizing undesirable side effects. A similar approach was applied in a clinical study conducted by Bishnoi et al. (2017) on recurrent, advanced, and metastatic cancer; lowering the chemotherapeutic dosage but prolonging therapy to achieve the desired outcomes extended patient survival while reducing undesirable side effects including nausea, neuropathy, fatigue, sepsis, and thrombocytopenia38. Similar results have been reported for advanced cancers, with improved quality of life achieved by application of this therapeutic approach39,40,41.

Cancer is defined as dysregulated cell proliferation combined with suppressed programmed cell death42. Anticancer agents that can induce programmed cell death in cancer cells, particularly apoptosis, while minimizing unwanted effects on the surrounding normal cells have received considerable attention because they offer a promising strategy for cancer prevention and treatment43. The majority of chemotherapy drugs currently in clinical use inhibit cell growth and induce apoptosis44. Cellular apoptosis is characterized by a series of typical morphological events, including membrane blebbing, cell shrinkage, chromatin condensation, nuclear fragmentation, fragmentation into membrane-bound apoptotic bodies, and translocation of membrane PS45,46. Microscopy analysis in the present study showed chromatin condensation, nuclear fragmentation, and binding of Annexin V dye to PS on the outer leaflet of the plasma membrane in the treated cells, suggesting that R-S5-C1-H1 was able to induce apoptotic cell death. However, further studies are needed to elucidate the underlying mechanism of apoptosis induced by E. indica.

According to the metabolome analysis in this study, several important non-polar lipid metabolites in the cancer cells were affected by treatment with R-S5-C1-H1. However, during polar metabolome extraction, residues from the culture medium (consisting of polar compounds, nutrients, and amino acids) become included during the process of snap freezing and extraction. Although concentrations of these residues may be low, they can introduce significant artifacts into polar metabolome profiling at the cellular level. Therefore, for polar metabolites, differences between control and E. indica-treated groups may be highly contributed by the culture medium instead of the responding metabolome. In negative ionization mode, co-clustering of the non-polar metabolomes of the treatment and control groups profiled was observed, indicating that the profiled metabolomes were not significantly different. Thus, in the current study, we focused on the non-polar metabolome in positive ionization mode. Lipids are small molecules with inherent hydrophobic properties. These molecules form key components of membranes, act as energy supplements, and play roles in inter- and intracellular signaling as well as in metabolic regulation. A mixture of complex lipids such as TGs, diacylglycerols (DGs), and monoacylglycerols (MGs), which consist of fatty acyl groups (varying in length and degree of unsaturation) esterified to glycerol backbones, is present in eukaryotic cells. During untargeted non-polar extract profiling using high-resolution mass spectrometry, TGs, DGs, and MGs can be distinguished from other lipid classes using positive electrospray ionization by the formation of ammoniated adduct ions, i.e., [M+ NH4]+ ions23,47. Collision-induced dissociation (CID) of the ammoniated ions yields an abundance of product ions corresponding to the loss of ammonia plus each of the fatty acyl groups as a free carboxylic acid47,48.

From our analysis, cellular energy homeostasis was suspected to be disrupted, as TGs were increased significantly. In vivo, it is possible to determine dysregulation of energy metabolism by analysis of glucose, pyruvate, and lactate concentrations in tissue. However, in vitro, monitoring of cultured cells based on these polar metabolites tends to be challenging as minute metabolite changes can be obscured by nutrients or supplements in the culture media, leading to false-positive results. Accumulation of TGs has been associated with cellular stress, including apoptosis49. Li et al. (2018) suggested that cells undergoing chemically-induced oxidative stress or apoptosis utilize TGs (with polyunsaturated fatty acyl chains) to protect against toxicity and limit cell death49. This phenomenon could explain the upregulation of TGs in all three cancer cell lines following treatment with R-S5-C1-H1 in the present study. Mitochondrial dysfunction and inhibition of lipid biosynthesis eventually lead cells to undergo apoptosis50.

Diacylglycerols play important roles as secondary messengers, activating proteins involved in various signaling cascades; thus, accumulation in biological systems is strictly regulated51. Phosphorylation of DGs to form phospholipids via diacylglycerol kinase leads to the attenuation of DG levels in the cell membrane, altering the availability of intracellular DG signaling protein and regulating cell growth, trafficking, differentiation, and migration51. During phospholipid synthesis, phosphocholine/phosphoethanolamine is transferred from cytidine diphosphate (CDP)-choline/ethanolamine to DGs by phosphotransferase, producing phosphatidylcholine (PC)/phosphatidylethanolamine (PE)52. Phosphatidylcholines account for nearly 50% of glycerophospholipid species in eukaryotic cells. The abundance of these phospholipids highlights their importance in biological and cellular systems as either membrane bilayer constituents or secondary messengers. Phosphorylcholine-containing lipid species are distinguished from other phospholipid classes using electrospray positive ionization, whereby CID of the protonated ions yields an abundant product ion, m/z 184, corresponding to a phosphorylcholine head group bonded at the sn-3 position of the glycerol backbone and fatty acyl substituents at sn-1 and/or sn-223,53. Meanwhile, PE, an aminophospholipid, is the second most abundant phospholipid in mammalian membranes after PC54. It is concentrated on the cytosol-facing leaflet of the cell membrane, whereas phosphorylcholine-containing lipids (including PC and SM) are localized in the outer plasma membrane55. Phosphatidylethanolamine plays an important role in post-translational modifications, influencing membrane topology and promoting cell and organelle membrane fusion, oxidative phosphorylation, mitochondrial biogenesis, and autophagy56. Phosphorylethanolamine species are distinguished from other phospholipid classes by CID of the protonated ions, yielding an abundant daughter ion, [M+ H-141]+, with neutral loss of m/z 141, corresponding to the phosphorylethanolamine head group bonded at the sn-3 position57. Inhibition of the CDP-choline pathway is a common feature of apoptosis58, while perturbation of PC levels leads to the induction of apoptotic cell death irrespective of the PC synthesis steps that are impaired59. From our analysis, levels of DG, PC, and PE were significantly perturbed in the treated cells, suggesting that R-S5-C1-H1 might induce apoptosis by modifying phospholipid synthesis.

Sphingolipids, amphipathic lipids with diverse molecular structures and functions, act as key regulators in cancer cell survival and death60. These molecules all contain a long-chain (sphingoid) base (aliphatic amino alcohol) consisting of a long hydrocarbon tail with hydroxyl groups at sn-1 and sn-3 positions and an amine group at sn-2. The functions attributed to bioactive sphingolipids include cell growth, cell death, senescence, inflammation, immune responses, nutrient uptake, metabolism, responses to stress stimuli, and autophagy61,62. Collision-induced dissociation of protonated Cer species produces a unique fragment, m/z 264, corresponding to the sphingosine backbone (d18:1), whereas dihydroceramide (d18:0) and sphingadienine (d18:2) contribute fragments of m/z 266 and 262, respectively.

Ceramides are central molecules in sphingolipid metabolism and are composed of a sphingosine backbone and a long fatty acid chain (varying in length from C14 to C26) amide-bonded to the second carbon of the backbone. They serve as metabolic and structural precursors for complex sphingolipids (which contain hydrophilic head groups) such as SM, ceramide-1-phosphate, and carbohydrate-conjugated Cers63. Regulation of Cer is maintained by complex and integrated metabolic pathways involving specialized enzymes63. Sphingomyelinase (SMase) is an enzyme that hydrolyzes SM to yield Cer in the cytoplasm and Golgi apparatus respectively64, while the choline head group becomes bonded to DG to form PC. Cer has been reported to be synthesized de novo from the common sphinganine (d18:0) backbone65. The unsaturated sphingosine backbone (d18:1 or d18:2) enhances intermolecular hydrogen bonding in the polar region, tightening Cer molecule packing in the cellular membrane62. Importantly, Cer has the capability to assemble channels in the outer mitochondrial membranes to promote the release of cytochrome C for caspase-9 activation, leading to apoptosis66,67. It has been established that Cer is proapoptotic, with its generation preceding the onset of apoptotic signaling68. Many apoptotic stimuli, including cellular stress and cytotoxic drugs, have been shown to increase the levels of endogenous Cer through multiple mechanisms68, including inhibition of Cer metabolism enzymes69. In contrast, SM is mainly found in the Golgi apparatus70 due to the transport of Cer from the endoplasmic reticulum71. Upon arrival in the Golgi apparatus, Cer is converted into SM via the enzyme sphingomyelin synthase (SMS), which can induce cell proliferation, migration, and survival72. Metabolism of SM plays a role in cancer progression, whereby the breakdown of SM directly results in Cer production73. In addition, glucosylceramide (GlcCer), a pro-survival and antiapoptotic molecule generated from the transfer of glucose to Cer by glucosylceramide synthase (GCS), has been shown to be associated with multi-drug resistance in cancer cells69. Cers can be also produced from the degradation of GlcCers by the enzyme glucocerebrosidase (GCase)74. It has been suggested that the intracellular pool of Cer affects the pro-proliferative and antiapoptotic effects of GlcCer, whereby an increase in the synthesis of GlcCer decreases the level of Cer in cells75. In this study, the level of Cer was generally increased in the R-S5-C1-H1-treated cancer cells compared to controls, whereas levels of SM and GlcCer were decreased following treatment. These results suggest that R-S5-C1-H1 may induce apoptosis in the studied cancer cells by perturbing the balance of key sphingolipids through effects on sphingolipid metabolism.

In conclusion, R-S5-C1-H1 affected the viability of H1299, MCF-7, and SK-HEP-1 cells in a dose- and time-dependent manner. It also induced apoptotic cell death in the treated cells. Metabolomics based on LC–MS demonstrated that R-S5-C1-H1 caused different metabolic responses related to damage and rejuvenation in the cancer cell lines tested. Upregulation of TGs suggests protection of cells against R-S5-C1-H1 exposure, while perturbation of PC and PE levels may indicate induction of apoptosis by R-S5-C1-H1 through modulation of phospholipid synthesis. Dysregulated SM and increased levels of Cer species also suggest initiation of apoptosis, as these amphipathic lipids are responsible for the release of proapoptotic proteins from the mitochondria. Therefore, through high-throughput, unbiased profiling, we were able to screen for differences in multiple metabolites simultaneously. This method has advantages over conventional methods of analyzing the effects of therapeutic interventions, as it can facilitate the overall understanding of molecular events.

Wild E. indica was collected from the lowlands of Kota Belud (latitude 6.290833°N and longitude 116.424722°E), Sabah, Malaysia, in accordance with the Sabah Wildlife Conservation Enactment 1997. The plant was identified by Mr. Julius Kulip from the Institute for Tropical Biology and Conservation, Universiti Malaysia Sabah, and a voucher specimen (BORH 2263) was deposited at the Borneensis Herbarium of the Institute for Tropical Biology and Conservation, Universiti Malaysia Sabah, Kota Kinabalu, Malaysia.

The collected plant samples were thoroughly cleaned and stored at -20 °C prior to being freeze-dried using a Labconco freeze-dryer (Labconco, Kansas City, MO, USA). Freeze-dried E. Indica was ground using a heavy-duty blender. The homogenized sample was then exhaustively extracted using a modified extraction method23, whereby an optimized ratio of 1:1:1 (v/v/v) of double-distilled water (Milli-Q system (Merck, Darmstadt, Germany) at a resistivity of > 18.2 MΩ-cm/methanol (MeOH)/chloroform was employed. Following this, the lower layer of the mixture was removed and vacuum-concentrated to obtain a semisolid crude extract.

Crude extracts of the different plant parts (seeds, stem, leaf, and roots) were screened by being used to treat H1299, MCF-7, and SK-HEP-1 cell lines. Based on the crude extract screening (Fig. S1a), the root extract, which exhibited the highest activity against all cell lines, was selected for solid-phase extraction (SPE). An OASIS HLB Cartridge SPE column was used to purify the extract according to the solvent mixture stated in Table S2. From the SPE screening, R-S5 was the most active fraction (Fig. S1b); this was then subjected to normal-phase open-column chromatography in columns packed with silica gel with mesh size ranging between 0.063 and 0.200 mm (Merck, Darmstadt, Germany). Elution was initiated with hexane:ethyl acetate (9:1, v/v) and completed with MeOH (Table S3). The most active fraction, R-S5-C1 (Fig. S1c), was further separated using the Agilent 1200 series HPLC system (Agilent Technologies, Waldbronn, Germany) coupled with a Thermo Scientific Accucore C18 column (2.1 mm × 150 mm  ×  2.6 µm; Sunnyvale, CA, USA). Mobile phase A consisted of a mixture of deionized water with 0.1% formic acid and 1% ammonium acetate (NH4CH3CO2), while mobile phase B consisted of a mixture of acetonitrile and MeOH (6:4 v/v) with 0.1% formic acid and 1% NH4CH3CO2. The elution gradient was programmed to increase linearly from 1 to 70% of solvent B in 7 min, followed by 100% of solvent B from 7.1 to 10 min, maintained for 3 min. Later, the column was conditioned with the initial gradient for 1 min before the next sample injection. According to the results of the final screening, R-S5-C1-H1 was the most active fraction against all cell lines (Fig. S1d) and was therefore selected for further analyses. R-S5-C1-H1 was fully suspended in dimethylsulfoxide (DMSO), with a final concentration of DMSO ≤ 0.1% (v/v) prior to all cell experiments.

Non-small cell lung carcinoma H1299 (ATCC CRL-5803™), hepatocellular carcinoma SK-HEP-1 (ATCC HTB-52™), and breast adenocarcinoma MCF-7 (ATCC HTB-22™) cell lines were purchased from the American Type Culture Collection (Manassas, VA, USA). SK-HEP-1 cells were maintained in Dulbecco’s Modified Eagle’s Medium, while the other cell lines were cultured in Roswell Park Memorial Institute 1640 medium. All media were supplemented with 10% fetal bovine serum (Tico Europe, Amstelveen, Netherlands) and 1% penicillin–streptomycin (Nacalai Tesque, Kyoto, Japan). The cells were grown in a humidified environment with 5% CO2 at 37 °C.

The viability of H1299, MCF-7, and SK-HEP-1 cells was assessed using the colorimetric [3-(4,5-Dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide] (MTT) assay (Nacalai Tesque, Kyoto, Japan). Briefly, cells were seeded at 5  ×  103 cells/well (100 μL) on a 96-well plate (Nunc, Thermo Fisher Scientific, Denmark) and incubated for 24 h to allow cell attachment. After 24 h, the cells were treated with varying concentrations of R-S5-C1-H1, ranging from 0.05 to 100 μg/mL, for 24 h, 48 h, and 72 h. Doxorubicin (0.01 to 100 μg/mL) served as a positive control. The formazan precipitate in each well was suspended in DMSO, and absorbance was measured at 570 nm using the Tecan Infinite 200 Pro microplate reader (Männedorf, Switzerland). Percentage cell viability (%) was calculated as \(\frac{Absorbance\, of\, sample}{Absorbance\, of\, control }\times 100\%\). All experiments were performed in triplicate, and data were presented as mean ± standard error of the mean (SEM). Statistically significant differences at p ≤ 0.05 significance level were determined using the t-test (*p ≤ 0.05; **p ≤ 0.001).

All cells were seeded at 2.5  ×  105 cells/well (500 µL) in a 4-well µ-slide (ibidi, Martinsried, Germany) and incubated overnight before being treated with the IC50 of R-S5-C1-H1 for 48 h. Treated cells were fixed and stained concurrently with Hoechst 33342 at 1 μg/mL (Sigma, USA) and Annexin V at 10 μg/mL (Nacalai Tesque, Kyoto, Japan) at 37 °C for 30 min in the dark. Morphological observations were performed using a fluorescence microscope, BX 53 (Olympus, Tokyo, Japan), to visualize nuclear condensation and PS externalization with ultraviolet and green fluorescence filters, respectively (Olympus, Tokyo, Japan).

Each cell line was seeded at 1  ×  106 cells/dish in 12 culture dishes. Six dishes of each cell line were treated with the IC50 of R-S5-C1-H1 for 48 h, while the remaining six dishes served as untreated controls. The culture medium from each culture dish was collected and centrifuged at 1500 g to acquire a cell pellet. Adhered cells were snap-frozen with pre-cooled HPLC grade MeOH, collected by scraping, and then transferred into clean disposable borosilicate tubes together with the cell pellets from the previous step. The combined cells were extracted using a modified Bligh and Dyer extraction protocol65. Briefly, cells were mixed thoroughly with 4 mL MeOH/chloroform (1:1 v/v), followed by 2 mL of 0.05 M sodium chloride, and centrifuged at 1500 g at 4 °C for 30 min. Both polar and non-polar of cell extracts were transferred separately and vaporized using a speed vacuum concentrator system (Eppendorf, Hamburg, Germany) before being stored at -80 °C. Prior to LC–MS/MS analysis, the cell extracts were redissolved in 1.2 mL MeOH. Each cell extract was profiled in triplicate to ensure the consistency of the instrument.

Cell extract profiling was carried out according to a previously developed method76,77. Briefly, 10 μL or 30 μL of extract were applied to a Vanquish UHPLC system (Thermo Scientific, Waltham, MA, USA) coupled to the ultra-high resolution Qq-Time-of-flight Impact II mass spectrometer (Bruker, Billerica, MA, USA) in positive and negative electrospray ionization modes, respectively. A pentafluorophenyl column, Kinetex F5 (2.1 mm × 100 mm × 2.6 µm; Phenomenex, Torrance, California, USA), was used for chromatographic separation at 35 °C and the mobile phase flow rate was maintained at 0.6 mL/min. Mobile phase A consisted of a mixture of deionized water with 0.1% formic acid and 1% NH4CH3CO2, and mobile phase B consisted of a mixture of acetonitrile and MeOH (6:4 v/v) with 0.1% formic acid and 1% NH4CH3CO2. The gradient elution was programmed to increase linearly from 1 to 70% of solvent B in 7 min, followed by 100% solvent B from 7.1 to 10 min, maintained for 3 min. Following this, the column was conditioned with the initial gradient for 1 min before the next sample was introduced. Data acquisition was set between m/z 50 and 1500, while positive and negative electrospray ionization voltages were set at 3.5 kV and − 3.5 kV, respectively. Ion source gas temperature was set at 325 °C, drying gas flow at 10 L/min, and nebulizer flow at 3 Bar. The mass spectrometer was calibrated with Tune Mix (Sigma-Aldrich, St Louis, MO, USA) before analysis of each batch. A mass calibrant, sodium formate, was introduced between 0.1 and 0.3 min during each acquisition. Post-acquisition analyte m/z values were calibrated against the introduced sodium formate.

Acquired raw data were pre-processed using MZmine 235 to compensate for the variation in retention times and m/z values in each analysis. These pre-process data were exported as a peak list in comma-separated values (CSV) format, with rows representing the integrated and normalized peak areas and columns representing the samples. Exported .csv files were used for multivariate analyses with Metaboanalyst 5.0 (version 5.0; http://www.metaboanalyst.ca)78. Metabolite features with missing values > 50% were removed, and missing values were replaced by 1/5 of minimum positive values of the corresponding variables by default. The acquired data sets were log-transformed and Pareto-scaled prior to multivariable analyses to examine data distribution and identify outliers. Multivariate analyses of data including PCA, OPLS-DA, and hierarchical cluster analysis (heat map) were obtained. Subsequently, for OPLS-DA, the model examined internal relationships in matrix X-variables and response matrix Y. Model quality was estimated by R2X or R2Y values and Q2 values. To avoid OPLS-DA model over-fitting, an additional cross-validation tool (permutation test with Y variables randomized 100 times) was performed for each model. Goodness-of-fit parameters for the OPLS-DA model, R2X, R2Y, and Q2 were calculated (varying from 0 to 1). The R2Y and Q2 in the permutated test described the fitness of the data and the predictability of the derived model, respectively. In addition, the metabolite features (non-parametric t-test, p < 0.05) and corresponding fold changes (fold change threshold ≥ 2) showed how the selected differential metabolites varied between the treatment and control groups; these metabolites were selected and underwent further compound matching and analysis. In the current study, all metabolomics data analyses focused on the non-polar layer only and were performed on the positive ionization mode data sets, as the negative ionization mode revealed no significant differences (data excluded).

Based on the normalized acquired data, statistical significance was determined by t-tests (non-parametric analysis). Identification of metabolites was based on ion fragmentation pattern recognition23,79 and fragmentation spectra (MS/MS) matching using Metfrag80, Human Metabolome Database (HMDB)81 and Lipid Metabolites and Pathway Strategy (LIPIDMAPS)82. Mass-to-charge ratios complemented with fragmentation spectra and acceptable mass tolerance (at 5 ppm) revealed the metabolome identity. Furthermore, determination of lipid molecular formulae was achieved via the description of fatty acid substituents through mathematical calculation of mass balance for molecular weight based on high-resolution mass accuracy.

All data generated or analysed during this study are included in this published article (and its Supplementary Information file). Additional raw data files can be available from the corresponding authors upon reasonable request.

Hajdu, S. I. A note from history: Landmarks in history of cancer, part 1. Cancer 117, 1097–1102. https://doi.org/10.1002/cncr.25553 (2011).

Article PubMed Google Scholar

Sung, H. et al. Global cancer statistics 2020: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 71, 209–249. https://doi.org/10.3322/caac.21660 (2021).

Article PubMed Google Scholar

Wild, C.P., Weiderpass, E., Stewart, B.W. (eds). World cancer report: Cancer Research for Cancer Prevention., (Lyon, France: International Agency for Research on Cancer., Available from: http://publications.iarc.fr/586. Licence: CC BY-NC-ND 3.0 IGO., 2020).

Ferlay J et al. Global Cancer Observatory: Cancer Tomorrow. Lyon, France: International Agency for Research on Cancer, <Available from: https://gco.iarc.fr/tomorrow, accessed [11 Sep 2021].> (2020).

Advancing Cancer Therapy. Nature Cancer 2, 245-246, doi:https://doi.org/10.1038/s43018-021-00192-x (2021).

World Health, O. WHO traditional medicine strategy: 2014–2023. 76 p. Vietnamese version published by WPRO: http://iris.wpro.who.int/handle/10665.1/11340 The Thai version is published by the National Health Commission Office, Thailand (World Health Organization, 2013).

Zhao, Z. et al. Prevention and treatment of COVID-19 using traditional Chinese medicine: A review. Phytomedicine 85, 153308. https://doi.org/10.1016/j.phymed.2020.153308 (2021).

Article CAS PubMed Google Scholar

Wu, X. V., Dong, Y., Chi, Y., Yu, M. & Wang, W. Traditional Chinese Medicine as a complementary therapy in combat with COVID-19—A review of evidence-based research and clinical practice. J. Adv. Nurs. 77, 1635–1644. https://doi.org/10.1111/jan.14673 (2021).

Article PubMed Google Scholar

Kuruppu, A. I., Paranagama, P. & Goonasekara, C. L. Medicinal plants commonly used against cancer in traditional medicine formulae in Sri Lanka. Saudi Pharmac. J. 27, 565–573. https://doi.org/10.1016/j.jsps.2019.02.004 (2019).

Article Google Scholar

Aggarwal, B. B. & Heber, D. Immunonutrition: Interactions of Diet, Genetics, and Inflammation (1st ed.). (CRC Press., 2014).

Newman, D. J. & Cragg, G. M. Natural products as sources of new drugs over the nearly four decades from 01/1981 to 09/2019. J. Nat. Prod. 83, 770–803. https://doi.org/10.1021/acs.jnatprod.9b01285 (2020).

Article CAS PubMed Google Scholar

Atanasov, A. G. et al. Natural products in drug discovery: advances and opportunities. Nat. Rev. Drug Discovery 20, 200–216. https://doi.org/10.1038/s41573-020-00114-z (2021).

Article CAS PubMed Google Scholar

Kinghorn, A. D. Review of anticancer agents from natural products. J. Nat. Prod. 78, 2315–2315. https://doi.org/10.1021/acs.jnatprod.5b00617 (2015).

Article CAS Google Scholar

Huang, M., Lu, J.-J. & Ding, J. Natural products in cancer therapy: Past, present and future. Nat. Prod. Bioprospect. 11, 5–13. https://doi.org/10.1007/s13659-020-00293-7 (2021).

Article PubMed PubMed Central Google Scholar

Kuehnbaum, N. L. & Britz-McKibbin, P. New advances in separation science for metabolomics: Resolving chemical diversity in a post-genomic era. Chem. Rev. 113, 2437–2468. https://doi.org/10.1021/cr300484s (2013).

Article CAS PubMed Google Scholar

Fiehn, O. Metabolomics—the link between genotypes and phenotypes. Plant Mol. Biol. 48, 155–171. https://doi.org/10.1023/A:1013713905833 (2002).

Article CAS PubMed Google Scholar

Gao, D. et al. Metabolomics study on the antitumor effect of marine natural compound flexibilide in HCT-116 colon cancer cell line. J. Chromatogr. B 1014, 17–23. https://doi.org/10.1016/j.jchromb.2016.01.003 (2016).

Article CAS Google Scholar

Fang, M. et al. Thermal degradation of small molecules: A global metabolomic investigation. Anal. Chem. 87, 10935–10941. https://doi.org/10.1021/acs.analchem.5b03003 (2015).

Article CAS PubMed PubMed Central Google Scholar

Dhakshinamoorthy, S., Dinh, N.-T., Skolnick, J. & Styczynski, M. P. Metabolomics identifies the intersection of phosphoethanolamine with menaquinone-triggered apoptosis in an in vitro model of leukemia. Mol. BioSyst. 11, 2406–2416. https://doi.org/10.1039/C5MB00237K (2015).

Article CAS PubMed PubMed Central Google Scholar

Johnson, C. H., Ivanisevic, J. & Siuzdak, G. Metabolomics: Beyond biomarkers and towards mechanisms. Nat. Rev. Mol. Cell Biol. 17, 451–459. https://doi.org/10.1038/nrm.2016.25 (2016).

Article CAS PubMed PubMed Central Google Scholar

Schmidt, D. R. et al. Metabolomics in cancer research and emerging applications in clinical oncology. CA Cancer J. Clin. 71, 333–358. https://doi.org/10.3322/caac.21670 (2021).

Article PubMed Google Scholar

Sakagami, H. et al. Metabolomic profiling of sodium fluoride-induced cytotoxicity in an oral squamous cell carcinoma cell line. Metabolomics 10, 270–279. https://doi.org/10.1007/s11306-013-0576-z (2014).

Article CAS Google Scholar

Ling, Y. S., Lim, L. R., Yong, Y. S., Tamin, O. & Puah, P. Y. MS-based metabolomics revealing Bornean Sinularia sp. extract dysregulated lipids triggering programmed cell death in Hepatocellular carcinoma. Nat. Prod. Res. https://doi.org/10.1080/14786419.2018.1531288 (2018).

Article PubMed Google Scholar

Liu, Q., Triplett, J. K., Wen, J. & Peterson, P. M. Allotetraploid origin and divergence in Eleusine (Chloridoideae, Poaceae): Evidence from low-copy nuclear gene phylogenies and a plastid gene chronogram. Ann. Bot. 108, 1287–1298. https://doi.org/10.1093/aob/mcr231 (2011).

Article PubMed PubMed Central Google Scholar

Ette, O. E., Peace, M. E. & UbulomObot, D. A systematic review on Eleucine indica (L.) Gaertn.): From ethnomedicinal uses to pharmacological activities. J. Med. Plants Stud. 8, 262–274 (2011).

Google Scholar

Al-Zubairi, A. S. et al. Eleucine indica possesses antioxidant, antibacterial and cytotoxic properties. Evid.-Based Complement. Alternat. Med. 2011, 965370. https://doi.org/10.1093/ecam/nep091 (2011).

Article PubMed PubMed Central Google Scholar

Haber, R. M. & Semaan, M. T. Two new records from Lebanon: Chamaesyce nutans (Lag.) small (Euphorbiaceae) and Eleusine indica (L.) Gaertner (Poaceae). Turk. J. Bot. 31, 341–343 (2007).

Google Scholar

De Melo, G. O. et al. C-Glycosylflavones from the aerial parts of Eleusine indica inhibit LPS-induced mouse lung inflammation. Planta Med. 71, 362–363 (2005).

Article PubMed Google Scholar

Sagnia, B. et al. Antioxidant and anti-inflammatory activities of extracts from Cassia alata, Eleusine indica, Eremomastax speciosa, Carica papaya and Polyscias fulva medicinal plants collected in Cameroon. PLoS ONE 9, e103999. https://doi.org/10.1371/journal.pone.0103999 (2014).

Article CAS ADS PubMed PubMed Central Google Scholar

Ettebong, E. O., Nwafor, P. A. & Okokon, J. E. In vivo antiplasmodial activities of ethanolic exract and fractions of Eleucine indica. Asian Pac. J. Trop. Med. 5, 673–676. https://doi.org/10.1016/S1995-7645(12)60105-9 (2012).

Article CAS PubMed Google Scholar

Ong, S. L., Nalamolu, K. R. & Lai, H. Y. Potential lipid-lowering effects of Eleusine indica (L) Gaertn. Extract on high-fat-diet-induced hyperlipidemic rats. Pharmacogn. Mag. 13, S1-s9. https://doi.org/10.4103/0973-1296.203986 (2017).

Article PubMed PubMed Central Google Scholar

Kulip, J. A preliminary survey of traditional medicinal plants in the west coast and interior of Sabah. J. Trop. For. Sci. 10, 271–274 (1997).

Google Scholar

Ogbole, O. O., Segun, P. A. & Adeniji, A. J. In vitro cytotoxic activity of medicinal plants from Nigeria ethnomedicine on Rhabdomyosarcoma cancer cell line and HPLC analysis of active extracts. BMC Complement. Altern. Med. 17, 494. https://doi.org/10.1186/s12906-017-2005-8 (2017).

Article CAS PubMed PubMed Central Google Scholar

Hansakul, P., Wongnoppavich, A., Ingkaninan, K., Seewaboon, S. & Watcharin, P. Apoptotic induction activity of Dactyloctenium aegyptium (L.) PB and Eleusine indica (L.) Gaerth extracts on human lung and cervical cancer cell lines. Songklanakarin J. Sci. Technol. 31, 2 (2009).

Google Scholar

Pluskal, T., Castillo, S., Villar-Briones, A. & Orešič, M. MZmine 2: Modular framework for processing, visualizing, and analyzing mass spectrometry-based molecular profile data. BMC Bioinform. 11, 395. https://doi.org/10.1186/1471-2105-11-395 (2010).

Article CAS Google Scholar

Ghagane, S. C. et al. In vitro antioxidant and anticancer activity of Leea indica leaf extracts on human prostate cancer cell lines. Integr. Med. Res. 6, 79–87. https://doi.org/10.1016/j.imr.2017.01.004 (2017).

Article PubMed PubMed Central Google Scholar

Niksic, H. et al. Cytotoxicity screening of Thymus vulgaris L. essential oil in brine shrimp nauplii and cancer cell lines. Sci. Rep. 11, 13178. https://doi.org/10.1038/s41598-021-92679-x (2021).

Article CAS ADS PubMed PubMed Central Google Scholar

Bishnoi, R., Shah, C., Bejjanki, H., Bennett, J. A. & Reisman, D. N. An alternative approach with a low dose and prolonged chemotherapy for palliative treatment of locally advanced, metastatic or recurrent squamous cell head and neck cancer. Appl. Cancer Res. 37, 43. https://doi.org/10.1186/s41241-017-0049-1 (2017).

Article Google Scholar

Montagna, E. et al. Metronomic chemotherapy combined with bevacizumab and erlotinib in patients with metastatic HER2-negative breast cancer: clinical and biological activity. Clin. Breast Cancer 12, 207–214. https://doi.org/10.1016/j.clbc.2012.03.008 (2012).

Article CAS PubMed Google Scholar

Orlando, L. et al. Trastuzumab in combination with metronomic cyclophosphamide and methotrexate in patients with HER-2 positive metastatic breast cancer. BMC Cancer 6, 225. https://doi.org/10.1186/1471-2407-6-225 (2006).

Article CAS PubMed PubMed Central Google Scholar

Wong, N. S. et al. Phase I/II trial of metronomic chemotherapy with daily dalteparin and cyclophosphamide, twice-weekly methotrexate, and daily prednisone as therapy for metastatic breast cancer using vascular endothelial growth factor and soluble vascular endothelial growth factor receptor levels as markers of response. J. Clin. Oncol. 28, 723–730. https://doi.org/10.1200/jco.2009.24.0143 (2010).

Article CAS PubMed Google Scholar

Yuan, B. Z., Jefferson, A. M., Millecchia, L., Popescu, N. C. & Reynolds, S. H. Morphological changes and nuclear translocation of DLC1 tumor suppressor protein precede apoptosis in human non-small cell lung carcinoma cells. Exp. Cell Res. 313, 3868–3880. https://doi.org/10.1016/j.yexcr.2007.08.009 (2007).

Article CAS PubMed Google Scholar

Lopez, L. et al. Cupressus lusitanica (Cupressaceae) leaf extract induces apoptosis in cancer cells. J. Ethnopharmacol. 80, 115–120 (2002).

Article CAS PubMed Google Scholar

Fribley, A. M. et al. Proteasome inhibitor PS-341 induces apoptosis in cisplatin-resistant squamous cell carcinoma cells by induction of Noxa. J. Biol. Chem. 281, 31440–31447. https://doi.org/10.1074/jbc.M604356200 (2006).

Article CAS PubMed Google Scholar

Elmore, S. Apoptosis: A review of programmed cell death. Toxicol. Pathol. 35, 495–516. https://doi.org/10.1080/01926230701320337 (2007).

Article CAS PubMed PubMed Central Google Scholar

Mukhopadhyay, S., Panda, P. K., Sinha, N., Das, D. N. & Bhutia, S. K. Autophagy and apoptosis: Where do they meet?. Apoptosis 19, 555–566. https://doi.org/10.1007/s10495-014-0967-2 (2014).

Article CAS PubMed Google Scholar

Krank, J., Murphy, R. C., Barkley, R. M., Duchoslav, E. & McAnoy, A. In: Methods in Enzymology Vol. 432, 1–20 (Academic Press, 2007).

Google Scholar

Li, X., Collins, E. J. & Evans, J. J. Examining the collision-induced decomposition spectra of ammoniated triglycerides as a function of fatty acid chain length and degree of unsaturation. II. The PXP/YPY series. Rapid Commun. Mass Spectrometry RCM 20, 171–177. https://doi.org/10.1002/rcm.2270 (2006).

Article CAS ADS PubMed Google Scholar

Li, N., Sancak, Y., Frasor, J. & Atilla-Gokcumen, G. E. A protective role for triacylglycerols during apoptosis. Biochemistry 57, 72–80. https://doi.org/10.1021/acs.biochem.7b00975 (2018).

Article CAS PubMed Google Scholar

Boren, J. & Brindle, K. M. Apoptosis-induced mitochondrial dysfunction causes cytoplasmic lipid droplet formation. Cell Death And Differentiation 19, 1561–1570, doi:https://doi.org/10.1038/cdd.2012.34. https://www.nature.com/articles/cdd201234#supplementary-information (2012).

Topham, M. K. & Prescott, S. M. in Handbook of Cell Signaling (Second Edition) (ed Edward A. Dennis) 1177–1182 (Academic Press, 2010).

Fagone, P. & Jackowski, S. Phosphatidylcholine and the CDP–choline cycle. Biochimica et Biophysica Acta (BBA) - Molecular and Cell Biology of Lipids 1831, 523–532, doi:https://doi.org/10.1016/j.bbalip.2012.09.009 (2013).

Ling, Y. S. et al. Two-dimensional LC-MS/MS to enhance ceramide and phosphatidylcholine species profiling in mouse liver. Biomed. Chromatogr. BMC 28, 1284–1293. https://doi.org/10.1002/bmc.3162 (2014).

Article CAS PubMed Google Scholar

Rockenfeller, P. et al. Phosphatidylethanolamine positively regulates autophagy and longevity. Cell Death And Differentiation 22, 499, doi:https://doi.org/10.1038/cdd.2014.219. https://www.nature.com/articles/cdd2014219#supplementary-information (2015).

van Meer, G., Voelker, D. R. & Feigenson, G. W. Membrane lipids: where they are and how they behave. Nat. Rev. Mol. Cell Biol. 9, 112–124. https://doi.org/10.1038/nrm2330 (2008).

Article CAS PubMed PubMed Central Google Scholar

Calzada, E., Onguka, O. & Claypool, S. M. Phosphatidylethanolamine metabolism in health and disease. Int. Rev. Cell Mol. Biol. 321, 29–88. https://doi.org/10.1016/bs.ircmb.2015.10.001 (2016).

Article CAS PubMed Google Scholar

Schwudke, D. et al. Lipid profiling by multiple precursor and neutral loss scanning driven by the data-dependent acquisition. Anal. Chem. 78, 585–595. https://doi.org/10.1021/ac051605m (2006).

Article CAS PubMed Google Scholar

Morton, C. C., Aitchison, A. J., Gehrig, K. & Ridgway, N. D. A mechanism for suppression of the CDP-choline pathway during apoptosis. J. Lipid Res. 54, 3373–3384. https://doi.org/10.1194/jlr.M041434 (2013).

Article CAS PubMed PubMed Central Google Scholar

Bladergroen, B. A. et al. Inhibition of phosphatidylcholine and phosphatidylethanolamine biosynthesis in rat-2 fibroblasts by cell-permeable ceramides. Eur. J. Biochem. 264, 152–160. https://doi.org/10.1046/j.1432-1327.1999.00589.x (1999).

Article CAS PubMed Google Scholar

Moro, K. et al. Ceramide species are elevated in human breast cancer and are associated with less aggressiveness. Oncotarget 9, 2. https://doi.org/10.18632/oncotarget.24903 (2018).

Article Google Scholar

Hannun, Y. A. & Obeid, L. M. Sphingolipids and their metabolism in physiology and disease. Nat. Rev. Mol. Cell Biol. 19, 175–191. https://doi.org/10.1038/nrm.2017.107 (2018).

Article CAS PubMed Google Scholar

Mencarelli, C. & Martinez-Martinez, P. Ceramide function in the brain: When a slight tilt is enough. Cell. Mol. Life Sci. CMLS 70, 181–203. https://doi.org/10.1007/s00018-012-1038-x (2013).

Article CAS PubMed Google Scholar

Ponnusamy, S. et al. Sphingolipids and cancer: Ceramide and sphingosine-1-phosphate in the regulation of cell death and drug resistance. Fut. Oncol. (Lond., Engl.) 6, 1603–1624. https://doi.org/10.2217/fon.10.116 (2010).

Article CAS Google Scholar

Andrieu-Abadie, N. & Levade, T. Sphingomyelin hydrolysis during apoptosis. Biochem. Biophys. Acta. 1585, 126–134. https://doi.org/10.1016/S1388-1981(02)00332-3 (2002).

Article CAS PubMed Google Scholar

Ling, Y. S., Liang, H.-J., Chung, M.-H., Lin, M.-H. & Lin, C.-Y. NMR- and MS-based metabolomics: various organ responses following naphthalene intervention. Mol. BioSyst. 10, 1918–1931. https://doi.org/10.1039/C4MB00090K (2014).

Article CAS PubMed Google Scholar

Young, M. M., Kester, M. & Wang, H. G. Sphingolipids: Regulators of crosstalk between apoptosis and autophagy. J. Lipid Res. 54, 5–19. https://doi.org/10.1194/jlr.R031278 (2013).

Article CAS PubMed PubMed Central Google Scholar

Schüll, S. et al. Cytochrome c oxidase deficiency accelerates mitochondrial apoptosis by activating ceramide synthase 6. Cell Death Dis. 6, e1691. https://doi.org/10.1038/cddis.2015.62 (2015).

Article CAS PubMed PubMed Central Google Scholar

Huang, W.-C., Chen, C.-L., Lin, Y.-S. & Lin, C.-F. Apoptotic sphingolipid ceramide in cancer therapy. J. Lipids 2011, 565316. https://doi.org/10.1155/2011/565316 (2011).

Article CAS PubMed PubMed Central Google Scholar

Gucluler, G., Piskin, O. & Baran, Y. The roles of antiapoptotic sphingosine kinase-1 and glucosylceramide genes in drug induced cell death of MCF-7 breast cancer cells. J. BUON 16, 646–651 (2011).

CAS PubMed Google Scholar

D'Avanzo, N. in Current Topics in Membranes Vol. 78 (eds Robert J. French & Sergei Yu Noskov) 353–407 (Academic Press, 2016).

Simons, K., Klose, C. & Surma, M. in Lipidomics 21–34 (2012).

Taniguchi, M. & Okazaki, T. Ceramide/sphingomyelin rheostat regulated by sphingomyelin synthases and chronic diseases in murine models. J. Lipid. Atheroscler. 9, 380–405. https://doi.org/10.12997/jla.2020.9.3.380 (2020).

Article PubMed PubMed Central Google Scholar

Tallima, H., Azzazy, H. M. E. & El Ridi, R. Cell surface sphingomyelin: Key role in cancer initiation, progression, and immune evasion. Lipids Health Dis. 20, 150. https://doi.org/10.1186/s12944-021-01581-y (2021).

Article CAS PubMed PubMed Central Google Scholar

Zigdon, H., Meshcheriakova, A. & Futerman, A. H. From sheep to mice to cells: Tools for the study of the sphingolipidoses. Biochim. Biophys. Acta Mol. Cell Biol. Lipids 1189–1199, 2014. https://doi.org/10.1016/j.bbalip.2014.02.008 (1841).

Article CAS Google Scholar

Reza, S., Ugorski, M. & Suchański, J. Glucosylceramide and galactosylceramide, small glycosphingolipids with significant impact on health and disease. Glycobiology https://doi.org/10.1093/glycob/cwab046 (2021).

Article PubMed PubMed Central Google Scholar

Yoong-Soon, Y., Jiann, C. E. T., Hsin-Chang, C., Ping-Chin, L. & Soon, L. Y. A comparative study of pentafluorophenyl and octadecylsilane columns in high-throughput profiling of biological fluids. J. Chin. Chem. Soc. 64, 699–710. https://doi.org/10.1002/jccs.201600873 (2017).

Article CAS Google Scholar

Puah, P. Y. et al. Extractable impurities from fluoropolymer-based membrane filters—interference in high-throughput, untargeted analysis. RSC Adv. 9, 31918–31927. https://doi.org/10.1039/c9ra06198c (2019).

Article CAS ADS Google Scholar

Pang, Z. et al. MetaboAnalyst 5.0: Narrowing the gap between raw spectra and functional insights. Nucleic Acids Res. 49, 388–396. https://doi.org/10.1093/nar/gkab382 (2021).

Article CAS Google Scholar

Wong, C., Ling, Y. S., Wee, J. L. S., Mujahid, A. & Müller, M. A comparative UHPLC-Q/TOF–MS-based eco-metabolomics approach reveals temperature adaptation of four Nepenthes species. Sci. Rep. 10, 21861. https://doi.org/10.1038/s41598-020-78873-3 (2020).

Article CAS ADS PubMed PubMed Central Google Scholar

Ruttkies, C., Schymanski, E. L., Wolf, S., Hollender, J. & Neumann, S. MetFrag relaunched: Incorporating strategies beyond in silico fragmentation. J. Cheminform. 8, 3. https://doi.org/10.1186/s13321-016-0115-9 (2016).

Article CAS PubMed PubMed Central Google Scholar

Wishart, D. S. et al. HMDB 4.0: the human metabolome database for 2018. Nucleic Acids Res. 46, D608-D617. https://doi.org/10.1093/nar/gkx1089 (2018).

Sud, M. et al. LMSD: LIPID MAPS structure database. Nucleic Acids Res. 35, D527–D532. https://doi.org/10.1093/nar/gkl838 (2007).

Download references

This work was supported by Universiti Malaysia Sabah (UMS) under the UMSGreat research grants GUG 0402-2/2019 and GUG 0278-2/2018. Many thanks to Mr. Samudi bin Surag and Miss Lim Leong Lui for their kind help throughout the experiments.

These authors contributed equally: Perng Yang Puah and Dexter Jiunn Herng Lee.

Department of Biomedical Sciences, Faculty of Medicine and Health Sciences, Universiti Malaysia Sabah, Jalan UMS, 88400, Kota Kinabalu, Sabah, Malaysia

Perng Yang Puah & Siat Yee Fong

Biotechnology Research Institute, Universiti Malaysia Sabah, Jalan UMS, 88400, Kota Kinabalu, Sabah, Malaysia

Dexter Jiunn Herng Lee & Yee Soon Ling

Medical Department, Sarawak General Hospital, Jalan Hospital, 93586, Kuching, Sarawak, Malaysia

Soo Huan Puah

Medical Department, Hospital Seberang Jaya, Jalan Tun Hussein Onn, Seberang Jaya, 13700, Permatang Pauh, Penang, Malaysia

Soo Huan Puah

Department of Surgery, Faculty of Medicine and Health Sciences, Universiti Malaysia Sabah, Jalan UMS, 88400, Kota Kinabalu, Sabah, Malaysia

Nik Amin Sahid Nik Lah

CAIQ Certification Sdn Bhd Kota Kinabalu, Sabah, Malaysia

Yee Soon Ling

Borneo Medical and Health Research Centre, Faculty of Medicine and Health Sciences, Universiti Malaysia Sabah, Jalan UMS, 88400, Kota Kinabalu, Sabah, Malaysia

Siat Yee Fong

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

You can also search for this author in PubMed Google Scholar

P.Y.P., D.J.H.L., Y.S.L., and S.Y.F. conceived the experiments, P.Y.P., D.J.H.L. and Y.S.L. conducted experiments and data analysis and wrote the manuscript, S.H.P. and N.A.S. contributed to statistical analysis, and Y.S.L and S.Y.F. supervised the work and designed the experiments. All authors discussed the data and reviewed the manuscript.

Correspondence to Yee Soon Ling or Siat Yee Fong.

The authors declare no competing interests.

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and Permissions

Puah, P.Y., Lee, D.J.H., Puah, S.H. et al. High-throughput metabolomics reveals dysregulation of hydrophobic metabolomes in cancer cell lines by Eleusine indica. Sci Rep 12, 9347 (2022). https://doi.org/10.1038/s41598-022-13575-6

Download citation

Received: 13 December 2021

Accepted: 18 May 2022

Published: 06 June 2022

DOI: https://doi.org/10.1038/s41598-022-13575-6

Anyone you share the following link with will be able to read this content:

Sorry, a shareable link is not currently available for this article.

Provided by the Springer Nature SharedIt content-sharing initiative

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.