Preliminary Study to Explore the Immune-Enhancement Mechanism of Platycodon grandiﬂorus Extract through Comparative Transcriptome Analysis

: Platycodon grandiﬂorus has huge potential medicinal applications in different diseases and immune enhancement (IE). In the past, limited studies have been conducted to decipher the molecular mechanism behind IE of Platycodon grandiﬂorus extract (PGE) despite the fact that IE can be one of the important factors for application of PGE in different diseases including cancer. In this study, whole transcriptome expression and subsequent analysis in PGE treated/untreated organisms were carried out to explore the molecular mechanism behind IE, and subsequently anticancer effect. Current study has the possibility of interaction of PGE treatment with diet, so further study is required to validate the results. The immunity-related genes were found to be up-regulated in the differential expression analysis, which is in line with the IE potential of PGE and the literature. Similarly, Siglec-15, the immune-suppressor gene that has recently been suggested as the normalization target in cancer, was one of the most down-regulated genes. In a set analysis with immune-suppressor genes, the high number of immune-suppressor genes were found to be down-regulated, which indicated that down-regulation of immune-suppressor genes can be the major mechanism behind the IE in PGE treatment. Identiﬁed genes, important in immunity and cancer, are highly recommended for anticancer and IE targets in future studies. down-regulated immune-suppressor genes.


