Naringin Exerts Therapeutic Effects on Mice Colitis: A Study Based on Transcriptomics Combined With Functional Experiments

Naringin has been shown to exert protective effects in an animal model of ulcerative colitis, but detailed mechanisms remain unclear. This study aimed to investigate function and signaling mechanisms underlying naringin-induced therapeutic effects on colitis. Two mouse models were established to mimic human Inflammatory bowel disease (IBD) by treating drinking water with dextran sodium sulphate or intra-colonic administration of 2, 4, 6-trinitrobenzene sulfonic acid. Transcriptomics combined with functional experiments were used to investigate underlying mechanisms. Colitis symptoms, including weight loss and high disease activity index were significantly reversed by naringin. The inflammatory response, oxidative reactions, and epithelial cell apoptosis that occur with colitis were also alleviated by naringin. After naringin treatment, transcriptomics results identified 753 differentially expressed mRNAs that were enriched in signaling pathways, including the neuroactive ligand-receptor interaction, calcium signaling, and peroxisome proliferator-activated receptor (PPAR) signaling. The naringin-induced alleviation of colitis was significantly inhibited by the PPAR-γ inhibitor BADGE. In IEC-6 and RAW264.7 cells incubated with lipopolysaccharide (LPS), NF-κB-p65, a downstream protein of PPAR-γ, was significantly increased. Naringin suppressed LPS-induced high expression of NF-κB-p65, which was inhibited by small interfering RNA targeting PPAR-γ. Our study clarifies detailed mechanisms underlying naringin-induced therapeutic effects on mice colitis, and PPAR-γ was found to be the main target of naringin by functional experiments both in vivo and in vitro. Our study supplies new scientific information for the use of naringin in colitis treatment.


INTRODUCTION
Inflammatory bowel disease (IBD) is an array of chronic inflammatory disorders within the gastrointestinal tract and includes primarily ulcerative colitis (UC) and Crohn's disease (CD). CD affects the entire gastrointestinal tract from mouth to anus, whereas UC mainly affects colon (Sairenji et al., 2017). Since 2000, incidence ranges of CD are 6-11/100,000 and that of UC are 6-15/100,000 in the Western world (Kaplan and Windsor, 2021). The main symptoms of IBD include diarrhea, abdominal pain, mucus stool and bloody stool, which can be accompanied by a variety of intestinal and parenteral complications, such as intestinal obstruction, stenosis, fibrosis, and joint lesions. Moderate and severe UC or colonic CD is a high-risk factor for colon cancer (Grivennikov, 2013). IBD has become one of the most complex and refractory intestinal diseases which bring a heavy economic and social burden to the health systems of various countries all over the world.
The IBD pathogenesis is still under investigation, recent studies have suggested that both genetic and environmental factors are involved. Current drugs include sulfasalazine, 5aminosalicylic acid, broad spectrum antibiotics, and corticosteroids are often applied in the treatment of IBD (Abraham et al., 2017). However, many side effects like nausea, anorexia, cytopenia, myalgia, and malfunctions of the kidney, liver, and lungs have been reported (Abraham et al., 2017). Although biopharmaceuticals have been shown to exert good therapeutic effects on IBD, a large part of patients do not respond or lose response to biopharmaceuticals over time (Yanai and Hanauer, 2011). In addition, biologic therapy may increase the risk of serious infections and malignancies (Holmer and Singh, 2019). Therefore, it is valuable to screen or develop new drugs to expand therapy options.
In recent years, herbal medicines with multi-target actions have been widely used in the treatment of complex inflammatory diseases. Citrus peel is a kind of herbal medicine and is used as an antiinflammatory drug in China. Naringin (4′,5,7-trihydroxyflavanone-7-rhamnoglucoside) is a major and active flavanone glycoside isolated from citrus fruit species (Ali and El Kader, 2004). Naringin has been found to exert anti-inflammatory, anti-oxidative, and anti-cancer effects (Jeon et al., 2001;Zhang et al., 2016). Cao et al. find that naringin exerts protective effects on dextran sodium sulphate (DSS)induced colitis via modulation of PPAR-γ activity (Cao et al., 2018). However, the detailed mechanisms of naringin on colitis are still not fully understood because multiple targets may be involved. In addition, naringin-induced modulation of PPAR-γ should also be confirmed in cells and other animal models of colitis.
In this study, transcriptomics and functional experiments are used to uncover the detailed mechanisms of naringin-induced therapeutic effects on colitis. RNA sequencing (RNA-seq) is a deep-sequencing approach in transcriptome profiling that provides an impartial and accurate method for measuring the levels of transcripts and their isoforms (Kukurba and Montgomery, 2015). Two mouse models are established to mimic human IBD by treating drinking water with DSS or intra-colonic administration of 2, 4, 6-trinitrobenzene sulfonic acid (TNBS).

