Next Article in Journal
Osteogenic Commitment of Human Periodontal Ligament Cells Is Predetermined by Methylation, Chromatin Accessibility and Expression of Key Transcription Factors
Next Article in Special Issue
Transcriptomic Profiling of DNA Damage Response in Patient-Derived Glioblastoma Cells before and after Radiation and Temozolomide Treatment
Previous Article in Journal
Pregnane X Receptor Mediates Atherosclerosis Induced by Dicyclohexyl Phthalate in LDL Receptor-Deficient Mice
Previous Article in Special Issue
Combining HDAC and MEK Inhibitors with Radiation against Glioblastoma-Derived Spheres
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cells
Volume 11
Issue 7
10.3390/cells11071127
cells-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Single-Cell Molecular Characterization to Partition the Human Glioblastoma Tumor Microenvironment Genetic Background

by
Francesca Lessi
1,*,
Sara Franceschi
1,
Mariangela Morelli
1,
Michele Menicagli
1,
Francesco Pasqualetti
2,3,
Orazio Santonocito
4,
Carlo Gambacciani
4,
Francesco Pieri
4,
Filippo Aquila
4,
Paolo Aretini
1 and
Chiara Maria Mazzanti
1
1
Section of Genomics and Transcriptomics, Fondazione Pisana per la Scienza, San Giuliano Terme, 56017 Pisa, Italy
2
Department of Radiation Oncology, Azienda Ospedaliera Universitaria Pisana, University of Pisa, 56126 Pisa, Italy
3
Department of Oncology, University of Oxford, Oxford OX1 2JD, UK
4
Division of Neurosurgery, Spedali Riuniti di Livorno—USL Toscana Nord-Ovest, 57124 Livorno, Italy
*
Author to whom correspondence should be addressed.
Cells 2022, 11(7), 1127; https://0-doi-org.brum.beds.ac.uk/10.3390/cells11071127
Submission received: 4 March 2022 / Revised: 17 March 2022 / Accepted: 24 March 2022 / Published: 26 March 2022
(This article belongs to the Special Issue Molecular and Cellular Mechanisms of Cancers: Glioblastoma II)
Browse Figures
Review Reports Versions Notes

Abstract

:
Background: Glioblastoma (GB) is a devastating primary brain malignancy. The recurrence of GB is inevitable despite the standard treatment of surgery, chemotherapy, and radiation, and the median survival is limited to around 15 months. The barriers to treatment include the complex interactions among the different cellular components inhabiting the tumor microenvironment. The complex heterogeneous nature of GB cells is helped by the local inflammatory tumor microenvironment, which mostly induces tumor aggressiveness and drug resistance. Methods: By using fluorescent multiple labeling and a DEPArray cell separator, we recovered several single cells or groups of single cells from populations of different origins from IDH-WT GB samples. From each GB sample, we collected astrocytes-like (GFAP+), microglia-like (IBA1+), stem-like cells (CD133+), and endothelial-like cells (CD105+) and performed Copy Number Aberration (CNA) analysis with a low sequencing depth. The same tumors were subjected to a bulk CNA analysis. Results: The tumor partition in its single components allowed single-cell molecular subtyping which revealed new aspects of the GB altered genetic background. Conclusions: Nowadays, single-cell approaches are leading to a new understanding of GB physiology and disease. Moreover, single-cell CNAs resource will permit new insights into genome heterogeneity, mutational processes, and clonal evolution in malignant tissues.

1. Introduction

Glioblastoma (GB) is the most aggressive and deadly primary tumor of the central nervous system in adults with an overall survival of fewer than 15 months [1]. The extremely poor prognosis of GB, despite the development in recent decades of new and innovative therapies, is enhanced by the resistance developed towards radio and chemotherapy [2]. In this tumor, as well as in other cancer types, the tumor microenvironment (TME) plays a pivotal role in treatment resistance [2]. The GB microenvironment is composed of a massive number of different cells, and besides malignant astrocytes and cancer stem cells, stromal cells, endothelial cells, pericytes, and a huge number of immune cells are present [3]. Moreover, intratumoral heterogeneity (ITH), which is one of the major features of GB tumors, is also hugely involved in anticancer treatment resistance [4,5] and is critical to promote tumoral growth and aggressiveness [6]. In support of this last remark, it has recently been demonstrated in GB that different sub-clones co-exist within the same tumor that respond differently to differing therapies [7]. These sub-populations of cells show distinct genomic profiles that reveal an individual behavior peculiar to the whole cell population [8]. Currently, the single-cell approach in GB is becoming increasingly popular. Reaching single-cell resolution enables avoiding the averaging of bulk analysis and capturing the heterogeneity of cells. Copy number aberration (CNA) is one of the most important somatic alterations in cancer [9,10], defined as somatic changes to the chromosome structure such as the gain and/or deletion of a particular DNA segment (>1 kb) [11]. The most common CNAs in GB include the loss, or partial loss, of chromosomes 9 and 10; the gain of chromosomes 7, 19, and 20; the focal deletion of the CDKN2A/B locus (9p21.3); and the focal high-level amplification of the EGFR locus (7p11.2) [12,13]. In particular, it is well known that CNAs targeting chromosomes 7 and 10 are some of the earliest events in GB tumor evolution [14]. The analysis of these aberrations is interesting because CNAs are detected with much greater accuracy than individual mutations and are associated with ITH in most cancers [15]. Moreover, the aggregation of cells sharing the same CNA profiles allows improving the phylogenetic analysis at the single nucleotide level [16].
In this work, we collected three human GB tumors and after dissociation, a certain number of single and groups of single cells were isolated through DEPArray technology, paying particular attention to four cell populations: astrocytes-like, microglia-like cells, endothelial-like cells, and stem-like cells. Afterwards, we investigated the genomic aberrations (CNA analysis) in these different types of tumor cells, thus performing a single-cell CNA analysis. The whole parental tumors were subjected to a bulk CNA analysis as well, to compare their molecular profiles with the single-cell results. The tumor partition in its single components allowed single-cell molecular profiling which revealed new aspects of the GB altered genetic background. Our work demonstrates that the single-cell approach is more representative and detailed than the bulk analysis, which contributes to a deeper insight into the basic molecular mechanisms of GB. Moreover, we presented an innovative approach to isolate and characterize different tumor populations of cells at the single-cell level.

2. Materials and Methods

2.1. Human Glioblastoma Tissue Collection

The study has been performed according to the Declaration of Helsinki and the samples’ collection protocol was approved by the Ethics Committee of the University Hospital of Pisa (787/2015). Tumor tissues were obtained from patients who underwent surgical resection of histologically confirmed GB after informed consent. Samples were obtained from the Unit of Neurosurgery of Livorno Civil Hospital. Three patient cases (GB01, GB02, and GB03) were included in the present study, the clinical and demographic data and the pathological and therapeutical information are summarized in Table 1. All cases had a diagnosis of GB with no previous history of any brain neoplasia and were not carrying R132 IDH1 or R172 IDH2 mutations. Surgically resected tumors were collected and stored in MACS tissue storage solution (Miltenyi Biotec, Bergisch Gladbach, Germany) at 4 °C for 2–4 h. Each tumor sample was washed with Dulbecco’s phosphate-buffered saline (DPBS) in a sterile dish and portioned with a scalpel into about 0.5–2 cm2 pieces under a biological hood. Afterward, they were vital frozen at −140 °C in 90% fetal bovine serum (FBS) and 1% dimethyl sulfoxide (DMSO) for further analyses.

