Biological Differentiation of Dampness-Heat Syndromes in Chronic Hepatitis B: From Comparative MicroRNA Microarray Profiling to Biomarker Identification

Wen, Li; Jiang, Cen; Wan, Ting-Jun; Wang, Dong; Yan, Di; Li, Gui-Yu; Su, Yue; Liu, Xi-Yang; Rong, Li-Jun; Ye, Hua; Li, Bai-Xue; Feng, Quan-Sheng

doi:https://doi.org/10.1155/2020/7234893

Evidence-Based Complementary and Alternative Medicine

On this page

Abstract Introduction Materials and Methods Results Discussion Conclusion Data Availability Conflicts of Interest Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2020 | Article ID 7234893 | https://doi.org/10.1155/2020/7234893

Biological Differentiation of Dampness-Heat Syndromes in Chronic Hepatitis B: From Comparative MicroRNA Microarray Profiling to Biomarker Identification

Li Wen,¹Cen Jiang,¹Ting-Jun Wan,¹Dong Wang,¹Di Yan,¹Gui-Yu Li,¹Yue Su,¹Xi-Yang Liu,¹Li-Jun Rong,²Hua Ye,³Bai-Xue Li,¹and Quan-Sheng Feng¹

Academic Editor: Juntra Karbwang

Received27 Sept 2019

Revised11 Dec 2019

Accepted18 Dec 2019

Published20 Jan 2020

Abstract

Increasing interest is aroused by traditional Chinese medicine (TCM) treatment of chronic hepatitis B (CHB) based on specific TCM syndrome. As the most common CHB syndromes, spleen-stomach dampness-heat (SSDH) syndrome and liver-gallbladder dampness-heat (LGDH) syndrome are still apt to be confused in TCM diagnosis, greatly hindering the stable exertion of TCM effectiveness. It is urgently needed to provide objective and biological evidences for differentiation and identification of the two significant syndromes. In this study, microRNA (miRNA) microarray analyses coupled with bioinformatics were employed for comparative miRNA profiling of SSDH and LGDH patients. It was found that the two syndromes had both the same and different significantly differentially expressed miRNAs (SDE-miRNAs). Commonness and specificity were also both found between their SDE-miRNA-based bioinformatics analyses, including Hierarchical Clustering, Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, and miRNA-GO/pathway networks. Furthermore, syndrome-specific SDE-miRNAs were identified as the potential biomarkers, including hsa-miR-1273g-3p and hsa-miR-4419b for SSDH as well as hsa-miR-129-1-3p and hsa-miR-129-2-3p for LGDH. All these laid biological and clinical bases for classification and diagnosis of the two significant CHB dampness-heat syndromes including SSDH and LGDH, providing more opportunities for better application of TCM efficacy and superiority in CHB treatment.

1. Introduction

Chronic hepatitis B (CHB) is a potentially life-threatening liver disease caused by hepatitis B virus (HBV) infection. It can progress to cirrhosis and hepatocellular carcinoma (HCC), leading to continuously increasing morbidity [1, 2]. It is estimated that 257 million people are living with HBV and Ca. 887 thousand people die annually from HBV-related disease worldwide [3]. These result in increasing healthcare cost and socioeconomic burdens [4]. CHB always represents a major global health challenge, and its prevention and control remains a significant issue.

Traditional Chinese medicine (TCM), with thousands of years of effective clinical practice, has become an important complementary and alternative medical system and aroused increased attention [5, 6]. It has been verified that TCM can relieve the clinical symptoms, reduce the liver injury, and slow the disease progression of CHB patients. [7–9]. TCM syndrome (also called “ZHENG”) is the basic concept of TCM theory. It describes special phenotype with comprehensive symptoms and signs of patients at a particular stage of disease [10]. The identification of TCM syndrome is the key to guide the specific TCM prescription [11, 12].

According to the standards of TCM syndrome differentiation for viral hepatitis [13], dampness-heat syndrome is the top popular CHB syndrome. In our previous study on CHB syndrome distribution (from 1260 CHB patients) [14], it was also found that dampness-heat syndrome occupied the largest proportion (46%) of the total CHB syndromes. All these indicated the significant importance of dampness-heat syndrome in CHB. Notably, spleen-stomach dampness-heat (SSDH) syndrome and liver-gallbladder dampness-heat (LGDH) syndrome are not only the most common dampness-heat syndromes but also the top two CHB syndromes, taking up 17% and 15%, respectively. However, it is still difficult to discriminate between SSDH and LGDH syndromes and apt to confuse their TCM treatments because they have some similar symptoms (also, see Supplementary Materials). Thus, the two dampness-heat syndromes were focused on in this study.

Currently, the classification of TCM syndromes is still debated [15] since the biological evidence remains lacking. In addition, syndrome differentiation currently depends on traditional observation, auscultation, interrogation, palpation, and the clinical experiences of TCM practitioners [16, 17]. Such diagnosis is usually accompanied with subjectivity, ambiguity, and nonrepeatability to some extent, which greatly impedes the stable exertion of TCM effectiveness and superiority [18]. This issue is particularly prominent for the identification of SSDH and LGDH syndromes. In traditional TCM diagnosis, SSDH and LGDH were classified by both the common and different symptoms (see the Supplementary Materials). Patients with SSDH syndrome usually present some special symptoms such as distending pain in lateral thorax, whereas LGDH patients are usually accompanied with ventosity, distention, and fullness. Meanwhile, SSDH and LGDH patients share the same symptoms such as yellow and slimy fur, nausea, and yellow urine, which may lead to doubt in their differences and confusion in discrimination. Thus, it is urgently needed to scientifically prove the syndrome classification and explore objective diagnosis for SSDH and LGDH syndromes.

