The Scientific World Journal

Volume 2014 (2014), Article ID 431792, 14 pages

http://dx.doi.org/10.1155/2014/431792

Research Article

Exploration of Potential Roles of a New LOXL2 Splicing Variant Using Network Knowledge in Esophageal Squamous Cell Carcinoma

Bing-Li Wu,¹ Guo-Qing Lv,¹ Hai-Ying Zou,¹ Ze-Peng Du,² Jian-Yi Wu,¹ Pi-Xian Zhang,¹ Li-Yan Xu,³ and En-Min Li¹

¹Department of Biochemistry and Molecular Biology, Shantou University Medical College, Shantou, Guangdong 515041, China

²Department of Pathology, Shantou Central Hospital, Affiliated Shantou Hospital of Sun Yat-sen University, Shantou, Guangdong 515041, China

³Institute of Oncologic Pathology, Shantou University Medical College, Shantou, Guangdong 515041, China

Received 28 March 2014; Revised 1 July 2014; Accepted 14 July 2014; Published 31 August 2014

Academic Editor: Rui Medeiros

Copyright © 2014 Bing-Li Wu 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.

Abstract.

LOXL2 (lysyl oxidase-like 2), an enzyme that catalyzes oxidative deamination of lysine residue, is upregulated in esophageal squamous cell carcinoma (ESCC). A LOXL2 splice variant LOXL2-e13 and its wild type were overexpressed in ESCC cells followed by microarray analyses. In this study, we explored the potential role and molecular mechanism of LOXL2-e13 based on known protein-protein interactions (PPIs), following microarray analysis of KYSE150 ESCC cells overexpressing a LOXL2 splice variant, denoted by LOXL2-e13, or its wild-type counterpart. The differentially expressed genes (DEGs) of LOXL2-WT and LOXL2-e13 were applied to generate individual PPI subnetworks in which hundreds of DEGs interacted with thousands of other proteins. These two DEG groups were annotated by Functional Annotation Chart analysis in the DAVID bioinformatics database and compared. These results found many specific annotations indicating the potential specific role or mechanism for LOXL2-e13. The DEGs of LOXL2-e13, comparing to its wild type, were prioritized by the Random Walk with Restart algorithm. Several tumor-related genes such as ERO1L, ITGA3, and MAPK8 were found closest to LOXL2-e13. These results provide helpful information for subsequent experimental identification of the specific biological roles and molecular mechanisms of LOXL2-e13. Our study also provides a work flow to identify potential roles of splice variants with large scale data.

1. Introduction

The lysyl oxidase (LOX) family, which is composed of five enzymes (LOX and LOXL1/2/3/4), catalyzes oxidative deamination of lysine residues in their protein substrates, generating highly reactive aldehyde residues that initiate inter- and intramolecular cross-linkages [1]. LOX family members are present in several human tissues, including the placenta, heart, lung, kidney, and pancreas [2–6], and are critical for multiple biological functions, such as growth, development, senescence, chemotaxis, and cell mobility [7]. LOXL2 has been emphasized in recent years because of its critical roles in carcinomas. Upregulation of LOXL2 has been detected in many tumor cell lines or clinical samples and also closely correlates with tumor invasion and metastasis [8–11]. LOXL2 protein distributes in either extracellularly or intracellularly [12]. Secreted LOXL2 is able to mediate extracellular matrix remodeling by upregulation of tissue inhibitor of metalloproteinase-1 (TIMP-1) and matrix metalloproteinase-9 (MMP-9) [13]. Intracellular LOXL2 is able to positively regulate the epithelial-mesenchymal transition (EMT) inducer Snail by enhancing Snail stability and functional activity and promoting EMT and tumor progression through downregulation of E-cadherin [14]. Moreover, LOXL2 modulates focal adhesions, tight junctions, and cell polarity complexes in basal breast carcinoma cells through activation of the FAK signaling pathway [15]. The mechanisms of intracellular LOXL2 action are not yet fully known. Recently, LOXL2 has been found to be associated with chromatin and reported to be involved in histone H3 deamination, a novel function that is dependent on the LOXL2 catalytic domain [16]. These analyses suggest the functions of LOXL2 in carcinoma are multifaceted and complicated. Therefore, delineation of LOXL2 function will provide a broad understanding of carcinogenesis.

In our previous study, LOXL2 was found to be overexpressed in esophageal squamous cell carcinoma (ESCC) cell lines and clinical samples and was significantly associated with lymph node metastasis [17]. Immunohistochemistry results showed the expression level of LOXL2 in ESCC is decreased in the nucleus but increased in the cytoplasm. Overall survival rates of ESCC patients with decreased nuclear expression or increased cytoplasmic expression of LOXL2 are significantly lower than those of the patients with the reverse expression pattern [17]. In a recent study, we identified a splice variant of LOXL2 lacking exon 13, denoted by LOXL2-e13, which is also expressed in ESCC cell lines and clinical samples [18]. To reveal the biological roles and molecular mechanisms of LOXL2 and its variants, we overexpressed wild-type LOXL2 (LOXL2-WT) and LOXL2-e13 in ESCC KYSE150 cell line and analyzed the mRNA profiles by the PrimeView Human Gene Expression Array (Affymetrix Corp., St Clara, CA, USA).

