Quality variation and biosynthesis of anti-inflammatory compounds for Capparis spinosa based on the metabolome and transcriptome analysis

Introduction Capparis spinosa L. fruits as edible and medicinal plant, has anti-inflammatory activities. The different morphological characteristics of C. spinosa fruits from Ili, Turpan, and Karamay may affect their anti-inflammatory components and functions. Methods The anti-inflammatory activity of C. spinosa fruit was assessed using an LPS-induced inflammatory cell model. Furthermore, the differences in anti-inflammatory compounds were analyzed by metabolome and RNA-seq. Additionally, the anti-inflammatory mechanism was elucidated using network pharmacology. Results In the study, we found that the 95% ethanol extracts (CSE) obtained from the three kinds of fruits showed remarkable anti-inflammatory effects both in vivo and in vitro. However, the CSE derived from Ili fruits significantly reduced CD86 levels on DCs. As a result of metabolomic analysis, the metabolic profiles of Ili fruits differed significantly from those of the other two habitats, which were consistent with transcriptome analysis. A total of 15 compounds exhibiting anti-inflammatory activity were subjected to screening, revealing a greater accumulation of flavonoids in the Turpan and Karamay districts. Notably, phenolic compounds were identified as the principal anti-inflammatory components in C. spinosa. Conclusion There were significant differences in the morphology, metabolites, transcriptional levels, and anti-inflammatory activity of C. spinosa from the three districts.

investigation, we found that more than 60% of the capers in Xinjiang came from Turpan. It has not been studied whether such morphological and habitat differences affect the antiinflammatory components and activity of C. spinosa.
The variation of habitat and morphology of traditional Chinese medicine (TCM) sometimes affects the content of active substances. The dry root and rhizome of Salvia miltiorrhiza Bunge (Danshen) are a TCM. Tanshinones are important active compounds, and their accumulation in the pericardium affects root color. The contents of Tanshinone IIA and Tanshinone I in the orange roots of the mutant were significantly reduced. However, the key enzyme genes involved in biosynthesis did not differ at the transcriptional level between the two kinds of S. miltiorrhiza (Zhan et al., 2019). Cistanche deserticola is an edible and medicinal plant, with phenylethanoid glycosides (PhGs) as its major active compounds. However, the content of PhGs in samples from three ecotypes grown in saline-alkali land, grassland, and sandy land was significantly different. Compared with other ecotypes, the content of PhGs was higher in saline-alkali soils, which may be due to the up-regulation of PhGs biosynthesis genes . Sixty samples of glycyrrhiza uralensis from different districts in Gansu, China, were analyzed, and the results showed significant differences in active compounds (Bai et al., 2020).
In this study, we conducted a comparative analysis of the antiinflammatory activity, metabolites, and transcriptome data of C. spinosa fruits from different habitats. We identified the habitats with the strongest anti-inflammatory activity and screened the major anti-inflammatory components in C. spinosa fruits. Through gene-metabolic network analysis, we identified several key metabolic pathways and genes that influence the synthesis of anti-inflammatory compounds. These findings provided new insights into the further development, cultivation, and genetic modification of C. spinosa.

Fruit morphology and sampling location information
C. spinosa fruits (6/habitat) were collected from three districts, including Turpan, Karamay, and Ili in Xinjiang Uygur Autonomous Region of China, which were named as groups A, B, and C. The fruits were identified by Dr. Jiang He from the Xinjiang Institute of Materia Medica. The fruits from Turpan were green and medium in size ( Figure 1A), the fruits from Karamay were the smallest and bottle-green or brown ( Figure 1B), and the fruits from Ili were the largest and brown ( Figure 1C). The GPS coordinates and climate information of the habitats can be found in Table 1. Turpan had the least precipitation, the most arid climate, and the highest annual average temperature, while Ili had the most humid climate and the lowest annual average temperature, and Karamay had the middle climate index. The contents of N, P, and K in Ili soil were the highest.

