Next Article in Journal
Available In Vitro Models for Human Satellite Cells from Skeletal Muscle
Next Article in Special Issue
A Novel Role of BIRC3 in Stemness Reprogramming of Glioblastoma
Previous Article in Journal
Restraint upon Embryonic Metatarsal Ex Vivo Growth by Hydrogel Reveals Interaction between Quasi-Static Load and the mTOR Pathway
Previous Article in Special Issue
Microglia-Derived Olfactomedin-like 3 Promotes Pro-Tumorigenic Microglial Function and Malignant Features of Glioma Cells
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Profiling Glioblastoma Cases with an Expression of DCX, OLIG2 and NES

1
Department of Biochemistry and Molecular Biology, Medical University of Lublin, 20-093 Lublin, Poland
2
Laboratory for Cytogenetics and Genome Research, Department of Human Genetics, KU Leuven, B-3000 Leuven, Belgium
3
Department of Neurosurgery and Pediatric Neurosurgery, Medical University of Lublin, 20-090 Lublin, Poland
4
Institute of Biological Bases of Animal Production, University of Life Sciences, 20-950 Lublin, Poland
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2021, 22(24), 13217; https://doi.org/10.3390/ijms222413217
Submission received: 17 October 2021 / Revised: 29 November 2021 / Accepted: 1 December 2021 / Published: 8 December 2021
(This article belongs to the Special Issue Advances in Molecular Genetics of Brain Tumors)

Abstract

:
Glioblastoma (GBM) remains the leading cause of cancer-related deaths with the lowest five-year survival rates among all of the human cancers. Multiple factors contribute to its poor outcome, including intratumor heterogeneity, along with migratory and invasive capacities of tumour cells. Over the last several years Doublecortin (DCX) has been one of the debatable factors influencing GBM cells’ migration. To resolve DCX’s ambiguous role in GBM cells’ migration, we set to analyse the expression patterns of DCX along with Nestin (NES) and Oligodendrocyte lineage transcription factor 2 (OLIG2) in 17 cases of GBM, using immunohistochemistry, followed by an analysis of single-cell RNA-seq data. Our results showed that only a small subset of DCX positive (DCX+) cells was present in the tumour. Moreover, no particular pattern emerged when analysing DCX+ cells relative position to the tumour margin. By looking into single-cell RNA-seq data, the majority of DCX+ cells were classified as non-cancerous, with a small subset of cells that could be regarded as glioma stem cells. In conclusion, our findings support the notion that glioma cells express DCX; however, there is no clear evidence to prove that DCX participates in GBM cell migration.

1. Introduction

Even though malignant brain tumours account only for a small percentage of all adult cancers, they lead to an extensive amount of cancer-related deaths [1]. Moreover, the five-year survival rates are among the lowest for all human cancers [2], regardless of treatment modality [3]. This remarkable resistance results mostly from tumour heterogeneity and its high propensity for malignant progression. One of the vital pathophysiologic features contributing to this dismal prognosis is their strong migrational capacity [4] for significant dispersal beyond the macroscopic tumour borders [5]. Interestingly, a similar migratory ability is one of the principal features of neuronal progenitor cells (NPC) during CNS development [6,7]. Given that, it is worth noticing that data generated by the Cancer Genome Atlas Research Network proved that one of the main gene profiles of GBM, a proneural profile, involves DCX expression [8].
DCX itself, when it is mutated, is responsible for an X-linked form of lissencephaly, affecting the organisation of neocortical layering in the cerebral cortex [9]. Subsequent studies have shown that DCX directly binds to microtubules, thereby regulating their polymerisation and stabilisation [10]. This process is crucial for a multipolar mode of neuroblast migration [11], a transient stage in neuronal progenitor migration where migrating cells search for environmental signals that will determine their mode of migration [12]. DCX expression has been restricted to migrating neuroblasts in developing and adult animals [12]. However, Daou et al. [13] proved that DCX expression might be found in various neuroepithelial origin tumours.
Interestingly, DCX was highly expressed in both high-grade and low-grade invasive tumours. Moreover, invasive tumours have been shown to express higher levels of DCX when compared to circumscribed tumours; no expression in normal brain tissue surrounding the tumour was found. A more recent follow-up study evaluated the sensitivity and specificity of DCX immunostaining to detect infiltrating glioma cells [14]. It confirmed that DCX is explicitly expressed in infiltrating gliomas but not in reactive astrocytes. Santra et al. presented a different hypothesis [15]. Their data indicate that DCX mRNA transcripts were not detected in primary glioma cells, while DCX expressing cells were revealed in tumour penumbra. DCX positive cells within glioma tumours, in their opinion, were either infiltrating neuroblasts or pre-existing neuronal cells.
According to Verhaak et al. [8], one of the proneural subtype’s signature genes is Oligodendrocyte lineage transcription factor 2 (OLIG2), a family member of basic helix–loop–helix transcription factors. It plays a crucial role in the early stages of brain development in oligodendrocyte precursor cells and neural progenitor cells by enhancing mitosis and limiting cell differentiation [16]. In gliomas, OLIG2 was expressed by glioma stem cells [17]. On top of that, cells expressing Rai (ShcC/N-Shc), a member of the family of Shc-like adaptor proteins, are involved in non-neoplastic cell migration, co-expressed DCX, and OLIG2 [18]. Accordingly, whenever DCX and OLIG2 is observed in migrating cells, it would imply that these cells have stem-cell capabilities.
On the other hand, Bott et al. [19] recently found that DCX function in complex with nestin (NES) is a marker of neural progenitors. Although initially detected in neuronal stem cells, its presence in multiple other tissues (including gliomas) has recently been described. NES is a type VI intermediate filament. It plays a role in several key aspects of primary cell functioning: self-renewal, proliferation, survival, differentiation, and migration. With regard to the later, Bott et al. established that cdk5/p35 selectively phosphorises DCX due to the DCX-NES complex’s presence. That could directly influence the growth cone during migration.
Based on the observations mentioned above, we contend that the role of DCX expression in glioma cells’ migration is still a matter of debate. Therefore, in order to further clarify this role, we decided to conduct an experiment focusing on the distribution of glioma cells, both inside and in the tumour’s margin. Simultaneously, in order to elucidate the possible role of DCX in the context of the proneuronal subtype pathogenesis with its potential interactions with microtubules and intermediate filaments, we have marked OLIG2 and NES.

