WO2007069882A1 - Genetic brain tumor markers - Google Patents
Genetic brain tumor markers Download PDFInfo
- Publication number
- WO2007069882A1 WO2007069882A1 PCT/NL2005/000855 NL2005000855W WO2007069882A1 WO 2007069882 A1 WO2007069882 A1 WO 2007069882A1 NL 2005000855 W NL2005000855 W NL 2005000855W WO 2007069882 A1 WO2007069882 A1 WO 2007069882A1
- Authority
- WO
- WIPO (PCT)
- Prior art keywords
- genes
- oligodendroglial
- tumor
- clustering
- subject
- Prior art date
Links
Classifications
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q1/00—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions
- C12Q1/68—Measuring or testing processes involving enzymes, nucleic acids or microorganisms; Compositions therefor; Processes of preparing such compositions involving nucleic acids
- C12Q1/6876—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes
- C12Q1/6883—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for diseases caused by alterations of genetic material
- C12Q1/6886—Nucleic acid products used in the analysis of nucleic acids, e.g. primers or probes for diseases caused by alterations of genetic material for cancer
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B20/00—ICT specially adapted for functional genomics or proteomics, e.g. genotype-phenotype associations
- G16B20/20—Allele or variant detection, e.g. single nucleotide polymorphism [SNP] detection
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B25/00—ICT specially adapted for hybridisation; ICT specially adapted for gene or protein expression
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B25/00—ICT specially adapted for hybridisation; ICT specially adapted for gene or protein expression
- G16B25/10—Gene or protein expression profiling; Expression-ratio estimation or normalisation
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q2600/00—Oligonucleotides characterized by their use
- C12Q2600/106—Pharmacogenomics, i.e. genetic variability in individual responses to drugs and drug metabolism
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q2600/00—Oligonucleotides characterized by their use
- C12Q2600/118—Prognosis of disease development
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q2600/00—Oligonucleotides characterized by their use
- C12Q2600/158—Expression markers
-
- C—CHEMISTRY; METALLURGY
- C12—BIOCHEMISTRY; BEER; SPIRITS; WINE; VINEGAR; MICROBIOLOGY; ENZYMOLOGY; MUTATION OR GENETIC ENGINEERING
- C12Q—MEASURING OR TESTING PROCESSES INVOLVING ENZYMES, NUCLEIC ACIDS OR MICROORGANISMS; COMPOSITIONS OR TEST PAPERS THEREFOR; PROCESSES OF PREPARING SUCH COMPOSITIONS; CONDITION-RESPONSIVE CONTROL IN MICROBIOLOGICAL OR ENZYMOLOGICAL PROCESSES
- C12Q2600/00—Oligonucleotides characterized by their use
- C12Q2600/172—Haplotypes
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B20/00—ICT specially adapted for functional genomics or proteomics, e.g. genotype-phenotype associations
Definitions
- the invention relates to the field of diagnosis of tumors, especially brain tumors, more especially oligodendroglial tumors, more particular to the prediction of susceptibility to treatment for patients with brain tumor.
- Diffuse gliomas are the most common primary central nervous system tumors in adults (Legler, J.M. et al., (1999) J. Natl. Cancer Inst. 91; 1382-1390; Macdonald, D.R. (2003) Semin. Oncol. 30: 72-76) and it is estimated that approximately 18,000 new patients per annum are diagnosed with a primary brain tumor in the USA (CBTRUS 2004-2005 statistical report). The worldwide standard for grading and classification of these tumors is at present the WHO classification (Kleihues, P and Cavenee, W.K., World Health Organization Classification of Tumours of the Nervous System, Lyon: WHO/IARC, 2000).
- gliomas Based on their histological appearance gliomas can be divided into astrocytic tumors, pure oligodendroglial tumors and mixed oligoastrocytic tumors. The latter two are grouped together as oligodendroglial tumors.
- the oligodendrogliomas comprise approximately 20% of all gliomas, and in comparison to most other gliomas, have a relatively long average survival time (5-12 years) after diagnosis (Okamoto, Y. et al., (2004) Acta Neuropathol. 28:28; Johannesen, T.B. et al. (2003) J. Neurosurg. 99: 854-862).
- oligodendroglial tumors are their sensitivity to therapy, especially radiotherapy and chemotherapy.
- the majority of oligodendroglial tumors respond favourably to chemotherapy with alkylating agents (either temolozomide or PCV, a combination therapy of procarbazine, CCNU, and vincristine), whereas other gliomas are often chemoresistant (Van den Bent, M.J. et al. (1998) Neurology 51: 1140-1145; Van den Bent, M.J. et al. (2003) J. Clin. Oncol. 21: 2525-2528).
- alkylating agents either temolozomide or PCV, a combination therapy of procarbazine, CCNU, and vincristine
- oligodendral tumors renders it therefore important to correctly identify this subtype of gliomas.
- histological classification and grading of gliomas has a significant subjective component.
- malignant gliomas can also be classified according to their gene expression profile (Nutt, CL. et al. (2003) Cancer Res. 63: 1602-1607).
- a common genomic aberration is a combined loss of the short arm of chromosome 1 (Ip) and the long arm of chromosome 19 (19q) (Okamoto, Y et al.,2004; Cairncross J.G. et al., (1998) J. Natl. Cancer Inst.90: 1473- 1479; Kros J.M. et al., (1999) J. Pathol. 188:282-288; Smith J.S. et al., (1999) Oncogene 18:4144-4152; Thiessen B.
- Loss of heterozygosity (LOH) on both chromosomal arms is correlated with a favourable response to therapy: A response to treatment is observed in 80-90% of oligodendroglial tumors with Ip LOH and in 25-30% without Ip LOH (Cairncross, J.G. et al, 1998; Thiessen, B. et al, 2003; van den Bent, M.J. et al., 2003).
- Other chromosomal aberrations observed at lower frequency include LOH on 1Oq and amplification of 7pll (Kitange G. et al. (2004) Genes
- Chromosomes Cancer Chromosomes Cancer. These aberrations are correlated with poor prognosis and are negatively correlated with LOH on Ip and 19q. This correlation between response to treatment and chromosomal aberrations can therefore help identify chemosensitive oligodendroglial tumors. However, predicting the tumors' response to treatment by its chromosomal status also incorrectly classifies a significant percentage of tumors.
- Expression profiling can be an alternative approach to identify oligodendroglial tumors that will benefit from therapeutic treatment. Although expression profiling has been performed on oligodendroglial tumors, mRNA expression has thus far not been correlated to treatment response.
- the current inventors have now surprisingly shown that gene expression can be used to be correlated with susceptibility to treatment and increased survival,independant of the (Ip and 19q) chromosomal status of the tumor. Further, also correlations have been found between gene expression and loss of Ip and 19q.
- the invention now comprise a method for producing a classification scheme for oligodendroglial tumors comprising the steps of: a) providing a plurality of reference samples, said reference samples comprising cell samples from a plurality of reference subjects suffering from oligodendroglial tumors, with known responsiveness to therapy and survival; b) providing reference profiles by establishing a gene expression profile, matched with parameters for sensitivity to treatment and survival for each of said reference samples individually; c) clustering said individual reference profiles according to a statistical procedure, comprising:
- the clustering of said gene expression profiles is performed based on the information of differentially-expressed genes and the sensitivity to treatment and/or survival of the subject, wherein, preferably, the clustering of said gene expression profiles with respect to treatment response is performed based on the information of the genes of Table 3, whereas the clustering of said gene expression profiles with respect to survival is performed based on the information of the genes of Table 4.
- Another embodiment of the invention is a method for classifying an oligodendroglial tumor of a subject suffering from an glial tumor, comprising the steps of: a) providing a classification scheme for oligodendroglial tumors according to the above described method; b) providing a subject profile by establishing a gene expression profile for said subject; c) clustering the subject profile together with reference profiles; d) determining in said scheme the clustered position of said subject profile among the reference profiles, and e) assigning to said glial tumor the oligodendroglial tumor class that corresponds to said clustered position.
- the gene expression profile with respect to treatment response comprises the expression parameters of a set of genes according to table 3, still more preferably 1 to 50 genes of the genes of table 3, whereas the gene expression profile with respect to survival comprises the expression parameters of a set of genes according to Table 4, more preferably 1 tot 50 genes of the genes of Table 4.
- a further embodiment of the invention is a method of determining the prognosis for a subject suffering from an oligodendroglial tumor, said method comprising the steps of: a) providing a classification scheme for oligodendroglial tumors by producing such a scheme according to the above described method; b) determining the prognosis for each olidendroglial tumor class in said scheme based on clinical records for the subjects comprised in said class; c) establishing the oligodendroglial class of a subject suffering from an oligodendroglial tumor by classifying the oligodendroglial tumor in said subject according to a method according to the invention, and d) assigning to said subject the prognosis corresponding to the established oligodendroglial tumor class of said subject.
- the invention provides for a method of determining the prognosis for a subject suffering from an oligodendroglial tumor, said method comprising the steps of: a) isolation of RNA from tumor cells of said subject; b) preparation of antisense, biotinylated RNA to the RNA of step a); c) hybridisation of said antisense, biotinylated DNA on Affymetrix U133A or U133 Plus2.0 GeneChips®; d) normalising the measured values for the gene set of Table 3; e) clustering the obtained data together with reference data, obtained from a reference set of patients with known prognoses; and f) determining the prognosis on basis of the subgroup/cluster to which the data of the subject are clustering.
- the invention provides for an oligonucleotide microarray of maximal 500 probesets, comprising at least 1, preferably at least 2, more preferably at least 25, still more preferably at least 100 oligonucleotide probes which each are capable of hybridizing under stringent conditions to different genes of the oligodendroglial tumor-associated genes selected from Table 3.
- the invention provides for an oligonucleotide microarray of maximal 500 probesets, comprising at least 1, preferably at least 2, more preferably at least 25, still more preferably at least 100 oligonucleotide probes which each are capable of hybridizing under stringent conditions to different genes of the oligodendroglial tumor-associated genes selected from Table 4.
- the invention provides for a method using an oligonucleotide microarray, which can be used for the determination of the presence of Ip LOH, 19q LOH or lp/19q LOH.
- the microarray for these determination should comprise the genesets of Table 5, 6 and 7, respectively.
- the invention also comprises an oligonucleotide microarray of maximal 500 probesets, comprising at least 1, preferably at least 2, more preferably at least 25, still more preferably at least 50 oligonuclaotide probes which each are capable of hybridizing under stringent conditions to different genes of the oligodendroglial tumor-associated genes selected form Table 5, 6 and 7, respectively.
- the invention also comprises a kit- of-parts comprising an oligonucleotide microarray as described above and means for comparing a gene expression profile determined by using said microarray with a database of oligodendroglial tumor reference expression profiles.
- EGFR ampl amplification of the EFGR chromosomal locus:
- Response response to therapy IS complete response, I — I partial response, I — I stable disease, I l progressive disease.
- Sur ⁇ tot survival (years) from time of diagnosis, jjj >10, Q 7-10, Q 3-7, Q] ⁇ 3.A: patient alive at time of analysis.
- PCA Principle components analysis
- the Ip and 19q status are color coded with
- BH no LOH on Ip and 19q
- LOH on Ip and 19q
- FIG. 3 PCA and hierarchical clustering based on 16 probesets differentially expressed between chemosensitive (CR+PR (complete response, partial response)) and chemoresistant (SD+PD, stable disease, progressive disease)) oligodendroglial tumors.
- B Hierarchical clustering based on 16 differentially expressed probesets. Relative expression levels of individual genes (columns) are plotted against individual tumor samples (rows). Gene expression levels are color coded with red and green indicating high (+1.8) and low green (-1.8) expression respectively.
- Dendrograms denote hierarchical clustering of samples (top) and genes (left) using Wards method. Hierarchical clustering separates tumors that fully respond to chemotherapy (CR) from tumors that do not respond (SD+PD). Furthermore, hierarchical clustering also clearly separates tumors with poor prognosis (subgroup III in figure 1) from other oligodendroglial tumors. Responses in oligodendroglial tumors are color coded with [33 complete response, ⁇ ZH partial response, ⁇ Hl stable disease,
- PCA hierarchical clustering based on 103 probesets associated with survival after diagnosis.
- A PCA identifies three main clusters of samples: oligodendroglial tumors with short survival, oligodendroglial tumors with long survival and control samples. Two low-grade samples (38 and 42, survival ⁇ 10 years I I ) cluster between control and tumor samples. PCA analysis separates short vs. long survivors on the first principle component axis (PCAl) whereas control and tumor samples are separated by the second PCA axis.
- PCAl first principle component axis
- B Hierarchical clustering based on 103 differentially expressed probesets. Relative expression levels of individual genes (columns) are plotted against individual tumor samples (rows).
- Gene expression levels are color coded with red and green indicating high (+2) and low green (-2) expression respectively.
- Dendrograms denote hierarchical clustering (Euclidian distance) of samples (top) and genes (left). Interestingly, the subgroups identified by hierarchical clustering are virtually identical to the subgroups that were identified by unsupervised clustering (figure 1). Survival after diagnosis is depicted astasia >10 years survival.) 1 ⁇ 10 years survival, 1 I ⁇ 7 years survival, I 1 ⁇ 4 years survival,
- the current inventors performed expression profiling on oligodendroglial tumors and correlated the results to response to treatment, survival after diagnosis and common chromosomal aberrations.
- One of the findings was that the chromosomal aberrations led to ⁇ 50% expression of some but not all of the genes that had been deleted, Thus, this means that it is not straightforward to use the expression data of the genes from the Ip and 19q loci for the determination of the presence of a loss of heterozygosity (LOH) in these areas.
- LHO loss of heterozygosity
- the present inventors have found that a subset of genes, which show a reduced expression when one of the chromosomal arms Ip and 19q are deleted can be used to detect these chromosomal aberrations.
- classifying is used in its art-recognized meaning and thus refers to arranging or ordering items, i.e. gene expression profiles, by classes or categories or dividing them into logically hierarchical classes, subclasses, and sub-subclasses based on the characteristics they have in common and/or that distinguish them.
- classifying refers to assigning, to a class or kind, an unclassified item.
- a "class” then being a grouping of items, based on one or more characteristics, attributes, properties, qualities, effects, parameters, etc., which they have in common, for the purpose of classifying them according to an established system or scheme.
- classification scheme is used in its art-recognized meaning and thus refers to a list of classes arranged according to a set of pre- established principles, for the purpose of organizing items in a collection or into groups based on their similarities and differences.
- clustering refers to the activity of collecting, assembling and/or uniting into a cluster or clusters items with the same or similar elements, a “cluster” referring to a group or number of the same or similar items, i.e. gene expression profiles, gathered or occurring closely together based on similarity of characteristics. "Clustered” indicates an item has been subjected to clustering.
- clustered position refers to the location of an individual item, i.e. a gene expression profile, in amongst a number of clusters, said location being determined by clustering said item with at least a number of items from known clusters.
- the process of clustering used in a method of the present invention may be any mathematical process known to compare items for similarity in characteristics, attributes, properties, qualities, effects, parameters, etc..
- Statistical analysis such as for instance multivariance analysis, or other methods of analysis may be used.
- methods of analysis such as self-organising maps, hierarchical clustering, multidimensional scaling, principle component analysis, supervised learning, k-nearest neighbours, support vector machines, discriminant analysis, partial least square methods and/or Pearson's correlation coefficient analysis are used.
- Pearson's correlation coefficient analysis significance analysis of microarrays (SAM) and/or prediction analysis of microarrays (PAM) are used to cluster gene expression profiles according to similarity.
- a highly preferred method of clustering comprises similarity clustering of gene expression profiles wherein the expression level of differentially-expressed genes, having markedly lower or higher expression than the geometric mean expression level determined for all genes in all profiles to be clustered, is log(2) transformed, and wherein the transformed expression levels of all differentially-expressed genes in all profiles to be clustered is clustered by using K-means.
- a numerical query may then be used to select a subset of genes used in the process of hierarchical clustering (Eisen et al., 1998), thus, numerical queries may be run to select differentially expressed genes relative to the calculated geometric mean to select a smaller group of genes for hierarchical clustering.
- Unsupervised sample clustering using genes obtained by numerical or threshold filtering is used to identify discrete clusters of samples as well as the gene-signatures associated with these clusters.
- gene signatures is used herein to refer to the set of genes that define the discrete position of the cluster apart from all other clusters, and includes cluster-specific genes.
- a numerical or threshold filtering is used to select genes for the analysis that are most likely of diagnostic relevance.
- Hierarchical clustering allows for visualization of large variation in gene expression across samples or present in most samples, and these genes could be used for unsupervised clustering so that clustering results are not affected by the noise from absent or non- changed genes.
- K-means clustering may be performed on all genes, the Pearson correlation is preferably calculated based on a subset of genes.
- the larger the threshold for accepting a deviation or change from the geometric mean the smaller the number of genes that is selected by this filtering procedure.
- Different cut-off or threshold values were used to prepare lists with different numbers of genes. The higher the number of genes selected and included on such lists, the more noise is generally encountered within the dataset, because there will be a relatively large contribution of non-tumor pathway related genes in such lists.
- the filtering and selection procedure is preferably optimized such that the analysis is performed on as many genes as possible, while minimizing the noise.
- All genes with changed expression values in at least one sample higher than or equal to 1.5 times the log(2) transformed expression values and genes with changed expression values lower than or equal to -1.5 times the log(2) transformed expression value means are selected for unsupervised clustering.
- the subset of genes showing a markedly higher or lower expression than the geometric mean may for instance be a value that is more than 1.5 times the geometric mean value, preferably more than 2 times the geometric mean value, Likewise, a markedly lower expression than the geometric mean expression level may for instance be a value that is less than 0.8 times the geometric mean value, preferably less than 0.6 times the geometric mean value.
- the present invention now provides several methods to accurately identify known as well as newly discovered diagnostically, prognostically and therapeutically relevant subgroups of oligodendroglial tumors, as well as methods that can predict if treatment is likely to be effective.
- the basis of these methods resides in the measurement of (oligodendroglial tumor-specific) gene expression in subjects suffering from brain tumors.
- the methods and compositions of the invention thus provide tools useful in choosing a therapy for brain tumor patients, including methods for assigning an brain tumor patient to a brain tumor class or cluster, methods of choosing a therapy for a brain tumor patient, and methods of determining the survival prognosis for a brain tumor patient.
- the methods of the invention comprise in various aspects the steps of establishing a gene expression profile of subject samples, for instance of reference subjects suffering from a brain tumor or of a subject diagnosed or classified as having a brain tumor.
- the expression profiles of the present invention are generated from samples from subjects having a brain tumor.
- the samples from the subject used to generate the expression profiles of the present invention can be derived from a tumor biopsy, wherein the sample comprises preferably more than 75% tumor cells.
- Gene expression profiling or “expression profiling” is used herein in its art-recognised meaning and refers to a method for measuring the transcriptional state (mRNA) or the translational state (protein) of a plurality of genes in a cell. Depending on the method used, such measurements may involve the genome-wide assessment of gene expression, but also the measurement of the expression level of selected genes, resulting in the establishment of a “gene expression profile” or “expression profile”, which terms are used in that meaning hereinbelow.
- an "expression profile” comprises one or more values corresponding to a measurement of the relative abundance of a gene expression product. Such values may include measurements of RNA levels or protein abundance.
- the expression profile can comprise values representing the measurement of the transcriptional state or the translational state of the gene.
- the transcriptional state of a sample includes the idensities and relative abundance of the RNA species, especially mRNAs present in the sample. Preferably, a substantial fraction of all constituent RNA species in the sample are measured, but at least a sufficient fraction to characterize the transcriptional state of the sample is measured.
- the transcriptional state can be conveniently determined by measuring transcript abundance by any of several existing gene expression technologies.
- Translational state includes the identities and relative abundance of the constituent protein species in the sample. As is known to those of skill in the art, the transcriptional state and translational state are often related.
- Each value in the expression profiles as determined and embodied in the present invention is a measurement representing the absolute or the relative expression level of a differentially-expressed gene.
- the expression levels of these genes may be determined by any method known in the art for assessing the expression level of an RNA or protein molecule in a sample.
- expression levels of RNA may be monitored using a membrane blot (such as used in hybridization analysis such as Northern, Southern, dot, and the like), or microwells, sample tubes, gels, beads or fibers (or any solid support comprising bound nucleic acids). See U.S. Patent Nos. 5,770,722, 5,874,219, 5,744,305, 5,677,195 and 5,445,934, to which explicit reference is made.
- the gene expression monitoring system may also comprise nucleic acid probes in solution.
- microarrays are used to measure the values to be included in the expression profiles. Microarrays are particularly well suited for this purpose because of the reproducibility between different experiments.
- DNA microarrays provide one method for the simultaneous measurement of the expression levels of large numbers of genes. Each array consists of a reproducible pattern of capture probes attached to a solid support. Labeled RNA or DNA is hybridized to complementary probes on the array and then detected by laser scanning. Hybridization intensities for each probe on the array are determined and converted to a quantitative value representing relative gene expression levels. See, the Experimental section. See also, U.S. Pat. Nos.
- High-density oligonucleotide arrays are particularly useful for determining the gene expression profile for a large number of RNA's in a sample.
- total RNA isolated from the sample is converted to labeled cRNA and then hybridized to an oligonucleotide array.
- Each sample is hybridized to a separate array.
- Relative transcript levels are calculated by reference to appropriate controls present on the array and in the sample. See, for example, the Experimental section.
- the values in the expression profile are obtained by measuring the abundance of the protein products of the differentially-expressed genes.
- the abundance of these protein products can be determined, for example, using antibodies specific for the protein products of the differentially-expressed genes.
- antibody refers to an immunoglobulin molecule or immunologically active portion thereof, i.e., an antigen-binding portion.
- immunologically active portions of immunoglobulin molecules include F(ab) and F(ab')2 fragments which can be generated by treating the antibody with an enzyme such as pepsin.
- the antibody can be a polyclonal, monoclonal, recombinant, e.g., a chimeric or humanized, fully human, non-human, e.g., murine, or single chain antibody. In a preferred embodiment it has effector function and can fix complement.
- the antibody can be coupled to a toxin or imaging agent.
- a full-length protein product from a differentially-expressed gene, or an antigenic peptide fragment of the protein product can be used as an immunogen.
- Preferred epitopes encompassed by the antigenic peptide are regions of the protein product of the differentially-expressed gene that are located on the surface of the protein, e.g., hydrophilic regions, as well as regions with high antigenicity.
- the antibody can be used to detect the protein product of the differentially- expressed gene in order to evaluate the abundance and pattern of expression of the protein.
- These antibodies can also be used diagnostically to monitor protein levels in tissue as part of a clinical testing procedure, e.g., to, for example, determine the efficacy of a given therapy. Detection can be facilitated by coupling (i.e., physically linking) the antibody to a detectable substance (i.e., antibody labeling).
- detectable substances include various enzymes, prosthetic groups, fluorescent materials, luminescent materials, bioluminescent materials, and radioactive materials.
- suitable enzymes include horseradish peroxidase, alkaline phosphatase, (3- galactosidase, or acetylcholinesterase;
- suitable prosthetic group complexes include streptavidin/biotin and avidin/biotin;
- suitable fluorescent materials include umbelliferone, fluorescein, fluorescein isothiocyanate, rhodamine, dichlorotriazinylamine fluorescein, dansyl chloride, quantum dots or phycoerythrin;
- an example of a luminescent material includes luminol;
- bioluminescent materials include luciferase, luciferin, and aequorin, and examples of suitable radioactive material include 125 I, 131 I, 35 S or 3 H.
- the subject profile is compared to the reference profile to determine whether the subject expression profile is sufficiently similar to the reference profile.
- the subject expression profile is compared to a plurality of reference expression profiles to select the reference expression profile that is most similar to the subject expression profile. Any method known in the art for comparing two or more data sets to detect similarity between them may be used to compare the subject expression profile to the reference expression profiles.
- the subject expression profile and the reference profile are compared using a supervised learning algorithm such as the support vector machine (SVM) algorithm, prediction by collective likelihood of emerging patterns (PCL) algorithm, the k-nearest neighbour algorithm, or the Artificial Neural Network algorithm.
- SVM support vector machine
- PCL collective likelihood of emerging patterns
- a subject expression profile shows "statistically significant similarity” or "sufficient similarity” to a reference profile
- statistical tests may be performed to determine whether the similarity between the subject expression profile and the reference expression profile is likely to have been achieved by a random event. Any statistical test that can calculate the likelihood that the similarity between the subject expression profile and the reference profile results from a random event can be used.
- the accuracy of assigning a subject to an oligodendroglial tumor class based on similarity between differentially-expressed genes is affected largely by the heterogeneity within the patient population, as is reflected by the deviation from the geometric mean. Therefore, when more accurate diagnoses are required, the stringency in evaluating the similarity between the subject and the reference profile should be increased by changing the numerical query.
- the method used for comparing a subject expression profile to one or more reference profiles is preferably carried out by re-running the subsequent analyses in a (n+1) modus by performing clustering methods as described herein. Also, in order to identify the oligodendroglial tumor class reference profile that is most similar to the subject expression profile, as performed in the methods for establishing the oligodendroglial tumor class of a subject having a brain tumor, i.e. by diagnosing presence of an oligodendroglial tumor in a subject or by classifying the oligodendroglial tumor in a subject, profiles are clustered according to similarity and it is determined whether the subject profile corresponds to a known class of reference profiles.
- this method is used wherein the clustered position of the subject profile, obtained after performing the clustering analysis of the present invention, is compared to any known oligodendroglial tumor class. If the clustered position of the subject profile is within a cluster of reference profiles, i.e. forms a cluster therewith after performing the similarity clustering method, it is said that the oligodendroglial tumor of the subject corresponds to the oligodendroglial tumor class of reference profiles.
- the expression profiles comprise values representing the expression levels of genes that are differentially-expressed in oligodendroglial tumor classes.
- the term "differentially-expressed” as used herein means that the measured expression level of a particular gene in the expression profile of one subject differs at least n-fold from the geometric mean calculated from all patient profiles.
- the expression level may be also be up-regulated or down-regulated in a sample from a subject in comparison with a sample from a normal brain sample , or in comparison with the mean of all oligodendroglial tumor patients. Examples of genes that are differentially expressed in brain tumor patients which respond to therapy and brain tumor patients which do not respond to therapy, short vs. long survivors and Ip and/or 19q LOH vs no loss are listed in Tables 3, 4, 5, 6 and 7.
- a differentially-expressed gene is not necessarily informative for determining the presence of different oligodendroglial tumor classes, nor is every differentially-expressed gene suitable for performing diagnostic tests.
- a cluster-specific differential gene expression is most likely to be informative only in a test among subjects having brain tumors. Therefore, a diagnostic test performed by using cluster-specific gene detection should preferably be performed on a subject in which the presence of an oligodendroglial tumor is confirmed. This confirmation may for instance be obtained by standard macroscopic and microscopic detection methods.
- the present invention provides groups of genes that are differentially-expressed in diagnostic oligodendroglial tumor biopsy and surgical resection samples of patients in different therapeutic groups (i.e. responders/non-responders, or short-survivors/long-survivors). Values representing the expression levels of the nucleic acid molecules detected by the probes were analyzed as described in the Experimental section using Omniviz and SAM analysis tools. Omniviz software was used to perform all clustering steps such as K-means, Hierarchical and Pearson correlation tests. SAM was used specifically to identify the genes underlying the clinically relevant groups identified in the Pearson correlation analysis. PAM is used to decide the minimum number of genes necessary to diagnose all individual patients within the given groups of the Pearson correlation.
- the present invention thus provides a method of classifying oligodendroglial tumors.
- a total of 28 brain tumor samples analysed on a DNA microarray consisting of 54675 probe sets, representing approximately 23000 genes could be classified.
- the classification into patient groups was performed on the basis of strong correlation between their individual differential expression profiles within a group for 1881 probe sets ( ⁇ 1413 genes).
- the methods used to analyze the expression level values to identify differentially-expressed genes were employed such that optimal results in clustering, i.e. unsupervised ordering, were obtained.
- the genes that defined the position or clustering of these patient groups could be determined and the minimal sets of genes required to accurately predict the prognostically important classes could be derived.
- the method for classifying oligodendroglial tumors according to the present invention may result in a distinct pattern and therefore in a different classification scheme when other (numbers of) subjects are used as reference, or when other types of oligonucleotide microarrays for establishing gene expression profiles are used.
- the present invention thus provides a comprehensive classification of oligodendroglial tumors covering previously identified therapeutically defined classes. Further analysis of classes by significance analysis of microarrays (SAM) to determine the minimum number of genes that defined or predicted these classes resulted in the establishment of cluster-specific genes or signature genes.
- SAM significance analysis of microarrays
- the methods of the present invention comprise in some aspects the step of defining cluster-specific genes by selecting those genes of which the expression level characterizes the clustered position of the corresponding oligodendroglial tumor class within a classification scheme of the present invention.
- Such cluster-specific genes are selected preferably on the basis of SAM analysis. This method of selection comprises the following.
- the methods of the present invention comprise in some aspects the step of establishing whether the level of expression of cluster-specific genes in a subject shares sufficient similarity to the level of expression that is characteristic for an individual oligodendroglial tumor class. This step is necessary in determining the presence of that particular oligodendroglial tumor class in a subject under investigation, in which case the expression of that gene is used as a prognostic marker. Whether the level of expression of cluster-specific genes in a subject shares sufficient similarity to the level of expression of that particular gene in an individual oligodendroglial tumor class may for instance be determined by setting a threshold value.
- the present invention also reveals genes with a high differential level of expression in specific oligodendroglial tumor classes compared to the geometric mean of all reference subjects. These highly differentially-expressed genes are selected from the genes shown in Tables 3-7, These genes and their expression products are useful as markers to predict the responsiveness to treatment, Ip and/or 19q loss of heterozygosity or survival chance in a patient. Antibodies or other reagents or tools may be used to detect the presence of these markers of brain tumor.
- the present invention also reveals gene expression profiles comprising values representing the expression levels of genes in the various identified oligodendroglial tumor classes. In a preferred embodiment, these expression profiles comprise the values representing the differential expression levels.
- the expression profiles of the invention comprise one or more values representing the expression level of a gene having differential expression in a defined oligodendroglial tumor class.
- Each expression profile contains a sufficient number of values such that the profile can be used to distinguish treatment response groups, to distinguish groups with different survival, an to distinguish groups with Ip and/or 19q LOH.
- the expression profile comprises more than one or two values corresponding to a differentially-expressed gene, for example at least 3 values, at least 4 values, at least 5 values, at least 6 values, at least 7 values, at least 8 values, at least 9 values, at least 10 values, at least 11 values, at least 12 values, at least 13 values, at least 14 values, at least 15 values, at least 16 values, at least 17 values, at least 18 values, at least 19 values, at least 20 values, at least 22 values, at least 25 values, at least 27 values, at least 30 values, at least 35 values , at least 40 values, at least 45 values, at least 50 values, at least 75 values, at least 100 values, at least 125 values, at least 150 values, at least 175 values, at least 200 values, at least 250 values, at least 300 values, at least 400 values, at least 500 values, at least 600 values, at least 700 values, at least 800 values, at least 900 values, at least 1000 values, at least 1200 values, at least 1500 values, or at least 2000 or more values.
- the diagnostic accuracy of assigning a subject to an oligodendroglial tumor class will vary based on the number of values contained in the expression profile. Generally, the number of values contained in the expression profile is selected such that the diagnostic accuracy is at least 85%, at least 87%, at least 90%, at least 91%, at least 92%, at least 93%, at least 94%, at least 95%, at least 96%, at least 97%, at least 98%, or at least 99%, as calculated using methods described elsewhere herein, with an obvious preference for higher percentages of diagnostic accuracy.
- the diagnostic accuracy of assigning a subject to an oligodendroglial tumor class will vary based on the strength of the correlation between the expression levels of the differentially-expressed genes within that specific oligodendroglial tumor class.
- the values in the expression profiles represent the expression levels of genes whose expression is strongly correlated with that specific oligodendroglial tumor class, it may be possible to use fewer number of values (genes) in the expression profile and still obtain an acceptable level of diagnostic or prognostic accuracy.
- the strength of the correlation between the expression level of a differentially-expressed gene and a specific oligodendroglial tumor class may be determined by a statistical test of significance. For example, the chi square test used to select genes in some embodiments of the present invention assigns a chi square value to each differentially-expressed gene, indicating the strength of the correlation of the expression of that gene to a specific oligodendroglial tumor class. Similarly, the T-statistics metric and the Wilkins 1 metric both provide a value or score indicative of the strength of the correlation between the expression of the gene and its specific oligodendroglial tumor class.
- These scores may be used to select the genes of which the expression levels have the greatest correlation with a particular oligodendroglial tumor class to increase the diagnostic or prognostic accuracy of the methods of the invention, or in order to reduce the number of values contained in the expression profile while maintaining the diagnostic or prognostic accuracy of the expression profile.
- a database is kept wherein the expression profiles of reference subjects are collected and to which database new profiles can be added and clustered with the already existing profiles such as to provide the clustered position of said new profile among the already present reference profiles.
- the addition of new profiles to the database will improve the diagnostic and prognostic accuracy of the methods of the invention.
- SAM or PAM analysis tools are used to determine the strength of such correlations.
- the methods of the invention comprise the steps of providing an expression profile from a sample from a subject affected by oligodendroglial tumor and comparing this subject expression profile to one or more reference profiles that are associated with a particular oligodendroglial tumor class with a known prognosis, or a class with a favourable response to therapy.
- the prognosis of a subject affected by an oligodendroglial tumor can be predicted by determining whether the expression profile from the subject is sufficiently similar to a reference profile associated with an established prognosis, such as a good prognosis or a bad prognosis.
- a preferred intervention strategy, or therapeutic treatment can then be proposed for said subject, and said subject can be treated according to said assigned strategy.
- treatment of a subject with an oligodendroglial can be optimized according to the identified cluster.
- the present invention provides a method of determining the prognosis for a brain tumor patient, said method comprising the steps of providing a classification scheme for oligodendroglial tumors by producing such a scheme according to a method of the invention for reference subjects having known post-therapy lifetimes.
- the present invention provides for the assignment of the various clinical data recorded to reference subjects affected by brain tumors. This assignment preferably occurs in a database. This has the advantage that once a new subject is identified as belonging to a particular oligodendroglial tumor class, then the prognosis that is assigned to that class may be assigned to that subject.
- compositions that are useful in determining the gene expression profile for a subject affected by an oligodendroglial tumor and selecting a reference profile that is similar to the subject expression profile.
- compositions include arrays comprising a substrate having capture probes that can bind specifically to nucleic acid molecules that are differentially-expressed in oligodendroglial tumor classes.
- a computer-readable medium having digitally encoded reference profiles useful in the methods of the claimed invention.
- arrays comprising capture probes for detection of polynucleotides (transcriptional state) or for detection of proteins (translational state) in order to detect differentially-expressed genes of the invention.
- array is intended a solid support or substrate with peptide or nucleic acid probes attached to said support or substrate.
- Arrays typically comprise a plurality of different nucleic acid or peptide capture probes that are coupled to a surface of a substrate in different, known locations.
- arrays may generally be produced using mechanical synthesis methods or light directed synthesis methods which incorporate a combination of photolithographic methods and solid phase synthesis methods.
- oligonucleotide microarrays will be used for determining the transcriptional state
- peptide microarrays will be used for determining the translational state of a cell.
- Nucleic acid or “oligonucleotide” or “polynucleotide” or grammatical equivalents used herein means at least two nucleotides covalently linked together. Oligonucleotides are typically from about 5, 6, 7, 8, 9, 10, 12, 15, 25, 30, 40, 50 or more nucleotides in length, up to about 100 nucleotides in length. Nucleic acids and polynucleotides are a polymers of any length, including longer lengths, e.g., 200, 300, 500, 1000, 2000, 3000, 5000, 7000, 10,000, etc.
- a nucleic acid of the present invention will generally contain phosphodiester bonds, although in some cases, nucleic acid analogs are included that may have alternate backbones, comprising, e.g., phosphoramidate, phosphorothioate, phosphorodithioate, or O- methylphophoroamidite linkages (see Eckstein, Oligonucleotides and Analogues: A Practical Approach, Oxford University Press); and peptide nucleic acid backbones and linkages.
- Other analog nucleic acids include those with positive backbones; non-ionic backbones, and non-ribose backbones, including those described in U.S. Pat. Nos.
- nucleic acids containing one or more carbocyclic sugars are also included within one definition of nucleic acids. Modifications of the ribose-phosphate backbone may be done for a variety of reasons, e.g. to increase the stability and half-life of such molecules in physiological environments or as probes on a biochip. Mixtures of naturally occurring nucleic acids and analogues can be made; alternatively, mixtures of different nucleic acid analogues, and mixtures of naturally occurring nucleic acids and analogues may be made.
- PNA peptide nucleic acids
- These backbones are substantially non-ionic under neutral conditions, in contrast to the highly charged phosphodiester backbone of naturally occurring nucleic acids. This results in two advantages.
- the PNA backbone exhibits improved hybridization kinetics. PNAs have larger changes in the melting temperature (T m ) for mismatched versus perfectly matched basepairs. DNA and RNA typically exhibit a 2-4 0 C drop in T m for an internal mismatch. With the non-ionic PNA backbone, the drop is closer to 7-9 0 C.
- PNAs are not degraded by cellular enzymes, and thus can be more stable.
- the nucleic acids may be single stranded or double stranded, as specified, or contain portions of both double stranded or single stranded sequence.
- the depiction of a single strand also defines the sequence of the complementary strand; thus the sequences described herein also provide the complement of the sequence.
- the nucleic acid may be DNA, both genomic and cDNA, RNA or a hybrid, where the nucleic acid may contain combinations of deoxyribo- and ribo-nucleotides, and combinations of bases, including uracil, adenine, thymine, cytosine, guanine, inosine, xanthine hypoxanthine, isocytosine, isoguanine, etc.
- Transcript typically refers to a naturally occurring RNA, e.g., a pre-mRNA, hnRNA, or mRNA.
- nucleoside includes nucleotides and nucleoside and nucleotide analogues, and modified nucleosides such as amino modified nucleosides.
- nucleoside includes non-natural occurring analogue structures. Thus, e.g. the individual units of a peptide nucleic acid, each containing a base, are referred to herein as a nucleoside.
- nucleic acid probe or oligonucleotide is defined as a nucleic acid capable of binding to a target nucleic acid of complementary sequence through one or more types of chemical bonds, usually through complementary base pairing, usually through hydrogen bond formation.
- a probe may include natural (i.e., A, G, C, or T) or modified bases (7-deazaguanosine, inosine, etc.).
- the bases in a probe may be joined by a linkage other than a phosphodiester bond, so long as it does not functionally interfere with hybridization.
- probes may be peptide nucleic acids in which the constituent bases are joined by peptide bonds rather than phosphodiester linkages. It will be understood by one of skill in the art that probes may bind target sequences lacking complete complementarity with the probe sequence depending upon the stringency of the hybridization conditions.
- the probes are preferably directly labeled such as with isotopes, chromophores, lumiphores, chromogens, or indirectly labeled such as with biotin to which a streptavidin complex may later bind or with enzymatic labels.
- By assaying for the hybridization of the probe to its target nucleic acid sequence one can detect the presence or absence of the select sequence or subsequence. Diagnosis or prognosis may be based at the genomic level, or at the level of RNA or protein expression.
- oligonucleotide probes that can be used in diagnostic methods of the present invention.
- such probes are immobilised on a solid surface as to form an oligonucleotide microarray of the invention.
- the oligonucleotide probes useful in methods of the present invention are capable of hybridizing under stringent conditions to oligodendroglial tumor- associated nucleic acids, such as to one or more of the genes selected from Table 2 or Table 3. Techniques for the synthesis of arrays using mechanical synthesis methods are described in, e.g., U.S. Patent No. 5,384,261, to which reference is made herein.
- arrays may be peptides or nucleic acids on beads, gels, polymeric surfaces, fibers such as fiber optics, glass or any other appropriate substrate, for the purpose of which reference is made to U.S. Pat. Nos. 5,770,358, 5,789,162, 5,708,153, 6,040,193 and 5,800,992. Arrays may be packaged in such a manner as to allow for diagnostics or other manipulation of an all-inclusive device. Reference is for example made to U.S. Pat. Nos. 5,856,174 and 5,922,591.
- the arrays provided by the present invention comprise capture probes that can specifically bind a nucleic acid molecule that is differentially- expressed in oligodendroglial tumor classes. These arrays can be used to measure the expression levels of nucleic acid molecules to thereby create an expression profile for use in methods of determining the therapeutic treatment and prognosis for oligodendroglial tumor patients.
- each capture probe in the array detects a nucleic acid molecule selected from the nucleic acid molecules designated in Tables 2 or Table 3.
- the designated nucleic acid molecules include those differentially-expressed in oligodendroglial tumor classes.
- the arrays of the invention comprise a substrate having a plurality of addresses, where each address has a capture probe that can specifically bind a target nucleic acid molecule.
- the number of addresses on the substrate varies with the purpose for which the array is intended.
- the arrays may be low-density arrays or high-density arrays and may contain 4 or more, 8 or more, 12 or more, 16 or more, 20 or more, 24 or more, 32 or more, 48 or more, 64 or more, 72 or more 80 or more, 96, or more addresses, or 192 or more, 288 or more, 384 or more, 768 or more, 1536 or more, 3072 or more, 6144 or more, 9216 or more, 12288 or more, 15360 or more, or 18432 or more addresses.
- the substrate has no more than 12, 24, 48, 96, or 192, or 384 addresses, no more than 500, 600, 700, 800, or 900 addresses, or no more than 1000, 1200, 1600, 2400, or 3600
- the invention also provides a computer-readable medium comprising one or more digitally encoded expression profiles, where each profile has one or more values representing the expression of a gene that is differentially-expressed in an oligodendroglial tumor class.
- the preparation and use of such profiles is well within the reach of the skilled person (see e.g. WO 03/083140).
- the digitally-encoded expression profiles are comprised in a database. See, for example, U.S. Patent No. 6,308,170.
- the present invention also provides kits useful for predicting the responsiveness to therapy in subjects affected by an oligodendroglial tumor. These kits comprise an array and a computer readable medium.
- the array comprises a substrate having addresses, where each address has a capture probe that can specifically bind a nucleic acid molecule (by using an oligonucleotide array) or a peptide (by using a peptide array) that is differentially-expressed in an oligodendroglial tumor class.
- the results are converted into a computer-readable medium that has digitally-encoded expression profiles containing values representing the expression level of a nucleic acid molecule detected by the array.
- mRNA in the sample is labeled, or labeled cDNA is prepared by using mRNA as a template, and the labeled mRNA or cDNA is subjected to hybridization with the array, so that mRNAs being expressed in the sample are simultaneously detected, whereby their expression levels can be determined.
- Genes each of which expression is altered due to an oligodendrioglial tumor can be found by determining expression levels of various genes in the tumor cells and classified into certain types as described above and comparing the expression levels with the expression level in a control tissue.
- the method for determining the expression levels of genes is not particularly limited, and any of techniques for confirming alterations of the gene expressions mentioned above can be suitably used. Among all, the method using the array is especially preferable because the expressions of a large number of genes can be simultaneously determined. Suitable arrays are commercially available, e.g., from Affymetrix.
- mRNA is prepared from tumor cells, and then reverse transcription is carried out with the resulting mRNA as a template.
- labeled cDNA can be obtained by using, for instance, any suitable labeled primers or labeled nucleotides.
- the labeling substance used for labeling there can be used substances such as radioisotopes, fluorescent substances, chemiluminescent substances and substances with fluophor, and the like.
- the fluorescent substance includes Cy2, FluorX, Cy3, Cy3.5, Cy5, Cy5.5, Cy7, fluorescein isothiocyanate (FITC), Texas Red, Rhodamine and the like.
- FITC fluorescein isothiocyanate
- samples to be tested cancer samples to be tested in the present selection method
- a sample to be used as a control are each labeled with different fluorescent substances, using two or more fluorescent substances, from the viewpoint of enabling simultaneous detection.
- labeling of the samples is carried out by labeling mRNA in the samples, cDNA derived from the mRNA, or nucleic acids produced by transcription or amplification from cDNA.
- the hybridization is carried out between the above-mentioned labeled cDNA and the array to which a nucleic acid corresponding to a suitable gene or its fragment is immobilized.
- the hybridization may be performed according to any known processes under conditions that are appropriate for the array and the labeled cDNA to be used. For instance, the hybridization can be performed under the conditions described in Molecular Cloning, A laboratory manual, 2nd ed., 9.52-9.55 (1989).
- the hybridization between the nucleic acids derived from the samples and the array is carried out, under the above-mentioned hybridization conditions.
- the degradation of mRNA may take place due to actions of ribonuclease.
- the mRNA levels in both of these samples are adjusted using a standard gene with relatively little alterations in expressions.
- genes exhibiting differential expression levels in both samples can be detected.
- a signal which is appropriate depending upon the method of labeling used is detected for the array which is subjected to hybridization with the nucleic acid sample labeled by the method as described above, whereby the expression levels in the samples to be tested can be compared with the expression level in the control sample for each of the genes on the array.
- the genes thus obtained which have a significant difference in signal intensities are genes each of which expression is altered specifically for certain oligodendroglial tumor classes.
- the present invention also provides a computer-readable medium comprising a plurality of digitally-encoded expression profiles wherein each profile of the plurality has a plurality of values, each value representing the expression of a gene that is differentially-expressed in an oligodendroglial tumor class.
- the invention also provides for the storage and retrieval of a collection of data relating to oligodendroglial tumor specific gene expression data of the present invention, including sequences and expression levels in a computer data storage apparatus, which can include magnetic disks, optical disks, magneto-optical disks, DRAM, SRAM, SGRAM, SDRAM, RDRAM, DDR RAM, magnetic bubble memory devices, and other data storage devices, including CPU registers and on-CPU data storage arrays.
- kits are also provided by the invention.
- such kits may include any or all of the following: assay reagents, buffers, oligodendroglial tumor class-specific nucleic acids or antibodies, hybridization probes and/or primers, antisense polynucleotides, ribozymes, arrays, antibodies, Fab fragments, capture peptides etc.
- kits may include instructional materials containing directions (i.e., protocols) for the practice of the methods of this invention.
- instructional materials typically comprise written or printed materials, they are not limited to such. Any medium capable of storing such instructions and communicating them to an end user is contemplated by this invention.
- Such media include, but are not limited to electronic storage media (e.g., magnetic discs, tapes, cartridges, chips), optical media (e.g., CD ROM), and the like.
- Such media may include addresses to internet sites that provide such instructional materials.
- One such internet site may provide a database of oligodendroglial tumor reference expression profiles useful for performing similarity clustering of a newly determined subject expression profiles with a large set of reference profiles of oligodendroglial subjects comprised in said database.
- the database includes clinically relevant data such as patient prognosis, effects of methods of treatment and other characteristics relating to the oligodendroglial tumor patient.
- kits comprising an array of the invention and a computer-readable medium having digitally-encoded reference profiles with values representing the expression of nucleic acid molecules detected by the arrays. These kits are useful for assigning a brain tumor patient subject to an oligodendroglial tumor class.
- kits-of-parts according to the invention comprises an oligonucleotide microarray according to the invention and means for comparing a gene expression profile determined by using said microarray with a database of oligodendroglial reference expression profiles.
- the present invention also comprises kits of parts suitable for performing a method of the invention as well as the use of the various products of the invention, including databases, microarrays, oligonucleotide probes and classification schemes in diagnostic or prognostic methods of the invention.
- the present invention discloses a number of genes that are differentially-expressed in oligodendroglial tumor classes. These differentially- expressed genes are shown in Tables 3-7. Many of the treatment sensitivity- associated transcripts (Table 3) are involved in transcriptional regulation, interaction with the extracellular matrix or affect cytoskeletal dynamics. For example genes involved in regulation of transcription include: i) PAX8, a member of the paired box gene family of transcription factors; ii) SpIlO, a protein that can function as an activator of transcription; iii) RENTl, a protein involved in mRNA nuclear export and nonsense-mediated mRNA decay; and iv) TNFSF13, a member of the tumor necrosis factor ligand family that activate transcription via e.g.
- NF- ⁇ B TNFSF13 transgenic mice develop lymphoid tumors (Planelles, L. et al., (2004) Cancer Cell 6:399-408).
- Transcripts involved in the cellular interaction with the extracellular matrix include: i) MANlCl, an ⁇ -mannosidase involved in the maturation of N-linked glycans; ii) CHSYl, which synthesizes chondroitin sulfate, a widely expressed glycosaminoglycan and iii) LGALS9, a member of the tandem-repeat type galectins that bind beta-galactoside. LGALS9 is expressed at high levels in distant metastasis of breast cancer (for review see (Hirashima, M.
- IQGAPl a scaffolding protein that interacts with components of the cytoskeleton. Overexpression of IQGAPl enhances cell migration (Mataraza, J.M. et al., (2003) J. Biol. Chem. 278: 41237-41245).
- genes expressed at high levels in chemoresistant oligodendroglial tumors include i) AQPl, a water channel often highly expressed in malignant gliomas that plays a role in migration and neovascularization of tumors; ii) TRIM56, a member of the tripartite motif family and iii) ARH, an adaptor protein that interacts with the LDL receptor.
- AQPl a water channel often highly expressed in malignant gliomas that plays a role in migration and neovascularization of tumors
- TRIM56 a member of the tripartite motif family
- ARH an adaptor protein that interacts with the LDL receptor.
- genes that are differentially expressed suggests that different molecular pathways are affected in the tumors of short and long survivors.
- the genetic background of the tumor therefore appears to be an important factor in determining the prognosis of the patient, although other factors also can contribute significantly to patient survival (e.g. tumor location). Therefore, genes that are differentially expressed between long and short survivors can help identify patient subgroups that are associated with favorable prognosis.
- Functional analysis reveals that many transcripts upregulated in short survivors are involved in the regulation of transcription.
- Examples include, i) BTEBl, a member of the SPl-like/KLF family of transcription regulators, ii) BCLlO, an activator NF- ⁇ B, iii) DRl, a transcriptional repressor, iv) JUN, part of the API transcription factor complex, v) PTPN12 and vi) PTP4A2, members of the protein tyrosine phosphatase family that regulate processes including cell growth, differentiation, mitotic cycle, and oncogenic transformation, vii) SFRS4, a member of the SR family of splicing factors, and viii) LMO4, a LIM domain containing protein that may play a role as a transcriptional regulator.
- transcripts encoding proteins involved in RNA translation are downregulated in short survivors. They include five ribosomal proteins (RPL24, RPL3, PRL7, RPLP2 and RPS3) and proteins involved in post- transcriptional modification like CUGBPl and RBM4.
- This invention shows that expression profiling can identify transcripts associated with chromosomal aberrations, therapeutic response and survival after diagnosis in patients suffering from oligodendroglial tumors.
- Tumor samples Patients were chosen with (anaplastic) oligodendroglioma or mixed oligoastrocytoma with enhancing disease at the time of chemotherapy. Patients were treated in a single institution (Erasmus MC) in clinical trials evaluating the efficacy of Temozolomide or PCV. Only patients with an evaluable for response to chemotherapy were included in this study. Treatment response was evaluated by MRI and scored according to
- Tumor size was defined as the product of the two largest perpendicular tumor diameters.
- Complete response (CR) was defined as disappearance of all contrast-enhancing tumor on two subsequent scans at least one month apart, the patient being off steroids and neurologically stable or improved.
- Partial response (PR) was defined as >50 % reduction in tumor area on two subsequent scans at least one-month apart, steroids stable or decreased and neurologically stable or improved.
- Progressive disease (PD) was defined as >25 % increase in tumor area, new tumor on MRI or neurological deterioration and steroids stable or increased. All other situations were considered stable disease (SD).
- Samples were collected immediately after surgical resection, snap frozen, and stored at -800C in the Erasmus MC brain tumor tissue bank. Samples were visually inspected on 10 ⁇ m H&E stained frozen sections by the neuropathologist (J.M.K). Samples with less than 80% tumor were omitted from this study. Tissue adjacent to the inspected sections was subsequently used for nucleic acid isolation. Using these criteria, 28 oligodendroglial tumors were selected (Table 1). Four additional tumor samples with insufficient RNA quantity for array analysis were selected for confirmation of differentially expressed genes using QPCR.
- RNA and genomic DNA were extracted using Trizol (Life-Technologies) according to the manufacturers instructions.
- Genomic DNA present in the organic phase was precipitated using ethanol, washed in 0.1M Na-citrate, 10% ethanol and dissolved in 8 mM NaOH whereafter the pH was adjusted to 8.4 using IM Hepes (free acid).
- RNA quality was assessed on agarose gel and Bioanalyser (Agilent).
- cDNA synthesis and cRNA labeling was performed using the alternative protocol for one-cycle cDNA synthesis.
- Biotin-labeled cRNA was generated using the ENZO Highyield RNA transcript labeling kit (ENZO life sciences inc, NY).
- Affymetrix (Santa Clara, CA) HG U133-plus2 microarrays were hybridized overnight with 15 ⁇ g biotin labeled cRNA.
- 54.675 probesets (a probeset is a set of oligonucleotide probes that examines the expression of a single transcript) are spotted on these arrays allowing expression profiling of virtually all human transcripts. Multiple probesets may be directed against the same transcript.
- Microarrays were then washed using fluidics stations according to standard Affymetrix protocols. Microsatellite analysis
- Microsatellites were amplified by PCR on 10 ng genomic DNA using forward and reversed primers and a fluorescently labeled M13 (-21) primer. Primers and cycling conditions are stated in supplementary table 1. PCR products were precipitated, dissolved in formamide and run on an ABI 3100 genetic analyzer (Applied Biosystems). Samples were analyzed using Genescan 3.7 software (Applied Biosystems) and scored by two independent researchers. Since nonneoplastic tissues were not available for most of the tumor samples, allelic losses were statistically determined as described (Harkes I.C., et al. (2003) Br. J. Cancer 89:2289-2292). Allelic loss was assumed when the tumor sample had a homozygous allele pattern for all microsatellites within the locus (P ⁇ 0.05 for each locus).
- Fluorescence In Situ Hybridization lp/19q status of samples with non-informative microsatellite analysis was determined using Fluorescence In Situ Hybridization (FISH) as previously described (Stege E.M. et al., (2005) Cancer 103:802-809).
- Locus-specific probes for Ip36 (D1S32), centromere 1 (pUCl.77), 19ql3.4 (Bac clone 426G3), and 19pl3 (Bac clones 95711, 153P24, and 95906) were labeled with either biotin- 16-dUTP, digoxigenin-16-dUTP (Roche Diagnostics, Mannheim, Germany) or Spectrum Orange (Vysis Illinois, USA) as previously described (23). Probes were detected using FITC-labeled sheep-anti-digoxigenin (Roche Diagnostics) and/or CY3-labeled avidin (Brunschwig Chemie, Amsterdam, The Netherlands).
- Nuclei were counterstained with 4',6-diamidino-2-phenylindole (DAPI). Sixty non-overlapping nuclei were enumerated per hybridization. Ratios were calculated as the number of signals of the marker divided by the number of signals of the reference. Ratio ⁇ 0.80 were considered allelic loss. Semi-Quantitative RT-PCR
- RT-PCR Semi-quantitative RT-PCR was performed using SYBR Green PCR master mix (Applied Biosystems) according to the manufacturers instructions. Expression levels were evaluated relative to HPRT and PDGB controls. Intron spanning primers were designed against 16 genes (supplementary table 2). All primers had an amplification efficiency >80% (determined by serial dilution) and generated a single amplification product at a temperature above 77 0 C (determined by melting point analysis). Cycling was performed on an ABI7700 sequence detection system (Applied Biosystems); cycling conditions are stated in supplementary table 2. Amplification of the EFGR receptor was determined by semi-quantitative PGR using identical conditions as described above. 20 ng genomic DNA was used for each reaction.
- the amount of product amplified using genomic EGFR primers was compared to the amount of product amplified using primers on different chromosomes lying within the F3 and the FGFR3 loci.
- Statistical analysis was performed using the Mann Whitney U test (eatworms.swmed.edu/ ⁇ leon/stats/utest.cgi), values are ⁇ SEM.
- SAM analysis statistical analysis of microarrays
- SAM calculates a score for each probeset on the basis of the change in expression relative to the SD of all measurements.
- analyses were performed using stringent statistical parameters with a false discovery rate (FDR) of less than 1 probeset.
- FDR false discovery rate
- PCA structures a dataset using as few variables as possible and is a mathematical way to reduce data dimensionality.
- PCA summarizes the most important variance in a dataset as principle components.
- Weighted average was used to perform most clustering analysis, in which the distance between two clusters is defined as the average of distances between all pairs of objects. Unlike clustering based on unweighted averages, the weighted average ascribes equal weight to the two branches of the dendrogram that are about to be fused. Ward's hierarchical clustering method forms groups in a manner that minimizes the loss associated with each grouping. At each step in this analysis, the two clusters whose fusion results in minimum increase in information loss are combined. Results
- EGFR amplification and LOH on 1Oq was identified in 4/28 (14%) oligodendroglial tumors, three of which showed combined EGFR amplification and 1Oq LOH.
- a response to chemotherapy was observed in 12/14 (86%) samples with Ip35.2 LOH and 6/14 (43%) without loss of Ip35.2.
- Similar results were obtained when comparing the response rate to LOH on 19q or to combined LOH on Ip and 19q (table 1). All four tumors in which the EGFR genomic region was amplified had retained both copies of Ip and 19q and showed no response to chemotherapy (progressive disease for all). 3/4 tumors with 1Oq LOH showed no response to treatment.
- Unsupervised clustering identifies a number of subgroups, summarized in Figure 1.
- a first subgroup consists mainly of control samples but also includes low-grade tumor samples. Because the amount of tumor present in all samples was high (determined by visual inspection of sections prior to the sample used for expression profiling), this close homology to control brain tissue is likely to reflect an intrinsic property of low-grade oligodendroglial tumors.
- the low- grade oligodendroglioma samples have a higher homology to samples from whole cortex than to samples from white matter.
- Group II consists of tumor samples that have LOH on Ip and 19q and has a relatively good prognosis: All but one sample respond favorably to chemotherapy and most (4/6) patients with CR are found in this group.
- Group III has the worst prognosis: None of the tumors respond to chemotherapy, the average time of survival after diagnosis was short (1.9 ⁇ 0.2 years) as was the average time after surgical resection (1.5 ⁇ 0.3 years). All tumors of this subgroup have retained both copies of Ip and 19q and are characterized by an amplification of the EGFR locus. The samples between groups II and III have a more mixed appearance, there is some degree of correlation with both groups I and group III. Many samples with PR and all samples with SD are found in this group. Survival after diagnosis and surgical resection is intermediate between groups II and III: 8.3 ⁇ 1.5 and 2.3 ⁇ 0.3 years respectively.
- oligodendroglial tumors include those that encode proteins expressed in mature oligodendrocytes: myelin associated oligodendrocyte basic protein (MOBP), myelin oligodendrocyte glycoprotein (MOG), myelin associated glycoprotein (MAG), claudin 11 (CLDNIl) and myelin basic protein (MBP).
- MOBP myelin associated oligodendrocyte basic protein
- MOG myelin oligodendrocyte glycoprotein
- MAG myelin associated glycoprotein
- CLDNIl myelin basic protein
- MBP myelin basic protein
- transcripts are expressed ( ⁇ SD) at 0.052 ⁇ 0.021 (4 probesets), 0.10 ⁇ 0.013 (4 probesets), 0.086 (1 probeset), 0.30 ⁇ 0.25 (2 probesets), and 0.21 ⁇ 0.17 (7 probesets) levels of control brain mRNA respectively.
- This downregulation was observed in each sample.
- the strong downregulation in low-grade samples confirms the hypothesis that their homology to control brain tissue (see figure 1) is a result of the genes expressed by the tumor.
- the downregulation of MOG was confirmed using RT-PCR (table 2).
- PDGFRa is often highly expressed in oligodendroglial tumors (Riemenschneider M.J. et al., (2004) Acta Neuropathol. (Berlin) 107:277-282.).
- this gene was not present in the set of tumor-associated genes identified by our screen. Closer inspection reveals that, although PDGFRa is on average upregulated 4.1 fold, the high variation of upregulation (4.1 ⁇ 4.7) indicates that this transcript is not a reliable marker for the amount of tumor present in the sample. In fact, we failed to observe any upregulation in 10/28 samples. The select upregulation of PDGFRa in a subset of samples was confirmed using RTPCR.
- the differentially expressed genes located on the lost chromosomal arm(s) 136/376 (36.1%) probesets are located on Ip, 25/64 (39.1%) on 19q and 49/60 (82%) on Ip or 19q.
- the ratio ( ⁇ SD) loss vs. no loss is 0.53 ⁇ 0.22 (Ip), 0.54 ⁇ 0.07 (19q) and 0.53 ⁇ 0.09 (Ip and 19q) indicating that loss of one allele reduces expression levels by ⁇ 50%.
- all but two of the differentially expressed probesets that are located on the lost chromosomal(s) are downregulated. This correlation between chromosomal loss and expression level therefore suggest that these genes have an allele-number dependent expression level.
- the differentially expressed genes can be identified across the entire chromosomal arms and suggests the entire arms have been lost.
- PCA Principle components analysis
- hierarchical clustering of genes associated with LOH on Ip and 19q is depicted in figure 2.
- All anaplastic oligodendrogliomas with combined loss/retention of Ip and 19q were correctly distributed by the first principal component axis, PCAl.
- This correct distribution includes 7 samples (2 samples that have retained both Ip and 19q copies and 5 samples with LOH on Ip and 19q) that were omitted from the clustering analysis.
- Further confirmation of a subset of differentially expressed genes by RT-PCR is shown in table 2 (including 4 additional oligodendroglial tumors).
- M male; F: female; ctr: normal brain; ctr/w: control brain white matter; OD oligodendroglioma (grade II); AOD: anaplastic oligodendroglioma, AOA anaplastic oligoastrocytoma; LOH: loss of heterozygosity; ampl: amplification of the EGFR locus; ther.: therapy: PCV: combination therapy of procarbazine, CCNU, and vincristine; ten ⁇ o: temozolomide. Treatment response was scored according to McDonald's criteria (20) CR: complete response; PR: partial response; SD: stable disease; PD: progressive disease. Surv tot: patient survival after diagnosis (years); Surv op: patient survival after surgical resection of the sample used in this study.
- QPCR ctr expression of the examined transcript in control samples (either no loss of Ip, 19q, no tumor or CR/PR) relative to PDGB expression levels
- QPCR marker expression of the examined transcript in test samples (either LOH on Ip, 19q, no tumor or CR/PR) relative to PDGB expression levels.
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Genetics & Genomics (AREA)
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Biophysics (AREA)
- Molecular Biology (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Biotechnology (AREA)
- General Health & Medical Sciences (AREA)
- Theoretical Computer Science (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Medical Informatics (AREA)
- Evolutionary Biology (AREA)
- Analytical Chemistry (AREA)
- Bioinformatics & Computational Biology (AREA)
- Organic Chemistry (AREA)
- Pathology (AREA)
- Immunology (AREA)
- Wood Science & Technology (AREA)
- Zoology (AREA)
- Oncology (AREA)
- Microbiology (AREA)
- Hospice & Palliative Care (AREA)
- Biochemistry (AREA)
- General Engineering & Computer Science (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
The present invention relates to methods of genetic analysis for the prediction of treatment sensitivity and survival prognosis of patients with brain tumors, especially oligodendroglial tumors. The invention provides a method for producing a classification scheme for oligodendroglial tumors comprising the steps of a) providing a plurality of reference samples, said reference samples comprising cell samples from a plurality of reference subjects suffering from oligodendroglial tumors; b) providing reference profiles by establishing a gene expression profile for each of said reference samples individually; c) clustering said individual reference profiles according to similarity, and d) assigning an oligodendroglial tumor class to each cluster.
Description
Title: Genetic brain tumor markers
The invention relates to the field of diagnosis of tumors, especially brain tumors, more especially oligodendroglial tumors, more particular to the prediction of susceptibility to treatment for patients with brain tumor.
Diffuse gliomas are the most common primary central nervous system tumors in adults (Legler, J.M. et al., (1999) J. Natl. Cancer Inst. 91; 1382-1390; Macdonald, D.R. (2003) Semin. Oncol. 30: 72-76) and it is estimated that approximately 18,000 new patients per annum are diagnosed with a primary brain tumor in the USA (CBTRUS 2004-2005 statistical report). The worldwide standard for grading and classification of these tumors is at present the WHO classification (Kleihues, P and Cavenee, W.K., World Health Organization Classification of Tumours of the Nervous System, Lyon: WHO/IARC, 2000). Based on their histological appearance gliomas can be divided into astrocytic tumors, pure oligodendroglial tumors and mixed oligoastrocytic tumors. The latter two are grouped together as oligodendroglial tumors. The oligodendrogliomas comprise approximately 20% of all gliomas, and in comparison to most other gliomas, have a relatively long average survival time (5-12 years) after diagnosis (Okamoto, Y. et al., (2004) Acta Neuropathol. 28:28; Johannesen, T.B. et al. (2003) J. Neurosurg. 99: 854-862). Two malignancy grades are recognized in oligodendrocytic tumors, Grades II (low-grade) and III (anaplastic) (Collins, V.P. (2004) J. Neurol. Neurosurg. Psych. 75 Suppl. 2: ii2-iill).
One of the striking differences between oligodendroglial tumors and other glioma subtypes is their sensitivity to therapy, especially radiotherapy and chemotherapy. The majority of oligodendroglial tumors respond favourably to chemotherapy with alkylating agents (either temolozomide or PCV, a combination therapy of procarbazine, CCNU, and vincristine), whereas other gliomas are often chemoresistant (Van den Bent, M.J. et al. (1998) Neurology 51: 1140-1145; Van den Bent, M.J. et al. (2003) J. Clin. Oncol. 21:
2525-2528). The most favourable clinical behaviour of oligodendral tumors renders it therefore important to correctly identify this subtype of gliomas. Unfortunately, histological classification and grading of gliomas has a significant subjective component. However, malignant gliomas can also be classified according to their gene expression profile (Nutt, CL. et al. (2003) Cancer Res. 63: 1602-1607).
In oligodendroglial tumors, there is a strong correlation between chromosomal aberrations and response to treatment (chemotherapy and/or radiotherapy). For example, a common genomic aberration is a combined loss of the short arm of chromosome 1 (Ip) and the long arm of chromosome 19 (19q) (Okamoto, Y et al.,2004; Cairncross J.G. et al., (1998) J. Natl. Cancer Inst.90: 1473- 1479; Kros J.M. et al., (1999) J. Pathol. 188:282-288; Smith J.S. et al., (1999) Oncogene 18:4144-4152; Thiessen B. et al., (2003) J. Neurooncol. 64:271-278; van den Bent, M.J. et al, (2003) Cancer 97:1276-1284.). Loss of heterozygosity (LOH) on both chromosomal arms is correlated with a favourable response to therapy: A response to treatment is observed in 80-90% of oligodendroglial tumors with Ip LOH and in 25-30% without Ip LOH (Cairncross, J.G. et al, 1998; Thiessen, B. et al, 2003; van den Bent, M.J. et al., 2003). Other chromosomal aberrations observed at lower frequency include LOH on 1Oq and amplification of 7pll (Kitange G. et al. (2004) Genes
Chromosomes Cancer). These aberrations are correlated with poor prognosis and are negatively correlated with LOH on Ip and 19q. This correlation between response to treatment and chromosomal aberrations can therefore help identify chemosensitive oligodendroglial tumors. However, predicting the tumors' response to treatment by its chromosomal status also incorrectly classifies a significant percentage of tumors.
Thus, there still is a need for a more accurate prediction whether a patient with oligodendroglial tumors will be responsive to treatment and/or to predict the survival of a brain tumor patient. Expression profiling can be an alternative approach to identify oligodendroglial tumors that will benefit from
therapeutic treatment. Although expression profiling has been performed on oligodendroglial tumors, mRNA expression has thus far not been correlated to treatment response.
The current inventors have now surprisingly shown that gene expression can be used to be correlated with susceptibility to treatment and increased survival,independant of the (Ip and 19q) chromosomal status of the tumor. Further, also correlations have been found between gene expression and loss of Ip and 19q.
Summary of the Invention
The invention now comprise a method for producing a classification scheme for oligodendroglial tumors comprising the steps of: a) providing a plurality of reference samples, said reference samples comprising cell samples from a plurality of reference subjects suffering from oligodendroglial tumors, with known responsiveness to therapy and survival; b) providing reference profiles by establishing a gene expression profile, matched with parameters for sensitivity to treatment and survival for each of said reference samples individually; c) clustering said individual reference profiles according to a statistical procedure, comprising:
(i) K-means clustering; (ii) hierarchical clustering; and (iii) Pearson correlation coefficient analysis; and d) assigning an oligodendroglial tumor class according to sensitivity to treatment and/or survival to each cluster.
Specifically in such a method the clustering of said gene expression profiles is performed based on the information of differentially-expressed genes and the sensitivity to treatment and/or survival of the subject, wherein, preferably, the
clustering of said gene expression profiles with respect to treatment response is performed based on the information of the genes of Table 3, whereas the clustering of said gene expression profiles with respect to survival is performed based on the information of the genes of Table 4. Another embodiment of the invention is a method for classifying an oligodendroglial tumor of a subject suffering from an glial tumor, comprising the steps of: a) providing a classification scheme for oligodendroglial tumors according to the above described method; b) providing a subject profile by establishing a gene expression profile for said subject; c) clustering the subject profile together with reference profiles; d) determining in said scheme the clustered position of said subject profile among the reference profiles, and e) assigning to said glial tumor the oligodendroglial tumor class that corresponds to said clustered position.
Preferably herein the gene expression profile with respect to treatment response comprises the expression parameters of a set of genes according to table 3, still more preferably 1 to 50 genes of the genes of table 3, whereas the gene expression profile with respect to survival comprises the expression parameters of a set of genes according to Table 4, more preferably 1 tot 50 genes of the genes of Table 4. A further embodiment of the invention is a method of determining the prognosis for a subject suffering from an oligodendroglial tumor, said method comprising the steps of: a) providing a classification scheme for oligodendroglial tumors by producing such a scheme according to the above described method; b) determining the prognosis for each olidendroglial tumor class in said scheme based on clinical records for the subjects comprised in said class; c) establishing the oligodendroglial class of a subject suffering from an oligodendroglial tumor by classifying the oligodendroglial tumor in said subject according to a method according to the invention, and
d) assigning to said subject the prognosis corresponding to the established oligodendroglial tumor class of said subject.
Alternatively, the invention provides for a method of determining the prognosis for a subject suffering from an oligodendroglial tumor, said method comprising the steps of: a) isolation of RNA from tumor cells of said subject; b) preparation of antisense, biotinylated RNA to the RNA of step a); c) hybridisation of said antisense, biotinylated DNA on Affymetrix U133A or U133 Plus2.0 GeneChips®; d) normalising the measured values for the gene set of Table 3; e) clustering the obtained data together with reference data, obtained from a reference set of patients with known prognoses; and f) determining the prognosis on basis of the subgroup/cluster to which the data of the subject are clustering. In another embodiment, the invention provides for an oligonucleotide microarray of maximal 500 probesets, comprising at least 1, preferably at least 2, more preferably at least 25, still more preferably at least 100 oligonucleotide probes which each are capable of hybridizing under stringent conditions to different genes of the oligodendroglial tumor-associated genes selected from Table 3. Alternatively, the invention provides for an oligonucleotide microarray of maximal 500 probesets, comprising at least 1, preferably at least 2, more preferably at least 25, still more preferably at least 100 oligonucleotide probes which each are capable of hybridizing under stringent conditions to different genes of the oligodendroglial tumor-associated genes selected from Table 4.
In oligodendrogliomas there is a strong correlation between LOH on lp/19q and response to treatment. In another embodiment, the invention provides for a method using an oligonucleotide microarray, which can be used for the determination of the presence of Ip LOH, 19q LOH or lp/19q LOH. Particularly, the microarray for these determination should comprise the
genesets of Table 5, 6 and 7, respectively. Accordingly, the invention also comprises an oligonucleotide microarray of maximal 500 probesets, comprising at least 1, preferably at least 2, more preferably at least 25, still more preferably at least 50 oligonuclaotide probes which each are capable of hybridizing under stringent conditions to different genes of the oligodendroglial tumor-associated genes selected form Table 5, 6 and 7, respectively.
For the above described methods, the invention also comprises a kit- of-parts comprising an oligonucleotide microarray as described above and means for comparing a gene expression profile determined by using said microarray with a database of oligodendroglial tumor reference expression profiles.
LEGENDS TO THE FIGURES Figure 1.
Correlation plot of all samples. Samples are plotted against each other to determine the degree of similarity based on expressed genes. Red and blue denote high and low similarity respectively (scale bar). Below the correlation plot is a graphic representation of histological and patient data. Tissue: origin of sample HiI control cortex, CU control white matter, I I low-grade oligodendroglioma.^] anaplastic oligodendroglioma. Ip, 19q, 1Oq L017:|||j no LOH, L_J LOH. (LOH: loss of heterozygosity). EGFR ampl: amplification of the EFGR chromosomal locus: |ϋf no amplification, | ] amplification, Response: response to therapy IS complete response, I — I partial response, I — I stable disease, I l progressive disease. Surυ tot: survival (years) from time of diagnosis, jjj >10, Q 7-10, Q 3-7, Q] <3.A: patient alive at time of analysis.
Figure 2.
Principle components analysis (PCA) and hierarchical clustering of 60 probesets differentially expressed between oligodendroglial tumors with combined Ip and 19q LOH and those that have retained both Ip and 19q arms. A: samples are separated on their Ip and 19q chromosomal status by the first principle component axis (PCAl) whereas PCA2 separates control brain from anaplastic oligodendroglial tumors. The Ip and 19q status are color coded with
BH = no LOH on Ip and 19q,[]^| = LOH on Ip and 19q, and [^J LOH on either Ip or 19q. B: Hierarchical clustering shows relative expression levels of individual genes (columns) plotted against individual tumor samples (rows). For clarity, control brain samples were omitted from the clustering analysis. Gene expression levels are color coded with red and green indicating high (+2) and low green (-2) expression respectively (on a Iog2 scale). Dendrograms denote hierarchical clustering (Euclidian distance) of samples (top) and genes (left). The Ip and 19q status in indicated below the hierarchical clustering ( mn = no LOH, [ J = LOH). As can be seen, hierarchical clustering clearly identifies two main subgroups associated with lp/19q LOH.
Figure 3. PCA and hierarchical clustering based on 16 probesets differentially expressed between chemosensitive (CR+PR (complete response, partial response)) and chemoresistant (SD+PD, stable disease, progressive disease)) oligodendroglial tumors. A: samples are separated on their response to chemotherapy by the first principle component axis (PCAl) whereas PCA2 separates control brain from anaplastic oligodendroglial tumors. B: Hierarchical clustering based on 16 differentially expressed probesets. Relative expression levels of individual genes (columns) are plotted against individual tumor samples (rows). Gene expression levels are color coded with red and green indicating high (+1.8) and low green (-1.8) expression respectively. Dendrograms denote hierarchical clustering of samples (top) and genes (left) using Wards method. Hierarchical
clustering separates tumors that fully respond to chemotherapy (CR) from tumors that do not respond (SD+PD). Furthermore, hierarchical clustering also clearly separates tumors with poor prognosis (subgroup III in figure 1) from other oligodendroglial tumors. Responses in oligodendroglial tumors are color coded with [33 complete response, \ZH partial response, ^Hl stable disease, | [ progressive disease, | | control brain. Ip chromosomal status is depicted as LU no loss of Ip and 1 I Ip LOH.
Figure 4. PCA hierarchical clustering based on 103 probesets associated with survival after diagnosis. A: PCA identifies three main clusters of samples: oligodendroglial tumors with short survival, oligodendroglial tumors with long survival and control samples. Two low-grade samples (38 and 42, survival < 10 years I I ) cluster between control and tumor samples. PCA analysis separates short vs. long survivors on the first principle component axis (PCAl) whereas control and tumor samples are separated by the second PCA axis. B: Hierarchical clustering based on 103 differentially expressed probesets. Relative expression levels of individual genes (columns) are plotted against individual tumor samples (rows). Gene expression levels are color coded with red and green indicating high (+2) and low green (-2) expression respectively. Dendrograms denote hierarchical clustering (Euclidian distance) of samples (top) and genes (left). Interestingly, the subgroups identified by hierarchical clustering are virtually identical to the subgroups that were identified by unsupervised clustering (figure 1). Survival after diagnosis is depicted astasia >10 years survival.) 1 <10 years survival, 1 I <7 years survival, I 1 <4 years survival, | | patient still alive or, I — I control brain.
Detailed Description of the Invention
The current inventors performed expression profiling on oligodendroglial tumors and correlated the results to response to treatment, survival after diagnosis and common chromosomal aberrations. One of the findings was that the chromosomal aberrations led to ~50% expression of some but not all of the genes that had been deleted, Thus, this means that it is not straightforward to use the expression data of the genes from the Ip and 19q loci for the determination of the presence of a loss of heterozygosity (LOH) in these areas. Yet, the present inventors have found that a subset of genes, which show a reduced expression when one of the chromosomal arms Ip and 19q are deleted can be used to detect these chromosomal aberrations. The genes, which can distinguish between the presence or absence of Ip have been listed in Table 5, for LOH of 19q the genes are listed in Table 6, and Table 7 gives the list of discriminating genes for combined Ip and 19q LOH. This means that gene expression data can be used for the determination of LOH of Ip and/or 19q. This is advantageous, since currently for said determination a FISH (Fluorescence In Situ Hybridisation) or LOH (loss of heterozygosity)-PCR is used, which are specialised tests, using labelled probes. Now it has been established that a similar determination can be achieved by using standard array technology.
Further, the present study shows that the currently used predictions, based on loss of Ip, were only correctly assigned to the correct treatment response group in 20/28 (71%) of the cases, both because of positive and negative misclassifications
The term "classifying" is used in its art-recognized meaning and thus refers to arranging or ordering items, i.e. gene expression profiles, by classes or categories or dividing them into logically hierarchical classes, subclasses, and sub-subclasses based on the characteristics they have in common and/or that distinguish them. In particular "classifying" refers to
assigning, to a class or kind, an unclassified item. A "class" then being a grouping of items, based on one or more characteristics, attributes, properties, qualities, effects, parameters, etc., which they have in common, for the purpose of classifying them according to an established system or scheme. The term "classification scheme" is used in its art-recognized meaning and thus refers to a list of classes arranged according to a set of pre- established principles, for the purpose of organizing items in a collection or into groups based on their similarities and differences.
The term "clustering" refers to the activity of collecting, assembling and/or uniting into a cluster or clusters items with the same or similar elements, a "cluster" referring to a group or number of the same or similar items, i.e. gene expression profiles, gathered or occurring closely together based on similarity of characteristics. "Clustered" indicates an item has been subjected to clustering. The term "clustered position" refers to the location of an individual item, i.e. a gene expression profile, in amongst a number of clusters, said location being determined by clustering said item with at least a number of items from known clusters.
The process of clustering used in a method of the present invention may be any mathematical process known to compare items for similarity in characteristics, attributes, properties, qualities, effects, parameters, etc.. Statistical analysis, such as for instance multivariance analysis, or other methods of analysis may be used. Preferably methods of analysis such as self- organising maps, hierarchical clustering, multidimensional scaling, principle component analysis, supervised learning, k-nearest neighbours, support vector machines, discriminant analysis, partial least square methods and/or Pearson's correlation coefficient analysis are used. In another preferred embodiment of a method of the present invention Pearson's correlation coefficient analysis, significance analysis of microarrays (SAM) and/or prediction analysis of microarrays (PAM) are used to cluster gene expression
profiles according to similarity. A highly preferred method of clustering comprises similarity clustering of gene expression profiles wherein the expression level of differentially-expressed genes, having markedly lower or higher expression than the geometric mean expression level determined for all genes in all profiles to be clustered, is log(2) transformed, and wherein the transformed expression levels of all differentially-expressed genes in all profiles to be clustered is clustered by using K-means. A numerical query may then be used to select a subset of genes used in the process of hierarchical clustering (Eisen et al., 1998), thus, numerical queries may be run to select differentially expressed genes relative to the calculated geometric mean to select a smaller group of genes for hierarchical clustering.
Unsupervised sample clustering using genes obtained by numerical or threshold filtering is used to identify discrete clusters of samples as well as the gene-signatures associated with these clusters. The term gene signatures is used herein to refer to the set of genes that define the discrete position of the cluster apart from all other clusters, and includes cluster-specific genes. A numerical or threshold filtering is used to select genes for the analysis that are most likely of diagnostic relevance. Hierarchical clustering allows for visualization of large variation in gene expression across samples or present in most samples, and these genes could be used for unsupervised clustering so that clustering results are not affected by the noise from absent or non- changed genes.
Thus, while K-means clustering may be performed on all genes, the Pearson correlation is preferably calculated based on a subset of genes. Generally speaking the larger the threshold for accepting a deviation or change from the geometric mean, the smaller the number of genes that is selected by this filtering procedure. Different cut-off or threshold values were used to prepare lists with different numbers of genes. The higher the number of genes selected and included on such lists, the more noise is generally encountered within the dataset, because there will be a relatively large
contribution of non-tumor pathway related genes in such lists. The filtering and selection procedure is preferably optimized such that the analysis is performed on as many genes as possible, while minimizing the noise.
All genes with changed expression values in at least one sample higher than or equal to 1.5 times the log(2) transformed expression values and genes with changed expression values lower than or equal to -1.5 times the log(2) transformed expression value means are selected for unsupervised clustering.
The subset of genes showing a markedly higher or lower expression than the geometric mean may for instance be a value that is more than 1.5 times the geometric mean value, preferably more than 2 times the geometric mean value, Likewise, a markedly lower expression than the geometric mean expression level may for instance be a value that is less than 0.8 times the geometric mean value, preferably less than 0.6 times the geometric mean value.
Independently (see Fig. 1) a Pearson correlation coefficient analysis was performed on the samples (1881 probesets), which showed that clustering of patients is feasible.
The present invention now provides several methods to accurately identify known as well as newly discovered diagnostically, prognostically and therapeutically relevant subgroups of oligodendroglial tumors, as well as methods that can predict if treatment is likely to be effective. The basis of these methods resides in the measurement of (oligodendroglial tumor-specific) gene expression in subjects suffering from brain tumors. The methods and compositions of the invention thus provide tools useful in choosing a therapy for brain tumor patients, including methods for assigning an brain tumor patient to a brain tumor class or cluster, methods of choosing a therapy for a brain tumor patient, and methods of determining the survival prognosis for a brain tumor patient.
The methods of the invention comprise in various aspects the steps of establishing a gene expression profile of subject samples, for instance of reference subjects suffering from a brain tumor or of a subject diagnosed or classified as having a brain tumor. The expression profiles of the present invention are generated from samples from subjects having a brain tumor. The samples from the subject used to generate the expression profiles of the present invention can be derived from a tumor biopsy, wherein the sample comprises preferably more than 75% tumor cells.
"Gene expression profiling" or "expression profiling" is used herein in its art-recognised meaning and refers to a method for measuring the transcriptional state (mRNA) or the translational state (protein) of a plurality of genes in a cell. Depending on the method used, such measurements may involve the genome-wide assessment of gene expression, but also the measurement of the expression level of selected genes, resulting in the establishment of a "gene expression profile" or "expression profile", which terms are used in that meaning hereinbelow. As used herein, an "expression profile" comprises one or more values corresponding to a measurement of the relative abundance of a gene expression product. Such values may include measurements of RNA levels or protein abundance. Thus, the expression profile can comprise values representing the measurement of the transcriptional state or the translational state of the gene. In relation thereto, reference is made to U.S. Pat. Nos. 6,040,138, 5,800,992, 6,020135, 6,344,316, and 6,033,860.
The transcriptional state of a sample includes the idensities and relative abundance of the RNA species, especially mRNAs present in the sample. Preferably, a substantial fraction of all constituent RNA species in the sample are measured, but at least a sufficient fraction to characterize the transcriptional state of the sample is measured. The transcriptional state can be conveniently determined by measuring transcript abundance by any of several existing gene expression technologies.
Translational state includes the identities and relative abundance of the constituent protein species in the sample. As is known to those of skill in the art, the transcriptional state and translational state are often related.
Each value in the expression profiles as determined and embodied in the present invention is a measurement representing the absolute or the relative expression level of a differentially-expressed gene. The expression levels of these genes may be determined by any method known in the art for assessing the expression level of an RNA or protein molecule in a sample. For example, expression levels of RNA may be monitored using a membrane blot (such as used in hybridization analysis such as Northern, Southern, dot, and the like), or microwells, sample tubes, gels, beads or fibers (or any solid support comprising bound nucleic acids). See U.S. Patent Nos. 5,770,722, 5,874,219, 5,744,305, 5,677,195 and 5,445,934, to which explicit reference is made. The gene expression monitoring system may also comprise nucleic acid probes in solution.
In one embodiment of the invention, microarrays are used to measure the values to be included in the expression profiles. Microarrays are particularly well suited for this purpose because of the reproducibility between different experiments. DNA microarrays provide one method for the simultaneous measurement of the expression levels of large numbers of genes. Each array consists of a reproducible pattern of capture probes attached to a solid support. Labeled RNA or DNA is hybridized to complementary probes on the array and then detected by laser scanning. Hybridization intensities for each probe on the array are determined and converted to a quantitative value representing relative gene expression levels. See, the Experimental section. See also, U.S. Pat. Nos. 6,040,138, 5,800,992 and 6,020,135, 6,033,860, and 6,344,316, to which explicit reference is made. High-density oligonucleotide arrays are particularly useful for determining the gene expression profile for a large number of RNA's in a sample.
In one approach, total RNA isolated from the sample is converted to labeled cRNA and then hybridized to an oligonucleotide array. Each sample is hybridized to a separate array. Relative transcript levels are calculated by reference to appropriate controls present on the array and in the sample. See, for example, the Experimental section.
In another embodiment, the values in the expression profile are obtained by measuring the abundance of the protein products of the differentially-expressed genes. The abundance of these protein products can be determined, for example, using antibodies specific for the protein products of the differentially-expressed genes. The term "antibody" as used herein refers to an immunoglobulin molecule or immunologically active portion thereof, i.e., an antigen-binding portion. Examples of immunologically active portions of immunoglobulin molecules include F(ab) and F(ab')2 fragments which can be generated by treating the antibody with an enzyme such as pepsin. The antibody can be a polyclonal, monoclonal, recombinant, e.g., a chimeric or humanized, fully human, non-human, e.g., murine, or single chain antibody. In a preferred embodiment it has effector function and can fix complement. The antibody can be coupled to a toxin or imaging agent. A full-length protein product from a differentially-expressed gene, or an antigenic peptide fragment of the protein product can be used as an immunogen. Preferred epitopes encompassed by the antigenic peptide are regions of the protein product of the differentially-expressed gene that are located on the surface of the protein, e.g., hydrophilic regions, as well as regions with high antigenicity. The antibody can be used to detect the protein product of the differentially- expressed gene in order to evaluate the abundance and pattern of expression of the protein. These antibodies can also be used diagnostically to monitor protein levels in tissue as part of a clinical testing procedure, e.g., to, for example, determine the efficacy of a given therapy. Detection can be facilitated by coupling (i.e., physically linking) the antibody to a detectable substance (i.e., antibody labeling). Examples of detectable substances include various
enzymes, prosthetic groups, fluorescent materials, luminescent materials, bioluminescent materials, and radioactive materials. Examples of suitable enzymes include horseradish peroxidase, alkaline phosphatase, (3- galactosidase, or acetylcholinesterase; examples of suitable prosthetic group complexes include streptavidin/biotin and avidin/biotin; examples of suitable fluorescent materials include umbelliferone, fluorescein, fluorescein isothiocyanate, rhodamine, dichlorotriazinylamine fluorescein, dansyl chloride, quantum dots or phycoerythrin; an example of a luminescent material includes luminol; examples of bioluminescent materials include luciferase, luciferin, and aequorin, and examples of suitable radioactive material include 125I, 131I, 35S or 3H.
Once the values comprised in the subject expression profile and the reference expression profile or expression profiles are established, the subject profile is compared to the reference profile to determine whether the subject expression profile is sufficiently similar to the reference profile. Alternatively, the subject expression profile is compared to a plurality of reference expression profiles to select the reference expression profile that is most similar to the subject expression profile. Any method known in the art for comparing two or more data sets to detect similarity between them may be used to compare the subject expression profile to the reference expression profiles. In some embodiments, the subject expression profile and the reference profile are compared using a supervised learning algorithm such as the support vector machine (SVM) algorithm, prediction by collective likelihood of emerging patterns (PCL) algorithm, the k-nearest neighbour algorithm, or the Artificial Neural Network algorithm. To determine whether a subject expression profile shows "statistically significant similarity" or "sufficient similarity" to a reference profile, statistical tests may be performed to determine whether the similarity between the subject expression profile and the reference expression profile is likely to have been achieved by a random event. Any statistical test that can calculate the likelihood that the similarity between the subject
expression profile and the reference profile results from a random event can be used. The accuracy of assigning a subject to an oligodendroglial tumor class based on similarity between differentially-expressed genes is affected largely by the heterogeneity within the patient population, as is reflected by the deviation from the geometric mean. Therefore, when more accurate diagnoses are required, the stringency in evaluating the similarity between the subject and the reference profile should be increased by changing the numerical query. The method used for comparing a subject expression profile to one or more reference profiles is preferably carried out by re-running the subsequent analyses in a (n+1) modus by performing clustering methods as described herein. Also, in order to identify the oligodendroglial tumor class reference profile that is most similar to the subject expression profile, as performed in the methods for establishing the oligodendroglial tumor class of a subject having a brain tumor, i.e. by diagnosing presence of an oligodendroglial tumor in a subject or by classifying the oligodendroglial tumor in a subject, profiles are clustered according to similarity and it is determined whether the subject profile corresponds to a known class of reference profiles. In assigning a subject oligodendroglial tumor to a specific oligodendroglial tumor class for instance, this method is used wherein the clustered position of the subject profile, obtained after performing the clustering analysis of the present invention, is compared to any known oligodendroglial tumor class. If the clustered position of the subject profile is within a cluster of reference profiles, i.e. forms a cluster therewith after performing the similarity clustering method, it is said that the oligodendroglial tumor of the subject corresponds to the oligodendroglial tumor class of reference profiles.
In some embodiments of the present invention, the expression profiles comprise values representing the expression levels of genes that are differentially-expressed in oligodendroglial tumor classes. The term "differentially-expressed" as used herein means that the measured expression level of a particular gene in the expression profile of one subject differs at least
n-fold from the geometric mean calculated from all patient profiles. The expression level may be also be up-regulated or down-regulated in a sample from a subject in comparison with a sample from a normal brain sample , or in comparison with the mean of all oligodendroglial tumor patients. Examples of genes that are differentially expressed in brain tumor patients which respond to therapy and brain tumor patients which do not respond to therapy, short vs. long survivors and Ip and/or 19q LOH vs no loss are listed in Tables 3, 4, 5, 6 and 7.
It should be noted that many genes will occur, of which the measured expression level differs at least n-fold from the geometric mean expression level for that gene of all reference profiles. This may for instance be due to the different physiological state of the measured cells, to biological variation or to the presence of other diseased states. Therefore, the presence of a differentially-expressed gene is not necessarily informative for determining the presence of different oligodendroglial tumor classes, nor is every differentially-expressed gene suitable for performing diagnostic tests. Moreover, a cluster-specific differential gene expression, as defined herein, is most likely to be informative only in a test among subjects having brain tumors. Therefore, a diagnostic test performed by using cluster-specific gene detection should preferably be performed on a subject in which the presence of an oligodendroglial tumor is confirmed. This confirmation may for instance be obtained by standard macroscopic and microscopic detection methods.
The present invention provides groups of genes that are differentially-expressed in diagnostic oligodendroglial tumor biopsy and surgical resection samples of patients in different therapeutic groups (i.e. responders/non-responders, or short-survivors/long-survivors). Values representing the expression levels of the nucleic acid molecules detected by the probes were analyzed as described in the Experimental section using Omniviz and SAM analysis tools. Omniviz software was used to perform all clustering steps such as K-means, Hierarchical and Pearson correlation tests. SAM was
used specifically to identify the genes underlying the clinically relevant groups identified in the Pearson correlation analysis. PAM is used to decide the minimum number of genes necessary to diagnose all individual patients within the given groups of the Pearson correlation. In short, expression profiling was carried out on biopsy material from 28 brain tumor patients. Unsupervised clustering was used to identify novel (sub)groups within the Pearson correlation following the hierarchical clustering. After running the SAM analysis the diagnostic gene-signatures (incl. cluster-specific genes) were obtained. It appeared that a clustering separating the different groups of patients could be performed on the basis of differential expression of a plurality of genes.
The present invention thus provides a method of classifying oligodendroglial tumors. Using this method, a total of 28 brain tumor samples analysed on a DNA microarray consisting of 54675 probe sets, representing approximately 23000 genes, could be classified. The classification into patient groups was performed on the basis of strong correlation between their individual differential expression profiles within a group for 1881 probe sets (~1413 genes). The methods used to analyze the expression level values to identify differentially-expressed genes were employed such that optimal results in clustering, i.e. unsupervised ordering, were obtained. The genes that defined the position or clustering of these patient groups could be determined and the minimal sets of genes required to accurately predict the prognostically important classes could be derived. It should be understood that the method for classifying oligodendroglial tumors according to the present invention may result in a distinct pattern and therefore in a different classification scheme when other (numbers of) subjects are used as reference, or when other types of oligonucleotide microarrays for establishing gene expression profiles are used. The present invention thus provides a comprehensive classification of oligodendroglial tumors covering previously identified therapeutically
defined classes. Further analysis of classes by significance analysis of microarrays (SAM) to determine the minimum number of genes that defined or predicted these classes resulted in the establishment of cluster-specific genes or signature genes. The methods of the present invention comprise in some aspects the step of defining cluster-specific genes by selecting those genes of which the expression level characterizes the clustered position of the corresponding oligodendroglial tumor class within a classification scheme of the present invention. Such cluster-specific genes are selected preferably on the basis of SAM analysis. This method of selection comprises the following.
The methods of the present invention comprise in some aspects the step of establishing whether the level of expression of cluster-specific genes in a subject shares sufficient similarity to the level of expression that is characteristic for an individual oligodendroglial tumor class. This step is necessary in determining the presence of that particular oligodendroglial tumor class in a subject under investigation, in which case the expression of that gene is used as a prognostic marker. Whether the level of expression of cluster-specific genes in a subject shares sufficient similarity to the level of expression of that particular gene in an individual oligodendroglial tumor class may for instance be determined by setting a threshold value.
The present invention also reveals genes with a high differential level of expression in specific oligodendroglial tumor classes compared to the geometric mean of all reference subjects. These highly differentially-expressed genes are selected from the genes shown in Tables 3-7, These genes and their expression products are useful as markers to predict the responsiveness to treatment, Ip and/or 19q loss of heterozygosity or survival chance in a patient. Antibodies or other reagents or tools may be used to detect the presence of these markers of brain tumor.
The present invention also reveals gene expression profiles comprising values representing the expression levels of genes in the various identified oligodendroglial tumor classes. In a preferred embodiment, these expression profiles comprise the values representing the differential expression levels. Thus, in one embodiment the expression profiles of the invention comprise one or more values representing the expression level of a gene having differential expression in a defined oligodendroglial tumor class. Each expression profile contains a sufficient number of values such that the profile can be used to distinguish treatment response groups, to distinguish groups with different survival, an to distinguish groups with Ip and/or 19q LOH. The expression profile comprises more than one or two values corresponding to a differentially-expressed gene, for example at least 3 values, at least 4 values, at least 5 values, at least 6 values, at least 7 values, at least 8 values, at least 9 values, at least 10 values, at least 11 values, at least 12 values, at least 13 values, at least 14 values, at least 15 values, at least 16 values, at least 17 values, at least 18 values, at least 19 values, at least 20 values, at least 22 values, at least 25 values, at least 27 values, at least 30 values, at least 35 values , at least 40 values, at least 45 values, at least 50 values, at least 75 values, at least 100 values, at least 125 values, at least 150 values, at least 175 values, at least 200 values, at least 250 values, at least 300 values, at least 400 values, at least 500 values, at least 600 values, at least 700 values, at least 800 values, at least 900 values, at least 1000 values, at least 1200 values, at least 1500 values, or at least 2000 or more values.
It is recognized that the diagnostic accuracy of assigning a subject to an oligodendroglial tumor class will vary based on the number of values contained in the expression profile. Generally, the number of values contained in the expression profile is selected such that the diagnostic accuracy is at least 85%, at least 87%, at least 90%, at least 91%, at least 92%, at least 93%, at least 94%, at least 95%, at least 96%, at least 97%, at least 98%, or at least
99%, as calculated using methods described elsewhere herein, with an obvious preference for higher percentages of diagnostic accuracy.
It is recognized that the diagnostic accuracy of assigning a subject to an oligodendroglial tumor class will vary based on the strength of the correlation between the expression levels of the differentially-expressed genes within that specific oligodendroglial tumor class. When the values in the expression profiles represent the expression levels of genes whose expression is strongly correlated with that specific oligodendroglial tumor class, it may be possible to use fewer number of values (genes) in the expression profile and still obtain an acceptable level of diagnostic or prognostic accuracy.
The strength of the correlation between the expression level of a differentially-expressed gene and a specific oligodendroglial tumor class may be determined by a statistical test of significance. For example, the chi square test used to select genes in some embodiments of the present invention assigns a chi square value to each differentially-expressed gene, indicating the strength of the correlation of the expression of that gene to a specific oligodendroglial tumor class. Similarly, the T-statistics metric and the Wilkins1 metric both provide a value or score indicative of the strength of the correlation between the expression of the gene and its specific oligodendroglial tumor class. These scores may be used to select the genes of which the expression levels have the greatest correlation with a particular oligodendroglial tumor class to increase the diagnostic or prognostic accuracy of the methods of the invention, or in order to reduce the number of values contained in the expression profile while maintaining the diagnostic or prognostic accuracy of the expression profile. Preferably, a database is kept wherein the expression profiles of reference subjects are collected and to which database new profiles can be added and clustered with the already existing profiles such as to provide the clustered position of said new profile among the already present reference profiles. Furthermore, the addition of new profiles to the database will improve the diagnostic and prognostic accuracy of the
methods of the invention. Preferably, in a method of the present invention SAM or PAM analysis tools are used to determine the strength of such correlations.
The methods of the invention comprise the steps of providing an expression profile from a sample from a subject affected by oligodendroglial tumor and comparing this subject expression profile to one or more reference profiles that are associated with a particular oligodendroglial tumor class with a known prognosis, or a class with a favourable response to therapy. By identifying the oligodendroglial tumor class reference profile that is most similar to the subject expression profile, e.g. when their clustered positions fall together, the subject can be assigned to an oligodendroglial tumor class. The oligodendroglial class assigned is that with which the reference profile(s) are associated. Similarly, the prognosis of a subject affected by an oligodendroglial tumor can be predicted by determining whether the expression profile from the subject is sufficiently similar to a reference profile associated with an established prognosis, such as a good prognosis or a bad prognosis. Whenever a subject's expression profile can be assigned to one of the reference profile(s), a preferred intervention strategy, or therapeutic treatment can then be proposed for said subject, and said subject can be treated according to said assigned strategy. As a result, treatment of a subject with an oligodendroglial can be optimized according to the identified cluster.
In one aspect, the present invention provides a method of determining the prognosis for a brain tumor patient, said method comprising the steps of providing a classification scheme for oligodendroglial tumors by producing such a scheme according to a method of the invention for reference subjects having known post-therapy lifetimes. The present invention provides for the assignment of the various clinical data recorded to reference subjects affected by brain tumors. This assignment preferably occurs in a database. This has the advantage that once a new subject is identified as belonging to a
particular oligodendroglial tumor class, then the prognosis that is assigned to that class may be assigned to that subject.
The present invention provides compositions that are useful in determining the gene expression profile for a subject affected by an oligodendroglial tumor and selecting a reference profile that is similar to the subject expression profile. These compositions include arrays comprising a substrate having capture probes that can bind specifically to nucleic acid molecules that are differentially-expressed in oligodendroglial tumor classes. Also provided is a computer-readable medium having digitally encoded reference profiles useful in the methods of the claimed invention.
The present invention provides arrays comprising capture probes for detection of polynucleotides (transcriptional state) or for detection of proteins (translational state) in order to detect differentially-expressed genes of the invention. By "array" is intended a solid support or substrate with peptide or nucleic acid probes attached to said support or substrate. Arrays typically comprise a plurality of different nucleic acid or peptide capture probes that are coupled to a surface of a substrate in different, known locations. These arrays, also described as "microarrays" or colloquially "chips" have been generally described in the art, and reference is made U.S. Patent. Nos. 5,143,854, 5,445,934, 5,744,305, 5,677,195, 6,040,193, 5,424,186,6,329,143, and 6,309,831 and Fodor et al. (1991) Science 251:767-77. These arrays may generally be produced using mechanical synthesis methods or light directed synthesis methods which incorporate a combination of photolithographic methods and solid phase synthesis methods. Typically, "oligonucleotide microarrays" will be used for determining the transcriptional state, whereas "peptide microarrays" will be used for determining the translational state of a cell.
"Nucleic acid" or "oligonucleotide" or "polynucleotide" or grammatical equivalents used herein means at least two nucleotides covalently linked together. Oligonucleotides are typically from about 5, 6, 7, 8, 9, 10, 12, 15, 25, 30, 40, 50 or more nucleotides in length, up to about 100
nucleotides in length. Nucleic acids and polynucleotides are a polymers of any length, including longer lengths, e.g., 200, 300, 500, 1000, 2000, 3000, 5000, 7000, 10,000, etc. A nucleic acid of the present invention will generally contain phosphodiester bonds, although in some cases, nucleic acid analogs are included that may have alternate backbones, comprising, e.g., phosphoramidate, phosphorothioate, phosphorodithioate, or O- methylphophoroamidite linkages (see Eckstein, Oligonucleotides and Analogues: A Practical Approach, Oxford University Press); and peptide nucleic acid backbones and linkages. Other analog nucleic acids include those with positive backbones; non-ionic backbones, and non-ribose backbones, including those described in U.S. Pat. Nos. 5,235,033 and 5,034,506, and Chapters 6 and 7, ASC Symposium Series 580, Carbohydrate Modifications in Antisense Research, Sanghui & Cook, eds. Nucleic acids containing one or more carbocyclic sugars are also included within one definition of nucleic acids. Modifications of the ribose-phosphate backbone may be done for a variety of reasons, e.g. to increase the stability and half-life of such molecules in physiological environments or as probes on a biochip. Mixtures of naturally occurring nucleic acids and analogues can be made; alternatively, mixtures of different nucleic acid analogues, and mixtures of naturally occurring nucleic acids and analogues may be made.
Particularly preferred are peptide nucleic acids (PNA) which includes peptide nucleic acid analogues. These backbones are substantially non-ionic under neutral conditions, in contrast to the highly charged phosphodiester backbone of naturally occurring nucleic acids. This results in two advantages. First, the PNA backbone exhibits improved hybridization kinetics. PNAs have larger changes in the melting temperature (Tm) for mismatched versus perfectly matched basepairs. DNA and RNA typically exhibit a 2-40C drop in Tm for an internal mismatch. With the non-ionic PNA backbone, the drop is closer to 7-90C. Similarly, due to their non-ionic nature, hybridization of the bases attached to these backbones is relatively insensitive
to salt concentration. In addition, PNAs are not degraded by cellular enzymes, and thus can be more stable.
The nucleic acids may be single stranded or double stranded, as specified, or contain portions of both double stranded or single stranded sequence. As will be appreciated by those in the art, the depiction of a single strand also defines the sequence of the complementary strand; thus the sequences described herein also provide the complement of the sequence. The nucleic acid may be DNA, both genomic and cDNA, RNA or a hybrid, where the nucleic acid may contain combinations of deoxyribo- and ribo-nucleotides, and combinations of bases, including uracil, adenine, thymine, cytosine, guanine, inosine, xanthine hypoxanthine, isocytosine, isoguanine, etc.
"Transcript" typically refers to a naturally occurring RNA, e.g., a pre-mRNA, hnRNA, or mRNA. As used herein, the term "nucleoside" includes nucleotides and nucleoside and nucleotide analogues, and modified nucleosides such as amino modified nucleosides. In addition, "nucleoside" includes non- naturally occurring analogue structures. Thus, e.g. the individual units of a peptide nucleic acid, each containing a base, are referred to herein as a nucleoside.
As used herein a "nucleic acid probe or oligonucleotide" is defined as a nucleic acid capable of binding to a target nucleic acid of complementary sequence through one or more types of chemical bonds, usually through complementary base pairing, usually through hydrogen bond formation. As used herein, a probe may include natural (i.e., A, G, C, or T) or modified bases (7-deazaguanosine, inosine, etc.). In addition, the bases in a probe may be joined by a linkage other than a phosphodiester bond, so long as it does not functionally interfere with hybridization. Thus, e.g., probes may be peptide nucleic acids in which the constituent bases are joined by peptide bonds rather than phosphodiester linkages. It will be understood by one of skill in the art that probes may bind target sequences lacking complete complementarity with the probe sequence depending upon the stringency of the hybridization
conditions. The probes are preferably directly labeled such as with isotopes, chromophores, lumiphores, chromogens, or indirectly labeled such as with biotin to which a streptavidin complex may later bind or with enzymatic labels. By assaying for the hybridization of the probe to its target nucleic acid sequence, one can detect the presence or absence of the select sequence or subsequence. Diagnosis or prognosis may be based at the genomic level, or at the level of RNA or protein expression.
The skilled person is capable of designing oligonucleotide probes that can be used in diagnostic methods of the present invention. Preferably, such probes are immobilised on a solid surface as to form an oligonucleotide microarray of the invention. The oligonucleotide probes useful in methods of the present invention are capable of hybridizing under stringent conditions to oligodendroglial tumor- associated nucleic acids, such as to one or more of the genes selected from Table 2 or Table 3. Techniques for the synthesis of arrays using mechanical synthesis methods are described in, e.g., U.S. Patent No. 5,384,261, to which reference is made herein. Although a planar array surface is preferred, the array may be fabricated on a surface of virtually any shape or even a multiplicity of surfaces. Arrays may be peptides or nucleic acids on beads, gels, polymeric surfaces, fibers such as fiber optics, glass or any other appropriate substrate, for the purpose of which reference is made to U.S. Pat. Nos. 5,770,358, 5,789,162, 5,708,153, 6,040,193 and 5,800,992. Arrays may be packaged in such a manner as to allow for diagnostics or other manipulation of an all-inclusive device. Reference is for example made to U.S. Pat. Nos. 5,856,174 and 5,922,591. The arrays provided by the present invention comprise capture probes that can specifically bind a nucleic acid molecule that is differentially- expressed in oligodendroglial tumor classes. These arrays can be used to measure the expression levels of nucleic acid molecules to thereby create an expression profile for use in methods of determining the therapeutic treatment and prognosis for oligodendroglial tumor patients.
In some embodiments, each capture probe in the array detects a nucleic acid molecule selected from the nucleic acid molecules designated in Tables 2 or Table 3. The designated nucleic acid molecules include those differentially-expressed in oligodendroglial tumor classes. The arrays of the invention comprise a substrate having a plurality of addresses, where each address has a capture probe that can specifically bind a target nucleic acid molecule. The number of addresses on the substrate varies with the purpose for which the array is intended. The arrays may be low-density arrays or high-density arrays and may contain 4 or more, 8 or more, 12 or more, 16 or more, 20 or more, 24 or more, 32 or more, 48 or more, 64 or more, 72 or more 80 or more, 96, or more addresses, or 192 or more, 288 or more, 384 or more, 768 or more, 1536 or more, 3072 or more, 6144 or more, 9216 or more, 12288 or more, 15360 or more, or 18432 or more addresses. In some embodiments, the substrate has no more than 12, 24, 48, 96, or 192, or 384 addresses, no more than 500, 600, 700, 800, or 900 addresses, or no more than 1000, 1200, 1600, 2400, or 3600 addresses.
The invention also provides a computer-readable medium comprising one or more digitally encoded expression profiles, where each profile has one or more values representing the expression of a gene that is differentially-expressed in an oligodendroglial tumor class. The preparation and use of such profiles is well within the reach of the skilled person (see e.g. WO 03/083140). In some embodiments, the digitally-encoded expression profiles are comprised in a database. See, for example, U.S. Patent No. 6,308,170. The present invention also provides kits useful for predicting the responsiveness to therapy in subjects affected by an oligodendroglial tumor. These kits comprise an array and a computer readable medium. The array comprises a substrate having addresses, where each address has a capture probe that can specifically bind a nucleic acid molecule (by using an oligonucleotide array) or a peptide (by using a peptide array) that is
differentially-expressed in an oligodendroglial tumor class. The results are converted into a computer-readable medium that has digitally-encoded expression profiles containing values representing the expression level of a nucleic acid molecule detected by the array. By using the array described above, the amounts of various kinds of nucleic acid molecules contained in a nucleic acid sample can be simultaneously determined. In addition, there is an advantage such that the determination can be carried out even with a small amount of the nucleic acid sample. For instance, mRNA in the sample is labeled, or labeled cDNA is prepared by using mRNA as a template, and the labeled mRNA or cDNA is subjected to hybridization with the array, so that mRNAs being expressed in the sample are simultaneously detected, whereby their expression levels can be determined.
Genes each of which expression is altered due to an oligodendrioglial tumor can be found by determining expression levels of various genes in the tumor cells and classified into certain types as described above and comparing the expression levels with the expression level in a control tissue.
The method for determining the expression levels of genes is not particularly limited, and any of techniques for confirming alterations of the gene expressions mentioned above can be suitably used. Among all, the method using the array is especially preferable because the expressions of a large number of genes can be simultaneously determined. Suitable arrays are commercially available, e.g., from Affymetrix.
For instance, mRNA is prepared from tumor cells, and then reverse transcription is carried out with the resulting mRNA as a template. During this process, labeled cDNA can be obtained by using, for instance, any suitable labeled primers or labeled nucleotides.
As to the labeling substance used for labeling, there can be used substances such as radioisotopes, fluorescent substances, chemiluminescent substances and substances with fluophor, and the like. For instance, the
fluorescent substance includes Cy2, FluorX, Cy3, Cy3.5, Cy5, Cy5.5, Cy7, fluorescein isothiocyanate (FITC), Texas Red, Rhodamine and the like. In addition, it is desired that samples to be tested (cancer samples to be tested in the present selection method) and a sample to be used as a control are each labeled with different fluorescent substances, using two or more fluorescent substances, from the viewpoint of enabling simultaneous detection. Here, labeling of the samples is carried out by labeling mRNA in the samples, cDNA derived from the mRNA, or nucleic acids produced by transcription or amplification from cDNA. Next, the hybridization is carried out between the above-mentioned labeled cDNA and the array to which a nucleic acid corresponding to a suitable gene or its fragment is immobilized. The hybridization may be performed according to any known processes under conditions that are appropriate for the array and the labeled cDNA to be used. For instance, the hybridization can be performed under the conditions described in Molecular Cloning, A laboratory manual, 2nd ed., 9.52-9.55 (1989).
The hybridization between the nucleic acids derived from the samples and the array is carried out, under the above-mentioned hybridization conditions. When much time is needed for the time period required for procedures from the collection of samples to the determination of expression levels of genes, the degradation of mRNA may take place due to actions of ribonuclease. In order to determine the difference in the gene expressions in the samples to be tested (i.e., tumor cells or biopsies from oligodendroglial tumor patients) and the gene expressions in a control sample, it is preferable that the mRNA levels in both of these samples are adjusted using a standard gene with relatively little alterations in expressions.
Thereafter, by comparing the hybridization results of the samples to be tested with those of the control sample, genes exhibiting differential expression levels in both samples can be detected. Concretely, a signal which is appropriate depending upon the method of labeling used is detected for the
array which is subjected to hybridization with the nucleic acid sample labeled by the method as described above, whereby the expression levels in the samples to be tested can be compared with the expression level in the control sample for each of the genes on the array. The genes thus obtained which have a significant difference in signal intensities are genes each of which expression is altered specifically for certain oligodendroglial tumor classes.
The present invention also provides a computer-readable medium comprising a plurality of digitally-encoded expression profiles wherein each profile of the plurality has a plurality of values, each value representing the expression of a gene that is differentially-expressed in an oligodendroglial tumor class. The invention also provides for the storage and retrieval of a collection of data relating to oligodendroglial tumor specific gene expression data of the present invention, including sequences and expression levels in a computer data storage apparatus, which can include magnetic disks, optical disks, magneto-optical disks, DRAM, SRAM, SGRAM, SDRAM, RDRAM, DDR RAM, magnetic bubble memory devices, and other data storage devices, including CPU registers and on-CPU data storage arrays. Typically, the data records are stored as a bit pattern in an array of magnetic domains on a magnetizable medium or as an array of charge states or transistor gate states, such as an array of cells in a DRAM device (e.g., each cell comprised of a transistor and a charge storage area, which may be on the transistor). For use in diagnostic, research, and therapeutic applications suggested above, kits are also provided by the invention. In the diagnostic and research applications such kits may include any or all of the following: assay reagents, buffers, oligodendroglial tumor class-specific nucleic acids or antibodies, hybridization probes and/or primers, antisense polynucleotides, ribozymes, arrays, antibodies, Fab fragments, capture peptides etc. In addition, the kits may include instructional materials containing directions (i.e., protocols) for the practice of the methods of this invention. While the
instructional materials typically comprise written or printed materials, they are not limited to such. Any medium capable of storing such instructions and communicating them to an end user is contemplated by this invention. Such media include, but are not limited to electronic storage media (e.g., magnetic discs, tapes, cartridges, chips), optical media (e.g., CD ROM), and the like. Such media may include addresses to internet sites that provide such instructional materials. One such internet site may provide a database of oligodendroglial tumor reference expression profiles useful for performing similarity clustering of a newly determined subject expression profiles with a large set of reference profiles of oligodendroglial subjects comprised in said database. Preferably the database includes clinically relevant data such as patient prognosis, effects of methods of treatment and other characteristics relating to the oligodendroglial tumor patient.
The invention encompasses for instance kits comprising an array of the invention and a computer-readable medium having digitally-encoded reference profiles with values representing the expression of nucleic acid molecules detected by the arrays. These kits are useful for assigning a brain tumor patient subject to an oligodendroglial tumor class.
In a preferred embodiment a kit-of-parts according to the invention comprises an oligonucleotide microarray according to the invention and means for comparing a gene expression profile determined by using said microarray with a database of oligodendroglial reference expression profiles. The present invention also comprises kits of parts suitable for performing a method of the invention as well as the use of the various products of the invention, including databases, microarrays, oligonucleotide probes and classification schemes in diagnostic or prognostic methods of the invention.
The present invention discloses a number of genes that are differentially-expressed in oligodendroglial tumor classes. These differentially- expressed genes are shown in Tables 3-7. Many of the treatment sensitivity- associated transcripts (Table 3) are involved in transcriptional regulation,
interaction with the extracellular matrix or affect cytoskeletal dynamics. For example genes involved in regulation of transcription include: i) PAX8, a member of the paired box gene family of transcription factors; ii) SpIlO, a protein that can function as an activator of transcription; iii) RENTl, a protein involved in mRNA nuclear export and nonsense-mediated mRNA decay; and iv) TNFSF13, a member of the tumor necrosis factor ligand family that activate transcription via e.g. NF-κB. TNFSF13 transgenic mice develop lymphoid tumors (Planelles, L. et al., (2004) Cancer Cell 6:399-408). Transcripts involved in the cellular interaction with the extracellular matrix include: i) MANlCl, an α-mannosidase involved in the maturation of N-linked glycans; ii) CHSYl, which synthesizes chondroitin sulfate, a widely expressed glycosaminoglycan and iii) LGALS9, a member of the tandem-repeat type galectins that bind beta-galactoside. LGALS9 is expressed at high levels in distant metastasis of breast cancer (for review see (Hirashima, M. et al., (2004) Glycoconj. J. ,19:593-600). Also two treatment sensitivity associated transcripts that are involved in regulation of cytoskeletal dynamics were identified: i) ARPClB, involved in the branching of actin filaments and downregulated in gastric cancers; and ii) IQGAPl, a scaffolding protein that interacts with components of the cytoskeleton. Overexpression of IQGAPl enhances cell migration (Mataraza, J.M. et al., (2003) J. Biol. Chem. 278: 41237-41245). Other genes expressed at high levels in chemoresistant oligodendroglial tumors include i) AQPl, a water channel often highly expressed in malignant gliomas that plays a role in migration and neovascularization of tumors; ii) TRIM56, a member of the tripartite motif family and iii) ARH, an adaptor protein that interacts with the LDL receptor. In summary, the genes identified in this invention that are associated with treatment sensitivity (Table 3) are involved in several discrete cellular processes and further study on these transcripts may help identify the molecular mechanisms that underlie treatment sensitivity.
Comparison of expression profiles to patient survival after diagnosis identified 103 differentially expressed probesets (Table 4). The observation that many genes are differentially expressed suggests that different molecular pathways are affected in the tumors of short and long survivors. The genetic background of the tumor therefore appears to be an important factor in determining the prognosis of the patient, although other factors also can contribute significantly to patient survival (e.g. tumor location). Therefore, genes that are differentially expressed between long and short survivors can help identify patient subgroups that are associated with favorable prognosis. Functional analysis reveals that many transcripts upregulated in short survivors are involved in the regulation of transcription. Examples include, i) BTEBl, a member of the SPl-like/KLF family of transcription regulators, ii) BCLlO, an activator NF-κB, iii) DRl, a transcriptional repressor, iv) JUN, part of the API transcription factor complex, v) PTPN12 and vi) PTP4A2, members of the protein tyrosine phosphatase family that regulate processes including cell growth, differentiation, mitotic cycle, and oncogenic transformation, vii) SFRS4, a member of the SR family of splicing factors, and viii) LMO4, a LIM domain containing protein that may play a role as a transcriptional regulator. In contrast, transcripts encoding proteins involved in RNA translation are downregulated in short survivors. They include five ribosomal proteins (RPL24, RPL3, PRL7, RPLP2 and RPS3) and proteins involved in post- transcriptional modification like CUGBPl and RBM4.
This invention shows that expression profiling can identify transcripts associated with chromosomal aberrations, therapeutic response and survival after diagnosis in patients suffering from oligodendroglial tumors.
As described above this knowledge can be used to identify patient classes with a high likelihood to respond to treatment and patient classes with favorable prognosis.
The following examples are offered by way of illustration and not by way of limitation.
EXAMPLE
Methods:
Tumor samples: Patients were chosen with (anaplastic) oligodendroglioma or mixed oligoastrocytoma with enhancing disease at the time of chemotherapy. Patients were treated in a single institution (Erasmus MC) in clinical trials evaluating the efficacy of Temozolomide or PCV. Only patients with an evaluable for response to chemotherapy were included in this study. Treatment response was evaluated by MRI and scored according to
McDonald's criteria (Macdonald D.R. et al., (1990) J. Clin. Oncol. 8:1277-1280). Tumor size was defined as the product of the two largest perpendicular tumor diameters. Complete response (CR) was defined as disappearance of all contrast-enhancing tumor on two subsequent scans at least one month apart, the patient being off steroids and neurologically stable or improved. Partial response (PR) was defined as >50 % reduction in tumor area on two subsequent scans at least one-month apart, steroids stable or decreased and neurologically stable or improved. Progressive disease (PD) was defined as >25 % increase in tumor area, new tumor on MRI or neurological deterioration and steroids stable or increased. All other situations were considered stable disease (SD). Samples were collected immediately after surgical resection, snap frozen, and stored at -800C in the Erasmus MC brain tumor tissue bank. Samples were visually inspected on 10 μm H&E stained frozen sections by the neuropathologist (J.M.K). Samples with less than 80% tumor were omitted from this study. Tissue adjacent to the inspected sections was subsequently
used for nucleic acid isolation. Using these criteria, 28 oligodendroglial tumors were selected (Table 1). Four additional tumor samples with insufficient RNA quantity for array analysis were selected for confirmation of differentially expressed genes using QPCR.
Nucleic acid isolation:
Tissues were homogenized using a polytron following which RNA and genomic DNA were extracted using Trizol (Life-Technologies) according to the manufacturers instructions. Total RNA, present in the aqueous phase after extraction, was precipitated in isopropanol, redissolved in diethyl- pyrocarbonate treated water and further purified on RNeasy mini columns (Qiagen). Genomic DNA present in the organic phase was precipitated using ethanol, washed in 0.1M Na-citrate, 10% ethanol and dissolved in 8 mM NaOH whereafter the pH was adjusted to 8.4 using IM Hepes (free acid).
cDNA synthesis and array hybridization
RNA quality was assessed on agarose gel and Bioanalyser (Agilent). cDNA synthesis and cRNA labeling was performed using the alternative protocol for one-cycle cDNA synthesis. Biotin-labeled cRNA was generated using the ENZO Highyield RNA transcript labeling kit (ENZO life sciences inc, NY). Affymetrix (Santa Clara, CA) HG U133-plus2 microarrays were hybridized overnight with 15 μg biotin labeled cRNA. 54.675 probesets (a probeset is a set of oligonucleotide probes that examines the expression of a single transcript) are spotted on these arrays allowing expression profiling of virtually all human transcripts. Multiple probesets may be directed against the same transcript. Microarrays were then washed using fluidics stations according to standard Affymetrix protocols.
Microsatellite analysis
Microsatellites were amplified by PCR on 10 ng genomic DNA using forward and reversed primers and a fluorescently labeled M13 (-21) primer. Primers and cycling conditions are stated in supplementary table 1. PCR products were precipitated, dissolved in formamide and run on an ABI 3100 genetic analyzer (Applied Biosystems). Samples were analyzed using Genescan 3.7 software (Applied Biosystems) and scored by two independent researchers. Since nonneoplastic tissues were not available for most of the tumor samples, allelic losses were statistically determined as described (Harkes I.C., et al. (2003) Br. J. Cancer 89:2289-2292). Allelic loss was assumed when the tumor sample had a homozygous allele pattern for all microsatellites within the locus (P<0.05 for each locus).
Fluorescence In Situ Hybridization lp/19q status of samples with non-informative microsatellite analysis was determined using Fluorescence In Situ Hybridization (FISH) as previously described (Stege E.M. et al., (2005) Cancer 103:802-809). Locus-specific probes for Ip36 (D1S32), centromere 1 (pUCl.77), 19ql3.4 (Bac clone 426G3), and 19pl3 (Bac clones 95711, 153P24, and 95906) were labeled with either biotin- 16-dUTP, digoxigenin-16-dUTP (Roche Diagnostics, Mannheim, Germany) or Spectrum Orange (Vysis Illinois, USA) as previously described (23). Probes were detected using FITC-labeled sheep-anti-digoxigenin (Roche Diagnostics) and/or CY3-labeled avidin (Brunschwig Chemie, Amsterdam, The Netherlands). Nuclei were counterstained with 4',6-diamidino-2-phenylindole (DAPI). Sixty non-overlapping nuclei were enumerated per hybridization. Ratios were calculated as the number of signals of the marker divided by the number of signals of the reference. Ratio <0.80 were considered allelic loss.
Semi-Quantitative RT-PCR
Semi-quantitative RT-PCR was performed using SYBR Green PCR master mix (Applied Biosystems) according to the manufacturers instructions. Expression levels were evaluated relative to HPRT and PDGB controls. Intron spanning primers were designed against 16 genes (supplementary table 2). All primers had an amplification efficiency >80% (determined by serial dilution) and generated a single amplification product at a temperature above 77 0C (determined by melting point analysis). Cycling was performed on an ABI7700 sequence detection system (Applied Biosystems); cycling conditions are stated in supplementary table 2. Amplification of the EFGR receptor was determined by semi-quantitative PGR using identical conditions as described above. 20 ng genomic DNA was used for each reaction. The amount of product amplified using genomic EGFR primers was compared to the amount of product amplified using primers on different chromosomes lying within the F3 and the FGFR3 loci. Statistical analysis was performed using the Mann Whitney U test (eatworms.swmed.edu/~leon/stats/utest.cgi), values are ± SEM.
Data analysis:
Arrays were omitted from the analysis when the number of present calls < 35% and when the 573' ratio of GAPDH controls >3. Probesets that were absent
(according to Affymetrix MAS5.0 software) in at least 33 of the 34 microarrays were omitted from further analysis. Raw intensities of the remaining probesets (36875) of each chip were Iog2 transformed and normalized using quantile normalization. For each probeset, the geometric mean of the hybridization intensities of all samples was calculated. The level of expression of each probeset was determined relative to this geometric mean and Iog2 transformed. The geometric mean of the hybridization signal of all samples was used to ascribe equal weight to gene-expression levels. Unsupervised clustering was performed using Omniviz version 3.6.0 (Omniviz, Maynard, MA) software. Probesets whose expression levels differed more than 2 fold
from the geometric mean in at least one sample were selected for the unsupervised clustering analysis. Similarities between samples is plotted using Omniviz software as Pearson's correlations.
Differentially expressed genes were identified using statistical analysis of microarrays (SAM analysis) (Tusher V.G. et al., Proc. Natl. Acad. Sci. U. S. A. 98:5116-5121). Such supervised analysis correlates gene expression with an external variable. SAM calculates a score for each probeset on the basis of the change in expression relative to the SD of all measurements. Unless otherwise indicated, analyses were performed using stringent statistical parameters with a false discovery rate (FDR) of less than 1 probeset. Differentially expressed probesets were imported into Spotfire DecisionSite (Spotfire, Somerville, MA) to perform principle components analysis (PCA) and hierarchical clustering. Data were Iog2 transformed followed by calculation of the z-score for each probeset. PCA structures a dataset using as few variables as possible and is a mathematical way to reduce data dimensionality. PCA summarizes the most important variance in a dataset as principle components. For more information on the use of PCA in microarray analysis microarrays see (Raychaudhuri S. et al. (2000), In: Hunter L, Altman B, Dunker AK, Klein TE, Lauderdale K, editors. Pacific Symposium on Biocomputing 1999. Honolulu, Hawaii: World Scientific Press; 2000) and references therein. Hierarchical clustering groups data based on their similarities in gene expression profiles. Weighted average was used to perform most clustering analysis, in which the distance between two clusters is defined as the average of distances between all pairs of objects. Unlike clustering based on unweighted averages, the weighted average ascribes equal weight to the two branches of the dendrogram that are about to be fused. Ward's hierarchical clustering method forms groups in a manner that minimizes the loss associated with each grouping. At each step in this analysis, the two clusters whose fusion results in minimum increase in information loss are combined.
Results
Samples:
Patient data, histological diagnosis, chromosomal aberrations, and response to chemotherapy are summarized in table 1. In total we performed expression analysis on 28 oligodendroglial tumors (2 lowgrade and 26 anaplastic oligodendrogliomas), and 6 control brain samples (4 samples from whole cortex, 2 from white matter only). We identified 14/28 samples (50%) with loss of most/all of the short arm of chromosome 1 (sample 18 had a predicted loss distal to Ip33) and 16/28 (57%) samples with loss of 19q (see Table 1). Most tumors showed combined loss or retention of Ip and 19q: only three tumors showed loss of 19q without loss of Ip, one showed LOH on Ip35.2 without loss of 19q. EGFR amplification and LOH on 1Oq was identified in 4/28 (14%) oligodendroglial tumors, three of which showed combined EGFR amplification and 1Oq LOH. When comparing the response rate (CR+PR vs. PD+SD) to loss of the telomeric end of chromosome 1, a response to chemotherapy was observed in 12/14 (86%) samples with Ip35.2 LOH and 6/14 (43%) without loss of Ip35.2. Similar results were obtained when comparing the response rate to LOH on 19q or to combined LOH on Ip and 19q (table 1). All four tumors in which the EGFR genomic region was amplified had retained both copies of Ip and 19q and showed no response to chemotherapy (progressive disease for all). 3/4 tumors with 1Oq LOH showed no response to treatment.
Unsuperυised clustering:
Unsupervised clustering identifies a number of subgroups, summarized in Figure 1. A first subgroup consists mainly of control samples but also includes low-grade tumor samples. Because the amount of tumor present in all samples was high (determined by visual inspection of sections prior to the sample used for expression profiling), this close homology to control brain tissue is likely to reflect an intrinsic property of low-grade oligodendroglial tumors. The low- grade oligodendroglioma samples have a higher homology to samples from
whole cortex than to samples from white matter. Group II consists of tumor samples that have LOH on Ip and 19q and has a relatively good prognosis: All but one sample respond favorably to chemotherapy and most (4/6) patients with CR are found in this group. Patients in this group also have a relatively long survival both after diagnosis (15.3±3.6 years) and after surgical resection of the tumor (4.8±1.5 years). Group III has the worst prognosis: None of the tumors respond to chemotherapy, the average time of survival after diagnosis was short (1.9±0.2 years) as was the average time after surgical resection (1.5±0.3 years). All tumors of this subgroup have retained both copies of Ip and 19q and are characterized by an amplification of the EGFR locus. The samples between groups II and III have a more mixed appearance, there is some degree of correlation with both groups I and group III. Many samples with PR and all samples with SD are found in this group. Survival after diagnosis and surgical resection is intermediate between groups II and III: 8.3±1.5 and 2.3±0.3 years respectively.
Supervised clustering: tumor vs. controls
We first performed supervised clustering to identify genes that are differentially expressed between control and tumors tissue. SAM analysis identified 1881 differentially expressed probesets (~1413 genes). Strongest downregulated transcripts in oligodendroglial tumors include those that encode proteins expressed in mature oligodendrocytes: myelin associated oligodendrocyte basic protein (MOBP), myelin oligodendrocyte glycoprotein (MOG), myelin associated glycoprotein (MAG), claudin 11 (CLDNIl) and myelin basic protein (MBP). These transcripts are expressed (±SD) at 0.052±0.021 (4 probesets), 0.10±0.013 (4 probesets), 0.086 (1 probeset), 0.30±0.25 (2 probesets), and 0.21±0.17 (7 probesets) levels of control brain mRNA respectively. This downregulation was observed in each sample. The strong downregulation in low-grade samples confirms the hypothesis that their homology to control brain tissue (see figure 1) is a result of the genes expressed
by the tumor. The downregulation of MOG was confirmed using RT-PCR (table 2).
It has been reported that PDGFRa is often highly expressed in oligodendroglial tumors (Riemenschneider M.J. et al., (2004) Acta Neuropathol. (Berlin) 107:277-282.). However, this gene was not present in the set of tumor-associated genes identified by our screen. Closer inspection reveals that, although PDGFRa is on average upregulated 4.1 fold, the high variation of upregulation (4.1±4.7) indicates that this transcript is not a reliable marker for the amount of tumor present in the sample. In fact, we failed to observe any upregulation in 10/28 samples. The select upregulation of PDGFRa in a subset of samples was confirmed using RTPCR.
Supervised clustering on chromosomal aberrations
Supervised clustering was performed to identify genes associated with specific chromosomal losses. For this we compared expression profiles of samples with i) Ip LOH (n=9) vs. no loss (n=9), ii) 19q LOH (n=ll) vs. no loss (n=7), and iii) combined Ip and 19q LOH (n=6) with no loss on either arm (n=6). SAM analysis identified 376, 64 and 60 probesets as being differentially expressed following loss of Ip, 19q or Ip and 19q respectively. Probesets are listed in supplementary table 3. Interestingly, many of the identified probesets are located on the lost chromosomal arm(s): 136/376 (36.1%) probesets are located on Ip, 25/64 (39.1%) on 19q and 49/60 (82%) on Ip or 19q. Of the differentially expressed genes located on the lost chromosomal arm(s), the ratio (±SD) loss vs. no loss is 0.53±0.22 (Ip), 0.54±0.07 (19q) and 0.53±0.09 (Ip and 19q) indicating that loss of one allele reduces expression levels by ~50%. In fact, all but two of the differentially expressed probesets that are located on the lost chromosomal(s) are downregulated. This correlation between chromosomal loss and expression level therefore suggest that these genes have an allele-number dependent expression level. Furthermore, the differentially expressed genes
can be identified across the entire chromosomal arms and suggests the entire arms have been lost.
Principle components analysis (PCA) and hierarchical clustering of genes associated with LOH on Ip and 19q is depicted in figure 2. All anaplastic oligodendrogliomas with combined loss/retention of Ip and 19q were correctly distributed by the first principal component axis, PCAl. This correct distribution includes 7 samples (2 samples that have retained both Ip and 19q copies and 5 samples with LOH on Ip and 19q) that were omitted from the clustering analysis. Further confirmation of a subset of differentially expressed genes by RT-PCR is shown in table 2 (including 4 additional oligodendroglial tumors).
Genes associated with chemosensitivity
We next performed supervised clustering to identify genes that are associated with response to chemotherapy. For this analysis we compared mRNA expression levels between tumors that show a response to chemotherapy (CR+PR), and those that do not (SD+PD). Such comparison using SAM (FDR<1 gene) identified 16 differentially expressed probesets that are listed in the supplementary table 3. 160 differentially expressed probesets (137 genes) were identified using less stringent statistical analysis (FDR=4.9%), of which 31 (27 genes) are located on chromosomes Ip or 19q (19%). Confirmation of differentially expressed genes was performed using RT-PCR on IQGAP, MANlCl, TRIM56 and AQPl transcripts (table 2). PCA based on the 16 genes associated with chemotherapeutic response identifies three main subgroups (figure 3): Samples with no response to chemotherapy (SD and PD, red), samples with response to treatment (CR and PR, green), and control samples (gray). Similarly, hierarchical clustering also separates the majority of oligodendroglial tumors with response to chemotherapy from those that show no or little response to treatment (figure 3). Similar results were obtained when clustering was performed on 160
differentially expressed probesets identified using FDR=4.9%. Most oligodendroglial tumors were correctly distributed on their response to treatment by the first principal component axis, PCAl: PCAl>0 in 14/18 samples that respond to treatment whereas PCAl<0 in 10/10 samples with no response to treatment. Only 4/28 samples were therefore incorrectly classified based on expression of genes associated with chemosensitivity. In comparison, 8/28 samples are incorrectly classified when predicting response to treatment based on the Ip chromosomal status: 6/14 tumors without LOH on Ip show response to treatment and 2/14 with LOH on Ip do not respond to treatment.
Genes associated with survival
We next performed supervised clustering to identify genes associated with overall survival after diagnosis. For this analysis we compared expression profiles of tumors from patients with the shortest survival time (2.0 ± 0.3 years, n=7) with those with the longest survival time (17.6 ± 4.4 years, n=8) after diagnosis. SAM analysis identified 103 probesets (92 genes, see supplementary data) associated with patient survival. 30 (29%) of these probesets are located on either Ip or 19q chromosomal arms. PCA of survival- associated genes identifies three main clusters of samples: oligodendroglial tumors with short survival, oligodendroglial tumors with long survival and control samples. Low-grade samples cluster between control and tumor samples. Similar subgroups were identified by hierarchical clustering using these probesets (figure 4). It is interesting to note that the subgroups identified by hierarchical clustering are virtually identical to the subgroups that were identified by unsupervised clustering (figure 1). Most oligodendroglial tumors were correctly distributed on survival after diagnosis by the first principal component axis, PCAl: PCAl>0 in 12/14 samples with favorable prognosis (i.e. survival time > 7 years after diagnosis) whereas PCAl<0 in 8/11 samples with relatively short survival after diagnosis (i.e. < 7 years).
Table 1. Summary of patient data, histological diagnosis and response to chemotherapy of samples used in this study. M: male; F: female; ctr: normal brain; ctr/w: control brain white matter; OD oligodendroglioma (grade II); AOD: anaplastic oligodendroglioma, AOA anaplastic oligoastrocytoma; LOH: loss of heterozygosity; ampl: amplification of the EGFR locus; ther.: therapy: PCV: combination therapy of procarbazine, CCNU, and vincristine; tenαo: temozolomide. Treatment response was scored according to McDonald's criteria (20) CR: complete response; PR: partial response; SD: stable disease; PD: progressive disease. Surv tot: patient survival after diagnosis (years); Surv op: patient survival after surgical resection of the sample used in this study.
47
Table 2. Confirmation of a subset of differentially expressed genes identified by expression profiling. Differential expression of most transcripts was reconfirmed by RT-PCR. The relative expression levels between control (either no loss of Ip, 19q, no tumor or CR/PR) and test set (either LOH on Ip, 19q, tumor or SD/PD) also remained similar on the array (rel expr array) and by RT-PCR (rel expr QPCR). QPCR ctr: expression of the examined transcript in control samples (either no loss of Ip, 19q, no tumor or CR/PR) relative to PDGB expression levels; QPCR marker: expression of the examined transcript in test samples (either LOH on Ip, 19q, no tumor or CR/PR) relative to PDGB expression levels. Statistical analysis was performed on QPCR ctr vs. marker using the Mann Whitney U test (two tailed), values are ±SE.
Table 3.: Differentially expressed probesets, which are able to discriminate on basis of response tot treatment
Table 4.: Differentially expressed probesets, which are able to discriminate on basis of patient survival
Table 5.: Differentially expressed probesets, which are able to discriminate on basis of loss of heterozygosity (LOH) on the Ip locus
Table 6.: Differentially expressed probesets, which are able to discriminate on basis of loss of heterozygosity (LOH) on the 19q locus
Table 7.: Differentially expressed probesets, which are able to discriminate on basis of loss of heterozygosity (LOH) on both the Ip and 19 q loci
Claims
1. A method for producing a classification scheme for oligodendroglial tumors comprising the steps of: a) providing a plurality of reference samples, said reference samples comprising cell samples from a plurality of reference subjects suffering from oligodendroglial tumors, with known responsiveness to therapy and survival or with known loss of heterozygosity of Ip and/or 19q; b) providing reference profiles by establishing a gene expression profile, matched with parameters for treatment sensitivity, survival and loss of heterozygosity for each of said reference samples individually; c) clustering said individual reference profiles according to a statistical procedure, comprising:
(i) K-means clustering; (ii) hierarchical clustering; and (iii) Pearson correlation coefficient analysis; and d) assigning an oligodendroglial tumor class according to treatment sensitivity and/or survival and/or loss of heterozygosity to each cluster.
2. Method according to claim 1, wherein the clustering of said gene expression profiles is performed based on the information of differentially - expressed genes and the treatment sensitivity and/or survival and/or loss of heterozygosity of the subject.
3. Method according to claim 1 or 2, wherein the clustering of said gene expression profiles with respect to treatment response is performed based on the information of the genes of Table 3.
4. Method according to claim 1 or 2, wherein the clustering of said gene expression profiles with respect to survival is performed based on the information of the genes of Table 4.
5. Method according to claim 1 or 2, wherein the clustering of said gene expression profiles with respect to loss of heterozygosity of Ip is performed based on the information of the genes of Table 5.
6. Method according to claim 1 or 2, wherein the clustering of said gene expression profiles with respect to loss of heterozygosity of 19q is based on the information of the genes of Table 6.
7. Method according to claim 1 or 2, wherein the clustering of said gene expression profiles with respect to loss of heterozygosity of Ip and 19q is performed based on the information of the genes of Table 7.
8. A method for classifying an oligodendroglial tumor of a subject suffering from oligodendroglioaml tumor, comprising the steps of: a) providing a classification scheme for oligodendroglial tumors according to the method of any one of claims 1-7; b) providing a subject profile by establishing a gene expression profile for said subject; c) clustering the subject profile together with reference profiles; d) determining in said scheme the clustered position of said subject profile among the reference profiles, and e) assigning to said oligodendroglial tumor the oligodendroglial tumor class that corresponds to said clustered position.
9. Method according to claim 8, wherein gene expression profile with respect to treatment response comprises the expression parameters of a set of genes according to table 3, still more preferably 1 to 50 genes of the genes of table 3.
10. Method according to claim 8, wherein gene expression profile with respect to survival comprises the expression parameters of a set of genes according to table 4, still more preferably 1 to 50 genes of the genes of table 4.
11. Method according to claim 8, wherein gene expression profile with respect to Ip loss of heterozygosity comprises the expression parameters of a set of genes according to table 5, still more preferably 1 to 50 genes of the genes of table 5.
12. Method according to claim 8, wherein gene expression profile with respect to 19q heterozygosity comprises the expression parameters of a set of genes according to table 6, still more preferably 1 to 50 genes of the genes of table 6.
13. Method according to claim 8, wherein gene expression profile with respect to Ip and 19q loss of heterozygosity comprises the expression parameters of a set of genes according to table 7, still more preferably 1 to 50 genes of the genes of table 7.
14. A method of determining the prognosis for a subject suffering from an oligodendroglial tumor, said method comprising the steps of: a) providing a classification scheme for oligodendroglial tumors by producing such a scheme according to the method of any one of claims 1-
7; b) determining the prognosis for each olidendroglial tumor class in said scheme based on clinical records for the subjects comprised in said class; c) establishing the oligodendroglial class of a subject suffering from an oligodendroglial tumor by classifying the oligodendroglial tumor in said subject according to a method of any of claims 8-13, and d) assigning to said subject the prognosis corresponding to the established oligodendroglial tumor class of said subject.
15. A method of determining the prognosis for a subject suffering from an oligodendroglial tumor, said method comprising the steps of: a) isolation of RNA from tumor cells of said subject; b) preparation of antisense, biotinylated RNA to the RNA of step a); c) hybridisation of said antisense, biotinylated DNA on Affymetrix U133A or U133 Plus2.0 GeneChips®; d) normalising the measured values for the gene set of Table 3; e) clustering the obtained data together with the reference data, obtained from a reference set of patient with known prognosis; and f) determining the prognosis on basis of the subgroup/cluster to which the data of the subject are clustering.
16. Oligonucleotide microarray of maximal 500 probesets, comprising at least 1, preferably at least 2, more preferably at least 25, still more preferably at least 100 oligonucleotide probes which each are capable of hybridizing under stringent conditions to different genes of the oligodendroglial tumor-associated genes selected from Table 3.
17. Oligonucleotide microarray of maximal 500 probesets, comprising at least 1, preferably at least 2, more preferably at least 25, still more preferably at least 100 oligonucleotide probes which each are capable of hybridizing under stringent conditions to different genes of the oligodendroglial tumor-associated genes selected from Table 4.
18. Oligonucleotide microarray of maximal 500 probesets, comprising at least 1, preferably at least 2, more preferably at least 25, still more preferably at least 100 oligonucleotide probes which each are capable of hybridizing under stringent conditions to different genes of the oligodendroglial tumor-associated genes selected from Table 5.
19. Oligonucleotide microarray of maximal 500 probesets, comprising at least 1, preferably at least 2, more preferably at least 25, still more preferably at least 100 oligonucleotide probes which each are capable of hybridizing under stringent conditions to different genes of the oligodendroglial tumor-associated genes selected from Table 6.
20. Oligonucleotide microarray of maximal 500 probesets, comprising at least 1, preferably at least 2, more preferably at least 25, still more preferably at least 100 oligonucleotide probes which each are capable of hybridizing under stringent conditions to different genes of the oligodendroglial tumor-associated genes selected from Table 7.
21. Kit-of-parts comprising an oligonucleotide microarray according to any of claims 16 to 20 and means for comparing a gene expression profile determined by using said microarray with a database of oligodendroglial tumor reference expression profiles.
Priority Applications (3)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US12/097,385 US20090215055A1 (en) | 2005-12-13 | 2005-12-13 | Genetic Brain Tumor Markers |
PCT/NL2005/000855 WO2007069882A1 (en) | 2005-12-13 | 2005-12-13 | Genetic brain tumor markers |
EP05822035A EP1969506A1 (en) | 2005-12-13 | 2005-12-13 | Genetic brain tumor markers |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
PCT/NL2005/000855 WO2007069882A1 (en) | 2005-12-13 | 2005-12-13 | Genetic brain tumor markers |
Publications (1)
Publication Number | Publication Date |
---|---|
WO2007069882A1 true WO2007069882A1 (en) | 2007-06-21 |
Family
ID=36046860
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
PCT/NL2005/000855 WO2007069882A1 (en) | 2005-12-13 | 2005-12-13 | Genetic brain tumor markers |
Country Status (3)
Country | Link |
---|---|
US (1) | US20090215055A1 (en) |
EP (1) | EP1969506A1 (en) |
WO (1) | WO2007069882A1 (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7794940B2 (en) | 2006-05-05 | 2010-09-14 | Mayo Foundation For Medical Education And Research | Detecting t(1;19) |
US10584144B2 (en) | 2015-03-09 | 2020-03-10 | University Of Kentucky Research Foundation | RNA nanoparticles for brain tumor treatment |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9852458B2 (en) | 2012-08-14 | 2017-12-26 | Ebay Inc. | Building containers of uncategorized items |
WO2014086949A2 (en) * | 2012-12-05 | 2014-06-12 | Royal College Of Surgeons In Ireland | System and method for the personalisation and optimization of anti-cancer treatments |
CN112634321B (en) * | 2020-10-21 | 2022-02-01 | 武汉大学 | Dam building particle material mechanical test system and method based on virtual reality combination |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005005601A2 (en) * | 2003-06-09 | 2005-01-20 | The Regents Of The University Of Michigan | Compositions and methods for treating and diagnosing cancer |
US20050209786A1 (en) * | 2003-12-11 | 2005-09-22 | Tzong-Hao Chen | Method of diagnosing biological states through the use of a centralized, adaptive model, and remote sample processing |
Family Cites Families (18)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5034506A (en) * | 1985-03-15 | 1991-07-23 | Anti-Gene Development Group | Uncharged morpholino-based polymers having achiral intersubunit linkages |
US5235033A (en) * | 1985-03-15 | 1993-08-10 | Anti-Gene Development Group | Alpha-morpholino ribonucleoside derivatives and polymers thereof |
US5143854A (en) * | 1989-06-07 | 1992-09-01 | Affymax Technologies N.V. | Large scale photolithographic solid phase synthesis of polypeptides and receptor binding screening thereof |
US6040138A (en) * | 1995-09-15 | 2000-03-21 | Affymetrix, Inc. | Expression monitoring by hybridization to high density oligonucleotide arrays |
US5744101A (en) * | 1989-06-07 | 1998-04-28 | Affymax Technologies N.V. | Photolabile nucleoside protecting groups |
US5424186A (en) * | 1989-06-07 | 1995-06-13 | Affymax Technologies N.V. | Very large scale immobilized polymer synthesis |
US5800992A (en) * | 1989-06-07 | 1998-09-01 | Fodor; Stephen P.A. | Method of detecting nucleic acids |
WO1993006121A1 (en) * | 1991-09-18 | 1993-04-01 | Affymax Technologies N.V. | Method of synthesizing diverse collections of oligomers |
EP0624059A4 (en) * | 1991-11-22 | 1994-12-21 | Affymax Technologies N.V. | Combinatorial strategies for polymer synthesis. |
US5556752A (en) * | 1994-10-24 | 1996-09-17 | Affymetrix, Inc. | Surface-bound, unimolecular, double-stranded DNA |
US5545531A (en) * | 1995-06-07 | 1996-08-13 | Affymax Technologies N.V. | Methods for making a device for concurrently processing multiple biological chip assays |
US5856174A (en) * | 1995-06-29 | 1999-01-05 | Affymetrix, Inc. | Integrated nucleic acid diagnostic device |
EP0880598A4 (en) * | 1996-01-23 | 2005-02-23 | Affymetrix Inc | Nucleic acid analysis techniques |
US5958905A (en) * | 1996-03-26 | 1999-09-28 | Texas Biotechnology Corporation | Phosphoramidates, phosphinic amides and related compounds and the use thereof to modulate the activity of endothelin |
WO1999005324A1 (en) * | 1997-07-25 | 1999-02-04 | Affymetrix, Inc. | System for providing a polymorphism database |
EP1027456B1 (en) * | 1997-10-31 | 2005-03-16 | Affymetrix, Inc. (a Delaware Corporation) | Expression profiles in adult and fetal organs |
JP2002502588A (en) * | 1998-02-06 | 2002-01-29 | アフィメトリックス インコーポレイテッド | Quality control methods in the manufacturing process |
US6020135A (en) * | 1998-03-27 | 2000-02-01 | Affymetrix, Inc. | P53-regulated genes |
-
2005
- 2005-12-13 EP EP05822035A patent/EP1969506A1/en not_active Withdrawn
- 2005-12-13 WO PCT/NL2005/000855 patent/WO2007069882A1/en active Application Filing
- 2005-12-13 US US12/097,385 patent/US20090215055A1/en not_active Abandoned
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2005005601A2 (en) * | 2003-06-09 | 2005-01-20 | The Regents Of The University Of Michigan | Compositions and methods for treating and diagnosing cancer |
US20050209786A1 (en) * | 2003-12-11 | 2005-09-22 | Tzong-Hao Chen | Method of diagnosing biological states through the use of a centralized, adaptive model, and remote sample processing |
Non-Patent Citations (3)
Title |
---|
KOHLMANN ALEXANDER ET AL: "Molecular characterization of acute leukemias by use of microarray technology.", GENES CHROMOSOMES AND CANCER, vol. 37, no. 4, August 2003 (2003-08-01), pages 396 - 405, XP002375048, ISSN: 1045-2257 * |
KUROKAWA YUKINORI ET AL: "Molecular prediction of response to 5-fluorouracil and interferon-alpha combination chemotherapy in advanced hepatocellular carcinoma.", CLINICAL CANCER RESEARCH, vol. 10, no. 18, 15 September 2004 (2004-09-15), pages 6029 - 6038, XP002374752, ISSN: 1078-0432 * |
MUKASA AKITAKE ET AL: "Distinction in gene expression profiles of oligodendrogliomas with and without allelic loss of 1p", ONCOGENE, vol. 21, no. 25, 6 June 2002 (2002-06-06), pages 3961 - 3968, XP002374750, ISSN: 0950-9232 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7794940B2 (en) | 2006-05-05 | 2010-09-14 | Mayo Foundation For Medical Education And Research | Detecting t(1;19) |
US10584144B2 (en) | 2015-03-09 | 2020-03-10 | University Of Kentucky Research Foundation | RNA nanoparticles for brain tumor treatment |
US11325939B2 (en) | 2015-03-09 | 2022-05-10 | University Of Kentucky Research Foundation | RNA nanoparticles for brain tumor treatment |
Also Published As
Publication number | Publication date |
---|---|
US20090215055A1 (en) | 2009-08-27 |
EP1969506A1 (en) | 2008-09-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP4980878B2 (en) | Classification, diagnosis, and prognosis of acute myeloid leukemia by gene expression profiling | |
JP6190434B2 (en) | Gene expression markers to predict response to chemotherapeutic agents | |
EP2768985B1 (en) | Colorectal cancer associated circulating nucleic acid biomarkers | |
EP2121988B1 (en) | Prostate cancer survival and recurrence | |
CN104364654A (en) | Methods for predicting and detecting cancer risk | |
JP2008521412A (en) | Lung cancer prognosis judging means | |
US20110159492A1 (en) | Diffuse Large B-Cell Lymphoma Markers and Uses Therefore | |
KR20140121524A (en) | Method for preparing prognosis prediction model of gastric cancer | |
WO2007069882A1 (en) | Genetic brain tumor markers | |
US20180051342A1 (en) | Prostate cancer survival and recurrence | |
WO2014110475A1 (en) | Biomarkers in cancer, methods, and systems related thereto | |
US7601532B2 (en) | Microarray for predicting the prognosis of neuroblastoma and method for predicting the prognosis of neuroblastoma | |
EP1683862B1 (en) | Microarray for assessing neuroblastoma prognosis and method of assessing neuroblastoma prognosis | |
AU2011204944B2 (en) | Gene expression markers for predicting response to chemotherapy | |
WO2010143941A1 (en) | Classification and risk-assignment of childhood acute myeloid leukaemia (aml) by gene expression signatures | |
US20110269636A1 (en) | Materials and methods for identifying patients at heighten risk for developing her2+ related brain tumors | |
CA3240150A1 (en) | Methods and materials for assessing homologous recombination deficiency in breast cancer subtypes | |
US20050142585A1 (en) | Determination of phenotype of cancer and of precancerous tissue | |
KR20140121523A (en) | Novel system for predicting prognosis of gastric cancer | |
JP2012511323A (en) | Method for identifying ERBB2 changes in tumors | |
WO2009151314A1 (en) | Classification and risk-assignment of childhood acute lymphoblastic leukaemia (all) by gene expression signatures |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
121 | Ep: the epo has been informed by wipo that ep was designated in this application | ||
NENP | Non-entry into the national phase |
Ref country code: DE |
|
WWE | Wipo information: entry into national phase |
Ref document number: 2005822035 Country of ref document: EP |
|
WWE | Wipo information: entry into national phase |
Ref document number: 12097385 Country of ref document: US |