- Research article
- Open Access
- Open Peer Review
Exploring serum and immunoglobulin G N-glycome as diagnostic biomarkers for early detection of breast cancer in Ethiopian women
BMC Cancer volume 19, Article number: 588 (2019)
Alterations in protein glycosylation patterns have potentially been targeted for biomarker discovery in a wide range of diseases including cancer. Although there have been improvements in patient diagnosis and survival for breast cancer (BC), there is no clinically validated serum biomarker for its early diagnosis. Here, we profiled whole serum and purified Immunoglobulin G (IgG) fraction N-glycome towards identification of non-invasive glycan markers of BC.
We employed a comprehensive glycomics approach by integrating glycoblotting-based glycan purification with MALDI-TOF/MS based quantitative analysis. Sera of BC patients belonging to stages I-IV and normal controls (NC) were collected from Ethiopian women during 2015–2016. IgG was purified by affinity chromatography using protein G spin plate and further subjected to glycoblotting for glycan release. Mass spectral data were further processed and evaluated rigorously, using various bioinformatics and statistical tools.
Out of 35 N-glycans that were significantly up-regulated in the sera of all BC patients compared to the NC, 17 complex type N-glycans showed profound expression abundance and diagnostic potential (AUC = 0.8–1) for the early stage (I and II) BC patients. Most of these glycans were core-fucosylated, multiply branched and sialylated structures, whose abundance has been strongly associated with greater invasive and metastatic potential of cancer. N-glycans quantified form IgG confirmed their abundance in BC patients, of which two core-fucosylated and agalactosylated glycans (m/z 1591, 1794) could specifically distinguish (AUC = 0.944 and 0.921, p ≤ 0.001) stage II patients from NC. Abundance of such structural features in IgG is associated with a decrease in its immunosuppressive potential towards tumor cells, which in part may correlate with the aggressive nature of BC commonly noticed in black population.
Our comprehensive study has addressed for the first time both whole serum and IgG N-glycosylation signatures of native black women suffering from BC and revealed novel glyco-biomarkers with marked overexpression and distinguishing ability at early stage patients. Further studies on direct identification of the intact glycoproteins using a glycoprteomics approach will provide a deeper understanding of specific biomarkers towards their clinical utility.
Based on a recent report, breast cancer (BC) has become the most common malignancy and the leading cause of cancer death among women globally . In clinical settings, the disease progression is classified into four major stages (stages I-IV) based on pathological assessments named as TNM system (‘T’ for the size of the tumor; ‘N’ for the involvement of lymph nodes; ‘M’ for metastasis to other parts of the body). Distant metastases account for over 90% of breast cancer deaths with bone, lung, liver, and brain ranked as the vital target organs for breast cancer metastasis [2, 3]. If diagnosed at an early stage, BC has a high rate of survival (5-year survival rate after diagnosis is 99% when the tumor is still localized, but 27% for distant-stage disease) . The current method for early-stage diagnosis of BC is mammography. Nevertheless, in addition to its inaccurate diagnosis, many patients do not receive mammograms regularly due to lack of access to the facility, high cost, procedural discomfort, and perceived risks of radiation exposure . Consequently, there has been an increasing interest in the development and validation of serum-based non-invasive biomarkers whose alteration in the serum precedes the appearance of a malignancy, of which there is currently none .
Glycosylation, a complex post-translational modification involving covalent attachment of sugars to proteins or lipids, has recently received attention as a key component of cancer progression. Hence, unique expression or quantitative alterations in tumor-associated glycans can provide novel diagnostic and therapeutic targets [6, 7]. The fact that most human serum proteins are heavily glycosylated makes them a potential reservoir of glycans released from tissues and cells, reflecting their physiological and pathological states [8, 9]. In cancer progression, serum N-glycan alterations have been associated with proliferation, invasion, metastasis, aggressiveness, angiogenesis, and immune regulation of tumor cells .
With recent remarkable improvements in analytical tools, mass spectrometry is a widely used method to analyze glycans, proteins, or glycoproteins for clinical biomarker discovery . In this regard, previous studies have suggested N-glycan biomarkers for a wide range of cancer types [12,13,14,15,16,17,18,19]. However, the glycomic studies performed to date on breast cancer were focusing either on animal models, cultured cell lines, a small number of human patients, or didn’t employ the current state of methods for glycan purification and quantitative analysis from multiple types of biological samples [19,20,21,22,23,24]. Apart from total serum N-glycome, recent studies have also illustrated the association of altered glycosylation patterns of immunoglobulin G (IgG) with the progression of several tumors [25,26,27]. Using a recently developed novel glycoblotting method , we could screen serum N-glycan biomarkers for different pathological states including cancer, neurodegenerative diseases, and other eye, kidney, and arthritis related disorders [16,17,18, 29,30,31,32,33,34].
Driven by our advanced glycoblotting technology and the unmet needs for better molecular biomarkers and drug targets, this study was aimed to explore the potentials of N-glycans released from patient serum glycoproteins and from purified IgG fraction towards discovery of clinical biomarker for breast cancer. For this purpose, our study population were from Ethiopia where breast cancer incidence has been increasing drastically in younger women, while inadequate and ineffective control measures are increasing the mortality rate.
Study population and sample collection
Serum samples were collected from 115 BC female patients clinically belonged to stages I-IV and as normal controls (NC), serum samples from 33 gender, ethnic, and age matched healthy volunteers were obtained. All study subjects were Ethiopians and serum collection was performed during 2015–2016 in the Oncology Center of Black Lion Specialized Teaching Hospital (the largest and the only hospital with cancer treatment center in the country) of Addis Ababa University, Ethiopia. Informed consent had been signed and obtained from all study participants. The study was performed in accordance with the ethical standards of the Declaration of Helsinki on approval by the ethical review boards of Addis Ababa University, School of Medicine, Ethiopia and Hokkaido University, Faculty of Advanced Life Sciences, Japan. After freezed in plastic vials at − 80 °C and carefully packed with dried CO2 in foam box, serum samples were immediately transported by airplane and reached Japan within 24 h.
Reagents and materials
Ammonium bicarbonate 99% (ABC), 1-propanesulfonic acid, 2-hydroxyl-3-myristamido (PHM), 3-Methyl-1-p-tolyltriazene (MTT), and disialyloctasaccharide were purchased from Tokyo Chemical Industry Co., Ltd. (Tokyo, Japan). Peptide N-glycosidase F (PNGase F) was purchased from New England BiolabsR Inc. (Ipswich, Massachusetts, USA), whereas BlotGlyco H beads were purchased from Sumitomo Bakelite, Co. Ltd. (Tokyo, Japan). Dithiothreitol (DTT), iodoacetamide (IAA), O-benzylhydroxylamine hydrochloride (BOA), α-cyano-4-hydroxycin-namic acid (CHCA) and trypsin were from Sigma-Aldrich, Inc. (St. Louis, MO, USA). Other reagents and solvents were obtained from Wako Pure Chemical Industries, Ltd. (Tokyo, Japan), unless otherwise stated. Protein G spin plate was purchased from Thermo Scientific. SweetBlot™ (automated glycan processing and incubating machine) was from System Instruments Co., Ltd. (Hachioji, Japan). Multi Screen SlvinertR filter plates were purchased from Millipore Co., Inc. (Tokyo, Japan). All mass quantifications were done using MALDI-TOF/MS (Ultraflex III, Brukers Daltonics, Germany).
Immunoglobulin G purification from whole serum
IgG was purified from serum by affinity chromatography using 96-well protein G spin plate (Thermo Scientific) that is applied by fixing on another wash/collection plate as per the manufacturer’s protocol. First, 50 μL of serum was diluted 4 times with binding buffer (PBS, pH 7.2) and applied to the wells. The plate assembly was placed on a plate shaker and incubated for 30 min with moderate agitation, followed by centrifugation at 1000×g for 1 min with discarding the flow-through. The resin on the purification plate was washed 4 times by adding 300 μL of PBS to each well and then centrifuging at 1000×g for 1 min, discarding the flow-through each time. Before elution, 20 μL of neutralization buffer (1 M Tris-HCl, pH 9) was added to each well of new collection plate to maintain IgG stability. Subsequently, IgG was eluted by adding 200 μL of 0.1 M formic acid, pH 2.7 and centrifuging at 1000×g for 1 min. This step was repeated twice for complete elution. After confirming the purity of the IgG fraction by SDS-PAGE, it was dried using speed vac and then reconstituted in 50 μL of pure water to be applied for glycan release and purification by glycoblotting.
Enzymatic release of N-glycans and glycoblotting
All the procedures for N-glycan release, purification, labeling, and spotting were performed in the SweetBlot 7 automated system (all-in-one approach) as previously described [28, 35,36,37]. Accordingly, whole serum or IgG fraction purified from each sample was pretreated enzymatically for producing oligosaccharides carrying reducing terminal that will enable them to chemically ligate with hydrazide-functionalized BlotGlycoH beads during the subsequent glycoblotting process. Hence, 10 μL of whole serum or IgG fraction was auto-transferred into a 96 well polymerase chain reaction (PCR) plate and was dissolved with freshly prepared 0.33 M Ammonium bicarbonate (ABC) containing 120 mM 1, 4-dithiothreitol (DTT) and 0.4% of 1-propanesulfonic acid, 2-hydroxyl-3-myristamido (PHM) in 10 mM ABC. As an internal standard, 12 μL of 60 μM disialyloctasaccharide was also added and mixed in each well, enabling us to perform eventual absolute quantification of individual glycan level. Solubilized proteins were reduced by incubation at 60 °C for 30 min, followed by alkylation with 10 μL of 123 mM iodoacetamide in the dark at room temperature for 1 h. The mixture was then digested with 400 U of trypsin in 1 mM HCl at 37 °C for 2 h, and subsequently heated at 90 °C for 10 min to inactivate the enzyme. Soon after being cooled to room temperature, free N-glycans were released using 2 U of Peptide N-Glycosidase F (PNGase F) and incubating at 37 °C for 6 h. In the glycoblotting-based strategy (diagrammed in Fig. 1), 250 μL of BlotGlyco H bead (Sumitomo Bakelite Co., Ltd., 10 mg/mL suspension with water) was placed into a well of a Multi Screen Solvinert filter plate (Millipore), by vacuuming to remove the water. 20 μL of PNGase F digested mixture containing released N-glycans was mixed with the bead in each well, followed by the addition of 180 μL of 2% acetic acid (AcOH) in acetonitrile (ACN). To capture the N-glycans in sample mixtures specifically onto beads via reversible hydrazone bonds, the plate was incubated at 80 °C for 45 min. Then, the beads were washed twice with each 200 μL of 2 M guanidine-HCl in ABC, water, and 1% triethyl amine in methanol. Unreacted hydrazide functional groups on beads were capped by incubation with 10% acetic anhydride in MeOH for 30 min at room temperature. After removing the solution by vacuum, the beads were successively washed twice with each 200 μL of 10 mM HCl, MeOH, and dioxane. To stabilize sialic acids, on-bead methyl esterification of carboxyl groups in sialic acids was carried out by incubation with fresh 100 mM 3-methyl-1-p-tolyltriazene (MTT) in dioxane at 60 °C for 90 min. Each well was then washed twice using 200 μL of dioxane, water, methanol, and water again. The glycans captured on beads were subjected to the trans-iminization reaction with 20 μL of 50 mM O-benzyloxyamine hydrochloride (BOA) and with mild acid hydrolysis of hydrazone bond using 180 μL of 2% AcOH in ACN at 80 °C for 45 min. Finally, the released and BOA-tagged N-glycans were eluted with 75 μL of water under vaccum and applied for subsequent MALDI-TOF/MS analysis.
MALDI-TOF/MS analysis of N-glycans
BOA-labeled N-glycans were directly dissolved with an equivalent volume of self-prepared ionic liquid matrix solution (100 mM α-cyano-4-hydroxycinnamic acid diethylamine salt in MeOH: H2O: DMSO: 10 mM NaOH: = 50:39:10:1), after which 2.5 μL of the sample-matrix mixture was spotted onto MTP 384 target plate (polished steel TF, Bruker Daltonics). Each sample was spotted into 4 replicates to ensure reproducibility of the experiment. After being dried to crystallize the spots, N-glycans were analyzed using Ultraflex III MALDI-TOF/MS (Bruker Daltonics, Germany) by operation in positive-ion reflector mode, typically summing 1, 000 laser shots for each spot.
Data processing and evaluation
The acquired MALD-TOF/MS spectra were further processed and evaluated using FlexAnalysis v. Three Software (Bruker Daltonics, Germany) and Microsoft Excel. N-glycans were selected based on their quantitative reproducibility after evaluation using the calibration curve of serially diluted human serum standards. The intensity of monoisotopic peak of each N-glycan was normalized using 60 μM of an internal standard (Disialyloctasaccharide, produced by Tokyo Chemical Industry) and this normalized data was used for further statistical analysis and quantitative comparison. Structural compositions of N-glycans were speculated using GlycoMod Tool (http://br.expasy.org/tools/glycomod/).
Data for N-glycan abundance were analyzed using SPSS software. Multiple comparisons among NC and clinical stage groups were performed using a bonferroni one-way analysis of variance (ANOVA) while comparison between two (NC vs whole BC) groups was performed using Independent-Samples T-Test. Mean value differences were considered to be significant at 95% confidence interval (p ≤ 0.05). The diagnostic potential of significantly differed individual glycans or glyco-subclasses was further analyzed by receiver operator characteristic (ROC) test. Area under the curve (AUC) value generated from ROC test was used to assess the diagnostic accuracy of potential glycan biomarkers. Accordingly, AUC values of 0.9–1, 0.8–0.9, 0.7–0.8, and < 0.7 indicate a “highly accurate”, “accurate”, “moderately accurate”, and uninformative test, respectively. We also used Graph Pad prism 5 software to show the data distribution in scatter dot plot.
Descriptive information on the selected study participants is presented in Table 1. All the data were collected during the time of diagnosis and individuals who were pregnant, were receiving any treatment, or were with a medical history of disease which can affect the glycan profile, were excluded from the study. Age and gender were matched among the clinical groups. However, with the majority of the young patients, still the control group was found to be slightly younger (Table 1, A and B). Overall, the relatively younger age of the study population (uncommon in prior cancer studies) can be taken as an advantage in order to rule out the impact of aging on glycan profile.
Quantitative reproducibility test
Since each serum or IgG sample was spotted in four replicates, four normalized spectral data for each sample were averaged before statistical analysis. Before selection of the detected N-glycans, each peak was evaluated for its quantitative reproducibility using serially diluted standard human serum samples that had been experimented in the same plate beside the patient samples. The peak area of each glycan detected in the dilution series (0.5×, 0.75×, 1×, 1.25×, 1.5×, 1.75×, 2×, and 2.25×) was normalized using known concentration of internal standard, after which a standard calibration curve was plotted for each glycan as shown in Fig. 2. Quantitative reliability was then judged based on minimum outlier scores, good slope linearity, and significant level of the correlation coefficient across the standard curve. Accordingly, glycan peaks that meet these criteria were selected and considered for further statistical analysis in the main study samples.
Total serum and IgG N-glycan profiles
Using an integrated protocol for glycan release and purification, coupled with MALDI-TOF/MS analysis, 46 N-glycans (Table 2) could be detected in the serum of entire study samples. 37 (80.43%) of the detected N-glycans belong to complex type whereas high-mannose and hybrid types comprise 5 (10.87%) and 4 (8.9%), respectively. About one-third of these serum glycans were detected in the IgG fraction as well (indicated with “+” sign in Table 2). The large-scale mass spectra of serum N-glycans (Additional file 1: Figure S1) have shown quite reproducible patterns with nonnegligible differences in the peak intensity of many glycans (shown in Fig. 3) between patients and controls. Similarly, the collective mass spectra for IgG N-glycome is shown in Additional file 2: Figure S2 from which recognizable variations in terms of peak intensity and detection patterns between BC patients and controls can be seen in the individual mass spectrum (Fig. 4).
Evaluation for potential serum glycan biomarker identification
To identify individual glycans that showed differential expression patterns between patient and control groups, we performed Independent-Samples T-Test using the quantified data of each N-glycan peak. The majority of N-glycans exhibited up-regulated expression in which the levels of 35 N-glycans were significantly higher in the sera of whole BC patients compared to the controls (Additional file 3: Table S1). On the contrary, only two glycans showed reduction in the patient group in a slightly significant (m/z 1915, p = 0.048) and a non-significant (m/z = 2220) manner. Further multiple comparisons among the BC stages and NC group demonstrated that a substantial increase in glycan abundance was predominant at an early stage of the disease (Fig. 5). The serum glycans that had shown significant (p ≤ 0.05) expression alteration during ANOVA or T-test analysis were further considered for receiver operating characteristic (ROC) test. Their ability to distinguish cancer patients from controls was then evaluated based on their area under the curve (AUC) value generated from ROC analysis. As a result, 17 serum N-glycans with the highest diagnostic performance (detailed in Table 3) were selected as biomarkers for different BC stages and/or entire cancer patients. All these glycan structures belong to complex type N-glycans that include 4 bisecting, 8 core fucosylated, and 12 sialylated to a different degree with multiply branched antennas. Each of these identified glyco-biomarkers could demonstrate significant abundance with “highly accurate” (0.9–1) or “accurate” (0.8–0.9) diagnostic performance at an early stage (stage I, II, or both) or all patients, whereas no powerful biomarker was obtained exclusively at the late stages (BC-III or BC-IV) of the disease. As shown in Table 3 and Fig. 6, only 2 glycans that have bisected and fucosylated structures (m/w 2261 and 2264) for BC-III, and 3 monosialylated glycans (m/w 2261, 2439, 3007) for BC-IV met the biomarker criteria and adequately distinguished the patients from the controls.
Dot plots and ROC curves for some of the identified candidate biomarkers of stage I, stage II, or all BC patients compared with NC are diagramed in Fig. 5. Strong correlation with the disease stages was observed in which a substantial increase in the glycans abundance was predominant at an early stage of the disease. Comparing with the control group, individual expression scores within the cancer groups have shown dispersed distribution. As clearly stated in Fig. 6, discriminating power of the candidate biomarkers that fulfilled both criteria (p ≤ 0.05 and AUC ≥0.8) was either specific to stage I and stage II or shared by multiple cancer stages. Among the 17 selected candidate biomarkers, 13 glycoforms for BC stage I patients, and 12 glycoforms for BC stage II patients could strongly predict and distinguish them from the control group. As an important observation, the most branched and sialylated structures have shown marked abundance in the serum of all cancer stages with a greater increase at an early stage (Fig. 7). Pair wise, the first (m/z 3049 vs 3195), the second (m/z 3414 vs 3560), and the third (m/z 3719 vs 3865) pairs are biosynthetically consistent, only differed by core-fucose moiety.
Serum glycotyping pattern towards predicting early stage BC
Apart from individual serum N-glycan alteration, group of glyco-subclasses sharing certain structural features were also compared among different BC stages and controls. With no glyco-subclass altered significantly at stage III, we found statistically significant abundance in core-fucosylated, bisecting, bi-, tri-, and tetra-antennary, bi-, tri-, and tetra-sialylated N-glycans in the serum of stages I and II patients. From ROC analysis, greater diagnostic accuracy to predict early-stage BC was noticed in bi-antennary, bi-sialylated, tri-antennary, and tri-sialylated glycans (AUC = 0.895, .89, 0.863, 0.884, respectively) for stage I patients at p < 0.001. Moreover, tri-antennary, and tri-sialylated glycans demonstrated an AUC value of 0.88 and 0.902 to distinguish stage II patients from NC (Fig. 8). Larger abundance of these glycotypes associated with the early stages of the disease agrees with the aforementioned expression patterns of individual N-glycan biomarkers.
IgG-derived N-glycans as BC biomarkers
To determine whether the altered glycosylation patterns of total serum N-glycome were associated with the immune system or not, we purposely focused on in-depth IgG N-glycome profiling of the study subjects. The result illustrated statistically significant elevation of 5 glycans (m/z 1445, 1591, 1753, 1794, and 2423) in the whole patient group (Fig. 9a), of which 3 glycans (m/z 1591, 1753, 1794) had substantial increase in the stage II patients comparing to the controls (Fig. 9b). Overall, 12 of the 15 detected IgG glycans belong to complex biantennary type within which core-fucosylation was noticed as a structural feature in 9 of them, including the 3 abundantly expressed in patients. Moreover, in the same way as in serum, only two IgG glycans (m/z 1915 and 2220) showed an overall slight reduction in the patient group.
ROC test analysis indicated that, two IgG glycoforms (m/z 151 and 1794) could distinguish all BC patients and more strongly stage-II patients from the controls (AUC = 0.944 and 0.921, respectively; Fig. 10a, b). It should be noted that these glycans have shown quite similar patterns of abundance across the IgG and serum of all cancer stages and controls (Fig. 10b, c), reflecting that their carrier glycoprotein in the patients’ serum is mainly IgG. To determine IgG as a quantitatively reliable carrier of which circulatory N-glycans, ratio of N-glycan levels quantified from IgG fraction to that of whole serum were calculated and summarized in Additional file 4: Table S2. Thus, among the glycans whose larger proportion originated from IgG in BC patients include the glycans with m/z 1445 (53.25%), 1591 (64.76%), 1607 (41.72%), 1753 (59.41%), 1794 (46.64%), 1871 (53.24%), and 1915 (53.44%).
Apart from the cancer genome, deciphering alterations in protein glycosylation has been of critical importance not only to understand the molecular mechanisms of cancer progression, but also as a promising target for discovering novel diagnostic and therapeutic agents . Despite significant advances in breast cancer health care system, the current diagnostic methods give procedural discomfort and suffer from lack of reliability and specificity to predict the disease at its early stage. Hence, novel serological biomarkers and drug targets are still needed to further improve breast cancer diagnosis and treatment . The present study quantitatively profiled whole serum and IgG N-glycosylation alterations in Ethiopian breast cancer patients and matched controls. To the best of our knowledge, this is the first comprehensive report simultaneously addressing whole serum and purified IgG N-glycomics of breast cancer patients and more importantly employing study samples from native black African population whose glycosylation profile has not been studied elsewhere. We used a sensitive and an efficient glycoblotting method for glycan enrichment, after which we could identify and quantify a wide range of N-glycan structures by MALDI-TOF/MS analysis.
Structurally, majority of them belong to core-fucosylated, multiply branched and sialylated N-glycan types. With an overall quantitative up-regulation in all the breast cancer stages compared to NC, these serum glycans showed the highest abundance in the early stage (stage I and II) patients, making them promising diagnostic targets. Such a drastic and unidirectional quantitative shift over a wide-range of glycans was an unpredicted observation. ANOVA and ROC test analysis independently confirmed that these 17 N-glycans had high discriminating power to differentiate cancer stages from the healthy controls. Eight of these putative glyco-biomarkers could mutually distinguish (AUC = 0.8–1) both stage I and II patients from NC within which one glycan (m/z 2261) demonstrated high diagnostic performance in all stages of I-IV (Fig. 6).
Quantitative analysis of glyco-subclasses sharing the same structural residue (fucose, bisecting GlcNAc, antenna, and sialic acid) further intensified the predominant branching and sialylation features associated with early stage patients compared to controls. Among such glycotypes, it was observed that total bi-antennary and bi-sialylation for stage I, total tri-antennary and tri-sialylation for stage II strongly predicted and classified them from NC (Fig. 8). These results illustrate that aberrant serum glycan alterations in the breast cancer pathogenesis begin at an early stage of the disease and might be considered for non-invasive diagnostic biomarker identification. The current results are different from a previous study  that suggested 8 sialylated serum N-glycans (none of them had bisecting GlcNAc) as biomarkers associated with breast cancer. In the previous study, the predominant change (down-regulation of lower m/z and up-regulation of higher m/z glycans) of the suggested markers was observed in the late stage (stage IV) patients compared to controls. The candidate biomarkers identified in the present study, however, comprise asialo, bisecting, and hperbranched/hypersialylated glycans including only three of their reported glycan markers (m/z 2744, 3560, and 3719). Notably, our candidate biomarkers showed only an up-regulation pattern that was primarily noticed at the earlier stages (stage I or II) compared with the controls. As they used C18 column for glycan purification and permethylation (classical strategy)-based MALDI-TOF/MS analysis of relative intensities, it should be clear that the methodology and patient background of the two studies are quite different which makes their comparisons complicated. Thus, we hypothesize that the N-glycosylation changes associated with breast cancer progression in young black subjects seem to be unique and needs further attention.
It is well recognized that hypersialylation is a crucial feature in the progression of many cancers, whereby the negative charges in the terminal sialic acids of sugars interfere with epithelial cadherin-mediated cell–cell adhesion, enhancing the migratory and metastatic capacity of tumor cells [6, 38]. Interestingly, higher expressions of the glycosyltransferases responsible for core fucosylation (FUT8), branching (GnT-V), and sialylation (both α-2, 3 and α-2, 6 sialytransferases) have been associated with a greater potential for motility, invasion and metastasis in breast cancer [39,40,41,42]. One unexpected funding of our study was the up-regulation of a bisecting tetraantennary monosialylated glycan (m/z 3007) in the cancer groups as the N-acetylglucosaminyltransferase III (GnT-III) and its bisecting GlcNAc structures have been supposed to inhibit further branching in the biosynthesis pathway and suppress cancer metastasis . Kinetic studies on GnT-V, another key enzyme for the synthesis of branched N-glycans, indicated that the enzyme can rarely use the bisecting GlcNAc as a substrate with a very low Vmax value and produce branched bisecting structures . Taking this into account, finding such a bisected and branched structure whose abundance was strongly associated with BC stages I, II, and IV may provide significant implications for a drastic up-regulation of the preceding steps in the biosynthetic machinery that could unusually satisfy the low affinity of the branching enzyme. Apart from the overall glycan alterations indicated earlier, such specific biochemical flux may in part reflect the aggressiveness and invasiveness of the disease. Additionally, this glycan seems to be unique to Ethiopian population as it was detected in neither of Japanese hepatocellular carcinoma patients up on simultaneous experimentation nor in any of our previous results of various diseases.
It is unclear and hardly reported what kind of carrier proteins are involved in breast cancer associated alterations of the serum N-glycan level. In a more detailed analysis, we purified the IgG fraction from the whole serum of all study participants and then integrated (Fig. 1) to our high-throughput glycoblotting method for subsequent glycan release and purification. MALDI-TOF/MS based quantitative analysis confirmed that many of the biantennary glycans including some of the suggested serum biomarkers (m/z 1591 and 2118) for early stage BC were originated from serum IgG (Additional file 4: Table S2). Two core-fucosylated and agalactosylated IgG glycans (m/z 1591 and 1794) were significantly up-regulated in the breast cancer patients, more specifically distinguishing patients with stage II breast cancer from NC. Their expression patterns in IgG and serum of controls and stages I-IV patients were found to be consistent (Fig. 10a-c), which proves that IgG is a reliable carrier for a larger proportion of their concentration in serum. Being one of the major N-glycosylated serum proteins and an essential part of the humoral immune system, IgG structure and function is importantly modulated by the abundance and types of glycoforms attached to its fragment crystallizable (Fc) region . Increased core fucose in the IgG Fc domain has been shown to suppress its antibody-dependent cell-mediated cytotoxicity (ADCC) by decreasing its affinity towards cellular Fc receptors [46, 47]. On the other hand, increased IgG Fc galactose has recently been reported to enhance its complement-dependent cytotoxicity (CDC) activity via C1q binding [48, 49]. Based on our IgG glycomics result and these reports, the increased fucosylated and agalactosylated glycan structures observed in the cancer group appears to be a synergistic mechanism, allowing the tumor cells to escape from the immune system. In line with these results, increased IgG fucosylation (of non-sialylated glycans) and agalactosylation features have been previously reported in colorectal and gastric cancers [26, 27].
On the contrary, a recent study by Kawaguchi-Sakita N, et al.,  has reported increased galactosylated IgG N-glycans associated with Japanese breast cancer patients. Accordingly, 7 galactosylated biantennary N-glycans when merged into a group showed significant abundance and predicted the probability of BC with a moderate accuracy whereas, unlike in the present study, no single IgG glycan could strongly distinguish BC patients from NC. Such dissimilarities can be caused not only by differences in analytical procedures but also differences in the patient background such as race, life style, age (as the majority of their study subjects were already in the menopausal state; over 50s y/o), all of which have been reported to affect the glycosylation pathway . To ensure the reliability of our analytical procedure, quadruplicate of each sample was subjected to MALDI-TOF/MS analysis and the data for glycan levels were obtained from absolute quantification unlike many prior reports that had measured relative abundances. As such, the versatility of our glycoblotting method has been evidenced by its suitability for glycan purification from diverse samples (serum, cell lines, tissues, and CSF) [16,17,18, 29,30,31,32,33,34, 51] using only 10 μL of sample aliquots and overnight reaction. Apart from nonnegligible distinguishing features associated with BC, our IgG N-glycomics results clearly show that other unknown carrier proteins are greatly involved in the serum N-glycosylation alteration during breast cancer progression. Particularly, serum proteins that carry the hyperbranched and hypersialylated candidate glyco-biomarkers of BC are yet to be identified. Alpha-1-acid glycoprotein (AGP), an acute phase protein synthesized in the liver and secreted into the circulation, is believed to be one potential carrier for the bi-, tri-, or tetra- antennary complex type glycoforms. However, little is known about the association between breast cancer progression and AGP glycosylation pattern .
Many of the aberrant alterations observed in serum and IgG glycosylations of the early stage BC patients provide somehow different insight from previous reports which is most probably due to the relatively younger and black nature of the present study participants. In this regard, findings from various studies highlighted that young and black women are more likely to face an aggressive type of breast cancer, as it tends to have a poorer prognosis, a higher proliferative rate, a triple negative hormone receptor status, and a higher chance of breast cancer susceptibility genes (BRCA1 or BRCA2) mutations [52,53,54,55]. Based on the observations of the present study, being diagnostic biomarkers primarily at the early stage BC patients may raise a practical reliability doubt, especially in the context of African countries where awareness on cancer is low and patients are often diagnosed at advanced stages. Considering the advantage of early stage cancer markers to allow increased treatment options, this challenge can be alleviated by strengthening awareness creation and regular check-up strategies for breast cancer. In addition to the case-control groups, inter-and intra-subject variability of glycan levels were observed in the current study (much dispersed distribution has been shown within the cancer group). This suggests the potentials of protein glycosylation to be targeted for personalized medicine, as recently reported . The current study, however, had limitation to reveal the carrier proteins of the identified biomarkers apart from IgG. Also, the suggested biomarkers from relatively small sample size demands further verification study on a larger sample size of well-matched case-control subjects before using them in the clinical area.
In conclusion, this study comprehensively evaluated the total serum and IgG N-glycosylation signatures of native black cancer patients for the first time and identified novel N-glycan biomarkers showing strong diagnostic potential mainly at early stages of BC. Apart from the whole serum, our IgG focused N-glycome profiling provides evidence on BC-associated glycosylation alterations from the immune system perspective. Further study on direct identification of the intact glycoproteins using a glycoprteomics approach is expected to provide deeper understanding of specific biomarkers towards their clinical utility. Additionally, specific N-glycan profiling among the IgG subclasses of BC patients may lead to a focused diagnosis and therapy.
Availability of data and materials
The datasets used and analyzed during the current study are available from the corresponding author on reasonable request.
Antibody-dependent cell-mediated cytotoxicity
Analysis of variance
Area under the curve
Body mass index
Breast cancer susceptibility gene
Matrix-assisted laser desorption/ionization mass spectrometry
- PNGase F:
Receiver operating characteristic
Sodium dodecyl sulfate polyacrylamide gel electrophoresis
Tumor, Node, Metastasis
Fitzmaurice C, Allen C, Barber RM, et al. Global, regional, and National Cancer Incidence, mortality, years of life lost, years lived with disability, and disability-adjusted life-years for 32 Cancer groups, 1990 to 2015: a systematic analysis for the global burden of disease study. JAMA Oncol. 2017;3:524–48.
American Cancer Society. Breast Cancer Facts & Figures 2017-2018. Atlanta: American Cancer Society, Inc; 2017.
Xin J, Ping M. Targeting breast Cancer metastasis. Breast Cancer (Auckl). 2015;9:23–34.
Kevin DS, Jennifer B, Gerardine MS, Anthony MM. Alpha-1-Acid Glycoprotein (AGP) as a Potential Biomarker for Breast Cancer: IntechOpen Limited (UK); 2012. https://doi.org/10.5772/48177.
Misek DE, Kim EH. Protein biomarkers for the early detection of breast Cancer. Int J Proteomics. 2011;2011. https://doi.org/10.1155/2011/343582.
Pinho SS, Reis CA. Glycosylation in cancer: mechanisms and clinical implications. Nat Rev Cancer. 2015;15:540–55.
Stowell SR, Ju T, Cummings RD. Protein glycosylation in cancer. Annu Rev Pathol. 2015;10:473–510.
Clerc F, Reiding KR, Jansen BC, Kammeijer GS, Bondt A, Wuhrer M. Human plasma protein N-glycosylation. Glycoconj J. 2016;33:309–43.
Anderson NL, Anderson NG. The human plasma proteome: history, character, and diagnostic prospects. Mol Cell Proteomics. 2002;1:845–67.
Pinho SS, Carvalho S, Marcos-Pinto R, et al. Gastric cancer: adding glycosylation to the equation. Trends Mol Med. 2013;19:664–76.
Christopher AC, Stefani NT, Lori JS, Daniel WC. Advances in mass spectrometry-based clinical biomarker discovery. Clin Proteomics. 2016;13:1.
An HJ, Miyamoto S, Lancaster KS, Kirmiz C, Li B, Lam KS, Leiserowitz GS, Lebrilla CB. Profiling of glycans in serum for the discovery of potential biomarkers for ovarian cancer. J Proteome Res. 2006;5:1626–35.
Snyder CM, Alley WR Jr, Campos MI, Svoboda M, Goetz JA, Vasseur JA, Jacobson SC, Novotny MV. Complementary Glycomic analyses of sera derived from colorectal Cancer patients by MALDI-TOF-MS and microchip electrophoresis. Anal Chem. 2016;88:9597–605.
Sethi MK, Hancock WS, Fanayan S. Identifying N-glycan biomarkers in colorectal Cancer by mass spectrometry. Acc Chem Res. 2016;49:2099–106.
Qin R, Zhao J, Qin W, et al. Discovery of non-invasive glycan biomarkers for detection and surveillance of gastric Cancer. J Cancer. 2017;8:1908–16.
Nouso K, Amano M, Ito YM, et al. Clinical utility of high-throughput glycome analysis in patients with pancreatic cancer. J Gastroenterol. 2013;48:1171–9.
Ishibashi Y, Tobisawa Y, Hatakeyama S, Ohashi T, Tanaka M, Narita S, Koie T, Habuchi T, Nishimura S, Ohyama C, Yoneyama T. Serum tri- and tetra-antennary N-glycan is a potential predictive biomarker for castration-resistant prostate cancer. Prostate. 2014;74:1521–9.
Kamiyama T, Yokoo H, Furukawa J, Kurogochi M, Togashi T, Miura N, Nakanishi K, Kamachi H, Kakisaka T, Tsuruga Y, Fujiyoshi M, Taketomi A, Nishimura S, Todo S. Identification of novel serum biomarkers of hepatocellular carcinoma using glycomic analysis. Hepatology. 2013;57:2314–25.
Kyselova Z, Mechref Y, Kang P, Goetz JA, Dobrolecki LE, Sledge GW, Schnaper L, Hickey RJ, Malkas LH, Novotny MV. Breast Cancer diagnosis and prognosis through quantitative measurements of serum glycan profiles. Clin Chem. 2008;54:1166–75.
Haakensen VD, Steinfeld I, Saldova R, Shehni AA, Kifer I, Naume B, Rudd PM, Børresen-Dale A-L, Yakhini Z. Serum N-glycan analysis in breast cancer patients – relation to tumour biology and clinical outcome. Mol Oncol. 2016;10:59–72.
De Leoz ML, Young LJ, An HJ, Kronewitter SR, Kim J, Miyamoto S, Borowsky AD, Chew HK, Lebrilla CB. High-mannose glycans are elevated during breast cancer progression. Mol Cell Proteomics. 2011;10:M110.002717.
Goetz JA, Mechref Y, Kang P, Novotny MV. Glycomic profiling of invasive and non-invasive breast cancer cells. Glycoconj J. 2009;26:117–31.
Kirmiz C, Li B, An HJ, Clowers BH, Chew HK, Lam KS, Ferrige A, Alecio R, Borowsky AD, Sulaimon S, Lebrilla CB, Miyamoto SA. A serum glycomics approach to breast cancer biomarkers. Mol Cell Proteomics. 2007;6:43–55.
Lee LY, Thaysen-Andersen M, Baker MS, Packer NH, Hancock WS, Fanayan S. Comprehensive N-Glycome profiling of cultured human epithelial breast cells identifies unique Secretome N-glycosylation signatures enabling tumorigenic subtype classification. J Proteome Res. 2014;13:4783–95.
Kawaguchi-Sakita N, Kaneshiro-Nakagawa K, Kawashima M, Sugimoto M, Tokiwa M, Suzuki E, Kajihara S, Fujita Y, Iwamoto S, Tanaka K, Toi M. Serum immunoglobulin G fc region N-glycosylation profiling by matrix-assisted laser desorption/ionization mass spectrometry can distinguish breast cancer patients from cancer-free controls. Biochem Biophys Res Commun. 2016;469:1140–5.
Vučković F, Theodoratou E, Thaçi K, et al. IgG Glycome in colorectal Cancer. Clin Cancer Res. 2016;22:3078–86.
Kodar K, Stadlmann J, Klaamas K, et al. Immunoglobulin G fc N-glycan profiling in patients with gastric cancer by LC-ESI-MS: relation to tumor progression and survival. Glycoconj J. 2012;29:57–66.
Nishimura S-I. Toward automated glycan analysis. Adv Carbohydr Chem Biochem. 2011;65:219–71.
Hatakeyama S, Amano M, Tobisawa Y, et al. Serum N-glycan alteration associated with renal cell carcinoma detected by high-throughput glycan analysis. J Urol. 2014;191:805–13.
Gizaw ST, Koda T, Amano M, Kamimura K, Ohashi T, Hinou H, Nishimura S-I. A comprehensive glycome profiling of Huntington's disease transgenic mice. BBA-Gen Subjects. 2015;1850:1704–18.
Gizaw ST, Ohashi T, Tanaka M, Hinou H, Nishimura S-I. Glycoblotting method allows for rapid and efficient glycome profiling of human Alzheimer’s disease brain, serum and cerebrospinal fluid towards potential biomarker discovery. BBA-Gen Subjects. 2016;1860:1716–27.
Inafuku S, Noda K, Amano M, Ohashi T, Yoshizawa C, Saito W, Murata M, Kanda A, Nishimura S-I, Ishida S. Alteration of N-glycan profiles in Diabetic Retinopathy. Invest Ophthalmol Vis Sci. 2015;56:5316–22.
Hatakeyama S, Amano M, Tobisawa Y, Yoneyama T, Tsushima M, Hirose K, Yoneyama T, Hashimoto Y, Koie T, Saitoh H, Yamaya K, Funyu T, Nisimura S-I, Ohyama C. Serum N-glycan profiling predicts prognosis in patients undergoing hemodialysis. Sci World J. 2013;2013:268407.
Urita A, Matsuhashi T, Onodera T, Nakagawa H, Hato M, Amano M, Seito N, Minami A, Nishimura S, Iwasaki N. Alterations of High-mannose Type N-Glycosylation in Human and Mouse Ostesarthritis Cartilage. Arthritis Rheum. 2011;63:3428–38.
Kita Y, Miura Y, Furukawa J, Nakano M, Shinohara Y, Ohno M, Takimoto A, Nishimura S-I. Quantitative glycomics of human whole serum glycoprotein based on the standardized protocol for liberating N-glycans. Mol Cell Proteomics. 2007;6:1437–45.
Furukawa J, Shinohara Y, Kuramoto H, Miura Y, Shimaoka H, Kurogochi M, Nakano M, Nishimura S-I. Comprehensive approach to structural and functional glycomics based on chemoselective glycoblotting and sequential tag conversion. Anal Chem. 2008;80:1094–101.
Nishimura S-I, Niikura K, Kurogochi M, Matsushita T, Fumoto M, Hinou H, Kamitani M, Nakagawa H, Deguchi K, Miura M, Monde K, Kondo H. High-throughput protein glycomics: combined use of chemoselective glycoblotting and MALDI-TOF/TOF mass spectroscopy. Angew Chem Int Ed. 2004;44:91–6.
Pearce OM, Läubli H. Sialic acids in cancer biology and immunity. Glycobiology. 2016;26:111–28.
Tu C-F, Wu M-Y, Lin Y-C, Kannagi R, Yang R-B. FUT8 promotes breast cancer cell invasiveness by remodeling TGF-β receptor core fucosylation. Breast Cancer Res. 2017;19:111.
Handerson T, Camp R, Harigopal M, Rimm D, Pawelek J. Beta1,6-branched oligosaccharides are increased in lymph node metastases and predict poor outcome in breast carcinoma. Clin Cancer Res. 2005;11:2969–73.
Recchi MA, Hebbar M, Hornez L, Harduin-Lepers A, Peyrat JP, Delannoy P. Multiplex reverse transcription polymerase chain reaction assessment of sialyltransferase expression in human breast cancer. Cancer Res. 1998;58:4066–70.
Cui H, Lin Y, Yue L, Zhao X, Liu J. Differential expression of the alpha 2,3-sialic acid residues in breast cancer is associated with metastatic potential. Oncol Rep. 2011;25:1365–71.
Zhao Y, Nakagawa T, Itoh S, et al. N-Acetylglucosaminyltransferase III antagonizes the effect of N-Acetylglucosaminyltransferase V on α3β1 integrin-mediated cell migration. J Biol Chem. 2006;281:32122–30.
Sasai K, Ikeda Y, Fujii T, Tsuda T, Taniguchi N. UDP-GlcNAc concentration is an important factor in the biosynthesis of beta1,6-branched oligosaccharides: regulation based on the kinetic properties of N-acetylglucosaminyltransferase V. Glycobiology. 2002;12:119–27.
Arnold JN, Wormald MR, Sim RB, Rudd PM, Dwek RA. The impact of glycosylation on the biological function and structure of human immunoglobulins. Annu Rev Immunol. 2007;25:21–50.
Lauc G, Pezer M, Rudan I, Campbell H. Mechanisms of disease: the human N-glycome. BBA-Gen subjects. 2016;1860:1574–82.
Li T, DiLillo DJ, Bournazos S, Giddens JP, Ravetch JV, Wang L-X. Modulating IgG effector function by fc glycan engineering. Proc Natl Acad Sci. 2017;114:3485–90.
Peschke B, Keller CW, Weber P, Quast I, Lünemann JD. Fc Galactosylation of human immunoglobulin gamma isotypes improves C1q binding and enhances complement-dependent cytotoxicity. Front Immunol. 2017;8:646.
Li W, Zhu Z, Chen W, Feng Y, Dimitrov DS. Crystallizable fragment Glycoengineering for therapeutic antibodies development. Front Immunol. 2017;8:1554.
Pucic M, Knezevic A, Vidic J, Adamczyk B, Novokmet M, Polasek O, Gornic O, Supraha-Goreta S, Wormald MR, Redzic I, Campbell H, Wright A, Hastie ND, Wilson JF, Rudan I, Wuhrer M, Rudd PM, Josic D, Lauc G. 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:M111.010090.
Terashima M, Amano M, Onodera T, Shin-Ichiro Nishimura S-I, Iwasaki N. Quantitative glycomics monitoring of induced pluripotent- and embryonic stem cells during neuronal differentiation. Stem Cell Res. 2014;13:454–64.
DeSantis CE, Fedewa SA, Goding SA, Kramer JL, Smith RA, Jemal A. Breast Cancer statistics, 2015: convergence of incidence rates between black and white women. CA Cancer J Clin. 2016;66:31–42.
Swanson GM, Haslam SZ, Azzouz F. Breast Cancer among Young African-American women. Cancer. 2002;97:273–9.
Assi HA, Khoury KE, Dbouk H, Khalil LE, Mouhieddine TH, El Saghir NS. Epidemiology and prognosis of breast cancer in young women. J Thorac Dis. 2013;5(Suppl 1):S2–8.
Morris GY, Mitchell EP. Higher incidence of aggressive breast cancers in African-American women: a review. J Natl Med Assoc. 2008;100:698–702.
Almeida A, Kolarich D. The promise of protein glycosylation for personalised medicine. BBA-Gen Subject. 2016;1860:1583–95.
We would like to acknowledge Mr. Masakazu Tanaka and Mr. Tetsu Ohashi for their kind help in designing the experiments, data processing and evaluation. We are thankful for the study participants for the generous donation of their time and efforts. We are also grateful to Ms. Maki Morita and Ms. Tomoko Takahashi for their valuable technical assistance in preparation of the manuscript.
This work was partly supported by a grant for “Innovation COE project for future Medicine and Medicinal Research” from MEXT and JSPS KAKENHI Grant Number 25220206, JSPS Core-to-Core Program B. Asia-Africa Science Platforms, and by Japan Science and Technology Agency (JST) through a grant for “Development of Systems and Technology for Advanced Measurement and Analysis (SENTAN)” and “The Matching Program for innovations in Future Drug Discovery and Medical Care”. These funding bodies were involved with the design of the study, and collection, analysis, and interpretation of data.
Ethics approval and consent to participate
The study was approved by the ethical review boards of Addis Ababa University, School of Medicine, Ethiopia (reference number DRERC 04/14) and Hokkaido University, Faculty of Advanced Life Science, Japan (reference number 07–01) and was performed in accordance with the ethical principles of the declaration of Helsinki. Written informed consent was obtained from all study participants. All the information obtained was kept confidential. The right to participate, to withdraw from the study, or not to participate was respected.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Stacked-view MALDI-TOF Mass Spectra from large-scale serum N-glycomics of normal controls and BC patients. The raw mass spectra were further subjected to compositional analysis and structural elucidation. (TIF 8303 kb)
Figure S2. Large-scale MALDI-TOF Mass Spectra of IgG N-glycans of normal controls and BC patients. The raw mass spectra were further subjected to compositional analysis and structural elucidation. (TIF 7665 kb)
Table S1. List of serum N-glycans statistically differed between whole cancer patients and normal controls. All glycan peaks except m/z 1915 showed increased abundance (p ≤ 0.05) in the BC patients compared to NC. Peak numbers (peak #) are according to the list of detected glycans as shown in Table 2. (DOCX 15 kb)
: Table S2. Percentage of IgG N-glycan level relative to their concentration in serum. NC = normal control group, BC = whole breast cancer group. Indicated in bold are those N-glycans whose quantity in IgG showed considerable contribution to their expression level in serum. (DOCX 11612 kb)