Skip to main content

Comparisons of gene coexpression network modules in breast cancer and ovarian cancer

Abstract

Background

Breast cancer and ovarian cancer are hormone driven and are known to have some predisposition genes in common such as the two well known cancer genes BRCA1 and BRCA2. The objective of this study is to compare the coexpression network modules of both cancers, so as to infer the potential cancer-related modules.

Methods

We applied the eigen-decomposition to the matrix that integrates the gene coexpression networks of both breast cancer and ovarian cancer. With hierarchical clustering of the related eigenvectors, we obtained the network modules of both cancers simultaneously. Enrichment analysis on Gene Ontology (GO), KEGG pathway, Disease Ontology (DO), and Gene Set Enrichment Analysis (GSEA) in the identified modules was performed.

Results

We identified 43 modules that are enriched by at least one of the four types of enrichments. 31, 25, and 18 modules are enriched by GO terms, KEGG pathways, and DO terms, respectively. The structure of 29 modules in both cancers is significantly different with p-values less than 0.05, of which 25 modules have larger densities in ovarian cancer. One module was found to be significantly enriched by the terms related to breast cancer from GO, KEGG and DO enrichment. One module was found to be significantly enriched by ovarian cancer related terms.

Conclusion

Breast cancer and ovarian cancer share some common properties on the module level. Integration of both cancers helps identifying the potential cancer associated modules.

Background

Despite of decades of intensive study and substantial progress in understanding breast cancer and ovarian cancer, these two diseases remain the most prevalent malignancy in women, and the important causes of death in women. Among solid gynaecological tumors, breast cancer is the most often diagnosed tumor while ovarian cancer is the most deadly gynaecological neoplasia. Both breast and ovarian cancer are hormone driven and are known to have some predisposition genes in common. The major genes associated with susceptibility to breast and ovarian cancer are the two well-known high-penetrance cancer genes: BRCA1 and BRCA2 [13]. However, mutations in these genes account for only a very small percent of all breast cancers and ovarian cancers. Other genes such as TP53, PTEN, and STK11/LKB1, are even less common causes of breast and ovarian cancer [4]. Despite tremendous efforts to conquer such malignant diseases, research on studying the mechanism of cancer development and developing effective preventive measures is still a hot topic.

The high speed development of high-throughput technologies such as next generation sequencing of the human genome, gene expression microarrays, identification of the changes of copy numbers has dramatically accelerated the study aiming at predicting and curing such diseases. Many works have been published to address the topics on associated susceptibilities, potential biomarkers, cancer predictions and so on [519]. Several of them put breast cancer and ovarian cancer together in their studies [510]. These works either borrowed information from each other with the assumption that both cancers have similar etiologies [5, 7, 9], or conducted research on the differences between the pathogenic mechanisms [6, 8]. Some review papers also analyzed the related research progresses of both cancers together [4, 10]. In this work, we also put both cancers together to study the complex gene coexpression patterns with the network tools.

Network has been widely applied to study the complex interactions between genes, proteins, and other small molecules. It is also a popular tool for studying breast cancers and ovarian cancers [1116, 2022]. When using networks to study the complex interactions, one typical concept is module, which is the densely connected subnetworks. With the network module analysis, we may infer the susceptibility genes in cancer [16], identify the biomarkers [14], and predict the prognosis of the cancer patients [11]. Current studies on cancer-related network modules are mainly based on one network, which is also true for breast cancer [14] and ovarian cancer [16]. Even in the study of more than one network, the module identification process is one network by one network, and then comparisons between different networks are performed [16]. Recently, several module identification methods on multiple networks are proposed [2224]. Among them, the method proposed in [22] introduced an algorithm to find the differential modules in different networks. It mainly concentrates on the differential part. While in the paper [24], the method not only can find the modules in each network, but also can align the modules at the same time. Thus both the common and the differential parts can be detected. In the following, we compared the modules that were identified from the gene coexpression networks of breast cancer and ovarian cancer using the method in [24]. We analyzed the basic properties of the modules including density, average degree, distribution difference etc., and we did enrichment analysis of Gene Ontology (GO), KEGG pathway, Disease Ontology (DO) and Gene Set Enrichment Analysis (GSEA). By comparing the modules, both the common properties and the differences between the two cancers are detected.

Methods

Data sets

The level 3 gene expression data for breast cancer (BRCA) and ovarian cancer (OV) were downloaded from The Cancer Genome Atlas (TCGA). The gene expression data were generated with UNC AgilentG4502A. We chose the samples from the solid tissues only. There are 526 and 572 samples for both BRCA and OV, respectively. The expression value of 17,814 genes was measured. The missing data for each specific gene was imputed with the mean value of the known samples. We also downloaded the most updated protein-protein interaction (PPI) data for humans from BioGrid https://thebiogrid.org . We chose the genes having PPIs in the gene expression data, and 9603 genes were selected.