MicroRNAs (miRNAs) are endogenous, noncoding, single-stranded RNAs of 19–25 nucleotides in length [19]. They play vital roles in regulating the global signaling networks and pathways involved in the pathological progression of liver diseases [20]. MiRNAs are attractive as potential biomarkers in recent years because they are specific to various disease states and reasonably stable under various environments [21]. Due to the advantages of high throughput, sensitivity, and accuracy [22], miRNA microarray profiling is regarded as a powerful strategy for demonstrating the expression levels of overall miRNAs [23]. It has become a predictive biomarker signature for detecting and distinguishing human diseases from healthy control (HC) over the past decade, such as colorectal cancer [24], chronic hepatitis C [25], and hepatocellular carcinoma [26]. Su et al. utilized circulating miRNAs to reveal the differences between two TCM syndromes in CHB [15, 27, 28], which inspired us to employ comparative miRNAs profiling to differentiate between the two CHB dampness-heat syndromes in this work.

We aim at biologically verifying the differentiation of SSDH and LGDH syndromes and identifying syndrome-specific miRNA biomarkers. Using miRNA microarray profiling, significantly differentially expressed miRNAs (SDE-miRNAs) of SSDH and LGDH patients were screened out firstly. These SDE-miRNAs were comparably analyzed by bioinformatics assays including Gene Ontology (GO) [29], Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways [30], and miRNA-GO/pathway networks [31]. Thus the biological similarities and differences between SSDH and LGDH can be dissected and proved. Quantitative real-time polymerase chain reaction (qRT-PCR) was further conducted to measure the level of syndrome-specific SDE-miRNAs, so that potential diagnostic biomarkers for each syndrome can be validated.

2. Materials and Methods

2.1. Clinical Evaluation

Participants were satisfied with inclusion criteria including CHB diagnostic criteria of western medicine derived from “Guidelines for Prevention and Treatment of Chronic Hepatitis B (December 2015)” issued by Hepatology Branch and Infectious Disease Branch of China Medical Association. The diagnostic criteria for SSDH syndromes of CHB referred to the Consensus of Experts in TCM Diagnosis and Treatment (2012) as well as the research results and related monographs of the National Major Project of Science and Technology (no. 2012ZX10005001). The LGDH patients were diagnosed according to the Guidelines for the Diagnosis and Treatment of CHB in Traditional Chinese Medicine (2018) issued by Hepatobiliary Diseases Professional Committee of China Association of Traditional Chinese Medicine. The detailed contents of the inclusion criteria and exclusion criteria can also be seen in the Supplementary Materials.

2.2. Patients and Blood Collection

CHB patients with SSDH or LGDH syndromes were aged from 22 to 51 years and came from Chengdu Public Health and Medical Center. For both the miRNA microarray profiling and qRT-PCR analyses, peripheral blood of 15 SSDH patients, 15 LGDH patients, and 15 HC participants were collected in an EDTA anticoagulant tube. The plasma was collected from peripheral blood (4 mL) following the manufacturer’s instructions. In brief, peripheral blood was centrifuged for 10 min at 1,700 g and 4°C. The supernatant was collected and centrifuged for 10 min at 2,000 g and 4°C. The supernatant was transferred into polypropylene tube as blood plasma sample and stored at −80°C.

2.3. RNA Extraction and miRNA Microarray Analysis

RNA isolation was performed on each blood plasma sample obtained as above. Total RNA was isolated by using Trizol, which was then followed by Qiagen miRNeasy Mini kit (Qiagen, Germany) according to the manufacturer’s instructions. The isolated RNA was quantified with the spectrophotometer. The quality of isolated RNA was inspected by formaldehyde gel electrophoresis.

Genome-wide microRNA microarray profiling was performed using a human miRNA microarray platform (Agilent, USA). In short, 200 nanogram of miRNA was labelled using miRNA Complete Labeling and Hyb kit. Dried samples were placed into the hybridization oven overnight. After hybridization and washing, signals were measured by Agilent microarray scanner (G2565CA). The picture analysis and data extraction were processed with Agilent Feature Extraction software. The data were normalized by using Agilent Gene Spring software.

2.4. Bioinformatics Analysis

TargetScanHuman and miRTarBase were used for target gene prediction of the SDE-miRNAs. Hierarchical clustering analysis was performed with Cluster 3.0 software, which was originally written by Michael Eisen at Stanford University. GO and signaling pathway analysis were performed based on the DAVID Bioinformatics Resources 6.7 and KEGG database. The miRNA-target gene interactions, miRNA-GO networks and miRNA-pathway networks were analyzed by using Cytoscape software.

2.5. qRT-PCR for miRNA Verification

Significantly dysregulated miRNAs were validated by qRT-PCR. The miRNAs were separated by using miRcute serum/plasma miRNA extraction and separation kit (TIANGEN, China). RNA was reverse-transcribed to cDNA by employing cDNA synthesis kit (Exiqon, Denmark). qRT-PCR system (Eppendorf, Germany) combining with the predesigned primers (ABM, Canada) was used for miRNA quantification. For the reaction conditions, polymerase activation/denaturation was performed for 10 min at 95°C. 40 amplification cycles at 95°C for 10 seconds, 63°C for 15 seconds, and 72°C for 32 seconds were performed for miRNA quantification, followed by signal detection. The relative amount of miRNA was normalized against U6 snRNA (the internal control), and the fold change (FC) in the amount of each miRNA compared with the HC group was calculated by using the 2^−ΔΔCT method.

2.6. Statistical Analysis

The calculation of mean ± standard deviation (SD) and Students’ t test was performed with GraphPad Prism 6.0 software. Only those miRNAs with FC >2 (or <0.5) and compared with the HC group were considered as the SDE-miRNAs, and only these with all FC >1.5 (or <0.7) and were considered as the potential biomarkers in the study.

