Positron emission tomography (PET) tumor segmentation and quantification: Development of new algorithms

Tumor functional volume (FV) and its mean activity concentration (mAC) are the quantities derived from positron emission tomography (PET). These quantities are used for estimating radiation dose for a therapy, evaluating the progression of a disease and also use it as a prognostic indicator for predicting outcome. PET images have low resolution, high noise and affected by partial volume effect (PVE). Manually segmenting each tumor is very cumbersome and very hard to reproduce. To solve the above problem I developed an algorithm, called iterative deconvolution thresholding segmentation (IDTS) algorithm; the algorithm segment the tumor, measures the FV, correct for the PVE and calculates mAC. The algorithm corrects for the PVE without the need to estimate camera's point spread function (PSF); also does not require optimizing for a specific camera. My algorithm was tested in physical phantom studies, where hollow spheres (0.5-16 ml) were used to represent tumors with a homogeneous activity distribution. It was also tested on irregular shaped tumors with a heterogeneous activity profile which were acquired using physical and simulated phantom. The physical phantom studies were performed with different signal to background ratios (SBR) and with different acquisition times (1-5 min). The algorithm was applied on ten clinical data where the results were compared with manual segmentation and fixed percentage thresholding method called T50 and T60 in which 50% and 60% of the maximum intensity respectively is used as threshold. The average error in FV and mAC calculation was 30% and -35% for 0.5 ml tumor. The average error FV and mAC calculation were ~5% for 16 ml tumor. The overall FV error was ∼10% for heterogeneous tumors in physical and simulated phantom data. The FV and mAC error for clinical image compared to manual segmentation was around -17% and 15% respectively. In summary my algorithm has potential to be applied on data acquired from different cameras as its not dependent on knowing the camera's PSF. The algorithm can also improve dose estimation and treatment planning.^

Synthesis and Biological Evaluation of Novel Folic Acid Receptor-Targeted, β-Cyclodextrin-Based Drug Complexes for Cancer Treatment

Drug targeting is an active area of research and nano-scaled drug delivery systems hold tremendous potential for the treatment of neoplasms. In this study, a novel cyclodextrin (CD)-based nanoparticle drug delivery system has been assembled and characterized for the therapy of folate receptor-positive [FR(+)] cancer. Water-soluble folic acid (FA)-conjugated CD carriers (FACDs) were successfully synthesized and their structures were confirmed by 1D/2D nuclear magnetic resonance (NMR), matrix-assisted laser desorption ionization time-of-flight mass spectrometer (MALDI-TOF-MS), high performance liquid chromatography (HPLC), Fourier transform infrared spectroscopy (FTIR), and circular dichroism. Drug complexes of adamatane (Ada) and cytotoxic doxorubicin (Dox) with FACD were readily obtained by mixed solvent precipitation. The average size of FACD-Ada-Dox was 1.5–2.5 nm. The host-guest association constant Ka was 1,639 M−1 as determined by induced circular dichroism and the hydrophilicity of the FACDs was greatly enhanced compared to unmodified CD. Cellular uptake and FR binding competitive experiments demonstrated an efficient and preferentially targeted delivery of Dox into FR-positive tumor cells and a sustained drug release profile was seen in vitro. The delivery of Dox into FR(+) cancer cells via endocytosis was observed by confocal microscopy and drug uptake of the targeted nanoparticles was 8-fold greater than that of non-targeted drug complexes. Our docking results suggest that FA, FACD and FACD-Ada-Dox could bind human hedgehog interacting protein that contains a FR domain. Mouse cardiomyocytes as well as fibroblast treated with FACD-Ada-Dox had significantly lower levels of reactive oxygen species, with increased content of glutathione and glutathione peroxidase activity, indicating a reduced potential for Dox-induced cardiotoxicity. These results indicate that the targeted drug complex possesses high drug association and sustained drug release properties with good biocompatibility and physiological stability. The novel FA-conjugated β-CD based drug complex might be promising as an anti-tumor treatment for FR(+) cancer.

The lived experience of healing touch with cancer patients

The recent use of complementary therapies by cancer patients has prompted the study of the use of Healing Touch, an energy based therapy, to learn the meaning of the experience. By using Ray's Caring Inquiry, a phenomenologic-hermeneutic process, the lived experience of receiving Healing Touch was elicited from three cancer patients. Through the interactions of the Healing Touch practitioners, the cancer patient participants, and the energy in and around them, specific themes were expressed: body-physical, emotion-feeling, mental-knowing, and spirit-essence. Further abstracting lead to the metathemes sensation and perception. Through a change in consciousness, a oneness/wholeness was experienced. The unity of meaning elicited was the Rhythm of Oneness Through Energy which is the connecting, opening, and cocreating through caring, the wholeness of each to become one through rhythms of energy. ^