2.2. Tumor Dissociation to Single-Cell Suspensions

Frozen GB tissues were defrosted in a water bath at 37 °C, washed with DPBS in a sterile dish and cut with a scalpel into small pieces. We used 0.11 gr, 0.16 gr, and 0.14 gr of GB01, GB02, and GB03 respectively. These finely minced tumor chunks were transferred in a C-tube (Miltenyi Biotec, Bergisch Gladbach, Germany) with the appropriate volume of buffer X following the protocol (Brain Tumor Dissociation Kit, Miltenyi Biotec, Bergisch Gladbach, Germany) for tumor dissociation with the gentleMACs Dissociator (Miltenyi Biotec, Bergisch Gladbach, Germany).

2.3. Immunofluorescence of Single-Cell Suspensions

The cell suspensions obtained were transferred to 1.5 mL LoBind tubes and washed three times with DPBS. After centrifugation at 300× g for 10 min at room temperature, the supernatant was removed, and the cells were resuspended in 400 µL of running buffer composed of MACS BSA stock solution (Miltenyi Biotec, Bergisch Gladbach, Germany) 1:20 with autoMACS Rinsing Solution (Miltenyi Biotec, Bergisch Gladbach, Germany). The cells were fixed by adding 400 µL of paraformaldehyde 4%; cells were incubated with fixation solution for 20 min at room temperature. To stop the reaction, the sample tubes were filled with DPBS and centrifuged at 400× g for 5 min at room temperature. Afterward, we performed two washes with DPBS to the sample tubes and then we incubated the pellet with blocking solution for 10 min at room temperature (BSA 3% in DPBS). The blocking reaction was stopped by filling the tube with DPBS, before centrifugation at 400× g for 5 min at room temperature. The cells were resuspended in running buffer and counted with a Luna Automated Cell Counter (Logos Biosystems, Anyang-si, Gyeonggi-do, Korea). For the immunofluorescence, a maximum of 100,000 fixed cells was used for the staining. The antibodies chosen for the staining were: anti-GFAP APC (130-124-040, Miltenyi Biotec, Bergisch Gladbach, Germany) for astrocytes, anti-IBA1 PE (ab209942, Abcam, Cambridge, UK) for microglia/macrophages cells, anti-CD105 PerCP/Cy5.5 (ab234265, Abcam, Cambridge, UK) for endothelial cells, anti-CD133 FITC (11-1339-42, eBioscience, San Diego, CA, USA) for stem cells, and Hoechst 33342 (62249, Thermofisher Scientific, Waltham, MA, USA) for nuclei. A total of 20 µL of anti-CD105 and 25 µL of anti-CD133 were added to the cell suspensions and mixed by gently pipetting. The samples were incubated for 15 min in the dark at 4 °C. The reaction was stopped by adding 1 mL of running buffer and mixed by gently pipetting. Then, the sample tubes were centrifuged at 400× g for 10 min at room temperature, the supernatant was removed, and the cells were resuspended with 100 µL of Inside Perm Buffer (Inside Stain Kit, Miltenyi Biotec, Bergisch Gladbach, Germany). A total of 8 µL of anti-GFAP and 2.5 µL of anti-IBA1 were added to the cell suspensions and mixed by gently pipetting. The samples were incubated for 20 min in the dark at room temperature. The reaction was stopped by adding 1 mL of Inside Perm Buffer, and mixed by gently pipetting. Then, the sample tubes were centrifuged at 400× g for 10 min at room temperature, the supernatant was removed, and the cells were resuspended with 1 mL of running buffer. Then, 1 µL of Hoechst (1 mg/mL) was added to the sample tubes and mixed by gently pipetting. The samples were incubated for 5 min in the dark at room temperature. Then, the sample tubes were centrifuged at 400× g for 10 min at room temperature and resuspended in 200 µL of running buffer.

2.4. Single-Cell Isolation by DEPArrayTM NxT

Single cells were isolated and sorted with DEPArray NxT (Menarini, Silicon Biosystems, Bologna, Italy). After the immunofluorescence of the single cell suspensions was measured, the cells were counted; we used a maximum of 24,000 cells to load the DEPArray NxT Cartridge. The samples were washed two times with 1 mL of SB115 Buffer (Menarini, Silicon Biosystems, Bologna, Italy) and the cells were loaded onto the DEPArray NxT cartridge following the protocol instructions. CellBrowser™ (Menarini, Silicon Biosystems, Bologna, Italy) analysis software, integrated into the DEPArray™ system, allows the user to view and select cells from the particle database according to multiple criteria, based on qualitative and quantitative marker evaluation and cell morphology. This software enables the user to create populations and sub-populations of cells using analysis tools such as scatter plots, histograms, and image panels. Cells become un-routable based on their positions; when these are out of the cage, it is no longer possible to move them and therefore complete the recovery. First of all, we excluded clusters of two or three cells, clumps, and spurious events and focused only on single cells with the desired fluorescence, analyzing only the “centered” DAPI cells in the cage. The single cells were selected manually based on fluorescence labeling and morphology. About 20 different single cells were recovered for each tumor patient and volume reduction was performed with a VRNxT-Volume Reduction Instrument (Menarini, Silicon Biosystems, Bologna, Italy) according to the instruction manual. The isolated cells were stored at −20 °C until later downstream analysis.

2.5. Immunofluorescence of GB Tissues

Formalin-fixed, paraffin-embedded (FFPE) tissue blocks, obtained from our GB samples, were cut into 2–4 μm thick sections. Antigen unmasking was achieved with Epitote Retrieval Solution (pH = 8) (Leica Microsystems, Wetzlar, Germany) in a microwave. GFAP monoclonal (ASTRO6) (MA5-12023, Thermofisher Scientific, Waltham, MA, USA) and IBA1 polyclonal (091-19741, Wako Chemicals, Neuss, Germany) primary antibodies were then applied at dilutions of 1:100 and 1:1000, respectively, overnight at 4 °C. The goat anti-rabbit Alexa Fluor 488 (Thermofisher Scientific, Waltham, MA, USA ) and goat anti-mouse Alexa Fluor 568 (Thermofisher Scientific, Waltham, MA, USA ) were diluted 1:500 and incubated for 1 h. Cells were counterstained with DAPI (Sigma Aldrich, St. Louis, MO, USA) and visualized using an Olympus Fluoview 3000 confocal microscope at a magnification of 60×.

2.6. DNA Extraction from Fresh Tissues

Genomic DNA was extracted directly from up to 50 mg of fresh tissue of GB01, GB02, and GB03 using the Maxwell® 16 Instrument with the Maxwell® 16 Tissue DNA Purification Kit (Promega, Madison, WI, USA). DNA concentration was determined using the Qubit Fluorometer (Life Technologies, Carlsbad, CA, USA) and the quality was assessed using the Agilent 2200 Tapestation (Agilent Technologies, Santa Clara, CA, USA) system.

2.7. Ampli1™ Whole Genome Amplification and Low Pass Analysis

