Next Article in Journal
Genetic Variations of the DPYD Gene and Its Relationship with Ancestry Proportions in Different Ecuadorian Trihybrid Populations
Next Article in Special Issue
Evaluation of Neutrophil-to-Lymphocyte Ratio (NLR), Platelet-to-Lymphocyte Ratio (PLR) and Systemic Immune–Inflammation Index (SII) as Potential Biomarkers in Patients with Sporadic Medullary Thyroid Cancer (MTC)
Previous Article in Journal
Atomic Force Microscopy Application for the Measurement of Infliximab Concentration in Healthy Donors and Pediatric Patients with Inflammatory Bowel Disease
Previous Article in Special Issue
Associations of Complete Blood Count Parameters with Disease-Free Survival in Right- and Left-Sided Colorectal Cancer Patients
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Machine Learning-Assisted FTIR Analysis of Circulating Extracellular Vesicles for Cancer Liquid Biopsy

1
Fondazione Policlinico Universitario “A. Gemelli” IRCCS, 00168 Rome, Italy
2
Dipartimento di Neuroscienze, Sezione di Fisica, Università Cattolica del Sacro Cuore, 00168 Rome, Italy
3
Sezione di Medicina Interna, Dipartimento di Medicina e Chirurgia Traslazionale, Università Cattolica del Sacro Cuore, 00168 Rome, Italy
4
Dipartimento di Scienze di Laboratorio e Infettivologiche, Fondazione Policlinico Universitario “A. Gemelli” IRCCS, 00168 Rome, Italy
*
Author to whom correspondence should be addressed.
These authors shared the seniorship of this work.
J. Pers. Med. 2022, 12(6), 949; https://doi.org/10.3390/jpm12060949
Submission received: 12 May 2022 / Revised: 1 June 2022 / Accepted: 5 June 2022 / Published: 10 June 2022
(This article belongs to the Special Issue Cancer Biomarker Research and Personalized Medicine 2.0)

Abstract

:
Extracellular vesicles (EVs) are abundantly released into the systemic circulation, where they show remarkable stability and harbor molecular constituents that provide biochemical information about their cells of origin. Due to this characteristic, EVs are attracting increasing attention as a source of circulating biomarkers for cancer liquid biopsy and personalized medicine. Despite this potential, none of the discovered biomarkers has entered the clinical practice so far, and novel approaches for the label-free characterization of EVs are highly demanded. In this regard, Fourier Transform Infrared Spectroscopy (FTIR) has great potential as it provides a quick, reproducible, and informative biomolecular fingerprint of EVs. In this pilot study, we investigated, for the first time in the literature, the capability of FTIR spectroscopy to distinguish between EVs extracted from sera of cancer patients and controls based on their mid-IR spectral response. For this purpose, EV-enriched suspensions were obtained from the serum of patients diagnosed with Hepatocellular Carcinoma (HCC) of nonviral origin and noncancer subjects. Our data point out the presence of statistically significant differences in the integrated intensities of major mid-IR absorption bands, including the carbohydrate and nucleic acids band, the protein amide I and II bands, and the lipid CH stretching band. Additionally, we used Principal Component Analysis combined with Linear Discriminant Analysis (PCA-LDA) for the automated classification of spectral data according to the shape of specific mid-IR spectral signatures. The diagnostic performances of the proposed spectral biomarkers, alone and combined, were evaluated using multivariate logistic regression followed by a Receiving Operator Curve analysis, obtaining large Areas Under the Curve (AUC = 0.91, 95% CI 0.81–1.0). Very interestingly, our analyses suggest that the discussed spectral biomarkers can outperform the classification ability of two widely used circulating HCC markers measured on the same groups of subjects, namely alpha-fetoprotein (AFP), and protein induced by the absence of vitamin K or antagonist-II (PIVKA-II).

Graphical Abstract

1. Introduction

In many clinical situations, cancer diagnosis requires single or repeated tissue biopsies of suspected cancerous regions. This procedure is invasive and sometimes associated with pain, discomfort, and risk for the patients. Additionally, the tissue region that needs to be sampled can be hardly or completely inaccessible by interventional procedures, and/or highly heterogeneous, thus leading to ambiguous conclusions. These drawbacks might limit the frequency with which a region can be sampled to check for cancer, thus hindering the possibility to perform an accurate diagnosis, especially in the early stages of the pathology. In this regard, liquid biopsy offers a promising diagnostic alternative, as it relies on the analysis of nonsolid biological tissues, such as blood. As such, it is noninvasive or minimally invasive and can be done more frequently by taking multiple liquid biopsies in a short period, being instrumental for diagnostic purposes and therapy monitoring [1].
Extracellular vesicles (EVs) are lipid-bound vesicles secreted by most cell types into the extracellular space [2,3,4,5,6,7], they abound in bio-fluids and harbor molecular constituents from their parent cells, including proteins, lipids, and genetic material. Due to their easy accessibility and capability of representing their cells of origin, EVs are attracting growing interest as cancer biomarkers in the field of liquid biopsy and personalized medicine. Despite this potential, EVs have still not been widely applied in diagnostics and novel approaches for their characterization are highly demanded.
In this regard, Fourier-transform infrared spectroscopy (FTIR), especially in the Attenuated Total Reflection (ATR) mode, is emerging as a promising tool for label-free molecular profiling of EVs, with several papers published on the subject [4,8,9,10,11,12,13,14,15,16,17,18,19,20]. Most of these works involve in vitro experiments aimed at assessing sample composition and purity [15,16], distinguishing different EV classes [18] or EVs derived from cells under different states/phenotypes/culture conditions [4,8,10,12,21,22,23]. So far, very few FTIR papers have focused on EVs purified from biofluids of subjects enrolled in clinical studies. These papers include Zlotogorski-Hurvitz et al., who studied salivary EVs of patients diagnosed with oral cancer [11], and Yap et al., who characterized EVs purified from the urine of prostate cancer patients [14,24]. To the best of our knowledge, the pioneering paper of Martins et al. on EVs obtained from neurological patients is the only article that uses serum as a biofluid of origin [19], and no studies have been yet published on cancer. Therefore, more data are needed to confirm the capability of FTIR spectroscopy to correctly classify EVs obtained from the serum of patients in different clinical conditions, with emphasis on cancer patients, as well as to test classification algorithms suitable for the high throughput analysis of FTIR spectra of EVs. The latter point is extremely important as EV spectral data can be automatically analyzed with machine learning (ML) methods, giving physicians a direct diagnosis instead of a biophysical parameter that needs further interpretation [8,13,19,20].
This pilot study aims to support the evidence that ATR-FTIR can be used as a fast, reproducible, and effective method for the classification of EV-enriched suspensions purified from sera of cancer patients. For this purpose, samples were extracted from a cohort of 20 patients diagnosed with Hepatocellular Carcinoma (HCC) of nonviral etiology and 19 noncancer subjects. This choice is motivated by the large clinical and research interest in serum biomarkers for HCC screening [25,26,27,28,29,30] and, also, by the growing number of papers reporting potential EV-based HCC biomarkers [31,32]. Here, multivariate statistical analysis and machine learning methods are used for the comparison of mid-IR spectra in the two groups, as schematically represented in Figure 1. Our data highlight the presence of statistically significant differences in the intensity and the shape of several mid-IR spectral bands that could be used as sensitive spectral biomarkers of the pathology. To further stress the relevance of the proposed approach for the development of novel liquid biopsy techniques, we compared the classification performances of our spectral biomarkers with those of serum alpha-fetoprotein (AFP), a glycoprotein produced in early fetal life by the liver and by a variety of tumors, which is currently the most widely used circulating marker in HCC with prognosis and risk of recurrence purposes. Additionally, we measured the serum levels of the protein induced by the absence of vitamin K or antagonist-II (PIVKA-II), a more sensitive circulating marker than AFP for differentiating HCC at all stages in patients with cirrhosis or chronic hepatitis. Notably, a Receiving Operator Curve (ROC) Analysis showed that our spectral biomarkers outperformed both AFP and PIVKA-II in distinguishing the two groups, a result that deserves further studies based on a larger sample size.

2. Materials and Methods

2.1. Clinical and Laboratory

All the subjects enrolled in this study have given written informed consent. We selected 20 patients with a recent diagnosis of nonviral HCC, according to EASL guidelines [33]. The exclusion criteria were as follows: history of previous treatment (i.e., hepatic resection, liver transplant, anti-angiogenic drugs, and radiofrequency); diagnosis of extra-hepatic neoplasia; HCV or HBV infection; Child-Pugh B or C; obstructive jaundice; and estimated creatinine clearance < 30 mL/min. The inclusion criteria were as follows: Hepatocellular carcinoma at diagnosis; no pharmacological treatment; Child-Pugh A; age > 18. Age- and gender-matched subjects were included in the control group. The exclusion criteria of controls were as follows: evidence of neoplasia; HCV or HBV infection; the presence of monoclonal components; anticoagulant therapy; presence of monoclonal component; and negative inflammatory markers.
Serum samples and clinical data of subjects were processed anonymously. A wide panel of biochemical parameters was measured on the collected blood samples, including Glucose, Cholesterol, Total bilirubin, Albumin, Alanine aminotransferase (ALT), Asand partate aminotransferase (AST), Alkaline phosphatase (ALP), and Gamma-Glutamyl Transferase (GGT). Patients’ sera were collected for quantitative measurement of AFP and PIVKA-II, performed using Lumipulse® G (Fujirebio, Tokyo, Japan), based on Chemiluminescent Enzyme [25,33]. Vascular invasion, portal trunk, and/or main portal branches, were diagnosed by MRI with hepato-specific contrast medium.

2.2. Sample Purification

