RNA expression patterns in serum microvesicles from patients with glioblastoma multiforme and controls
© Noerholm et al; licensee BioMed Central Ltd. 2012
Received: 30 October 2011
Accepted: 17 January 2012
Published: 17 January 2012
RNA from exosomes and other microvesicles contain transcripts of tumour origin. In this study we sought to identify biomarkers of glioblastoma multiforme in microvesicle RNA from serum of affected patients.
Microvesicle RNA from serum from patients with de-novo primary glioblastoma multiforme (N = 9) and normal controls (N = 7) were analyzed by microarray analysis. Samples were collected according to protocols approved by the Institutional Review Board. Differential expressions were validated by qRT-PCR in a separate set of samples (N = 10 in both groups).
Expression profiles of microvesicle RNA correctly separated individuals in two groups by unsupervised clustering. The most significant differences pertained to down-regulated genes (121 genes > 2-fold down) in the glioblastoma multiforme patient microvesicle RNA, validated by qRT-PCR on several genes. Overall, yields of microvesicle RNA from patients was higher than from normal controls, but the additional RNA was primarily of size < 500 nt. Gene ontology of the down-regulated genes indicated these are coding for ribosomal proteins and genes related to ribosome production.
Serum microvesicle RNA from patients with glioblastoma multiforme has significantly down-regulated levels of RNAs coding for ribosome production, compared to normal healthy controls, but a large overabundance of RNA of unknown origin with size < 500 nt.
KeywordsExosomes Microvesicles Microarray Biomarkers Serum Glioma
Cancer molecular diagnostics is becoming increasingly important with the accumulating knowledge of the molecular mechanisms underlying various types of cancers and the implications for treatment option selection and prognosis. For patients with glioblastoma multiforme (GBM), treatment planning currently takes into account radiographic imaging, which documents volume and location of disease , and in some cases mutational analysis , methylation status of genomic DNA with particular emphasis on the DNA repair gene for methyl guanidine methyl transferase (MGMT ) and gene expression patterns of the tumour, which allows the broad categorization of tumours that are histologically similar into molecular subtypes . To date most molecular studies have utilized primary explant cultures or frozen, formalin fixed tumour tissue derived at the time of surgical resection. These methods have the disadvantage that the part of the tumour specimen chosen for analysis may not represent the rest of the tumour, and the molecular profile of the recurrent tumour may be very different from the original biopsy. It would be very useful to have a way to monitor and evaluate the tumour gene expression pattern over time in a non-invasive assay, such as through a blood sample test. Over the past few years a growing list of studies has reported on the ability to use expression profiling tests on RNA derived from blood samples to differentiate between healthy controls and patients with certain types of cancer [5–8], to classify separate patient populations  or to predict clinical outcome . The ability to conduct nucleic acid expression profiling assays on a blood sample rather than on tumours has a wide range of implications for patient welfare, including the ability to conduct longitudinal disease monitoring in situations where tumour tissue is not easily accessible or one is trying to sample metastatic cancer. Because the blood harbors nucleic acid of both tumour and non-tumour origin, it is possible that this approach may capture not only direct nucleic acid changes seen in the tumour cells, but also a component of the host response to the presence of tumour. For example, studies so far have reported on RNA extracted from Peripheral Blood Mononuclear Cells (PBMC) or other fractions of circulating blood cells where changes in the cell RNA profile appears to represent the host's response to the malignancy  rather than the tumour itself. Different groups have isolated RNA from circulating tumour cells  and from cell-free body fluids . Given the short half-life of unprotected RNA in serum , it is likely that most of the cell-free RNA is protected in the exosomes/microvesicle fraction or in the case of microRNAs (miRNAs) by protein complexes in the blood [14, 15].
Microvesicles are very stable and can protect cell-free RNA stored in the freezer for many years. This is a great advantage compared to analyzing circulating tumour cells where the blood needs to be processed within hours of collection. In addition, circulating tumour cells have not yet been described in glioma patients . In this study, full microarray analysis was carried out on mRNA isolated from serum microvesicles (including exosomes and other shedding microvesicles ) from GBM patients and controls to test the hypothesis that this mRNA could be used to reflect tumour-associated changes in the exosomal/microvesicle fraction of serum RNA. RNA species showing differential expression were chosen for quantitative reverse transcriptase (qRT-PCR) validation. This study is the first to report the ability to differentiate GBM patients from normal controls based on a gene expression blood test and the first to report differential expression analysis using RNA extracted from exosomes/microvesicles isolated from clinical patient serum samples, as compared to controls.
Blood samples from patients diagnosed with de-novo primary GBM were collected immediately prior to surgery (before opening of the dura mater) into a BD Vacutainer SST (#367985) at Massachusetts General Hospital (MGH). Patients were following standard of care at MGH, including fasting prior to surgery and most of the patients were treated with steroids to alleviate vasogenic edema and pain. Blood from normal healthy controls was collected from de-identified volunteers recruited at the MGH blood bank. All samples were collected with informed consent according to the appropriate protocols approved by the Institutional Review Board at MGH. The blood was left to clot for 30 min at room temperature (r.t.) and serum was isolated, according to manufacturer's recommendations within 2 h of collection. Serum was filtered by slowly passing it through a 0.8 μm syringe filter (Millopore, Billerica, MA, USA) and aliquoted into 1.8 mL cryotubes (Fisher Scientific, Waltham, MA, USA) and stored at -80°C until use.
Isolation of microvesicle RNA
Isolation of RNA from microvesicles was performed as previously described  with a few modifications. Briefly, 1 mL of serum was transferred to an ultracentrifuge tube, diluted 1:3 with cold PBS and centrifuged at 120,000 g for 80 min at 8°C and the supernatant was carefully aspirated off without disturbing the microvesicle pellet. The pellet was resuspended, treated for 15 min with 4 U of DNase I (Ambion, Austin, TX, USA) (in 25 μL of the accompanying buffer), 700 μL miRNeasy lysis buffer (Qiazol Reagent) (Qiagen, Valencia, CA, USA) was then added to the tube and the RNA was isolated following the manufacturer's recommendations. After elution of the RNA from the column in 30 μL nuclease-free water (Ambion), the RNA was precipitated by adding 2.5 volumes 100% EtOH, 1/10 3 M sodium acetate (pH 5.2) and incubated at -20°C for 1 h. Samples were then centrifuged for 20 min at 16,000 g and the supernatant was removed. The pellet was left to dry at r.t. and dissolved in 14 μL nuclease-free water and stored at -80°C until needed. RNA quality and concentration was assessed with the Agilent Bioanalyzer RNA Pico Chip and the Nanodrop 2000 (Thermo Scientific, Wilmington, DE, USA).
Linear amplification and array hybridization
Linear amplification and hybridization to Agilent microarrays was carried out by Miltenyi, according to manufacturer's recommendations. Briefly, exoRNA was linearly amplified and fluorescently labeled with Cy3 using Low Input Quick Amp Labeling Kit (Agilent, Santa Clara, CA, USA) and 1.4 μg amplified RNA was hybridized to Agilent 4×44K Human Microarrays, washed and scanned. Raw data was generated by image analysis using Feature Extraction (Agilent).
Microarray data analysis
The raw data exported from Agilent Feature Extraction v9.1 was pre-processed and normalized using R/Bioconductor and the packages limma, Agi4×44PreProcess and vsn (see R script in the Additional file 1 information for details). Where not mentioned in the text, the data from quartile normalization after background subtraction were used.
To reduce the risk that the normalization procedure introduced unintended biases or artifacts, we normalized the data in three different ways using: 1) variance stabilized normalization (VSN); 2) quartile normalization with background subtraction; and 3) quartile normalization without background subtraction. Although there were small differences between the three methods many of the same genes turned out to be dysregulated regardless of the applied normalization.
Clustering analysis, heat maps and Linear Discriminant Analysis of the normalized data was done using dChip (http://biosun1.harvard.edu/complab/dchip). Normalized data was transferred to Excel and filtered with various criteria, as described in the text. Gene lists of interest were uploaded and analyzed with the online Gene Ontology Tool DAVID 6.7 http://david.abcc.ncifcrf.gov/. The raw data and the quartile normalized mean data with background subtraction has been deposited in GEO with accession# GSE24084.
Reverse transcription and qPCR analysis
Twelve μL of the RNA isolated from 1 mL of serum were reverse transcribed using Superscript VILO cDNA synthesis kit (Invitrogen, Carlsbad, CA, USA), according to manufacturer's recommendations. Samples were then preamplified using the TaqMan® PreAmp Master Mix (Applied Biosystems, Carlsbad, CA, USA). Briefly, 12.5 μL of the cDNA was added to the PreAmp Master Mix together with all the genes of interest and pre- amplified for 14 cycles, according to the manufacturer's recommendations. The samples were then diluted 1:10 and TaqMan qRT-PCR was performed on all samples for all the selected genes. The amplification was performed using ABI PRISM 7500 with the following program: 50°C, 2 min; 95°C, 10 min; 40 cycles of 95°C, 15 s, 60°C, 1 min on standard mode. Additional file 2: Table S1 contains a list of all the commercially available and custom made probes used.
Unsupervised clustering separates GBM patients from controls in expression array analysis
Down-and up-regulated genes passing the criteria described in "Results".
Homo sapiens thymosin-like 3 (TMSL3), mRNA [NM_183049]
Homo sapiens guanine nucleotide binding protein (G protein), gamma 11 (GNG11), mRNA [NM_004126]
Homo sapiens ribosomal protein S4, Y-linked 2 (RPS4Y2), mRNA [NM_001039567]
Homo sapiens regulator of G-protein signalling 10 (RGS10), transcript variant 1, mRNA [NM_001005339]
Homo sapiens ribosomal protein S4, Y-linked 1 (RPS4Y1), mRNA [NM_001008]
Homo sapiens beta-2-microglobulin (B2M), mRNA [NM_004048]
Homo sapiens chemokine (C-C motif) ligand 5 (CCL5), mRNA [NM_002985]
Homo sapiens mRNA; cDNA DKFZp586A0722 (from clone DKFZp586A0722). [AL049447]
Homo sapiens glutathione peroxidase 1 (GPX1), transcript variant 2, mRNA [NM_201397]
Homo sapiens ribosomal protein S29 (RPS29), transcript variant 1, mRNA [NM_001032]
Homo sapiens immediate early response 2 (IER2), mRNA [NM_004907]
Homo sapiens RAP1B, member of RAS oncogene family (RAP1B), transcript variant 1, mRNA [NM_015646]
Homo sapiens monocyte to macrophage differentiation-associated (MMD), mRNA [NM_012329]
Homo sapiens MYC associated factor × (MAX), transcript variant 3, mRNA [NM_145113]
Homo sapiens transmembrane protein 111 (TMEM111), mRNA [NM_018447]
Homo sapiens ribosomal protein S29 (RPS29), transcript variant 1, mRNA [NM_001032]
Homo sapiens ribosomal protein L19 (RPL19), mRNA [NM_000981]
Homo sapiens ribosomal protein L13 (RPL13), transcript variant 2, mRNA [NM_033251]
PREDICTED: Homo sapiens similar to 40S ribosomal protein S6 (LOC392497), mRNA [XR_018138]
Homo sapiens myosin regulatory light chain MRCL3 (MRCL3), mRNA [NM_006471]
Homo sapiens ribosomal protein L30 (RPL30), mRNA [NM_000989]
Homo sapiens ribosomal protein S27 (metallopanstimulin 1) (RPS27), mRNA [NM_001030]
40S ribosomal protein S21. [Source:Uniprot/SWISSPROT;Acc:P63220] [ENST00000337102]
Homo sapiens replication protein A1, 70 kDa (RPA1), mRNA [NM_002945]
oe37f10.y1 Human keratoconus cornea, unamplified, od [CV575560]
Homo sapiens ring finger and KH domain containing 1 (RKHD1), mRNA [NM_203304]
Homo sapiens zinc finger protein 784 (ZNF784), mRNA [NM_203374]
Leukocyte elastase inhibitor (LEI) (Serpin B1) (Monocyte/neutrophil elastase inhibitor) (M/NEI) (EI).
PREDICTED: Homo sapiens similar to TBP-associated factor 15 isoform 1 (LOC390427), mRNA [XM_372498]
Homo sapiens forkhead box D3 (FOXD3), mRNA [NM_012183]
Homo sapiens chromosome 17 open reading frame 74 (C17orf74), mRNA [NM_175734]
Homo sapiens SLIT and NTRK-like family, member 4 (SLITRK4), mRNA [NM_173078]
Homo sapiens agouti related protein homolog (mouse) (AGRP), transcript variant 1, mRNA [NM_001138]
chromosome × open reading frame 18 (CXorf18), misc RNA [RefSeq_dna;Acc:XR_018001] [ENST00000359466]
Homo sapiens tenascin XB (TNXB), transcript variant XB, mRNA [NM_019105]
ALU1_HUMAN (P39188) Alu subfamily J sequence contamination warning entry, partial (13%) [THC2755576]
Homo sapiens SIN3 homolog B, transcription regulator (yeast), complete cds. [BC063531]
Homo sapiens cytokine receptor-like factor 1 (CRLF1), mRNA [NM_004750]
Homo sapiens zinc finger protein 219 (ZNF219), mRNA [NM_016423]
ALU2_HUMAN (P39189) Alu subfamily SB sequence contamination warning entry, partial (4%) [THC2654039]
AA418814 Soares_NhHMPu_S1 Homo sapiens cDNA clone IMAGE:767978 3', mRNA sequence [AA418814]
Homo sapiens trinucleotide repeat containing 4 (TNRC4), mRNA [NM_007185]
Homo sapiens cDNA FLJ25506 fis, clone CBR05185. [AK098372]
Homo sapiens trinucleotide repeat containing 6B (TNRC6B), transcript variant 1, mRNA [NM_015088]
Homo sapiens homeobox A4 (HOXA4), mRNA [NM_002141]
Homo sapiens interleukin 26 (IL26), mRNA [NM_018402]
Gene ontology analysis
Normalization of array data can be done in a number of different ways, each having its own advantages and disadvantages. To reduce the risk that the normalization procedure introduced unintended biases or artifacts, we normalized the data in three different ways, as described in Methods. For gene ontology analysis we used the top two hundred genes that passed the criteria p < 0.05 after FDR-correction in all three normalization procedures mentioned above (Additional file 2: Table S2 and Table S3). All genes listed in the tables of this paper met the same selection criteria.
Down-regulated genes in serum microvesicles from GBM patients
We analyzed the up- and down-regulated genes independently to try to identify brain- or GBM-specific genes up-regulated in the GBM patients due to the microvesicle shedding activity of the tumour, as well as genes down-regulated in normal cells, such as platelets, lymphocytes and endothelial cells  that also shed microvesicles into the blood. Gene ontology analysis of the 200 most down-regulated genes (from -4.21 to -1.92-fold, Additional file 2: Table S2) using the online tool, DAVID  revealed that an overwhelming majority of the resulting GO terms are related to ribosome functions, e.g. "Ribosome", "Translational Elongation" and "RNA Binding". The most significant cluster had an enrichment score of 54.82 and contained 115 of the 200 input genes, most of which were mRNAs coding for ribosomal proteins in both the large and small subunits of the ribosome (e.g. RPL11, RPS29, RPLP1, RPS27A, etc.). These genes have previously been shown to be very highly expressed in lymphocytes relative to other blood cells [23, 24]. Evaluation of the percentage of lymphocytes in the white blood cell (WBC) counts in the GBM patients used in this study, showed that 7 out of the 9 patients used for array analysis and 9 out of the 10 patients used for qPCR validation had values of 3-19% lymphocytes, which is well below the reference interval for normal healthy individuals of 22-28%  and is consistent with immunosuppression observed in most GBM patients . Unfortunately, the de-identified normal control samples used for this study were obtained through a blood bank so no detailed information is available on their blood cell counts.
Up-regulated genes in microvesicles from GBM patients
When we left out the down-regulated genes and attempted unsupervised clustering using only genes that were up-regulated in microvesicles from GBM patients (Table 1/Additional file 2: Table S3), no clear separation of the GBMs from the normal controls was evident. Further, gene ontology analysis of the up-regulated genes in Table 1/Additional file 2: Table S3 resulted in clusters and GO terms with enrichment scores substantially lower (< 3) than for the down-regulated genes in Table 1/Additional file 2: Table S2 (> 50).
Validation analysis by qRT-PCR
We moved on to see if some of the dysregulated genes from the microarray could be validated as dysregulated by qRT-PCR in an independent set of 10 patients and 10 controls. We selected individual up-regulated genes from the lists in Table 1/Additional file 2: Table S3 in search of markers associated with GBM tumours. We sought to identify mRNAs up-regulated in the microvesicles in peripheral blood from GBM patients, which were expressed at low to undetectable levels in blood cells, and elevated in either GBM or normal brain cells. We reasoned that since the majority of microvesicles in blood are believed to arise from platelets  and other blood cells, the RNA expression of these cells could serve as an approximation of the normal blood microvesicle RNA profile. We compared public datasets in Gene Expression Omnibus (GEO) for GBM (GSE15824) to those of platelets (GSE11524) and PBMCs (GSE22224) from normal controls to derive a list of "GBM-like" genes highly expressed in GBM, but low or absent in platelets and PBMCs and similarly a list of "blood cell-like" genes highly expressed in platelets and PBMCs. When we compared the "GBM-like" genes to the genes found to be up-regulated in GBM patient blood microvesicle RNA in our own dataset, we found no overlap to support GBM tumour origin of these RNAs. On the other hand, many of the "blood cell-like" genes were the same as those found in our study to be down-regulated in GBM patients relative to normal controls (data not shown).
Our findings demonstrate the feasibility of genome wide gene expression analysis of RNA extracted from the cell-free microvesicle fraction (exoRNA) of frozen biobanked serum samples, and that clear differences can be observed between the exoRNA expression profiles from GBM patients and normal controls. The two major differences between exoRNA from GBM patients vs. normal controls were found to be: 1) a significantly reduced level of mRNAs encoding ribosomal protein genes; and 2) an overall significant up-regulation of RNA amounts in the serum of GBM patients, which seems to be primarily due to a larger fraction < 300 nt and could not be attributed to any of the mRNAs investigated by qRT-PCR.
ExoRNA from tumour microvesicles in serum will always be diluted into the background of exoRNA coming from normal non-malignant cells. Depending on the purpose of the study, normal exoRNA may also generate useful information. Previous studies have indicated that RNA expression patterns in blood cells can change as a response to tumours . To look specifically at the genes dysregulated in the tumour, one would need to enrich the tumour specific microvesicle fraction before extracting the RNA.
Microvesicles are shed by many different cell types into the blood stream, including circulating WBCs, platelets, endothelial cells, and dendritic cells . Shedding has been shown to be higher from tumour cells in culture compared to normal fibroblasts  and tumour-derived microvesicle concentrations in serum increase as a function of increased malignancy [29, 30]. However, microvesicles of tumour origin circulating in the blood are mixed with microvesicles from other cell sources and it is likely that tumour microvesicles in many cases constitute a relatively small fraction of the total population.
There is a large contribution of microvesicles from normal cells and platelets in blood. In control plasma, for example, platelet-derived microvesicles have been reported to constitute the majority of all microvesicles  and although tumour-derived microvesicles increase in serum as a function of degree of malignancy in ovarian cancer  and that the total number of microvesicles in blood increase with disease progression , it is not clear what proportion of the increased number of vesicles is coming from the tumour vs. from normal cells as a response to the tumour. Thus, down-regulated genes in GBM exoRNA could either be caused by altered biological processes in normal cells as a response to or consequence of the disease, or by large numbers of tumour microvesicles with low levels of particular RNAs.
A major fraction of the genes we observed to be significantly down-regulated in serum microvesicles from GBM patients in this study code for various ribosomal proteins (e.g. RPL11, RPS29, RPLP1, RPS27A, etc.). However, these genes were not found to be among the most abundant transcripts in platelets in a study by , which we also confirmed in our analysis of the public dataset from GEO (GSE11524). We speculate that these RNAs would therefore also not be abundant in the platelet-derived microvesicles, which constitute a substantial fraction of the plasma microvesicles . However, these genes have been shown to be very highly expressed in lymphocytes relative to other blood cells [23, 24] and were also found to be up-regulated in our analysis of the public dataset for PBMC (GSE22224) relative to GBM. This led us to speculate that the lower abundance of these mRNAs in circulating microvesicles from GBM patients may be a consequence of a reduced level of lymphocyte-derived microvesicles in patients as compared to normal controls, since cancer patients are known to often be immune-compromised . Evaluation of the GBM patients in this study confirmed that their lymphocyte counts were lower than the normal reference level interval, but because we were unable to measure lymphocyte counts in the normal controls, it is impossible to completely correlate transcript down-regulation with lymphocyte depletion in the GBM patients. None of the GBM patients had received chemotherapy prior to diagnosis, but several of them were on steroids and seizure medications at the time of blood draw, which are two potential sources of lymphocytopenia. Evaluation of gene expression differences between lymphocytopenia caused by cancer and by medication would require a much larger study of more carefully chosen subjects.
The observation that many ribosomal protein genes are down-regulated in GBM serum microvesicles may prove to be a valuable contributory marker for GBM and other pathological states. Sharma et al.  also observed a reduction in many of the same genes in PBMCs including lymphocytes, from breast cancer patients compared to normal controls. The down-regulation of ribosomal protein genes is therefore unlikely to provide a specific GBM component for a diagnostic classification, but may serve as a more generic pathological indicator and can provide confirmatory support when combined with other exoRNA profiles. Interestingly, many ribosomal proteins have extra-ribosomal functions that go beyond the function of protein biosynthesis. Ribosomal protein genes have been shown to have important regulatory functions also in cancer cells. Some ribosomal proteins are tumour suppressors , oncogenes or have regulatory functions in tumour progression, invasion and metastasis .
The "ribosome" GO-clusters from down-regulated genes had very high enrichment scores from DAVID (> 50) compared to any of the clusters obtained with up-regulated genes (< 2.5), confirming that the patterns of down-regulated genes were much more significant. For gene ontology analysis it is important to analyze a relatively large number of genes (e.g. 100 ) in order to avoid stochastic errors, e.g. if a list of 10 genes is analyzed, the presence of a single gene (10%) of a certain class or belonging to a certain pathway will appear overrepresented compared to the prevalence in the genome, when in fact this is just a sampling artifact. The lack of strong GO associations between the up-regulated genes in our study does not mean that individual genes might not be strongly dysregulated, significant and predictive in GBM exoRNA, but simply that the up-regulated genes do not appear to be associated by already known relationships. Up-regulation of specific genes in exoRNA from GBM patients could be a biological response from normal cells to the presence of the tumour, but they could also be specifically derived from the tumour cells and as such be useful as markers of the tumour. However, identifying uniquely elevated levels of tumour-derived exoRNAs in a high background of normal exoRNA from other serum microvesicle sources is more challenging than identifying systemic changes, such as the depressed lymphocyte count and ribosomal protein RNA levels.
The observation that the up-regulated genes chosen for qRT-PCR validation was not verified may be attributed to the fact that we used two different sets of serum samples for the microarray analysis and the qRT-PCR validation, making borderline significant increases difficult to confirm. Similarly, the positive validation of the down-regulated genes in a separate set of samples makes this observation stronger.
We observed that, coinciding with the significantly increased amounts of microvesicle RNA in GBM patient serum compared to normal controls, the expression levels of all the genes we investigated by qRT-PCR were lower than in normal controls when normalized to total RNA amount (Figure 5). It is likely that this is a general trend across most mRNAs, and since all these genes are down-regulated other transcripts must be up-regulated to account for the overall increase in RNA amounts.
The Bioanalyzer profiles of the different RNA serum samples suggest that this up-regulation of RNA in the GBM patients mainly falls in the range of < 300 nt. This could explain why the majority of the tumour exoRNA transcripts appeared down-regulated on the array. If the < 300 nt RNAs do not contain sequences belonging to the coding genes, they would not have been picked up by the microarray and the mRNA-fraction recognized by the capture probes would be relatively smaller since the same amount of input RNA was used for hybridization. There are many RNA species that are not covered by the microarray, including miRNA (mature and precursors), repetitive elements and other non-coding RNAs, as well as single stranded DNA [28, 34].
Balaj et al.  showed that human endogenous retroviruses and other transposable elements, as well as single stranded DNA fragments are very abundant in microvesicles from cancer cells, as compared to normal fibroblasts, including Alu, LINE and HERV sequences and it is likely that these and other non-coding RNAs are contributing to the increased amounts of RNA we observe in GBM patient serum (single stranded DNA is also detected on the Bioanalyzer RNA chip). We tested a single transposable element Alu-Y by qRT-PCR and found it to be extremely abundant in the microvesicles, but did not observe any differential expression between the GBMs and normal controls. The human genome consists of about 40% retrotransposon sequences , and there is an increasing number of publications showing the dysregulation of transposable elements and other non-coding RNAs in cancer [36–38].
The very existence of the peak of < 500 nt RNAs in GBM serum exoRNA make it a fertile ground for biomarker discovery, and it warrants further investigation to establish the exact nature and distribution of the nucleic acid species contained in this fraction. A number of regulatory non-coding RNAs are transcribed off coding elements . Collectively, down-regulation of specific RNAs and up-regulation of RNA levels in exosomes/microvesicles from serum of GBM patients, as compared to controls provides promising biomarkers. In addition, in other studies, tumour mutant mRNAs  and elevated oncogene mRNA  have been detected in serum exosomes/microvesicles from human GBM patients and mice bearing medulloblastoma tumours, respectively. Tumour mutant RNA from prostate cancer patients can also be found in urine exosomes/microvesicles . More detailed analysis of exoRNA released by tumour cells into serum should be possible as isolation methods are developed with tumour specific surface markers for different types of cancer, e.g. with magnetic activated sorting  and microfluidic capture .
Currently, there are no biomarkers for gliomas to distinguish tumor recurrence from radiation necrosis or to monitor tumor response to therapy. Here we demonstrate that frozen, biobanked glioblastoma serum microvesicles contain RNA in amounts and qualities sufficient to perform qPCR and microarray analysis, and we show that glioma exoRNA profiles are unique compared to tumor free individuals. The most significant expression differences pertained to down-regulated genes in the GBM patient exoRNA, which we were able to validate by qRT-PCR. However, overall yields of exoRNA from GBM patient serum microvesicles was higher than yields from normal controls, but the additional RNA was primarily of size < 500 nt, evident by a peak of small RNAs on the bioanalyzer profiles. Gene ontology analysis of the down-regulated genes indicated these are primarily mRNAs coding for ribosomal proteins and other genes related to ribosome production. The observation may be explained by a reduced contribution of exoRNA from lymphocytes, which have previously been reported to have high expression levels of these transcripts. This theory was supported by the lymphocyte counts in the blood of the patients investigated in this study, which was substantially lower than the normal reference interval.
Methyl guanidine methyl transferase
Peripheral blood mononuclear cells
Variance stabilized normalization
White blood cell
Gene expression omnibus
False discovery rate
Quantitative reverse transcriptase PCR.
We would like to thank Lori LoGuidice for kind assistance with RNA extractions and Suzanne McDavitt for review and formatting of the manuscript. This work was funded by American Brain Tumor Association Translational grant (JS, LB), NIH/NCI grant CA141226 (XOB), CA156009 (XOB) and NIH/NCI grant CA069246 (BC, AS and FHH). Brain Tumor Funders Collaborative and Harvard Catalyst (BSC). We gratefully acknowledge the Neuroscience Center PCR Core facility, which is funded by PHS grant P30NS045776.
- Dean BL, Drayer BP, Bird CR, Flom RA, Hodak JA, Coons SW, Carey RG: Gliomas: classification with MR imaging. Radiology. 1990, 174: 411-415.View ArticlePubMedGoogle Scholar
- Hartmann C, Hentschel B, Tatagiba M, Schramm J, Schnell O, Seidel C, Stein R, Reifenberger G, Pietsch T, von Deimling A, et al: Molecular Markers in Low-Grade Gliomas: Predictive or Prognostic?. Clin Cancer Res. 2011, 17: 4588-4599. 10.1158/1078-0432.CCR-10-3194.View ArticlePubMedGoogle Scholar
- Kreth S, Thon N, Eigenbrod S, Lutz J, Ledderose C, Egensperger R, Tonn JC, Kretzschmar HA, Hinske LC, Kreth FW: O-methylguanine-DNA methyltransferase (MGMT) mRNA expression predicts outcome in malignant glioma independent of MGMT promoter methylation. PLoS One. 2011, 6: e17156-10.1371/journal.pone.0017156.View ArticlePubMedPubMed CentralGoogle Scholar
- Cancer Genome Atlas Research Network: Comprehensive genomic characterization defines human glioblastoma genes and core pathways. Nature. 2008, 455: 1061-1068. 10.1038/nature07385.View ArticleGoogle Scholar
- Sharma P, Sahni NS, Tibshirani R, Skaane P, Urdal P, Berghagen H, Jensen M, Kristiansen L, Moen C, Sharma P, et al: Early detection of breast cancer based on gene-expression patterns in peripheral blood cells. Breast Cancer Res. 2005, 7: R634-R644. 10.1186/bcr1203.View ArticlePubMedPubMed CentralGoogle Scholar
- Osman I, Bajorin DF, Sun TT, Zhong H, Douglas D, Scattergood J, Zheng R, Han M, Marshall KW, Liew CC: Novel blood biomarkers of human urinary bladder cancer. Clin Cancer Res. 2006, 12: 3374-3380. 10.1158/1078-0432.CCR-05-2081.View ArticlePubMedGoogle Scholar
- Han M, Liew CT, Zhang HW, Chao S, Zheng R, Yip KT, Song ZY, Li HM, Geng XP, Zhu LX, et al: Novel blood-based, five-gene biomarker set for the detection of colorectal cancer. Clin Cancer Res. 2008, 14: 455-460. 10.1158/1078-0432.CCR-07-1801.View ArticlePubMedGoogle Scholar
- Aarøe J, Lindahl T, Dumeaux V, Saebø S, Tobin D, Hagen N, Skaane P, Lönneborg A, Sharma P, Børresen-Dale AL: Gene expression profiling of peripheral blood cells for early detection of breast cancer. Breast Cancer Res. 2010, 12: R7-10.1186/bcr2472.View ArticlePubMedPubMed CentralGoogle Scholar
- Showe MK, Vachani A, Kossenkov AV, Yousef M, Nichols C, Nikonova EV, Chang C, Kucharczuk J, Tran B, Wakeam E, et al: Gene expression profiles in peripheral blood mononuclear cells can distinguish patients with non-small cell lung cancer from patients with nonmalignant lung disease. Cancer Res. 2009, 69: 9202-9210. 10.1158/0008-5472.CAN-09-1378.View ArticlePubMedPubMed CentralGoogle Scholar
- Burczynski ME, Twine NC, Dukart G, Marshall B, Hidalgo M, Stadler WM, Logan T, Dutcher J, Hudes G, Trepicchio WL, et al: Transcriptional profiles in peripheral blood mononuclear cells prognostic of clinical outcomes in patients with advanced renal cell carcinoma. Clin Cancer Res. 2005, 11: 1181-1189.PubMedGoogle Scholar
- Maheswaran S, Haber DA: Circulating tumor cells: a window into cancer biology and metastasis. Curr Opin Genet Dev. 2010, 20: 96-99. 10.1016/j.gde.2009.12.002.View ArticlePubMedPubMed CentralGoogle Scholar
- Schmidt B, Engel E, Carstensen T, Weickmann S, John M, Witt C, Fleischhacker M: Quantification of free RNA in serum and bronchial lavage: a new diagnostic tool in lung cancer detection?. Lung Cancer. 2005, 48: 145-147. 10.1016/j.lungcan.2004.09.013.View ArticlePubMedGoogle Scholar
- Tsui NB, Ng EK, Lo YM: Stability of endogenous and added RNA in blood specimens, serum, and plasma. Clin Chem. 2002, 48: 1647-1653.PubMedGoogle Scholar
- Arroyo JD, Chevillet JR, Kroh EM, Ruf IK, Pritchard CC, Gibson DF, Mitchell PS, Bennett CF, Pogosova-Agadjanyan EL, Stirewalt DL, et al: Argonaute2 complexes carry a population of circulating microRNAs independent of vesicles in human plasma. Proc Natl Acad Sci USA. 2011, 108: 5003-5008. 10.1073/pnas.1019055108.View ArticlePubMedPubMed CentralGoogle Scholar
- Turchinovich A, Weiz L, Langheinz A, Burwinkel B: Characterization of extracellular circulating microRNA. Nucleic Acids Res. 2011, 39: 7223-7233. 10.1093/nar/gkr254.View ArticlePubMedPubMed CentralGoogle Scholar
- Bernstein JJ, Woodard CA: Glioblastoma cells do not intravasate into blood vessels. Neurosurgery. 1995, 36: 124-132. 10.1227/00006123-199501000-00016.View ArticlePubMedGoogle Scholar
- Cocucci E, Racchetti G, Meldolesi J: Shedding microvesicles: artefacts no more. Trends Cell Biol. 2009, 19: 43-51. 10.1016/j.tcb.2008.11.003.View ArticlePubMedGoogle Scholar
- Skog J, Würdinger T, van Rijn S, Meijer D, Gainche L, Curry WTJ, Carter BS, Krichevsky AM, Breakefield XO: Glioblastoma microvesicles transport RNA and protein that promote tumor growth and provide diagnostic biomarkers. Nat Cell Biol. 2008, 10: 1470-1476. 10.1038/ncb1800.View ArticlePubMedPubMed CentralGoogle Scholar
- Dobbin KK, Zhao Y, Simon RM: How large a training set is needed to develop a classifier for microarray data?. Clin Cancer Res. 2008, 14: 108-114. 10.1158/1078-0432.CCR-07-0443.View ArticlePubMedGoogle Scholar
- Benjamini Y, Hochberg Y: Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. J Royal Stat Soc, Series B (Methodolofical). 1995, 1: 289-300.Google Scholar
- Théry C, Ostrowski M, Segura E: Membrane vesicles as conveyors of immune responses. Nat Rev Immunol. 2009, 9: 581-593. 10.1038/nri2567.View ArticlePubMedGoogle Scholar
- da Huang W, Sherman BT, Lempicki RA: Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009, 4: 44-57.View ArticlePubMedGoogle Scholar
- Whitney AR, Diehn M, Popper SJ, Alizadeh AA, Boldrick JC, Relman DA, Brown PO: Individuality and variation in gene expression patterns in human blood. Proc Natl Acad Sci USA. 2003, 100: 1896-1901. 10.1073/pnas.252784499.View ArticlePubMedPubMed CentralGoogle Scholar
- Palmer C, Diehn M, Alizadeh AA, Brown PO: Cell-type specific gene expression profiles of leukocytes in human peripheral blood. BMC Genomics. 2006, 7: 115-10.1186/1471-2164-7-115.View ArticlePubMedPubMed CentralGoogle Scholar
- Stedman's Medical Dictionary. 2000, Lippincott Williams & Wilkins, 27
- Gomez GG, Kruse CA: Mechanisms of malignant glioma immune resistance and sources of immunosuppression. Gene Ther Mol Biol. 2006, 10: 133-146.PubMedPubMed CentralGoogle Scholar
- Baran J, Baj-Krzyworzeka M, Weglarczyk K, Szatanek R, Zembala M, Barbasz J, Czupryna A, Szczepanik A, Zembala M: Circulating tumour-derived microvesicles in plasma of gastric cancer patients. Cancer Immunol Immunother. 2010, 59: 841-850. 10.1007/s00262-009-0808-2.View ArticlePubMedGoogle Scholar
- Balaj L, Lessard R, Dai L, Cho Y-J, Pomeroy SL, Breakefield XO, Skog J: Tumour microvesicles contain retrotransposon elements and amplified oncogene sequences. Nat Commun. 2011, 2: 180-View ArticlePubMedPubMed CentralGoogle Scholar
- Logozzi M, De Milito A, Lugini L, Borghi M, Calabrò L, Spada M, Perdicchio M, Marino ML, Federici C, Iessi E, et al: High levels of exosomes expressing CD63 and caveolin-1 in plasma of melanoma patients. PLoS One. 2009, 4: e5219-10.1371/journal.pone.0005219.View ArticlePubMedPubMed CentralGoogle Scholar
- Taylor DD, Gercel-Taylor C: MicroRNA signatures of tumor-derived exosomes as diagnostic biomarkers of ovarian cancer. Gynecol Oncol. 2008, 110 (1): 13-21. 10.1016/j.ygyno.2008.04.033.View ArticlePubMedGoogle Scholar
- Gnatenko DV, Dunn JJ, Schwedes J, Bahou WF: Transcript profiling of human platelets using microarray and serial analysis of gene expression (SAGE). Methods Mol Biol. 2009, 496: 245-272. 10.1007/978-1-59745-553-4_16.View ArticlePubMedGoogle Scholar
- Amsterdam A, Sadler KC, Lai K, Farrington S, Bronson RT, Lees JA, Hopkins N: Many ribosomal protein genes are cancer genes in zebrafish. PLoS Biol. 2004, 2: E139-10.1371/journal.pbio.0020139.View ArticlePubMedPubMed CentralGoogle Scholar
- Lai MD, Xu J: Ribosomal proteins and colorectal cancer. Curr Genomics. 2007, 8: 43-49. 10.2174/138920207780076938.View ArticlePubMedPubMed CentralGoogle Scholar
- van der Vos KE, Balaj L, Skog J, Breakefield XO: Brain tumor microvesicles: Insights into intercellular communication in the nervous system. Cell Mol Neurobiol. 2011, 31: 949-959. 10.1007/s10571-011-9697-y.View ArticlePubMedPubMed CentralGoogle Scholar
- Goodier JL, Kazazian HHJ: Retrotransposons revisited: the restraint and rehabilitation of parasites. Cell. 2008, 135: 23-35. 10.1016/j.cell.2008.09.022.View ArticlePubMedGoogle Scholar
- White RJ: RNA polymerase III transcription and cancer. Oncogene. 2004, 23: 3208-3216. 10.1038/sj.onc.1207547.View ArticlePubMedGoogle Scholar
- Ting DT, Lipson D, Paul S, Brannigan BW, Akhavanfard S, Coffman EJ, Contino G, Deshpande V, Iafrate AJ, Letovsky S, et al: Aberrant overexpression of satellite repeats in pancreatic and other epithelial cancers. Science. 2011, 331: 593-596. 10.1126/science.1200801.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen W, Böcker W, Brosius J, Tiedge H: Expression of neural BC200 RNA in human tumours. J Pathol. 1997, 183: 345-351. 10.1002/(SICI)1096-9896(199711)183:3<345::AID-PATH930>3.0.CO;2-8.View ArticlePubMedGoogle Scholar
- Mercer TR, Dinger ME, Mattick JS: Long non-coding RNAs: insights into functions. Nat Rev Genet. 2009, 10: 155-159. 10.1038/nrg2521.View ArticlePubMedGoogle Scholar
- Nilsson J, Skog J, Nordstrand A, Baranov V, Mincheva-Nilsson L, Breakefield XO, Widmark A: Prostate cancer-derived urine exosomes: A novel approach to biomarkers for prostate cancer. Brit J Cancer. 2009, 100: 1603-1607. 10.1038/sj.bjc.6605058.View ArticlePubMedPubMed CentralGoogle Scholar
- Chen C, Skog J, Hsu C-H, Lessard R, Breakefield XO, Toner M, Irimia D: Microfluidic isolation and transcriptome analysis of serum microvesicles. Lab Chip. 2010, 5: 505-511.View ArticleGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2407/12/22/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.