Search tips
Search criteria 


Logo of bmcamBioMed Centralsearchsubmit a manuscriptregisterthis articleBMC Complementary and Alternative Medicine
BMC Complement Altern Med. 2014; 14: 14.
Published online 2014 January 11. doi:  10.1186/1472-6882-14-14
PMCID: PMC3893506

Sho-saiko-to, a traditional herbal medicine, regulates gene expression and biological function by way of microRNAs in primary mouse hepatocytes



Sho-saiko-to (SST) (also known as so-shi-ho-tang or xiao-chai-hu-tang) has been widely prescribed for chronic liver diseases in traditional Oriental medicine. Despite the substantial amount of clinical evidence for SST, its molecular mechanism has not been clearly identified at a genome-wide level.


By using a microarray, we analyzed the temporal changes of messenger RNA (mRNA) and microRNA expression in primary mouse hepatocytes after SST treatment. The pattern of genes regulated by SST was identified by using time-series microarray analysis. The biological function of genes was measured by pathway analysis. For the identification of the exact targets of the microRNAs, a permutation-based correlation method was implemented in which the temporal expression of mRNAs and microRNAs were integrated. The similarity of the promoter structure between temporally regulated genes was measured by analyzing the transcription factor binding sites in the promoter region.


The SST-regulated gene expression had two major patterns: (1) a temporally up-regulated pattern (463 genes) and (2) a temporally down-regulated pattern (177 genes). The integration of the genes and microRNA demonstrated that 155 genes could be the targets of microRNAs from the temporally up-regulated pattern and 19 genes could be the targets of microRNAs from the temporally down-regulated pattern. The temporally up-regulated pattern by SST was associated with signaling pathways such as the cell cycle pathway, whereas the temporally down-regulated pattern included drug metabolism-related pathways and immune-related pathways. All these pathways could be possibly associated with liver regenerative activity of SST. Genes targeted by microRNA were moreover associated with different biological pathways from the genes not targeted by microRNA. An analysis of promoter similarity indicated that co-expressed genes after SST treatment were clustered into subgroups, depending on the temporal expression patterns.


We are the first to identify that SST regulates temporal gene expression by way of microRNA. MicroRNA targets and non-microRNA targets moreover have different biological roles. This functional segregation by microRNA would be critical for the elucidation of the molecular activities of SST.

Keywords: Sho-saiko-to, Microarray, MicroRNA, Hepatocyte, Pathway, Promoter


Sho-saiko-to (SST) (also known as so-shi-ho-tang or xiao-chai-hu-tang) is a botanical formulation composed of seven herbal materials (see Additional file 1: Table S1) and is widely used for the treatment of chronic hepatitis and liver cirrhosis in Korea, Japan, and China [1]. SST and its major components (e.g., baicalin, baicalein, glycyrrhizin, and saikosaponin-D) have marked antiproliferative activity on hepatocellular carcinoma [2-4], prevent liver injury [5], and promote liver regeneration in animal models [6,7]. These pharmacologic effects of SST involve the immunomodulation of diverse immune cells and immune molecules [8,9]. However, because of the complex nature of the chemical components of SST, focusing only on specific components or on a few target genes is inadequate to understand the diverse biological activities of SST. Therefore, it is necessary to apply a multiple target-based approach to elucidate the molecular mechanisms mediated by the multiple components of SST.

Recent advances in high-throughput technology such as the microarray has made it possible to investigate the effects of drugs at the whole-genome level [10]. One high-throughput technology is the microRNA array, which can detect the expression level of whole microRNAs that have been discovered to date [11]. MicroRNA is a small noncoding RNA molecule composed of approximately 22 nucleotides that pair to sites in messenger RNA (mRNA) and directly repress post-transcription in eukaryotic cells [12]. Many reports suggest that microRNAs are involved in diverse biological functions such as proliferation, differentiation, and development. The search for targets of microRNA shows that many mammalian mRNAs are the conserved targets of microRNA [13]. This suggests an important role of microRNA in regulating gene expression. Therefore, using the information of mRNA and microRNA is important to elucidate the precise mechanism of gene expression. The integrated multi-omics approach actually reveals a novel regulatory network of gene expression in diverse biological situations such as disease research [14-16], genome research [17], and herbal research [18,19]. We also previously reported the usefulness of a genome-wide approach in elucidating the molecular effects of herbal extracts [20,21].

By using an integrated genomic analysis of genes and microRNAs in the present study, we attempted to identify SST-induced gene expression changes in primary mouse hepatocytes. The results indicated that SST regulated gene expression through microRNA in a functionally coordinated manner. Our approach could give perspective on the role of microRNAs in the pharmacological effects of SST.


Primary mouse hepatocyte isolation and culture

Six-week-old male ICR mice were purchased from Samtako Bio Inc. (Seoul, Korea). Primary mouse hepatocytes were prepared by using the collagenase perfusion method. In brief, the 6-week-old male mice were anesthetized by an intraperitoneal injection of Zoletil-50 and 2% Rompun, which were cannulated through the right ventricle. The livers were perfused with ethylene glycol tetra-acetic acid (0.5 mM) in Hepes-buffered Hank’s balanced salt solution (HBSS; pH 7.4) for 5–6 min (flow rate 5 mL/min). The livers were then perfused for another 20 min with Hepes-buffered HBSS containing collagenase (Sigma, USA) (flow rate 5 mL/min). The hepatocytes were dispersed, washed, and purified on a Percoll density gradient (Sigma). Hepatocyte preparations with viability greater than 85%, as determined by the trypan blue exclusion protocol, were used. The isolated hepatocytes were suspended, and then transferred to gelatin-coated culture dishes or plates at a density of approximately 5 × 105 cells/mL. The hepatocytes were allowed to attach onto culture dishes or plates coated with gelatin for 4–6 hours in William’s Media E (Sigma) containing 1% penicillin/streptomycin, 2 mM of L-glutamine and 10% fetal bovine serum. After the attachment, the hepatocytes were washed with HBSS and provided fresh medium. They were incubated overnight at 37°C, 95% air, and 5% carbon dioxide. The hepatocytes were then deprived of the serum and used for experiments. All animal experimental procedures were approved by Institutional Animal Care and Use Committee of the Korea Institute of Oriental Medicine (Permit Number: KIOM 12–024) and performed in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals at the Korea Institute of Oriental Medicine.