Introduction
Platycodon grandiflorus is a perennial flowering plant commonly known as balloon flower found in East Asia, including Korea, China, Japan and some parts of Russia. P. grandiflorus is the only member of the genus Platycodon and belongs to the family Campanulaceae. P. grandiflorus is the medicine food homology species that has been consumed as both food and medicine from ancient time mostly in Korea, China and Japan [1]. Although P. grandiflorus also has high ornamental value as the flower buds are balloon shaped and found in different colors such as blue, purple, red, violet, white or pink, their extraordinary medicinal properties make the plant highly valuable for anti-hypertension, anti-tumor, anti-obesity, anti-atherogenic, anti-inflammatory, immune enhancement, and reduction in blood sugar levels and blood pressure [1]. Generally, saponins (includes Platycodin A, Platycodin C, Platycodin D, Platycodin D2, Platycodin D3, Deapioplatycodin D, 16-oxo-platycodin D, Polygalin D), flavonoids such as (Luteolin-7-O-glucoside Apigenin-7-O-glucoside), Lobetyol palmitic acid, oleic acid, and gamma-aminobutyric acid are considered as the bioactive component in P. grandiflorus extract (PGE) [1]. The anticancer property of PGE against many types of cancer and metastasis can be considered as one of the most prospective properties of PGE that has potential for cancer treatment [2][3][4]. Diverse medicinal properties of PGE can be primarily attributed to immune-enhancement properties of PGE, which result in recovery against a number of diseases, including different types of cancers [1]. Immune-enhancement effects of PGE were observed in vitro and in vivo studies, which showed activation of both innate and adaptive immune responses in immunosuppressive conditions [5][6][7][8]. The effect of PGE treatment on expression of genes related to immune response can help to understand of the molecular mechanism of PGE in immune enhancement and subsequently other important medicinal properties of PGE i.e., anticancer activity against various cancers and metastasis. Concerning the importance of the immune-enhancing property of PGE, our previous studies found that cell proliferation of macrophages and splenocytes correlates positively with the PGE treatment. It also increased the phagocytic activity of macrophage cells, increased interleukin levels and induced phosphorylation of the p38 mitogen-activated protein kinase, which activated the nuclear factor-kappa B, as well as production of nitric oxide and cytokines [6]. In another study with treatment of fermented products of PGE, we also observed similar effects as the activation and phosphorylation levels of nuclear factor kappa B (NF-κB) were found to be increased [5,9]. In recent experiments with the rat model, the researchers found that PGE treatment increased the cell viability, activity of natural killer cells and cytotoxic T lymphocytes. They also observed the increased production of cyclophosphamide-induced cytokines and different immunoglobulins in splenocytes [5,9]. Overall the study found that PGE recovers the damage in the spleen and thymus caused by cyclophosphamide.
Although several studies were carried out on the role of PGE in immune enhancement, these studies were mainly focused on exploring the macrophage and other immune cell activity, levels of cytokines and interleukins [7,8] to confirm the immune-enhancement property. Limited studies have been performed to explore the molecular mechanism of immune enhancement, whereas no study was carried out to analyze the whole transcriptome expression in PGE treatment for a broader view on the mechanism of immune enhancement. Studying the expression profile of the different genes in the whole transcriptome can help to elucidate the mechanism of immune enhancement and the effects on different pathways in PGE treatment. The molecular mechanism of immune enhancement in PGE treatment can be of great importance in optimal use of PGE in immune enhancement and subsequently its use as an anticancer treatment and other pharmacological applications. High-throughput technology such as next-generation sequencing made it possible to study the whole transcriptome; RNA-Seq technology has especially proven its potential. RNA-Seq has been successfully used in several whole transcriptome studies to explore the molecular mechanism in the immune response in different organisms [10][11][12][13] The RNA-Seq study coupled with computational analysis can provide a comprehensive understanding of molecular mechanisms involved in different and specific biological phenomena and disease conditions with the information on gene expression and function [12,[14][15][16][17]. Pre-existing RNA-Seq datasets in the public repository can also be helpful to unveil and/or validate the molecular mechanism reported from related studies [12]. Hence, the available RNA-Seq dataset from PGE-treated organisms was considered in the current research to study expression of various genes. The differentially expressed gene playing a role in immunity-related functions can help to understand the mechanism of immune enhancement in PGE treatment for optimal use of PGE in different diseases and discovery of new targets for therapeutics.

Materials and Methods
The methodology for the current research was designed and followed in a stepwise manner to understand the transcriptomic changes important in immune modulation occurring in PGE treatment ( Figure 1).

Materials and Methods
The methodology for the current research was designed and followed in a stepwise manner to understand the transcriptomic changes important in immune modulation occurring in PGE treatment ( Figure 1).

RNA-Seq Data Sets Used in the Study
The only available RNA-Seq dataset of the PGE-treated organisms was extracted from the NCBI-SRA database [18]. The dataset was in three groups, i.e., RNA-Seq datasets from epididymal white adipose tissue (eWAT) of mice fed with Normal Diet (ND), High Fat Diet (HFD) and PGE-treated HFD (PTFD). Therefore, 3 samples from every group were used in the study for the RNA-Seq analysis as the possible triplicate in the experiment, i.e., SRR8581331, SRR8581334 and SRR8581337 in group PTFD; SRR8581333, SRR8581335 and SRR8581336 in ND group; and SRR8581332, SRR8581338 and SRR8581339 in HFD group. Finally, two comparisons (comparison-I and comparison-II) were considered for the study according to PGE treatment and diet provided to the mice. In comparison-I, the PGE-treated HFD (PTFD) mice were grouped as PGE Treated group (PTG) and HFD mice as PGE Not-Treated group (PNTG). In comparison-II, again two groups were considered-the first group was the same as comparison-I PTFD as PTG, and the second group, i.e., PNTG was created through merging the ND and HFD mice to make comparison-II more challenging, as PNTG group in comparison-II has both normal-diet and high-fat diet mice.

RNA-Seq Data Sets Used in the Study
The only available RNA-Seq dataset of the PGE-treated organisms was extracted from the NCBI-SRA database [18]. The dataset was in three groups, i.e., RNA-Seq datasets from epididymal white adipose tissue (eWAT) of mice fed with Normal Diet (ND), High Fat Diet (HFD) and PGE-treated HFD (PTFD). Therefore, 3 samples from every group were used in the study for the RNA-Seq analysis as the possible triplicate in the experiment, i.e., SRR8581331, SRR8581334 and SRR8581337 in group PTFD; SRR8581333, SRR8581335 and SRR8581336 in ND group; and SRR8581332, SRR8581338 and SRR8581339 in HFD group. Finally, two comparisons (comparison-I and comparison-II) were considered for the study according to PGE treatment and diet provided to the mice. In comparison-I, the PGE-treated HFD (PTFD) mice were grouped as PGE Treated group (PTG) and HFD mice as PGE Not-Treated group (PNTG). In comparison-II, again two groups were consideredthe first group was the same as comparison-I PTFD as PTG, and the second group, i.e., PNTG was created through merging the ND and HFD mice to make comparison-II more challenging, as PNTG group in comparison-II has both normal-diet and high-fat diet mice.

RNA-Seq Analysis for Assembly and Differential Expression
All pair-end reads from nine different samples of both (PNTG and PTG) groups were converted into fastq reads through an SRA tool kit [19]. Furthermore, the fastq reads were analyzed through AfterQC for the quality control [20]. The suggested parameters were used in AfterQC analysis, which includes filtering, trimming, error removing and quality control of the reads [20]. Finally, only the good quality reads according to suggested parameters from quality-control analysis were subjected to alignment with genome using an HISAT2 alignment tool [21]. GRCm38, the mouse (Mus musculus) reference genome assembly released by the Genome Reference Consortium, was used for alignment of reads for all the samples that resulted in Sequence Alignment Map (SAM) files. SAM files were sorted and converted to bam files as per requirement for further assembly process using SAMtools [22]. The assembly was done using StringTie [23] with the -e option, which is required to combine the results of all samples for differential expression analysis. To combine the assembly results for all selected nine samples from both groups (PTG and PNTG) python script (prepDE.py) was used. Further, iDEP 9.1 was used for different expression-related analysis through EdgeR and DESeq2 [24][25][26]. Suggested parameters for differential expression (i.e., FDR < 0.1 and minimum 2-fold change) were used to identify the DEG in both comparisons (comparison-I and comparison-II) [24,25].

Function Enrichment of DEG
All differentially expressed genes (DEG) from both comparisons were subjected to functional enrichment analysis using PANTHER (protein analysis through evolutionary relationship), which is based on biological process, cellular component, molecular function, protein class and pathway enrichment [27,28]. Each analysis was carried out with suggested statistical cutoff, and the results were stored in tab-separated variable format (TSV) files [27,28].

Identification of DEG with Immunity Related Function in Subset Analysis
Two different immunity-related gene sets were prepared and compared with identified DEG between PTG and PNTG groups in both the comparisons. Information of all 871 immunity-related mouse genes present in IMGT/GENE-DB [29] corresponding to 1321 allele was extracted in tabular format. Furthermore, the gene names were converted to gene accession numbers using g:Profiler server [28]. All these genes were used to create immunity-related gene set-1. Immunity-related gene set-2, which is the collection of immune-suppressive genes, was created with all genes of mouse homologous to immunosuppressive genes reported in HisgAtlas 1.0 database (a human immunosuppression gene database) [30]. All immunosuppressive genes present in HisgAtlas database were mapped to homologous genes in mouse using g:Profiler server [28]. Names of these immune-suppressive genes in mouse were converted to gene accession numbers to create immunity-related gene set-2. These immunity-related gene set-1 and gene set-2 were used to identify common up-and down-regulated DEG in both the comparisons considered in the study from set analysis [31]. The list of ensembl ids for up-and down-regulated genes in both comparison and prepared gene set-1 and gene set-2 was used in set analysis through InteractiVenn [31].

Protein-Protein Interaction Network
All down-regulated common immune-suppressive genes identified in both comparisons were subjected to protein-protein interaction analysis using STRING version 11.0 [32]. Protein-protein interaction and KEGG network enrichment analysis were carried out using Mus musculus, which was selected as the target organism [33]. The ensemble ids of genes were used as input, and default parameters were utilized for the analysis. According to output, the results were stored in TSV and image format files.

Preprocessing and Alignment of Reads to Reference Genome
Different analyses were carried out on all nine samples to obtain the whole transcription data in the form of assembled genes/transcripts for further differential expression and other downstream analysis ( Figure 1). The preprocessing and quality control were carried out to select the good quality reads for alignment of the reads with reference genome, and finally more than 98.7% reads (with average of more than 99.5%) were obtained as good reads from all the selected samples (Table 1). After preprocessing, only the good quality reads were subjected to alignment with efficient method HISAT2 using the mouse reference genome GRCm38. Finally, high-alignment rate was achieved in all samples used in the study, i.e., more than 98% overall alignment rate in the alignment (Table 1).

Assembly and Expression Analysis
The alignment was used to assemble the transcripts through StringTie, which also calculates the read counts and fragment per kilobase of transcript per million map reads (FPKM) value corresponding to every assembled gene/transcript. Finally, the count matrix was made from the StringTie output of every sample as count matrix is the recommended input for the differential-expression analysis [23]. The read count matrix of comparison-I and comparison-II had six samples and nine samples, respectively, which were used for further analysis in iDEP 9.1. In expression analysis, genes with 0.5 minimal counts per million (CPM) in at least one library were selected and taken for the clustering and heatmap analysis. Heatmaps were generated using 1000 most differentially expressed genes in both the comparisons (comparison-I and comparison-II). Hierarchical clustering within the heatmaps clustered the PTG and PNTG groups in separate clusters in both the comparisons (Figures 2 and 3). In differential expression analysis, a total of 1951 and 130 genes were found to be up-regulated in comparison-I and comparison-II, respectively (Table D1 and Table D2). Similarly, 2396 and 859 genes were found to be down-regulated in comparison-I and comparison-II, respectively (Table D1 and Table D2). Different plots such as MA plot, volcano plot and scattered plot were also constructed from the DEGs in both the comparisons (Figures 2 and 3).

Functional Enrichment Analysis
Identified DEGs were subjected to functional enrichment analysis using PANTHER (protein analysis through evolutionary relationship) on the basis of biological process, cellular component, molecular function, protein class and pathway enrichment. In comparison-I, the DEGs were enriched in 8 (transcription regulator activity, molecular transducer activity, binding, structural molecule activity, molecular function regulator, catalytic activity, translation regulator activity and transporter activity) different molecular functions, 21 different biological processes, 14 different cellular components, 21 diverse protein classes and 137 various pathways (Tables S1a-S5a). Similarly, in comparison-II the DEGs were enriched in 7 (transcription regulator activity, molecular transducer activity, binding, structural molecule activity, molecular function regulator, catalytic activity and transporter activity) different molecular functions, 20 different biological processes, 14 different cellular components, 19 different protein classes and 85 different pathways (Tables S1b-S5b).

DEG with Immunity Related Function
In set analysis, among the up-regulated genes the 61 and 8 genes from comparison-I and comparison-II, respectively, were found to be present in gene set-1, i.e., listed in mouse immunity-related genes (IMGT/GENE-DB) [29]; no gene was listed in gene set-1 within down-regulated genes from both the comparisons (comparison-I and comparison-II). Similarly, 219 and 125 genes among down-regulated genes from comparison-I and comparison-II, respectively, were found to be listed in gene set-2, which consisted of known immunosuppressive genes. All down-regulated immune-suppressive genes (125 genes) in comparison-II were also found to be common in 219 down-regulated immune-suppressive genes in comparison-I. However, only 98 and 6 genes from comparison-I and comparison-II, respectively, were found to be present ingene set-2 (Figures 4 and 5). II). Similarly, 219 and 125 genes among down-regulated genes from comparison-I and comparison-II, respectively, were found to be listed in gene set-2, which consisted of known immunosuppressive genes. All down-regulated immune-suppressive genes (125 genes) in comparison-II were also found to be common in 219 down-regulated immunesuppressive genes in comparison-I. However, only 98 and 6 genes from comparison-I and comparison-II, respectively, were found to be present ingene set-2 (Figures 4 and 5).