The subjects enrolled in this study underwent peripheral blood sampling and then serum samples were obtained and treated as reported in previous studies [8,13]. EVs were isolated from serum by using the ExoQuick ULTRA precipitation kit (System Biosciences, Palo Alto, CA, USA). The isolation was performed according to the ExoQuick ULTRA manual and within 24 h of the serum sampling from patients. Briefly, blood cell debris was removed from the samples by two sequential centrifugations (3000× g for 15 min and 12,000× g for 10 min, respectively) and the ExoQuick reagent (67 µL) was added to 250 µL of each clarified biofluid. After 30 min of incubation, final centrifugation was performed and the pellet, which contains precipitated EVs, was resuspended in the ExoQuick buffers and finally purified using the provided columns. The samples were analyzed without any further manipulations. Regarding sample purity, a caveat was necessary. As reported in [34,35], the described purification protocol is likely to introduce contaminants from serum, such as lipoproteins. For this reason, our preparation is more adequately referred to as EV-enriched suspension, rather than a pure EV sample. Due to the presence of such contaminants, we decided not to assign the observed spectral signatures to specific molecular species. Conversely, we preferred to use robust statistical methods to highlight clinically relevant differences between the two groups that can be instrumental to develop novel HCC diagnostic tools. A structural and morphological characterization of the samples extracted with this method is reported in Figure S1.

2.3. FTIR Measurements and Data Analysis

The EV-enriched suspensions isolated from sera were analyzed using a Bruker ALPHA II compact FTIR Spectrometer, equipped with an attenuated total reflection (ATR) module (Eco-ATR). According to several papers investigating the spectral characteristics of EVs using ATR-FTIR, sample solution droplets were deposited on the ATR crystal and a thin film was obtained by slow evaporation of the solvent under ambient conditions. Specifically for each measurement, 5 μL of the EV solution was air-dried for at least 10–15 min on a high-throughput ZnSe crystal [4,8,9,10,11,12,13,14,15,16,17,18,19,20]. This step was repeated three times, resulting in a total volume of 15 μL of dried solution. All IR spectra were acquired in the wavelength range included from 4000 to 650 cm−1 and, for each spectrum 54 scans, at a resolution of 2 cm−1 was averaged. The background was acquired before the measurements and then subtracted from each sample spectra. The spectra were registered and preprocessed using the commercial OPUS 8.5 SP1 software, dedicated to the analysis of IR spectral data. Data were then exported and further analyzed with the programming software R. Spectra were normalized after linear baseline subtraction. Average Spectra were computed together with the corresponding 95% confidence intervals. Data were visualized with the ggplot2 software package. Principal Component Analysis (PCA) and Linear Discriminant Analysis (LDA) of spectra were carried out after subtracting a linear baseline from the analyzed absorption band. This analysis protocol was chosen as PCA alone, as well as PCA-LDA in combination, were successfully exploited in similar papers for discriminating EVs of different origins according to their vibrational fingerprints [8,9,10,11,12,13,14,15,16]. Additionally, PCA-LDA has been widely employed in other FTIR-based classification problems in different fields, including diagnostics and forensic sciences [17,18,19,20,21]. More specifically, PCA is an unsupervised statistical learning technique and was used here to find a low-dimensional representation of spectra that retains as much as possible of the variation in the original dataset [22]. Mathematically, this is achieved by computing the covariance matrix of the whole dataset and finding the matrix eigenvalues with the corresponding eigenvectors. The supervised LDA statistical learning algorithm was then used to maximize the separation between the two groups, a task that is performed by solving an eigenvalue problem similar to PCA but starting from the matrix obtained by dividing the between- and the within-class scatter matrices.
Differences between the two groups were assessed using a Wilcoxon non-parametric test for independent samples. Sensitivity and Specificity values were computed from the confusion matrix and validated using the Leave-One-Out-Cross-Validation (LOOCV) method. In this procedure, the database with the 39 samples is divided into 39 groups. At each algorithm iteration, one subset is selected as the testing set, while the remaining ones are used as the training set. The model is verified by the testing set, thus obtaining a classification rate. The mean of 39 classification rates is calculated and used as a robust statistical estimator [20].

3. Results

3.1. The Ratio of Different Molecular Classes within EVs Are Different in Cancer Patients and Controls

In this study, we compare the mid-IR spectral response of EV-enriched suspensions obtained from sera of HCC patients and noncancer subjects to search for novel spectral circulating biomarkers. The demographic, clinical, biochemical, and spectral characteristics of the recruited subjects are summarized in Table 1. In Figure 2A,B, two mid-IR absorbance spectra measured on the two groups of subjects are reported. Spectra are computed by averaging the data acquired on EV-enriched suspensions obtained from each of the recruited subjects, separately. For this purpose, absorbance was first normalized after linear baseline subtraction (see Material and Methods), 95% confidence bands (shaded area) are superimposed on the average curves (continuous lines).
At the level of the whole mid-IR range, the average curves look qualitatively resemblant, showing similar spectral signatures. A more in-depth analysis of Figure 2 shows the presence of changes in the spectral shape and in the relative intensity of several absorption bands, namely the carbohydrates and nucleic acids band at 1000–1200 cm−1 [16,19], the protein Amide I and II bands [13,36,37], the lipid C-H stretching band at 2800–3000 cm−1 [23,38], and the C=O stretching at approximately 1735 cm−1, which has been assigned to purine base and ester groups of lipids in EV samples [9,36,37]. These bands are highlighted in green and red color in controls and patients, respectively, and are considered a major source of spectral biomarkers in EV research [4,8,9,10,11,12,13,14,15,16,17,18,19,20,23]. To investigate more in-depth the difference between the two groups, an enlarged detail of the four mentioned IR bands for controls (green) and oncologic patients (red) is reported in Figure 3. Linear baseline subtraction, followed by area normalization, was performed before averaging curves.
A box-plot analysis of the computed areas is reported in Figure 4 for the amide I and II bands at 1470–1700 cm−1 (A), the lipid CH stretching band at 2800–3000 cm−1 (B), the carbohydrates and nucleic acids band at 1000–1200 cm−1 (C), and the C=O stretching peak at 1720–1740 cm−1 (D). The result of a Wilcoxon rank–sum test for independent samples was superimposed on each plot. Data are summarized in Table 1. Statistically significant differences in the computed areas were detected for all the studied bands, with emphasis on the CH stretching band that displays a remarkably low p-value. Interestingly, highly significant differences are detected by computing the ratios between the CH-stretching band at 2800–3000 cm−1 and the Amide I-II bands at 1470–1700 cm−1 (Figure 4E), which is known as the spectral Lipid-to-Protein Ratio (LPR) [8,18], and provide a quantitative index correlated to the relative amount of lipids and proteins in the analyzed samples. Remarkable differences are also observed in Figure 4F that reports the ratio between the CH-stretching band and the absorption peak at 1000–1200 cm−1, which can be defined—in analogy to Figure 4E—as the spectral Lipid-to-Nucleic acids Ratio (LNR). Taken together, Figure 4E,F hint at the occurrence of an alteration in the ratio between different molecular classes in EVs obtained from cancer patients compared to controls.

3.2. Machine Learning-Assisted Classification of HCC Patients and Controls

Aside from intensity values, differences in the average line shape between the two groups can be observed for all the IR bands reported in Figure 3, except for the carbohydrates and nucleic acid band. In Figure 5A–I, we investigate the possibility to use multivariate statistical and Machine Learning methods to automatically classify patients according to their group membership based on the spectral shape of the analyzed bands. For this purpose, we show the projections of the spectral data in the plane of the two principal components (PC1 and PC2) for the amide I–II region (A), the CH-stretching band of lipids (D), and the CO stretching peak at 1720–1760 cm−1 (G). We then used Linear Discriminant Analysis (LDA) to point out the direction (LD1) that maximizes the separation between the two groups, which is represented as a continuous black line together with the corresponding datapoints projection. The LD1 discriminant scores for the three analyzed bands are reported in Figure 5B,E,H (lower panels), respectively. A sufficient degree of clustering among data points belonging to the two groups can be observed, as confirmed by the corresponding box-plot analysis (upper panels), which points out the presence of statistically significant differences. In this context, remarkably low p-values are measured for the IR band analyzed in Figure 5H (p = 4.9 × 10−5). The PCA-LDA classification performances were tested with the leave-one-out-cross validation method (LOOCV) and the results are visualized in Figure 5C,F,I, respectively. These results are extremely interesting if compared with those obtained on conventional HCC circulating markers, such as AFP, which shows a sensitivity of 0.61 (95% CI 0.60–0.62) and a specificity of 0.86 (95% CI 0.85–0.87) in the 20–100 ng/mL concentration range, as demonstrated by a recent metanalysis [28].

3.3. A Combined Spectral Biomarker for HCC Diagnosis

In Figure 6, we evaluate the classification performance of each of the spectral biomarkers highlighted in Figure 4 and Figure 5 using receiving operator characteristics (ROC) curves and computing the corresponding Areas Under the Curve (AUC). The same analysis is carried out for two widely studied HCC circulating markers, namely AFP and PIVKA-II. The corresponding AUC values are summarized in Table 2, together with the corresponding confidence intervals. Interestingly enough, all the measured AUCs are statistically different from the random classifier (AUC = 0.5), which is represented by the continuous black line in Figure 6A. In Figure 6B, we investigate the possibility to use biomarkers in combination to improve the accuracy of the classification. For this purpose, we performed a multivariate logistic regression including all the markers in Table 2, followed by a stepwise logistic regression aimed at selecting the most informative set of parameters. This procedure selects two spectral biomarkers that can be used in combination, namely the integrated area of the CH stretching band at 2800–3000 cm−1 and the LD1 coefficients computed for the spectral signature at 1720–1760 cm−1. The result of the logistic regression performed on the stepwise model shows that the two selected markers possess statistically significant regression coefficients. The ROC curve of the stepwise model is reported in red in Figure 6B, showing a very large area under the curve (0.91, 95% CI: 0.8–1), which confirms the diagnostic potential of the combined spectral biomarker. Youden’s method was used to compute the sensitivity and specificity starting from the combined ROC curve and obtaining 1.00 and 0.76, respectively. ROC curves obtained on the same groups of subjects for AFP (gold curve) and PIVKA (violet curve) are reported for comparison in Figure 6B, obtaining lower—albeit comparable—AUC values, namely 0.81 and 0.86, respectively discussed.

