SURVIVAL PREDICTION FOR BRAIN TUMOR PATIENTS USING GENE EXPRESSION DATA

Brain tumor is one of the most aggressive types of cancer in humans, with an estimated median survival time of 12 months and only 4% of the patients surviving more than 5 years after disease diagnosis. Until recently, brain tumor prognosis has been based only on clinical information such as tumor grade and patient age, but there are reports indicating that molecular profiling of gliomas can reveal subgroups of patients with distinct survival rates. We hypothesize that coupling molecular profiling of brain tumors with clinical information might improve predictions of patient survival time and, consequently, better guide future treatment decisions. In order to evaluate this hypothesis, the general goal of this research is to build models for survival prediction of glioma patients using DNA molecular profiles (U133 Affymetrix gene expression microarrays) along with clinical information. First, a predictive Random Forest model is built for binary outcomes (i.e. short vs. long-term survival) and a small subset of genes whose expression values can be used to predict survival time is selected. Following, a new statistical methodology is developed for predicting time-to-death outcomes using Bayesian ensemble trees. Due to a large heterogeneity observed within prognostic classes obtained by the Random Forest model, prediction can be improved by relating time-to-death with gene expression profile directly. We propose a Bayesian ensemble model for survival prediction which is appropriate for high-dimensional data such as gene expression data. Our approach is based on the ensemble "sum-of-trees" model which is flexible to incorporate additive and interaction effects between genes. We specify a fully Bayesian hierarchical approach and illustrate our methodology for the CPH, Weibull, and AFT survival models. We overcome the lack of conjugacy using a latent variable formulation to model the covariate effects which decreases computation time for model fitting. Also, our proposed models provides a model-free way to select important predictive prognostic markers based on controlling false discovery rates. We compare the performance of our methods with baseline reference survival methods and apply our methodology to an unpublished data set of brain tumor survival times and gene expression data, selecting genes potentially related to the development of the disease under study. A closing discussion compares results obtained by Random Forest and Bayesian ensemble methods under the biological/clinical perspectives and highlights the statistical advantages and disadvantages of the new methodology in the context of DNA microarray data analysis.

THE CONSEQUENCES OF DISRUPTING THE MDM2-P53 BALANCE IN HEMATOPOIESIS

The bone marrow accommodates hematopoietic stem cells and progenitors. These cells provide an indispensible resource for replenishing the blood constituents throughout an organism’s life. A tissue with such a high turn-over rate mandates intact cycling checkpoint and apoptotic pathways to avoid inappropriate cell proliferation and ultimately the development of leukemias. p53, a major tumor suppressor, is a transcription factor that regulates cell cycle, and induces apoptosis and senescence. Mice inheriting a hypomorphic p53 allele in the absence of Mdm2, a p53 inhibitor, have elevated p53 cell cycle activity and die by postnatal day 13 due to hematopoietic failure. Hematopoiesis progresses normally during embryogenesis until it moves to the bone marrow in late development. Increased oxidative stress in the bone marrow compartment postnatally is the impediment for normal hematopoiesis via activation of p53. p53 in turn stimulates the generation of more reactive oxygen species and depletes bone marrow cellularity. Also, p53 exerts various defects on the hematopoietic niche by increasing mesenchymal lineage populations and their differentiation. Hematopoietic defects are rescued with antioxidants or when cells are cultured at low oxygen levels. Deletion of p16 partially rescues bone marrow cellularity and progenitors via a p53-independent pathway. Thus, although p53 is required to inhibit tumorigenesis, Mdm2 is required to control ROS-induced p53 levels for sustainable hematopoiesis and survival during homeostasis.

Interruptions in Workflow for RNs in a Level One Trauma Center