Animals
Eighty male C57BL/6 mice (6-8 weeks old, weighing 18-20 g) were obtained from the laboratory animal center, Dalian Medical University, Dalian city, China [Certificate of Conformity: No. SYXK (Liao) 2018-0007]. The experimental protocol was approved by Dalian Medical University Animal Care and Ethics Committee (No. AEE20046). The animals were acclimatized to laboratory conditions (23°C, 12 h/12 h light/dark, 50% humidity, ad libitum access to food, and water) for 2 weeks prior to experiments. The animal protocol was designed to minimize pain and discomfort to the animals. Mice were housed one per cage and were deprived of food for 12 h before experiments. All mice were euthanized by barbiturate overdose (intravenous injection, 150 mg/kg pentobarbital sodium) for intestinal tissue collection. All animal experiments were carried out in accordance with the National Institutes of Health guide for the care and use of laboratory animals. Animal studies are reported in compliance with the ARRIVE guidelines (Kilkenny et al., 2010). All applicable institutional and/or national guidelines for the care and use of animals were followed.

Experimental Design
For the DSS-colitis model, mice were divided randomly into five groups of six mice each. The mice were treated as follows: group I was a sham-operated control with gavage administration of saline; group II was untreated colitis; and groups III, IV, and V were treated with the following: SASP (500 mg/kg body weight, gavage administration, dissolved in saline), low-dose naringin (20 mg/kg body weight, gavage administration, dissolved in saline), and high-dose naringin (40 mg/kg body weight, gavage administration, dissolved in saline), respectively. All treatments occurred 1 day after colitis induction with DSS. SASP is an antiinflammatory drug that is, widely used for the clinical treatment of diseases such as IBD and, therefore, it was used as a positive control for the effects of naringin on colitis. Mice in group III, IV, and V were administered SASP or naringin by gavage once a day for seven successive days. Mice in groups II-V were given drinking water containing DSS (4% w/v) dissolved in autoclaved distilled water, to induce colitis symptoms (Feng et al., 2020). Control group I mice were given autoclaved distilled water and otherwise, treated in the same way. No abnormal moribund mice were found during the study. Food intake and bodyweight of mice were recorded once a day. On the 8th day, mice in group I-V were given autoclaved distilled water as drinking water. On the 9th day, distal colon samples in groups I-V were harvested for biochemical studies. Mice in the PPAR-γ inhibition group were injected intraperitoneally with 30 mg/kg BADGE solution per day for seven consecutive days, the experimental designs are the same as above.
For the TNBS-colitis model study, mice were divided randomly into five groups of six mice each. The mice were treated as follows: group I was a sham-operated control with gavage administration of saline; group II was untreated colitis; and groups III, IV, and V were treated with the following: SASP (500 mg/kg body weight, gavage administration, dissolved in saline), low-dose naringin (20 mg/kg body weight, gavage administration, dissolved in saline), and high-dose naringin (40 mg/kg body weight, gavage administration, dissolved in saline), respectively. All treatments occurred 1 day after colitis induction. Mice in group III, IV, and V were administered SASP (Xiong et al., 2017), naringin by gavage once a day for seven successive days. Colitis was induced with TNBS, as described previously (Dong et al., 2020). A catheter was inserted through the anus to approximately the level of the splenic flexure under urethane anesthesia. The colon was then infused with 0.1 ml of TNBS dissolved in ethanol (50% v/v) at a dose of 125 mg/kg. The mice could eat and drink ad libitum starting 1 h after the operation. On the 8th day, mice in groups I-V were given autoclaved distilled water as drinking water. On the 9th day, distal colon samples in groups I-V were harvested for biochemical studies. Protein extraction, western blotting, and ELISA were performed as previously described (Dong et al., 2020).