Whole-genome amplification on all recovered single cells was performed using the Ampli1™ WGA Kit version 02 (Menarini, Silicon Biosystems, Bologna, Italy) following the manufacturer’s instructions. The same procedure was adjusted for the DNA obtained from fresh tissues starting from 1 µL of 1 ng/µL. Afterward, the WGA product was cleaned up with SPRIselect Beads (Beckman Coulter, Brea, CA, USA) and sequencing-ready libraries were prepared with an Ampli ™ LowPass Kit (Menarini, Silicon Biosystems, Bologna, Italy) to detect chromosomal aneuploidies and copy number aberrations (CNAs) with a low sequencing depth. To sequence our libraries, we used an Ion 520/530-OT2 kit (Ion Torrent, Life Technologies, Grand Island, NY, USA) with the Ion 530 Chip (Ion Torrent, Life Technologies, Grand Island, NY, USA). The runs were conducted on the Ion S5 system (Ion Torrent, Life Technologies, Grand Island, NY, USA).

2.8. CNA Calling

The data obtained from low-pass whole genome sequencing were processed with the IchorCNA tool [17]. The CNA segmented number profiles obtained from IchorCNA were processed with the CNApp tool [18] with default cutoffs.

3. Results

3.1. Isolation of Single-Cells from GB Fresh Tissues with DEPArrayTM NxT

Three GB fresh tissues obtained from the Unit of Neurosurgery of Livorno Civil Hospital were analyzed with DEPArrayTM NxT, the overview of the procedure is shown in Figure 1 in which H&E images for each tumor tissue are also present. After DEPArray NxT Cartridge loading, we selected the routable cells using the CellBrowser™ analysis software. In detail, for GB01, 2880 routable cells, for GB02, 17,378 routable cells, and for GB03, 4788 routable cells, were observed. After that, we performed the exclusion of cell clusters obtaining single and routable cells: 2654, 9535, and 4278 cells respectively for GB01, GB02, and GB03.

Cell Populations in GB01, GB02, and GB03

We chose four different conventional markers to identify the four most representative subpopulations of GB (astrocytes, microglia, endothelial cells, and stem cells): GFAP, IBA1, CD105, and CD133. We decided to call them: astrocyte-like, microglia-like, endothelial-like, and stem-like cells because of their similarity to these particular cells. Moreover, we found some cells with double fluorescence staining.
An example of the main populations is shown in Figure 2. In Figure 3, percentages of the main populations found in the three samples are summarized, while in the Supplementary Materials (Figure S1), double fluorescence stained cells and unlabeled cells are shown. We recovered both single cells and groups of a maximum of five single cells with the same characteristics. The recovered cells for the three samples are summarized in Figure 4. In particular, for GB01 we selected 20 cells: 3 single astrocytes-like, 3 groups of astrocytes-like, 4 microglia-like single cells, 2 groups of microglia-like cells, 1 group of endothelial-like cells, 1 single stem-like cell, 2 single astrocytes/microglia-like cells (positive for both GFAP and IBA1), and 3 single cells and 1 group of single cells without labeling (positive to Hoechst 33342 only). For GB02, 26 cells were recovered: 6 single astrocytes-like, 5 microglia-like single cells, 5 single endothelial-like cells, 3 groups of endothelial-like cells, 5 single stem/endothelial-like cells (positive both for CD133 and CD105), 1 group of stem/endothelial-like cells (positive both for CD133 and CD105), and 1 single cell without labeling (Hoechst 33342 signal only). Finally, for GB03, 17 cells were selected: 6 single astrocytes-like, 5 single microglia-like cells, and 6 single endothelial-like cells.

3.2. Copy Number Aberrations (CNAs) Analysis

3.2.1. GB Bulk Tissues

Cellular genomic profiling was performed on the selected cells using the Ampli1™ LowPass kit to identify genome-wide CNAs at the single-cell level and to obtain information on ITH. The same analysis was also carried out on the DNA obtained from fresh tumor tissues (GB01, GB02, and GB03), to compare the bulk molecular profile to the one derived from single cells. In Figure 5, the CNA pattern of the fresh GB tissues is shown: as expected, each sample has a different CNA configuration due to GB ITH. However, all three samples presented chromosome 10q deletion, and GB01 and GB02 also presented chromosome 7 amplification, which represent typical GB alterations. Consequently, for each sample, tumors in bulk and single-cell CNAs were compared.

3.2.2. GB Single Cells

The summarized results obtained from CNA analysis on single cells are described in Table 2, Table 3 and Table 4. In GB01, we found a group of wild type endothelial-like cells; of these, there were six microglia-like cells (four single cells and two groups of cells), two were wild type (one single cell and one group of cells), one cell showed a chr 19 deletion only, and the other cells showed different alterations, sharing a chr 10 deletion, and chr 7, 9q, and 17q amplifications; six astrocytes-like (three single cells and three groups of cells) were altered, sharing a chr 10 deletion, and chr 7, 9q, and 17q amplifications; one stem-like cell with a chr 1p and 10 deletion and chr 7, 9, 17q, and 19q amplifications. In GB01, moreover, two cells with double staining (GFAP and IBA1) were found with the same alterations, chr 1p and 10 deletions and chr 7, 9, 17q, and 19q amplifications. Finally, three out of four unstained cells (three single and one group of cells) were wild type and one exhibited chr 1p, 10, and 17p deletion and chr 7, 9q, 17q, and 19q amplifications. In GB02, we found eight endothelial-like cells (five single and three groups of cells), two were wild type and the others carried a chr 19 deletion except for only one having chr 9p, 10, 13q, 14q, and 22q deletions. Then, of six single astrocytes-like cells, one was wild type, one had a chr 19 deletion, and the others shared chr 9p, 10, 13q, 14q, and 22q deletions. Indeed, of five single microglia-like cells, two were wild type and three had a chr 19 deletion. In GB02, we selected six double staining cells (CD133 and CD105 positive), of these, four were wild type and the others shared chr 10, 13q, 14q, and 22q deletions. Finally, one unstained cell was wild type. GB03 counted six single endothelial-like cells, four of which were wild type and the other two presented different alterations sharing in particular chr 9p, 10, and 22q deletion and chr 7, 9q, and 20 amplifications. Five single microglia-like cells were all wild type. Finally, six single astrocytes-like were selected, one was wild type while the other cells showed all the same alterations: chr 9p, 10, and 22q deletions and chr 7, 9q, and 20 amplifications.

3.2.3. Comparison between Bulk Tissues and Single Cells

In Figure 6, the comparison between bulk fresh tumor CNAs and single cell CNAs obtained with the CNApp tool is shown. Cells with CNAs have similar alterations to those found in the bulk tissues and also show additional alterations. The molecular alteration profiles in single cells are more strongly highlighted, as in bulk tumors many alterations may be hidden since many cells are analyzed at the same time.

3.3. Double Staining Cells Immunofluorescence

To confirm the presence of double staining cells in our tissues, we performed immunofluorescence with anti-GFAP (red) and anti-IBA1 (green) on our tissues’ slides. We observed some cells with these characteristics in GB01 tissue slides, as shown in Figure 7.

4. Discussion

