System Review about Function Role of ESCC Driver Gene KDM6A by Network Biology Approach

Ran, Jihua; Li, Hui; Li, Huiwu

doi:https://doi.org/10.1155/2016/1970904

The Scientific World Journal

On this page

Abstract Introduction Discussion Conclusion Acknowledgments Supplementary Materials References Copyright Related Articles

Research Article | Open Access

Volume 2016 | Article ID 1970904 | https://doi.org/10.1155/2016/1970904

System Review about Function Role of ESCC Driver Gene KDM6A by Network Biology Approach

Jihua Ran,^1,2Hui Li,³and Huiwu Li^1,4

Academic Editor: Manoj K. Kashyap

Received27 Jan 2016

Accepted29 Mar 2016

Published17 May 2016

Abstract

Background. KDM6A (Lysine (K)-Specific Demethylase 6A) is the driver gene related to esophageal squamous cell carcinoma (ESCC). In order to provide more biological insights into KDM6A, in this paper, we treat PPI (protein-protein interaction) network derived from KDM6A as a conceptual framework and follow it to review its biological function. Method. We constructed a PPI network with Cytoscape software and performed clustering of network with Clust&See. Then, we evaluate the pathways, which are statistically involved in the network derived from KDM6A. Lastly, gene ontology analysis of clusters of genes in the network was conducted. Result. The network includes three clusters that consist of 74 nodes connected via 453 edges. Fifty-five pathways are statistically involved in the network and most of them are functionally related to the processes of cell cycle, gene expression, and carcinogenesis. The biology themes of clusters 1, 2, and 3 are chromatin modification, regulation of gene expression by transcription factor complex, and control of cell cycle, respectively. Conclusion. The PPI network presents a panoramic view which can facilitate for us to understand the function role of KDM6A. It is a helpful way by network approach to perform system review on a certain gene.

1. Introduction

Esophageal cancer is the eighth most common malignancy and the sixth cause of cancer deaths worldwide [1, 2]. Esophageal squamous cell carcinoma (ESCC) and adenocancer are the two histologic types that make up for greater than 90 percent of the diagnoses of esophageal cancers [3]. In East Asia, the majority of esophageal cancers are ESCC [2]. As the most of cancers, the pathomechanism of ESCC remains elusive. There is an agreement that cancers arise owing to mutations in a subset of genes that confer growth advantage. Not all mutations can contribute to the development of the cancers. The genes which harbor mutation conferring benefit to survival of tumor cell are called driver genes. The counterpart of driver gene is passenger gene, whose mutation does not contribute to oncogenesis [4]. Defining a driver gene in physiologic terms is easy but identifying which mutations are drivers and which are passengers is more difficult, so numerous statistical methods to identify driver genes have been described. Vogelstein and colleagues identified about 140 driver genes from Catalogue of Somatic Mutations in Cancer database by their unique ruler [5]. Of them, KDM6A is the driver gene related to ESCC.

In eukaryotic cells, DNA is packaged into chromatin whose functional unit is nucleosome. A nucleosome is an octameric structure composed of two histones each of H2A, H2B, H3, and H4 encircled by 147 bp of DNA [6]. Modifications of histones such as methylation, acetylation, phosphorylation, ubiquitination, and sumoylation regulate the structure of chromatin and determine the active or repressive chromatin states. Of them, histone methylation can function in gene activation or repression, depending on which residues are targeted. Methylation of histone H3 on Lysine 4 (H3K4me) is an active chromatin modification, while methylation of histone H3 on Lysine 27 (H3K27me) is associated with repression of gene activity [7]. The PRC2 (polycomb repressive complex 2) mediates transcriptional repression by methylation of histone on H3K27 [8–11]. KDM6A counteracts PRC2 and activates chromatin transcriptionally by demethylation of H3K27 [12]. This is the main biological function of KDM6A. However, KDM6A is vital in a wide array of functions including cell cycle regulation, cell differentiation, and stem cell specification [13–15]. The mutation of KDM6A is involved in various types of cancer across both solid and liquid tumors [16–18].