Discussion
P. grandiflorus has great potential for medicinal use in various diseases due to different pharmacological properties. Immune enhancement is one of the major outcomes of PGE treatment, which is considered to be one of the important factors behind the positive conclusion in different cancers and other diseases such as cough, phlegm and tuberculosis [1]. The immune-enhancement mechanism of PGE is encouraging the researchers to elucidate the PGE and/or its major constituents (platycodin D, etc.) as a possible treatment in many chronic diseases [1]. As immune enhancement is a complex process, its mechanism is currently attributed to the combined effect of several genes and different pathways. Different studies have reported immune enhancement in both innate and adaptive immune responses through increased activity of immune-system cells and molecules [5][6][7]; no study has been carried out to study immune enhancement through whole transcriptome analysis in PGE-treated organisms, which can elucidate the molecular mechanism and genes behind the immune enhancement.
Hence, to study the role of PGE on the whole transcriptome for immune enhancement, we analyzed the PGE-treated RNA-Seq datasets from PGE-treated mice in two comparisons. In comparison-II, the mice fed on PGE and HFD were compared with mice fed on HFD as well as normal diet. Having both normal and high-fat diet mice in PNTG group in comparison-II made this comparison more challenging, but it was necessary to minimize the effect of HFD in PGE-treated mice. The possibility of interaction between the HFD and PGE treatment cannot be ruled out, but unavailability of transcriptome data from normal PGE-treated organisms limits the current study. Thus, this preliminary study must be further extended according to availability of datasets in the future. The animal study with normal diet and PGE treatment can be used to generate whole transcriptome