Assessment of Disease Activity Index
Mice body weights were monitored daily. The DAI scores were assessed according to our previous study (Melgar et al., 2005). Briefly, an observer blinded to the treatments combined the scores for weight loss, stool bleeding, and stool consistency on the last day.

Isolation of Colonic Epithelial Cells
Colonic epithelial cells were isolated according to a previous study (Dong et al., 2020). Briefly, distal colons were isolated from euthanized mice and immediately rinsed with ice-cold phosphate-buffered saline (PBS) to clear luminal contents. The distal colon was then opened longitudinally along the mesenteric border. Tissue was cut into approximately 2 mm long pieces and submerged in 40 ml of ice-cold PBS with 5 mM EDTA in a 50 ml Falcon tube. The pieces of tissue in PBS-EDTA were then incubated at 37°C with gentle rocking for 30 min. Following incubation, colonic tissue was vigorously shaken to disperse colonic crypts and surface epithelium in solution. Supernatant was then loaded into 1.5 ml microcentrifuge tubes and spun at 1 × 10 3 g for 5 min to pellet suspended cells. The isolated epithelial cells were used in downstream applications.

RNA-Seq and Bioinformatics
Three naringin treated DSS-induced C57BL/6 mice, three DSSinduced C57BL/6 mice and three untreated mice were used for mRNA profile detection. Total RNA was extracted using TRIzol reagent (Invitrogen), and quality was evaluated using the RNA Nano 6000 Assay Kit of the Bioanalyzer 2,100 system (Agilent Technologies, CA, United States).
A total of 3 µg RNA per sample was used to generate an mRNA sequencing library using NEBNext ® Ultra ™ RNA Library Prep Kit for Illumina ® (NEB, United States) following manufacturer's recommendations and index codes were added to each sample for identification (Gao et al., 2012). Briefly, mRNA was purified and fragmented. Then, first and second strand cDNA was synthesized, ends were repaired and adenylated, adapters were ligated, and fragments were enriched with PCR amplification. Library quality was assessed on the Agilent Bioanalyzer 2,100 system. Libraries were sequenced on an Illumina platform. Raw data (raw reads) of fastqc format were first processed through Trimmomatic (Bolger et al., 2014). In this step, clean data (clean reads) were obtained by removing reads containing adapter, reads containing poly-N, and low-quality reads from raw data. Reference genome and gene model annotation files were downloaded from the genome website directly. Index of the reference genome was built and paired-end clean reads were aligned to the reference genome using Hisat2 (v2.0.5), the selected mapping tool. The mapped reads of each sample were assembled by StringTie (v1.3.3b) in a reference-based approach. FeatureCounts v1.5.0-p3 was used to count the reads mapped to each gene. Then, FPKM of each gene was calculated based on the length of the gene and read counts mapped to this gene (Conesa et al., 2016). FPKM, expected number of Fragments Per Kilobase of transcript sequence per Million base pairs sequenced, considers the effect of sequencing depth and gene length for the read counts at the same time. Differential expression analysis of two conditions/groups (two biological replicates per condition) was performed using the DESeq2 R package (1.16.1). The resulting p-values were adjusted using the Benjamini and Hochberg approach. Genes with an adjusted p-value < 0.05 found by DESeq2 were assigned as differentially expressed. We used the clusterProfiler R package to test the statistical enrichment of differential expression genes in KEGG pathways. The raw sequencing data sets were deposited in the Sequence Read Archive of the NCBI (https://www.ncbi.nlm. nih. gov/sra) under accession number PRJNA741857.