2. Results

Our heuristic approach to quantify the number of cells with each of the antigens allowed us to obtain data that can be compared with the results from single-cell RNA-seq. The fundamental component of this method is the determination of cell nuclei in the examined ROI. Cellular nuclei (Figure 1A) act as seeds to assess adjacent fluorescent signals: OLIG2 was expected to overlap with nuclei (Figure 1B,F,H), whereas NES (Figure 1C,E,H) and DCX (Figure 1D,G,H) signal should tightly adhere to it.
As expected, cell counts were significantly different in various parts of the tumour. More cells were tumour-adjacent to the margin or in non-specific tumour sites than in the margin (p = 1 × 10−4 and p = 1.9 × 10−7, respectively)—Figure 2H. No significant cell count change was detected while comparing tumours adjacent to margin and non-specific tumour sites (p = 0.91).
The DCX+ cells were only a small subpopulation of cells in total, although they were present in all types of tissues—Figure 2A. The median frequency of DCX+ cells was 3.1% in margin tissue and 1.1% in adjacent tumour tissue, without significant difference (p = 0.82). The majority of ROIs did not have a context of margin tissue. Thus, it was impossible to assess the exact placement within the tumour. They were grouped and called other/non-specific tumour tissue. A batch comparison of DCX+ cells in other tumour tissue versus tumour-adjacent to margin did not show statistical importance (p = 0.24).
Contrary to the frequency of DCX+ cells, NES+ cells were much more abundant in studied cases (Figure 2B). The median frequency was 17.7% in margin tissue, 20.6% in adjacent tumour tissue, and 25.4% in the rest of ROIs tumours, yet without significant difference (the lowest p = 0.78 in the batch comparison, and p = 0.30 when accounting pairs of margins vs. adjacent tumour).
OLIG2+ cells were more frequent than DCX+ cells (Figure 2C). Non-specific tumour sites and tumours close to the margin had a similar median (8.2% and 5.1%, respectively). The median of the tumour margin, on the other hand, was 19.8%. Pairwise comparison of the margin and the tumour-adjacent to the margin shows a significant change in OLIG2+ cells frequency (p = 4.9 × 10−4).
No significant changes in frequency in cells expressing more than one of the analysed markers were detected. The majority of measured frequencies were next to 0% (Figure 2D–G).
Furthermore, we tested if subsequent tumour resections impacted cell count with specific staining (Figure 3A–G). We have focused on the first three subsequent resections, as they have a representative number of cases. There were no significant trends in all stainings (DCX+: p-value = 0.90, OLIG2+: 0.72, NES+: p-value = 0.89; df = 2).
Then, we tested if there are any correlations between the frequency of cells with different antigens and sites (Supplementary Figure S1). The frequency of cells expressing either NES or OLIG2 had a strong positive correlation between margins and tumour sites adjacent to margins. Although cells expressing DCX lacked this correlation, there was a strong negative correlation between the frequency of DCX+ cells in the margin and the frequency of NES+ cells both in tumour and margin sites.
Finally, we did not observe any change in overall survival regarding tested antigens (Supplementary Figure S2).
To better understand the function of DCX, we opted for the single-cell approach, which allowed annotating the cells as “cancer” and “normal” according to the number of GBM specific copy number variations (CNV) prediction (Supplementary Figure S3) based on transcriptomics of single cells (Figure 4A). To avoid false-positive identification of malignant cells, we considered only the canonical aberrations for GBM, which are a duplication of chromosome 7 or chromosome 10 loss. The cells defined as “cancer” and “normal” were 39.77% and 60.23%. The cells identified as ‘normal’ were clustered together, mainly in the two big clusters—one of them was specific for the cells derived from the foetal origin, whereas the second big cluster contained the non-cancer brain cells (Figure 4B–D). The ‘cancer’ cells, which were defined as GBM specific CNVs containing, had been clustered into the number of overlapping and sample-specific clusters, which shows their high heterogeneity (Figure 4A–C).
Based on the single-cell transcriptomics, we also annotated the cluster cells to define the analysed cells’ phenotype (Figure 5A–C). The complete list of cluster-specific genes are attached to the Supplementary Data chapter as a CSV file, whereas the distribution of the cells among the clusters by their “cancer” and “normal” status are shown in Supplementary Figure S4. These annotations were used to predict the role of DCX, OLIG2, and NES in tumour biology. We visualised the expressions of these genes that were shown among defined cell clusters in Figure 6. The cells defined as “normal” showed the expression of specific markers for neurons, astrocytes, dendritic cells, macrophages, plasmacytoid dendritic cells, circulating foetal cells, undifferentiated cells, and oligodendrocytes. These cells that had been labelled as “cancer” expressed the markers specific for circulating foetal cells, basophils, astrocytes, dendritic cells, neurons, proliferative cells, undifferentiated cells, and oligodendrocytes. The cells’ distribution by their “cancer” or “normal” status is shown in Supplementary Figure S4.
OLIG2 level in the cells annotated as “cancer” was found only in this cluster of cells which also expressed the markers specific for astrocytes, neurons, and oligodendrocytes; however, it also was found in the cells with undifferentiated and proliferative properties. NES was the most specific for the cells annotated as ‘cancer’, for which basophil markers are characteristic. DCX in the “cancer” identified cells was expressed in these clusters, identified according to the transcriptomics signatures specific for oligodendrocytes, undifferentiated, and proliferative cells (Figure 6A). In normal cells, OLIG2 was found mainly in the undifferentiated cells and in lower levels in clusters of proliferative cells and this group of cells with the gene signature that is specific for neurons, lymphocytes, and macrophages. NES expression was mostly seen in the clusters of cells that expressed the markers of neurons, lymphocytes, and macrophages and those that harboured the transcripts specific for proliferative and undifferentiated cells. DCX in normal cells was specific for neurons and astrocytes (Figure 6B).
We also evaluated the correlations between the expression of DCX, OLIG2, and NES genes at a single-cell level in the population of cells derived from the tumour tissue to show the potential relation of these genes. We found significant (p < 0.01) but rather weak positive correlations between transcript levels of these genes conducting Spearman’s rank test: DCX vs. OLIG2 (r = 0.2179733), NES vs. OLIG2 (r = 0.3019327), NES vs. DCX (r = 0.1535854). When we look closer at the populations of these cells that had been annotated as ‘normal’, simultaneous expression of OLIG2 and DCX was observed in the clusters of these cells that were expressing markers specific for highly proliferative or undifferentiated cells as well as in the cluster of cells expressing markers that are common for astrocytes, oligodendrocytes and neurons. On the other hand, in the population of cancerous cells, NES and DCX were simultaneously expressed by the cells that were positive mostly for markers of: proliferative, undifferentiated or astrocytic cells.

