Plasma Extracellular Vesicle Long RNAs Have Potential as Biomarkers in Early Detection of Colorectal Cancer

Background: Early detection of colorectal cancer (CRC) is crucial to the treatment and prognosis of patients. Traditional screening methods have disadvantages.

Methods: 231 blood samples were collected from 86 CRC, 56 colorectal adenoma (CRA), and 89 healthy individuals, from which extracellular vesicle long RNAs (exLRs) were isolated and sequenced. An CRC diagnostic signature (d-signature) was established, and prognosis-associated cell components were evaluated.

Results: The exLR d-signature for CRC was established based on 17 of the differentially expressed exLRs. The d-signature showed high diagnostic efficiency of CRC and control (CRA and healthy) samples with an area under the curve (AUC) of 0.938 in the training cohort, 0.943 in the validation cohort, and 0.947 in an independent cohort. The d-signature could effectively differentiate early-stage (stage I–II) CRC from healthy individuals (AUC 0.990), as well as differentiating CEA-negative CRC from healthy individuals (AUC 0.988). A CRA d-signature was also generated and could differentiate CRA from healthy individuals both in the training (AUC 0.993) and validation (AUC 0.978) cohorts. The enrichment of class-switched memory B-cells, B-cells, naive B-cells, and mast cells showed increasing trends between CRC, CRA, and healthy cohorts. Class-switched memory B-cells, mast cells, and basophils were positively associated with CRC prognosis while natural killer T-cells, naive B-cells, immature dendritic cells, and lymphatic endothelial cells were negatively associated with prognosis.

Conclusions: Our study identified that the exLR d-signature could differentiate CRC from CRA and healthy individuals with high efficiency and exLR profiling also has potential in CRA screening and CRC prognosis prediction.

Introduction

Colorectal cancer (CRC) ranks the third common cancer in men and the second in women, as well as the second cause of cancer death worldwide, which remains an enormous socioeconomic burden on society (1, 2). Meanwhile, colorectal adenoma (CRA) usually take years to develop to invasive or metastatic CRC, which makes CRC one of the cancers most suitable for early detection (3).

Early detection of CRC is the key to reducing invasive treatment, morbidity, mortality, and treatment cost (3). CRC screening methods include invasive and non-invasive tests. Colonoscopy is widely known as the golden standard but limited by invasiveness and low compliance rate (4). The guaiac fecal occult blood test (gFOBT) and fecal immunochemical test for hemoglobin (FIT) are most widely used because they are convenient, cheap, and non-invasive. However, these fecal tests have limitations of low sensitivity or specificity (3). CT colonography, anther non-invasive test, is costly and not sensitive to tumors less than 10 mm (3, 5). From the above, blood tests tend to be more acceptable for CRC screening, but no reliable detecting method or markers have been widely acknowledged (6).

Extracellular vesicles (EVs) are extracellular membrane vesicles originated and released from endocytosis and exocytosis, containing proteins, DNA, RNA, and lipids (7). Due to the protection of the lipid membrane, EV RNAs are likely to be more stable than other free plasma RNA. Long RNAs have been identified in human blood EVs, including messenger RNA (mRNA), long non-coding RNA (lncRNA), and circular RNA (circRNA), which have emerged as promising markers for cancer diagnosis recently and have already been evaluated in some cancers (8–10). However, difficulties in EV research lie on the lack of efficient and stable methods for plasma EVs isolating and purifying. Fortunately, an optimized strategy for plasma EV long RNA (exLR) sequencing (exLR-seq) has been developed and reliable positive data have been obtained in our recent studies (11, 12).

In this study, a CRC diagnostic signature (d-signature) based on plasma exLR profiling was identified and validated, which could differentiate CRC from control (CRA and healthy) individuals efficiently. We also evaluated cell components and signaling pathways between CRC, CRA, and healthy groups, and associated prognostic significance were revealed.

Patients and Methods

Patients

From February 2018 to January 2019, 194 blood samples were collected from 72 CRC patients, 42 CRA patients, and 80 age- and sex-matched healthy participants receiving routine medical examination. The diagnoses of all CRC and CRA patients were pathologically confirmed, and these participants did not have a history of other malignant tumors. All enrolled CRC patients underwent surgical treatment without preoperative chemotherapy or radiotherapy at the Colorectal Surgery Department of Fudan University Shanghai Cancer Center. 37 blood samples (14 CRC, 14 CRA, 9 healthy) were collected in an independent center from Fujian Medical University Union Hospital.

EVs Identification and exLR-seq Analysis

The optimized strategy for plasma exLR-seq included several steps as follows: plasma sample collection, EV purification, transmission electron microscopy (TEM), size distribution measurement, RNA isolation, and RNA-seq library preparation (11). To be brief, the blood samples of CRC and CRA patients were collected before the excision of tumor and centrifuged twice at 3,000 and 13,000 rpm, respectively. The EV RNAs were isolated using the exoRNeasy Serum/Plasma Kit, and the EVs were photographed using a TEM. The size distribution was analyzed using Flow NanoAnalyzer. EV markers TSG101 and CD63 were estimated by Western blots. The RNA-seq library was prepared using SMART technology and sequenced by the Illumina sequencing platform. Details of these steps are found in Supplementary Materials.

ExLR-Seq Analysis for Quantifying Gene Expression

The qualified FASTQ files generated from RNA-seq were aligned to the human genome (hg38) using STAR v2.5.3 with default parameters (13). The mapped sequencing reads in the resulting BAM files were then assigned to genes by featureCounts v1.6.3 (14). Considering that the transcriptome library was reversely stranded, “-s” was set as 2 for performing strand-specific read counting. Genes were annotated with GENCODE v.29. The read count of each gene was converted to transcripts per million (TPM) as follows:

$$TP{M}_i=\frac{R{C}_i}{L_i}\ast \left(\frac{1}{{\sum }_{j=1}^{LR}\frac{R{C}_j}{L_j}}\right)\ast {10}^6$$

Where RC_i stands for the count of reads mapped to the gene and L_i is the length of the gene. LR is the number of long RNA genes including protein coding and long non-coding genes.

Differential Expression Analysis and Pathway Enrichment Analysis

We calculated the correlation coefficient between each two samples based on TPM expression profiles and filtered poor samples with the median of correlation coefficients smaller than 0.9. The final dataset analyzed in our study contained 72 CRC samples and 122 control (42 CRA and 80 healthy) samples. To explore differentially expressed genes (DEGs) between these two cohorts, we applied R package “limma” on TPM expression profiles (15). The Benjamini–Hochberg approach was used to adjust the p values for multiple testing. A gene with a fold change (FC) bigger than 1.5 and adjusted p value smaller than 0.05 was defined as a DEG. To investigate the differential pathways between CRC and control samples, R package “clusterProfiler” was used for KEGG pathway enrichment analysis based on the DEGs (16).

Selecting Effective Feature Genes and Building CRC/CRA-Identification Model

The whole dataset was randomly divided into training cohort (48 CRC and 82 control) and validation cohort (24 CRC and 40 control). With respect to the training cohort, we firstly conducted DEG analysis. To elect informative and functional signature genes for effectively distinguishing CRC samples from control samples, we focused on these upregulated protein coding or long non-coding genes in CRC samples. Then, we employed the minimum redundancy maximum relevance (mRMR) algorithm to rank these candidate genes. This was implemented using the mRMR package with the “MIQ” feature selection scheme (http://home.penglab.com/proj/mRMR/) (17). Next, we applied the incremental feature selection (IFS) strategy to determine the optimal subset of feature genes based on the support vector machine (SVM) (18). The first feature set was constructed with the top one gene. The remaining ranked feature genes were added one by one incrementally for producing new feature sets. Each new feature set was composed of the previous set adding with a new feature gene. Each feature gene set was evaluated with the area under the curve (AUC) value derived from the SVM model using leave one out cross-validation (LOOCV). Finally, the optimal CRC-identification model was built using the feature gene set with the highest AUC value. This model was then applied to classify the validation cohort for further assessing the prediction performance of these feature genes. SVM models were constructed using the LibSVM software package downloaded from https://www.csie.ntu.edu.tw/~cjlin/libsvm/ (19). The CRA-identification model was built in the same way.

Cell Type and Pathway Estimation

To infer the cell types of EV origins, we performed xCell analysis on TPM expression profiles using R package “xCell,” a gene signature-based method that integrates the advantages of gene set enrichment with deconvolution approaches (20). We obtained the enrichment scores of 64 immune and stromal cell types and further investigated the influence of each cell type on the overall survival (OS) and disease-free survival (DFS) of CRC samples. The survival analysis and Kaplan–Meier plotting were implemented by R package “survminer.” The single sample gene set enrichment analysis (ssGSEA) algorithm was used to calculate the enrichment scores of the canonical MSigDB pathways (C2, KEGG) (21). This was carried out on R package “GSVA” with the method of “ssGSEA” (22). To explore the significant different cell types and pathways among CRC, CRA, and normal cohorts, the Wilcoxon-rank sum test was used for comparison between any two cohorts and the one-way analysis of variance (ANOVA) test was used for comparisons among the three cohorts.

Results

Patient Characteristics

In general, 194 participants were involved in our center, consisting of 72 CRC patients, 42 CRA patients, and 80 healthy individuals. The clinicopathological information is listed in Table 1. No obvious difference was seen in age, gender, or tumor site between the three groups. We included more early-stage CRC (stage I–II, 53 cases) than advanced CRC (stage III–IV, 19 cases) because this study was designed to mainly focus on the early detection of CRC. All the CRC patients were followed up for at least 24 months. Death events were observed in 13 stage IV CRC patients, and tumor recurrence or metastasis events were observed in 8 stage II/III CRC patients.

TABLE 1

Table 1 Clinicopathological information of 194 participants.

EVs Isolation and exLR-seq

The isolated EVs observed by TEM were round capsule bubbles. The scanning electron microscope images of EVs are shown in Figure 1A. Since types of EVs (exosomes, microvesicles, and apoptotic bodies) should be distinguished by diameter, we analyzed the size distribution by flow cytometry (10). The size distribution result revealed abundant peaks ranging from 50 to 200 nm and a mean diameter of 103.9 ± 38.6 nm (Figure 1B), indicating that morphologically most of the isolated EVs were exosomes with definition of 40–200 nm in diameter (10). Western blot analysis confirmed that the EV markers CD63 and TSG101 were enriched in EVs but not peripheral blood mononuclear cells (PMBCs), while the negative-control protein marker calnexin was enriched in PMBCs but not EVs (Figure 1C). Afterward, exLR-seq was conducted and no obvious difference of detected mRNA, lncRNA, and pseudogene amount was observed between the three groups (Figure 1D). Unsupervised hierarchical clustering revealed clear separations of CRC and control (CRA and healthy) samples, as well as CRC, CRA, and healthy samples (Figure 1E). The differentially expressed exLRs were enriched for some cancer-associated pathways, such as transcriptional misregulation in cancer and NF-kappa B signaling pathway (Figure 1F). Therefore, we hypothesized that exLRs have potential as diagnostic biomarkers of CRC.

FIGURE 1

Figure 1 Plasma EVs and exLR-seq. (A) Photograph of EVs using a TEM. (B) Size distribution of EVs. (C) Western blot analysis of EV markers TSG101 and CD63 in PMBC and EVs. (D) Amount of exLRs for each sample among CRC, CRA, and healthy individuals. (E) Unsupervised hierarchical clustering of the exLRs differentially expressed between CRC and control (class I); CRC, healthy, and CRA (class II). (F) KEGG pathway enrichment analysis for differentially expressed exLRs.

Establishment of an exLR d-Signature for CRC

To identify the diagnostic potential of exLRs, we developed an exLR-based d-signature for CRC. The flowchart of the establishment of the d-signature is presented in Figure 2A. By random sampling, the cohort was divided into a training cohort (48 CRC, 82 control) and a validation cohort (24 CRC, 40 control). Next, we selected 66 long RNA genes upregulated in CRC samples compared with control samples by DEG analysis (expression frequency >0.5, log₂(FC) >0.59 and adjusted p value < 0.05). MRMR and SVM were used to select the optimal feature gene set among the training cohort. The top 17 genes of the ranked 66 genes were selected to build the SVM prediction model (Table 2). Unsupervised hierarchical clustering of the 17 genes showed relatively high consistency between predicting CRC and true CRC individuals in both training and validation cohorts (Figures 2B, C). The d-signature was applied in the training cohort and validation cohort to assess the diagnostic efficiency. We generated receiver operating characteristic (ROC) plots, displaying the performance of the d-signature in the training cohort, the validation cohort, and the independent cohort (Figures 2D–F). The training sensitivity, specificity, and accuracy were 82.93%, 93.75%, and 86.15%, respectively (Figure 2D and Table 3). The validation sensitivity, specificity, and accuracy were 87.50%, 91.67%, and 87.50%, respectively (Figure 2E and Table 3). The sensitivity, specificity, and accuracy of the independent cohort were 71.43%, 95.65%, and 86.49% (Figure 2F and Table 3). The CRC d-signature showed high diagnostic efficiency both in the training cohort and the validation cohort, as well as the independent cohort.

FIGURE 2

Figure 2 Establishment of the exLR d-signature. (A) Flowchart of establishment of the d-signature. (B, C) Unsupervised hierarchical clustering of the 17 genes in training cohort (B) and validation cohort (C). (D–F) ROC curve for the exLR d-signature in the training (D), validation (E), and independent (F) cohorts. ^aSelection of lncRNA or protein-coding genes with (1) expression frequency >0.5; (2) log2(FC) >0.59, adjusted p value < 0.05.

TABLE 2

Table 2 Basic information and expression of the 17 feature genes.

TABLE 3

Table 3 The exLR d-signature in diagnosis of CRC, CRA, and healthy participants.

The exLR d-Signature for Early Detection of CRC

We further evaluated the performance of the exLR d-signature in subgroups. The d-signature could differentiate between healthy, CRA, and CRC cohorts, and an increasing trend of the diagnostic probability was shown among the three cohorts, which is consistent with the process of CRC carcinogenesis (Figure 3A). Performance of the d-signature was then assessed among different stages of the CRC and control cohorts. As shown in Figure 3B, the d-signature had diagnostic ability for CRC of stages I, II, III, and IV. The sensitivity, specificity, and accuracy of the d-signature to differentiate CRC from CRA were 76.19%, 84.72%, and 79.83% (Figure 3C and Table 3). The diagnostic efficiency was higher for the d-signature to differentiate between CRC and healthy cohorts (sensitivity 92.50%, specificity 94.44%, accuracy 89.47%, Figure 3D and Table 3). As for the early-stage (stage I–II) CRC versus CRA subgroup, the sensitivity, specificity, and accuracy were 85.71%, 81.13%, and 82.11% (Figure 3E and Table 3). The sensitivity, specificity, and accuracy for the d-signature to differentiate between early-stage (stage I–II) CRC and healthy cohorts were 95.00%, 96.23%, and 92.48%, respectively (Figure 3F and Table 3).

FIGURE 3

Figure 3 Prediction performance of the exLR d-signature in subgroups. (A) The d-signature in distinguishing healthy, CRA, and CRC individuals. (B) The d-signature in control and stage I–IV CRC participants. The ROC curve for the d-signature in CRC and CRA (C), CRC and healthy (D), early-stage (stage I–II) CRC and CRA (E), early-stage (stage I–II) CRC and healthy (F), CEA-negative CRC and CRA (G), and CEA-negative CRC and healthy (H) cohorts.

Carcinoembryonic antigen (CEA) is one of the most common cancer markers but limited by low diagnostic efficiency when used along for CRC diagnosis (23). The performance of the d-signature in distinguishing CEA-negative CRC from CRA cohorts is presented in Figure 3G and Table 3 (sensitivity 76.19%, specificity 87.81%, accuracy 80.72%). High performance was observed of the d-signature to differentiate CEA-negative CRC from healthy cohorts (sensitivity 92.50%, specificity 97.56%, accuracy 92.56%, Figure 3H and Table 3). The diagnostic ability of the d-signature to differentiate between CRA and CRC, especially early-stage (stage I–II) and CEA-negative CRC, was of great significance to determine whether the tumor should be resected endoscopically or surgically in clinical practice. Meanwhile, the high efficiency of the d-signature to differentiate between healthy and CRC individuals, including early-stage and CEA-negative CRC individuals, was supposed to have an important potential role in CRC screening.

Potential of the exLR d-Signature in Detecting CRA

In addition to the diagnosis of CRC, detection of CRA is also a very important link in the management of CRC, considering CRA as precancerous lesions of CRC. In this part, we developed another exLR-based d-signature for CRA in the same way as building the CRC d-signature. Unsupervised hierarchical clustering revealed a clear separation of CRA and healthy samples (Figure 4A). KEGG analysis showed that the differentially expressed exLRs were enriched for some tumor-associated pathways (Figure 4B). Unsupervised hierarchical clustering of the top 7 genes selected to build the CRA-identification model showed high consistency between predicting CRA and true CRA individuals in both the training and validation cohorts (Figures 4C, D). Encouraging results of the CRA d-signature were observed both in the training (sensitivity 89.29%, specificity 98.15%, accuracy 95.12%) and validation (sensitivity 71.43%, specificity 96.15%, accuracy 87.50%) cohorts (Figures 4E, F and Table 3).

FIGURE 4

Figure 4 Potential of the exLR d-signature in differentiating CRA and healthy participants. (A) Unsupervised hierarchical clustering of the differentially expressed exLRs between CRA and healthy cohorts. (B) KEGG pathway enrichment analysis for the differentially expressed exLRs between CRA and healthy cohorts. (C, D) Unsupervised hierarchical clustering of the 7 genes selected for d-signature establishment in the training cohort (C) and validation cohort (D). (E, F) ROC curve for the exLR d-signature in the training (D) and validation (E) cohorts.

Estimation of Cell Populations and Prognostic Prediction

EVs are produced by many cell types including immune cells, serving as communicators of immune-modulatory activities that affect the tumor microenvironment and antitumor immune responses (24). We used xCell to infer cell populations represented in EVs. Abundances of 64 immune and stromal cell types based on gene expression profile were estimated, and 21 of them showed statistical differences, including epithelial, lymphoid, myeloid, stem, and stroma cells (Figure 5A). Low enrichment of class-switched memory B-cells, B-cells, naive B-cells, and mast cells was observed in the CRC group compared with CRA and healthy groups, and there was a slight increasing trend among CRC, CRA, and healthy cohorts, implying that the tumor-immune microenvironment had been affected in the CRC group (Figure 5B). In the analysis of prognostic significance, a positive correlation was observed between longer OS and the abundance of class-switched memory B-cells and mast cells, while a negative correlation was observed between OS and the abundance of natural killer T-cells (NKT cells) and naive B-cells (Figure 5C). A high basophil level was associated with longer DFS, while a high level of immature dendritic cells and lymphatic endothelial cells predicted shorter DFS (Figure 5D). These prognosis-associated cell populations were supposed to play a role in CRC prognostic prediction. Besides, we assessed the pathway enrichment of differentially expressed transcriptomes between the CRC, CRA, and healthy groups by performing gene set enrichment via KEGG analysis, showing that the differentially expressed exLRs were enriched in the intestinal immune network for the IgA production pathway and the circadian rhythm mammal pathway with a gradual rising trend between the three groups (Figure 5E). These results presented the potential applications of the exLR profiling.

FIGURE 5

Figure 5 Analyses of cell components, survival, and signaling pathways. (A) Heatmap of unsupervised hierarchical clustering of the 21 cell types in different groups. (B) Box plots of selected cell-type abundance between CRC, CRA, and healthy groups. Prognostic significance of selected cell types by (C) OS and (D) DFS. (E) ssGSEA score and statistical significance for selected KEGG pathways differing between CRC, CRA, and healthy groups. ***p < 0.001; **p < 0.01; *p < 0.05; NS, not significant.

Discussion

In this study, exLR-seq expression profiles were gained from 231 CRC, CRA, and healthy blood samples. To our knowledge, this is the first study focusing on the early detection potential of exLRs between CRC, CRA, and healthy individuals. The preliminary findings seem to be inspiring as certain diagnostic and prognosis prediction efficiency was achieved.

Extracellular vesicles, known as small membranous vesicles released by cells, have recently been identified to contain long RNAs, which may serve as biomarkers in the diagnosis, therapeutic sensitivity prediction, and prognostic prediction of tumors (8, 9, 12, 25). Although the clinical application of EVs is still in its infancy, EVs are increasingly recognized as promising biomarkers for tumor diagnosis and prognosis (10). However, previous studies are mainly focused on protein and miRNAs in EVs. In reviewing the literature, no published study was found to in-depth analyze the diagnostic or prognostic value of exLRs in CRC due to the limitation of methodology and size of samples.

Nowadays, the incidence and mortality of colorectal cancer remain high in both developed and developing countries. Early detection is a key to reducing morbidity and the socioeconomic burden. Traditional detection methods, including colonoscopy, gFOBT, FIT, and CT colonography, all have some limitations of invasiveness, high expense, or low efficiency (2, 3). Emerging screening strategies, such as ctDNA, circulating tumor cells, and septin-9, have been studied widely. Nonetheless, results in relevant studies have shown much lower diagnostic efficiency of CRA and early-stage CRC than that of advanced-stage CRC (6, 26).

A diagnostic signature based on plasma exLR profiling was developed in this study. We first verified EVs from TEM morphology, size distribution analysis, and Western blot analysis. These all corresponded to the characteristics of EVs (7). ExLR profiling of plasma samples from 194 participants was successfully performed using an optimized exLR-seq strategy we recently developed (11). We established a d-signature of 17 exLRs for CRC detection, which could efficiently differentiate CRC from control (CRA and healthy) cohorts (training AUC = 0.938, validation AUC = 0.943, independent cohort AUC = 0.947). In clinical practice, people with positive testing results are supposed to take colonoscopy examination to identify the results. The d-signature makes it possible to screen high-risk patients efficiently and reliably, standing a good chance of easing the suffering of the screened people and improving screening compliance.

High sensitivity and specificity were identified for the d-signature to differentiate CRC from CRA, which was of great significance in clinical practice, especially when it comes to early-stage (stage I–II) CRC or CEA-negative CRC. In clinical practice, CRA patients need no additional surgery if the polyp has been completely endoscopically resected with favorable histologic features, while radical surgery plays a vital role in the treatment of most early-stage CRC patients (27, 28). Different diagnoses of CRC or CRA lead to different treatment strategies, and this d-signature is supposed to provide reference for clinicians and patients to make decisions. Compared with differentiating between CRA and CRC cohorts, the d-signature had higher diagnostic efficiency to differentiate between healthy and CRC cohorts, including early-stage (stage I–II) CRC and CEA-negative CRC. This is of great significance for improving the efficiency of CRC screening, considering the limitations of traditional non-invasive CRC screening methods (3, 5).

The 17 genes used to establish the d-signature comprised 16 protein-coding genes and one lncRNA gene, all of which were upregulated in CRC samples. The H2BFS expression level in lung cancer tissue has been reported to be higher than that in normal lung tissue (29). However, its expression in CRC remains unknown. In a previous study, a high expression level of XCL2 was revealed to be associated with NK cells in tumor-immune activities (30). DMC1, short for “downregulated in multiple cancers-1,” plays an important role in DNA binding and repairing, with loss expression identified in multiple human cancers (31). The different expression levels in this study might be explained by using peripheral blood samples but not tumor tissue samples. KLHDC8B is suggested to have a role in the formation of Hodgkin/Reed–Sternberg cells in familial Hodgkin lymphoma (32). CA3 expression is reported to promote the transformation and invasive ability of hepatocellular carcinoma cells (33). Overexpressed CYP20A1 is observed in some pathological types of lung cancer and associated with prognosis according to a previous study (34). The expression of HIST1H2BB is reduced in ovarian cancer cells and might have growth-suppressing roles (35). STK3 is a critical molecule of the Hippo pathway that controls cell development, proliferation, and apoptosis (36). The expression level of CBWD1 has been reported to be associated with melanoma (37). The tumor-associated significance of the other seven genes (HIST2H2AA4, UQCRHL, AC008269.1, RAB6D, APOL4, HIST1H2AI, ANKAR, SGMS1) remains unclear.

This study was mainly designed to build a d-signature for CRC screening, and we were surprised to find that a similar model might be very efficient in CRA diagnosis. However, due to the limitation of CRA cohort size, we believe that the encouraging initial results need to be reconfirmed in further study with larger cohorts.

In this study, statistical differences of 21 immune cell types estimated based on the gene expression profile were observed between CRC, CRA, and healthy cohorts. Actually, the relationship between systemic immune cells and CRC still remains poorly understood, even though some studies with a small sample size have yielded some preliminary conclusions (38, 39). In this study, differences in immune cell subset distribution were observed between CRC, CRA, and healthy cohorts, such as reduced percentage of class-switched memory B-cells, B-cells, naive B-cells, and mast cells in the CRC cohort. This study also showed correlations between survival and these cells. A decreased percentage of peripheral blood B-cells and naive-B cells in the CRC cohort compared with the healthy cohort has been reported previously, whereas the percentage of peripheral blood memory B cells was increased in the CRC cohort in that study (39). Contrary prognostic implications of class-switched memory B-cells and naive B-cells were revealed in this study, and both the tumor progression-enhancing and -suppressing effects of B-cells have been reported in previous literature (40, 41). Activation or suppression of B cells may play an important role in CRC carcinogenesis, which needs to be identified in further studies. The difference of peripheral blood mast cell count between CRC and healthy cohorts has not been reported, and its relationship with survival remains controversial (42, 43). High levels of NKT cells were related to poor prognosis in this study; a similar result has been reported previously (38). In a recent study, a decreased level of circulating basophils was found linked to aggressive biology and poor survival, which is similar to the result of this study (44). In this study, a high level of immature dendritic cells predicted poor survival. Actually, a dendritic cell-infiltrating level has been reported to be positively correlated with layilin and a high layilin level was linked to poor survival in colorectal cancer patients (45). A lymphatic endothelial cell level was associated with poor survival in this study. Lymphatic vessel invasion has been identified as an independent prognostic factor for poor survival in colorectal cancer, and CRC-associated intestinal lymphatic endothelial cells were revealed to be able to regulate tumor progression (46). Further studies are needed to evaluate the role of peripheral blood immune cells in CRC progression and the potential of EVs estimating peripheral blood immune cells.

Furthermore, differentially expressed exLRs between CRC, CRA, and healthy cohorts were enriched in two pathways, the intestinal immune network for the IgA production pathway and the pathway of circadian rhythm of mammal. IgA deficiency is associated with a number of immune-mediated diseases, and it has also been proved to be associated with increased risk of gastrointestinal cancer in a nationwide population-based cohort study (47). Circadian rhythms of cell cycle–related molecule expression have been extensively reported (48). In a recently published study, circadian disruption was revealed to be associated with tumor-associated immune cell remodeling, resulting in facilitation of tumor growth (49).

Limitations and prospects of this study are listed as follows. First, the independent cohort size was limited and the diagnostic performance of the CRC d-signature needs to be validated in more independent centers. Second, we are continuing to recruit participants to identify the efficiency of the CRA d-signature. Third, the potential of EVs in predicting chemotherapy resistance is under study.

In summary, our study evaluated the value of exLRs serving as markers in the detection of CRC. The d-signature we have established can differentiate CRC from control (CRA and healthy) cohorts efficiently, which is supposed to improve CRC early detection efficiency in clinical practice. The exLR profiling can also indicate immune cell distribution and associated prognostic significance. We believe that this d-signature can contribute to the early detection of CRC and improve CRC prognosis in the near future.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Ethics Statement

The studies involving human participants were reviewed and approved by the Ethics Committee of the Fudan University Shanghai Cancer Center. The patients/participants provided their written informed consent to participate in this study.

Author Contributions

T-AG, H-YL, Z-ZZ, H-BH, S-LH, and YX were responsible for the study concept and study design. T-AG, H-YL, and CL performed the data acquisition. H-YL, Y-TJ, YL, and Y-CL were responsible for the methodology, software, formal analysis, and visualization. T-AG and H-YL wrote the original draft. YX, S-LH, and Z-ZZ edited and revised the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the National Natural Science Foundation of China (82072694, 81872294), the Shanghai Science and Technology Innovation Action Plan (20JC1419000), and the Shanghai Committee of Science and Technology (20DZ1100101, 19511121202).

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 all the participants included in this study.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2022.829230/full#supplementary-material

Abbreviations

CRC, colorectal cancer; exLR, extracellular vesicle long RNA; CRA, colorectal adenoma; AUC, area under the curve; gFOBT, guaiac fecal occult blood test; FIT, fecal immunochemical test for hemoglobin; CTC, CT colonography; EV, extracellular vesicle; mRNA, messenger RNA; lncRNA, long non-coding RNA; circRNA, circular RNA; exLR-seq, extracellular vesicle long RNA sequencing; d-signature, diagnostic signature; TEM, transmission electron microscopy; TPM, transcripts per million; DEG, differentially expressed gene; FC, fold change; mRMR, minimum redundancy maximum relevance; IFS, incremental feature selection; SVM, support vector machine; OS, overall survival; DFS, disease-free survival; ssGSEA, single sample gene set enrichment analysis; ANOVA, analysis of variance; NA, not available; PMBC, peripheral blood mononuclear cell; ROC, receiver operating characteristic; CEA, carcinoembryonic antigen; NKT cell, natural killer T-cell.

References

1. Sung H, Ferlay J, Siegel RL, Laversanne M, Soerjomataram I, Jemal A, et al. Global Cancer Statistics 2020: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA Cancer J Clin (2021) 71(3):209–49. doi: 10.3322/caac.21660

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Favoriti P, Carbone G, Greco M, Pirozzi F, Pirozzi REM, Corcione F. Worldwide Burden of Colorectal Cancer: A Review. Updates Surg (2016) 68(1):7–11. doi: 10.1007/s13304-016-0359-y

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Schreuders EH, Ruco A, Rabeneck L, Schoen RE, Sung JJY, Young GP, et al. Colorectal Cancer Screening: A Global Overview of Existing Programmes. Gut (2015) 64(10):1637–49. doi: 10.1136/gutjnl-2014-309086

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Taylor DP, Cannon-Albright LA, Sweeney C, Williams MS, Haug PJ, Mitchell JA, et al. Comparison of Compliance for Colorectal Cancer Screening and Surveillance by Colonoscopy Based on Risk. Genet Med (2011) 13(8):737–43. doi: 10.1097/GIM.0b013e3182180c71

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Pickhardt PJ, Hassan C, Halligan S, Marmo R. Colorectal Cancer: CT Colonography and Colonoscopy for Detection–Systematic Review and Meta-Analysis. Radiology (2011) 259(2):393–405. doi: 10.1148/radiol.11101887

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Ladabaum U, Dominitz JA, Kahi C, Schoen RE. Strategies for Colorectal Cancer Screening. Gastroenterology (2020) 158(2):418–32. doi: 10.1053/j.gastro.2019.06.043

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Yáñez-Mó M, Siljander PRM, Andreu Z, Bedina Zavec A, Borràs FE, Buzas EI, et al. Biological Properties of Extracellular Vesicles and Their Physiological Functions. J Extracellular Vesicles (2015) 4(1):27066. doi: 10.3402/jev.v4.27066

CrossRef Full Text | Google Scholar

8. Del Re M, Biasco E, Crucitta S, Derosa L, Rofi E, Orlandini C, et al. The Detection of Androgen Receptor Splice Variant 7 in Plasma-Derived Exosomal RNA Strongly Predicts Resistance to Hormonal Therapy in Metastatic Prostate Cancer Patients. Eur Urol (2017) 71(4):680–7. doi: 10.1016/j.eururo.2016.08.012

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Zhao R, Zhang Y, Zhang X, Yang Y, Zheng X, Li X, et al. Exosomal Long Noncoding RNA HOTTIP as Potential Novel Diagnostic and Prognostic Biomarker Test for Gastric Cancer. Mol Cancer (2018) 17(1):68. doi: 10.1186/s12943-018-0817-x

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Shao H, Im H, Castro CM, Breakefield X, Weissleder R, Lee H. New Technologies for Analysis of Extracellular Vesicles. Chem Rev (2018) 118(4):1917–50. doi: 10.1021/acs.chemrev.7b00534

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Li Y, Zhao J, Yu S, Wang Z, He X, Su Y, et al. Extracellular Vesicles Long RNA Sequencing Reveals Abundant mRNA, circRNA, and lncRNA in Human Blood as Potential Biomarkers for Cancer Diagnosis. Clin Chem (2019) 65(6):798–808. doi: 10.1373/clinchem.2018.301291

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Yu S, Li Y, Liao Z, Wang Z, Wang Z, Li Y, et al. Plasma Extracellular Vesicle Long RNA Profiling Identifies a Diagnostic Signature for the Detection of Pancreatic Ductal Adenocarcinoma. Gut (2020) 69(3):540–50. doi: 10.1136/gutjnl-2019-318860

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Dobin A, Davis CA, Schlesinger F, Drenkow J, Zaleski C, Jha S, et al. STAR: Ultrafast Universal RNA-Seq Aligner. Bioinformatics (2013) 29(1):15–21. doi: 10.1093/bioinformatics/bts635

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Liao Y, Smyth GK, Shi W. FeatureCounts: An Efficient General-Purpose Program for Assigning Sequence Reads to Genomic Features. Bioinformatics (2014) 30(7):923–30. doi: 10.1093/bioinformatics/btt656

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. Limma Powers Differential Expression Analyses for RNA-Sequencing and Microarray Studies. Nucleic Acids Res (2015) 43(7):e47. doi: 10.1093/nar/gkv007

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Yu G, Wang LG, Han Y, He QY. ClusterProfiler: An R Package for Comparing Biological Themes Among Gene Clusters. Omics (2012) 16(5):284–7. doi: 10.1089/omi.2011.0118

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Radovic M, Ghalwash M, Filipovic N, Obradovic Z. Minimum Redundancy Maximum Relevance Feature Selection Approach for Temporal Gene Expression Data. BMC Bioinf (2017) 18(1):9. doi: 10.1186/s12859-016-1423-9

CrossRef Full Text | Google Scholar

18. Sayed S, Nassef M, Badr A, Farag I. A Nested Genetic Algorithm for Feature Selection in High-Dimensional Cancer Microarray Datasets. Expert Syst Appl (2019) 121:233–43. doi: 10.1016/j.eswa.2018.12.022

CrossRef Full Text | Google Scholar

19. Chang C, Lin C. LIBSVM: A Library for Support Vector Machines. ACM Trans Intelligent Syst Technol (2011) 2(3):27. doi: 10.1145/1961189.1961199

CrossRef Full Text | Google Scholar

20. Aran D, Hu Z, Butte AJ. XCell: Digitally Portraying the Tissue Cellular Heterogeneity Landscape. Genome Biol (2017) 18(1):220. doi: 10.1186/s13059-017-1349-1

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Liberzon A, Subramanian A, Pinchback R, Thorvaldsdottir H, Tamayo P, Mesirov JP. Molecular Signatures Database (MSigDB) 3.0. Bioinformatics (2011) 27(12):1739–40. doi: 10.1093/bioinformatics/btr260

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Hanzelmann S, Castelo R, Guinney J. GSVA: Gene Set Variation Analysis for Microarray and RNA-Seq Data. BMC Bioinf (2013) 14(1):7. doi: 10.1186/1471-2105-14-7

CrossRef Full Text | Google Scholar

23. Moertel CG, O’Fallon JR, Go VL, O’Connell MJ, Thynne GS. The Preoperative Carcinoembryonic Antigen Test in the Diagnosis, Staging, and Prognosis of Colorectal Cancer. Cancer-Am Cancer Soc (1986) 58(3):603–10. doi: 10.1002/1097-0142(19860801)58:3<603::aid-cncr2820580302>3.0.co;2-k

CrossRef Full Text | Google Scholar

24. Alipoor SD, Mortaz E, Varahram M, Movassaghi M, Kraneveld AD, Garssen J, et al. The Potential Biomarkers and Immunological Effects of Tumor-Derived Exosomes in Lung Cancer. Front Immunol (2018) 9:819. doi: 10.3389/fimmu.2018.00819

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Del Re M, Marconcini R, Pasquini G, Rofi E, Vivaldi C, Bloise F, et al. PD-L1 mRNA Expression in Plasma-Derived Exosomes is Associated With Response to Anti-PD-1 Antibodies in Melanoma and NSCLC. Brit J Cancer (2018) 118(6):820–4. doi: 10.1038/bjc.2018.9

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Church TR, Wandell M, Lofton-Day C, Mongin SJ, Burger M, Payne SR, et al. Prospective Evaluation of Methylated SEPT9 in Plasma for Detection of Asymptomatic Colorectal Cancer. Gut (2013) 63(2):317–25. doi: 10.1136/gutjnl-2012-304149

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Tanaka S, Kashida H, Saito Y, Yahagi N, Yamano H, Saito S, et al. Japan Gastroenterological Endoscopy Society Guidelines for Colorectal Endoscopic Submucosal Dissection/Endoscopic Mucosal Resection. Dig Endosc (2020) 32(2):219–39. doi: 10.1111/den.13545

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Benson AB, Venook AP, Al-Hawary MM, Arain MA, Chen YJ, Ciombor KK, et al. NCCN Guidelines Insights: Rectal Cancer, Version 6.2020. J Natl Compr Canc Netw (2020) 18(7):806–15. doi: 10.6004/jnccn.2020.0032

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Zeng Z, Lu J, Wu D, Zuo R, Li Y, Huang H, et al. Poly(ADP-Ribose) Glycohydrolase Silencing-Mediated H2B Expression Inhibits Benzo(a)Pyrene-Induced Carcinogenesis. Environ Toxicol (2021) 36(3):291–7. doi: 10.1002/tox.23034

PubMed Abstract | CrossRef Full Text | Google Scholar

30. de Andrade LF, Lu Y, Luoma A, Ito Y, Pan D, Pyrdol JW, et al. Discovery of Specialized NK Cell Populations Infiltrating Human Melanoma Metastases. JCI Insight (2019) 4(23):e133103. doi: 10.1172/jci.insight.133103

CrossRef Full Text | Google Scholar

31. Harada H, Nagai H, Tsuneizumi M, Mikami I, Sugano S, Emi M. Identification of DMC1, a Novel Gene in the TOC Region on 17q25.1 That Shows Loss of Expression in Multiple Human Cancers. J Hum Genet (2001) 46(2):90–5. doi: 10.1007/s100380170115

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Salipante SJ, Mealiffe ME, Wechsler J, Krem MM, Liu Y, Namkoong S, et al. Mutations in a Gene Encoding a Midbody Kelch Protein in Familial and Sporadic Classical Hodgkin Lymphoma Lead to Binucleated Cells. Proc Natl Acad Sci USA (2009) 106(35):14920–5. doi: 10.1073/pnas.0904231106

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Dai HY, Hong CC, Liang SC, Yan MD, Lai GM, Cheng AL, et al. Carbonic Anhydrase III Promotes Transformation and Invasion Capability in Hepatoma Cells Through FAK Signaling Pathway. Mol Carcinog (2008) 47(12):956–63. doi: 10.1002/mc.20448

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Li M, Li A, He R, Dang W, Liu X, Yang T, et al. Gene Polymorphism of Cytochrome P450 Significantly Affects Lung Cancer Susceptibility. Cancer Med (2019) 8(10):4892–905. doi: 10.1002/cam4.2367

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Valle BL, Rodriguez-Torres S, Kuhn E, Diaz-Montes T, Parrilla-Castellar E, Lawson FP, et al. HIST1H2BB and MAGI2 Methylation and Somatic Mutations as Precision Medicine Biomarkers for Diagnosis and Prognosis of High-Grade Serous Ovarian Cancer. Cancer Prev Res (Phila) (2020) 13(9):783–94. doi: 10.1158/1940-6207.CAPR-19-0412

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Thompson BJ, Sahai E. MST Kinases in Development and Disease. J Cell Biol (2015) 210(6):871–82. doi: 10.1083/jcb.201507005

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Zhang T, Choi J, Kovacs MA, Shi J, Xu M, Goldstein AM, et al. Cell-Type-Specific eQTL of Primary Melanocytes Facilitates Identification of Melanoma Susceptibility Genes. Genome Res (2018) 28(11):1621–35. doi: 10.1101/gr.233304.117

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Krijgsman D, de Vries NL, Skovbo A, Andersen MN, Swets M, Bastiaannet E, et al. Characterization of Circulating T-, NK-, and NKT Cell Subsets in Patients With Colorectal Cancer: The Peripheral Blood Immune Cell Profile. Cancer Immunol Immunother (2019) 68(6):1011–24. doi: 10.1007/s00262-019-02343-7

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Shimabukuro-Vornhagen A, Schlosser HA, Gryschok L, Malcher J, Wennhold K, Garcia-Marquez M, et al. Characterization of Tumor-Associated B-Cell Subsets in Patients With Colorectal Cancer. Oncotarget (2014) 5(13):4651–64. doi: 10.18632/oncotarget.1701

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Punt CJ, Barbuto JA, Zhang H, Grimes WJ, Hatch KD, Hersh EM. Anti-Tumor Antibody Produced by Human Tumor-Infiltrating and Peripheral Blood B Lymphocytes. Cancer Immunol Immunother (1994) 38(4):225–32. doi: 10.1007/BF01533513

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Barbera-Guillem E, Nelson MB, Barr B, Nyhus JK, May KJ, Feng L, et al. B Lymphocyte Pathology in Human Colorectal Cancer. Experimental and Clinical Therapeutic Effects of Partial B Cell Depletion. Cancer Immunol Immunother (2000) 48:541–9. doi: 10.1007/pl00006672

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Mao Y, Feng Q, Zheng P, Yang L, Zhu D, Chang W, et al. Low Tumor Infiltrating Mast Cell Density Confers Prognostic Benefit and Reflects Immunoactivation in Colorectal Cancer. Int J Cancer (2018) 143(9):2271–80. doi: 10.1002/ijc.31613

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Mehdawi L, Osman J, Topi G, Sjolander A. High Tumor Mast Cell Density is Associated With Longer Survival of Colon Cancer Patients. Acta Oncol (2016) 55(12):1434–42. doi: 10.1080/0284186X.2016.1198493

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Liu Q, Luo D, Cai S, Li Q, Li X. Circulating Basophil Count as a Prognostic Marker of Tumor Aggressiveness and Survival Outcomes in Colorectal Cancer. Clin Transl Med (2020) 9(1):6. doi: 10.1186/s40169-019-0255-4

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Pan JH, Zhou H, Cooper L, Huang JL, Zhu SB, Zhao XX, et al. LAYN is a Prognostic Biomarker and Correlated With Immune Infiltrates in Gastric and Colon Cancers. Front Immunol (2019) 10:6. doi: 10.3389/fimmu.2019.00006

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Ungaro F, Colombo P, Massimino L, Ugolini GS, Correale C, Rasponi M, et al. Lymphatic Endothelium Contributes to Colorectal Cancer Growth via the Soluble Matrisome Component GDF11. Int J Cancer (2019) 145(7):1913–20. doi: 10.1002/ijc.32286

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Ludvigsson JF, Neovius M, Ye W, Hammarstrom L. IgA Deficiency and Risk of Cancer: A Population-Based Matched Cohort Study. J Clin Immunol (2015) 35(2):182–8. doi: 10.1007/s10875-014-0124-2

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Masri S, Cervantes M, Sassone-Corsi P. The Circadian Clock and Cell Cycle: Interconnected Biological Circuits. Curr Opin Cell Biol (2013) 25(6):730–4. doi: 10.1016/j.ceb.2013.07.013

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Aiello I, Fedele M, Roman F, Marpegan L, Caldart C, Chiesa JJ, et al. Circadian Disruption Promotes Tumor-Immune Microenvironment Remodeling Favoring Tumor Cell Proliferation. Sci Adv (2020) 6(42):eaaz4530. doi: 10.1126/sciadv.aaz4530

PubMed Abstract | CrossRef Full Text | Google Scholar