To choose the differential expressed genes for BRCA and OV, t-test and Kolmogorov-Smirnov test were applied. The genes having p-values less 0.01 in both tests were chosen, where the p-values were adjusted by controlling the false discovery rate. We then got 7742 genes.

Coexpression network construction

We first computed the pairwise Pearson correlation coefficient r ij to measure the coexpression levels between gene i and gene j. To make the correlations across the two networks comparable, we ‘normalized’ the Pearson correlation coefficients with the same method as shown in [25]. Fisher’s z transformation score for each r ij was first computed as \(z_{ij}=0.5\log \frac {1+r_{ij}}{1-r_{ij}}\). The z-scores were then normalized to make z ij follow normal distribution. After this, the ‘normalized’ correlations were obtained by transforming back to r ij . With these steps, the correlations in the two networks will be on the same level. We then did hard thresholding to make the network be unweighted. We took 99.5% quantile of the transformed correlation coefficients as the threshold. If the absolute value of the correlation coefficient is greater than the threshold, we assigned an edge between the corresponding genes, otherwise, there is no edge. With this method, the average degree is about 37 in both networks.

Module identification in BRCA and OV coexpression networks

Before we did module identification, we first removed the genes having no links with any other genes in both networks. The left gene number became 6779. Consider the constructed gene coexpression network G1 (BRCA), G2 (OV) consisting of 6779 genes. We let the adjacency matrices for both networks be A1,A2, where A k (i,j)=1 represents there is an edge between gene i and gene j. We first applied the model proposed in [24] to cluster the genes in lower subspace for both networks. This model aims at finding the clusters in multiple networks and aligning the clusters at the same time. The main idea is to use spectral clustering to find the cluster in each network, and align them by maximizing the cluster similarity of multiple networks. Here, we only have two networks. We assume the putative number of clusters M in both networks is given first.

Let Sk be the assignment of the 6779 vertices into M clusters for the network G k , where \({S}_{im}^{k}=1\), if iG k belongs to the m-th cluster, otherwise \({S}_{im}^{k}=0\) for i=1,2,,6779,m=1,2,,M,k=1,2. The optimization model is formulated as:

$$\begin{array}{@{}rcl@{}} \min &&\,\, \sum\limits_{k=1,2}\sum\limits_{m=1}^{M}\frac{{S_{,m}^{k}}^{T}(D_{k}-A_{k}){S}_{,m}^{k}}{{S_{,m}^{k}}^{T}{S}_{,m}^{k}}-\beta\sum\limits_{k,l=1,2}\sum\limits_{m=1}^{M}\frac{{S_{,m}^{k}}^{T}S_{,m}^{l}}{\|S_{,m}^{k}\|_{2}\|S_{,m}^{l}\|_{2}}\\ \text{s.t.} &&\,\, S_{i,m}^{k}\in \{0,1\},i=1,2,\cdots,6779; m=1,2,\cdots,M; k=1,2;\\ &&\,\,\sum\limits_{m=1}^{M}{S}_{,m}^{k}=\mathbf{1},\text{for}\,\,k=1,2. \end{array} $$
(1)

The first term in the objective function is to do clustering in both networks separately, and the second term defines the similarity of the clusters in both networks measured by cosine function. β balances the contributions from inter-network and intra-network.

To solve the optimization problem (1), the variable \(S^{k}_{im}\) is relaxed. Using the same technique as spectral clustering, the optimization problem is transformed to:

$$\begin{array}{@{}rcl@{}} \min \,\, {\tilde\Psi}(\tilde S)=\text{Tr}(\tilde S^{T}C\tilde S),\,\text{s.t.} \,\, \tilde S^{T}\tilde S=2I_{M}, \end{array} $$

where \( C=\left (\begin {array}{cc} L_{1}&{}\mathbf {0}\\ \mathbf {0}&L_{2}\\ \end {array} \right) -\beta \left (\begin {array}{cccc} \mathbf {0}&I_{n}\\ I_{n}&0\\ \end {array} \right), \tilde S = \left (\begin {array}{c} \tilde S^{1}\\ \tilde S^{2}\\ \end {array} \right), \)\(\tilde S_{,m}^{k}=\frac {S_{,m}^{k}}{\|S_{,m}^{k}\|_{2}}\), and L k =D k A k , 0 is an n×n matrix with all entries being zero.

By computing the eigenvectors corresponding to the M smallest eigenvalues of matrix C, the original vertices in the networks are projected to a space of dimension M. To get the clusters, we may use k-means clustering to cluster the data points similar to spectral clustering. Due to the large size of the network, k-means does not work well. Instead, we applied hierarchical clustering with complete linkage to cluster the vertices. The distance is chosen to be the spearman distance. This is because when the size of matrix C is large, the range of the eigenvector entries is large, but their order is comparatively stable. The algorithm is summarized in ‘Algorithm’.

