Gene expression profiling of human glioblastomas for the identification of predictive factors and characterisation of molecular mechanism implicated in response to therapy and outcome

ABSTRACT Poor outcome for glioblastoma patients is largely due to resistance to chemoradiation therapy. While epigenetic inactivation of MGMT mediated DNA repair is highly predictive for benefit from the alkylating agent therapy Temozolomide, additional mechanisms for resistance associated with molecular alterations exist. Furthermore, new concepts in cancer suggest that resistance to treatment may be linked to cancer stem cells that escape therapy and act as source for tumour recurrence. We determined gene expression signatures associated with outcome in glioblastoma patients enrolled in a phase II and phase III clinical trial establishing the new combination therapy of radiation plus concomitant and adjuvant Temozolomide. Correlating stable gene clusters emerging from unsupervised analysis with survival of 42 treated patients identified a number of biological processes associated with outcome. Most prominent, a gene cluster dominated by HOX genes and comprising PROM1, was associated with resistance. PROM1 encodes CD133, a marker for a subpopulation of tumour cells enriched for glioblastoma stem- like cells. The core of this correlated HOX cluster was comprised in the top genes of a "self-renewal signature" defined in a mouse model for MLL-AF9 initiated leukaemia. The association of the HOX gene cluster with tumour resistance was confirmed in two external data sets of 146 malignant glioma As additional resistance factors we identified over-expression of the epidermal growth factor receptor gene, EGFR, while increased gene expression related to biological features of tumour host interaction, including markers for tumour vascular and cell adhesion, and innate immune response, were associated with better outcome. The "self-renewal" signature associated with resistance to the new combination chemoradiation therapy provides first clinical evidence that glioma stem like cells may implicated in resistance in a uniformly treated cohort of glioblastoma patients. This study underlines the need to target the tumour stem cell compartment, and provides some testable hypothesis for biological mechanisms relevant for malignant behaviour of glioblastoma that may be targeted in new treatment approaches. Résumé Le glioblastome, tumeur cérébrale primaire maligne la plus fréquente, est connue pour son mauvais pronostique. Des avancées chimiothérapeutiques récentes avec des agents alkylants comme le témozolomide (TMZ), ont permis une amélioration notable dans la survie de certains patients. Les bénéficiaires ont la caractéristique commune de présenter une particularité génétique, la methylation du MGMT (methylguanine methyltransferase). Néanmoins, d'autres mécanismes de résistance en fonction des aberrations moléculaires existent. Nous avons établi les profils d'expressions génétiques des patients traités par irradiation et TMZ dans des études cliniques de phase II et III. En combinant des méthodes non-supervisées et supervisées, de l'étude de la cohorte des patients traités nous avons découvert des groupes de gènes associés à la survie. Un ensemble de gènes contenant les gènes Hox semble lié au mécanisme de résistance au traitement. Récemment, les gènes Hox ont été décrits comme faisant partie d"une signature d'autorenouvellement (self-renewal) des cellules souches cancéreuses de la leucémie. L'autorenouvellement est un processus grâce auquel les cellules souches se maintiennent tout au long de la vie. Cette association à la résistance est confirmée dans deux autres études indépendantes. Un autre facteur de résistance au traitement est la surexpression du gène EGFR. D'autre part, deux groupes de gènes associés à la relation entre hôte-tumeur tels que les marqueurs des vaisseaux tumoraux et de la réponse immunitaire innée s'avèrent avoir un effet positif sur la survie des patients traités. La découverte de la signature d'autorenouvellement comme facteur de résistance à la nouvelle chimio-radiothérapie offre une preuve clinique que les cellules souches cancéreuses sont impliquées dans la résistance au traitement. If est donc logique de penser que le traitement ciblé contre des cellules souches cancéreuses va dans l'avenir permettre des thérapies anticancéreuses plus performantes.

Testing Spatial Autocorrelation in Weighted Networks : the Modes Permutation Test

In a weighted spatial network, as specified by an exchange matrix, the variances of the spatial values are inversely proportional to the size of the regions. Spatial values are no more exchangeable under independence, thus weakening the rationale for ordinary permutation and bootstrap tests of spatial autocorrelation. We propose an alternative permutation test for spatial autocorrelation, based upon exchangeable spatial modes, constructed as linear orthogonal combinations of spatial values. The coefficients obtain as eigenvectors of the standardised exchange matrix appearing in spectral clustering, and generalise to the weighted case the concept of spatial filtering for connectivity matrices. Also, two proposals aimed at transforming an acessibility matrix into a exchange matrix with with a priori fixed margins are presented. Two examples (inter-regional migratory flows and binary adjacency networks) illustrate the formalism, rooted in the theory of spectral decomposition for reversible Markov chains.