Secretion of tumor necrosis factor by human leukocytes stimulated with shark cartilage

To assess the role of shark cartilage as an immune modulator, acid, salt-soluble, and phosphate-buffered saline extracts were prepared from three different commercial sources (SL, TL, FDC) of cartilage and used to stimulate human leukocytes in vitro. Duplicate leukocyte cultures were set up, each containing 50 $\mu$l of endotoxin-free extract, 200 $\mu$l of cell suspension (2.4-2.5 $\times$ 10$\sp5$ cells) and 100 $\mu$l of medium and incubated at 37$\sp\circ$C. Cultures stimulated with LPS (5 $\mu$g/ml) or medium served as the positive and negative controls, respectively. Culture supernatants were assayed for TNF$\alpha$ by ELISA. Cartilage extracts stimulated cells to release significant levels of TNF$\alpha$ (p $<$.005); the highest response was obtained with the acid extract of SL cartilage. In comparison, response to corresponding extracts of bovine cartilage was lower (p $<$.05). The stimulatory activity was reduced (85%) following proteolytic digestion, and lost when extract was heated (60$\sp\circ$C, 20 min) or treated with urea (6M), suggesting that the active component(s) is a protein. ^

Fas-FasL expression and interactions in mouse tumor cell lines: Implications for tumor immune escape

The Fas system, comprising the Fas receptor (Fas/Apo-1/CD95) and its ligand, Fas ligand (FasL), is a central mediator of programmed cell death in various physiological and pathological processes. FasL exists as transmembrane and soluble forms and induces apoptosis on crosslinking with Fas receptor. Recent evidence indicated that tumor cells exploit this system for their immunologic escape that includes the loss of Fas and the gain of FasL expression. In the present study, nine mouse tumor cell lines of diverse origin were examined immunocytochemically for the expression of Fas and FasL. Nine of nine cell lines expressed FasL, and five of nine cell lines expressed Fas. FasL expression in these tumor cell lines was demonstrated to be functional by its induction of apoptosis in Fas-sensitive target cells in coculture experiments. These results suggest that FasL may be a prevalent mediator of immune privilege in mouse malignancies, and support the recently proposed "counterattack model" for local elimination of tumor-reactive immune cells by tumor cell-derived FasL.^ Culture supernatant of four cell lines expressing FasL showed cytotoxic effect on Fas-sensitive target cells, indicating the possibility of secreted FasL in the medium. The Fas-expressing cell lines were sensitized to anti-Fas antibody cytotoxicity following treatment with IL-2 and IFN-$\gamma$, suggesting cytokine stimulation as an effective target for future immunotherapeutic strategies. ^

Evaluation of the cytotoxic and cytostatic activities of medicinal plants used by Peruvian healers against cancer-related symptoms

Twelve plants used medicinally in Callejon de Huaylas, Department of Ancash, northeastern Peru were selected and screened in vitro for cytotoxic and cytostatic activities. Traditional preparations, aqueous extracts and organic extracts (methanol:dimethyl chloride) were tested against murine leukemia P388 cells using flow cytometry. Seventy-five percent or more of the traditional and aqueous extracts were cytostatic at concentrations of 1mg/ml. For organic extracts, cytostatic activity ranged from 8.3% (at 6.25 μg/ml) to 58.3% (at 100 μg/ml). Quinchamalium procumbens, Ophryosporus chilca and Baccharis genistelloides showed strong activity. Extracts of Brachyotum rostratum, Monnina salicifolia, and Orthrosanthus chimboracensis were particularly interesting, since they were cytostatic but not cytotoxic at concentrations of 0.5 mg/ml. These Andean plants merit further analysis. The high percentage of activity found among the traditional preparations suggests that the traditional medical knowledge of Callejon de Huaylas healers deserves respect and merits further research. ^

Thioredoxin and Jab1 control estrogen- and antiestrogen-mediated progression of the cell cycle through p27 interactions

