Single-cell and WGCNA uncover a prognostic model and potential oncogenes in colorectal cancer

Abstract

Background Colorectal cancer (CRC) is one of the leading causes of cancer-related death worldwide. Single-cell transcriptome sequencing (scRNA-seq) can provide accurate gene expression data for individual cells. In this study, a new prognostic model was constructed by scRNA-seq and bulk transcriptome sequencing (bulk RNA-seq) data of CRC samples to develop a new understanding of CRC. Methods CRC scRNA-seq data were downloaded from the GSE161277 database, and CRC bulk RNA-seq data were downloaded from the TCGA and GSE17537 databases. The cells were clustered by the FindNeighbors and FindClusters functions in scRNA-seq data. CIBERSORTx was applied to detect the abundance of cell clusters in the bulk RNA-seq expression matrix. WGCNA was performed with the expression profiles to construct the gene coexpression networks of TCGA-CRC. Next, we used a tenfold cross test to construct the model and a nomogram to assess the independence of the model for clinical application. Finally, we examined the expression of the unreported model genes by qPCR and immunohistochemistry. A clone formation assay and orthotopic colorectal tumour model were applied to detect the regulatory roles of unreported model genes. Results A total of 43,851 cells were included after quality control, and 20 cell clusters were classified by the FindCluster () function. We found that the abundances of C1, C2, C4, C5, C15, C16 and C19 were high and the abundances of C7, C10, C11, C13, C14 and C17 were low in CRC tumour tissues. Meanwhile, the results of survival analysis showed that high abundances of C4, C11 and C13 and low abundances of C5 and C14 were associated with better survival. The WGCNA results showed that the red module was most related to the tumour and the C14 cluster, which contains 615 genes. Lasso Cox regression analysis revealed 8 genes (PBXIP1, MPMZ, SCARA3, INA, ILK, MPP2, L1CAM and FLNA), which were chosen to construct a risk model. In the model, the risk score features had the greatest impact on survival prediction, indicating that the 8-gene risk model can better predict prognosis. qPCR and immunohistochemistry analysis showed that the expression levels of MPZ, SCARA3, MPP2 and PBXIP1 were high in CRC tissues. The functional experiment results indicated that MPZ, SCARA3, MPP2 and PBXIP1 could promote the colony formation ability of CRC cells in vitro and tumorigenicity in vivo. Conclusions We constructed a risk model to predict the prognosis of CRC patients based on scRNA-seq and bulk RNA-seq data, which could be used for clinical application. We also identified 4 previously unreported model genes (MPZ, SCARA3, MPP2 and PBXIP1) as novel oncogenes in CRC. These results suggest that this model could potentially be used to evaluate the prognostic risk and provide potential therapeutic targets for CRC patients.

Single cells and TRUST4 reveal immunological features of the HFRS transcriptome

“…In our investigation, we employed the CellChat (version 1.6.1) R package to delve into cell–cell communication intricacies. This tool facilitated the exploration of ligand-receptor interactions on the cell surface, offering insights into intercellular information transmission ( 12 ). Leveraging gene expression data, we deduced protein expression and established a comprehensive cell interaction network.…”

Study of Archaeal Diversity in the Arctic Meltwater Lake Region

“…WGCNA is an R package, and its basic process is to first cluster highly correlated genes and form a module, summarize these clusters with hub genes, correlate the modules with each other and with out-of-sample features and calculate the correlation, and finally find the hub genes [66]. It is widely used in the biomedical field to analyze hub genes for diseases, such as aortic dissection [67], colorectal cancer [68], and endometrial carcinoma [69]. In this study, we used WGCNA to explore the correlations between highly correlated archaeal modules and physicochemical factors to identify significantly influential factors and hub archaea.…”

Identification of driving genes of familial adenomatous polyposis by differential gene expression analysis and weighted gene co-expression network analysis

“…Co-expression networks of filtered genes were constructed through the usage of “WGCNA” R package [ 36 , 37 ]. After assessment of the expression matrix assessed via the average method with the “hclust” function, the clustered genes of gene chips comprising GSM3024204, GSM3024195 and GSM3024196 were defined as biased and were therefore excluded from the further analysis (Fig.…”