An understanding of interruptions in healthcare is important for the design, implementation, and evaluation of health information systems and for the management of clinical workflow and medical errors. The purpose of this study is to identify and classify the types of interruptions experienced by Emergency Department(ED) nurses working in a Level One Trauma Center. This was an observational field study of Registered Nurses (RNs) employed in a Level One Trauma Center using the shadowing method. Results of the study indicate that nurses were both recipients and initiators of interruptions. Telephones, pagers, and face-to-face conversations were the most common sources of interruptions. Unlike other industries, the healthcare community has not systematically studied interruptions in clinical settings to determine and weigh the necessity of the interruption against their sometimes negative results such as medical errors, decreased efficiency, and increased costs. Our study presented here is an initial step to understand the nature, causes, and effects of interruptions, thereby improving both the quality of healthcare and patient safety. We developed an ethnographic data collection technique and a data coding method for the capturing and analysis of interruptions. The interruption data we collected are systematic, comprehensive, and close to exhaustive. They confirmed the findings from earlier studies by other researchers that interruptions are frequent events in critical care and other healthcare settings. We are currently using these data to analyze the workflow dynamics of ED clinicians, to identify the bottlenecks of information flow, and to develop interventions to improve the efficiency of emergency care through the management of interruptions.

Repression of tumor necrosis factor-related apoptosis-inducing ligand (TRAIL) but not its receptors during oral cancer progression.

BACKGROUND: TRAIL plays an important role in host immunosurveillance against tumor progression, as it induces apoptosis of tumor cells but not normal cells, and thus has great therapeutic potential for cancer treatment. TRAIL binds to two cell-death-inducing (DR4 and DR5) and two decoy (DcR1, and DcR2) receptors. Here, we compare the expression levels of TRAIL and its receptors in normal oral mucosa (NOM), oral premalignancies (OPM), and primary and metastatic oral squamous cell carcinomas (OSCC) in order to characterize the changes in their expression patterns during OSCC initiation and progression. METHODS: DNA microarray, immunoblotting and immunohistochemical analyses were used to examine the expression levels of TRAIL and its receptors in oral epithelial cell lines and in archival tissues of NOM, OPM, primary and metastatic OSCC. Apoptotic rates of tumor cells and tumor-infiltrating lymphocytes (TIL) in OSCC specimens were determined by cleaved caspase 3 immunohistochemistry. RESULTS: Normal oral epithelia constitutively expressed TRAIL, but expression was progressively lost in OPM and OSCC. Reduction in DcR2 expression levels was noted frequently in OPM and OSCC compared to respective patient-matched uninvolved oral mucosa. OSCC frequently expressed DR4, DR5 and DcR1 but less frequently DcR2. Expression levels of DR4, DR5 and DcR1 receptors were not significantly altered in OPM, primary OSCC and metastatic OSCC compared to patient-matched normal oral mucosa. Expression of proapoptotic TRAIL-receptors DR4 and DR5 in OSCC seemed to depend, at least in part, on whether or not these receptors were expressed in their parental oral epithelia. High DR5 expression in primary OSCC correlated significantly with larger tumor size. There was no significant association between TRAIL-R expression and OSSC histology grade, nodal status or apoptosis rates of tumor cells and TIL. CONCLUSION: Loss of TRAIL expression is an early event during oral carcinogenesis and may be involved in dysregulation of apoptosis and contribute to the molecular carcinogenesis of OSCC. Differential expressions of TRAIL receptors in OSCC do not appear to play a crucial role in their apoptotic rate or metastatic progression.

Gene Therapy for Pediatric Cancer: State of the Art and Future Perspectives.