3. Results

3.1. Patient Characteristics and miRNA Microarray Profiling

In order to gain biological insights into the similarities and differences between SSDH and LGDH syndromes, and further discover their miRNA biomarkers, 45 participants were included in our study, including 15 SSDH patients, 15 LGDH patients, and 15 HC volunteers. As shown in Supplementary Table S1, the sex and age were not significantly different between the three groups (). The distribution of ALT, AST, and TBIL were also found no significant differences between SSDH and LGDH groups (). The virus load in SSDH was higher than that in LGDH, but there was no significant difference (). The above clinical baseline data indicated that the participants were available for the next comparative miRNA microarray profiling.

According to the quantitative data of the identified miRNAs in microarray analysis, those with the FC > 2 (or<0.5) and compared with the HC group were screened as the SDE-miRNAs, which were considered to be more valuable for further bioinformatics analysis [20]. Of these, 7 SDE-miRNAs were upregulated, and 3 SDE-miRNAs were downregulated in the SSDH group (Supplementary Table S2), while 12 SDE-miRNAs were upregulated in the LGDH group (Supplementary Table S3). Moreover, 4 of these SDE-miRNAs, namely, hsa-miR-122-5p, hsa-miR-320e, hsa-miR-1260a, and hsa-miR-483-3p, were shared by both SSDH and LGDH groups.

3.2. Hierarchical Clustering Analysis and Target Gene Prediction

To deeply dissect the biological commonalities and differences between SSDH and LGDH syndromes, intensive bioinformatics analysis of the SDE-miRNAs were performed. Firstly, the hierarchical clustering analysis [32] of SDE-miRNAs is shown in Figure 1. Then, the target genes of the SDE-miRNAs were predicted based on agreement between the databases of TargetScanHuman and miRTarBase [33]. A total of 1000 and 1200 target genes were predicted for the 10 SDE-miRNAs in SSDH and the 12 SDE-miRNAs in LGDH, respectively. To further show the regulatory relationships between SDE-miRNAs and their target genes, the biological interaction networks [34] were built (Supplementary Figure S1).

(a)

(b)

3.3. Comparative GO Annotation for SSDH and LGDH Syndromes

GO enrichment analysis were performed to determine the biological functions of the target genes of SDE-miRNAs. GO analysis mainly consists of three components: biological processes, cellular components, and molecular functions [29]. As shown in Figure 2, the top 10 terms of each GO component were plotted and compared between SSDH and LGDH syndromes. Notably, although the two dampness-heat syndromes shared some common terms in each component, there were some syndrome-specific terms. For biological processes (red histograms), positive regulation of cellular biosynthetic process (GO:0031328) and positive regulation of biosynthetic process (GO:0009891) were specific for SSDH, and whole phosphorus metabolic process (GO:0006793) and phosphate metabolic process (GO:0006796) were specifically related to LGDH syndrome. For cellular components (green histograms), the targets genes of the SDE-miRNAs in SSDH and LGDH syndromes were specifically involved with insoluble fraction (GO:0005626) and Golgi apparatus (GO:0005794), respectively. As for molecular functions (blue histograms), 4 terms were specific for each TCM syndromes, including ion binding (GO:0043167), metal ion binding (GO:0046872), transcription activator activity (GO:0016563) and transcription repressor activity (GO:0016564) for SSDH syndrome, as well as transcription factor binding (GO:0008134), protein kinase activity (GO:0004672), protein dimerization activity (GO:0046983), and identical protein binding (GO:0042802) for LGDH syndrome.

(a)

(b)

3.4. Comparative KEGG Pathway Analysis for SSDH and LGDH Syndromes

KEGG pathway analysis [30] was carried out to further understand the functions and signaling pathways of the target genes. Twenty most commonly observed pathways were showed to compare the two syndromes (Figure 3). The result showed that the target genes of SDE-miRNAs in SSDH and LGDH group both mainly functioned in the MAPK signaling pathway (hsa04010), neurotrophin signaling pathway (hsa04722), chemokine signaling pathway (hsa04062), and endocytosis (hsa04144). Notably, the target genes in SSDH group specifically linked to 8 signaling pathways including ubiquitin-mediated proteolysis (hsa04120), p53 signaling pathway (hsa04115), regulation of actin cytoskeleton (hsa04810), cell adhesion molecules (CAMs) (hsa04514), purine metabolism (hsa00230), Wnt signaling pathway (hsa04310), glioma (hsa05214) and epithelial cell signaling in Helicobacter pylori infection (hsa05120). For LGDH syndrome, the target genes were also particularly involved in 8 signaling pathways including cytokine-cytokine receptor interaction (hsa04060), focal adhesion (hsa04510), Jak-STAT signaling pathway (hsa04630), TGF-beta signaling pathway (hsa04350), renal cell carcinoma (hsa05211), pancreatic cancer (hsa05212), colorectal cancer (hsa05210) and ErbB signaling pathway (hsa04012).

(a)

(b)

3.5. Comparative miRNA-GO/Pathway Network Analysis for SSDH and LGDH Syndromes