As we all know, proteins do not function alone while they orchestrate all biological processes by interacting with other proteins. Analysis of PPI network consisting of KDM6A and its interactors facilitates providing more biological insights into KDM6A for us. In this paper, we treat PPI network derived from KDM6A as a conceptual framework and follow it to review its biological function.

2. Method

2.1. Construction of PPI Network Derived from KDM6A

We obtained PPIs from STRING database, a precomputed database for the exploration of protein-protein interactions. The prediction methods in STRING include neighborhood gene fusion, cooccurrence, coexpression, experiments, databases, and text mining. A confidence score can be assigned to each prediction method and every interaction has a final aggregate score. The newest version of STRING 9.1 covers approximately 2.5 million proteins from 630 different organisms [19]. In this study, the interactions restricted to Homo sapiens were downloaded. In order to avoid false link in greatest extent, we set the score of inclusion criteria larger than 0.9. We constructed a network that consists of not only the direct PPI neighbors of KDM6A but also their secondary neighbors. The network was constructed using Cytoscape, a popular and highly versatile software platform for the analysis, operation, and visualization of large networks [20].

2.2. Construction of Cellular Pathways Database

All of the pathways with their gene members were downloaded from an integrated pathway database, Molecular Signatures Database (MSigDB) [21], which is a large collection of annotated functional gene sets. There are 880 canonical pathways with 6804 proteins members in the database, including the metabolic and signaling pathways collected from BioCarta (http://www.biocarta.com/), KEGG [22], and Reactome [23].

2.3. Identification of Pathways Involved in PPI Network Derived from KDM6A

To examine the evidence of association of a given pathway with PPI network derived from KDM6A, Fisher’s exact test based on the cumulative hypergeometric distribution was employed. The value was calculated to evaluate statistical significance of a given pathway by the formula as follows: In this formula, represents the total number of proteins in the background population, represents the number of proteins in the PPI network derived from KDM6A, and denotes the number of proteins within the given pathways. The number of proteins that overlapped with both proteins in the PPI network and this pathway is denoted as . In this study, a pathway is considered statistically associating with the PPI network derived from KDM6A under the condition of its value being less than 0.05.

2.4. Ontology Analysis of Proteins in PPI Network Derived from KDM6A

Usually, gene ontology analysis is mainly conducted in individual gene. However, in this study, we perform ontology enrichment analysis based on a group of genes. As a web tool, NOA not only is able to perform analysis based on a group of genes but also takes interactions between genes into consideration [24]. Every group of genes corresponds to the proteins in the clusters of PPI network derived from KDM6A. There are many clustering methods that have been used to highlight groups of densely connected nodes [25]. Clusters in networks can offer mechanistic hypotheses of disease because they are highly interconnected molecular complexes or signaling pathways [26]. In this study, clusters were found with Cytoscape plugin called Clust&See dedicated to the identification, visualization, and analysis of clusters extracted from PPI network [27].

3. Result and Discussion

3.1. General Description of the PPI Network Derived from KDM6A

The network consists of 74 nodes connected via 453 edges. KDM6A has seven direct neighbors, which are RBBP5, WDR5, ASH2L, MLL2, PAXIP1, NCOA6, and RBL2, respectively (Figure 1). KDM6A contains a tetratricopeptide motif predicted to mediate protein-protein interactions [28] and is not only a member of a stable multiprotein complex that demethylates H3K27me3 but also a member of the MLL2 H3K4 methyltransferase complex which can facilitate gene expression [29, 30]. Its catalytic activity has been linked to regulation of homeobox (HOX) and RB transcriptional networks [31]. Different protein partners modulate the recruitment of KDM6A to specific chromatin regions to target specific genes [32, 33]. There are three clusters whose nodes are distinguished with different colors in the network (Figure 1). Every cluster has its biological theme detailed later and there are links among these clusters.

3.2. Pathways Involved in PPI Network Derived from KDM6A

Of eight hundred and eighty (880) pathways, fifty-five (55) pathways are statistically involved in PPI network derived from KDM6A (, see Table in Supplementary Material available online at http://dx.doi.org/10.1155/2016/1970904). These fifty-five pathways can be catalogue to eight classes including cell cycle, gene expression, lipid metabolism, cancer, apoptosis, signal transduction, development, and DNA repair. Most of pathways are related to cell cycle, gene expression, and cancer (Figure 2).

There are two points in current mainstream discipline about KDM6A: (i) as a general factor, KDM6A activates gene transcription during development. HOX genes encode transcription factors that regulate embryogenesis and guide tissue differentiation [34]. KDM6A removes H3K27me3 from HOX genes to restore their activity and control HOX gene expression [35]. (ii) KDM6A acts with methyltransferase complex to facilitate gene expression for regulating transcriptional networks of RB [36].

The evidence indicates that KDM6A evolved in pathways such as inflammation, apoptosis, cell cycle, and DNA repair related. KDM6A is implicated in IL-4 mediated transcriptional activation of the arachidonate 15-lipoxygenase-1 (ALOX15) gene. ALOX15 oxygenates polyunsaturated fatty acids and biomembranes, which generate multiple lipid signaling mediators involved in inflammation [37]. KDM6A is required for hormone-mediated transcriptional regulation of apoptosis and autophagy genes in Drosophila salivary glands [38]. DNA methylation affects the expression of genes involved in cell cycle checkpoint, apoptosis, and DNA repair [39].

3.3. Ontology Analysis of the Protein Interactions in the Three Clusters

We manually choose top three terms ranked by value generated by analysis tool NOA (Table 1). The theme of cluster 1 is chromatin modification by methyl transfer of histone. It is the main function of KDM6A as previously mentioned. As part of an H3K4-methyltransferase complex containing MLL2, PTIP, ASC2, ASH2, RBQ3, and WDR5, KDM6A can promote H3K4 methylation [30]. It indicates that KDM6A and methyltransferase complex concert with active chromatin modification by demethylation of histone on H3K27 and methylation on K3K4. Clinically, KDM6A or MLL4 was associated with poor prognosis in patients with breast cancer. KDM6A interacts with a c-terminal region of MLL4 and coordinates regulation of cotarget gene expression [40].

The proteins in cluster 2 mainly compose the transcription factor complex. The direct neighbor of KDM6A, NCOA6 (nuclear receptor coactivator 6, NCOA6), binds nuclear receptors and stimulates the transcriptional activities in a hormone-dependent fashion. Mass spectrometry analysis demonstrated that PTIP associated with ASH2L, RBBP5, WDR5, NCOA6, and KDM6A [29]. It suggests that there exist interactions among NCOA6, KDM6A, and other methyltransferases.

The nodes in cluster 3 are RBL2 (retinoblastoma-like 2, RBL2) and its interactors which control cell cycle. RBL2 is a key regulator of entry into cell division, directly involved in heterochromatin formation by maintaining overall chromatin structure through stabilizing histone methylation. Even though there is no direct experimental evidence that KDM6A interacts with RBL2, a study found that ectopic expression of KDM6A enhanced the expression of retinoblastoma tumor suppressor gene RB and its related gene RBL2 [36].

4. Conclusion and Prospect

PPI network derived from KDM6A presents a panoramic view which can facilitate the understanding of the function role of KDM6A for us. The results from both sides of identification of pathway and ontology analysis of network confirm and complement one another. It is a helpful way by network approach to perform system review on a certain gene. Considering that these could be nonspecific interactions, the network can only be regarded as a guide map. Some links between the proteins in the network need to be validated by experiment to explore the biology role of KDM6A in the progression of carcinogenesis. For example, we need to investigate the expression pattern of these genes at mRNA or protein level and the possible regulatory relationships between these genes in clinical ESCC specimen or esophageal cancer cell line. Similar studies in other cancers were completed and reported [40, 41].

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

The present study was supported by a grant from the National Natural Science Foundation of China (Grant no. 81460359).

Supplementary Materials

List of pathways involved in PPI network derived from KDM6A.

Supplementary Material

References

F. Kamangar, G. M. Dores, and W. F. Anderson, “Patterns of cancer incidence, mortality, and prevalence across five continents: defining priorities to reduce cancer disparities in different geographic regions of the world,” Journal of Clinical Oncology, vol. 24, no. 14, pp. 2137–2150, 2006.
View at: Publisher Site | Google Scholar
E. A. Mathé, H. N. Giang, E. D. Bowman et al., “MicroRNA expression in squamous cell carcinoma and adenocarcinoma of the esophagus: associations with survival,” Clinical Cancer Research, vol. 15, no. 19, pp. 6192–6200, 2009.
View at: Publisher Site | Google Scholar
J. M. Daly, W. A. Fry, A. G. Little et al., “Esophageal cancer: results of an American College of Surgeons patient care evaluation study,” Journal of the American College of Surgeons, vol. 190, no. 5, pp. 562–572, 2000.
View at: Publisher Site | Google Scholar
C. Greenman, P. Stephens, R. Smith et al., “Patterns of somatic mutation in human cancer genomes,” Nature, vol. 446, no. 7132, pp. 153–158, 2007.
View at: Publisher Site | Google Scholar
B. Vogelstein, N. Papadopoulos, V. E. Velculescu, S. Zhou, L. A. Diaz Jr., and K. W. Kinzler, “Cancer genome landscapes,” Science, vol. 339, no. 6127, pp. 1546–1558, 2013.
View at: Publisher Site | Google Scholar
R. Margueron and D. Reinberg, “Chromatin structure and the inheritance of epigenetic information,” Nature Reviews Genetics, vol. 11, no. 4, pp. 285–296, 2010.
View at: Publisher Site | Google Scholar
C. Martin and Y. Zhang, “The diverse functions of histone lysine methylation,” Nature Reviews Molecular Cell Biology, vol. 6, no. 11, pp. 838–849, 2005.
View at: Publisher Site | Google Scholar
B. Czermin, R. Melfi, D. McCabe, V. Seitz, A. Imhof, and V. Pirrotta, “Drosophila enhancer of Zeste/ESC complexes have a histone H3 methyltransferase activity that marks chromosomal Polycomb sites,” Cell, vol. 111, no. 2, pp. 185–196, 2002.
View at: Publisher Site | Google Scholar
R. Cao, L. Wang, H. Wang et al., “Role of histone H3 lysine 27 methylation in polycomb-group silencing,” Science, vol. 298, no. 5595, pp. 1039–1043, 2002.
View at: Publisher Site | Google Scholar
J. Muller, C. M. Hart, N. J. Francis et al., “Histone methyltransferase activity of a Drosophila Polycomb group repressor complex,” Cell, vol. 111, no. 2, pp. 197–208, 2002.
View at: Publisher Site | Google Scholar
A. Kuzmichev, K. Nishioka, H. Erdjument-Bromage, P. Tempst, and D. Reinberg, “Histone methyltransferase activity associated with a human multiprotein complex containing the enhancer of zeste protein,” Genes and Development, vol. 16, no. 22, pp. 2893–2905, 2002.
View at: Publisher Site | Google Scholar
B. Schuettengruber, D. Chourrout, M. Vervoort, B. Leblanc, and G. Cavalli, “Genome regulation by polycomb and trithorax proteins,” Cell, vol. 128, no. 4, pp. 735–745, 2007.
View at: Publisher Site | Google Scholar
G. L. Sen, D. E. Webster, D. I. Barragan, H. Y. Chang, and P. A. Khavari, “Control of differentiation in a self-renewing mammalian tissue by the histone demethylase JMJD3,” Genes and Development, vol. 22, no. 14, pp. 1865–1870, 2008.
View at: Publisher Site | Google Scholar
J. K. Wang, M.-C. Tsai, G. Poulin et al., “The histone demethylase UTX enables RB-dependent cell fate control,” Genes and Development, vol. 24, no. 4, pp. 327–332, 2010.
View at: Publisher Site | Google Scholar
T. Satoh, O. Takeuchi, A. Vandenbon et al., “The Jmjd3-Irf4 axis regulates M2 macrophage polarization and host responses against helminth infection,” Nature Immunology, vol. 11, no. 10, pp. 936–944, 2010.
View at: Publisher Site | Google Scholar
G. L. Dalgliesh, K. Furge, C. Greenman et al., “Systematic sequencing of renal carcinoma reveals inactivation of histone modifying genes,” Nature, vol. 463, no. 7279, pp. 360–363, 2010.
View at: Publisher Site | Google Scholar
Y. Gui, G. Guo, Y. Huang et al., “Frequent mutations of chromatin remodeling genes in transitional cell carcinoma of the bladder,” Nature Genetics, vol. 43, no. 9, pp. 875–878, 2011.
View at: Publisher Site | Google Scholar
B. G. Mar, L. Bullinger, E. Basu et al., “Sequencing histone-modifying enzymes identifies UTX mutations in acute lymphoblastic leukemia,” Leukemia, vol. 26, no. 8, pp. 1881–1883, 2012.
View at: Publisher Site | Google Scholar
A. Franceschini, D. Szklarczyk, S. Frankild et al., “STRING v9.1: protein-protein interaction networks, with increased coverage and integration,” Nucleic Acids Research, vol. 41, no. 1, pp. D808–D815, 2013.
View at: Publisher Site | Google Scholar
P. Shannon, A. Markiel, O. Ozier et al., “Cytoscape: a software environment for integrated models of biomolecular interaction networks,” Genome Research, vol. 13, no. 11, pp. 2498–2504, 2003.
View at: Publisher Site | Google Scholar
A. Liberzon, A. Subramanian, R. Pinchback, H. Thorvaldsdóttir, P. Tamayo, and J. P. Mesirov, “Molecular signatures database (MSigDB) 3.0,” Bioinformatics, vol. 27, no. 12, pp. 1739–1740, 2011.
View at: Publisher Site | Google Scholar
M. Kanehisa, S. Goto, M. Furumichi, M. Tanabe, and M. Hirakawa, “KEGG for representation and analysis of molecular networks involving diseases and drugs,” Nucleic Acids Research, vol. 38, no. 1, pp. D355–D360, 2009.
View at: Publisher Site | Google Scholar
D. Croft, G. O'Kelly, G. Wu et al., “Reactome: a database of reactions, pathways and biological processes,” Nucleic Acids Research, vol. 39, no. 1, pp. D691–D697, 2011.
View at: Publisher Site | Google Scholar
J. Wang, Q. Huang, Z.-P. Liu et al., “NOA: a novel Network Ontology Analysis method,” Nucleic Acids Research, vol. 39, article e87, 2011.
View at: Publisher Site | Google Scholar
T. Aittokallio and B. Schwikowski, “Graph-based methods for analysing networks in cell biology,” Briefings in Bioinformatics, vol. 7, no. 3, pp. 243–255, 2006.
View at: Publisher Site | Google Scholar
T. Ideker and R. Sharan, “Protein networks in disease,” Genome Research, vol. 18, no. 4, pp. 644–652, 2008.
View at: Publisher Site | Google Scholar
L. Spinelli, P. Gambette, C. E. Chapple et al., “Clust&See: a Cytoscape plugin for the identification, visualization and manipulation of network clusters,” BioSystems, vol. 113, no. 2, pp. 91–95, 2013.
View at: Publisher Site | Google Scholar
T. Swigut and J. Wysocka, “H3K27 Demethylases, at Long Last,” Cell, vol. 131, no. 1, pp. 29–32, 2007.
View at: Publisher Site | Google Scholar
Y.-W. Cho, T. Hong, S. Hong et al., “PTIP associates with MLL3- and MLL4-containing histone H3 lysine 4 methyltransferase complex,” The Journal of Biological Chemistry, vol. 282, no. 28, pp. 20395–20406, 2007.
View at: Publisher Site | Google Scholar
I. Issaeva, Y. Zonis, T. Rozovskaia et al., “Knockdown of ALR (MLL2) reveals ALR target genes and leads to alterations in cell adhesion and growth,” Molecular and Cellular Biology, vol. 27, no. 5, pp. 1889–1903, 2007.
View at: Publisher Site | Google Scholar
J. Van der Meulen, F. Speleman, and P. Van Vlierberghe, “The H3K27me3 demethylase UTX in normal development and disease,” Epigenetics, vol. 9, no. 5, pp. 658–668, 2014.
View at: Publisher Site | Google Scholar
K. Agger, P. A. C. Cloos, J. Christensen et al., “UTX and JMJD3 are histone H3K27 demethylases involved in HOX gene regulation and development,” Nature, vol. 449, no. 7163, pp. 731–734, 2007.
View at: Publisher Site | Google Scholar
G. G. Welstead, M. P. Creyghton, S. Bilodeau et al., “X-linked H3K27me3 demethylase Utx is required for embryonic development in a sex-specific manner,” Proceedings of the National Academy of Sciences of the United States of America, vol. 109, no. 32, pp. 13004–13009, 2012.
View at: Publisher Site | Google Scholar
S. D. Weatherbee, G. Halder, J. Kim, A. Hudson, and S. Carroll, “Ultrabithorax regulates genes at several levels of the wing-patterning hierarchy to shape the development of the Drosophila haltere,” Genes and Development, vol. 12, no. 10, pp. 1474–1482, 1998.
View at: Publisher Site | Google Scholar
F. Lan, P. E. Bayliss, J. L. Rinn et al., “A histone H3 lysine 27 demethylase regulates animal posterior development,” Nature, vol. 449, no. 7163, pp. 689–694, 2007.
View at: Publisher Site | Google Scholar
M. Terashima, A. Ishimura, M. Yoshida, Y. Suzuki, S. Sugano, and T. Suzuki, “The tumor suppressor Rb and its related Rbl2 genes are regulated by Utx histone demethylase,” Biochemical and Biophysical Research Communications, vol. 399, no. 2, pp. 238–244, 2010.
View at: Publisher Site | Google Scholar
H. Han, D. Xu, C. Liu, H.-E. Claesson, M. Björkholm, and J. Sjöberg, “Interleukin-4-mediated 15-lipoxygenase-1 trans-activation requires UTX recruitment and H3K27me3 demethylation at the promoter in A549 cells,” PLoS ONE, vol. 9, no. 1, article e85085, 2014.
View at: Publisher Site | Google Scholar
D. Denton, M. T. Aung-Htut, N. Lorensuhewa et al., “UTX coordinates steroid hormone-mediated autophagy and cell death,” Nature Communications, vol. 4, article 2916, 2013.
View at: Publisher Site | Google Scholar
M. Toyota, H. Suzuki, T. Nishizaka, A. Sato, and T. Yamashita, “Molecular mechanisms involved in epigenetic alterations in cancer,” Gan To Kagaku Ryoho, vol. 37, no. 9, pp. 1650–1653, 2010.
View at: Google Scholar
J.-H. Kim, A. Sharma, S. S. Dhar et al., “UTX and MLL4 coordinately regulate transcriptional programs for cell proliferation and invasiveness in breast cancer cells,” Cancer Research, vol. 74, no. 6, pp. 1705–1717, 2014.
View at: Publisher Site | Google Scholar
X. L. Wang, C. L. Li, L. H. Yang, Z. M. Jiang, Y. T. Gui, and Z. M. Cai, “Clinical significance of high expression with UTX in renal cell carcinoma,” Beijing Da Xue Xue Bao, vol. 46, no. 6, pp. 926–930, 2014.
View at: Google Scholar

Copyright

Copyright © 2016 Jihua Ran 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

1539

Downloads

904

Citations