Hundreds and thousands of interactions between either extracellular or intracellular proteins compose a network. With recent advances in high-throughput technologies in protein-protein interactions (PPIs), network knowledge can give rise to understanding the biological function and dynamic behavior of cellular systems, generating new biological hypotheses and providing important clues for experimental verification [19–21]. In this study, two PPI subnetworks were generated by mapping DEGs of LOXL2-WT and LOXL2-e13 to the human PPI dataset. These DEGs were annotated by Functional Annotation Chart in the DAVID bioinformatics database. Annotations were compared to reveal the potentially specific roles or mechanisms of LOXL2-e13. This analysis can provide important clues for the future identification of specific roles of LOXL2-e13 from the viewpoint of system analysis.

2. Materials and Methods

2.1. Differentially Expressed Genes

The detailed manipulations of overexpression and microarray were prepared in another manuscript [18]. Briefly, LOXL2-WT and LOXL2-e13 were cloned into the pcDNA3.0 plasmid and overexpressed by transfection into ESCC KYSE150 cells (a generous gift from Professor Mingzhou Guo, Department of Gastroenterology and Hepatology, Chinese PLA General Hospital, Beijing, China). Transfection with an empty plasmid was used as a control. The mRNA expression profiles were analyzed by PrimeView Human Gene Expression Array (Affymetrix Corp., St Clara, CA, USA), and the raw data was normalized by the RMA algorithm. The mRNA expression profiles have been submitted to the GEO database (http://www.ncbi.nlm.nih.gov/geo) under accession number of GSE53645. The DEGs of LOXL2-WT versus empty plasmid (LOXL2-WT-DEGs) and LOXL2-e13 versus empty plasmid (LOXL2-e13-DEGs) were differentiated using a threshold of 1.5-fold change. Moreover, e13-WT-DEGs were also obtained by comparing the expression profile of LOXL2-e13 to that of LOXL2-WT using the same threshold.

2.2. Protein-Protein Interaction Subnetwork Generation

Human Protein Reference Database (HPRD) (http://www.hprd.org) is a database of human proteins, including protein-protein interactions, posttranslational modifications, enzyme-substrate relationships, and disease associations [22]. The Biological General Repository for Interaction Datasets (BioGRID) is a freely accessible database of physical and genetic interactions for various species, which is available at http://www.thebiogrid.org [23]. All PPI data were collected manually from the published literature. These two datasets are popular in PPI network research or other high-throughput data analyses because of their reliability and up-to-date release [24, 25]. The newest PPI data version of Homo sapiens species from these two databases was downloaded for our study, and we integrated them by removing the redundant interactions. There were 18595 unique proteins and 174552 interactions in this integrated PPI dataset, and we regarded it as the parent PPI network.

First, we used Cytoscape to map the LOXL2-WT-DEGs and LOXL2-e13-DEGs PPI datasets to the HPRD&BioGRID parent PPI network to generate the respective PPI subnetworks. Cytoscape (http://www.cytoscape.org) is a free software for visualizing, modeling, and analyzing molecular and genetic interaction networks and has been widely applied in the field of “omics” research (e.g., genomics, transcriptomics, proteomics, and metabolomics) [26]. In Cytoscape, the PPI network is presented by a graph in which each protein is represented by a node, and each interaction between two nodes is represented by an edge. To increase the reliability and limit the protein perturbation to a certain level, the DEG PPI subnetwork construction was limited to the first neighbor proteins of each DEG.

2.3. Network Topological Parameter Analysis

In order to gain insight into the organization and structure of the PPI network, NetworkAnalyzer (http://med.bioinf.mpi-inf.mpg.de/netanalyzer) was applied to analyze the network topological parameters [27]. In the study of networks, the degree of a node is the number of its connections to other nodes, and the degree distribution is the probability distribution of these degrees over the whole network. The node-degree distribution of PPI subnetworks follows a power law, one of the most important network topological characteristics, and was analyzed as we previously described [28]. Briefly, the edges in all networks were treated as undirected. Node-degree distribution is the number of nodes with a degree . By fitting a line on datasets, the pattern of their dependencies can be visualized. NetworkAnalyzer considers only data points with positive coordinate values for fitting the line where the power law curve of the form is transformed into a linear model and the value (coefficient of determination) provides a measure of how well the data points fit to the curve. An value closer to 1 indicates a better fit.

2.4. Functional Annotation of Differentially Expressed Genes and Comparison

Functional Annotation Chart in the DAVID bioinformatics database (http://david.abcc.ncifcrf.gov/) is a tool to identify the enriched biological terms associated with a large gene list [29]. Compared to other similar enrichment analysis tools, Functional Annotation Chart has an extended annotation coverage, increasing from only GO to over 40 annotation categories, including protein sequence features, protein-protein interactions, protein functional domains, disease associations, pathways, homology, gene functional summaries, gene tissue expression, and literature. For example, INTERPRO annotation provides functional analysis of protein sequences by classifying them into families and predicting the presence of domains and important sites [30]. SMART annotation (Simple Modular Architecture Research Tool) allows the identification and annotation of genetically mobile domains and the analysis of domain architecture [31]. The LOXL2-e13-DEGs and LOXL2-WT-DEGs were submitted to the DAVID bioinformatics database and their Functional Annotation Chart results obtained. The terms from the Functional Annotation Chart results with were visualized by the Enrichment Map plugin in Cytoscape, which is a network-based visualization method for gene-set enrichment results [32]. The two Functional Annotation Chart results of LOXL2-e13-DEGs and LOXL2-WT-DEGs were compared and displayed.

2.5. Prioritization Analyses of e13-WT-DEGs

We applied the Random Walk with Restart (RWR) algorithm to analyze the prioritization of e13-WT-DEGs when considering their relationships with LOXL2-e13. RWR is an algorithm for graph analysis, for example, a PPI network, in which, at every tick time, a random walker has a chance to get back to one or more start nodes, from any current node, with a fixed, common, and constant probability [33, 34]. For example, when node(A) is a start node, an RWR result with node(A) provides the stationary probabilities from node(A) to all nodes of the given graph. These probabilities can be considered to be the affinity or proximity from node(A) to individual nodes. RWR is formally defined as the following equation:

where is the restart probability, is the column-normalized adjacency matrix of the network graph, and is a vector of size equal to the number of nodes in the graph where the ith element holds the probability of being at node at time step . In this study, RWR was executed with our custom program in the e13-WT-DEGs PPI subnetwork with LOXL2 set as the seed node.

3. Results

3.1. Protein-Protein Interaction Subnetwork of Differentially Expressed Genes

Using a 1.5-fold change as the threshold, we obtained 349 DEGs, including 217 upregulated genes and 132 downregulated genes, from comparing LOXL2-e13 versus empty plasmid, while there were 606 upregulated genes and 351 downregulated genes from comparison of LOXL2-WT versusempty plasmid. A single protein could not work alone to achieve so many biological effects but must cooperate with other proteins through their physical interactions. To find how many and what kinds of proteins are connected with the DEGs, two PPI subnetworks were constructed by mapping the LOXL2-WT-DEGs and LOXL2-e13-DEGs to the HPRD&BioGRID parent PPI network to generate their respective PPI subnetworks.

The LOXL2-e13-DEGs PPI subnetwork contained 4768 nodes and 67769 edges, including 289 DEGs (Figure 1(a)). We also found, among the interactions between DEGs, a small protein cluster containing 84 nodes and 103 edges (Figure 1(b)). The LOXL2-WT-DEGs PPI subnetwork was composed of 6429 nodes and 97233 edges with 821 DEGs (Figure 1(c)). The interactions between LOXL2- and WT-DEGs were also found. To our surprise, there was a huge protein cluster containing 318 DEGs. These two PPI subnetworks indicated that overexpression of LOXL2-WT or LOXL2-e13 greatly perturbed the PPI network in ESCC cells due to DEGs interacting with hundreds and thousands of proteins to achieve the biological consequences of the LOX2 protein itself. To our great interest, a large number of interactions between DEGs were found, indicating these DEGs might form protein complexes in specific temporal or spatial patterns, and their upregulation or downregulation might favor their counterpart to conduct its biological functions.

(a)

(b)

(c)

(d)

Figure 1: DEGs were mapped to the HPRD&BioGRID parent PPI network to generate PPI subnetworks. (a) PPI subnetworks of LOXL2-e13-DEGs. (b) PPI subnetworks of LOXL2-WT-DEGs. (c) Interactions between LOXL2-e13-DEGs. (d) Interactions between LOXL2-WT- DEGs. Nodes are labeled by different colors to indicate the expression trend of proteins. Green nodes represent proteins encoded by downregulated genes, while red nodes represent proteins encoded by upregulated genes. Interacting proteins without significantly different expression are represented as blue nodes.

3.2. Topological Parameters of PPI Subnetworks

Network topological parameters are important characteristics for evaluating the networks. In this study, PPI subnetwork topological parameters were analyzed by NetworkAnalyzer. The distributions of node degree approximately followed power law distributions with and 0.852, respectively (Figure 2). Thus, these two PPI subnetworks were characterized to be scale-free, one of the most important parameters to recognize a true complex biological network [35]. The other topological parameters of these two subnetworks, such as clustering coefficient, network centralization, and network density, are shown in Table 1.

Table 1: Network parameters of the LOXL2-e13-DEGs and LOXL2-WT-DEGs PPI subnetwork.

PPI subnetwork

Correlation
Clustering coefficient
Network centralization
Network density

LOXL2-e13

0.852
0.754
0.261
0.680
0.006

LOXL2-WT

0.863
0.618
0.242
0.681
0.005

(a)

(b)

Figure 2: Power law fit of node-degree distribution for PPI subnetworks. The node degree is represented on the -axis and the number of nodes with is represented on the -axis. The graph displays a decreasing trend of degree distribution with an increase in the number of links, indicating scale-free topology. (a) Node-degree distribution of the LOXL2-e13-DEGs PPI subnetwork. (b) Node-degree distribution of LOXL2-WT-DEGs PPI subnetwork.

3.3. Annotations and Comparison of DEGs

To understand the function of these DEGs beyond the traditional gene listing and GO annotation, we submitted the LOXL2-e13-DEGs and LOXL2-WT-DEGs to Functional Annotation Chart analysis in the DAVID bioinformatics database. The Functional Annotation Chart results with were visualized and displayed as a network by Enrichment Map.

In Figure 3, the functional annotation terms were represented as nodes. The node with more significance was larger. Edge width was defined by the overlap coefficient between these categories (overlap coefficient cut-off was set as 0.6). The more genes shared by two nodes, the wider their edge width. For LOXL2-WT-DEGs, there were 166 terms from GO categories (e.g., biological process, cellular component, and molecular function) (Figure 3(a)). Except for these, another 96 terms were included from the following annotation categories: 13 INTERPRO, 3 SMART, 31 SP_PIR_KEYWORDS, 38 UP_SEQ_FEATURE, and 5 KEGG_PATHWAY. For the LOXL2-e13-DEGs, in addition to the 71 terms from GO categories, the Functional Annotation Chart results also contained another 51 terms from the following annotation categories: 12 INTERPRO, 2 SMART, 17 SP_PIR_KEYWORDS, 8 UP_SEQ_FEATURE, and 5 KEGG_PATHWAY (Figure 3(b)). These results provide more information than merely using GO enrichment alone, including protein sequence features, protein domains, and pathways.

(a)

(b)

(c)

(d)

Figure 3: Visualization of Functional Annotation Chart analysis for DEGs and their comparison. (a) Visualization of Functional Annotation Chart analysis of LOXL2-e13-DEGs. (b) Visualization of Functional Annotation Chart analysis of LOXL2-WT-DEGs. (c) The unique Functional Annotation Chart of LOXL2-e13-DEGs. (d) The unique Functional Annotation Chart of LOXL2-WT-DEGs.

To illustrate the specific characters or functions of LOXL2-e13 in comparison to wild type, we compared and displayed the differences of the two Functional Annotation Chart results. The number of unique annotation terms for LOXL2-e13 was 89, while that for LOXL2-WT was 229 (Figures 3(c) and 3(d)). From INTERPROT annotation, the protein sequences of LOXL2-e13-DEGs were mainly comprised of domains for “IPR007125:Histone core,” “IPR009072:Histone-fold,” and “IPR000558:Histone H2B.” LOXL2-e13-DEGs were only annotated by the terms for “IPR007872:Zinc finger, DPH-type.” To our surprise, LOXL2-WT-DEGs proteins showed several different zinc annotations, including “IPR007087:Zinc finger, C2H2-type,” “IPR015880:Zinc finger, C2H2-like,” and “IPR013087:Zinc finger, C2H2-type/integrase, DNA-binding” from INTERPROT; “GO:0008270~zinc ion binding” from GO “molecular function”; “zinc-finger,” “zinc finger region : C2H2-type (1–17), C2H2-type 19, and C2H2-type 21” from SP_PIR_KEYWORDS annotation. Other sequence feature annotations from PIR_SUPERFAMILY indicated LOXL2-e13-DEGs were characterized by “histone H2B,” “chaperone HSP70,” and “serpin,” while LOXL2-WT-DEGs were annotated by “zinc finger protein ZFP-36.” A specific KEGG pathway potentially correlated to carcinoma, “hsa04620:Toll-like receptor signaling pathway,” was found for the LOXL2-e13-DEGs network. As for LOXL2-WT-DEGs, four carcinoma-related KEGG pathways were found: “hsa04010:MAPK signaling pathway,” “hsa04110:Cell cycle,” “hsa04115:p53 signaling pathway,” and “hsa05200:Pathways in cancer.”

Both LOXL2-WT and LOXL2-e13 were involved in regulation of gene expression, but they might play different roles or at different levels. It showed LOXL2-e13-DEGs might affect chromosome remodeling, such as “GO:0006334~nucleosome assembly,” “GO:0031497~chromatin assembly,” “GO:0006333~chromatin assembly or disassembly,” “GO:0006323~DNA packaging,” “GO:0000786~nucleosome,” and “GO:0065004~protein-DNA complex assembly” from GO annotation. However, LOXL2-WT-DEGs might regulate expression mainly at a transcriptional level, as indicated by “GO:0045449~regulation of transcription,” “GO:0006350~transcription,” “GO:0006355~regulation of transcription, DNA-dependent,” “GO:0010629~negative regulation of gene expression,” “GO:0016604~nuclear body,” and “GO:0000775~chromosome, centromeric region.”

To our great interest, these two charts also provided clues revealing the different roles of LOXL2-WT and LOXL-e13 in malignancy. The GO “biological process” annotation showed LOXL2-WT-DEGs were involved in cell motility (GO:0048870), cell migration (GO:0016477), and regulation of cell adhesion (GO:0030155). Though LOXL-e13-DEG was not annotated by such terms related to the malignant degree of tumors, it could not be excluded that LOXL2-e13 also plays important roles in carcinoma development. Significantly, LOXL2-e13-DEGs could regulate the activities of several enzymes or modify proteins, including enzyme inhibitor activity (GO:0004857), Rab-protein geranylgeranyltransferase complex (GO:0005968), endopeptidase inhibitor activity (GO:0004866), glutaminase activity (GO:0004359), peptidase inhibitor activity (GO:0030414), serine-type endopeptidase inhibitor activity (GO:0004867), protein serine/threonine kinase activity (GO:0004674), and GO:0016310~phosphorylation.

3.4. Prioritization of e13-WT-DEGs

To provide a deep view of the different influence on ESCC cells between LOXL2-e13 and its wild type, we selected the DEGs between LOXL2-e13 and its wild-type counterpart for future analyses. We obtained 275 e13-WT-DEGs, including 118 downregulated genes and 157 upregulated genes. To obtain a full view of the influence of the e13-WT-DEGs, its PPI subnetwork was also constructed and contained 2741 nodes and 43507 edges, including 100 upregulated genes and 54 downregulated genes (Figure 4(a)). These results indicated that overexpression of LOXL2-e13 had a widely different impact on the mRNA expression profile compared to wild-type LOXL2, providing important clues to reveal the specific functions and potential molecular mechanisms of LOXL2-e13.

(a)

(b)

(c)

(d)

Figure 4: Prioritization analyses of e13-WT-DEGs based on the PPI subnetwork. (a) PPI subnetwork of e13-WT-DEGs. (b) Random Walk with Restart algorithm was used to score all proteins in the PPI network with LOXL2 set as the seed node. The size of each node in the PPI subnetwork was designed as a gradient based on the scores. (c) The e13-WT-DEGs were extracted from (b) to illustrate their sizes. (d) The DEGs were rearranged according to their closeness to LOXL2-e13 protein.

We applied the RWR algorithm to identify how the e13-WT-DEGs would be ranked by their closeness in the context of LOXL2-e13 overexpression. The RWR algorithm gave each protein member in the e13-WT-DEGs PPI subnetwork a probability score, which ranged from −2.56 to −8.74 (the more negative the score, the less significant the protein) (Figure 4(b)). Subsequently, only the nodes of e13-WT-DEGs were extracted from the modified PPI subnetwork for a better illustration (Figure 4(c)). Since some of the e13-WT-DEGs scored closely, we rearranged e13-WT-DEGs after their scores were log10-transformed and classified by the score ranges. For example, only the seed node LOXL2 was classed as the A layer; DEGs with a log-transformed score of −2.0~−2.99 were classified as the B layer, and DEGs with a log-transformed score of −3.0~−3.99 were classified as the C layer (Figure 4(d)). We found upregulated ERO1L and ITGA3 and downregulated PLAA ranked in the first class close to LOXL2-e13-DEGs and upregulated MAPK8 ranked in the second class. These results provided the priorities of e13-WT-DEGs when considering their relationships with LOXL2-e13, which could be applied as candidate genes for subsequent experimental identification.

4. Discussion

Esophageal squamous cell carcinoma (ESCC), the major histopathologic form of esophageal cancer, is one of the most prevalent cancers in Asia and is the fourth leading cause of cancer death in China [36, 37]. The functions of a splice variant might be different and vary compared to wild type [38]. Much evidence has indicated that splicing abnormalities are a hallmark of cancer [39–41]. The potential roles for splice variants in cancer might involve cell migration, cell growth, hormone responsiveness, apoptosis, and response to chemotherapy [42]. For example, two splice variants of K-Ras, K-Ras 4A, and 4B, which arise from two alternative versions of exon 4, have been found [43]. These two variants show antagonistic biological effects; K-Ras 4A exerts proapoptotic effects, while K-Ras 4B is an antiapoptotic protein. Both variants are coexpressed in many tissues, but their ratio is altered in sporadic colorectal cancer, favoring the antiapoptotic 4B isoform [44]. This indicates that the expression of tumour-specific splice variants significantly affects many cellular events, critical for cancer biology, which are still far from being illustrated.

Inspired by this evidence, we overexpressed LOXL2-e13 and determined the mRNA expression profile, in comparison with wild-type LOXL2 overexpression, to explore its specific biological roles. PPI subnetworks were generated by mapping the DEGs to a public PPI dataset to gain a full view of their influence. We show that LOXL2-WT-DEGs and LOXL2-e13-DEGs interact with thousands of other proteins, which suggests LOXL2-WT and LOXL2-e13 can greatly impact cellular activity through the cascades of interactions of the DEGs.

To better illustrate their potential molecular functions, Functional Annotation Chart analysis was applied to classify the DEGs based on related multiple gene function annotations, which allows investigators to analyze genes from many different biological aspects in a single space. Our results indicate that there are distinguishing differences between LOXL2-e13 and LOXL2 in regard to GO enrichment, protein sequence features, protein domains, and pathways. LOXL2-e13-DEGs contain a list of upregulated histone proteins, including HIST1H2AC, HIST1H2BC, HIST1H2BD, HIST2H2BE, HIST1H2BJ, HIST1H2AK, HIST2H4A, and HIST2H4B, which mediate DNA organization and play a dominant role in regulating eukaryotic transcription. On the other hand, there are as many as 78 proteins, containing zinc finger, which can interact with LOXL2-WT-DEGs, such as ZNF19, ZNF292, KLF11, CREB5, SNAI2, and SP2. It has been suggested that zinc-finger-containing proteins function in gene transcription, translation, mRNA trafficking, cytoskeleton organization, epithelial development, cell adhesion, protein folding, chromatin remodeling, and zinc sensing [45].

The other evidence for regulation of transcription by both DEG networks comes from the GO annotations, which show that LOXL2-e13-DEGs might affect chromosome remodeling, based on annotations “nucleosome assembly,” “DNA packaging,” and “protein-DNA complex assembly.” LOXL2-WT-DEGs might regulate expression mainly at the transcriptional level, such as “regulation of transcription” and “regulation of transcription, DNA-dependent.” These results suggest that both LOXL2-WT and LOXL2-e13 are involved in gene expression regulation at different levels. This could be one of the specific molecular mechanisms for LOXL-e13 compared to wild type.

The other distinguishing property for LOXL2-e13 is that its DEGs involve enzyme activities or modifications, such as “Rab-protein geranylgeranyltransferase complex,” “serine-type endopeptidase inhibitor activity,” and “peptidase inhibitor activity.” The term “serine-type endopeptidase inhibitor activity” contained 7 genes, including APP, SERPINB2, SERPINB1, SLPI, SERPINB3, SERPINB13, and OVOS2. SERPINB1 (serine protease inhibitor, clade B, member 1) has been shown to be downregulated in hepatocellular carcinoma and is correlatively related with cancer metastasis and intravasation [46]. Interference of SERPINB1 promotes migration and invasion of HCC cells, with an apparent increase in the level of active matrix metalloproteinase-2 (MMP2) [46]. It is interesting to note that SERPINB1 is also decreased 1.58-fold upon LOXL2-e13 overexpression. The “peptidase inhibitor activity” annotation contains 10 genes, including APP, CARD16, SERPINB2, BIRC6, SERPINB1, SLPI, SERPINB3, HSPA5, SERPINB13, and OVOS2. BIRC6 is upregulated 1.56-fold in the LOXL2-e13 expression profile. BIRC6 is a member of the inhibitors of apoptosis protein (IAP) family, which is thought to protect a variety of cancer cells from apoptosis [47]. BIRC6 has been found to be overexpressed in various carcinomas, such as prostate cancer, melanoma, and non-small-cell lung cancer [47–49]. The relationship between LOXL2 and metabolism was not found. We assumed that LOXL2-e13 might influence metabolism in ESCC cells and that this might be a specific biological role or molecular mechanism for LOXL2-e13.

It was important to identify which e13-WT-DEG responds to explain differing functions for LOXL2-e13. We prioritized the e13-WT-DEGs by application of the RWR algorithm based on current knowledge of PPI networks. For example, ERO1L, ITGA3, PLAA, and MAPK8 were DEGs that ranked the closest to LOXL2-e13. It was reported that the upregulation of ERO1L (ERO1-like (S. cerevisiae)) is specifically induced in hypoxic microenvironments coinciding with that of upregulated VEGF expression in human tumors, such as in hepatocellular carcinoma and glioblastoma [50]. In support of this, we find VEGF is upregulated 1.68-fold in the e13-WT-DEGs. ITGA3 (integrin, alpha 3) plays an important role in tumors. Loss of ITGA3 prevents skin tumor formation by promoting epidermal turnover and depletion of slow-cycling cells [51]. ITGA3 also serves as a biomarker for various conditions, such as a biomarker for estimation of the risks of locoregional and hematogenous dissemination of oral squamous cell carcinoma, a marker for NMuMG cells undergoing epithelial-mesenchymal transition and for cancer cells with aggressive phenotypes, and a diagnostic marker for the clinical outcome of tongue squamous cell carcinoma [52, 53]. MAPK8, also named JNK, is involved in a wide variety of cellular processes, such as proliferation, differentiation, transcription regulation, and development in both normal tissues and tumors [54]. So it is expected that LOXL2-e13 causes broad changes in mRNA expression profile, including some critical tumor-related genes, enabling LOXL2-e13 to play new and specific roles in ESCC compared to its wild-type counterpart.

5. Conclusions

In summary, we provide evidence, through generation, annotation, and comparison of PPI networks, to indicate that LOXL2-e13 might play different roles compared to wild-type LOXL2. These results provide helpful information in the experimental identification of its biological roles and explanation of its molecular mechanisms. With the development of high-throughput techniques for protein-protein interactions, the public PPI network database and its use in subsequent analyses will become more acute and reliable. Our analyses also provide a work flow to test the different roles of a splicing variant with large scale data.

Conflict of Interests

The authors have declared that no competing interests exist.

Authors’ Contribution

Bing-Li Wu and Guo-Qing Lv contributed equally to this work.

Acknowledgments

This work was supported by Grants from the National Natural Science Foundation of China (no. 81000900), the National Basic Research Program (no. 2012CB526600), the National High Technology Research and Development Program of China (no. 2012AA02A503 and no. 2012AA02A209), the Natural Science Foundation of China-GuangDong Joint Fund (no. U1301227), the Foundation for Distinguished Young Talents in Higher Education of Guangdong (LYM09081), and Shantou University Medical Research Fund. The authors thank Dr. Stanley Lin, Department of Pathophysiology, The Key Immunopathology Laboratory of Guangdong Province, Shantou University Medical College, for the assistance in revising the paper.

References

H. A. Lucero and H. M. Kagan, “Lysyl oxidase: an oxidative enzyme and effector of cell function,” Cellular and Molecular Life Sciences, vol. 63, no. 19-20, pp. 2304–2316, 2006.

Y. Kim, C. D. Boyd, and K. Csiszar, “A new gene with sequence and structural similarity to the gene encoding human lysyl oxidase,” The Journal of Biological Chemistry, vol. 270, no. 13, pp. 7176–7182, 1995.

K. Kenyon, W. S. Modi, S. Contente, and R. M. Friedman, “A novel human cDNA with a predicted protein similar to lysyl oxidase maps to chromosome 15q24-q25,” The Journal of Biological Chemistry, vol. 268, no. 25, pp. 18435–18437, 1993.

C. Jourdan-le Saux, O. le Saux, C. Gleyzal, P. Sommer, and K. Csiszar, “The mouse lysyl oxidase-like 2 gene (mLOXL2) maps to chromosome 14 and is highly expressed in skin, lung and thymus,” Matrix Biology, vol. 19, no. 2, pp. 179–183, 2000.

C. Jourdan-Le Saux, A. Tomsche, A. Ujfalusi, L. Jia, and K. Csiszar, “Central nervous system, uterus, heart, and leukocyte expression of the LOXL3 gene, encoding a novel lysyl oxidase-like protein,” Genomics, vol. 74, no. 2, pp. 211–218, 2001.

L. Asuncion, B. Fogelgren, K. S. K. Fong, S. F. T. Fong, Y. Kim, and K. Csiszar, “A novel human lysyl oxidase-like gene (LOXL4) on chromosome 10q24 has an altered scavenger receptor cysteine rich domain,” Matrix Biology, vol. 20, no. 7, pp. 487–491, 2001.

L. I. Smith-Mungo and H. M. Kagan, “Lysyl oxidase: properties, regulation and multiple functions in biology,” Matrix Biology, vol. 16, no. 7, pp. 387–398, 1998.

H. Peinado, G. Moreno-Bueno, D. Hardisson et al., “Lysyl oxidase-like 2 as a new poor prognosis marker of squamous cell carcinomas,” Cancer Research, vol. 68, no. 12, pp. 4541–4550, 2008.

D. P. Macartney-Coxson, K. A. Hood, H. J. Shi et al., “Metastatic susceptibility locus, an 8p hot-spot for tumour progression disrupted in colorectal liver metastases: 13 candidate genes examined at the DNA, mRNA and protein level,” BMC Cancer, vol. 8, article 187, 2008.

H. Offenberg, N. Brünner, F. Mansilla, F. Ørntoft Torben, and K. Birkenkamp-Demtroder, “TIMP-1 expression in human colorectal cancer is associated with TGF-B1, LOXL2, INHBA1, TNF-AIP6 and TIMP-2 transcript profiles,” Molecular Oncology, vol. 2, no. 3, pp. 233–240, 2008.

S. F. T. Fong, E. Dietzsch, K. S. K. Fong et al., “Lysyl oxidase-like 2 expression is increased in colon and esophageal tumors and associated with less differentiated colon tumors,” Genes Chromosomes and Cancer, vol. 46, no. 7, pp. 644–655, 2007.

H. J. Moon, J. Finney, and L. Xu, “MCF-7 cells expressing nuclear associated lysyl oxidase-like 2 (LOXL2) exhibit an epithelial-to-mesenchymal transition (EMT) phenotype and are highly invasive in vitro,” The Journal of Biological Chemistry, vol. 288, pp. 30000–30008, 2013.

H. E. Barker, J. Chang, T. R. Cox et al., “LOXL2-mediated matrix remodeling in metastasis and mammary gland involution,” Cancer Research, vol. 71, no. 5, pp. 1561–1572, 2011.

H. Peinado, M. D. C. Iglesias-de La Cruz, D. Olmeda et al., “A molecular role for lysyl oxidase-like 2 enzyme in Snail regulation and tumor progression,” EMBO Journal, vol. 24, no. 19, pp. 3446–3458, 2005.

G. Moreno-Bueno, F. Salvador, A. Martín et al., “Lysyl oxidase-like 2 (LOXL2), a new regulator of cell polarity required for metastatic dissemination of basal-like breast carcinomas,” EMBO Molecular Medicine, vol. 3, no. 9, pp. 528–544, 2011.

N. Herranz, N. Dave, A. Millanes-Romero et al., “Lysyl oxidase-like 2 deaminates lysine 4 in histone H3,” Molecular Cell, vol. 46, no. 3, pp. 369–376, 2012.

T. Li, L. Xu, Z. Wu et al., “Reduced nuclear and ectopic cytoplasmic expression of lysyl oxidase-like 2 is associated with lymph node metastasis and poor prognosis in esophageal squamous cell carcinoma,” Human Pathology, vol. 43, no. 7, pp. 1068–1076, 2012.

G. Q. Lv, H. Y. Zou, L. D. Liao, et al., “Identification of a novel lysyl oxidase-like 2 alternative splicing isoform, LOXL2 Δe13, in esophageal squamous cell carcinoma,” Biochemistry and Cell Biology, 2014.

M. G. Kann, “Protein interactions and disease: computational approaches to uncover the etiology of diseases,” Briefings in Bioinformatics, vol. 8, no. 5, pp. 333–346, 2007.

Z. Wu, X. Zhao, and L. Chen, “Identifying responsive functional modules from protein-protein interaction network,” Molecules and Cells, vol. 27, no. 3, pp. 271–277, 2009.

I. W. Taylor, R. Linding, D. Warde-Farley et al., “Dynamic modularity in protein interaction networks predicts breast cancer outcome,” Nature Biotechnology, vol. 27, no. 2, pp. 199–204, 2009.

T. S. Keshava Prasad, R. Goel, K. Kandasamy et al., “Human protein reference database—2009 update,” Nucleic Acids Research, vol. 37, no. 1, pp. D767–D772, 2009.

A. Chatr-Aryamontri, B. Breitkreutz, S. Heinicke et al., “The BioGRID interaction database: 2013 update,” Nucleic Acids Research, vol. 41, no. 1, pp. D816–D823, 2013.

G. C. Koh, P. Porras, B. Aranda, H. Hermjakob, and S. E. Orchard, “Analyzing protein-protein interaction networks,” Journal of Proteome Research, vol. 11, no. 4, pp. 2014–2031, 2012.

M. E. Cusick, H. Yu, A. Smolyar, et al., “Literature-curated protein interaction datasets,” Nature Methods, vol. 6, pp. 39–46, 2009.

M. S. Cline, M. Smoot, E. Cerami et al., “Integration of biological networks and gene expression data using Cytoscape,” Nature Protocols, vol. 2, no. 10, pp. 2366–2382, 2007.

Y. Assenov, F. Ramírez, S. E. Schelhorn, T. Lengauer, and M. Albrecht, “Computing topological parameters of biological networks,” Bioinformatics, vol. 24, no. 2, pp. 282–284, 2008.

B. Wu, C. Li, P. Zhang, et al., “Dissection of miRNA-miRNA interaction in esophageal squamous cell carcinoma,” PLoS One, vol. 8, Article ID e73191, 2013.

D. W. Huang, B. T. Sherman, Q. Tan et al., “DAVID Bioinformatics Resources: expanded annotation database and novel algorithms to better extract biology from large gene lists,” Nucleic Acids Research, vol. 35, no. 2, pp. W169–W175, 2007.

S. Hunter, R. Apweiler, T. K. Attwood et al., “InterPro: the integrative protein signature database,” Nucleic Acids Research, vol. 37, no. 1, pp. D211–D215, 2009.

I. Letunic, T. Doerks, and P. Bork, “SMART 7: recent updates to the protein domain annotation resource,” Nucleic Acids Research, vol. 40, no. 1, pp. D302–D305, 2012.

D. Merico, R. Isserlin, O. Stueker, A. Emili, and G. D. Bader, “Enrichment map: a network-based method for gene-set enrichment visualization and interpretation,” PLoS ONE, vol. 5, no. 11, Article ID e13984, 2010.

K. C. Chipman and A. K. Singh, “Predicting genetic interactions with random walks on biological networks,” BMC Bioinformatics, vol. 10, article 17, 2009.

K. Macropol, T. Can, and A. K. Singh, “RRW: repeated random walks on genome-scale protein networks for local cluster discovery,” BMC Bioinformatics, vol. 10, article 1471, p. 283, 2009.

A. L. Barabási and Z. N. Oltvai, “Network biology: understanding the cell's functional organization,” Nature Reviews Genetics, vol. 5, pp. 101–113, 2004.

D. M. Parkin, F. Bray, J. Ferlay, and P. Pisani, “Global cancer statistics, 2002,” CA: A Cancer Journal for Clinicians, vol. 55, no. 2, pp. 74–108, 2005.

L. Ke, “Mortality and incidence trends from esophagus cancer in selected geographic areas of China Circa 1970–90,” International Journal of Cancer, vol. 102, no. 3, pp. 271–274, 2002.

B. J. Blencowe, “Alternative Splicing: new insights from global analyses,” Cell, vol. 126, no. 1, pp. 37–47, 2006.

J. P. Venables, “Unbalanced alternative splicing and its significance in cancer,” BioEssays, vol. 28, no. 4, pp. 378–386, 2006.

R. I. Skotheim and M. Nees, “Alternative splicing in cancer: noise, functional, or systematic?” International Journal of Biochemistry and Cell Biology, vol. 39, no. 7-8, pp. 1432–1449, 2007.

M. Ladomery, “Aberrant alternative splicing is another hallmark of cancer,” International Journal of Cell Biology, vol. 2013, Article ID 463786, 6 pages, 2013.

G. S. Wang and T. A. Cooper, “Splicing in disease: disruption of the splicing code and the decoding machinery,” Nature Reviews Genetics, vol. 8, no. 10, pp. 749–761, 2007.

M. Malumbres and M. Barbacid, “RAS oncogenes: the first 30 years,” Nature Reviews Cancer, vol. 3, no. 6, pp. 459–465, 2003.

S. J. Plowman, R. L. Berry, S. A. Bader et al., “K-ras 4A and 4B are co-expressed widely in human tissues, and their ratio is altered in sporadic colorectal cancer,” Journal of Experimental and Clinical Cancer Research, vol. 25, no. 2, pp. 259–267, 2006.

J. H. Laity, B. M. Lee, and P. E. Wright, “Zinc finger proteins: new insights into structural and functional diversity,” Current Opinion in Structural Biology, vol. 11, no. 1, pp. 39–46, 2001.

X. Cui, Y. Liu, and C. Wan, “Decreased expression of SERPINB1 correlates with tumor invasion and poor prognosis in hepatocellular carcinoma,” Journal of Molecular Histology, vol. 45, pp. 59–68, 2014.

C. G. Low, I. S. U. Luk, D. Lin et al., “BIRC6 protein, an inhibitor of apoptosis: role in survival of human prostate cancer cells,” PLoS ONE, vol. 8, no. 2, Article ID e55837, 2013.

E. Tassi, M. Zanon, C. Vegetti et al., “Role of Apollon in human melanoma resistance to antitumor agents that activate the intrinsic or the extrinsic apoptosis pathways,” Clinical Cancer Research, vol. 18, no. 12, pp. 3316–3327, 2012.

X. Dong, D. Lin, C. Low et al., “Elevated expression of birc6 protein in non-small-cell lung cancers is associated with cancer recurrence and chemoresistance,” Journal of Thoracic Oncology, vol. 8, no. 2, pp. 161–170, 2013.

D. May, A. Itin, O. Gal, H. Kalinski, E. Feinstein, and E. Keshet, “Ero1-Lα plays a key role in a HIF-1-mediated pathway to improve disulfide bond formation and VEGF secretion under hypoxia: implication for cancer,” Oncogene, vol. 24, no. 6, pp. 1011–1020, 2005.

N. Sachs, P. Secades, L. Van Hulst, M. Kreft, J. Song, and A. Sonnenberg, “Loss of integrin α3 prevents skin tumor formation by promoting epidermal turnover and depletion of slow-cycling cells,” Proceedings of the National Academy of Sciences of the United States of America, vol. 109, no. 52, pp. 21468–21473, 2012.

M. Nagata, A. A. Noman, K. Suzuki, et al., “ITGA3 and ITGB4 expression biomarkers estimate the risks of locoregional and hematogenous dissemination of oral squamous cell carcinoma,” BMC Cancer, vol. 13, article 410, 2013.

T. Shirakihara, T. Kawasaki, A. Fukagawa et al., “Identification of integrin α3 as a molecular marker of cells undergoing epithelial-mesenchymal transition and of cancer cells with aggressive phenotypes,” Cancer Science, vol. 104, no. 9, pp. 1189–1197, 2013.

J. Liu and A. Lin, “Wiring the cell signaling circuitry by the NF-κB and JNK1 crosstalk and its applications in human diseases,” Oncogene, vol. 26, no. 22, pp. 3267–3278, 2007.