- 1Department of Gastrointestinal Surgery, The First Affiliated Hospital, Sun Yat-sen University, Guangzhou, China
- 2Laboratory of General Surgery, The First Affiliated Hospital, Sun Yat-sen University, Guangzhou, China
- 3Center for Digestive Disease, The Seventh Affiliated Hospital of Sun Yat-sen University, Shenzhen, China
Colorectal cancer (CRC) is a leading cause of cancer-related deaths worldwide. Exosomes have great potential as liquid biopsy specimens due to their presence and stability in body fluids. However, the function and diagnostic values of exosomal genes in CRC are poorly understood. In the present study, exosomal data of CRC and healthy samples from the exoRBase 2.0 and Gene Expression Omnibus (GEO) databases were used, and 38 common exosomal genes were identified. Through the least absolute shrinkage and selection operator (Lasso) analysis, support vector machine recursive feature elimination (SVM-RFE) analysis, and logistic regression analysis, a diagnostic model of the training set was constructed based on 6 exosomal genes. The diagnostic model was internally validated in the test and exoRBase 2.0 database and externally validated in the GEO database. In addition, the co-expression analysis was used to cluster co-expression modules, and the enrichment analysis was performed on module genes. Then a protein–protein interaction and competing endogenous RNA network were constructed and 10 hub genes were identified using module genes. In conclusion, the results provided a comprehensive understanding of the functions of exosomal genes in CRC as well as a diagnostic model related to exosomal genes.
Introduction
Colorectal cancer (CRC) is one of the most common malignant tumors of the digestive tract and presents the second-highest mortality rate (9.2%) among all cancers (Arnold et al., 2017; Bray et al., 2018). Tumor metastases and invasion are associated with poor prognosis and lead to a low five-year survival rate in CRC patients (Siegel et al., 2018). CRC is typically detected by colonoscopy, measuring carcinoembryonic antigen (CEA) levels, multitarget stool DNA testing, and the septin 9 gene methylation blood test (Ahluwalia et al., 2021). Although colonoscopy is a highly sensitive method, it is invasive and uncomfortable for patients and its accuracy depends on the skill level and experience of the endoscopist (Schreuders et al., 2015). CEA levels have been widely used as tumor markers for the detection of CRC. However, it is still limited in terms of sensitivity and specificity. Therefore, there is an urgent need to identify more effective and less invasive surrogate CRC-specific diagnostic markers for the rapid, noninvasive, and high sensitivity screening of patients.
The use of exosomes as a new noninvasive approach for diagnosing diseases has attracted growing attention (Mousavi et al., 2019). Exosomes are extracellular vesicles containing messenger RNAs (mRNA), microRNAs (miRNA), long noncoding RNAs (lncRNA), circular RNAs (circRNA), DNA, lipids, and proteins with sizes between 40 and 150 nm and density between 1.13 and 1.19 g/ml (Yáñez-Mó et al., 2015; Shi et al., 2021). Notably, exosomal biomarkers and molecule information remain relatively stable in most body fluids because they are protected from degradation and external impact (Xiao et al., 2019). Therefore, exosomes have great potential as liquid biopsy specimens for various diseases (Wang et al., 2017; Liu et al., 2019). In particular, cancer cells secrete significantly more exosomes than normal cells (Mousavi et al., 2019; Nabariya et al., 2020). Cancer-derived exosomes likely serve as new circulating biomarkers for the early detection of cancer as they carry cargo reflective of genetic or signaling alterations in the cancer cells of origin (Li et al., 2015; Melo et al., 2015). Therefore, exosomes may be an ideal candidate to act as a biomarker for CRC.
In the present study, we identified several differentially expressed exosomal genes from public databases to understand the underlying molecular changes and biological mechanisms. For the diagnosis of CRC patients, we established a 6–exosomal gene diagnosis model using the least absolute shrinkage and selection operator (LASSO), support vector machine recursive feature elimination (SVM-RFE), and logistic regression analyses. This model was verified using a receiver operating characteristic (ROC) curve in internal and external sets. In addition, the co-expression analysis was used to cluster co-expression modules, and the enrichment analysis was performed on module genes. We also established protein–protein interaction (PPI) networks to investigate hub genes and constructed competing endogenous RNA (ceRNA) networks related to serum exosomal genes in CRC. The purpose of the present study was to obtain further insight into the underlying functions of exosomal genes and to identify any potential diagnostic exosomal genes using the bioinformatics analysis in CRC.
Materials and Methods
Data Source and Identification of Differentially Expressed Exosomal Genes
The flowchart of this study is presented in Supplementary Figure S1. We extracted serum exosomal data relative to mRNA, lncRNA, and circRNA from 35 CRC patients and 118 healthy persons from the exoRBase 2.0 database (Li S. et al., 2018) (http://www.exorbase.org/). We also downloaded serum exosomal data from the Gene Expression Omnibus (GEO) database (http://www.ncbi.nlm.nih.gov/geo/). The GSE100063 and GSE100206 datasets contain exosomal data (mRNA) from 12 CRC patients and 32 normal persons. All CRC and healthy exosomal data were available and included in the study. After data normalization and log base 2 transformations, differentially expressed genes between CRC patients and healthy individuals were identified using the limma R package (Ritchie et al., 2015). Differentially expressed genes were defined as those whose expression differences were associated with an adjusted p <0.05. This study was approved by the Ethics Committee of the First Affiliated Hospital of Sun Yat-sen University [Approval number (2021)687].
Identification of the Diagnosis-Related Exosomal Gene Signature Associated With Colorectal Cancer
Genes that were differentially expressed (adjusted p <0.05) in the exoRBase 2.0, GSE100063, and GSE100206 datasets were selected. First, 70% of the samples in exoRBase 2.0 were randomly selected as a training set. The LASSO regression analysis was used to obtain diagnostic exosomal genes from the training set. The SVM-RFE analysis was used simultaneously to screen exosomal genes in the training set for CRC diagnosis. Then, we combined the LASSO and SVM-RFE analyses to obtain the common exosomal genes. Next, we selected common exosomal genes that were regulated (up- or downregulated) in the same direction to build a diagnosis-related exosomal gene signature by the multivariate logistic regression analysis. Finally, the ROC curve analysis and the area under the curve (AUC) were used to estimate the diagnostic value of the diagnosis-related exosomal gene signature using the pROC package in R (Robin et al., 2011).
Validation of the Diagnosis-Related Exosomal Gene Signature
As a validation set, 30% of the samples in exoRBase 2.0, GSE100063, and GSE100206 datasets were selected. To validate whether candidate exosomal genes had an important diagnostic value in patients with CRC, we also measured the ROC curve and the AUC value in the validation datasets. p <0.05 is considered statistically significant.
Cell Culture, Human Plasma Samples, and Isolation of Exosomes
HCoEpiC, HCT116, and SW480 cell lines (purchased from ATCC) were cultured in a DMEM medium (Cellmax) containing 10% fetal calf serum, and 100 U/ml each of penicillin and streptomycin at 37°C with 5% CO2. HCoEpiC, HCT116, and SW480 cells were cultured with a full medium at 80% confluency and replaced with a fresh medium without fetal bovine serum. After a 48-h culture, the cell medium was harvested. In addition, thirty-two CRC patient serums and seventeen healthy human serums were collected in the First Affiliated Hospital of Sun Yat-sen University in December 2021. The samples were stored at −80 before exosome extraction. The exosomes were collected from the cell culture supernatant by differential centrifugations. In addition, exosome morphology was identified by transmission electron microscopy (TEM), the nanoparticle tracking analysis (NTA), and the expression of exosome surface markers CD9, TSG101, and HSP70 were evaluated by the Western blotting analysis.
Validation of Exosomal Gene Expression Levels
Total RNA from exosomes was extracted using the TRIzol reagent (Invitrogen, NYC, United States) and reverse transcribed with the PrimeScript RT kit (Takara, China). Real-time PCR was carried out using the SYBR PreMix Ex Taq II kit (Takara, China). GAPDH was used as the normalized control of mRNA. The relative expression levels of mRNA in exosomes were calculated by the 2−ΔΔCT method. The primer sequences used in this study are listed in Supplementary Table S1.
The Co-Expression Analysis
The weighted gene co-expression network analysis (WGCNA) was used to identify the co-expression network of differentially expressed exosomal genes in exoRBase 2.0 using the WGCNA package in R (Langfelder and Horvath, 2008). A weighted adjacency was constructed by calculating using the Pearson correlations of all gene pairs. Soft power β = 7 was selected to construct a standard scale-free network. The similarity matrix, which was constructed using the Pearson’s correlation coefficients of all gene pairs, was transformed into a topological overlap matrix (TOM) as well as the corresponding dissimilarity (1-TOM). Then, a hierarchical clustering dendrogram of the 1-TOM matrix was used to classify similar gene expressions into different gene co-expression modules. Afterward, a module-clinical trait association was calculated to identify functional modules in a co-expression network. The brown and gray modules were selected for further analysis.
The Gene Ontology Term and Kyoto Encyclopedia of Genes and Genomes Functional Enrichment Analyses
To further clarify the potential biological functions of the exosomal genes in the modules, we performed the Gene Ontology (GO) term and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses using the clusterProfiler package in R (Yu et al., 2012). The GO terms and KEGG pathways with p <0.05 were considered significantly enriched.
Protein–Protein Interaction Network Construction and Hub Gene Identification
In this study, the exosomal genes in the brown and gray modules were analyzed. The PPI network was constructed by the online STRING database (https://string-db.org), and interaction with a combined score >0.7 was considered statistically significant. Cytoscape, an open-source bioinformatics software platform, was used to visualize molecular interaction networks. Hub genes were identified using the Degree algorithm of the cytoHubba plugin in Cytoscape.
Exosomal ceRNA Network Construction
The TargetScan database (www.targetscan.org) and miRanda (http://cbio.mskcc.org/microrna_data/miRanda-aug2010.tar.gz) statistical models that predict the effects of miRNAs binding to canonical sites of mRNA were used to accomplish miRNA prediction. miRNA/circRNA and miRNA/lncRNA interactions were predicted using ENCORI (http://starbase.sysu.edu.cn/) and miRcode (http://www.mircode.org/), respectively. The ceRNA regulatory network was constructed according to the ceRNA regulatory mechanism and the differentially expressed lncRNAs, circRNAs, and mRNAs in exosomes, and the network were visualized by Cytoscape.
Results
Identification of Differential Expression of Exosomal Genes
The RNA sequencing exosomal data for CRC and healthy samples were downloaded from the exoRBase 2.0 and GEO databases. A total of 2839 differentially expressed exosomal genes were obtained from the exoRBase 2.0 dataset (Figure 1A) and 475 differentially expressed exosomal genes in the GSE100063 and GSE100206 datasets (Figure 1B). A total of 38 differentially expressed exosomal genes were common in the two databases (Figure 1C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/863747/fgene-13-863747-HTML/image_m/fgene-13-863747-g001.jpg)
FIGURE 1. Differentially expressed exosomal genes and identification of diagnostic exosomal genes in CRC. (A) Differentially expressed exosomal genes between CRC patients and controls in the exoRBase 2.0 database. (B) Exosomal genes differentially expressed between CRC patients and controls in the GSE100063 and GSE100206 datasets. (C) The intersection of differentially expressed exosomal genes in the exoRBase 2.0, the GSE100063, and the GSE100206 dataset. (D) The LASSO method identified 13 diagnostic exosomal genes. (E) The SVM-RFE method identified 19 diagnostic exosomal genes. (F) The intersection of diagnostic exosomal genes in the two analyses. CRC, colorectal cancer; LASSO, least absolute shrinkage and selection operator; SVM-RFE, support vector machine recursive feature elimination.
Development and Verification of the Diagnostic Model
We randomly split all the samples in the exoRBase 2.0 dataset into a training set (70%) and a validation set (30%). We combined the LASSO and SVM-RFE methods to obtain 9 common exosomal genes (Figures 1D–F). Next, in the training and validation set, we selected common exosomal genes that were in the same regulated direction to build a diagnosis-related exosomal gene signature using the multivariate logistic regression analysis (Figures 2A,B). Therefore, we obtained 6 exosomal genes including H3F3A, MYL6, FBXO7, TUBA1C, MEF2C, and BANK1. Furthermore, these 6 exosomal genes were identified based on the model according to the following formula: index = H3F3A × (3.67971578945741) + MYL6 × (2.01995321634272) + FBXO7 × (−0.86536340517915) + TUBA1C × (0.184108184598825) + MEF2C × (−3.86524121779742) + BANK1 × (−5.35388054514927). The AUC for the gene signature was 0.981 in the training set (Figure 2C) and 0.923 in the internal test set (Figure 2D). Importantly, we also used external datasets to validate the gene signature: it showed good diagnostic ability, with an AUC of 0.995 (Figure 2E).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/863747/fgene-13-863747-HTML/image_m/fgene-13-863747-g002.jpg)
FIGURE 2. Potential exosomal genes for the diagnosis of CRC. The relative expression level of diagnostic exosomal genes (A) between CRC serum and healthy samples in the training set and (B) in the external validation dataset. ROC curves of the exosomal gene signature in the (C) training set and (D) in the internal validation set of the exoRBase 2.0 database. (E) ROC curves of the exosomal gene signature in the external validation set of the GSE100063 and GSE100206 database. CRC, colorectal cancer; ROC, receiver operating characteristic.
Validation of Exosomal Gene Expression Levels
To further validate the expression levels of these model exosomal genes, we extracted exosomes from HCoEpiC, HCT116, SW480 cells, and human plasmas. The purification of the exosomes was validated by TEM, NTA, and Western blotting analysis. TEM detected double-layer spherical vesicles ranging from 30 to 160 nm in size, which confirmed the presence of exosomes (Supplementary Figure S2A). NTA characterized the size and concentration of exosomes (Supplementary Figure S2B). Western blotting analysis estimated the quantity and purity of exosomes by detecting exosomal marker proteins (CD9, TSG101, and HSP70) (Supplementary Figure S2C). Next, we validated the expression of signature exosomal genes in the exosomes obtained by RT-qPCR. The results showed that H3F3A, MYL6, and TUBA1C were significantly upregulated in CRC cells, and MEF2C and FBXO7 were significantly downregulated in CRC cells. However, no significant differences were observed in BANK1 expression between HCT116 and SW480 cells and HCoEpic cells (Supplementary Figure S2D). Furthermore, we validated the expression levels of the model exosomal genes in the serums from the CRC patients and healthy humans. The results showed that the mRNA levels of three genes (H3F3A, TUBA1C, and MYL6) were significantly elevated in the CRC exosomes, whereas BANK1, MEF2C, and FBXO7 were downregulated in the CRC exosomes when compared with those in the healthy human exosomes (Figure 3). These results were consistent with the results of the exoRBase 2.0 and GEO databases.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/863747/fgene-13-863747-HTML/image_m/fgene-13-863747-g003.jpg)
FIGURE 3. Validation of exosomal gene expression levels in CRC patient serums and controls. ∗∗∗∗p < 0.0001; ∗∗p < 0.01; ∗p < 0.05; CRC, colorectal cancer.
The Weighted Gene Co-Expression Network Analysis and Key Module Identification
To identify exosomal genes associated with CRC, we analyzed differentially expressed exosomal genes between CRC and healthy samples in the exoRBase 2.0 database using WGCNA. The included samples were clustered with the hierarchical average linkage clustering method (Figure 4A). The optimal soft power threshold for WGCNA was set to 7 to preserve the scale-free topology and effective connectivity (Figure 4B). Four co-expression modules of differentially expressed exosomal genes were established (Figure 4C). The brown and gray modules were found to have a positive correlation with CRC.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/863747/fgene-13-863747-HTML/image_m/fgene-13-863747-g004.jpg)
FIGURE 4. Identifying the functions of CRC-associated exosomal genes. (A) Clustering dendrogram. (B) Determination of soft-thresholding power in the weighted gene co-expression network analysis. (C) Module–trait associations evaluated by correlations between CRC and clinical traits. (D) The GO enrichment analysis of the exosomal genes in the brown and gray modules. (E) The KEGG pathway of the exosomal genes in the brown and gray modules. (F) The PPI network of genes in the brown and gray modules and hub genes screening. CRC, colorectal cancer; PPI, protein–protein interaction
Enrichment of Module Genes
To explore the potential function of exosomal genes in the brown and gray modules, the enrichment analysis of GO and KEGG was performed. The GO analysis of these differentially expressed exosomal genes revealed that “small molecule catabolic process,” “contractile fiber,” and “structural constituent of muscle” were the most frequent biological terms for biological process, cellular components, and molecular functions, respectively (Figure 4D). The KEGG analysis revealed that these exosomal genes were mainly enriched in the “apelin signaling pathway,” “ubiquitin mediated proteolysis,” “spinocerebellar ataxia,” and “tight junction” (Figure 4E).
Protein–Protein Interaction Network Construction and Hub Genes Screening
The PPI network of the exosomal genes in the brown and gray modules was constructed through the STRING database and visualized with Cytoscape. The PPI network and hub genes identified from the network were obtained through the degree algorithm of the CytoHubba plugin. According to degree scores, the top-scoring genes, including HIST1H3E, HIST1H3J, HIST1H3A, HIST1H2BC, SEH1L, H3F3A, HIST1H2BJ, HIST1H2BF, HIST1H2BB, and RHEB, were considered the hub genes (Figure 4F).
Construction of the Exosomal ceRNA Network
After the differentially expressed exosomal genes were identified in the brown and gray co-expression modules, these exosomal genes were used to construct ceRNA networks. Exosomal lncRNA and circRNA that were differentially expressed (adjusted p <0.05) in exoRBase 2.0 were selected. Based on mRNA, lncRNA, circRNA, and the predicted corresponding miRNA, we constructed a ceRNA network. The ceRNA network consisted of 5 circRNA nodes, 2 lncRNA nodes, 40 miRNA nodes, and 72 mRNA nodes (Figure 5).
Discussion
CRC is a highly malignant cancer with a poor prognosis. CRC patients are at substantial risk of recurrence and metastasis. Therefore, early diagnosis is very important to improve the clinical prognosis of CRC patients. Liquid biopsy, a recent and hot topic in cancer detection, has been considered for the early diagnosis of cancer (Chang et al., 2019). The functional states of cancer cells could be assessed by testing the exosome they secreted, which provides the basis for exosome-mediated noninvasive cancer liquid biopsy (Yang et al., 2020). Currently, various studies show that the increased or decreased expression of exosomes plays an important role in different kinds of cancer, including CRC (Hon et al., 2017; Galamb et al., 2019; Xiao et al., 2020). Thus, exosomes are crucial potential candidates for the early detection of CRC.
To date, exosomal gene–based diagnostic models have not been described for CRC, although considerable efforts have been made to develop prognostic signatures based on differentially expressed genes (Chen et al., 2019; Dai et al., 2020; Wang et al., 2020). The present study analyzed differences in exosomal gene expression between CRC patients and healthy humans. Importantly, we identified diagnostic exosomal genes based on a comprehensive analysis that could serve as valuable biomarkers in the clinical setting.
Potential exosomal gene modules related to CRC were identified with the WGCNA analysis. The brown and gray modules were found to have a positive correlation with CRC. To better understand the potential function of exosomal genes among the brown and gray modules, the GO and KEGG enrichment analysis and the PPI network were conducted. The results of the functional and pathway enrichment analyses showed that the exosomal genes in the modules were mainly enriched in the structural components of the muscle and apelin signaling pathway. The PPI network of exosomal genes in the brown and gray modules was constructed and 10 hub genes were selected using the CytoHubba plugin in Cytoscape. Among these genes, H3F3A, as a diagnostic biomarker in the present study, was also identified as a hub gene, and maybe as a promoter of CRC progression and metastasis. The ceRNA network played a critical role in the initiation and progression of CRC (Ke et al., 2019; Ma et al., 2020). In the present study, we further constructed ceRNA networks based on those key exosomal genes. This approach provided a novel view of the RNA–RNA crosstalk in the exosome and indicated the potential diagnostic and therapeutic functions of exosomal ceRNA networks in CRC.
Our exosomal gene–based model highlighted 6 exosomal genes, that is, H3F3A, MYL6, FBXO7, TUBA1C, MEF2C, and BANK1. These genes and their biological functions have been studied in some tumors. H3F3A is one of two genes encoding histone H3.3, a noncanonical histone variant, and has been established as a major driver gene of malignant gliomas (Schwartzentruber et al., 2012; Sturm et al., 2012; Wu et al., 2012). MYL6, a gene that encodes hexameric ATPase cellular motor protein, is upregulated in circulating tumor cells of many cancers (Yadavalli et al., 2017). FBXO7 may have a proto-oncogenic role in epithelial tumors (Laman et al., 2005). TUBA1C, as a type of tubulin, is associated with tumor cell death and cell proliferation, and TUBA1C overexpression is predicted to have a poor prognosis (Li C.-W. et al., 2018; Li et al., 2020). MEF2C was traditionally considered a development-associated factor and can inhibit tumor growth in vitro and in vivo (Bai et al., 2015). BANK1, which encodes a protein adaptor that is predominantly expressed in B cells, is a putative tumor suppressor gene in B-cell lymphomagenesis (Yan et al., 2014). Our results showed that diagnostic exosomal genes may participate in the development of human CRC; however, the underlying molecular mechanism of these genes in the prognosis of CRC requires further investigation. Our experimental results showed that the differential expression levels of diagnostic exosomal genes in cell lines were approximately in agreement with the serum exosomal data in these public databases.
To the best of our knowledge, this is the first reported exosomal gene–based model for CRC. However, the potential limitations of this study should also be considered when interpreting the findings. First, we could not explore the association between these exosomal genes and the prognosis of CRC patients due to the lack of therapeutic and prognostic information. Second, because the data we analyzed were obtained from public databases, further experimental studies are necessary to validate our findings.
To conclude, we investigated the potential functions and diagnostic values of exosomal genes in CRC through a comprehensive bioinformatics analysis. A diagnostic exosomal gene model was constructed. This model can assess the value of exosomal genes to diagnose CRC using a noninvasive method and might be useful for the development of individualized treatment for CRC patients, but the feasibility of its use in the population needs to be further validated.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.
Ethics Statement
The studies involving human participants were reviewed and approved by the Ethics Committee of the First Affiliated Hospital of Sun Yat-sen University [Approval number (2021)687]. The ethics committee waived the requirement of written informed consent for participation.
Author Contributions
TL and WS conceived the concept, instructed data analysis, and revised the manuscript. YZ and XW conducted most data analysis, prepared figures and tables, and wrote the manuscript draft. WL and WF helped with some analysis and interpretation of data. All authors contributed to the article and approved the submitted version.
Funding
This project was supported by the National Natural Science Foundation of China (grant number: 81871908). Guangdong Provincial Natural Science Foundation (grant number: 2018A030313715). Guangzhou Science and Technology Plan-General Project (grant number: 201904010036). General projects of Guangdong Provincial Natural Science Foundation (grant number: 82070529).
Conflict 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.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors, and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
The authors would like to thank the exoRBase 2.0 and GEO databases for the availability of the data and thank Dr. ZeBing Song for his help and guidance.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.863747/full#supplementary-material
Supplementary Figure S1 | Overall workflow of the present study.
Supplementary Figure S2 | Identification of exosomes and validation of signature exosomal gene expression levels. (A) Morphologies of exosomes in HCoEpic, HCT116, and SW480 cell lines were observed by transmission electron microscopy. (B) The size and concentration of exosomes by the nanoparticle tracking analysis. (C) Exosomes were analyzed by Western blotting using antibodies against exosomal markers. (D) The relative expression fold of diagnostic exosomal genes was determined. ****P<0.0001; ***P<0.001; **P<0.01.
References
Ahluwalia, P., Kolhe, R., and Gahlay, G. K. (2021). The Clinical Relevance of Gene Expression Based Prognostic Signatures in Colorectal Cancer. Biochim. Biophys. Acta (BBA) - Rev. Cancer 1875 (2), 188513. doi:10.1016/j.bbcan.2021.188513
Arnold, M., Sierra, M. S., Laversanne, M., Soerjomataram, I., Jemal, A., and Bray, F. (2017). Global Patterns and Trends in Colorectal Cancer Incidence and Mortality. Gut 66 (4), 683–691. doi:10.1136/gutjnl-2015-310912
Bai, X. L., Zhang, Q., Ye, L. Y., Liang, F., Sun, X., Chen, Y., et al. (2015). Myocyte Enhancer Factor 2C Regulation of Hepatocellular Carcinoma via Vascular Endothelial Growth Factor and Wnt/β-Catenin Signaling. Oncogene 34 (31), 4089–4097. doi:10.1038/onc.2014.337
Bray, F., Ferlay, J., Soerjomataram, I., Siegel, R. L., Torre, L. A., and Jemal, A. (2018). Global Cancer Statistics 2018: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA A Cancer J. Clin. 68 (6), 394–424. doi:10.3322/caac.21492
Chang, L., Ni, J., Zhu, Y., Pang, B., Graham, P., Zhang, H., et al. (2019). Liquid Biopsy in Ovarian Cancer: Recent Advances in Circulating Extracellular Vesicle Detection for Early Diagnosis and Monitoring Progression. Theranostics 9 (14), 4130–4140. doi:10.7150/thno.34692
Chen, L., Lu, D., Sun, K., Xu, Y., Hu, P., Li, X., et al. (2019). Identification of Biomarkers Associated with Diagnosis and Prognosis of Colorectal Cancer Patients Based on Integrated Bioinformatics Analysis. Gene 692, 119–125. doi:10.1016/j.gene.2019.01.001
Dai, G. P., Wang, L. P., Wen, Y. Q., Ren, X. Q., and Zuo, S. G. (2020). Identification of Key Genes for Predicting Colorectal Cancer Prognosis by Integrated Bioinformatics Analysis. Oncol. Lett. 19 (1), 388–398. doi:10.3892/ol.2019.11068
Galamb, O., Barták, B. K., Kalmár, A., Nagy, Z. B., Szigeti, K. A., Tulassay, Z., et al. (2019). Diagnostic and Prognostic Potential of Tissue and Circulating Long Non-Coding RNAs in Colorectal Tumors. World J. Gastroenterol. 25 (34), 5026–5048. doi:10.3748/wjg.v25.i34.5026
Hon, K. W., Abu, N., Ab Mutalib, N.-S., and Jamal, R. (2017). Exosomes as Potential Biomarkers and Targeted Therapy in Colorectal Cancer: A Mini-Review. Front. Pharmacol. 8, 583. doi:10.3389/fphar.2017.00583
Ke, M. J., Ji, L. D., and Li, Y. X. (2019). Explore Prognostic Marker of Colorectal Cancer Based on ceRNA Network. J. Cell. Biochem. 120 (12), 19358–19370. doi:10.1002/jcb.28860
Laman, H., Funes, J. M., Ye, H., Henderson, S., Galinanes-Garcia, L., Hara, E., et al. (2005). Transforming Activity of Fbxo7 Is Mediated Specifically Through Regulation of Cyclin D/cdk6. Embo J. 24 (17), 3104–3116. doi:10.1038/sj.emboj.7600775
Langfelder, P., and Horvath, S. (2008). WGCNA: An R Package for Weighted Correlation Network Analysis. BMC Bioinf. 9, 559. doi:10.1186/1471-2105-9-559
Li, Y., Zheng, Q., Bao, C., Li, S., Guo, W., Zhao, J., et al. (2015). Circular RNA Is Enriched and Stable in Exosomes: A Promising Biomarker for Cancer Diagnosis. Cell Res. 25 (8), 981–984. doi:10.1038/cr.2015.82
Li, C.-W., Chiu, Y.-K., and Chen, B.-S. (2018a). Investigating Pathogenic and Hepatocarcinogenic Mechanisms from Normal Liver to HCC by Constructing Genetic and Epigenetic Networks via Big Genetic and Epigenetic Data Mining and Genome-Wide NGS Data Identification. Dis. Markers 2018, 1–22. doi:10.1155/2018/8635329
Li, S., Li, Y., Chen, B., Zhao, J., Yu, S., Tang, Y., et al. (2018b). exoRBase: A Database of circRNA, lncRNA and mRNA in Human Blood Exosomes. Nucleic Acids Res. 46 (D1), D106–d112. doi:10.1093/nar/gkx891
Li, J., Zhou, Y., Yan, Y., Zheng, Z., Hu, Y., and Wu, W. (2020). Sulforaphane-Cysteine Downregulates CDK4/CDK6 and Inhibits Tubulin Polymerization Contributing to Cell Cycle Arrest and Apoptosis in Human Glioblastoma Cells. Aging 12 (17), 16837–16851. doi:10.18632/aging.103537
Liu, H., Li, P.-w., Yang, W.-q., Mi, H., Pan, J.-l., Huang, Y.-c., et al. (2019). Identification of Non-Invasive Biomarkers for Chronic Atrophic Gastritis from Serum Exosomal microRNAs. BMC Cancer 19 (1), 129. doi:10.1186/s12885-019-5328-7
Ma, Z., Han, C., Xia, W., Wang, S., Li, X., Fang, P., et al. (2020). circ5615 Functions as a ceRNA to Promote Colorectal Cancer Progression by Upregulating TNKS. Cell Death Dis. 11 (5), 356. doi:10.1038/s41419-020-2514-0
Melo, S. A., Luecke, L. B., Kahlert, C., Fernandez, A. F., Gammon, S. T., Kaye, J., et al. (2015). Glypican-1 Identifies Cancer Exosomes and Detects Early Pancreatic Cancer. Nature 523 (7559), 177–182. doi:10.1038/nature14581
Mousavi, S., Moallem, R., Hassanian, S. M., Sadeghzade, M., Mardani, R., Ferns, G. A., et al. (2019). Tumor-Derived Exosomes: Potential Biomarkers and Therapeutic Target in the Treatment of Colorectal Cancer. J. Cell. Physiol. 234 (8), 12422–12432. doi:10.1002/jcp.28080
Nabariya, D. K., Pallu, R., and Yenuganti, V. R. (2020). Exosomes: The Protagonists in the Tale of Colorectal Cancer? Biochim. Biophys. Acta (BBA) - Rev. Cancer 1874 (2), 188426. doi:10.1016/j.bbcan.2020.188426
Ritchie, M. E., Phipson, B., Wu, D., Hu, Y., Law, C. W., Shi, W., et al. (2015). Limma Powers Differential Expression Analyses for RNA-Sequencing and Microarray Studies. Nucleic Acids Res. 43 (7), e47. doi:10.1093/nar/gkv007
Robin, X., Turck, N., Hainard, A., Tiberti, N., Lisacek, F., Sanchez, J.-C., et al. (2011). pROC: An Open-Source Package for R and S+ to Analyze and Compare ROC Curves. BMC Bioinf. 12, 77. doi:10.1186/1471-2105-12-77
Schreuders, E. H., Ruco, A., Rabeneck, L., Schoen, R. E., Sung, J. J. Y., Young, G. P., et al. (2015). Colorectal Cancer Screening: A Global Overview of Existing Programmes. Gut 64 (10), 1637–1649. doi:10.1136/gutjnl-2014-309086
Schwartzentruber, J., Korshunov, A., Liu, X.-Y., Jones, D. T. W., Pfaff, E., Jacob, K., et al. (2012). Driver Mutations in Histone H3.3 and Chromatin Remodelling Genes in Paediatric Glioblastoma. Nature 482 (7384), 226–231. doi:10.1038/nature10833
Shi, H., Wang, M., Sun, Y., Yang, D., Xu, W., and Qian, H. (2021). Exosomes: Emerging Cell-Free Based Therapeutics in Dermatologic Diseases. Front. Cell Dev. Biol. 9, 736022. doi:10.3389/fcell.2021.736022
Siegel, R. L., Miller, K. D., and Jemal, A. (2018). Cancer Statistics, 2018. CA A Cancer J. Clin. 68 (1), 7–30. doi:10.3322/caac.21442
Sturm, D., Witt, H., Hovestadt, V., Khuong-Quang, D.-A., Jones, D. T. W., Konermann, C., et al. (2012). Hotspot Mutations in H3F3A and IDH1 Define Distinct Epigenetic and Biological Subgroups of Glioblastoma. Cancer Cell 22 (4), 425–437. doi:10.1016/j.ccr.2012.08.024
Wang, J., Yan, F., Zhao, Q., Zhan, F., Wang, R., Wang, L., et al. (2017). Circulating Exosomal miR-125a-3p as a Novel Biomarker for Early-Stage Colon Cancer. Sci. Rep. 7 (1), 4150. doi:10.1038/s41598-017-04386-1
Wang, Y., Chen, Y., Xiao, S., and Fu, K. (2020). Integrated Analysis of the Functions and Prognostic Values of RNA-Binding Proteins in Colorectal Cancer. Front. Cell Dev. Biol. 8, 595605. doi:10.3389/fcell.2020.595605
Wu, G., Broniscer, A., McEachron, T. A., Lu, C., Paugh, B. S., Becksfort, J., et al. (2012). Somatic Histone H3 Alterations in Pediatric Diffuse Intrinsic Pontine Gliomas and Non-Brainstem Glioblastomas. Nat. Genet. 44 (3), 251–253. doi:10.1038/ng.1102
Xiao, Y., Li, Y., Yuan, Y., Liu, B., Pan, S., Liu, Q., et al. (2019). The Potential of Exosomes Derived from Colorectal Cancer as a Biomarker. Clin. Chim. Acta 490, 186–193. doi:10.1016/j.cca.2018.09.007
Xiao, Y., Zhong, J., Zhong, B., Huang, J., Jiang, L., Jiang, Y., et al. (2020). Exosomes as Potential Sources of Biomarkers in Colorectal Cancer. Cancer Lett. 476, 13–22. doi:10.1016/j.canlet.2020.01.033
Yadavalli, S., Jayaram, S., Manda, S. S., Madugundu, A. K., Nayakanti, D. S., Tan, T. Z., et al. (2017). Data-Driven Discovery of Extravasation Pathway in Circulating Tumor Cells. Sci. Rep. 7, 43710. doi:10.1038/srep43710
Yan, J., Nie, K., Mathew, S., Tam, Y., Cheng, S., Knowles, D. M., et al. (2014). Inactivation of BANK1 in a Novel IGH-Associated Translocation T(4;14)(q24;q32) Suggests a Tumor Suppressor Role in B-Cell Lymphoma. Blood Cancer J. 4 (5), e215. doi:10.1038/bcj.2014.36
Yáñez-Mó, M., Siljander, P. R.-M., Andreu, Z., Bedina Zavec, A., Borràs, F. E., Buzas, E. I., et al. (2015). Biological Properties of Extracellular Vesicles and Their Physiological Functions. J. Extracell. Vesicles 4, 27066. doi:10.3402/jev.v4.27066
Yang, D., Zhang, W., Zhang, H., Zhang, F., Chen, L., Ma, L., et al. (2020). Progress, Opportunity, and Perspective on Exosome Isolation - Efforts for Efficient Exosome-Based Theranostics. Theranostics 10 (8), 3684–3707. doi:10.7150/thno.41580
Keywords: exosome, diagnostic model, functions, colorectal cancer, bioinformatics analysis
Citation: Lei T, Zhang Y, Wang X, Liu W, Feng W and Song W (2022) A Diagnostic Model Using Exosomal Genes for Colorectal Cancer. Front. Genet. 13:863747. doi: 10.3389/fgene.2022.863747
Received: 27 January 2022; Accepted: 19 May 2022;
Published: 15 July 2022.
Edited by:
Liang Zhao, University of São Paulo, BrazilReviewed by:
Ehsan Nazemalhosseini-Mojarad, Shahid Beheshti University of Medical Sciences, IranLaura La Paglia, Institute for High Performance Computing and Networking (ICAR), Italy
Copyright © 2022 Lei, Zhang, Wang, Liu, Feng and Song. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Wu Song, songwu@mail.sysu.edu.cn
†These authors have contributed equally to this work