3. Discussion

DCX status in GBM has been a subject of debate since Rich [20] and Daou [13] reported the expression of DCX on mRNA and protein levels in 2005. While Rich et al. correlated the mRNA level of DCX with poor diagnosis, Daou et al. proved that DCX had more intense staining towards the margin of the tumour using immunostaining. Not all samples of GBM had the DCX expression, however. Interestingly, our results did not confirm most of those observations. Although DCX was detected in most samples, DCX+ cells were only a small subset of the GBM’s cell population, both in tumours and margins. Only in a few cases were there many DCX+ cells. Moreover, there was no shift in DCX+ cell frequency towards one of the sites. Although inconclusive, this may imply that DCX is not directly connected to GBM cell migration.
On the other hand, Santra et al. [15,21] reported DCX as a marker for a favourable patient outcome. They found that cells with DCX overexpression had lower invasion abilities, thus, the authors concluded that DCX positive cells in glioma either infiltrate neuroblasts or pre-existing neuronal cells. Importantly, our data partially supports this idea. Although most of the cells expressing DCX within the tumour were marked as non-cancerous, four subpopulations were marked as GBM origin. They followed expression patterns similar to oligodendrocytes, astrocytes, as well as proliferative and undifferentiated cells. This aligns with the idea that DCX is expressed by GSC.
In general, NES+ and OLIG2+ cells were more frequent than DCX+ ones, which is consistent with the bulk values reported in The Human Protein Atlas [22]. Interestingly, cells expressing NES or OLIG2 were present both in the margin and in tumours, while being highly correlated. The lack of a significant drop in NES+ cells is in contradiction with a study by Smith et al. [23]; this might be explained by a different definition of margin adopted in our study. Smith et al. defined margin as a region with 5-aminolevulinic acid (5ALA) fluorescence during surgery that is beyond the T1 enhancing region on magnetic resonance imaging (MRI) [23]. However, when considering regions by their relative position to tumour sites, both studies are consistent. We also found that OLIG2+ cells frequency differed among sites: more cells with OLIG2 expression were in margin than in tumour.
On the other hand, the frequency of NES+ cells on the margin and within the tumour were reversely correlated with DCX+ cells in the margin. A negative correlation seems not to support the existence of the DCX-NES complex in the tumour. It is worth keeping in mind that there is a gross difference between the frequency of DCX- and NES-positive cells, and sample size does not compensate for that. When we look closer at the populations of these cells that had been annotated as ‘normal’, simultaneous expression of OLIG2 and DCX was observed in the clusters of these cells that were expressing markers specific for highly proliferative or undifferentiated cells as well as in the cluster of cells expressing markers that are common for astrocytes, oligodendrocytes and neurons. On the other hand, in the population of cancerous cells, NES and DCX were simultaneously expressed by the cells that were positive mostly for markers of: proliferative, undifferentiated or astrocytic cells.
Another interesting finding was that subsequent resections did not significantly influence assessed markers’ frequency. In fact, all three of them seemed to be associated with some form of undifferentiated cells. It was reported previously that glioma stem cells have higher invasion capabilities [24], suggesting that stemness might be a phenotypic response to changing the tumour’s microenvironment [25]. Our data also support the notion that there is a set proportion of cells expressing DCX/NES/OLIG2, regardless of treatment, guided by tumour plasticity.
The GBM cells revealed high heterogeneity, which we saw as clustering into several different groups of cells derived by transcriptomics. Others also reported a similar relationship [26], which confirms the value of our analytical approach. However, annotating the cells to “cancer” or “normal” subpopulations using our quite simple approach may not be entirely accurate as we might omit other possible genetic aberrations than those used canonically for GMB. Nonetheless, our approach seems to be valuable, for it limits the false positive detection of malignant cells. Moreover, the algorithms used in our investigation allowed us to properly annotate the neurons and the immune cells (i.e., macrophages, dendritic and T cells) as non-cancer cells. However, we also identified a distinct cluster of cells that we recognised as “cancer” because of genetic aberrations parallel to CD63 basophil marker expression, while not being positive for other canonical basophil markers (e.g., CD123). That suggests that these cells should not be classified as immune cells [27].
On the technical side, our primary consideration was to distinguish individual cells. Fluorescent signals from cell nuclei and OLIG2 were easy to partition between cells, as in the vast majority of cases, it was single point luminescence. On the other hand, DCX fluorescence was more challenging to evaluate because of the branched structure of a cell’s cytoskeleton. Finally, the hardest to assess was NES: not only were the signal figures branched, but NES+ cells were also more prevalent, with a tendency to be clumped.
In conclusion, our findings support the notion that DCX is indeed expressed by glioma cells, but there is no clear evidence to prove that it may participate in GBM cell migration. Other GSC markers: NES and OLIG2 are in much higher abundance and are present both in tumours and their margin.