Despite the new therapies developed in the last few years, GB still remains an incurable and devastating disease [19]. The adjective “multiforme”, often used to define GB, was coined in 1926 by Percival Bailey and Harvey Cushing [20] to describe the various appearances of necrosis, cysts, and hemorrhage. As a matter of fact, this definition also fits from a molecular point of view to explain the high degree of heterogeneity in GB. The poor prognosis of GB patients is mainly associated with ITH, which represents the presence in the tumor mass of multiple sub-clones, each characterized by different molecular and genomic alterations [21]. The sub-clones’ alterations are certainly masked during bulk tumor analysis [22]. There are several approaches to assess the degree of ITH, such as flow cytometry or more innovative methods such as single-cell sequencing and DEPArray analysis. These are certainly three technologies used to decipher ITH, but none can replace the others; rather, they aim to be complementary. Recently, in some single-cell sequencing studies, to investigate the ITH, CNAs investigations were conducted instead of the identification of individual mutations with a gain in sensitivity and accuracy [4,15,23]. Regarding these different techniques, single-cell RNA seq is mainly a discovery analysis: recent single-cell transcriptome studies in GB have made it possible to identify tumor cell populations and to highlight tumor plasticity and hierarchy [24,25]. DEPArray analysis, instead, allows us to select, isolate, and analyze specific cells or groups of cells providing a higher level of precision and accuracy in cell selection than traditional flow cytometry, with a high transfer efficiency and unprecedented purity for molecular analysis. In this work, we decided to focus our attention on some of the most representative GB populations: astrocytes, microglia, stem cells, and endothelial cells. We have assumed that we have isolated the above-mentioned cells, based on the positivity of the chosen markers. However, as only one marker is used per cell, we cannot be sure that we have exactly the hypothesized cell, so we use the suffix “–like” to describe the cells isolated. Through selection and isolation with DEPArray, we investigated the molecular alterations of the isolated cells by comparing them with whole tumor tissue, in terms of CNAs. Astrocytes are star-shaped cells of the brain with different active roles in both healthy people and in brain pathological conditions [26]. For example, they regulate neural signaling and give support in blood-brain barrier (BBB) formation [26]. Regarding GB tumorigenesis, a much-debated topic concerns the cell-of-origin in the cancer stem cell (CSC) or hierarchical hypothesis: GB stem cells (GSCs) or glioma initiating cells seem to be responsible for tumor formation [3]. They are a small population of stem cells characterized by self-renewal and differentiation properties [27]. GSCs are involved in tumor growth, invasion, and recurrence development [28]. Based on this theory, GSCs can arise from neural stem cells [29] but also from already differentiated astrocytes transformed through genetic and epigenetic mutations [30,31]. Therefore, based on this hypothesis, the cell population initiating GB is composed of a mixture of cells including astrocytes and stem cells. In our work, most of the astrocytes-like cells in all three tumors, were altered with a CNA pattern identical to the bulk tumor. In some cells, more alterations were observed than in the bulk, in support of the concept of the higher sensitivity and accuracy of the single-cell analysis approach. Indeed, the only stem-like cell collected in GB02 showed a CNA pattern typical of a transformed tumor cell. This suggests that the cumulative acquisition of mutations in the stem cells can be responsible for invasive cancer generation.
In the brain, microglial cells, a specialized population of macrophage-like cells, represent resident innate immune cells and are involved in many crucial physiological processes [32]. Microglia have been ignored for a long time but by now it is common knowledge that these cells are an integral part of the tumor, constituting approximately 30% of tumor mass [33] and participating in tumor progression and anti-cancer treatment resistance [34]. Indeed, microglial cells have a key role in many brain diseases [35]. From our results, we observed some microglia-like cells with normal chromosomes sets, as we expected, but we also found some cells presenting CNAs, indicating that within the tumor there are also microglia cells with potential tumoral behavior. From a transcriptional point of view, some alterations have been described in GB microglia [36]. In 2020, Maas and colleagues defined a particular type of transformed microglial cells. In this context, tumoral GB cells hijack microglial gene expression to enhance tumor proliferation, suppressing the immune response [37].
Endothelial cells (ECs) represent the principal components of the BBB [38]. Different brain pathologies, including GB, show molecular alterations of ECs [39]. In GB, vessels are necessary for cancer cell spreading and it has been demonstrated that ECs regulate tumor invasion through crosstalk with GB cells [40]. Our results illustrate the presence of wild type endothelial-like cells also carrying CNAs, confirming that the tumor mass can contain tumor-ECs (also defined tumor-associated ECs) as has been highlighted in some recent publications [41,42,43]. In these papers, the tumor-associated ECs showed different phenotypic and functional characteristics concerning normal ECs. Moreover, the relationship between ECs and GB tumor cells was demonstrated in two recent studies, in particular it was observed that tumor-derived ECs and GB stem cells shared the same genomic mutations and that CD144 and VDGFR2 genes are expressed by the emerging endothelium [44,45].
Moreover, in our study, we observed and then recovered some cells with a double signal of labeling: astrocytes/microglia-like cells in GB01 and stem/endothelial-like cells in GB02. Indeed, in the literature, the detection of dual positive cells has been reported in experiments using our same technology, especially in the circulating tumor cell studies [46]. The presence of these double stained cells, in particular GFAP+/IBA1+, was also confirmed by immunofluorescence experiments (Figure S2) to strengthen our findings. Fais et al. in 2007 introduced the concept of cannibalism as an exclusive property of malignant tumor cells [47]. Moreover, Coopman et al. assumed that phagocytosis is the mechanism used by invasive tumor cells to allow migration into the surrounding tissues [48]. In this regard, in malignant gliomas, phagocytic tumor cells were detected, particularly in GB [49,50]. A different hypothesis could be the cell fusion formation, for example, Huysentruyt et al. observed fusion between macrophages and tumor cells [51].
A further aspect that emerged from our results is the detection both in GB01 and in GB02 of some unstained cells with CNAs. We observed, in fact, that not all the astrocytes are positive for GFAP and it has also been demonstrated in the literature that GFAP is not an astrocytes-exclusive marker, as GFAP expression in GB varies significantly [52].
The use of CNAs as a method of evaluating tumor cells is more popular lately. The CNA burden is assessed in different tumors, such as in prostate cancer [53], meaning as the analysis of the variable amounts of amplifications or deletions in different patients. In particular, Hieronymus et al. [53] observed that patients with a high CNA burden showed a greater risk of relapse after treatment. For this reason, CNA analysis can also be considered as a useful marker. Therefore, the tumor CNA burden, rather than individual CNAs, can be associated with cancer outcomes. Recently, CNA analysis has been evaluated as more advantageous than mutational analysis for diagnostic reasons in particular in association with survival [54]: CNAs and miRNA analysis had a better performance than mutational data for poorly predicted survival. In addition, in melanoma, Roh et al. demonstrated that the association of CNAs and the mutational burden can be very useful for prognosis and the response to therapy [55].
To the best of our knowledge, this is the first time that our approach is used to partition a GB tumor tissue into its cellular components and provide its molecular profile. Single-cell CNA analysis has the potential to yield new insights into the molecular dynamics of cellular populations. Measuring single-cell genome alterations in tissues and cell populations will greatly advance the clonal decomposition of malignant tissues, resolving rare cell population genotypes and identifying DNA amplification and the deletion states of individual cells, which are difficult to establish when cellular information is destroyed in bulk sequencing. A novel feature of our approach is also the capture, by brightfield and immunofluorescence imaging, of the morphologic features of cells, permitting analytical integration with genomic properties.