While modern treatments have led to a dramatic improvement in survival for pediatric malignancy, toxicities are high and a significant proportion of patients remain resistant. Gene transfer offers the prospect of highly specific therapies for childhood cancer. "Corrective" genes may be transferred to overcome the genetic abnormalities present in the precancerous cell. Alternatively, genes can be introduced to render the malignant cell sensitive to therapeutic drugs. The tumor can also be attacked by decreasing its blood supply with genes that inhibit vascular growth. Another possible approach is to modify normal tissues with genes that make them more resistant to conventional drugs and/or radiation, thereby increasing the therapeutic index. Finally, it may be possible to attack the tumor indirectly by using genes that modify the behavior of the immune system, either by making the tumor more immunogenic, or by rendering host effector cells more efficient. Several gene therapy applications have already been reported for pediatric cancer patients in preliminary Phase 1 studies. Although no major clinical success has yet been achieved, improvements in gene delivery technologies and a better understanding of mechanisms of tumor progression and immune escape have opened new perspectives for the cure of pediatric cancer by combining gene therapy with standard therapeutic available treatments.

Interruptions in workflow for RNs in a Level One Trauma Center.

An understanding of interruptions in healthcare is important for the design, implementation, and evaluation of health information systems and for the management of clinical workflow and medical errors. The purpose of this study is to identify and classify the types of interruptions experienced by ED nurses working in a Level One Trauma Center. This was an observational field study of Registered Nurses employed in a Level One Trauma Center using the shadowing method. Results of the study indicate that nurses were both recipients and initiators of interruptions. Telephone, pagers, and face-to-face conversations were the most common sources of interruptions. Unlike other industries, the outcomes caused by interruptions resulting in medical errors, decreased efficiency and increased cost have not been systematically studied in healthcare. Our study presented here is an initial step to understand the nature, causes, and effects of interruptions, and to develop interventions to manage interruptions to improve healthcare quality and patient safety. We developed an ethnographic data collection technique and a data coding method for the capturing and analysis of interruptions. The interruption data we collected are systematic, comprehensive, and close to exhaustive. They confirmed the findings from early studies by other researchers that interruptions are frequent events in critical care and other healthcare settings. We are currently using these data to analyze the workflow dynamics of ED clinicians, identify the bottlenecks of information flow, and develop interventions to improve the efficiency of emergency care through the management of interruptions.

Genetics and epidemiology of gallbladder disease in New World native peoples.

Native peoples of the New World, including Amerindians and admixed Latin Americans such as Mexican-Americans, are highly susceptible to diseases of the gallbladder. These include cholesterol cholelithiasis (gallstones) and its complications, as well as cancer of the gallbladder. Although there is clearly some necessary dietary or other environmental risk factor involved, the pattern of disease prevalence is geographically associated with the distribution of genes of aboriginal Amerindian origin, and levels of risk generally correspond to the degree of Amerindian admixture. This pattern differs from that generally associated with Westernization, which suggests a gene-environment interaction, and that within an admixed population there is a subset whose risk is underestimated when admixture is ignored. The risk that an individual of a susceptible New World genotype will undergo a cholecystectomy by age 85 can approach 40% in Mexican-American females, and their risk of gallbladder cancer can reach several percent. These are heretofore unrecognized levels of risk, especially of the latter, because previous studies have not accounted for admixture or for the loss of at-risk individuals due to cholecystectomy. A genetic susceptibility may, thus, be as "carcinogenic" in New World peoples as any known major environmental exposure; yet, while the risk has a genetic basis, its expression as gallbladder cancer is so delayed as to lead only very rarely to multiply-affected families. Estimates in this paper are derived in part from two studies of Mexican-Americans in Starr County and Laredo, Texas.

LOSS OF GPRC5A ENHANCES SURVIVAL IN NORMAL AND MALIGNANT LUNG EPITHELIAL CELLS BY ELICITING PERSISTENT STAT3 ACTIVATION INDUCED BY AUTOCRINE LIF