4. Discussion

In this paper, we discuss the proof of concept of a novel liquid biopsy approach for cancer diagnosis. Our method is based on the molecular profiling of EVs with FTIR spectroscopy. For this purpose, EV-enriched suspensions were obtained from the serum of patients diagnosed with HCC of nonviral origin and noncancer subjects. We have chosen to investigate EVs because of their crucial role in the emerging liquid biopsies for cancer diagnosis [39,40,41,42,43]. Such a pivotal role is due to their large availability in biofluids and specific molecular content, which reflects the state of the cells of origin. Additionally, EV-specific molecular cargoes change during cancer evolution, thus being also instrumental for disease staging [44]. As far as HCC is concerned, EVs play a particularly relevant role in the disease progression, as cancer cells are capable of influencing many biological pathways through the release of extracellular vesicles in the tumor microenvironment [32]. These effects include a local regulation of the Epithelial to Mesenchymal Transition (EMT) [45], and the activation of the inflammatory microenvironment to increase cancer cell invasiveness, converting fibroblasts and macrophages to CAFs and TAMs [32,46,47]. Moreover, HCC-related EVs are believed to regulate the functions of immune and endothelial cells, inducing immune escape and angiogenesis [48,49].
HCC diagnosis is currently carried out by combining different types of information, which include the levels of selected circulating markers, medical imaging (US/MRI/CT), and histopathological biopsies [25,33,50]. Unfortunately, all these methods present some limitations, which negatively affect patient prognosis. Circulating markers, AFP included, have usually low sensitivity and specificity, although improved performance can be obtained through the combination of multiple markers [25,26,27,29,33]. Despite a remarkable specificity, medical imaging has a poor capability of detecting small tumors, being plagued by a large number of false negatives (FNs) and, consequently, a poor sensitivity. Histopathological biopsies also present several problems such as their invasiveness and a large number of FNs, which are partly due to the intrinsic variability of the bioptic tissue that complicates the sampling of the suspected region [25,32,33,51,52]. EV-based liquid biopsy can contribute to overcoming some of these shortcomings, as it is non- or minimally invasive and can provide a homogeneous sampling. Here, we used multivariate statistical analyses and Machine Learning (ML) methods to search for IR spectral biomarkers of the disease. We highlighted relevant differences in, both, the integrated intensity (Figure 4) and the line shape (Figure 5) of major mid-IR absorption bands of biomolecules within EVs. Specifically, HCC patients presented with a decreased intensity of the amide I and II bands (1470–1700 cm−1) and the nucleic acid and carbohydrate band (1000–1200 cm−1) compared to the controls, together with an increased intensity of the lipid CH stretching band (2800–3000 cm−1). These findings hint at a difference in the ratio of molecular classes within EVs in the two groups, as confirmed by the measured spectral Lipid to Protein ratio (LPR, Figure 4E) and the Lipid to carbohydrates and Nucleic acid Ratio (Figure 4F). This point deserves a more in-depth discussion as it helps point out the complementarity of FTIR with respect to conventional omics techniques for the biochemical characterization of EVs. Genomics, proteomics, and lipidomics, indeed, provide detailed quantitative information on the EV molecular cargos but involve complex sampling procedures that can alter the ratio among different molecular classes. On the contrary, FTIR spectroscopy appears to be perfectly suited to provide semiquantitative information on the relative amount of lipids, proteins, DNA, and RNA in EVs, also highlighting possible biochemical changes that depend on the clinical conditions of patients [8,9]. Investigating quantitative parameters, such as the shape or intensity of spectra, is a reasonable approach, as long as each sample has been uniformly processed. With this idea in mind, several works applied ML algorithms for the analysis of several circulating biomarkers and the development of novel cancer detection methods [11,53,54,55,56,57,58]
In this context, the same LPR here exploited has been also used in the pioneering paper of Mihàly et al. for discriminating among different classes of EVs [18] and in Romanò et al. to classify EVs obtained from cancer cells in different stages of the epithelial to mesenchymal transition [8].
Aside from the relative bands’ intensity, we found statistically significant differences in the shape of three specific mid-IR bands encompassing the spectral ranges 1470–1700 cm−1 (Figure 5A), 1720–1760 cm−1 (Figure 5B), and 2800–3000 cm−1 (Figure 5C). These differences were pointed out using ML, thus allowing us to automatically classify patients according to their clinical conditions. We believe that this approach is highly interesting in a clinical setting as it is designed to provide physicians with a direct diagnostic response instead of a set of spectral parameters that need to be interpreted by specialized personnel. Additionally, compared to IR peak assignment, the use of ML for spectral classification has the potential to reduce the detrimental effects due to the presence of non-EV contaminants in purified samples, a well-known problem for vibrational spectroscopy of EV samples and associated with the lack of a reproducible gold standard purification method to isolate EVs from serum with high yield and purity [4,59,60,61]. For instance, Lee et al. investigated SERS spectra from EVs purified with several isolation techniques and concluded that commercial kit reagents had a strong affinity to EVs leading to relevant coprecipitation [62]. Accordingly, Pereira et al. showed that a significant portion of the IR spectral signals of EVs isolated from cell culture media is ascribed to kit reagents [10]. Therefore, coprecipitants from the isolation kits along with serum proteins, such as albumin, lipoproteins, or immunoglobulins, are potential contaminants that may affect the downstream analyses of EVs. In such a complex condition, the IR peak assignment is extremely challenging and potentially misleading, while the results of an ML-based classification are likely to be more robust, as the interaction between EV and contaminants, which may lead to coprecipitation, may be different in different types of samples. This is especially true if one considers recent findings demonstrating the formation of a protein corona around extracellular vesicles dissolved in complex media [63].
EV-based HCC biomarkers were reviewed in the excellent papers of Wang and colleagues [32] and Szabo & Momen-Heravi [64], which showed that microRNAs, long noncoding RNAs (lncRNAs), mRNAs, circular RNAs, and proteins within EVs are promising sources of HCC biomarkers for diagnosis and staging. MicroRNAs are highly abundant in EVs and difficult to be degraded, thus being an ideal biomarker candidate. Altered microRNA levels in HCC patients have been measured for EV miR-9-3p[65], miR-21[66], miR-93[67], miR-92b [68], miR-718 [69,70], miR-122 [71,72], miR-638 [73,74,75], and miR-125b [76], to mention a few. Altered EV lncRNA (LINC00161[77], lncRNA-HEIH [78,79], ENSG00000258332.1, and LINC000635[80], lncRNA-ATB[81]), mRNA (mRNA-hnRNPH1), circular RNA (circPTGR1[82], and circDB [83]) levels were also detected in HCC patients compared to controls. EV proteins are also important markers in the diagnosis of HCC and alterations in the levels of several proteins have been detected including, G3BP, PIGR [84], hepcidin [85], and SMAD3 protein[86]. Despite the promising diagnostic role of the mentioned RNA and protein structures, clinical trials have pointed out that a single biomarker at a time display low sensitivity and specificity. Comprehensive detection and evaluation using multiple molecules in combination may be more effective in diagnostics [25,32,33]. In this context, the spectral differences highlighted in the intensity of the nucleic acid and carbohydrates band and the shape and intensity of the amide I and II bands are highly interesting, as the former is associated with nucleic acids [19,87], while the latter with the vibrational modes of proteins, but also with the in-plane vibrations of nitrogenous bases of RNA [9,18,88,89]. Therefore, an alteration in the levels of the mentioned protein and RNA markers is likely to affect the shape and the intensities of the analyzed mid-IR bands, producing detectable spectral changes (Figure 5). Additionally, we found highly significant differences in the shape of the lipid-ester peak at about 1730 cm−1 and in the lipid stretching band at 2800–3000 cm−1. This is an interesting finding and confirms that—aside from RNA and protein markers—the EV lipid content is an effective source of HCC markers, as supported by the fact that HCC development and progression have been linked to fatty acid metabolism dysregulation, in which aberrant activation of oncogenic signaling pathways alters the expression and activity of lipid-metabolizing enzymes [90].
Finally, in Figure 5 and Figure 6, we evaluated the classification performances of our spectral biomarkers alone and in combination using the LOOCV method and a ROC curve analysis followed by AUC calculation. The obtained results are interesting if compared with those obtained on conventional HCC circulating markers, i.e., AFP and PIVKA-II, in the same groups of subjects (Figure 6B,C). Despite that AFP is currently the most widely used serum HCC biomarker, its sensitivity at low cut-off values (less than 20 ng/mL) is approximately 60%, and the specificity is still inadequate [91,92]. Furthermore, serum AFP levels remain normal in 15–30% of advanced HCC [93], as, in our case, stressing the importance of finding other circulating markers capable of detecting AFP-negative cancer patients. In this context, PIVKA-II, Protein Induced by Vitamin K deficiency or antagonist-II (PIVKA-II), is believed to be a favorable biomarker to detect AFP-negative HCC [27,29,94,95]. Figure 6B,C show that our combined spectral biomarker possesses a larger, albeit comparable, AUC value than AFP and PIVKA-II measured on the same subjects. Very interestingly, if the sensitivity is computed from the ROC curves using the Youden’s method, we obtained a 70% sensitivity for AFP, 85% for PIVKA-II and 100% for our combined spectral biomarker. A similar outcome was obtained with the “closest top-left” method for computing sensitivity. Additionally, since no significant correlation was found between spectral markers and the levels of PIVKA-II and AFP (data not shown), it would be interesting to investigate whether the formers can contribute to the detection of HCC-negative cancer patients.
A limitation of this study is the reduced sample size, which is a common bottleneck affecting the vast majority of the published FTIR study on this subject. Therefore, the results here presented need to be validated on a larger sample size, which could serve also as a platform to test other machine learning algorithm than PCA-LDA, such as PLS-DA, which have been shown to outperform PCA-LDA on large datasets [96,97].

