The sEVs exist in various body fluids, which is convenient for non-invasive detection (38). CircRNAs are stable, conservative, and specific expression of cells and tissues, which suggests that they have the potential to be used as molecular diagnostic and prognostic markers (39). sEV derived circRNAs combine the advantages of using sEVs with the specificity of circRNAs, enhancing their potential application as early non-invasive biomarkers. sEVs derived from pathological cells can carry their disease-specific circRNA into the peripheral blood. Therefore, the detection of sEV-derived circRNAs in serum may be feasible in the diagnosis of tumor disease. CircRNAs have also been considered as EV biomarkers to monitor the progression and chemoresistance of some type of cancers. In addition, it has been discovered that circRNAs are stably expressed in sEVs and these circRNAs are suggested to be a promising candidate for biomarkers in cancer (40).
CircRNA expression profile in mesothelioma
According to results of our previous study (37), 290 circRNAs derived from host genes PHKB, SLC45A4, ARHGEF28, FBXW4, TAF15, PLEKHM1, RALGPS1, STIL, L3MBTL4, ANKRD27, NHS, ILKAP, and PTK2 in PM cell lines were upregulated using high throughput human lncRNA microarrays through fold change (FC) (37). For this paper, we selected hsa_circ_0007386 (the one most highly expressed in four mesothelioma cells) as a representative circRNA for the PHKB gene (37) and investigated for changes in its expression levels to determine whether it could be employed as a potential biomarker for PM diagnosis.
sEV characterization
For sEV characterization, particle size and concentration were evaluated by NTA, the morphology was evaluated with Cryo-EM and the expression of the sEV common protein markers (CD63, CD81, and CD9) was assessed using Western blotting (Fig.2). The tetraspanins (CD63, CD81, and CD9) were detected in the sEVs derived from studied cell line using Western Blotting (Fig. 2 A & Fig. 1A-E Supplementary information). The size distribution of sEVs, measured by NTA and the concentration of sEVs enriched from cells is shown in Fig. 2B. As shown in Fig. 2C, most of the particles had a mean size of 100-200 nm. For more and precise characterisation of the sEVs, the Cryo-EM imaging was performed. Under cryo-EM, the specimens were imaged under extremely low temperature (below − 175 °C) so that sEVs retained its original spherical shape. The results of cryo-EM for sEVs confirmed their expected size and morphology (Fig. 2D).
The hsa_circ_0007386 RNA binding sites.
Numerous databases are available for RNA binding sites to circRNA targets. In this study, we have used the Circular RNA interactome database (https://circinteractome.nia.nih.gov/) to find out the RNA-binding sites matching. The results are listed in Fig. 3.
The hsa_circ_0007386 highly expressed in PM cell lines and sEVs.
For the digital PCR analysis, the negative control for the hsa_circ_0007386 was used to adjust the threshold for achieving the correct signal from the positive samples and the signal detected from the top right corner of the quadruplet, regarded as a positive control. In all samples, the background signal was used as the threshold above which signals were considered positive. Results of digital PCR revealed that among the studied cell lines, the hsa_circ_0007386 was overexpressed in NCI H-28 cells, with 1218.3 copies per µL input sample, followed by the H2052, H226, H2452, and MSTO-211H, with 696.7, 555.6, 412.5, and 187 copies per µL, respectively. The normal mesothelial cell line (MeT5A) showed 448 copies per µL. By normalizing the copy numbers using the ratio of hsa_circ_0007386 to hsa_circ_0000284, ratios were 0.183, 0.168, 0.098, 0.0943, and 0.0908 respectively (Fig. 4A). Based on findings by Zhong at al.,(41) hsa_circ_0000284 was selected as the internal control due to its demonstrated superior stability. This characteristic renders it an optimal candidate not only for circRNAs but also for broader RNA applications, serving as a reliable reference gene (41). The results also indicate that the copy number ratio of hsa_circ_0007386 to the circ RNA reference was significantly lower in non-PM cell lines including melanoma (Colo794; 0.0340, Colo679;0.0042, A375; 0.00755), gastric (Kato III; 0.0467, MKN45;0.0066), lung (H460; 0.031, H3122; 0.0167), colon (HCT 115;0.0216, HCT 116;0.0323), breast (MCF-7; 0.0235), liver (HepG2;0.068), and prostate (LnCAP; 0.00433) cancer cell lines (p<0.0001).
the hsa_circ_0007386 was overexpressed in H226 derived sEVs, with 36.23 copy per µL, followed by MSTO-211H, H2452, H2052, and NCI-H28, with 23, 11.71, 9.54, and 6.77 copies per µL of mixture, respectively. By normalizing the copy numbers found in sEVs using the hsa_circ_0007386 to hsa_circ_0000284 ratio, results differed from those attained using cells RNA extracts, and so did the ranking across cell lines as follows; H2452, NCI-H28, MSTO-211H, H2052, and H226 with the following ratios 0.187, 0.171, 0.121, 0.1045, and 0.087, respectively (Fig. 4B). Interestingly, the copy number ratio of this circRNA biomarker to the circ RNA reference was significantly lower in non-PM cell lines including melanoma (Colo 794; 0.0580, Colo679;0.0053, A375; 0.0084), gastric (Kato III; 0.055, MKN45;0.0054), lung (H460; 0.036, H3122; 0.032,), colon (HCT 116;0.0452, HCT 115; 0.018), breast (MCF-7; 0.0233), and prostate (LnCAP; 0.0066) cancer cell lines (p<0.0001). Indicating that the sorting of specific circRNA species to sEVs compared to the cells may be actively regulated.
Results of our study revealed the specificity of the has_circ_0007386 derived from both cells and sEVs as a specific biomarker for early diagnosis of PM, as we have observed significant expression of this biomarker in the PM cell derived sEVs (Fig. 5A) and PM cell lines (Fig. 5B) compared to the normal mesothelial cell line (Met5A) and the microglial cell line (BV-2) (p<0.0001).
The results of our study on PM and non-PM cells derived RNAs confirm the 90% sensitivity and 93.3% specificity of the hsa_circ_0007386 as a biomarker in PM diagnosis. However, this number was 81.81% and 87.5% in the sEV respectively. The specificity was calculated by dividing the number of true negative samples which are here non-PM cell lines to the sum of true negative and false positive which are here sum pf non-PM cell lines and the cell lines that shows higher ratio of has_circ_0007386 to reference circRNA compared to the PM cell lines (42). Similarly for the sensitivity calculation, the number of true positive samples which are here PM cell lines to the sum of true positive and false negative which are here sum of PM cell lines and the cell lines that shows higher ratio of has_circ_0007386 to reference circRNA compared to the PM cell lines (42). While our results demonstrate consistencies in the copy numbers or copy number ratios of hsa_circ_0007386 and hsa_circ_0000284 between PM, and non-PM derived sEV samples and their corresponding cellular counterparts, the findings also reveal a higher copy number ratio of hsa_circ_0007386 to the reference hsa_circ_0000284 in sEV-derived samples compared to the cells. The heightened presence of hsa_circ_0007386 within PM-derived sEVs, as compared to its cellular counterpart, underscores the propensity of circRNAs for encapsulation within sEVs. By scrutinizing the contents of sEVs, we can attain more precise and sensitive findings than those derived from the analyses of cellular extracts alone. Notably, our investigation revealed the maximal copy number of sEV-derived hsa_circ_0007386 within MSTO-211H, representative of the biphasic variant of PM. Early identification of this PM subtype holds substantial promise for advancing our understanding of disease progression. Consequently, the potential clinical value of sEV-derived hsa_circ_0007386 in less-invasive liquid biopsy approaches deserves consideration, with the outlook of replacing traditional tissue biopsies in the future. Our findings underscore the diagnostic specificity of hsa_circ_0007386 as a viable biomarker for early PM detection, with significantly elevated expression observed in PM cell lines compared to non-PM, normal mesothelial, and other cancer cell lines (p-value < 0.0001). This highlights the promising diagnostic potential of hsa_circ_0007386 across both sEV and cellular contexts. To ensure precision and reliability, further validation utilizing patient samples is imperative for consolidating these outcomes. In different stages of different development diseases, disease-related circRNA can be sorted into sEVs to be enriched and transported to target cells or target organs for release. Many studies have shown that differential expression of sEV-derived circRNAs in the body fluid was associated with the pathological characteristics of tumor vascular invasion, lymph node metastasis, poor survival and TNM stage (43-45). Studies have shown that circRNAs can be packaged and function in sEVs (46). However, the mechanism behind the selective packaging of specific circRNAs into sEVs is not yet clear and requires further investigation. In a study by Zhang et al (47), it was showed that circRNA polo-like kinase 1 (circPLK1) was upregulated in Malignant Pleural Mesothelioma (MPM) tumor tissues and cell lines. CircPLK1 knockdown suppressed the proliferation, migration, invasion and stemness of MPM cells in MPM progression. The studies mentioned have indicated the diagnostic efficacy of different circ-RNAs as a cancer marker.
In the PM cell derived sEVs, a notable observation arose when comparing the copy number ratio of hsa_circ_0007386 to the circRNA reference across various PM cell derived sEVs, particularly in the context of one lung cancer cell line (H1975), and one liver cancer (HepG2). Despite an apparently higher ratio in H1975 (0.1630), and HepG2 (0.086) compared to other PM cell lines (H20525, MSTO-211H, and H226) a deeper analysis utilizing raw data and absolute copy number calculations revealed that the observed copy number of this putative circRNA biomarker in the above-mentioned cell lines did exhibit statistical significance (p<0.0001) (Fig. 6A). The decision to employ exact copy number calculations was driven by the criterion of fewer than 100 copies per 40 µL of digital PCR sample, aiming to enhance result accuracy. Similarly, in the PM cell lines, comparable to our observations in the cell derived sEVs, a notable observation arose when comparing the copy number ratio of hsa_circ_0007386 to the circRNA reference across various PM cell lines, particularly in the context of one lung cancer cell line (H1975). Despite an apparently higher ratio in H1975 (0.12) compared to other PM cell lines (MSTO-211H, H226, and H2452), a deeper analysis utilizing raw data and absolute copy number calculations revealed that the observed copy number of this putative circRNA biomarker in the H1975 cell line did exhibit statistical significance (p<0.0001) (Fig. 6B). The decision to employ exact copy number calculations was driven by the criterion of fewer than 100 copies per 40 µL of digital PCR sample, aiming to enhance result accuracy (48).
Results of a systematic review and meta-analysis revealed that circRNAs have the potential to be biomarkers for diagnosis and prognosis of cancers (49). It has been described by Stella et al. that two circRNAs that localised in serum derived sEVs including circSMARCA5 (hsa_circ_0001445) and circHIPK3 (hsa_circ_0000284) could be potential biomarkers for glioblastoma and could distinguish glioblastoma patients from healthy controls with high accuracy (50). In contrast, another study has reported that the continued high expression of sEV derived circRNA-100338 in the serum of HCC HCC (Hepatocellular Carcinoma) patients undergoing therapeutic hepatectomy may be related to lung metastasis and poor survival (51).
There are some advantages in analyzing sEV-encapsulated non-coding RNAs compared to whole plasma/serum. Firstly, as extracellular vesicles can be secreted by a variety of cells, the contents of sEVs can be used as biomarkers for diagnosis or prognosis in various diseases (52). Secondly, it is easier to sort circRNA into sEVs than linear RNAs (35). In addition, sEVs derived from cancers contain highly specific RNA, and they can also prevent the nucleic acid molecules from degradation by RNase in the blood (53). However, there are still many issues to be resolved before sEVs-derived circRNAs can be employed as reliable biomarkers, such as preservation of specimens, cell source of sEVs, sEVs isolation methods, etc.
Future investigations focusing on sEV circRNAs in various biological contexts, such as the hematopoietic system, immune response, nervous disorders, cancer development, and other diseases, will provide further insights into the enigmatic nature of sEV circRNAs. Consequently, elucidating the mechanisms of cancer pathogenesis and identifying potential novel diagnostic biomarkers or therapeutic targets are expected to be prominent areas of research in the future.