Quantitative Real-Time Polymerase Chain Reaction
Validation of differentially expressed genes was performed by qRT-PCR using SuperReal PreMix Plus (SYBR Green, FP205, Tiangen) with an Applied Biosystems StepOnePlus Real-Time PCR System, according to the manufacturer's instructions. Total RNA was isolated from colonic epithelial cells using RNAiso Plus (9108, TaKaRa) and reverse-transcribed to single-stranded cDNA using a reverse transcription system (KR116, Tiangen). Primer sequences are listed in Table 1. A separate qRT-PCR using a primer for the detection of GAPDH was used as a control.

Statistical Analysis
The animal experiments, in vitro experiments, and data analysis were conducted according to a single-blind study design. Data were compared between three or more groups using one-way ANOVA, and between two groups using Student's t-test. Data were expressed as the mean ± standard deviation. Data were normally distributed, and each group showed similar variances. Further evaluation was carried out using Kruskal-Wallis rank sum tests. All experiments were repeated at least three times and a p value < 0.05 was considered statistically significant.

Effects of Naringin on Normal Mice and Cells
The chemical structure of naringin is shown in Figure 1A. Naringin had no significant effect on body weight or food intake in normal mice ( Figures 1B,C). Cytokine profiles [TNF-α, interleukins (IL)-1β, and interferon (INF)-γ] in colon tissue were not affected by naringin (Figures 1D-F). There was no significant effect on expression of PPAR-γ in colon tissue after gavage administration of 10 and 40 mg/kg naringin for 7 days, while the expression of PPAR-γ was decreased by 160 mg/kg naringin ( Figure 1G). Naringin (5-160 μM, 12 h) did not significantly affect cell viability in RAW264.7 cells. Naringin (5-80 μM, 12 h) exerted no significant effects on cell viability in IEC6 cells, while 160 μM significantly inhibited cell viability ( Figure 1H). According to our preliminary experiments and previous reports (Ahmad et al., 2015;Cao et al., 2018), the maximum doses of naringin were 40 mg/kg in vivo and 20 μM in vitro.

Naringin Ameliorates DSS-Colitis Symptoms
The significant weight loss, high colon weight-to-length ratio, decreased food intake and high DAI scores in the untreated model group showed that instillation of DSS to the colon resulted in reproducible colitis in mice (Figures 2A-D). The overall morphology of the colon is shown in Figure 2E.
Histopathological damage was characteristically induced by DSS ( Figures 2F,G). A variety of colitis symptoms, including damaged villi structure (green ring), infiltration of inflammatory cells (black arrow), submucosal edema (blue ring), and muscle fiber separation (blue arrow) were visible in the H and E staining ( Figure 2G). Naringin and SASP reversed the pathological changes in colitis.

Effects of Naringin on Inflammation and Apoptosis in DSS-Colitis
As shown in Figures 3A-D, MPO activity and expression of proinflammatory factors (TNF-α, IL-1β, and INF-γ) were significantly increased in DSS-colitis, and naringin significantly reversed these changes. DSS-induced colitis up-regulated inflammation-related proteins (iNOS and p-NF-κB-p65) and an apoptosis-related protein cl-caspase3 expression. Finally, an attenuated trend was observed in groups treated with naringin ( Figures 3E-G). These results confirmed that naringin has antiinflammatory and anti-apoptotic activities in DSS-colitis.