Selection of parameter β and M The parameter β controls the connections between the vertices in both networks. When β=0, it is equivalent to finding the clusters in two networks separately. When β becomes larger, the corresponding vertices in both networks tend to belong to the same cluster. We note that even when a group of vertices are densely connected in the first network, while their corresponding parts are isolated in the second network, the method will put all the isolated vertices in the same cluster as in the first network. Here, since both networks were controlled to have a close number of total connections, we directly set β=1, which means the connection weight between two networks is the same as that within both networks. We note that when β>1 and it is within a reasonable range, the results do not change much.

The number of clusters M was chosen according to the eigenvalues of matrix C. M corresponds to the first big eigengap [26]. We note that here M is not the number of clusters in either of the two networks because by choosing β, the isolated vertices in one network can also be clustered together depending on the other network. It should be the number of the union clusters in both networks. Thus the method can find both the consistent clusters and the differential clusters having quite different connection probabilities.

After the above clustering procedures, we can get a cluster label for each gene. The corresponding clusters in both networks may include different genes. We take the union of the genes as the cluster. Due to the large number of genes, the clustering method may have some bias. Some unconnected subnetworks may be clustered together. Before going to further analysis, we need to check each cluster such that it cannot include unconnected subnetworks. These resulted subnetworks are defined as the modules we identified.

Enrichment analysis

To see the associations between the identified gene modules and gene functions, we did Gene Ontology (GO)[27], KEGG pathway [28], and Disease Ontology (DO) [29] enrichment analysis for each module. GO annotates genes to biological processes (BP), molecular functions (MF), and cellular components (CC) in a directed acyclic graph structure. We only considered BP terms here. KEGG annotates genes to pathways, and Disease Ontology (DO) annotates genes with human disease associations. To see whether the identified modules show statistically significant, concordant differences between the two diseases, we also did Gene Set Enrichment Analysis (GSEA) [30]. GSEA was also done for the three types of enrichment analysis: GO, KEGG, and DO. We implemented the enrichment analysis with ‘clusterProfiler’ [31]. For all the cases, we let the cutoff be the Benjamini-Hochberg adjusted p-value 0.05, and recorded all the enriched terms with p-value less than 0.05.

Results

We did module identification in the gene coexpression networks of both breast cancer and ovarian cancer simultaneously. Figure 1 shows the first 350 eigenvalues of the matrix C. We chose M to be 215. After we clustered the genes into 215 clusters with hierarchical clustering, we removed those with size less than 5 and greater than 800. With the method described, finally, we got 62 modules.

Fig. 1
figure 1

Eigenvalues of matrix C

To look at the module structures in both networks, we first computed the average degree (\({\bar d}_{BRCA}\), \({\bar d}_{OV}\)) and density(D BRCA ,D OV ) for each module. Besides, we did statistical test for each module to see whether the modules in the two networks have the same connection distribution. We assume the connection between any two vertices is randomly generated following Bernoulli distribution with a given probability. We applied t-test to see whether the probability in the two networks is the same. The p-values are recorded. For all the 62 modules, after we did all the enrichment analysis, we removed those modules that have no enriched terms. Finally, 43 modules were found to have at least one type of enrichment. We put the module size, the average degree, the density, the t-test p-value, and the number of enriched terms for GO, KEGG pathway, DO, GSEA in Table 1. For GSEA, we only found GSEA GO enrichment terms.

Table 1 Statistics of the identified modules

From the t-test p-values of the modules in both coexpression networks, the structures of 29 modules are significantly different with p-values less than 0.05, of which 25 modules have larger densities in OV. Thirty one modules are enriched by GO terms, 25 modules are enriched by KEGG pathways, 18 modules are enriched by DO terms, and one module is enriched by GSEA GO terms. One module (module 14) is enriched by all the four terms. Nine modules (module 7,15,17,19,21,29,33,34,37) are enriched by GO, KEGG, and DO. We checked the details of each module and the enriched terms. All the enrichment results are put in the Additional files 1, 2, 3, 4 and 5.

GO enrichment analysis

We listed the modules that have an enriched p-value less than 10−5 in Table 2. Three modules including module 6, 36, and 39 have the same connection distribution in both cancers. These modules have small sizes, and very significant enrichments. Six among the 7 genes in module 6 are involved in the regionalization and pattern specification process. Module 36 is mainly related to glutathione metabolic process and xenobiotic metabolic process. Five of the 6 genes are involved in these processes. Module 39 is mainly related to skeletal system development. Five among 6 genes are involved in this process. In other three modules 5, 14, and 27, the connections in the OV coexpression network are much denser. There are several isolated genes in the BRCA coexpresson network. These three modules are involved in many complex biological processes significantly. One typical example is the enriched term ‘GO:0016259 selenocysteine metabolic process’. There are 89 genes in the background 16655 genes involved in this process. In module 5, 44 among the 195 (overlapping with the background) genes are involved in this process. For the term ‘GO:0006614 SRP-dependent cotranslational protein targeting to membrane’, 45 genes among the 195 genes are involved in the process compared to the background 108 genes of 16655 genes in this process. These show the high correlations among the genes involved in the same process. Compared to OV, the genes involving in the same biological processes in BRCA have much less correlations, which leads to the sparser module structures.

