Abstract

The functions of genes involved in idiopathic portal hypertension (IPH) remain unidentified. The present study was undertaken to identify the functions of genes expressed in blood samples from patients with IPH through comprehensive analysis of gene expression using DNA microarrays. The data were compared with data from healthy individuals to explore the functions of genes showing increased or decreased expression in patients with IPH. In cluster analysis, no dominant probe group was shown to differ between patients with IPH and healthy controls. In functional annotation analysis using the Database for Annotation Visualization and Integrated Discovery tool, clusters showing dysfunction in patients with IPH involved gene terms related to the immune system. Analysis using network-based pathways revealed decreased expression of adenosine deaminase, ectonucleoside triphosphate diphosphohydrolase 4, ATP-binding cassette, subfamily C, member 1, transforming growth factor-β, and prostaglandin E receptor 2; increased expression of cytochrome P450, family 4, subfamily F, polypeptide 3, and glutathione peroxidase 3; and abnormalities in the immune system, nucleic acid metabolism, arachidonic acid/leukotriene pathways, and biological processes. These results suggested that IPH involved compromised function of immunocompetent cells and that such dysfunction may be associated with abnormalities in nucleic acid metabolism and arachidonic acid/leukotriene-related synthesis/metabolism.

1. Introduction

Idiopathic portal hypertension (IPH) is characterized by portal hypertension due to obstruction or stenosis of the intrahepatic peripheral portal branches [1, 2]. Patients with IPH present with splenomegaly, anemia, and portal hypertension and are free of obstruction of the extrahepatic portal vein or hepatic vein, hematological disease, parasitic disease, granulomatous liver disease, or congenital hepatic fibrosis. IPH does not usually lead to liver cirrhosis [35]. IPH resembles a disease called noncirrhotic portal fibrosis in India or hepatoportal sclerosis/noncirrhotic portal hypertension in Western countries [5, 6]. The etiology of IPH remains unclear. To date, researchers have suggested that IPH may be attributed to intrahepatic peripheral portal vein thrombosis, splenic factors, abnormal autoimmunity, and related factors [711].

Our research team has shown that the connective tissue growth factor (CTGF) gene is expressed specifically in patients with IPH. Additionally, we previously reported that patients with IPH exhibit high levels of CTGF in the blood and overexpression of CTGF mRNA in liver tissues [12]. However, in rats with transient overexpression of CTGF induced by recombinant adenovirus, no changes were observed in liver tissues, despite the expression of fibrosis-related genes [13]. Thus, we aimed to examine the specific expression of other genes in patients with IPH.

Recently, comprehensive gene expression analysis with DNA microarrays has been used in a variety of diseases, such as cancer and immunological disease, and has been employed for pathophysiological analysis of these diseases [1417]. Such analysis enables detection of genes or biological pathways showing significant increases or decreases in expression in the target specimens; this is useful for estimation of the features of the target disease. To the best of our knowledge, no studies have reported comprehensive analysis of gene expression using DNA microarrays in patients with IPH. Thus, in this study, we aimed to perform comprehensive gene expression analysis using DNA microarrays in blood samples from patients with IPH to identify genes showing significant changes in expression. Moreover, we performed biological pathway analysis using the data obtained from DNA microarrays.

2. Materials and Methods

2.1. Patients

Analysis was conducted using blood samples from four patients with IPH satisfying the IPH diagnostic criteria (prepared under the Intractable Hepatobiliary Disease Program of the Ministry of Health, Labour and Welfare, Japan) and four healthy volunteers as controls [18, 19]. The diagnosis of IPH was based on the following criteria: (1) general findings: one or more cytopenia, normal function or mild dysfunction of the liver, or collateral circulation such as upper gastrointestinal varix, (2) imaging analysis: splenomegaly, atrophy of peripheral parenchyma and enlargement of central parenchyma in the liver, irregular hepatic surface, increase in splenic venous flow, abnormality of intrahepatic peripheral portal branches, mutual binding of hepatic veins, or increase of hepatic venous pressure, (3) pathological findings: collapse or stenosis of peripheral branches of intrahepatic portal veins, fibrosis of splenic sinus, or Gamma-Gandy bodies formation in splenic trabeculae, (4) exclusion of other diseases causing portal hypertension such as cirrhosis, extrahepatic occlusion, Budd-Chiari syndrome, blood disease, parasitic disease, granulomatous liver disease, congenital liver fibrosis, chronic viral hepatitis, and primary biliary cirrhosis, among others. These blood samples had been stored without identifying information at the Sample Storage Center of Kyushu University installed by the Portal Hemodynamics Abnormalities Study Group within the framework of the Intractable Hepatobiliary Disease Program of the Ministry of Health, Labour and Welfare, Japan. The samples were randomly chosen for use in this study. This study was conducted in compliance with the ethical principles of  the Declaration of  Helsinki (1964) and was approved by the ethics committees of Osaka City University and Kyushu University. Informed consent was obtained from all participants.

