Skip to main content

Advertisement

A new method for mining information of co-expression network based on multi-cancers integrated data

Abstract

Background

Gene co-expression network is a favorable method to reveal the nature of disease. With the development of cancer, the way to build gene co-expression networks based on cancer data has been become a hot spot. However, there are still a limited number of current node measurement methods and node mining strategies for multi-cancers network construction.

Methods

In this paper, we introduce a new method for mining information of co-expression network based on multi-cancers integrated data, named PMN. We construct the network by combining the different types of relevant measures (linear and nonlinear rules) for different nodes based on integrated gene expression data of multi-cancers from The Cancer Genome Atlas (TCGA). For mining genes, we combine different properties (local and global characteristics) of the nodes.

Results

We uncover more suspicious abnormally expressed genes and shared pathways of different cancers. And we have also found some proven genes and pathways; of course, there are some suspicious factors and molecules that need clinical validation.

Conclusions

The results demonstrate that our method is very effective in excavating gene co-expression genes of multi-cancers.

Background

In recent years, the study of molecular biology networks has provided a powerful tool for exploring complex life activities. With the diversity and high dimension of biological sequencing data, the advantages of network pattern mining become more and more noticeable. The network can macroscopically construct the connection between biomolecules, which can enrich the amount of excavation information. Gene expression data reflect the abundance of the directly or indirectly measured gene transcription product by mRNA in cells. These data can be utilized to analyze genes that have changed their expression and the relationship between them. And gene expression data also imply how the activities of genes are affected under different conditions. Constructing co-expression networks based on gene expression data is a universal approach to molecular biology networks. There are two main types of patterns based on gene expression data: regard genes as nodes or samples as nodes. Typical representatives such as Weighted Gene Co-expression Network Analysis (WGCNA) [1] and Similarity Network Fusion (SNF) [2]. In this paper, our improvements are based on nodes for network construction.

The similarity measures based on gene expression profiles during the process of network construction are mainly obtained by calculating the linear or non-linear correlation coefficient of different gene expression profiles. Among them, universal linear correlation measurements include Pearson Correlation Coefficient (PCC), Spearman Correlation Coefficient (SCC), and Partial Correlation Coefficient which exclude other variables and so on. Measurements of non-linear correlations include, for example, Mutual Information (MI) and so on. Among them, PCC and MI are more commonly utilized. WGCNA adds an index to the PCC so that the distribution of correlation coefficient values gradually accords with the scale-free distribution. WGCNA is a method of mining module information from sequencing data. WGCNA is the current effective method with higher recognition. But it is based on the module method, lack of directional mining for the node. Butte and Kohane firstly used MI as a measure of correlation between the genes and constructed an MI genes-related network [3]. Since then, the application of MI has been widely up [4]. There are no specific standards for the association measure, and some scholars have compared above two types of measures and found that each has its own advantages and disadvantages [5].

