- Open Access
Exploration of anti-inflammatory mechanism of forsythiaside A and forsythiaside B in CuSO4-induced inflammation in zebrafish by metabolomic and proteomic analyses
Journal of Neuroinflammation volume 17, Article number: 173 (2020)
Inflammation is a general pathological phenomenon during severe disturbances to the homeostasis. Forsythiaside A (FA) and forsythiaside B (FB), isolated from the dried fruit of Forsythia suspensa (Thunb.) Vahl, are phenylethanoid compounds that show a significant anti-inflammatory effect. However, the properties and therapeutic mechanisms of this effect have not yet been systematically elucidated.
In this study, the anti-inflammatory effects of FA and FB were investigated in CuSO4-induced inflammation in zebrafish larvae. Intracellular generation of reactive oxygen species (ROS) and nitric oxide (NO) was investigated using fluorescence probes. Metabolomic and proteomic analyses using liquid chromatography-mass spectrometry were carried out to identify the expressions of metabolites and proteins associated with the anti-inflammatory mechanism of FA and FB. Quantitative polymerase chain reaction (PCR) was performed to detect the progressive changes in gene expression.
FA and FB inhibited neutrophils migration to the damaged neuromasts and remarkably reduced CuSO4-induced ROS and NO generation in zebrafish larvae. Metabolomic analysis pointed to the involvement of nicotinate and nicotinamide metabolism, energy metabolism, pyrimidine metabolism, and purine metabolism. Proteomic analysis identified 146 differentially expressed proteins between the control and model groups. These included collagen [collagen type II alpha 1b precursor (col2a1b), collagen alpha-2(IX) chain precursor (col9a2), collagen type IX alpha I precursor (col9a1b)], nucleoside diphosphate kinase 3 isoform X1 (Nme3), WD repeat-containing protein 3 (Wdr3), and 28S ribosomal protein S7 mitochondrial precursor (Mrps7). FA and FB were shown to reverse the abnormal expressions of potential metabolite and protein biomarkers and alleviate CuSO4-induced damage to the neuromasts in the zebrafish lateral line.
Our results indicate that FA and FB possess remarkable anti-inflammatory properties, protecting against CuSO4-induced neuromasts damage in zebrafish larvae. The results also suggest a multi-component and multi-regulatory therapeutic mechanism for FA and FB.
The inflammatory response is a key component in normal homeostasis that protects the body from irritation and restores damaged tissue structure and function. Generally, inflammatory reactions are beneficial to the body. However, excessive and uncontrolled inflammation can cause chronic diseases such as cancer and neurodegenerative diseases, including Parkinson's disease and Alzheimer's disease [1, 2]. During this process, over-activation of the macrophages and neutrophils can induce the secretion of interleukin-6 (IL-6), tumor necrosis factor-α (TNF-α), and interleukin-1β (IL-1β), which are significant mediators of the inflammatory response .
Forsythiae Fructus is the dried fruit of Forsythia suspensa (Thunb.) Vahl. As an antipyretic and anti-inflammatory agent used in traditional Chinese medicine, Forsythiae Fructus has been used to treat various infectious diseases, such as acute nephritis and ulcers . Forsythiaside A (FA, Fig. 1a) and forsythiaside B (FB, Fig. 1b) are the major bioactive components extracted from Forsythiae Fructus. They were reported to possess anti-inflammatory and anti-bacterial properties [5, 6]. However, few studies have been conducted to systematically investigate their anti-inflammatory properties and therapeutic mechanisms.
Zebrafish (Danio rerio), with morphological and physiological functions similar to humans, are being widely used in pharmacology research. It has unique advantages such as high fecundity, easy breeding, small size, availability of various transgenic lines, simple in vivo operation, and the possibility to perform real-time imaging [7,8,9]. The transparency of zebrafish larvae and the availability of transgenic zebrafish lines make it possible to monitor the inflammatory processes and observe cells behavior in vivo . Studies have also shown that the zebrafish immune system shares significant similarities with humans, and almost all the human immune system cells have counterparts in zebrafish . The neuromasts in the zebrafish lateral line system are composed of mechanosensory hair cells and can be damaged by physical or chemical stimulation . CuSO4 is a metal chemical that induces neuromasts damage in zebrafish, followed by an infiltration of inflammatory cells, which in turn results in a progressive disruption of the neuromasts structure . Unfortunately, the specific role of FA and FB in CuSO4-induced inflammation in zebrafish and the associated molecular mechanisms remain to be further explored.
The biological interpretation of common pharmacological research might be a great challenge due to the complicated biochemical regulatory mechanism at multiple levels. An integrated analysis of omics data is promising and can provide an improved comprehension of biological mechanisms by identifying the potential biomarkers and interpreting their intrinsic relations. The combination of metabolomics and proteomics is a powerful tool that frequently applied in complex mechanism exploration. Therefore, we investigated the anti-inflammatory mechanisms of FA and FB through an untargeted metabolic and proteomic profiling of zebrafish, using liquid chromatography-mass spectrometry. Based on the results, we propose a potential mechanism to explain the inflammatory behavior in zebrafish neuromasts during CuSO4 exposure, and systematically elucidate the anti-inflammatory effects and mechanisms of FA and FB. We thus expand the understanding of FA and FB in the treatment of inflammation and their potential clinical applications.
FA (MUST-18010303) and FB (MUST-18081202) were purchased from Chengdu MUST Bio-technology Co., Ltd. (Chengdu, China). Tricaine was obtained from Sigma-Aldrich (Shanghai, China). Acetonitrile, methanol, and formic acid of HPLC grade were collected from Merck Chemicals (Shanghai, China), Wokai Chemical Technology Co., Ltd. (Shanghai, China) and TCI Chemical Industry Development Co., Ltd. (Shanghai, China), respectively. DAF-FMDA and DCF-DA were purchased from Yeasen Bio-technology Co., Ltd. (Shanghai, China). TRIzol reagents were purchased from Ambion Life Technologies (Carlsbad, CA, USA). 5 × All-In-one MasterMix and Eva Green 2 × RT-qPCR MasterMix-Low RoX were purchased from Applied Biological Materials Inc. (Richomnd, BC, Canada). PCR primer sequences were synthesized in TSINGKE Biological Technology (Chengdu, China). Other chemicals and reagents used in this study were obtained from Kelong Chemical Reagent Factory (Chengdu, China).
Wild-type AB strain and neutrophil-specific transgenic (mpx: EGFP) zebrafish, in which neutrophils were labelled with an enhanced green fluorescent protein, were obtained from China Zebrafish Resource Center (Wuhan, China), and raised at 28.5 ± 1.0 °C on a 14-h light/10-h dark cycle. Zebrafish embryos were collected by natural spawning and maintained in the embryonic medium (5.00 mmol L−1 NaCl, 0.44 mmol L−1 CaCl2, 0.33 mmol L−1 MgSO4, 0.17 mmol L−1 KCl, and 0.1% methylene blue, equilibrated to pH 7.0) till drug administration [14, 15].
Anti-inflammatory effects of FA and FB against CuSO4-induced inflammation in zebrafish
In vivo neutrophils recruitment assay
All experiments were carried out on 3 days post-fertilization (dpf) zebrafish larvae in the present study. Transgenic (mpx: EGFP) zebrafish larvae at 3 dpf were randomly transferred to a 12-well plate (15 larvae per well) and assigned to 8 groups for neutrophils recruitment investigation. In the control group, larvae were incubated in the embryonic medium. In the model group, larvae were exposed to 10 μmol L−1 CuSO4 for 40 min. In FA group, larvae were exposed to 1 h FA pre-treatment at 120, 60, 30 μmol L−1, followed by a mixture incubation of 10 μmol L−1 CuSO4 and FA (120, 60, 30 μmol L−1) for 40 min. In FB group, larvae were exposed to 1 h FB pre-treatment at 150, 75, 37.5 μmol L−1, followed by a mixture incubation of 10 μmol L−1 CuSO4 and FB (150, 75, 37.5 μmol L−1) for 40 min. All drugs were diluted in the embryonic medium. After the treatments, zebrafish larvae were washed with fresh medium and subsequently anesthetized with tricaine. Then, zebrafish larvae were photographed under a Leica M165Fic fluorescence microscope (Leica Microsystems, Germany). Finally, Image Pro Plus 6.0 software (Media Cybernetics, USA) was applied to quantify the recruitment of neutrophils to the zebrafish neuromasts.
Intracellular production of ROS and NO in zebrafish
2′,7′-dichlorodihydrofluorescein diacetate (DCF-DA) and diaminofluorophore 4-amino-5-methylamino-2′,7′-difluorofluorescein diacetate (DAF-FMDA) were used as fluorescence probes, respectively, to investigate intracellular reactive oxygen species (ROS) and nitric oxide (NO) accumulation in CuSO4-induced inflammation in zebrafish . Wild-type AB strain zebrafish larvae were randomly transferred to a twelve-well plate (15 larvae per well), and assigned to eight groups following a previously-described method for modeling and drug administration. After the treatments, the zebrafish were moved to another twelve-well plate (15 larvae per well) and treated with the DCF-DA (0.05 μmol L−1) or DAF-FMDA (5 μmol L−1) solution. After 1 h incubation in the dark, the zebrafish larvae were washed with fresh medium and subsequently anesthetized with tricaine. The zebrafish larvae were photographed under a Leica M165Fic fluorescence microscope. Image Pro Plus 6.0 software was used to analyze the fluorescence intensity in individual zebrafish larvae to quantify the accumulation of ROS and NO.
Statistical analysis of neutrophils migration and ROS, NO accumulation
Zebrafish larvae were randomly assigned to each treatment group, and pharmacodynamic experiments were repeated three times. Statistical analyses were conducted by SPSS 25.0 (SPSS Inc., Chicago, IL, USA), and graphs were generated by GraphPad Prism 6.0 (GraphPad, San Diego, CA, USA). All statistics regarding neutrophils migration, ROS, and NO accumulation were evaluated by independent sample t test or Mann-Whitney U test (two groups) and one-way ANOVA or Kruskal-Wallis test followed by pairwise comparisons (three or more groups) depending on whether the data were normally distributed. Data were shown as mean ± SD and p < 0.05 were assumed for statistical significance indication.
Sample processing and metabolites detection
Zebrafish larvae from the control, model, FA (120 μmol L−1), and FB (150 μmol L−1) groups were transferred to 2-mL centrifuge tubes, respectively. We then added 1 mL 80% methanol and steel balls to the tubes. Samples were ground for 1 min in a Tissue Grinding Device (SCIENTZ-48, Xinzhi Biotechnology Co., Ltd.) at 70 Hz. Samples were then placed in an Ultrasonic Machine (KW-100TDV, Shumei Ultrasonic Instrument Co., Ltd.) at room temperature for 30 min and subsequently incubated on ice for 30 min. After being centrifuged for 10 min (14,000 rpm, 4 °C), the supernatant was transferred to a new tube. Samples were dried by vacuum concentration and then dissolved with 400 μL 2-chlorobenzalanine methanol aqueous solution (1:1, 4 °C). Samples were filtered through a 0.22-μm microfilter membrane and then prepared for liquid chromatography-mass spectrometry (LC-MS) detection. Quality control samples were obtained by mixing a 20-μL aliquot from each sample together.
Chromatographic separation was performed on a High-performance liquid chromatography (HPLC, Ultimate 3000 system, Thermo), equipped with a Waters column (1.8 μm, 150 × 2.1 mm) at 40 °C. The autosampler temperature was set at 8 °C. Gradient elution of analytes was carried out at a flow rate of 0.25 mL min−1 with 5 mmol L−1 ammonium formate in water (A) and acetonitrile (B) in the negative model or with 0.1% formic acid in water (C) and 0.1% formic acid in acetonitrile (D) in the positive model. Two microliter samples were injected after equilibration. A linear gradient of solvent B (negative model) or solvent D (positive model) was programmed as follows: 0–1 min, 2% B or 0.1% D; 1–9 min, 2–50% B/D; 9–12 min, 50–98% B/D; 12–13.5 min, 98% B/D; 13.5–14 min, 98–2% B/D; 14–20 min, 2% D or 14–17 min, 2% B.
The electrospray ionization mass detection (ESI-MSn) in positive and negative ion modes was performed on a Mass spectrometer (Q Exactive Focus, Thermo) with a spray voltage of 3.8 kV and − 2.5 kV, respectively. Sheath gas and auxiliary gas were set at 45 and 15 arbitrary units, respectively. The capillary temperature was set at 325 °C . The Orbitrap Analyzer (Orbitrap Fusion Lumos mass spectrometer, CA, USA) performed a full scan with a mass to charge ratio (m/z) range of 81–1000 and a mass resolution of 70,000. Data-dependent experiments were performed on mass spectrometry/mass spectrometry (MS/MS) with higher-energy collisional dissociation (HCD) scanning mode. The standard collision energy was 30 eV. Dynamic exclusion was implemented to remove some unnecessary information from the MS/MS spectrum.
Data processing and multivariate analyses of metabolites
Data from the ultra-performance liquid chromatography-mass spectrometry/mass spectrometry (UPLC-MS/MS) were displayed in mzXML format using the ProteoWizard software (v3.0.8789). Data pre-processing was performed by the XCMS package of R language (v3.3.2) using autoscaling, mean-centering, and scaling to unit variance. Parameters were set as follows: bw = 2, ppm = 15, mzwid = 0.015, peak width = c (5, 30), mzdiff = 0.01, method = cent wave. The area of metabolite peaks in each sample was normalized by the summation method applied in Metaboanalyst (http://www.metaboanalyst.ca/). Multivariate analyses of metabolites, including standard peak areas and retention time, were performed using the R language to further locate the specific metabolites in the dataset. A biomarker was generated and queried for accurate molecular weights, with mass errors of less than 30 part per million (ppm). The fracture patterns of potential biomarkers were analyzed, and their identification was performed using the Human Metabolome Database (HMDB) (http://www.hmdb.ca) and METLIN (http://metlin.scripps.edu/). The Kyoto Encyclopedia of Genes and Genomes Database (KEGG) (http://www.kegg.jp/) was used to plot a metabolomic pathway network diagram. Using Metaboanalyst, we performed a pathway enrichment analysis to screen metabolomic pathways of the identified metabolites in the current study.
Statistical analysis of metabolomics
In metabolomic analysis, 10 biological replications were made in each group in the present study. Orthogonal projections to latent structures discriminant analysis (OPLS-DA) was used to filter metabolites, and variable importance in the projection (VIP) values were used to promote group discrimination. The S-Plot (R language: v3.3.2) was applied to detect the differential metabolites with values of VIP ≥ 1. One-way ANOVA of SPSS 25.0 was conducted to determine whether the potential metabolites were significantly different (p ≤ 0.05) among the control, model, FA, and FB groups. Finally, the metabolites with VIP ≥ 1 and p ≤ 0.05 were considered as potential metabolite biomarkers [17, 18].
Protein extraction and digestion
Proteins were extracted from zebrafish larvae of the control, model, FA (120 μmol L−1), and FB (150 μmol L−1) groups, as previously described . Briefly, 500 μL lysis buffer (2% sodium deoxycholate, 50 mM ammonium bicarbonate, 75 mM sodium chloride) was added to the samples, which were then ground on ice in an Ultrasonic Crushing Machine (Scientz-JY92, Ningbo Xinzhi Biotechnology Co., Ltd.) for 10 min, with cycles of 2 s on and 4 s off, at 15% power. The ground material was centrifuged for 10 min (10,000×g, 4 °C). We then added 10 mM dithiothreitol to the supernatant to precipitate the proteins at − 20 °C . The precipitation procedure was repeated with acetone until the supernatant became colorless, as previously described . Then, the supernatant was incubated at 56 °C for 1 h and subsequently alkylated with 55 mmol L−1 iodoacetamide for 45 min in the dark at room temperature. Proteins were then resuspended in the lysis buffer, and the concentration of protein was detected using the BCA assay. The protein solution (100 μg) was digested by Trypsin Gold (40:1, protein:trypsin) at 37 °C overnight. Peptides desalination was carried out using a Strata X C18 column, and the samples were then vacuum-dried following the manufacturer's instructions.
TMT labeling and fractionation
The peptides were labeled using Tandem Mass Tags (TMT) six plex Isobaric Label Reagent Set (Thermo Scientific, 90061) following the manufacturer’s instructions. Samples of the control and model groups were labeled with tags 126 and 127, respectively, while samples of FA and FB groups were labeled with tags 128 and 129, respectively. The labeled peptides were separated using a Shimadzu LC-20AB HPLC Pump system (Shimadzu, Kyoto, Japan), equipped with a high-pH reversed-phase (RP) column (5 μm, Phenomenex, CA, USA). The labeled peptides were firstly reconstituted to 2 mL with buffer A (5% ACN and 95% H2O, pH adjusted to 9.8 with ammonia) and then loaded on the column for separation with a gradient of 5% buffer B (5% H2O and 95% ACN, pH adjusted to 9.8 with ammonia) for 10 min, 5–35% buffer B for 40 min, and 35–95% buffer B for 1 min at a flow rate of 1 mL min−1. Finally, the system was maintained for 3 min at this condition, followed by a decrease to 5% buffer B within 1 min and equilibration in 5% buffer B for 10 min. The elution process was monitored by measuring absorbance at 214 nm. Fractions were collected every minute and pooled as 20 fractions. The eluted peptides were subsequently concentrated by vacuum centrifugation.
The fractions were resuspended in buffer A (2% acetonitrile and 0.1% formic acid) and then centrifuged for 10 min (20,000×g). The supernatants were loaded onto a C18 trap column on an LC-20 AD nano-HPLC instrument. An internally packed analytical C18 column was used to elute and separate the peptides. The elution was performed at a flow rate of 600 nL min−1 under the following conditions: 8–35% buffer B (5% H2O and 95% ACN, pH adjusted to 9.8 with ammonia), 0–35.00 min; 60% buffer B, 35.01–40.00 min; 80% buffer B, 40.01–45.00 min; 5% buffer B, 45.01–45.10 min. The equilibration time was 10 min. The eluted peptides were subjected to nano-electrospray ionization, followed by MS/MS (Orbitrap Fusion Lumos mass spectrometer, CA, USA) and nano HPLC analyses. The mass spectrometry analyses were performed with a scan range of 350–1800 m/z, and the survey scans were obtained at 120,000 mass resolution of 400 m/z using an Orbitrap analyzer. Dynamic exclusion parameters were a repeat count of 2 and a repeat duration of 30 s.
Data processing and proteins identification
The Proteome Discoverer software v2.1 (Thermo Fisher Scientific, MA, USA) was applied to process and quantify the raw data files. A protein search was performed in the RefSeq human protein database (24078 sequences, released in 2017) based on the SEQUEST algorithm. Parameters were set as follows: fixed modifications, including carbamidomethylation of cysteine, and modification at N-terminus and K of TMT six-plex. Methionine oxidation was applied as variable modification. Trypsin was set as the enzyme, allowing two missing cleavages. Verification and identification were performed using the Percolator software (Proteome Discoverer sequest (v2.1)). The protein and peptide profile matching the false discovery rate (FDR) was set at 0.01 . Total protein intensity was generated by summing all reported ion intensities of unique peptides that matched each protein. Corrections of the bias and background were performed by checking protein quantification and normalization. Proteins containing at least two unique peptides were further quantified.
Statistical analysis of proteomics
In proteomic analysis, proteins were extracted from three biological replications in each group. The fold changes of proteins in the control, model, FA, and FB groups were calculated as mean value according to the relative and absolute quantification ratio of the protein isobaric tags. Student’s t tests of SPSS 25.0 were performed to further determine whether the differential proteins were significantly different (p < 0.05) between groups. Accordingly, proteins with p < 0.05 and fold change > 1.2 were identified as differentially expressed biomarkers.
Interactive network construction of metabolites and proteins
For pathway and network analyses, the correlation coefficient and p value of differentially expressed metabolites and proteins (with NCBI IDs) were calculated by the R language. The metabolites and proteins with p < 0.05 and absolute correlation coefficient > 0.9 were screened out, and further imported into Cytoscape software (v3.4.0) to generate the final association analysis network diagram.
Reverse transcription quantitative real-time PCR (RT-qPCR) analysis
We used RT-qPCR to detect and quantify the mRNA expressions of nucleoside diphosphate kinase 3 isoform X1 (Nme3), WD repeat-containing protein 3 (Wdr3), 28S ribosomal protein S7 mitochondrial precursor (Mrps7), and Collagen. We also quantified the inflammatory mediators such as IL-6, IL-1β, TNF-α, and genes involved in nuclear transcription factor-kappa B (NF-κB), mitogen-activated protein kinases (MAPK), and Janus kinase/signal transducer and activator of transcription (JAK-STAT) signaling pathways. In brief, zebrafish larvae were washed three times with RNase-free water and homogenized with the TRIzol reagent. Homogenized tissues were extracted by chloroform and centrifuged for 15 min (12,000×g, 4 °C). The supernatant was transferred to a new tube, and isopropyl alcohol of equal volume was added for RNA precipitation. After being centrifuged for 10 min (12,000×g, 4 °C), 75% ethanol was used to wash the pellet, which was then dried and then suspended in 50 μL RNase-free water. The optical density (OD) at 260/280 nm was measured for RNA purity detection. The RT-qPCR reaction conditions were set as follows: 95 °C for 10 min, 40 cycles of 95 °C for 15 s and 60 °C for 30 s. The relative mRNA expressions levels were calculated by the 2−ΔΔCT method. All primers used for gene amplification were designed using the Primer-BLAST (NCBI), and their sequences were listed in Table S1. For gene expression analysis, one-way ANOVA of SPSS 25.0 was used for statistical comparisons. Data were shown as mean ± S.D and p < 0.05 were assumed for the probability level for statistical significance.
Inhibitory effects of FA and FB against CuSO4-induced inflammation in zebrafish
The primary lateral line system of zebrafish was established within 3 days after fertilization. Studies have previously reported that the addition of CuSO4 can swiftly destroy the hair cells in zebrafish lateral line neuromasts through oxidation and cell death . The zebrafish larvae of each treatment group were observed using a fluorescence microscope, and the number of fluorescent cells in the region of about ten cell diameters within the horizontal muscle was counted (Fig. 2a). In the control group, most of the neutrophils were assembled in the posterior blood island or caudal hematopoietic tissue, in which most leukocytes were distributed in this development stage (Fig. 2b). In contrast, zebrafish larvae exposed to CuSO4 for 40 min formed a typical neutrophils cluster in the horizontal muscles, which indicated that neutrophils could migrate to the lateral line neuromasts, with the induction of CuSO4. On the contrary, co-treatment with FA or FB inhibited neutrophils migration and decreased the assembled number of neutrophils, as shown in Fig. 2c (the red dotted lines represented the location of the horizontal muscles). These results indicated that FA and FB could dose-dependently inhibit the migration of neutrophils to the zebrafish lateral line neuromasts (Fig. 2d).
Inhibitory effects of FA and FB against CuSO4-induced intracellular ROS and NO generation in zebrafish
The overproduction of ROS by polymorphonuclear neutrophils at inflammatory sites might cause endothelial dysfunction and tissue damage . Therefore, DCF-DA was used to detect ROS production in CuSO4-induced inflammation in zebrafish. Fig. 3a was a representative image of intracellular ROS production. The control group, not treated with CuSO4, FA, or FB, exhibited a dark and weak fluorescent image. However, the model group exposed to CuSO4 presented a brighter and stronger fluorescent image. Thus, the addition of CuSO4 could lead to intracellular ROS generation in zebrafish. However, FA and FB dose-dependently decreased intracellular ROS accumulation in zebrafish (Fig. 3b).
Besides, NO is also a common conduction molecular and plays an essential role in both the initiation and development of inflammation . Therefore, the fluorescence probe DAF-FMDA was performed to measure the intracellular NO production induced by CuSO4, for further evaluating the inhibitory effects of FA and FB in zebrafish. As a result, zebrafish larvae exposed to CuSO4 exhibited a high generation of NO (Fig. 3c). On the contrary, zebrafish larvae treated with FA or FB showed a significant decrease of NO accumulation in a dose-dependent manner (Fig. 3d). These results indicated that FA and FB could relieve CuSO4-induced inflammation in zebrafish by inhibiting the intracellular generations of ROS and NO.
Multivariate analyses of UPLC-MS/MS data
UPLC-MS/MS was used to analyze all zebrafish larvae samples in both positive and negative ion modes. Representative base peak chromatograms (BPC) of zebrafish larvae samples from the control, model, FA, and FB groups were obtained under optimal conditions (Fig. 4a, b). Fifteen thousand four hundred eighty-nine variables (ESI+) and 26308 variables (ESI−) could simultaneously be detected in 20 min.
Multivariate analyses suggested clear separations from the control, model, FA, and FB groups. The principal component analysis (PCA) model and the profiles of the control, model, FA, and FB groups showed a tendency towards separation (Fig. 4c, d). The partial least squares-discriminant analysis (PLS-DA) was further performed to distinguish the differences of metabolites in the four groups, and the scores plot of PLS-DA indicated that these groups could be obviously distinguished (Fig. 4e, f). It was found that the control and model groups were separated in both positive mode and negative mode, indicating that a CuSO4-induced acute inflammation model in zebrafish larvae was successfully established. FA (120 μmol L−1) and FB (150 μmol L−1) groups were separated from the model group and were closer to the control group, which highlighted the anti-inflammatory effects of drugs.
Identification and analysis of metabolites
The heat-maps of metabolites changes of control vs. model, model vs. FA, and model vs. FB were shown in Fig. 4g–i. Based on the outcome of MS/MS and the information from online database, we carried out the identification of these metabolites. As a result, a total of 88 metabolites were found to be significantly altered in the model group relative to the control group. Treatment with FA and FB reversed the changes of 14 and 35 metabolic biomarkers, respectively. As listed in Table 1, 36 metabolic biomarkers were found to be involved in CuSO4-induced inflammation after treatment with FA or FB. Both FA and FB could reverse the abnormal expressions of the following metabolites involved in inflammatory response: uridine 5-diphosphate, carnitine, Val-Abu-OH, ciliatine, cytidine, benzaldehyde, isobutyryl carnitine, N-acetyl-L-phenylalanine, gamma-glutamyl-leucine, dihydrojasmonic acid, cyclic GMP, 12-hydroxydodecanoic acid, and N-arachidonoyl dihydroxypropylamine. FA and FB might exert anti-inflammatory effects by regulating the expressions of these metabolites [25,26,27,28,29,30].
Based on the information of the KEGG database, we performed the topology and pathway enrichment analyses to identify the involved metabolomic pathways. We also applied the metaboanalyst to calculate the −log(p) value and the pathway impact value based on pathway enrichment and topology analyses, respectively. According to the values of −log(p) and pathway impact, we finally identified the potential pathways associated with the effect of FA and FB in CuSO4-induced inflammation, summarized in Table S2. The identified biomarkers and related pathways were shown in Fig. 4j. Our findings pointed mainly to the involvement of nicotinate and nicotinamide metabolism, energy metabolism, pyrimidine metabolism, and purine metabolism.
UPLC-MS/MS was applied to analyze the proteomic profiles of CuSO4-induced inflammation in zebrafish after FA and FB treatment. As a result, 5212 proteins were detected, and 2976 proteins were further identified and quantified with at least two unique peptides and the false discovery rate (FDR) < 1%. According to the criteria of fold change > 1.2 and p < 0.05, there were 146 differentially expressed proteins (DEPs) in the model group, including 60 upregulated proteins and 86 downregulated proteins, compared with the control group. Additionally, FA and FB treatment reversed the expressions of 51 and 18 DEPs, respectively (Table 2).
Analyses of the DEPs
The heat-maps of changes in differential proteins of control vs. model, model vs. FA, and model vs. FB were presented in Fig. 5a–c. The DEPs were categorized according to the following Gene Ontology (GO) classes: biological process, molecular function, and cellular components (Fig. S1). KEGG analysis was further carried out to identify the biological pathways associated with CuSO4-induced inflammation in zebrafish, so as to clarify the therapeutic mechanism of FA and FB. As shown in Fig. 5d–f, the identified proteins were involved in a variety of pathways, including cellular processes, environmental information processing, genetic information processing, human diseases, metabolism, and organismal systems, covering a wide range of biological pathways in the inflammatory response and neuromasts damage.
The identified DEPs were imported to the STRING database (https://string-db.org/) for protein-protein interaction (PPI) network constructions (Fig. S2). Compared with the control group, 146 DEPs were identified in the model group, among these DEPs, 101 proteins were found to be involved in the network, 40 proteins were related to each other, and 5 proteins did not display any linkage at a confidence level of string score = 0.4. This network reflected complex functional relationships among the identified proteins in the present study.
Interactive network construction
The differentially expressed metabolites and proteins from control vs. model, model vs. FA, and model vs. FB were imported into Cytoscape software to conduct the significant network analyses (Fig. S3). Further integrated analyses of metabolomic and proteomic studies exhibited the altered pathways in response to CuSO4-induced acute inflammation, including lipid metabolism, amino acid metabolism, and Nucleotide metabolism (Fig. 6a–c). FA and FB were likely to reverse expressions of metabolites and proteins involved in these metabolic pathways, thus alleviating the injury of hair cells in zebrafish.
We used RT-qPCR to detect and quantify the mRNA levels of Wdr3, Nme3, collagen, and Mrps7. We also quantified the inflammatory mediators IL-6, IL-1β, and TNF-α that contribute to the early inflammatory phase, and mediate the recruitment of other immune cells to the damaged areas. Finally, we quantified the genes involved in the NF-κB, MAPK, and JAK-STAT signaling pathways. As a result, gene expressions of Wdr3, Nme3, collagen, and Mrps7 in FA and FB groups were decreased, compared with the model group, which was consistent with proteomic results (Fig. 7a). Moreover, compared with the model group, FA and FB significantly downregulated the expression of IL-6, IL-1β, and TNF-α (Fig. 7B). FB could reverse mRNA expressions of genes involved in NF-κB, MAPK, and JAK/STAT signaling pathways. However, FA mainly exerted an anti-inflammatory effect through NF-κB and MAPK signaling pathways (Fig. 7c).
The inflammatory response is a common and highly regulated biological process in the body that helps to clear harmful irritants and effectively initiate reparation . In the process of inflammation, excessive recruitment and uncontrolled migration of neutrophils to injured sites could result in tissue damage and disease . In this study, we used CuSO4 to induce inflammation in zebrafish so that we could monitor the behavior of the neutrophils during the critical transitions phase. CuSO4 is a metal chemical that induces neuromasts damage in the zebrafish lateral line system, and can easily induce the characteristics of oxidative stress and inflammation. In addition to their regenerative capacity, the structure, function, and molecular activity in the zebrafish lateral line hair cells are similar to hair cells in mammals . In the present study, neutrophils have significantly accumulated in the damaged neuromasts following exposure to CuSO4. However, the number of neutrophils in the neuromasts was clearly reduced when zebrafish were treated with FA or FB, which indicates that by exhibiting anti-inflammatory activities, FA and FB might exert excellent neuroprotective properties against CuSO4-induced neuroinflammation.
Excessive ROS and NO production in the body might lead to cell and tissue damage. Therefore, inhibiting this excessive production can effectively impede the progress of the inflammatory response. Such inhibition could be considered as a potential anti-inflammatory drug development target. Exposure to CuSO4 significantly increases ROS and NO production in zebrafish. As previously reported, the oxidative stress caused by the overproduction of ROS can continuously stimulate the neurons, leading to their damage and dysfunction [33, 34]. However, FA and FB have markedly inhibited ROS and NO accumulation in the CuSO4-induced damaged neuromasts, which further suggests the remarkable neuroprotective properties of FA and FB, preventing the inflammation induced by CuSO4 exposure.
Metabolomic analysis has identified multiple metabolites that differed between the control and model groups. Among them, treatment with FA and FB reversed the abnormal expressions of uridine 5-diphosphate, carnitine, Val-Abu-OH, ciliatine, cytidine, benzaldehyde, isobutyryl carnitine, N-acetyl-L-phenylalanine, gamma-glutamyl-leucine, dihydrojasmonic acid, cyclic GMP, 12-hydroxydodecanoic acid, and N-arachidonoyl dihydroxypropylamine. Some metabolites have been reported to participate in the modulation of particular inflammation-related molecules or mediators and play a crucial role in the inflammatory cascade reaction. For example, uridine 5-diphosphate can induce the expression of chemokines such as monocyte chemotactic protein-1 (MCP-1) and macrophage inflammatory protein-1α (MIP-1α) in microglia and astrocytes. These chemokines are the primary effector molecules that mediate the recruitment of inflammatory cells to the damaged tissues [25, 35]. Therefore, the potential role of uridine 5-diphosphate in the inflammatory response is likely to be related to the involvement of chemokines in inflammation. The depletion of carnitine was previously reported to be prevalent in hemodialysis patients, which was assumed to be associated with their inflammatory response in patients . Moreover, exogenous carnitine was required to satisfy energy consumption, which might reflect on carnitine’s indirect modulation of the inflammatory response by regulating energy metabolism. The specific role of energy metabolism in inflammation will be more thoroughly discussed later. Cytidine can be metabolized to uridine through the action of cytidine deaminase. This action activates and increases nucleoside-sensing of the Toll-like receptor 8, which is supported by the overexpression of TNF-α in monocytes or myeloid cells [37, 38]. Thus, the involvement of cytidine in inflammation might be related to the inflammatory cascade reaction induced by the activation of Toll-like receptors.
Further metabolomic pathway enrichment analysis pointed mainly to the involvement of nicotinate and nicotinamide metabolism, energy metabolism, pyrimidine metabolism, and purine metabolism. Most of these pathways are directly or partly involved in inflammation [39,40,41,42,43]. The involvement of nicotinate and nicotinamide metabolism in the inflammatory response might be related to some metabolic enzymes as well as secondary metabolites produced during the metabolic processes. For example, nicotinate phosphoribosyltransferase is an important inflammatory mediator that can bind to Toll-like receptor 4, which in turn induces the activation of inflammasomes and NF-κB. This activation leads to the secretion of the inflammatory mediators IL-1β, IL-8, and TNF-α . Moreover, the identified metabolite nicotinamide ribotide participates in nicotinate and nicotinamide metabolism. This metabolite, also known as nicotinamide mononucleotide (NMN), can be converted into the principal precursor of nicotinamide adenine dinucleotide (NAD+) through transamination [44, 45]. To activate Sirtuin 1 (SIRT1), NAD+ regulates the ratio between NAD+ and nicotinamide adenine dinucleotide (NADH). SIRT1, in turn, positively regulates NF-κB expression, and NF-κB is a crucial transcription factor involved in pro-inflammatory cytokines generation [46, 47]. Presently, it can be speculated that FB might decrease NAD+ synthesis by inhibiting the expression of nicotinamide ribotide. This cascade may further result in an inactivation of SIRT1 and the NF-кB signaling pathway, along with a decrease in the release of inflammatory factors.
Differentially expressed metabolites such as cytidine, cytosine, deoxyuridine, uracil, XMP, dIMP, cyclic GMP, and guanosine were identified in the present study. They are all involved in pyrimidine and purine metabolism. To reach a reasonable comprehension of the involvement of pyrimidine and purine metabolism in the inflammatory response, the participation of pyrimidine and purine nucleotide receptors (P2Y receptors) had to be well explained . That is, P2Y receptors are present practically in all immune cells, and can specifically mediate the inflammatory response to various cell injuries by recognizing pyrimidine and purine nucleotides . For example, the P2Y6 receptor shows the highest expression in macrophages, dendritic cells, and microglia. The activation of this receptor could lead to immune phenotypes such as macrophage and microglial dysfunction, and inactivation of effector T cells. It also plays a significant role in vascular inflammation [48,49,50]. Moreover, the inhibition of P2Y6 receptor was shown to be accompanied by a decrease in the release of IL-8, IL-1α, and IL-6, a release that was induced by monosodium urate crystals in human keratinocytes . Therefore, the significant role of pyrimidine and purine metabolism in the inflammatory response might be closely related to the immune functions played by various P2Y receptors. Consequently, it can be speculated that FA and FB might exert their anti-inflammatory effects by regulating these metabolites, and through them modulating cellular immune functions and the release of various inflammatory mediators.
Energy consumption also plays an essential role in the development of inflammation. In the present study, Starch and sucrose metabolism, amino sugar and nucleotide sugar metabolism, galactose metabolism, as well as the pentose phosphate pathway (PPP) and the glycolysis, have all been shown to be related to energy metabolism. During the inflammatory response, there is an increased demand for energy by the immune cells. This extra energy is needed to satisfy the demand by housekeeping functions and multiple immune cells-specific functions. Such functions include cell migration, cytokinesis, antigen processing and presentation, and lymphocyte activation [38, 52]. The inflammatory response can increase energy consumption by as much as 25–60% [53, 54]. As reported, the use of drugs to regulate the energy metabolism of immune cells could be a potential treatment strategy [55, 56]. For example, glucose-6-phosphate (G6P) is at the intersection of multiple energy metabolism pathways, through which it can be converted into PPP for nicotinamide adenine dinucleotide phosphate (NADPH) and ribose-5-phosphate. The activation of PPP contributes to the inhibition of P53. This inhibition leads to increased glucose consumption, NADPH production, and biosynthesis . In our study, G6P was found to be significantly upregulated in the model group relative to the control group. Treatment with FB remarkably downregulated G6P, suggesting that FB might serve as an anti-inflammatory agent by regulating the expression of G6P, as well as the reduction in energy consumption. Meanwhile, the relieved damage in the FB group might also contribute to the decreased production of G6P and the reduced requirement of energy consumption. That is, in CuSO4-induced acute inflammation, zebrafish need more energy to conduct a series of inflammation-related activities and immune cells activation. They also need it to repair the damaged tissues. Zebrafish larvae in the FA and FB groups tended to require less energy, which was basically used to maintain housekeeping functions. In other words, the reduced energy consumption might be secondary to the reduced inflammatory response in the FA and FB groups.
The proteomic analysis has identified 146 DEPs when comparing the model group and control groups. Among these DEPs, 51 and 18 proteins were reversed in the FA and FB groups, respectively. The KEGG pathway classification showed that the DEPs covered a variety of cellular processes and metabolic pathways involved in the inflammatory response and neuromast damage. The identified DEPs participate in multiple biological processes and categories such as immune system and infectious disease, nervous system and neurodegenerative disease, signal transduction, and cell growth and death. Our data suggest that the exposure to CuSO4 could have damaged both the immune and the nervous systems of the zebrafish. Such damage involves a variety of signal transduction as well as cell death and regeneration. The identified proteins also participate in various metabolomic pathways, including amino acid metabolism, lipid metabolism, nucleotide metabolism, and energy metabolism. These pathways could influence multiple inflammation-related processes, as will be further elaborated in the discussion below on integrating the metabolomic and proteomic analyses. The changes in expression of the identified DEPs that are actively or passively involved in the inflammatory response, including collagen (col2a1b, col9a2, col9a1b), Wdr3, Mrps7, and Nme3, have all been reversed both in the FA and FB groups [58,59,60,61]. Previous studies have shown that the organization of ROS, NF-κB, and collagen at the injury margins is involved in promoting wound healing . High levels of TNF-α might lead to an excessive inflammatory response, along with increased deposition of collagen in fibroblast-treated wounds . In addition to the generation of collagen, remodeling its structure, as induced by tissue damage, also plays a crucial role in the repair process. This remodeling might be achieved by proteases such as the matrix metalloproteinase (MMP) [64, 65]. The critical role of MMP9 in collagen reorganization and regeneration has been previously reported. One report suggested that MMP9 could aggravate wound damage in chronic injury, but promote wound healing and regeneration in acute injury . This implies a critical role for collagen in tissue damage and repair. Our findings suggested that both FA and FB might guide treatments focused on prevention (neuroprotection) or repair (anti-inflammatory) of CuSO4-induced hair cells damage in the zebrafish larvae by modulating the productions of collagen.
Mitochondrial ribosomal protein (MRP) subunits are nuclear-encoded and were crucial for mitochondrial functions and mitochondrial-encoded protein synthesis. These subunits were also reported to be involved in mitochondrial diseases [61, 67]. Mrps7, a 12S ribosomal RNA binding subunit, is necessary for the assembly of small ribosomal subunits. Studies have confirmed that Mrps7 mutations could cause mitochondrial respiratory chain dysfunction and congenital sensorineural deafness . Mitochondria are an essential component of the innate inflammatory response. Their dysfunction was shown to be a vital trigger of inflammation and could result in metabolic alteration and mitochondrial damage-induced release of inflammatory mediators . Thus, we speculated that alterations in Mrps7 expression induced by CuSO4 exposure might be related to mitochondrial dysfunction in the zebrafish. This dysfunction could be the main source of ROS following the increase in energy demands and efficient mitochondrial function-dependent excessive inflammatory activities .
As a member of the WD repeat-containing proteins, Wdr3 participates in several cellular processes such as cell cycle progression and signal transduction [71, 72]. Studies have confirmed that some Wdr proteins participate in the MAPK and STAT3 signaling pathways . Moreover, the expression of misfolded Wdr proteins in dendritic cells could activate the inflammasome and result in a substantial release of IL-18 . Here, we can speculate that FA and FB may have reversed the abnormal expression of Wdr3, which in turn regulate Wdr proteins-mediated cell cycle progression, signal transduction, and the release of inflammatory mediators from the inflammasomes. Nme3, also known as NM-23, participates in a variety of physiological and pathological cellular processes such as differentiation, development, cellular signaling, and cellular function. It was reported that Nm23-H1 could interact with p53 and positively mediate the apoptosis and cell cycle arrest that was induced by p53 . Besides, Nme3 was confirmed to modulate pro-inflammatory transcription by activating NF-κB . This modulation activates the Toll-like receptor 5-mediated NF-κB signaling pathway in a MyD88-dependent manner . In the present study, both FA and FB downregulated Nme3 expression, suggesting that FA and FB might possess a regulatory effect on cellular signaling and function.
Integration of the metabolomic and proteomic analyses indicates that the lipid metabolism, amino acid metabolism, and nucleotide metabolism are significantly altered as a result of CuSO4 exposure. The analyses helped to form a more comprehensive understanding of the metabolic pathways involved in CuSO4-induced inflammation in zebrafish. Moreover, pathways identified through the metabolomic analysis, including the pyrimidine and purine metabolism, are indispensable during the metabolism of lipids, amino acid, and nucleotides. Through these metabolic pathways, they perform their crucial role in lipids activation, energy conservation, and nucleotide synthesis . Presently, lipids are considered as an effective signaling molecule that regulates multiple cellular responses. A variety of lipid-derived mediators are involved in the regulation of inflammation and the response to infection [77, 78]. For example, the lipid-derived mediators, prostaglandin E2 (PGE2) and prostacyclin (PGI2), function as vasodilatation agents during the initial stage of inflammation. These mediators help to recruit various immune cells from the blood, including neutrophils and macrophages, and mobilize them to the infected or damaged tissues . Therefore, it can be speculated that FA and FB might exert therapeutic effects by regulating lipid metabolism. In turn, this regulation might affect the role of lipid-derived mediators in the inflammatory response. Amino acid metabolism is also considered to mediate the immune response in mammals. Their metabolism might be necessary for pivotal enzymes, and their degradation could be crucial immune checkpoints in autoimmunity . The activation of enzymes such as indoleamine 2,3-dioxygenase 1 (IDO1) and arginase 1 (ARG1) could induce an immunoregulatory effect on dendritic cells and inflammatory cells. This effect could coordinate the immune response and resistance to inflammation [81, 82]. Additionally, nucleotide metabolism could occur under inflammatory conditions and following tissue injury . The various nucleotides released to the extracellular domain, including ATP, ADP, and UDP, could specifically bind with individual P2Y receptors. These receptors, in turn, participate in handling tissue damage and the recruitment of immune cells [83,84,85]. Therefore, FA and FB might exert their anti-inflammatory activities by modulating the amino acid metabolism and nucleotide metabolism, and thus reducing the immune response and alleviating the damage caused to the neuromasts.
Moreover, in the present study, we also detected changes in mRNA expression of genes involved in the NF-κB, MAPK, and JAK-STAT inflammatory signaling pathways, and in the production or release of inflammatory cytokines such as IL-6, IL-1β, and TNF-α. NF-κB, MAPK, and JAK-STAT pathways are involved in the signal transduction of the inflammatory response, including regulation of the secretion of various chemokines and inflammatory cytokines. These contribute to an amplified inflammatory response and worsening of the tissue damage [86,87,88]. Our findings suggested that FA and FB might ameliorate inflammation by negatively interfering with these signaling pathways. The inflammatory cytokines IL-6, IL-1β, and TNF-α are the most studied pro-inflammatory cytokines. They could contribute to the early inflammatory phase and participate in multiple inflammation-related processes. These effects are partly mediated by their crucial role in the recruitment of additional inflammatory cells to the infected or injured tissues . Furthermore, IL-6, TNF-α, and IL-1β were also reported to play a significant role in neurocyte damage-induced neuropathic pain . Immediately after nerve injury, neutrophils and macrophages are recruited to the injured area, resulting in a substantial release of cytokines such as IL-6, TNF-α, and IL-1β. These cytokines can then directly regulate the activity of the neurons [89, 90]. In the present study, CuSO4 was used to specifically induce neuromasts damage in the lateral line system of zebrafish. This damage initiated the infiltration of immune cells such as neutrophils, and possibly also the release of inflammatory mediators such as IL-6, TNF-α, and IL-1β. However, both FA and FB were found to inhibit the recruitment of neutrophils and the secretion of inflammatory mediators induced by CuSO4. Thus, FA and FB were able to reduce inflammation and alleviate neuromasts damage. It was further shown that through anti-inflammatory activities, FA and FB exhibit neuroprotection against the CuSO4-induced neuroinflammatory response.
Presently, we applied the integration of metabolomics and proteomics to explore the anti-inflammatory properties and mechanisms of FA and FB against CuSO4-induced neuromasts damage in zebrafish larvae. Our results indicated that both FA and FB significantly inhibited the injury of CuSO4-induced hair cells in zebrafish neuromasts. Based on our study, we presumed that FA and FB might exhibit neuroprotective properties in the damage of CuSO4-induced zebrafish neuromasts, evidenced by the reduced recruitment of neutrophils and inhibited expressions of ROS, NO, and inflammatory mediators including IL-6, IL-1β, and TNF-α. Metabolomic and proteomic analyses indicated that FA and FB might share similar anti-inflammatory mechanisms to some extent because of the same maternal nucleus structure. The different substituents of hydroxyl groups of FA and FB might indicate a minor difference in anti-inflammatory activities: FA exhibited no effect on JAK/STAT signaling pathway. On the contrary, FB could reverse the overexpression of STAT3 and JAK3.
In summary, we proposed a potential mechanism of CuSO4-induced inflammation and the action diagram of FA and FB (Fig. 8). FA and FB administration resulted in reductions of ROS and NO, gene expressions of Wdr3, collagen, Nme3 and Mrps7, and other inflammatory mediators induced by CuSO4, thus reducing inflammation. As discussed above, a series of metabolites, metabolic pathways, proteins, and gene products were altered in response to CuSO4 exposure. However, FA and FB treatments significantly reversed these changes through a complicated regulatory mechanism. Our findings suggested that these identified metabolites and proteins may serve as potential anti-inflammatory biomarkers. FA and FB exhibited excellent anti-inflammatory activities against CuSO4-induced inflammation in zebrafish. Due to the multi-component and multi-regulatory process in response to CuSO4 exposure, FA and FB might influence downstream conduction process by affecting some of the upstream inflammatory processes. For instance, reductions in NO and ROS may lead to secondary effects on their downstream products, such as the decreased expressions of downstream inflammatory mediators IL-6, IL-1β, and TNF-α via modulating specific inflammatory signaling pathways. Accordingly, the reversed expressions of Collagen, Nme3, Wdr3, and Mrps7 in FA and FB group against CuSO4 exposure may also be secondary to the regulation of some metabolic pathways upstream. Also, the fact that zebrafish were both pre- and post-treated with FA or FB may provide another potential therapeutic mechanism, that is, pretreatment may attenuate the initiation of inflammation wherein post-treatment may inhibit the inflammatory cascade as well. Both administration methods will inhibit the inflammatory response through different mechanisms and have different therapeutic significance, which can be further studied in the future.
In conclusion, we analyzed the inhibitory effect and therapeutic mechanism of FA and FB against CuSO4-induced inflammation in zebrafish using UPLC-MS/MS based metabolomics and TMT-6 Plex based proteomics for the first time. The alterations of identified metabolic intermediates and protein biomarkers might help to elucidate the anti-inflammatory mechanisms of FA and FB against CuSO4-induced inflammation in zebrafish. Our results also provided a potential novel and sequential mechanism for the initiation of inflammation induced by CuSO4 and suggested a multi-component and multi-regulatory therapeutic mechanism for FA and FB.
Availability of data and materials
The datasets generated and/or analyzed in the current study are available from the corresponding author on reasonable request.
Tumor necrosis factor-α
Diamino fluorophore 4-amino-5-methylamino-2′,7′-difluorofluorescein diacetate
Reactive oxygen species
Liquid chromatography-mass spectrometry
High-performance liquid chromatography
Electrospray ionization mass detection
Mass to charge ratio
Mass spectrometry/mass spectrometry
Higher energy collisional dissociation
Ultra-performance liquid chromatography-mass spectrometry/mass spectrometry
Orthogonal projections to latent structures discriminant analysis
Variable importance in the projection
Part per million
Human metabolome database
Kyoto Encyclopedia of Genes and Genomes Database
Tandem mass tags
Nuclear transcription factor-kappa B
Mitogen-activated protein kinases
Janus kinase/signal transducer and activator of transcription
Base peak chromatogram
Principal component analysis
Partial least squares-discriminant analysis
Differentially expressed proteins
False discovery rate
Reverse transcription quantitative real-time PCR
Monocyte chemotactic protein-1
Macrophage inflammatory protein-1α
Nicotinamide single nucleotide
Nicotinamide adenine dinucleotide
Nicotinamide adenine dinucleotide
Pyrimidine and purine nucleotide
Pentose phosphate pathway
Nicotinamide adenine dinucleotide phosphate
Mitochondrial ribosomal protein
Indoleamine 2,3-dioxygenase 1
- PGE2 :
- PGI2 :
Amor S, Peferoen LA, Vogel DY, Breur M, Valk P, Baker D, et al. Inflammation in neurodegenerative diseases--an update. Immunology. 2014;142(2):151–66.
Freire MO, Van Dyke TE. Natural resolution of inflammation. Periodontol. 2013;63(1):149–64.
Muralidharan S, Mandrekar P. Cellular stress response and innate immune signaling: integrating pathways in host defense and inflammation. J Leukocyte Biol. 2013;94(6):1167–84.
Chen HY, Lin YH, Huang JW, Chen YC. Chinese herbal medicine network and core treatments for allergic skin diseases: implications from a nationwide database. J Ethnopharmacol. 2015;168:260–7.
Pan CW, Zhou GY, Chen WL. Protective effect of forsythiaside a on lipopolysaccharide/d-galactosamine-induced liver injury. Int Immunopharmacol. 2015;26(1):80–5.
Jiang WL, Fu FH, Xu BM, Tian JW, Zhu HB, Hou J. Cardioprotection with forsythoside b in rat myocardial ischemia-reperfusion injury: relation to inflammation response. Phytomedicine. 2010;17(8–9):635–9.
MacRae CA, Peterson RT. Zebrafish as tools for drug discovery. Nat Rev Drug Discov. 2015;14(10):721–31.
Dooley K, Zon LI. Zebrafish: a model system for the study of human disease. Curr Opin Genet Dev. 2000;10(3):252–6.
García-Moreno D, Tyrkalska SD, Valera-Pérez A, Gómez-Abenza E, Pérez-Oliva AB, Mulero V. The zebrafish: A research model to understand the evolution of vertebrate immunity. Fish Shellfish Immunol. 2019;90:215–22.
Lee SH, Ko CI, Jee Y, Jeong Y, Jeon YJ. Anti-inflammatory effect of fucoidan extracted from ecklonia cava in zebrafish model. Carbohyd Polym. 2013;92(1):84–9.
Meeker ND, Trede NS. Immunology and zebrafish: spawning new models of human disease. Dev Comp Immunol. 2008;32(7):0–757.
Ghysen A, Dambly-Chaudière C. Development of the zebrafish lateral line. Curr Opin Neurobiol. 2004;14(1):67–73.
d'Alençon CA, Peña OA, Wittmann C, Gallardo VE, Jones RA, Loosli F, et al. A high-throughput chemically induced inflammation assay in zebrafish. BMC Biol. 2010;8:1–16.
Kimmel CB, Ballard WW, Kimmel SR, Ullmann B, Schilling TF. Stages of embryonic development of the zebrafish. Dev Dyn. 1995;203(3):253–310.
Haffter P, Granato M, Brand M, Mullins MC, Hammerschmidt M, Kane DA, et al. The identification of genes with unique and essential functions in the development of the zebrafish, Danio rerio. Development. 1996;123:1–36.
Yang Y, Zheng K, Mei W, Wang Y, Hu J. Anti-inflammatory and proresolution activities of bergapten isolated from the roots of ficus hirta in an in vivo zebrafish model. Biochem Biophys Res Commun. 2018;496(2):763–9.
Trygg J, Wold S. Orthogonal projections to latent structures (O-PLS). J Chemometrics. 2002;16:119–28.
Wang JB, Pu SB, Sun Y, Li ZF, Niu M, Yan XZ, et al. Metabolomic Profiling of Autoimmune Hepatitis: The Diagnostic Utility of Nuclear Magnetic Resonance Spectroscopy. J Proteome Res. 2014;13(8):3792–801.
Unwin RD, Griffiths JR, Whetton AD. Simultaneous analysis of relative protein expression levels across multiple samples using itraq isobaric tags with 2d nano lc–ms/ms. Nat Protoc. 2010;5(9):1574–82.
Crowell AM, Wall MJ, Doucette AA. Maximizing recovery of water-soluble proteins through acetone precipitation. Anal Chim Acta. 2013;796:48–54.
Käll L, Canterbury JD, Weston J, Noble WS, Maccoss MJ. Semi-supervised learning for peptide identification from shotgun proteomics datasets. Nat Methods. 2007;4(11):923–5.
Olivari FA, Hernández PP, Allende ML. Acute copper exposure induces oxidative stress and cell death in lateral line hair cells of zebrafish larvae. Brain Res. 2008;1244:1–12.
Mills EL, O'Neill LA. Reprogramming mitochondrial metabolism in macrophages as an anti-inflammatory signal. Eur J Immunol. 2016;46(1):13–21.
Tripathi P, Tripathi P, Kashyap L, Singh V. The role of nitric oxide in inflammatory reactions. Fems Immunol Med Mic. 2007;51(3):443–52.
Kim B, Jeong HK, Kim JH, Lee SY, Jou I, Joe EH. Uridine 5'-diphosphate induces chemokine expression in microglia and astrocytes through activation of the P2Y6 receptor. J Immunol. 2011;186(6):3701–9.
Shimizu T, Marusawa H, Endo Y, Chiba T. Inflammation-mediated genomic instability: roles of activation-induced cytidine deaminase in carcinogenesis. Cancer Sci. 2012;103(7):1201–6.
Khalatbari-Soltani S, Tabibi H. Inflammation and L-carnitine therapy in hemodialysis patients: a review. Clin Exp Nephrol. 2015;19(3):331–5.
Millet P, Opiekun M, Martin T, Beauchamp GK, Kimball BA. Cytokine contributions to alterations of the volatile metabolome induced by inflammation. Brain Behav Immun. 2018;69:312–20.
Li QJ, Wang ZG, Xie Y, Liu Q, Hu HL, Gao YX. Mechanistic evaluation of gastro-protective effects of KangFuXinYe on indomethacin-induced gastric damage in rats. Chin J Nat Med. 2020;18(1):47–56.
Wu J, Sun L, Chen X, Du F, Shi H, Chen C, et al. Cyclic GMP-AMP is an endogenous second messenger in innate immune signaling by cytosolic DNA. Science. 2013;339(6121):826–30.
Gernez Y, Tirouvanziam R, Chanez P. Neutrophils in chronic inflammatory airway diseases: can we target them and how? Eur Respir J. 2010;35(3):467–9.
Christine DC, Dora S, Fabien S, Kelly D, Nicolas G, Alain G. The lateral line of zebrafish: a model system for the analysis of morphogenesis and neural development in vertebrates. Biol Cell. 2003;95(9):579–87.
de Oliveira DM, Ferreira Lima RM, El-Bachá RS. Brain rust: recent discoveries on the role of oxidative stress in neurodegenerative diseases. Nutr Neurosci. 2012;15(3):94–102.
Wang J, Song Y, Chen Z, Leng SX. Connection between Systemic Inflammation and Neuroinflammation Underlies Neuroprotective Mechanism of Several Phytochemicals in Neurodegenerative Diseases. Oxid Med Cell Longev. 2018;2018:1972714.
Fernandez EJ, Lolis E. Structure, function, and inhibition of chemokines. Annu Rev Pharmacol Toxicol. 2002;42:469–99.
Csiky B, Bene J, Wittmann I, Sulyok E, Melegh B. Effect of hemodialysis session on the dynamics of carnitine ester profile changes in L-carnitine pretreated end-stage renal disease patients. Int Urol Nephrol. 2013;45(3):847–55.
Micozzi D, Carpi FM, Pucciarelli S, Polzonetti V, Polidori P, Vilar S, et al. Human cytidine deaminase: a biochemical characterization of its naturally occurring variants. Int J Biol Macromol. 2014;63:64–74.
Furusho K, Shibata T, Sato R, Fukui R, Motoi Y, Zhang Y, et al. Cytidine deaminase enables Toll-like receptor 8 activation by cytidine or its analogs. Int Immunol. 2019;31(3):167–73.
Ma Y, Bao Y, Wang S, Li T, Chang X, Yang G, et al. Anti-Inflammation Effects and Potential Mechanism of Saikosaponins by Regulating Nicotinate and Nicotinamide Metabolism and Arachidonic Acid Metabolism. Inflammation. 2016;39(4):1453–61.
Managò A, Audrito V, Mazzola F, Sorci L, Gaudino F, Gizzi K, et al. Extracellular nicotinate phosphoribosyltransferase binds Toll like receptor 4 and mediates inflammation. Nat Commun. 2019;10(1):4116.
Chokesuwattanaskul S, Phelan MM, Edwards SW, Wright HL. A robust intracellular metabolite extraction protocol for human neutrophil metabolic profiling. PLoS One. 2018;13(12):e0209270.
Buttgereit F, Burmester GR, Brand MD. Bioenergetics of immune functions: fundamental and therapeutic aspects. Immunol Today. 2000;21(4):192–9.
Le Duc D, Schulz A, Lede V, Schulze A, Thor D, Brüser A, et al. P2Y Receptors in Immune Response and Inflammation. Adv Immunol. 2017;136:85–121.
Penberthy WT. Editorial [hot topic: nicotinamide adenine dinucleotide biology and disease (executive editor: w. todd penberthy)]. Curr Pharm Design. 2009;15(1):1–2.
Zhai RG, Rizzi M, Garavaglia S. Nicotinamide/nicotinic acid mononucleotide adenylyltransferase, new insights into an ancient enzyme. Cell Mol Life Sci. 2009;66(17):2805–18.
Matsushita T, Sasaki H, Takayama K, Ishida K, Kuroda R. The overexpression of sirt1 inhibited osteoarthritic gene expression changes induced by interleukin-1 in human chondrocytes. J Orthop Res. 2013;31(4):531–7.
Imai SI. The nad world: a new systemic regulatory network for metabolism and aging—sirt1, systemic nad biosynthesis, and their importance. Cell Biochem Biophys. 2009;53(2):65–74.
Bar I, Guns PJ, Metallo J, Cammarata D, Wilkin F, Boeynams JM, et al. Knockout mice reveal a role for P2Y6 receptor in macrophages, endothelial cells, and vascular smooth muscle cells. Mol Pharmacol. 2008;74(3):777–84.
Giannattasio G, Ohta S, Boyce JR, Xing W, Balestrieri B, Boyce JA. The purinergic G protein-coupled receptor 6 inhibits effector T cell activation in allergic pulmonary inflammation. J Immunol. 2011;187(3):1486–95.
Riegel AK, Faigle M, Zug S, Rosenberger P, Robaye B, Boeynaems JM, et al. Selective induction of endothelial P2Y6 nucleotide receptor promotes vascular inflammation. Blood. 2011;117(8):2548–55.
Uratsuji H, Tada Y, Kawashima T, Kamata M, Hau CS, Asano Y, et al. P2Y6 receptor signaling pathway mediates inflammatory responses induced by monosodium urate crystals. J Immunol. 2012;188(1):436–44.
Besedovsky HO, del Rey A. Immune-neuro-endocrine interactions: facts and hypotheses. Endocr Rev. 1996;17(1):64–102.
Straub RH, Cutolo M, Buttgereit F, Pongratz G. Energy regulation and neuroendocrine-immune control in chronic inflammatory diseases. J Intern Med. 2010;267(6):543–60.
Torine IJ, Denne SC, Wright-Coltart S, Leitch C. Effect of late-onset sepsis on energy expenditure in extremely premature infants. Pediatr Res. 2007;61(5 Pt 1):600–3.
Seki M, Ushiyama C, Seta N, Abe K, Fukazawa T, Asakawa J, et al. Apoptosis of lymphocytes induced by glucocorticoids and relationship to therapeutic efficacy in patients with systemic lupus erythematosus. Arthritis Rheum. 1998;41(5):823–30.
Buttgereit F, Brink I, Thiele B, Burmester GR, Hiepe F, Hall ED. Effects of methylprednisolone and 21-aminosteroids on mitogen-induced interleukin-6 and tumor necrosis factor-alpha production in human peripheral blood mononuclear cells. J Pharmacol Exp Ther. 1995;275(2):850–3.
Jiang P, Du W, Wang X, Mancuso A, Gao X, Wu M, et al. P53 regulates biosynthesis through direct inactivation of glucose-6-phosphate dehydrogenase. Nat Cell Biol. 2011;13(3):310–6.
Cui SJ, Fu Y, Liu Y, Kou XX, Zhang JN, Gan YH, et al. Chronic inflammation deteriorates structure and function of collagen fibril in rat temporomandibular joint disc. Int J Oral Sci. 2019;11(1):2.
Flentie K, Gonzalez C, Kocher B, Wang Y, Zhu H, Marasa J, et al. Nucleoside Diphosphate Kinase-3 (NME3) Enhances TLR5-Induced NFκB Activation. Mol Cancer Res. 2018;16(6):986–99.
Pfajfer L, Mair NK, Heredia RJ, Genel F, Gulez N. Ardeniz Ömür, et al. Mutations affecting the actin regulator WD repeat-containing protein 1 lead to aberrant lymphoid immunity. J Allergy Clin Immunol. 2018;142(5):1589–604.
Li HB, Wang RX, Jiang HB, Zhang ED, Tan JQ, Xu HZ, et al. Mitochondrial Ribosomal Protein L10 Associates with Cyclin B1/Cdk1 Activity and Mitochondrial Function. DNA Cell Biol. 2016;35(11):680–90.
LeBert D, Squirrell JM, Freisinger C, Rindy J, Golenberg N, Frecentese G, et al. Damage-induced reactive oxygen species regulate vimentin and dynamic collagen-based projections to mediate wound repair. Elife. 2018;7:e30703.
Bodnar RJ, Yang T, Rigatti LH, Liu F, Evdokiou A, Kathju S, et al. Pericytes reduce inflammation and collagen deposition in acute wounds. Cytotherapy. 2018;20(8):1046–60.
Cheng F, Shen Y, Mohanasundaram P, Lindström M, Ivaska J, Ny T, Eriksson JE. Vimentin coordinates fibroblast proliferation and keratinocyte differentiation in wound healing via tgf-β–slug signaling. Proc Natl Acad Sci USA. 2016;113(30):4320–7.
Eckes B, Colucciguyon E, Smola H, Nodder S, Babinet C, Krieg T, et al. Impaired wound healing in embryonic and adult mice lacking vimentin. J Cell Sci. 2000;113(13):2455–62.
Lebert DC, Squirrell JM, Rindy J, Broadbridge E, Lui Y, Zakrzewska A, et al. Matrix metalloproteinase 9 modulates collagen matrices and wound repair. Development. 2015;142(12):2136–46.
Brodersen DE, Nissen P. The social life of ribosomal proteins. FEBS J. 2005;272(9):2098–108.
Menezes MJ, Guo Y, Zhang J, Riley LG, Cooper ST, Thorburn DR, et al. Mutation in mitochondrial ribosomal protein s7 (mrps7) causes congenital sensorineural deafness, progressive hepatic and renal failure and lactic acidemia. Hum Mol Genet. 2015;24(8):2297–307.
West AP. Mitochondrial dysfunction as a trigger of innate immune responses and inflammation. Toxicology. 2017;391:54–63.
de Mello AH, Costa AB, Engel JDG, Rezin GT. Mitochondrial dysfunction in obesity. Life Sci. 2018;192:26–32.
Van der Voorn L, Ploegh HL. The wd-40 repeat. Febs Lett. 1992;307(2):131–4.
Smith TF, Gaitatzes C, Saxena K, Never EJ. The wd repeat: a common architecture for diverse functions. Trends Biochem Sci. 1999;24(5):181.
Zhu Y, Wang Y, Xia C, Li D, Li Y, Zeng W, et al. Wdr26: a novel gβ-like protein, suppresses mapk signaling pathway. J Cell Biochem. 2004;93(3):579–87.
Boissan M, Dabernat S, Peuchant E, Schlattner U, Lascu I, Lacombe ML. The mammalian nm23/ndpk family: from metastasis control to cilia movement. Mol Cell Biochem. 2009;329(1-2):51–62.
Amendola R, Martinez R, Negroni A, Tanno B, Calabretta B, Raschellà G. DR-nm23 gene expression in neuroblastoma cells: relationship to integrin expression, adhesion characteristics, and differentiation. J Natl Cancer Inst. 1997;89(17):1300–10.
Chitrakar I, Kim-Holzapfel DM, Zhou W, French JB. Higher order structures in purine and pyrimidine metabolism. J Struct Biol. 2017;197(3):354–64.
Zhang C, Wang K, Yang L, Liu R, Chu Y, Qin X, et al. Lipid metabolism in inflammation-related diseases. Analyst. 2018;143(19):4526–36.
Bennett M, Gilroy DW. Lipid Mediators in Inflammation. Microbiol Spectr. 2016;4(6):10.
Wallace JL. Prostaglandin biology in inflammatory bowel disease. Gastroenterol Clin North Am. 2001;30(4):971–80.
Mondanelli G, Iacono A, Carvalho A, Orabona C, Volpi C, Pallotta MT, et al. Amino acid metabolism as drug target in autoimmune diseases. Autoimmun Rev. 2019;18(4):334–48.
Mondanelli G, Bianchi R, Pallotta MT, Orabona C, Albini E, Iacono A, et al. A Relay Pathway between Arginine and Tryptophan Metabolism Confers Immunosuppressive Properties on Dendritic Cells. Immunity. 2017;46(2):233–44.
Banchereau J, Steinman RM. Dendritic cells and the control of immunity. Nature. 1998;392:245–52.
Idzko M, Ferrari D, Eltzschig HK. Nucleotide signalling during inflammation. Nature. 2014;509(7500):310–7.
Myrtek D, Idzko M. Chemotactic activity of extracellular nucleotideson human immune cells. Purinergic Signal. 2007;3(1–2):5–11.
Ferrari D, la Sala A, Panther E, Norgauer J, Di Virgilio F, Idzko M. Activation of human eosinophils via P2 receptors: novel findings and future perspectives. J Leukoc Biol. 2006;79(1):7–15.
Hayden MS, Ghosh S. NF-κB in immunobiology. Cell Res. 2011;21(2):223–44.
Thalhamer T, McGrath MA, Harnett MM. MAPKs and their relevance to arthritis and inflammation. Rheumatology (Oxford). 2008;47(4):409–14.
Banerjee S, Biehl A, Gadina M, Hasni S, Schwartz DM. JAK-STAT Signaling as a Target for Inflammatory and Autoimmune Diseases: Current and Future Prospects. Drugs. 2017;77(5):521–46.
Scholz J, Woolf CJ. The neuropathic pain triad: neurons, immune cells and glia. Nat Neurosci. 2007;10(11):1361–8.
Wolf G, Gabay E, Tal M, Yirmiya R, Shavit Y. Genetic impairment of interleukin-1 signaling attenuates neuropathic pain, autotomy, and spontaneous ectopic neuronal activity, following nerve injury in mice. Pain. 2006;120(3):315–24.
We thank Dr. Mei-chen Liu for her critical revision of the manuscript for important metabolic pathway analysis content.
Studies involving human participants, human data, or human tissue
The study was supported by National Natural Science Foundation of China (nos. 81373943, 81573583, 81891012, 81630101, and U19A2010).
Ethics approval and consent to participate
All experiments were performed consistent with animal use in science in European Communities Council Directive (2010/63/EU) and followed by Chengdu University of Traditional Chinese Medicine of Laboratory Animal Ethics Committee.
Consent for publication
The authors declare no conflicts of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Primers used for RT-qPCR.
Metabolic pathways associated to varied influences of FA and FB in CuSO4-induced zebrafish inflammation.
GO analysis of differentially-expressed proteins.
Protein-protein interactions of differentially expressed proteins.
Interactive network construction by integrated metabolomics and proteomics analysis.
About this article
Cite this article
Gong, L., Yu, L., Gong, X. et al. Exploration of anti-inflammatory mechanism of forsythiaside A and forsythiaside B in CuSO4-induced inflammation in zebrafish by metabolomic and proteomic analyses. J Neuroinflammation 17, 173 (2020). https://doi.org/10.1186/s12974-020-01855-9
- Forsythiaside A
- Forsythiaside B