5. Conclusions

In conclusion, in this work, we were able to isolate single cells from fresh GB tissues based on markers that assigned them to the four cell subpopulations: astrocytes, microglia, endothelial cells, and stem cells. CNA analysis allows us to distinguish the tumor cells inside the tumor microenvironment. This is a preliminary work based on an innovative technique, single-cell CNA analysis with DEPArray, to select single tumor cells and study their molecular alterations in depth. This new type of experimental approach is proposed as a complementary procedure to conventional methodologies and provides a baseline for further analyses that aim to explore in depth the different subpopulations in the GB microenvironment. Moreover, the single-cell approach allows a very sensitive analysis rather than bulk analysis, obtaining molecular profiles more accurately. In such an aggressive and lethal tumor, any kind of information is crucial and can be useful to better understand the mechanisms underlying the development of the tumor and its propagation.

Supplementary Materials

The following supporting information can be downloaded at: https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/cells11071127/s1, Figure S1: Pie charts of the double-stained cells; Figure 7: Immunofluorescence assay.

Author Contributions

Conceptualization, C.M.M. and F.L.; methodology, F.L.; software, P.A.; validation, M.M. (Michele Menicagli), M.M. (Mariangela Morelli) and. S.F.; investigation, F.L.; resources, O.S., C.G., F.P. (Francesco Pieri), and F.A.; data curation, F.P. (Francesco Pasqualetti); writing—original draft preparation, F.L.; writing—review and editing, F.L. and S.F.; visualization, M.M. (Mariangela Morelli); supervision, C.M.M.; project administration, F.P. (Francesco Pasqualetti). All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

This study was conducted in accordance with the Declaration of Helsinki, and approved by the Ethics Committee of the University Hospital of Pisa (787/2015).

Informed Consent Statement

