BioMed Research International
Volume 2016 (2016), Article ID 4241293, 8 pages
http://dx.doi.org/10.1155/2016/4241293

Review Article

Differential Regulatory Analysis Based on Coexpression Network in Cancer Research

Junyi Li,1,2 Yi-Xue Li,1,2,3,4 and Yuan-Yuan Li2,3,4

¹Key Lab of Computational Biology, CAS-MPG Partner Institute for Computational Biology, Shanghai Institutes for Biological Sciences, Chinese Academy of Sciences, Shanghai 200031, China
²Shanghai Center for Bioinformation Technology, 1278 Keyuan Road, Shanghai 201203, China
³Shanghai Industrial Technology Institute, 1278 Keyuan Road, Shanghai 201203, China
⁴Shanghai Engineering Research Center of Pharmaceutical Translation, 1278 Keyuan Road, Shanghai 201203, China

Received 14 April 2016; Revised 9 June 2016; Accepted 12 June 2016

Academic Editor: Zhenguo Zhang

Copyright © 2016 Junyi Li 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.

With rapid development of high-throughput techniques and accumulation of big transcriptomic data, plenty of computational methods and algorithms such as differential analysis and network analysis have been proposed to explore genome-wide gene expression characteristics. These efforts are aiming to transform underlying genomic information into valuable knowledges in biological and medical research fields. Recently, tremendous integrative research methods are dedicated to interpret the development and progress of neoplastic diseases, whereas differential regulatory analysis (DRA) based on gene coexpression network (GCN) increasingly plays a robust complement to regular differential expression analysis in revealing regulatory functions of cancer related genes such as evading growth suppressors and resisting cell death. Differential regulatory analysis based on GCN is prospective and shows its essential role in discovering the system properties of carcinogenesis features. Here we briefly review the paradigm of differential regulatory analysis based on GCN. We also focus on the applications of differential regulatory analysis based on GCN in cancer research and point out that DRA is necessary and extraordinary to reveal underlying molecular mechanism in large-scale carcinogenesis studies.

1. Introduction

In the past decade, plenty of computational methods and algorithms such as differential analysis and network analysis [1, 2] are proposed to explore genome-wide gene expression characteristics with rapid development of high-throughput technologies and accumulation of big transcriptomic data. These efforts in computational genomic area are dedicated to transform underlying genomic information into valuable knowledges in biological and medical research fields [3, 4]. Recently, tremendous integrative research aims to interpret the development and progress of cancers because elucidating molecular regulatory mechanisms, especially the dysregulation mechanisms, of neoplastic diseases makes great benefit in medical and pharmaceutical aspects. Although partial different regulatory functions of cancer hallmarks such as evading growth suppressors and resisting cell death [5] have been revealed, the whole dysregulation mechanisms are far from clear.

Cancer is a complex disease and an effective way to study regulatory role of genes involved in cancer is to summarize them into network [6]. It is suggested that genes having similar or correlated expression patterns might contribute to the same regulatory function and gene coexpression patterns revealed by coexpression network analysis may lead to more insightful discovery on the underlying regulatory mechanisms [2, 7]. By comparing the difference of the regulatory networks between cancer and normal status, specific differential network of genes can be identified as dysfunctional in cancer. A large number of reverse engineering approaches have been developed to construct regulatory network from gene expression data. For examples, Xiao suggested Boolean model to analyze and stimulate the gene regulatory network [8]. Some methods based on Bayesian model lead to Bayesian networks and they are widely applied [9–11]. Nonlinear differential equation model is also developed to construct the regulatory network [12]. Prior biological knowledge such as transcription factor- (TF-) target regulatory relationships or miRNA-target regulatory relationships can also be integrated into modelling framework [11, 13, 14]. These reverse and forward integrated approaches are supposed to have smaller false positive rate to extract informative insights of transcriptomic behaviors.

Although network analysis provides the possibility to comprehensively understand biological processes, it does increase the computational complexity. Decreasing the searching space before network analysis is necessary in high dimension data analysis. An obvious strategy of reducing the computational burden is to build a subnetwork around a given set of genes such as previously reported disease-related genes [15] or around differentially expressed genes [16–18]. Differential expression analysis (DEA) compares the mean expression value of genes between case and control samples and identifies significantly differentially expressed genes by statistical tests. In current transcriptomic analysis procedure, DEA has become the basic and the very first analysis step.

Recently, differential coexpression analysis (DCEA) increasingly plays a robust complement to DEA [2] and is widely used in discovering the system properties of carcinogenesis features. By calculating the change of correlations between gene pairs instead of mean expression level, DCEA provides more information about phenotypic change-related regulatory network [19–24]. Therefore, differential regulatory analysis based on coexpression network may detect more insights into regulatory mechanisms.