5. Conclusions

In this work, we present, for the first time in the literature, the proof-of-concept of a liquid biopsy approach for cancer diagnosis based on the molecular profiling of serum EVs with FTIR spectroscopy. A cohort of HCC patients of nonviral etiology was selected as a model system and compared to a cohort of noncancer subjects. As deeper argued in the discussion section, the alteration of the molecular content in EVs isolated from HCC patients has been reported by several excellent works. While the study of these biochemical changes at the single-molecule level is needed to have a better understanding of tumor biology, the combination of multiple circulating markers is highly demanded to improve diagnostic performances in a clinical setting. In this context, FTIR spectroscopy appears to be extremely promising as it can provide a label-free comprehensive molecular fingerprint of EVs through the analysis of specific mid-IR absorption bands. Here, we analyze the FTIR spectra of serum EVs by applying multivariate statistics and Machine Learning to highlight clinically relevant variations between the two groups. Our data point out the presence of statistically significant differences in the integrated intensities of major mid-IR absorption bands, including the carbohydrate and nucleic acids band, the protein amide I and II bands, and the lipid CH stretching band. These findings are efficiently captured by the lipid-to-protein spectral ratio (LPR, Figure 4E) and the Lipid-to-Nucleic Acids spectral ratio (LNR, Figure 4F), which show that HCC-derived EVs have an altered balance among the levels of different molecular classes compared to the controls. Additionally, we used PCA-LDA for the automated classification of spectral data according to the line shape of the mentioned mid-IR spectral bands. This analysis identified three potential spectral biomarkers for HCC diagnosis, namely the shape of the Amide I and II bands, the lipid CH stretching band, and the lipid-ester peak at about 1740 cm−1. The diagnostic performances of all the proposed spectral biomarkers, alone and in combination, were evaluated by ROC curve analysis. A stepwise multivariate logistic regression allowed us to select the most informative subset of spectral biomarkers suitable for patient classification. The combined biomarker displayed good classification potential (Figure 6), as confirmed by a large AUC value (0.91) and a sensitivity and specificity of 1 and 0.76, respectively. Very interestingly, our combined spectral biomarker appears to outperform the classification performance of two widely used circulating HCC markers—namely, AFP and PIVKA-II—measured on the same group of subjects, a result that deserves a dedicated study based on a larger sample size.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/jpm12060949/s1. Figure S1: (A) microscopic characterization of a representative purified sample performed with atomic force microscopy. The image is acquired on an EV-enriched suspension extracted from a pool of sera of 5 out of 19 controls recruited in the study. Several similar images were acquired, showing reproducible features. The image shows a large number of rounded particles with the expected size distribution for an EV sample measured with AFM on mica. Particles also display the typical cup shape morphology as shown by the selected (B) line profile. Additionally, smaller and partly aggregated particles can be observed, which could be due to the presence of contaminants, such as lipoproteins or albumin. (C) An Intensity-weighted diameter distribution measured with DLS. Measurements are acquired on 5 out of the 19 EV-enriched suspensions obtained from control subjects. Measures are reported as mean together with the 95% confidence bands. Two peaks can be observed, one cantered at approximately 32 nm and the other at approximately 190 nm. The first peak is compatible with the average size observed in the AFM image. Moreover, one can observe that microscopic images do not show particles with sizes comparable to the larger peak. This was expected as DLS Intensity data are weighted with the six power of the particle radius, which im-plies that approximately one over 60.000 particle is as great as 200 nm. Similar size comparison can be detected in the (D) TEM image (scale bar = 200 nm). TEM measurements were not acquired on the same samples recruited in the study, but they were purified using the same protocol, thus providing a further confirmation that despite contaminants, the kit allows a substantial enrichment in EVs.

Author Contributions

Conceptualization: G.C., R.D.S. and U.B.; Sample preparation: R.D.S., M.V., S.R., U.B. and L.M.; methodology: G.C., R.D.S., L.M., G.L.R. and U.B.; data curation: G.C., F.D.G. and R.D.S.; original draft preparation: G.C. and R.D.S.; review and editing: G.L.R., L.M., M.P. and M.D.S.; supervision: G.C., M.P., M.D.S., L.M., G.L.R. and U.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Italian Ministry of Health, “Progetto Giovani Ricercatori 2015-2016”, grant number GR-2016-02363310.

Institutional Review Board Statement

The study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Ethics Committee of Università Cattolica del Sacro Cuore (Protocol ID: 2078).

Informed Consent Statement

All participants provided written informed consent.

Data Availability Statement