To further understand the association relationships between SDE-miRNAs and their corresponding GO terms and signaling pathways, miRNA-GO networks (Figure 4) and miRNA-pathway networks (Figure 5) were analyzed and demonstrated. As was shown in Figures 4 and 5, the SDE-miRNAs (the red squares) of SSDH or LGDH were all involved in the same network, suggesting that there were close biological correlations among the SDE-miRNAs of the syndrome. Furthermore, according to the study on biology network [35], the size of the red square indicates the complexity of regulatory relationship and the importance of the miRNA. For miRNA-GO networks, hsa-miR-122-5p, hsa-miR-1260a, hsa-miR-3196, and hsa-miR-15b-5p played crucial roles in modulating the molecular networks in SSDH syndrome (Figure 4(a)), whereas hsa-miR-483-3p, hsa-miR-22-3p, hsa-miR-21-5p, and hsa-miR-129-1-3p might be key regulators of pathogenesis in LGDH syndrome (Figure 4(b)). Moreover, the visualized miRNA-pathway networks indicated that hsa-miR-483-3p, hsa-miR-122-5p, hsa-miR-3196, and hsa-miR-15b-5p played prominent roles in signaling pathways related to SSDH syndrome (Figure 5(a)). Comparatively, hsa-miR-483-3p, hsa-miR-21-5p, hsa-miR-129-2-3p, and hsa-miR-22-3p carried considerable weight in LGDH groups (Figure 5(b)).

(a)

(b)

(a)

(b)

3.6. Validation of the Biomarkers for SSDH and LGDH Syndromes

To validate the microarray results and further identify potential biomarkers for confirming and distinguishing between SSDH and LGDH syndromes in CHB, syndrome-specific SDE-miRNAs were randomly selected for qRT-PCR analysis. As shown in Figure 6, the expression level of 6 SDE-miRNAs was measured. To gain more reliability for syndrome-specificity, the expression level of each SDE-miRNA in one syndrome must suffer the comparison with that in the other syndrome and that in the HC group simultaneously. Only the SDE-miRNAs with both FC values >1.5 (or <0.7) and both values <0.05 can be considered as biomarkers. For example, the expression level of hsa-miR-1273g-3p in SSDH group was significantly higher than both that in the HC group (FC = 2.12, ) and that in the LGDH group (FC = 2.45, ) (Figure 6(a)). Likewise, hsa-miR-4419b was significantly highly expressed in the SSDH group when compared with both HC and LGDH groups, with FC values of 1.88, 1.85 and value was both lower than 0.05 (Figure 6(b)). For LGDH group, a significant hsa-miR-129-1-3p increase was shown when compared with the HC group (FC = 2.71, ) and SSDH group (FC = 3.12, ) (Figure 6(d)). Hsa-miR-129-2-3p was significantly overexpressed in LGDH group when compared with HC and SSDH groups, with FC values of 2.79, 2,49 and both values lower than 0.001 (Figure 6(e)). These results are also in accordance with the microarray data. Thus, it was believed that hsa-miR-1273g-3p and hsa-miR-4419b could be employed as SSDH biomarkers and that hsa-miR-129-1-3p and hsa-miR-129-2-3p were available as LGDH biomarkers. However, hsa-miR-3196 (Figure 6(c)) and hsa-miR-21-5p (Figure 6(f)) with were excluded as biomarkers.

(a)

(b)

(c)

(d)

(e)

(f)

4. Discussion

TCM has been proved effective for CHB treatment in practice for a long time. However, biological basis for classification and diagnosis of TCM syndromes of CHB are still lacking. Dampness-heat syndrome occupies the largest proportion of the total CHB syndromes. The top two dampness-heat syndromes, SSDH and LGDH are apt to be confused since they have some similar symptoms. Thus, the two syndromes were investigated in this study. Based on the microarray profiling data, the SDE-miRNAs including hsa-miR-122-5p, hsa-miR-320e, hsa-miR-1260a, and hsa-miR-483-3p might be related with the progression of CHB dampness-heat syndrome, as they were identified to be significant in both SSDH and LGDH groups. Hsa-miR-1273g-3p, hsa-miR-4419b, hsa-miR-451a hsa-miR-3196, hsa-miR-223-3p, and hsa-miR-15b-5p might be specifically linked to SSDH syndrome. Meanwhile, hsa-miR-129-1-3p, hsa-miR-129-2-3p, hsa-miR-21-5p, hsa-miR-1304-3p hsa-miR-30d-5p, hsa-miR-762, hsa-miR-4532, and hsa-miR-22-3p might play an important and special role in LGDH syndrome. These results indicated that both syndrome-common and syndrome-specific mechanisms may exist between SSDH and LGDH.

Most of the SDE-miRNAs have already been reported in previous virus hepatopathy-related studies. In general, miRNA could act as a cellular antiviral defense or help viruses establish a favorable environment for their replication and survival [20]. Hsa-miR-122 is a liver-specific miRNA and associated with the immune control of chronic HBV infection [36]. It has been found that hsa-miR-122-5p was upregulated during different phases of chronic HBV infection [37], which is in accordance with our microarray results. Hsa-miR-483-3p was differentially expressed in peripheral blood mononuclear cell from the chronic asymptomatic carriers [38] and upregulated in HBV-associated HCC [39]. Hsa-miR-320 was identified as the miRNA whose expression levels were altered by hepatitis virus infection [40] and showed downregulation in hepatitis B patients [41]. It has been demonstrated that hsa-miR-451a was differentially expressed in CHB and might play a crucial role in global signaling networks and pathways involved in CHB pathogenesis [20]. Hsa-miR-15b has been reported to be important during HBV infection because it could promote HBV replication by augmenting HBV enhancer I activity via directly targeting hepatocyte nuclear factor 1α [42]. Hsa-miR-1273g-3p might affect the activation and apoptosis of HSCs by targeting PTEN in hepatitis virus [43]. Hsa-miR-223-3p was found to be significantly dysregulated in HBV-positive patients [44]. It was further confirmed as a novel noninvasive biomarker of HBV-positive HCC at a very early stage of liver disease [45]. Hsa-miR-22 played a prominent role in HBV-related diseases [46], as it might be involved in HBV infection [47]. Hsa-miR-22-3p was found differentially expressed during different pathologic processes of CHB [20]. It has been indicated that hsa-miR-21-5p participated in inflammatory responses and hepatocyte proliferation and was closely related with liver disease [48–50]. Serum hsa-miR-21-5p was significantly elevated in CHB patients [51], which is consistent with our results. Hsa-miR-30d was significantly overexpressed in HBV-associated HCC patients [52], which was in accordance with the upregulation of hsa-miR-30d in our study. Khairy et al. have demonstrated that there was significant fold change of hsa-miR-129 in HCC patients compared with the HC group [53]. Hsa-miR-129-2 has been evaluated as a potential early diagnostic biomarker for HBV-related HCC since 85% of HCC patients at stage I could be distinguished by their miR-129-2 methylation levels [54]. Hsa-miR-4532 has been identified as a vital miRNA demonstrating a strong expressional response to HBV and could be used as an early diagnostic biomarker of hepatitis B [41]. All these above suggested that most of the identified SDE-miRNAs were biologically related to virus hepatopathy, and thus the microarray results were further validated. Some of the SDE-miRNAs (such as hsa-miR-1260a, hsa-miR-3196, and hsa-miR-762) have not been reported in virus hepatopathy-related studies, indicating that the finding in our research is novel and needs further investigation.