Cell viability assay

SST was kindly provided by Dr. Hyeun Kyoo Shin (Basic Herbal Medicine Research Group, Korea Institute of Oriental Medicine, Republic of Korea). Preparation of SST was described as previously [22]. In brief, crude seven herbal medicines were extracted in distilled water at 100°C for 2 hours, filtered, and then freeze-dried. We confirmed the safety of SST by using an in vitro colorimetric cell proliferation kit (methyl thiazolyl tetrazoliym [MTT]) (Roche Applied Science, Germany) as described previously [23]. In brief, hepatocytes were first cultured in 48-well plates at a density of 1.0 × 105 cells/well for 24 hours. After incubation, the cells were washed with phosphate-buffered saline and treated with different concentrations of SST (0.1–1.0 mg/mL) for 24 hours. The cells were hereafter washed and incubated for 1 hour with MTT (500 μg/mL). Formazan crystals were dissolved by using dimethyl sulfoxide (100 μL/well). The absorbance was measured colorimetrically at 570 nm.

Microarray experiment and quantitative real-time polymerase chain reaction

Mouse primary hepatocytes were treated with 500 μg/mL of SST at a density of 1.0 × 106 cells per 60-mm dish for 1–24 hours in triplication. The total RNA from hepatocytes was isolated with Tri-reagent (Sigma) in accordance with the manufacturer’s instructions. The quality of purified RNA was measured by using the Agilent 2100 Bioanalyzer (Agilent Technologies, USA); only samples with a RNA integrity number (RIN) greater than 7.0 were included in the microarray analysis. RNAs from the triplication of experiments at each time point were pooled to exclude experimental bias. For the gene expression microarray, isolated RNA was amplified and labeled by using the low RNA input linear amplification Kit PLUS and then hybridized to a microarray (Agilent Mouse Whole Genome 44 K; Agilent Technologies, USA) that contained approximately 44,000 probes (approximately 26,600 unique genes) in accordance with the manufacturer’s instructions. For microRNA expression microarray, the microRNA was labeled and hybridized to Agilent Mouse miRNA Microarray (Release 17.0) by using the Agilent miRNA Labeling and Hyb Kit (Agilent Technologies, USA). Approximately 1100 microRNAs, based on the annotation of miRBase Release 17.0, were presented in microarray. The arrays were then scanned with the Agilent Microarray Scanner (Agilent Technologies, USA). For quantitative real-time polymerase chain reaction (Q-PCR) analysis, mRNA and microRNA were reverse-transcribed, amplified, and detected by using Taqman probes (ABI, USA) in triple time, as described previously [24].

Microarray data analysis