Validation of anti-inflammatory function
C. spinosa fruits were extracted with 95% ethanol, then the extract was concentrated and lyophilized, named CSE. The yield of the extract was about 24%. CSEs prepared with C. spinosa fruits from Turpan, Karamay, and Ili were named CSEA, CSEB, and CSEC, respectively. To detect the anti-inflammatory effects of CSE in vitro, different concentrations (0.5, 1.0 and, 1.5 mg/mL) of CSEs were used to treat DCs in the presence of LPS. Results of flow cytometry showed that there was no significant change in cell proportion among all groups, suggesting that the selected doses of CSEs did not affect the viability of DCs (Figure 2A). Compared with the LPS group, CSEs did not inhibit the expression of surface molecule CD40 induced by LPS, whereas CSEC (1.5 mg/mL) significantly suppressed the expression of CD86 induced by LPS ( Figure 2B). Three pro-inflammatory cytokines, including TNF-a, IL-6, and IL-12p40, were significantly decreased ( Figure 2C).
In vivo, CSE treatment significantly suppressed the expression of pro-inflammatory cytokines TNF-a and IL-6 and relieved the degree of ear edema in the TPA-induced mouse ear edema model ( Figure 2D).

FIGURE 1
The fruits from three districts. (A) Turpan; (B) Karamay; (C) Ili. These results showed that C. spinosa from all three different districts had anti-inflammatory effects, but the effect of fruits from Ili was better than the other two groups in inhibiting DC maturation.