A major problem with breast cancer treatment is the prevalence of antiestrogen resistance, be it de novo or acquired after continued use. Many of the underlying mechanisms of antiestrogen resistance are not clear, although estrogen receptor-mediated actions have been identified as a pathway that is blocked by antiestrogens. Selective estrogen receptor modulators (SERMs), such as tamoxifen, are capable of producing reactive oxygen species (ROS) through metabolic activation, and these ROS, at high levels, can induce irreversible growth arrest that is similar to the growth arrest incurred by SERMs. This suggests that SERM-mediated growth arrest may also be through ROS accumulation. Breast cancer receiving long-term antiestrogen treatment appears to adapt to this increased, persistent level of ROS. This, in turn, leads to the disruption of reversible redox signaling that involves redox-sensitive phosphatases and protein kinases and transcription factors. This has downstream consequences for apoptosis, cell cycle progression, and cell metabolism. For this dissertation, we explored if altering the ROS formed by tamoxifen also alters sensitivity of the drug in resistant cells. We explored an association with a thioredoxin/Jab1/p27 pathway, and a possible role of dysregulation of thioredoxin-mediated redox regulation contributing to the development of antiestrogen resistance in breast cancer. We used standard laboratory techniques to perform proteomic assays that showed cell proliferation, protein concentrations, redox states, and protein-protein interactions. We found that increasing thioredoxin reductase levels, and thus increasing the amount of reduced thioredoxin, increased tamoxifen sensitivity in previously resistant cells, as well as altered estrogen and tamoxifen-induced ROS. We also found that decreasing levels of Jab1 protein also increased tamoxifen sensitivity, and that the downstream effects showed a decrease p27 phosphorylation in both cases. We conclude that the chronic use of tamoxifen can lead to an increase in ROS that alters cell signaling and causing cell growth in the presence of tamoxifen, and that this resistant cell growth can be reversed with an alteration to the thioredoxin/Jab1 pathway.

The epidemiology of prostate cancer among multiethnic men

The research goal was to document differences in the epidemiology of prostate cancer among multicultural men [non-Hispanic White (NHW), Hispanic (H), non-Hispanic Black (NHB)], and Black subgroups, particularly among NHB subgroups [US-born (USB) and Caribbean-born (CBB)]. Study findings will be useful in supporting further research into Black subgroups. Aim 1 explored changes over time in reported prostate cancer prevalence, by race/ethnicity and by birthplace (within the Black subgroups). Aim 2 investigated relationships between observed and latent variables. The analytical approaches included confirmatory factor analysis (CFA for measurement models) and structural equation modeling (SEM for regression models). National Center for Health Statistics, National Health Interview Survey (NHIS) data from 1999–2008 were used. The study sample included men aged 18 and older, grouped by race/ethnicity. Among the CBB group, survey respondents were limited to the English-speaking Caribbean. Prostate cancer prevalence, by race showed a higher trend among NHB men than NHW men overall, however differences over time were not significant. CBB men reported a higher proportion of prostate cancer among cancers diagnosed than USB men overall. Due to small sample sizes, stable prostate cancer prevalence trends could not be assessed over time nor could trends in the receipt of a PSA exam among NHB men when stratified by birthplace. USB and CBB men differ significantly in their screening behavior. The effect of SES on PSA screening adjusted for risk factors was statistically significant while latent variable lifestyle was not. Among risk factors, family history of cancer exhibited a consistent positive effect on PSA screening for both USB and CBB men. Among the CBB men, the number of years lived in the US did not significantly affect PSA screening behavior. When NHB men are stratified by birthplace, CBB men had a higher overall prevalence of prostate cancer diagnoses than USB men although not statistically significant. USB men were 2 to 3 times more likely to have had a PSA exam compared to CBB men, but among CBB men birthplace did not make a significant difference in screening behavior. Latent variable SES, but not lifestyle, significantly affected the likelihood of a PSA exam.

The potential role of environmental exposures and genomic signaling in development of central nervous system tumors