Informed consent was obtained from all subjects involved in this study.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Stupp, R.; Hegi, M.E.; Mason, W.P.; van den Bent, M.J.; Taphoorn, M.J.; Janzer, R.C.; Ludwin, S.K.; Allgeier, A.; Fisher, B.; Belanger, K.; et al. Effects of radiotherapy with concomitant and adjuvant temozolomide versus radiotherapy alone on survival in glioblastoma in a randomised phase III study: 5-year analysis of the EORTC-NCIC trial. Lancet Oncol. 2009, 10, 459–466. [Google Scholar] [CrossRef]
  2. Zhang, X.; Ding, K.; Wang, J.; Li, X.; Zhao, P. Chemoresistance caused by the microenvironment of glioblastoma and the corresponding solutions. Biomed. Pharmacother. 2019, 109, 39–46. [Google Scholar] [CrossRef] [PubMed]
  3. Hambardzumyan, D.; Bergers, G. Glioblastoma: Defining Tumor Niches. Trends Cancer 2015, 1, 252. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Patel, A.P.; Tirosh, I.; Trombetta, J.J.; Shalek, A.K.; Gillespie, S.M.; Wakimoto, H.; Cahill, D.P.; Nahed, B.V.; Curry, W.T.; Martuza, R.L.; et al. Single-cell RNA-seq highlights intratumoral heterogeneity in primary glioblastoma. Science 2014, 344, 1396–1401. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Zhu, Y.; Parada, L.F. The molecular and genetic basis of neurological tumours. Nat. Rev. Cancer 2002, 2, 616–626. [Google Scholar] [CrossRef]
  6. Yap, T.A.; Gerlinger, M.; Futreal, P.A.; Pusztai, L.; Swanton, C. Intratumor Heterogeneity: Seeing the Wood for the Trees. Sci. Transl. Med. 2012, 4, 127ps10. [Google Scholar] [CrossRef] [Green Version]
  7. Meyer, M.; Reimand, J.; Lan, X.; Head, R.; Zhu, X.; Kushida, M.; Bayani, J.; Pressey, J.C.; Lionel, A.C.; Clarke, I.D.; et al. Single cell-derived clonal analysis of human glioblastoma links functional and genomic heterogeneity. Proc. Natl. Acad. Sci. USA 2015, 112, 851–856. [Google Scholar] [CrossRef] [Green Version]
  8. Shlush, L.I.; Hershkovitz, D. Clonal evolution models of tumor heterogeneity. Am. Soc. Clin. Oncol. Educ. Book. Am. Soc. Clin. Oncol. Annu. Meet. 2015, 35, e662–e665. [Google Scholar] [CrossRef]
  9. Albertson, D.G.; Collins, C.; McCormick, F.; Gray, J.W. Chromosome aberrations in solid tumors. Nat. Genet. 2003, 34, 369–376. [Google Scholar] [CrossRef]
  10. Shlien, A.; Malkin, D. Copy number variations and cancer. Genome Med. 2009, 1, 1–9. [Google Scholar] [CrossRef] [Green Version]
  11. Stratton, M.R.; Campbell, P.J.; Futreal, P.A. The cancer genome. Nature 2009, 458, 719–724. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. McLendon, R.; Friedman, A.; Bigner, D.; Van Meir, E.G.; Brat, D.J.; Mastrogianakis, G.M.; Olson, J.J.; Mikkelsen, T.; Lehman, N.; Aldape, K.; et al. Comprehensive genomic characterization defines human glioblastoma genes and core pathways. Nature 2008, 455, 1061–1068. [Google Scholar] [CrossRef]
  13. Assem, M.; Sibenaller, Z.; Agarwal, S.; Al-Keilani, M.S.; Alqudah, M.A.Y.; Ryken, T.C. Enhancing Diagnosis, Prognosis, and Therapeutic Outcome Prediction of Gliomas Using Genomics. Omics A J. Integr. Biol. 2012, 16, 113–122. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Stichel, D.; Ebrahimi, A.; Reuss, D.; Schrimpf, D.; Ono, T.; Shirahata, M.; Reifenberger, G.; Weller, M.; Hänggi, D.; Wick, W.; et al. Distribution of EGFR amplification, combined chromosome 7 gain and chromosome 10 loss, and TERT promoter mutation in brain tumors and their potential for the reclassification of IDHwt astrocytoma to glioblastoma. Acta Neuropathol. 2018, 136, 793–803. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Taylor, A.M.; Shih, J.; Ha, G.; Gao, G.F.; Zhang, X.; Berger, A.C.; Schumacher, S.E.; Wang, C.; Hu, H.; Liu, J.; et al. Genomic and Functional Approaches to Understanding Cancer Aneuploidy. Cancer Cell 2018, 33, 676. [Google Scholar] [CrossRef] [Green Version]
  16. Laks, E.; McPherson, A.; Zahn, H.; Lai, D.; Steif, A.; Brimhall, J.; Biele, J.; Wang, B.; Masud, T.; Ting, J.; et al. Clonal Decomposition and DNA Replication States Defined by Scaled Single-Cell Genome Sequencing. Cell 2019, 179, 1207. [Google Scholar] [CrossRef]
  17. Adalsteinsson, V.A.; Ha, G.; Freeman, S.S.; Choudhury, A.D.; Stover, D.G.; Parsons, H.A.; Gydush, G.; Reed, S.C.; Rotem, D.; Rhoades, J.; et al. Scalable whole-exome sequencing of cell-free DNA reveals high concordance with metastatic tumors. Nat. Commun. 2017, 8, 1–13. [Google Scholar] [CrossRef] [Green Version]
  18. Franch-Expósito, S.; Bassaganyas, L.; Vila-Casadesús, M.; Hernández-Illán, E.; Esteban-Fabró, R.; Díaz-Gay, M.; Lozano, J.J.; Castells, A.; Llovet, J.M.; Castellví-Bel, S.; et al. CNApp, a tool for the quantification of copy number alterations and integrative analysis revealing clinical implications. Elife 2020, 9, e50267. [Google Scholar] [CrossRef]
  19. Fan, F.; Zhang, H.; Dai, Z.; Zhang, Y.; Xia, Z.; Cao, H.; Yang, K.; Hu, S.; Guo, Y.; Ding, F.; et al. A comprehensive prognostic signature for glioblastoma patients based on transcriptomics and single cell sequencing. Cell. Oncol. (Dordr) 2021, 44, 917–935. [Google Scholar] [CrossRef]
  20. Bailey, P.; Cushing, H. A Classification of the Tumors of the Glioma Group on a Histogenetic Basis with a Correlated Study of Prognosis. J. Am. Med. Assoc. 1926, 87, 268. [Google Scholar] [CrossRef]
  21. Burrell, R.A.; McGranahan, N.; Bartek, J.; Swanton, C. The causes and consequences of genetic heterogeneity in cancer evolution. Nature 2013, 501, 338–345. [Google Scholar] [CrossRef] [PubMed]
  22. Castro, L.N.G.; Tirosh, I.; Suvà, M.L. Decoding Cancer Biology One Cell at a Time. Cancer Discov. 2021, 11, 960. [Google Scholar] [CrossRef] [PubMed]
  23. Tirosh, I.; Izar, B.; Prakadan, S.M.; Wadsworth, M.H.; Treacy, D.; Trombetta, J.J.; Rotem, A.; Rodman, C.; Lian, C.; Murphy, G.; et al. Dissecting the multicellular ecosystem of metastatic melanoma by single-cell RNA-seq. Science 2016, 352, 189–196. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. 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, 1–19. [Google Scholar] [CrossRef] [PubMed]
  25. Neftel, C.; Laffy, J.; Filbin, M.G.; Hara, T.; Shore, M.E.; Rahme, G.J.; Richman, A.R.; Silverbush, D.; Shaw, M.L.; Hebert, C.M.; et al. An integrative model of cellular states, plasticity and genetics for glioblastoma. Cell 2019, 178, 835. [Google Scholar] [CrossRef] [PubMed]
  26. Reemst, K.; Noctor, S.C.; Lucassen, P.J.; Hol, E.M. The Indispensable Roles of Microglia and Astrocytes during Brain Development. Front. Hum. Neurosci. 2016, 10, 566. [Google Scholar] [CrossRef] [Green Version]
  27. Singh, S.K.; Clarke, I.D.; Terasaki, M.; Bonn, V.E.; Hawkins, C.; Squire, J.; Dirks, P.B. Identification of a Cancer Stem Cell in Human Brain Tumors. Cancer Res. 2003, 63, 5821–5828. [Google Scholar]
  28. Bao, S.; Wu, Q.; McLendon, R.E.; Hao, Y.; Shi, Q.; Hjelmeland, A.B.; Dewhirst, M.W.; Bigner, D.D.; Rich, J.N. Glioma stem cells promote radioresistance by preferential activation of the DNA damage response. Nature 2006, 444, 756–760. [Google Scholar] [CrossRef]
  29. Lee, J.H.; Lee, J.E.; Kahng, J.Y.; Kim, S.H.; Park, J.S.; Yoon, S.J.; Um, J.Y.; Kim, W.K.; Lee, J.K.; Park, J.; et al. Human glioblastoma arises from subventricular zone cells with low-level driver mutations. Nature 2018, 560, 243–247. [Google Scholar] [CrossRef]
  30. Matias, D.; Balça-Silva, J.; da Graça, G.C.; Wanjiru, C.M.; Macharia, L.W.; Nascimento, C.P.; Roque, N.R.; Coelho-Aguiar, J.M.; Pereira, C.M.; Dos Santos, M.F.; et al. Microglia/astrocytes–glioblastoma crosstalk: Crucial molecular mechanisms and microenvironmental factors. Front. Cell. Neurosci. 2018, 12, 235. [Google Scholar] [CrossRef] [Green Version]
  31. Gimple, R.C.; Bhargava, S.; Dixit, D.; Rich, J.N. Glioblastoma stem cells: Lessons from the tumor hierarchy in a lethal cancer. Genes Dev. 2019, 33, 591–609. [Google Scholar] [CrossRef] [PubMed]
  32. Keane, L.; Cheray, M.; Blomgren, K.; Joseph, B. Multifaceted microglia-key players in primary brain tumour heterogeneity. Nat. Rev. Neurol. 2021, 17, 243–259. [Google Scholar] [CrossRef] [PubMed]
  33. Graeber, M.B.; Scheithauer, B.W.; Kreutzberg, G.W. Microglia in brain tumors. Glia 2002, 40, 252–259. [Google Scholar] [CrossRef] [PubMed]
  34. Leite, D.M.; Zvar Baskovic, B.; Civita, P.; Neto, C.; Gumbleton, M.; Pilkington, G.J. A human co-culture cell model incorporating microglia supports glioblastoma growth and migration, and confers resistance to cytotoxics. FASEB J. 2020, 34, 1710–1727. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Wolf, S.A.; Boddeke, H.W.G.M.; Kettenmann, H. Microglia in Physiology and Disease. Annu. Rev. Physiol. 2017, 79, 619–643. [Google Scholar] [CrossRef] [PubMed]
  36. Bowman, R.L.; Klemm, F.; Akkari, L.; Pyonteck, S.M.; Sevenich, L.; Quail, D.F.; Dhara, S.; Simpson, K.; Gardner, E.E.; Iacobuzio-Donahue, C.A.; et al. Macrophage Ontogeny Underlies Differences in Tumor-Specific Education in Brain Malignancies. Cell Rep. 2016, 17, 2445–2459. [Google Scholar] [CrossRef] [Green Version]
  37. Maas, S.L.N.; Abels, E.R.; Van De Haar, L.L.; Zhang, X.; Morsett, L.; Sil, S.; Guedes, J.; Sen, P.; Prabhakar, S.; Hickman, S.E.; et al. Glioblastoma hijacks microglial gene expression to support tumor growth. J. Neuroinflamm. 2020, 17, 1–18. [Google Scholar] [CrossRef] [PubMed]
  38. Daneman, R.; Prat, A. The blood–brain barrier. Cold Spring Harb. Perspect. Biol. 2015, 7, a020412. [Google Scholar] [CrossRef] [Green Version]
  39. Xie, Y.; He, L.; Lugano, R.; Zhang, Y.; Cao, H.; He, Q.; Chao, M.; Liu, B.; Cao, Q.; Wang, J.; et al. Key molecular alterations in endothelial cells in human glioblastoma uncovered through single-cell RNA sequencing. JCI Insight 2021, 6, e150861. [Google Scholar] [CrossRef]
  40. Chouleur, T.; Tremblay, M.L.; Bikfalvi, A. Mechanisms of invasion in glioblastoma. Curr. Opin. Oncol. 2020, 32, 631–639. [Google Scholar] [CrossRef]
  41. Charalambous, C.; Hofman, F.M.; Chen, T.C. Functional and phenotypic differences between glioblastoma multiforme-derived and normal human brain endothelial cells. J. Neurosurg. 2005, 102, 699–705. [Google Scholar] [CrossRef] [PubMed]
  42. Miebach, S.; Grau, S.; Hummel, V.; Rieckmann, P.; Tonn, J.C.; Goldbrunner, R.H. Isolation and culture of microvascular endothelial cells from gliomas of different WHO grades. J. Neurooncol. 2006, 76, 39–48. [Google Scholar] [CrossRef] [PubMed]
  43. Luissint, A.C.; Artus, C.; Glacial, F.; Ganeshamoorthy, K.; Couraud, P.O. Tight junctions at the blood brain barrier: Physiological architecture and disease-associated dysregulation. Fluids Barriers CNS 2012, 9, 1–12. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Wang, R.; Chadalavada, K.; Wilshire, J.; Kowalik, U.; Hovinga, K.E.; Geber, A.; Fligelman, B.; Leversha, M.; Brennan, C.; Tabar, V. Glioblastoma stem-like cells give rise to tumour endothelium. Nature 2010, 468, 829–835. [Google Scholar] [CrossRef] [PubMed]
  45. Ricci-Vitiani, L.; Pallini, R.; Biffoni, M.; Todaro, M.; Invernici, G.; Cenci, T.; Maira, G.; Parati, E.A.; Stassi, G.; Larocca, L.M.; et al. Tumour vascularization via endothelial differentiation of glioblastoma stem-like cells. Nature 2010, 468, 824–830. [Google Scholar] [CrossRef] [PubMed]
  46. Reduzzi, C.; Vismara, M.; Gerratana, L.; Silvestri, M.; De Braud, F.; Raspagliesi, F.; Verzoni, E.; Di Cosimo, S.; Locati, L.D.; Cristofanilli, M.; et al. The curious phenomenon of dual-positive circulating cells: Longtime overlooked tumor cells. Semin. Cancer Biol. 2020, 60, 344–350. [Google Scholar] [CrossRef]
  47. Fais, S. Cannibalism: A way to feed on metastatic tumors. Cancer Lett. 2007, 258, 155–164. [Google Scholar] [CrossRef]
  48. Coopman, P.J.; Do, M.T.; Thompson, E.W.; Mueller, S.C. Phagocytosis of cross-linked gelatin matrix by human breast carcinoma cells correlates with their invasive capacity. Clin. Cancer Res. 1998, 4, 507–515. [Google Scholar]
  49. Persson, A.; Englund, E. The glioma cell edge—Winning by engulfing the enemy? Med. Hypotheses 2009, 73, 336–337. [Google Scholar] [CrossRef]
  50. Chang, G.H.F.; Barbaro, N.M.; Pieper, R.O. Phosphatidylserine-dependent phagocytosis of apoptotic glioma cells by normal human microglia, astrocytes, and glioma cells. Neuro Oncol. 2000, 2, 174–183. [Google Scholar] [CrossRef]
  51. Huysentruyt, L.C.; Seyfried, T.N. Perspectives on the mesenchymal origin of metastatic cancer. Cancer Metastasis Rev. 2010, 29, 695–707. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Wilhelmsson, U.; Eliasson, C.; Bjerkvig, R.; Pekny, M. Loss of GFAP expression in high-grade astrocytomas does not contribute to tumor development or progression. Oncogene 2003, 22, 3407–3411. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Hieronymus, H.; Murali, R.; Tin, A.; Yadav, K.; Abida, W.; Moller, H.; Berney, D.; Scher, H.; Carver, B.; Scardino, P.; et al. Tumor copy number alteration burden is a pan-cancer prognostic factor associated with recurrence and death. Elife 2018, 7, e37294. [Google Scholar] [CrossRef] [PubMed]
  54. Gómez-Rueda, H.; Martínez-Ledesma, E.; Martínez-Torteya, A.; Palacios-Corona, R.; Trevino, V. Integration and comparison of different genomic data for outcome prediction in cancer. BioData Min. 2015, 8, 1–12. [Google Scholar] [CrossRef] [Green Version]
  55. Roh, W.; Chen, P.L.; Reuben, A.; Spencer, C.N.; Prieto, P.A.; Miller, J.P.; Gopalakrishnan, V.; Wang, F.; Cooper, Z.A.; Reddy, S.M.; et al. Integrated molecular analysis of tumor biopsies on sequential CTLA-4 and PD-1 blockade reveals markers of response and resistance. Sci. Transl. Med. 2017, 9, eaah3560. [Google Scholar] [CrossRef] [Green Version]