Learning ability and longevity: a symmetrical evolutionary trade-off in Drosophila.

Learning ability can be substantially improved by artificial selection in animals ranging from Drosophila to rats. Thus these species have not used their evolutionary potential with respect to learning ability, despite intuitively expected and experimentally demonstrated adaptive advantages of learning. This suggests that learning is costly, but this notion has rarely been tested. Here we report correlated responses of life-history traits to selection for improved learning in Drosophila melanogaster. Replicate populations selected for improved learning lived on average 15% shorter than the corresponding unselected control populations. They also showed a minor reduction in fecundity late in life and possibly a minor increase in dry adult mass. Selection for improved learning had no effect on egg-to-adult viability, development rate, or desiccation resistance. Because shortened longevity was the strongest correlated response to selection for improved learning, we also measured learning ability in another set of replicate populations that had been selected for extended longevity. In a classical olfactory conditioning assay, these long-lived flies showed an almost 40% reduction in learning ability early in life. This effect disappeared with age. Our results suggest a symmetrical evolutionary trade-off between learning ability and longevity in Drosophila.

Evaluating predictive densities of U.S. output growth and inflation in a large macroeconomic data set

We evaluate conditional predictive densities for U.S. output growth and inflationusing a number of commonly used forecasting models that rely on a large number ofmacroeconomic predictors. More specifically, we evaluate how well conditional predictive densities based on the commonly used normality assumption fit actual realizationsout-of-sample. Our focus on predictive densities acknowledges the possibility that, although some predictors can improve or deteriorate point forecasts, they might have theopposite effect on higher moments. We find that normality is rejected for most modelsin some dimension according to at least one of the tests we use. Interestingly, however,combinations of predictive densities appear to be correctly approximated by a normaldensity: the simple, equal average when predicting output growth and Bayesian modelaverage when predicting inflation.

LIF and sIL-2R plasma concentrations in IVF patients on the day of embryo transfer: predictive markers of IVF outcome.

Successful implantation is still the limiting step in IVF. We hypothesized that maternal plasma concentrations of certain cytokines at the time of embryo transfer could predict the likelihood of successful implantation and pregnancy. sIL-2R, IL-6, LIF, and MMP2 concentrations were measured in plasma from 160 IVF patients (natural and stimulated IVF cycles) on the morning of the embryo transfer (ET0) and 14days later (ET+14). Patients were ultimately subdivided into four groups depending on the IVF treatment outcome (pregnancy failure, biochemical pregnancy, first-trimester miscarriage and normal term delivery). In natural and stimulated IVF cycles at ET0, sIL-2R concentrations were threefold higher in biochemical pregnancies than in pregnancy failures (P=0.020), and in natural cycles only, 2.5-fold higher in normal term deliveries than in pregnancy failures (P=0.023). Conversely, in natural and stimulated IVF cycles at ET0, LIF concentrations were one third lower in biochemical pregnancies/first-trimester miscarriages compared with pregnancy failures (P=0.042). We suggest that high sIL-2R and low LIF concentrations in maternal plasma on the morning of the embryo transfer might be associated with increased risks of early pregnancy loss, while a basal level of sIL-2R is necessary for normal term delivery outcome. Both cytokine measurements might therefore be useful in the management of IVF patients, and modulation of their concentrations could be investigated as a therapeutic alternative for women with abnormal concentrations at the time of embryo transfer.

Testing I(1) against I(d) alternatives in the presence of deteministic components

This paper discusses the role of deterministic components in the DGP and in the auxiliaryregression model which underlies the implementation of the Fractional Dickey-Fuller (FDF) test for I(1) against I(d) processes with d [0, 1). This is an important test in many economic applications because I(d) processess with d < 1 are mean-reverting although, when 0.5 = d < 1, like I(1) processes, they are nonstationary. We show how simple is the implementation of the FDF in these situations, and argue that it has better properties than LM tests. A simple testing strategy entailing only asymptotically normally distributedtests is also proposed. Finally, an empirical application is provided where the FDF test allowing for deterministic components is used to test for long-memory in the per capita GDP of several OECD countries, an issue that has important consequences to discriminate between growth theories, and on which there is some controversy.

Bringing game theory to hypothesis testing: Establishing finite sample bounds on inference

