- Research article
- Open Access
- Open Peer Review
Study of ZHENG differentiation in Hepatitis B-caused cirrhosis: a transcriptional profiling analysis
BMC Complementary and Alternative Medicinevolume 14, Article number: 371 (2014)
In traditional Chinese medicine (TCM) clinical practice, ZHENG (also known as TCM syndrome) helps to understand the human homeostasis and guide individualized treatment. However, the scientific basis of ZHENG remains unclear due to limitations of current reductionist approaches.
We collected the leukocyte samples of three hepatitis B-caused cirrhosis (HBC) patients with dampness-heat accumulation syndrome (DHAS) and three HBC patients with liver depression and spleen deficiency syndrome (LDSDS) for microarray analysis. We generated Gene-Regulatory-Networks (GeneRelNet) from the differentially expressed genes (DEGs) of microarray date. Core genes were validated using anther independent cohort of 40 HBC patients (20 DHAS, 20 LDSDS) with RT-PCR.
There were 2457 mapped genes were differentially expressed between DHAS and LDSDS (Fold change ≥ 2.0, P < 0.05). There were markedly different genes co-expression patterns in DHAS and LDSDS. Furthermore, three differential co-expression genes including purine nucleoside phosphorylase (PNP); aquaporin 7 (AQP7) and proteasome 26S subunit, non-ATPase 2 (PSMD2) were screened by GeneRelNets, and their mRNA expressions were further validated by real time RT-PCR. The results were consistent with microarray. The PNP (P = 0.007), AQP7 (P = 0.038) and PSMD2 (P = 0.009) mRNA expression is significant difference between DHAS and LDSDS using the non-parametric test. Furthermore, we constructed an mRNA panel of PNP, AQP7 and PSMD2 (PAP panel) by logistic regression model, and evaluated the PAP panel to distinguish DHAS from LDSDS by area under the receiver operating characteristic curve (AUC) analysis, which showed a higher accuracy (AUC = 0.835). Gene ontology (GO) analysis indicated that the DHAS is most likely related to system process while the functions overrepresented by LDSDS most related to the response to stimulus.
This study suggested that there are particular transcriptional profiles, genes co-expressions patterns and functional properties of DHAS and LDSDS, and PNP, AQP7, and PSMD2 may be involved in ZHENG differentiation of DHAS and LDSDS in HBC.
Chronic hepatitis B virus (HBV) infection is a serious health problem because its worldwide distribution and it has potential adverse sequelae such as cirrhosis and hepatocellular carcinoma (HCC) [1, 2]. It was estimated that more than 200,000 chronic HBV carriers die of liver cirrhosis each year globally . Since HBV replication may persist after the development of chronic hepatitis B (CHB) , hepatitis B- caused cirrhosis (HBC) is easily happen. Moreover, HBC increases the risk of developing HCC by at least 40 times over the risk of an average person .
Traditional Chinese medicine (TCM) is extensively used for treatment of liver disease in China, with the great advantages in early intervention, combination therapies, personalized medicine, etc. TCM treatment is based on ZHENG (also known as TCM syndrome) differentiation, namely, discerning patterns of imbalances within the body and between the body and the environment by the analysis of symptoms and signs of patients. It is the key principle of TCM. Since ZHENG differentiation depended on clinical observation and TCM practitioner’s experience, it would be subjective and unrepeatable. The phenotype-oriented diagnosis and therapy in TCM needs to be objective and quantified evaluation for the inheritance .
TCM has long been practiced as an empirical and holistic medicine, which happens to share the same concept with the systems biological medicine. The advent of high-throughput technologies, such as transcript microarray chips, has allowed simultaneous interrogation of multiple molecular components at any given time. Advances in gene expression profiling have been used to identify key differentially expressed genes (DEGs) and pathways in HBV-related patients with different ZHENGs. In the previous study, we reported that the molecular mechanism of “same TCM syndrome for different diseases and different TCM syndrome for same disease” . Besides, we integrated the CHB gene expression with topological features of different ZHENGs resulted in enhanced diagnosis of CHB. The miRNA-target network analysis elucidated dysfunctional interactions among genes and miRNAs-target genes of three different ZHENGs . Above all, the transcriptional profiles of ZHENG differentiation in HBC could be interpreted using the systems-omics approach.
The dampness-heat accumulation syndrome (DHAS) and the liver depression and spleen deficiency syndrome (LDSDS) are the major ZHENGs in HBC . In this study, we demonstrated that the transcriptional profiles and the functional properties are significantly different in DHAS and LDSDS. Besides, the differential expression or co-expression genes could be the potential biomarkers for ZHENG differentiation with DHAS and LDSDS in HBC patients.
Samples and RNA extraction
Morning fasting venous blood samples of total 46 HBC patients (23 DHAS, 23 LDSDS) were obtained from Shuguang Hospital and Longhua Hospital in Shanghai, China. Six HBC patients (3 DHAS, 3 LDSDS) were used for microarray analysis. Another independent cohort of 40 HBC patients (20 DHAS, 20 LDSDS) was used as validation group. The research protocol was approved by the respective institutional review boards, and informed consent was obtained for all study participants. The study was approved by the Official Ethics Committee of the Shanghai University of TCM and written informed consent to participate in the study was obtained from all subjects included. Diagnosis standard of cirrhosis is referred to “Chronic hepatitis B prevention and treatment guidelines” . The TCM ZHENG types were identified by three chief or deputy physicians, according to “evaluation criteria of the clinical diagnosis, drug efficacy, and ZHENG differentiation for cirrhosis (pilot program)” . All the outpatients were diagnosed by attending TCM physicians at the first time and then identified by three chief TCM physicians, who were consistently diagnosed as DHAS or LDSDS by all of the physicians were enrolled in this study [10, 11]. The patients with other hepatotropic virus hepatitis, chronic severe hepatitis, serious primary disease or pregnancy were excluded. The exclusion criterion of HBC: (1) cases complicated with other hepatotropic virus; (2) age < 18 or age > 65; (3) associated with serious primary disease of heart, kidney, lung, endocrine, blood, metabolic and gastrointestinal; (4) psychotic patients; (5) pregnant or lactating women. Samples of three DHAS, three LDSDS and three normal control were used to microarray detection. Another independent twenty DHAS and twenty LDSDS in HBC patients were used to verify the accuracy of the results. The characteristics of the study participants are presented in Table 1. The leukocyte from blood samples were isolated with Ficoll optimized density gradient separation  and saved at -80°C.
Total RNA was extracted using a “two-step” protocol as described previously . Total RNA of leukocyte from the whole blood was extracted using TRIzol® Reagent (Invitrogen, Carlsbad, CA, US) and saved at -80°C. A Quality Control was carried out with NanoDrop ND-1000.
Microarrays and data analysis
The biotinylated cDNAs were hybridized to NimbleGen Homo sapiens 12 × 135K Array (Roche, CAT No. A6484-00-01). This product analyzes the expression level for over 135,000 transcripts and variants, including 44049 well-substantiated human genes.
Microarray data preprocessing was performed using the Gene-Pix software. Raw expression data were log2 transformed and normalized by quantile normalization. Probes were considered robustly expressed if Signal/Noise (SNR) > 2.
The t-test function in R software was used to select difference expressed gene (threshold: P value < 0.05) in ZHENGs of DHAS and LDSDS. They were both compared with normal control. Heatmap analysis which executed in Cluster 3.0 was computing the hierarchical clustering in both rows and columns according to the set of gene values and drawing a color image as a visible result. Gene Ontology (GO) was analyzed using a DAVID online analysis tool (http://david.abcc.ncifcrf.gov).
We used a stepwise logistic regression model to combine diagnostic mRNA markers based on the data obtained by qRT-PCR. The predicted probability of distinguishing of DHAS and LDSDS was used as a surrogate marker to construct a receiver operating characteristic (ROC) curve. The AUC was used as an accuracy index for evaluating the diagnostic performance of the selected mRNA panel. All tests were two-tailed and P < 0.05 was considered statistically significant.
Gene-Regulatory-Networks (GeneRelNet) construction
GeneRelNets  were presented to find the interactions among genes, indicating the gene co-expression relationship. GeneRelNets was constructed by DEGs that obtained by comparing the DHAS or LDSDS. They were built according to the normalized signal intensity of specific expression genes. For each pair of genes, we calculate the Pearson correlation and choose the significant correlation pairs to construct the network.
In network analysis, Degree centrality, Clustering coefficient and Core genes were analyzed. Degree centrality, defined as the link numbers one node has to the other, is the most simplest and important measures of the centrality of a gene within a network that determine the relative importance . Clustering coefficient is a measure of degree to which nodes in a graph tend to cluster together. The higher clustering coefficient of an individual gene, the higher is the affection among its neighbor genes . The purpose of network structure analysis of each ZHENG is to locate differential co-expression genes, which connect most adjacent genes and have biggest degrees with high Clustering coefficient. While considering different networks, the Core genes were determined by the degree differences between two class samples . They always own the biggest degree differences and selected from the located differential co-expression genes.
The cDNAs were synthesized by the Invitrogen First-Strand cDNA Synthesis kits (Invitrogen, Carlsbad, CA, US). Samples will be qualified for further processing if the A260/A280 spectrophotometeric ratio is between 1.8 and 2.1. One μg of total RNA was transcribed into cDNA in a 20 μL reaction volume.
Real-time RT-PCR was used to verify the differential expression of three genes that detected by the NimbleGen GeneChip. The primers used are listed in Table 2. Each real-time RT-PCR reaction (in 25 μL) contained 2 × SYBR Green Real-time RT-PCR Master Mix, 0.4 μM primers, and 0.5 μL of template cDNA. The cycling conditions consisted of an initial, single cycle of 5 min at 95°C, followed by 40 cycles of 30 sec at 95°C, 30 sec at 54°C, 15 sec at 72°C, and fluorescence acquisition at 83°C for 1 sec. The cDNA was synthesized using reverse transcriptase (TOYOBO, Osaka, Japan), oligo (dT) and random primers with 5 μg RNA from the same samples as those used in the microarray. The PCR amplifications were performed in duplicate for each sample. The gene expression levels were quantified relative to the expression of β-actin by employing an optimized comparative Ct (ΔΔCt) value method. The differences in gene expression levels between groups were compared using the Student's t-test. A p value < 0.05 was considered significant.
Identification of core genes of DHAS and LDSDS in HBC
We analyze the DEGs of DHAS and LDSDS from the microarrays data. Both samples from DHAS and LDSDS patients were compared with normal controls and then compared to each other for the array analysis. DEGs were defined as |foldchange | ≥ 2, p < 0.05. 3133 genes were identified after comparison of DHAS and normal control. 3627 genes were identified after comparison of LDSDS and normal control. Then the DHAS group compared with LDSDS group. After the two steps of DEGs selection, a total of 2457 mapped genes were identified to be differentially expressed in DHAS compared to LDSDS with |foldchange (DHAS/LDSDS) | ≥ 2, p < 0.05. (Figure 1), among which 1855 genes were up-regulated and 602 genes were down-regulated.
Based on the DEGs of DHAS and LDSDS (Additional file 1), GeneRelNets were constructed, respectively (Figure 2). The structure of GeneRelNet reflects the situation of genes co-expressions, which is interaction among genes. There were markedly different genes co-expression patterns in DHAS and LDSDS (Figure 2). Furthermore, the genes co-expressions analysis was used to screen the candidate genes that may involve in the differentiation of DHAS and LDSDS.
In network analysis, the differential co-expressed genes with higher Degree Centrality (degree ≥ 10) and higher Clustering Coefficient (Clustering Coefficient ≥ 0. 2) were considered to have important regulation and control ability. They were located from GeneRelNet of DHAS or LDSDS for ZHENG differentiation. Core genes in DHAS and LDSDS had higher degree differences (ΔDegree) (Figure 3). In this study, the differential co-expression genes were defined with ΔDegree ≥ 8. Core genes were PNP (ΔDegree, 16), AQP7 (ΔDegree, 9), NFE2 (ΔDegree, 9) and PSMD2 (ΔDegree, 8) (Table 3).
Since PNP and AQP7 were the top differential co-expression genes, and PSMD2 had the highest Degree in both DHAS and LDSDS. They were chosen to verify the transcript expressions and compare the difference between DHAS and LDSDS in HBC.
Validation of core genes in DHAS and LDSDS GeneRelNets
Following the computational analyses of GeneRelNets, quantitative RT-PCR was carried out using a cohort of independent samples of DHAS (n = 20) and LDSDS (n = 20) in HBC. The expression result of PNP, AQP7 and PSMD2 was consistent with the results from microarrays. At the transcript expression level, statistical significance was noted for AQP7 and PSMD2 (Figure 4). PNP transcript levels were elevated significantly (P = 0.007) in DHAS group compared with LDSDS group. AQP7 (P = 0.038) and PSMD2 (P = 0.009) expression levels were significantly elevated in DHAS group compared to LDSDS group (Figure 4).
Sensitivity and specificity of core genes of DHAS and LDSDS for ZHENG differentiation
The levels of PNP, AQP7 and PSMD2 mRNA expressions between DHAS and LDSDS were compared. As shown in Figure 5, AUC of PNP, AQP7 and PSMD2 were 0.566 (P = 0.473, 95% CI, 0.372 to 0.761), 0.702 (P = 0.234, 95% CI, 0.426 to 0.794) and 0.767(P = 0.008, 95% CI, 0.587 to 0.900) respectively (Figure 5A).
Furthermore, we applied a stepwise logistic regression model to combine the three core genes to distinguishing DHAS from LDSDS. The predicted probability from the logit model based on the three mRNA panel (PAP panel). The logit model (P = LDSDS) = -0.634 + 13.701 PNP-230.17 AQP7-512.468 PSMD2 was used to construct the ROC-curve. The diagnostic performance for the established mRNA panel was evaluated by using ROC analysis. The AUC for the PAP panel was 0.835 (P = 0.000, 95% CI, 0.708 to 0.962, Figure 5B).
Functional properties of DEGs between DHAS and LDSDS in HBC
To further determine differential co-expression genes in DHAS and LDSDS involved in the cellular behavior, we conducted a GO enrichment analysis and searched for functional categories.The results showed that these differential co-expression genes were enriched for neurological system process and sensory perception in DHAS. By contrast, the functions overrepresented by LDSDS samples related to the sensory perception of smell, sensory perception of chemical stimulus, response to stimulus, neurological system process, and cognition, etc. Illustrative gene function diagrams corresponding to DHAS and LDSDS are given (Figure 6), respectively.
TCM ZHENG differs significantly in HBC patients . DHAS and LDSDS are two typical TCM ZHENGs with different clinical phenotype in HBC . In TCM, DHAS is recognized to dampness-heat accumulation. It represents accumulation of dampness-heat in the liver and gallbladder resulting in impaired bile flow and downward pouring of dampness-heat; LDSDS is recognized to the dampness and heat accumulated in liver and gallbladder. It represents a syndrome marked by hypochondriac and abdominal painful distension, depressed mood, frequent sighing, anorexia, uncomfortable loose bowels or alleviation of abdominal pain after defecation, borborygmi with flatus and white slimy tongue coating, the same as the syndrome of liver stagnation and spleen deficiency, also known as the liver-spleen disharmony syndrome . It was the first time to stratify the DHAS and LDSDS of HBC patients with the transcriptional technology.
Elucidating which genetic factors contribute to the ZHENGs differentiation may have important therapeutic implications in TCM practice. High-throughput, genome-wide analytical technologies, especial expression microarrays could be both accurate and precise when properly implemented. It has been used for disease classification and identification of causal mechanisms in different human liver diseases, such as alcoholic liver disease , fatty liver disease , HCV-associated liver disease , liver cancer , and even liver tumor metastases . However, it is lack of microarray-based diagnostic algorithm for ZHENG differentiation in HBC.
Since TCM is holistic and systems medicine , TCM ZHENG differentiation may be demonstrated through the investigation in gene regulated network. In this study, we identified the whole gene expressions of DHAS and LDSDS in HBC, and we found the significant difference in gene expression patterns between DHAS and LDSDS (Figure 1). Moreover, we propose using co-expression-based gene networks (GeneRelNets) to identify the biomarkers that are associated with DHAS and LDSDS syndromes. In GeneRelNets, the links between genes (nodes) were determined by the extent of their correlated pattern of expression across multiple genes and a variety of biological relationships. Since the DHAS and LDSDS represent the different ZHENGs of HBC, the GeneRelNets analysis showed the markedly different genes co-expressions and their patterns (Figure 2), and there were differential co-expression genes (Figure 3). The degree centrality difference of two GeneRelNets in DHAS and LDSDS indicated a possibility of distinguishing DHAS and LDSDS.
The differential co-expression genes of PNP, AQP7 and PSMD2 between DHAS and LDSDS were verified (Figure 4) in another independent patient cohorts. To observe the sensitivity and specificity of these mRNAs, ROC curve analysis was conducted to differentiating DHAS and LDSDS. The AUC of the three genes was 0.566, 0.372 and 0.702, respectively (Figure 5-A). The mRNA panel with the three mRNAs from the multivariate logistic regression model demonstrated high accuracy (AUC = 0.835) in the differentiation of DHAS and LDSDS (Figure 5-B). The result indicated that the three genes may be involved in DHAS and LDSDS differentiation in HBC. PNP, a key enzyme in the purine salvage pathway, is located primarily in the cytoplasm of endothelial cells, Kupffer cells, and hepatocytes and is released into hepatic sinusoids during necrosis  there are studies showed the value of PNP as a marker of liver damage in rodents . AQP7 are transmembrane proteins that belong to the subset of aquaglyceroporins. The lack of AQP7 in both cases resulted in an increased accumulation of glycerol within the adipocytes leading to adipocyte hypertrophy . A large-scale screening study  has been reported that PSMD2 offers a growth advantage to NIH3T3 cells under certain conditions, which was known as a collagen-producing fibroblast cell. Therefore, PSMD2 might be involved in the generation of liver fibrosis. Epidemiological studies [29, 30] have reported that HBC develop to HCC often following with TCM ZHENG changes from excess to deficiency. In this study, AQP7 and PSMD2 expression in LDSDS was significantly higher than that in DHAS. We speculated that both HBC prognosis and HCC morbidity in LDSDS patients may be more than that in DHAS patients. A future study to elucidate how such molecular consequences are mediated in PNP, AQP7 and PSMD2-involed relationship between HBC and TCM ZHENG is warranted. It indicated a possibility of distinguishing DHAS and LDSDS by the expressions of PNP, AQP7 and PSMD2 their combination.
Furthermore, we performed GO function enrichment analysis to understand the different functions of DHAS and LDSDS. The results indicated that there are similarities of gene functions in DHAS and LDSDS, they all mainly involved in system development, anatomical structure development and so on. Otherwise, the DHAS was most likely related to cellular process when the functions overrepresented by LDSDS samples most related to the response to stimulus (Figure 6). Previous studies [31, 32] showed that the central neurobiological mechanism of LDSDS closely correlates to the hypothalamic-pituitary-adrenal axis, brain-gut axis, myriad central neurotropic factors, neurotransmitters and neuropeptides involving in many encephalic regions such as the hypothalamus, hippocampus, cortex, amygdale, etc. This funding was consisting with our GO analysis, which related to the chronic stress.
The limitations of this study are as follows. (1) A small sample size limited the value and generalization of the results; a larger population study should be further considered in both test group and validation group. (2) Only two typical ZHENGs in HBC were discussed in this study. Other ZHENGs types of HBC required for further research. (3) Since the gene regulation network is complicate and remaining unclear in HBV-caused chronic diseases, the mechanism of selected genes needs further study.
In this study, we demonstrated the transcriptional profiles, differential co-expression genes and their functional properties for ZHENG differentiation in HBC patients. The results suggested that there are particular transcriptional profiles, genes co-expressions patterns and their functional properties of DHAS and LDSDS, and several interesting genes including PNP, AQP7, PSMD2 and their combination may be involved in ZHENG differentiation of DHAS and LDSDS in HBC. Although further investigation with larger number of samples and other type ZHENGs is requested, these results supported the hypothesis that transcriptional profiles may be used to the ZHENG differentiation in HBC.
traditional Chinese medicine
hepatitis B-caused cirrhosis
dampness-heat accumulation syndrome
liver depression and spleen deficiency syndrome
differentially expressed genes
purine nucleoside phosphorylase
proteasome 26S subunit, non-ATPase 2
Chronic hepatitis B virus
chronic hepatitis B
hepatitis B-caused cirrhosis
Alpha Fetal Protein
gamma glutamyl transferase
Lavanchy D: Hepatitis B virus epidemiology, disease burden, treatment, and current and emerging prevention and control measures. J Viral Hepat. 2004, 11: 97-107. 10.1046/j.1365-2893.2003.00487.x.
Liaw YF, Chu CM: Hepatitis B virus infection. Lancet. 2009, 373: 582-592. 10.1016/S0140-6736(09)60207-5.
Perz JF, Armstrong GL, Farrington LA, Hutin YJ, Bell BP: The contributions of hepatitis B virus and hepatitis C virus infections to cirrhosis and primary liver cancer worldwide. J Hepatol. 2006, 45: 529-538. 10.1016/j.jhep.2006.05.013.
Schütte K, Bornschein J, Malfertheiner P: Hepatocellular carcinoma–epidemiological trends and risk factors. Dig Dis. 2009, 27: 80-92.
Su SB, Lu A, Li S, Jia W: Evidence-based ZHENG: a traditional Chinese medicine syndrome. Evid base Compl Alternative Med. 2012, 2012: 246538-http://www.hindawi.com/journals/ecam/2012/246538/,
Guo ZZ, Yu SH, Guan Y, Li YY, Lu YY, Zhang H, Su SB: Molecular mechanisms of same TCM syndrome for different diseases and different TCM syndrome for same disease in chronic hepatitis B and liver cirrhosis. Evid base Compl Alternative Med. 2012, 2012: 120350-http://www.hindawi.com/journals/ecam/2012/120350/,
Chen QL, Lu YY, Zhang GB, Song YN, Zhou Q-M, Zhang H, Zhang W, Su S-B: Progression from excessive to deficient syndromes in chronic Hepatitis B: a dynamical network analysis of miRNA array data. Evid base Compl Alternative Med. 2013, 2013: 120350-
Zhang YX, Wei BH: Standard of clinic diagnosis, syndrome differentiation and assessing curative effect on hepatocirrhosis. Chin J Integrated Trad West Med Gastro-spleen. 1994, 14: 237-238.
Association CM: Chronic hepatitis B prevention and treatment guidelines. Chinese J Infect Dis. 2001, 19: 56-62.
Zhao Y, Gou XJ, Dai JY, Peng JH, Feng Q, Sun SJ, Cao HJ, Zheng NN, Fang JW, Jiang J, Su SB, Liu P, Hu YY, Zhang YY: Differences in metabolites of different tongue coatings in patients with chronic hepatitis B. Evid base Compl Alternative Med. 2013, 2013: 204908-
Sun SJ, Dai JY, Fang JW, Gou XJ, Cao HJ, Zheng NN, Wang Y, Zhang W, Zhang YY, Jia W: Differences of excess and deficiency Zheng in patients with chronic hepatitis B by urinary metabonomics. Evid base Compl Alternative Med. 2013, 2013: 738245-http://www.hindawi.com/journals/ecam/2013/738245/,
Gertler R, Rosenberg R, Fuehrer K, Dahm M, Nekarda H, Siewert JR: Detection of circulating tumor cells in blood using an optimized density gradient centrifugation. Molecular Staging of Cancer. 2003, Springer, 149-155.
Pujana MA, Han JDJ, Starita LM, Stevens KN, Tewari M, Ahn JS, Rennert G, Moreno V, Kirchhoff T, Gold B: Network modeling links breast cancer susceptibility and centrosome dysfunction. Nat Genet. 2007, 39: 1338-1349. 10.1038/ng.2007.2.
Barabási AL, Oltvai ZN: Network biology: understanding the cell's functional organization. Nat Rev Genet. 2004, 5: 101-113. 10.1038/nrg1272.
Ravasz E, Somera AL, Mongru DA, Oltvai ZN, Barabási A-L: Hierarchical organization of modularity in metabolic networks. Science. 2002, 297: 1551-1555. 10.1126/science.1073374.
Carlson MR, Zhang B, Fang ZX, Mischel PS, Horvath S, Nelson SF: Gene connectivity, function, and sequence conservation: predictions from modular yeast co-expression networks. BMC Genomics. 2006, 7: 40-10.1186/1471-2164-7-40.
Zhu LL, Meng H, Jiang J: Study on TCM syndrome typing of chronic hepatitis B. Chin J Integrated Tradit West Med. 2008, 28: 20-
Organization WH: WHO international standard terminologies on traditional medicine in the western pacific region. Book WHO international standard terminologies on traditional medicine in the western pacific region (Editor ed.^eds.). 2007, City: City
Deaciuc IV, Song ZY, McClain CJ: Lessons from large-scale gene profiling of the liver in alcoholic liver disease. Hepatol Res. 2005, 31: 187-192. 10.1016/j.hepres.2005.02.004.
Hoekstra M, Li Z, Kruijt JK, Eck MV, Berkel TJV, Kuiper J: The expression level of non-alcoholic fatty liver disease-related gene PNPLA3 in hepatocytes is highly influenced by hepatic lipid status. J Hepatol. 2010, 52: 244-251. 10.1016/j.jhep.2009.11.004.
Branch AD, Walewski JL: The coming impact of gene expression profiling on the diagnosis and treatment of HCV-associated liver disease. Antiviral Res. 2001, 52: 173-179. 10.1016/S0166-3542(01)00182-6.
Kurokawa Y, Matoba R, Takemasa I, Nakamori S, Tsujie M, Nagano H, Dono K, Umeshita K, Sakon M, Ueno N: Molecular features of non-B, non-C hepatocellular carcinoma: a PCR-array gene expression profiling study. J Hepatol. 2003, 39: 1004-1012. 10.1016/S0168-8278(03)00473-2.
Trolet J, Hupé P, Huon I, Lebigot I, Decraene C, Delattre O, Sastre-Garau X, Saule S, Thiéry J-P, Plancher C: Genomic profiling and identification of high-risk uveal melanoma by array CGH analysis of primary tumors and liver metastases. Invest Ophthalmol Vis Sci. 2009, 50: 2572-2580. 10.1167/iovs.08-2296.
Buriani A, Garcia Bermejo ML, Bosisio E, Xu Q, Li H, Dong X, Simmonds MSJ, Carrara M, Tejedor N, Lucio-Cazana J, Hylands PJ: Omic techniques in systems biology approaches to traditional Chinese medicine research: Present and future. J Ethnopharmacol. 2012, 140: 535-544. 10.1016/j.jep.2012.01.055.
Ohuchi T, Tada K, Akamatsu K: Endogenous ET-1 contributes to liver injury induced by galactosamine and endotoxin in isolated perfused rat liver. Am J Physiol Gastrointest Liver Physiol. 1995, 31: G997-
Ozer J, Ratner M, Shaw M, Bailey W, Schomaker S: The current state of serum biomarkers of hepatotoxicity. Toxicology. 2008, 245: 194-205. 10.1016/j.tox.2007.11.021.
Hara-Chikuma M, Sohara E, Rai T, Ikawa M, Okabe M, Sasaki S, Uchida S, Verkman A: Progressive adipocyte hypertrophy in aquaporin-7-deficient mice adipocyte glycerol permeability as a novel regulator of fat accumulation. J Biol Chem. 2005, 280: 15493-15496. 10.1074/jbc.C500028200.
Wan DF, Gong Y, Qin WX, Zhang PP, Li JJ, Wei L, Zhou XM, Li HN, Qiu XK, Zhong F: Large-scale cDNA transfection screening for genes related to cancer development and progression. Proc Natl Acad Sci U S A. 2004, 101: 15724-15729. 10.1073/pnas.0404089101.
Yang CJ, Liu HW, Wang LC: Study on the differential gene expressions of chronic hepatitis B patients of gan depression pi deficiency syndrome and pi-wei damp-heat syndrome. Chin J Integrated Tradit West Med. 2012, 32: 1032-1037.
Liu SN, Tao XP, Wang RB: Study on Chinese syndromes developing rules of chronic hepatitis B. Chin J Trad Med Sci Technol. 2008, 3: 002-
Li XH, Li JJ, Liu YY, Chen JX: Central neurobiological mechanism of liver depression and spleen deficiency syndrome based on chronic stress: a review. J Chinese Int Medicine. 2012, 10: 1-10.3736/jcim20120101.
Yue LF, Ding J, Chen JX, Yue GX, Liang Y, Huo SK, Li JJ: Establishment and review of rat model of syndrome of liver depression with spleen insufficiency. J Beijing University of Trad Chinese Med. 2008, 6: 396-400.
The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1472-6882/14/371/prepub
This study was supported by the National Science and Technology Major Project of China (no. 2012ZX10005001-004 and no. 2009ZX09311-003). Key Program of National Natural Science Fundation of China (81330084), and E-Institutes of Shanghai Municipal Education Commission (no. E03008). The “Yang Fan” Program of Shanghai Committee of science and Technology Fund annotation (no. 14YF1411400).
The authors declare that they have no competing interests.
YYL, QLC and YG collected the samples, carried out the molecular genetic studies, and drafted the manuscript. YYL and ZZG carried out the RT-PCR. YYL, HZ and WZ participated in the design of the study and performed the statistical analysis. YYH and SBS conceived of the study, and participated in its design and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.