Next Article in Journal
Improved Production and Biophysical Analysis of Recombinant Silicatein-α
Next Article in Special Issue
Fully-Connected Neural Networks with Reduced Parameterization for Predicting Histological Types of Lung Cancer from Somatic Mutations
Previous Article in Journal
Caveolin-3: A Causative Process of Chicken Muscular Dystrophy
Previous Article in Special Issue
The Development of a Skin Cancer Classification System for Pigmented Skin Lesions Using Deep Learning
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Colorectal Cancer Prediction Based on Weighted Gene Co-Expression Network Analysis and Variational Auto-Encoder

1
Basic Experimental Center of Natural Science, University of Science and Technology Beijing, Beijing 100083, China
2
School of Mathematics and Physics, University of Science and Technology Beijing, Beijing 100083, China
*
Author to whom correspondence should be addressed.
Biomolecules 2020, 10(9), 1207; https://doi.org/10.3390/biom10091207
Submission received: 11 July 2020 / Revised: 30 July 2020 / Accepted: 1 August 2020 / Published: 20 August 2020
(This article belongs to the Special Issue Application of Artificial Intelligence for Medical Research)

Abstract

:
An effective feature extraction method is key to improving the accuracy of a prediction model. From the Gene Expression Omnibus (GEO) database, which includes 13,487 genes, we obtained microarray gene expression data for 238 samples from colorectal cancer (CRC) samples and normal samples. Twelve gene modules were obtained by weighted gene co-expression network analysis (WGCNA) on 173 samples. By calculating the Pearson correlation coefficient (PCC) between the characteristic genes of each module and colorectal cancer, we obtained a key module that was highly correlated with CRC. We screened hub genes from the key module by considering module membership, gene significance, and intramodular connectivity. We selected 10 hub genes as a type of feature for the classifier. We used the variational autoencoder (VAE) for 1159 genes with significantly different expressions and mapped the data into a 10-dimensional representation, as another type of feature for the cancer classifier. The two types of features were applied to the support vector machines (SVM) classifier for CRC. The accuracy was 0.9692 with an AUC of 0.9981. The result shows a high accuracy of the two-step feature extraction method, which includes obtaining hub genes by WGCNA and a 10-dimensional representation by variational autoencoder (VAE).

1. Introduction

Colorectal cancer (CRC) is a malignant tumor that ranks third for incidence and second for mortality worldwide [1]. Despite advances in medical technology, most CRC patients go to the hospital because of pain when the cancer is in the middle and late stage. These statistics call for an effective early diagnosis and novel prognostic markers for CRC. Classification models for an accurate screening and diagnosis are required for precision medicine [2]. An effective feature extraction method can avoid the problem of overfitting [3]. The feature extraction method is a key factor which affects the performance of the classifier so as to improve the accuracy of microarray gene expression datasets in predicting CRC.
Researchers have presented many methods for predicting CRC, such as SVM, logistic regression model, discriminative deep belief networks (DDBN), and so on. Zhao et al. selected factors such as the weight, tumor type, and tumor grade as the classification features of a classifier based on logistic regression (LR) and a support vector machine (SVM) [4]. Agesen et al. selected the gene expression data of 13 genes closely related to CRC as classification features to construct a ColoGuideEx classifier, so as to achieve the effective classification of stage II colorectal cancer [5]. Gabere et al. used the gene expression data of 30 genes closely related to CRC as classification features of SVM for predicting CRC, the accuracy of classifier meeting 0.95 [6]. Cubiella et al. applied 11 factors, such as age, gender, faecal haemoglobin, and carcinoembryonic antigen to a logistic regression model for screening CRC [7]. Karabulut et al. constructed discriminative deep belief networks (DDBN) with gene expression profile data for colorectal cancer prediction [8]. Yong et al. applied the miRNA expression profile data of differential genes in cancer tissue and blood samples to a logistic regression model for predicting CRC [9].
To screen and diagnose cancer, researchers usually focus on genes whose expression are significantly different between cancer samples and normal samples. Pathogenic cancer genes are usually correlated, and the occurrence of cancer is the result of an abnormal expression of multiple genes [10]. If the expression of certain genes in a physiological process is similar between different tissues, then these genes are functionally related, and thus they can be treated as a whole group [11]. Gene modules are made by grouping genes according to relevance, and genes with a higher correlation were, accordingly, classified as one type by weighted gene co-expression network analysis (WGCNA). WGCNA is similar to cluster analysis, but WGCNA is a more biologically significant statistical method. WGCNA can be used for finding clusters (modules) of highly correlated genes. It can also be used to find modules that are closely related to disease. WGCNA has been used to find genes related to lung squamous cell carcinoma and colorectal cancer [12,13].
Microarray gene expression data have the characteristics of a high dimension and small sample size, and there are also a large number of redundant genes that will cause a curse of dimensionality, interfere with the diagnosis, and reduce the accuracy of a classifier. Thus, it is important to reduce the dimensions of gene expression data. Common methods to reduce the dimensionality are a principal component analysis [14], Fisher’s linear discriminant analysis [15], and autoencoder [16]. Autoencoder is a type of unsupervised learning technology that uses neural networks to learn to express low-dimension features from high-dimensional features. In many tasks, the autoencoder has a better effect than principal component analysis. There are many improved autoencoder algorithms, such as stack autoencoder [17], sparse autoencoder [18], denoising autoencoder [19], and variational autoencoder (VAE). AE just maps the data to a function. VAE map the data information to a distribution, and the distribution can effectively summarize the information on this kind of data [20]. VAE successfully reduced the dimension of microarray gene expression data, and the effective features closely related to cancer were extracted from microarray gene expression data [21,22].
We used WGCNA to construct co-expressed gene modules and identified a key gene module that was most closely related to cancer. We analyzed the hub genes in this key module. Hub genes may serve as effective early diagnostic and prognostic markers for CRC. We used the hub genes as a type of feature classifier. We applied the VAE to genes with significantly different expression data, and we mapped the results to a 10-dimensional representation as another type of cancer feature classifier.