Discussion
P. grandiflorus has great potential for medicinal use in various diseases due to different pharmacological properties. Immune enhancement is one of the major outcomes of PGE treatment, which is considered to be one of the important factors behind the positive conclusion in different cancers and other diseases such as cough, phlegm and tuberculosis [1]. The immune-enhancement mechanism of PGE is encouraging the researchers to elucidate the PGE and/or its major constituents (platycodin D, etc.) as a possible treatment in many chronic diseases [1]. As immune enhancement is a complex process, its mechanism is currently attributed to the combined effect of several genes and different pathways. Different studies have reported immune enhancement in both innate and adaptive immune responses through increased activity of immune-system cells and molecules [5][6][7]; no study has been carried out to study immune enhancement through whole transcriptome analysis in PGE-treated organisms, which can elucidate the molecular mechanism and genes behind the immune enhancement.
Hence, to study the role of PGE on the whole transcriptome for immune enhancement, we analyzed the PGE-treated RNA-Seq datasets from PGE-treated mice in two comparisons. In comparison-II, the mice fed on PGE and HFD were compared with mice fed on HFD as well as normal diet. Having both normal and high-fat diet mice in PNTG group in comparison-II made this comparison more challenging, but it was necessary to minimize the effect of HFD in PGE-treated mice. The possibility of interaction between the HFD and PGE treatment cannot be ruled out, but unavailability of transcriptome data from normal PGE-treated organisms limits the current study. Thus, this preliminary study must be further extended according to availability of datasets in the future. The animal study with normal diet and PGE treatment can be used to generate whole transcriptome datasets that would be most appropriate for the immune-enhancement study. Furthermore, we included Appl. Sci. 2021, 11, 226 9 of 13 all of the nine samples in comparison-II, which increased the statistical significance of this study. High alignment rate was achieved in the alignment of reads with the reference genome (>98% overall alignment rate) in all nine samples, which signposted the good quality of samples used in the study ( Table 1). The state-of-art RNA-seq analysis tools such as HISAT2, StringTie, iDEP9.1, DESeq2 and EdgeR were used [21,23] for different assembly and expressions analysis (Figure 1). Almost all (127 out of 130) up-and (853 out of 859) down-regulated genes in comparison-II were common with up-and down-regulated genes in comparison-I, which suggest the common pattern in both comparisons (Table D1 and  Table D2). Siglec-15 was one of most down-regulated genes in both comparisons and was recently discovered as the therapeutic target in normalization cancer immunotherapy [34]. Down-regulation of Siglec-15 may be the main reason behind both immune enhancement and anticancer activity of PE. In a very recent landmark study, down-regulation of this immune-suppressor gene was found to suppress antigen-specific T cell responses in both in vitro and in vivo conditions [34]. Blockage of Siglec-15 through specific monoclonal antibodies reversed the suppression of T cells, enhanced tumor immunity and inhibited the growth of tumors in multiple established tumor models [34]. Down-regulation of Siglec-15 protein molecules can also be responsible for immune-enhancement activity of other plant extracts or compounds of unknown mechanism. Further research is highly required to study this candidate gene/protein.
All DEGs were used for the functional enrichment analysis through prediction of five important biological categories to predict genes important in immunity-related functions. In comparison-I, 150 DEGs were enriched in immune-system processes in the biological function category (Table S1a) which can be candidate genes in the immune-modulation effect of PGE. Similarly, 82 DEGs were enriched in defense/immunity protein (PC00090) in the protein-class category (Table S4b). In pathway enrichment analysis 70, 26,21,18,12,4 DEGs were found to be from inflammation mediated by the chemokine-and-cytokine signaling pathway (P00031), B cell activation (P00010), T cell activation (P00053), the interleukin signaling pathway (P00036), the toll receptor signaling pathway (P00054), and the interferon-gamma signaling pathway (P00035), respectively ( Table S5b). All of these enriched immune-system-related pathways directly modulate the immune function and response [35,36]. Similarly, in comparison-II the functional enrichment 51 DEGs were enriched in immune-system processes in the biological-function category (Table S1b). In addition, 29 DEGs were enriched in defense/immunity protein (PC00090) in the proteinclass category (Table S4b). In pathway-enrichment analysis 33,13,9,5,4,2 DEGs were found to be from inflammation mediated by the chemokine-and-cytokine signaling pathway (P00031), B cell activation (P00010), T cell activation (P00053), toll-receptor signaling pathway (P00054), interleukin-signaling pathway (P00036) and interferon-gamma signaling pathway (P00035), respectively (Table S5b). Molecular function of these candidate DEGs can be studied in future studies, which may unveil their precise role in immune enhancement [35,36].
To further explore the role of DEGs in context to immune modulation, two gene sets, i.e., gene set-1 and set-2 were created and used in set analysis. Gene set-1 was composed of immunity genes present in IMGT/GENE-DB [29] that are involved in immune function, and gene set-2 comprised genes that suppress the immune function [30]. In set analysis, it was found that no genes were down-regulated from gene set-1 in both the comparisons. However, 61 and 8 genes were found to be up-regulated in comparison-I and comparison-II, respectively, which was consistent with the immune-enhancement properties PGE observed in different studies [5,6]. These genes were mostly composed of immunoglobulin genes that directly act as receptors in immune response [29]. These immunoglobulins were also found to be up-regulated in earlier studies carried out to study immune enhancement in PGE treatments [5][6][7]. Similarly, when gene set-2 was used in set analysis, a high number of immune suppressive genes were found to be down-regulated compared to up-regulated genes in both the comparisons (Figures 4 and 5). The mechanism that downregulates immune-suppressive genes is expected to increase immunity-related function and immune enhancement [30]. Therefore, a higher number of down-regulated immunesuppressive genes again justifies the immune-enhancing effect of PGE treatment [30]. All (125) down-regulated immune-suppressive genes in comparison I were also found to be down-regulated in comparison II. Expression of immunosuppressive genes can be used to study immune modulation in future high-throughput studies as the current finding is in line with the immune-enhancement effect of PGE. Furthermore, the protein-protein interaction (PPI) analysis was carried out in addition to the function enrichment and pathway analysis for common down-regulated immune-suppressor genes. The protein-protein interaction network of these common down-regulated immune-suppressor genes resulted in a highly connected network compared to random possible connections (Figure 6), which suggests that there must be some pattern in expression of genes in our study that may be attributed to the PGE treatment and immune-suppressive nature of the genes. Furthermore, KEGG enrichment analysis through the STRING PPI server revealed that three cancer-related pathways were associated with 16 immune-suppressive genes (Table S7) [32,33]. Most of these genes such as Cd44, Csf1r, Cxcr4, Fn1, Il7r, Lpar5, Mmp2, Pak1, Rac2, and Tlr2 were found to be studied in cancer and their expressions were associated with cancer progression and metastasis as expected, and these genes were proposed as the target in cancer [37][38][39][40][41][42][43][44][45][46]. Therefore, these immune-suppressor genes are proposed as the important candidates for further study in cancer application along with Siglec-15 (Table S7). Although our current study is based on only available RNA-Seq datasets from the eWAT of mice [18], there is a possibility of extending our study to further detailed research on other tissues or cells through deep sequencing and/or single cell sequencing to consider precise transcriptomic changes occurring in immune cells such as B-cell, T-cell or other immunity-related cells [47]. However, the preliminary nature of the current study further requires the validation of proposed genes in functional assays, which may include inhibition and/or overexpression studies. The current research has provided the strong basis of the immune-enhancement effect of PGE in terms of down-regulation of immune-suppressive and up-regulation of immune-enhancing genes that warrants the further research on Siglec-15 and other 16 immunosuppressive genes involved in cancer.

