Patent Grant
10422005
This invention relates to methods and compositions for determining the prognosis of cancer, particularly colorectal cancer, in a patient. Specifically, this invention relates to the use of genetic markers for determining the prognosis of cancer, such as colorectal cancer, based on prognostic signatures.
Colorectal cancer (CRC) is one of the most common cancers in the developed world, and its incidence is continuing to increase. Although the progression of colorectal cancer from benign polyp to adenoma to carcinoma is well studied (1), the molecular events influencing the transition and establishment of metastasis are less well understood. The prognosis and treatment of CRC currently depends on the clinico-pathological stage of disease at the time of diagnosis, and primary surgical treatment. Unfortunately disease stage alone does not allow accurate prediction of outcome for individual patients. If patient outcomes could be predicted more accurately treatments could be tailored to avoid under-treating patients destined to relapse, or over-treating patients who would be helped by surgery alone.
Many attempts have been made to identify markers that predict clinical outcome in CRC. Until recently most studies focused on single proteins or gene mutations with limited success in terms of prognostic information (2). Microarray technology enables the identification of sets of genes, called classifiers or signatures that correlate with cancer outcome. This approach has been applied to a variety of cancers, including CRC (3-5), but methodological problems and a lack of independent validation has cast doubt over the findings (6, 7). Furthermore, doubts about the ability of classifiers/signatures to predict outcome have arisen due to poor concordance of identified by different researchers using different array platforms and methodologies (8).
There is a need for further tools to predict the prognosis of colorectal cancer. This invention provides further methods, compositions, kits, and devices based on prognostic cancer markers, specifically colorectal cancer prognostic markers, to aid in the prognosis and treatment of cancer.
In certain embodiments there is provided a set of markers genes identified to be differentially expressed in recurrent and non-recurrent colorectal tumours. This set of genes can be used to generate prognostics signatures, comprising two or more markers, capable of predicting the progression of colorectal tumour in a patient.
The individual markers can differentially expressed depending on whether the tumour is recurrent or not. The accuracy of prediction can be enhanced by combining the markers together into a prognostic signature for, providing for much more effective individual tests than single-gene assays. Also provided for is the application of techniques, such as statistics, machine learning, artificial intelligence, and data mining to the prognostics signatures to generate prediction models. In another embodiment, expression levels of the markers of a particular prognostic signature in the tumour of a patient can then be applied to the prediction model to determine the prognosis.
In certain embodiments, the expression level of the markers can be established using microarray methods, quantitative polymerase chain reaction (qPCR), or immunoassays.
This invention is described with reference to specific embodiments thereof and with reference to the figures, in which:
Definitions
Before describing embodiments of the invention in detail, it will be useful to provide some definitions of terms used herein.
The term “marker” refers to a molecule that is associated quantitatively or qualitatively with the presence of a biological phenomenon. Examples of “markers” include a polynucleotide, such as a gene or gene fragment, RNA or RNA fragment; or a gene product, including a polypeptide such as a peptide, oligopeptide, protein, or protein fragment; or any related metabolites, by products, or any other identifying molecules, such as antibodies or antibody fragments, whether related directly or indirectly to a mechanism underlying the phenomenon. The markers of the invention include the nucleotide sequences (e.g., GenBank sequences) as disclosed herein, in particular, the full-length sequences, any coding sequences, any fragments, or any complements thereof, and any measurable marker thereof as defined above.
The terms “CCPM” or “colorectal cancer prognostic marker” or “CCPM family member” refer to a marker with altered expression that is associated with a particular prognosis, e.g., a higher or lower likelihood of recurrence of cancer, as described herein, but can exclude molecules that are known in the prior art to be associated with prognosis of colorectal cancer. It is to be understood that the term CCPM does not require that the marker be specific only for colorectal tumours. Rather, expression of CCPM can be altered in other types of tumours, including malignant tumours.
The terms “prognostic signature,” “signature,” and the like refer to a set of two or more markers, for example CCPMs, that when analysed together as a set allow for the determination of or prediction of an event, for example the prognostic outcome of colorectal cancer. The use of a signature comprising two or more markers reduces the effect of individual variation and allows for a more robust prediction. Non-limiting examples of CCPMs are set forth in Tables 1, 2, 5, and 9, while non-limiting examples of prognostic signatures are set forth in Tables 3, 4, 8A, 8B, and 9, herein. In the context of the present invention, reference to “at least one,” “at least two,” “at least five,” etc., of the markers listed in any particular set (e.g., any signature) means any one or any and all combinations of the markers listed.
The term “prediction method” is defined to cover the broader genus of methods from the fields of statistics, machine learning, artificial intelligence, and data mining, which can be used to specify a prediction model. These are discussed further in the Detailed Description section.
The term “prediction model” refers to the specific mathematical model obtained by applying a prediction method to a collection of data. In the examples detailed herein, such data sets consist of measurements of gene activity in tissue samples taken from recurrent and non-recurrent colorectal cancer patients, for which the class (recurrent or non-recurrent) of each sample is known. Such models can be used to (1) classify a sample of unknown recurrence status as being one of recurrent or non-recurrent, or (2) make a probabilistic prediction (i.e., produce either a proportion or percentage to be interpreted as a probability) which represents the likelihood that the unknown sample is recurrent, based on the measurement of mRNA expression levels or expression products, of a specified collection of genes, in the unknown sample. The exact details of how these gene-specific measurements are combined to produce classifications and probabilistic predictions are dependent on the specific mechanisms of the prediction method used to construct the model.
“Sensitivity”, “specificity” (or “selectivity”), and “classification rate”, when applied to the describing the effectiveness of prediction models mean the following:
“Sensitivity” means the proportion of truly positive samples that are also predicted (by the model) to be positive. In a test for CRC recurrence, that would be the proportion of recurrent tumours predicted by the model to be recurrent. “Specificity” or “selectivity” means the proportion of truly negative samples that are also predicted (by the model) to be negative. In a test for CRC recurrence, this equates to the proportion of non-recurrent samples that are predicted to by non-recurrent by the model. “Classification Rate” is the proportion of all samples that are correctly classified by the prediction model (be that as positive or negative).
As used herein “antibodies” and like terms refer to immunoglobulin molecules and immunologically active portions of immunoglobulin (Ig) molecules, i.e., molecules that contain an antigen binding site that specifically binds (immunoreacts with) an antigen. These include, but are not limited to, polyclonal, monoclonal, chimeric, single chain, Fc, Fab, Fab′, and Fab2 fragments, and a Fab expression library. Antibody molecules relate to any of the classes IgG, IgM, IgA, IgE, and IgD, which differ from one another by the nature of heavy chain present in the molecule. These include subclasses as well, such as IgG1, IgG2, and others. The light chain may be a kappa chain or a lambda chain. Reference herein to antibodies includes a reference to all classes, subclasses, and types. Also included are chimeric antibodies, for example, monoclonal antibodies or fragments thereof that are specific to more than one source, e.g., a mouse or human sequence. Further included are camelid antibodies, shark antibodies or nanobodies.
The terms “cancer” and “cancerous” refer to or describe the physiological condition in mammals that is typically characterized by abnormal or unregulated cell growth. Cancer and cancer pathology can be associated, for example, with metastasis, interference with the normal functioning of neighbouring cells, release of cytokines or other secretory products at abnormal levels, suppression or aggravation of inflammatory or immunological response, neoplasia, premalignancy, malignancy, invasion of surrounding or distant tissues or organs, such as lymph nodes, etc. Specifically included are colorectal cancers, such as, bowel (e.g., large bowel), anal, and rectal cancers.
The term “colorectal cancer” includes cancer of the colon, rectum, and/or anus, and especially, adenocarcinomas, and may also include carcinomas (e.g., squamous cloacogenic carcinomas), melanomas, lymphomas, and sarcomas. Epidermoid (nonkeratinizing squamous cell or basaloid) carcinomas are also included. The cancer may be associated with particular types of polyps or other lesions, for example, tubular adenomas, tubulovillous adenomas (e.g., villoglandular polyps), villous (e.g., papillary) adenomas (with or without adenocarcinoma), hyperplastic polyps, hamartomas, juvenile polyps, polypoid carcinomas, pseudopolyps, lipomas, or leiomyomas. The cancer may be associated with familial polyposis and related conditions such as Gardner's syndrome or Peutz-Jeghers syndrome. The cancer may be associated, for example, with chronic fistulas, irradiated anal skin, leukoplakia, lymphogranuloma venereum, Bowen's disease (intraepithelial carcinoma), condyloma acuminatum, or human papillomavirus. In other aspects, the cancer may be associated with basal cell carcinoma, extramammary Paget's disease, cloacogenic carcinoma, or malignant melanoma.
The terms “differentially expressed,” “differential expression,” and like phrases, refer to a gene marker whose expression is activated to a higher or lower level in a subject (e.g., test sample) having a condition, specifically cancer, such as colorectal cancer, relative to its expression in a control subject (e.g., reference sample). The terms also include markers whose expression is activated to a higher or lower level at different stages of the same condition; in recurrent or non-recurrent disease; or in cells with higher or lower levels of proliferation. A differentially expressed marker may be either activated or inhibited at the polynucleotide level or polypeptide level, or may be subject to alternative splicing to result in a different polypeptide product. Such differences may be evidenced by a change in mRNA levels, surface expression, secretion or other partitioning of a polypeptide, for example.
Differential expression may include a comparison of expression between two or more markers (e.g., genes or their gene products); or a comparison of the ratios of the expression between two or more markers (e.g., genes or their gene products); or a comparison of two differently processed products (e.g., transcripts or polypeptides) of the same marker, which differ between normal subjects and diseased subjects; or between various stages of the same disease; or between recurring and non-recurring disease; or between cells with higher and lower levels of proliferation; or between normal tissue and diseased tissue, specifically cancer, or colorectal cancer. Differential expression includes both quantitative, as well as qualitative, differences in the temporal or cellular expression pattern in a gene or its expression products among, for example, normal and diseased cells, or among cells which have undergone different disease events or disease stages, or cells with different levels of proliferation.
The term “expression” includes production of polynucleotides and polypeptides, in particular, the production of RNA (e.g., mRNA) from a gene or portion of a gene, and includes the production of a polypeptide encoded by an RNA or gene or portion of a gene, and the appearance of a detectable material associated with expression. For example, the formation of a complex, for example, from a polypeptide-polypeptide interaction, polypeptide-nucleotide interaction, or the like, is included within the scope of the term “expression”. Another example is the binding of a binding ligand, such as a hybridization probe or antibody, to a gene or other polynucleotide or oligonucleotide, a polypeptide or a protein fragment, and the visualization of the binding ligand. Thus, the intensity of a spot on a microarray, on a hybridization blot such as a Northern blot, or on an immunoblot such as a Western blot, or on a bead array, or by PCR analysis, is included within the term “expression” of the underlying biological molecule.
The terms “expression threshold,” and “defined expression threshold” are used interchangeably and refer to the level of a marker in question outside which the polynucleotide or polypeptide serves as a predictive marker for patient survival without cancer recurrence. The threshold will be dependent on the predictive model established are derived experimentally from clinical studies such as those described in the Examples below. Depending on the prediction model used, the expression threshold may be set to achieve maximum sensitivity, or for maximum specificity, or for minimum error (maximum classification rate). For example a higher threshold may be set to achieve minimum errors, but this may result in a lower sensitivity. Therefore, for any given predictive model, clinical studies will be used to set an expression threshold that generally achieves the highest sensitivity while having a minimal error rate. The determination of the expression threshold for any situation is well within the knowledge of those skilled in the art.
The term “long-term survival” is used herein to refer to survival for at least 5 years, more preferably for at least 8 years, most preferably for at least 10 years following surgery or other treatment.
The term “microarray” refers to an ordered or unordered arrangement of capture agents, preferably polynucleotides (e.g., probes) or polypeptides on a substrate. See, e.g., Microarray Analysis, M. Schena, John Wiley & Sons, 2002; Microarray Biochip Technology, M. Schena, ed., Eaton Publishing, 2000; Guide to Analysis of DNA Microarray Data, S. Knudsen, John Wiley & Sons, 2004; and Protein Microarray Technology, D. Kambhampati, ed., John Wiley & Sons, 2004.
The term “oligonucleotide” refers to a polynucleotide, typically a probe or primer, including, without limitation, single-stranded deoxyribonucleotides, single- or double-stranded ribonucleotides, RNA: DNA hybrids, and double-stranded DNAs. Oligonucleotides, such as single-stranded DNA probe oligonucleotides, are often synthesized by chemical methods, for example using automated oligonucleotide synthesizers that are commercially available, or by a variety of other methods, including in vitro expression systems, recombinant techniques, and expression in cells and organisms.
The term “polynucleotide,” when used in the singular or plural, generally refers to any polyribonucleotide or polydeoxribonucleotide, which may be unmodified RNA or DNA or modified RNA or DNA. This includes, without limitation, single- and double-stranded DNA, DNA including single- and double-stranded regions, single- and double-stranded RNA, and RNA including single- and double-stranded regions, hybrid molecules comprising DNA and RNA that may be single-stranded or, more typically, double-stranded or include single- and double-stranded regions. Also included are triple-stranded regions comprising RNA or DNA or both RNA and DNA. Specifically included are mRNAs, cDNAs, and genomic DNAs, and any fragments thereof. The term includes DNAs and RNAs that contain one or more modified bases, such as tritiated bases, or unusual bases, such as inosine. The polynucleotides of the invention can encompass coding or non-coding sequences, or sense or antisense sequences. It will be understood that each reference to a “polynucleotide” or like term, herein, will include the full-length sequences as well as any fragments, derivatives, or variants thereof.
“Polypeptide,” as used herein, refers to an oligopeptide, peptide, or protein sequence, or fragment thereof, and to naturally occurring, recombinant, synthetic, or semi-synthetic molecules. Where “polypeptide” is recited herein to refer to an amino acid sequence of a naturally occurring protein molecule, “polypeptide” and like terms, are not meant to limit the amino acid sequence to the complete, native amino acid sequence for the full-length molecule. It will be understood that each reference to a “polypeptide” or like term, herein, will include the full-length sequence, as well as any fragments, derivatives, or variants thereof.
The term “prognosis” refers to a prediction of medical outcome, for example, a poor or good outcome (e.g., likelihood of long-term survival); a negative prognosis, or poor outcome, includes a prediction of relapse, disease progression (e.g., tumour growth or metastasis, or drug resistance), or mortality; a positive prognosis, or good outcome, includes a prediction of disease remission, (e.g., disease-free status), amelioration (e.g., tumour regression), or stabilization.
The term “proliferation” refers to the processes leading to increased cell size or cell number, and can include one or more of: tumour or cell growth, angiogenesis, innervation, and metastasis.
The term “qPCR” or “QPCR” refers to quantative polymerase chain reaction as described, for example, in PCR Technique: Quantitative PCR, J. W. Larrick, ed., Eaton Publishing, 1997, and A-Z of Quantitative PCR, S. Bustin, ed., IUL Press, 2004.
The term “tumour” refers to all neoplastic cell growth and proliferation, whether malignant or benign, and all pre-cancerous and cancerous cells and tissues.
“Stringency” of hybridization reactions is readily determinable by one of ordinary skill in the art, and generally is an empirical calculation dependent upon probe length, washing temperature, and salt concentration. In general, longer probes require higher temperatures for proper annealing, while shorter probes need lower temperatures. Hybridization generally depends on the ability of denatured DNA to reanneal when complementary strands are present in an environment below their melting temperature. The higher the degree of desired homology between the probe and hybridisable sequence, the higher the relative temperature which can be used. As a result, it follows that higher relative temperatures would tend to make the reaction conditions more stringent, while lower temperatures less so. Additional details and explanation of stringency of hybridization reactions, are found e.g., in Ausubel et al., Current Protocols in Molecular Biology, Wiley Interscience Publishers, (1995).
“Stringent conditions” or “high stringency conditions”, as defined herein, typically: (1) employ low ionic strength and high temperature for washing, for example 0.015 M sodium chloride/0.0015 M sodium citrate/0.1% sodium dodecyl sulfate at 50° C.; (2) employ a denaturing agent during hybridization, such as formamide, for example, 50% (v/v) formamide with 0.1% bovine serum albumin/0.1% Ficoll/0.1% polyvinylpyrrolidone/50 mM sodium phosphate buffer at pH 6.5 with 750 mM sodium chloride, 75 mM sodium citrate at 42° C.; or (3) employ 50% formamide, 5×SSC (0.75 M NaCl, 0.075 M sodium citrate), 50 mM sodium phosphate (pH 6.8), 0.1% sodium pyrophosphate, 5×, Denhardt's solution, sonicated salmon sperm DNA (50 μg/ml), 0.1% SDS, and 10% dextran sulfate at 42° C., with washes at 42° C. in 0.2×SSC (sodium chloride/sodium citrate) and 50% formamide at 55° C., followed by a high-stringency wash comprising 0.1×SSC containing EDTA at 55° C.
“Moderately stringent conditions” may be identified as described by Sambrook et al., Molecular Cloning: A Laboratory Manual, New York: Cold Spring Harbor Press, 1989, and include the use of washing solution and hybridization conditions (e.g., temperature, ionic strength, and % SDS) less stringent that those described above. An example of moderately stringent conditions is overnight incubation at 37° C. in a solution comprising: 20% formamide, 5×SSC (150 mM NaCl, 15 mM trisodium citrate), 50 mM sodium phosphate (pH 7.6), 5×Denhardt's solution, 10% dextran sulfate, and 20 mg/ml denatured sheared salmon sperm DNA, followed by washing the filters in 1×SSC at about 37-50° C. The skilled artisan will recognize how to adjust the temperature, ionic strength, etc. as necessary to accommodate factors such as probe length and the like.
The practice of the present invention will employ, unless otherwise indicated, conventional techniques of molecular biology (including recombinant techniques), microbiology, cell biology, and biochemistry, which are within the skill of the art. Such techniques are explained fully in the literature, such as, Molecular Cloning: A Laboratory Manual, 2nd edition, Sambrook et al., 1989; Oligonucleotide Synthesis, MJ Gait, ed., 1984; Animal Cell Culture, R. I. Freshney, ed., 1987; Methods in Enzymology, Academic Press, Inc.; Handbook of Experimental Immunology, 4th edition, D. M. Weir & C C. Blackwell, eds., Blackwell Science Inc., 1987; Gene Transfer Vectors for Mammalian Cells, J. M. Miller & M. P. Calos, eds., 1987; Current Protocols in Molecular Biology, F. M. Ausubel et al., eds., 1987; and PCR: The Polymerase Chain Reaction, Mullis et al., eds., 1994.
Description Of Embodiments Of The Invention
In colorectal cancer, discordant results have been reported for prognostic markers. The present invention discloses the use of microarrays to reach a firmer conclusion, and to determine the prognostic role of specific prognostic signatures in colorectal cancer. The microarray-based studies shown herein indicate that particular prognostic signatures in colorectal cancer are associated with a prognosis. The invention can therefore be used to identify patients at high risk of recurrence of cancer, or patients with a high likelihood of recovery.
The present invention provides for markers for the determination of disease prognosis, for example, the likelihood of recurrence of tumours, including colorectal tumours. Using the methods of the invention, it has been found that numerous markers are associated with the prognosis of colorectal cancer, and can be used to predict disease outcome. Microarray analysis of samples taken from patients with various stages of colorectal tumours has led to the surprising discovery that specific patterns of marker expression are associated with prognosis of the cancer. The present invention therefore provides for a set of genes, outlined in Table 1 and Table 2, that are differentially expressed in recurrent and non-recurrent colorectal cancers. The genes outlined in Table 1 and Table 2 provide for a set of colorectal cancer prognostic makers (CCPMs).
A decrease in certain colorectal cancer prognostic markers (CCPMs), for example, markers associated with immune responses, is indicative of a particular prognosis. This can include increased likelihood of cancer recurrence after standard treatment, especially for colorectal cancer. Conversely, an increase in other CCPMs is indicative of a particular prognosis. This can include disease progression or the increased likelihood of cancer recurrence, especially for colorectal cancer. A decrease or increase in expression can be determined, for example, by comparison of a test sample, e.g., patient's tumour sample, to a reference sample, e.g., a sample associated with a known prognosis. In particular, one or more samples from patient(s) with non-recurrent cancer could be used as a reference sample.
For example, to obtain a prognosis, expression levels in a patient's sample (e.g., tumour sample) can be compared to samples from patients with a known outcome. If the patient's sample shows increased or decreased expression of one or more CCPMs that compares to samples with good outcome (no recurrence), then a positive prognosis, or recurrence is unlikely, is implicated. If the patient's sample shows expression of one or more CCPMs that is comparable to samples with poor outcome (recurrence), then a positive prognosis, or recurrence of the tumour is likely, is implicated.
As further examples, the expression levels of a prognostic signature comprising two or more CCPMS from a patient's sample (e.g., tumour sample) can be compared to samples of recurrent/non-recurrent cancer. If the patient's sample shows increased or decreased expression of CCPMs by comparison to samples of non-recurrent cancer, and/or comparable expression to samples of recurrent cancer, then a negative prognosis is implicated. If the patient's sample shows expression of CCPMs that is comparable to samples of non-recurrent cancer, and/or lower or higher expression than samples of recurrent cancer, then a positive prognosis is implicated.
As one approach, a prediction method can be applied to a panel of markers, for example the panel of CCPMs outlined in Table 1 and Table 2, in order to generate a predictive model. This involves the generation of a prognostic signature, comprising two or more CCPMs.
The disclosed CCPMs in Table 1 and Table 2 therefore provide a useful set of markers to generate prediction signatures for determining the prognosis of cancer, and establishing a treatment regime, or treatment modality, specific for that tumour. In particular, a positive prognosis can be used by a patient to decide to pursue standard or less invasive treatment options. A negative prognosis can be used by a patient to decide to terminate treatment or to pursue highly aggressive or experimental treatments. In addition, a patient can chose treatments based on their impact on the expression of prognostic markers (e.g., CCPMs).
Levels of CCPMs can be detected in tumour tissue, tissue proximal to the tumour, lymph node samples, blood samples, serum samples, urine samples, or faecal samples, using any suitable technique, and can include, but is not limited to, oligonucleotide probes, quantitative PCR, or antibodies raised against the markers. It will be appreciated that by analyzing the presence and amounts of expression of a plurality of CCPMs in the form of prediction signatures, and constructing a prognostic signature (e.g., as set forth in Tables 3, 4, 8A, 8B, and 9), the sensitivity and accuracy of prognosis will be increased. Therefore, multiple markers according to the present invention can be used to determine the prognosis of a cancer.
The invention includes the use of archived paraffin-embedded biopsy material for assay of the markers in the set, and therefore is compatible with the most widely available type of biopsy material. It is also compatible with several different methods of tumour tissue harvest, for example, via core biopsy or fine needle aspiration. In certain aspects, RNA is isolated from a fixed, wax-embedded cancer tissue specimen of the patient. Isolation may be performed by any technique known in the art, for example from core biopsy tissue or fine needle aspirate cells.
In one aspect, the invention relates to a method of predicting a prognosis, e.g., the likelihood of long-term survival of a cancer patient without the recurrence of cancer, comprising determining the expression level of one or more prognostic markers or their expression products in a sample obtained from the patient, normalized against the expression level of other RNA transcripts or their products in the sample, or of a reference set of RNA transcripts or their expression products. In specific aspects, the prognostic marker is one or more markers listed in Tables 1, 2, or 5, or is included as one or more of the prognostic signatures derived from the markers listed in Tables 1, 2, and 5, or the prognostic signatures listed in Tables 3, 4, 8A, 8B, or 9.
In further aspects, the expression levels of the prognostic markers or their expression products are determined, e.g., for the markers listed in Tables 1, 2, or 5, a prognostic signature derived from the markers listed in Tables 1, 2, and 5, e.g., for the prognostic signatures listed in Tables 3, 4, 8A, 8B, or 9. In another aspect, the method comprises the determination of the expression levels of a full set of prognosis markers or their expression products, e.g., for the markers listed in Tables 1, 2, or 5, or, a prognostic signature derived from the markers listed in Tables 1, 2, and 5, e.g., for the prognostic signatures listed in Tables 3, 4, 8A, 8B, or 9.
In an additional aspect, the invention relates to an array (e.g., microarray) comprising polynucleotides hybridizing to two or more markers, e.g., for the markers listed in Tables 1, 2, and 5, or a prognostic signature derived from the markers listed in Tables 1, 2, and 5, e.g., the prognostic signatures listed in Tables 3, 4, 8A, 8B, and 9. In particular aspects, the array comprises polynucleotides hybridizing to prognostic signature derived from the markers listed in Tables 1, 2, and 5, or e.g., for the prognostic signatures listed in Tables 3, 4, 8A, 8B, or 9. In another specific aspect, the array comprises polynucleotides hybridizing to the full set of markers, e.g., for the markers listed in Tables 1, 2, or 5, or, e.g., for the prognostic signatures listed in Tables 3, 4, 8A, 8B, or 9.
For these arrays, the polynucleotides can be cDNAs, or oligonucleotides, and the solid surface on which they are displayed can be glass, for example. The polynucleotides can hybridize to one or more of the markers as disclosed herein, for example, to the full-length sequences, any coding sequences, any fragments, or any complements thereof. In particular aspects, an increase or decrease in expression levels of one or more CCPM indicates a decreased likelihood of long-term survival, e.g., due to cancer recurrence, while a lack of an increase or decrease in expression levels of one or more CCPM indicates an increased likelihood of long-term survival without cancer recurrence.
Drosophila)
General Approaches to Prognostic Marker Detection
The following approaches are non-limiting methods that can be used to detect the proliferation markers, including CCPM family members: microarray approaches using oligonucleotide probes selective for a CCPM; real-time qPCR on tumour samples using CCPM specific primers and probes; real-time qPCR on lymph node, blood, serum, faecal, or urine samples using CCPM specific primers and probes; enzyme-linked immunological assays (ELISA); immunohistochemistry using anti-marker antibodies; and analysis of array or qPCR data using computers.
Other useful methods include northern blotting and in situ hybridization (Parker and Barnes, Methods in Molecular Biology 106: 247-283 (1999)); RNase protection assays (Hod, BioTechniques 13: 852-854 (1992)); reverse transcription polymerase chain reaction (RT-PCR; Weis et al., Trends in Genetics 8: 263-264 (1992)); serial analysis of gene expression (SAGE; Velculescu et al., Science 270: 484-487 (1995); and Velculescu et al., Cell 88: 243-51 (1997)), MassARRAY technology (Sequenom, San Diego, Calif.), and gene expression analysis by massively parallel signature sequencing (MPSS; Brenner et al., Nature Biotechnology 18: 630-634 (2000)). Alternatively, antibodies may be employed that can recognize specific complexes, including DNA duplexes, RNA duplexes, and DNA-RNA hybrid duplexes or DNA-polypeptide duplexes.
Primary data can be collected and fold change analysis can be performed, for example, by comparison of marker expression levels in tumour tissue and non-tumour tissue; by comparison of marker expression levels to levels determined in recurring tumours and non-recurring tumours; by comparison of marker expression levels to levels determined in tumours with or without metastasis; by comparison of marker expression levels to levels determined in differently staged tumours; or by comparison of marker expression levels to levels determined in cells with different levels of proliferation. A negative or positive prognosis is determined based on this analysis. Further analysis of tumour marker expression includes matching those markers exhibiting increased or decreased expression with expression profiles of known colorectal tumours to provide a prognosis.
A threshold for concluding that expression is increased will be dependent on the particular marker and also the particular predictive model that is to be applied. The threshold is generally set to achieve the highest sensitivity and selectivity with the lowest error rate, although variations may be desirable for a particular clinical situation. The desired threshold is determined by analysing a population of sufficient size taking into account the statistical variability of any predictive model and is calculated from the size of the sample used to produce the predictive model. The same applies for the determination of a threshold for concluding that expression is decreased. It can be appreciated that other thresholds, or methods for establishing a threshold, for concluding that increased or decreased expression has occurred can be selected without departing from the scope of this invention.
It is also possible that a prediction model may produce as it's output a numerical value, for example a score, likelihood value or probability. In these instances, it is possible to apply thresholds to the results produced by prediction models, and in these cases similar principles apply as those used to set thresholds for expression values.
Once the expression level, or output of a prediction model, of a predictive signature in a tumour sample has been obtained, the likelihood of the cancer recurring can then be determined.
From the markers identified, prognostic signatures comprising one or more CCPMs can be used to determine the prognosis of a cancer, by comparing the expression level of the one or more markers to the disclosed prognostic signature. By comparing the expression of one or more of the CCPMs in a tumour sample with the disclosed prognostic signature, the likelihood of the cancer recurring can be determined. The comparison of expression levels of the prognostic signature to establish a prognosis can be done by applying a predictive model as described previously.
Determining the likelihood of the cancer recurring is of great value to the medical practitioner. A high likelihood of re-occurrence means that a longer or higher dose treatment should be given, and the patient should be more closely monitored for signs of recurrence of the cancer. An accurate prognosis is also of benefit to the patient. It allows the patient, along with their partners, family, and friends to also make decisions about treatment, as well as decisions about their future and lifestyle changes. Therefore, the invention also provides for a method establishing a treatment regime for a particular cancer based on the prognosis established by matching the expression of the markers in a tumour sample with the differential expression signature.
It will be appreciated that the marker selection, or construction of a prognostic signature, does not have to be restricted to the CCPMs disclosed in Tables 1, 2, or 5, herein, or the prognostic signatures disclosed in Tables 3, 4, 8A, 8B, and 9, but could involve the use of one or more CCPMs from the disclosed signatures, or a new signature may be established using CCPMs selected from the disclosed marker lists. The requirement of any signature is that it predicts the likelihood of recurrence with enough accuracy to assist a medical practitioner to establish a treatment regime.
Reverse Transcription PCR (RT-PCR)
Of the techniques listed above, the most sensitive and most flexible quantitative method is RT-PCR, which can be used to compare RNA levels in different sample populations, in normal and tumour tissues, with or without drug treatment, to characterize patterns of expression, to discriminate between closely related RNAs, and to analyze RNA structure.
For RT-PCR, the first step is the isolation of RNA from a target sample. The starting material is typically total RNA isolated from human tumours or tumour cell lines, and corresponding normal tissues or cell lines, respectively. RNA can be isolated from a variety of samples, such as tumour samples from breast, lung, colon (e.g., large bowel or small bowel), colorectal, gastric, esophageal, anal, rectal, prostate, brain, liver, kidney, pancreas, spleen, thymus, testis, ovary, uterus, etc., tissues, from primary tumours, or tumour cell lines, and from pooled samples from healthy donors. If the source of RNA is a tumour, RNA can be extracted, for example, from frozen or archived paraffin-embedded and fixed (e.g., formalin-fixed) tissue samples.
The first step in gene expression profiling by RT-PCR is the reverse transcription of the RNA template into cDNA, followed by its exponential amplification in a PCR reaction. The two most commonly used reverse transcriptases are avian myeloblastosis virus reverse transcriptase (AMV-RT) and Moloney murine leukaemia virus reverse transcriptase (MMLV-RT). The reverse transcription step is typically primed using specific primers, random hexamers, or oligo-dT primers, depending on the circumstances and the goal of expression profiling. For example, extracted RNA can be reverse-transcribed using a GeneAmp RNA PCR kit (Perkin Elmer, CA, USA), following the manufacturer's instructions. The derived cDNA can then be used as a template in the subsequent PCR reaction.
Although the PCR step can use a variety of thermostable DNA-dependent DNA polymerases, it typically employs the Taq DNA polymerase, which has a 5′-3′ nuclease activity but lacks a 3′-5′ proofreading endonuclease activity. Thus, TaqMan (q) PCR typically utilizes the 5′ nuclease activity of Taq or Tth polymerase to hydrolyze a hybridization probe bound to its target amplicon, but any enzyme with equivalent 5′ nuclease activity can be used.
Two oligonucleotide primers are used to generate an amplicon typical of a PCR reaction. A third oligonucleotide, or probe, is designed to detect nucleotide sequence located between the two PCR primers. The probe is non-extendible by Taq DNA polymerase enzyme, and is labeled with a reporter fluorescent dye and a quencher fluorescent dye. Any laser-induced emission from the reporter dye is quenched by the quenching dye when the two dyes are located close together as they are on the probe. During the amplification reaction, the Taq DNA polymerase enzyme cleaves the probe in a template-dependent manner. The resultant probe fragments disassociate in solution, and signal from the released reporter dye is free from the quenching effect of the second fluorophore. One molecule of reporter dye is liberated for each new molecule synthesized, and detection of the unquenched reporter dye provides the basis for quantitative interpretation of the data.
TaqMan RT-PCR can be performed using commercially available equipment, such as, for example, ABI PRISM 7700 Sequence Detection System (Perkin-Elmer-Applied Biosystems, Foster City, Calif., USA), or Lightcycler (Roche Molecular Biochemicals, Mannheim, Germany). In a preferred embodiment, the 5′ nuclease procedure is run on a real-time quantitative PCR device such as the ABI PRISM 7700tam Sequence Detection System. The system consists of a thermocycler, laser, charge-coupled device (CCD), camera, and computer. The system amplifies samples in a 96-well format on a thermocycler. During amplification, laser-induced fluorescent signal is collected in real-time through fibre optics cables for all 96 wells, and detected at the CCD. The system includes software for running the instrument and for analyzing the data.
5′ nuclease assay data are initially expressed as Ct, or the threshold cycle. As discussed above, fluorescence values are recorded during every cycle and represent the amount of product amplified to that point in the amplification reaction. The point when the fluorescent signal is first recorded as statistically significant is the threshold cycle.
To minimize errors and the effect of sample-to-sample variation, RT-PCR is usually performed using an internal standard. The ideal internal standard is expressed at a constant level among different tissues, and is unaffected by the experimental treatment. RNAs most frequently used to normalize patterns of gene expression are mRNAs for the housekeeping genes glyceraldehyde-3-phosphate-dehydrogenase (GAPDH) and-actin.
Real-time Quantitative PCR (qPCR)
A more recent variation of the RT-PCR technique is the real time quantitative PCR, which measures PCR product accumulation through a dual-labeled fluorigenic probe (i.e., TaqMan probe). Real time PCR is compatible both with quantitative competitive PCR and with quantitative comparative PCR. The former uses an internal competitor for each target sequence for normalization, while the latter uses a normalization gene contained within the sample, or a housekeeping gene for RT-PCR. Further details are provided, e.g., by Held et al., Genome Research 6: 986-994 (1996).
Expression levels can be determined using fixed, paraffin-embedded tissues as the RNA source. According to one aspect of the present invention, PCR primers and probes are designed based upon intron sequences present in the gene to be amplified. In this embodiment, the first step in the primer/probe design is the delineation of intron sequences within the genes. This can be done by publicly available software, such as the DNA BLAT software developed by Kent, W. J., Genome Res. 12 (4): 656-64 (2002), or by the BLAST software including its variations. Subsequent steps follow well established methods of PCR primer and probe design.
In order to avoid non-specific signals, it is useful to mask repetitive sequences within the introns when designing the primers and probes. This can be easily accomplished by using the Repeat Masker program available on-line through the Baylor College of Medicine, which screens DNA sequences against a library of repetitive elements and returns a query sequence in which the repetitive elements are masked. The masked sequences can then be used to design primer and probe sequences using any commercially or otherwise publicly available primer/probe design packages, such as Primer Express (Applied Biosystems); MGB assay-by-design (Applied Biosystems); Primer3 (Steve Rozen and Helen J. Skaletsky (2000) Primer3 on the WWW for general users and for biologist programmers in: Krawetz S, Misener S (eds) Bioinformatics Methods and Protocols: Methods in Molecular Biology. Humana Press, Totowa, N.J., pp 365-386).
The most important factors considered in PCR primer design include primer length, melting temperature (Tm), and G/C content, specificity, complementary primer sequences, and 3′ end sequence. In general, optimal PCR primers are generally 17-30 bases in length, and contain about 20-80%, such as, for example, about 50-60% G+C bases. Melting temperatures between 50 and 80° C., e.g., about 50 to 70° C., are typically preferred. For further guidelines for PCR primer and probe design see, e.g., Dieffenbach, C. W. et al., General Concepts for PCR Primer Design in: PCR Primer, A Laboratory Manual, Cold Spring Harbor Laboratory Press, New York, 1995, pp. 133-155; Innis and Gelfand, Optimization of PCRs in: PCR Protocols, A Guide to Methods and Applications, CRC Press, London, 1994, pp. 5-11; and Plasterer, T. N. Primerselect: Primer and probe design. Methods Mol. Biol. 70: 520-527 (1997), the entire disclosures of which are hereby expressly incorporated by reference.
Microarray Analysis
Differential expression can also be identified, or confirmed using the microarray technique. Thus, the expression profile of CCPMs can be measured in either fresh or paraffin-embedded tumour tissue, using microarray technology. In this method, polynucleotide sequences of interest (including cDNAs and oligonucleotides) are plated, or arrayed, on a microchip substrate. The arrayed sequences (i.e., capture probes) are then hybridized with specific polynucleotides from cells or tissues of interest (i.e., targets). Just as in the RT-PCR method, the source of RNA typically is total RNA isolated from human tumours or tumour cell lines, and corresponding normal tissues or cell lines. Thus RNA can be isolated from a variety of primary tumours or tumour cell lines. If the source of RNA is a primary tumour, RNA can be extracted, for example, from frozen or archived formalin fixed paraffin-embedded (FFPE) tissue samples and fixed (e.g., formalin-fixed) tissue samples, which are routinely prepared and preserved in everyday clinical practice.
In a specific embodiment of the microarray technique, PCR amplified inserts of cDNA clones are applied to a substrate. The substrate can include up to 1, 2, 5, 10, 15, 20, 25, 30, 35, 40, 45, 50, or 75 nucleotide sequences. In other aspects, the substrate can include at least 10,000 nucleotide sequences. The microarrayed sequences, immobilized on the microchip, are suitable for hybridization under stringent conditions. As other embodiments, the targets for the microarrays can be at least 50, 100, 200, 400, 500, 1000, or 2000 bases in length; or 50-100, 100-200, 100-500, 100-1000, 100-2000, or 500-5000 bases in length. As further embodiments, the capture probes for the microarrays can be at least 10, 15, 20, 25, 50, 75, 80, or 100 bases in length; or 10-15, 10-20, 10-25, 10-50, 10-75, 10-80, or 20-80 bases in length.
Fluorescently labeled cDNA probes may be generated through incorporation of fluorescent nucleotides by reverse transcription of RNA extracted from tissues of interest. Labeled cDNA probes applied to the chip hybridize with specificity to each spot of DNA on the array. After stringent washing to remove non-specifically bound probes, the chip is scanned by confocal laser microscopy or by another detection method, such as a CCD camera. Quantitation of hybridization of each arrayed element allows for assessment of corresponding mRNA abundance. With dual colour fluorescence, separately labeled cDNA probes generated from two sources of RNA are hybridized pairwise to the array. The relative abundance of the transcripts from the two sources corresponding to each specified gene is thus determined simultaneously. An exemplary protocol for this is described in detail in Example 4.
The miniaturized scale of the hybridization affords a convenient and rapid evaluation of the expression pattern for large numbers of genes. Such methods have been shown to have the sensitivity required to detect rare transcripts, which are expressed at a few copies per cell, and to reproducibly detect at least approximately two-fold differences in the expression levels (Schena et al., Proc. Natl. Acad. Sci. USA 93 (2): 106-149 (1996)). Microarray analysis can be performed by commercially available equipment, following manufacturer's protocols, such as by using the Affymetrix GenChip technology, Illumina microarray technology or Incyte's microarray technology. The development of microarray methods for large-scale analysis of gene expression makes it possible to search systematically for molecular markers of cancer classification and outcome prediction in a variety of tumour types.
RNA Isolation, Purification, and Amplification
General methods for mRNA extraction are well known in the art and are disclosed in standard textbooks of molecular biology, including Ausubel et al., Current Protocols of Molecular Biology, John Wiley and Sons (1997). Methods for RNA extraction from paraffin embedded tissues are disclosed, for example, in Rupp and Locker, Lab Invest. 56: A67 (1987), and De Sandres et al., BioTechniques 18: 42044 (1995). In particular, RNA isolation can be performed using purification kit, buffer set, and protease from commercial manufacturers, such as Qiagen, according to the manufacturer's instructions. For example, total RNA from cells in culture can be isolated using Qiagen RNeasy mini-columns. Other commercially available RNA isolation kits include MasterPure Complete DNA and RNA Purification Kit (EPICENTRE (D, Madison, Wis.), and Paraffin Block RNA Isolation Kit (Ambion, Inc.). Total RNA from tissue samples can be isolated using RNA Stat-60 (Tel-Test). RNA prepared from tumour can be isolated, for example, by cesium chloride density gradient centrifugation.
The steps of a representative protocol for profiling gene expression using fixed, paraffin-embedded tissues as the RNA source, including mRNA isolation, purification, primer extension and amplification are given in various published journal articles (for example: T. E. Godfrey et al. J. Molec. Diagnostics 2: 84-91 (2000); K. Specht et al., Am. J. Pathol. 158: 419-29 (2001)). Briefly, a representative process starts with cutting about 10 μm thick sections of paraffin-embedded tumour tissue samples. The RNA is then extracted, and protein and DNA are removed. After analysis of the RNA concentration, RNA repair and/or amplification steps may be included, if necessary, and RNA is reverse transcribed using gene specific promoters followed by RT-PCR. Finally, the data are analyzed to identify the best treatment option(s) available to the patient on the basis of the characteristic gene expression pattern identified in the tumour sample examined.
Immunohistochemistry and Proteomics
Immunohistochemistry methods are also suitable for detecting the expression levels of the proliferation markers of the present invention. Thus, antibodies or antisera, preferably polyclonal antisera, and most preferably monoclonal antibodies specific for each marker, are used to detect expression. The antibodies can be detected by direct labeling of the antibodies themselves, for example, with radioactive labels, fluorescent labels, hapten labels such as, biotin, or an enzyme such as horse radish peroxidase or alkaline phosphatase. Alternatively, unlabeled primary antibody is used in conjunction with a labeled secondary antibody, comprising antisera, polyclonal antisera or a monoclonal antibody specific for the primary antibody. Immunohistochemistry protocols and kits are well known in the art and are commercially available.
Proteomics can be used to analyze the polypeptides present in a sample (e.g., tissue, organism, or cell culture) at a certain point of time. In particular, proteomic techniques can be used to assess the global changes of polypeptide expression in a sample (also referred to as expression proteomics). Proteomic analysis typically includes: (1) separation of individual polypeptides in a sample by 2-D gel electrophoresis (2-D PAGE); (2) identification of the individual polypeptides recovered from the gel, e.g., by mass spectrometry or N-terminal sequencing, and (3) analysis of the data using bioinformatics. Proteomics methods are valuable supplements to other methods of gene expression profiling, and can be used, alone or in combination with other methods, to detect the products of the proliferation markers of the present invention.
Once the expression level of one or more prognostic markers in a tumour sample has been assessed the likelihood of the cancer recurring can then be determined. The inventors have identified a number of markers that are differentially expressed in non-recurring colorectal cancers compared to recurring colorectal cancers in patient data sets. The markers are set out in Tables 1, 2, and 9, in the examples below.
Selection of Differentially Expressed Genes.
An early approach to the selection of genes deemed significant involved simply looking at the “fold change” of a given gene between the two groups of interest. While this approach hones in on genes that seem to change the most spectacularly, consideration of basic statistics leads one to realize that if the variance (or noise level) is quite high (as is often seen in microarray experiments), then seemingly large fold-change can happen frequently by chance alone.
Microarray experiments, such as those described here, typically involve the simultaneous measurement of thousands of genes. If one is comparing the expression levels for a particular gene between two groups (for example recurrent and non-recurrent tumours), the typical tests for significance (such as the t-test) are not adequate. This is because, in an ensemble of thousands of experiments (in this context each gene constitutes an “experiment”), the probability of at least one experiment passing the usual criteria for significance by chance alone is essentially unity. In a test for significance, one typically calculates the probability that the “null hypothesis” is correct. In the case of comparing two groups, the null hypothesis is that there is no difference between the two groups. If a statistical test produces a probability for the null hypothesis below some threshold (usually 0.05 or 0.01), it is stated that we can reject the null hypothesis, and accept the hypothesis that the two groups are significantly different. Clearly, in such a test, a rejection of the null hypothesis by chance alone could be expected 1 in 20 times (or 1 in 100). The use of t-tests, or other similar statistical tests for significance, fail in the context of microarrays, producing far too many false positives (or type 1 errors)
In this type of situation, where one is testing multiple hypotheses at the same time, one applies typical multiple comparison procedures, such as the Bonferroni Method (43). However such tests are too conservative for most microarray experiments, resulting in too many false negative (type II) errors.
A more recent approach is to do away with attempting to apply a probability for a given test being significant, and establish a means for selecting a subset of experiments, such that the expected proportion of Type I errors (or false discovery rate; 47) is controlled for. It is this approach that has been used in this investigation, through various implementations, namely the methods provided with BRB Array Tools (48), and the limma (11, 42) package of Bioconductor (that uses the R statistical environment; 10, 39).
General Methodology for Data Mining: Generation of Prognostic Signatures
Data Mining is the term used to describe the extraction of “knowledge”, in other words the “know-how”, or predictive ability from (usually) large volumes of data (the dataset). This is the approach used in this study to generate prognostic signatures. In the case of this study the “know-how” is the ability to accurately predict prognosis from a given set of gene expression measurements, or “signature” (as described generally in this section and in more detail in the examples section).
The specific details used for the methods used in this study are described in Examples 17-20. However, application of any of the data mining methods (both those described in the Examples, and those described here) can follow this general protocol.
Data mining (49), and the related topic machine learning (40) is a complex, repetitive mathematical task that involves the use of one or more appropriate computer software packages (see below). The use of software is advantageous on the one hand, in that one does not need to be completely familiar with the intricacies of the theory behind each technique in order to successfully use data mining techniques, provided that one adheres to the correct methodology. The disadvantage is that the application of data mining can often be viewed as a “black box”: one inserts the data and receives the answer. How this is achieved is often masked from the end-user (this is the case for many of the techniques described, and can often influence the statistical method chosen for data mining. For example, neural networks and support vector machines have a particularly complex implementation that makes it very difficult for the end user to extract out the “rules” used to produce the decision. On the other hand, k-nearest neighbours and linear discriminant analysis have a very transparent process for decision making that is not hidden from the user.
There are two types of approach used in data mining: supervised and unsupervised approaches. In the supervised approach, the information that is being linked to the data is known, such as categorical data (e.g. recurrent vs. non recurrent tumours). What is required is the ability to link the observed response (e.g. recurrence vs. non-recurrence) to the input variables. In the unsupervised approach, the classes within the dataset are not known in advance, and data mining methodology is employed to attempt to find the classes or structure within the dataset.
In the present example the supervised approach was used and is discussed in detail here, although it will be appreciated that any of the other techniques could be used.
The overall protocol involves the following steps:
The methods can be by first performing the step of data mining process (above), and then applying the appropriate known software packages. Further description of the process of data mining is described in detail in many extremely well-written texts. (49)
There are many variations of these methodologies that can be explored (49), and many new methodologies are constantly being defined and developed. It will be appreciated that any one of these methodologies can be applied in order to obtain an acceptable result. Particular care must be taken to avoid overfitting, by ensuring that all results are tested via a comprehensive validation scheme.
Validation
Application of any of the prediction methods described involves both training and cross-validation (43, 55) before the method can be applied to new datasets (such as data from a clinical trial). Training involves taking a subset of the dataset of interest (in this case gene expression measurements from colorectal tumours), such that it is stratified across the classes that are being tested for (in this case recurrent and non-recurrent tumours). This training set is used to generate a prediction model (defined above), which is tested on the remainder of the data (the testing set).
It is possible to alter the parameters of the prediction model so as to obtain better performance in the testing set, however, this can lead to the situation known as overfitting, where the prediction model works on the training dataset but not on any external dataset. In order to circumvent this, the process of validation is followed. There are two major types of validation typically applied, the first (hold-out validation) involves partitioning the dataset into three groups: testing, training, and validation. The validation set has no input into the training process whatsoever, so that any adjustment of parameters or other refinements must take place during application to the testing set (but not the validation set). The second major type is cross-validation, which can be applied in several different ways, described below.
There are two main sub-types of cross-validation: K-fold cross-validation, and leave-one-out cross-validation
K-fold cross-validation: The dataset is divided into K subsamples, each subsample containing approximately the same proportions of the class groups as the original. In each round of validation, one of the K subsamples is set aside, and training is accomplished using the remainder of the dataset. The effectiveness of the training for that round is gauged by how correctly the classification of the left-out group is. This procedure is repeated K-times, and the overall effectiveness ascertained by comparison of the predicted class with the known class.
Leave-one-out cross-validation: A commonly used variation of K-fold cross validation, in which K=n, where n is the number of samples.
Combinations of CCPMS, such as those described above in Tables 1 and 2, can be used to construct predictive models for prognosis.
Prognostic Signatures
Prognostic signatures, comprising one or more of these markers, can be used to determine the outcome of a patient, through application of one or more predictive models derived from the signature. In particular, a clinician or researcher can determine the differential expression (e.g., increased or decreased expression) of the one or more markers in the signature, apply a predictive model, and thereby predict the negative prognosis, e.g., likelihood of disease relapse, of a patient, or alternatively the likelihood of a positive prognosis (continued remission).
A set of prognostic signatures have been developed. In the first instance, there are two signatures developed by cross-comparison of predictive ability between two datasets: the set of microarray experiments encompassing the German colorectal cancer samples, and the set of microarray experiments encompassing the New Zealand samples (discussed in example 6). In the second instance there has been an exhaustive statistical search for effective signatures based solely on the German dataset (discussed in example 17).
As described in Example 6 below, a prognostic signature comprising 19 genes has been established from a set of colorectal samples from Germany (Table 4). Another prognostic signature, of 22 genes, has also been established from samples of colorectal tumours from patients in New Zealand (Table 3). By obtaining a patient sample (e.g., tumour sample), and matching the expression levels of one or more markers in the sample to the differential expression profile, the likelihood of the cancer recurring can be determined.
In certain aspects, this invention provides methods for determining the prognosis of a cancer, comprising: (a) providing a sample of the cancer; (b) detecting the expression level of a CCPM family member in said sample; and (c) determining the prognosis of the cancer. In one aspect, the cancer is colorectal cancer.
In other aspects, the invention includes a step of detecting the expression level of a CCPM mRNA. In other aspects, the invention includes a step of detecting the expression level of a CCPM polypeptide. In yet a further aspect, the invention includes a step of detecting the level of a CCPM peptide. In yet another aspect, the invention includes detecting the expression level of more than one CCPM family member in said sample. In a further aspect the CCPM is a gene associated with an immune response. In a further aspect the CCPM is selected from the markers set forth in Tables 3, 4, 8A, 8B, or 9. In a still further aspect, the CCPM is included in a signature selected from the signatures set forth in Tables 3, 4, 8A, 8B, or 9.
In a further aspect the invention comprises detecting the expression level of; WDR44, RBMS1, SACM1L, SOAT1, PBK, G3BP2, ZBTB20, ZNF410, COMMD2, PSMC1, COX10, GTF3C5, HMMR, UBE2L3, GNAS, PPP2R2A, RNASE2, SCOC PSMD9, EIF3S7, ATP2B4, and ABCC9. In a further aspect the invention comprises detecting the expression level of; CXCL10, FAS, CXCL0, TLK1, CXCL11, PBK, PSAT1, MAD2L1, CA2, GZMB, SLC4A4, DLG7, TNFRSF11A, KITLG, INDO, GBP1, CXCL13, CLCA4, and PCP4.
In still further aspects, the invention includes a method of determining a treatment regime for a cancer comprising: (a) providing a sample of the cancer; (b) detecting the expression level of a CCPM family member in said sample; (c) determining the prognosis of the cancer based on the expression level of a CCPM family member; and (d) determining the treatment regime according to the prognosis.
In still further aspects, the invention includes a device for detecting a CCPM, comprising: a substrate having a CCPM capture reagent thereon; and a detector associated with said substrate, said detector capable of detecting a CCPM associated with said capture reagent. Additional aspects include kits for detecting cancer, comprising: a substrate; a CCPM capture reagent; and instructions for use. Yet further aspects of the invention include method for detecting a CCPM using qPCR, comprising: a forward primer specific for said CCPM; a reverse primer specific for said CCPM; PCR reagents; a reaction vial; and instructions for use.
Additional aspects of this invention comprise a kit for detecting the presence of a CCPM polypeptide or peptide, comprising: a substrate having a capture agent for said CCPM polypeptide or peptide; an antibody specific for said CCPM polypeptide or peptide; a reagent capable of labeling bound antibody for said CCPM polypeptide or peptide; and instructions for use.
In yet further aspects, this invention includes a method for determining the prognosis of colorectal cancer, comprising the steps of: providing a tumour sample from a patient suspected of having colorectal cancer; measuring the presence of a CCPM polypeptide using an ELISA method. In specific aspects of this invention the CCPM of the invention is selected from the markers set forth in Tables 1, 2, 5, or 9. In still further aspects, the CCPM is included in a prognostic signature selected from the signatures set forth in Tables 3, 4, 8A, 8B, or 10.
The examples described herein are for purposes of illustrating embodiments of the invention. Other embodiments, methods, and types of analyses are within the scope of persons of ordinary skill in the molecular diagnostic arts and need not be described in detail hereon. Other embodiments within the scope of the art are considered to be part of this invention.
Two cohorts of patients were included in this study, one set from New Zealand (NZ) and the second from Germany (DE). The NZ patients were part of a prospective cohort study that included all disease stages, whereas the DE samples were selected from a tumour bank. Clinical information is shown in Table 6, while
Primary colorectal tumor samples from 149 NZ patients were obtained from patients undergoing surgery at Dunedin Hospital and Auckland Hospital between 1995-2000. Tumor samples were snap frozen in liquid nitrogen. All surgical specimens were reviewed by a single pathologist (H-S Y) and were estimated to contain an average of 85% tumor cells. Among the 149 CRC patients, 12 had metastatic disease at presentation, 35 developed recurrent disease, and 102 were disease-free after a minimum of 5-year follow up.
Primary colorectal tumor samples from DE patients were obtained from patients undergoing surgery at the Surgical Department of the Technical University of Munich between 1995-2001. A group of 55 colorectal carcinoma samples was selected from banked tumours which had been obtained fresh from surgery, snap frozen in liquid nitrogen. The samples were obtained from 11 patients with stage I cancer and 44 patients with stage II cancer. Twenty nine patients were recurrence-free and 26 patients had experienced disease recurrence after a minimum of 5-year follow up.
Tumor content ranged between 70 and 100% with an average of 87%.
121
NZ tumours: Tumours were homogenized and RNA was extracted using Tri-Reagent (Progenz, Auckland, New Zealand). The RNA was then further purified using RNeasy mini column (Qiagen, Victoria, Australia). Ten micrograms of RNA was labelled with Cy5 dUTP using the indirect amino-allyl cDNA labelling protocol.
A reference RNA from 12 different cell lines was labelled with Cy3 dUTP. The fluorescently labelled cDNA were purified using a QiaQuick PCR purification kit (Qiagen, Victoria, Australia) according to the manufacturer's protocol.
DE tumours: Tumours were homogenized and RNA was isolated using RNeasy Mini Kit (Qiagen, Hilden, Germany). cRNA preparation was performed as described previously (9), purified on RNeasy Columns (Qiagen, Hilden, Germany), and eluted in 55 μl of water. Fifteen micrograms of cRNA was fragmented for 35 minutes at 95° C. and double stranded cDNA was synthesized with a oligo-dT-T7 primer (Eurogentec, Koln, Germany) and transcribed using the Promega RiboMax T7-kit (Promega, Madison, Wis.) and Biotin-NTP labelling mix (Loxo, Dossenheim, Germany).
NZ tumours: Hybridisation of the labelled target cDNA was performed using MWG Human 30K Array oligonucleotides printed on epoxy coated slides. Slides were blocked with 1% BSA and the hybridisation was done in pre-hybridisation buffer at 42° C. for at least 12 hours followed by a high stringency wash. Slides were scanned with a GenePix Microarray Scanner and data was analyzed using GenePix Pro 4.1 Microarray Acquisition and Analysis Software (Axon, Calif.).
DE tumours: cRNA was mixed with B2-control oligonucleotide (Affymetrix, Santa Clara, Calif.), eukaryotic hybridization controls (Affymetrix, Santa Clara, Calif.), herring sperm (Promega, Madison, Wis.), buffer and BSA to a final volume of 300 μl and hybridized to one microarray chip (Affymetrix, Santa Clara, Calif.) for 16 hours at 45° C. Washing steps and incubation with streptavidin (Roche, Mannheim, Germany), biotinylated goat-anti streptavidin antibody (Serva, Heidelberg, Germany), goat-IgG (Sigma, Taufkirchen, Germany), and streptavidin-phycoerythrin (Molecular Probes, Leiden, Netherlands) was performed in an Affymetrix Fluidics Station according to the manufacturer's protocol. The arrays were then scanned with a HP-argon-ion laser confocal microscope and the digitized image data were processed using the Affymetrix® Microarray Suite 5.0 Software.
NZ data: Data pre-processing and normalization was performed in the R computing environment (10). A log2 transformation was applied to the foreground intensities from each channel of each array. Data from each spot was used on a per array basis to perform print-tip loss normalization via the limma package (11) from the Bioconductor suite of analysis tools (12). Scale normalization (13) was then used to standardize the distribution of log intensity ratios across arrays. Post-normalization cluster analysis revealed the presence of a gene-specific print-run effect present in the data. Analysis of variance (ANOVA) normalization was used to estimate and remove print run effects from the data for each gene. Replicate array data was available for 46 of the 149 samples. Cluster analysis of the entire data set indicated that the duplicate arrays clustered well with each other suggesting internal consistency of the array platform. Genes with low intensity, large differences between replicates (mean log2 difference between duplicates higher than 0.5), and unknown proteins were removed from the data set. After the initial normalization procedure, a subset of 10,318 genes was chosen for further analysis.
DE data: All Affymetrix U133A GeneChips passed quality control to eliminate scans with abnormal characteristics, that is, abnormal low or high dynamic range, high perfect match saturation, high pixel noise, grid misalignment problems, and low mean signal to noise ratio. Background correction and normalization were performed in the R computing environment (10, 40). Background corrected and normalized expression measures from probe level data (cel-files) were obtained using the robust multi-array average function (14) implemented in the Bioconductor package affy.
Data analysis was performed using the BRB Array-Tools package. Gene selection was performed using a random variance model t-test. In the DE data, 318 genes were found to be differentially expressed when using a significance threshold of 0.001. As most of the differentially expressed genes exhibited relatively small changes in expression, a condition requiring the mean log2 fold change between the two classes to be higher than 1.1 was added to the gene selection process for the DE data. Gene-based prognostic signatures were produced using leave one out cross validation (LOOCV) in each of the NZ and DE data sets. To avoid the problem of over-fitting, both the gene selection and signature construction were performed during each LOOCV iteration. After LOOCV, the prediction rate was estimated by the fraction of samples correctly predicted. In order to find a gene set that could make the best prediction for unknown samples, different t-test thresholds using a random variance model were investigated in conjunction with six classification methods: compound covariate classifier (CCP), diagonal linear discriminant analysis (DLD), 3-nearest neighbours (3-NN), 1-nearest neighbours (1-NN), nearest centroid (NC), and support vector machines (SVM).
To establish the validity of the NZ and DE prognosis signatures, reciprocal validation was performed, with the NZ signature validated using the DE data set, and vice versa. To test the NZ genes, probes relating to the 22 genes from the NZ signature were identified in the DE data, and LOOCV was used to assess the performance of a signature for the DE samples, based only on these probes. Similarly, probes relating to the 19 genes in the DE signature were identified in the NZ data and LOOCV was used to assess the performance of a signature for the NZ samples. In both cases a significance threshold of 0.999 was used to ensure that all genes were used in each LOOCV iteration. Differences between the platforms (in particular, log-ratio data versus log-intensity data) meant that direct application of a prediction rule across data sets was not feasible. The consequence of this is that only the gene sets, and not the prediction rules used, can be generalized to new samples. The significance of the LOOCV prediction results was calculated by permuting the class labels of the samples and finding the proportion of times that the permuted data resulted in a higher LOOCV prediction rate than that obtained for the unpermuted data. All permutation analysis involved 2000 permutations, with small P-values indicating that prediction results were unlikely to be due to chance.
Kaplan-Meier survival analysis for censored data was performed using the survival package within the R computing environment. Survival was defined to be “disease free survival” post surgery. For each analysis, survival curves were constructed, and the log-rank test (15) was used to assess the presence of significant differences between the curves for the two groups in question. Censoring was taken into account for both the NZ and DE data sets. For the disease-free survival data, right censoring prior to five years could only occur for non-recurrent patients as a result of either death, or the last clinical follow-up occurring at less than five years. Odds ratios and confidence intervals were produced using the epitools package for R.
Genes in the DE data which had a Pearson correlation coefficient greater than 0.75 with at least one of the four chemokines appearing in the predictor in the non-relapse group were selected for ontology analysis. Ontology was performed using DAVID (hypertext transfer protocol://apps1.niaid.nih.gov/david/).
To identify robust prognostic signatures to predict disease relapse for CRC, two independent sets of samples from NZ and DE were used to generate array expression data sets from separate series of primary tumours with clinical follow-up of five or more years. After normalization, each data set was analyzed using the same statistical methods to generate a prognostic signature, which was then validated on the alternate series of patients. As such, the DE prognostic signature was validated on the NZ data set and the NZ prognostic signature was validated on the DE data set.
DE Data Set: The BRB Array Tools class comparison procedure was used to detect probes exhibiting statistically significant differences in average intensity between relapse and non-relapse samples. The RVM (random variance model) was again used to produce p-values for each probe in the data set. In this second round, a total of 325 probes were found to be significantly differentially expressed between the two sample classes using an arbitrary significance threshold of 0.05. Note this selection of genes did not apply any fold-change threshold, and used a significance cut off of 0.05, rather than the threshold of 0.001 that was used in Example 6. The purpose of this less stringent threshold (p=0.05 instead of p=0.001) was to put forward a larger number of genes for construction of the second round of signatures (see example 17) These probes represent 270 unique genes (Table 1 and Table 2).
Explicitly, the test for significance (random variance model) comprises the following: generating a test statistic for each gene which was identical to that of a standard two sample t-test (45) except that the estimate of the pooled variance was obtained by representing the variance structure across all genes as an F-distribution, and then using the parameters, a and b, of this distribution (obtained via maximization of the empirical likelihood function) to form the following estimate of the pooled variance (see next page),
where S2 is the new estimate of the pooled variance, S2pooled is the standard estimate of pooled variance (45), n is the number of samples, and a and b are the parameters of the F-distribution (46). Based on the t-statistic formed, a t-distribution with n−2+2a degrees of freedom was used to obtain a p-value for each gene. To adjust for multiple hypothesis testing, the False Discovery Rate controlling procedure of Benjamini and Hochberg (7) was used to produce adjusted p-values for each gene. A gene was considered to have undergone significant differential expression if its adjusted p-value was less than 0.05.
In order to identify additional genes that can be used as prognostic predictors, correlation analysis was carried out using the R statistical computing software package. This analysis revealed 167 probes that had a Pearson correlation coefficient (40, 44, 45) of at least 0.8. Of these probes, 51 were already present in the set of 325 significantly differentially expressed probes, while the remaining 116 were reported as non-significant (using a 0.05 threshold for the FDR, or “false-discovery rate” (47) controlling procedure, the RVM, or rando variance model). These 116 probes represent 111 distinct genes (Table 2).
The NZ data set was generated using oligonucleotide printed microarrays. Six different signatures were constructed, with a support vector machine (SVM) using a gene selection threshold of 0.0008 yielding the highest LOOCV prediction rate, and producing a 22-gene signature (77% prediction rate, 53% sensitivity, 88% specificity; p=0.002, Tables 7, 8A, and 8B). For Tables 8A and 8B, the gene descriptions are shown in Tables 3 and 4, respectively.
§95% confidence interval
The NZ signature had an odds ratio for disease recurrence in the NZ patients of 8.4 (95% CI 3.5-21.4).
The DE data set was generated using Affymetrix arrays resulting in a 19-gene (22-probe) and 3-nearest neighbour (3-NN) signature (selection threshold 0.002, log2 fold change>1.1, 84% classification rate, 85% sensitivity, 83% specificity, p<0.0001, Tables 3, 4, 7). The DE signature had an odds ratio for recurrence in the DE patients of 24.1 (95% CI 5.3-144.7). Using Kaplan-Meier analysis, disease-free survival in NZ and DE patients was significantly different for those predicted to recur or not recur (NZ signature, p<0.0001,
To validate the NZ signature, the 22 genes were used to construct a SVM signature in the DE data set by LOOCV. A prediction rate of 71% was achieved, which was highly significant (p=0.002; Table 7). The odds ratio for recurrence in DE patients, using the NZ signature, was 5.9 (95% CI 1.6-24.5). We surmise that the reduction in prediction rate, from 77% in NZ patients to 71% in DE patients (Table 7), was due to four genes from the NZ signature not being present in the DE data. Disease-free survival for DE patients predicted to relapse, according to the NZ signature, was significantly lower than disease-free survival for patients predicted not to relapse (p=0.0049,
The DE signature was next validated by using the 19 genes to construct a 3-NN signature in the NZ data set by LOOCV. The prediction rate of 67% was again significant (p=0.046; Table 7), confirming the validity of the DE signature. The odds ratio for recurrence in NZ patients, using the DE signature, was 2.6 (95% CI 1.2-6.0). We consider that the reduction of the prediction rate was due to five genes from DE signature not being present in the NZ data set. This was confirmed when removal of these five genes from the DE data set resulted in a reduction of the LOOCV prediction rate from 84% to 67% (Table 7). Disease-free survival for NZ patients predicted to relapse, according to the DE signature, was significantly lower than disease-free survival for patients predicted not to relapse (p=0.029;
Significant differences in disease-free survival between patients predicted to relapse or not relapse were also observed within the same clinico-pathological stage (
A number of genes in the NZ signature (Table 3) including G3BP2 (16), RBMS1 (17), HMMR (18), UBE2L3 (19), GNAS (20), RNASE2 (21) and ABCC9 (22) have all been reported to be involved in cancer progression, while RBMS1 (23), EIF3S7 (24) and GTF3C5 (25) are involved in transcription or translation. PBK is a protein kinase, which is involved in the process of mitosis (26), and the only gene common to the NZ and DE signatures. Eleven of 19 genes in the DE signature (Table 4) are involved in the immune response including 4 chemokine ligands (CXCL9, CXCL10, CXCL11, CXCL13; (27)), PBK (28), INDO (29), GBP1 (30), GZMB (31), KITLG (32), and two receptors of the tumor necrosis factor family (TNFRSF11A, FAS; 33)).
Eighty six genes were found to be moderately correlated (Pearson correlation coefficient >0.75) with at least one of the four chemokine ligands in the DE data. Ontology analysis found that 39 of these 65 genes were in the category of immune response (p<10-26). This result suggests a key role for the host immune response in determining CRC recurrence.
It has been shown that the two different prognostic signatures can be used to improve the current prognosis of colorectal cancer.
For the DE signature, it was surprising and unexpected that the stage I/II samples could be used to predict stage III outcome. It was also surprising that many genes associated with recurrent disease are related to the immune response. The immune response has an important role in the progression of different cancers and T-lymphocyte infiltration in CRC patients is an indicator of good prognosis (36-38). All of the eleven immune response (Table 5) genes were down-regulated in recurrent patients which would be unexpected based on known biological mechanisms.
To further confirm these results, 4 chemokine genes were chosen for further analysis. Chemokine ligands not only reflect the activity of the immune system and mediate leukocyte recruitment but also are involved in chemotaxis, cell adhesion and motility, and angiogenesis (36). To investigate the role of the immune response genes, 86 genes co-expressed with the chemokine ligands were identified. Almost half of these genes had a Gene Ontology classification within the “immune response” category suggesting that the primary function of these genes in the recurrence process is the modulation of the immune response. Furthermore, CD4+ and CD8+ T cell antigens (CD8A, CD3, PRF1, TRA@, TRB@) or functionally related antigens, for example, major histocompatibility molecules, interferon gamma induced proteins, and IL2RB, were found in the co-expressed gene list. The activation of tumor specific CD4+ T cells and CD8+ T cells has been shown to result in tumour rejection in a mouse colorectal cancer model (37). Collectively, these findings suggest that the lymphocytes form part of a tumor-specific host response involved in minimising the spread of cells from the primary tumour.
The performance of the two prognostic signatures described above was excellent in terms of cross-validation between the two data sets. Further studies were carried out, using a purely statistical approach, to develop a range of signatures, in addition to the aforementioned, that would also predict prognosis for other data sets. One of the additional goals of these studies was to ensure that the method used to normalize the microarray data (robust multi-array average) was not exerting undue influence on the choice of genes.
The effect of normalization on feature selection was thoroughly investigated by generating gene lists from 1000 stratified sub-samples of the original set of tumours, each time removing 5 samples (i.e. 1/11 of the total number of samples) from the data set. (This is effectively the same as performing 11-fold cross-validation). A tally was made of the number of times each gene appeared in the “top-n” gene lists (i.e., top 10, top 20, top 100, and top 325). This value was termed the “top count”. Top counts were generated using three different normalization methods (40) (
Genes from the gene lists (see Table 1 and Table 2), were used to generate signatures by random sampling. The generation of samples was weighted, such that genes with higher “top count” were more likely to be selected. A range of signatures was generated, using between 2 and 55 Affymetrix probes. Signatures were selected if they exhibited >80% median classification rate, using three methods of classifiers: k-nearest neighbours, with k=1; k-nearest neighbours, with k=3; and support vector machines, with a linear kernel function, and using leave-one-out cross-validation.
On average, longer prognostic signatures were preferred over shorter signatures in terms of ability to predict prognosis for new data sets (
In addition, because, on average, longer signatures (>10 genes/signature) tended to perform better, we selected signatures with 20 or more genes/signatures from a pool of signatures with 30 or more probes/signature. It is expected that these signatures (Table 10) will perform with a classification rate of around 70% when applied to other data sets, on the basis of the results shown in
In selection of the gene signatures described here, two different statistical methods were used to characterise the signatures: k-nearest neighbours, and support vector machines. These methods are provided as packages to the R statistical software system (ref), through the packages class (ref) and e1071 (ref).
The signatures described in this document were tested as follows. In both cases, the data used to develop the prediction models for a given signature were the gene expression values (raw normalised intensities from the Affymetrix array data) for the probes corresponding to genes that comprise that signature, across both recurrent and non-recurrent samples:
Wherein in the description reference has been made to integers or components having known equivalents, such equivalents are herein incorporated as if individually set fourth.
Although the invention has been described by way of example and with reference to possible embodiments thereof, it is to be appreciated that improvements and/or modifications may be made without departing from the scope thereof.
| Number | Date | Country | Kind |
|---|---|---|---|
| 544432 | Dec 2005 | NZ | national |
This application is a Continuation under 35 U.S.C. § 1.111(a) of PCT/NZ2006/000343, International Filing Date 22 Dec. 2006, which claims the benefit of New Zealand Provisional Patent Application No. 544432 filed Dec. 23, 2005, both of which are incorporated by reference herein in their entirety.
| Number | Date | Country |
|---|---|---|
| 2008046195 | Aug 2009 | SG |
| WO2002068677 | Sep 2002 | WO |
| WO2002068677 | Sep 2002 | WO |
| WO2004090550 | Oct 2004 | WO |
| WO2005064009 | Jul 2005 | WO |
| Entry |
|---|
| Wang et al. (J Clin Oncol., May 2004, 22:1564-1571). |
| Butler et al. (Br J Cancer. Jan. 2000; 82(1): 131-135). |
| Dieffenbach et al. (Genome Res.,1993, 3: S30-S37). |
| “Affymatrix GeneChip Human Genome U133 Array Set HG-U133A” Geo Host, Mar. 11, 2002. |
| Eschrich, S et al., “Molecular staging for survival prediction of colorectal cancer patients”; Journal of Clinical Oncology (May 2005) vol. 23, No. 15, pp. 3526-3535. |
| Arango, D., et al.; “Gene expression profiling predicts recurrence in Dukes' C. colorectal cancer”; Gastroenterology (Sep. 2005) vol. 129, pp. 874-884. |
| Int'l Preliminary Exam Rp, Jul. 2, 2004, PCT/NZ2003/000045. |
| Database Medline Abstract, Database Accession No. PMID:16078572; Sichuan Da Xue Xue Bao Yi Xue Ban Jul. 2005; 36(4):503-5; Lin Man-Hua et al: “Expression and clinical significance of nm23-H1, Fas and FasL in colorectal carcinoma tissues”. |
| Goo, Jian-Jun, “Evaluative significance of the expression of Fas antigen and P53 protein in colorectal cancer for the biological characters and prognosis,” Basic Research, Reference Number XP008082293. (Best copy available.) |
| Eschrich, Steven; Molecular Staging or Survival Prediction of Colorectal Cancer Patients; Journal of Clinical Oncology; May 20, 2005; vol. 23, No. 15, pp. 3526-3535. |
| English Translation of Korean office action dated May 10, 2013; Korean Patent Appl. No. 7018020/2008. |
| Abstract with reference No. XP-002512619. |
| European Search Report dated Nov. 2, 2009; Application No. 06835771.4-1222/1977237; Pacific Edge Biotechnology Limited; reference FBP20410. |
| Number | Date | Country | |
|---|---|---|---|
| 20090181384 A1 | Jul 2009 | US |
| Number | Date | Country | |
|---|---|---|---|
| Parent | PCT/NZ2006/000343 | Dec 2006 | US |
| Child | 12214782 | US |