Table 2 Enriched GO terms with p-value <10−5 in all the modules

KEGG enrichment analysis

The KEGG pathways that enrich the modules having a p-value less than 10−4 are listed in Table 3. Module 5 is enriched by the pathway ‘hsa03010 Ribosome’. Forty-six out of 130 (overlapping with the background) genes in this module belong to this pathway compared to 154 out of 7274 in the background genes. From the GO enriched terms of this module, it is clear that this module is mainly involved in the translation process. The pathways related to module 14 are mainly related to diseases. ‘IL-17 signaling pathway’ [32], ‘TNF signaling pathway’ [33], ‘Rheumatoid arthritis’ [34], and ‘MAPK signaling pathway’ [35] were shown to have relations with BRCA. For these enriched terms, no existing literatures addressed their associations with OV to the best of our knowledge. Module 23 is also related to cancers. Eight genes in this module belong to ‘Pathways in cancer’. It is also enriched by ‘breast cancer’ with a p-value 0.049, with 3 genes associated with breast cancer in this module. This module has no enrichments related to OV. The pathway ‘hsa05418 Chemical carcinogenesis’ that enriches module 36 is also related to cancers. Chemical carcinogens may contribute significantly to the causation of a sizable fraction, perhaps a majority, of human cancers [36].

Table 3 Enriched KEGG terms with p-value <10−4 in all the modules

DO enrichment analysis

Table 4 lists the DO enriched terms in all the modules with p-value less than 0.01, and one term that is related to OV. Module 14 is enriched by BRCA, with 7 among 20 (overlapping with the background) genes in this module associated with BRCA. Module 38 is enriched by OV with 4 among 11 (overlapping with the background) genes in this module being associated with OV. There are no OV related enriched terms in module 14, and no BRCA related terms in module 38. This is mainly due to the different known genes associated with the two cancers. We note that module 34 is also enriched by the female organ cancer.

Table 4 Enriched DO terms with p-value <0.01 in all the modules

From all the above analysis, we found that most modules are the general modules that may not be associated with BRCA and OV. The validated enriched terms related to BRCA are much more than that of OV. One reason may be there are more researches conducted on BRCA. Among all the enriched modules, module 14 is enriched by BRCA, and module 38 is enriched by OV. Module 36 is not enriched by these two diseases, but the enriched terms are related to cancer treatment. Among these three modules, the structure of module 14 is shown to be different with a p-value 1.76E-11 in these two cancers. In the following, we give some details of these three modules.

Enrichment analysis for module 14

Module 14 consists of 25 genes. The module structure in both networks is significantly different with a p-value 1.76E-11. Figure 2 shows the module structures. It is much denser in OV compared to that in BRCA. Figure 3 shows the dotplot of GO, KEGG, DO enrichment results. For GO and DO, we plotted the first 30 enriched terms with the smallest p-values. We plotted all the 27 enriched terms for KEGG. Figure 4 shows the associations between the genes and the enriched terms. We selected the most enriched 12 GO terms with p-value less than 10−4, 10 KEGG terms with p-value less than 0.01, and 15 DO terms with p-value less than 0.01.

Fig. 2
figure 2

Network structure of module 14. a BRCA; b OV

Fig. 3
figure 3

Enrichment results for module 14. a GO, 30 terms with minimal p-values; b KEGG, all enriched 27 terms; c DO, 30 terms with minimal p-values

Fig. 4
figure 4

Enrichment results for module 14. a GO terms with p-value <10−4; b KEGG terms with p-value <0.01; c DO terms with p-value <0.01

The enriched GO terms are mainly related to different responses, such at inflammatory response, response to molecule of bacterial origin, responses to wounding, immune response, etc.. Several responses have been shown to be associated with cancers [37, 38]. For example, the inflammation as the seventh hallmark of cancer plays important roles in cancer development. Inflammatory cells may facilitate angiogenesis and promote the growth, invasion, and metastasis of tumor cells, which may change the genetic instability in cancer cells. Controlling the regulation of inflammatory response has a potential in both prevention and treatment of cancer [37]. There are 9 genes included in this module associated with the enriched term ‘regulation of inflammatory response’, which achieves a p-value of 5.43E-08. This module also includes 7 genes that are associated with the enriched term ‘adaptive immune response’, which achieves a p-value of 8.78E-05. The immune response can result in the proliferation of antigen-specific lymphocytes. When antibodies and T-cell receptors are expressed and up-regulated, immunity is acquired. Then the immune systems will initiate antigenic responses against carcinomas. A new approach to the treatment of cancer is immunotherapy, which aims to up-regulate the immune system in order that it may better control carcinogenesis [38]. Figure 4a shows the relations between the 12 GO terms with p-value less than 10−4 and the 25 genes in this module. The genes ‘CXCL1’, ‘CCL20’, ‘ANXA1’, ‘S100A9’, ‘S100A12’, ‘TNF’, ‘DUSP10’, ‘TNFAIP2’ are included in more than 6 enriched terms. There are 10 genes associated with the most enriched terms ‘response to lipopolysaccharide’ and ‘response to molecule of bacterial origin’. For these two terms, we have not found the related literature that addresses their relations to cancers.