The etiology of central nervous system tumors (CNSTs) is mainly unknown. Aside from extremely rare genetic conditions, such as neurofibromatosis and tuberous sclerosis, the only unequivocally identified risk factor is exposure to ionizing radiation, and this explains only a very small fraction of cases. Using meta-analysis, gene networking and bioinformatics methods, this dissertation explored the hypothesis that environmental exposures produce genetic and epigenetic alterations that may be involved in the etiology of CNSTs. A meta-analysis of epidemiological studies of pesticides and pediatric brain tumors revealed a significantly increased risk of brain tumors among children whose mothers had farm-related exposures during pregnancy. A dose response was recognized when this risk estimate was compared to those for risk of brain tumors from maternal exposure to non-agricultural pesticides during pregnancy, and risk of brain tumors among children exposed to agricultural activities. Through meta-analysis of several microarray studies which compared normal tissue to astrocytomas, we were able to identify a list of 554 genes which were differentially expressed in the majority of astrocytomas. Many of these genes have in fact been implicated in development of astrocytoma, including EGFR, HIF-1α, c-Myc, WNT5A, and IDH3A. Reverse engineering of these 554 genes using Bayesian network analysis produced a gene network for each grade of astrocytoma (Grade I-IV), and ‘key genes’ within each grade were identified. Genes found to be most influential to development of the highest grade of astrocytoma, Glioblastoma multiforme (GBM) were: COL4A1, EGFR, BTF3, MPP2, RAB31, CDK4, CD99, ANXA2, TOP2A, and SERBP1. Lastly, bioinformatics analysis of environmental databases and curated published results on GBM was able to identify numerous potential pathways and geneenvironment interactions that may play key roles in astrocytoma development. Findings from this research have strong potential to advance our understanding of the etiology and susceptibility to CNSTs. Validation of our ‘key genes’ and pathways could potentially lead to useful tools for early detection and novel therapeutic options for these tumors.

4-hydroxy estradiol-induced oxidant-mediated signaling is involved in the development of breast cancer

Breast cancer is a disease associated with excess exposures to estrogens. While the mode of cancer causation is unknown, others have shown that oxidative stress induced by prolonged exposure to estrogens mediates renal, liver, endometrial and mammary tumorigenesis though the mechanism(s) underling this process is unknown. In this study, we show that 4-hydroxyl 17β-estradiol (4-OHE2), a catechol metabolite of estrogen, induces mammary tumorigenesis in a redox dependent manner. We found that the mechanism of tumorigenesis involves redox activations of nuclear respiratory factor-1 (NRF1); a transcriptions factor associated with regulation of mitochondria biogenesis and oxidative phosphorylation (OXPHOS), as well as mediation of cell survival and growth of cells during periods of oxidative stress. Key findings from our study are as follows: (i) Prolonged treatments of normal mammary epithelial cells with 4-OHE2, increased the formation of intracellular reactive oxygen species (ROS). (ii) Estrogen-induced ROS activates redox sensitive transcription factors NRF1. (iii) 4-OHE2 through activation of serine-threonine kinase and histone acetyl transferase, phosphorylates and acetylate NRF1 respectively. (iv) Redox mediated epigenetic modifications of NRF1 facilitates mammary tumorigenesis and invasive phenotypes of breast cancer cells via modulations of genes involved in proliferation, growth and metastasis of exposed cells. (v) Animal engraftment of transformed clones formed invasive tumors. (vi) Treatment of cells or tumors with biological or chemical antioxidants, as well as silencing of NRF1 expressions, prevented 4-OHE2 induced mammary tumorigenesis and invasive phenotypes of MCF-10A cells. Based on these observations, we hypothesize that 4-OHE2 induced ROS epigenetically activate NRF1 through its phosphorylation and acylation. This, in turn, through NRF1-mediated transcriptional activation of the cell cycle genes, controls 4-OHE2 induced cell transformation and tumorigenesis.^

The role of redox signaling in the molecular mechanism of tamoxifen resistance in breast cancer.

The emergence of tamoxifen or aromatase inhibitor resistance is a major problem in the treatment of breast cancer. The molecular signaling mechanism of antiestrogen resistance is not clear. Understanding the mechanisms by which resistance to these agents arise could have major clinical implications for preventing or circumventing it. Therefore, in this dissertation we have investigated the molecular mechanisms underlying antiestrogen resistance by studying the contributions of reactive oxygen species (ROS)-induced redox signaling pathways in antiestrogen resistant breast cancer cells. Our hypothesis is that the conversion of breast tumors to a tamoxifen-resistant phenotype is associated with a progressive shift towards a pro-oxidant environment of cells as a result of oxidative stress. The hypothesis of this dissertation was tested in an in vitro 2-D cell culture model employing state of the art biochemical and molecular techniques, including gene overexpression, immunoprecipitation, Western blotting, confocal imaging, ChIP, Real-Time RT-PCR, and anchorage-independent cell growth assays. We observed that tamoxifen (TAM) acts like both an oxidant and an antioxidant. Exposure of tamoxifen resistant LCC2 cell to TAM or 17 beta-estradiol (E2) induced the formation of reactive oxidant species (ROS). The formation of E2-induced ROS was inhibited by co-treatment with TAM, similar to cells pretreated with antioxidants. In LCC2 cells, treatments with either E2 or TAM were capable of inducing cell proliferation which was then inhibited by biological and chemical antioxidants. Exposure of LCC2 cells to tamoxifen resulted in a decrease in p27 expression. The LCC2 cells exposed to TAM showed an increase in p27 phosphorylation on T157 and T187. Conversely, antioxidant treatment showed an increase in p27 expression and a decrease in p27 phosphorylation on T157 and T187 in TAM exposed cells which were similar to the effects of Fulvestrant. In line with previous studies, we showed an increase in the binding of cyclin E-Cdk2 and in the level of p27 in TAM exposed cells that overexpressed biological antioxidants. Together these findings highly suggest that lowering the oxidant state of antiestrogen resistant LCC2 cells, increases LCC2 susceptibility to tamoxifen via the cyclin dependent kinase inhibitor p27.

