Patent Grant
8065093
The application contains a sequence listing which has been submitted via CD-R in lieu of a printed paper copy, and is hereby incorporated by reference in its entirety. The CD-R, recorded on Mar. 9, 2005, are labeled “CRF,” “Copy 1,” and “Copy 2” and each contains only one identical 1.95 MB file (38271767.APP).
The present invention relates generally to systems, compositions, and methods for predicting disease susceptibility in a patient.
Mutations in p53 are thought to occur in more than 50% of human cancers and are most frequently observed in the DNA binding and transactivation domains, underscoring the importance of its transcriptional activity in suppressing tumor development. In sporadic breast cancers, unlike most cancer types, p53 mutations are only observed in approximately 20% of cases. However, that breast cancer is frequently observed in individuals with germline mutations of p53 (i.e., Li-Fraumeni syndrome) suggests a particularly important role for p53 inactivation in breast carcinogenesis, and perhaps a similarly important role for other factors capable of compromising p53 function.
For example, the reduced transcriptional activation of p53 following hypermethylation and subsequent inhibition of the HOXA5 transcription factor has recently been implicated as a possible epigenetic mechanism in reducing p53 expression in breast cancers. In both breast tumors and other cancer types, amplification and overexpression of the MDM2 gene, whose product promotes p53 degradation, has been implicated in oncogenesis. Moreover, both deletion and epigenetic silencing of the p14ARF gene, a negative regulator of MDM2, has been observed in various cancer types. Thus, p53 deficiency in breast carcinogenesis can potentially arise from a number of mechanisms other than p53 gene mutation.
There is evidence that the p53 status has prognostic significance in a number of cancer types and in particular breast cancer. In breast cancer, p53 mutations confer worse overall and disease-free survival, and a higher incidence of tumor recurrence, independent of other risk factors. Recent evidence suggests that p53 inactivation renders breast tumors resistant to certain DNA-damaging chemotherapies and endocrine therapies presumably through loss of p53-dependent apoptosis.
However, in all of these studies, the prognostic capability and degree of therapeutic resistance of the p53 mutants was found to depend largely on mutant-specific attributes, such as the type of mutations or the precise domain in which the mutation occurs. Importantly, this latter observation is consistent with findings from previous studies showing that not all p53 mutations have equal effects: some simply confer loss of function, while others have a dominant negative effect (such as trans-dominant suppression of wildtype p53 or oncogenic gain of function), while still others show only a partial loss of function where, for example, only a small subset of p53 downstream transcriptional target genes are dysregulated. For these reasons, no single molecular assessment of p53 status appears to provide an absolute indication of the complete p53 function.
There is a need for methods that better assess the effects of different p53 mutations on cell function in general and gene expression in particular, in an effort to enable better cancer prognosis and diagnosis.
Accordingly, the present invention provides methods, systems, and compositions that provide a more useful measure of in vivo p53 functionality. These methods, systems, and compositions may be employed for the classification, prognosis, and diagnosis of cancers.
In one aspect of the present invention there is provided a method for predicting disease outcome in a patient, the method comprising the steps of: obtaining gene expression profiles from a plurality of genes from tumor samples, wherein said tumor samples may be mutant or wildtype for the p53 gene; comparing said gene expression profiles to determine which genes are differentially expressed in the mutant or wildtype tumors; deriving from said differentially expressed genes a set of genes to predict p53 mutational status; and using the set of genes to predict disease outcome in the patient.
In another aspect of the present invention there is provided a method for predicting disease outcome in a late-stage breast cancer patient, the method comprising the steps of: obtaining gene expression profiles from a plurality of genes from tumor samples, wherein said tumor samples may be mutant or wildtype for the p53 gene; comparing said gene expression profiles to determine which genes are differentially expressed in the mutant or wildtype tumors; deriving from said differentially expressed genes a set of genes to predict p53 mutational status; and using the set of genes to predict disease outcome in the late-stage breast cancer patient wherein the set of genes are selected from the group consisting of GenBank accession numbers: BG271923 (SEQ ID NO: 22), NM—002466 (SEQ ID NO: 31), D38553 (SEQ ID NO: 11), NM—000909 (SEQ ID NO: 9), NM—024843 (SEQ ID NO: 1), R73030 (SEQ ID NO: 29), NM—003226 (SEQ ID NO: 28), AW299538 (SEQ ID NO: 5) and AI990465 (SEQ ID NO: 25).
In yet another aspect of the present invention there is provided a method for predicting clinical outcome in an early-stage, locally-treated breast cancer patient, the method comprising the steps of: obtaining gene expression profiles from a plurality of genes from tumor samples, wherein said tumor samples may be mutant or wildtype for the p53 gene; comparing said gene expression profiles to determine which genes are differentially expressed in the mutant or wildtype tumors; deriving from said differentially expressed genes a set of genes to predict p53 mutational status; and using the set of genes to predict disease outcome in the early-stage, locally-treated breast cancer patient wherein the set of genes are selected from the group consisting of GenBank accession numbers: AI961235 (SEQ ID NO-23), BG271923 (SEQ ID NO: 22), NM—002466 (SEQ ID NO: 31), BC001651 (SEQ ID NO: 14), D38553 (SEQ ID NO: 11), AK000345 (SEQ ID NO: 26), BC004504 (SEQ ID NO: 8), NM—000909 (SEQ ID NO: 9), NM—024843 (SEQ ID NO: 1), R73030 (SEQ ID NO: 29), AI435828 (SEQ ID NO: 20), AI810764 (SEQ ID NO: 24), AI922323 (SEQ ID NO: 10), NM—003225 (SEQ ID NO: 32), NM—003226 (SEQ ID NO: 28), AW299538 (SEQ ID NO: 5), NM—003462 (SEQ ID NO: 16), AI990465 (SEQ ID NO: 25), NM—004392 (SEQ ID NO: 15), NM—001267 (SEQ ID NO: 7) and AI826437 (SEQ ID NO: 3).
In a further aspect of the present invention there is provided a method for predicting clinical outcome in a liver cancer patient, the method comprising the steps of: obtaining gene expression profiles from a plurality of genes from tumor samples, wherein said tumor samples may be mutant or wildtype for the p53 gene; comparing said gene expression profiles to determine which genes are differentially expressed in the mutant or wildtype tumors; deriving from said differentially expressed genes a set of genes to predict p53 mutational status; and using the set of genes to predict disease outcome in the liver cancer patient wherein the set of genes are selected from the group consisting of GenBank accession numbers: NM—002466 (SEQ ID NO: 31), BC001651 (SEQ ID NO: 14), D38553 (SEQ ID NO: 11), NM—024843 (SEQ ID NO: 1), AI435828 (SEQ ID NO: 20), AI810764 (SEQ ID NO: 24), NM—003226 (SEQ ID NO: 28) and AW299538 (SEQ ID NO: 5).
In a still further aspect of the present invention there is provided a method of identifying a group of genes for predicting disease outcome in a patient, the method comprising the steps of: obtaining gene expression profiles from a plurality of genes from tumor samples, wherein said tumor samples may be mutant or wildtype for the p53 gene; comparing said gene expression profiles to determine which genes are differentially expressed in the mutant or wildtype tumors; ranking the differentially expressed genes according to their ability to predict p53 mutational status; training the ranked genes to distinguish between mutant and wildtype p53 gene expression profiles; obtaining a p53 classifier including a set of genes capable of predicting p53 mutational status; validating the p53 classifier in independent datasets; and assessing the ability of the p53 classifier to predict disease outcome in the patient.
In another aspect of the present invention there is provided a computer system for predicting disease outcome in a patient, the computer system comprising: a computer having a processor and a memory, the memory having executable code stored thereon for execution by the processor for performing the steps of: obtaining gene expression profiles from a plurality of genes from tumor samples, wherein said tumor samples may be mutant or wildtype for the p53 gene; comparing said gene expression profiles to determine which genes are differentially expressed in the mutant or wildtype tumors; deriving from said differentially expressed genes a set of genes to predict p53 mutational status; and using the set of genes to predict disease outcome in the patient.
In yet another aspect of the present invention there is provided a diagnostic tool for predicting disease susceptibility in a patient comprising a plurality of genes capable of predicting p53 mutational status immobilized on a solid support.
In a still further aspect of the present invention there is provided a nucleic acid array for predicting disease susceptibility in a patient comprising a solid support and displayed thereon nucleic acid probes corresponding to genes capable of predicting p53 mutational status in the patient.
These aspects and embodiments are described in greater detail below.
Definitions
As used in this application, the singular form “a,” “an,” and “the” include plural references unless the context clearly dictates otherwise. For example, the term “an agent” includes a plurality of agents, including mixtures thereof.
An individual is not limited to a human being but may also be other organisms including but not limited to a mammal, invertebrate, plant, fungus, virus, bacteria, or one or more cells derived from any of the above.
As used herein the term “comprising” means “including”. Variations of the word “comprising”, such as “comprise” and “comprises”, have correspondingly varied meanings.
Throughout this disclosure, various aspects of this invention can be presented in a range format. It should be understood that the description in range format is merely for convenience and brevity and should not be construed as an inflexible limitation on the scope of the invention. Accordingly, the description of a range should be considered to have specifically disclosed all the possible subranges as well as individual numerical values within that range. For example, description of a range such as from 1 to 6 should be considered to have specifically disclosed subranges such as from 1 to 3, from 1 to 4, from 1 to 5, from 2 to 4, from 2 to 6, from 3 to 6 etc., as well as individual numbers within that range, for example, 1, 2, 3, 4, 5, and 6. This applies regardless of the breadth of the range.
As used herein, the term “histologic grade” or “tumor grade” refers to characteristics of tumors classified according to the Elston-Ellis system of grading tumors.
As used herein, “p53 status” refers to the mutational status of the p53 gene. A p53 mutant tumor contains a mutation in the p53 gene that alters the function of the protein. A p53 wildtype tumor contains no detectable mutation in the p53 gene.
As used herein “Disease-specific survival” or DSS is a survival assessment where the end point being examined is death because of a disease, for example, breast cancer.
As used herein, “Disease-free survival” or DFS is a survival assessment where the end points are either tumor recurrence (i.e., the cancer comes back as the consequence of distant metastasis to other sites in the body) or death because of breast cancer without evidence of distant metastasis.
As used herein, an “array” is an intentionally created collection of molecules which can be prepared either synthetically or biosynthetically. The molecules in the array can be identical or different from each other. The array can assume a variety of formats, e.g., libraries of soluble molecules; libraries of compounds tethered to resin beads, silica chips, or other solid supports.
As used herein, a “nucleic acid library or array” is an intentionally created collection of nucleic acids which can be prepared either synthetically or biosynthetically in a variety of different formats (e.g., libraries of soluble molecules; and libraries of oligonucleotides tethered to resin beads, silica chips, or other solid supports). Additionally, the term “array” is meant to include those libraries of nucleic acids which can be prepared by spotting nucleic acids of essentially any length (e.g., from 1 to about 1000 nucleotide monomers in length) onto a substrate. The term “nucleic acid” as used herein refers to a polymeric form of nucleotides of any length, either ribonucleotides, deoxyribonucleotides or peptide nucleic acids (PNAs), that comprise purine and pyrimidine bases, or other natural, chemically or biochemically modified, non-natural, or derivatized nucleotide bases. The backbone of the polynucleotide can comprise sugars and phosphate groups, as may typically be found in RNA or DNA, or modified or substituted sugar or phosphate groups. A polynucleotide may comprise modified nucleotides, such as methylated nucleotides and nucleotide analogs. The sequence of nucleotides may be interrupted by non-nucleotide components. Thus the terms nucleoside, nucleotide, deoxynucleoside and deoxynucleotide generally include analogs such as those described herein. These analogs are those molecules having some structural features in common with a naturally occurring nucleoside or nucleotide such that when incorporated into a nucleic acid or oligonucleotide sequence, they allow hybridization with a naturally occurring nucleic acid sequence in solution. Typically, these analogs are derived from naturally occurring nucleosides and nucleotides by replacing and/or modifying the base, the ribose or the phosphodiester moiety. The changes can be tailor made to stabilize or destabilize hybrid formation or enhance the specificity of hybridization with a complementary nucleic acid sequence as desired.
As used herein, the term “complementary” refers to the hybridization or base pairing between nucleotides or nucleic acids, such as, for instance, between the two strands of a double stranded DNA molecule or between an oligonucleotide primer and a primer binding site on a single stranded nucleic acid to be sequenced or amplified. Complementary nucleotides are, generally, A and T (or A and U), or C and G. Two single stranded RNA or DNA molecules are said to be complementary when the nucleotides of one strand, optimally aligned and compared and with appropriate nucleotide insertions or deletions, pair with at least about 80% of the nucleotides of the other strand, usually at least about 90% to 95%, and more preferably from about 98 to 100% of the nucleotides of the other strand. Alternatively, complementarity exists when an RNA or DNA strand will hybridize under selective hybridization conditions to its complement. Typically, selective hybridization will occur when there is at least about 65% complementarity over a stretch of at least 14 to 25 nucleotides, preferably at least about 75%, and more preferably at least about 90% complementarity.
As used herein, a “fragment,” “segment,” or “DNA segment” refers to a portion of a larger DNA polynucleotide or DNA. A polynucleotide, for example, can be broken up, or fragmented into, a plurality of segments. Various methods of fragmenting nucleic acids are well known in the art. These methods may be, for example, either chemical or physical in nature. Chemical fragmentation may include partial degradation with a DNase; partial depurination with acid; the use of restriction enzymes; intron-encoded endonucleases; DNA-based cleavage methods, such as triplex and hybrid formation methods, that rely on the specific hybridization of a nucleic acid segment to localize a cleavage agent to a specific location in the nucleic acid molecule; or other enzymes or compounds which cleave DNA at known or unknown locations. Physical fragmentation methods may involve subjecting the DNA to a high shear rate. High shear rates may be produced, for example, by moving DNA through a chamber or channel with pits or spikes, or forcing the DNA sample through a restricted size flow passage, e.g., an aperture having a cross sectional dimension in the micron or submicron scale. Other physical methods include sonication and nebulization. Combinations of physical and chemical fragmentation methods may likewise be employed such as fragmentation by heat and ion-mediated hydrolysis. See for example, Sambrook et al., “Molecular Cloning: A Laboratory Manual,” 3rd Ed. Cold Spring Harbor Laboratory Press, Cold Spring Harbor, N.Y. (2001) (“Sambrook et al.) which is incorporated herein by reference for all purposes. These methods can be optimized to digest a nucleic acid into fragments of a selected size range. Useful size ranges may be from 100, 200, 400, 700 or 1000 to 500, 800, 1500, 2000, 4000 or 10,000 base pairs. However, larger size ranges such as 4000, 10,000 or 20,000 to 10,000, 20,000 or 500,000 base pairs may also be useful.
As used herein, the term “hybridization” refers to the process in which two single-stranded polynucleotides bind non-covalently to form a stable double-stranded polynucleotide. The term “hybridization” may also refer to triple-stranded hybridization. The resulting (usually) double-stranded polynucleotide is a “hybrid.” The proportion of the population of polynucleotides that forms stable hybrids is referred to herein as the “degree of hybridization”. Hybridization conditions will typically include salt concentrations of less than about 1M, more usually less than about 500 mM and less than about 200 mM. Hybridization temperatures can be as low as 5° C., but are typically greater than 22° C., more typically greater than about 30° C., and preferably in excess of about 37° C. Hybridizations are usually performed under stringent conditions, i.e. conditions under which a probe will hybridize to its target subsequence. Stringent conditions are sequence-dependent and are different in different circumstances. Longer fragments may require higher hybridization temperatures for specific hybridization. As other factors may affect the stringency of hybridization, including base composition and length of the complementary strands, presence of organic solvents and extent of base mismatching, the combination of parameters is more important than the absolute measure of any one alone. Generally, stringent conditions are selected to be about 5° C. lower than the thermal melting point (Tm) for the specific sequence at a defined ionic strength and pH. The Tm is the temperature (under defined ionic strength, pH and nucleic acid composition) at which 50% of the probes complementary to the target sequence hybridize to the target sequence at equilibrium.
Typically, stringent conditions include salt concentration of at least 0.01 M to no more than 1 M Na ion concentration (or other salts) at a pH 7.0 to 8.3 and a temperature of at least 25° C. For example, conditions of 5× SSPE (750 mM NaCl, 50 mM NaPhosphate, 5 mM EDTA, pH 7.4) and a temperature of 25-30° C. are suitable for allele-specific probe hybridizations. For stringent conditions, see for example, Sambrook, Fritsche and Maniatis. “Molecular Cloning A Laboratory Manual” 2nd Ed. Cold Spring Harbor Press (1989) and Anderson “Nucleic Acid Hybridization” 1st Ed., BIOS Scientific Publishers Limited (1999), which are hereby incorporated by reference in their entireties for all purposes above.
As used herein, “hybridization probes” are nucleic acids (such as oligonucleotides) capable of binding in a base-specific manner to a complementary strand of nucleic acid. Such probes include peptide nucleic acids, as described in Nielsen et al., Science 254:1497-1500 (1991), Nielsen Curr. Opin. Biotechnol., 10:71-75 (1999) and other nucleic acid analogs and nucleic acid mimetics.
As used herein, “mRNA” or “mRNA transcripts” include, but are not limited to pre-mRNA transcript(s), transcript processing intermediates, mature mRNA(s) ready for translation and transcripts of the gene or genes, or nucleic acids derived from the mRNA transcript(s). Transcript processing may include splicing, editing and degradation. As used herein, a nucleic acid derived from an mRNA transcript refers to a nucleic acid for whose synthesis the mRNA transcript or a subsequence thereof has ultimately served as a template. Thus, a cDNA reverse transcribed from an mRNA, a cRNA transcribed from that cDNA, a DNA amplified from the cDNA, an RNA transcribed from the amplified DNA, etc., are all derived from the mRNA transcript and detection of such derived products is indicative of the presence and/or abundance of the original transcript in a sample. Thus, mRNA derived samples include, but are not limited to, mRNA transcripts of the gene or genes, cDNA reverse transcribed from the mRNA, cRNA transcribed from the cDNA, DNA amplified from the genes, RNA transcribed from amplified DNA, and the like.
As used herein, a “probe” is a molecule that can be recognized by a particular target. In some embodiments, a probe can be surface immobilized. Examples of probes that can be investigated by this invention include, but are not restricted to, agonists and antagonists for cell membrane receptors, toxins and venoms, viral epitopes, hormones (e.g. opioid peptides, steroids, etc.), hormone receptors, peptides, enzymes, enzyme substrates, cofactors, drugs, lectins, sugars, oligonucleotides, nucleic acids, oligosaccharides, proteins, and monoclonal antibodies.
As used herein, a “target” is a molecule that has an affinity for a given probe. Targets may be naturally-occurring or man-made molecules. Also, they can be employed in their unaltered state or as aggregates with other species. Targets may be attached, covalently or noncovalently, to a binding member, either directly or via a specific binding substance. Examples of targets which can be employed by this invention include, but are not restricted to, antibodies, cell membrane receptors, monoclonal antibodies and antisera reactive with specific antigenic determinants (such as on viruses, cells or other materials), drugs, oligonucleotides, nucleic acids, peptides, cofactors, lectins, sugars, polysaccharides, cells, cellular membranes, and organelles. Targets are sometimes referred to in the art as anti-probes. As the term targets is used herein, no difference in meaning is intended. A “Probe Target Pair” is formed when two macromolecules have combined through molecular recognition to form a complex.
The patent or application file contains at least one drawing executed in color. Copes of this patent or patent application publication with color drawing(s) will be provided by the Office upon request and payment of the necessary fee.
The accompanying drawings, which are incorporated in and form a part of this specification, illustrate embodiments of the invention and, together with the description, serve to explain the disclosed principles of the invention:
The present invention has many preferred embodiments and relies on many patents, applications and other references for details known to those of the art. Therefore, when a patent, application, or other reference is cited or repeated below, it should be understood that it is incorporated by reference in its entirety for all purposes as well as for the proposition that is recited.
Embodiments of the disclosed methods, systems, and compositions for classification, prognosis, and diagnosis of cancers will now be described. These methods, systems, and compositions provide a more useful measure of in vivo p53 functionality and thereby provide a better prognostic indicator of patient outcome as compared to p53 mutation status alone. Other advantages inherent in the disclosed embodiments of the methods, systems, and compositions will be apparent from the following description.
p53 mutations in cancer development and progression can result in trans-dominant suppression of the wild-type p53 allele conferring loss of p53 activity or an oncogenic gain of function independent of wildtype p53. Additionally, the altered activity of some effectors of p53 function, including those that directly influence p53 expression, may contribute to p53 deficiency recapitulating the p53-mutant phenotype. In breast cancer, these effects manifest in more aggressive tumors, therapeutic resistance, and poor clinical outcome.
In accordance with providing a more useful measure of in vivo p53 functionality, disclosed herein is a “p53 classifier”, an expression signature deduced from differences in the molecular configurations of p53 wildtype and mutant tumors. The classifier may comprise a defined number of genes, for example, at least 3 genes. In other embodiments, the classifier may comprise from about 3 genes to about 500 genes. Table 1 provides a listing of the 500 genes. In some embodiments, an optimized p53 classifier comprises 32 genes (Table 2). The optimized 32-gene classifier could distinguish p53 mutant and wildtype tumors with significant accuracy and could predict recurrence and survival in populations representing all therapeutic groups. Moreover, the p53 classifier was a more significant predictor of survival than p53 mutation status alone and remained significant by multivariate analysis independent of other clinical predictors where p53 mutation status did not. Furthermore, downregulation of p53 expression in the absence of mutations was sufficient to induce a mutant (mt) phenotype tumor behaviour in both transcriptional activity and clinical outcome.
In independent datasets of both breast and liver cancers, and regardless of other clinical features, subsets of the optimized p53 classifier could predict p53 status with significant accuracy. As a predictor of disease-specific survival (DSS), the classifier significantly outperformed p53 mutational status alone in both a large patient cohort with heterogeneous treatment, as well as in a set of patients who received postoperative adjuvant endocrine therapy alone.
Moreover, in an independent cDNA microarray study comprised mostly of stage 3 patients who received chemotherapy in the neoadjuvant setting, a 9-gene subset of the p53 classifier was a highly significant predictor of both disease-specific and disease-free survival. The genes of the p53 classifier could accurately discern not only which patients would relapse and die following chemotherapy, but also which late stage patients would survive their cancer.
A 21-gene subset of the classifier could also significantly distinguish molecular subgroups of early-stage radiation-treated patients who would go on to develop a distant metastasis within 5 years from those who would not.
Therefore, by defining among other aspects, a p53 classifier described herein, the methods, systems and compositions of the present invention demonstrate a much greater impact of p53 on human tumor behaviour than previously appreciated and thereby provide a better approach for clinically assessing p53 function.
One aspect of the present invention provides a method for predicting disease outcome in a patient, the method comprising the steps of: obtaining gene expression profiles from a plurality of genes from tumor samples, wherein said tumor samples may be mutant or wildtype for the p53 gene; comparing said gene expression profiles to determine which genes are differentially expressed in the mutant or wildtype tumors; deriving from said differentially expressed genes a set of genes to predict p53 mutational status; and using the set of genes to predict disease outcome in the patient. The disease outcome may be selected from the group consisting of disease-specific survival, disease-free survival, tumor recurrence and therapeutic response. The disease may be any cancer but is preferably breast cancer or liver cancer.
The predicted p53 mutational status may be obtained by ranking the differentially expressed genes according to their association with p53 mutational status, ER (estrogen receptor) status and histologic grade of the tumor. A multivariate ranking procedure such as a Linear Model Fit may be employed to rank the genes. The ranked genes may be subjected to supervised learning to enable them to distinguish between mutant and wildtype gene expression profiles. An example of a supervised learning method that may be employed is Diagonal Linear Discriminant Analysis (DLDA).
In some embodiments, the set of genes with the ability to predict p53 mutational status may comprise at least 3 genes, preferably about 3-500 genes and most preferably about 32 genes. The 32 genes making up the optimized p53 classifier may be selected from the group comprising the list of genes in Table 1. In some embodiments, the 32 genes may include GenBank accession numbers: AI961235 (SEQ ID NO: 23), BG271923 (SEQ ID NO: 22), NM—002466 (SEQ ID NO: 31), BC001651 (SEQ ID NO: 14), D38553 (SEQ ID NO: 11), AK000345 (SEQ ID NO: 26), AA742697 (SEQ ID NO: 21), AL080170 (SEQ ID NO: 30), BF245284 (SEQ ID NO: 27), BC004504 (SEQ ID NO: 8), H15261 (SEQ ID NO: 2), NM—000909 (SEQ ID NO: 9), NM—024843 (SEQ ID NO: 1), R73030 (SEQ ID NO: 29), NM—030896 (SEQ ID NO: 17), AI435828 (SEQ ID NO: 20), AL512727 (SEQ ID NO: 6), AW242997 (SEQ ID NO: 18), AI810764 (SEQ ID NO: 24), AI922323 (SEQ ID NO: 10), AL360204 (SEQ ID NO: 13), NM—003225 (SEQ ID NO: 32), NM—003226 (SEQ ID NO: 28), AW299538 (SEQ ID NO: 5), NM—003462 (SEQ ID NO: 16), AI990465 (SEQ ID NO: 25), NM—004392 (SEQ ID NO: 15), NM—001267 (SEQ ID NO: 7), AF269087 (SEQ ID NO: 4), AI826437 (SEQ ID NO: 3), AL355392 (SEQ ID NO: 12), and AU156421 (SEQ ID NO: 19).
The present invention also provides a method for predicting disease outcome in a late-stage breast cancer patient, the method comprising the steps of obtaining gene expression profiles from a plurality of genes from tumor samples, wherein said tumor samples may be mutant or wildtype for the p53 gene; comparing said gene expression profiles to determine which genes are differentially expressed in the mutant or wildtype tumors; deriving from said differentially expressed genes a set of genes to predict p53 mutational status; and using the set of genes to predict disease outcome in the late-stage breast cancer patient wherein the set of genes are selected from the group consisting of GenBank accession numbers: BG271923, NM—002466, D38553, NM—000909, NM—024843, R73030, NM—003226, AW299538 and AI990465. All GenBank accession numbers are associated with a sequence and a SEQ ID NO. as shown in
The present invention also provides a method for predicting clinical outcome in an early-stage, locally-treated breast cancer patient, the method comprising the steps of: obtaining gene expression profiles from a plurality of genes from tumor samples, wherein said tumor samples may be mutant or wildtype for the p53 gene; comparing said gene expression profiles to determine which genes are differentially expressed in the mutant or wildtype tumors; deriving from said differentially expressed genes a set of genes to predict p53 mutational status; and using the set of genes to predict disease outcome in the early-stage, locally-treated breast cancer patient wherein the set of genes are selected from the group consisting of GenBank accession numbers: AI961235, BG271923, NM—002466, BC001651, D38553, AK000345, BC004504, NM—000909, NM—024843, R73030, AI435828, AI810764, AI922323, NM—003225, NM—003226, AW299538, NM—003462, AI990465, NM—004392, NM—001267 and AI826437.
The present invention also provides a method for predicting clinical outcome in a liver cancer patient, the method comprising the steps of: obtaining gene expression profiles from a plurality of genes from tumor samples, wherein said tumor samples may be mutant or wildtype for the p53 gene; comparing said gene expression profiles to determine which genes are differentially expressed in the mutant or wildtype tumors; deriving from said differentially expressed genes a set of genes to predict p53 mutational status; and using the set of genes to predict disease outcome in the liver cancer patient wherein the set of genes are selected from the group consisting of GenBank accession numbers: NM—002466, BC001651, D38553, NM—024843, AI435828, AI810764, NM—003226 and AW299538.
The present invention also provides a method of identifying a group of genes for predicting disease outcome in a patient, the method comprising the steps of: obtaining gene expression profiles from a plurality of genes from tumor samples, wherein said tumor samples may be mutant or wildtype for the p53 gene; comparing said gene expression profiles to determine which genes are differentially expressed in the mutant or wildtype tumors; ranking the differentially expressed genes according to their ability to predict p53 mutational status; training the ranked genes to distinguish between mutant and wildtype p53 gene expression profiles; obtaining a p53 classifier including a set of genes capable of predicting p53 mutational status; validating the p53 classifier in independent datasets; and assessing the ability of the p53 classifier to predict disease outcome in the patient.
In the above-disclosed method of identifying a group of genes for predicting disease outcome in a patient, the differentially expressed genes may be ranked by a multivariate ranking procedure according to their association with p53 status, ER (estrogen receptor) status and histologic grade of the tumor. The multivariate ranking procedure may be a Linear Model-Fit method or any other method known to one of skill in the art. The step of training may comprise employing a supervised learning method, such as Diagonal Linear Discriminant Analysis (DLDA) or any other supervised learning method known to one of skill in the art.
The p53 classifier disclosed above may comprise at least 3 genes, preferably between about 3-500 genes and more preferably about 32 genes. This 32-gene p53 classifier is an “optimized classifier” which may include genes selected from the group consisting of GenBank accession numbers: AI961235, BG271923, NM—002466, BC001651, D38553, AK000345, AA742697, AL080170, BF245284, BC004504, H15261, NM—000909, NM—024843, R73030, NM—030896, AI435828, AL512727, AW242997, AI810764, AI922323, AL360204, NM—003225, NM—003226, AW299538, NM—003462, AI990465, NM—004392, NM—001267, AF269087, AI826437, AL355392 and AU156421.
The disease outcome may be selected from the group consisting of disease-specific survival, disease-free survival, tumor recurrence and therapeutic response. In one disclosed embodiment, a 9-gene partial classifier may predict clinical outcome in a late-stage breast cancer patient. The 9-gene partial classifier may include genes selected from the group consisting of GenBank accession numbers: BG271923, NM—002466, D38553, NM—000909, NM—024843, R73030, NM—003226, AW299538 and AI990465.
In another disclosed embodiment, a 21-gene partial classifier may predict clinical outcome in an early-stage, locally-treated breast cancer patient. The 21-gene partial classifier may include genes selected from the group consisting of GenBank accession numbers: AI961235, BG271923, NM—002466, BC001651, D38553, AK000345, BC004504, NM—000909, NM—024843, R73030, AI435828, AI810764, AI922323, NM—003225, NM—003226, AW299538, NM—003462, AI990465, NM—004392, NM—001267 and AI826437.
In yet another disclosed embodiment, a 8-gene partial classifier may predict clinical outcome in a liver cancer patient. The 8-gene partial classifier may include genes selected from the group consisting of GenBank accession numbers: NM—002466, BC001651, D38553, NM—024843, AI435828, AI810764, NM—003226 and AW299538.
The present invention also provides a computer system for predicting disease outcome in a patient, the computer system comprising: a computer having a processor and a memory, the memory having executable code stored thereon for execution by the processor for performing the steps of: obtaining gene expression profiles from a plurality of genes from tumor samples, wherein said tumor samples may be mutant or wildtype for the p53 gene; comparing said gene expression profiles to determine which genes are differentially expressed in the mutant or wildtype tumors; deriving from said differentially expressed genes a set of genes to predict p53 mutational status; and using the set of genes to predict disease outcome in the patient.
The present invention also provides a diagnostic tool for predicting disease susceptibility in a patient comprising a plurality of genes capable of predicting p53 mutational status immobilized on a solid support. The solid support may be a microarray, for example. In one embodiment, the plurality of genes immobilized on the solid support may include genes selected from the group consisting of GenBank accession numbers: AI961235, BG271923, NM—002466, BC001651, D38553, AK000345, AA742697, AL080170, BF245284, BC004504, H15261, NM—000909, NM—024843, R73030, NM—030896, AI435828, AL512727, AW242997, AI810764, AI922323, AL360204, NM—003225, NM—003226, AW299538, NM—003462, AI990465, NM—004392, NM—001267, AF269087, AI826437, AL355392 and AU156421. In another embodiment, the plurality of genes immobilized on the solid support may include genes selected from the group consisting of GenBank accession numbers: BG271923, NM—002466, D38553, NM—000909, NM—024843, R73030, NM—003226, AW299538 and AI990465. In yet another embodiment, the plurality of genes immobilized on the solid support may include genes selected from the group consisting of GenBank accession numbers: AI961235, BG271923, NM—002466, BC001651, D38553, AK000345, BC004504, NM—000909, NM—024843, R73030, AI435828, AI810764, AI1922323, NM—003225, NM—003226, AW299538, NM—003462, AI990465, NM—004392, NM—001267 and AI826437. In a still further embodiment, the plurality of genes immobilized on the solid support may include genes selected from the group consisting of GenBank accession numbers: NM—002466, BC001651, D38553, NM—024843, AI435828, AI810764, NM—003226 and AW299538.
The present invention also provides a nucleic acid array for predicting disease susceptibility in a patient comprising a solid support and displayed thereon nucleic acid probes corresponding to genes capable of predicting p53 mutational status in the patient. The nucleic acid array may comprise at least 8, 32, 100, 250 or 500 nucleic acid probes.
Thus, the disclosed methods, systems and compositions are capable of discerning p53-deficient from p53-enabled breast tumors and may be effective in gauging p53 activity in other cancer types. As much as 14% of breast tumors that are otherwise p53 wildtype at the DNA sequence level may be deficient for p53 by other means. Moreover, the classifier is a significant predictor of disease-specific survival and recurrence in various breast cancer populations and therefore will have clinical utility in predicting these endpoints, particularly in the context of therapeutic agents that function predominantly through p53-dependent cell death pathways.
To gain insight into the molecular variation between p53 mutant (mt) and p53 wildtype (wt) breast tumors, high-density oligonucleotide microarrays were utilized to analyze a population-based series of 257 biopsies, all of which were previously sequenced for mutations in the p53 coding regions (Bergh, J., Norberg, T., Sjogren, S., Lindgren, A. & Holmberg, L. Complete sequencing of the p53 gene provides prognostic information in breast cancer patients, particularly in relation to adjuvant systemic therapy and radiotherapy. Nat Med 1, 1029-34 (1995), incorporated herein by reference).
The original patient material consisted of freshly frozen breast tumors from a population-based cohort of 315 women representing 65% of all breast cancers resected in Uppsala County during the time period Jan. 1, 1987 to Dec. 31, 1989 (Bergh et al., previously incorporated by reference). After surgery, the viable part of the fresh tumor was cut in two; one part was immediately frozen in isopentane and stored at −70° C. until analysis, and the other was fixed in 10% formalin and prepared for histopathologic examination. Frozen tumor tissue was available from 299 of the original 315 patients. Out of these, 270 had RNA of sufficient quantity and quality for microarray experiments, and after Affymetrix quality control, expression profiles of 260 tumors were further analysed. The present study was approved by the ethical committee at the Karolinska Institute.
Mutational analysis of the p53 gene (TP53) was carried out in the original 315 tumors as described previously in Bergh et al. (previously incorporated by reference). Among the 260 tumors included in the present study, 59 had p53 mutations found by cDNA sequence analysis of exons 2 to 11 (Bergh et al., previously incorporated by reference). In three samples p53 status could not be evaluated. Clinico-pathological characteristics were derived from the patient records and from routine clinical measurements at the time of diagnosis. Estrogen receptor status was determined by ligand binding assay as part of the routine clinical procedure. An experienced pathologist determined the Elston-Ellis grades of the tumors, classifying the tumors into low, medium and high-grade tumors (Elston, C. W. & Ellis, I. O. Pathological prognostic factors in breast cancer. I. The value of histological grade in breast cancer: experience from a large study with long-term follow-up. Histopathology 19, 403-10 (1991), incorporated herein by reference). Axillary lymph node metastases were found in 84 of these 260 patients while 166 were node-negative. Ten patients had unknown node status, as no axillary examination was performed due to advanced age or concomitant serious disease. Systemic adjuvant therapy was offered to all node-positive patients. In general, premenopausal women were offered chemotherapy and postmenopausal women received endocrine treatment. Out of the 260 patients included in the present study, 149 did not receive adjuvant therapy. Overall survival of the patients was based on information from the Swedish population registry, and date and cause of death were obtained from a review of the patient records in late 1999.
RNA from 59 tumors known to contain p53 mutations resulting in amino acid-level alterations, and from 198 tumors known to have wildtype p53 were analyzed on Affymetrix U133A and U133B arrays.
Extraction of total RNA was carried out using the Qiagen RNeasy Mini Kit (Qiagen, Germany). Frozen tumors were cut into small pieces and homogenized for around 30-40 seconds in test tubes (maximum 40 mg/tube) containing RLT buffer (RNeasy lysis buffer) with mercaptoethanol. The mixtures were then treated with Proteinase K for 10 minutes at 55° C., which in previous RNA extractions demonstrated improved RNA yield (Egyhazi, S. et al. Proteinase K added to the extraction procedure markedly increases RNA yield from primary breast tumors for use in microarray studies. Clin Chem 50, 975-6 (2004), incorporated herein by reference). In the following centrifugation steps on RNeasy columns, DNase treatment was also included to increase the RNA quality. The integrity of the RNA extracts was tested on an Agilent 2100 Bioanalyzer (Agilent Technologies, Rockville, Md., U.S.A), measuring the 28S:18S ribosomal RNA ratio. RNA extracts of high quality were stored at −70° C. until microarray analysis.
Preparation of in vitro transcription (IVT) products (i.e., target) and oligonucleotide array hybridization and scanning were performed according to the Affymetrix protocol (Affymetrix Inc., Santa Clara, Calif., U.S.A). First-strand cDNA was synthesized from a starting amount of 2-5 μg total RNA using a T7-linked oligo-dT primer, followed by second-strand synthesis. Double-stranded cDNA was purified using phenol/chloroform extraction and phase lock gel. Biotinylated cRNA targets were prepared from the cDNA templates in IVT reactions. The labeled cRNA targets were purified using Qiagen RNeasy Mini Kit and subsequently chemically fragmented. Ten μg of the fragmented, biotinylated cRNA was hybridized to the Affymetrix oligonucleotide human array set, HG-U133A&B, which contains 45,000 probe sets representing more than 39,000 transcripts derived from approximately 33,000 well-substantiated human genes. Hybridization was carried out in a hybridization oven at 45° C. and rotation was set at 60 rpm for 16 h. The arrays were washed and stained in the Fluidics Station 400 (Affymetrix Inc., Santa Clara, Calif., U.S.A) in accordance with the Affymetrix protocol. Staining was carried out using streptavidin-phycoerythrin (SAPE, final concentration of 10 μg/ml) and signal amplification with a biotinylated anti-streptavidin antibody and a second SAPE staining. The arrays were washed and scanned according to the manufacturer's instructions.
The raw expression data was processed using Microarray Suite 5.0 software (Affymetrix Inc., Santa Clara, Calif., U.S.A) and normalized using the global mean method. For each microarray, probeset signal values were scaled by adjusting the mean log intensity to a target signal value of 500. Samples with suboptimal average signal intensities were re-labeled and re-hybridized on new arrays. If microarray artifacts were visible, the samples were re-hybridized on new chips using the same fragmented probe, or alternatively, if the defective areas were small, the affected probes were censored from further analysis. The normalized expression data from both U133A and B chips were combined and natural log transformed.
The extent to which gene expression patterns could distinguish p53 mt and wt tumors was first investigated. By Wilcoxon rank-sum test 3,330 Affymetrix probe-sets representing ˜2,770 distinct genes (according to UniGene build #167) were identified whose expression patterns distinguished p53 mt and wt tumors with a false discovery rate (FDR)-adjusted p value of p<0.001. A number of these genes were found to be known transcriptional targets of p53 including PERP, RRM2, SEMA3B, TAP1, GTSE1, CHECK1, and CHEK2. Shown in
The hierarchical structure of the gene expression profiles was next investigated. As in the tumors, two predominant clusters were observed: one consisting of ˜200 genes more highly expressed in the mutant-like tumor cluster, and the other representing ˜50 genes more highly expressed in the wildtype-like cluster. Within the former, the genes most highly correlated with p53 mutant status were associated with cell cycle progression including, CDC2, CDC20, CCNB1, CCNB2, CKS2, CDCA1, CDCA3, CDCA8, CENPA, TOP2A, PTTG1 and MCM6. This finding is consistent with the observation that wt p53 has a negative regulatory effect on cell cycle genes. Of the genes more highly expressed in the wildtype-like cluster, the presence of several estrogen-regulated and ER status-associated genes including STC2, NCOR1, and ADRA2A was observed.
Further examination of the tumors revealed that in addition to p53 status, the predominant tumor clusters were also correlated with other clinical features, namely estrogen receptor (ER) status and tumor grade. The estrogen receptor status of a cell has been found to be correlated with cancer in several instances. Normal breast cells usually have receptors for estrogen. However, cancer cells arising in the breast do not always have receptors for estrogen. Breast cancers that have estrogen receptors are said to be “estrogen receptor-positive,” while those breast cancers that do not possess estrogen receptors are “estrogen receptor-negative.” In estrogen receptor-positive cancers, cancer cell growth is under the control of estrogen. In contrast, the growth of estrogen receptor-negative cancer cells is not governed by estrogen.
Segregating with the mutant-like cluster were observed 86% of estrogen receptor-negative (ER−) tumors (pcs=1 7×10−10), 96% of grade III tumors (pcs=2.5×10−19) and only 3% of grade I tumors (pfe=6.9×10−15). This result owes, in part, to the fact that the p53 mutants in this study are positively correlated with ER negativity (pcs=1.7×10−6) and grade III status (pcs=1.2×10−11), and is consistent with previous reports demonstrating that p53 mutant breast cancers are significantly correlated with negative ER status and higher tumor grade. See for example, Cattoretti, G., Rilke, F., Andreola, S., D'Amato, L. & Delia, D. P53 expression in breast cancer. Int J Cancer 41, 178-83 (1988); Isola, J., Visakorpi, T., Holli, K. & Kallioniemi, O. P. Association of overexpression of tumor suppressor protein p53 with rapid cell proliferation and poor prognosis in node-negative breast cancer patients. J Natl Cancer Inst 84, 1109-14 (1992); Andersen, T. I. et al. Prognostic significance of TP53 alterations in breast carcinoma. Br J Cancer 68, 540-8 (1993) and Bhargava, V. et al. The association of p53 immunopositivity with tumor proliferation and other prognostic indicators in breast cancer. Mod Pathol 7, 361-8 (1994), all of which are incorporated herein by reference.
However, it was also observed that among the p53 wt tumors within the mutant-like cluster, there, too, was a significant over-representation of ER-(pcs=2.0×10−6) and grade III tumors (pfe=7.1×10−11). Thus, by univariate statistical analysis, a large number of genes highly associated with p53 status have been identified that are capable of segregating tumors in a manner correlated with p53 status, but also histologic grade and ER status.
The finding that a fraction of p53 wt tumors were found to cluster together with the majority of p53 mutants suggests the possibility that these tumors may in fact be p53 deficient through mechanisms other than p53 mutation. Conversely, the discovery of p53 mutants with molecular configurations reminiscent of most wt tumors suggests that these tumors might in fact express functionally intact p53. However, the tumor group assignments in this case were based on genes selected by a univariate ranking procedure that did not account for the association of p53 status with ER and grade status. This raised the possibility that, to some extent, the selected genes included those that are mostly grade and/or ER associated, which may have biased the clustering of the tumors towards these properties rather than p53 status, per se.
Therefore, a robust gene expression-based classifier for predicting p53 status was developed by designing a predictive model including a multivariate linear regression method known as linear model-fit (LMF) for ranking p53 status-correlated genes independent of histologic grade and ER status.
For gene selection, a linear model was fitted to the gene expression data with expression level as the response, and p53 status, ER status and grade status as the predictor variables. As an initial filter for removing genes not well correlated with the predictor variables, all genes with a p-value fit greater than 0.001 were excluded. Using ER and grade as additional predictors allowed for filtering out genes whose expression patterns could be mostly explained by either ER or grade status. When applied, the LMF ranking procedure markedly reduced the rank of many known cell cycle-regulated genes compared to the univariate Wilcoxon rank-sum (WR) method, indicating that these genes are best explained by high grade rather than p53 status (
For class prediction purposes, the genes were ranked in decreasing order of the absolute value of the p53 status coefficient. For building the classifier, a variant of the maximum likelihood method, DLDA (diagonal linear discriminant analysis) was employed. This had previously been applied to class determination problems using microarray data, described for example, in Dudoit, S., Frilyand, J. & Speed, T. P. Comparison of discrimination methods for the classification of tumors using gene expression data. Journal of the American Statistical Association 97, 77-87 (2002), incorporated herein by reference. The set of predictor genes with greatest classification accuracy was chosen by leave-one-out cross validation.
The accuracy of the classifier as a function of the number of genes it comprised is plotted in
C. elegans)
Drosophila)
Drosophila)
Drosophila);
The 500-gene classifier: The genes are ranked according to their correlation with p53 status. The genes are identified by their GenBank Accession Nos., Affymetrix Probeset IDs, Unigene IDs, Unigene Names and Unigene Symbols.
For sequences and SEQ ID NOs for the genes described in Table 1, see
The performance of the p53 classifier in the context of independent datasets was then evaluated.
Two publicly available microarray datasets where p53 status was known, were therefore accessed: a breast cancer study by Sorlie et al (Sorlie, T. et al. Repeated observation of breast tumor subtypes in independent gene expression data sets. Proc Natl Acad Sci USA 100, 8418-23 (2003), incorporated herein by reference) and a liver cancer study by Chen et al (Chen, X. et al. Gene expression patterns in human liver cancers. Mol Biol Cell 13, 1929-39 (2002), incorporated herein by reference). Both studies were conducted on cDNA microarray platforms.
In the Sorlie dataset, 69 breast tumors were sequenced for p53 mutations. This subset of tumors was queried for the availability of expression data corresponding to the genes of the classifier. Twenty-eight genes in the classifier mapped to UniGene IDs (build #167). Though over half of these genes mapped to the Sorlie et. al. microarray, few were expressed in the majority of the tumors, and a number of tumors possessed measurements for less than half of the genes. Only 9 genes in the classifier were found to correspond to cDNA probes (representing 9 different genes) having expression measurements present in >50% of the tumors, where the tumors possessed measurements for >50% of the genes (resulting in a subset of 44 well-sampled tumors). Using this 9-gene subset of the classifier to hierarchically cluster the tumors (
A cDNA-microarray based liver cancer dataset where p53 status was ascertained by immunohistochemistry, IHC (Chen, X. et al. Gene expression patterns in human liver cancers. Mol Biol Cell 13, 1929-39 (2002), incorporated herein by reference) was next analyzed. In this study, p53 protein levels were ascertained by IHC. Here, 8 classifier genes could be mapped to all 59 tumors assayed for p53 status (with each gene having data present in 90% or more of all tumors, and where each tumor contained data for >50% of the genes). With similar statistical significance as that seen in the breast cancer dataset (i.e, pfe=3.5×10−4), this 8-gene subset of the classifier was able to cluster the HCC samples into two predominant clusters correlated with p53 status: 87% of the mutants in one cluster, and 61% of the wildtypes in the other (
Optimized 32-gene p53 Classifier: The genes are identified by their GenBank Accession Nos., Affymetrix Probeset IDs, Unigene IDs, Unigene Names and Unigene Symbols.
It is widely accepted that in breast cancer and other tumor types p53 status is prognostic of clinical outcomes such as tumor recurrence, patient survival, and therapeutic response. The hypothesis that a classifier based on p53 activity would out-perform p53 mutation status alone as a prognostic indicator of clinical outcomes was tested.
The classifier and sequence-level p53 mutation status were compared with respect to their abilities to predict disease-specific survival (DSS) in all 257 patients of the Uppsala cohort regardless of treatment type or clinical stage.
The significance of the hazard ratio generated using the p53 classifier to segregate patients was an order of magnitude greater than that obtained using p53 mutation status alone (pw=0.00057 versus pw=0.012, respectively) (
Next, the prognostic significance of the classifier on the Sorlie et al cDNA microarray dataset was examined (Sorlie, T. et al. Repeated observation of breast tumor subtypes in independent gene expression data sets. Proc Natl Acad Sci U S A 100, 8418-23 (2003), incorporated herein by reference).
Here, the 9-gene partial classifier that could distinguish mt and wt tumors both with 77% accuracy, was used to hierarchically cluster 76 well-sampled tumor specimens with associated patient survival information (
For hierarchical cluster analysis, log expression values were mean centered and normalized, and genes and tumors were clustered using the Pearson correlation metric and average linkage (Cluster and TreeView software courtesy Dr. Michael Eisen; software available on Lawrence Berkeley National Laboratory, UC Berkeley's website). For survival analysis, patients were stratified according to the p53 classifier output or, as in one case, according to p53 mutation status. The Kaplan Meier estimate was used to compute survival curves for the different patient groups and the Wald Test was used to assess the statistical significance of the resultant hazard ratio. The
For association tests (i.e., to ascertain the significance of the number of observed events in two or more groups), the Chi-square test was employed. When the number of events was sufficiently small (<5) in any category, Fisher's Exact test was applied instead of Chi-square test.
For the statistical analysis of expression levels for p53 downstream target genes and upstream effectors, two-tailed two-group T tests were employed to determine differentially expressed genes between the p53 wt and mt tumors (
It would be evident to one of skill in the art that the method embodiments of the present invention are not limited to the statistical methods disclosed herein. Embodiments of the present invention encompass equivalent analytical methods. The p-value abbreviations used herein include:
pwr=Wilcoxon rank-sum test
pt=T test
pcs=Chi-square test
pfe=Fisher's Exact test
pw=Wald test
pmc=Monte Carlo estimate
Promoter analysis for p53 binding sites was performed on each of the classifier genes with a known transcription start site (TSS). BEARR (Vega, V. B., Bangarusamy, D. K., Miller, L. D., Liu, E. T. & Lin, C. Y. BEARR: Batch Extraction and Analysis of cis-Regulatory Regions. Nucleic Acids Res 32, W257-60 (2004), incorporated herein by reference) was used to extract promoter sequences (3000 bp upstream to 500 bp downstream of the TSS) and predict putative binding sites using the P53 position weight matrix obtained from TRANSFAC (Kel, A. E. et al. MATCH: A tool for searching transcription factor binding sites in DNA sequences. Nucleic Acids Res 31, 3576-9 (2003), incorporated herein by reference) version 6.0 (Matrix accession: M00272) as well as simple pattern search based on the canonical p53 binding site consensus 5′-RRRCWWGYYYN(0-13)RRRCWWGYYY-3′ (el-Deiry, W. S., Kern, S. E., Pietenpol, J. A., Kinzler, K. W. & Vogelstein, B. Definition of a consensus binding site for p53. Nat Genet 1, 45-9 (1992), incorporated herein by reference.
To further test the robustness of the classifier in predicting patient outcome, its performance in other relevant therapeutic treatment groups was analyzed. Recently, it has been observed that p53 mt breast tumors show greater resistance to endocrine therapy than p53 wt tumors, and this has been explained, in part, by the uncoupling of p53-dependent apoptosis in the resistant tumors (Berns, E. M. et al. Complete sequencing of TP53 predicts poor response to systemic therapy of advanced breast cancer. Cancer Res 60, 2155-62 (2000), incorporated herein by reference). To test the ability of the classifier to predict outcome in a hormone therapy-specific patient cohort, a subpopulation of the Uppsala cohort consisting of 68 ER+ patients who received only adjuvant tamoxifen treatment following surgery, was examined.
Next, the prognostic performance of the classifier on a set of 97 breast tumors published by van't Veer et al (van't Veer, L. J. et al. Gene expression profiling predicts clinical outcome of breast cancer. Nature 415, 530-6 (2002), incorporated herein by reference) was examined.
Here, all of the samples were controlled for clinical uniformity, i.e., <5 cm in size (T1/T2), with no advanced disease (pN0), from patients less than 55 years of age at diagnosis, treated by surgery and subsequent radiotherapy only (with the exception of 5 patients who received adjuvant systemic therapy). From the 32-gene classifier, 24 probes corresponding to 21 genes could be mapped to all 97 tumors with survival information. Upon clustering the tumors, approximately 4 clusters with similar average distance correlations were observed that significantly distinguished patients who would develop a distant metastasis within 5 years (pfe=2.2×10−4) (
To gain some mechanistic insights, the functional annotations of the classifier genes were analysed for clues to explain the correlation between their expression levels and p53 status and patient outcome. Surprisingly, it was found that none of the classifier genes are known transcriptional targets of p53, nor have they been previously implicated in the p53 pathway. Promoter analysis of the 21 genes with defined promoter regions revealed no evidence of the canonical p53 binding site, or recently described novel p53 binding sites, within any of the promoters.
Twelve of the genes are of unknown function. However, of the characterized genes, a number are associated with cell growth and proliferation (MYBL2, TFF1, BRRN1, CHAD, SCGB3A1, DACH, CDCA8), transcription (LAF4, NY-BR-1, DACH, MYBL2), ion transport (CACNG4, CYBRD1, LRP2), and breast cancer biology (SCGB3A1, TFF1, STC2, NY-BR-1, AGR2). Speculatively, some of these genes may contribute mechanistically to the poor prognosis of the p53 mutant-like tumors. For example, MYBL2, which was observed to be upregulated in the p53 mutant-like tumors, is a growth-promoting transcription factor closely related to the c-MYB oncogene. It maps to a chromosomal region frequently amplified in breast cancer (20q13) and has previously been reported to be overexpressed in breast cancer cell lines and sporadic ovarian carcinomas (Forozan, F. et al. Comparative genomic hybridization analysis of 38 breast cancer cell lines: a basis for interpreting complementary DNA microarray data. Cancer Res 60, 4519-25 (2000) and Tanner, M. M. et al. Frequent amplification of chromosomal region 20q12-q13 in ovarian cancer. Clin Cancer Res 6, 1833-9 (2000), both of which are incorporated herein by reference. SCGB3A1 (HIN1), which was observed to be downregulated in the p53 mutant-like tumors, is a putative tumor suppressor gene that can inhibit breast cancer cell growth when overexpressed and has been found to be transcriptionally silenced by hypermethylation of its promoter in early stages of breast tumorigenesis (Krop, I. E. et al. HIN-1, a putative cytokine highly expressed in normal but not cancerous mammary epithelial cells. Proc Natl Acad Sci USA 98, 9796-801 (2001), incorporated herein by reference).
It was observed that a number of cancers with wild type p53 sequence status were classified as p53 mutant by expression profiling using the 32-gene classifier. If the “misclassified” p53 wt tumors were in fact p53 deficient, they would possess certain molecular characteristics reflective of perturbations of the p53 pathway, and these characteristics would be found in the majority of p53 mutant tumors. First, the possibility that p53 deficiency could result from reduced transcript levels either by transcriptional repression of the p53 gene (TP53) or by the shortening of its mRNA half-life, was considered. The t test was used to compare the relative expression levels of TP53 (using the TP53 probe-sets present on the microarray) among the different tumor classes (
Table 3 shows a comparative analysis of p53 mutations. (I) Severe mutations were defined as insertions, deletions, or stop codons. Of the remaining missense point mutations (mpms; 11 in the wt-like group, 27 in the mt-like group) we determined the frequency of occurrence of (II) the most common missense point mutations in p53 as defined by the IARC TP53 Mutation Database (available online on the website of the International Agency for Research on Cancer, IARC), and (III) mutants previously shown, in vitro, to possess dominant negative activity were determined. P-values were calculated using Fisher's Exact test.
This strategy was applied to known transcriptional targets of p53, which were hypothesized to show altered transcription in p53-deficient tumors to some extent. Indeed, a number of p53 target genes demonstrated altered patterns of expression (
CHEK1 and CHEK2, both positive upstream effectors of p53 that phosphorylate p53 and thereby promote its stabilization, are known to be transcriptionally repressed by p53. A significant increase in the mRNA levels of these genes in both the p53 wt and mt tumors of the mutant-like class was observed. It was also observed that the 12 mt tumors misclassified as wildtype-like displayed significantly lower expression of these genes compared to the other 47 p53 mutants. Notably, no differential expression of the p53-regulated genes CDKN1A (p21), GADD45, PPM1D (WIP1), TP5313 (PIG3), TNFRSF6, BBC3 (PUMA), APAF1 or BCL2 was observed in these breast tumor specimens.
Taken together, these data suggest that the classifier can distinguish tumors based on some aspects of p53 transcriptional activity that are inhibited in both the p53 mutant and wildtype tumors of the mutant-like class, yet operative in the p53 wildtype tumors (and to some extent the 12 p53 mutant tumors) of the wildtype-like class.
Perhaps paradoxically, it was observed that the p53-inducible genes PERP, BAX and SFN (14-3-3 sigma) were all expressed at significantly higher levels in the 28 misclassified wt tumors, rather than at lower levels like their inducible gene counterparts described above. However, the significant overexpression of these genes in the p53 mt tumors classified as mutant-like was also observed, suggesting that in breast cancer, these genes may be induced by alternate regulatory mechanisms in the context of mutant or deficient p53.
Intriguingly, another positive upstream effector of p53, ATR, which is thought to enhance p53 activity in a manner similar to that of CHEK1 and CHEK2, was also found expressed at significantly higher levels in the p53 mutants and p53 wt tumors of the mutant-like class, even though this gene is not known to be modulated in a p53-dependent manner. Of note, no significant differences in the expression levels of the upstream effectors, ATM or PRKDC (DNA-PK) were observed.
The expression levels of other upstream modulators of p53 activity were then examined in order to ascertain possible alternate mechanisms by which p53 expression and activity might be reduced in the mutant-like p53 wt tumors. First, it was observed that several known positive regulators of p53 transactivation were significantly reduced in both the wildtypes and mutants of the mutant-like class including HOXA5, USF1, EGR1 and TP53BP1. HOXA5, USF1, and EGR1 are all transcription factors known to bind the p53 promoter and enhance its expression. Interestingly, deficiencies in all three have previously been implicated in breast carcinogenesis. Recently the coordinate loss of both p53 and HOXA5 mRNA and protein expression was observed in a panel of human breast cancer cell lines, and the HOXA5 promoter was found to be methylated in 16 of 20 p53-negative human breast tumors. USF1, which is structurally related to the c-Myc oncoprotein, has been found to have reduced transcriptional activity in breast cancer cell lines, and has recently been shown to activate the expression of estrogen receptor alpha. EGR1, a DNA damage-responsive gene with antiproliferative and apoptotic functions, can inhibit tumorigenicity when exogenously expressed in human breast cancer cells, and has been observed to have reduced expression in human and mouse breast cancer cell lines and tumors. TP53BP1 is not thought to be a transcription factor, but rather a BRCT domain-containing substrate of ATM that is phosphorylated in response to DNA damage. This gene product is known to bind the central DNA-binding domain of p53 and thus enhance the transcriptional activation of p53 target genes. A significantly reduced expression of all four genes in the 28 p53 wt tumors classified as mutant-like was found, and in the cases of USF1 and TP53BP1, significantly higher expression in the p53 mutants classified as wildtype-like. Interestingly, it was also observed that their expression levels are also significantly lower in the 47 p53 mt tumors classified as mutant-like, suggesting a possible positive feedback loop whereby wildtype p53 can enhance expression of these genes and impaired p53 cannot. Together, these observations suggest the possibility that either acting separately or in combination, these genes may be important for intact p53 activity in the breast, and when transcriptionally silenced, contribute to p53 deficiency.
Finally, the expression of several known negative regulators of p53 activity were examined. Notably, MDM2, which negatively regulates p53 through phosphorylation-mediated degradation of the p53 protein, and whose overexpression at the protein level has been implicated in a variety of cancers, was not found to be differentially expressed at the transcript level in the experiments described herein. However, both PLK1 and GTSE1 were. The M-phase regulator PLK1 has recently been shown to bind to the DNA-binding domain of p53 and thus inhibit its transcriptional activity in vitro. GTSE1 (B99) binds the C-terminal regulatory domain of p53 causing the inhibition of p53 transactivation function as well as a reduction of intracellular levels of p53 protein. Intriguingly, the transcript levels of both genes were among the most highly significantly overexpressed in both p53 wt and mt tumors of the mt-like class, suggesting a possible role for these gene products in suppression of p53 function in breast carcinogenesis.
The spectrum of p53 mutations for correlations that might explain the misclassification of the 12 p53-mutant tumors as wildtype-like was next analyzed. First, it was observed that only one mutation was common to the wildtype-like and the mutant-like tumors: a Tyr>Cys at amino acid 220 in the DNA-binding domain. Of the 47 p53 mt tumors correctly classified as mutants, it was observed that 42% (20/47) possessed “severe” mutations defined as insertions (n=2), deletions (n=11) and stop codons (n=7) (Table 3-I) resulting in frameshifts and subsequent trunctation, whereas in the 12 mutants classified as wildtype-like, only 1 (8%) contained a severe mutation: a 3-bp insertion in the DNA-binding domain resulting in the inframe addition of a glycine residue (pfe=0.025). Using the IARC TP53 Mutation Database (available online on the website of the International Agency for Research on Cancer, IARC), which, as of June 2003, has indexed 18,585 somatic and 225 germline mutations of p53, the frequencies of occurrence of the most common p53 mutations in human cancer (representing ˜20% of all p53 mutations; Table 1-II) in the 12 wt-like mutants and the 47 mt-like mutants were compared. None of the common mutations were found to overlap with the subset of 11 missense point mutations (mpms) in the wt-like group, compared to 9 of 27 in the mt-like group (pfe=0.029). The mpms in each tumor group was then cross-compared with the IARC TP53 Mutation Database's comprehensive listing of 418 mutants previously analyzed for dominant negative function in at least one of 44 previously published studies. As Table 2-III shows, it was found that only one of the 11 mpms among the 12 wt-like mutants had been demonstrated previously to have dominant negative activity, compared to 12 of 27 within the mt-like group (pfe=0.039). Together, these data suggest that at the sequence level, the 12 p53 mutants classified as wildtype-like may in fact comprise of mostly “benign” p53 mutant forms compared to those 47 classified as mutant-like, in agreement with their molecular consistencies with the majority of p53 wt tumors in our expression analyses.
Comparative analysis of p53 mutations. (I) Severe mutations were defined as insertions, deletions, or stop codons. Of the remaining missense point mutations (mpms; 11 in the wt-like group, 27 in the mt-like group) we determined the frequency of occurrence of (II) the most common missense point mutations in p53 as defined by the IARC TP53 Mutation Database (http://www.iarc.fr/p53/index.html), and (III) mutants previously shown, in vitro, to possess dominant negative activity. P-values were calculated using Fisher's Exact test.
The practice of the present invention may employ conventional biology methods known to the skilled artisan, software and systems. The foregoing examples have described methods for predicting disease outcome in a patient. In another aspect, there is also provided a computer system for predicting disease outcome in a patient. The computer system may comprise a computer having a processor and a memory, the memory having executable code stored thereon for execution by the processor for performing the steps of obtaining gene expression profiles from a plurality of genes from tumor samples, wherein said tumor samples may be mutant or wildtype for the p53 gene; comparing said gene expression profiles to determine which genes are differentially expressed in the mutant or wildtype tumors; deriving from said differentially expressed genes a set of genes to predict p53 mutational status; and using the set of genes to predict disease outcome in the patient.
A suitable computer system may be a general purpose computer such as a PC or a Macintosh, for example. Computer software products of the invention typically include a computer readable medium having computer-executable instructions for performing the logic steps of the method of the invention. Suitable computer readable media include floppy disk, CD-ROM/DVD/DVD-ROM, hard-disk drive, flash memory, ROM/RAM, magnetic tapes etc. The computer executable instructions may be written in a suitable computer language or a combination of several languages. Basic computational biology methods are described in, e.g. Setubal and Meidanis et al., Introduction to Computational Biology Methods (PWS Publishing Company, Boston, 1997); Salzberg, Searles, Kasif, (Ed.), Computational Methods in Molecular Biology, (Elsevier, Amsterdam, 1998); Rashidi and Buehler, Bioinformatics Basics: Application in Biological Science and Medicine (CRC Press, London, 2000) and Ouelette and Bzevanis Bioinformatics: A Practical Guide for Analysis of Gene and Proteins (Wiley & Sons, Inc., 2nd Ed., 2001).
Additionally, the present invention may have preferred embodiments that include methods for providing genetic information over networks such as the Internet.
Additionally, some embodiments of the present invention may provide a plurality of pharmaceutical targets for designing chemotherapeutic drugs for a variety of cancers. For example, the 32 genes most correlated with p53 mutational status could serve as potential molecular targets for chemotherapy. Chemotherapy drugs (cytotoxics) and antihormonal treatments are commonly used to treat cancers. In several patients however, treatment regimens involving cytotoxics and antihormonals have been known to cause mild to severe side effects. In breast cancer for example, these side effects include vomiting, nausea, alopecia and fatigue. The future of effective treatment for cancer thus resides with drugs that are more specific for their targets. According to some studies, about 68% of breast cancer drugs in the clinical developmental pipeline are of the targeted class. Therefore, molecular signatures such as those embodied in certain aspects of the present invention will provide important leads or will prove to be targets in their own right for targeted chemotherapeutic drugs.
In conclusion, the disclosed embodiments of the present invention define a gene expression signature a gene expression signature that can predict p53 status and survival in human breast tumours (the p53 signature or classifier). In independent datasets of both breast and liver cancers, and regardless of other clinical features, subsets of the p53 signature can predict p53 status with significant accuracy. As a predictor of disease-specific survival (DSS), the signature significantly outperformed p53 mutation status alone in a large patient cohort with heterogeneous treatment. The p53 signature could significantly distinguish patients having more or less benefit from systemic adjuvant therapies and loco-regional radiotherapy. Though the p53 pathway may be compromised at some level in most human cancers, analysis of transcripts involved in the p53 pathway suggests that the p53 expression signature defines an operational configuration of this pathway in breast tumors (more so than p53 mutation status alone) that impacts patient survival, and therapeutic response. In cancer, it is clear that not all p53 mutations have equal effects: some simply confer loss of function, while others have a dominant negative effect (such as trans-dominant suppression of wildtype p53 or oncogenic gain of function), while still others show only a partial loss of function where, for example, only a small subset of p53 downstream transcriptional target genes are dysregulated. For these reasons, no single molecular assessment of p53 status appears to provide an absolute indication of the complete p53 function. The embodiments disclosed herein suggest that by looking at the downstream indicators of p53 function, the functional status of p53 may be ascertained more precisely than using sequencing or biochemical means.
It is to be understood that the above description in intended to be illustrative and not restrictive. Many variations of the invention will be apparent to those of skill in the art upon reviewing the above description. The scope of the invention should be determined with reference to the appended claims, along with the full scope of equivalents to which such claims are entitled. All cited references, including patent and non-patent literature, are incorporated herewith by reference in their entireties for all purposes.
| Number | Name | Date | Kind |
|---|---|---|---|
| 6208983 | Parra et al. | Mar 2001 | B1 |
| 6306087 | Barnhill et al. | Oct 2001 | B1 |
| 6468476 | Friend et al. | Oct 2002 | B1 |
| 6714925 | Barnhill et al. | Mar 2004 | B1 |
| 6757412 | Parsons et al. | Jun 2004 | B1 |
| 20040058340 | Dai et al. | Mar 2004 | A1 |
| Number | Date | Country |
|---|---|---|
| WO 2004065583 | Aug 2004 | WO |
| Number | Date | Country | |
|---|---|---|---|
| 20060074565 A1 | Apr 2006 | US |