In this review, we will introduce the paradigm of differential regulatory analysis (DRA) based on gene coexpression network (GCN). We also focus on the applications of DRA based on GCN in cancer research and point out that DRA is necessary and extraordinary to reveal underlying regulatory mechanism in large-scale carcinogenesis studies.

2. Paradigm of Differential Regulatory Analysis Based on Gene Coexpression Network

Differential regulatory analysis based on gene coexpression network has been widely used in carcinogenesis regulation research and basically includes three procedures as shown in Figure 1: constructing gene coexpression network based on transcriptomic data, regulatory analysis according to gene coexpression network, and differential regulatory comparison between different conditions.

Figure 1: Paradigm of differential regulatory analysis based on gene coexpression network. The paradigm of differential regulatory analysis based on gene coexpression network includes but is not limited to three procedures. (a) Constructing gene coexpression network based on genomic transcriptomic data. (b) Regulatory analysis according to gene coexpression network. (c) Differential regulatory comparison between different conditions.

2.1. Construction of Gene Coexpression Network

In a gene coexpression network, genes are nodes and their correlations are represented by the edges of network. Pearson correlation coefficient (PCC) is the mostly used score to measure the tendency of gene expression correlation [25–28]. The value of PCC ranges from −1 to 1 and higher absolute value of PCC means higher correlation between gene pairs. When constructing gene coexpression network, a correlation threshold is selected. After removing the nonsignificant edges or negligible coordinated gene pairs by the threshold, gene coexpression network is constructed by the significantly correlated gene pairs remained as shown in Figure 1(a) [29].

Weighted correlation network analysis (WGCNA) [27, 28] is widely used for constructing coexpression network based on gene expression data and implementing network analysis [19, 21, 30]. It summarizes clusters of highly correlated genes by defining a continuous network adjacency which is a power of initial one to reduce the low-adjacency gene pairs. WGCNA analyzes the cluster structure and explores the relationships between modules or that between modules and genes. GCN topological characters can be well studied by WGCNA and this great advantage makes WGCNA one of most used GCN construction methods in research. There are many coexpression module detection methods provided by WGCNA for users to choose their own preferred one. Meanwhile, choosing a suitable threshold is required for GCN construction in WGCNA. One limitation of WGCNA is that its GCN construction is undirected. Other prior knowledge is needed if further regulatory analysis based on GCN is designed.

Link-based quantitative methods in DCGL [26, 31] employ a half-thresholding strategy to construct specific GCNs. That is, if at least one of the two coexpression values of a specific link exceeds the threshold, the link in both coexpression networks from two different conditions is kept [26, 31]. In this way, minute variations are ignored by filtering out those noninformative links whose correlation values in both networks are insignificant. GCNs constructed by DCGL are also without directions.

Gaussian graphical model (GGM) is another approach to construct gene coexpression network [32–34]. Based on the assumption that the covariance of gene pairs follows a multivariate Gaussian distribution, partial correlation between gene pairs is calculated as the degree of correlation after the effects of other genes are removed. Unlike the fact that correlations in PCC-based method are calculated by gene pairs themselves, correlations of gene pair in GGM-based methods take into account information of other genes, which makes GGM-based GCN more similar to real biological network.

Some algorithms are proposed with focus on how to infer the structure of gene correlation relationships. Algorithm for the Reconstruction of Accurate Cellular Networks (ARACNE) [35, 36] is a method of GCN construction, which also pays attention to partial network properties of GCN by counting gene triplets. The Context Likelihood of Relatedness (CLR) [37] algorithm calculates the relative correlation based on the empirical correlations over surrounding genes. MRNET is an iterative feature selection algorithm and uses a maximum relevance and minimum redundancy criterion [38].

DECODE (differential coexpression and differential expression) combines the information from both differential coexpression and differential expression to set up the thresholds systematically based on a chi-square maximization [39].

A recent GCN construction algorithm is proposed by Planar Filtered Network Analysis (PFNA) and Multiscale Embedded Gene Coexpression Network Analysis (MEGENA) [40, 41]. According to this algorithm, GCN is constructed on a spherical surface so that links between gene pairs do not cross the others. They have advantages in extracting most relevant information from similarity matrix from complex network system based on topological sphere.

2.2. Regulatory Analysis according to Gene Coexpression Network