High-Throughput Transcriptome Sequencing and KEGG Pathway Enrichment Analysis
We applied transcriptome sequencing to investigate the possible mechanisms of naringin in the treatment of colitis. A total of 753 differentially expressed mRNAs before and after naringin Frontiers in Pharmacology | www.frontiersin.org August 2021 | Volume 12 | Article 729414 treatment were identified with fold changes > 2 and p < 0.05 and some representative genes were shown in heatmap ( Figures  4A,B). In addition, compared with the control group, 834 differentially expressed mRNAs were observed in the model group ( Figure 5A). The KEGG enrichment results showed the top 20 enriched KEGG pathways in the model group vs. naringin treatment group ( Figure 4C) and the control group vs. model group ( Figure 5B). There are three common enriched KEGG pathways including the neuroactive ligand-receptor interaction, calcium signaling, and PPAR signaling among the top 20 enriched KEGG pathways between the model group vs. naringin treatment group and the control group vs. model group ( Figure 4C, Figure 5B). Based on the genes contained in the top 11 pathways, we constructed a KEGG-target network of naringin-induced treatment including 133 targets ( Figure 4D). Most targets were associated with different pathways, suggesting synergistic effects are involved in naringin-induced treatment. Therefore, naringin may target different pathways and targets.

Effects of Naringin on PPAR-γ
According to reference validation, we subsequently studied proteins in the PPAR signaling pathways to uncover naringininduced treatment effects on colitis. We used qRT-PCR to detect the roles of some mRNAs in PPAR signaling pathways, including PPAR-α, PPAR-γ, PPAR-δ, Scd1, and Hmgcs2 in naringin- induced treatment effects ( Figure 6). There were significant differences in the expressions of PPAR-α and PPAR-γ between the naringin treatment and DSS-colitis groups. The expressions of PPAR-δ, Scd1, and Hmgcs2 were not significantly affected by naringin in the colitis model. We further studied the protein expressions of PPAR-α and PPAR-γ. Compared with the sham group, PPAR-γ expression was significantly decreased in DSS-colitis, but naringin treatment significantly stimulated PPAR-γ expression ( Figures 7B-D). However, naringin treatment did not significantly affect PPAR-α protein expression ( Figure 7A). The results suggest that PPAR-α and PPAR-γ may be involved in naringininduced treatment effects, however, PPAR-γ plays a more important role.
Next, we examined the role of PPAR-γ for naringin-induced treatment effects. In vivo, the PPAR-γ inhibitor BADGE was used to investigate the role of PPAR-γ in naringin-induced improvement of colitis. As shown in Figure 8A, PPAR-γ expression was significantly inhibited by BADGE, and BADGE significantly exacerbated colitis symptoms including weight loss ( Figure 8B), high DAI, increased colon weight-to-length ratio, and histological score ( Figures 8C,D,F). The overall morphology of the colon is shown in Figure 8E. Every colitis symptom was reversed by naringin, but this affect was significantly inhibited by BADGE.
Both RAW264.7 and IEC-6 cells were used to confirm the stimulatory effects of naringin on PPAR-γ. In RAW264.7 cells, pro-inflammatory cytokines, including TNF-α, IL-1β, and INF-γ, were significantly increased by LPS, and the increase was reversed by naringin ( Figures 9A-C). P-NF-κB-p65, a target gene of PPAR-γ, was also significantly increased by LPS ( Figures  9D-F). The LPS induced increase in p-NF-κB-p65 was significantly reversed by naringin, but was significantly inhibited by siRNA targeting PPAR-γ. Similar results were also obtained using IEC-6 cells ( Figures 9G-I). These results suggested that PPAR-γ is a potential target of naringin.

Effects of Naringin on TNBS-Colitis
A TNBS-colitis model was also established to study the effects of naringin. Compared with the sham group, significant weight loss, high colon weight-to-length ratio, lower food intake, high DAI scores, and increased histological score were observed in the untreated colitis model group ( Figures 10A-D,F). The overall morphology of the colon is shown in Figure 10E. Naringin significantly reversed these colitis symptoms. Higher levels of MPO activity and the expression of pro-inflammatory cytokines  were also reversed by naringin treatment (Figures 10H-K). Naringin ameliorated the above colitis symptoms as evidenced by H and E staining ( Figure 10G). Damaged villi structure (green ring), infiltration of inflammatory cells (black arrow), and submucosal edema (blue ring) shown in H and E staining were alleviated by naringin. Naringin can reverse the increasing expression of inflammation-related proteins (iNOS and p-NF-κB-p65), and apoptosis-related protein cl-caspase3 in TNBS-colitis ( Figures  11D-F). Furthermore, decreased PPAR-γ expression in TNBScolitis was also stimulated by gastric administration of naringin ( Figures 11A-C). In accordance with our DSS-colitis results, naringin also significantly alleviated TNBS-colitis.