Conclusions
Comparative analysis of the whole transcriptome of PGE-treated mice was carried out successfully with the help of available transcriptomic resources. Up-regulation of immunity-related genes/proteins such as immunoglobulins, was in line with earlier reported immune-enhancement effects of PGE and justified our current study. Downregulation of immune-suppression genes was identified as one of the major mechanisms in immune enhancement in PGE treatment. One of the most down-regulated genes (Siglec-15) was an immunosuppressive gene that was recently discovered as a therapeutic candidate in cancer. The results of the current study can have limitations from the possible interaction of HFD and PGE treatment, so further in vitro and in vivo studies are also required to generate the appropriate transcriptome datasets and further validation. Our finding can be considered an active pathway of mechanisms important for both immune enhancement and anticancer application of PGE.
Supplementary Materials: The following are available online at https://www.mdpi.com/2076-3 417/11/1/226/s1, Table S1: title. Table S1a: Results of DEGs in panther enrichment analysis for biological function in comparison-I, Table S2a: Results of DEGs in panther enrichment analysis for molecular function in comparison-I, Table S3a: Results of DEGs in panther enrichment analysis for cellular components in comparison-I, Table S4a: Results of DEGs in panther enrichment analysis for protein class in comparison-I, Table S5a: Results of DEGs in panther enrichment analysis for pathways in comparison-I, Table S1b: Results of DEGs in panther enrichment analysis for biological function in comparison-II, Table S2b: Results of DEGs in panther enrichment analysis for molecular function in comparison-II, Table S3b: Results of DEGs in panther enrichment analysis for cellular components in comparison-II, Table S4b: Results of DEGs in panther enrichment analysis for protein class in comparison-II, Table S5b: Results of DEGs in panther enrichment analysis for pathways in