A novel 3-D segmentation algorithm for anatomic liver and tumor volume calculations for liver cancer treatment planning

Three-Dimensional (3-D) imaging is vital in computer-assisted surgical planning including minimal invasive surgery, targeted drug delivery, and tumor resection. Selective Internal Radiation Therapy (SIRT) is a liver directed radiation therapy for the treatment of liver cancer. Accurate calculation of anatomical liver and tumor volumes are essential for the determination of the tumor to normal liver ratio and for the calculation of the dose of Y-90 microspheres that will result in high concentration of the radiation in the tumor region as compared to nearby healthy tissue. Present manual techniques for segmentation of the liver from Computed Tomography (CT) tend to be tedious and greatly dependent on the skill of the technician/doctor performing the task. ^ This dissertation presents the development and implementation of a fully integrated algorithm for 3-D liver and tumor segmentation from tri-phase CT that yield highly accurate estimations of the respective volumes of the liver and tumor(s). The algorithm as designed requires minimal human intervention without compromising the accuracy of the segmentation results. Embedded within this algorithm is an effective method for extracting blood vessels that feed the tumor(s) in order to plan effectively the appropriate treatment. ^ Segmentation of the liver led to an accuracy in excess of 95% in estimating liver volumes in 20 datasets in comparison to the manual gold standard volumes. In a similar comparison, tumor segmentation exhibited an accuracy of 86% in estimating tumor(s) volume(s). Qualitative results of the blood vessel segmentation algorithm demonstrated the effectiveness of the algorithm in extracting and rendering the vasculature structure of the liver. Results of the parallel computing process, using a single workstation, showed a 78% gain. Also, statistical analysis carried out to determine if the manual initialization has any impact on the accuracy showed user initialization independence in the results. ^ The dissertation thus provides a complete 3-D solution towards liver cancer treatment planning with the opportunity to extract, visualize and quantify the needed statistics for liver cancer treatment. Since SIRT requires highly accurate calculation of the liver and tumor volumes, this new method provides an effective and computationally efficient process required of such challenging clinical requirements.^

Multifunctional nanoparticles in cancer: In vitro characterization, in vivo distribution, and cellular response after laser-NIR dye-induced heating

A novel biocompatible and biodegradable polymer, termed poly(Glycerol malate co-dodecanedioate) (PGMD), was prepared by thermal condensation method and used for fabrication of nanoparticles (NPs). PGMD NPs were prepared using the single oil emulsion technique and loaded with an imaging/hyperthermia agent (IR820) and a chemotherapeutic agent (doxorubicin, DOX). The size of the void PGMD NPs, IR820-PGMD NPs and DOX-IR820-PGMD NPs were approximately 90 nm, 110 nm, and 125 nm respectively. An acidic environment (pH=5.0) induced higher DOX and IR820 release compared to pH=7.4. DOX release was also enhanced by exposure to laser, which increased the temperature to 42°C. Cytotoxicity of DOX-IR820-PGMD NPs was comparable in MES-SA but was higher in Dx5 cells compared to free DOX plus IR820 (p<0.05). The combination of hyperthermia (HT) and chemotherapy improved cytotoxicity in both cell lines. We also explored the cellular response after rapid, short-term and low thermal dose (laser/Dye/NP) induced-heating, and compared it to slow, long-term and high thermal dose cell incubator heating by investigating the reactive oxygen species (ROS) level, hypoxia-inducible factor-1&agr; (HIF-1&agr;) and vascular endothelial growth factor (VEGF) expression. The cytotoxicity of IR820-PGMD NPs after laser/Dye/NP HT resulted in higher cancer cell killing compared to incubator HT. ROS level, HIF-1&agr; and VEGF expression were elevated under incubator HT, while maintained at the baseline level under the laser/Dye/NP HT. In vivo mouse studies showed that NP formulation significantly improved the plasma half-life of IR820 after tail vein injection. Significant lower IR820 content was observed in kidney in DOX-IR820-PGMD NP treatment as compared to free IR820 treatment in our biodistribution studies (p<0.05). In conclusion, both IR820-PGMD NPs and DOX-IR820-PGMD NPs were successfully developed and used for both imaging and therapeutic purposes. Rapid and short-term laser/Dye/NP HT, with a low thermal dose, did not up-regulate HIF-1&agr; and VEGF expression, whereas slow and long-term incubator HT, with a high thermal dose, can enhance expression of both HIF-1&agr; and VEGF.^

