Integrative Analysis of mRNA,microRNA,and Protein Correlates of Relative Cerebral Blood Volume Values in GBM Reveals the Role for Modulators of Angiogenesis and Tumor Proliferation

Abstract

Dynamic susceptibility contrast-enhanced magnetic resonance imaging is routinely used to provide hemodynamic assessment of brain tumors as a diagnostic as well as a prognostic tool. Recently, it was shown that the relative cerebral blood volume (rCBV), obtained from the contrast-enhancing as well as -nonenhancing portion of glioblastoma (GBM), is strongly associated with overall survival. In this study, we aim to characterize the genomic correlates (microRNA, messenger RNA, and protein) of this vascular parameter. This study aims to provide a comprehensive radiogenomic and radioproteomic characterization of the hemodynamic phenotype of GBM using publicly available imaging and genomic data from the Cancer Genome Atlas GBM cohort. Based on this analysis, we identified pathways associated with angiogenesis and tumor proliferation underlying this hemodynamic parameter in GBM.

Keywords

imaging-genomics perfusion imaging rCBV angiogenesis pathway analysis data integration

Introduction

Glioblastoma (GBM) is a WHO grade IV tumor with an extremely poor prognosis and a median survival of 12–15 months. Many imaging features are used to characterize GBM, such as structural magnetic resonance imaging (MRI), diffusion MRI, and perfusion MRI.¹ Several of these derived radiology-based image features (denoted radiophenotypes) have shown association with the underlying genetics and biology of the tumor.²

Identification of such noninvasive surrogates of clinical outcome and the associated molecular differences might point to a deeper understanding of the underlying disease state.^3,4 Several studies^5–7 have demonstrated that imaging phenotypes derived from MRI studies might have the potential to serve as noninvasive proxies of cancer-associated clinical outcomes. In this regard, we follow up on recent studies based on dynamic susceptibility contrast-enhanced MRI (DSC-MRI) that identifies a hemodynamic parameter associated with overall survival. Specifically, these studies^5,8 pointed to the association of rCBV_max (maximum relative cerebral blood volume of the contrast-enhancing component of GBM) and rCBV_NER (rCBV in the nonenhancing region of the tumor) with patient survival.

The Cancer Genome Atlas (TCGA)⁹ contains comprehensive genomic data (spanning expression, copy number, methylation, microRNA (miRNA), etc.). The Cancer Imaging Archive (TCIA)¹⁰ is a publicly available repository of patient-derived images across multiple modalities (single-photon emission computed tomography, computed tomography, and MRI) and multiple tumor types in many of the patients who have corresponding molecular data in TCGA. Using genomic and imaging data from the TCGA-GBM archive, we carried out a bioinformatics analysis to identify molecular correlates of imaging features associated with survival.