Signal transduction and activator of transcription 3 (Stat3) is activated by cytokines and growth factors in many cancers. Persistent activation of Stat3 plays important role in cell growth, survival, and transformation through regulating its targeted genes. Previously, we found that mice with a deletion of the G protein-coupled receptor, family C, group 5, member a (Gprc5a) gene develop lung tumors indicating that Gprc5a is a tumor suppressor. In the present study, we examined he mechanism of Gprc5a-mediated tumor suppression. We found that epithelial cells from Gprc5a knockout mouse lung (Gprc5a-/- cells) survive better in vitro in medium deprived of exogenous growth factors and form more colonies in semi-solid medium than their counterparts from wildtype mice (Gprc5a+/+ cells). The phosphorylation of tyrosine 705 on Stat3 and the expression of Stat3-regulated anti-apoptotic genes Bcl-XL, Cryab, Hapa1a, and Mcl1 were higher in the Gprc5a-/- than in Gprc5a+/+ cells. In addition, their responses to Lif were different; Stat3 activation was persistent by Lif treatment in the Gprc5a-/- cells, but was transient in the Gprc5a+/+ cells. The persistent activation of Stat3 by Lif in Gprc5a-/- cells is due to a decreased level of Socs3 protein, a negative inhibitor of the Lif-Stat3 signaling. Restoration of Socs3 inhibited the persistent Stat3 activation in Gprc5a-/- cells. Lung adenocarcinoma cells isolated from Gprc5a-/- mice also exhibited autocrine Lif-mediated Stat3 activation. Treatment of Gprc5a-/- cells isolated from normal and tumor tissue with AG490, a Stat3 signaling inhibitor, or with dominant negative Stat3(Y705F) increased starvation-induced apoptosis and inhibited anchorage-independent growth. These results suggest that persistent Stat3 activation increased the survival and transformation of Gprc5a-/- lung cells. Thus, the tumor suppressive effects of Gprc5a are mediated, at least in part, by inhibition of Stat3 signaling through regulating the stability of the Socs3 protein.

BENCHMARKING AND IMPLEMENTATION OF A NEW INDEPENDENT MONTE CARLO DOSE CALCULATION QUALITY ASSURANCE AUDIT TOOL FOR CLINICAL TRIALS

Introduction Commercial treatment planning systems employ a variety of dose calculation algorithms to plan and predict the dose distributions a patient receives during external beam radiation therapy. Traditionally, the Radiological Physics Center has relied on measurements to assure that institutions participating in the National Cancer Institute sponsored clinical trials administer radiation in doses that are clinically comparable to those of other participating institutions. To complement the effort of the RPC, an independent dose calculation tool needs to be developed that will enable a generic method to determine patient dose distributions in three dimensions and to perform retrospective analysis of radiation delivered to patients who enrolled in past clinical trials. Methods A multi-source model representing output for Varian 6 MV and 10 MV photon beams was developed and evaluated. The Monte Carlo algorithm, know as the Dose Planning Method (DPM), was used to perform the dose calculations. The dose calculations were compared to measurements made in a water phantom and in anthropomorphic phantoms. Intensity modulated radiation therapy and stereotactic body radiation therapy techniques were used with the anthropomorphic phantoms. Finally, past patient treatment plans were selected and recalculated using DPM and contrasted against a commercial dose calculation algorithm. Results The multi-source model was validated for the Varian 6 MV and 10 MV photon beams. The benchmark evaluations demonstrated the ability of the model to accurately calculate dose for the Varian 6 MV and the Varian 10 MV source models. The patient calculations proved that the model was reproducible in determining dose under similar conditions described by the benchmark tests. Conclusions The dose calculation tool that relied on a multi-source model approach and used the DPM code to calculate dose was developed, validated, and benchmarked for the Varian 6 MV and 10 MV photon beams. Several patient dose distributions were contrasted against a commercial algorithm to provide a proof of principal to use as an application in monitoring clinical trial activity.

An Innovative Phase I Trial Design Allowing for the Identification of Multiple Potential Maximum Tolerated Doses with Combination Therapy of Targeted Agents