Reduced Androgen Receptor Expression Accelerates the Onset of ERBB2 Induced Breast Tumors in Female Mice

Androgen receptor (AR) is commonly expressed in both the epithelium of normal mammary glands and in breast cancers. AR expression in breast cancers is independent of estrogen receptor alpha (ERα) status and is frequently associated with overexpression of the ERBB2 oncogene. AR signaling effects on breast cancer progression may depend on ERα and ERBB2 status. Up to 30% of human breast cancers are driven by overactive ERBB2 signaling and it is not clear whether AR expression affects any steps of tumor progression in this cohort of patients. To test this, we generated mammary specific Ar depleted mice (MARKO) by combining the floxed allele of Ar with the MMTV-cre transgene on an MMTV-NeuNT background and compared them to littermate MMTV-NeuNT, Arfl/+ control females. Heterozygous MARKO females displayed reduced levels of AR in mammary glands with mosaic AR expression in ductal epithelium. The loss of AR dramatically accelerated the onset of MMTV-NeuNT tumors in female MARKO mice. In this report we show that accelerated MMTV-NeuNT-dependent tumorigenesis is due specifically to the loss of AR, as hormonal levels, estrogen and progesterone receptors expression, and MMTV-NeuNT expression were similar between MARKO and control groups. MMTV-NeuNT induced tumors in both cohorts displayed distinct loss of AR in addition to ERα, PR, and the pioneer factor FOXA1. Erbb3 mRNA levels were significantly elevated in tumors in comparison to normal mammary glands. Thus the loss of AR in mouse mammary epithelium accelerates malignant transformation rather than the rate of tumorigenesis.

Reverse Engineering of Modified Genes by Bayesian Network Analysis Defines Molecular Determinants Critical to the Development of Glioblastoma

In this study we have identified key genes that are critical in development of astrocytic tumors. Meta-analysis of microarray studies which compared normal tissue to astrocytoma revealed a set of 646 differentially expressed genes in the majority of astrocytoma. Reverse engineering of these 646 genes using Bayesian network analysis produced a gene network for each grade of astrocytoma (Grade I–IV), and ‘key genes’ within each grade were identified. Genes found to be most influential to development of the highest grade of astrocytoma, Glioblastoma multiforme were: COL4A1, EGFR, BTF3, MPP2, RAB31, CDK4, CD99, ANXA2, TOP2A, and SERBP1. All of these genes were up-regulated, except MPP2 (down regulated). These 10 genes were able to predict tumor status with 96–100% confidence when using logistic regression, cross validation, and the support vector machine analysis. Markov genes interact with NFkβ, ERK, MAPK, VEGF, growth hormone and collagen to produce a network whose top biological functions are cancer, neurological disease, and cellular movement. Three of the 10 genes - EGFR, COL4A1, and CDK4, in particular, seemed to be potential ‘hubs of activity’. Modified expression of these 10 Markov Blanket genes increases lifetime risk of developing glioblastoma compared to the normal population. The glioblastoma risk estimates were dramatically increased with joint effects of 4 or more than 4 Markov Blanket genes. Joint interaction effects of 4, 5, 6, 7, 8, 9 or 10 Markov Blanket genes produced 9, 13, 20.9, 26.7, 52.8, 53.2, 78.1 or 85.9%, respectively, increase in lifetime risk of developing glioblastoma compared to normal population. In summary, it appears that modified expression of several ‘key genes’ may be required for the development of glioblastoma. Further studies are needed to validate these ‘key genes’ as useful tools for early detection and novel therapeutic options for these tumors.