Small sample properties are of fundamental interest when only limited data is avail-able. Exact inference is limited by constraints imposed by speci.c nonrandomizedtests and of course also by lack of more data. These e¤ects can be separated as we propose to evaluate a test by comparing its type II error to the minimal type II error among all tests for the given sample. Game theory is used to establish this minimal type II error, the associated randomized test is characterized as part of a Nash equilibrium of a .ctitious game against nature.We use this method to investigate sequential tests for the di¤erence between twomeans when outcomes are constrained to belong to a given bounded set. Tests ofinequality and of noninferiority are included. We .nd that inference in terms oftype II error based on a balanced sample cannot be improved by sequential sampling or even by observing counter factual evidence providing there is a reasonable gap between the hypotheses.

MGMT promoter methylation is prognostic but not predictive for outcome to adjuvant PCV chemotherapy in anaplastic oligodendroglial tumors: a report from EORTC Brain Tumor Group Study 26951.

PURPOSE: O6-methylguanine-methyltransferase (MGMT) promoter methylation has been shown to predict survival of patients with glioblastomas if temozolomide is added to radiotherapy (RT). It is unknown if MGMT promoter methylation is also predictive to outcome to RT followed by adjuvant procarbazine, lomustine, and vincristine (PCV) chemotherapy in patients with anaplastic oligodendroglial tumors (AOT). PATIENTS AND METHODS: In the European Organisation for the Research and Treatment of Cancer study 26951, 368 patients with AOT were randomly assigned to either RT alone or to RT followed by adjuvant PCV. From 165 patients of this study, formalin-fixed, paraffin-embedded tumor tissue was available for MGMT promoter methylation analysis. This was investigated with methylation specific multiplex ligation-dependent probe amplification. RESULTS: In 152 cases, an MGMT result was obtained, in 121 (80%) cases MGMT promoter methylation was observed. Methylation strongly correlated with combined loss of chromosome 1p and 19q loss (P = .00043). In multivariate analysis, MGMT promoter methylation, 1p/19q codeletion, tumor necrosis, and extent of resection were independent prognostic factors. The prognostic significance of MGMT promoter methylation was equally strong in the RT arm and the RT/PCV arm for both progression-free survival and overall survival. In tumors diagnosed at central pathology review as glioblastoma, no prognostic effect of MGMT promoter methylation was observed. CONCLUSION: In this study, on patients with AOT MGMT promoter methylation was of prognostic significance and did not have predictive significance for outcome to adjuvant PCV chemotherapy. The biologic effect of MGMT promoter methylation or pathogenetic features associated with MGMT promoter methylation may be different for AOT compared with glioblastoma.

Testing demographic models of effective population size.

Due to practical difficulties in obtaining direct genetic estimates of effective sizes, conservation biologists have to rely on so-called 'demographic models' which combine life-history and mating-system parameters with F-statistics in order to produce indirect estimates of effective sizes. However, for the same practical reasons that prevent direct genetic estimates, the accuracy of demographic models is difficult to evaluate. Here we use individual-based, genetically explicit computer simulations in order to investigate the accuracy of two such demographic models aimed at investigating the hierarchical structure of populations. We show that, by and large, these models provide good estimates under a wide range of mating systems and dispersal patterns. However, one of the models should be avoided whenever the focal species' breeding system approaches monogamy with no sex bias in dispersal or when a substructure within social groups is suspected because effective sizes may then be strongly overestimated. The timing during the life cycle at which F-statistics are evaluated is also of crucial importance and attention should be paid to it when designing field sampling since different demographic models assume different timings. Our study shows that individual-based, genetically explicit models provide a promising way of evaluating the accuracy of demographic models of effective size and delineate their field of applicability.

Learning of courtship components in Drosophila mercatorum (Paterson & Wheller) (Diptera, Drosophilidae)

In Drosophila, courtship is an elaborate sequence of behavioural patterns that enables the flies to identify conspecific mates from those of closely related species. This is important because drosophilids usually gather in feeding sites, where males of various species court females vigorously. We investigated the effects of previous experience on D. mercatorum courtship, by testing if virgin males learn to improve their courtship by observing other flies (social learning), or by adjusting their pre-existent behaviour based on previous experiences (facilitation). Behaviours recorded in a controlled environment were courtship latency, courtship (orientation, tapping and wing vibration), mating and other behaviours not related to sexual activities. This study demonstrated that males of D. mercatorum were capable of improving their mating ability based on prior experiences, but they had no social learning on the development of courtship.