4. Materials and Methods

4.1. Samples Collection and Selection of Regions of Interest (ROIs)

Paraffin-embedded samples were obtained from the Department of Neurosurgery and Paediatric Neurosurgery archive, Medical University of Lublin (Poland). We selected only those samples subjected to at least two resections, and at least one sample from each patient was diagnosed with GBM. Moreover, we have focused on case studies since 2011. For each of the samples, a trained neuropathologist performed second-hand diagnosis, along with marking ROIs. The goal was to mark regions either with representative tumour tissue or with the border between tumour and margin tissue. In total, we collected and assessed 17 cases over 46 paraffin-embedded tissues, marking 60 ROIs.
Additionally, all ROIs were classified using two types of categories: tissue origin and surgery sequence number. Tissue origin contains four subcategories: GBM cases (1) when a tumour was adjacent to margin, (2) tissue adjacent to margin, (3) tumour-only when tumour ROI was far from margin/margin was not detected, and (4) tumour-only non-GBM, with ROIs containing a lower-grade astrocytoma. The second category showed surgery sequence number corresponding to samples taken during the first, second or third surgery.

4.2. Multiplex Immunofluorescence Staining

Paraffin-embedded tissues were cut into three µm sections and placed on glass slides (Thermo Fisher Scientific, Waltham, MA, USA: 10149870). Slides were baked overnight at 60 °C. The next day, slides were deparaffinised and hydrated in a series of xylene and ethyl alcohol. Antigens were retrieved by microwave-HIER in a citrate buffer (pH 6.0) over 20 min, followed by another 20 min in RT to cool down. Endogenous peroxide activity was quenched by a mixture of 1% H2O2 and 1% sodium azide for 20 min. All non-specific binding sites were blocked by incubation with blocking solution: PBS (Merck, Darmstadt, Germany: P4417-50TAB) with 1% BSA (Jackson ImmunoResearch, Cambridgeshire, UK: 001-000-161), 5% NDS (Cambridgeshire, UK: 017-000-121) and 1% Triton X-100, for 1 h. OLIG2, NES, and DCX were marked sequentially. Firstly, sections were incubated with appropriate primary antibodies (anti-NES: 1:100, Novus Biologicals, Centennial, CO, USA, MAB1259; anti-OLIG2: 1:100, AF2418; anti-DCX: 1:100, Cambridge, United Kingdom, ab18723; diluted in PBS with 1% BSA, 1% Triton X-100). After washing, corresponding secondary antibodies conjugated with HRP were used (anti-Mouse–Rabbit–Goat: 1:500, Cambridgeshire, UK, 715-036-150, 711-036-152, 705-036-147; diluted in PBS with 1% BSA, 1% Triton X-100). The visualisation was carried out using the Tyramide Signal Amplification (Thermo Fisher Scientific, Waltham, MA, USA: B40953, B40958, B40955) setup. In between fluorescent stainings, HRP activity was quenched by a mixture of 1% H2O2 and 1% sodium azide for 20 min to not interfere with subsequent staining. Finally, nuclei were detected with Hoechst-33342 (Thermo Fisher Scientific, Waltham, MA, USA: H3570), and slides were mounted (Thermo Fisher Scientific, Waltham, MA, USA: P36982). During incubations, slides were covered with a hand-cut parafilm piece, and all washing between steps was conducted with PBS 3 × 5 min. Incubations with primary and secondary antibodies lasted 1 h each.

4.3. Image Acquisition

All images were captured using a Nikon Ti Confocal microscope, using four lasers for fluorescence: 405, 488, 563, and 647 nm. The NISelements (ver 3.22.08, Melville, NY, USA) software was used to set up analyses. The corresponding primary antibody’s negative control was used to control non-specific staining for each batch of analysed slides. In this regard, specific laser intensity and gain parameters were chosen to remove any signal corresponding to 488, 563 and 647 nm. In addition, epifluorescence site conformation was done before confocal imaging to ensure the highest quality of the images.