The data presented in this study are available upon reasonable request to the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ilié, M.; Hofman, P. Pros: Can tissue biopsy be replaced by liquid biopsy? Transl. Lung Cancer Res. 2016, 5, 420. [Google Scholar] [CrossRef] [Green Version]
  2. Petrovčíková, E.; Vičíková, K.; Leksa, V. Extracellular vesicles—biogenesis, composition, function, uptake and therapeutic applications. Biologia 2018, 73, 437–448. [Google Scholar] [CrossRef]
  3. György, B.; Szabó, T.G.; Pásztói, M.; Pál, Z.; Misják, P.; Aradi, B.; László, V.; Pállinger, E.; Pap, E.; Kittel, A. Membrane vesicles, current state-of-the-art: Emerging role of extracellular vesicles. Cell. Mol. Life Sci. 2011, 68, 2667–2688. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Di Santo, R.; Romanò, S.; Mazzini, A.; Jovanović, S.; Nocca, G.; Campi, G.; Papi, M.; De Spirito, M.; Di Giacinto, F.; Ciasca, G. Recent Advances in the Label-Free Characterization of Exosomes for Cancer Liquid Biopsy: From Scattering and Spectroscopy to Nanoindentation and Nanodevices. Nanomaterials 2021, 11, 1476. [Google Scholar] [CrossRef] [PubMed]
  5. Bebelman, M.P.; Smit, M.J.; Pegtel, D.M.; Baglio, S.R. Biogenesis and function of extracellular vesicles in cancer. Pharmacol. Ther. 2018, 188, 1–11. [Google Scholar] [CrossRef] [PubMed]
  6. Akers, J.C.; Gonda, D.; Kim, R.; Carter, B.S.; Chen, C.C. Biogenesis of extracellular vesicles (EV): Exosomes, microvesicles, retrovirus-like vesicles, and apoptotic bodies. J. Neurooncol. 2013, 113, 1–11. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Margolis, L.; Sadovsky, Y. The biology of extracellular vesicles: The known unknowns. PLoS Biol. 2019, 17, 1–12. [Google Scholar] [CrossRef] [PubMed]
  8. Romanò, S.; Di Giacinto, F.; Primiano, A.; Gervasoni, J.; Mazzini, A.; Papi, M.; Urbani, A.; Serafino, A.; De Spirito, M.; Krasnowskad, E.K. Label-free spectroscopic characterization of exosomes reveals cancer cell differentiation. Anal. Chim. Acta 2021, 1192, 339359. [Google Scholar] [CrossRef] [PubMed]
  9. Lee, J.; Wen, B.; Carter, E.A.; Combes, V.; Grau, G.E.R.; Lay, P.A. Infrared spectroscopic characterization of monocytic microvesicles (microparticles) released upon lipopolysaccharide stimulation. FASEB J. 2017, 31, 2817–2827. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Fernandez, A.; João Canha Pereira, M.; Ramalhete, L.; Aleixo, S.; Lobato da Silva, C.; Cabral, J.; Calado, C.R. Evaluating the Impact of Culture Conditions on Human Mesenchymal Stem/Stromal Cell-derived Exosomes through FTIR Spectroscopy. J. Extracell. Vesicles 2020, 9, 366. [Google Scholar]
  11. Zlotogorski-Hurvitz, A.; Dekel, B.Z.; Malonek, D.; Yahalom, R.; Vered, M. FTIR-based spectrum of salivary exosomes coupled with computational-aided discriminating analysis in the diagnosis of oral cancer. J. Cancer Res. Clin. Oncol. 2019, 145, 685–694. [Google Scholar] [CrossRef]
  12. Perut, F.; Graziani, G.; Roncuzzi, L.; Zini, N.; Avnet, S.; Baldini, N. FT-IR Spectral Signature of Sensitive and Multidrug-Resistant Osteosarcoma Cell-Derived Extracellular Nanovesicles. Cells 2022, 11, 778. [Google Scholar] [CrossRef]
  13. Romanò, S.; Di Giacinto, F.; Primiano, A.; Mazzini, A.; Panzetta, C.; Papi, M.; Di Gaspare, A.; Ortolani, M.; Gervasoni, J.; De Spirito, M.; et al. Fourier Transform Infrared Spectroscopy as a useful tool for the automated classification of cancer cell-derived exosomes obtained under different culture conditions. Anal. Chim. Acta 2020, 1140, 219–227. [Google Scholar] [CrossRef] [PubMed]
  14. Yap, X.-L.; Ong, T.-A.; Lim, J.; Wood, B.; Lee, W.-L. Study of prostate cancer-derived extracellular vesicles in urine using IR spectroscopy. Prog. Drug Discov. Biomed. Sci. 2019, 2, 4–7. [Google Scholar] [CrossRef]
  15. Pascucci, L.; Coccè, V.; Bonomi, A.; Ami, D.; Ceccarelli, P.; Ciusani, E.; Viganò, L.; Locatelli, A.; Sisto, F.; Doglia, S.M.; et al. Paclitaxel is incorporated by mesenchymal stromal cells and released in exosomes that inhibit in vitro tumor growth: A new approach for drug delivery. J. Control. Release 2014, 192, 262–270. [Google Scholar] [CrossRef]
  16. Baddela, V.S.; Nayan, V.; Rani, P.; Onteru, S.K.; Singh, D. Physicochemical Biomolecular Insights into Buffalo Milk-Derived Nanovesicles. Appl. Biochem. Biotechnol. 2016, 178, 544–557. [Google Scholar] [CrossRef] [PubMed]
  17. Merdalimova, A.; Chernyshev, V.; Nozdriukhin, D.; Rudakovskaya, P.; Gorin, D.; Yashchenok, A. Identification and analysis of exosomes by surface-enhanced Raman spectroscopy. Appl. Sci. 2019, 9, 1135. [Google Scholar] [CrossRef] [Green Version]
  18. Mihály, J.; Deák, R.; Szigyártó, I.C.; Bóta, A.; Beke-Somfai, T.; Varga, Z. Characterization of extracellular vesicles by IR spectroscopy: Fast and simple classification based on amide and C[sbnd]H stretching vibrations. Biochim. Biophys. Acta Biomembr. 2017, 1859, 459–466. [Google Scholar] [CrossRef] [PubMed]
  19. Martins, T.S.; Magalhães, S.; Rosa, I.M.; Vogelgsang, J.; Wiltfang, J.; Delgadillo, I.; Catita, J.; da Cruz e Silva, O.A.B.; Nunes, A.; Henriques, A.G. Potential of FTIR Spectroscopy Applied to Exosomes for Alzheimer’s Disease Discrimination: A Pilot Study. J. Alzheimer’s Dis. 2020, 74, 391–405. [Google Scholar] [CrossRef]
  20. Uthamacumaran, A.; Elouatik, S.; Abdouh, M.; Berteau-Rainville, M.; Gao, Z.-H.; Arena, G. Machine Learning Characterization of Cancer Patients-Derived Extracellular Vesicles using Vibrational Spectroscopies. arXiv 2021, arXiv:2107.10332. [Google Scholar]
  21. Li, J.; Lee, Y.; Johansson, H.J.; Mäger, I.; Vader, P.; Nordin, J.Z.; Wiklander, O.P.B.; Lehtiö, J.; Wood, M.J.A.; El Andaloussi, S. Serum-free culture alters the quantity and protein composition of neuroblastoma-derived extracellular vesicles. J. Extracell. Vesicles 2015, 9, 26883. [Google Scholar] [CrossRef] [PubMed]
  22. Di Giacinto, F.; Tartaglione, L.; Nardini, M.; Mazzini, A.; Romanò, S.; Rizzo, G.E.; Papi, M.; De Spirito, M.; Pitocco, D.; Ciasca, G. Searching for the Mechanical Fingerprint of Pre-diabetes in T1DM: A Case Report Study. Front. Bioeng. Biotechnol. 2020, 8, 1099. [Google Scholar] [CrossRef] [PubMed]
  23. Stępień, E.; Kamińska, A.; Surman, M.; Karbowska, D.; Wróbel, A.; Przybyło, M. Fourier-Transform InfraRed (FT-IR) spectroscopy to show alterations in molecular composition of EV subpopulations from melanoma cell lines in different malignancy. Biochem. Biophys. Rep. 2021, 25, 100888. [Google Scholar] [CrossRef]
  24. Zhang, P.; Wang, L.; Fang, Y.; Zheng, D.; Lin, T.; Wang, H. Label-free exosomal detection and classification in rapid discriminating different cancer types based on specific Raman phenotypes and multivariate statistical analysis. Molecules 2019, 24, 2947. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Basile, U.; Miele, L.; Napodano, C.; Ciasca, G.; Gulli, F.; Pocino, K.; De Matthaeis, N.; Liguori, A.; De Magistris, A.; Marrone, G.; et al. The diagnostic performance of PIVKA-II in metabolic and viral hepatocellular carcinoma: A pilot study. Eur. Rev. Med. Pharmacol. Sci. 2020, 24, 12675–12685. [Google Scholar] [PubMed]
  26. Wan, H.-G.; Xu, H.; Gu, Y.-M.; Wang, H.; Xu, W.; Zu, M.-H. Comparison osteopontin vs AFP for the diagnosis of HCC: A meta-analysis. Clin. Res. Hepatol. Gastroenterol. 2014, 38, 706–714. [Google Scholar] [CrossRef]
  27. Ricco, G.; Cosma, C.; Bedogni, G.; Biasiolo, A.; Guarino, M.; Pontisso, P.; Morisco, F.; Oliveri, F.; Cavallone, D.; Bonino, F. Modeling the time-related fluctuations of AFP and PIVKA-II serum levels in patients with cirrhosis undergoing surveillance for hepatocellular carcinoma. Cancer Biomark. 2020, 29, 189–196. [Google Scholar] [CrossRef]
  28. Zhang, J.; Chen, G.; Zhang, P.; Zhang, J.; Li, X.; Gan, D.; Cao, X.; Han, M.; Du, H.; Ye, Y. The threshold of alpha-fetoprotein (AFP) for the diagnosis of hepatocellular carcinoma: A systematic review and meta-analysis. PLoS ONE 2020, 15, e0228857. [Google Scholar] [CrossRef] [PubMed]
  29. Yu, R.; Tan, Z.; Xiang, X.; Dan, Y.; Deng, G. Effectiveness of PIVKA-II in the detection of hepatocellular carcinoma based on real-world clinical data. BMC Cancer 2017, 17, 608. [Google Scholar] [CrossRef] [PubMed]
  30. Pelizzaro, F.; Cardin, R.; Penzo, B.; Pinto, E.; Vitale, A.; Cillo, U.; Russo, F.P.; Farinati, F. Liquid Biopsy in Hepatocellular Carcinoma: Where Are We Now? Cancers 2021, 13, 2274. [Google Scholar] [CrossRef] [PubMed]
  31. Chen, R.; Xu, X.; Tao, Y.; Qian, Z.; Yu, Y. Exosomes in hepatocellular carcinoma: A new horizon. Cell Commun. Signal. 2019, 17, 1. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Wang, H.; Lu, Z.; Zhao, X. Tumorigenesis, diagnosis, and therapeutic potential of exosomes in liver cancer. J. Hematol. Oncol. 2019, 12, 133. [Google Scholar] [CrossRef] [Green Version]
  33. Pocino, K.; Napodano, C.; Marino, M.; Di Santo, R.; Miele, L.; De Matthaeis, N.; Gulli, F.; Saporito, R.; Rapaccini, G.L.; Ciasca, G. A Comparative Study of Serum Angiogenic Biomarkers in Cirrhosis and Hepatocellular Carcinoma. Cancers 2022, 14, 11. [Google Scholar] [CrossRef]
  34. Van Deun, J.; Mestdagh, P.; Sormunen, R.; Cocquyt, V.; Vermaelen, K.; Vandesompele, J.; Bracke, M.; De Wever, O.; Hendrix, A. The impact of disparate isolation methods for extracellular vesicles on downstream RNA profiling. J. Extracell. Vesicles 2014, 3, 24858. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Huang, C.; Fisher, K.P.; Hammer, S.S.; Navitskaya, S.; Blanchard, G.J.; Busik, J. V Plasma exosomes contribute to microvascular damage in diabetic retinopathy by activating the classical complement pathway. Diabetes 2018, 67, 1639–1649. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Naumann, D. FT-Infrared and FT-Raman Spectroscopy in Biomedical Research. Appl. Spectrosc. Rev. 2001, 36, 239–298. [Google Scholar] [CrossRef]
  37. Kim, S.Y.; Khanal, D.; Tharkar, P.; Kalionis, B.; Chrzanowski, W. None of us is the same as all of us: Resolving the heterogeneity of extracellular vesicles using single-vesicle, nanoscale characterization with resonance enhanced atomic force microscope infrared spectroscopy (AFM-IR). Nanoscale Horiz. 2018, 3, 430–438. [Google Scholar] [CrossRef] [PubMed]
  38. Arrondo, J.L.R.; Goñi, F.M. Infrared Spectroscopic Studies of Membrane Lipids. In Biomolecular Structure and Dynamics; Springer: Dordrecht, The Netherlands, 1997; pp. 229–242. [Google Scholar]
  39. Wang, J.-C.; Bégin, L.R.; Bérubé, N.G.; Chevalier, S.; Aprikian, A.G.; Gourdeau, H.; Chevrette, M. Down-regulation of CD9 expression during prostate carcinoma progression is associated with CD9 mRNA modifications. Clin. Cancer Res. 2007, 13, 2354–2361. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  40. Peinado, H.; Alečković, M.; Lavotshkin, S.; Matei, I.; Costa-Silva, B.; Moreno-Bueno, G.; Hergueta-Redondo, M.; Williams, C.; García-Santos, G.; Ghajar, C.M. Melanoma exosomes educate bone marrow progenitor cells toward a pro-metastatic phenotype through MET. Nat. Med. 2012, 18, 883–891. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Jin, H.; Liu, P.; Wu, Y.; Meng, X.; Wu, M.; Han, J.; Tan, X. Exosomal zinc transporter ZIP4 promotes cancer growth and is a novel diagnostic biomarker for pancreatic cancer. Cancer Sci. 2018, 109, 2946–2956. [Google Scholar] [CrossRef] [Green Version]
  42. Sandfeld-Paulsen, B.; Jakobsen, K.R.; Bæk, R.; Folkersen, B.H.; Rasmussen, T.R.; Meldgaard, P.; Varming, K.; Jørgensen, M.M.; Sorensen, B.S. Exosomal proteins as diagnostic biomarkers in lung cancer. J. Thorac. Oncol. 2016, 11, 1701–1710. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Kimura, H.; Yamamoto, H.; Harada, T.; Fumoto, K.; Osugi, Y.; Sada, R.; Maehara, N.; Hikita, H.; Mori, S.; Eguchi, H. CKAP4, a DKK1 receptor, is a biomarker in exosomes derived from pancreatic cancer and a molecular target for therapy. Clin. Cancer Res. 2019, 25, 1936–1947. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Jafari, A.; Babajani, A.; Abdollahpour-Alitappeh, M.; Ahmadi, N.; Rezaei-Tavirani, M. Exosomes and cancer: From molecular mechanisms to clinical applications. Med. Oncol. 2021, 38, 45. [Google Scholar] [CrossRef] [PubMed]
  45. Greening, D.W.; Gopal, S.K.; Mathias, R.A.; Liu, L.; Sheng, J.; Zhu, H.-J.; Simpson, R.J. Emerging roles of exosomes during epithelial–mesenchymal transition and cancer progression. In Proceedings of the Seminars in Cell & Developmental Biology; Elsevier: Amsterdam, The Netherlands, 2015; Volume 40, pp. 60–71. [Google Scholar]
  46. Fang, T.; Lv, H.; Lv, G.; Li, T.; Wang, C.; Han, Q.; Yu, L.; Su, B.; Guo, L.; Huang, S. Tumor-derived exosomal miR-1247-3p induces cancer-associated fibroblast activation to foster lung metastasis of liver cancer. Nat. Commun. 2018, 9, 191. [Google Scholar] [CrossRef] [Green Version]
  47. Hsieh, C.-H.; Tai, S.-K.; Yang, M.-H. Snail-overexpressing cancer cells promote M2-like polarization of tumor-associated macrophages by delivering MiR-21-abundant exosomes. Neoplasia 2018, 20, 775–788. [Google Scholar] [CrossRef]
  48. Ludwig, N.; Whiteside, T.L. Potential roles of tumor-derived exosomes in angiogenesis. Expert Opin. Ther. Targets 2018, 22, 409–417. [Google Scholar] [CrossRef]
  49. Shao, C.; Yang, F.; Miao, S.; Liu, W.; Wang, C.; Shu, Y.; Shen, H. Role of hypoxia-induced exosomes in tumor biology. Mol. Cancer 2018, 17, 120. [Google Scholar] [CrossRef]
  50. Bruix, J. EASL Panel of Experts t on HCC. Clinical management of hepatocellular carcinoma. Conclusions of the Barcelona-2000 EASL conference. J. Hepatol. 2001, 35, 421–430. [Google Scholar] [CrossRef]
  51. Ciasca, G.; Sassun, T.E.; Minelli, E.; Antonelli, M.; Papi, M.; Santoro, A.; Giangaspero, F.; Delfini, R.; De Spirito, M. Nano-mechanical signature of brain tumours. Nanoscale 2016, 8, 19629–19643. [Google Scholar] [CrossRef]
  52. Minelli, E.; Ciasca, G.; Sassun, T.E.; Antonelli, M.; Palmieri, V.; Papi, M.; Maulucci, G.; Santoro, A.; Giangaspero, F.; Delfini, R.; et al. A fully-automated neural network analysis of AFM force-distance curves for cancer tissue diagnosis. Appl. Phys. Lett. 2017, 111, 143701. [Google Scholar] [CrossRef] [Green Version]
  53. Xiong, H.; Huang, Z.; Yang, Z.; Lin, Q.; Yang, B.; Fang, X.; Liu, B.; Chen, H.; Kong, J. Recent Progress in Detection and Profiling of Cancer Cell-Derived Exosomes. Small 2021, 17, 2007971. [Google Scholar] [CrossRef]
  54. Al-Sowayan, B.S.; Al-Shareeda, A.T. Nanogenomics and Artificial Intelligence: A Dynamic Duo for the Fight Against Breast Cancer. Front. Mol. Biosci. 2021, 8, 219. [Google Scholar] [CrossRef] [PubMed]
  55. Quagliarini, E.; Di Santo, R.; Pozzi, D.; Caracciolo, G. Protein corona-enabled serological tests for early stage cancer detection. Sensors Int. 2020, 1, 100025. [Google Scholar] [CrossRef]
  56. Saba, T. Recent advancement in cancer detection using machine learning: Systematic survey of decades, comparisons and challenges. J. Infect. Public Health 2020, 13, 1274–1289. [Google Scholar] [CrossRef] [PubMed]
  57. Caputo, D.; Digiacomo, L.; Cascone, C.; Pozzi, D.; Palchetti, S.; Di Santo, R.; Quagliarini, E.; Coppola, R.; Mahmoudi, M.; Caracciolo, G. Synergistic Analysis of Protein Corona and Haemoglobin Levels Detects Pancreatic Cancer. Cancers 2021, 13, 93. [Google Scholar] [CrossRef] [PubMed]
  58. Di Santo, R.; Quagliarini, E.; Digiacomo, L.; Pozzi, D.; Di Carlo, A.; Caputo, D.; Cerrato, A.; Montone, C.M.; Mahmoudi, M.; Caracciolo, G. Protein corona profile of graphene oxide allows detection of glioblastoma multiforme using a simple one-dimensional gel electrophoresis technique: A proof-of-concept study. Biomater. Sci. 2021, 9, 4671–4678. [Google Scholar] [CrossRef]
  59. Li, P.; Kaslan, M.; Lee, S.H.; Yao, J.; Gao, Z. Progress in exosome isolation techniques. Theranostics 2017, 7, 789–804. [Google Scholar] [CrossRef]
  60. Im, H.; Shao, H.; Park, Y.I.; Peterson, V.M.; Castro, C.M.; Weissleder, R.; Lee, H. Label-free detection and molecular profiling of exosomes with a nano-plasmonic sensor. Nat. Biotechnol. 2014, 32, 490–495. [Google Scholar] [CrossRef] [Green Version]
  61. Liga, A.; Vliegenthart, A.D.B.; Oosthuyzen, W.; Dear, J.W.; Kersaudy-Kerhoas, M. Exosome isolation: A microfluidic road-map. Lab. Chip 2015, 15, 2388–2394. [Google Scholar] [CrossRef] [Green Version]
  62. Lee, C.; Carney, R.P.; Hazari, S.; Smith, Z.J.; Knudson, A.; Robertson, C.S.; Lam, K.S.; Wachsmann-Hogiu, S. 3D plasmonic nanobowl platform for the study of exosomes in solution. Nanoscale 2015, 7, 9290–9297. [Google Scholar] [CrossRef]
  63. Tóth, E.Á.; Turiák, L.; Visnovitz, T.; Cserép, C.; Mázló, A.; Sódar, B.W.; Försönits, A.I.; Petővári, G.; Sebestyén, A.; Komlósi, Z. Formation of a protein corona on the surface of extracellular vesicles in blood plasma. J. Extracell. Vesicles 2021, 10, e12140. [Google Scholar] [CrossRef] [PubMed]
  64. Szabo, G.; Momen-Heravi, F. Extracellular vesicles in liver disease and potential as biomarkers and therapeutic targets. Nat. Rev. Gastroenterol. Hepatol. 2017, 14, 455–466. [Google Scholar] [CrossRef] [PubMed]
  65. Tang, J.; Li, Y.; Liu, K.; Zhu, Q.; Yang, W.-H.; Xiong, L.-K.; Guo, D.-L. Exosomal miR-9-3p suppresses HBGF-5 expression and is a functional biomarker in hepatocellular carcinoma. Minerva Med. 2017, 109, 15–23. [Google Scholar] [CrossRef] [PubMed]
  66. Meng, F.; Henson, R.; Wehbe–Janek, H.; Ghoshal, K.; Jacob, S.T.; Patel, T. MicroRNA-21 regulates expression of the PTEN tumor suppressor gene in human hepatocellular cancer. Gastroenterology 2007, 133, 647–658. [Google Scholar] [CrossRef] [Green Version]
  67. Xue, X.; Wang, X.; Zhao, Y.; Hu, R.; Qin, L. Exosomal miR-93 promotes proliferation and invasion in hepatocellular carcinoma by directly inhibiting TIMP2/TP53INP1/CDKN1A. Biochem. Biophys. Res. Commun. 2018, 502, 515–521. [Google Scholar] [CrossRef]
  68. Nakano, T.; Chen, I.; Wang, C.; Chen, P.; Tseng, H.; Huang, K.; Hu, T.; Li, L.; Goto, S.; Cheng, Y. Circulating exosomal miR-92b: Its role for cancer immunoediting and clinical value for prediction of posttransplant hepatocellular carcinoma recurrence. Am. J. Transplant. 2019, 19, 3250–3262. [Google Scholar] [CrossRef]
  69. Sugimachi, K.; Matsumura, T.; Hirata, H.; Uchi, R.; Ueda, M.; Ueo, H.; Shinden, Y.; Iguchi, T.; Eguchi, H.; Shirabe, K. Identification of a bona fide microRNA biomarker in serum exosomes that predicts hepatocellular carcinoma recurrence after liver transplantation. Br. J. Cancer 2015, 112, 532–538. [Google Scholar] [CrossRef]
  70. Zhong, Y.; Li, Y.; Song, T.; Zhang, D. MiR-718 mediates the indirect interaction between lncRNA SEMA3B-AS1 and PTEN to regulate the proliferation of hepatocellular carcinoma cells. Physiol. Genom. 2019, 51, 500–505. [Google Scholar] [CrossRef]
  71. Ahsani, Z.; Mohammadi-Yeganeh, S.; Kia, V.; Karimkhanloo, H.; Zarghami, N.; Paryan, M. WNT1 gene from WNT signaling pathway is a direct target of miR-122 in hepatocellular carcinoma. Appl. Biochem. Biotechnol. 2017, 181, 884–897. [Google Scholar] [CrossRef]
  72. Nassirpour, R.; Mehta, P.P.; Yin, M.-J. miR-122 regulates tumorigenesis in hepatocellular carcinoma by targeting AKT3. PLoS ONE 2013, 8, e79655. [Google Scholar] [CrossRef] [Green Version]
  73. Cheng, J.; Chen, Y.; Zhao, P.; Li, N.; Lu, J.; Li, J.; Liu, Z.; Lv, Y.; Huang, C. Dysregulation of miR-638 in hepatocellular carcinoma and its clinical significance. Oncol. Lett. 2017, 13, 3859–3865. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Zhang, Y.; Zhang, D.; Jiang, J.; Dong, L. Loss of miR-638 promotes invasion and epithelial-mesenchymal transition by targeting SOX2 in hepatocellular carcinoma. Oncol. Rep. 2017, 37, 323–332. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Wan, Y.; Cui, R.; Gu, J.; Zhang, X.; Xiang, X.; Liu, C.; Qu, K.; Lin, T. Identification of four oxidative stress-responsive microRNAs, miR-34a-5p, miR-1915-3p, miR-638, and miR-150-3p, in hepatocellular carcinoma. Oxid. Med. Cell. Longev. 2017, 2017, 5189138. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  76. Xue, X.; Zhao, Y.; Wang, X.; Qin, L.; Hu, R. Development and validation of serum exosomal microRNAs as diagnostic and prognostic biomarkers for hepatocellular carcinoma. J. Cell. Biochem. 2019, 120, 135–142. [Google Scholar] [CrossRef] [Green Version]
  77. Sun, L.; Su, Y.; Liu, X.; Xu, M.; Chen, X.; Zhu, Y.; Guo, Z.; Bai, T.; Dong, L.; Wei, C. Serum and exosome long non coding RNAs as potential biomarkers for hepatocellular carcinoma. J. Cancer 2018, 9, 2631. [Google Scholar] [CrossRef]
  78. Yang, F.U.; Zhang, L.; Huo, X.; Yuan, J.; Xu, D.; Yuan, S.; Zhu, N.; Zhou, W.; Yang, G.; Wang, Y. Long noncoding RNA high expression in hepatocellular carcinoma facilitates tumor growth through enhancer of zeste homolog 2 in humans. Hepatology 2011, 54, 1679–1689. [Google Scholar] [CrossRef]
  79. Zhang, C.; Yang, X.; Qi, Q.; Gao, Y.; Wei, Q.; Han, S. lncRNA-HEIH in serum and exosomes as a potential biomarker in the HCV-related hepatocellular carcinoma. Cancer Biomark. 2018, 21, 651–659. [Google Scholar] [CrossRef]
  80. Xu, H.; Chen, Y.; Dong, X.; Wang, X. Serum exosomal long noncoding RNAs ENSG00000258332. 1 and LINC00635 for the diagnosis and prognosis of hepatocellular carcinoma. Cancer Epidemiol. Prev. Biomark. 2018, 27, 710–716. [Google Scholar] [CrossRef] [Green Version]
  81. Lee, Y.R.; Kim, G.; Tak, W.Y.; Jang, S.Y.; Kweon, Y.O.; Park, J.G.; Lee, H.W.; Han, Y.S.; Chun, J.M.; Park, S.Y. Circulating exosomal noncoding RNAs as prognostic biomarkers in human hepatocellular carcinoma. Int. J. Cancer 2019, 144, 1444–1452. [Google Scholar] [CrossRef] [Green Version]
  82. Wang, G.; Liu, W.; Zou, Y.; Wang, G.; Deng, Y.; Luo, J.; Zhang, Y.; Li, H.; Zhang, Q.; Yang, Y. Three isoforms of exosomal circPTGR1 promote hepatocellular carcinoma metastasis via the miR449a–MET pathway. EBioMedicine 2019, 40, 432–445. [Google Scholar] [CrossRef] [Green Version]
  83. Zeng, X.; Chen, Q.; Gong, G.; Yuan, F.; Wang, T.; Zhang, G.; Li, X.; Wang, D.; Wang, Q. A simple formula to correct for the effects of storage time and temperature on the insulin concentration. J. Clin. Lab. Anal. 2020, 34, e23255. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  84. Arbelaiz, A.; Azkargorta, M.; Krawczyk, M.; Santos-Laso, A.; Lapitz, A.; Perugorria, M.J.; Erice, O.; Gonzalez, E.; Jimenez-Agüero, R.; Lacasta, A. Serum extracellular vesicles contain protein biomarkers for primary sclerosing cholangitis and cholangiocarcinoma. Hepatology 2017, 66, 1125–1143. [Google Scholar] [CrossRef] [PubMed]
  85. Tseng, H.-H.; Chang, J.-G.; Hwang, Y.-H.; Yeh, K.-T.; Chen, Y.-L.; Yu, H.-S. Expression of hepcidin and other iron-regulatory genes in human hepatocellular carcinoma and its clinical implications. J. Cancer Res. Clin. Oncol. 2009, 135, 1413–1420. [Google Scholar] [CrossRef] [PubMed]
  86. Fu, Q.; Zhang, Q.; Lou, Y.; Yang, J.; Nie, G.; Chen, Q.; Chen, Y.; Zhang, J.; Wang, J.; Wei, T. Primary tumor-derived exosomes facilitate metastasis by regulating adhesion of circulating tumor cells via SMAD3 in liver cancer. Oncogene 2018, 37, 6105–6118. [Google Scholar] [CrossRef] [Green Version]
  87. Rai, V.; Mukherjee, R.; Routray, A.; Ghosh, A.K.; Roy, S.; Ghosh, B.P.; Mandal, P.B.; Bose, S.; Chakraborty, C. Serum-based diagnostic prediction of oral submucous fibrosis using FTIR spectrometry. Spectrochim. Acta Part A Mol. Biomol. Spectrosc. 2018, 189, 322–329. [Google Scholar] [CrossRef]
  88. Movasaghi, Z.; Rehman, S.; Rehman, I.U. Fourier transform infrared (FTIR) spectroscopy of biological tissues. Appl. Spectrosc. Rev. 2008, 43, 134–179. [Google Scholar] [CrossRef]
  89. Baker, M.J.; Gazi, E.; Brown, M.D.; Shanks, J.H.; Gardner, P.; Clarke, N.W. FTIR-based spectroscopic analysis in the identification of clinically aggressive prostate cancer. Br. J. Cancer 2008, 99, 1859–1866. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  90. Wang, M.; Han, J.; Xing, H.; Zhang, H.; Li, Z.; Liang, L.; Li, C.; Dai, S.; Wu, M.; Shen, F. Dysregulated fatty acid metabolism in hepatocellular carcinoma. Hepatic Oncol. 2016, 3, 241–251. [Google Scholar] [CrossRef]
  91. Luo, P.; Wu, S.; Yu, Y.; Ming, X.; Li, S.; Zuo, X.; Tu, J. Current status and perspective biomarkers in AFP negative HCC: Towards screening for and diagnosing hepatocellular carcinoma at an earlier stage. Pathol. Oncol. Res. 2020, 26, 599–603. [Google Scholar] [CrossRef]
  92. Trevisani, F.; D’Intino, P.E.; Morselli-Labate, A.M.; Mazzella, G.; Accogli, E.; Caraceni, P.; Domenicali, M.; De Notariis, S.; Roda, E.; Bernardi, M. Serum α-fetoprotein for diagnosis of hepatocellular carcinoma in patients with chronic liver disease: Influence of HBsAg and anti-HCV status. J. Hepatol. 2001, 34, 570–575. [Google Scholar] [CrossRef]
  93. Han, L.-L.; Lv, Y.; Guo, H.; Ruan, Z.-P.; Nan, K.-J. Implications of biomarkers in human hepatocellular carcinoma pathogenesis and therapy. World J. Gastroenterol. WJG 2014, 20, 10249. [Google Scholar] [CrossRef] [PubMed]
  94. Wang, X.; Zhang, W.; Liu, Y.; Gong, W.; Sun, P.; Kong, X.; Yang, M.; Wang, Z. Diagnostic value of prothrombin induced by the absence of vitamin K or antagonist-II (PIVKA-II) for early stage HBV related hepatocellular carcinoma. Infect. Agent. Cancer 2017, 12, 47. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  95. Park, H.; Park, J.Y. Clinical significance of AFP and PIVKA-II responses for monitoring treatment outcomes and predicting prognosis in patients with hepatocellular carcinoma. Biomed. Res. Int. 2013, 2013, 310427. [Google Scholar] [CrossRef] [PubMed]
  96. Rashid, N.A.; Hussain, W.; Abdullah, A. Performance of classification analysis: A comparative, study between PLS-DA and integrating PCA+ LDA. Math Stat 2019, 7, 24–28. [Google Scholar] [CrossRef]
  97. Lee, L.C.; Liong, C.-Y.; Jemain, A.A. Partial least squares-discriminant analysis (PLS-DA) for classification of high-dimensional (HD) data: A review of contemporary practice strategies and knowledge gaps. Analyst 2018, 143, 3526–3539. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Schematic workflow of the presented EV-based liquid biopsy. Briefly, serum samples from HCC patients and controls were collected. EVs were isolated from serum samples and ATR-FTIR measurements were performed. Finally, automated classification of the two groups was achieved by machine learning analysis of the acquired spectra.