After the gene coexpression network is constructed based on the transcriptomic data, regulatory information can be extracted by various regulatory analysis methods from the GCN according to research desires as shown in Figure 1(b). The most common way is to use the prior knowledge of TF-target regulatory relationships or miRNA-target regulatory relationships to highlight the specific regulatory subnetwork [25].

Another important method looking for regulatory elements is clustering method. A major goal of coexpression analysis is discovering biologically related modules or gene groups. In WGCNA, hierarchical clustering method is used to identify highly correlated gene subnetworks [27]. Significantly compact subclusters are verified by the average shortest path distance within each cluster over the cluster size. Nonnegative matrix factorization clustering method [42–44] is also employed to clustering coexpressed genes with features of interests.

According to topological structure of gene coexpression network, some interesting characters such as hub genes can help to explain regulatory function contained in the GCN. Connectivity in GCN presents how a gene connects to other genes and hub genes are the ones with very high level of connectivity. Hub genes are normally connections between different gene modules and should be a specific research focus for investigations into cancer-correlated gene modules. For instance, Yang et al. build gene coexpression networks based on transcriptomic and clinical data of four cancer types and discovered that prognostic mRNA genes tended not to be hub genes [24]. They suggested that hubs genes coordinate genes over different pathways to participate in the regulatory processes. Chou et al. investigated endometrial cancers (ECs) hub genes by constructing WGCNA coexpression network and these hub genes are involved in antigen processing, cell adhesion, and cell-cycle regulation [21]. On the other hand, loss of connectivity in coexpression network is a common topological trait among the different kinds of cancer [45].

2.3. Differential Regulatory Comparison between Different Conditions

Distinguishing different regulatory elements between different conditions such as tumor and normal tissue, or different cancer types, or even cross-species (human and mouse) [46] help to understand the dysfunctional regulation.

There are two ways to perform differential regulatory comparison between different conditions. The first way is to construct gene coexpression network based on each condition and compare the difference between constructed GCNs to extract different regulation elements as shown in Figure 1(c) [24, 30]. The other way is to calculate the significantly different correlation between various conditions and build a network based on these selected gene pairs [47, 48]. Differential regulatory comparison between different conditions is able to find the differential genes or gene modules across different conditions, providing useful information as well. In DCGL v2, differentially coexpressed TFs are defined as differential regulated genes (DRGs), and DRGs are ranked for prioritizing regulators that are putatively causative to the phenotype of interests in DRrank function in the R package of DCGL2 [25]. For example, RIF algorithm in DRrank function combines three types of transcriptomic information and assigns a high score to those TFs that are “cumulatively most differentially wired to the abundant most differentially expressed genes” [25, 49].

3. Applications of Differential Regulatory Analysis Based on GCN in Cancer Research

Recently, differential regulatory analysis based on GCN is applied in more and more cancer studies. In the following, we give some examples of its applications and summarize advantages of this integrative method.

3.1. Revealing Dysfunctional Regulatory Genes and Subnetworks in Cancer Research

The direct advantage of differential regulatory analysis (DRA) is that DRA is able to distinguish dysfunctional regulatory subnetworks or pathways in cancer status. For example, Jiang et al. constructed highly preserved gene ontology biological process (GO_BP) gene coexpression network and prostate cancer coexpression network by using WGCNA approach. With regulatory analysis they discovered 548 GO_BP coexpression modules and 294 prostate cancer coexpression modules. By comparing the difference of these modules, they identified 55 conserved prostate cancer coexpression modules [30]. And there are five modules which are significantly enriched with prostate cancer candidate genes. These five modules are featured with regulation of apoptosis, response to stress, cellular localization, and protein localization [30]. Udyavar et al. performed coexpression network construction based on a dataset of combined normal, adenocarcinoma, squamous cell carcinoma, and small-cell lung cancer (SCLC) tissue specimens by WGCNA. They compared the distribution of significant modules across four types of samples and derived an SCLC-specific hub network classifier and identified spleen tyrosine kinase as candidate biomarker and therapeutic target for SCLC [19].