In the KEGG enriched terms, ‘TNF signaling pathway’ [33],‘Rheumatoid arthritis’ [34], ‘MAPK signaling pathway’ [35], and ‘IL-17 signaling pathway’ [32] have shown to be associated with BRCA. TNF is a major inflammatory cytokine shown to be highly expressed in breast carcinomas. It induces a wide range of intracellular signal pathways including apoptosis and cell survival as well as inflammation and immunity [33]. ‘MARK signaling pathway’ is involved in various cellular functions, including cell proliferation, differentiation and migration. Research on signaling pathway switch in breast cancer shows that in a large proportion of breast cancer, MARK signaling pathway is repressed, while another important pathway is activated. This mechanism may have impacts on the balance between self-renewal, proliferation, and differentiation of the tumor-initiating cells [35]. IL-17 plays crucial roles in both acute and chronic inflammatory responses. It is shown to have a direct association with breast cancer invasion in human breast tumors. IL-17 directly induced breast cancer cell invasion. There should be a potential mechanism for breast cancer invasion and tumor progression [32]. ‘Rheumatoid arthritis’ is mainly related to immune systems. Research shows that the risk of breast cancer is increased in non-Caucasians patients with rheumatoid arthritis while it decreased in Caucasian population [34]. Figure 4b shows the associations between the enriched KEGG terms and the genes. Four genes including ‘CXCL1’, ‘TNF’, ‘FOS’, and ‘IL1A’ connect to at least 6 of the 10 terms. There are at least 6 genes associated with ‘TNF signaling pathway’, ‘MAPK signaling pathway’, and ‘IL-17 signaling pathway’. For these enriched pathways, we have not found their associations with OV.

In the DO enriched terms, ‘breast carcinoma’ reaches the p-value 0.001. Seven genes including ‘S100A9’, ‘TNF’, ‘SOC53’, ‘CD55’, ‘IL1A’, ‘ANXA1’, ‘GADD45A’ among the 25 genes in this module are associated with BRCA. Several other diseases also enrich module 14, including arteriosclerosis disease, nutrition disease, etc.. However, OV is not on this list. Figure 4c shows the associations between the genes and the diseases. ‘S100A9’, ‘TNF’, ‘SOCS3’, and ‘IL1A’ connect to at least 10 diseases among the 15 enriched diseases.

In the GSEA study, we ordered the genes according to the t-test p-value between the two diseases and did the analysis. Finally, 17 GO terms enrich this module. Table 5 shows the enriched terms. The 8 sequential genes having the largest t-test p-values are all in the enriched biological processes. They are ‘TNFAIP3’, ‘S100A9’, ‘BCL3’, ‘MAFF’, ‘TMEM173’, ‘JUNB’, ‘CEBPD’, ‘NFKBIZ’. Figure 5 shows the patterns for the running enrichment score. All the enriched terms have a similar pattern for these 8 genes.

Fig. 5
figure 5

GSEA GO enrichment for module 14

Table 5 Enriched GO terms with GSEA for module 14

By comparison of this module structure between BRCA and OV, we found module 14 is closely related to BRCA from the above enrichment analysis. From the gene coexpression network of BRCA (Fig. 2), we see some breast cancer associated genes are isolated, such as ‘TNF’, ‘CD55’, ‘IL1A’. In OV network, these genes are highly correlated and clustered into one module. This may show that due to the tumor, the correlations of some cancer related genes decrease in BRCA.

Enrichment analysis for module 36

This module includes 6 densely connected genes in both diseases, which has a p-value 0.46 in t-test for the structure difference. Of the 6 genes, 5 genes GSTM1, GSTM2, GSTM3, GSTM4, GSTM5 encode the glutathione S-transferase that belongs to the mu class. These genes function in the detoxification of electrophilic compounds such as carcinogens, therapeutic drugs, environmental toxins and products of oxidative stress, by conjugation with glutathione. Thus this module is enriched by the related biological processes such as: ‘glutathione derivative metabolic process’, ‘xenobiotic metabolic process’, ‘sulfur compound biosynthetic process’, etc.. It is also enriched by the related pathways such as ‘glutathione metabolism’, ‘drug metabolism-cytochrome P450’, ‘chemical carcinogenesis’, and so on. Another gene is BCAR3, which is associated with estrogen resistance and breast cancer. It is translated to the breast cancer anti-estrogen resistance protein 3. Although this module is not enriched by BRCA and OV in DO significantly, it is related to the treatment of breast cancer [39].

Enrichment analysis for module 38