4.4. Image Analysis

If possible, three different random square spots with a side of 600 pixels from each tissue image were selected to perform further analysis. In total, 152 random square spots were generated across 60 ROIs.
First, stack images were divided by channels corresponding to the wavelength used to detect each antigen: NES, OLIG2, DCX, and nuclei. Next, all features were marked, applying the following three heuristic assumptions:
-
Nucleus—any oval and coherent figure with a signal pattern corresponding to Hoechst-33342;
-
DCX/NES—nucleus with an adjacent signal corresponding to DCX or NES; in case of a signal adjacent to more than one nucleus, all were counted as positive;
-
OLIG2—nucleus overlapping with a signal corresponding to OLIG2.
A nuclei drove Voronoi’s diagram was created to find all positive cells for each of the fluorescent signals. Finally, DCX/NES/OLIG2 layers with marks were used to find and count all overlaps with cells in Voronoi’s diagram. The analysis was carried out in ImageJ (Fiji) [28].

4.5. Image Deposition and Sharing

We set up a local data share service based on Digital Slide Archive (DSA) [29]. All images in TIFF format were converted to pyramidal TIFF format using the ImageMagick tool [30] and uploaded to DSA along with all necessary metadata.

4.6. Single-Cell RNA-Seq and Data Processing

The data used here were generated and published previously by Couturier et al. [26]. It was further processed using the CellRanger pipeline. We included data for 12 CRC patients (originally named as OPK333B, OPK338B, OPK346B, OPK363, OPK364B, OPK368B, OPK389B, OPK390, OPK397, OPK402B, OPK407, OPK409) and three samples derived from the foetal brain (HFA567, HFA570, HFA571). First, the raw gene expression matrix was filtered and normalised using the Seurat R package. Then, the dataset was filtered according to the following criteria: cells with >1000 unique molecular identifier (UMI) counts; >500 genes and <5000 genes; and <10% of mitochondrial gene expression in UMI counts. The gene expression matrices from 15,000 randomly down-sampled cells were log-normalised to the total UMI counts per cell, scaled and finally clustered and visualised using t-SNE projection. The major cell types were characterised by comparing the canonical marker genes found in tissue-specific cell taxonomy reference database CellMatch and the differentially expressed genes for each cluster using the scCATCH automatic annotation algorithm [31].
To identify evidence for somatic large-scale chromosomal copy number alterations, we used inferCNV of the Trinity CTAT Project algorithm [32] in the reference with non-cancer, foetal brain cells. The approach used to annotate the cells as malignant (‘cancer’), was based on the presence of the most common canonical aberration for GBM (duplication of chromosome 7 or chromosome 10 loss), and it is line with the analytical strategy used previously by Couturier et al. [26].

4.7. Statistics

Data files from image analysis were imported and analysed in Rstudio, along with ggplot, corrplot, and dplyr packages. Pairwise comparisons using the Wilcoxon rank-sum test was used to check the cell count in different regions of samples. The Wilcoxon signed-rank test was used to find any consequential difference between margin and adjacent tumour tissue in terms of any staining. Mann–Whitney-Wilcoxon Test compared any other tumour tissue to the ones that were adjacent to the margin. The Kruskal–Wallis rank-sum test allowed to check any change in staining after subsequent tumour resections. Spearman’s rank correlation coefficient measured correlations. To check overall survival, Cox proportional-hazards models were used. In all tests, p < 0.05 was considered significant. All subsamples of the same image were averaged.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/ijms222413217/s1.

Author Contributions

Conceptualisation, R.R., B.J. and A.O.; Methodology, A.O., M.K. and I.T.; Software, A.O.; Validation, R.R., B.J. and A.O.; Formal Analysis, A.O. and M.K.; Investigation, A.O., B.J., M.K., I.T., and S.K.; Resources, B.J., and D.Z.; Data Curation, R.R. and A.O.; Writing—Original Draft Preparation, A.O., and M.K.; Writing—Review and Editing, R.R., and B.J.; Visualization, A.O., M.K.; Supervision, R.R.; Project Administration, R.R.; Funding Acquisition, A.O. and R.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Polish National Science Centre, grant number 2015/17/N/NZ4/02839.

Institutional Review Board Statement

The study material was collected for histopathological diagnosis and came from the archives of Independent Public Clinical Hospital No. 4 in Lublin, Poland. In accordance with article 26(4) of the Act on Patient’s Rights and Patient’s Ombudsman of 2008, November the 6th (i.e., Journal of Laws of 2016, item 186 as amended) medical records can also be shared with universities or research institutes to be used for scientific purposes without revealing surnames or other data enabling identification of the person to which it relates. According to the judgment of the Supreme Court number VCSK 256/10, tissue sections from human organism and histopathological specimen constitute medical records. Thereby, tissue sections and specimen can be shared with medical colleges according to the law under conditions other than patient’s consent. Therefore, no Institutional Review Board Statement is required under Polish Law.

Informed Consent Statement

The study material was collected for histopathological diagnosis and came from the archives of Independent Public Clinical Hospital No. 4 in Lublin, Poland. In accordance with article 26(4) of the Act on Patient’s Rights and Patient’s Ombudsman of 2008, November the 6th (i.e., Journal of Laws of 2016, item 186 as amended) medical records can also be shared with universities or research institutes to be used for scientific purposes without revealing surnames or other data enabling identification of the person to which it relates. According to the judgment of the Supreme Court number VCSK 256/10, tissue sections from human organism and histopathological specimen constitute medical records. Thereby, tissue sections and specimen can be shared with medical colleges according to the law under conditions other than patient’s consent. Therefore, no Informed Consent from the patients is required under Polish Law.