Cancer is considered as a complex disease with multilevel progressing process. DRA based on GCN is able to bring light to dynamic regulatory relationships of cancer in its different progress levels. For instance, Cao et al. first constructed gene coexpression networks of normal, adenoma, and carcinoma-specific gastric carcinogenesis to decrease the searching space for potential regulatory genes [48]. After these potential regulatory genes are acquired, three differential networks are constructed. By comparing constructed differential networking information and signaling pathway information of three developing stages, the regulation roles of GATA6 and ESRRG and their signaling pathways in gastric carcinogenesis were suggested [48]. This work frame is prospective and extendable to other cancer researches. Wu et al. also performed a system-level study of gastric cancer by constructing five phenotype-specific coexpression networks [50]. Their comparison analysis of connectivity reveals that hub genes which only exit in the normal networks play important roles in gastric tumorigenesis and hub genes only related to tumor networks are enriched in specific biological terms. Ruan et al. identified specific pathways associated with renal cell carcinoma (RCC) based on differentially coexpressed links which were detected by three methods: Pearson’s correlation, Bayesian network, and WGCNA. These RCC-related pathways help to explain underlying regulatory mechanisms of RCC [29]. Khosravi et al. built independent gene regulatory networks from each prostate cancer and found critical transcription factors involved in prostate cancer based on hub type variation [51]. These dynamic network studies across cancer stages well reveal the change of regulatory patterns during the progression of cancers.

3.2. Practical and Beneficial in Multilevel Network Analysis When Integrated with miRNAs and lncRNAs Data

With the emerging roles of microRNAs (miRNAs) and long noncoding RNAs (lncRNAs) in gene regulatory networks, more and more genomic studies on miRNAs and lncRNAs are performed in cancer research aspects [52–58]. By integrating genomic miRNAs and lncRNAs with mRNAs data to construct multi-level co-expression network and analyze differential regulatory mechanism, the understanding of carcinogenesis functions of miRNAs and lncRNAs can be greatly improved from the views of system level.

For example, Lin et al. constructed a cross-cancer miRNA differential coexpression network and identified two potential miRNA-regulated oncomodules associated with poor survival outcomes in patients [59]. This study suggested that disruption of miRNA positive coexpression in cancer might contribute to cancer development. There are many efforts made to discover the regulatory action of lncRNAs in cancer [56, 57]. DRA based on coexpression network is also extendable to lncRNAs-cancer gene network analysis. InCaNet is a data resource which contains precalculated significant coexpression pairs of 9641 lncRNAs and 2544 well-classified cancer genes in 2922 matched TCGA samples [60]. And InCaNet helps to explore regulatory functions of particular lncRNA-cancer gene interaction in cancer studies. Most lncRNAs’ regulatory functions are unknown and DRA based on lncRNAs-cancer gene network has the exact ability to perform the predictions. According to the assumption that genes or nodes in a subgroup may execute similar functions, Cogill and Wang identified a list of previously uncharacterized lncRNAs coexpressed with key cancer genes in their study [20] and Hao et al. inferred lncRNAs related to esophageal squamous cell carcinoma ESCC from constructed coexpression network and differential regulatory analysis [61]. Moreover, an integrated miRNA-mRNA-lncRNA coexpression network analysis was performed by Wu et al. to study the oestrogen receptor-regulated transcriptome in breast cancer [62]. All these multilevel studies expand our understanding of regulatory mechanism in cancer biology.

3.3. Applicable in Medical and Pharmaceutical Aspects

Since DRA based on GCN is a system-level analysis method and explores the regulatory mechanisms of diseases, it has been widely applied in medical and pharmaceutical aspects. By integrating other pieces of information such as clinical information or drug-target genes information, DRA based on GCN has more potential to contribute to theoretical base of medical and pharmaceutical researches.

For example, prognostic genes are very important for cancer prognosis and treatment. By integrating survival information, Yang et al. studied the system-level prognostic genes across four cancer types by DRA based on GCN [24]. Discovering new biomarkers or molecular subtypes of cancer is also valuable for stratification in clinical studies. DRA-based signature has the ability to classify patients into different subtypes with different clinical results. Meanwhile, DRA-based signature is featured with different regulatory patterns of each subtype. Wu et al. revealed a novel three-transcription-factor signature including AHR, NFIL3, and ZNF423 for glioma molecular subtypes by DRA based on GCN. This three-gene DRA-based signature clusters glioma patients into three major subtypes which are significantly different in patient survival as well as transcriptomic patterns [47]. Jin et al. captured a 12-gene network module of ovarian cancer by constructing weighted survival and differential coexpression network and this module shows a close correlation with cell death [63]. All these prognostic DRA studies based on GCN help to provide a more accurate survival prediction. Moreover, DRA-based prognostic signature has more potential to explore carcinogenesis mechanisms which lead to a better precision medicine in cancer diagnosis and treatment.

4. Discussion