Figure 1. Schematic workflow of the presented EV-based liquid biopsy. Briefly, serum samples from HCC patients and controls were collected. EVs were isolated from serum samples and ATR-FTIR measurements were performed. Finally, automated classification of the two groups was achieved by machine learning analysis of the acquired spectra.
Jpm 12 00949 g001
Figure 2. Average absorbance spectra of control subjects (A) and HCC patients (B). Data are reported as mean ± 95% confidence bands. The analysis performed in this study is principally focused on the highlighted green (controls) and red (patients) bands.
Figure 2. Average absorbance spectra of control subjects (A) and HCC patients (B). Data are reported as mean ± 95% confidence bands. The analysis performed in this study is principally focused on the highlighted green (controls) and red (patients) bands.
Jpm 12 00949 g002
Figure 3. Enlarged details of characteristics absorption bands of biomolecules in EVs in control (green) and HCC patients (red), namely the carbohydrates and nucleic acids band at 1000–1200 cm−1 (A,B), the amide I and II bands at 1470–1700 cm−1 (C,D), the C=O stretching band at 1720–1740 cm−1 (E,F), and the lipid CH stretching band at 2800–3000 cm−1 (G,H).
Figure 3. Enlarged details of characteristics absorption bands of biomolecules in EVs in control (green) and HCC patients (red), namely the carbohydrates and nucleic acids band at 1000–1200 cm−1 (A,B), the amide I and II bands at 1470–1700 cm−1 (C,D), the C=O stretching band at 1720–1740 cm−1 (E,F), and the lipid CH stretching band at 2800–3000 cm−1 (G,H).
Jpm 12 00949 g003
Figure 4. Box-plot analysis of areas computed for the amide I and II bands at 1470–1700 cm−1 (A), the lipid CH stretching band at 2800–3000 cm−1 (B), the carbohydrates and nucleic acids band at 1000–1200 cm−1 (C), and the C=O stretching band at 1720–1740 cm−1 (D). The ratio between the integrated area computed in panels B and A (E) and B and C (F).
Figure 4. Box-plot analysis of areas computed for the amide I and II bands at 1470–1700 cm−1 (A), the lipid CH stretching band at 2800–3000 cm−1 (B), the carbohydrates and nucleic acids band at 1000–1200 cm−1 (C), and the C=O stretching band at 1720–1740 cm−1 (D). The ratio between the integrated area computed in panels B and A (E) and B and C (F).
Jpm 12 00949 g004
Figure 5. (A) Projections of the spectral data in the PC1–PC2 plane for the region enclosed between 1470–1700 cm−1, together with the LD1 line (black continuous line), which maximizes the separation between controls (green points) and HCC patients (red points). (B) LD1 discriminant scores from the region enclosed between 1470 and 1700 cm−1 are reported in the lower panel. The presence of statistically significant difference between the LD1 scores of the two groups is supported by a box-plot analysis (upper panel). (C) The PCA-LDA (1470–1700 cm−1) classification performances are reported in terms of sensitivity, sensibility, and accuracy. (DF) The same analyses were performed for the region enclosed between 2800–3000 cm−1 and (GI) the region enclosed between 1720–1760 cm−1.
Figure 5. (A) Projections of the spectral data in the PC1–PC2 plane for the region enclosed between 1470–1700 cm−1, together with the LD1 line (black continuous line), which maximizes the separation between controls (green points) and HCC patients (red points). (B) LD1 discriminant scores from the region enclosed between 1470 and 1700 cm−1 are reported in the lower panel. The presence of statistically significant difference between the LD1 scores of the two groups is supported by a box-plot analysis (upper panel). (C) The PCA-LDA (1470–1700 cm−1) classification performances are reported in terms of sensitivity, sensibility, and accuracy. (DF) The same analyses were performed for the region enclosed between 2800–3000 cm−1 and (GI) the region enclosed between 1720–1760 cm−1.
Jpm 12 00949 g005
Figure 6. ROC analysis of the all the analyzed spectral biomarkers (A) and of AFP (gold), PIVKA-II (violet), and of the combined spectral biomarkers obtained with stepwise logistic regression (red) (B); AUC values for each of the investigated markers, alone and in combination, together with the corresponding 95% confidence interval (C). The same color scale is used in the three panels.
Figure 6. ROC analysis of the all the analyzed spectral biomarkers (A) and of AFP (gold), PIVKA-II (violet), and of the combined spectral biomarkers obtained with stepwise logistic regression (red) (B); AUC values for each of the investigated markers, alone and in combination, together with the corresponding 95% confidence interval (C). The same color scale is used in the three panels.
Jpm 12 00949 g006
Table 1. Demographical, clinical, biochemical, and spectral characteristics of the recruited subjects.
Table 1. Demographical, clinical, biochemical, and spectral characteristics of the recruited subjects.
CharacteristicCTRL, N = 19 1HCC, N = 20 1p-Value 2
Age (years)65.8 (4.8)68.7 (5.9)0.14
Gender 0.7
F26%20%
M74%80%
Cholesterol (mg/dL)180 (19)170 (46)0.5
Triglycerides (mg/dL)122 (38)101 (30)0.023
GPT-ALT (UI/L)15 (3)45 (23)<0.001
GOT-AST (UI/L)20 (5)47 (21)<0.001
AST/ALT1.43 (0.53)3.86 (11.37)0.3
GGT (UI/L)31 (5)179 (147)<0.001
ALP (UI/L)171 (17)272 (137)0.049
Hb (mmol/L)14.32 (0.96)12.78 (2.78)0.15
Creatinine (mg/dL)0.79 (0.18)1.42 (1.29)0.071
Azotemia (mg/dL)15 (3)21 (10)0.067
Bilirubin (mg/dL)0.69 (0.23)1.63 (1.04)<0.001
Area (1470–1700 cm−1)12.0 (3.1)9.7 (3.6)0.033
Area (1000–1200 cm−1)4.16 (1.22)3.32 (0.99)0.033
Area (2800–3000 cm−1)0.023 (0.006)0.030 (0.008)0.006
Area (1720–1760 cm−1)0.0012 (0.0007)0.0017 (0.0008)0.041
LD1 (1470–1700 cm−1)−0.54 (0.77)0.52 (1.18)0.002
LD1 (2800–3000 cm−1)−0.48 (1.24)0.46 (0.70)0.0075
LD1 (1720–1760 cm−1)−0.77 (1.10)0.73 (0.90)<0.001
PIVKA14 (4)3205 (9987)<0.001
AFP2 (1)193 (417)0.001
1 Mean (SD); %. 2 Wilcoxon rank–sum test; Fisher’s exact test.
Table 2. ROC-AUC Table analysis of the significant spectral biomarkers highlighted in Figure 4 and Figure 5. Conventional HCC circulating biomarkers—namely AFP and PIVKA-II—are also included in the analysis together with the results of a stepwise combined model selecting the most informative biomarkers for the classification of the two groups.
Table 2. ROC-AUC Table analysis of the significant spectral biomarkers highlighted in Figure 4 and Figure 5. Conventional HCC circulating biomarkers—namely AFP and PIVKA-II—are also included in the analysis together with the results of a stepwise combined model selecting the most informative biomarkers for the classification of the two groups.
VariableAUC95% CI
Area (1470–1700 cm−1)0.7000.524–0.876
Area (1000–1200 cm−1)0.7000.534–0.866
Area (2800–3000 cm−1)0.7550.596–0.915
Area (1720–1760 cm−1)0.6920.518–0.866
LD1 (1470–1700 cm−1)0.7760.627–0.926
LD1 (2800–3000 cm−1)0.7080.534–0.881
LD1 (1720–1760 cm−1)0.8260.694–0.959
PIVKA–II0.8050.649–0.961
AFP0.8620.728–0.996
Stepwise model *0.9100.803–1
* LD1 (1720–1760 cm−1) + Area (2800–3000 cm−1).
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Di Santo, R.; Vaccaro, M.; Romanò, S.; Di Giacinto, F.; Papi, M.; Rapaccini, G.L.; De Spirito, M.; Miele, L.; Basile, U.; Ciasca, G. Machine Learning-Assisted FTIR Analysis of Circulating Extracellular Vesicles for Cancer Liquid Biopsy. J. Pers. Med. 2022, 12, 949. https://doi.org/10.3390/jpm12060949