2.2. Sample Adjustment and Microarray Data

All blood samples were stored at −80°C. Total RNA was extracted by ISOGEN-LS (NIPPON GENE Co., Ltd., Tokyo, Japan). RNA levels were measured with a spectrophotometer SmartSpec 3000 (Bio-Rad Laboratories, Inc., Hercules, CA, USA), and the quality of each RNA sample was checked with an Agilent 2100 Bioanalyzer (Agilent Technologies, Inc., Santa Clara, CA, USA). Three samples from patients with IPH and three samples from healthy volunteers (with RNA integrity numbers of ≥ 6) were subjected to microarray analysis. The WT-Ovation Pico RNA Amplification System (NuGEN Technologies, Inc., San Carlos, CA, USA) was employed for cDNA synthesis, conversion into sense chain, fragmentation, and biotin labeling, and a Human Gene ST 1.0 Array (Affymetrix, Inc., Santa Clara, CA, USA) was then used for hybridization of the sense chain cDNA. Array scanning was conducted with a GeneChip 3000 Scanner (Affymetrix, Inc.) to yield image data. The array image data from each sample were converted into files in a format enabling numerical extraction using GeneChip Operating Software (GCOS), a standard data analysis system included with the GeneChip system. Microarray data in this paper have been submitted to National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO) and are accessible through the GEO series accession number GSE69601.

2.3. Cluster Analysis

The signal intensity values for data from the three patients with IPH and the three healthy controls were divided by the mean signal intensity values for all six samples, and the quotient was converted into a logarithm. The calculated relative signal intensity values were presented on a heat map and subjected to cluster analysis using Cluster 3.0 software with the following parameters: k-means; number of clusters, 12; correlation (uncentered) selected for similarity metric; and number of runs, 1000. The results were visualized with Java TreeView software [20].

2.4. Functional Annotation Analysis for Classification of Gene Function

