In this study, we have identified PGD2S as a candidate biomarker for pediatric medulloblastoma after studying the extent of variability that exists in the CSF proteome. Protein abundance ratios and associated CV values were calculated to estimate the degree of CSF protein heterogeneity. In addition, the abundance ratios of proteins varied largely in both control and tumor groups in agreement with the high dynamic range reported in other studies [
22,
23]. We identified altered levels of apo
E, apo
J and PGD2S in the CSF of medulloblastoma patients. Alterations in the levels of apo
E have been reported in glioblastoma and various other cancers including pancreatic and lung [
24–
26]. On the other hand, up-regulation of apo
J, also known as clusterin, has been suggested to play a role in cancer progression by conferring a protective role against apoptosis in cancer cells [
27,
28]. Similar to our finding, Pucci
et al reported over expression of apo
J in the cytoplasm of primary colon cancer cells that is released into the extracellular space [
29].
The most intriguing finding in our study is the alteration of prostaglandin D2 synthase isoforms (PGD2S), a glycoprotein that is abundant in CSF. PGD2S, also known as β- trace, is one of the most abundant glycoprotein in human CSF [
30]; it is synthesized and secreted by both glial cells and the choroid plexus [
31]. An ideal biomarker is one that is at least expressed in higher amounts (if not exclusively ) in the CSF of children with tumors compared to controls, in which case, it is logical to speculate that the tumor cells are the source of that protein. Conversely, since PGD2S is reduced in tumor CSF samples, we speculate that the reduction is a host response to the presence of the tumor. Consistent with this hypothesis is the fact that PGD2S abundance did not vary depending upon disease burden (e.g. between tumors with metastases and completely resected tumors) even when measured in samples obtained 10–14 days after tumor resection. One potential advantage of a biomarker derived from a host response is that it can be amplified, analogous to the way that serum antibody levels created by an immune response can indicate previous exposure to an infectious agent. Although PGD2S has been investigated in various diseases, only two studies quantitated PGD2S levels in CSF [
32–
36]. Saso
et al reported reduced levels of PGD2S in the CSF of brain tumor patients; however, the study did not specify the types of tumors that were investigated [
35]. Our study differs in its specific focus on medulloblastoma and the pediatric population. The other report by Huang
et al showed a reduced level of PGD2S in the CSF of patients with acute inflammatory demyelinating polyneuropathy [
36]. There is still very limited knowledge on the mechanisms involved in the alteration of PGD2S levels in CSF. The known established function of PGD2S is the synthesis of prostaglandin PGD2 from PGH2. Therefore, we quantitated levels of PGD2 in our CSF samples using ELISA but found no correlation to PGD2S quantities in CSF samples from tumor and control subjects (data not shown) most likely due to its instability in biological fluids and rapid degradation to the 15-deoxy prostaglandin J2 metabolite [
37].
We performed statistical analysis to investigate confounding variables that could affect the proteome. Age was not found to be a significant covariate in the comparative analysis of PGD2S spot intensities between control and tumor samples. Likewise, there was no significant variation based upon the presence of metastases or histological subtypes (see
supporting information). None of the PGD2S isoform levels were different between CSF samples associated with gross residual disease and those associated with completely resected disease. This observation could suggest that either these spots are sensitive enough to indicate the presence of minimal residual disease after surgical removal of the tumor or that the turnover of PGD2S in the CSF requires more than 2 weeks. The remaining, though unlikely, possibility is that PGD2S levels never return to normal levels. Further work with serial samples will be required to determine the true explanation.
Because neither age, histological subtype, effect of surgery nor presence of metastases correlated with the intensity of the PGD2S spots, the observed alteration in PGD2S was not attributed to any factor other than disease status. Furthermore, the consideration of blood contamination was rendered moot by the fact that PGD2S is a protein. Proteins that exhibit great inter-individual variability make poor candidate biomarkers for disease. As depicted in the
supporting information, although the CV values of individual spots attributable to inter-individual variation showed no dependency on spot abundance ratios, average CV was found to be 0.7. In other words, a typical CSF protein varies between different individuals by about 70%. Knowledge about inter-individual proteome variations in human samples (both tissue and biological fluid) is currently limited with only two published reports on this topic. The first study reported an average CV of 0.18 within the platelet proteome [
38]. In the second, Hu
et al conducted both intra and inter-individual variability studies on the CSF proteome and though they reported that the inter-individual variability is much more than the intra-individual variation, they did not report the average CV [
10]. Given the high variability that exists in the CSF proteome, an understanding of the contribution of technical variation is also critical for differential proteomic investigations. In our hands, the contribution of technical variability to the observed total variability was 9%, well within the anticipated values for 2-DE [
39].
In addition to the decrease in total PGD2S levels, there was a marked difference in the levels of various isoforms between patients and controls. In particular, there was at least a 2 fold reduction of the acidic isoforms of PGD2S in medulloblastoma samples (see ). This prompted us to identify the glycosylation patterns in the various isoforms of PGD2S. Since the linear MS spectra for all three spots exhibited pairs of peaks with masses of 4163 and 3870 as well as 3960 and 3667, which differ in mass by 292 Da (the mass of a sialic acid molecule), we hypothesized that all three spots contained sialylated glycopeptides but different number of sialic molecules. This was confirmed by neuraminidase treatment (see ). We were able to detect only two of the three predicted glycopeptides from PGD2S, one spanning amino acid sequence 67 to 86 and the other starting at 87 and ending at 92. We speculate that these PGD2S isoforms are a result of differential glycosylation patterns of the three predicted sites in the protein. Further analysis of purified PGD2S using accurate mass measurements is required to fully characterize its glycosylation pattern in CSF but such an analysis is beyond the scope of this clinical study. In conclusion, the observations described herein indicate that PGD2S is a CSF biomarker for medulloblastoma, which could be useful in detecting the efficacy of treatment and the recurrence of disease.