Module 38 includes 11 genes. Figure 6 shows the module structure in both cancers. The connection probability in these two networks is statistically the same with t-test, although the detailed connections are different. This module is enriched by 36 GO terms with adjusted p-value less than 0.05. The gene number involved in the related biological processes is at most 4. It is enriched by 3 diseases including ovarian cancer. Four of the 14 genes are associated with OV using DO enrichment. One typical gene is ‘WT1’, which connects to several other genes in the OV network, while it has no connections in the BRCA network. This gene is necessary for the development of the ovaries in females, and thus is associated with ovarian cancer. The ‘WT1’ protein has been found to bind a host of cellular factors such as p53. It has been ranked as the No.1 target for cancer immunotherapy by the National Cancer Institute. However, it has no associations with BRCA to the best of our knowledge. A densely connected subnetwork in this module is Kallikrein-related peptidases (KLK5, KLK6, KLK7, KLK8, KLK10). This gene family can be taken as the novel cancer biomarkers as shown in [40]. The potential of KLKs as diagnostic, prognostic, and treatment monitoring biomarkers for many types of malignancies has been extensively investigated including breast cancer and ovarian cancer. Overall, this module is closely related to cancers including breast cancer and ovarian cancer, and is more associated with ovarian cancer.

Fig. 6
figure 6

Network structure of module 38. a BRCA; b OV

Discussion and conclusion

Breast cancer and ovarian cancer are important causes of death for women. Both of them are harmone driven and are known to have some common susceptible genes such as BRCA1 and BRCA2. Several published works have studied both cancers together with the assumption that they have the same etiologies. Coexpression network modules for both breast cancer and ovarian cancer have been studied separately by several researchers. However, there are no comparisions between the coexpression network modules between breast cancer and ovarian cancer. By comparing the modules in both cancers, we aim at finding more relations between both cancers including both the similarities and the differences.

In this work, we compared the coexpression network modules of both cancers by simultaneously identifying the modules of both cancers. By projecting the genes into a lower space representation, we did hierarchical clustering of all the genes with complete linkage of spearman distance to get the modules. 43 modules were identified to be enriched by at least one of GO, KEGG pathway, and DO terms. In most of these modules, density in OV network is larger than that in BRCA network. There are 31, 25 and 18 modules enriched by GO, KEGG pathway, and DO terms, respectively. By checking the details of each enriched term, we found that one module (module 14) is enriched by GO, KEGG and DO terms related to breast cancer. Among all GO and KEGG enriched terms, several have been validated. And there are 7 genes associated with breast cancer in this module. However, there is no enrichment information related to ovarian cancer. Another module 38, which is enriched by ovarian cancer in DO, is closely related to several cancers including breast cancer and ovarian cancer, and is more associated with ovarian cancer. From the analysis of the genes in module 36, we found that it is related to the treatment of breast cancer, although it is not enriched by breast cancer terms in DO.

Comparison of the identified modules shows the differences of breast cancer and ovarian cancer on the module level. Different modules are enriched by the two cancers significantly. It also shows some common properties of both cancers such as the KLKs family in module 36. More importantly, by simultaneous clustering of both cancers, the potential cancer related modules can be identified. For example, module 14 is significantly enriched by breast cancer related terms, but much connection information is borrowed from the ovarian cancer network. This may imply the associations between this module and breast cancer.