To extract information on the functions of genes identified as differentially expressing using the GeneChip numerical data, the Database for Annotation Visualization and Integrated Discovery (DAVID) tool was used; this database includes the Gene Ontology Database (http://geneontology.org/) and the Kyoto Encyclopedia of Genes and Genomes (http://www.genome.jp/kegg/) [2123]. For each gene, the ratio between the mean values of the IPH group and the healthy control group, that is, the expression ratio (fold-change), and value obtained using -test were calculated. We used the volcano plot method, which reflects both biological significance and statistical significance. By changing the cut-off value of the expression ratio and value to include gene function terms that may be important, four gene groups were defined, including three groups showing a significant reduction in expression in the presence of IPH (Group A: healthy control/IPH > 2.0, ; Group B: healthy control/IPH > 1.5, ; Group C: healthy control/IPH > 1.5, ) and one group showing significant elevation of expression in the presence of IPH (Group D: IPH/healthy control > 1.5, ). For each group, functional annotation chart analysis and functional annotation clustering analysis were performed.

2.5. Biological Interpretation Using Gene Ontology and Network-Based Pathway Analysis

Gene network analysis was performed using Ingenuity iReport (Ingenuity Systems, http://www.ingenuity.com/, Redwood City, CA, USA). This analysis enabled us to identify differentially expressed genes and molecular interactions for the target disease. On the basis of iReport data, differentially expressed genes were focused using -tests and fold changes, followed by biological interpretation with Ingenuity ontology and canonical pathways. For each pathway, biological process, and disease classification, categories were arranged in the order of involvement rate ranking. Furthermore, relevant genes and networks were also explored by Ingenuity Pathways Analysis (IPA; Ingenuity Systems), with reference to the top 10 significant pathways selected by IPA. Biological interpretation was made from these results.

3. Results

3.1. Cluster Analysis

Figure 1 represents the results of cluster analysis. Samples from all three healthy controls showed similar trends (increases or decreases) in gene expression, while samples from the three patients with IPH showed differing trends. No dominant probe group showing a clear difference between the IPH group and the healthy control group was obtained.

3.2. Functional Annotation Analysis of Differentially Expressed Genes

Table 1 shows the number of charts extracted with the functional annotation chart and the number of clusters extracted by functional annotation clustering. The number of charts was smallest (three) in Group D. Among all groups, clusters for Group C (exhibiting the highest enrichment score (3.37)) contained many gene terms related to the immune system, including “lymphocyte activation” and “leukocyte activation” (Figure 2(a)). This finding suggested that immunological abnormalities were involved in IPH. In Group D, the enrichment score was low (0.94), and no gene function terms were found to be specifically increased in the presence of IPH (Figure 2(b)).

3.3. Biological Interpretation Using Network-Based Pathway Analysis

Supplementary Table in Supplementary Material available online at http://dx.doi.org/10.1155/2015/349215 lists all genes selected on the basis of the iReport. For each pathway, biological process, and disease classification, we extracted categories in the top 25 positions of the involvement rate ranking (Tables 24). The highest ranked target in pathway classification was reduced adenosine deaminase (ADA) expression. Additionally, the target at rank 21 in the pathway classification was abnormal purine metabolism. Targets 8 and 15 in pathway classification and targets 1, 2, and 11 in biological process classification were associated with abnormalities in the synthesis/metabolism of arachidonic acid- (AA-) prostaglandin- (PG-) leukotriene (LT). Expression levels of cytochrome P450, family 4, subfamily F, polypeptide 3 (CYP4F3), and glutathione peroxidase 3 (GPX3), which are involved in leukotriene B4 (LTB4) metabolism, were increased, while the expression levels of the PG E receptor 2 (PTGER2) were reduced. Furthermore, target 22 in pathway classification and target 13 in disease classification were associated with abnormal endothelin (ET) signals. The expression of the ET receptor type A (EDNRA) was reduced, while the expression of natriuretic peptide receptor 3 (NPR3) was increased. Additionally, the expression levels of cluster of differentiation 44 (CD44) and transforming growth factor (TGF)-β were reduced. Figure 3(a) shows combination of the first pathway with the second pathway, and Figure 3(b) shows combination of the first pathway with the third and fifth pathways; each of these plots was prepared with reference to the top 10 significant pathways selected by IPA. Reduced expression of the ectonucleoside triphosphate diphosphohydrolase 4 (ENTPD4), ATP-binding cassette, subfamily C, member 1 (ABCC1), ADA, and TGF-β and increased expression of CYP4F3 were noted. Thus, the results from the iReport were also observed within the IPA network.

4. Discussion

IPH is a rare disease. In Japan, IPH was officially added to the list of specific intractable diseases covered by the medical expense subsidy program in January of 2015. IPH has been reported to be associated with autoimmune abnormalities [10, 11, 24]; however, the exact pathophysiology of IPH remains unclear. In this study, we attempted, for the first time, to explore the features of IPH through comprehensive gene analysis with DNA microarrays.

The results of cluster analysis showed that there were no dominant increases/decreases in gene expression in the IPH group as compared to the healthy control group. IPH has been conventionally viewed as a syndrome observed in individuals free of any other disease that can elevate portal pressure. The diagnosis of IPH generally involves ruling out other potential diagnoses, and the pathophysiological features of IPH are not uniform. Thus, these characteristics of IPH may explain the above-mentioned results.

In the analysis of functional annotation charts, many genes showed significantly reduced expression in the IPH group, but only three charts in the IPH group showed significantly increased expression. In functional annotation clustering analysis, the enrichment score of the Group C cluster was the highest (3.37), and this cluster contained many gene terms related to the immune system. Coupled with the observation that some particular trends are found when the enrichment score rises to 2.0 or higher, our data suggested that immunological abnormalities are involved in the pathophysiology of IPH. We also found that the enrichment score in Group D was relatively low (0.94), with no gene function terms showing increased expression in the presence of IPH.

We attempted further clarification of the IPH pathophysiology through network analysis. First, canonical pathway analysis was conducted using the Ingenuity iReport. In pathway classification, reduced expression of ADA was noted. Severe immunodeficiency has been shown to be induced by accumulation of hazardous adenosine metabolites resulting from ADA defects [25, 26]. Thus, reduced expression of ADA may cause induction of autoimmune disease from immunodeficiency, resulting in attack of specific organs. This theory is consistent with the view that immune system abnormalities, particularly suppressor T-cell dysfunction and accompanying autoimmunity, are candidate factors responsible for IPH [24, 27, 28]. ADA dysfunction has also been shown to affect liver differentiation and function [2931]. It is likely that reduced ADA activity serves as a core factor, and combination of ADA with autoimmunity-dependent injury and direct liver injury leads to manifestation of IPH. In this analysis, abnormal purine metabolism was selected. Based on this connection, ADA deficiency induces accumulation of intermediate metabolites due to abnormal nucleic acid metabolism, thereby exerting cytotoxic effects and disrupting lymphocyte differentiation, leading to onset of severe immunodeficiency [32]. Similar to the effects of ADA deficiency, patients with IPH also have abnormal systemic nucleic acid metabolism as a primary feature; this may lead to induction of abnormal lymphocyte differentiation and abnormal function.

The results of the current study showed that the expression levels of CYP4F3 and GPX3, which are known to be involved in degradation of LTB4 and the LTB4 precursor 5-HPETE, were increased. In contrast, the expression of PTGER2 was decreased. LTB4 is known to activate neutrophils [33], and PGE2 is known to have vasodilatory effects [34]. These results may be viewed as compensatory reactions, corresponding to the increased production of LTB4 and PGE2; elevation of these components in the blood may compensate for leukocyte dysfunction and increased blood pressure caused by other factors. Furthermore, analysis of ET signals revealed reduced EDNRA expression and increased NPR3 expression. Decreased expression of the ET receptor leads to vascular contraction [35], while increased expression of NPR3 leads to vascular dilation [36]. Taken together, these findings suggest that NPR3 expression increases as a compensatory reaction to vascular contraction, possibly causing dilation of blood vessels.

When viewed as a whole, the expression of CD44 and TGF-β was reduced. CD44 is known to induce T-cell homing and activate T-helper 1 (Th1) cells, T-helper 2 (Th2) cells, cytotoxic T lymphocytes (CTLs), and natural killer (NK) cells and is thus involved in stem cell proliferation and differentiation [3739]. TGF-β suppresses the proliferation of blood cells and induces the proliferation and activation of mesenchymal cells involved in blood cell differentiation [40, 41]. Thus, in patients with IPH, immune system responses and the differentiation-inducing environment appear to be disturbed.

Furthermore, we explored the network pathway involved in IPH using IPA network analysis. Although the IPA network results did not reveal any additional finding, many genes selected by iReport were also observed in the high-ranked IPA pathways, thus strongly supporting the iReport results.

Taken together, these findings suggest that abnormal systemic nucleic acid metabolism is the first factor involved in the onset of IPH and that this abnormality induces abnormal differentiation and function of immunocompetent cells, although our results are preliminary and require further analysis. In this context, abnormalities in AA-related synthesis/metabolism and ET signals can be viewed as compensatory reactions of blood cells corresponding to changes in the external environment of blood cells, for example, increases/decreases in the levels of certain factors in the serum.

The present study enrolled four IPH patients and three of four samples were used for gene analyses. We cannot deny the possibility that our results did not cover the patterns of whole IPH patients because the sample number was small. In addition, we conducted comprehensive gene expression analysis of blood samples. Splenomegaly is a characteristic symptom in patients with IPH and is often accompanied by histological abnormalities, such as splenic sinus hyperplasia and increased collagen fibers. Thus, analysis of splenic tissue specimens may be useful for clarifying the pathophysiology of IPH. However, unlike the liver, in which percutaneous biopsy is relatively easy, the spleen is difficult to be biopsied or resected in healthy individuals and in patients with IPH during routine clinical practice. To date, no reports have described comprehensive gene expression analysis in patients with IPH. Therefore, additional genetic studies are needed to confirm the biological roles of the molecular pathology of IPH determined in the present study. The identification of genetic abnormalities in IPH may lead to the development of genetic therapy.

5. Conclusion

Blood samples from patients with IPH and healthy individuals as controls were subjected to comprehensive gene analysis with DNA microarrays. In cluster analysis, we did not observe any dominant probe groups showing clear differences between the IPH group and the healthy control group. In functional annotation analysis, the cluster of compromised function among patients with IPH included many gene terms primarily related to the immune system. Network analysis revealed suppression of immune function, abnormal nucleic acid metabolism, and abnormal AA-related synthesis/metabolism. These results suggested that abnormal systemic nucleic acid metabolism was a key factor involved in the onset of IPH and that this abnormality induced abnormal differentiation and function in immunocompetent cells. The abnormalities in AA-related synthesis/metabolism and ET signals may represent compensatory reactions in blood cells, corresponding to changes in the external environment of blood cells, such as increases/decreases in the levels of certain factors in the serum.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors thank Dr. Yoshiyuki Takahara for his assistance of constructive discussion. This study is partially supported by the Health and Labour Sciences Research Grant from the Ministry of Health, Labour and Welfare, Japan, regarding Research on Intractable Diseases, Portal Hemodynamics Abnormalities.

Supplementary Materials

Supplementary Table 1: All genes showing significantly increased or decreased expression in patients with IPH selected on the basis of the iReport.

  1. Supplementary Materials