In terms of data, the current gene network is primarily focused on single-disease dataset. With the development of next generation sequencing technology, the production of a large number of multi-omics data provides a new direction for disease research scholars to study [6, 7]. Multi-omics data can be better carried out examination research by integrating variety datasets. Integrated data include the integration of different data types of the same disease, as well as data integration of diverse diseases. These data have important biological significance for discovering the pathogenesis of diseases and providing rich and implicit information. The data from The Cancer Genome Atlas (TCGA, https://cancergenome.nih.gov/) are more interested in the changes in gene expression of human cells after birth (acquired variation). Most of the datasets provided by TCGA are integrated data, multi-group data. In addition, there are few ways to build a network based on integrated data of multi-cancers to further exploit causative agents, aiming to find the common disease-causing pathways or biotic factors in multi-cancers. In the most process of network constructions, the definition of measure between nodes is single, linear or non-linear. In the node mining aspect, the node evaluation has yet to be improved.

Based on the integrated gene expression data of the three types of cancers from TCGA, our paper provides a novel idea of constructing a co-expression network and attempts to discover common disease-causing modules of three kinds of cancers so as to provide some references for the research and development of cancers. The procedure of constructing a network draws on the combination of linear (PCC) and nonlinear (MI) similarity between genes at the same time which is called PMN, eliminating the single association of network nodes and enriching the information contained in gene networks. Considering the local and global effects in node evaluation, it is helpful for the comprehensive exploitation of nodes.

Methods

In this section, we will introduce the steps of PMN, including network construction and strategies of evaluating nodes.

The PCC between two genes X and Y is defined as follows:

$$ {\mathrm{PCC}}_{XY}=\frac{\sum_i^n\left({X}_i-\overline{X}\right)\left({Y}_i-\overline{Y}\right)}{\sqrt{\sum_{i=1}^n{\left({X}_i-\overline{X}\right)}^2{\sum}_i^n{\left({Y}_i-\overline{Y}\right)}^2}}, $$
(1)

where Xi and Yi represent the observations of the gene X and Y on the i ‐ th sample, respectively. n is the number of samples, \( \overline{X} \) and \( \overline{Y} \) are the mean observations on the gene X and Y in all samples.

The entropy of discrete random variables is defined as follows:

$$ H(X)=-\sum \limits_xp(x)\log p(x), $$
(2)

where p(x) = Pr {X = x} is the probability density function, the mutual information between two genes X and Y can be described as:

$$ MI\left(X,Y\right)=H(X)+H(Y)-H\left(X,Y\right), $$
(3)

where H(X, Y) =  − ∑xyp(x, y) log p(x, y) is joint entropy. Many studies have shown that MI is better for discrete values [8]. But gene expression data contain all continuous values. In order to make better use of the MI measure, we learn from the treatment in Minet [8] method to discretize continuous values.

The relationship between genes is not unique and the existing single measure is difficult to accurately describe the relationship between genes. We consider the different relationships between genes, so that the stronger relationships can be retained in the networks. The definition of adjacency matrix A = [aij] is based on the construction of two measures to retain the linear or non-linear relationships of genes:

$$ {a}_{ij}=\left\{\begin{array}{l}1\kern1.8em if\kern0.4em PCC>{t}_1\kern0.5em or\kern0.5em MI>{t}_2,\\ {}0\kern1.6em else,\end{array}\right. $$
(4)

where t1 and t2 are the filter thresholds of the matrices about PCC and MI. Theoretically, there is no strict requirements to formulate t1 and t2 .

For node mining, on the one hand, degree centrality of node Dx is taken into consideration: degree centrality is the most basic topological property that describes a single node in the network. The degree centrality of node v refers to the number of edges connected directly to v in the network, implying the local characteristics of a node. On the other hand, the betweenness centrality of a node is also excavable for mining nodes. The betweenness of a node is a measure of the sum of its proportions appearing in the shortest path between other nodes. The betweenness of the node is defined as follows:

$$ {B}_x={\sum}_{i\ne x\ne j\in X}\frac{\sigma_{ixj}}{\sigma_{ij}}, $$
(5)

where σij represents the counts of shortest path between node i and j, but σixj is the number of shortest path via node x. The betweenness indicates the role a node plays in connecting other nodes to each other. The higher the betweenness, the more important the node is in maintaining tight connectivity for networks, which reflecting the global characteristics of a node.

In this paper, we define a weighted value Wx for each gene node in the network by combining above two aspects to select the abnormally expressed genes:

$$ {W}_x={D}_x\times {B}_x. $$
(6)

Taking the type of values into account, multiplication can better reflect the effect of two characteristics. The measured values of the degree are greater than 0, and the values of the betweenness are less than 1. Considering the type of values and conducting extensive experimental tests, multiplication can better reflect the effect of two characteristics. It is also a combination of local and global features for nodes. Detailed schematic diagram of the overall process flow shown in Fig. 1, the following is a detailed process and details.

Fig. 1
figure1

Diagram of PMN for mining information on multi-cancers integrated data

Experiments

Datasets

The experimental data are cancer integrated gene expression data of three types of cancer from the TCGA (RNA-seq data at level 3) database: Pancreatic adenocarcinoma (PAAD), Esophageal carcinoma (ESCA) and Head and Neck squamous cell carcinoma (HNSC). There are more than 20,000 genes in each dataset, and the samples are different from each other. Among them, HNSC contains 398 tumor samples, ESCA 183 contains tumor samples, and 176 tumor samples are from PAAD. After removing the missing values, we integrated them into a matrix which contains 757 samples with 15,442 features for gene correlation calculation.

Threshold selection

Values of PCC correlation matrix are in the range of [0, 1], and correlation values of the MI matrix have no specific scope. Taking into account the overall situation, after all the correlation values of the two matrices are sorted, curve fitting is performed, and then the first inflection point is selected separately as thresholds (t1 = 0.9363 and t2 = 0.7271) of the MI and PCC, as is shown in Fig. 2 (red cross symbol is inflection point). Since there are plenty of correlation values and in order to facilitate the display, we only intercepted the previous points to exhibit.

Fig. 2
figure2

Inflection point selection

Results

After threshold filtering, we can get 242 gene nodes and 6 modules and each interaction group covers at least 18 nodes, which are visualized by Cytoscape [9], as is illustrated in Fig. 3. Due to too few genes in other modules, we cannot calculate the relative role of nodes, and we will not consider them here. In the Fig. 3, the size of a node indicates the degree of a node, and the color depth indicates the size of a messenger. Through the network of nodes in the interrelationship, we can obtain the enrichment of pathways for each module. The connections between the nodes indicate the correlation between genes, strictly dependent on PCC and MI. The discovered node interaction group in the Fig. 3 contains fewer nodes, mainly due to the threshold selection method. Threshold inflection point selection method is also a truncated threshold selection strategy; this strategy will lead to the network more serious block. And in the networks, related and unrelated nodes become straightforward, interactive modules bound not too much. There are also levels of correlations (differences in connectivity for each node) in our networks. Although in many genes, only 1% of the information about gene interactions is retained, it is from another perspective that these retainable modules must be highly relevant modules. It should be more effective to analyze the role of nodes in these highly correlated and well-structured modules.

Fig. 3
figure3

Networks of 6 Modules under PCC and MI

Finally, genes with larger weighted values implied a higher importance in the interaction group, and we list the top 10 for analysis, genes are listed in Table 1. The corresponding icon can also be found in Fig. 3. The corresponding notes are available on the GeneCards (http://www.genecards.org/).

Table 1 The top 10 genes with higher weights

According to the proven material, ASXL2 is related to cancer in several tissue types, such as breast, bladder, pancreas, ovary, prostate, and blood. This gene plays a significant role in neurodevelopment, cardiac function, adipogenesis, and osteoclastogenesis; the latest studies have shown that reproductive mutation gene ASXL2 in the maintenance of normal hematopoietic stem cell function of Acute Myeloid Leukemia (AML) has antitumor effect [10]. Abnormal expression of BUB1B has turned out to be associated with the development of lung adenocarcinoma, breast cancer, colorectal cancer and prostate cancer [11,12,13]. CLEC14A is considered as a candidate for tumor vascular targeting due to its presence at higher levels in tumor endothelium than in normal tissue endothelium. KIF11 abnormalities in prostate cancer, ovarian cancer were found [14, 15]. Some studies have shown that MYOZ1 has correlation with oral squamous cell carcinoma and it could be regarded as a molecular target for cancer applications. The over-expression of detection for TPX2 can be used as a prognostic indicator and therapeutic target for gastric cancer, prostate cancer, cervical cancer and so on [16,17,18]. Some scholars detected a decrease in CD48 expression in hepatocellular carcinoma cells [19].

The relevant pathways are found in 6 modules and we take the top 10 by p-values as is listed in the Table 2. Pathway enrichment databases involved: Kyoto Encyclopedia of Genes and Genomes (KEGG) and Reactome. Cell adhesion molecules (CAMs) play an important role in the tumorigenesis of cervical cancer [20]. The up-regulation of cyclin B1/Cdc2 plays a pivotal role in treated breast cancer cells in human in stagnation of Mitotic Prometaphase [21]. In the study of Staphylococcus aureus infection, patients with malignant tumors have a high rate of S. aureus infection [22]. The changes in the cancer cells of MHC class II antigen presentation have been noticed. Correspondingly, the development of efficient vaccines to activate the immune system against cancer has been investigated [23]. Over-expression of FOXM1 transcription factor network is linked to the development and pathogenesis of cervical cancer [24].

Table 2 Cancer related suspicious pathways

In order to further prove the advantages of PMN, we compare it with WGCNA. Firstly, we compare the number of pathways that found in the networks built by PMN and WGCNA under same scale (the number of nodes is almost, WGCNA: 230, PMN: 242), as shown in Table 3. It is clear that the network constructed by PMN which can find more pathways. And in order to understand more about these pathways, we contrast the p-values of the common pathways of the more relevant network modules obtained by the two methods. In the common pathways between the two, the difference about p-values is still quite significant, as is listed in Table 4.

Table 3 Comparison for pathways of PMN and WGCNA
Table 4 Comparison for p-values of common pathways found by PMN and WGCNA

The genes found by our method can match the pathways more accurately and with a lower error rate. WGCNA differs from PMN in that WGCNA is a method that focuses on module analysis, whose all modules together cover all genes, and focusing on the overall contribution of each module. In each module of WGCNA, the interactions of genes are very close and the ways of node mining is single. It is not conducive to distinguish the role of nodes in the single module. In principle, our method only retains some of the more relevant modules, focusing on the role of genes within the modules.

Discussion

The abnormally expressed genes have different statuses in cell networks and may have different effects on the network. Different network structures also have different tolerances for gene mutations. In this sense, the tumor is actually a kind of molecular network disease. So, it will be one-sided to ignore the whole effects and emphasize the role of single gene. PMN utilizes network structure to simulate gene interactions of cancers and discovers some abnormally expressed genes in the node relationships of networks, which have a significant effect in digging out co-expressed abnormally genes of multi-cancers.

As a kind of biomolecule network, the co-expression network has the nature of scale-free: most nodes have a small number of connections, and a small number of nodes have a large number of connections. This characteristic shows that in the dynamic process possessed by the biomolecular network, a few molecules represented by the nodes play a key role. PMN also follows this rule and takes the mining of key pathogenic nodes as the primary goal. Based on the global and local characteristics of the nodes, more suspicious abnormal genes are detected.

Among the top genes found by the PMN, ASXL2, BUB1B, KIF11 and TPX2 have been clinically proven to be associated with multi-cancers, which are a variety of genes commonly identified as mutated in cancers. The genes that abnormally and simultaneously expressed in multi-cancers are naturally considered to be key genes for the development of malignant tumors. These genes are of great significance for the prognosis of other cancers. MYO1F, SASH3 and KLHL41 are the candidate pathogenic genes that shared with multi-cancers in this paper. Until now, no research has shown that the improper expression of these genes may cause consequences for canceration, but it can still provide clues for clinical medicine. The progress of the three cancer research is not particularly impressive; we found some common anomalies of cancers that went well with the research by building networks based on the integrated dataset. To some extent, it will be beneficial for the study of three kinds of cancer.

Although the linear association rules can describe the relationship between most of the nodes, it cannot deny that the non-linear relationship exists in the gene expression data in various forms. Therefore, MI is a supplement to PCC, which makes the structure of network more complete. Degree centrality of a node is a very important symbol in the network, but simply considering the degree is one-sided. Betweenness centrality measures the effect of a node on the tightness of the entire network from a global perspective. The combination of the two can detect more critical nodes to some extent. For example, both KIF11 and CD48 are nodes that are less prominent in connectivity and betweenness, but the weighted values are prominent, indicating that the weighted is reasonable.

On the network scale, our approach only preserves a handful of the key modules of gene interaction, also gives the strategy of node mining in the module. It eliminates the influence of irrelevant information in high-dimensional data. And to a certain extent, it has some advantages for extracting key information.

Conclusions

In this paper, we introduce a method of constructing a network that combines the linear (PCC) and non-linear (MI) measurements between genes, which is based on gene expression integrated data of three kinds of cancers from TCGA. Network construction eliminates the single association between network gene nodes and enriches the information contained in the gene network. At the same time, considering the degree and betweenness of the node from the local and global perspective on the evaluation of nodes, it is conducive to mining more critical abnormally expressed node in the entire network. By obtaining genes and pathways shared by several common cancers, the effectiveness of this approach is demonstrated. Besides, the availability of additional genes and pathways may provide evidence for cancer research. Of course, the method of this article also has a lot to be improved. For example, gene expression data itself carry some noise, which will bring some interference to the latter network. Whether the weighting values of nodes can be obtained by a better combination of their characteristics, which is also the focus of our next study.

Availability of data and materials

The datasets that support the findings of this study are available in https://cancergenome.nih.gov/.

Abbreviations

ESCA:

Esophageal carcinoma

HNSC:

Head and Neck squamous cell carcinoma

MI:

Mutual Information

PAAD:

Pancreatic adenocarcinoma

PCC:

Pearson Correlation Coefficient

TCGA:

The Cancer Genome Atlas

WGCNA:

Weighted Gene Co-expression Network Analysis

References

  1. 1.

    Langfelder P, Horvath S. WGCNA: an R package for weighted correlation network analysis. BMC Bioinformatics. 2008;9:559.

  2. 2.

    Wang B, Mezlini AM, Demir F, Fiume M, Tu Z, Brudno M, Haibekains B, Goldenberg A. Similarity network fusion for aggregating data types on a genomic scale. Nat Methods. 2014;11:333–7.

  3. 3.

    Butte AJ, Kohane IS. Mutual information relevance networks: functional genomic clustering using pairwise entropy measurements. Biocomputing. 1999;2000:418–29.

  4. 4.

    Steuer R, Kurths J, Daub CO, Weise J, Selbig J. The mutual information: Detecting and evaluating dependencies between variables. Bioinformatics. 2002;18(Suppl 2):S231–40.

  5. 5.

    Song L, Langfelder P, Horvath S. Comparison of co-expression measures: mutual information, correlation, and model based indices. BMC Bioinformatics. 2012;13(1):328.

  6. 6.

    Ge SG, Xia JF, Wen S, Zheng CH. Cancer subtype discovery based on integrative model of multigenomic data. IEEE/ACM Trans Comput Biol Bioinform. 2017;14(5):1115–21.

  7. 7.

    Yang C, Ge SG, Zheng CH. ndmaSNF: cancer subtype discovery based on integrative framework assisted by network diffusion model. Oncotarget. 2017;8(51):89021–32.

  8. 8.

    Meyer PE, Lafitte F, Bontempi G. minet: a R/bioconductor package for inferring large transcriptional networks using mutual information. BMC Bioinformatics. 2008;9(1):461.

  9. 9.

    Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.

  10. 10.

    Micol JB, Pastore A, Inoue D, Duployez N, Kim E, Lee SC, Durham BH, Chung YR, Cho H, Zhang XJ. ASXL2 is essential for haematopoiesis and acts as a haploinsufficient tumour suppressor in leukemia. Nat Commun. 2017;8:15429.

  11. 11.

    Myrie KA, Percy MJ, Azim JN, Neeley CK, Petty EM. Mutation and expression analysis of human BUB1 and BUB1B in aneuploid breast cancer cell lines. Cancer Lett. 2000;152(2):193–9.

  12. 12.

    Burum-Auensen E, Deangelis PM, Schjølberg AR, Røislien J, Mjåland O, Clausen OPF. Reduced level of the spindle checkpoint protein BUB1B is associated with aneuploidy in colorectal cancers. Cell Prolif. 2008;41(4):645–59.

  13. 13.

    Fu X, Chen G, Cai Z, Wang C, Liu Z, Lin Z, Wu Y, Liang Y, Han Z, Liu J. Overexpression of BUB1B contributes to progression of prostate cancer and predicts poor outcome in patients with prostate cancer. Onco Targets Ther. 2016;9:2211–20.

  14. 14.

    Xia Q, Shan M, Ding T, Al E. Screening for differential genes of the prostate cancer and bioinformatics analysis of their interaction. China Oncol. 2017;27(3):169–76.

  15. 15.

    Xu Z, Yu Z, Cao Y, Dinh TLA, Jing W, Min Z. Identification of candidate biomarkers and analysis of prognostic values in ovarian cancer by integrated bioinformatics analysis. Med Oncol. 2016;33(11):130.

  16. 16.

    Pan HW, Su HH, Hsu CW, Huang GJ, Wu TL. Targeted TPX2 increases chromosome missegregation and suppresses tumor cell growth in human prostate cancer. Onco Targets Ther. 2017;10:3531–43.

  17. 17.

    Jiang P, Shen K, Wang X, Song H, Yue Y, Liu T. TPX2 regulates tumor growth in human cervical carcinoma cells. Mol Med Rep. 2014;9(6):2347–51.

  18. 18.

    Li B, Palma MD, Nathanson KL, Xu G, Smalley K, Herlyn M. Abstract 659: overexpression of TPX2 accelerates cell cycle progression in melanoma and is an indicator of response to Aurora a kinase inhibitors. Cancer Res. 2010;70(8 Supplement):659.

  19. 19.

    Midorikawa Y, Tsutsumi S, Taniguchi H, Ishii M, Kobune Y, Kodama T, Makuuchi M, Aburatani H. Identification of genes associated with dedifferentiation of hepatocellular carcinoma with expression profiling analysis. Jpn J Cancer Res. 2002;93(6):636–43.

  20. 20.

    Douville MA, Grognec PL. Abnormal immunoexpression of cell adhesion molecules (CAMs) in cervical cancer. Int J Surg Pathol. 2011;19(6):733–42.

  21. 21.

    Choi HJ, Fukui M, Zhu BT. Role of cyclin B1/Cdc2 up-regulation in the development of mitotic prometaphase arrest in human breast cancer cells treated with nocodazole. PLoS One. 2011;6(8):e24312.

  22. 22.

    Young LS, Perdreauremington F, Winston LG. Clinical, epidemiologic, and molecular evaluation of a clonal outbreak of methicillin-resistant Staphylococcus aureus infection. Clin Infect Dis. 2004;38(8):1075–83.

  23. 23.

    Thibodeau J, Bourgeoisdaigneault MC, Lapointe R. Targeting the MHC class II antigen presentation pathway in cancer immunotherapy. Oncoimmunology. 2012;1(6):908–16.

  24. 24.

    Chan DW, Yu S, Chiu PM, Yao KM, Liu V, Cheung A, Ngan H. Over-expression of FOXM1 transcription factor is associated with cervical cancer progression and pathogenesis. J Pathol. 2008;215(3):245–52.

Download references

Acknowledgments

Thanks for all the reviewers’ comments concerning our manuscript, those comments are all valuable and very helpful for revising and improving our work.

About this supplement

This article has been published as part of BMC Medical Genomics, Volume 12 Supplement 7, 2019: Selected articles from the 14th International Symposium on Bioinformatics Research and Applications (ISBRA-18): medical genomics. The full contents of the supplement are available at https://bmcmedgenomics.biomedcentral.com/articles/supplements/volume-12-supplement-7

Funding

This study was supported by the National Natural Science Foundation of China under grant Nos. 61872220, 61572284, 61702299, and 61701279; by the Co-Innovation Center for Information Supply & Assurance Technology, Anhui University under grant No.ADXXBZ201704.). The funding bodies had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Author information

MXH and GYL designed and implemented the method, performed the experiments, and drafted the manuscript. JXL and SSY participated in the design of the study and performed experimental examination. JS and RZ contributed to the data analysis. All authors read and approved the final manuscript.

Correspondence to Ying-Lian Gao or Jin-Xing Liu.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Hou, M., Gao, Y., Liu, J. et al. A new method for mining information of co-expression network based on multi-cancers integrated data. BMC Med Genomics 12, 155 (2019). https://doi.org/10.1186/s12920-019-0608-2

Download citation

Keywords

  • Gene co-expression network
  • TCGA
  • Abnormally expressed genes
  • Mutual information
  • Pearson correlation coefficient