Qualitative and quantitative metabolites
Base peak ion chromatograms of the three groups were shown as fingerprints. There were significant differences among the samples of Ili and the other two districts ( Figure 3A). A total of 82 compounds were identified, of which 47 were in positive ion mode (Tables S1, S3) and 43 were in negative ion mode (Tables S2,  S4). Eight compounds were found in both modes. Principal component analysis (PCA) was performed on the 18 samples (6/ habitat) ( Figure 3B) or together with quality control (QC) ( Figure  S1). All samples were within Hotelling's T-squared ellipse, indicating that the detection system had good stability. In positive ion mode, PC1, PC2, and PC3 were 66.95%, 19.14%, and 8.98%, respectively, while in negative ion mode, PC1, PC2, and PC3 were 87.59%, 7.16% and, 1.29%, respectively ( Figure 3B). The model of orthogonal partial least squares discriminant analysis (OPLS-DA) The anti-inflammatory activity of capers from different districts in vitro and in vivo. Therapeutic effect of CSEs on mouse edema model. Data are means ± SE (n = 5). ##P < 0.01, ###P < 0.001 compared to DMSO; *P < 0.05, **P < 0.01, ***P < 0.001 compared to LPS/Model. was stable and reliable ( Figure S2). The score plots of PCA and OPLS-DA exhibited an obvious separation among the samples.
The identified compounds were analyzed by hierarchical clustering analysis. The heatmap of hierarchical cluster analysis (HCA) showed that group C was separated from the other two groups ( Figure 3C). Differential metabolites were screened by the Venn analysis in positive and negative modes, respectively ( Figure 3D; Table S5).
Cosine similarity analysis showed a significant difference among the three groups in negative ion mode, and the similarity value in group C was 86% ( Figure 3E). The clustering tree revealed that the distance between groups A and B was close, but group C was far from the other two groups ( Figure 3F). Stachydrine, ferulic acid, and chrysoeriol were selected to be quantified, and their contents in groups A, B, and C were consistent with the metabolomics ( Figure S3). These results suggested that the metabolic level of group C was considerably different from that of the other two groups. The metabolomics data were deposited in the OMIX repository, accession number OMIX004470 (https://ngdc. cncb.ac.cn/omix/preview/QoEYJ1oa).

Screening and analyzing of antiinflammatory compounds
By retrieving the database, 15 anti-inflammatory compounds were selected. They included one alkaloid, four lipids, three flavones, four phenols, one glycoside, and two other compounds ( Figure 4A).
Considering the quantitative value ( Figure 4B) and log2 (fold change) value of the compound ( Figure 4C), the contents of almost compounds were similar between groups A and C except for oleoyl ethyl amid, 3,4-dihydroxybenzaldehyde, 2,3,5,6-tetramethylpyraz and vanillic acid, but the contents of almost compounds in group A and B were different with that in group C except for 6-gingerol and ferulic acid. The contents of three flavonoids including diosmetin, chrysoeriol, and astragalin in groups A and B were significantly higher than that in group C. The contents of linoleic acid, monoolein, and cappariloside A in group C were significantly higher than that in groups A and B, whereas the content of 2,3,5,6tetramethylpyrazine in group B was highest. The contents of stachydrine, 6-gingerol, ferulic acid, and b-sitosterol in the three different districts showed no significant difference in the three groups. The correlation analysis of the quantitative values of metabolites showed that three flavones had a significantly negative correlation with stachydrine, linoleic acid, monoolein, cappariloside A, and b-sitosterol. There was a significantly positive correlation between the three flavones and linoleoyl ethanolamide. Linoleoyl ethanolamide was negatively correlated with stachydrine, linoleic acid, monoolein, and cappariloside A. Stachydrine had a significantly high positive correlation with linoleic acid, monoolein, and cappariloside A. Stachydrine was significantly correlated with 3,4-dihydroxybenzaldehyde, ferulic acid and bsitosterol. Among the four phenol compounds, ferulic acid, vanillic acid and, 3,4-dihydroxybenzaldehyde had a significant correlation, but the correlation among 6-gingerol and other compounds was weak. The contents of compounds in similar synthetic pathways existed a positive correlation ( Figure 4D).
The contents of most anti-inflammatory compounds were strongly correlated with environmental factors. High evaporation, high temperature, and high Olsen-P contributed to the accumulation of the three flavones compounds. The trend of linoleoyl ethanolamide was similar to that of flavones. The trends of stachydrine, linoleic acid, monoolein, 6-gingerol, ferulic acid, cappariloside A, and b-sitosterol were opposite to that of flavone. On the other hand, 2,3,5,6-tetramethylpyrazine had a weak correlation with environmental factors ( Figure 4E).
There was no difference in abscisic acid and tryptophan contents among the three groups, indicating that all the materials were in the same period of fruit development ( Figure S4). Compounds with consistent response trends to environmental factors also had a strong positive correlation in their content, indicating that the quality variation of C. spinosa was caused mainly by environmental factors.

Prediction of anti-inflammatory mechanisms by network pharmacology
To explore the anti-inflammatory mechanism of C. spinosa fruit, the above 15 anti-inflammatory compounds were used for network pharmacological analysis. 564 targets of compounds were collected, and 1,196 human targets of rheumatoid arthritis (RA) were obtained. Thus, after intersection analysis, 167 common targets were acquired ( Figure S5). The common targets were imputed into STRING to build the protein-protein interaction (PPI) network. Subsequently, the interaction network was constructed for illustrating the interactive relationship between compounds and targets ( Figure 5A). Based on degree value, the top 20 were the core targets including MAPK1, STAT3, IL-6, TNF, VEGFA (vascular endothelial growth factor A), SRC, MMP9, and JUN ( Figure 5B).
In addition, the common targets were subjected to GO and KEGG enrichment analysis, and the terms associated with RA of GO mainly focused on the regulation of interleukin-6 production, regulation of tumor necrosis factor production, regulation of interleukin-12 production, dendrite, immune receptor activity, cytokine activity, etc. ( Figure 5C). KEGG pathways mainly focused on the TNF signaling pathway, NF-kappa B signaling pathway, Cytokine-cytokine receptor interaction, etc. ( Figure 5D).
These results suggested that C. spinosa plays an anti-inflammatory role by affecting cytokine expression and lymphocyte differentiation.

Screening of anti-inflammatory compounds in vitro
To identify key active compounds in CSE that inhibit the maturation of LPS-induced DC, we examined the effect of stachydrine (ST), diosmetin (DIO), and 3,4-dihydroxybenzaldehyde (DHB) to treat DCs in the presence of LPS. Sachydrine and diosmetin did not inhibit the expression of costimulatory molecules (CD40 and CD86) and the secretion of pro-inflammatory cytokines IL-6 and TNF-a ( Figures 6A, B). 3,4-dihydroxybenzaldehyde significantly suppressed the expression of CD40 and IL-6 induced by LPS at the concentration of 1 mM, and it significantly suppressed the expression of CD86 and TNF-a at the concentration of 0.5 mM ( Figure 6C). The results showed that phenols were the main compounds that inhibited DCs maturation and secretion of pro-inflammatory cytokines in capers.

Identifying key genes in phenols biosynthetic pathway
Transcriptome data of C. spinosa fruits were analyzed. PCA showed PC1, PC2, and PC3 were 49.31%, 17.50%, and 9.04%, respectively, and all samples were within the Hotelling's Tsquared ellipse. Group C was significantly separated from the other two groups ( Figure S6A). Volcano Analysis showed that 2,432 genes were up-regulated, and 1,705 genes were downregulated in group A and group B, 9,104 genes were up-regulated, and 4,209 genes were down-regulated in group A and group C, 4,872 genes were up-regulated, and 2.320 genes were downregulated in group B and group C ( Figure S6B). HCA and Venn's analysis showed that group C was significantly different from the other two groups (Figures S6C, D). The enrichment differential pathways were concentrated in biosynthesis and metabolism pathways between groups A and C ( Figure S7). These data indicated that group C was significantly different from the other two groups at the transcriptional level. These transcript data were deposited into the NCBI database with identifier number PRJNA778809 (https://www.ncbi.nlm.nih.gov/sra/PRJNA778809).  Figure 7B). The qRT−PCR results showed good consistency with the RNA-Seq data. The qRT−PCR results indicated that the expression of CYP73A (Cluster-22225.2180) in groups A and B was significantly higher than that in group C, while the expression of CYP73A (Cluster-22225.4290) in group C was significantly higher than that in the other two groups. REF1 (Cluster-22225.7764) was highly expressed in all three groups ( Figure S8).
The expression level of REF1 (Cluster-22225.7764) indicated a strong positive correlation with ferulic acid content (r = 0.7667, P = 0.0214), which might be the key gene promoting ferulic acid biosynthesis in C. spinosa ( Figure 7C). These results indicated that C. spinosa fruits from different districts affected the biosynthesis of secondary metabolites by regulating gene transcription.

Materials and regents
The fruits of C. spinosa were collected from June to July 2019 from Turpan, Karamay, and Ili in Xinjiang Uygur Autonomous Region of China. The fruits were immediately frozen in liquid nitrogen and then stored at -80°C. The identification was carried out by Dr. Jiang He of Xinjiang Institute of Materia Medica.
BALB/c mice (20 ± 5 g) aged six weeks old were obtained from the experimental animal center of the Xinjiang Medical University (Urumqi, China). Mice were allowed free access to distilled water and standard food, mice room under conditions of maintained temperature (24~26°C) and humidity (69~71%) as well as a 12:12 h light-dark cycle. All experimental animal protocols were reviewed and approved by the Ethics Committee of Xinjiang University for the use of Laboratory Animals.

Preparation of C. spinosa fruits ethanol extracts
First, the dried fruits of C. spinosa were ground into a powder and passed through a 40-mesh sieve. Then, 50 g of the powder was mixed with 95% ethanol in a ratio of 1:6 (w/v). The mixture was ultra-sonicated for 20 minutes, followed by extraction for 2 hours in a water bath at 60°C . This extraction process was repeated three times, and the extracts were mixed. The extract was concentrated using rotary evaporation. Then, the extract was subjected to lyophilization and named 'CSE'.

Detection of anti-inflammatory effect in vitro
Immature DCs were induced from bone marrow cells of BALB/ c mice by Granulocyte-macrophage colony-stimulating factor (GM-CSF, Peprotech, USA) (Aipire et al., 2017). The cells were cultured in RPMI-1640 medium containing 10% heat-inactivated fetal bovine serum (FBS), 100 units/ml penicillin-streptomycin, 50mM b-mercaptoethanol and 20 ng/ml GM-CSF. Subsequently, they were placed in a 5% CO 2 , 37°C incubator. The cells were collected on the B C A FIGURE 6 Anti-inflammatory effects of three compounds in vitro. Data are means ± SE (n = 3). ## P < 0.01, ### P < 0.001 compared to Untreated; *P < 0.05, **P < 0.01, ***P < 0.001 compared to LPS. seventh day, cultured at a concentration of 1 × 10 6 /mL on a 24-well plate, and treated with different concentrations of CSEs or standards in the presence of 80 ng/mL LPS for 12h. The cells were collected and stained using the mAbs, including anti-mouse CD40, and anti-mouse APC-CD86. All samples were detected on FACSCalibur (BD Biosciences, USA), and the data were analyzed by the FlowJo software (v7.6). The cell supernatant was collected, and the cytokines (TNF-a, IL-6, and IL-12p40) were measured by ELISA kit. OD 450 was measured using a microplate reader (Bio-Rad, USA).

Treatment of TPA-induced mouse ear edema model
TPA was used to induce the ear edema model in mice. The 45 mice were randomly divided into a control group (DMSO), model group (TPA, 75 mg/kg), positive control group (dexamethasone, DXM, 2.5 mg/kg), low-dose group (200 mg/kg) and high-dose group (500 mg/kg), with five mice in each group. The medicine was applied to the right ear for 6 h, ear edema degree of mice was calculated by right ear mass (mg)/left ear mass (mg), and the cytokines (TNF-a and IL-6) were measured by ELISA kit.

Differential metabolite analysis
The raw data were imported into Progenesis QI software for retention time correction, peak identification, peak extraction, peak integration, and peak alignment, and the corresponding TCM metabolic database was established. Substance identification of the unknown components was performed using the self-built mass spectrometry database. Multivariate statistical analysis, including PCA, and OPLS-DA, was performed using SIMCA software (V15.0.2, Sartorius Stedim Data Analytics AB, Umea, Sweden). Also, R software or GraphPad Prism software (V8.0, La Jolla, CA, USA) was used HCA, Venn analysis, cosine similarity analysis, radar plot, clustering tree, pie chart, histogram, and Pearson correlation analysis were performed in combination with univariate analysis for metabolite accumulation patterns between different samples.

Validation of metabolomic data by UHPLC-QTOF-MS
Metabolites were extracted according to 2.2. Stachydrine, ferulic acid, and chrysinol (HPLC ≥ 98%, Yuanye, China) were dissolved in with methanolic solution, respectively. The prepared extracts were each taken 400 mL, mixed to prepare mixed standard solutions and a series of mixed standard solutions of different concentrations were prepared by stepwise dilution with methanol. The mass spectrometer apparatus and elution conditions were the same as in 3.4.

Analysis of network pharmacology
The compounds with anti-inflammatory functions were screened from metabolites for network pharmacological analysis. The chemical targets were screened by TCMSP databases (http:// lsp.nwu.edu.cn/tcmsp.php), Swiss Target Prediction Database (http://www.swisstargetprediction.ch/), SEA databases (https:// sea.bkslab.org/). The disease targets were screened by GeneCards (https://www.genecards.org), and OMIM databases (https:// www.omim.org/). Furthermore, the keyword 'rheumatoid arthritis' was used to retrieve all targets under the condition of 'Homo sapiens'. Venny software (v2.1.0) was used to intersect antiinflammatory component targets and disease targets, which were named 'the common targets'. Then, the common targets of RA and active ingredients were used to construct the PPI network by the STRING database (https://string-db.org/cgi/input.pl). The results were visualized using Cytoscape software (v3.6.1). The top 20 targets with degree value in the PPI network were named 'core targets'. The R Project was used to visualize core targets. GO enrichment analysis of common targets was performed by Metascape (http://metascape.org).

Extraction RNA and analysis of transcriptomics
Total RNA was extracted using TRIzol reagent, and genomic DNA was removed using DNase I (Takara) according to the manufacturer's instructions. RNA quality and quantification were measured using a 2100 Bioanalyser (Agilent) and an ND-2000 (Nano-Drop Technologies) RNA-seq transcriptome libraries were then prepared from total RNA (1 mg) using the TruSeqTM RNA Sample Preparation Kit (Illumina, San Diego, CA, USA). The expression levels of each transcript were calculated according to the mapping reads per million exon kilobase fragment method to identify these DEGs in different samples. EdgeR (Empirical Analysis of Digital Gene Expression in R) was used for differential expression analysis and NR (NCBI non-redundant protein sequences) was used to annotate the assembled unigenes. In addition, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analyses were performed for these DEGs using Goatools (https:// github.com/tanghaibao/Goatools) and KOBAS (http:// kobas.cbi.pku.edu.cn/home.do), respectively.

Correlation analysis mining antiinflammatory compounds and genes
Based on the KEGG pathway analysis (https://www.kegg.jp/ kegg/pathway.html), genes related to the synthesis of antiinflammatory compounds were screened. The key genes were identified by Spearman correlation analysis using the GraphPad Prism 8 software.

Validation of RNA−seq data by qRT−PCR
Total RNA was reverse transcribed using the reverse transcription kit. Then, qPCR amplification was performed on the qTOWER3G qPCR instrument. The primer sequences are shown in Table S1.

Statistical analysis
Statistical significance analysis was performed using one-way ANOVA or t-test and computed using GraphPad Prism software. A value of P < 0.05 was statistically significant.

Discussion
In this study, we compared the fruit quality differences of capers from different districts in Xinjiang through metabolome and transcriptome analysis and detected their anti-inflammatory activity in vitro and in vivo. The fruits from Ili were significantly different from Turpan and Karamay in metabolism and transcription levels. In addition, fruits from all three different districts inhibited the expression of pro-inflammatory cytokines, but ethanol extract of caper from Ili significantly inhibited the expression of CD86 on DCs induced by LPS. There were significant differences in the different districts, morphology, metabolites, gene expression level, and antiinflammatory activity of capers from the three districts.
C. spinosa fruit extract rich in flavonoids, indoles, and phenolic acids, effectively inhibited carrageenan-induced paw edema in mice (Zhou et al., 2010). C. spinosa fruit ethanol extracts significantly inhibited the secretions of IL-12p40, IL-6, IL-1b, and TNF-a in LPSinduced DC (Azeguli et al., 2017). C. spinosa leaf alcoholic extract significantly decreased immune cell infiltration, vasodilatation, and dermis thickness in the inflammatory site by inhibiting cytokine gene expression, including IFN-g, IL-17, and IL-4 in the contact hypersensitivity mice. The extract contained saponins, flavonoids, and alkaloids (Azhary et al., 2017). 70% alcohol extract of C. spinosa aerial parts significantly reduced the levels of TNF-a, COX-2, IL-1b, and IL-6 in the LPS-induced inflammation in microglia (Rahimi et al., 2020). Pro-inflammatory cytokines IL-6, TNF-a, and IL-1 played an important pathogenic role in inflammatory processes such as RA (Kaneshiro et al., 2019). The results of this study also confirmed that C. spinosa extract could treat inflammation by inhibiting the expression of pro-inflammatory cytokines both in vitro and in vivo. Thus, fruits from all three districts can be used in the treatment of inflammatory diseases.
At present, most studies on the anti-inflammatory effects and mechanisms of C. spinosa fruit are based on experimental animal models. However, it remains uncertain whether comparable positive effects can be achieved in the treatment of human diseases. We will further verify the anti-inflammatory effect of C. spinosa fruit on human cells. These results will lay the foundation for clinical trials of the fruit for the treatment of inflammatory diseases.
However, different geographical locations greatly influence the accumulation of active compounds and their function in C. spinosa. In this study, we found that the content of secondary metabolites of C. spinosa fruits from three different districts was significantly different. Consistently, the analysis of different samples from various regions of Sardinia (Italy) revealed qualitative and quantitative differences in the content of flavonoids, glucosinolates, anthocyanins, and phenolic acids (Maldini et al., 2016). The C. spinosa bud samples of Morocco, Turkey, and Italy were analyzed, and the results showed that the content of phenol in the Morocco sample was the highest, and that of rutin in the Italian decoction was the highest. Significantly, the highest acetylcholinesterase inhibitory activity was observed in Turkish and Moroccan samples, and the best butyrylcholinesterase inhibitory effects were detected in the Italian samples (Stefanucci et al., 2018). Therefore, the chemical composition of C. spinosa with different geographical distributions changed significantly, which endowed them with various biological activities. Among the antiinflammatory compounds, the flavonoids in fruits from Turpan were the highest, while those in fruits from Ili were the lowest. The biosynthesis and accumulation of flavonoids were related to drought stress (Nakabayashi et al., 2014;Yang et al., 2020). The humid climate in the Ili region was not conducive to the accumulation of flavonoids.
The extracts from all three districts reduced TNF-a and IL-6 both in vitro and in vivo. However, only a few authors reported the inhibition of TNF-a and IL-6 by flavonoids (Peluso et al., 2015). DCs are professional antigen-presenting cells, link innate and adaptive immune responses. They regulate T cell homeostasis and inflammatory response by expressing costimulatory molecules and releasing cytokines (Waisman et al., 2017). However, ethanol extract of caper from Ili had the best inhibitory effect on the expression of DC surface molecule CD86. Therefore, flavonoids may not be the main component of C. spinosa fruits in inhibiting DC maturation. Our experimental results showed that diosmetin had no anti-inflammatory activity at the concentration of 4 mM, which supported the above conclusions.
The contents of stachydrine were higher in 15 antiinflammatory compounds, and there was little difference in the contents among the three districts. But stachydrine also had no anti-inflammatory activity. Therefore, the anti-inflammatory components in capers need to be further verified by subsequent experiments. Ferulic acid reduced TNF-a, IL-1b, IL-6, IL-23, and MMP9, and increased TGF-b in inflammatory conditions (Doss et al., 2018;Ganesan and Rasool, 2019;Yin et al., 2019;Zhu et al., 2020). On the other hand, 6-gingerol inhibited COX-2, IL-17, IL-6, IL-1b, IL-8, and TNF-a, and increased IL-10 to treat inflammation (Kim et al., 2005;Saha et al., 2016;Zhang and Zheng, 2018;Liu et al., 2020;Sheng et al., 2020). 3,4-Dihydroxybenzaldehyde was a phenolic compound with high content in capers. Our results showed that it significantly inhibited DCs maturation and secretion of pro-inflammatory cytokines. These studies suggested that phenols had good antiinflammatory activity.
Plants respond to environmental changes by regulating gene expression and thus affect metabolic processes. Cistanche deserticola was investigated and the 12 key genes involved in PhGs biosynthesis were found to be differentially expressed among three ecotypes . The content of starting compound phenylalanine was the highest in group C. However, due to the down-regulated expression of most genes in the synthetic pathway, there was no difference in the content of the final product 6-gingerol among the three groups. As an intermediate product, ferulic acid did not accumulate in the three groups. According to the content of compounds, the biosynthetic pathway of 3,4-dihydroxybenzaldehyde played a dominant role.

Conclusion
In this study, we found there were significant differences in the morphology, metabolites, gene expression level, and antiinflammatory activity of capers from the three different districts. All the 95% ethanol extracts from the three districts showed remarkable anti-inflammatory effects, but the anti-inflammatory effects of Capparis spinosa L. vary by origin. The extract from Ili significantly inhibited the expression levels of CD86 on DCs induced by LPS. Phenols were the main compounds that inhibited DCs maturation and secretion of pro-inflammatory cytokines in capers. The gene of REF1 (Cluster-22225.7764) might be the key gene promoting ferulic acid biosynthesis in C. spinosa.

Data availability statement
The metabolomics data presented in the study are deposited in the OMIX repository, accession number OMIX004470. The transcript data presented in the study are deposited in the NCBI repository, accession number PRJNA778809.

Ethics statement
The animal study was reviewed and approved by Ethics Committee of Xinjiang University, Xinjiang University (XJUAE-2019-017).

Author contributions
XL, AA, XW, and JL designed the research and analyzed the data. XL, AA, YD, PH, and LM performed the experiments. JL, XW, and XL wrote the manuscript. All authors contributed to the article and approved the submitted version.