Testing the descriptive performance of the rank-dependent utility in the domain of health profiles

Expected utility theory (EUT) has been challenged as a descriptive theoryin many contexts. The medical decision analysis context is not an exception.Several researchers have suggested that rank dependent utility theory (RDUT)may accurately describe how people evaluate alternative medical treatments.Recent research in this domain has addressed a relevant feature of RDU models-probability weighting-but to date no direct test of this theoryhas been made. This paper provides a test of the main axiomatic differencebetween EUT and RDUT when health profiles are used as outcomes of riskytreatments. Overall, EU best described the data. However, evidence on theediting and cancellation operation hypothesized in Prospect Theory andCumulative Prospect Theory was apparent in our study. we found that RDUoutperformed EU in the presentation of the risky treatment pairs in whichthe common outcome was not obvious. The influence of framing effects onthe performance of RDU and their importance as a topic for future researchis discussed.

Exact and approximate stepdown methods for multiple hypothesis testing

Consider the problem of testing k hypotheses simultaneously. In this paper,we discuss finite and large sample theory of stepdown methods that providecontrol of the familywise error rate (FWE). In order to improve upon theBonferroni method or Holm's (1979) stepdown method, Westfall and Young(1993) make eective use of resampling to construct stepdown methods thatimplicitly estimate the dependence structure of the test statistics. However,their methods depend on an assumption called subset pivotality. The goalof this paper is to construct general stepdown methods that do not requiresuch an assumption. In order to accomplish this, we take a close look atwhat makes stepdown procedures work, and a key component is a monotonicityrequirement of critical values. By imposing such monotonicity on estimatedcritical values (which is not an assumption on the model but an assumptionon the method), it is demonstrated that the problem of constructing a validmultiple test procedure which controls the FWE can be reduced to the problemof contructing a single test which controls the usual probability of a Type 1error. This reduction allows us to draw upon an enormous resamplingliterature as a general means of test contruction.

Combining multivariate density forecasts using predictive criteria

This paper combines multivariate density forecasts of output growth, inflationand interest rates from a suite of models. An out-of-sample weighting scheme based onthe predictive likelihood as proposed by Eklund and Karlsson (2005) and Andersson andKarlsson (2007) is used to combine the models. Three classes of models are considered: aBayesian vector autoregression (BVAR), a factor-augmented vector autoregression (FAVAR)and a medium-scale dynamic stochastic general equilibrium (DSGE) model. Using Australiandata, we find that, at short forecast horizons, the Bayesian VAR model is assignedthe most weight, while at intermediate and longer horizons the factor model is preferred.The DSGE model is assigned little weight at all horizons, a result that can be attributedto the DSGE model producing density forecasts that are very wide when compared withthe actual distribution of observations. While a density forecast evaluation exercise revealslittle formal evidence that the optimally combined densities are superior to those from thebest-performing individual model, or a simple equal-weighting scheme, this may be a resultof the short sample available.

Testing financing constraints on firm investment using variable capital

We consider a dynamic multifactor model of investment with financing imperfections,adjustment costs and fixed and variable capital. We use the model to derive a test offinancing constraints based on a reduced form variable capital equation. Simulation resultsshow that this test correctly identifies financially constrained firms even when the estimationof firms investment opportunities is very noisy. In addition, the test is well specified inthe presence of both concave and convex adjustment costs of fixed capital. We confirmempirically the validity of this test on a sample of small Italian manufacturing companies.

The organization of research corporations and researcher ability

This paper analyzes the formation of Research Corporations as an alternative governance structure for performing R&D compared to pursuing in-house R&D projects. Research Corporations are privatefor-profit research centers that bring together several firms with similar research goals. In a Research Corporation formal authority over the choice of projects is jointly exercised by the top management of the member firms. A private for-profit organization cannot commit not to interfere with the project choice of the researchers. However, increasing the number of member firms of the Research Corporation reduces the incentive of member firms to meddle with the research projects of researchers because exercising formal authority over the choice of research projects is a public good. The Research Corporation thus offers researchers greater autonomy than a single firm pursuing an identical research program in its in-house R&D department. This attracts higher ability researchers to the Research Corporation compared to the internal R&D department. The paper uses the theoretical model to analyze the organization of the Microelectronics and Computer Technology Corporation (MCC). The facts of this case confirm the existence of a tension between control over the choice of research projects and the ability of researchers that the organization is able to attract or hold onto.