AMA Style

Di Santo R, Vaccaro M, Romanò S, Di Giacinto F, Papi M, Rapaccini GL, De Spirito M, Miele L, Basile U, Ciasca G. Machine Learning-Assisted FTIR Analysis of Circulating Extracellular Vesicles for Cancer Liquid Biopsy. Journal of Personalized Medicine. 2022; 12(6):949. https://doi.org/10.3390/jpm12060949

Chicago/Turabian Style

Di Santo, Riccardo, Maria Vaccaro, Sabrina Romanò, Flavio Di Giacinto, Massimiliano Papi, Gian Ludovico Rapaccini, Marco De Spirito, Luca Miele, Umberto Basile, and Gabriele Ciasca. 2022. "Machine Learning-Assisted FTIR Analysis of Circulating Extracellular Vesicles for Cancer Liquid Biopsy" Journal of Personalized Medicine 12, no. 6: 949. https://doi.org/10.3390/jpm12060949

APA Style

Di Santo, R., Vaccaro, M., Romanò, S., Di Giacinto, F., Papi, M., Rapaccini, G. L., De Spirito, M., Miele, L., Basile, U., & Ciasca, G. (2022). Machine Learning-Assisted FTIR Analysis of Circulating Extracellular Vesicles for Cancer Liquid Biopsy. Journal of Personalized Medicine, 12(6), 949. https://doi.org/10.3390/jpm12060949

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

Article Metrics

Back to TopTop