This article has Open Peer Review reports available.
Identification of a novel set of genes reflecting different in vivo invasive patterns of human GBM cells
© Monticone et al.; licensee BioMed Central Ltd. 2012
Received: 2 December 2011
Accepted: 15 July 2012
Published: 17 August 2012
Most patients affected by Glioblastoma multiforme (GBM, grade IV glioma) experience a recurrence of the disease because of the spreading of tumor cells beyond surgical boundaries. Unveiling mechanisms causing this process is a logic goal to impair the killing capacity of GBM cells by molecular targeting.
We noticed that our long-term GBM cultures, established from different patients, may display two categories/types of growth behavior in an orthotopic xenograft model: expansion of the tumor mass and formation of tumor branches/nodules (nodular like, NL-type) or highly diffuse single tumor cell infiltration (HD-type).
We determined by DNA microarrays the gene expression profiles of three NL-type and three HD-type long-term GBM cultures. Subsequently, individual genes with different expression levels between the two groups were identified using Significance Analysis of Microarrays (SAM). Real time RT-PCR, immunofluorescence and immunoblot analyses, were performed for a selected subgroup of regulated gene products to confirm the results obtained by the expression analysis.
Here, we report the identification of a set of 34 differentially expressed genes in the two types of GBM cultures. Twenty-three of these genes encode for proteins localized to the plasma membrane and 9 of these for proteins are involved in the process of cell adhesion.
This study suggests the participation in the diffuse infiltrative/invasive process of GBM cells within the CNS of a novel set of genes coding for membrane-associated proteins, which should be thus susceptible to an inhibition strategy by specific targeting.
Massimiliano Monticone and Antonio Daga contributed equally to this work
The Glioblastoma multiforme (GBM, stage IV Glioma) arise from neuroglial cells or their progenitors and represents the most aggressive brain tumor, with 15 months median survival after diagnosis, causing 4% of all cancer-related death despite recent improvement of diagnostic and treatment procedures. Surgery represents the standard treatment procedure. However, the vast majority of the patients affected by GBM experience a recurrence of the disease because of the spreading of cells beyond the limits of the resection . The identification of the affected region of the central nervous system (CNS) to be resected is a major challenge. Neither advanced imaging techniques nor histological examination warrant against leaving some tumor cells in adjacent normal-looking brain tissue. Histologically normal brain tissue acquired at a distance greater than 4 cm from the GBM/Oligodendroglioma tumor was shown to give rise to tumor colonies in soft agar culture . Therefore, the ability of GBM cells to invade the host tissue is one of the biological features of this disease that eventually has the most detrimental impact on the life expectancy of the patient . In addition, these cells are difficult to eradicate since they invade areas of the CNS with an intact blood–brain barrier. As a consequence, the targeting of the GBM invasion process is a major topic of interest [3–5].
In the past years, some reports have focused on the inverse relationship between growth/apoptosis sensitivity and migration of glioma cells  and production by glioma cells of factors able to enhance invasion in an autocrine fashion . Other research groups have shown the ability of the microenvironment to influence migration properties via cell-extracellular matrix interactions and paracrine stimuli [3, 4, 7–9].
Key to the study of GBM invasion is the availability of a reliable culture system, in order to preserve the tumorigenic potential of cells derived from patients, and of “in vivo” models suitable to address questions and test hypothesis concerning this process.
To these aims, we have successfully established long-term cell cultures from surgical tumor samples obtained from several GBM patients and demonstrated their ability to generate GBM xenografts by serial transplantation[10, 11]. In particular, we observed that these cultures displayed two types of “in vivo” growth behavior in these transplants. The first one was mainly expansive while the second, causing the host’s white and gray matters substitution by tumor cells, was highly diffusive. The aim of the present study was to identify by microarray analysis if the two GBM culture types were characterized by differential gene expression. We actually identified a set of differentially expressed genes. Some of these were known to be involved directly or indirectly in promoting glioma invasion, which supported our results [12–14]. Other genes, however, were not previously described in association with glioma invasion. This study provides, therefore, a novel set of potential target genes for future research and development of treatment strategies intended to inhibit the invasion by GBM cells of healthy brain tissue.
Patients, tumor characteristics and survival of mice inoculated with patient GBM cells
Patient and GBM cells code
Newly diagnosed/ recurrence
Survival (days) of mice bearing secondary tumors
Early subependymal involvement
Cells isolation was described in detail elsewhere [11, 16]. Briefly, bioptic samples were plated in regular plastic dishes using proliferation medium. One to two weeks after plating, cellular aggregates, resembling neurospheres derived from normal neural precursors, were detectable in all glioblastoma cultures. Neurospheres-like aggregates were collected, dissociated to obtain a single cell suspension and seeded on Matrigel (BD Biosciences, San Jose, CA, USA) coated flasks. Under these conditions, cells grown as a monolayer gave rise to tumor when injected orthotopically in nude mice [11, 16]. In addition, these long-term cultures of GBM tumorigenic cells maintained the ability to generate neurospheres-like aggregates repeatedly when transferred into flasks lacking Matrigel coating. In particular, the long-term cultures of GBM tumorigenic cells used in the present study did not show differences in growth pattern in vitro because they were all able to attach, spread and proliferate on Matrigel. Similarly, when seeded in culture flasks lacking the Matrigel coat, they were all able to generate neurosphere-like aggregates and proliferate with no major phenotype differences in neurosphere dimension and appearance (data not shown).
The proliferation medium was DMEM-F12/Neurobasal additioned with 1% v/v B27 supplement, (Gibco Ltd, Paisley, Scotland), 2 mM L-glutamine (Gibco Ltd), recombinant human fibroblast growth factor (FGF-2, 10 ng/ml Peprotech, London, UK), recombinant human epidermal growth factor (EGF, 20 ng/ml Peprotech). The medium was changed twice a week. Normal human astrocytes were purchased from ScienCell Research Laboratories (Carlsbad, CA) and cultured following the manufacturer’s instructions. Normalized expression levels for selected stem cell markers, nervous system markers, PDGF receptors and IDH genes in the tumorigenic long-term GBM cell cultures used in this study, are shown in Additional File 1: Table S1.
Intracranial tumorigenicity assays
Human GBM cells in vivo tumorigenicity was tested by cell intracranial inoculation in 6 – 8-week-old NOD/SCID mice (Charles River Laboratories, Lecco, Italy).
NOD/SCID mice were housed in pathogen-free conditions, according to the National Regulation on Animal Research Resources. For intracranial inoculation, at least three mice for each GBM cell culture were used. Mice were anesthetized with i.m. ketamine and xylazine. Thereafter, the animal was positioned into a stereotactic frame (David Kopf Instruments, Tujunga, CA, USA) and a hole was made 2 mm lateral and 1 mm anterior from the bregma. Cells (105) were injected using a Hamilton syringe (Sigma-Aldrich, Milan, Italy) at a depth of 3.5 mm in a vol of 2 μl. Mice were monitored for about 6 months for disease symptoms and were sacrificed by CO2 asphyxiation when they showed weight loss or any severe sign of disease. The brains of all sacrificed mice were collected and processed for histology. The survival in days of mice inoculated with the long-term GBM cultures used in this study is reported in Table1.
All experiments were performed in compliance with guidelines approved by the Ethical Committee for Animal Use in Cancer Research at the Istituto Nazionale per la Ricerca sul Cancro in Genoa, Italy. Under our experimental conditions, the minimum number of GBM cells required to give rise to tumor in mice was 104.
RNA extraction and quality control analysis
Total RNA was isolated from cultured GBM cells using miRNeasy® mini kit (Qiagen, Milan, Italy) with DNase treatment. RNA concentration and purity were determined by measuring absorbance at 260 and 280 nm; 2 μg total RNA was run on a 1% denaturing gel and 100 ng were loaded on the 2100 Bioanalyzer (Agilent, Palo Alto, CA) to verify RNA integrity.
Amplification of RNA and array hybridization
According to the recommendations of the manufacturer, 100 ng of total RNA was used in the first-round synthesis of double-stranded cDNA. The RNA was reverse-transcribed using a Whole Transcript cDNA synthesis and amplification kit (Affymetrix UK Ltd., High Wycombe, UK). The resulting biotin-labeled cRNA was purified using an IVT clean-up kit (Affymetrix) and quantified using a UV spectrophotometer (A260/280; Beckman, Palo Alto, CA). An aliquot (15 μg) of cRNA was fragmented by heat and ion-mediated hydrolysis at 94°C for 35 minutes. The fragmented cRNA, run on the Bioanalyzer (Agilent Technologies, Santa Clara, CA) to verify the correct elettropherogram, was hybridized in a hybridization oven (16 hours, 45 °C) to a Human Gene 1.0 ST array (Affymetrix) representing whole-transcript coverage. Each one of the 28869 genes was represented on the array by approximately 26 probes spread across the full length of the gene, providing a more complete and more accurate picture of gene expression than the 3’ based expression array design. The washing and staining procedures of the arrays with phycoerythrin-conjugated streptavidin (Invitrogen, Monza, Italy) was completed in the Fluidics Station 450 (Affymetrix). The arrays were subsequently scanned using a confocal laser GeneChip Scanner 3000 7 G and the GeneChip Command Console (Affymetrix).
GeneChip microarray analysis and data normalization
Affymetrix raw data files [cell intensity (CEL) files] were used as input files in expression console environment (Affymetrix). Briefly, CEL files were processed using the Robust Multi-Array Analysis (RMA) procedure , an algorithm that is publicly available at http://www.bioconductor.org. The RMA method was used to convert the intensities from the multiple probes of a probe set into a single expression value with greater precision and reduced background noise (relying on the perfect match probes only and thus ignoring the mismatch probes) and then to normalize by sketch quantile normalization. Quality assessments were also performed in the expression console environment. This procedure, based on various metrics, allowed us to identify a chip as an outlier (see for details Quality assessment of exon and gene arrays http://www.affymetrix.com/support/technical/whitepapers/exon_gene_arrays_qa_whitepaper.pdf). Significance Analysis of Microarrays (SAM), Principal Component Analysis (PCA) of variance and Hierarchical Clustering (HCL), after mean scaling and log2 transformation, were performed with the software tool from The Institute for Genomic Research (TIGR) MeV (multiple experimental viewer) (http://www.tigr.org/software/tm4/mev.html) .
Individual genes with different expression levels, among the two groups, were identified using SAM . The false discovery rate expressed as q-value was used to evaluate statistical significance. For comparison purposes, an arbitrary filter was applied excluding all genes that did not exhibit a difference in expression of at least 2-fold. Genes differentially expressed were investigated using a two-class analysis.
We used PCA to reduce the complexity of high-dimensional data and to simplify the task of identifying patterns and sources of variability in these large data sets. The results from SAM were visualized using HCL . All the microarray information has been submitted to the National Center for Biotechnology Information Gene Expression Omnibus web site (http://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?token=zfutlaeiauqmmbo&acc=GSE16805).
Pathways identification by Expression Analysis Systemic Explore (EASE)
Gene lists from Affymetrix results were examined using the EASE program, accessible via http://david.abcc.ncifcrf.gov/. EASE is a customized stand-alone software application with statistical functions for discovering biological themes within gene lists. This software assigns genes of interest into functional categories based on the Gene Ontology database (GO, http://www.geneontology.org/index.shtml) and uses the Fisher's exact test statistics to determine the probability of observing the number of genes within a list of interest versus the number of genes in each category on the array. A more detailed analysis of the genes' association with physiological pathways was performed using the Kyoto Encyclopedia of Genes and Genomes (KEGG, http://www.genome.jp/kegg/pathway.html). Each identified process was confirmed through PubMed/Medline (http://www.ncbi.nlm.nih.gov/sites/entrez?db=pubmed).
Starting from about 1 μg of total RNA, cDNA was synthesized by using an Oligo(dT)20, random hexamers mix and a Superscript III first-strand synthesis system supermix for RT-PCR (Invitrogen). cDNAs, diluted 5–20 times, were then subjected to PCR analysis.
Relative quantification was performed by real-time quantitative RT-PCR (qPCR). Briefly, qPCR was performed and analyzed using the Mastercycler ep Realplex (Eppendorf AG, Hamburg, Germany). Primers were designed across a common exon–exon splice junction by using the tool available at https://www.roche-applied-science.com/sis/rtpcr/upl/index.jsp?id=UP030000 (Roche Applied Science, Monza, Italy)(see Additional File 2: Table S2). Reactions were carried out in triplicates and amplicons were measured by SYBR Green fluorescence (5 Prime, Hamburg, Germany) according to the manufacturer's recommendations. The dissociation curve analysis was used to define the specificity of the products by the presence of a single dissociation peak on the thermal melting curve.
The gene coding for the peptidyl-prolyl cis-trans isomerase A (PPIA) was used as the endogenous control for normalization because, in the microarray data, it showed in all conditions the steadiest expression in our experimental setting when compared with other housekeeping genes.
Multiplex ligation-dependent probe amplification (MLPA) analysis
Genomic DNA was isolated from long-term GBM cell cultures and normal human astrocytes using QIAamp DNA microkit (Qiagen). The MLPA analysis (SALSA MLPA KITs P175-A1 Tumour-Gain and P294-A1 Tumour-loss, (MRC Holland, Amsterdam, the Netherlands) was performed using 100 ng of genomic DNA, diluted in 5 μl of TE buffer, following the manufacturer’s instructions. The resulting DNA fragments were identified and quantified by using capillary electrophoresis on an ABI XL3130 genetic analyzer (Applied Biosystems, Foster City, CA, USA) and the Genemapper program (version 4.0 - Applied Biosystems). Data were analyzed with the Coffalyser software (MRC-Holland). For each patient long-term culture of tumorigenic GBM cells, gains and losses were assigned by comparing the peaks between the patient and the reference samples (DNA from normal human astrocytes). A value =0 was considered as a biallelic loss, a value <0.7 corresponded to a DNA loss, a value between 1.3 and 2.0 corresponded to a DNA gain. The examined region was defined as amplified for values >2. All experiments were done in triplicate.
Fluorescence in situ hybridization (FISH)
FISH was performed using a whole chromosome painting for chromosome 10 (WCP-10), obtained from Mariano Rocchi, Resources for Human Molecular Cytogenetics Project, by Telethon and the Italian Association for Cancer Research (AIRC). WCP-10 was amplified and labeled with Cy-dUTP (Amersham Pharmacia Biotech, Piscataway, NJ, USA) by degenerate oligonucleotide-primed polymerase chain reaction (DOP-PCR) . Hybridization was performed as described , with minor modification. Briefly, slides carrying metaphase spreads were denatured in 70% (vol/vol) formamide/2x SSC, pH 7, at 70°C for 2 min], and dehydrated in a 4°C ethanol series. The hybridization mix consisted of 50% formamide, 2x SSC, 10% dextran sulfate, carrier DNA (sonicated salmon sperm DNA) at 500 μg/ml, human DNA C0t1 and Cy3-labeled WCP-10 at 2 μg/ml. This mixture was denaturated at 76°C for 6 min, incubated at 37°C for 30 min, applied to the slides under a glass coverslip, and sealed with rubber cement. After overnight incubation at 37°C in Hybrite™ (Abbott Molecular/Vysis, Abbott Park, IL, USA), the slides were washed for 2 min at 72°C in 0.4x SSC/0.3% NP40, pH 7 and for 1 min at R.T. in 2x SSC/0,1% NP40, pH7. Slides were counterstained with 4’,6- Diamidino-2-phenylindole, DAPI (Sigma-Aldrich) and mounted in Mowiol. Hybridization signals were evaluated using a digital image analysis system based on an epifluorescence (Provis AX70, Olympus, Milan, Italy), and images were acquired with a digital CCD camera C4742 Orca II (Hamamatsu, Japan) driven by Cytovision (Applied Imaging Corp., Santa Clara, CA, USA). The DAPI and Cy3 images were acquired with selective single-bandpass filters at 1000× optical magnifications.
Antibodies for Sema5A (ab51957) were purchased from Abcam (Abcam, Cambridge, UK). The Nestin antibody (MAB353) was from Chemicon-Millipore (Millipore SPA, Milan, Italy). The BCAN antibody (HPA007865) was from Sigma-Aldrich. The Pmp2 antibody (12717-AP) was purchased from Protein Tech (Protein Tech Group, Inc., Chicago, IL, USA).
For xenograft tumor analysis, brains were cryopreserved and 10-μm cryostat (Leica Microsystems, Milan, Italy) sections were cut. Sections bearing tumors were identified by H&E. Cryosections containing tumors were permeabilized in PBS containing 0.2% Triton X-100 and blocked in 5% normal FCS-PBS. After incubation with primary antibodies (anti-Nestin), sections were stained with the appropriate secondary antibodies.
Cells plated on Matrigel-coated glass coverslips were fixed in 3% paraformaldehyde in PBS, pH 7.6 containing 2% sucrose for 5 min at room temperature. After PBS rinse, cells were permeabilized by a solution containing 20 mM Hepes pH 7.4, 300 mM sucrose, 50 mM NaCl, 3 mM MgCl2 and 0.5% Triton X-100 for 5 min on ice. Non-specific binding was prevented by incubation with pure goat serum for 30 min on ice. Slides were incubated with primary antibody in PBS supplemented with 10% goat serum or fetal bovine serum (antibody dilution buffer) for 2 hours on ice. After extensive PBS washes, slides were incubated for 30 min with the appropriate Alexa594 or Alexa488-conjugated secondary antibody in antibody dilution buffer. Nuclei were stained with Hoechst 33258 (Invitrogen).
Image acquisition were performed at 23°C with an Axiovert 200 M microscope equipped with the following filter sets: Zeiss 49, Zeiss 10 and Omega XF102-2 (Omega Optical, Brattleboro, VT, USA), which were used to detect Hoechst, Alexa 488 and Alexa 594 respectively.
SDS PAGE and immunoblot analysis
Cell lysates were obtained as described previously . SDS PAGE and immunoblotting were performed using precasted 4-12% polyacrylamide gels (Invitrogen) using manufacturer’s instructions. The blotting membranes were from Millipore (Millipore S.p.A., Vimodrone, Italy). The Pmp2 antibody (12717-AP) was purchased from Protein Tech (Protein Tech Group, Inc., Chicago, IL, USA). The antibodies for Bcan, Megf10, Pcdh10, Pcdh15 and Actin were all from Sigma-Aldrich (HPA007865, HPA026876, HPA011220, AV50153 and A 2066, respectively). The antibody for CD109 (Sc-271085) was purchased from Santa Cruz Biotechnology (Santa Cruz Biotechnology, Santa Cruz, CA, USA). The Gria2 antibody was acquired from Abcam ab40878.
Human GBM cells in mouse orthotopic transplantation show two distinct patterns of tumor growth and tissue invasion
Cultures of human tumorigenic GBM cells display unique chromosomal copy number aberrations unrelated to distinct in vivo invasion patterns
To investigate whether the two different in vivo invasion patterns of tumorigenic GBM cells used in the present study could be related to differences in gains or losses at genomic loci corresponding to known proto-oncogenes or tumor suppressor genes, we performed a Multiplex ligation-dependent probe amplification (MLPA) analysis. We investigated a total of 51 loci and found a number of alterations with respect to genomic DNA from primary cultures of human astrocytes, used a normal control reference. We observed that the cell cultures from each patient displayed a unique profile of losses and gains with the only exception for the CDKN2A biallelic loss, which was common to all six cultures (see Additional File 3: Table S3). However, we also observed that PT1, PT2 and PT3 (NL-type) cells shared a combined DNA gain at the EGFR locus and a DNA loss at the PTEN locus, which was not observed in the HD-type of GBM cells (see Additional File 3: Table S3). In fact, among HD-type of GBM cells we found either six copies of the EGFR locus or biallelic loss at the PTEN locus or no change at these two loci (PT4, PT6 and PT5, respectively, Additional File 3: Table S3). To investigate whether the two types of tumorigenic GBM cell cultures differed for chromosome 10 copy number, we performed a Fluorescence in situ hybridization (FISH) analysis using a whole chromosome 10-specific painting probe. The result of this analysis showed that chromosome 10 copy number cannot differentiate the two types of tumorigenic GBM cell cultures (Additional File 4: Table S4).
Cultures of human tumorigenic GBM cells displaying distinct in vivo invasion patterns display a set of differentially expressed genes
We hypothesized that the different in vivo invasion pattern observed could be related to gene differentially expressed in the two groups of long-term tumorigenic GBM cultures. In order to identify these genes, we compared the gene expression patterns, obtained by microarray analysis performed on the Affymetrix platform. In particular, we used the TIGR MEV program and mRNAs extracted from GBM cultures previously isolated from patients 1–3 (PT1-3), belonging to the NL-type and from GBM cultures previously isolated from patients 4–6 (PT 4–6) belonging to the HD-type.
It should be noticed that the principal component analysis showed a scattering of the samples belonging to the NL-type (PT1-3). Based on gene expression, a molecular classification of glioblastoma into four different profiles called ProNeural, Neural, Classical and Mesenchymal was recently proposed . For each of the 4 profiles, the corresponding centroid on the basis of the expression of 840 genes was identified. We extracted the expression profiles of these 840 genes from our dataset to classify our samples according to the correlation with each centroid. In this way we classified PT1, PT2 and PT3 as Neural, PT4 and PT6 as Mesenchymal (data not shown). Unfortunately, we were unable to similarly classify the profile of PT5. It is likely that PT5 was not similar to those previously published. Therefore, because PT1-3 were all Neural, their scattering in the principal component analysis appeared to result from the expression of genes different from those 840 used for the classification.
Probe sets ID, with relative mRNA Accession or Ensembl Transcript ID, gene symbols of genes regulated in the comparison between cultures of human GBM tumorigenic cells with the highly diffusive (HD-type) and the nodular-like (NL-type) invasive patterns, as determined by using SAM software with two-class unpaired analysis and the additional requirement of at least a 2-fold change in gene expression (Ratio)
Probe Set ID
mRNA Accession / ENSEMBL transcript ID
Ratio HD / NL expression
To gain a more mechanistic understanding of the processes associated to the HD phenotype, the EASE score  was used to identify Gene Ontology (GO) functional categories, which were significantly over-represented. After filtering the results, to avoid redundant and/or generic categories, two statistically significant GO term Biological Process and Cell Compartment were found associated with the HD-type of invasive phenotype: cell adhesion and intrinsic to membrane having a P-value of 1.1×10-4 and 1.9×10-5 respectively. In particular, 23 out of 34 genes were coding for intrinsic membrane protein and among these 23 genes, 9 were associated to cell adhesion processes (Table2).
Validation at the protein level of the differential gene expression by NL-type and HD-type of human GBM tumorigenic cultures
Targeting the GBM cell ability to invade the surrounding healthy CNS tissue is a goal to be obtained at the earliest stage of the disease in order to reduce the occurrence of relapses after surgery.
In order to contribute to the identification of novel target genes having a potential role in this process, we undertook a screening study by gene expression analysis.
We first observed that cells from our collection of GBM cultures, when transplanted in mice brains, were able to grow either as a nodule-like mass with a tendency to remain confined in the subcortical regions or to pervasively infiltrate as single cells the entire host’s CNS and eventually substitute the host’s tissue with the GBM tissue.
With this observation in mind, we hypothesized that specific genes able to drive or to inhibit the diffuse infiltrative behavior could be found up- or down-regulated respectively in GBM cultures displaying this in vivo phenotype.
Our gene expression analysis, by using principal component, first demonstrated that GBM cultures belonging to the diffuse infiltrating type clustered together separately from those belonging to the less infiltrating class. Secondly, cell adhesion appeared predominant among the functional processes associated to the 34 differentially expressed genes. Moreover, most genes belonging to this regulated set (23) appeared to encode for proteins intrinsic to the plasma membrane.
Real time PCR, immunofluorescence and immunoblot analyses, performed by using probes against eight arbitrarily selected gene products out of thirty-four, confirmed the differential regulation determined by microarray gene expression analysis.
Interestingly, among the genes that we found to be up-regulated in the HD-type, a few like GRIA2, BCAN and LPAR4 were already shown to be directly or indirectly implicated with glioma cell invasion [12–14, 26].
It is worth of note that because we did not observe a differential regulation of genes downstream of the EGFR, PI3K, P53 and AKT pathways, the combination of gains at the EGFR locus and losses at the PTEN locus, found in the NL-type and not present in the HD-type of cultures, was likely not responsible for the different in vivo invasion behavior. One of the possible explanations, for this lack of differential regulation of genes downstream of the mentioned pathways, is the finding of amplification at the EGFR locus (six copies) in PT4 and the biallelic loss at the PTEN locus in PT6 (both HD-type cells), which may result in stimulation of these pathways similar to the combination of monoallelic EGFR gain and PTEN loss.
Our data also showed that, although the long-term tumorigenic GBM cultures used in this study displayed loss of heterozygosis (LOH) at specific loci on chromosome 10q (i.e., PTEN and RET), they contained 2 copies (PT2, PT3, PT5, PT6) or 3 copies (PT1, PT4) of chromosome 10, independently from the in vivo invasion behavior. Therefore, monosomy for chromosome 10, which is frequent but not always observed in GBMs , could not be responsible for the observed different invasive property of these cells.
The difference in invasive behavior among HD-type and NL-typecould reflect in part the fact that two of the HD-type cultures were Mesenchymal (PT4 and PT6) whereas the NL-types were all Neural according to a GBM signature previously reported . Interestingly, Tchoghandjian and collaborators, observed similar infiltrative and multifocal clusters in vivo invasion patterns displayed by Mesenchymal and Neural GBM-derived stem-like cells, respectively . The result that BCAN was found expressed at the highest level in PT5 and PT6, which we classified as Mesenchymal, is in apparent disagreement with previous studies [23, 29]. However, several reports showed that BCAN is associated to invasive glioma and promotes glioma invasion after proteolytic cleavage and fibronectin binding [12, 14, 30–34], which support our data. Moreover, the signature published by Verhaak and collaborators , which included BCAN and classified GBM cells, was based on the global expression of 840 genes. Therefore, in our opinion, some variation in the expression of a single gene may still occur without affecting GBM cell classification.
Previous studies have linked GBM cell lines in vitro and in vivo transplant growth patterns as well as gene signatures with cortical and deep tumor location in patient brains [28, 29, 35, 36]. The analysis of the patient tumor characteristics, from which the long-term GBM cultures were established in our study, showed overall similarity in tumor location and lobe oforigin between the two invasive behavior types. A possible exception might be the frontal lobe involvement, which was present in two out of three patientsfrom which were derived the HD-type GBM cells and absent from those patients from which were established the NL-type of GBM cells. In our opinion, our series of cases is small to confirm or disprove a significant correlation between in vitro and in vivo growth patterns and invasive phenotype and tumor or patient characteristics, which was, however, not the aim of this study.
In conclusion, we think that the present study has identified for the first time a set of genes that are likely to be implicated in the diffuse infiltration ability of GBM cells. In addition, most of the genes in this set encode for membrane proteins, which are, therefore, amenable to in vivo targeting approaches with conventional or recombinant antibodies.
Acknowledgements and Funding
We thank Daniela Marubbi for expert technical assistance. This work was supported by grants from Ministero della Salute and from Fondazione S. Paolo.
- Demuth T, Berens ME: Molecular mechanisms of glioma cell migration and invasion. J Neurooncol. 2004, 70 (2): 217-228. 10.1007/s11060-004-2751-6.View ArticlePubMedGoogle Scholar
- Silbergeld DL, Chicoine MR: Isolation and characterization of human malignant glioma cells from histologically normal brain. J Neurosurg. 1997, 86 (3): 525-531. 10.3171/jns.1997.86.3.0525.View ArticlePubMedGoogle Scholar
- D'Abaco GM, Kaye AH: Integrins: molecular determinants of glioma invasion. J Clin Neurosci. 2007, 14 (11): 1041-1048. 10.1016/j.jocn.2007.06.019.View ArticlePubMedGoogle Scholar
- Harpold HL, Alvord EC, Swanson KR: The evolution of mathematical modeling of glioma proliferation and invasion. J Neuropathol Exp Neurol. 2007, 66 (1): 1-9. 10.1097/nen.0b013e31802d9000.View ArticlePubMedGoogle Scholar
- Hoelzinger DB, Demuth T, Berens ME: Autocrine factors that sustain glioma invasion and paracrine biology in the brain microenvironment. J Natl Cancer Inst. 2007, 99 (21): 1583-1593. 10.1093/jnci/djm187.View ArticlePubMedGoogle Scholar
- Giese A, Bjerkvig R, Berens ME, Westphal M: Cost of migration: invasion of malignant gliomas and implications for treatment. J Clin Oncol. 2003, 21 (8): 1624-1636. 10.1200/JCO.2003.05.063.View ArticlePubMedGoogle Scholar
- Nakada M, Nakada S, Demuth T, Tran NL, Hoelzinger DB, Berens ME: Molecular targets of glioma invasion. Cell Mol Life Sci. 2007, 64 (4): 458-478. 10.1007/s00018-007-6342-5.View ArticlePubMedGoogle Scholar
- Park JB, Kwak HJ, Lee SH: Role of hyaluronan in glioma invasion. Cell Adh Migr. 2008, 2 (3): 202-207. 10.4161/cam.2.3.6320.View ArticlePubMedPubMed CentralGoogle Scholar
- Salhia B, Tran NL, Symons M, Winkles JA, Rutka JT, Berens ME: Molecular pathways triggering glioma cell invasion. Expert Rev Mol Diagn. 2006, 6 (4): 613-626. 10.1586/14737188.8.131.523.View ArticlePubMedGoogle Scholar
- Castriconi R, Daga A, Dondero A, Zona G, Poliani PL, Melotti A, Griffero F, Marubbi D, Spaziante R, Bellora F, et al: NK cells recognize and kill human glioblastoma cells with stem cell-like properties. J Immunol. 2009, 182 (6): 3530-3539. 10.4049/jimmunol.0802845.View ArticlePubMedGoogle Scholar
- Gangemi RM, Griffero F, Marubbi D, Perera M, Capra MC, Malatesta P, Ravetti GL, Zona GL, Daga A, Corte G: SOX2 silencing in glioblastoma tumor-initiating cells causes stop of proliferation and loss of tumorigenicity. Stem Cells. 2009, 27 (1): 40-48. 10.1634/stemcells.2008-0493.View ArticlePubMedGoogle Scholar
- Gary SC, Hockfield S: BEHAB/brevican: an extracellular matrix component associated with invasive glioma. Clin Neurosurg. 2000, 47: 72-82.PubMedGoogle Scholar
- Lyons SA, Chung WJ, Weaver AK, Ogunrinu T, Sontheimer H: Autocrine glutamate signaling promotes glioma cell invasion. Cancer Res. 2007, 67 (19): 9463-9471. 10.1158/0008-5472.CAN-07-2034.View ArticlePubMedPubMed CentralGoogle Scholar
- Viapiano MS, Hockfield S, Matthews RT: BEHAB/brevican requires ADAMTS-mediated proteolytic cleavage to promote glioma invasion. J Neurooncol. 2008, 88 (3): 261-272. 10.1007/s11060-008-9575-8.View ArticlePubMedPubMed CentralGoogle Scholar
- Griffero F, Daga A, Marubbi D, Capra MC, Melotti A, Pattarozzi A, Gatti M, Bajetto A, Porcile C, Barbieri F, et al: Different response of human glioma tumor-initiating cells to epidermal growth factor receptor kinase inhibitors. J Biol Chem. 2009, 284 (11): 7138-7148.View ArticlePubMedPubMed CentralGoogle Scholar
- Monticone M, Biollo E, Fabiano A, Melotti A, Corte G, Fabbi M, Daga A, Romeo F, Maffei M, Giaretti W, et al: Cell Cultures Used in Studies Focused on Targeting Glioblastoma Tumor-Initiating Cells - Response. Molecular Cancer Research. 2010, 8 (2): 291-View ArticleGoogle Scholar
- Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, Speed TP: Exploration, normalization, and summaries of high density oligonucleotide array probe level data. Biostatistics. 2003, 4 (2): 249-264. 10.1093/biostatistics/4.2.249.View ArticlePubMedGoogle Scholar
- Saeed AI, Sharov V, White J, Li J, Liang W, Bhagabati N, Braisted J, Klapa M, Currier T, Thiagarajan M, et al: TM4: a free, open-source system for microarray data management and analysis. Biotechniques. 2003, 34 (2): 374-378.PubMedGoogle Scholar
- Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci U S A. 2001, 98 (9): 5116-5121. 10.1073/pnas.091062498.View ArticlePubMedPubMed CentralGoogle Scholar
- Eisen MB, Spellman PT, Brown PO, Botstein D: Cluster analysis and display of genome-wide expression patterns. Proc Natl Acad Sci U S A. 1998, 95 (25): 14863-14868. 10.1073/pnas.95.25.14863.View ArticlePubMedPubMed CentralGoogle Scholar
- Telenius H, Pelmear AH, Tunnacliffe A, Carter NP, Behmel A, Ferguson-Smith MA, Nordenskjold M, Pfragner R, Ponder BA: Cytogenetic analysis by chromosome painting using DOP-PCR amplified flow-sorted chromosomes. Genes Chromosomes Cancer. 1992, 4 (3): 257-263. 10.1002/gcc.2870040311.View ArticlePubMedGoogle Scholar
- Pinkel D, Straume T, Gray JW: Cytogenetic analysis using quantitative, high-sensitivity, fluorescence hybridization. Proc Natl Acad Sci U S A. 1986, 83 (9): 2934-2938. 10.1073/pnas.83.9.2934.View ArticlePubMedPubMed CentralGoogle Scholar
- Verhaak RG, Hoadley KA, Purdom E, Wang V, Qi Y, Wilkerson MD, Miller CR, Ding L, Golub T, Mesirov JP, et al: Integrated genomic analysis identifies clinically relevant subtypes of glioblastoma characterized by abnormalities in PDGFRA, IDH1, EGFR, and NF1. Cancer Cell. 2010, 17 (1): 98-110. 10.1016/j.ccr.2009.12.020.View ArticlePubMedPubMed CentralGoogle Scholar
- Rebhan M, Chalifa-Caspi V, Prilusky J, Lancet D: GeneCards: integrating information about genes, proteins and diseases. Trends in genetics : TIG. 1997, 13 (4): 163-10.1016/S0168-9525(97)01103-7.View ArticlePubMedGoogle Scholar
- Hosack DA, Dennis G, Sherman BT, Lane HC, Lempicki RA: Identifying biological themes within lists of genes with EASE. Genome Biol. 2003, 4 (10): R70-10.1186/gb-2003-4-10-r70.View ArticlePubMedPubMed CentralGoogle Scholar
- Lange K, Kammerer M, Saupe F, Hegi ME, Grotegut S, Fluri E, Orend G: Combined lysophosphatidic acid/platelet-derived growth factor signaling triggers glioma cell migration in a tenascin-C microenvironment. Cancer Res. 2008, 68 (17): 6942-6952. 10.1158/0008-5472.CAN-08-0347.View ArticlePubMedGoogle Scholar
- Horiguchi H, Hirose T, Sano T, Nagahiro S: Loss of chromosome 10 in glioblastoma: relation to proliferation and angiogenesis. Pathol Int. 1999, 49 (8): 681-686. 10.1046/j.1440-1827.1999.00934.x.View ArticlePubMedGoogle Scholar
- Tchoghandjian A, Baeza-Kallee N, Beclin C, Metellus P, Colin C, Ducray F, Adelaide J: Rougon G. 2011, Figarella-Branger D: Cortical and Subventricular ZoneGlioblastoma-Derived Stem-Like Cells Display Different Molecular Profilesand Differential In Vitro and In Vivo Properties. Annals of surgical oncology,Google Scholar
- Gunther HS, Schmidt NO, Phillips HS, Kemming D, Kharbanda S, Soriano R, Modrusan Z, Meissner H, Westphal M, Lamszus K: Glioblastoma-derived stem cell-enriched cultures form distinct subgroups according to molecular and phenotypic criteria. Oncogene. 2008, 27 (20): 2897-2909. 10.1038/sj.onc.1210949.View ArticlePubMedGoogle Scholar
- Gary SC, Zerillo CA, Chiang VL, Gaw JU, Gray G, Hockfield S: cDNA cloning, chromosomal localization, and expression analysis of human BEHAB/brevican, a brain specific proteoglycan regulated during cortical development and in glioma. Gene. 2000, 256 (1–2): 139-147.View ArticlePubMedGoogle Scholar
- Hu B, Kong LL, Matthews RT, Viapiano MS: The proteoglycan brevican binds to fibronectin after proteolytic cleavage and promotes glioma cell motility. J Biol Chem. 2008, 283 (36): 24848-24859. 10.1074/jbc.M801433200.View ArticlePubMedPubMed CentralGoogle Scholar
- Matthews RT, Gary SC, Zerillo C, Pratta M, Solomon K, Arner EC, Hockfield S: Brain-enriched hyaluronan binding (BEHAB)/brevican cleavage in a glioma cell line is mediated by a disintegrin and metalloproteinase with thrombospondin motifs (ADAMTS) family member. J Biol Chem. 2000, 275 (30): 22695-22703. 10.1074/jbc.M909764199.View ArticlePubMedGoogle Scholar
- Nakada M, Miyamori H, Kita D, Takahashi T, Yamashita J, Sato H, Miura R, Yamaguchi Y, Okada Y: Human glioblastomas overexpress ADAMTS-5 that degrades brevican. Acta Neuropathol. 2005, 110 (3): 239-246. 10.1007/s00401-005-1032-6.View ArticlePubMedGoogle Scholar
- Viapiano MS, Matthews RT, Hockfield S: A novel membrane-associated glycovariant of BEHAB/brevican is up-regulated during rat brain development and in a rat model of invasive glioma. J Biol Chem. 2003, 278 (35): 33239-33247. 10.1074/jbc.M303480200.View ArticlePubMedGoogle Scholar
- Joo KM, Kim SY, Jin X, Song SY, Kong DS, Lee JI, Jeon JW, Kim MH, Kang BG, Jung Y, et al: Clinical and biological implications of CD133-positive and CD133-negative cells in glioblastomas. Laboratory investigation; a journal of technical methods and pathology. 2008, 88 (8): 808-815. 10.1038/labinvest.2008.57.View ArticlePubMedGoogle Scholar
- Lottaz C, Beier D, Meyer K, Kumar P, Hermann A, Schwarz J, Junker M, Oefner PJ, Bogdahn U, Wischhusen J, et al: Transcriptional profiles of CD133+ and CD133- glioblastoma-derived cancer stem cell lines suggest different cells of origin. Cancer Res. 2010, 70 (5): 2030-2040. 10.1158/0008-5472.CAN-09-1707.View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/12/358/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.