Cells 11 01127 g001 550
Figure 1. Histological images of GB01, GB02, and GB03. Experimental design starting from tumor shredding to DEPArray analysis.
Figure 1. Histological images of GB01, GB02, and GB03. Experimental design starting from tumor shredding to DEPArray analysis.
Cells 11 01127 g001
Cells 11 01127 g002 550
Figure 2. Example of DEPArray images of single cells belonging to the main GB populations, stained in yellow with GFAP (astrocytes-like), in red with IBA1 (microglia-like), in purple with CD105 (endothelial-like cells), in green with CD133 (stem-like cells) and in blue with Hoechst. BF: Brightfield.
Figure 2. Example of DEPArray images of single cells belonging to the main GB populations, stained in yellow with GFAP (astrocytes-like), in red with IBA1 (microglia-like), in purple with CD105 (endothelial-like cells), in green with CD133 (stem-like cells) and in blue with Hoechst. BF: Brightfield.
Cells 11 01127 g002
Cells 11 01127 g003 550
Figure 3. Pie charts of the percentages of the main cell type populations found in GB01, GB02, and GB03.
Figure 3. Pie charts of the percentages of the main cell type populations found in GB01, GB02, and GB03.
Cells 11 01127 g003
Cells 11 01127 g004 550
Figure 4. Summary of all the recovered cells after DEPArray analysis from GB01, GB02, and GB03. Astrocytes-like, microglia-like, endothelial-like, and stem-like cells were collected. Double staining cells and only Hoechst positive cells are shown.
Figure 4. Summary of all the recovered cells after DEPArray analysis from GB01, GB02, and GB03. Astrocytes-like, microglia-like, endothelial-like, and stem-like cells were collected. Double staining cells and only Hoechst positive cells are shown.
Cells 11 01127 g004
Cells 11 01127 g005 550
Figure 5. CNA pattern of the fresh GB tissues in bulk. The chromosomal amplifications are shown in red, and in blue the deletions. The intensity of the red and blue color components correlates with the gain and loss values based on the results obtained from the CNApp tool.
Figure 5. CNA pattern of the fresh GB tissues in bulk. The chromosomal amplifications are shown in red, and in blue the deletions. The intensity of the red and blue color components correlates with the gain and loss values based on the results obtained from the CNApp tool.
Cells 11 01127 g005
Cells 11 01127 g006 550
Figure 6. Genome-wide chromosome arm CNA profile heatmap for GB01, GB02, and GB03. For each sample, the CNA profile of the single cells collected is shown and on the right the CNA profile of the tumor tissue in bulk.
Figure 6. Genome-wide chromosome arm CNA profile heatmap for GB01, GB02, and GB03. For each sample, the CNA profile of the single cells collected is shown and on the right the CNA profile of the tumor tissue in bulk.
Cells 11 01127 g006
Cells 11 01127 g007 550
Figure 7. Immunofluorescence assays of GB01 tissue: (a–c) three different portions of GB01 tissue with GFAP (red) and IBA1 (green) markers shown. Arrows indicate cells with the double-positive signal.
Figure 7. Immunofluorescence assays of GB01 tissue: (a–c) three different portions of GB01 tissue with GFAP (red) and IBA1 (green) markers shown. Arrows indicate cells with the double-positive signal.
Cells 11 01127 g007
Table
Table 1. Patient clinical, demographic, pathological, and therapeutical data.
Table 1. Patient clinical, demographic, pathological, and therapeutical data.
CasesAgeSexPrimary or RecurrenceBrain LocationIDH1/IDH2Pathology ReportTherapy Administered
GB0130MPrimaryparietal lobeWTGlioblastoma (Grade IV WHO) (GFAP+, MKI67-20%)Levetiracetam, Soldesam, Lansoprazole
GB0247MPrimaryright temporal lobeWTGlioblastoma (Grade IV WHO) (GFAP+, MKI67-30%)Levetiracetam, Dexamethasone, Omeprazole
GB0365MPrimaryright frontal lobeWTGlioblastoma (Grade IV WHO) (GFAP+, MKI67-20%)Levetiracetam, Lansoprazole, Dexamethasone (Mannitol pre-op)
Table
Table 2. CNAs results obtained after CNApp processing for single cells and groups of single cells collected in the GB01 sample.
Table 2. CNAs results obtained after CNApp processing for single cells and groups of single cells collected in the GB01 sample.
GB01
Single Cells CollectedCNA
Group of endothelial-like cellsWT
Microglia-like cell19-
Microglia-like cellWT
Microglia-like cell1p-, 7+, 9q+, 10-, 17q+, 19+
Microglia-like cell1q+, 2-, 5+, 7+, 9p-, 9q+, 10-, 17q+, 19+
Group of microglia-like cellsWT
Group of microglia-like cells 7+, 9q+, 10-, 17q+
Astrocyte-like1p-, 7+, 9q+, 10-, 17q+, 19+
Astrocyte-like1p-, 7+, 9q+, 10-, 17q+, 19+
Astrocyte-like1q+, 5+, 7+, 9p-, 9q+, 10-, 17p-, 17q+, 19q+
Group of astrocytes-like1p-, 7+, 9q+, 10-, 17p-, 17q+, 19+
Group of astrocytes-like1p-, 7+, 9+, 10-, 17p-, 17q+, 19+
Group of astrocytes-like1p-, 7+, 9q+, 10-, 17q+
Stem-like cell1p-, 7+, 9+, 10-, 17q+, 19q+
Astrocyte/microglia-like cell1p-, 7+, 9+, 10-, 17q+, 19q+
Astrocyte/microglia-like cell1p-, 7+, 9+, 10-, 17q+, 19+
Not stained cellWT
Not stained cellWT
Not stained cell1p-, 7+, 9q+, 10-, 17p-, 17q+, 19q+
Group of not stained cellsWT
Table
Table 3. CNAs results obtained after CNApp processing for single cells and groups of single cells collected in the GB02 sample.
Table 3. CNAs results obtained after CNApp processing for single cells and groups of single cells collected in the GB02 sample.
GB02
Single Cells CollectedCNA
Endothelial-like cellWT
Endothelial-like cell19-
Endothelial-like cell19-
Endothelial-like cell19-
Endothelial-like cell9p-, 10-, 13q-, 14q-, 22q-
Group of endothelial-like cellsWT
Group of endothelial-like cells19-
Group of endothelial-like cells19-
Astrocyte-likeWT
Astrocyte-like19-
Astrocyte-like9p-, 10-, 13q-, 14q-, 22q-
Astrocyte-like9p-, 10-, 11- 13q-, 14q-, 19p-, 22q-
Astrocyte-like9p-, 10-, 11- 13q-, 14q-, 20+, 22q-
Astrocyte-like9p-, 10-, 13q-, 14q-, 19p-, 20+ 22q-
Microglia-like cellWT
Microglia-like cellWT
Microglia-like cell19-
Microglia-like cell19-
Microglia-like cell19-
Endothelial/stem-like cellWT
Endothelial/stem-like cellWT
Endothelial/stem-like cellWT
Endothelial/stem-like cell9p-, 10-, 13q-, 14q-, 22q-
Endothelial/stem-like cell 10-, 11-, 13q-, 14q-, 16+, 22q-
Group of endothelial/stem-like cellsWT
Not stained cell9p-, 10-, 13q-, 14q-, 22q-
Table
Table 4. CNAs results obtained after CNApp processing for single cells and groups of single cells collected in the GB03 sample.
Table 4. CNAs results obtained after CNApp processing for single cells and groups of single cells collected in the GB03 sample.
GB03
Single Cells CollectedCNA
Endothelial-like cellWT
Endothelial-like cellWT
Endothelial-like cellWT
Endothelial-like cellWT
Endothelial-like cell7+, 9p-, 9q+, 10-, 20+, 22q-
Endothelial-like cell3q-, 7+, 9p-, 9q+, 10-, 20+, 22q-
Microglia-like cellWT
Microglia-like cellWT
Microglia-like cellWT
Microglia-like cellWT
Microglia-like cellWT
Astrocyte-likeWT
Astrocyte-like7+, 9p-, 9q+, 10-, 20+, 22q-
Astrocyte-like7+, 9p-, 9q+, 10-, 20+, 22q-
Astrocyte-like7+, 9p-, 9q+, 10-, 20+, 22q-
Astrocyte-like7+, 9p-, 9q+, 10-, 20+, 22q-
Astrocyte-like7+, 9p-, 9q+, 10-, 20+, 22q-
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lessi, F.; Franceschi, S.; Morelli, M.; Menicagli, M.; Pasqualetti, F.; Santonocito, O.; Gambacciani, C.; Pieri, F.; Aquila, F.; Aretini, P.; et al. Single-Cell Molecular Characterization to Partition the Human Glioblastoma Tumor Microenvironment Genetic Background. Cells 2022, 11, 1127. https://0-doi-org.brum.beds.ac.uk/10.3390/cells11071127

AMA Style

Lessi F, Franceschi S, Morelli M, Menicagli M, Pasqualetti F, Santonocito O, Gambacciani C, Pieri F, Aquila F, Aretini P, et al. Single-Cell Molecular Characterization to Partition the Human Glioblastoma Tumor Microenvironment Genetic Background. Cells. 2022; 11(7):1127. https://0-doi-org.brum.beds.ac.uk/10.3390/cells11071127

Chicago/Turabian Style

Lessi, Francesca, Sara Franceschi, Mariangela Morelli, Michele Menicagli, Francesco Pasqualetti, Orazio Santonocito, Carlo Gambacciani, Francesco Pieri, Filippo Aquila, Paolo Aretini, and et al. 2022. "Single-Cell Molecular Characterization to Partition the Human Glioblastoma Tumor Microenvironment Genetic Background" Cells 11, no. 7: 1127. https://0-doi-org.brum.beds.ac.uk/10.3390/cells11071127

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