Data Availability Statement

Single-cell RNA-seq data (EGAD00001006206) can be accessed via European Genome-Phenome Archive. A.O. acquired written permission to use mentioned data.

Acknowledgments

We would like to thank Ewa Hordyjewska-Kowalczyk for allowing us to use the workstation funded by the Polish National Science Centre [2017/27/N/NZ5/02940 to EHK].

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Louis, D.N.; Perry, A.; Reifenberger, G.; von Deimling, A.; Figarella-Branger, D.; Cavenee, W.K.; Ohgaki, H.; Wiestler, O.D.; Kleihues, P.; Ellison, D.W. The 2016 World Health Organization Classification of Tumors of the Central Nervous System: A summary. Acta Neuropathol. 2016, 131, 803–820. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Park, J.; de Lomana, A.; Marzese, D.; Juarez, T.; Feroze, A.; Hothi, P.; Cobbs, C.; Patel, A.; Kesari, S.; Huang, S.; et al. A Systems Approach to Brain Tumor Treatment. Cancers 2021, 13, 3152. [Google Scholar] [CrossRef] [PubMed]
  3. Zhang, H.; Wang, R.; Yu, Y.; Liu, J.; Luo, T.; Fan, F. Glioblastoma Treatment Modalities besides Surgery. J. Cancer 2019, 10, 4793–4806. [Google Scholar] [CrossRef]
  4. Vollmann-Zwerenz, A.; Leidgens, V.; Feliciello, G.; Klein, C.A.; Hau, P. Tumor Cell Invasion in Glioblastoma. Int. J. Mol. Sci. 2020, 21, 1932. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Petrecca, K.; Guiot, M.-C.; Panet-Raymond, V.; Souhami, L. Failure pattern following complete resection plus radiotherapy and temozolomide is at the resection margin in patients with Glioblastoma. J. Neuro-Oncol. 2013, 111, 19–23. [Google Scholar] [CrossRef]
  6. Alvarez-Buylla, A.; Temple, S. Stem cells in the developing and adult nervous system. J. Neurobiol. 1998, 36, 105–110. [Google Scholar] [CrossRef]
  7. Gage, F.H. Mammalian Neural Stem Cells. Science 2000, 287, 1433–1438. [Google Scholar] [CrossRef] [PubMed]
  8. Verhaak, R.G.W.; Hoadley, K.A.; Purdom, E.; Wang, V.; Qi, Y.; Wilkerson, M.D.; Miller, C.R.; Ding, L.; Golub, T.; Mesirov, J.P.; et al. Integrated Genomic Analysis Identifies Clinically Relevant Subtypes of Glioblastoma Characterized by Abnormalities in PDGFRA, IDH1, EGFR, and NF1. Cancer Cell 2010, 17, 98–110. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  9. Gleeson, J.G.; Allen, K.M.; Fox, J.W.; Lamperti, E.D.; Berkovic, S.; Scheffer, I.; Cooper, E.C.; Dobyns, W.; Minnerath, S.R.; Ross, M.; et al. Doublecortin, A Brain-Specific Gene Mutated in Human X-Linked Lissencephaly and Double Cortex Syndrome, Encodes a Putative Signaling Protein. Cell 1998, 92, 63–72. [Google Scholar] [CrossRef] [Green Version]
  10. Gleeson, J.G.; Lin, P.T.; Flanagan, L.A.; Walsh, C.A. Doublecortin Is a Microtubule-Associated Protein and Is Expressed Widely by Migrating Neurons. Neuron 1999, 23, 257–271. [Google Scholar] [CrossRef] [Green Version]
  11. Tabata, H.; Nakajima, K. Multipolar migration: The third mode of radial neuronal migration in the developing cerebral cortex. J. Neurosci. 2003, 23, 9996–10001. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. LoTurco, J. Doublecortin and a Tale of Two Serines. Neuron 2004, 41, 175–177. [Google Scholar] [CrossRef] [Green Version]
  13. Daou, M.-C.; Smith, T.W.; Litofsky, N.S.; Hsieh, C.C.; Ross, A.H. Doublecortin is preferentially expressed in invasive human brain tumors. Acta Neuropathol. 2005, 110, 472–480. [Google Scholar] [CrossRef] [PubMed]
  14. Masui, K.; Mawatari, S.Y.; Suzuki, S.O.; Iwaki, T. Evaluation of sensitivity and specificity of Doublecortin immunostatining for the detection of infiltrating glioma cells. Brain Tumor Pathol. 2008, 25, 1–7. [Google Scholar] [CrossRef] [PubMed]
  15. Santra, M.; Santra, S.; Roberts, C.; Zhang, R.L.; Chopp, M. Doublecortin induces mitotic microtubule catastrophe and inhibits glioma cell invasion. J. Neurochem. 2009, 108, 231–245. [Google Scholar] [CrossRef] [Green Version]
  16. Mitew, S.; Hay, C.; Peckham, H.; Xiao, J.; Koenning, M.; Emery, B. Mechanisms regulating the development of oligodendrocytes and central nervous system myelin. Neuroscience 2014, 276, 29–47. [Google Scholar] [CrossRef] [PubMed]
  17. Trépant, A.-L.; Bouchart, C.; Rorive, S.; Sauvage, S.; Decaestecker, C.; Demetter, P.; Salmon, I. Identification of OLIG2 as the most specific Glioblastoma stem cell marker starting from comparative analysis of data from similar DNA chip microarray platforms. Tumor Biol. 2014, 36, 1943–1953. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Ortensi, B.; Setti, M.; Osti, D.; Pelicci, G. Cancer stem cell contribution to Glioblastoma invasiveness. Stem Cell Res. Ther. 2013, 4, 18. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Bott, C.J.; McMahon, L.P.; Keil, J.M.; Yap, C.C.; Kwan, K.Y.; Winckler, B. Nestin selectively facilitates the phosphorylation of the Lissencephaly-linked protein Doublecortin (DCX) by cdk5/p35 to regulate growth cone morphology and Sema3a sensitivity in developing neurons. J. Neurosci. 2020, 40, 3720–3740. [Google Scholar] [CrossRef]
  20. Rich, J.N.; Hans, C.; Jones, B.; Iversen, E.S.; McLendon, R.E.; Rasheed, B.K.A.; Dobra, A.; Dressman, H.K.; Bigner, D.D.; Nevins, J.R.; et al. Gene Expression Profiling and Genetic Markers in Glioblastoma Survival. Cancer Res. 2005, 65, 4051–4058. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  21. Santra, M.; Zheng, X.; Roberts, C.; Santra, S.; Lu, M.; Panda, S.; Jiang, F.; Chopp, M. Single Doublecortin gene therapy significantly reduces glioma tumor volume. J. Neurosci. Res. 2010, 88, 304–314. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Uhlén, M.; Zhang, C.; Lee, S.; Sjöstedt, E.; Fagerberg, L.; Bidkhori, G.; Benfeitas, R.; Arif, M.; Liu, Z.; Edfors, F.; et al. A pathology atlas of the human cancer transcriptome. Science 2017, 357, eaan2507. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Smith, S.J.; Diksin, M.; Chhaya, S.; Sairam, S.; Estevez-Cebrero, M.A.; Rahman, R. The Invasive Region of Glioblastoma Defined by 5ALA Guided Surgery Has an Altered Cancer Stem Cell Marker Profile Compared to Central Tumour. Int. J. Mol. Sci. 2017, 18, 2452. [Google Scholar] [CrossRef] [Green Version]
  24. Cheng, L.; Wu, Q.; Guryanova, O.A.; Huang, Z.; Huang, Q.; Rich, J.N.; Bao, S. Elevated invasive potential of Glioblastoma stem cells. Biochem. Biophys. Res. Commun. 2011, 406, 643–648. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Dirkse, A.; Golebiewska, A.; Buder, T.; Nazarov, P.V.; Muller, A.; Poovathingal, S.; Brons, N.H.C.; Leite, S.; Sauvageot, N.; Sarkisjan, D.; et al. Stem cell-associated heterogeneity in Glioblastoma results from intrinsic tumor plasticity shaped by the microenvironment. Nat. Commun. 2019, 10, 1787. [Google Scholar] [CrossRef]
  26. Couturier, C.P.; Ayyadhury, S.; Le, P.U.; Nadaf, J.; Monlong, J.; Riva, G.; Allache, R.; Baig, S.; Yan, X.; Bourgey, M.; et al. Single-cell RNA-seq reveals that Glioblastoma recapitulates a normal neurodevelopmental hierarchy. Nat. Commun. 2020, 11, 3406. [Google Scholar] [CrossRef]
  27. Kim, Z.; Choi, B.S.; Kim, J.K.; Won, D.I. Basophil Markers for Identification and Activation in the Indirect Basophil Activation Test by Flow Cytometry for Diagnosis of Autoimmune Urticaria. Ann. Lab. Med. 2016, 36, 28–35. [Google Scholar] [CrossRef] [Green Version]
  28. Schindelin, J.; Arganda-Carreras, I.; Frise, E.; Kaynig, V.; Longair, M.; Pietzsch, T.; Preibisch, S.; Rueden, C.; Saalfeld, S.; Schmid, B.; et al. Fiji: An open-source platform for biological-image analysis. Nat. Methods 2012, 9, 676–682. [Google Scholar] [CrossRef] [Green Version]
  29. Gutman, D.A.; Khalilia, M.; Lee, S.; Nalisnik, M.; Mullen, Z.; Beezley, J.; Chittajallu, D.R.; Manthey, D.; Cooper, L.A.D. The Digital Slide Archive: A Software Platform for Management, Integration, and Analysis of Histology for Cancer Research. Cancer Res. 2017, 77, e75–e78. [Google Scholar] [CrossRef] [Green Version]
  30. The ImageMagick Development Team, ImageMagick. 2021. Available online: https://imagemagick.org (accessed on 15 October 2021).
  31. Shao, X.; Liao, J.; Lu, X.; Xue, R.; Ai, N.; Fan, X. scCATCH: Automatic Annotation on Cell Types of Clusters from Single-Cell RNA Sequencing Data. iScience 2020, 23, 100882. [Google Scholar] [CrossRef] [Green Version]
  32. Tickle, T.; Tirosh, I.; Georgescu, C.; Brown, M.; Haas, B. inferCNV of the Trinity CTAT Project; Klarman Cell Observatory, Broad Institute of MIT and Harvard: Cambridge, MA, USA, 2019; Available online: https://github.com/broadinstitute/inferCNV (accessed on 15 October 2021).