2. Materials and Methods

2.1. CRC Microarray Datasets

We downloaded the gene expression datasets GSE8671 [23], GSE9348 [24], GSE23878 [25], and GSE37364 [26] from the GPL570 (Affymetrix Human Genome U133 Plus 2.0 Array) platform, and these included 106 normal samples and 132 colorectal cancer samples (Supplementary Information Table S1).

2.2. Datasets Preprocessing

For repeated probe values, we calculated their median as the new value of the probe to ensure the uniqueness of the probe. The uploaders of the four datasets normalized each chip. Two datasets, GSE23878 and GSE37364, were logarithmically processed. We first processed the dataset without logarithm; all data were the result of l o g 2 ( c o u n t ) , with the count representing the amount of gene expression [27]. Each dataset in the GEO database has a small sample size, and statistical analysis performed on small samples is usually not robust [28]. Thus, it is necessary to perform a joint analysis of the gene expression data from many different sources. The four expression datasets were combined into one dataset according to the same gene symbols, and the expression data of 13,487 genes from 238 samples were obtained.
However, different conditions can lead to batch effects [29]. A number of approaches can be used to remove batch effects between different datasets, such as the Empirical Bayes method (ComBat) [30], singular value decomposition (SVD) [31], and distance weighted discrimination (DWD) [32]. The empirical Bayesian algorithm can adjust the batch effect in microarray expression data [33].
After combining the four sets of GSE data, we used the sva package in R language based on the empirical Bayesian framework [34] to remove the batch effect for the subsequent statistical analysis.
To remove the batch effect, we assume that the measured expression value of gene i in sample j of batch g can be expressed in a general form as follows:
x i j = x i j + b i j x + ε i j x
Then, we standardize it:
x i j ^ = x i j x i σ x i
In the above formula, x i is the average value of the expression value from gene i in all samples; σ x i is the standard deviation for x i of the expression value from gene i in all samples.

2.3. Screening of Differentially Expressed Genes

An important task of gene expression profiling analysis is to perform a variable screening based on known data in order to find genes with significantly different expressions between samples. The typical practice is to use a combination of fold change [35] and T test [36]. The most commonly used tools for differential gene expression analysis are DESeq2 [37], EdgeR [38], and limma [39]. We used the limma package, the P-values were corrected using the FDR correction toolkit in R language. We selected genes with significant differences in choice (adj. p < 0.05) and more than double the difference in gene expression between normal and cancer tissues ( | log 2 F C | > 1 ) to serve as differentially expressed genes.

2.4. Construction of Co-Expression Network and Identification of Key Modules

The scale-free gene network condition should be satisfied before conducting WGCNA. A hierarchical clustering tree (module) is built based on the correlation matrix of the gene co-expression, adjacency function, and the dissimilarity measurements of different nodes. The dynamic tree cut method can identify more accurately and more biologically significant co-expression modules [40]. A key module with a high correlation with a specific phenotype or disease can be extracted from these co-expression modules. The hub genes are identified according to the internal connectivity of genes in the key module and the correlation between genes and feature vectors of the key module [41,42].

2.5. Mining Hub Genes from the Key Module

The module membership (MM) was calculated via the WGCNA function signedKME, which correlated the module eigengene (ME) with gene expression values, thereby quantifying how close a gene was to a given module. The correlation between individual genes and a biological trait (cancer or normal) was defined as the gene significance (GS). The summation of adjacency performed for all genes in a particular network was calculated as the intramodular connectivity (K.in).
Hub genes in a module usually have high GS, high MM, and high K.in [43,44]. Therefore, we first selected genes with GS > 0.5 and MM > 0.8. Then, we screened for genes with significantly different expressions in normal and colorectal cancer samples from the genes that met the above conditions. Finally, we sorted the selected genes into descending order according to the K.in value and considered the top ten genes as the hub genes for this analysis. This approach provided a new method for reducing the dimension issues that often accompany, and potentially interfere with, cancer prediction and biological significance.

2.6. Dimension Reduction with VAE