Regulatory analysis is always a focal point in biological research. Understanding the function of each regulatory element in biological process is fundamental and challenging. Large-scale and multilevel sequencing data provide more opportunities to reveal molecular regulatory mechanism from the systematic viewpoint. Differential regulatory analysis is designed for distinguishing the differential regulatory elements in different conditions or dysfunctional regulation specific for an abnormal condition. For example, cancer is considered as a complex genetic disease and different phonotypes in cancer embody regulatory level mechanisms. Dysfunctional regulatory elements take priority in carcinogenesis studies because of their important roles in regulation as well as their potential in cancer treatments.

In the early past decade, gene expression data for specific phenotype is very limited and researchers had to use various cell lines data to construct conceptual gene regulatory networks [14]. Then it is very difficult to explain the differential regulatory relationships between tumor types. Recently, with the accumulation of large-scale and multiscale data, researchers are able to apply differential regulatory analysis to identify specific regulatory patterns in a given cancer type. Since differential regulatory analysis based on coexpression network has a system-level property, it has great strength to discover underlying molecular mechanism and dysfunctional regulatory elements from large-scale data of complex system. With the development of differential regulatory analysis based on coexpression network itself and its applications in more and more genomic and big data research, it presents its essential and prospective role in cancer research.

Technically, construction of coexpression networks as the basic starting point of differential regulatory analysis plays important and critical role in the whole investigation process. Therefore, both the data sets and the way to construct networks must be carefully examined. For example, Ballouz et al. suggest minimal experimental criteria to obtain useful functional connectivity and topology information of coexpression with microarrays greater than 20 samples and read depth greater than 10 M per sample [64]. Selection of differential regulatory genes relies on the method or algorithm chosen in DRA based on GCN [22, 65, 66]. For instance, differential analysis between two conditions, which is followed by regulatory analysis, can be performed after two gene coexpression networks are constructed. Sometimes, regulatory analysis based on coexpression networks between two conditions is conducted first and comparison of differential regulatory elements is performed later. All these steps in current differential regulatory analysis methods are relatively flexible depending on the aim of specific research and information available. The DRA based on GCN method might need further standardization and refinement to better serve more carcinogenesis research in medical and pharmaceutical fields.

5. Conclusion

In this review, we summarize the paradigm of differential regulatory analysis based on coexpression network of transcriptomic data and the applications of differential regulatory analysis based on GCN in cancer research. Differential regulatory analysis based on GCN is demonstrated as a necessary and potential tool to reveal underlying molecular mechanism in basic functional genomic research as well as practical carcinogenesis studies.

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

This work is supported by the grants from the National Key Scientific Instrument and Equipment Development Project of China (2012YQ03026108), National Natural Science Foundation of China (31501077), the National “973” Key Basic Research Development Program (2012CB316501 and 2013CB910801), and the Fundamental Research Program of Shanghai Municipal Commission of Science and Technology (14DZ1951300 and 14DZ2252000).

References

E. E. Schadt, “Molecular networks as sensors and drivers of common human diseases,” Nature, vol. 461, no. 7261, pp. 218–223, 2009.

A. de la Fuente, “From ‘differential expression’ to ‘differential networking’—identification of dysfunctional regulatory networks in diseases,” Trends in Genetics, vol. 26, no. 7, pp. 326–333, 2010.

E. E. Schadt, S. H. Friend, and D. A. Shaywitz, “A network view of disease and compound screening,” Nature Reviews Drug Discovery, vol. 8, no. 4, pp. 286–295, 2009.

E. E. Schadt, B. Zhang, and J. Zhu, “Advances in systems biology are enhancing our understanding of disease and moving us closer to novel disease treatments,” Genetica, vol. 136, no. 2, pp. 259–269, 2009.

D. Hanahan and R. A. Weinberg, “Hallmarks of cancer: the next generation,” Cell, vol. 144, no. 5, pp. 646–674, 2011.

P. K. Kreeger and D. A. Lauffenburger, “Cancer systems biology: a network modeling perspective,” Carcinogenesis, vol. 31, no. 1, pp. 2–8, 2009.

M. R. J. Carlson, B. Zhang, Z. Fang, P. S. Mischel, S. Horvath, and S. F. Nelson, “Gene connectivity, function, and sequence conservation: predictions from modular yeast co-expression networks,” BMC Genomics, vol. 7, article 40, 2006.

Y. Xiao, “A tutorial on analysis and simulation of boolean gene regulatory network models,” Current Genomics, vol. 10, no. 7, pp. 511–525, 2009.

D. Husmeier, “Sensitivity and specificity of inferring genetic regulatory interactions from microarray experiments with dynamic Bayesian networks,” Bioinformatics, vol. 19, no. 17, pp. 2271–2282, 2003.

