Screening Immunoactive Compounds of Ganoderma lucidum Spores by Mass Spectrometry Molecular Networking Combined With in vivo Zebrafish Assays

Ganoderma lucidum is a well-known herbal remedy widely used for treating various chronic diseases. Traditionally, the fruiting body is regarded as the medicinal part of this fungus, while recently, the therapeutic potentials of Ganoderma lucidum spore (GLS) is gaining increasing interests. However, detailed knowledge of chemical compositions and biological activities of the spore is still lacking. In this study, high-resolution mass spectrometry and molecular networking were employed for in-depth chemical profiling of GLS, sporoderm-broken GLS (BGLS) and sporoderm-removed GLS (RGLS), leading to the characterization of 109 constituents. The result also showed that RGLS contained more triterpenoids with much higher contents than BGLS and GLS. Moreover, the immunomodulatory activities of BGLS and RGLS were investigated in the zebrafish models of neutropenia or macrophage deficiency. RGLS exhibited more potent activities in alleviating vinorelbine-induced neutropenia or macrophage deficiency, and significantly enhanced phagocytic function of macrophages, which indicated the immunomodulatory activity of GLS was positively correlated with the content of triterpenoids. Further correlation analysis of chemical profiles of GLS and corresponding bioactivities by partial least squares regression identified the potential immunoactive compounds of GLS, including 20-hydroxylganoderic acid G, elfvingic acid A and ganohainanic acid C. Our findings suggest that combining mass spectrometry molecular networking with zebrafish-based bioassays and chemometrics is a feasible strategy to reveal complex chemical compositions of herbal medicines, as well as to discover their potential active constituents.

Introduction

Ganoderma lucidum, commonly known as Lingzhi or Reishi, has been used as an herbal remedy in China and many Asian countries for over 2000 years (Zhou et al., 2015; Yuan et al., 2018). According to traditional Chinese medicine (TCM) theory, G. lucidum can tonify “Qi,” and has been revered for its miracle cures and general health promoting benefits (Bishop et al., 2015). Modern scientific studies have proven that this medical macrofungus possesses various bioactivities, including immunomodulation, liver protection, diabetic treatment, anti-tumor and neuroprotective effects (Ahmad, 2018; Cao et al., 2018). Traditionally, the fruiting body of G. lucidum is used as the medicinal part and regarded as the source for many reported activities (Russell and Paterson, 2006; Hsu and Cheng, 2018). Less mature, but potentially even more valuable to therapeutic agent development, is the Ganoderma lucidum spore (GLS), the tiny reproduction unit of the fungus. Recently, GLS is gaining increasing acceptance and popularity as a functional food and nutraceutical, whose efficacy and safety have been suggested by multiple clinical studies in the treatment of cancers (Zhao et al., 2012; Hsu and Cheng, 2018), chronic periodontitis (Nayak et al., 2015) and Alzheimer disease (Wang et al., 2018). Although the use of GLS becomes popular, detailed knowledge of its chemical composition and biological activity is often lacking, as are data on the pharmacodynamics and clinical effects. Additionally, as GLS has outer bilayers of sporoderm, which is mainly composed of chitin and glucan (Lin and Wang, 2006), a variety of sporoderm-breaking techniques have been developed to release the components from the hard and resilient spores (Liu et al., 2005; Soccol et al., 2016). However, only a limited number of studies have been performed to investigate changes in chemical and biological properties of GLS after breaking the spore walls (Chen et al., 2012; Fu et al., 2012; Gao et al., 2013; Xu et al., 2014; Yang et al., 2017), and active constituents of GLS remain elusive (Liu et al., 2011; Yan et al., 2013).