Figure 1. Staining patterns. (A)—Nuclei (yellow), (B)—OLIG2 signal (red), (C)—NES signal (green), (D)—DCX signal (magenta), (E)—NES signal with nuclei, (F)—OLIG2 signal with nuclei, (G)—DCX signal with nuclei (H)—composition of (AD).
Figure 1. Staining patterns. (A)—Nuclei (yellow), (B)—OLIG2 signal (red), (C)—NES signal (green), (D)—DCX signal (magenta), (E)—NES signal with nuclei, (F)—OLIG2 signal with nuclei, (G)—DCX signal with nuclei (H)—composition of (AD).
Ijms 22 13217 g001
Figure 2. Overview of samples’ staining and staining overlap. Frequency of cells with different stainings with respect to tissues: (A)—DCX+ cells, (B)—NES+ cells, (C)—OLIG2+ cells. Frequency of cells with mixtures of all stainings with respect to tissues: (D)—OLIG2+DCX+ cells, (E)—OLIG2+NES+ cells, (F)—NES+DCX+ cells, (G)—DCX+OLIG2+NES+ cells. (H)—Total number of cells irrespective of antigen staining.
Figure 2. Overview of samples’ staining and staining overlap. Frequency of cells with different stainings with respect to tissues: (A)—DCX+ cells, (B)—NES+ cells, (C)—OLIG2+ cells. Frequency of cells with mixtures of all stainings with respect to tissues: (D)—OLIG2+DCX+ cells, (E)—OLIG2+NES+ cells, (F)—NES+DCX+ cells, (G)—DCX+OLIG2+NES+ cells. (H)—Total number of cells irrespective of antigen staining.
Ijms 22 13217 g002
Figure 3. Overview of samples underwent subsequent resections. Frequency of cells with different stainings for tissues: (A)—DCX+ cells, (B)—NES+ cells, (C)—OLIG2+ cells. Frequency of cells with mixtures of all stainings for tissues: (D)—OLIG2+DCX+ cells, (E)—OLIG2+NES+ cells, (F)—NES+DCX+ cells, (G)—DCX+OLIG2+NES+ cells.
Figure 3. Overview of samples underwent subsequent resections. Frequency of cells with different stainings for tissues: (A)—DCX+ cells, (B)—NES+ cells, (C)—OLIG2+ cells. Frequency of cells with mixtures of all stainings for tissues: (D)—OLIG2+DCX+ cells, (E)—OLIG2+NES+ cells, (F)—NES+DCX+ cells, (G)—DCX+OLIG2+NES+ cells.
Ijms 22 13217 g003
Figure 4. tSNE projection of single cells. (A)—Cancer and normal cells, (B)—All the cells coloured by sample name, (C)—“Cancer” cells only coloured by sample name, (D)—“Normal” cells only coloured by sample name.
Figure 4. tSNE projection of single cells. (A)—Cancer and normal cells, (B)—All the cells coloured by sample name, (C)—“Cancer” cells only coloured by sample name, (D)—“Normal” cells only coloured by sample name.
Ijms 22 13217 g004
Figure 5. tSNE projection of annotated single cells by expression of markers specific for various cell types. (A)—Cancer and normal cells, (B)—Cancer cells, (C)—Normal cells.
Figure 5. tSNE projection of annotated single cells by expression of markers specific for various cell types. (A)—Cancer and normal cells, (B)—Cancer cells, (C)—Normal cells.
Ijms 22 13217 g005
Figure 6. The mean expression of OLIG2, NES, and DCX transcripts among the annotated clusters of single cells that express markers for each cell type. (A) the cells annotated as “cancer”, (B) the cells annotated as “normal”.
Figure 6. The mean expression of OLIG2, NES, and DCX transcripts among the annotated clusters of single cells that express markers for each cell type. (A) the cells annotated as “cancer”, (B) the cells annotated as “normal”.
Ijms 22 13217 g006
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Odrzywolski, A.; Jarosz, B.; Kiełbus, M.; Telejko, I.; Ziemianek, D.; Knaga, S.; Rola, R. Profiling Glioblastoma Cases with an Expression of DCX, OLIG2 and NES. Int. J. Mol. Sci. 2021, 22, 13217. https://doi.org/10.3390/ijms222413217

AMA Style

Odrzywolski A, Jarosz B, Kiełbus M, Telejko I, Ziemianek D, Knaga S, Rola R. Profiling Glioblastoma Cases with an Expression of DCX, OLIG2 and NES. International Journal of Molecular Sciences. 2021; 22(24):13217. https://doi.org/10.3390/ijms222413217

Chicago/Turabian Style

Odrzywolski, Adrian, Bożena Jarosz, Michał Kiełbus, Ilona Telejko, Dominik Ziemianek, Sebastian Knaga, and Radosław Rola. 2021. "Profiling Glioblastoma Cases with an Expression of DCX, OLIG2 and NES" International Journal of Molecular Sciences 22, no. 24: 13217. https://doi.org/10.3390/ijms222413217

APA Style

Odrzywolski, A., Jarosz, B., Kiełbus, M., Telejko, I., Ziemianek, D., Knaga, S., & Rola, R. (2021). Profiling Glioblastoma Cases with an Expression of DCX, OLIG2 and NES. International Journal of Molecular Sciences, 22(24), 13217. https://doi.org/10.3390/ijms222413217

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

Article Metrics

Back to TopTop