Z. S. H. Chan, L. Collins, and N. Kasabov, “Bayesian learning of sparse gene regulatory networks,” BioSystems, vol. 87, no. 2-3, pp. 299–306, 2007.

S. Gao and X. Wang, “Quantitative utilization of prior biological knowledge in the Bayesian network modeling of gene expression data,” BMC Bioinformatics, vol. 12, article 359, 2011.

T. T. Vu and J. Vohradsky, “Nonlinear differential equation model for quantification of transcriptional regulation applied to microarray data of Saccharomyces cerevisiae,” Nucleic Acids Research, vol. 35, no. 1, pp. 279–287, 2007.

K. Tu, H. Yu, Y.-J. Hua et al., “Combinatorial network of primary and secondary microRNA-driven regulatory mechanisms,” Nucleic Acids Research, vol. 37, no. 18, pp. 5969–5980, 2009.

H. Yu, K. Tu, Y.-J. Wang et al., “Combinatorial network of transcriptional regulation and microRNA regulation in human cancer,” BMC Systems Biology, vol. 6, article 61, 2012.

B. Zhang, H. Li, R. B. Riggins et al., “Differential dependency network analysis to identify condition-specific topological changes in biological networks,” Bioinformatics, vol. 25, no. 4, pp. 526–532, 2009.

D. Nitsch, L.-C. Tranchevent, B. Thienpont et al., “Network analysis of differential expression for the identification of disease-causing genes,” PLoS ONE, vol. 4, no. 5, Article ID e5526, 2009.

R. Gill, S. Datta, and S. Datta, “A statistical framework for differential network analysis from microarray data,” BMC Bioinformatics, vol. 11, article 95, 2010.

D. Nitsch, J. P. Gonçalves, F. Ojeda, B. de Moor, and Y. Moreau, “Candidate gene prioritization by network analysis of differential expression using machine learning approaches,” BMC Bioinformatics, vol. 11, article 460, 2010.

A. R. Udyavar, M. D. Hoeksema, J. E. Clark et al., “Co-expression network analysis identifies Spleen Tyrosine Kinase (SYK) as a candidate oncogenic driver in a subset of small-cell lung cancer,” BMC Systems Biology, vol. 7, supplement 5, article S1, 2013.

S. B. Cogill and L. Wang, “Co-expression network analysis of human lncRNAs and cancer genes,” Cancer Informatics, vol. 13, supplement 5, pp. 49–59, 2014.

W.-C. Chou, A.-L. Cheng, M. Brotto, and C.-Y. Chuang, “Visual gene-network analysis reveals the cancer gene co-expression in human endometrial cancer,” BMC Genomics, vol. 15, article 300, 2014.

M. M. Bourdakou, E. I. Athanasiadis, and G. M. Spyrou, “Discovering gene re-ranking efficiency and conserved gene-gene relationships derived from gene co-expression network analysis on breast cancer data,” Scientific Reports, vol. 6, Article ID 20518, 2016.

S. Hong, H. Dong, L. Jin, and M. Xiong, “Gene co-expression network and functional module analysis of ovarian cancer,” International Journal of Computational Biology and Drug Design, vol. 4, no. 2, pp. 147–164, 2011.

Y. Yang, L. Han, Y. Yuan, J. Li, N. Hei, and H. Liang, “Gene co-expression network analysis reveals common system-level properties of prognostic genes across cancer types,” Nature Communications, vol. 5, article 3231, 2014.

J. Yang, H. Yu, B.-H. Liu et al., “DCGL v2.0: an R package for unveiling differential regulation from differential co-expression,” PLoS ONE, vol. 8, no. 11, Article ID e79729, 2013.

B.-H. Liu, H. Yu, K. Tu, C. Li, Y.-X. Li, and Y.-Y. Li, “DCGL: an R package for identifying differentially coexpressed genes and links from gene expression microarray data,” Bioinformatics, vol. 26, no. 20, Article ID btq471, pp. 2637–2638, 2010.

P. Langfelder and S. Horvath, “WGCNA: an R package for weighted correlation network analysis,” BMC Bioinformatics, vol. 9, article 559, 2008.

B. Zhang and S. Horvath, “A general framework for weighted gene co-expression network analysis,” Statistical Applications in Genetics and Molecular Biology, vol. 4, no. 1, article 17, Article ID 20518, 2005.

X. Ruan, H. Li, B. O. Liu et al., “A novel method to identify pathways associated with renal cell carcinoma based on a gene co-expression network,” Oncology Reports, vol. 34, no. 2, pp. 567–576, 2015.

