Tomas Bertok1Aniko Bertokova1Eduard Jane1Michal Hires1Juvissan Aguedo1Maria Potocarova2
Ludovit Lukac2Alica Vikartovska1
Peter Kasak3
Lubor Borsig 4,5and 
Jan Tkac1*- 1Institute of Chemistry, Slovak Academy of Sciences, Bratislava, Slovakia
- 2University Hospital Bratislava, Bratislava, Slovakia
- 3Center for Advanced Materials, Qatar University, Doha, Qatar
- 4Department of Physiology, University of Zurich, Zurich, Switzerland
- 5Comprehensive Cancer Center, Zurich, Switzerland
Colorectal cancer (CRC) is one of the most common types of cancer among men and women worldwide. Efforts are currently underway to find novel and more cancer-specific biomarkers that could be detected in a non-invasive way. The analysis of aberrant glycosylation of serum glycoproteins is a way to discover novel diagnostic and prognostic CRC biomarkers. The present study investigated a whole-serum glycome with a panel of 16 different lectins in search for age-independent and CRC-specific glycomarkers using receiver operating characteristic (ROC) curve analyses and glycan heat matrices. Glycosylation changes present in the whole serum were identified, which could lead to the discovery of novel biomarkers for CRC diagnostics. In particular, the change in the bisecting glycans (recognized by Phaseolus vulgaris erythroagglutinin) had the highest discrimination potential for CRC diagnostics in combination with human L selectin providing area under the ROC curve (AUC) of 0.989 (95% CI 0.950–1.000), specificity of 1.000, sensitivity of 0.900, and accuracy of 0.960. We also implemented novel tools for identification of lectins with strong discrimination power.
Graphical Abstract
Introduction
Cancer is regarded as the most important obstacle to increasing life expectancy in the 21st century. In both sexes, colorectal cancer (CRC) is one of the leading cancers, with incidence of ~1.8 M and mortality of 860,000 people annually (1), along with lung and prostate/breast cancers (2, 3). Moreover, it is estimated that, by 2030, there will be 2.2 M new CRC cases with an associated 1.1 M deaths (4, 5). Incidence increases rapidly with higher age and western lifestyle, while early diagnosis of localized tumors using sigmoidoscopy/colonoscopy with polypectomy increases the 5-year survival rate up to 90% (6). Point mutations in specific genes result in sporadic CRC in around 70% of all cases, while there is also a familial and inherited aspect behind CRC development and progression (7). A cancer biomarker is defined as a tumor characteristic that can be objectively measured and associated with the pathogenic process. Biomarkers are used for diagnostic, prognostic, predictive, and therapeutic purposes (8). In some cases, risk/predisposition biomarkers are used to identify people at significant risk of developing a disease. In addition to tissue biomarkers, such as cytokeratins, cadherin 17, or upregulation of telomerase expression (9), stool and blood biomarkers are of importance due to the less invasive collection of the sample.
CRC screening is performed using colonoscopy, flexible sigmoidoscopy, fecal occult blood testing, fecal DNA testing and by measuring the blood level of the carcinoembryonic antigen (CEA) (1, 10). Out of all these methods, colonoscopy is the gold standard due to the high sensitivity of the method, but the examination is uncomfortable for the patient and the accuracy of the examination depends on the skill level and experience of the operator (1). Although CEA is used for CRC diagnostics in combination with other examination methods, due to its low sensitivity and specificity, this biomarker is usually used in monitoring the effectiveness of the treatment and CRC recurrence (1). Additionally, CEA is also produced by other tissues in the human body, such as the breasts and pancreas (11), or is elevated under normal conditions, e.g., in smokers (12, 13). Moreover, besides CEA, other CRC biomarkers could be used in clinical practice with drawbacks similar to that of CEA (10, 11) and thus cannot be used in CRC screening but rather in monitoring the disease progression or CRC prognosis (9). Hence, there is a continuing need for biomarkers for early-stage CRC diagnosis, which can be detected in blood, serum, or plasma (liquid biopsy) (14, 15). Liquid biopsy-based assays were not yet recommended for CRC screening in the 2018 Guideline Updates from the American Cancer Society (16), but such tests may have a greater role in CRC screening in the future (16).
Very recently, novel blood-based biomarkers have been discovered, such as circulating tumor cells, circulating cell-free DNA, noncoding RNAs and microRNAs, and extracellular vesicles (mainly exosomes and oncosomes) (17). Besides several alternative CRC biomarkers, exosomes are regarded as a rich source of CRC biomarkers (proteins, glycoproteins, and various forms of RNAs), which can be effectively used for CRC diagnostics and CRC prognosis, as summarized in recent review papers (1, 18). Proteins expressed by exosomes were proven effective for the identification of CRC patients resistant to drug treatment or CRC patients with metastases (19).
A novel method for the discovery of CRC biomarkers is to identify glycosylation changes associated with CRC. Below, we describe the application of reverse-phase lectin microarrays for the identification of changes in the serum glycome as potential CRC biomarkers. A fluorescent microarray in combination with lectins is frequently used in the discovery of novel glycan-based biomarkers (20–24). In order to identify only those glycosylation changes that are associated with CRC and not with aging, serum samples in this study were analyzed as two distinct matrices, i.e., as the age matrix [healthy young (hY) vs. healthy old (hO)] and the CRC matrix (hO vs. CRC). Only those changes that are associated with the CRC matrix and not with the age matrix are deemed to be prospective novel CRC biomarkers.
Experimental
Serum samples were taken from 34 individuals (6 healthy individuals with no CRC confirmed and with no comorbidities diagnosed at that time with an average age of 33.0 ± 6.1 years (hY cohort), 10 individuals with no malignancy with an average age of 67.0 ± 8.6 years (hO cohort), and 18 colorectal cancer patients with histologically proven CRC with an average age of 73.0 ± 7.3 (CRC cohort)). The Ethics Committee of the Faculty of Medicine, Comenius University and University Hospital in Bratislava, Old Town Hospital, Bratislava, Slovakia, approved the use of the samples, and all participants signed an informed consent document prior to sample collection. The procedure was performed under the ethical guidelines of the last revision of the Helsinki Declaration. Untreated serum samples were taken during the morning fasted state using a gel and clot activator tube (Vacutest Kima, Piove di Sacco, Italy). After 30 min, the tubes were centrifuged at 25°C for 10 min at 2,500 g. The sera were transferred into sterile plastic vials and were stored in the form of aliquots at -80°C until use and used within 1 year.
Chemicals
All common chemicals [e.g., buffer components, bovine serum albumin (BSA), etc.] were purchased from Sigma Aldrich (USA). All solutions were freshly prepared prior to experiments in 0.055 μS deionized water (DW) and filtered using 0.2-μm sterile filters. Biotin conjugation kits for the biotinylation of unconjugated lectins were purchased from Abcam (UK). Lectins RPL-Fuc1 and RPL-Sia2 were obtained from GlycoSelect (Ireland). HPyL was purchased from GlycoDiag (France). P-selectin (P sel), L-selectin (HL sel), and E-selectin (HE sel) were used in the form of chimera proteins fused to IgG1 tail and obtained from Prof. Borsig with details provided in the paper (25). All the other lectins used were purchased in their biotinylated form from Vector Labs (USA). Conjugate streptavidin-CF647 was provided by Biotium (USA).
Reverse-Phase Lectin Microarrays
Reverse-phase lectin microarray experiments were performed with a phosphate buffer solution (PBS) (0.01 M, pH 7.4) as a printing buffer. Spotting temperature was set at 10°C and humidity at 60%. Subsequently, the slide was placed in a humidity chamber for 1 h at ambient temperature (AT) with humidity of 80%–90%, blocked using a blocking buffer at ambient temperature for 1 h and with slow shaking, rinsed under a gentle stream of printing buffer in a Petri dish, and drained. For blocking purposes, 70 μl of 3% bovine serum albumin (BSA) was used, as we observed lower background fluorescence intensity than with a Carbo-free blocking solution (Vector Labs, USA) for some of the lectins used. Samples diluted 50× were spotted in two different wells in triplicates using SpotBot3 Microarray Protein edition (Arrayit, USA) on epoxide-coated slides Nexterion E (Schott, Germany) using a previously optimized protocol. Subsequently, after spotting and blocking the slides (1 h at AT with shaking), 70 μl of biotinylated lectins (c = 5 μg/ml in PBS) was added and incubated at AT for 1 h. The slides were washed gently three times with PBS, and then 70 μl of streptavidin-CF647 conjugate (c = 0.1 μg/ml in PBS) was added for 15 min. After a washing step and additional wash with DW, fluorescence intensity was read at 635 nm using an InnoScan microarray reader (Arrayit, USA). The signal evaluated and ascribed to individual samples using Mapix software was an average value of at least three spots after background fluorescence subtraction.
Data Evaluation
All computations were performed using R software (version 3.6.3) (26) with CARET (Classification and Regression Training) and GLM packages (27). Receiver operating characteristic (ROC) curve, sensitivities, and specificities were estimated using the pROC package (28). The confidence intervals for area under the ROC curve (AUC) were computed by the bootstrap method with 2,000 stratified bootstrap replicates (29). Glycan heat matrices were prepared using OriginPro 2020.
Net reclassification improvement (NRI) and integrated discrimination improvement (IDI) were calculated according to Pencina et al. (30) using R package Hmics (https://CRAN.R-project.org/package=Hmisc). We computed continuous [category-free, NRI (>0)] NRI according to Pencina et al. (31). According to Pencina et al. (32), if NRI (>0) >0.6, there is strong discrimination, if NRI (>0) ~0.4, it is intermediate, and for NRI (>0) <0.2, discrimination between the models is considered weak. The heat maps are the NRI (>0) values. Only a combination of two lectins was considered in the evaluation. In the NRI (>0) and IDI evaluation, two models were compared. In our case, the old model of a single lectin was used and the updated model is the combination of two lectins.
The multicollinearity was tested using variance inflation factor (VIF) using R package car (33) with a VIF value above 10 indicating a multicollinearity problem (34).
Results and Discussion
Numerous studies have described glycosylation changes to be strongly associated with age (35–38). To clearly identify which glycan changes are strongly associated with age, we evaluated two healthy cohorts, i.e., hY (33.0 ± 6.1 years old) and hO (67.0 ± 8.6 years old). This discrimination is provided in the form of column graphs (Figure 1) for single lectins and a heat map as a submatrix 1 (hY vs. hO) in Figure 2 for the combination of two lectins.
Figure 1 Values of AUC for discrimination of the age matrix [healthy young (hY) vs. healthy old (hO) individuals] and the CRC matrix [hO vs. CRC patients] using single lectins (A). The ratio of AUCage vs. AUCCRC as determined for single lectins with the ratio below 0.8 is shown in green and the ratio above 1.2 is shown in brown (B). Standard deviation (SD) values are not shown for the better visual clarity of both figures.
Figure 2 Heat map generated using 16 lectins and their performance according to AUC value for two submatrices, i.e., hY vs. hO (left upper triangle) and hO vs. CRC (right lower triangle).
To evaluate glycosylation changes associated only with CRC, we also evaluated the CRC matrix, i.e., hO (67.0 ± 8.6 years old) vs. CRC patients (73.0 ± 7.3 years old). The results are shown in Figure 1 for single lectins and Figure 2 as a heat matrix for double lectins.
Single Lectins
Figure 1 shows the AUC values obtained when discrimination hY vs. hO (the age matrix) and hO vs. CRC (the CRC matrix) was investigated. More details of the clinical performance of single lectins to discriminate the age matrix vs. the CRC matrix are given in Table S1.
Weak Discrimination Power for Age Matrix vs. Colorectal Cancer Matrix
Discrimination hY vs. hO (the age matrix) and hO vs. CRC (the CRC matrix) for single lectins is shown in Figures 1A, B. The ratio AUCage/AUCCRC was introduced here to investigate whether a single lectin has a discrimination power for the age matrix (high AUCage/AUCCRC with a value higher than 1.2), a discrimination power for the CRC matrix (low AUCage/AUCCRC with a value lower than 0.8), or only a moderate discrimination potential (AUCage/AUCCRC within values 0.8–1.2). There are several lectins (i.e., AAL, DBA, MAA, PHAL, WGA, RPL-Fuc1, RPL-Sia2, P sel, and HE sel) with a really minor difference in the ratio AUCage/AUCCRC (i.e., in the range of 0.8–1.2) (Table S1 and Figure 1B). This means that such lectins can discriminate hY vs. hO and hO vs. CRC with similar performance, i.e., similar AUC. Hence, such lectins have only limited potential for CRC diagnostics.
Strong Discrimination Power for the Age Matrix
There are three lectins that can be used in selective discrimination of the age matrix over the CRC matrix, such as Concanavalin A (ConA) (AUCage/AUCCRC = 1.76); WFL (AUCage/AUCCRC = 1.32), and SNA I (AUCage/AUCCRC = 1.21) (Figure 1B and Table S1). Serum proteins are present in the blood in the following descending order: IgG (40.4%), IgA (9.0%), serotransferrin (8.5%), IgM (5.0%), haptoglobin (4.5%), α-1-acid glycoprotein (2.6%), ceruloplasmin (1.2%), β-2-glycoprotein I (0.7%), apolipoprotein B-100 (0.5%), apolipoprotein D (0.3%), and IgD (0.1%) (39). Hence, the most abundant protein in the blood is IgG, which might be the most likely carrier of glycans determined in this study. Several studies describe that, at early adulthood, there is a high abundance of digalactosylated IgG forms and that, with increasing age, a decrease in galactosylation and sialylation can be observed (37). Thus, the SNA I lectin should be a positive discriminator when used for the age matrix (i.e., hY vs. hO), as was experimentally confirmed in this study (Figure 1A). In the other study using the reverse-phase lectin microarray, it was found that the SNA I binding to transferrin isolated from human serum decreased with age (40), suggesting that there might be other protein carriers besides IgG-carrying glycans recognized by SNA I.
Another lectin, which is more suitable for discrimination of the age matrix rather than for discrimination of the CRC matrix, is WFL lectin. Since WFL especially recognizes GalNAc and LacdiNAc glycan structures, this might indicate that such glycan structures are present on other proteins in the blood and not on IgG, which has a rather conserved biantennary structure.
The significantly increased AUC value for ConA to discriminate the age matrix might indicate a decrease in the overall glycosylation of proteins present in the serum with age or a decrease in oligo-mannose-containing glycans. Changes in the glycosylation of proteins can also be an indicator of other processes in the body, including inflammation and other diseases (37).
There is one negative correlation for discriminating hY vs. hO, i.e., application of HPyL. This lectin is produced by GlycoDiag with the source of Human Polyomavirus 9 VP1 (Table 1). This lectin preferentially binds short, linear glycan sequences terminating in N-glycolylneuraminic acid (Neu5Gc) (46). Neu5Gc is a non-human derivative of sialic acid delivered to the human body in the form of red meat. The only explanation for the result obtained using HPyL is that the serum of older healthy individuals contains a larger amount of Neu5Gc as a result of a higher intake of red meat for elderly people.
Table 1 Lectin specificity for the lectins applied in this study.
Another glycan change referred to in the literature is a decrease in the level of Neu5Ac attached to Gal via α2,3-linkage (38) with age, which was also observed in this work with an AUC value of 0.767 for the MAA lectin for the age matrix.
Strong Discrimination Power for the Colorectal Cancer Matrix
There are, however, several lectins with a good discrimination power for the CRC matrix rather than for the age matrix, which is a feature applicable to CRC diagnostics, i.e., PHAE (AUCage/AUCCRC = 0.66), HPyL (AUCage/AUCCRC = 0.68), HL sel (AUCage/AUCCRC = 0.68), and RCA I (AUCage/AUCCRC = 0.74).
Three out of those four lectins have AUC exceeding a value of 0.8, i.e., PHAE, RCA I, and HL sel, with a potential application for CRC diagnostics.
Changes in the glycosylation of proteins in the sera of CRC patients were recently extensively reviewed (47). The CRC patients were observed to have a decreased occurrence in bisecting N-acetyl-lactosamine (LacNAc) structures and an increase in truncated paucimannosidic (high-mannose) structures for N-glycans and a decrease in Core 3 and Core 4 structures for O-glycans (47). Moreover, a decrease in extended O-glycans was observed with CRC development and progression (48).
PHAE lectin specifically recognizes bisecting glycans (49), and the present study showed that CRC patients possessed a lower amount of bisecting N-glycans, which is in agreement with previous conclusions (47). Moreover, a lectin microarray study focused on the analysis of glycosylation changes of transferrin isolated from human serum showed a decreased response toward transferrin’s glycans recognized by PHAE (40). This might indicate that transferrin could also be a potential carrier of the glycans recognized by PHAE in our study, since transferrin is an abundant protein in the serum (8.5% of total serum proteins).
The RCA I-recognizing lectin indicates that Gal and Galβ4GlcNAc-containing glycans are present at a lower level in CRC patients in comparison with healthy individuals. This finding is in agreement with the results from a lectin microarray study focused on a change of transferrin’s glycans related to CRC (40). This might indicate transferrin to be a potential carrier of glycans recognized by RCA I.
It is remarkable that HPyL-recognizing Neu5Gc in our study has a good negative discrimination power (i.e., increased level for healthy individuals in comparison with the CRC patients) (Figure 1A), since there is quite a strong correlation between red meat uptake and CRC (50, 51). It is well known that red meat contains sialic acid in the form of Neu5Gc, and it was recently shown that the level of Neu5Gc is higher in CRC patients than in healthy individuals (52). However, there is an explanation for this apparent discrepancy. We found an increased level of Neu5Gc in the healthy old cohort in comparison with the healthy young cohort. After CRC diagnosis, such CRC patients might be on a low red meat diet, which might result in a reduced amount of Neu5Gc in the CRC patients.
From all three selectins used in the study, it was found that HL sel had a strong discrimination potential for the CRC matrix (Figure 1A). A lectin-binding preference (Table 1) shows that HL sel has quite a strong binding preference for sulfated glycans. The results show a greater abundance of sulfated glycans in the healthy old cohort than in the CRC patient cohort. It is not yet clear what kind of glycosylation change to expect regarding the acetylation and sulfation of glycans during CRC development and/or progression (47).
The results related to discrimination of the ConA lectin are in agreement with the literature (40, 47) showing an increased level of mannosylated (truncated paucimannosidic) N-glycan structures associated with CRC.
Combination of Two Lectins
Evaluation of double biomarkers for discrimination of the age matrix and the CRC matrix resulted in 120 combinations for each matrix, i.e., in total 240 combinations (Table S1 and Figure 2). Accordingly, in the following section, we discuss only the double lectin combinations with lectins identified in the previous sections as lectins showing the best discrimination performance for the particular matrix.
Since SNA I showed AUC above 0.8 (i.e., 0.95) and the ratio of AUCage/AUCCRC above 1.2 (i.e., 1.21), only the double lectin biomarkers with SNA I for the age matrix will be evaluated here. As for the CRC matrix, only the double lectin combinations with PHAE (AUC = 0.939, AUCage/AUCCRC = 0.657), RCA I (AUC = 0.811; AUCage/AUCCRC = 0.74), or HL sel (AUC = 0.889, AUCage/AUCCRC = 0.675) will be evaluated.
Strong Discrimination Power for the Age Matrix
From all the 15 combinations of double lectin biomarkers with SNA I in Table 2, there is a selection of only those combinations with an AUC value that exceeds that of the AUC value of the single SNA I lectin, i.e., the AUC value of 0.95. Thus, in Table 2, clinical performance parameters are shown for two such combinations having a significant power to discriminate between the age matrix over the CRC matrix, as judging from the AUCage/AUCCRC ratio above 1.2, i.e., SNA I+RPL-Fuc1 (1.319) and SNA I+WFL (1.285). This means that two such combinations of lectins have the potential to be used in the selective discrimination of age. A literature review identified a paper where the ratio of two biantennary glycans present in the serum, i.e., G0F/G2F (G0F denotes an agalactosylated bianntennary glycan with core fucose, and G2F denotes a digalactosylated biannennary glycan with core fucose) increased with age (53). This ratio was then applied as a GlycoAgeTest for the detection of biological age (53). The other study confirmed a decrease in G2F and an increase in G0F of IgG (54) in agreement with the previous study (53). Moreover, Pučić et al. (54) identified a decrease in the IgG sialylation with aging. The results in the present study indicate that SNA I can detect the change in the sialylation of N-glycans associated with age (Figure 1), which confirms the observation in the study by Pučić et al. (54), and that most probably the protein carrier associated with such glycosylation changes is IgG. For example, changes in glycosylation including galactosylation, sialylation, and bisecting glycans were used in the calculation of the GlycanAge index (55).
Table 2 Clinical performance characteristics of double lectins for the age matrix with the best combination showed in red.
Our results indicate that a combination of SNA I with two other lectins, i.e., RPL-Fuc1 and WFL, has the potential to effectively discriminate age, in particular a decreased fucosylation (recognized by RPL-Fuc1 but not by AAL) and a decrease in the level of GalNAc- and LacdiNAc-containing glycans (recognized by WFL). It is obvious that there is an age-related change in the glycosylation of IgG and other serological proteins in the process called inflammaging (chronic and low-grade inflammation progressing with age) (56). Table 2 shows that a combination of lectins SNA I + WFL displayed the highest accuracy in discriminating the age matrix, as was confirmed by an additional statistical evaluation (Figure 3 left; Figure 4 left). Furthermore, the statistical evaluation showed CIF well below a threshold value of 10 for each lectin combination including SNA I + WFL (Figure 5 left), not indicating a multicollinearity problem. The only exception is the combination of SNA I with HL sel with a VUF value of 30.4, which can be explained by the fact that both lectins recognize the same glycan epitope (i.e., α2,6-linked sialic acid).
Figure 3 Heat maps for net reclassification improvement (NRI) (>0) analysis showing hY vs. hO (left) and hO vs. CRC (right).
Figure 4 Heat maps for integrated discrimination improvement (IDI) analysis showing hY vs. hO (left) and hO vs. CRC (right).
Figure 5 Heat maps for variance inflation factor (VIF) analysis showing hY vs. hO (left) and hO vs. CRC (right).
Figure 3 left (hY vs. hO) shows a strong intensity in the horizontal for lectin SNA and slightly less intensive for lectin WGA. These lectins provide strong discrimination in case of hY vs. hO (only these two lectins have AUC over 0.8 in single lectin evaluation as shown in Figure 1). Figure 3 right (hO vs. CRC) shows a high intensity for the lectins DBA, PHAE, RCA I, WGA, and selectins, which reached AUC above the value of 0.8 in a single lectin evaluation. Lectin PHA-E has the highest AUC, and PHAE lectin also shows a horizontal row with the highest values of NRI. Clearly, the NRI values confirm the results presented in Figure 1, which are the best lectins to improve discrimination hY vs. hO or hO vs. CRC, respectively.
The IDI values show similar conclusions (Figure 4) as made for the NRI values (Figure 3). In case of the hY vs. hO, the highest values of IDI were observed for lectins SNA and WGA (Figure 4 left), while for the discrimination of hO vs. CRC, the following lectins showed high values: PHAE, WGA, DBA, RCA I, and selectins (Figure 4 right).
We can conclude that it is easier to identify lectins with strong discrimination power hY vs. hO or hO vs. CRC using the IDI matrix (Figure 4) compared to the NRI matrix (Figure 3).
Strong Discrimination Power for the Colorectal Cancer Matrix
In the CRC matrix, only the double lectin combinations with PHAE (AUC = 0.939, AUCage/AUCCRC = 0.657), RCA I (AUC = 0.811; AUCage/AUCCRC = 0.74), or HL sel (AUC = 0.889, AUCage/AUCCRC = 0.675) are evaluated.
All combinations of PHAE with any other lectin shown in Table 3 provide AUC higher than for a single PHAE (AUC of 0.939) with the discrimination AUCage/AUCCRC ratio below 0.8, as was confirmed for most combinations of PHAE with other lectins by an additional statistical evaluation (Figure 3 right and Figure 4 right). The lectin combination PHAE + HL sel provided the highest discrimination accuracy based on the AUC value (0.989).
Table 3 Clinical performance characteristics of double lectins for the CRC matrix with the best combinations showed in red.
Three combinations of RCA I with other lectins provided significantly higher AUC than a single RCA I lectin with a higher discrimination ratio AUCage/AUCCRC (Table 3). A combination of RCA I + HL sel provided the highest discrimination accuracy for the CRC matrix based on AUC value (0.961).
With regard to double lectin combinations with HL sel, a high discrimination power is a result of the combination with RCA I and PHAE lectins, which is the anticipated outcome, since RCA I and PHAE lectins have a strong discrimination power for the CRC matrix as single lectins. The best discrimination power based on the AUC value (0.989) for the CRC matrix was obtained using the combination of HL sel + PHAE (Table 3). Furthermore, statistical evaluation showed CIF well below a threshold value of 10 for every lectin combination (Figure 5 right) with no multicollinearity problem.
From Table 3, it can be concluded that the best double lectin biomarkers are a combination of PHAE + HL sel with AUC of 0.989 (Table 3) with high discrimination power for the CRC matrix (AUCage/AUCCRC = 0.741).
Several studies describe the use of glycan analysis for CRC diagnostics. All the papers discussed below are instrument-based approaches applied to glycan analysis.
Detection of IgG Fc N-glycopeptides by nanomaterial enrichment with subsequent MALDI-TOF mass spectrometry (MS) analysis of plasma samples from 46 CRC patients and 67 healthy individuals was evaluated in the form of an ROC curve. Machine learning algorithm using analysis of 11 N-glycopeptides was able to distinguish CRC patients from healthy controls with the average AUC value of the ROC of 0.893 (57).
An analysis of glycopeptides isolated from the serum of 80 CRC patients and 50 healthy individuals using an instrument-based approach revealed that especially leucine-rich α-2-glycoprotein with fucosylated triantennary N-glycan was a prospective diagnostic CRC biomarker (with high AUC of 0.86, sensitivity of 0.80, and specificity of 0.74) especially in combination with CEA (12).
Analysis of the ratio of two glycans with an AUC value above 0.95 was able to discriminate between CRC patients (n = 20) and healthy individuals (n = 20) (58).
The changed glycosylation of IgG isolated from the serum of CRC patients (n = 36), people with benign disease (n = 23), and a healthy control (n = 19) was also studied (59). The results indicate that glycans can only moderately discriminate the benign disease from the early-stage CRC (AUC = 0.75), and the early-stage CRC vs. the late-stage CRC (AUC = 0.75) with significant discrimination observed only between the late-stage CRC vs. benign disease (AUC = 0.85) (59). Accordingly, such an approach is not suitable for early-stage CRC diagnostics.
Another study revealed the CRC to be associated with a decrease in IgG galactosylation and IgG sialylation and an increase in core-fucosylation of neutral glycans with a concurrent decrease in the core-fucosylation of sialylated glycans. Glycan analysis rendered it possible to discriminate CRC patients from the control cohort with the AUC value of 0.755 (60).
Significant differences in the total serum N-glycome between CRC patients and the control cohort were used for CRC diagnostics (61). The authors observed an increased branching and sialylation for CRC patients, while the control cohort showed predominantly biantennary glycans. Glycan analysis could discriminate CRC patients from the control cohort with the AUC value of 0.81, sensitivity of 0.72, and specificity of 0.79. The 5-year survival rate largely varied between CRC patients with an altered serum N-glycome (46%) and an N-glycome similar to controls (87%). Importantly, the total serum N-glycome showed a prognostic value beyond age and stage (61).
It is worth mentioning that instrument-based methods are regarded as low-throughput and high-cost approaches to glycan analysis (12). The present study is the first using a reverse-phase lectin microarray for evaluation of the potential of glycans for CRC diagnostics with important clinical performance parameters obtained. It is difficult to identify protein carriers carrying aberrant glycans detected by lectins in this study, but some clues could be found in our recent study (41). A clear advantage of the approach presented here is that the most promising lectins (PHAE, RCA I, and HL sel) having a diagnostic potential for CRC could be used in the selective enrichment of serum proteins. In that case, lectin columns (PHAE, RCA I, and HL sel) will be used for the fractionation of serum proteins with subsequent protein identification using MS techniques. Hence, it would be possible to identify serum glycoproteins carrying aberrant glycans recognized by PHAE, RCA I, or HL sel, which could be more reliable CRC biomarkers based on glycans than, for example, CEA.
Such a lectin combination has the potential for use in the selective fractionation of serum proteins using lectin columns or lectins attached to magnetic particles. In that way, we will pre-enrich glycoproteins with specific aberrant glycosylation profiles. Then, we will identify such glycoproteins by MS using peptide mapping. Finally, such identified glycoproteins can be tested as potential novel glycan-based biomarkers for CRC diagnostics using specific glycoprofiling, i.e., application of an antibody selective against such a protein for selective capture from serum samples with subsequent glycoprofiling using lectins in a sandwich configuration (antibody/glycoprotein/lectin) using ELISA format of analysis, which is fully compatible with clinical practice (62, 63).
Conclusions
This study used novel recombinant lectins and lectins of human origin in combination with reverse-phase lectin microarrays to investigate glycosylation changes in the whole serum associated with age and with CRC. To the best of our knowledge, the lectin-recognizing Neu5Gc (HPyL) was used for the first time on reverse-phase lectin microarrays. The study identified lectins (SNA I + WFL) that can be specifically used for the discrimination of age and thus for determining biological human age. Then, we identified the combination of lectins (PHAE + HL sel), which can detect only glycosylation changes associated with CRC and not with age.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material. Further inquiries can be directed to the corresponding author.
Ethics Statement
The studies involving human participants were reviewed and approved by the ethical committee of University Hospital in Bratislava. The patients/participants provided their written informed consent to participate in this study.
Author Contributions
TB, AB, and JT contributed with concept formulation and the study design. TB, AB, MH, JA, and AV conducted experimental work. MP and LL provided samples and associated data and prepared the first draft of the article. TB and EJ performed statistical data evaluation and data presentation. TB, PK, LB, and JT wrote the article. All authors contributed to the article and approved the submitted version.
Funding
The financial support received from VEGA 2/0130/20 and from the Slovak Research and Development Agency APVV 17-0300 is acknowledged. The authors would like to acknowledge the support received from the Ministry of Health of the Slovak Republic under the project registration number 2018/23-SAV-1. The authors wish to acknowledge the financial support received from the Innovative Training Network (No. 813120). This publication was jointly supported by Qatar University and Chemical Institute, Slovak Academy of Sciences-IRCC-2020-004. The finding achieved herein is solely the responsibility of the authors. The statements made herein are solely the responsibility of the authors. This publication was created with the support of the Operational Programme Integrated Infrastructure for the project: Centre for Biomedical Research-BIOMEDIRES-II stage, ITMS: 313011W428, co-financed by the European Regional Development Fund.
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 wish to thank Simona Adamusová for optimization of the reverse-phase lectin microarray conditions.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fonc.2021.735338/full#supplementary-material
References
1. Xiao Y, Zhong J, Zhong B, Huang J, Jiang L, Jiang Y, et al. Exosomes as Potential Sources of Biomarkers in Colorectal Cancer. Cancer Lett (2020) 476:13–22. doi: 10.1016/j.canlet.2020.01.033
2. Bray F, Ferlay J, Soerjomataram I, Siegel RL, Torre LA, Jemal A. Global Cancer Statistics 2018: GLOBOCAN Estimates of Incidence and Mortality Worldwide for 36 Cancers in 185 Countries. CA: A Cancer J Clin (2018) 68:394–424. doi: 10.3322/caac.21492
4. Ogunwobi OO, Mahmood F, Akingboye A. Biomarkers in Colorectal Cancer: Current Research and Future Prospects. Int J Mol Sci (2020) 21:5311. doi: 10.3390/ijms21155311
5. Arnold M, Sierra MS, Laversanne M, Soerjomataram I, Jemal A, Bray F. Global Patterns and Trends in Colorectal Cancer Incidence and Mortality. Gut (2017) 66:683–91. doi: 10.1136/gutjnl-2015-310912
6. Brenner H, Kloor M, Pox CP. Colorectal Cancer. Lancet (2014) 383:1490–502. doi: 10.1016/S0140-6736(13)61649-9
7. Mármol I, Sánchez-de-Diego C, Pradilla Dieste A, Cerrada E, Rodriguez Yoldi MJ. Colorectal Carcinoma: A General Overview and Future Perspectives in Colorectal Cancer. Int J Mol Sci (2017) 18:197. doi: 10.3390/ijms18010197
8. Chand M, Keller DS, Mirnezami R, Bullock M, Bhangu A, Moran B, et al. Novel Biomarkers for Patient Stratification in Colorectal Cancer: A Review of Definitions, Emerging Concepts, and Data. World J Gastrointestinal Oncol (2018) 10:145. doi: 10.4251/wjgo.v10.i7.145
9. Oh H-H, Joo Y-E. Novel Biomarkers for the Diagnosis and Prognosis of Colorectal Cancer. Intestinal Res (2020) 18:168–83. doi: 10.5217/ir.2019.00080
10. Zhang DM, Xie Q, Wang Q, Wang YP, Miao JS, Li L, et al. Mass Spectrometry Analysis Reveals Aberrant N-Glycans in Colorectal Cancer Tissues. Glycobiology (2019) 29:372–84. doi: 10.1093/glycob/cwz005
11. Tang Y, Cui YD, Zhang SF, Zhang LJ. “The Sensitivity and Specificity of Serum Glycan-Based Biomarkers for Cancer Detection”. In: Zhang L, editor. Glycans and Glycosaminoglycans as Clinical Biomarkers and Therapeutics, Pt a. Cambridge, MA, USA. (2019). p. 121–40.
12. Shinozaki E, Tanabe K, Akiyoshi T, Tsuchida T, Miyazaki Y, Kojima N, et al. Serum Leucine-Rich Alpha-2-Glycoprotein-1 With Fucosylated Triantennary N-Glycan: A Novel Colorectal Cancer Marker. BMC Cancer (2018) 18(1):406. doi: 10.1186/s12885-018-4252-6
13. Huang C-S, Chen C-Y, Huang L-K, Wang W-S, Yang S-H. Prognostic Value of Postoperative Serum Carcinoembryonic Antigen Levels in Colorectal Cancer Patients Who Smoke. PLoS One (2020) 15:e0233687. doi: 10.1371/journal.pone.0233687
14. Norcic G. Liquid Biopsy in Colorectal Cancer-Current Status and Potential Clinical Applications. Micromachines (2018) 9:300. doi: 10.3390/mi9060300
15. Bhardwaj M, Gies A, Weigl K, Tikk K, Benner A, Schrotz-King P, et al. Evaluation and Validation of Plasma Proteins Using Two Different Protein Detection Methods for Early Detection of Colorectal Cancer. Cancers (2019) 11(10):1426. doi: 10.3390/cancers11101426
16. Jacobson RA, Munding E, Hayden DM, Levy M, Kuzel TM, Pappas SG, et al. Evolving Clinical Utility of Liquid Biopsy in Gastrointestinal Cancers. Cancers (2019) 11:1164. doi: 10.3390/cancers11081164
17. Namee NM, O’Driscoll L. Extracellular Vesicles and Anti-Cancer Drug Resistance. Biochim Biophys Acta Rev Cancer (2018) 1870:123–36. doi: 10.1016/j.bbcan.2018.07.003
18. Mannavola F, Salerno T, Passarelli A, Tucci M, Internò V, Silvestris F. Revisiting the Role of Exosomes in Colorectal Cancer: Where Are We Now? Front Oncol (2019) 9. doi: 10.3389/fonc.2019.00521
19. Xiao Y, Li Y, Yuan Y, Liu B, Pan S, Liu Q, et al. The Potential of Exosomes Derived From Colorectal Cancer as a Biomarker. Clin Chim Acta (2019) 490:186–93. doi: 10.1016/j.cca.2018.09.007
20. Yokose T, Kabe Y, Matsuda A, Kitago M, Matsuda S, Hirai M, et al. O-Glycan-Altered Extracellular Vesicles: A Specific Serum Marker Elevated in Pancreatic Cancer. Cancers (2020) 12:2469. doi: 10.3390/cancers12092469
21. Wagatsuma T, Nagai-Okatani C, Matsuda A, Masugi Y, Imaoka M, Yamazaki K, et al. Discovery of Pancreatic Ductal Adenocarcinoma-Related Aberrant Glycosylations: A Multilateral Approach of Lectin Microarray-Based Tissue Glycomic Profiling With Public Transcriptomic Datasets. Front Oncol (2020) 10. doi: 10.3389/fonc.2020.00338
22. Hirabayashi J, Yamada M, Kuno A, Tateno H. Lectin Microarrays: Concept, Principle and Applications. Chem Soc Rev (2013) 42:4443–58. doi: 10.1039/c3cs35419a
23. Cheng LM, Cao L, Wu YR, Xie WJ, Li JQ, Guan F, et al. Bisecting N-Acetylglucosamine on EGFR Inhibits Malignant Phenotype of Breast Cancer via Down-Regulation of EGFR/Erk Signaling. Front Oncol (2020) 10:14. doi: 10.3389/fonc.2020.00929
24. Yu HJ, Li XJ, Chen MT, Zhang F, Liu XW, Yu JM, et al. Integrated Glycome Strategy for Characterization of Aberrant LacNAc Contained N-Glycans Associated With Gastric Carcinoma. Front Oncol (2019) 9:11. doi: 10.3389/fonc.2019.00636
25. Beauharnois ME, Lindquist KC, Marathe D, Vanderslice P, Xia J, Matta KL, et al. Affinity and Kinetics of Sialyl Lewis-X and Core-2 Based Oligosaccharides Binding to L-And P-Selectin. Biochemistry (2005) 44:9507–19. doi: 10.1021/bi0507130
26. R.C. Team. R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing (2017).
27. Kuhn M, Wing J, Weston S, Williams A, Keefer C, Engelhardt A, et al. Package ‘Caret’. The R Journal (2020).
28. Robin X, Turck N, Hainard A, Tiberti N, Lisacek F, Sanchez J-C, et al. pROC: An Open-Source Package for R and S+ to Analyze and Compare ROC Curves. BMC Bioinf (2011) 12:77. doi: 10.1186/1471-2105-12-77
29. Carpenter J, Bithell J. Bootstrap Confidence Intervals: When, Which, What? A Practical Guide for Medical Statisticians. Stat Med (2000) 19:1141–64. doi: 10.1002/(SICI)1097-0258(20000515)19:9<1141::AID-SIM479>3.0.CO;2-F
30. Pencina MJ, D’Agostino RB Sr, D’Agostino RB Jr, Vasan RS. Evaluating the Added Predictive Ability of a New Marker: From Area Under the ROC Curve to Reclassification and Beyond. Stat Med (2008) 27:157–72. doi: 10.1002/sim.2929
31. Pencina MJ, D’Agostino RB Sr, Steyerberg EW. Extensions of Net Reclassification Improvement Calculations to Measure Usefulness of New Biomarkers. Stat Med (2011) 30:11–21. doi: 10.1002/sim.4085
32. Pencina MJ, D’agostino RB, Pencina KM, Janssens ACJ, Greenland P. Interpreting Incremental Value of Markers Added to Risk Prediction Models. Am J Epidemiol (2012) 176:473–81. doi: 10.1093/aje/kws207
34. O’Brien RM. A Caution Regarding Rules of Thumb for Variance Inflation Factors. Qual Quantity (2007) 41:673–90. doi: 10.1007/s11135-006-9018-6
35. Miura Y, Endo T. Glycomics and Glycoproteomics Focused on Aging and Age-Related Diseases — Glycans as a Potential Biomarker for Physiological Alterations. Biochim Biophys Acta (BBA) (2016) 1860:1608–14. doi: 10.1007/s10719-020-09958-4
36. Miura Y, Hashii N, Ohta Y, Itakura Y, Tsumoto H, Suzuki J, et al. Characteristic Glycopeptides Associated With Extreme Human Longevity Identified Through Plasma Glycoproteomics. Biochim Biophys Acta (BBA) (2018) 1862:1462–71. doi: 10.1016/j.bbagen.2018.03.025
37. Gudelj I, Lauc G, Pezer M. Immunoglobulin G Glycosylation in Aging and Diseases. Cell Immunol (2018) 333:65–79. doi: 10.1371/journal.pone.0083077
38. Dotz V, Wuhrer M. N-Glycome Signatures in Human Plasma: 5. FEBS Lett (2019) 593:2966–76. doi: 10.1128/iai.64.12.5421-5424.1996
39. Doherty M, Theodoratou E, Walsh I, Adamczyk B, Stockmann H, Agakov F, et al. Plasma N-Glycans in Colorectal Cancer Risk. Sci Rep (2018) 8. doi: 10.1016/j.semcancer.2010.04.005
40. Penezic A, Krizakova M, Miljus G, Katrlik J, Nedic O. Diagnostic Potential of Transferrin Glycoforms-A Lectin-Based Protein Microarray Approach. Proteomics Clin Appl (2019) 13. doi: 10.1128/JVI.03455-13
41. Bertok T, Jane E, Chrenekova N, Hroncekova S, Bertokova A, Hires M, et al. Analysis of Serum Glycome by Lectin Microarrays for Prostate Cancer Patients-a Search for Aberrant Glycoforms. Glycoconjugate J (2020) 37:703–11. doi: 10.1016/j.bbagen.2016.01.013
42. Cummings RD, Darvill AG, Etzler ME, Hahn MG. “Glycan-Recognizing Probes as Tools”. In: Varki A, Cummings R, Esko J, editors. Essentials of Glycobiology, 3rd edition. Cold Spring Harbor (NY: Cold Spring Harbor Laboratory Press (2015-2017).
43. Houser J, Komarek J, Kostlanova N, Cioci G, Varrot A, Kerr SC, et al. A Soluble Fucose-Specific Lectin From Aspergillus Fumigatus Conidia - Structure, Specificity and Possible Role in Fungal Pathogenicity. PLoS One (2013) 8:e83077/1–e83077/15, 15 pp. doi: 10.1016/j.cellimm.2018.07.009
44. Ruhl S, Sandberg AL, Cole MF, Cisar JO. Recognition of Immunoglobulin A1 by Oral Actinomyces and Streptococcal Lectins. Infection Immun (1996) 64:5421–4. doi: 10.1002/1873-3468.13598
45. Läubli H, Borsig L. Selectins Promote Tumor Metastasis. Semin Cancer Biol (2010) 20(3):169–77. doi: 10.1038/s41598-018-26805-7
46. Khan ZM, Liu Y, Neu U, Gilbert M, Ehlers B, Feizi T, et al. Crystallographic and Glycan Microarray Analysis of Human Polyomavirus 9 VP1 Identifies N-Glycolyl Neuraminic Acid as a Receptor Candidate. J Virol (2014) 88(11):6100–11. doi: 10.1002/prca.201800185
47. Holst S, Wuhrer M, Rombouts Y. Chapter Six - Glycosylation Characteristics of Colorectal Cancer. In: Drake RR, Ball LE, editors. Advances in Cancer Research. Academic Press (2015). p. 203–56.
48. Gizaw ST, Gaunitz S, Novotny MV. Highly Sensitive O-Glycan Profiling for Human Serum Proteins Reveals Gender-Dependent Changes in Colorectal Cancer Patients. Analytical Chem. San Diego, CA, USA. (2019) 91:6180–9. doi: 10.1021/acs.analchem.9b00822
49. Wong MY, Chen K, Antonopoulos A, Kasper BT, Dewal MB, Taylor RJ, et al. XBP1s Activation can Globally Remodel N-Glycan Structure Distribution Patterns. Proc Natl Acad Sci (2018) 115:E10089–98. doi: 10.1073/pnas.1805425115
50. Zhao Z, Feng Q, Yin Z, Shuang J, Bai B, Yu P, et al. Red and Processed Meat Consumption and Colorectal Cancer Risk: A Systematic Review and Meta-Analysis. Oncotarget (2017) 8(18):30563–75. doi: 10.18632/oncotarget.20667
51. Hofseth LJ, Hebert JR, Chanda A, Chen H, Love BL, Pena MM, et al. Early-Onset Colorectal Cancer: Initial Clues and Current Views. Nat Rev Gastroenterol Hepatol (2020) 17:352–64. doi: 10.1038/s41575-019-0253-4
52. Ahn S, Lee H, Lim H, Kim Y. Comparison of Serum Levels of Neu5Gc Between Normal and Colorectal Cancer Patients With HPLC-MS/MS. Clin Chim Acta (2019) 493:S70. doi: 10.1016/j.cca.2019.03.155
53. Vanhooren V, Dewaele S, Libert C, Engelborghs S, De Deyn PP, Toussaint O, et al. Serum N-Glycan Profile Shift During Human Ageing. Exp Gerontol (2010) 45:738–43. doi: 10.1016/j.exger.2010.08.009
54. Pučić M, Knežević A, Vidič J, Adamczyk B, Novokmet M, Polašek O, et al. High Throughput Isolation and Glycosylation Analysis of IgG–Variability and Heritability of the IgG Glycome in Three Isolated Human Populations. Mol Cell Proteomics (2011) 10(10):M111.010090. doi: 10.1074/mcp.M111.010090
55. Krištić J, Vučković F, Menni C, Klarić L, Keser T, Beceheli I, et al. Glycans Are a Novel Biomarker of Chronological and Biological Ages. J Gerontol: Ser A (2013) 69:779–89. doi: 10.1093/gerona/glt190
56. Franceschi C, Garagnani P, Parini P, Giuliani C, Santoro A. Inflammaging: A New Immune–Metabolic Viewpoint for Age-Related Diseases. Nat Rev Endocrinol (2018) 14:576–90. doi: 10.1038/s41574-018-0059-4
57. Zou Y, Hu JX, Jie JZ, Lai JY, Li MN, Liu ZM, et al. Comprehensive Analysis of Human IgG Fc N-Glycopeptides and Construction of a Screening Model for Colorectal Cancer. J Proteomics (2020) 213:103616. doi: 10.1016/j.jprot.2019.103616
58. Peng Y, Lv J, Yang LJ, Wang DQ, Zhang Y, Lu HJ. A Streamlined Strategy for Rapid and Selective Analysis of Serum N-Glycome. Analytica Chim Acta (2019) 1050:80–7. doi: 10.1016/j.aca.2018.11.002
59. Liu S, Cheng LM, Fu Y, Liu BF, Liu X. Characterization of IgG N-Glycome Profile in Colorectal Cancer Progression by MALDI-TOF-Ms. J Proteomics (2018) 181:225–37. doi: 10.1016/j.jprot.2018.04.026
60. Vučković F, Theodoratou E, Thaçi K, Timofeeva M, Vojta A, Štambuk J, et al. IgG Glycome in Colorectal Cancer. Clin Cancer Res (2016) 22:3078–86. doi: 10.1158/1078-0432.CCR-15-1867
61. de Vroome SW, Holst S, Girondo MR, van der Burgt YEM, Mesker WE, Tollenaar RAEM, et al. Serum N -Glycome Alterations in Colorectal Cancer Associate With Survival. Oncotarget (2018) 9(55):30610–23. doi: 10.18632/oncotarget.25753
62. Bertokova A, Bertok T, Jane E, Hires M, Ďubjaková P, Novotná O, et al. Detection of N, N-Diacetyllactosamine (LacdiNAc) Containing Free Prostate-Specific Antigen for Early Stage Prostate Cancer Diagnostics and for Identification of Castration-Resistant Prostate Cancer Patients. Bioorg Med Chem (2021) 39:116156. doi: 10.1016/j.bmc.2021.116156
Keywords: colorectal cancer, lectin, microarray, glycosylation, biomarker
Citation: Bertok T, Bertokova A, Jane E, Hires M, Aguedo J, Potocarova M, Lukac L, Vikartovska A, Kasak P, Borsig L and Tkac J (2021) Identification of Whole-Serum Glycobiomarkers for Colorectal Carcinoma Using Reverse-Phase Lectin Microarray. Front. Oncol. 11:735338. doi: 10.3389/fonc.2021.735338
Received: 02 July 2021; Accepted: 02 November 2021;
Published: 09 December 2021.
Edited by:
Kwong Tsang, Precision Biologics, Inc., United StatesReviewed by:
Pierpaolo Correale, Azienda ospedaliera ‘Bianchi-Melacrino-Morelli’, ItalyShu Zhang, Fudan University, China
Copyright © 2021 Bertok, Bertokova, Jane, Hires, Aguedo, Potocarova, Lukac, Vikartovska, Kasak, Borsig and Tkac. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jan Tkac, amFuLnRrYWNAc2F2YmEuc2s=