DISCUSSION
In this study, we examined the effects and detailed mechanisms of naringin on mouse colitis. We showed that naringin had therapeutic effects on both DSS-and TNBS-induced colitis. Colitis symptoms, including neutrophil infiltration, cytokine profiles, and epithelial barrier dysfunction were significantly ameliorated by naringin. In this study, PPAR-γ was found to be an important target for naringin-induced improvement of colitis.
Based on RNA-seq analysis, 753 mRNAs were identified that differed between untreated colitis and naringin treated groups. We identified 11 important pathways involved using a KEGG-target network. Some of these pathways have previously been shown to be involved in colitis. For example, the disruption of Ca 2+ homeostasis is associated with colonic dysmotility in UC (Wang et al., 2016a). cAMP/ PKA-dependent and independent pathways can reduce the expression of inflammatory mediators and attenuate p38 phosphorylation to treat TNBS-colitis (Sun et al., 2017). The activation of PPAR-α, PPAR-δ, and Scd1 has been shown to play a role in DSS-colitis, and Hmgcs2 has an effect in TNBS-colitis (Bassaganya-Riera et al., 2004;Wang et al., 2016b;Wang et al., 2017;Liu et al., 2018). Cnr1, part of the neuroactive ligand-receptor interaction, is reported to be FIGURE 4 | High-throughput transcriptome sequencing and KEGG pathway enrichment analysis. Naringin (40 mg/kg) was administered to mice by gavage once a day for seven consecutive days and total RNA was extracted from colonic epithelial cells. Rich factor is the ratio of differentially expressed gene numbers annotated in this pathway term to all gene number annotated in this pathway term. Greater rich factor means greater intensiveness. Q-value is corrected p value ranging from 0 to 1. (D) Target-pathway network including candidate 133 targets and 11 KEGG pathways.
Frontiers in Pharmacology | www.frontiersin.org August 2021 | Volume 12 | Article 729414 associated with IBD symptoms (Storr et al., 2010). The expression of Mylk in the calcium signaling pathway is a molecular marker of intestinal fibrosis, a critical complication of CD (Rodansky et al., 2015). In our study, PPAR-α, PPAR-γ, PPAR-δ, Scd1, and Hmgcs2 involved in PPAR signaling pathways were selected for preliminary research. Although we found that PPAR-γ is the main target for naringin-induced improvement of colitis, there are multiple targets involved. The other potential targets of naringin will be studied in the future. PPAR-γ is a transcription factor that plays an important role in anti-inflammation, antioxidant, and phagocytemediated cleanup processes and is highly expressed in colonic epithelium and adipose tissue (Lefebvre et al., 1999; FIGURE 5 | High-throughput transcriptome sequencing and KEGG pathway enrichment analysis. Naringin (40 mg/kg) was administered to mice by gavage once a day for seven consecutive days and total RNA was extracted from colonic epithelial cells. (A) Volcano plot of differentially expressed mRNAs in model group vs. control group. (B) The scatter plots showing the top 20 enriched KEGG pathways in model group vs. control group. The Rich factor is the ratio of differentially expressed gene numbers annotated in this pathway term to all gene number annotated in this pathway term. Greater rich factor means greater intensiveness. Q-value is corrected p value ranging from 0 to 1.
Frontiers in Pharmacology | www.frontiersin.org August 2021 | Volume 12 | Article 729414 Dubuquoy et al., 2002). Several studies have demonstrated that expression of PPAR-γ is reduced in UC patients (Dubuquoy et al., 2003;Pedersen and Brynskov, 2010). In the present study, PPAR-γ expression was significantly decreased in both DSS-and TNBS-induced colitis. Reduced PPAR-γ expression may lead to an increase in expression of nuclear transcription factor NF-κB. PPAR-γ can directly bind to NF-κB p50/NF-κB p65 dimer and inhibit the degradation of IκBα, thereby blocking the activation of NF-κB and its nuclear translocation (Chen et al., 2014). Gavage administration of naringin significantly alleviated colitis symptoms in two animal models. Increased MPO activity and higher levels of pro-inflammatory factors (TNF-α, IL-1β, and INF-γ) in colitis were significantly reversed by naringin. High expression of the inflammation-related protein iNOS in colitis was significantly decreased by naringin. Naringin also significantly reduced cl-caspase3 expression to maintain epithelial barrier. Naringin induced alleviation of colitis symptoms was abolished by PPAR-γ inhibitor BADGE, which suggests that PPAR-γ may be a target for naringin induced therapeutic effects on colitis.
In RAW264.7 cells, higher levels of pro-inflammatory factors (TNF-α, IL-1β, and INF-γ) were significantly reversed by naringin. In IEC-6 cells and RAW264.7 cells, LPS-induced high expression of NF- κB-p65 was significantly reduced by naringin, and the reduction was inhibited by siRNA targeting PPAR-γ. All these results suggests that naringin exert treatment effects on colitis through inducing PPAR-γ activation. It has been reported that some PPAR-γ agonists, such as rosiglitazone, pioglitazone, exert inhibitory effects on inflammatory responses (Nitta et al., 2013;Zhou et al., 2021). In this study, we aimed to search and identify natural products which may serve as PPAR-γ agonists. Compared with chemical drugs, naringin is a monomer isolated from citrus fruit species with better safety. The two animal models of colitis were used in this study to confirm the treatment effects of naringin. DSS-colitis is a type of acute injury and it is caused by disrupting the integrity of the mucosal barrier. TNBS-colitis is characterized by T-cell-mediated immunity against haptenized proteins and luminal antigens (Gadaleta et al., 2017). More or less, DSS-induced intestinal inflammation is mainly mediated by an excessive Th1 T-cell response (CD-like IBD), while TNBS-induced intestinal inflammation is a Th2-mediated inflammation (UC-like IBD) (Strober et al., 2002). PPAR-γ activation can inhibit the inflammatory response by inhibiting NF-κB. It is also known that PPAR-γ upregulating tight junction improves the integrity of intestinal mucosal barrier (Huang et al., 2020). In this study, naringin had a good therapeutic effect on both types of colitis through activating PPAR-γ. We used transcriptome sequencing to detect the expression of genes only in naringin treated DSS-induced C57BL/6 mice. Whether the mechanisms of naringin treatment of TNBS-colitis is the same as that of DSS-colitis still needs further investigation.
In conclusion, we uncovered the drug target net of naringin for alleviating DSS-and TNBS-induced colitis. We established a KEGGtarget network to explore detailed mechanisms involved in naringin-induced amelioration of colitis symptoms. PPAR-γ may be a main potential target of naringin. There are still some limitations for our study, for example, the interaction between naringin and PPAR-γ is not clearly studied. The future plan is to extensively explore the interaction between naringin and PPAR-γ, including the exact binding site of PPAR-γ for naringin. The chemical structure of monomers isolated from Chinese medicinal plants is complex, and they can exert pharmacological effects through multiple targets. RNAseq may be a good method to help researchers to construct the drug target net for monomers isolated from Chinese medicinal plants.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: BioProject accession: PRJNA741857. Data are expressed as mean ± SD. Values in the sham group are set to 100% and other values are given relative to those in the sham group. **p < 0.01 compared with sham group; # p < 0.05, ## p < 0.01 compared with TNBS-colitis group; n 6 samples.

ETHICS STATEMENT
The animal study was reviewed and approved by the Laboratory Animal Ethics Committee of Dalian Medical University. Written informed consent was obtained from the owners for the participation of their animals in this study.