Detailed dissection of the identified SDE-miRNAs and the regulation mechanisms is of great significance for further understanding the pathogenesis and biological differences of SSDH and LGDH syndromes in CHB. Thus, intensive bioinformatics analysis based on the predicted target genes of SDE-miRNAs was carried out. According to GO data (Figure 2), the two dampness-heat syndromes shared some common GO terms, and they also owned syndrome-specific terms. KEGG pathway analysis (Figure 3) demonstrated that the target genes of SDE-miRNAs in the two syndromes both mainly functioned in the MAPK signaling pathway, neurotrophin signaling pathway, chemokine signaling pathway, and endocytosis. Meanwhile, there were also some syndrome-specific signaling pathways, such as ubiquitin-mediated proteolysis for SSDH syndrome and cytokine-cytokine receptor interaction for LGDH syndrome. Finally, the association relationship and the importance of the SDE-miRNAs in each TCM syndrome were determined by setting up miRNA-GO and miRNA-pathway networks (Figures 4 and 5). For SSDH and LGDH groups, although there were the same SDE-miRNAs (such as hsa-miR-483-3p and hsa-miR-122-5p) that played crucial roles in the networks, most of the key SDE-miRNAs in the networks of the two syndromes were different. All these above further validated that both biological similarity and differences existed between SSDH and LGDH syndromes. The results provided more evidence for the TCM practices, in which the syndromes are based on the same disease and discriminated by different symptoms of patients [11], and in which the prescription of TCM medications is determined by both the disease and the specific symptoms.

Furthermore, in order to provide more biological bases for objective diagnosis and accurate treatment of the two dampness-heat syndromes, qRT-PCR were experimented to identify potential biomarkers based on the syndrome-specific SDE-miRNAs. Comparisons among the three groups (SSDH, LGDH, and HC) were carried out. Only those with all FC >1.5 (or <0.7) and were considered, and thus more reliability for syndrome-specificity was provided. Taken together, 4 SDE-miRNAs were validated as the diagnostic biomarkers, including hsa-miR-1273g-3p and hsa-miR-4419b for SSDH as well as hsa-miR-129-1-3p and hsa-miR-129-2-3p for LGDH.

However, some limitations also exist in this work. In order to enhance the value of this study in future, a large number of CHB patients with each syndrome should be further enrolled and tested. The importance of the potential diagnostic biomarkers of SSDH and LGDH syndromes should be further confirmed. Studies on molecular regulating mechanisms of the identified SDE-miRNAs in each syndrome need to be carried out, so that the pathogenesis of dampness-heat syndrome in CHB and therapeutic targets for treatment can be deeply elucidated. Moreover, other dampness-heat syndromes should also be investigated to comprehensively understand the dampness-heat syndrome of CHB.

5. Conclusion

By utilizing miRNA array profiling, 10 and 12 SDE-miRNAs were identified in SSDH and LGDH syndromes of CHB, respectively. Among these SDE-miRNAs, 4 were found in both dampness-heat syndromes, but 6 and 8 specifically linked to SSDH and LGDH, respectively. These results indicated that both biological similarity and syndrome-specificity existed in SSDH and LGDH patients. This was further confirmed by bioinformatics analyses, in which both same and different GO terms, KEGG pathways, and miRNA-GO/pathway networks were found between the two syndromes. Furthermore, syndrome-specific SDE-miRNAs were identified as the potential biomarkers, including hsa-miR-1273g-3p and hsa-miR-4419b for SSDH and hsa-miR-129-1-3p and hsa-miR-129-2-3p for LGDH, respectively. All these laid scientific basis for the differentiation and diagnosis of the two significant dampness-heat syndromes in CHB, providing more opportunities for stable exertion and better application of the efficacy and superiority of TCM in CHB treatment.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

The authors thank all the blood donors from the Chengdu Military Region General Hospital (Sichuan, China) and Dr. Shan-Hong Tang for helping us collect the blood samples. This research was supported by the Traditional Chinese Medicine Bureau of Sichuan Province (no. 2018JC003), Postdoctoral Research Foundation of China (no. 2017M622984), Education Department of Sichuan Province (no. 18ZB0199), Chengdu University of Traditional Chinese Medicine Foundation (nos. ZRQN1715, ZRYY1737, BSH2018001, and YXRC2018002), National Major Science and Technology Program of China (no. 2017ZX10205501), Foreign High-End Talent Introduction Project of Sichuan Province (no. SYZ201983), Technology Bureau of Chengdu City (no. 2016-XT00-00033-GX), and National Natural Science Foundation of China (no. 81803976).