References

  1. Miki Y, Swensen J, Shattuckeidens D, Futreal PA, Harshman K, Tavtigian SV, Liu Q, Cochran C, Bennett L, Ding W, et al. A strong candidate for the breast and ovarian cancer susceptibility gene brca1. Science. 1994; 266(5182):66–71.

    Article  CAS  PubMed  Google Scholar 

  2. Wooster R, Neuhausen SL, Mangion J, Quirk Y, Ford D, Collins N, Nguyen K, Seal S, Tran T, Averill D, et al. Localization of a breast cancer susceptibility gene, brca2, to chromosome 13q12-13. Science. 1994; 265(5181):2088–90.

    Article  CAS  PubMed  Google Scholar 

  3. Claus EB, Schildkraut JM, Thompson WD, Risch N. The genetic attributable risk of breast and ovarian cancer. Cancer. 1996; 77(11):2318–24.

    Article  CAS  PubMed  Google Scholar 

  4. Wooster R, Weber BL. Breast and ovarian cancer. N Engl J Med. 2009; 348(23):2339–47.

    Article  Google Scholar 

  5. M’́a’́att’́a KM, Nurminen R, Kankuri-Tammilehto M, Kallioniemi A, Laasanen S-L, Johanna S. Germline emsy sequence alterations in hereditary breast cancer and ovarian cancer families. BMC Cancer. 2017. 17(496).

  6. Pearce CL, Hirschhorn JN, Wu AH, Burtt NP, Stram DO, Young S, Kolonel LN, Henderson BE, Altshuler D, Pike MC. Clarifying the progins allele association in ovarian and breast cancer risk: A haplotype-based analysis. J Natl Cancer Inst. 2005; 97(1):51–9.

    Article  CAS  PubMed  Google Scholar 

  7. Song H, Ramus SJ, Kjaer SK, Dicioccio RA, Chenevixtrench G, Pearce CL, Hogdall E, Whittemore AS, Mcguire V, Hogdall C, et al. Association between invasive ovarian cancer susceptibility and 11 best candidate snps from breast cancer genome-wide association study. Hum Mol Genet. 2009; 18(12):2297–304.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Eckstein N. Platinum resistance in breast and ovarian cancer cell lines. J Exp Clin Cancer Res. 2011; 30(1):91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  9. Thompson ER, Dragovic RL, Stephenson S, Eccles D, Campbell IG, Dobrovic A. A novel duplication polymorphism in the fanca promoter and its association with breast and ovarian cancer. BMC Cancer. 2005; 5(1):43.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Newman B, Millikan RC, King M. Genetic epidemiology of breast and ovarian cancers. Epidemiol Rev. 1997; 19(1):69–79.

    Article  CAS  PubMed  Google Scholar 

  11. Wang X, Wang S, Zhou L, Yu L, Zhang L. A network-pathway based module identification for predicting the prognosis of ovarian cancer patients. J Ovarian Res. 2016; 9(1):73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Wu G, Stein L. A network module-based method for identifying cancer prognostic signatures. Genome Biol. 2012; 13(12):1–14.

    Article  Google Scholar 

  13. Hofree M, Shen JP, Carter H, Gross AM, Ideker T. Network-based stratification of tumor mutations. Nat Methods. 2013; 10(11):1108–15.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. Yang R, Daigle BJ, Petzold LR, Doyle FJ. Core module biomarker identification with network exploration for breast cancer metastasis. BMC Bioinformatics. 2012; 13(1):12.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Chuang H, Lee E, Liu Y, Lee D, Ideker T. Network-based classification of breast cancer metastasis. Mol Syst Biol. 2007; 3(1):140.

    PubMed  PubMed Central  Google Scholar 

  16. Hong S, Dong H, Jin L, Xiong M. Gene co-expression network and functional module analysis of ovarian cancer. Int J Comput Biol Drug Des. 2011; 4(2):147.

    Article  PubMed  Google Scholar 

  17. Lee J-H, Zhao X-M, Yoon I, et al. Integrative analysis of mutational and transcriptional profiles reveals driver mutations of metastatic breast cancers. Cell Discov. 2016; 2:16025.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Liu K-Q, Liu Z-P, Hao J-K, et al. Identifying dysregulated pathways in cancers from pathway interaction networks. BMC Bioinformatics. 2012;13(126).

  19. Zhao X-M, Liu K-Q, Zhu G, et al. Identifying cancer-related micrornas based on gene expression data. Bioinformatics. 2015; 31(8):1226–34.

    Article  PubMed  Google Scholar 

  20. Song W-MS, Bin Z. Multiscale embedded gene co-expression network analysis. PLoS Comput Biol. 2015; 11(11):1004574.

    Article  Google Scholar 

  21. Gao C, McDowell IC, Zhao S, Brown CD, Engelhardt BE. Context specific and differential gene co-expression networks via bayesian biclustering. PLoS Comput Biol. 2016; 12(7):1004791.

    Article  Google Scholar 

  22. Zhu L, Ding Y, Chen C-Y, et al. Metadcn: meta-analysis framework for differential co-expression network detection with an application in breast cancer. Bioinformatics. 2017; 33(8):1121–9.

    PubMed  Google Scholar 

  23. Huang H, Chuang Y, Chen C. Affinity aggregation for spectral clustering. In: IEEE conference on Computer Vision and Pattern Recognition (CVPR) 2012: 2012. p. 773–80.

  24. Zhang S, Zhao H, Ng MK. Functional module analysis for gene coexpression networks with network integration. IEEE/ACM Trans Comput Biol Bioinforma. 2015; 12(5):1146–60.

    Article  CAS  Google Scholar 

  25. Li W, Liu C, Zhang T, Li H, Waterman MS, Zhou XJ. Integrative analysis of many weighted co-expression networks using tensor computation. PLoS Comput Biol. 2011;7(6).

  26. Von Luxburg U. A tutorial on spectral clustering. Stat Comput. 2007; 17(4):395–416.

    Article  Google Scholar 

  27. Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al. Gene ontology: tool for the unification of biology. Nat Genet. 2000; 25(1):25–9.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Kanehisa M, Goto S, Furumichi M, Tanabe M, Hirakawa M. Kegg for representation and analysis of molecular networks involving diseases and drugs. Nucleic Acids Res. 2010; 38:D355–D360.

    Article  CAS  PubMed  Google Scholar 

  29. Osborne JD, Flatow J, Holko M, Lin S, Kibbe WA, Zhu LJ, Danila M, Feng G, Chisholm RL. Annotating the human genome with disease ontology. BMC Genomics. 2009; 10(1):1–8.

    Article  Google Scholar 

  30. Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich AG, Pomeroy SL, Golub TR, Lander ES, et al. Gene set enrichment analysis: A knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci U S A. 2005; 102(43):15545–50.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  31. Yu G, Wang L, Han Y, He Q. clusterprofiler: an r package for comparing biological themes among gene clusters. Omics A J Int Biol. 2012; 16(5):284–7.

    Article  CAS  Google Scholar 

  32. Zhu X, Mulcahy LA, Mohammed R, Lee AHS, Franks HA, Kilpatrick L, Yilmazer A, Paish EC, Ellis IO, Patel PM, et al. Il-17 expression by breast-cancer-associated macrophages: Il-17 promotes invasiveness of breast cancer cell lines. Breast Cancer Res. 2008; 10(6):1–11.

    Article  Google Scholar 

  33. Kamel MW, Shouman S, Elmerzebany M, Kilic GS, Veenstra TD, Saeed M, Wagih ME, Diazarrastia CR, Patel D, Salama SA. Effect of tumour necrosis factor-alpha on estrogen metabolic pathways in breast cancer cells. J Cancer. 2012; 3(1):310–21.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  34. Tian G, Liang J, Wang Z, Zhou D. Breast cancer risk in rheumatoid arthritis: An update meta-analysis. BioMed Res Int. 2014; 2014:453012.

    PubMed  PubMed Central  Google Scholar 

  35. Guille A, Chaffanet M, Birnbaum D. Signaling pathway switch in breast cancer. Cancer Cell Int. 2013; 13(1):66.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Luch A. Nature and nurture lessons from chemical carcinogenesis. Nat Rev Cancer. 2005; 5(2):113–25.

    Article  CAS  PubMed  Google Scholar 

  37. Colotta F, Allavena P, Sica A, Garlanda C, Mantovani A. Cancer-related inflammation, the seventh hallmark of cancer: links to genetic instability. Carcinogenesis. 2009; 30(7):1073–81.

    Article  CAS  PubMed  Google Scholar 

  38. Adam JK, Odhav B, Bhoola KD. Immune responses in cancer. Pharmacol Ther. 2003; 99(1):113–32.

    Article  CAS  PubMed  Google Scholar 

  39. Van Agthoven T, Van Agthoven TLA, Dekker A, Der Spek PJV, Vreede L, Dorssers LCJ. Identification of bcar3 by a random search for genes involved in antiestrogen resistance of human breast cancer cells. EMBO J. 1998; 17(10):2799–808.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  40. Kontos CK, Scorilas A. Kallikrein-related peptidases (klks): a gene family of novel cancer biomarkers. Clin Chem Lab Med. 2012; 50(11):1877–91.

    Article  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