Autoencoder is a type of artificial neural network used to learn the compressed representation for a set of data. This neural network has the same number of neuron nodes in the input layer and output layer, but fewer in the hidden layer. The hidden layer variable of the VAE can learn the distribution of the original data to ensure that the hidden layer can better abstract the characteristics of the input data. Variable inference and reparameterization are used to reduce the dimension of gene expression data. We extracted low-dimensional information from gene expression data with VAE and used this low-dimensional information as a type of feature classifier for screening cancer.
Our pipeline for using WGCNA in order to identify hub genes of colorectal cancer and for using VAE in order to identify a 10-dimensional representation can be freely accessed from GitHub at https://github.com/gutmicrobes/WGCNA-VAE.git.

2.7. SVM

Support Vector Machine (SVM) is a class of generalized linear classifier, which has good robustness, with supervised learning that carries out a binary classification of data [45]. The advantages of SVM in predicting are quantitative and qualitative [46]. Researchers have used linear SVM to predict CRC [4]. We analyzed three data sets, GSE8671, GSE9348, and GSE23878, and we selected 10 hub genes and the 10-dimensional representation from the variational autoencoder (VAE) for 1166 genes with significantly different expressions as the classifier for CRC.

3. Results and Discussion

3.1. Extraction of Differential Genes

We integrated four datasets as the sample set for the final calculations. To test for the existence of batch effects and whether the four dataset samples were completely separated or whether they had some continuity, we performed a principal components analysis (PCA) on all gene expression data (Figure 1A). The four sample sets were not completely continuous; in particular, dataset GSE23878 was significantly different from the other datasets. We then performed a batch effect correction, and a PCA analysis was performed again on the corrected data. The results in Figure 1B show that the separation among the four datasets after the correction was essentially eliminated.
We divided gene expression data that came from the datasets GSE8671, GSE9348, and GSE23878 into normal and cancer groups for a differential analysis. A volcano plot shown in Figure S1 of the Supplementary Information shows the distribution of genes expressed between cancer and normal tissue.
An important task in gene expression profiling is to perform variable screening based on known data in order to find meaningful differential genes. Accordingly, we performed a differential analysis on the gene expression profiling data after removing the batch effect, as described above. We identified 1159 significantly altered genes from four GEO datasets, of which 419 genes were upregulated and 740 genes were downregulated.

3.2. Extracting Hub Gene by Weighted Correlation Network Analysis

3.2.1. Soft Threshold Screening

We used the pickSoftThreshold grid search function to calculate the β value of the scale-free network. Table S2 in the Supplementary Information summarizes the parameters corresponding to β . Included are the square of the correlation coefficient between log ( p ( k ) ) and log ( k ) ( S E T . R . s q ) , slope and correlation coefficient of the truncation index model (truncated.R.sq), average degree of connectivity (mean.k), median degree of connectivity (median.k), and maximum degree of connectivity (max.k). The main parameters we considered in the table were SET.R.sq and mean.k.
As the soft threshold increased, the SET.R.sq and mean.k showed corresponding changes (Supplementary Information, Figure S2). Generally, SET.R.sq must be greater than 0.85; that is, the network was considered as meeting the requirement of a no scale distribution. When β = 12 , the network met the requirements of a scale-free distribution (Figure S2). Therefore, we set the threshold β to 12.

3.2.2. Gene Module Construction

We calculated the topological matrix and the adjacency matrix according to the soft threshold and then transformed the topological matrix into a dissimilarity matrix. We clustered according to the degree of dissimilarity in order to obtain the system clustering tree. In the clustering tree, each “twig” corresponded to a module, and each piece of “leaf” represented a gene in the module. All samples were clustered so as to form 17 modules, the smallest of which contained 36 genes (Figure 2A). At most, the modules contained 6377 genes, with an average of 814 genes per module.
Each co-expressed gene module was a set of genes with similar expression patterns, and all gene expression amounts in the module were then combined in order to obtain a “characteristic value,” which we called a module characteristic value. The module eigenvalue was the result of the operation of the common eigenvalues of all genes in the module. We used the module eigenvalue to calculate the correlation, and the thermal map is plotted in Figure 2B.
Figure 2B shows a certain correlation between different modules. Modules with a high correlation can be further merged into the same module. Accordingly, we used the dynamic pruning method to further merge modules with a significant correlation into the same module. To form 12 modules, we defined a minimum of 30 genes per module and combined modules with a correlation coefficient greater than 0.8. Table 1 shows the 12 modules and the number of genes included.
To analyze the relationship between gene modules and cancer, we quantified cancer patients as 1, normal people as 0, and we calculated the ME for each module; all genes in the module are represented by ME. We calculated the Pearson correlation coefficient (PCC) between ME and cancer (Table 2) [47]. With the exception of the midnight blue module, all other modules were significantly correlated with cancer. The results of quantitatively analyzing the correlation between cancer and the module eigengene indicated that the gene module division after pruning had a biological significance.
In Table 2, a large PCC and a small p-value indicated a stronger relationship between the module and cancer. The PCC between the MEturquoise module and cancer was the highest, and the p-value was the lowest, making the correlation biologically significant. These genes were the most important in the module that distinguished cancer from normal tissue.

3.2.3. Hub Gene Identification