J. Jiang, P. Jia, Z. Zhao, and B. Shen, “Key regulators in prostate cancer identified by co-expression module analysis,” BMC Genomics, vol. 15, no. 1, article 1015, 2014.

H. Yu, B.-H. Liu, Z.-Q. Ye, C. Li, Y.-X. Li, and Y.-Y. Li, “Link-based quantitative methods to identify differentially coexpressed genes and gene Pairs,” BMC Bioinformatics, vol. 12, article 315, 2011.

T. Wang, Z. Ren, Y. Ding et al., “FastGGM: an efficient algorithm for the inference of gaussian graphical model in biological networks,” PLOS Computational Biology, vol. 12, no. 2, Article ID e1004755, 2016.

P. Ingkasuwan, S. Netrphan, S. Prasitwattanaseree et al., “Inferring transcriptional gene regulation network of starch metabolism in Arabidopsis thaliana leaves using graphical Gaussian model,” BMC Systems Biology, vol. 6, article 100, 2012.

S. Aburatani, F. Sun, S. Saito, M. Honda, S.-I. Kaneko, and K. Horimoto, “Gene systems network inferred from expression profiles in hepatocellular carcinogenesis by graphical gaussian model,” Eurasip Journal on Bioinformatics & Systems Biology, vol. 2007, Article ID 47214, 2007.

A. A. Margolin, I. Nemenman, K. Basso et al., “ARACNE: an algorithm for the reconstruction of gene regulatory networks in a mammalian cellular context,” BMC Bioinformatics, vol. 7, no. 1, article S7, 2006.

R. A. Chávez Montes, G. Coello, K. L. González-Aguilera, N. Marsch-Martínez, S. de Folter, and E. R. Alvarez-Buylla, “ARACNe-based inference, using curated microarray data, of Arabidopsis thaliana root transcriptional regulatory networks,” BMC Plant Biology, vol. 14, no. 1, article 97, 2014.

J. J. Faith, B. Hayete, J. T. Thaden et al., “Large-scale mapping and validation of Escherichia coli transcriptional regulation from a compendium of expression profiles,” PLoS biology, vol. 5, no. 1, article e8, 2007.

P. E. Meyer, K. Kontos, F. Lafitte, and G. Bontempi, “Information-theoretic inference of large transcriptional regulatory networks,” EURASIP Journal on Bioinformatics and Systems Biology, vol. 2007, Article ID 79879, 2007.

T. W. H. Lui, N. B. Y. Tsui, L. W. C. Chan, C. S. C. Wong, P. M. F. Siu, and B. Y. M. Yung, “DECODE: an integrated differential co-expression and differential expression analysis of gene expression data,” BMC Bioinformatics, vol. 16, no. 1, article 182, 2015.

W.-M. Song and B. Zhang, “Multiscale embedded gene co-expression network analysis,” PLoS Computational Biology, vol. 11, no. 11, Article ID e1004574, 2015.

M. Tumminello, T. Aste, T. Di Matteo, and R. N. Mantegna, “A tool for filtering information in complex systems,” Proceedings of the National Academy of Sciences of the United States of America, vol. 102, no. 30, pp. 10421–10426, 2005.

D. D. Lee and H. S. Seung, “Learning the parts of objects by non-negative matrix factorization,” Nature, vol. 401, no. 6755, pp. 788–791, 1999.

J. J.-Y. Wang, X. Wang, and X. Gao, “Non-negative matrix factorization by maximizing correntropy for cancer clustering,” BMC Bioinformatics, vol. 14, article 107, 2013.

K. Inamura, T. Fujiwara, Y. Hoshida et al., “Two subclasses of lung squamous cell carcinoma with different gene expression profiles and prognosis identified by hierarchical clustering and non-negative matrix factorization,” Oncogene, vol. 24, no. 47, pp. 7105–7113, 2005.

R. Anglani, T. M. Creanza, V. C. Liuzzi et al., “Loss of connectivity in cancer co-expression networks,” PLoS ONE, vol. 9, no. 1, Article ID e87075, 2014.

A. Aytes, A. Mitrofanova, C. Lefebvre et al., “Cross-species regulatory network analysis identifies a synergistic interaction between FOXM1 and CENPF that drives prostate cancer malignancy,” Cancer Cell, vol. 25, no. 5, pp. 638–651, 2014.