Since its emergence, mass spectrometry (MS) is increasingly perceived as an essential tool in nearly all phases of drug discovery and development, including lead identification, metabolism, pharmacokinetics, and assessment of drug quality and safety (Hofstadler and Sannes-Lowery, 2006; Pacholarz et al., 2012). The hyphenated techniques, such as liquid chromatography-MS (LC-MS), and tandem MS (MS²), which represent the most widely used tools in MS arsenal, have shown many unique strengths in the drug discovery process. Recently, this cutting-edge technique has also been introduced into the realm of natural products and herbal medicines, which have been the source for new pharmaceutical drugs (Newman and Cragg, 2016). Different from synthetic or highly purified drugs, herbal medicines are complex mixtures, which usually contain hundreds of different phytochemicals. These herbal constituents generate thousands of molecular ions and fragment ions in MS analysis, rendering it challenging to annotate the detected chemical signatures. To address this issue, many MS data processing strategies have been developed to accelerate the dereplication and discovery process (Wang et al., 2016b, 2019; Li et al., 2017, 2019). Among these approaches, molecular networking (Watrous et al., 2012) is an emerging tool well suited to this task. Molecular networking visualizes all the ions and their chemical relationships based on MS² spectra similarity, which is calculated by using a cosine score (Guthals et al., 2012; Quinn et al., 2017). In the generated molecular networks, each node represents a consensus MS² spectrum (merged spectra with the same precursor ion and similar MS² spectra), while the edges between the nodes indicate the degree of cosine similarity. Based on the established Global Natural Product Social Molecular Networking (GNPS, https://gnps.ucsd.edu) web platform (Wang et al., 2016a), researchers can elucidate the structure of analogs and structurally related molecules based on their MS spectra (Yang et al., 2013; Allard et al., 2016). Despite the great potential to aid dereplication and structure elucidation, the utility of molecular networking in herbal medicine discovery just begins (Ge et al., 2017; Li et al., 2018; Qiang et al., 2019), and improvements in method performance (e.g., in terms of algorithm, and bioactivity relevance) are still necessary.

Besides its applications in chemical identification, MS has been employed to screen active constituents from complex herbal medicines. Many strategies for active (constituents) identification based on MS have been proposed, which can be separated into two primary categories: the affinity ultrafiltration based strategies which directly detect target-ligand complexes or “free” ligands released from the noncovalent complexes (Chen et al., 2018); and the chemometrics based strategies which investigate the active constituents by correlating chemical profiles of herbal medicines to their bioactive effects (Chang et al., 2018). For example, Yang et al. developed an ultrafiltration high-performance liquid chromatography coupled with diode array detector and mass spectrometry (UF-HPLC-DAD-MS) method to screen tyrosinase inhibitors from mulberry leaves (Yang et al., 2012). Besides, in one of our previous studies, we explored the active constituents of a Chinese medicine (Wenxin Keli) by combining LC-MS, bioassays and an active index approaches (Liu et al., 2018). It is noteworthy that both types of methods own their distinctive pros and cons, and meticulous validation (e.g., in silico docking, dose-response tests, and in vivo pharmacological assays) is necessary when candidates are obtained.

In this work, the active constituents of GLS were investigated by correlating chemical profiles of GLS produced by different manufacturing processes to their respective in vivo activities using partial least squares regression. Molecular networking was employed for structure elucidation and in-depth chemical profiling of GLS. Moreover, the immunomodulatory activities of GLS were evaluated by the zebrafish models of neutrophil or macrophage deficiency, as well as the phagocytic capability of macrophages. The efficiency rate of each constituent was then calculated and ranked based on its peak areas in different GLS samples and the corresponding bioactivities. The experimental workflow is outlined schematically in Figure 1.

FIGURE 1

Figure 1. Scheme of identifying immunoactive compounds of Ganoderma lucidum spores by molecular networking, zebrafish assays and chemometrics.

Materials and Methods

Materials and Reagents

Ganoderma lucidum spore (GLS, batch No. 14072201), sporoderm-broken Ganoderma lucidum spore (BGLS, batch No. 15121401) and sporoderm-removed Ganoderma lucidum spore (RGLS, batch No. 16042301) were provided by Zhejiang Shouxiangu Institute of Rare Medicine Plant (Zhejiang, China) and the species origin of the samples were authenticated as Ganoderma lucidum (Curtis) P. Karst. by Prof. Mingyan Li in the institute. The samples were kept in the sample room in Pharmaceutical Informatics Institute of Zhejiang University. Vinorelbine (batch No. 140501) was purchased from Jiangsu Haosen Pharmaceutical Co., Ltd (Jiangsu, China). Reference standards including ganoderic acid A, B, C₂, D, DM, F, G, H, ganoderenic acid A, B, C, D, lucidenic acid A were purchased from Nature Standard (Shanghai, China).

HPLC-grade acetonitrile and methanol were purchased from Merck (Darmstadt, Germany). Formic acid (HPLC grade) was purchased from Roe Scientific (Newark, DE, USA). Ethanol and Methylcellulose were acquired from Shanghai Aladdin Bio-Chem Technology (Shanghai, China). Neutral red was purchased from Sigma Aldrich (Saint-Louis, MO, USA). Isopropyl alcohol was purchased from Hangzhou Changzheng Chemical Reagent (Hangzhou, China). Deionized water was prepared with an Elga PURELAB flex system (ELGA LabWater, UK).

Sample Preparation

2 mg of GLS, BGLS, and RGLS were dissolved in 1 mL methanol respectively, and ultrasonically extracted for 20 min, and then centrifuged at 10,000 rpm for 10 min. The supernatants were collected for LC-MS analysis.

BGLS and RGLS were dissolved in system fish water at 10 mg/mL respectively, which were used for zebrafish assays. The system fish water was composed of 200 mg Instant Ocean Salt in per litter of reverse osmosis water with final pH 6.9–7.2, conductivity 480–510 μS/cm, and hardness of 53.7–71.6 mg/L CaCO₃.

Instrumentation

An Acquity UPLC system (Waters, Milford, MA, USA) coupled with a Triple TOF 5600plus MS (AB SCIEX, Framingham, MA, USA) was employed for chemical identification. Analysis was performed in negative electrospray ionization (ESI-) mode under following parameters: scan range, m/z 100–1500; source voltage, −4.5 kV; source temperature, 550°C; curtain gas, 35 psi; gas 1 (N₂), 50 psi; gas 2 (N₂), 50 psi. Declustering potential (DP), collision energy (CE) and collision energy spread (CES) of information dependent acquisition (IDA)-mediated MS² were 100 V, 40 eV, and 20 eV, respectively.

Chromatographic separation was carried out on an Acquity BEH C18 column (100 mm × 2.1 mm, 1.7 μm, Waters) at 30°C with mobile phase A (0.1% v/v formic acid-water) and mobile phase B (acetonitrile). The flow rate was 0.4 mL/min and a linear gradient elution was programmed: 0–20 min, 20–40% B; 20–28 min, 40–50% B; 28–35 min, 50–70% B; 35–40 min, 70–80% B; 40–45 min, 80–90% B; 45–50 min, 90–95% B; 50–55 min, 95% B. The injection volume was 2 μL.

Construction of MS/MS Based Molecular Network

Tandem mass spectrometry molecular networks were generated using the GNPS platform (https://gnps.ucsd.edu/). Raw MS data were first converted to mzXML format with MSConvert (Kessner et al., 2008) and then uploaded to GNPS to create the molecular networks. The precursor ion mass tolerance was set to 0.5 Da and to a product ion tolerance of 0.1 Da. A network was constructed using 6 minimum matched peaks and a cosine score above 0.7. The spectra in the network were searched against the spectral libraries on GNPS. Results were open and visualized in Cytoscape 3.7.1.

Zebrafish Husbandry and Management

Tg (mpx:GFP) transgenic zebrafish that expressed GFP exclusively in neutrophils and Albino zebrafish were provided by Hunter Biotechnology, which is accredited by the International Association for Assessment and Accreditation of Laboratory Animal Care (AAALAC). Embryos were generated by natural pair-wise mating, and anesthetized in 0.016% (w/v) tricaine prior to observations.

Zebrafish Model of Neutropenia

2-dpf Tg (mpx:GFP) transgenic zebrafish embryos were distributed into 6-well plates, with 30 larvae in 3 mL system fish water for each well. Three groups, i.e., the control group, the model group and the treatment group, were set up. Vinorelbine was administered at 1 ng per larva by intravenous microinjection to generate the zebrafish model of neutropenia. The treatment group was incubated in BGLS or RGLS supplemented fish water after microinjection. The final concentrations of BGLS (22, 67, and 200 μg/mL) and RGLS (33, 100, and 300 μg/mL) were set according to the maximum tolerated concentrations (MTCs) assay (Supplementary Methods). All the groups were incubated in a 28°C incubator for 24 h. Ten larvae were randomly selected from each group and the numbers of neutrophils in the zebrafish were counted and recorded with a Nikon Multi-purpose Zoom Microscope AZ100. The neutrophil recovery rate was calculated using the following formula:

$\begin{array}{l}\text{neutrophil recovery rate}=\frac{N_{treatment}-N_{model}}{N_{model}}\times 100\%\end{array}$

where N_treatment and N_model were the numbers of neutrophils of the larvae in the treatment group and model group, respectively.

Zebrafish Model of Macrophage Deficiency

2-dpf Albino zebrafish embryos were distributed into 6-well plates, with 30 larvae in each well with 3 mL system fish water. Three groups were designed, including the control group, the model group and the treatment group. Vinorelbine was administered at 0.25 ng per larva by intravenous microinjection to generate the zebrafish model of macrophage deficiency. The treatment group was incubated in BGLS or RGLS supplemented fish water after microinjection. The final concentrations of BGLS (22, 67 and 200 μg/mL) and RGLS (111, 333 and 1000 μg/mL) were set according to the MTC assay in the Albino zebrafish (Supplementary Material). After 48 h incubation, the embryos were stained with 3 mL neutral red (2.5 μg/mL). Subsequently, the zebrafish embryos were immobilized in 3% methylcellulose, and the number of macrophages were counted and recorded with a Nikon dissecting microscope SMZ645. The macrophage formation efficiency was calculated using the following formula:

$\begin{array}{l}\text{macrophage formation efficiency}=\frac{N_{treatment}-N_{model}}{N_{model}}\times 100\%\end{array}$

where N_treatment and N_model were the numbers of macrophages of the larvae in the treatment group and model group, respectively.

Zebrafish Model of Macrophage Phagocytosis

A zebrafish model of PM2.5 phagocytosis was used to assess the phagocytic function of macrophages under the effect of BGLS or RGLS. The model was created by injecting 10 nL active carbon nanoparticles (ACNP, 2.3 mg/mL) to 3-dpf Albino zebrafish. The group design and the tested concentrations were identical to the macrophage deficiency assay. After 24 h incubation, embryos were stained with 3 mL neutral red (2.5 μg/mL) and immobilized in 3% methylcellulose. The number of macrophages that phagocytized ACNP was then counted with a Nikon dissecting microscope SMZ645. The macrophage phagocytosis efficiency was calculated using the following formula:

$\begin{array}{l}\text{macrophage phagocytosis efficiency}\\ =\frac{N_{treatment}-N_{model}}{N_{model}}\times 100\%\end{array}$

where N_treatment and N_model were the numbers of macrophages that phagocytized ACNP in the treatment group and model group, respectively.

Statistical Analysis

The data obtained from zebrafish assays was analyzed with GraphPad prism 7 software (GraphPad Software, USA). Parameter comparisons between groups were made with one-way ANOVA analysis of variance. The result was considered statistically significant when P-value < 0.05.

Results and Discussion

Molecular Networking to Profile Ganoderma lucidum Spores

Representative UPLC–Q–TOF/MS chromatograms of GLS (raw material), sporoderm-broken Ganoderma lucidum spores (BGLS) and sporoderm-removed Ganoderma lucidum spores (RGLS) are shown in Figure 2. Apparently, the chemical profiles of these samples varied considerably. Only a few peaks were detected in the chromatogram of GLS, which suggested that the intact sporoderm acted as a barrier against the release of constituents inside the spores. In comparison, both BGLS and RGLS had much more peaks than GLS. The peak intensities of RGLS were higher than those of BGLS, which could be ascribed to the removal of the sporoderm in RGLS.

FIGURE 2

Figure 2. Representative base peak chromatograms of different Ganoderma lucidum spore samples obtained by LC-QTOF-MS in negative ion mode. BGLS, Ganoderma lucidum spore; BGLS, sporoderm-broken Ganoderma lucidum spore; RGLS, sporoderm-removed Ganoderma lucidum spore.

Given the large quantity of constituents with diverse chemical structures in the spore, we next employed molecular networking for chemical identification, and focused on RGLS, which had the most peaks with higher intensities. The molecular network of RGLS based on MS² spectra similarity was created on GNPS, which contained 501 distinguishable precursor ions, visualized as nodes in the network with 68 clusters (node ≥2), and 186 single nodes. Constituents were then identified and dereplicated through automatic searching in the spectral libraries on GNPS. Besides molecular networking, targeted LC-MS analysis was also conducted to assist the identification based on several strategies developed in our previous studies (Xiao et al., 2014; Li et al., 2015). Briefly, molecular formulae were first generated according to the high-resolution MS data, then the putative identification of the peaks was assigned based on literature and database matching, and was further confirmed via MS² fragmentations. In addition, 13 constituents were unambiguously confirmed by comparisons with chemical standards in terms of retention time and mass spectra (Table 1). Nodes in the network corresponding to precursor ions of the 13 constituents were positioned and used to propagate molecular annotations, which accelerated dereplication of structurally related molecules. By applying these approaches, a total of 20, 96, and 109 constituents were identified or tentatively characterized from GLS, BGLS, and RGLS, respectively, including 99 triterpenoids, one linoleic acid, and 9 potentially new compounds (Table 1). Taking the cluster in Figure 3 as an example, node A showed a quasi-molecular ion [M-H]⁻ at m/z 513.2850, giving the formula C₃₀H₄₂O₇. Fragment ions at m/z 495, 451, and 436 corresponded to successive losses of H₂O, CO₂, and CH₃. Fragment ions at m/z 301 and 285 were characteristic ions formed by the cleavage of the D-ring ([M-H-C₁₁H₁₆O₄]⁻) and the subsequent loss of CH₄. In addition, this node showed identical retention time and similar mass spectra with those of the reference standard ganoderic acid D. Thus, this compound was unambiguously assigned as ganoderic acid D. Node B, which was adjacent to node A with high MS² spectral similarity, showed 2 Da difference in the quasi-molecular ion and many fragment ions, indicating the two compounds shared similar structures. Moreover, node B exhibited identical fragments at m/z 301 and 285 with node A, which represented their A-, B-, and C-ring might be identical. Therefore, this compound was rapidly identified as ganoderenic acid G, with only one double-bond difference in the side chain. Likewise, as the neighbor of node B, node C exhibited a similar fragmentation pattern with node A and B, thus was identified as lucidenic acid C. Therefore, combining molecular networking with targeted MS analysis, large-scale MS dataset can be explored rapidly without any prior knowledge regarding the chemical compositions of the samples, and greatly facilitated the discovery of novel analogs. Detailed MS information is displayed in Table 1.

TABLE 1

Table 1. Characterization of constituents in Ganoderma lucidum spore by LC-QTOF-MS.

FIGURE 3

Figure 3. A representative MS/MS similarity network of the sporoderm-removed Ganoderma lucidum spore. Orange nodes represent reference triterpenoids (i.e., ganoderic acid D, H, and F), while green nodes represent identified triterpenoids.

Zebrafish Assays to Assess Immunomodulatory Activities

The large difference of chemical profiles of BGLS and RGLS indicated their bioactivities might vary. Due to its convenience and optical accessibility, the zebrafish (D. rerio) has been widely adopted as a model for understanding the mechanisms of development and recently, there has been increasing use of the this organism in varied fields including immunology (Novoa and Figueras, 2012). Therefore, as GLS is considered as a potential immunotherapy agent (Cao et al., 2018), the zebrafish models of neutropenia and macrophage deficiency were employed to evaluate the immunomodulatory activities in terms of neutropenia recovery, macrophage formation, and macrophage phagocytosis. Vinorelbine was intravenously injected at 1 and 0.25 ng per larva to generate the neutropenia and macrophage deficiency models, respectively (Figures 4A, 5A). As shown in Figure 4B, the fluorescence intensity of the model group decreased significantly compared with the control, indicating the model was successfully established. After exposing to different concentrations of BGLS and RGLS for 24 h, the number of neutrophils in the larvae recovered with different degrees. BGLS of 22 μg/mL and RGLS of 33 μg/mL significantly improved the neutrophils compared with the model (P < 0.05 or 0.01), while RGLS exhibited more potent effects than BGLS (Figures 4C,D).

FIGURE 4

Figure 4. Assessment of immunoactivity of Ganoderma lucidum spores on zebrafish models of neutropenia. (A) Scheme of zebrafish husbandry and treatment. (B) Representative fluorescent images of Tg (mpx:GFP) transgenic larvae of control and model groups. (C) Neutrophils count in zebrafish of different groups. (D) Neutrophil recovery rate (%) of different groups. BGLS, sporoderm-broken Ganoderma lucidum spores; RGLS, sporoderm-removed Ganoderma lucidum spores. Data are expressed as mean ± SEM, n = 10. *P < 0.05, **P < 0.01 vs. Model.

FIGURE 5

Figure 5. Assessment of immunoactivity of Ganoderma lucidum spores on zebrafish models of macrophage deficiency and macrophage phagocytosis. (A) Overview of zebrafish development and experimental protocol. (B) Represent images of macrophages in healthy fish and zebrafish model of macrophage deficiency. (C) Macrophages count in zebrafish of different groups. (D) The number of macrophages that phagocytize active carbon nanoparticles of different groups. (E) Macrophage formation efficiency (%) of different groups. (F) Macrophage phagocytosis efficiency (%) of different groups. BGLS, sporoderm-broken Ganoderma lucidum spores; RGLS, sporoderm-removed Ganoderma lucidum spores. Data are expressed as mean ± SEM, n = 10. **P < 0.01, ***P < 0.001 vs. Model.

In the macrophage deficiency model created by vinorelbine (Figure 5B), however, only RGLS of higher concentrations (i.e., 333 and 1,000 μg/mL) could significantly promote the formation of macrophages compared with the model (P < 0.01), while BGLS showed very weak effects (Figures 5C,D). Likewise, in the macrophage phagocytosis model, RGLS of 333 and 1,000 μg/mL significantly increased the number of macrophages that phagocytized ACNP, while BGLS of tested concentrations exhibited slight improvement (Figures 5E,F).

Screening Active Compounds by PLSR-Based Activity Ranking

The above results suggested the immunomodulatory activity of GLS was positively correlated with the content of triterpenoids. However, more than 100 triterpenoids were detected in the GLS, and their contents varied considerably, making it difficult to assign the activity to a single molecule. In addition, as most of these triterpenoids were not commercially available, it would be rather labor-intensive and costly to purify and test the compounds individually. To rapidly screen potential active compounds, we employed the partial least squares regression (PLSR) algorithm to find the correlation between the peak area information (X variables) and the activity (Y variable, i.e., neutrophil recovery rate, macrophage formation efficiency or macrophage phagocytosis efficiency). The activity index was proposed to evaluate the contribution of each constituent to the activities, which was calculated using the following formula:

$\begin{array}{l}\text{Y}={\displaystyle \sum _{i=1}^{n=109}}{a}_i{x}_i{C}_i\end{array}$

where Y was the activity of the tested sample; a_i was the activity index of constituent i; x_i was the peak area of constituent i in the tested sample; C_i was the relative concentration of the tested sample, which was defined as c_i/c_max, where c_i was the concentration of the tested sample and c_max was the max concentration in each group.

The PLSR model was created according to PLS-W2A (Wegelin, 2000). The dataset obtained was transferred to the center of the multidimensional coordinate system to identify potential active compounds that had high activity index in all the three bioassays. All programs were operated by Spyder with Python 3.6 (Anaconda, Inc., Austin, Texas, USA) software.

As shown in Figure 6, 11 compounds, namely 20-hydroxylganoderic acid G, elfvingic acid A, elfvingic C, lucidenic acid I, ganohainanic acid C, ganoderic acid β, methyl lucidenate E2, dehydrolucidenic acid N, applanoxidic acid G, ganolucidic acid D and ganoderenic acid H were identified as potential immunoactive compounds that ranked high in neutropenia recovery, macrophage formation, and macrophage phagocytosis. However, these candidates needed meticulous validation in vivo and in vitro to confirm their pharmacological effects.

FIGURE 6

Figure 6. Identification of immunoactive compounds of Ganoderma lucidum spore by activity index. anf, amf, apf are the activity indices of the compound, representing its effects on neutrophil recovery, macrophage formation and macrophages phagocytosis, respectively.

Discussion

It has been reported that G. lucidum contains over 400 bioactive compounds, including triterpenoids, polysaccharides, nucleotides, sterols, steroids, fatty acids and proteins/peptides, which have various medicinal effects. In the present study, we identified chemical constitutes of GLS by LC-MS and the main constitutes of GLS were identified as triterpenoids. We further employed zebrafish as the animal model to investigate the immunomodulatory activities of GLS, due to its advantages in immune-related research, which included: (1) having complete (innate and adaptative) immune systems, and possessing neutrophils and macrophages in adults and larvae; (2) having a relatively rapid life cycle and was easy to maintain; (3) quickly producing large numbers of offspring that could be assayed in multi-well plates and treated with different chemicals. Finally, activity index was proposed to evaluate the contribution of each constituent to the activity to screen potential immunoactive ones.

In a previous study, Ahmadi and Riazipour revealed that G. lucidum could improve macrophage function through cytokine and NO release (Ahmadi and Riazipour, 2007). Besides, several pharmacological studies have shown that G. lucidum can play an antitumor role through the regulation of the immune system (Boh et al., 2007). A finding suggested that the triterpenes of G. lucidum inhibited anti-lung cancer in vitro and in vivo via enhancement of immunomodulation and induction of cell apoptosis (Liang et al., 2013). Therefore, it can be predicted that GLS may affect the formation and phagocytic function of macrophages to exert the immunomodulatory effect.

Among the compounds identified in G. lucidum, triterpenoids are most widely investigated. Besides the anticancer activities which have been reported by many studies, triterpenoids also showed good immunomodulatory effects. It has been reported that the triterpenoids extract from G. lucidum have anti-inflammatory and anti-proliferative effects, which are mediated through the inhibition of NF-κB and AP-1 signaling pathways in macrophages (Dudhgaonkar et al., 2009). As a major triterpene of GL and GLS, previous studies suggested that Ganoderic acid A played a role in immunomodulation through inhibiting the release of proinflammatory mediators like IL-1β, IL-6 and TNF-α (Chi et al., 2018). A study also suggested that Ganoderic acid A-treatment not only enhanced cell-mediated immune responses, but also potentiated antitumor immune responses by activating IFN-γ producing CD8+ T cells (Radwan et al., 2015). Liu et al. found that ganoderic acid C₁ had an effect on immune system and significantly suppressed murine macrophage TNF-α production, which was associated with suppression of NF-κB (Liu et al., 2015). These studies gave possible direction of our further mechanism study.

Conclusion

This work demonstrates the feasibility of mass spectrometry molecular networking coupled with zebrafish-based bioassays and chemometrics for active constituent identification of complex herbal medicine. An efficient MS data processing method integrating molecular networking and targeted LC-MS analysis was established to comprehensively profile GLS, leading to the characterization of 109 constituents. Immunomodulatory activities of different GLS samples were evaluated by zebrafish models of neutrophil or macrophage deficiency, in which RGLS showed better therapeutic effects than BGLS. Moreover, a three-dimensional activity index approach based on PLSR was developed to identify active constituents in GLS that were effective on all the three bioassays, which included 20-hydroxylganoderic acid G, elfvingic acid A, elfvingic C, and etc. Future works are needed to validate the pharmacological effects of these compounds using purified substances on various in vivo assays. Also possible for further work is employing ions in the molecular network rather than identified peaks for correlation analysis, to explore the potential relationship between bioactivities and certain ions/structures.

Data Availability Statement

The datasets generated for this study are available on request to the corresponding author.

Ethics Statement

The animal study was reviewed and approved by the International Association for Assessment and Accreditation of Laboratory Animal Care (AAALAC).

Author Contributions

ZL, JX, and YW designed the research. ZL, YS, and XZ performed the experimental work and analyzed the data. ZL, JX, and HW prepared the GLP extracts. ZL, YS, LZ, and YW participated in the preparation of manuscript. All authors proofread the paper and approved the final version of the manuscript.

Funding

This work was supported by grants from the National Natural Science Foundation of China (Grant Nos. 81803714, 81822047), the National Major Science and Technology Project of China (Grant No. 2017ZX09301012), the Science Technology Department of Zhejiang Province (2016C33060), and the Key R&D Foundation of Zhejiang Province (2019C02100).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2020.00287/full#supplementary-material

References

Ahmad, M. F. (2018). Ganoderma lucidum: persuasive biologically active constituents and their health endorsement. Biomed. Pharmacother. 107, 507–519. doi: 10.1016/j.biopha.2018.08.036

PubMed Abstract | CrossRef Full Text | Google Scholar

Ahmadi, K., and Riazipour, M. (2007). Effect of Ganoderma lucidum on cytokine release by peritoneal macrophages. Iran. J. Immunol. 4, 220–226.

PubMed Abstract | Google Scholar

Allard, P.-M., Peresse, T., Bisson, J., Gindro, K., Marcourt, L., Van Cuong, P., et al. (2016). Integration of molecular networking and in-silico MS/MS fragmentation for natural products dereplication. Anal. Chem. 88, 3317–3323. doi: 10.1021/acs.analchem.5b04804

PubMed Abstract | CrossRef Full Text | Google Scholar

Bishop, K. S., Kao, C. H. J., Xu, Y., Glucina, M. P., Paterson, R. R. M., and Ferguson, L. R. (2015). From 2000 years of Ganoderma lucidum to recent developments in nutraceuticals. Phytochemistry 114, 56–65. doi: 10.1016/j.phytochem.2015.02.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Boh, B., Berovic, M., Zhang, J., and Lin, Z. B. (2007). Ganoderma lucidum and its pharmaceutically active compounds. Biotechnol. Annu. Rev. 13, 265–301. doi: 10.1016/S1387-2656(07)13010-6

PubMed Abstract | CrossRef Full Text | Google Scholar

Cao, Y., Xu, X., Liu, S., Huang, L., and Gu, J. (2018). Ganoderma: a cancer immunotherapy review. Front. Pharmacol. 9:1217. doi: 10.3389/fphar.2018.01217

PubMed Abstract | CrossRef Full Text | Google Scholar

Chang, J. B., Lane, M. E., Yang, M., and Heinrich, M. (2018). Disentangling the complexity of a hexa-herbal Chinese medicine used for inflammatory skin conditions-predicting the active components by combining LC-MS-based metabolite profiles and in vitro pharmacology. Front. Pharmacol 9:1091. doi: 10.3389/fphar.2018.01091

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, G., Huang, B. X., and Guo, M. (2018). Current advances in screening for bioactive components from medicinal plants by affinity ultrafiltration mass spectrometry. Phytochem. Analysis 29, 375–386. doi: 10.1002/pca.2769

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, X., Liu, X., Sheng, D., Huang, D., Li, W., and Wang, X. (2012). Distinction of broken cellular wall Ganoderma lucidum spores and G. lucidum spores using FTIR microspectroscopy. Spectrochim. Acta A 97, 667–672. doi: 10.1016/j.saa.2012.07.046

PubMed Abstract | CrossRef Full Text | Google Scholar

Chi, B., Wang, S., Bi, S., Qin, W., Wu, D., Luo, Z., et al. (2018). Effects of ganoderic acid A on lipopolysaccharide-induced proinflammatory cytokine release from primary mouse microglia cultures. Exp. Ther. Med. 15, 847–853. doi: 10.3892/etm.2017.5472

PubMed Abstract | CrossRef Full Text | Google Scholar

Dudhgaonkar, S., Thyagarajan, A., and Sliva, D. (2009). Suppression of the inflammatory response by triterpenes isolated from the mushroom Ganoderma lucidum. Int. Immunopharmacol. 9, 1272–1280. doi: 10.1016/j.intimp.2009.07.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Fu, H., Jia, T., Zhang, Y., Fumimaru, M., Eri, O., Li, H., et al. (2012). The role of three kinds of ganoderma lucidum powder on regulation of immune functions of immunosuppressive mice model. Chinese J. Immunol. 28, 712–716. doi: 10.3969/j.issn.1000-484X.2012.08.008

CrossRef Full Text

Gao, G., Bao, H., and Bau, T. (2013). Anti-tumor activity of the nano-scale fruiting body powder and the wall-broken spore powder of Ganoderma lucidum in vitro. Mygosystema 32, 114–127. doi: 10.13346/j.mycosystema.2013.01.008

CrossRef Full Text

Ge, Y. W., Zhu, S., Yoshimatsu, K., and Komatsu, K. (2017). MS/MS similarity networking accelerated target profiling of triterpene saponins in Eleutherococcus senticosus leaves. Food Chem. 227, 444–452. doi: 10.1016/j.foodchem.2017.01.119

PubMed Abstract | CrossRef Full Text | Google Scholar

Guthals, A., Watrous, J. D., Dorrestein, P. C., and Bandeira, N. (2012). The spectral networks paradigm in high throughput mass spectrometry. Mol. Biosyst. 8, 2535–2544. doi: 10.1039/c2mb25085c

PubMed Abstract | CrossRef Full Text | Google Scholar

Hofstadler, S. A., and Sannes-Lowery, K. A. (2006). Applications of ESI-MS in drug discovery: interrogation of noncovalent complexes. Nat. Rev. Drug Disco. 5, 585–595. doi: 10.1038/nrd2083

PubMed Abstract | CrossRef Full Text | Google Scholar

Hsu, K.-D., and Cheng, K.-C. (2018). From nutraceutical to clinical trial: frontiers in Ganoderma development. Appl. Microbiol. Biot. 102, 9037–9051. doi: 10.1007/s00253-018-9326-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Kessner, D., Chambers, M., Burke, R., Agusand, D., and Mallick, P. (2008). ProteoWizard: open source software for rapid proteomics tools development. Bioinformatics 24, 2534–2536. doi: 10.1093/bioinformatics/btn323

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Z., Guo, X., Cao, Z., Liu, X., Liao, X., Huang, C., et al. (2018). New MS network analysis pattern for the rapid identification of constituents from traditional Chinese medicine prescription Lishukang capsules in vitro and in vivo based on UHPLC/Q-TOF-MS. Talanta 189, 606–621. doi: 10.1016/j.talanta.2018.07.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Z., Liu, T., Liao, J., Ai, N., Fan, X., and Cheng, Y. (2017). Deciphering chemical interactions between glycyrrhizae radix and coptidis rhizoma by liquid chromatography with transformed multiple reaction monitoring mass spectrometry. J. Sep. Sci. 40, 1254–1265. doi: 10.1002/jssc.201601054

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Z., Wang, Y., and Cheng, Y. (2019). Mass spectrometry-sensitive probes coupled with direct analysis in real time for simultaneous sensing of chemical and biological properties of botanical drugs. Anal. Chem. 91, 9001–9009. doi: 10.1021/acs.analchem.9b01251

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, Z., Xiao, S., Ai, N., Luo, K., Fan, X., and Cheng, Y. (2015). Derivative multiple reaction monitoring and single herb calibration approach for multiple components quantification of traditional Chinese medicine analogous formulae. J. Chromatogr. A 1376, 126–142. doi: 10.1016/j.chroma.2014.12.024

PubMed Abstract | CrossRef Full Text | Google Scholar

Liang, F., Ling, Y., Meng, D., Yan, C., Ming-Hua, Z., Jun-Fei, G., et al. (2013). Anti-lung cancer activity through enhancement of immunomodulation and induction of cell apoptosis of total triterpenes extracted from Ganoderma luncidum (Leyss. ex Fr.) Karst. Molecules 18, 9966–9981. doi: 10.3390/molecules18089966

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, Z. B., and Wang, P. Y. (2006). The pharmacological study of Ganoderma spores and their active components. J. Peking Univ. Health Sci. 38, 541–547.

PubMed Abstract | Google Scholar

Liu, C., Yang, N., Song, Y., Wang, L., Zi, J., Zhang, S., et al. (2015). Ganoderic acid C1 isolated from the anti-asthma formula, ASHMI™ suppresses TNF-α production by mouse macrophages and peripheral blood mononuclear cells from asthma patients. Int. Immunopharmacol. 27, 224–231. doi: 10.1016/j.intimp.2015.05.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, H., Chen, X., Zhao, X., Zhao, B., Qian, K., Shi, Y., et al. (2018). Screening and identification of cardioprotective compounds from wenxin keli by activity index approach and in vivo zebrafish model. Front. Pharmacol. 9:1288. doi: 10.3389/fphar.2018.01288

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, X., Wang, J. H., and Yuan, J. P. (2005). Pharmacological and anti-tumor activities of Ganoderma spores processed by top-down approaches. J. Nanosci. Nanotec. 5, 2001–2013. doi: 10.1166/jnn.2005.448

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, Y., Liu, Y., Qiu, F., and Di, X. (2011). Sensitive and selective liquid chromatography-tandem mass spectrometry method for the determination of five ganoderic acids in Ganoderma lucidum and its related species. J. Pharmaceut. Biomed. 54, 717–721. doi: 10.1016/j.jpba.2010.11.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Nayak, R. N., Dixitraj, P. T., Nayak, A., and Bhat, K. (2015). Evaluation of anti-microbial activity of spore powder of Ganoderma lucidum on clinical isolates of Prevotella intermedia: a pilot study. Contemp. Clin. Dent. 6, S248–252. doi: 10.4103/0976-237X.166834

PubMed Abstract | CrossRef Full Text | Google Scholar

Newman, D. J., and Cragg, G. M. (2016). Natural products as sources of new drugs from 1981 to 2014. J. Nat. Prod. 79, 629–661. doi: 10.1021/acs.jnatprod.5b01055

PubMed Abstract | CrossRef Full Text | Google Scholar

Novoa, B., and Figueras, A. (2012). Zebrafish: model for the study of inflammation and the innate immune response to infectious diseases. Adv. Exp. Med. Biol. 946, 253–275. doi: 10.1007/978-1-4614-0106-3_15

PubMed Abstract | CrossRef Full Text | Google Scholar

Pacholarz, K. J., Garlish, R. A., Taylor, R. J., and Barran, P. E. (2012). Mass spectrometry based tools to investigate protein-ligand interactions for drug discovery. Chem. Soc. Rev. 41, 4335–4355. doi: 10.1039/c2cs35035a

PubMed Abstract | CrossRef Full Text | Google Scholar

Qiang, L., Ting-Hao, K., Chongde, S., Kunsong, C., Cheng-Chih, H., and Xian, L. (2019). Comprehensive structural characterization of phenolics in litchi pulp using tandem mass spectral molecular networking. Food Chem. 282, 9–17. doi: 10.1016/j.foodchem.2019.01.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Quinn, R. A., Nothias, L. F., Vining, O., Meehan, M., Esquenazi, E., and Dorrestein, P. C. (2017). Molecular networking as a drug discovery, drug metabolism, and precision medicine strategy. Trends Pharmacol. Sci. 38, 143–154. doi: 10.1016/j.tips.2016.10.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Radwan, F. F. Y., Hossain, A., God, J. M., Leaphart, N., Elvington, M., Nagarkatti, M., et al. (2015). Reduction of myeloid-derived suppressor cells and lymphoma growth by a natural triterpenoid. J. Cell. Biochem. 116, 102–114. doi: 10.1002/jcb.24946

PubMed Abstract | CrossRef Full Text | Google Scholar

Russell, R., and Paterson, M. (2006). Ganoderma – a therapeutic fungal biofactory. Phytochemistry 67, 1985–2001. doi: 10.1016/j.phytochem.2006.07.004

PubMed Abstract | CrossRef Full Text | Google Scholar

Soccol, C. R., Bissoqui, L. Y., Rodrigues, C., Rubel, R., Sella, S., Leifa, F., et al. (2016). Pharmacological properties of biocompounds from spores of the lingzhi or reishi medicinal mushroom Ganoderma lucidum (Agaricomycetes): a review. Int. J. Med. Mushrooms 18, 757–767. doi: 10.1615/IntJMedMushrooms.v18.i9.10

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, G. H., Wang, L. H., Wang, C., and Qin, L. H. (2018). Spore powder of Ganoderma lucidum for the treatment of alzheimer disease: a pilot study. Medicine 97:e0636. doi: 10.1097/MD.0000000000010636

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, M., Carver, J. J., Phelan, V. V., Sanchez, L. M., Garg, N., Peng, Y., et al. (2016a). Sharing and community curation of mass spectrometry data with global natural products social molecular networking. Nat. Biotechnol. 34, 828–837. doi: 10.1038/nbt.3597

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, X., Zhang, A., Sun, H., Han, Y., and Yan, G. (2016b). Discovery and development of innovative drug from traditional medicine by integrated chinmedomics strategies in the post-genomic era. Trac-Trends Anal. Chem. 76, 86–94. doi: 10.1016/j.trac.2015.11.010

CrossRef Full Text | Google Scholar

Wang, X.-J., Ren, J.-L., Zhang, A.-H., Sun, H., Yan, G.-L., Han, Y., et al. (2019). Novel applications of mass spectrometry-based metabolomics in herbal medicines and its active ingredients: current evidence. Mass Spectrom. Rev. 38, 380–402. doi: 10.1002/mas.21589

PubMed Abstract | CrossRef Full Text | Google Scholar

Watrous, J., Roach, P., Alexandrov, T., Heath, B. S., Yang, J. Y., Kersten, R. D., et al. (2012). Mass spectral molecular networking of living microbial colonies. Proc. Natl. Acad. Sci. U.S.A. 109, E1743–E1752. doi: 10.1073/pnas.1203689109

PubMed Abstract | CrossRef Full Text | Google Scholar

Wegelin, J. A. (2000). A survey of Partial Least Squares (PLS) Methods, With Emphasis on the Two-Block Case. Seattle, WA: University of Washington, Tech. Rep.

Google Scholar

Xiao, S., Luo, K. D., Wen, X. X., Fan, X. H., and Cheng, Y. Y. (2014). A pre-classification strategy for identification of compounds in traditional Chinese medicine analogous formulas by high-performance liquid chromatography-mass spectrometry. J. Pharmaceut. Biomed. 92, 82–89. doi: 10.1016/j.jpba.2013.12.042

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, J., Liu, Z., Wang, Y., Zhu, H., Zheng, H., Li, M., et al. (2014). Determination of inorganic elements in Ganoderma lucidum spore powder with different treatments of wall-broken methods by microwave digestion-ICP-MS. Chin. J. Mod. Appl. Pharm. 31, 813–817.

Yan, Z., Xia, B., Qiu, M. H., Sheng, D. L., and Xu, H. X. (2013). Fast analysis of triterpenoids in Ganoderma lucidum spores by ultra-performance liquid chromatography coupled with triple quadrupole mass spectrometry. Biomed. Chromatogr. 27, 1560–1567. doi: 10.1002/bmc.2960

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, D., Limin, H., Yuqi, L., Liming, Z., Jike, L., Qingwu, H., et al. (2017). Effect of different cell wall disruption methods on in vitro antioxidant activity of crude polysaccharides Ganoderma lucidum spores. Food Sci. Chin. 38, 101–106. doi: 10.7506/spkx1002-6630-201717017

CrossRef Full Text

Yang, J. Y., Sanchez, L. M., Rath, C. M., Liu, X., Boudreau, P. D., Bruns, N., et al. (2013). Molecular networking as a dereplication strategy. J. Nat. Prod. 76, 1686–1699. doi: 10.1021/np400413s

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, Z., Zhang, Y., Sun, L., Wang, Y., Gao, X., and Cheng, Y. (2012). An ultrafiltration high-performance liquid chromatography coupled with diode array detector and mass spectrometry approach for screening and characterising tyrosinase inhibitors from mulberry leaves. Analy. Chim. Acta 719, 87–95. doi: 10.1016/j.aca.2012.01.018

PubMed Abstract | CrossRef Full Text | Google Scholar

Yuan, Y., Wang, Y., Sun, G., Wang, Y., Cao, L., Shen, Y., et al. (2018). Archaeological evidence suggests earlier use of Ganoderma in Neolithic China. Chinese Sci. Bull. 63, 1180–1188. doi: 10.1360/N972018-00188

CrossRef Full Text | Google Scholar

Zhao, H., Zhang, Q., Zhao, L., Huang, X., Wang, J., and Kang, X. (2012). Spore powder of Ganoderma lucidum improves cancer-related fatigue in breast cancer patients undergoing endocrine therapy: a pilot clinical trial. Evid.-Based Complement. Alternat. Med. 2012:809614. doi: 10.1155/2012/809614

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, L.-W., Cao, Y., Wu, S.-H., Vlasak, J., Li, D.-W., Li, M.-J., et al. (2015). Global diversity of the Ganoderma lucidum complex (Ganodermataceae, Polyporales) inferred from morphology and multilocus phylogeny. Phytochemistry 114, 7–15. doi: 10.1016/j.phytochem.2014.09.023

PubMed Abstract | CrossRef Full Text | Google Scholar