Based on the prognostic significance of the hemodynamic parameter (rCBV_max value, we sought to examine the molecular correlates of rCBV via differential expression max rCBV_max via differential expression analysis of messenger RNA (mRNA), miRNA, and protein levels. The intent was to identify molecular entities that are likely regulated across different layers of biological complexity: each modality (mRNA, miRNA, and protein) provides complementary information about the biology of the process of GBM tumorigenesis. Thus, we decided to combine the significantly associated genes (transcriptional correlate), miRNAs (post-transcriptional correlate), and proteins (signal transduction correlate) to identify if there was a concordant biological process (or processes) that could be inferred when these three–omic modalities were combined. Combining entities found associated with tumor hemodynamic response across different biological layers has the potential to illuminate important regulatory processes.

For this purpose, patients were dichotomized based on median rCBV_max value and then examined for differential expression of mRNA, miRNA, and protein levels. A molecular examination of the two phenotypic groups reveals that genes and miRNAs related to the molecular process of angiogenesis underlie this phenotypic transition between low and high rCBV_max groups.

Materials and Methods

Ethics Statement

All patients in this retrospective study had been previously de-identified by TCGA, a publicly available data set that contains no linkage to patient identifiers and is compliant with the Health Insurance Portability and Accountability Act; therefore, a waiver was obtained from the Henry Ford Hospital Institutional Review Board of Henry Ford Hospital, Detroit, MI. This research complied with the principles of the Declaration of Helsinki.

Patient Population

Based on the availability of perfusion (DSC) MRI image data from the TCIA and the corresponding genomic expression from TCGA, 50 GBM cases were identified (19 females, 30 males, and 1 unknown). These cases were selected to ensure consistency in glioma history and treatment profile. Specifically, patients with untreated primary GBM were selected. Clinical data for this population were obtained from the cBioPortal for Cancer Genomics.¹¹

Image Features

Radiology annotations of the 50 GBM cases and their preoperative images (from TCIA) were obtained from the TCGA Glioma Phenotype Research Group.¹⁰ This Health Insurance Portability and Accountability Act-compliant retrospective study was approved by the Henry Ford Hospital Institutional Review Board (IRB #6381). Patients underwent DSC T2*-weighted MR perfusion studies and had gene expression data available from TCGA. Fifteen cases were imaged in a 3 T scanner, while others were imaged in a 1.5 T scanner. rCBV_max was measured from a region of interest of 10 x 10 voxels⁵ placed on the hottest-appearing tumor region based on qualitative perfusion maps.

Genomic Features

Batch-corrected, Level 3 expression data for both mRNA (Affymetrix U133Av2 BI Platform: Affymetrix HT HG U133A) and miRNA (UNC Agilent Human miRNA 8x15K) were downloaded from the UT MD Anderson Genome Data Analysis Center Portal¹² version dated June 18, 2013. Protein expression from reverse-phase protein array (RPPA) was downloaded from the Cancer Proteome Portal.¹³ Forty-seven cases had miRNA and mRNA data available, whereas 15 cases had data for all three platforms. Instead of restricting our study to just the 15 common cases for differential expression analysis, we performed analysis of each platform separately and integrated the differential expression results subsequently via pathway analysis.

Defining Classes Based on Median rCBV_Max

Based on the prior finding that rCBV_max is associated with survival in a statistical significant manner,⁵ we sought to examine the molecular correlates of rCBV_max via differential expression analysis. Median rCBV_max across the entire dataset was used to divide the study population into two groups.

Differential Expression Analysis between Cutpoint-Induced Phenotype Classes

Differential expression analysis for mRNA, miRNA, and protein was done to investigate molecular differences between the median rCBV_max -induced phenotype groups. This was done using a two-sided t-test to identify genes, miRNAs, and proteins that are differentially expressed between the two classes. Because of the exploratory nature of this study, no correction was performed for multiple hypothesis testing.

Functional Analysis and Integrative Analysis

For relating differentially expressed miRNAs with differentially expressed mRNAs, the miRNA target filter feature (experimentally validated as well as high confidence sequence matches) is used within Ingenuity Pathway Analysis (IPA) software to find miRNAs that target the differentially expressed genes. Such integrative miRNA:mRNA analysis looks for target relationship (based on sequence affinity, using databases like TargetScan) between miRNAs and the mRNAs derived from differential expression analysis. In addition, the miRNA target filter in IPA looks for concordant changes in expression, ie, anticorrelated expression changes in miRNA:mRNA abundance.

Integrative Network Analysis of miRNA, mRNA and Protein 141 Entities via Ingenuity Pathway Network Analysis

Core analysis and functional analysis were performed on the differentially expressed gene, protein, and miRNA lists. We explored both “Direct and Indirect Interactions” using “Experimentally Observed OR High Confidence Predictions” in the IPA knowledgebase.

Results

Differential Expression Programs between Phenotypic Classes Reveal the Role for Cell Proliferation and Angiogenesis-Associated Pathways

The differential expression analysis procedure on the mRNA, protein, and miRNA data yields 326 genes, 76 miRNAs, and 8 proteins that are differentially expressed between these two phenotype classes. For completeness, the ontological analysis of these genes, miRNAs, and proteins (using the IPA tool) is presented (Table 1). The key pathways activated in this set are related to cellular development, cellular proliferation, cell death, interleukin signaling, and inflammatory response.^14,15

Table 1

IPA network results (integrating miRNA, mRNA, and protein expression).

TOP CANONICAL PATHWAYS
NAME	p-value	RATIO
LPS/IL-1 mediated inhibition of RXR function	3.74E-04	13/245 (0.053)
1L-17A signaling in gastric cells	1.14E-03	4/28 (0.143)
Tight junction signaling	2.5E-03	9/167 (0.054)
p53 signaling	2.5E-03	7/113 (0.062)
Pathogenesis of multiple sclerosis	1.17E-03	2/10 (0.2)
DISEASES AND DISORDERS
NAME	p-value	# MOLECULES
Cancer	1.79E-13–7.55E-03	311
Organismal injury and abnormalities	1.79E-13–7.55E-03	200
Reproductive system disease	1.79E-13–7.554E-03	163
Gastrointestinal disease	3.52E-11–6.97E-03	159
Dermatological diseases and conditions	3.59E-11–6.97E-03	41
PHYSIOLOGICAL SYSTEM DEVELOPMENT AND FUNCTION
NAME	p-value	# MOLECULES
Reproductive system development and function	1.80E-06–6.97E-03	25
Tissue morphology	1.80E-06–7.38E-03	80
Organismal development	6.06E-05–7.45E-03	109
Organismal development	2.09E-05–7.4E-03	65
Tumor morphology	2.21E-05–6.97E-03	48
MOLECULAR AND CELLULAR FUNCTIONS
NAME	p-value	# MOLECULES
Cellular development	1.70E-06–7.23E-03	124
Cellular growth and proliferation	1.80E-06–7.38E-03	130
Cell death and survival	6.06E-06–7.45E-03	130
Cell cycle	2.09E-05–7.45E-03	40
Cell morphology	2.21E-04–6.97E-03	77

Integrated miRNA:mRNA Analysis Reveals Multiple Molecules Underlying the Transition of this Perfusion-Associated Phenotype

We reasoned that combining the differentially expressed miRNAs and mRNAs might reveal a robust mechanism underlying the transition between the two phenotypic classes. Using the IPA tool, we combined the miRNA expression results with the mRNA expression results and found that eight differentially expressed miRNAs (miR-29b-3p, miR495–3p, miR30c/30d, miR-26a-5p, miR296-5p, miR128-3p, miR144-3p, and miR214-3p) target several differentially expressed genes (including PTEN, COL15A1, SPARC, ANPEP, CBFB, STRN, and TMED10) while exhibiting concordant expression correlation (the expression levels of the mRNAs are anticorrelated with the miRNA expression levels).

Several of these miRNAs have been shown to play a key role in GBM. miR-26a regulates PTEN and is observed to be amplified in high-grade gliomas in addition to promoting tumorigenesis.¹⁶ miR-29b has been shown to regulate invasion in GBM.¹⁷

Examination of the combined set of mRNAs, proteins, and miRNAs via IPA reveals several networks of potential interest. The primary biological processes (based on Fisher's exact tests for over-representation analysis) are LPS/IL-1-mediated signaling, IL-17A signaling, tight junction signaling, p53 signaling, pathogenesis of multiple sclerosis, etc. (presumably related to angiogenic activity [Table 1]).

Discussion

Perfusion MRI is a promising tool that is now being routinely used to interrogate tumor behavior, such as inflammation and vascularization. It is already being used for treatment monitoring (eg, assessment of changes due to treatment with antiangiogenic therapy). Since it is a much newer tool than conventional MRI, a study of the specific biological processes associated with perfusion imaging is useful to inform the radiologist about what specifically is being assessed using this modality.

Tumor blood volume (rCBV_max) is a measure of total tumor vasculature and angiogenesis. It has been shown to correlate with microvascular density¹⁸ as well as genes related to angiogenesis regulation¹⁹ and hence is a useful prognostic marker.^5,8 Our analysis reveals the role of several angiogenesis-related molecules (at the miRNA, mRNA, and protein levels) associated with this vascular parameter. This functional analysis that spans miRNA, mRNA, and protein entities produces several insights that could perhaps clarify the etiology of GBMs as measured via perfusion imaging in addition to possibly generating new hypotheses that could be validated in a laboratory. Specifically, animal models of GBM with overexpressed gene/miRNA levels might inform on the relationships of these molecules with tumor-associated angiogenic activity.

An examination of the IPA-derived pathways reveals roles for several pathways involved in tumorigenesis, tumor growth, and angiogenesis – specifically, IL-17 signaling is key to tumor survival, inflammation, and angiogenic activity.^20,21 Further IL-1 signaling has a role in tumor-mediated angiogenesis.^22,23 An examination of the molecular and cellular activities in IPA reveals a role for several tumor hallmarks, such as modulation of tumor morphology, cell proliferation, cell cycle, and apoptosis. Another interesting point is that the context of angiogenesis ties in the observations from our gene set analysis with the pathogenesis of multiple sclerosis, again suggesting a role for angiogenesis modulation.^24,25

Angiogenesis is one of the key hallmarks of cancer.²⁶ HER2 signaling and angiogenesis have been shown to be closely linked.^27,28 The role of EGFR in GBM is well documented.²⁹ Moreover, its role in modulation of invasion, angiogenesis, and metastasis has been reported; specifically, EGFR signaling contributes to angiogenesis via upregulation of MMP levels and cytokines such as VEGFA and IL8.³⁰ BRCA2 is a tumor suppressor,³¹ so it is expected that its level goes down in the more extreme rCBV phenotype (suggesting a more aggressive disease). However, its exact relationship with angiogenesis is not clearly characterized, especially in the context of glioma. Several of these proteins have been shown to play a key role in GBM. For example, JNK2 activation correlates strongly with EGFR expression as well as histological grade in GBM.³² JNK2 has been shown to control cell proliferation and apoptosis and thus been implicated in tumorigenesis.^33,34 Furthermore, JNK2 has been shown to an upstream regulator of NFAT activity, which is a critical regulator of angiogenesis.^35,36

Among the miRNAs listed in Table 2, miR495 has been shown to inhibit proliferation of GBM cells by downregulating CDK6, and thus is heavily downregulated in higher grade gliomas.³⁷ miRNA-26a regulates angiogenesis via modulation of the BMP/SMAD1 signaling.^38,39 miR-29b suppresses angiogenesis, migration, and invasion by regulating MMP2 levels; thus, it is expected that it would be down-regulated in the high rCBV_max group, suggesting an activation of these tumorigenic hallmarks. miR-30d has been shown to modulate tumor cell proliferation, angiogenesis, and metastasis.^40,41 Their specific role in glioma needs to be investigated more thoroughly.

Table 2

List of differentially expressed miRNA, mRNA, and protein entities.

ENTITY	miRNA^a	mRNA^b	PROTEIN
1	miR26a	PTEN	JNK2
2	miR29b	COL15A1	HER2
3	miR29b	SPARC	ACC1
4	miR30d	ANPEP	EGFR pY1068, EGFR pY1173
5	miR30d	CBFB	P90RSK
6	miR30d	STRN	BRCA2
7	miR30d	TMED10	Dvl3

Only the ones that have an experimentally proven relationship of sequence affinity are shown.

Conclusions and Future Work

This study reveals the utility of integrating genomic data to dissect the molecular programs underlying the vascular phenotype of the tumor, which is assessed with perfusion imaging. Our hypothesis is that instead of focusing solely on one -omic modality, combining -omic modalities is more likely to illuminate the molecular mechanisms underlying the disease processes. Thus, an integration of mRNA, miRNA, and RPPA data permits the combination of multiple layers of biological complexity: transcriptional, posttranscriptional, and signaling, to infer biological mechanisms underlying this hemodynamic characteristic of GBM.

There are several limitations to address in our study. First, this was a retrospective study performed on a publicly available patient subset, consisting of data acquired on multiple MRI systems with varying protocols. A study that evaluates the robustness of these integrative molecular programs across scanning protocols, scanner resolutions, and a larger sample size is essential to establishing their biological role. A separate dataset with uniformity of treatment regimen is the next step in validating the biological role of these genomic programs. Furthermore, incorporating molecular markers such as IDH mutation status and molecular subtype in a multivariate regression can help to deconvolve the underlying molecular programs underlying the perfusion phenotype in GBM. Also, the top ranking pathways and molecules identified via integrative analysis need to be validated through laboratory studies, possibly via in vivo mouse models of GBM. Finally, an experimental analysis of drugs targeting these signaling molecules and implicated pathways might permit the examination of novel or rational therapeutic combinations for GBM. Methodologically, this work infers differentially expressed entities from median-dichotomized patient groups from rCBV_max. However, examining associated genes/miRNAs/proteins using the original, continuous-valued rCBV_max values (via measures such as Spearman and Kendall correlation) is another possible avenue for examination. Future work will also involve the examination of computational features derived from perfusion data^42–44 and the assessment of molecular alterations (mRNA, miRNA, and protein) underlying those characteristics. This will permit the molecular examination of hemodynamic characteristics other than rCBV, thus illuminating other aspects beyond those routinely obtained in current clinical practice.

Author Contributions

Conception and design of research: AR, RJ. Data preparation, analysis and reporting: GM, GR, RJ, AR. All authors reviewed and approved of the final manuscript.

Footnotes

Acknowledgments

The authors would like to thank Dr. Laila Poisson for helpful discussions on this manuscript and Tomiko Fisher for assistance with manuscript editing.

References

Pope

W.B.

, Young

J.R.

, Ellingson

B.M.

Advances in MRI assessment of gliomas and response to anti-VEGF therapy.

Curr Neurol Neurosci Rep. 2011; 11(3): 336–44.

Gutman

D.A.

, Cooper

L.A.

, Hwang

S.N.

MR imaging predictors of molecular profile and survival: multi-institutional study of the TCGA glioblastoma data set.

Radiology. 2013; 267(2): 560–9.

Diehn

, Nardini

, Wang

D.S.

Identification of noninvasive imaging surrogates for brain tumor gene-expression modules.

Proc Natl Acad Sci U S A. 2008; 105(13): 5213–8.

Segal

, Sirlin

C.B.

, Ooi

Decoding global gene expression programs in liver cancer by noninvasive imaging.

Nat Biotechnol. 2007; 25(6): 675–80.

Jain

, Poisson

, Narang

Genomic mapping and survival prediction in glioblastoma: molecular subclassification strengthened by hemodynamic imaging biomarkers.

Radiology. 2013; 267(1): 212–20.

Zinn

P.O.

, Mahajan

, Sathyan

Radiogenomic mapping of edema/cellular invasion MRI-phenotypes in glioblastoma multiforme.

PLoS One. 2011; 6(10): e25451.

Pope

W.B.

, Mirsadraei

, Lai

Differential gene expression in glioblastoma defined by ADC histogram analysis: relationship to extracellular matrix molecules and survival.

AJNR Am J Neuroradiol. 2012; 33(6): 1059–64.

Jain

, Poisson

L.M.

, Gutman

Outcome prediction in patients with glioblastoma by using imaging, clinical, and genomic biomarkers: focus on the nonenhancing component of the tumor.

Radiology. 2014; 272(2): 484–93.

Haralick

R.M.

, Shanmugam

, Dinstein

I.H.

Textural features for image classification.

IEEE Trans Syst Man Cybern B Cybern. 1973; 3(6): 610–21.

10.

Gatrell

A.C.

, Bailey

T.C.

, Diggle

P.J.

, Rowlingson

B.S.

Spatial point pattern analysis and its application in geographical epidemiology.

Trans Inst Br Geogr. 1996; 21(1): 256–74.

11.

Tesar

, Shimizu

, Smutek

, Kobatake

, Nawano

Medical image analysis of 3D CT images based on extension of Haralick texture features.

Comput Med Imag Grap. 2008; 32(6): 513–20.

12.

, Ganeshan

, Kozarski

, Miles

K.A.

, Goh

Assessment of primary colorectal cancer heterogeneity by using whole-tumor texture analysis: contrast-enhanced CT texture as a biomarker of 5-year survival.

Radiology. 2013; 266(1): 177–84.

13.

, Lu

, Akbani

The Cancer Proteome Atlas; 2014. Available at: http://app1.bioinformatics.mdanderson.org/tcpa/_design/basic/index.html.

14.

Angelo

L.S.

, Kurzrock

Vascular endothelial growth factor and its relationship to inflammatory mediators.

Clin Cancer Res. 2007; 13(10): 2825–30.

15.

Imhof

B.A.

, Aurrand-Lions

Angiogenesis and inflammation face off.

Nat Med. 2006; 12(2): 171–2.

16.

Huse

J.T.

, Brennan

, Hambardzumyan

The PTEN-regulating microRNA miR-26a is amplified in high-grade glioma and facilitates gliomagenesis in vivo.

Genes Dev. 2009; 23(11): 1327–37.

17.

Cortez

M.A.

, Nicoloso

M.S.

, Shimizu

miR-29b and miR-125a regulate podoplanin and suppress invasion in glioblastoma.

Genes Chromosomes Cancer. 2010; 49(11): 981–90.

18.

Jain

, Gutierrez

, Narang

In vivo correlation of tumor blood volume and permeability with histologic and molecular angiogenic markers in gliomas.

AJNR Am J Neuroradiol. 2011; 32(2): 388–94.

19.

Jain

, Poisson

, Narang

Correlation of perfusion parameters with genes related to angiogenesis regulation in glioblastoma: a feasibility study.

AJNR Am J Neuroradiol. 2012; 33(7): 1343–8.

20.

Hayata

, Iwahashi

, Ojima

Inhibition of IL-17A in tumor microenvironment augments cytotoxicity of tumor-infiltrating lymphocytes in tumor- bearing mice.

PLoS One. 2013; 8(1): e53131.

21.

Numasaki

, Fukushi

, Ono

Interleukin-17 promotes angiogenesis and tumor growth.

Blood. 2003; 101(7): 2620–7.

22.

Carmi

, Voronov

, Dotan

The role of macrophage-derived IL-1 in induction and maintenance of angiogenesis.

J Immunol. 2009; 183(7): 4705–14.

23.

Voronov

, Carmi

, Apte

R.N.

The role IL-1 in tumor-mediated angiogenesis.

Front Physiol. 2014; 5: 114.

24.

Kirk

, Frank

J.A.

, Karlik

Angiogenesis in multiple sclerosis: is it good, bad or an epiphenomenon?

J Neurol Sci. 2004; 217(2): 125–30.

25.

Seabrook

T.J.

, Littlewood-Evans

, Brinkmann

, Pöllinger

, Schnell

, Hiestand

P.C.

Angiogenesis is present in experimental autoimmune encephalomyelitis and proangiogenic factors are increased in multiple sclerosis lesions.

J Neuroinflammation. 2010; 7: 95.

26.

Hanahan

, Weinberg

R.A.

Hallmarks of cancer: the next generation.

Cell. 2011; 144(5): 646–74.

27.

Alameddine

R.S.

, Otrock

Z.K.

, Awada

, Shamseddine

Crosstalk between HER2 signaling and angiogenesis in breast cancer: molecular basis, clinical applications and challenges.

Curr Opin Oncol. 2013; 25(3): 313–24.

28.

Kumar

, Yarmand-Bagheri

The role of HER2 in angiogenesis.

Semin Oncol. 2001; 28(5 suppl 16): 27–32.

29.

Hatanpaa

K.J.

, Burma

, Zhao

, Habib

A.A.

Epidermal growth factor receptor in glioma: signal transduction, neuropathology, imaging, and radioresistance.

Neoplasia. 2010; 12(9): 675–84.

30.

Box

, Peak

, Rogers

, Eccles

EGFR signaling in invasion, angiogenesis and metastasis. In: Haley

J.D.

, Gullick

W.J.

, eds. EGFR Signaling Networks in Cancer Therapy. Cancer Drug Discovery and Development. New York, NY: Humana Press; 2008: 245–64.

31.

Shahid

, Soroka

, Kong

E.H.

Structure and mechanism of action of the BRCA2 breast cancer tumor suppressor.

Nat Struct Mol Biol. 2014; 21(11): 962–8.

32.

Yeung

Y.T.

, McDonald

K.L.

, Grewal

, Munoz

Interleukins in glioblastoma pathophysiology: implications for therapy.

Br J Pharmacol. 2013; 168(3): 591–606.

33.

Ennis

B.W.

, Fultz

K.E.

, Smith

K.A.

Inhibition of tumor growth, angiogenesis, and tumor cell proliferation by a small molecule inhibitor of c-Jun N-terminal kinase.

J Pharmacol Exp Ther. 2005; 313(1): 325–32.

34.

Sabapathy

, Wagner

E.F.

JNK2: a negative regulator of cellular proliferation.

Cell Cycle. 2004; 3(12): 1520–3.

35.

Yao

, Cho

Y.Y.

, Bergen

III . Nuclear factor of activated T3 is a negative regulator of Ras-JNK1/2-AP-1 induced cell transformation. Cancer Res. 2007; 67(18): 8725–35.

36.

Zaichuk

T.A.

, Shroff

E.H.

, Emmanuel

, Filleur

, Nelius

, Volpert

O.V.

Nuclear factor of activated T cells balances angiogenesis activation and inhibition.

J Exp Med. 2004; 199(11): 1513–22.

37.

Adlakha

Y.K.

, Saini

Brain microRNAs and insights into biological functions and therapeutic potential of brain enriched miRNA-128.

Mol Cancer. 2014; 13: 33.

38.

Chai

Z.T.

, Kong

, Zhu

X.D.

MicroRNA-26a inhibits angiogenesis by down-regulating VEGFA through the PIK3C2alpha/Akt/HIF-1alpha pathway in hepatocellular carcinoma.

PLoS One. 2013; 8(10): e77957.

39.

Icli

, Wara

A.K.

, Moslehi

MicroRNA-26a regulates pathological and physiological angiogenesis by targeting BMP/SMAD1 signaling.

Circ Res. 2013; 113(11): 1231–41.

40.

Abba

, Patil

, Allgayer

MicroRNAs in the regulation of MMPs and metastasis.

Cancers. 2014; 6(2): 625–45.

41.

S.F.

, Chang

Y.W.

, Andreu-Vieyra

miR-30d, miR-181a and miR- 199a-5p cooperatively suppress the endoplasmic reticulum chaperone and signaling regulator GRP78 in cancer.

Oncogene. 2013; 32(39): 4694–701.

42.

Zacharaki

E.I.

, Morita

, Bhatt

, O'Rourke

D.M.

, Melhem

E.R.

, Davatzikos

Survival analysis of patients with high-grade gliomas based on data mining of imaging variables.

AJNR Am J Neuroradiol. 2012; 33(6): 1065–71.

43.

Ryoo

, Choi

S.H.

, Kim

J.H.

Cerebral blood volume calculated by dynamic susceptibility contrast-enhanced perfusion MR imaging: preliminary correlation study with glioblastoma genetic profiles.

PLoS One. 2013; 8(8): e71704.

44.

Lee

, Jain

, Khalil

Texture feature ratios from relative CBV maps of perfusion MRI are associated with patient survival in glioblastoma.

AJNR Am J Neuroradiol. 2015; 37(1): 37–43.

Integrative Analysis of mRNA,microRNA,and Protein Correlates of Relative Cerebral Blood Volume Values in GBM Reveals the Role for Modulators of Angiogenesis and Tumor Proliferation

Abstract

Keywords

Introduction

Materials and Methods

Ethics Statement

Patient Population

Image Features

Genomic Features

Defining Classes Based on Median rCBVMax

Differential Expression Analysis between Cutpoint-Induced Phenotype Classes

Functional Analysis and Integrative Analysis

Integrative Network Analysis of miRNA, mRNA and Protein 141 Entities via Ingenuity Pathway Network Analysis

Results

Differential Expression Programs between Phenotypic Classes Reveal the Role for Cell Proliferation and Angiogenesis-Associated Pathways

Integrated miRNA:mRNA Analysis Reveals Multiple Molecules Underlying the Transition of this Perfusion-Associated Phenotype

Discussion

Conclusions and Future Work

Author Contributions

Footnotes

Acknowledgments

References

Defining Classes Based on Median rCBV_Max