S. Wu, J. Li, M. Cao, Y.-X. Li, and Y.-Y. Li, “A novel integrated gene coexpression analysis approach reveals a prognostic three-transcription-factor signature for glioma molecular subtypes,” in Proceedings of the International Conference on Intelligent Biology and Medicine ((ICIBM '15), Indianapolis, Ind, USA, November 2015.

M. S. Cao, B. Y. Liu, W. T. Dai, W. X. Zhou, Y. X. Li, and Y. Y. Li, “Differential network analysis reveals dysfunctional regulatory networks in gastric carcinogenesis,” American Journal of Cancer Research, vol. 5, no. 9, pp. 2605–2625, 2015.

N. J. Hudson, A. Reverter, and B. P. Dalrymple, “A differential wiring analysis of expression data correctly identifies the gene containing the causal mutation,” PLoS Computational Biology, vol. 5, no. 5, Article ID e1000382, 2009.

J. Wu, X. Zhao, Z. Lin, and Z. Shao, “A system level analysis of gastric cancer across tumor stages with RNA-seq data,” Molecular BioSystems, vol. 11, no. 7, pp. 1925–1932, 2015.

P. Khosravi, V. H. Gazestani, M. Akbarzadeh, S. Mirkhalaf, M. Sadeghi, and B. Goliaei, “Comparative analysis of prostate cancer gene regulatory networks via hub type variation,” Avicenna Journal of Medical Biotechnology, vol. 7, no. 1, pp. 8–15, 2015.

C.-C. Lin, Y.-J. Chen, C.-Y. Chen, Y.-J. Oyang, H.-F. Juan, and H.-C. Huang, “Crosstalk between transcription factors and microRNAs in human protein interaction network,” BMC Systems Biology, vol. 6, article 18, 2012.

J. Xu, C.-X. Li, Y.-S. Li et al., “MiRNA-miRNA synergistic network: construction via co-regulating functional modules and disease miRNA topological features,” Nucleic Acids Research, vol. 39, no. 3, pp. 825–836, 2011.

G. A. Calin and C. M. Croce, “MicroRNA signatures in human cancers,” Nature Reviews Cancer, vol. 6, no. 11, pp. 857–866, 2006.

V. N. Kim and J.-W. Nam, “Genomics of microRNA,” Trends in Genetics, vol. 22, no. 3, pp. 165–173, 2006.

S. W. Cheetham, F. Gruhl, J. S. Mattick, and M. E. Dinger, “Long noncoding RNAs and the genetics of cancer,” British Journal of Cancer, vol. 108, no. 12, pp. 2419–2425, 2013.

S. Wang and E. J. Tran, “Unexpected functions of lncRNAs in gene regulation,” Communicative and Integrative Biology, vol. 6, no. 6, Article ID e27610, 2013.

A. Quitadamo, L. Tian, B. Hall, and X. Shi, “An integrated network of microRNA and gene expression in ovarian cancer,” BMC Bioinformatics, vol. 16, supplement 5, article S5, 2015.

C.-C. Lin, R. Mitra, F. Cheng, and Z. Zhao, “A cross-cancer differential co-expression network reveals microRNA-regulated oncogenic functional modules,” Molecular BioSystems, vol. 11, no. 12, pp. 3244–3252, 2015.

Y. Liu and M. Zhao, “lnCaNet: pan-cancer co-expression network for human lncRNA and cancer genes,” Bioinformatics, vol. 32, no. 10, pp. 1595–1597, 2016.

Y. Hao, W. Wu, F. Shi et al., “Prediction of long noncoding RNA functions with co-expression network in esophageal squamous cell carcinoma,” BMC Cancer, vol. 15, article 168, 2015.

Q. Wu, L. Guo, F. Jiang, L. Li, Z. Li, and F. Chen, “Analysis of the miRNA-mRNA-lncRNA networks in ER+ and ER- breast cancer cell lines,” Journal of Cellular and Molecular Medicine, vol. 19, no. 12, pp. 2874–2887, 2015.

N. Jin, H. Wu, Z. Miao et al., “Network-based survival-associated module biomarker and its crosstalk with cell death genes in ovarian cancer,” Scientific Reports, vol. 5, Article ID 11566, 2015.

S. Ballouz, W. Verleyen, and J. Gillis, “Guidance for RNA-seq co-expression network construction and analysis: safety in numbers,” Bioinformatics, vol. 31, no. 13, pp. 2123–2130, 2015.

O. Odibat and C. K. Reddy, “Ranking differential hubs in gene co-expression networks,” Journal of Bioinformatics and Computational Biology, vol. 10, no. 1, Article ID 1240002, 2012.

H. Yu, R. Mitra, J. Yang, Y. Li, and Z. Zhao, “Algorithms for network-based identification of differential regulators from transcriptome data: a systematic evaluation,” Science China Life Sciences, vol. 57, no. 11, pp. 1090–1102, 2014.