Treatment for cancer often involves combination therapies used both in medical practice and clinical trials. Korn and Simon listed three reasons for the utility of combinations: 1) biochemical synergism, 2) differential susceptibility of tumor cells to different agents, and 3) higher achievable dose intensity by exploiting non-overlapping toxicities to the host. Even if the toxicity profile of each agent of a given combination is known, the toxicity profile of the agents used in combination must be established. Thus, caution is required when designing and evaluating trials with combination therapies. Traditional clinical design is based on the consideration of a single drug. However, a trial of drugs in combination requires a dose-selection procedure that is vastly different than that needed for a single-drug trial. When two drugs are combined in a phase I trial, an important trial objective is to determine the maximum tolerated dose (MTD). The MTD is defined as the dose level below the dose at which two of six patients experience drug-related dose-limiting toxicity (DLT). In phase I trials that combine two agents, more than one MTD generally exists, although all are rarely determined. For example, there may be an MTD that includes high doses of drug A with lower doses of drug B, another one for high doses of drug B with lower doses of drug A, and yet another for intermediate doses of both drugs administered together. With classic phase I trial designs, only one MTD is identified. Our new trial design allows identification of more than one MTD efficiently, within the context of a single protocol. The two drugs combined in our phase I trial are temsirolimus and bevacizumab. Bevacizumab is a monoclonal antibody targeting the vascular endothelial growth factor (VEGF) pathway which is fundamental for tumor growth and metastasis. One mechanism of tumor resistance to antiangiogenic therapy is upregulation of hypoxia inducible factor 1α (HIF-1α) which mediates responses to hypoxic conditions. Temsirolimus has resulted in reduced levels of HIF-1α making this an ideal combination therapy. Dr. Donald Berry developed a trial design schema for evaluating low, intermediate and high dose levels of two drugs given in combination as illustrated in a recently published paper in Biometrics entitled “A Parallel Phase I/II Clinical Trial Design for Combination Therapies.” His trial design utilized cytotoxic chemotherapy. We adapted this design schema by incorporating greater numbers of dose levels for each drug. Additional dose levels are being examined because it has been the experience of phase I trials that targeted agents, when given in combination, are often effective at dosing levels lower than the FDA-approved dose of said drugs. A total of thirteen dose levels including representative high, intermediate and low dose levels of temsirolimus with representative high, intermediate, and low dose levels of bevacizumab will be evaluated. We hypothesize that our new trial design will facilitate identification of more than one MTD, if they exist, efficiently and within the context of a single protocol. Doses gleaned from this approach could potentially allow for a more personalized approach in dose selection from among the MTDs obtained that can be based upon a patient’s specific co-morbid conditions or anticipated toxicities.

MODELING SPORADIC TUMOR FORMATION DRIVEN BY TELOMERE DYSFUNCTION IN THE GASTROINTESTINAL TRACT

Colorectal cancer is a complex disease that is thought to arise when cells accumulate mutations that allow for uncontrolled growth. There are several recognized mechanisms for generating such mutations in sporadic colon cancer; one of which is chromosomal instability (CIN). One hypothesized driver of CIN in cancer is the improper repair of dysfunctional telomeres. Telomeres comprise the linear ends of chromosomes and play a dual role in cancer. Its length is maintained by the ribonucleoprotein, telomerase, which is not a normally expressed in somatic cells and as cells divide, telomeres continuously shorten. Critically shortened telomeres are considered dysfunctional as they are recognized as sites of DNA damage and cells respond by entering into replicative senescence or apoptosis, a process that is p53-dependent and the mechanism for telomere-induced tumor suppression. Loss of this checkpoint and improper repair of dysfunctional telomeres can initiate a cycle of fusion, bridge and breakage that can lead to chromosomal changes and genomic instability, a process that can lead to transformation of normal cells to cancer cells. Mouse models of telomere dysfunction are currently based on knocking out the telomerase protein or RNA component; however, the naturally long telomeres of mice require multiple generational crosses of telomerase null mice to achieve critically short telomeres. Shelterin is a complex of six core proteins that bind to telomeres specifically. Pot1a is a highly conserved member of this complex that specifically binds to the telomeric single-stranded 3’ G-rich overhang. Previous work in our lab has shown that Pot1a is essential for chromosomal end protection as deletion of Pot1a in murine embryonic fibroblasts (MEFs) leads to open telomere ends that initiate a DNA damage response mediated by ATR, resulting in p53-dependent cellular senescence. Loss of Pot1a in the background of p53 deficiency results in increased aberrant homologous recombination at telomeres and elevated genomic instability, which allows Pot1a-/-, p53-/- MEFs to form tumors when injected into SCID mice. These phenotypes are similar to those seen in cells with critically shortened telomeres. In this work, we created a mouse model of telomere ysfunction in the gastrointestinal tract through the conditional deletion of Pot1a that recapitulates the microscopic features seen in severe telomere attrition. Combined intestinal loss of Pot1a and p53 lead to formation of invasive adenocarcinomas in the small and large intestines. The tumors formed with long latency, low multiplicity and had complex genomes due to chromosomal instability, features similar to those seen in sporadic human colorectal cancers. Taken together, we have developed a novel mouse model of intestinal tumorigenesis based on genomic instability driven by telomere dysfunction.