According to the coefficient between ME and cancer, MEturquoise became our module of interest, and it contained 6865 genes. Therefore, we further mined the hub genes from this module.
A hub gene in a module usually has a high GS, high MM, and high K.in. Therefore, we first selected genes with GS > 0.5 and MM > 0.8. These conditions were met by 191 genes, of which 165 had significant differences between cancer and normal samples.
Then, to explore the functional association between gene modules and colorectal cancer, we performed a Gene Ontology and KEGG pathway enrichment analysis on the aforesaid 165 genes (see Supplementary Information, Figure S3 and Figure S4).
Finally, we selected the top 10 genes with the largest K.in as the hub genes for this study (Table 3).
GUCA2B (Guanylate Cyclase Activator 2B), a protein coding gene, is a physiological regulator of intestinal fluid and electrolyte transport. Nagaraj et al. considered GUCA2B as a noninvasive biomarker for the early detection of colorectal cancer [48]. Gene Ontology annotations for this gene include calcium-sensitive guanylate cyclase activator activity and guanylate cyclase activator activity. The GUCA2A gene is an important paralog of GUCA2B.
CDK4 (cyclin-dependent kinase 4) is important for cell cycle G1 phase progression. It is a regulator of the cell cycle. Lee et al. use CDK4/6 inhibitors to inhibit the proliferation of colorectal cancer cells effectively [49].
TRIP13 (thyroid hormone receptor interactor 13) messenger RNA was highly expressed in multiple CRC tissues. The depletion of TRIP13 in CRC cells suppressed cell proliferation, migration, and invasion [50].
As for EIF3B (eukaryotic translation initiation factor 3 subunit B), the researchers found that the silencing of EIF3B gene expression could significantly inhibit the colorectal cancer cell proliferation rate and clonability [51].
The CDH3 gene encodes cadherin, a type of cell adhesion molecule that is important in the formation of adherens junctions that bind cells [52]. CDH3 is located in a gene cluster of chromosome 16 that is involved in breast cancer and prostate cancer. An abnormal expression of CDH3 protein is also observed in cervical cancer. This gene is also associated with adolescent macular dystrophy and poor ectodermal release. CDH3 is overexpressed in colorectal tumors, and it has potential as a serum marker for colorectal cancer surveillance [53].
GUCA2A (Guanylate Cyclase Activator 2A) is a protein-coding gene. Among its related pathways are the Metabolism and Myometrial Relaxation and Contraction Pathways. Gene Ontology annotations for this gene include hormonal activity and guanylate cyclase activator activity. Zhang et al. revealed that GUCA2A was downregulated in CRC tissues [54].

3.3. Realization of VAE

We designed a VAE model to extract low-dimensional information from gene expression data as a type of feature for the analysis and prediction of cancer. There were more nodes between the input layer and the first hidden layer, with 680,000 parameters needing training. If we had extracted features directly from more than 10,000 genes with expression data, we would have needed to reduce the number of hidden layers of the model and the number of neuron nodes in the hidden layer, thereby reducing accuracy. Therefore, the input data were genes with significant differences. We adopted a neural network model with an input-output layer and five hidden layers. In the VAE, the numbers of hidden layer nodes in the model were 584, 100, 10, 100, and 584, and the input layer and output layer nodes were 1159. The learning rate was 0.0005, batch size was 20, and epoch was 6. To measure the dimensionality reduction effect of the model, the last 65 samples of the gene expression data were used as the test set, and all the other samples were used as the training set. We determined that the appropriate epoch to ensure the stable convergence of the model was six by calculating the reconstruction error and the loss function of the model on the test set. Finally, we used a VAE model to successfully scale the gene expression data into 10 dimensions and use them as one of the characteristics of the classifier.

3.4. Analysis Results from the Model of Cancer Prediction

We used the classifier SVM to identify a classification method suitable for the dataset described here. We used three datasets as the training set: GSE8671, GSE9348, and GSE23878; meanwhile, the dataset GSE37364 was the test set. The training set obtained 68 healthy samples and 105 cancer samples, and the test set obtained 38 healthy samples and 27 cancer samples. The training set and the test set basically balanced the positive and negative samples. SVMs are more flexible in binary tasks and have a variety of kernel functions to choose from. A new taxonomic data set was constructed by combining the hub genes with the expression of the low-dimensional features of the VAE. SVM was used to classify the new classification dataset, with a test accuracy of 0.9692 and AUC of 0.9981. The SVM classified the 65 samples in the test set exactly. The two-step dimensionality reduction method greatly improved the classification accuracy.

4. Conclusions

We obtained 1166 genes with significant differences in gene expression from 238 normal and cancer samples. The number of genes obtained by differential analysis was relatively large, and many of these genes interact. We used WGCNA to construct a scale-free network of gene co-expression data. WGCNA divided genes into different modules and found that the genes in each module were similar in expression. On the basis of a correlation analysis between the module eigengene and the samples, we identified the MEturquoise module, which may have an important function in carcinogenesis. Our results showed that the genes in this module were different from those in the other modules. That is, GUCA2A, GUCA2B, CDH3, and other hub genes were highly correlated with the occurrence of CRC. We used 10 hub genes as one type of feature for CRC prediction. We also used VAE to reduce the gene spectrum data to 10 dimensions and derive another type of feature for the CRC prediction model. The accuracy of the CRC prediction model was 0.9692.