S. Zhang’s research is supported in part by NSFC grant No.11471082, Science and Technology Commission of Shanghai Municipality 16JC1402600.

Funding

The publication of this work is supported by NSFC grant No. 11471082.

Availability of data and materials

The level 3 gene expression data for breast cancer (BRCA) and ovarian cancer (OV) were downloaded from The Cancer Genome Atlas (TCGA). The gene expression data are generated with UNC AgilentG4502A. The most updated protein-protein interaction (PPI) data for humans were downloaded from BioGrid https://thebiogrid.org.

About this supplement

This article has been published as part of BMC Systems Biology Volume 12 Supplement 1, 2018: Selected articles from the 16th Asia Pacific Bioinformatics Conference (APBC 2018): systems biology. The full contents of the supplement are available online at https://bmcsystbiol.biomedcentral.com/articles/supplements/volume-12-supplement-1.

Author information

Authors and Affiliations

Authors

Contributions

Not applicable.

Corresponding author

Correspondence to Shuqin Zhang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The author declares that she has no competing interests.

Publisher’s Note

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

Additional files

Additional file 1

The identified modules. We listed all the identified modules. The second column in each sheet lists the gene names. (XLSX 337 kb)

Additional file 2

Statistics of the modules. The module size, average degree, density, number of enriched GO, KEGG, DO terms are listed. The module names are the sheet name in Additional file 1. There is a map between the module names and the module number in the paper. (CSV 5 kb)

Additional file 3

GO enrichment results. This file lists all the GO-BP enriched terms for each module. ‘geneID’ is the Entrez gene id. (XLSX 193 kb)

Additional file 4

KEGG enrichment results. This file lists all the KEGG enriched terms for each module. ‘geneID’ is the Entrez gene id. (XLSX 59 kb)

Additional file 5

DO enrichment results. This file lists all the DO enriched terms for each module. ‘geneID’ is the Entrez gene id. (XLSX 54 kb)

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

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Zhang, S. Comparisons of gene coexpression network modules in breast cancer and ovarian cancer. BMC Syst Biol 12 (Suppl 1), 8 (2018). https://doi.org/10.1186/s12918-018-0530-9

Download citation

  • Published:

  • DOI: https://doi.org/10.1186/s12918-018-0530-9

Keywords