THE ROLE OF TYROSINE PHOSPHORYLATION IN THE FUNCTIONS OF THE TUMOR SUPPRESSOR GPRC5A

The retinoic acid inducible G protein coupled receptor family C group 5 type A (GPRC5A) is expressed preferentially in normal lung tissue but its expression is suppressed in the majority of human non-small cell lung cancer cell lines and tissues. This differential expression has led to the idea that GPRC5A is a potential tumor suppressor. This notion was supported by the finding that mice with a deletion of the Gprc5a gene develop spontaneous lung tumors. However, there are various tumor cell lines and tissue samples, including lung, that exhibit higher GPRC5A expression than normal tissues and some reports by other groups that GPRC5A transfection increased cell growth and colony formation. Obviously, GPRC5A has failed to suppress the development of the tumors and the growth of the cell lines where its expression is not suppressed. Since no mutations were detected in the coding sequence of GPRC5A in 20 NSCLC cell lines, it’s possible that GPRC5A acts as a tumor suppressor in the context of some cells but not in others. Alternatively, we raised the hypothesis that the GPRC5A protein may be inactivated by posttranslational modification(s) such as phosphorylation. It is well established that Serine/Threonine phosphorylation of G protein coupled receptors leads to their desensitization and in a few cases Tyrosine phosphorylation of GPCRs has been linked to internalization. Others reported that GPRC5A can undergo tyrosine phosphorylation in the cytoplasmic domain after treatment of normal human mammary epithelial cells (HMECs) with epidermal growth factor (EGF) or Heregulin. This suggested that GPRC5A is a substrate of EGFR. Therefore, we hypothesized that tyrosine phosphorylation of GPRC5A by activation of EGFR signaling may lead to its inactivation. To test this hypothesis, we transfected human embryo kidney (HEK) 293 cells with GPRC5A and EGFR expression vectors and confirmed that GPRC5A can be tyrosine phosphorylated after activation of EGFR by EGF. Further, we found that EGFR and GPRC5A can interact either directly or through other proteins and that inhibition of the EGFR kinase activity decreased the phosphorylation of GPRA5A and the interaction between GPRC5A and EGFR. In c-terminal of GPRC5A, There are four tyrosine residues Y317, Y320, Y347, Y350. We prepared GPRC5A mutants in which all four tyrosine residues had been replaced by phenylalanine (mutant 4F) or each individual Tyr residue was replaced by Phe and found that Y317 is the major site for EGFR mediated phosphorylation in the HEK293T cell line. We also found that EGF can induce GPRC5A internalization both in H1792 transient and stable cell lines. EGF also partially inactivates the suppressive function of GPRC5A on cell invasion activity and anchorage-independent growth ability of H1792 stable cell lines. These finding support our hypothesis that GPRC5A may be inactivated by posttranslational modification- tyrosine phosphorylation.