Supplementary Materials

Table S1: clinical diagnosis of patients who participated in this study. Table S2: the SDE-miRNAs in SSDH patients. Table S3: the SDE-miRNAs in LGDH patients. Figure S1: the miRNA-target gene networks built using the target genes of the SDE-miRNAs in SSDH (a) and LGDH (b). The red squares (the central nodes) represent SDE-miRNAs, and the blue spots represent the target genes. The lines indicate the interactions between SDE-miRNAs and the target genes. (Supplementary Materials)

References

J. Y. Dai, S. J. Sun, J. H. Peng et al., “Exploration of macro-micro biomarkers for dampness-heat syndrome differentiation in different diseases,” Evidence-Based Complementary and Alternative Medicine, vol. 2013, Article ID 706762, 9 pages, 2013.
View at: Publisher Site | Google Scholar
E. A. Bard-Chapeau, A.-T. Nguyen, A. G. Rust et al., “Transposon mutagenesis identifies genes driving hepatocellular carcinoma in a chronic hepatitis B mouse model,” Nature Genetics, vol. 46, no. 1, pp. 24–32, 2014.
View at: Publisher Site | Google Scholar
World Health Organization, Hepatitis B, WHO, Geneva, Switzerland, 2018, http://www.who.int/en/news-room/fact-sheets/detail/hepatitis-b.
S. T. Goldstein, F. Zhou, S. C. Hadler, B. P. Bell, E. E. Mast, and H. S. Margolis, “A mathematical model to estimate global hepatitis B disease burden and vaccination impact,” International Journal of Epidemiology, vol. 34, no. 6, pp. 1329–1339, 2005.
View at: Publisher Site | Google Scholar
M. Jiang, J. Yang, C. Zhang et al., “Clinical studies with traditional Chinese medicine in the past decade and future research and development,” Planta Medica, vol. 76, no. 17, pp. 2048–2064, 2010.
View at: Publisher Site | Google Scholar
W.-T. Liao, J.-H. Chiang, C.-J. Li, M.-T. Lee, C.-C. Su, and H.-R. Yen, “Investigation on the use of traditional Chinese medicine for polycystic ovary syndrome in a nationwide prescription database in taiwan,” Journal of Clinical Medicine, vol. 7, no. 7, p. 179, 2018.
View at: Publisher Site | Google Scholar
L. Zhang, G. Wang, W. Hou, P. Li, A. Dulin, and H. L. Bonkovsky, “Contemporary clinical research of traditional Chinese medicines for chronic hepatitis B in China: an analytical review,” Hepatology, vol. 51, no. 2, pp. 690–698, 2010.
View at: Publisher Site | Google Scholar
H. M. Li, Z. H. Ye, X. Gao et al., “Diwu Yanggan capsule improving liver histological response for patients with HBeAg-negative chronic hepatitis B: a randomized controlled clinical trial,” American Journal of Translational Research, vol. 10, no. 5, pp. 1511–1521, 2018.
View at: Google Scholar
Y. N. Song, J. J. Sun, Y. Y. Lu et al., “Therapeutic efficacy of fuzheng-huayu tablet based traditional Chinese medicine syndrome differentiation on hepatitis-B-caused cirrhosis: a multicenter double-blind randomized controlled trail,” Evidence-Based Complementary and Alternative Medicine, vol. 2013, Article ID 709305, 8 pages, 2013.
View at: Publisher Site | Google Scholar
J. Y. Dai, S. J. Sun, J. M. Cao et al., “Similar connotation in chronic hepatitis B and nonalcoholic fatty liver patients with dampness-heat syndrome,” Evidence-Based Complementary and Alternative Medicine, vol. 2013, Article ID 793820, 12 pages, 2013.
View at: Publisher Site | Google Scholar
H. Kang, Y. Zhao, C. Li et al., “Integrating clinical indexes into four-diagnostic information contributes to the Traditional Chinese Medicine (TCM) syndrome diagnosis of chronic hepatitis B,” Scientific Reports, vol. 5, p. 9395, 2015.
View at: Publisher Site | Google Scholar
H.-P. Xie, H.-Z. Yang, W.-K. Wu et al., “Chinese medicine syndrome distribution of chronic hepatitis B virus carriers in immunotolerant phase,” Chinese Journal of Integrative Medicine, vol. 20, no. 2, pp. 94–100, 2014.
View at: Publisher Site | Google Scholar
Branch of Hepatobiliary Diseases and Chinese Association of Chinese Medicine, “The standards of traditional Chinese medicine syndrome differentiation for viral hepatitis,” Journal of Clinical Hepatology, vol. 33, pp. 1839–1846, 2017, in Chinese.
View at: Google Scholar
Z. Wang, Research on Correlations Among TCM Common Type of Syndromes and Laboratory Index in 1260 Patients with CHB, Chengdu University of traditional Chinese Medicine, Chengdu, China, 2011, in Chinese.
H. Zhang, Y. Guan, Y. Y. Lu, Y. Y. Hu, S. Huang, and S. B. Su, “Circulating miR-583 and miR-663 Refer to ZHENG Differentiation in Chronic Hepatitis B,” Evidence-Based Complementary and Alternative Medicine, vol. 2013, Article ID 751341, 8 pages, 2013.
View at: Publisher Site | Google Scholar
F. Cheung, “TCM: made in China,” Nature, vol. 480, no. 7378, pp. S82–S83, 2011.
View at: Publisher Site | Google Scholar
Y. Zhao, H. Kang, J.-H. Peng, L. Xu, Z.-W. Cao, and Y.-Y. Hu, “Key symptoms selection for two major syndromes diagnosis of Chinese medicine in chronic hepatitis B,” Chinese Journal of Integrative Medicine, vol. 23, no. 4, pp. 253–260, 2017.
View at: Publisher Site | Google Scholar
L. Wen, Y. F. Liu, C. Jiang et al., “Comparative proteomic profiling and biomarker identification of traditional Chinese medicine-based HIV/AIDS syndromes,” Scientific Reports, vol. 8, no. 1, p. 4187, 2018.
View at: Publisher Site | Google Scholar
B. Yang, H. Lin, J. Xiao et al., “The muscle-specific microRNA miR-1 regulates cardiac arrhythmogenic potential by targeting GJA1 and KCNJ2,” Nature Medicine, vol. 13, no. 4, pp. 486–491, 2007.
View at: Publisher Site | Google Scholar
X. L. Hou, Y. Liang, J. N. Chen et al., “Expression profiling of cellular MicroRNA in asymptomatic HBsAg carriers and chronic hepatitis B patients,” BioMed Research International, vol. 2017, Article ID 6484835, 17 pages, 2017.
View at: Publisher Site | Google Scholar
I. Alevizos and G. G. Illei, “MicroRNAs as biomarkers in rheumatic diseases,” Nature Reviews Rheumatology, vol. 6, no. 7, pp. 391–398, 2010.
View at: Publisher Site | Google Scholar
P. T. Nelson, D. A. Baldwin, L. M. Scearce, J. C. Oberholtzer, J. W. Tobias, and Z. Mourelatos, “Microarray-based, high-throughput gene expression profiling of microRNAs,” Nature Methods, vol. 1, no. 2, pp. 155–161, 2004.
View at: Publisher Site | Google Scholar
J. Lange, “microRNA profiling on automated biochip platform reveals biomarker signatures from blood samples,” Nature Methods, vol. 7, p. 17, 2010.
View at: Publisher Site | Google Scholar
Y. Ma, P. Zhang, J. Yang, Z. Liu, Z. Yang, and H. Qin, “Candidate microRNA biomarkers in human colorectal cancer: systematic review profiling studies and experimental validation,” International Journal of Cancer, vol. 130, no. 9, pp. 2077–2087, 2012.
View at: Publisher Site | Google Scholar
B. Liu, Y. Xiang, and H.-S. Zhang, “Circulating microRNA-196a as a candidate diagnostic biomarker for chronic hepatitis C,” Molecular Medicine Reports, vol. 12, no. 1, pp. 105–110, 2015.
View at: Publisher Site | Google Scholar
L. Wei, B. Lian, Y. Zhang et al., “Application of microRNA and mRNA expression profiling on prognostic biomarker discovery for hepatocellular carcinoma,” BMC Genomics, vol. 15, no. Suppl 1, p. S13, 2014.
View at: Publisher Site | Google Scholar
Y. Guan, H. Zhang, W. Zhang, and S. B. Su, “Analysis of differential gene expression profile in peripheral blood of patients with chronic hepatitis B and syndromes of dual deficiency of liver and kidney yin and accumulation of dampness heat,” Journal of Chinese Integrative Medicine, vol. 10, no. 7, pp. 751–756, 2012.
View at: Publisher Site | Google Scholar
Z. Z. Guo, S. H. Yu, Y. Guan et al., “Molecular mechanisms of same TCM syndrome for different diseases and different TCM syndrome for same disease in chronic hepatitis B and liver cirrhosis,” Evidence-Based Complementray and Alternative Medicine, vol. 2012, Article ID 120350, 9 pages, 2012.
View at: Publisher Site | Google Scholar
J. Dutkowski, M. Kramer, M. A. Surma et al., “A gene ontology inferred from molecular networks,” Nature Biotechnology, vol. 31, no. 1, pp. 38–45, 2013.
View at: Publisher Site | Google Scholar
M. J. Moore, T. K. Scheel, J. M. Luna et al., “miRNA-target chimeras reveal miRNA 3’-end pairing as a major determinant of Argonaute target specificity,” Nature Communication, vol. 6, no. 1, p. 8864, 2015.
View at: Publisher Site | Google Scholar
Z. L. Tang, Y. L. Yang, Z. S. Wang, S. P. Zhao, Y. L. Mu, and K. Li, “Integrated analysis of miRNA and mRNA paired expression profiling of prenatal skeletal muscle development in three genotype pigs,” Scientific Reports, vol. 5, no. 1, Article ID 15544, 2015.
View at: Publisher Site | Google Scholar
K. Murakami, U. Günesdogan, J. J. Zylicz et al., “NANOG alone induces germ cells in primed epiblast in vitro by activation of enhancers,” Nature, vol. 529, no. 7586, pp. 403–407, 2016.
View at: Publisher Site | Google Scholar
T. Ishimoto, H. Sugihara, M. Watanabe et al., “Macrophage-derived reactive oxygen species suppress miR-328 targeting CD44 in cancer cells and promote redox adaptation,” Carcinogenesis, vol. 35, no. 5, pp. 1003–1011, 2014.
View at: Publisher Site | Google Scholar
M. E. Smoot, K. Ono, J. Ruscheinski, P.-L. Wang, and T. Ideker, “Cytoscape 2.8: new features for data integration and network visualization,” Bioinformatics, vol. 27, no. 3, pp. 431-432, 2011.
View at: Publisher Site | Google Scholar
Y. Sun, S. Shen, H. Tang et al., “miR-429 identified by dynamic transcriptome analysis is a new candidate biomarker for colorectal cancer prognosis,” OMICS: A Journal of Integrative Biology, vol. 18, no. 1, pp. 54–64, 2014.
View at: Publisher Site | Google Scholar
C. L. Jopling, “Modulation of hepatitis C virus RNA abundance by a liver-specific MicroRNA,” Science, vol. 309, no. 5740, pp. 1577–1581, 2005.
View at: Publisher Site | Google Scholar
M. R. Brunetto, D. Cavallone, F. Oliveri et al., “A serum microRNA signature is associated with the immune control of chronic hepatitis B virus infection,” PLoS One, vol. 9, no. 10, Article ID e110782, 2014.
View at: Publisher Site | Google Scholar
L. Chen, C. Li, Z. Peng, J. Zhao, G. Gong, and D. Tan, “miR-197 expression in peripheral blood mononuclear cells from hepatitis B virus-infected patients,” Gut and Liver, vol. 7, no. 3, pp. 335–342, 2013.
View at: Publisher Site | Google Scholar
W. Wang, L. J. Zhao, Y. X. Tan, H. Ren, and Z. T. Qi, “Identification of deregulated miRNAs and their targets in hepatitis B virus-associated hepatocellular carcinoma,” World Journal of Gastroenterology, vol. 18, no. 38, pp. 5442–5453, 2012.
View at: Publisher Site | Google Scholar
H. Ishida, T. Tatsumi, A. Hosui et al., “Alterations in microRNA expression profile in HCV-infected hepatoma cells: involvement of miR-491 in regulation of HCV replication via the PI3 kinase/Akt pathway,” Biochemical and Biophysical Research Communications, vol. 412, no. 1, pp. 92–97, 2011.
View at: Publisher Site | Google Scholar
S. Jahan, “Expressional studies of microRNAs in hepatitis B patients of quetta, Pakistan,” Pure and Applied Biology, vol. 6, no. 3, pp. 1044–1052, 2017.
View at: Publisher Site | Google Scholar
X. Dai, W. Zhang, H. Zhang et al., “Modulation of HBV replication by microRNA-15b through targeting hepatocyte nuclear factor 1α,” Nucleic Acids Research, vol. 42, no. 10, pp. 6578–6590, 2014.
View at: Publisher Site | Google Scholar
X. Niu, N. Fu, J. Du et al., “miR-1273g-3p modulates activation and apoptosis of hepatic stellate cells by directly targeting PTEN in HCV-related liver fibrosis,” FEBS Letters, vol. 590, no. 16, pp. 2709–2724, 2016.
View at: Publisher Site | Google Scholar
H. Zheng, A. E. Zou, M. A. Saad et al., “Alcohol-dysregulated microRNAs in hepatitis B virus-related hepatocellular carcinoma,” PLoS One, vol. 12, no. 5, Article ID e0178547, 2017.
View at: Publisher Site | Google Scholar
B. G. Giray, G. Emekdas, S. Tezcan et al., “Profiles of serum microRNAs; miR-125b-5p and miR223-3p serve as novel biomarkers for HBV-positive hepatocellular carcinoma,” Molecular Biology Reports, vol. 41, no. 7, pp. 4513–4519, 2014.
View at: Publisher Site | Google Scholar
C. Shi and X. Xu, “MicroRNA-22 is down-regulated in hepatitis B virus-related hepatocellular carcinoma,” Biomedicine & Pharmacotherapy, vol. 67, no. 5, pp. 375–380, 2013.
View at: Publisher Site | Google Scholar
G. Y. Wang, F. L. Dong, Z. Y. Xu et al., “MicroRNA profile in HBV-induced infection and hepatocellular carcinoma,” BMC Cancer, vol. 17, no. 1, p. 805, 2017.
View at: Publisher Site | Google Scholar
Q. Ruan, P. Wang, T. Wang et al., “MicroRNA-21 regulates T-cell apoptosis by directly targeting the tumor suppressor gene Tipe2,” Cell Death & Disease, vol. 5, no. 2, Article ID e1095, 2014.
View at: Publisher Site | Google Scholar
T. X. Lu, J. Hartner, E. J. Lim et al., “MicroRNA-21 limits in vivo immune response-mediated activation of the IL-12/IFN-gamma pathway, Th1 polarization, and the severity of delayed-type hypersensitivity,” Journal of Immunology, vol. 187, no. 6, pp. 3362–3373, 2011.
View at: Publisher Site | Google Scholar
F. J. Sheedy, E. Palsson-McDermott, E. J. Hennessy et al., “Negative regulation of TLR4 via targeting of the proinflammatory tumor suppressor PDCD4 by the microRNA miR-21,” Nature Immunology, vol. 11, no. 2, pp. 141–147, 2010.
View at: Publisher Site | Google Scholar
J. Xu, C. Wu, X. Che et al., “Circulating microRNAs, miR-21, miR-122, and miR-223, in patients with hepatocellular carcinoma or chronic hepatitis,” Molecular Carcinogenesis, vol. 50, no. 2, pp. 136–142, 2011.
View at: Publisher Site | Google Scholar
A. M. Liu, C. Zhang, J. Burchard et al., “Global regulation on microRNA in hepatitis B virus-associated hepatocellular carcinoma,” OMICS: a Journal of Integrative Biology, vol. 15, no. 3, pp. 187–191, 2011.
View at: Publisher Site | Google Scholar
A. Khairy, I. Hamza, O. Shaker, and A. Yosry, “Serum miRNA panel in Egyptian patients with chronic hepatitis C related hepatocellular carcinoma,” Asian Pacific Journal of Cancer Prevention, vol. 17, no. 5, pp. 2699–2703, 2016.
View at: Google Scholar
W. Yin, Y. Zhao, Y. J. Ji et al., “Serum/plasma microRNAs as biomarkers for HBV-related hepatocellular carcinoma in China,” BioMed Research International, vol. 2015, Article ID 965185, 8 pages, 2015.
View at: Publisher Site | Google Scholar

Copyright

Copyright © 2020 Li Wen et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

423

Downloads

936

Citations