Supplementary Materials

The following are available online at https://www.mdpi.com/2218-273X/10/9/1207/s1, Figure S1: Volcano map of differential genes between colorectal cancer sample and normal sample. Figure S2: Soft threshold Screening. Figure S3: GO Functional enrichment analysis. Figure S4: KEGG function enrichment analysis. Table S1: Sample distribution of four datasets. Table S2: Soft Threshold Screening for Network Construction.

Author Contributions

D.A. and Y.W. conceived and designed experiments. Y.W. and X.L. performed the experiments and analyzed data. D.A. supervised the study and manuscript writing with inputs from X.L., and H.P. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by grants from the National Natural Science Foundation of China (61873027, 61370131).

Acknowledgments

DA thanks Fengzhu Sun at the University of Southern California and Li C. Xia at Stanford University for their support and helpful discussions.

Conflicts of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

  1. Bray, F.F.J.; Soerjomataram, I.; Siegel, R.L.; Torre, L.A.; Jemal, A. Global cancer statistics 2018: GLOBOCAN estimates of incidence and mortality worldwide for 36 cancers in 185 countries. CA Cancer J. Clin. 2018, 68, 394–424. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Ai, L.; Tian, H.; Chen, Z.; Chen, H.; Xu, J.; Fang, J.Y. Systematic evaluation of supervised classifiers for fecal microbiota-based prediction of colorectal cancer. Oncotarget 2017, 8, 9546–9556. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Martin, M.L.; Nevado, A.; Carro, B. Detection of early stages of Alzheimer’s disease based on MEG activity with a randomized convolutional neural network. Artif. Intell. Med. 2020, 107, 101924. [Google Scholar]
  4. Zhao, D.D.; Liu, H.; Zheng, Y.J.; He, Y.L.; Lu, D.J.; Chen, L. A reliable method for colorectal cancer prediction based on feature selection and support vector machine. Med. Biol. Eng. Comput. 2019, 577, 901–912. [Google Scholar] [CrossRef] [PubMed]
  5. Agesen, T.H.; Sevvn, A.; Lind, G.E.; Nesbakken, A.; Skotheim, R.L.; Lothe, R.A. ColoGuideEx: A robust gene classifier specific for stage II colorectal cancer prognosis. Gut 2012, 61, 1560–1567. [Google Scholar] [CrossRef]
  6. Gabere, M.N.; Hussein, N.A.; Aziz, M.A. Filtered selection coupled with support vector machines generate a functionally relevant prediction model for colorectal cancer. Oncotargets Ther. 2016, 9, 3313–3325. [Google Scholar]
  7. Cubiella, J.; Vega, P.; Salve, M.; Ondina, M.D.; Alves, M.T.; Quintero, E.; Victoria, Á.S.; Fernando, F.B.; Boadas, J.; Campo, R. Development and external validation of a faecal immunochemical test-based prediction model for colorectal cancer detection in symptomatic patients. BMC Med. 2016, 14, 1–13. [Google Scholar] [CrossRef] [Green Version]
  8. Karabulut, E.M.; Ibrikci, T. Discriminative deep belief networks for microarray based cancer classification. Biomed. Res. 2017, 28, 1016–1024. [Google Scholar]
  9. Yong, F.L.; Law, C.W.; Wang, C.W. Potentiality of a triple microRNA classifier: miR-193a-3p, miR-23a and miR-338-5p for early detection of colorectal cancer. BMC Cancer 2013, 13, 280. [Google Scholar] [CrossRef] [Green Version]
  10. Bärlund, M.; Monni, O.; Kononen, J.; Cornelison, R.; Torhorst, J.; Sauter, G.; Kallioniemi, O.P.; Kallioniemi, A. Multiple genes at 17q23 undergo amplification and overexpression in breast cancer. Cancer Res. 2000, 60, 5340–5344. [Google Scholar]
  11. Carlson, M.R.; Zhang, B.; Fang, Z.; Mischel, P.S.; Horvath, S.; Nelson, S.F. Gene connectivity, function, and sequence conservation: Predictions from modular yeast co-expression networks. BMC Genom. 2006, 7, 40. [Google Scholar] [CrossRef] [PubMed]
  12. Tian, F.; Zhao, J.; Fan, X.; Kang, Z. Weighted gene co-expression network analysis in identification of metastasis-related genes of lung squamous cell carcinoma based on the Cancer Genome Atlas database. J. Thorac. Dis. 2017, 9, 42. [Google Scholar] [CrossRef] [Green Version]
  13. Qin, L.; Zeng, J.; Shi, N.; Chen, L.; Wang, L. Application of Weighted Gene co-expression Network Analysis to Explore the Potential Diagnostic Biomarkers for Colorectal Cancer. Mol. Med. Rep. 2020, 21, 2533–2543. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Lenz, M.; Müller, F.-J.; Zenke, M.; Schuppert, A. Principal components analysis and the reported low intrinsic dimensionality of gene expression microarray data. Sci. Rep. 2016, 6, 25696. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Huerta, E.B.; Duval, B.; Hao, J.K. A hybrid LDA and genetic algorithm for gene selection and classification of microarray data. Neurocomputing 2010, 73, 2375–2383. [Google Scholar] [CrossRef] [Green Version]
  16. Wang, Y.; Yao, H.; Zhao, S. Auto-encoder based dimensionality reduction. Neurocomputing 2016, 184, 232–242. [Google Scholar] [CrossRef]
  17. Shin, H.C.; Orton, M.R.; Collins, D.J.; Doran, S.J.; Leach, M.O. Stacked autoencoders for unsupervised feature learning and multiple organ detection in a pilot study using 4D patient data. IEEE Trans. Pattern Anal. Mach. Intell. 2012, 35, 1930–1943. [Google Scholar] [CrossRef]
  18. Ng, A. Sparse autoencoder. CS294A Lect. Notes 2011, 72, 1–19. [Google Scholar]
  19. Vincent, P.; Larochelle, H.; Bengio, Y.; Manzagol, P.A. Extracting and Composing Robust Features with Denoising Autoencoders. In Proceedings of the 25th International Conference on Machine Learning, Association for Computing Machinery, New York, NY, USA, 2014. [Google Scholar]
  20. Kingma, D.P.; Welling, M. Auto-Encoding Variational Bayes. In Proceedings of the International Conference on Learning Representations, Banff, AB, Canada, 14–16 April 2014. [Google Scholar]
  21. Chatrian, A.; Sirinukunwattana, K.; Verrill, C.; Rittscher, J. Towards the Identification of Histology Based Subtypes in Prostate Cancer. In Proceedings of the International Symposium on Biomedical Imaging, Venice, Italy, 24–27 April 2019. [Google Scholar]
  22. Wang, Z.X.; Wang, Y.D. Extracting a biologically latent space of lung cancer epigenetics variational autoencoders. BMC Bioinform. 2019, 20, 568. [Google Scholar] [CrossRef]
  23. Sabates-Bellver, J.; Van der Flier, L.G.; de Palo, M.; Cattaneo, E.; Maake, C.; Rehrauer, H.; Laczko, E.; Kurowski, M.A.; Bujnicki, J.M.; Menigatti, M.; et al. Transcriptome profile of human colorectal adenomas. Mol. Cancer Res. 2007, 5, 1263–1275. [Google Scholar] [CrossRef] [Green Version]
  24. Hong, Y.; Downey, T.; Eu, K.W.; Koh, P.K.; Cheah, P.Y. A ‘metastasis-prone’signature for early-stage mismatch-repair proficient sporadic colorectal cancer patients and its implications for possible therapeutics. Clin. Exp. Metastasis 2010, 27, 83–90. [Google Scholar] [CrossRef] [PubMed]
  25. Uddin, S.; Ahmed, M.; Hussain, A.; Abubaker, J.; Al-Sanea, N.; AbdulJabbar, A.; Ashari, L.H.; Alhomoud, S.; Al-Dayel, F.; Jehan, Z.; et al. Genome-wide expression analysis of Middle Eastern colorectal cancer reveals FOXM1 as a novel target for cancer therapy. Am. J. Pathol. 2011, 178, 537–547. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Valcz, G.; Patai, Á.V.; Kalmár, A.; Péterfia, B.; Fűri, I.; Wichmann, B.; Műzes, G.; Sipos, F.; Krenács, T.; Mihály, E. Myofibroblast-derived SFRP1 as potential inhibitor of colorectal carcinoma field effect. PloS ONE 2014, 9, E106143. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Cui, X.; Churchill, G.A. Statistical tests for differential expression in cDNA microarray experiments. Genome Biol. 2003, 4, 210. [Google Scholar] [CrossRef] [Green Version]
  28. Bevilacqua, V.; Pannarale, P.; Abbrescia, M.; Cava, C.; Paradiso, A.; Tommasi, S. Comparison of data-merging methods with SVM attribute selection and classification in breast cancer gene expression. BMC Bioinform. 2012, 13, S9. [Google Scholar] [CrossRef] [Green Version]
  29. Luo, J.; Schumacher, M.; Scherer, A.; Sanoudou, D.; Megherbi, D.; Davison, T.; Shi, T.; Tong, W.; Shi, L.; Hong, H.; et al. A comparison of batch effect removal methods for enhancement of prediction performance using MAQC-II microarray gene expression data. Pharm. J. 2010, 10, 278–291. [Google Scholar] [CrossRef] [Green Version]
  30. Johnson, W.E.; Li, C.; Rabinovic, A. Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics 2007, 8, 118–127. [Google Scholar] [CrossRef]
  31. Alter, O.; Brown, P.O.; Botstein, D. Singular value decomposition for genome-wide expression data processing and modeling. Proc. Natl. Acad. Sci. USA 2000, 97, 10101–10106. [Google Scholar] [CrossRef] [Green Version]
  32. Benito, M.; Parker, J.; Du, Q.; Wu, J.; Xiang, D.; Perou, C.M.; Marron, J.S. Adjustment of systematic microarray data biases. Bioinformatics 2004, 20, 105–114. [Google Scholar] [CrossRef] [Green Version]
  33. Stein, C.K.; Qu, P.; Epstein, J.; Buros, A.; Rosenthal, A.; Crowley, J.; Morgan, G.; Barlogie, B. Removing batch effects from purified plasma cell gene expression microarrays with modified ComBat. BMC Bioinform. 2015, 16, 63. [Google Scholar] [CrossRef] [Green Version]
  34. Leek, J.T.; Johnson, W.E.; Parker, H.S.; Jaffe, A.E.; Storey, J.D. The sva package for removing batch effects and other unwanted variation in high-throughput experiments. Bioinformatics 2012, 28, 882–883. [Google Scholar] [CrossRef] [PubMed]
  35. Gerhold, D.; Lu, M.; Xu, J.; Austin, C.; Caskey, C.T.; Rushmore, T. Monitoring expression of genes involved in drug metabolism and toxicology using DNA microarrays. Physiol. Genom. 2001, 5, 161–170. [Google Scholar] [CrossRef] [PubMed]
  36. Baldi, P.; Long, A.D. A Bayesian framework for the analysis of microarray expression data: Regularized t-test and statistical inferences of gene changes. Bioinformatics 2001, 17, 509–519. [Google Scholar] [CrossRef] [PubMed]
  37. Love, M.I.; Huber, W.; Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014, 15, 550. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Robinson, M.D.; McCarthy, D.J.; Smyth, G.K. EdgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics 2010, 26, 139–140. [Google Scholar] [CrossRef] [Green Version]
  39. Ritchie, M.E.; Phipson, B.; Wu, D.; Hu, Y.; Law, C.W.; Shi, W.; Smyth, G. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015, 43, E47. [Google Scholar] [CrossRef]
  40. Langfelder, P.; Zhang, B.; Horvath, S. Defining clusters from a hierarchical cluster tree: The Dynamic Tree Cut package for R. Bioinformatics 2007, 24, 719–720. [Google Scholar] [CrossRef]
  41. Zhang, B.; Horvath, S. A general framework for weighted gene co-expression network analysis. Stat. Appl. Genet. Mol. Biol. 2005, 4, 1. [Google Scholar] [CrossRef]
  42. Langfelder, P.; Horvath, S. WGCNA: An R package for weighted correlation network analysis. BMC Bioinform. 2008, 9, 559. [Google Scholar] [CrossRef] [Green Version]
  43. Lou, Y.; Tian, G.Y.; Song, Y.; Liu, Y.L.; Chen, Y.D.; Shi, J.P.; Yang, J. Characterization of transcriptional modules related to fibrosing-NAFLD progression. Sci. Rep. 2017, 7, 4748. [Google Scholar] [CrossRef] [Green Version]
  44. Hu, Y.; Pan, J.; Xin, Y.; Mi, X.; Wang, J.; Gao, Q.; Luo, H. Gene Expression Analysis Reveals Novel Gene Signatures Between Young and Old Adults in Human Prefrontal Cortex. Front. Aging Neurosci. 2018, 10, 259. [Google Scholar] [CrossRef] [PubMed]
  45. Foody, G.M.; Mathur, A. Toward intelligent training of supervised image classifications: Directing training data acquisition for SVM classification. Remote Sens. Environ. 2004, 1, 107–117. [Google Scholar] [CrossRef]
  46. Meeh, C.L.; Croshaw, R.; Crimm, H.; Miller, S.K.; Oroian, D.; Kowli, S.; Zhu, J.; Carver, W.; Wu, W.; Pena, E.A.; et al. A Gene Expression Classifier of Node-Positive Colorectal Cancer. Neoplasia 2009, 11, 1074–1083. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Pearson, K. Determination of the coefficient of correlation. Science 1909, 30, 23–25. [Google Scholar] [CrossRef]
  48. Nagaraj, S.H.; Reverter, A. A Boolean-based systems biology approach to predict novel genes associated with cancer: Application to colorectal cancer. BMC Syst. Biol. 2011, 5, 35. [Google Scholar] [CrossRef] [Green Version]
  49. Lee, S.M.; Helms, T.L.; Feng, N.; Chang, Q.E.; Tian, F.; Wu, J.Y.; Toniatti, C.; Heffernan, T.P.; Powis, G.; Kwong, L.N.; et al. Efficacy of the combination of MEK and CDK4/6 inhibitors in vitro and in vivo in KRAS mutant colorectal cancer models. Oncotarget 2016, 26, 39595–39608. [Google Scholar] [CrossRef] [Green Version]
  50. Kurita, K.; Maeda, M.; Mansour, M.A.; Kokuryo, T.; Uehara, K.; Yokoyama, Y.; Nagino, M.; Hamaguchi, M.; Senga, T. TRIP13 is expressed in colorectal cancer and promotes cancer cell invasion. Oncol. Lett. 2016, 12, 5240–5246. [Google Scholar] [CrossRef] [Green Version]
  51. Wang, Z.; Chen, J.; Sun, J.; Cui, Z.; Wu, H. RNA interference-mediated silencing of eukaryotic translation initiation factor 3, subunit B (EIF3B) gene expression inhibits proliferation of colon cancer cells. World J. Surg. Oncol. 2012, 10, 119–127. [Google Scholar] [CrossRef] [Green Version]
  52. Alimperti, S.; Andreadis, S.T. CDH2 and CDH11 act as regulators of stem cell fate decisions. Stem Cell Res. 2015, 14, 270–282. [Google Scholar] [CrossRef] [Green Version]
  53. Kumara, H.S.; Bellini, G.A.; Caballero, O.L.; Herath, S.A.; Su, T.; Ahmed, A.; Njoh, L.; Ckkic, V.; Whelan, R.L. P-Cadherin (CDH3) is overexpressed in colorectal tumors and has potential as a serum marker for colorectal cancer monitoring. Oncoscience 2017, 4, 139. [Google Scholar] [CrossRef] [Green Version]
  54. Zhang, H.; Du, Y.; Wang, Z.; Lou, R.; Wu, J.; Feng, J. Integrated Analysis of Oncogenic Networks in Colorectal Cancer Identifies GUCA2A as a Molecular Marker. Biochem. Res. Int. 2019, 2019, 1–13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. PCA analysis of four datasets before and after the batch effect correction. (A) PCA analysis of four datasets before batch effect correction. Dataset GSE23878 was significantly different from the other three datasets; (B) PCA analysis of four datasets after batch effect correction. The batch effects of the four datasets were basically eliminated after the correction.