The raw intensity of the probe signals was obtained by using Feature Extraction Software (Agilent Technologies, USA). Only array elements showing a signal intensity greater than 1.4-fold of the local background were considered well measured. The remaining elements were normalized using the quantile method [25]. The intensities for duplicated spots were averaged. The expression ratio of genes (or microRNAs) in the experimental samples was then determined by comparing them with genes (or microRNAs) in the control sample. The expression profile was hierarchically clustered by using the Cluster program and visualized using the TreeView program (both can be obtained from Figure 1 shows a schematic diagram of the overall analytical process.

Figure 1
A schematic illustration of the analysis strategy. Temporally expressed genes from the gene expression microarray dataset have been integrated with a microRNA expression microarray dataset in which putative target gene information of microRNA was obtained ...

Temporal expression of genes and microRNAs

The short time-series expression miner (STEM) program—which was originally developed for the temporal analysis of microarray experiments [26] —was used to identify genes or expression patterns that were changed temporally. Only genes with a fold ratio greater than 2 or less than 0.5 for at least one time point were included in the analysis. The statistical significance of the temporal patterns was calculated by using a permutation test (n = 1000) corrected by the false discovery rate (FDR).

Integration of mRNA and microRNA expression

The relationship between gene expression and microRNA expression was measured by using a permutation-based correlation method. First, a list of the predicted target genes of microRNAs, calculated by bioinformatic analysis of large public microarray datasets, was obtained from the MicroRNA Database (miRDB version 4.0) website ( [27,28]. Second, the Pearson correlation coefficient was measured between each microRNA expression in the microRNA microarray and each predicted target gene expression in the mRNA microarray. Only gene and microRNA pairs that showed a negative correlation coefficient were selected to form a correlation coefficient matrix between the predicted target genes and the microRNAs. The statistical significance of the resultant correlation coefficient matrix was estimated by using a random sampling-based permutation [29] in which the coefficient values from the original dataset were compared with the values from 1000 times randomly permuted datasets. Only target genes and microRNAs with a FDR less than 0.01 were selected as significant.

Pathway enrichment

The simple enriched pathways were estimated by the DAVID program [30] in which the p values of each pathway were calculated, based on Fisher’s exact test, from an input list of genes. For adjustment by multiple comparisons, the DAVID program used the FDR by the Benjamini procedure. For another pathway analysis, the Signaling Pathway Impact Analysis (SPIA) program [31] was implemented by using a subgroup of differentially expressed genes. The SPIA program calculated a global pathway significance p value (PG) that combines the enrichment p values and the perturbation p values by considering pathway topology with a random bootstrap iteration number of 3000. The FDR of the pathways was measured by applying the Benjamini algorithm in SPIA. The pathway information was obtained from the database of the Kyoto Encyclopedia of Genes and Genomes (KEGG,

Pathway activity

The activity of the pathways was measured by linearly combining the logarithmic expression value of all genes in each pathway to account for the accumulative effect of small changes by many genes [32]. Statistical significance was measured by the FDR in which the original pathway’s activity was compared with the randomly permutated activity values (1000 times). Pathways with a FDR less than 0.01 were selected as significant and then hierarchically clustered on the basis of similarity of activity values.

Core microRNA targets from multiple pathways

Core nodes (i.e., core genes) among multiple pathways were measured by implementing KEGGgraph R package (version 2.10) [33]. In brief, the core nodes were determined by calculating the relative betweenness centrality of nodes in which the number of ingoing and outgoing edges for each node was computed in the network structure of the multiple pathways. Nodes with a relative betweenness centrality greater than 0.01 were selected as the core microRNA targets.

Transcription factors binding sites analysis

Candidate binding sites for transcription factors in the promoter region were identified through sequence matching of the position weight matrix by implementing MotifDb R package (version 1.2.2, [34]. A total of 329 position weight matrices for mouse transcription factors were used. Of these, 47 matrices were from the JASPAR database ( [35,36] and 282 matrices were from the Universal PBM Resource for Oligonucleotide-Binding Evaluation (UniPROBE) database [37]. The nucleotide sequence of the promoter region of the gene (-2000 bp to +500 bp from the transcription start site) was obtained from the Mus musculus full genome, which was provided by the University of California, Santa Cruz (UCSC mm10 version). The presence of the transcription factor binding site (TFBS) within the promoter region of each gene was predicted by using the matchPWM algorithm in which a minimum score for counting a match was set at 90% [36]. Based on the resultant frequency of the matrices of the TFBS, the similarity of genes was determined by using Jaccard’s algorithm, which does not consider the absence of binding sites in two promoters as an indication of similarity [38]. Jaccard’s algorithm is effective in the promoter clustering of genes, as we previously reported [21].


Temporal pattern of genes and microRNA expression

The cytotoxic effect of SST on primary hepatocytes was not significant under the experimental condition (0.1–1.0 mg/mL) as shows (see Additional file 1: Figure S1). The concentration of SST therefore chosen for the study was 500 μg/mL because of its solubility and cytotoxicity in the microarray analysis. The expression profiles of genes and microRNAs, regulated by the treatment of SST, were measured by using microarray analysis in primary mouse hepatocytes. Figure 1 depicts the overall analysis. The expression pattern of genes shows that 1166 genes were dramatically changed in their expression levels at the time of SST treatment (Figure 2A). Among these patterns of gene expression, Sub-cluster 1 was composed of genes that temporally increased expression, whereas Sub-cluster 2 was composed of genes with temporally decreased expression. For a more systematic approach, we tried to isolate genes showing a specific temporal pattern by using a time-series analysis of the microarray. Figure 2B presents two representative statistically significant temporal patterns: the temporal up-regulated pattern (temporal up-pattern) and the temporal down-regulated (temporal down-pattern); the FDR was less than 0.001, which included most temporally expressed genes that were changed by SST. The temporal up-pattern included 463 temporally up-regulated genes and the temporal down-pattern included 177 genes down-regulated by SST. However, the expression of microRNAs did not show a clear temporal pattern after treatment with SST (Figure 2C) (see Additional file 1: Table S2) shows the full list of temporal pattern genes.

Figure 2
Temporal expression of genes and microRNAs after SST treatment in primary mouse hepatocytes. (A) Approximately 1100 differentially expressed genes with a fold ratio greater than 2 or less than 0.5 for at least one sample are clustered hierarchically. ...

Integration of gene and microRNA expression

To determine the putative targets of microRNA among the genes in the two temporal expression patterns, we developed an algorithm integrating the temporal expressions of the genes and the microRNAs. By using the predicted target genes from the miRDB ( [27,28], correlation coefficients were measured between the temporal expression of the predicted target genes and the microRNAs. After the permutation-based adjustment of the measured correlation coefficients, 174 genes with a FDR less than 0.01 were finally selected as the putative targets of microRNA from the two temporal patterns: 155 genes were identified from the temporal up-pattern and 19 genes were identified from the temporal down-pattern. For experimental confirmation of the expression pattern based on microarray testing, Q-PCR was performed for some genes (ABCC4 from temporal up-pattern, CYP3A11 from temporal down-pattern, and FOXA1 from non-pattern) and microRNAs (miR-23a-3p and miR-466b-3p). CYP3A11 was specifically the target of miR-23a-3p (Table 1). As shows (see Additional file 1: Figure S2), the overall patterns of gene expressions between microarray and Q-PCR were similar. Figure 3 shows the connection map between microRNA molecules and its target genes from temporal up-pattern (Figure 3A) and from temporal down-pattern (Figure 3B). Table 1 lists the microRNA targets. The number of putative microRNA targets (19 of 177 genes) in the temporal down-pattern was significantly lower than the number of targets (155 of 463 genes) in the temporal up-pattern (p value < 0.001). This unbalanced distribution of the microRNA target genes imply microRNAs have a specific biological role induced by SST. Therefore, we measured the functional involvement of microRNA targets via pathway analysis.

Table 1
The microRNA targets regulated by SST
Figure 3
The connection map between microRNAs and target genes altered by SST. Genes showing a negative correlation with microRNA expression were selected as microRNA targets by implementing permutation-based correlation method (i.e., FDR less than 0.01). The ...

Pathway analysis of putative microRNA targets

The pathways involved in the two temporal patterns of the genes were measured by applying two different approaches (Table 2), by simple enrichment pathway analysis, and by topology-based signaling pathways analysis. Simple enrichment analysis of the pathways, which measures enriched pathways from Fisher’s exact test, showed that the temporal up-pattern induced by SST was involved in the cell cycle pathway (i.e., KEGG ID mmu04110) and that the temporal down-pattern included drug metabolism-related pathways (e.g., mmu00982, mmu00983, mmu00980) and immune-related pathways such as the systematic lupus erythematosus pathway (mmu05322) and the complement and coagulation cascade pathways (mmu04610). Topology-based signaling pathway analysis, which calculates the enrichment score by taking into account the topology of each signaling pathway, also showed that the cell cycle pathway (i.e., mmu04110) was significantly enriched from the temporal up-pattern, whereas diverse signaling pathways (e.g., immune-related pathways and metabolism-related pathways) were selected as significant pathways from the temporal down-pattern. In addition, the cytokine-cytokine receptor interaction pathway (mmu04060) and the osteoclast differentiation pathway (mmu04380) were also significant pathways that were associated with the temporal up-pattern.

Table 2
Pathways enriched in temporal patterns by SST

We measured temporal changes in pathway activity by using the expression levels of all genes in each pathway. Figure 4 shows that many diverse pathways were temporally activated or repressed, according to the SST treatment. Pathways enriched from the temporal up-pattern and down-pattern showed temporally increased and decreased activity, respectively.

Figure 4
Pathway activities altered by SST in primary mouse hepatocytes. The temporal change of pathway activity is measured by linearly combining the logarithmic expression value of all genes in each pathway and then clustering them hierarchically. The columns ...

The functional association of microRNA target genes shows that only one pathway—the cell cycle pathway (mmu04110)—was statistically significant from the temporal up-pattern (the FDR was less than 0.01). From the temporal down-pattern, we measured one pathway that was also statistically significant: the xenobiotics metabolism pathway (mmu00980). On the other hand, non-microRNA targets from the temporal up-pattern were associated with signaling pathways such as the cytokine-cytokine receptor interaction pathway (mmu04060), the NF-κB signaling pathway (mmu04064), the mitogen-activated protein kinase (MAPK) pathway (mmu04010), and the osteoclast differentiation pathway (mmu04380). However, non-microRNA targets from the temporal down-pattern were associated with diverse pathways such as immune-related pathways and metabolism-related pathways. (see Additional file 1: Figure S3) shows the positions of the temporally regulated genes in each significant pathway.

By comparing pathways involved in the microRNA targets and microRNA non-targets, we speculated that microRNA was specific for the regulation of the cell cycle pathway from temporal up-pattern and the xenobiotics metabolism pathway from the temporal down-pattern.

Integration of multiple pathways

We found that only a few pathways (e.g., cell cycle pathway and xenobiotics metabolism pathway) were associated with microRNA target genes regulated by SST. However, as an individual gene, the microRNA target could play critical roles in diverse pathways. Therefore, we integrated all pathways that were significantly enriched by SST to identify key microRNA targets. From multiple pathways associated with the temporal up-pattern, the core microRNA targets selected were CCNA2, PTTG1, CDK1, CCNB2, CDC25B, CCL7, MAPK12 and ESPL1 (Figure 5A). From the temporal down-pattern, CYP2F2, CYP3A11, and CYP2C50 were selected as nodes with multiple roles (Figure 5B). The pathways containing these core targets of microRNA are shown below each network structure.

Figure 5
The core microRNA target genes in multiple pathways regulated by SST. (A) Nodes with high relative betweenness centrality were selected as the core microRNA targets in multiple pathways enriched in (A) the temporal up-pattern genes and (B) the temporal ...

TFBS analysis

The functional segregation of genes, based on the expression pattern, suggests that the gene transcription process would be the direct regulatory target of SST. Therefore, we investigated the possible association of the TFBS structure on the gene expression after SST treatment. By using the promoter region (-2000 bp to +500 bp from the transcription start site) of genes included in the temporal patterns, the correlation matrix of genes based on TFBS similarity was measured. The resultant clustering profile shows that genes in the temporal up-pattern are clearly distinguished from genes in the temporal down-pattern. As Figure 6A shows, two subgroups of genes were tightly clustered (i.e., Up-cluster and Down-cluster), which were primarily composed of genes from the temporal up-pattern and down-pattern, respectively. In addition to the main subgroups, there were other subgroups that also consisted exclusively of temporal up-pattern or down-pattern genes. The putative target genes of the microRNAs were interestingly also clustered into small subgroups (Figure 6A). This segregation of microRNA targets was more clearly observed in the temporal up-pattern genes (Figure 6B). One subgroup of microRNA targets was closely correlated with the similar TFBS structure (depicted as MicroRNA cluster in Figure 6B). MicroRNA target genes from the temporal down-pattern were also primarily concentrated on one cluster, although the number of target genes was small (Figure 6C). This separation of genes based on TFBS similarity indicates the presence of common cis-elements in the SST-regulated gene expression.

Figure 6
The clustering profile of temporally co-expressed genes by SST, based on the similarity of the transcription factor binding site (TFBS). Genes each from (A) both temporal up- and down-patterns, (B) from the temporal up-pattern, and (C) from the temporal ...


Despite the clinical usefulness of traditional herbal medicine, the complex nature of herbal chemical components prevents the elucidation of their exact molecular mechanisms. The herbal preparation of SST is also widely prescribed for the treatment of diverse liver diseases, but without clear understanding of its molecular mechanism [1]. What further complicates the situation is that SST is composed of seven different herbal plants (see Additional file 1: Table S1). Therefore, understanding the molecular activity of SST is limited when focusing on only a few major components or certain kinds of genes.

In the present study, we measured the global changes of genes and microRNAs expression induced by SST in cultured primary mouse hepatocytes, because the liver is a primary target organ of SST and is responsible for metabolizing xenobiotics. The expression profile shows two temporal expression patterns of genes after SST treatment, but no clear temporal pattern in microRNA expression (Figure 2). The microRNA expression levels after SST treatment were lower than the expression levels of genes. This suggests that a small number of microRNAs can regulate many genes. Therefore, it is critical to identify accurately the microRNA target genes. We used a correlation-based permutation approach to exclude possible false-positive links between microRNA and its putative target gene expression. The resultant 174 microRNA target genes were obtained from 463 temporal up-pattern genes and 19 targets were obtained from 177 temporal down-pattern genes (Figure 3 and Table 1). This indicates that microRNA is especially concentrated in the regulation of temporal up-pattern genes (p value < 0.001).

In addition to this unbalanced distribution of microRNA target genes, different biological functions were associated with microRNA targets in the two temporal patterns. For example, cell cycle pathway (mmu04110) was specifically involved in microRNA targets from the temporal up-pattern genes. On the other hand, non-microRNA target genes from temporal up-pattern genes were significantly associated (the FDR was less than 0.01) with cell signaling pathways such as the cytokine-cytokine receptor interaction pathway (mmu04060), the NF-κB signaling pathway (mmu04064), the MAPK signaling pathway (mmu04010), and the osteoclast differentiation pathway (mmu04380) (Table 2). Unlike the temporal up-pattern, microRNA targets from the temporal down-pattern were associated only with the xenobiotics metabolism pathway (mmu00980). Non-microRNA targets from the temporal down-pattern were involved in diverse pathways, among which were two primary categories: the immune-related pathway and the metabolism-related pathway. However, the number of microRNA targets from the temporal down-pattern was small. The SST-enriched peroxisome proliferator-activated receptor (PPAR) pathway is critical in regulating metabolism and proliferation by modulating E2F and AKT signaling in the liver regeneration process [39].

The temporal change of activity plot (Figure 4) indicated that many other pathways in addition to pathways listed in Table 2 were also activated or suppressed, reflecting the fact that diverse biological functions were influenced by the SST treatment. As expected, the cell cycle pathway (mmu04110) from the temporal up-pattern showed increased activity, whereas the immune-related pathways and drug metabolism pathways from the temporal down-pattern showed decreased activity. The regulatory role of SST on cell proliferation has interestingly been previously reported in studies indicating that SST has an antiproliferative effect on hepatocarcinoma primarily because of anticarcinogenic components such as baicalein, baicalin, and saikosaponin [2,40]. However, clinical evidence and recent reports also suggest that SST enhances liver function by promoting the regeneration of the liver in animal models [6,7]. Therefore, activation of cell cycle pathway and MAPK pathway in the present study could be explained by this liver-regenerative effect of SST.

Another major clinical effect of SST is immuno-modulatory activity in diverse diseases [41,42]. As evidenced in previous reports, SST can activate or repress immune processes, depending on the cell type and the clinical situation [9,43]. In our results, SST activated immune pathways such as the cytokine receptor pathway (mmu04060), the TNF signaling pathway (mmu04668), rheumatoid arthritis pathway (mmu05323), NOD-like receptor signaling pathway (mmu04621) but it also repressed other immune-related pathways such as the systemic lupus erythematosus pathway (mmu05322), the complement and coagulation pathway (mmu04610), and the Staphylococcus aureus infection pathway (mmu05150) (Figure 4).

This coordinated change, induced by SST on the activity of multiple pathways, implicates a common regulatory mechanism controlling the multiple pathways. We interestingly observed that some microRNA targets (e.g., CCNA2, PTTG1, CDK1, CCNB2, CDC25B, CCL7, MAPK12, and ESPL1 from the temporal up-pattern genes and CYP2F2, CYP3A11, and CYP2C50 from the temporal down-pattern genes) can act as core targets connected with multiple significant pathways from non-microRNA targets (Figure 5).

We mentioned in the previous paragraph that signal pathways regulated by SST (e.g., the cell cycle pathway, PPAR pathway, and MAPK pathway) could be associated with the liver regenerative activity of SST. This can be also confirmed by using individual core node genes. For example, CCNA2 and CCL7, main elements of cell cycle pathways and the cytokine receptor pathway, respectively, are associated with liver regeneration in the rat liver [44,45]. Also CDC25B can regulate mouse liver regeneration in association with FOXM1 by promoting hepatocyte proliferation [46,47]. CDK1, another key element in the cell cycle pathway, plays an essential role in the control of DNA replication in liver regeneration [48]. These previous reports suggest that core microRNA target genes in temporal up-pattern could be associated with the liver regeneration function of SST by enhancing cell proliferation function. On the other hand, core microRNA target genes in the temporal down-pattern (e.g., CYP2F2, CYP3A11, and CYP2C50) are exclusively associated with cytochrome P450 metabolism. However, there is interesting evidence that genes included in the cytochrome P450 family are also associated with liver regeneration. For example, early reduction of CYP activity has been observed in the regenerating rat liver, although the exact mechanism has not been elucidated [49]. The transcription of cytochrome P450 genes, including CYP3A11, moreover is reportedly suppressed by immune responses such as TNF-α in primary hepatocytes and hepatoma cells [50-52]. In consistent with the findings of previous reports, we observed the down-regulation of cytochrome P450 metabolism pathways and the activation of the cytokine pathway (mmu04060) and TNF signaling pathway (mmu04668) by SST (Table 2 and Figure 4), which imply the involvement of drug metabolism pathway and immune-pathways on liver regeneration process. To conclude, pathways identified in present study such as cell cycle pathway, drug metabolism-cytochrome P450 pathway and immune-related pathways, and individual core node genes could be possible molecular targets involved in liver regenerative process induced by SST. However, considering that SST has diverse pharmacological activities on various pathological conditions, the roles of these pathways and core node genes should be more precisely measured in a variety of physiological models.

We also observed that this coordinated regulation of gene expression by SST was predisposed in the genomic structure. As Figure 6A shows, the similarity in measurements of the TFBS clearly distinguished temporal up-pattern genes from temporal down-pattern genes. The present results imply that common cis-elements present in the promoter region of the genes could determine the temporal co-expression of genes induced by SST. Moreover, considering functions associated with each temporal pattern, the difference in TFBS structure between the two temporal patterns may be related to biological functions associated with each temporal pattern. For a clearer conclusion, a TFBS analysis should be performed of all genes at a genome level. It should also be elucidated whether resultant genes with a similar TFBS structure may be co-expressed by SST. What was more intriguing was that putative microRNA target genes also were clustered into separate subgroups, especially in the temporal up-pattern genes (Figure 6B). Recent research reveals that microRNA is involved in the promoter methylation of target genes to regulate the transcription level in association with transcription factors [53] and that this mechanism of gene expression would form the global regulatory network [12,54-56]; however, we do not know whether methylation-based regulation by microRNA is also involved in the present study. Moreover, there is no report on the role of the TFBS structure on the regulation of gene expression by microRNA. Therefore, we expect that our finding could give an important clue about the novel mechanism of gene expression by microRNA.


The present study is the first to indicate that SST systematically regulates gene expression by microRNA. We demonstrated that temporally up-regulated pattern by SST was associated with signaling pathways, including the cell cycle pathway, whereas the temporally down-regulated pattern included drug metabolism-related pathways and immune-related pathways, all of which could possibly contribute to the liver regenerative activity of SST. Also, this complex gene expression demonstrates that the effects of SST would be exerted from a delicately regulated mechanism on a genome-wide scale.


SST: Sho-saiko-to; HBSS: Hank’s balanced salt solution; STEM: Short time-series expression miner; FDR: False discovery rate; SPIA: Signaling pathway impact analysis; KEGG: Kyoto encyclopedia of genes and genomes; TFBS: Transcription factor binding site; UniPROBE: Universal PBM Resource for Oligonucleotide-Binding Evaluation; MAPK: Mitogen-activated protein kinase.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

KHS and YHK: conception, design of the experiment and preparation of the manuscript. BYK: conception, design of the experiment, analysis of data and preparation of the manuscript. All authors have read and approved the final manuscript.

Pre-publication history

The pre-publication history for this paper can be accessed here:

Supplementary Material

Additional file 1: Table S1:

Constituents of sho-saiko-to (SST). Table S2. List of genes in temporal patterns. Figure S1. The cytotoxic effect of Sho-saiko-to (SST) on primary hepatocytes. Hepatocytes are first cultured in 48-well plates at a density of 1.0 × 105 cells/well for 24 hours. After incubation, the cells are washed with phosphate-buffered saline and treated with different concentrations of SST (0.1–1.0 mg/mL) for 24 hours. Viability is measured in triplicate by using an in vitro colorimetric method (i.e., methyl thiazolyl tetrazoliym [MTT] assay). The viability is presented as the mean standard deviation (S.D.). Figure S2. Quantitative real-time polymerase chain reaction (Q-PCR). Mouse primary hepatocytes are treated with 500 μg/mL of SST at a density of 1.0 × 106 cells/60 mm dish for 1–24 hours in triplicate. The mRNA and microRNA are then reversetranscribed amplified, and detected by using Taqman probes (ABI, USA). The Q-PCR results are presented as the mean standard deviation (S.D.). Figure S3. Pathways enriched in the temporal up-pattern and temporal down-pattern. The position of each gene is denoted by red for the temporal up-pattern or blue for the temporal down-pattern in the pathways.


The authors would like to thank Dr. Hyeun Kyoo Shin (Basic Herbal Medicine Research Group, Korea Institute of Oriental Medicine, Republic of Korea) for supporting SST. This research was supported by the “study of high frequency TKM prescription based on microRNA (C13020)” funded by SME Partnership Center of Korea Institute of Oriental Medicine (KIOM) and in part by a grant (KIOM-2010-2) from the Inter-Institutional Collaboration Research Program under the Korea Research Council of Fundamental Science & Technology (KRCF).


  • Ohtake N, Nakai Y, Yamamoto M, Sakakibara I, Takeda S, Amagaya S, Aburada MJ, Chromatogr B. Separation and isolation methods for analysis of the active principles of Sho-saiko-to (SST) oriental medicine. Analyt Technol Biomed Life Sci. 2004;812:135–148. [PubMed]
  • Shimizu I. Sho-saiko-to: Japanese herbal medicine for protection against hepatic fibrosis and carcinoma. J Gastroenterol Hepatol. 2000;15(Suppl):D84–D90. [PubMed]
  • Morgan TR. Chemoprevention of hepatocellular carcinoma in chronic hepatitis C. Recent Results Cancer Res. 2011;188:85–99. [PubMed]
  • Shiota G, Maeta Y, Mukoyama T, Yanagidani A, Udagawa A, Oyama K, Yashima K, Kishimoto Y, Nakai Y, Miura T, Ito H, Murawaki Y, Kawasaki H. Effects of Sho-Saiko-to on hepatocarcinogenesis and 8-hydroxy-2′-deoxyguanosine formation. Hepatology. 2002;35:1125–1133. doi: 10.1053/jhep.2002.33066. [PubMed] [Cross Ref]
  • Chen MH, Chen JC, Tsai CC, Wang WC, Chang DC, Lin CC, Hsieh HY. Sho-saiko-to prevents liver fibrosis induced by bile duct ligation in rats. Am J Chin Med. 2004;32:195–207. doi: 10.1142/S0192415X04001862. [PubMed] [Cross Ref]
  • Miyamura M, Ono M, Kyotani S, Nishioka Y. Effects of sho-saiko-to extract on fibrosis and regeneration of the liver in rats. J Pharm Pharmacol. 1998;50:97–105. [PubMed]
  • Nakata R, Tsukamoto I, Miyoshi M, Kojo S. Liver regeneration after carbon tetrachloride intoxication in the rat. Biochem Pharmacol. 1985;34:586–588. doi: 10.1016/0006-2952(85)90195-9. [PubMed] [Cross Ref]
  • Ohtake N, Yamamoto M, Takeda S, Aburada M, Ishige A, Watanabe K, Inoue M. The herbal medicine Sho-saiko-to selectively inhibits CD8+ T-cell proliferation. Eur J Pharmacol. 2005;507:301–310. doi: 10.1016/j.ejphar.2004.11.037. [PubMed] [Cross Ref]
  • Chen MH, Chen JC, Tsai CC, Wang WC, Chang DC, Tu DG, Hsieh HY. The role of TGF-beta 1 and cytokines in the modulation of liver fibrosis by Sho-saiko-to in rat’s bile duct ligated model. J Ethnopharmacol. 2005;97:7–13. doi: 10.1016/j.jep.2004.09.040. [PubMed] [Cross Ref]
  • Sanoudou D, Mountzios G, Arvanitis DA, Pectasides D. Array-based pharmacogenomics of molecular-targeted therapies in oncology. Pharmacogenomics J. 2012;12:185–196. doi: 10.1038/tpj.2011.53. [PubMed] [Cross Ref]
  • Rukov JL, Shomron N. MicroRNA pharmacogenomics: post-transcriptional regulation of drug response. Trends Mol Med. 2011;17:412–423. doi: 10.1016/j.molmed.2011.04.003. [PubMed] [Cross Ref]
  • Chen K, Rajewsky N. The evolution of gene regulation by transcription factors and microRNAs. Nat Rev Genet. 2007;8:93–103. [PubMed]
  • Friedman RC, Farh KK, Burge CB, Bartel DP. Most mammalian mRNAs are conserved targets of microRNAs. Genome Res. 2009;19:92–105. [PubMed]
  • Creighton CJ, Hernandez-Herrera A, Jacobsen A, Levine DA, Mankoo P, Schultz N, Du Y, Zhang Y, Larsson E, Sheridan R, Xiao W, Spellman PT, Getz G, Wheeler DA, Perou CM, Gibbs RA, Sander C, Hayes DN, Gunaratne PH. Integrated analyses of microRNAs demonstrate their widespread influence on gene expression in high-grade serous ovarian carcinoma. Cancer genome atlas research network. PLoS One. 2012;7:e34546. doi: 10.1371/journal.pone.0034546. [PMC free article] [PubMed] [Cross Ref]
  • Genovesi LA, Carter KW, Gottardo NG, Giles KM, Dallas PB. Integrated analysis of miRNA and mRNA expression in childhood medulloblastoma compared with neural stem cells. PLoS One. 2011;6:e23935. doi: 10.1371/journal.pone.0023935. [PMC free article] [PubMed] [Cross Ref]
  • Cho JH, Gelinas R, Wang K, Etheridge A, Piper MG, Batte K, Dakhallah D, Price J, Bornman D, Zhang S, Marsh C, Galas D. Systems biology of interstitial lung diseases: integration of mRNA and microRNA expression changes. BMC Med Genomics. 2011;4:8. doi: 10.1186/1755-8794-4-8. [PMC free article] [PubMed] [Cross Ref]
  • Su Z, Xia J, Zhao Z. Functional complementation between transcriptional methylation regulation and post-transcriptional microRNA regulation in the human genome. BMC Genomics. 2011;12(Suppl 5):S15. doi: 10.1186/1471-2164-12-S5-S15. [PMC free article] [PubMed] [Cross Ref]
  • Wiench B, Chen YR, Paulsen M, Hamm R, Schröder S, Yang NS, Efferth T. Integration of different “-omics” technologies identifies inhibition of the IGF1R-Akt-mTOR signaling cascade involved in the cytotoxic effect of shikonin against leukemia cells. Evid Based Complement Alternat Med. 2013;2013:818709. [PMC free article] [PubMed]
  • Wen Z, Wang Z, Wang S, Ravula R, Yang L, Xu J, Wang C, Zuo Z, Chow MS, Shi L, Huang Y. Discovery of molecular mechanisms of traditional Chinese medicinal formula Si-Wu-Tang using gene expression microarray and connectivity map. PLoS One. 2011;6:e18278. doi: 10.1371/journal.pone.0018278. [PMC free article] [PubMed] [Cross Ref]
  • Kim BY, Lee J, Park SJ, Bang OS, Kim NS. Gene expression profile of the A549 human Non-small cell lung carcinoma cell line following treatment with the seeds of descurainia Sophia, a potential anticancer drug. Evid Based Complement Alternat Med. 2013;2013:584604. [PMC free article] [PubMed]
  • Kim BY, Cao LH, Kim JY. Common responses in gene expression by Ephedra herba in brain and heart of mouse. Phytother Res. 2011;25:1440–1446. doi: 10.1002/ptr.3434. [PubMed] [Cross Ref]
  • Shin IS, Lee MY, Kim Y, Seo CS, Kim JH, Shin HK. Subacute toxicity and stability of Soshiho-tang, a traditional herbal formula, in Sprague–Dawley rats. BMC Complement Altern Med. 2012;12:266. doi: 10.1186/1472-6882-12-266. [PMC free article] [PubMed] [Cross Ref]
  • Park H, Song KH, Jung PM, Kim JE, Ro H, Kim MY, Ma JY. Inhibitory effect of arctigenin from fructus arctii extract on melanin synthesis via repression of tyrosinase expression. Evid Based Complement Alternat Med. 2013;2013:965312. [PMC free article] [PubMed]
  • Song KH, Li T, Owsley E, Chiang JY. A putative role of micro RNA in regulation of cholesterol 7alpha-hydroxylase expression in human hepatocytes. J Lipid Res. 2010;51:2223–2233. doi: 10.1194/jlr.M004531. [PMC free article] [PubMed] [Cross Ref]
  • Bolstad BM, Irizarry RA, Astrand M, Speed TP. A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003;19:185–193. doi: 10.1093/bioinformatics/19.2.185. [PubMed] [Cross Ref]
  • Ernst J, Bar‒Joseph Z. STEM: a tool for the analysis of short time series gene expression data. BMC Bioinforma. 2006;7:191. doi: 10.1186/1471-2105-7-191. [PMC free article] [PubMed] [Cross Ref]
  • Wang X, el Naqa IM. Prediction of both conserved and nonconserved microRNA targets in animals. Bioinformatics. 2008;24:325–332. doi: 10.1093/bioinformatics/btm595. [PubMed] [Cross Ref]
  • Wang X. miRDB: a microRNA target prediction and functional annotation database with a wiki interface. RNA. 2008;14:1012–1017. doi: 10.1261/rna.965408. [PubMed] [Cross Ref]
  • Lee SY, Song KH, Koo I, Lee KH, Suh KS, Kim BY. Comparison of pathways associated with hepatitis B- and C-infected hepatocellular carcinoma using pathway-based class discrimination method. Genomics. 2012;99:347–354. doi: 10.1016/j.ygeno.2012.04.004. [PubMed] [Cross Ref]
  • Dennis G Jr, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, Lempicki RA. DAVID: Database for Annotation, Visualization, and Integrated Discovery. Genome Biol. 2003;4:P3. doi: 10.1186/gb-2003-4-5-p3. [PubMed] [Cross Ref]
  • Tarca AL, Draghici S, Khatri P, Hassan SS, Mittal P, Kim JS, Kim CJ, Kusanovic JP, Romero R. A novel signaling pathway impact analysis. Bioinformatics. 2009;25:75–82. doi: 10.1093/bioinformatics/btn577. [PMC free article] [PubMed] [Cross Ref]
  • Bang J, Jeon WK, Lee IS, Han JS, Kim BY. Biphasic functional regulation in hippocampus of Rat with chronic cerebral hypoperfusion induced by permanent occlusion of bilateral common carotid artery. PLoS One. 2013;8:e70093. doi: 10.1371/journal.pone.0070093. [PMC free article] [PubMed] [Cross Ref]
  • Zhang JD, Wiemann S. KEGGgraph: a graph approach to KEGG PATHWAY in R and bioconductor. Bioinformatics. 2009;25:1470–1471. doi: 10.1093/bioinformatics/btp167. [PMC free article] [PubMed] [Cross Ref]
  • Shannon P. MotifDb: an annotated collection of protein-DNA binding sequence motifs. R package version 1.2.2. []
  • Portales-Casamar E, Thongjuea S, Kwon AT, Arenillas D, Zhao X, Valen E, Yusuf D, Lenhard B, Wasserman WW, Sandelin A. JASPAR 2010: the greatly expanded open-access database of transcription factor binding profiles. Nucleic Acids Res. 2010;38:D105–D110. doi: 10.1093/nar/gkp950. [PMC free article] [PubMed] [Cross Ref]
  • Wasserman WW, Sandelin A. Applied bioinformatics for the identification of regulatory elements. Nat Rev Genet. 2004;5:276–587. doi: 10.1038/nrg1315. [PubMed] [Cross Ref]
  • Newburger DE, Bulyk ML. UniPROBE: an online database of protein binding microarray data on protein-DNA interactions. Nucleic Acids Res. 2009;37:D77–D82. doi: 10.1093/nar/gkn660. [PMC free article] [PubMed] [Cross Ref]
  • Veerla S, Hoglund M. Analysis of promoter regions of coexpressed genes identified by microarray analysis. BMC Bioinforma. 2006;7:384. doi: 10.1186/1471-2105-7-384. [PMC free article] [PubMed] [Cross Ref]
  • Liu HX, Fang Y, Hu Y, Gonzalez FJ, Fang J, Wan YJ. PPARβ Regulates Liver Regeneration by Modulating Akt and E2f Signaling. PLoS One. 2013;8:e65644. doi: 10.1371/journal.pone.0065644. [PMC free article] [PubMed] [Cross Ref]
  • Inoue T, Jackson EK. Strong antiproliferative effects of baicalein in cultured rat hepatic stellate cells. Eur J Pharmacol. 1999;378:129–135. doi: 10.1016/S0014-2999(99)00418-5. [PubMed] [Cross Ref]
  • Borchers AT, Sakai S, Henderson GL, Harkey MR, Keen CL, Stern JS, Terasawa K, Gershwin ME. Shosaiko-to and other Kampo (Japanese herbal) medicines: a review of their immunomodulatory activities. J Ethnopharmacol. 2000;73:1–13. doi: 10.1016/S0378-8741(00)00334-2. [PubMed] [Cross Ref]
  • Ohtake N, Nakai Y, Yamamoto M, Ishige A, Sasaki H, Fukuda K, Hayashi S, Hayakawa S. The herbal medicine Shosaiko-to exerts different modulating effects on lung local immune responses among mouse strains. Int Immunopharmacol. 2002;2:357–366. doi: 10.1016/S1567-5769(01)00161-8. [PubMed] [Cross Ref]
  • Kang H, Choi TW, Ahn KS, Lee JY, Ham IH, Choi HY, Shim ES, Sohn NW. Upregulation of interferon-gamma and interleukin-4, Th cell-derived cytokines by So-Shi-Ho-Tang (Sho-Saiko-To) occurs at the level of antigen presenting cells, but not CD4 T cells. J Ethnopharmacol. 2009;123:6–14. doi: 10.1016/j.jep.2009.02.045. [PubMed] [Cross Ref]
  • Wang WB, Fan JM, Zhang XL, Xu J, Yao W. Serial expression analysis of liver regeneration-related genes in rat regenerating liver. Mol Biotechnol. 2009;43:221–231. doi: 10.1007/s12033-009-9199-z. [PubMed] [Cross Ref]
  • Chen X, Xu C, Zhang F, Ma J. Comparative analysis of expression profiles of chemokines, chemokine receptors, and components of signaling pathways mediated by chemokines in eight cell types during rat liver regeneration. Genome. 2010;53:608–618. doi: 10.1139/G10-040. [PubMed] [Cross Ref]
  • Wang X, Kiyokawa H, Dennewitz MB, Costa RH. The Forkhead Box m1b transcription factor is essential for hepatocyte DNA replication and mitosis during mouse liver regeneration. Proc Natl Acad Sci U S A. 2002;99:16881–16886. doi: 10.1073/pnas.252570299. [PubMed] [Cross Ref]
  • Wang X, Krupczak-Hollis K, Tan Y, Dennewitz MB, Adami GR, Costa RH. Increased hepatic Forkhead Box M1B (FoxM1B) levels in old-aged mice stimulated liver regeneration through diminished p27Kip1 protein levels and increased Cdc25B expression. J Biol Chem. 2002;277:44310–44316. doi: 10.1074/jbc.M207510200. [PubMed] [Cross Ref]
  • Garnier D, Loyer P, Ribault C, Guguen-Guillouzo C, Corlu A. Cyclin-dependent kinase 1 plays a critical role in DNA replication control during rat liver regeneration. Hepatology. 2009;50:1946–1956. doi: 10.1002/hep.23225. [PubMed] [Cross Ref]
  • Starkel P, Laurent S, Petit M, van den Berge V, Lambotte L, Horsmans Y. Early down-regulation of cytochrome P450 3A and 2E1 in the regenerating rat liver is not related to the loss of liver mass or the process of cellular proliferation. Liver. 2000;20:405–410. doi: 10.1034/j.1600-0676.2000.020005405.x. [PubMed] [Cross Ref]
  • Morgan ET, Goralski KB, Piquette-Miller M, Renton KW, Robertson GR, Chaluvadi MR, Charles KA, Clarke SJ, Kacevska M, Liddle C, Richardson TA, Sharma R, Sinal CJ. Regulation of drug-metabolizing enzymes and transporters in infection, inflammation, and cancer. Drug Metab Dispos. 2008;36:205–216. doi: 10.1124/dmd.107.018747. [PubMed] [Cross Ref]
  • Aitken AE, Richardson TA, Morgan ET. Regulation of drug-metabolizing enzymes and transporters in inflammation. Annu Rev Pharmacol Toxicol. 2006;46:123–149. doi: 10.1146/annurev.pharmtox.46.120604.141059. [PubMed] [Cross Ref]
  • Morgan ET. Regulation of cytochromes P450 during inflammation and infection. Drug Metab Rev. 1997;29:1129–1188. doi: 10.3109/03602539709002246. [PubMed] [Cross Ref]
  • Taguchi YH. MicroRNA-mediated regulation of target genes in several brain regions is correlated to both microRNA-targeting-specific promoter methylation and differential microRNA expression. Bio Data Min. 2013;6:11. [PMC free article] [PubMed]
  • Shalgi R, Lieber D, Oren M, Pilpel Y. Global and local architecture of the mammalian microRNA-transcription factor regulatory network. PLoS Comput Biol. 2007;6:e131. [PubMed]
  • Hobert O. Gene regulation by transcription factors and microRNAs. Science. 2008;6:1785–1786. [PubMed]
  • Zhou Y, Ferguson J, Chang JT, Kluger Y. Inter- and intra-combinatorial regulation by transcription factors and microRNAs. BMC Genomics. 2007;6:396. [PMC free article] [PubMed]

Articles from BMC Complementary and Alternative Medicine are provided here courtesy of BioMed Central