Figure 1. PCA analysis of four datasets before and after the batch effect correction. (A) PCA analysis of four datasets before batch effect correction. Dataset GSE23878 was significantly different from the other three datasets; (B) PCA analysis of four datasets after batch effect correction. The batch effects of the four datasets were basically eliminated after the correction.
Biomolecules 10 01207 g001
Figure 2. (A) System clustering tree of all samples after classification. Part a is the clustering tree constructed by genes, part b is the gene module obtained by clustering, and part c is the gene module obtained by combining similar expression patterns; (B) The thermal map of the relationship between the eigenvalues of different modules.
Figure 2. (A) System clustering tree of all samples after classification. Part a is the clustering tree constructed by genes, part b is the gene module obtained by clustering, and part c is the gene module obtained by combining similar expression patterns; (B) The thermal map of the relationship between the eigenvalues of different modules.
Biomolecules 10 01207 g002
Table 1. Gene module and number of corresponding genes.
Table 1. Gene module and number of corresponding genes.
ColorTanBrownTurquoiseBlueGreenPurple
Number85279963772636400153
ColorBlackMagentaMidnight BlueRedCyanGrey60
Number3301756031811836
Table 2. The PCC between each module and cancer after pruning.
Table 2. The PCC between each module and cancer after pruning.
MEtanMEbrownMEturquoiseMEblueMEgreenMEpurple
PCC−0.1285−0.3052−0.9251−0.7075−0.20170.3457
p-value0.09200.00000.00000.00000.00780.0000
MEblackMEmagentaMEblueMEredMEcyanMEgrey60
PCC−0.4263−0.57530.06090.51270.29440.3082
p-value0.00000.00000.42600.00000.00000.0000
Table 3. Top 10 genes with the highest intramodular connectivity.
Table 3. Top 10 genes with the highest intramodular connectivity.
GENE NAMElogFCadj.P.ValGSMM.TurquoiseK.in
CDK41.30422.40 × 10−390.8128−0.9076933.8071
CDH36.39182.92 × 10−800.9417−0.9206922.1908
DKC11.35173.34 × 10−410.8226−0.9055918.2304
UBE2S1.90857.08 × 10−400.8118−0.8997906.0183
GUCA2B−6.44446.12 × 10−530.87230.9075895.3242
UBE2C2.14359.16 × 10−520.8706−0.9122895.1551
EIF3B1.31731.02 × 10−430.8368−0.8970894.6524
TRIP131.94102.07 × 10−460.8474−0.8905890.4061
GUCA2A−5.32012.29 × 10−500.86280.8628887.2683
GTF3A1.41265.21 × 10−390.8085−0.8173883.3120

Share and Cite

MDPI and ACS Style

Ai, D.; Wang, Y.; Li, X.; Pan, H. Colorectal Cancer Prediction Based on Weighted Gene Co-Expression Network Analysis and Variational Auto-Encoder. Biomolecules 2020, 10, 1207. https://doi.org/10.3390/biom10091207

AMA Style

Ai D, Wang Y, Li X, Pan H. Colorectal Cancer Prediction Based on Weighted Gene Co-Expression Network Analysis and Variational Auto-Encoder. Biomolecules. 2020; 10(9):1207. https://doi.org/10.3390/biom10091207

Chicago/Turabian Style

Ai, Dongmei, Yuduo Wang, Xiaoxin Li, and Hongfei Pan. 2020. "Colorectal Cancer Prediction Based on Weighted Gene Co-Expression Network Analysis and Variational Auto-Encoder" Biomolecules 10, no. 9: 1207. https://doi.org/10.3390/biom10091207

APA Style

Ai, D., Wang, Y., Li, X., & Pan, H. (2020). Colorectal Cancer Prediction Based on Weighted Gene Co-Expression Network Analysis and Variational Auto-Encoder. Biomolecules, 10(9), 1207. https://doi.org/10.3390/biom10091207

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop