Breast cancer survival and recurrence

Information

  • Patent Application
  • 20050100933
  • Publication Number
    20050100933
  • Date Filed
    June 18, 2004
    20 years ago
  • Date Published
    May 12, 2005
    19 years ago
Abstract
The invention provides for the identification and use of gene expression profiles, or patterns, with clinical relevance to breast cancer. In particular, the invention provides the identities of genes that are correlated with patient survival and breast cancer recurrence. The gene expression profiles may be embodied in nucleic acid expression, protein expression, or other expression formats and used to predict the survival of subjects afflicted with breast cancer and to predict breast cancer recurrence and. The profiles may also be used in the study and/or diagnosis of breast cancer cells and tissue, including the grading of invasive breast cancer, as well as for the study and/or determination of prognosis of a patient. When used for diagnosis or prognosis, the profiles may be used to determine the treatment of breast cancer based upon the likelihood of life expectancy and recurrence.
Description
FIELD OF THE INVENTION

The invention relates to the identification and use of gene expression profiles, or patterns, with clinical relevance to breast cancer. In particular, the invention provides the identities of genes that are correlated with patient survival and breast cancer recurrence. The gene expression profiles, whether embodied in nucleic acid expression, protein expression, or other expression formats, may be used to predict the survival of subjects afflicted with breast cancer and to predict breast cancer recurrence and. The profiles may also be used in the study and/or diagnosis of breast cancer cells and tissue, including the grading of invasive breast cancer, as well as for the study and/or determination of prognosis of a patient. When used for diagnosis or prognosis, the profiles are used to determine the treatment of breast cancer based upon the likelihood of life expectancy and recurrence.


BACKGROUND OF THE INVENTION

Breast cancer is by far the most common cancer among women. Each year, more than 180,000 and 1 million women in the U.S. and worldwide, respectively, are diagnosed with breast cancer. Breast cancer is the leading cause of death for women between ages 50-55, and is the most common non-preventable malignancy in women in the Western Hemisphere. An estimated 2,167,000 women in the United States are currently living with the disease (National Cancer Institute, Surveillance Epidemiology and End Results (NCI SEER) program, Cancer Statistics Review (CSR), www-seer.ims.nci.nih.gov/Publications/CSR1973 (1998)). Based on cancer rates from 1995 through 1997, a report from the National Cancer Institute (NCI) estimates that about 1 in 8 women in the United States (approximately 12.8 percent) will develop breast cancer during her lifetime (NCI's Surveillance, Epidemiology, and End Results Program (SEER) publication SEER Cancer Statistics Review 1973-1997). Breast cancer is the second most common form of cancer, after skin cancer, among women in the United States. An estimated 250,100 new cases of breast cancer are expected to be diagnosed in the United States in 2001. Of these, 192,200 new cases of more advanced (invasive) breast cancer are expected to occur among women (an increase of 5% over last year), 46,400 new cases of early stage (in situ) breast cancer are expected to occur among women (up 9% from last year), and about 1,500 new cases of breast cancer are expected to be diagnosed in men (Cancer Facts & Figures 2001 American Cancer Society). An estimated 40,600 deaths (40,300 women, 400 men) from breast cancer are expected in 2001. Breast cancer ranks second only to lung cancer among causes of cancer deaths in women. Nearly 86% of women who are diagnosed with breast cancer are likely to still be alive five years later, though 24% of them will die of breast cancer after 10 years, and nearly half (47%) will die of breast cancer after 20 years.


Every woman is at risk for breast cancer. Over 70 percent of breast cancers occur in women who have no identifiable risk factors other than age (U.S. General Accounting Office. Breast Cancer, 1971-1991: Prevention, Treatment and Research. GAO/PEMD-92-12; 1991). Only 5 to 10% of breast cancers are linked to a family history of breast cancer (Henderson I C, Breast Cancer. In: Murphy G P, Lawrence W L, Lenhard R E (eds). Clinical Oncology. Atlanta, Ga.: American Cancer Society; 1995:198-219).


Each breast has 15 to 20 sections called lobes. Within each lobe are many smaller lobules. Lobules end in dozens of tiny bulbs that can produce milk. The lobes, lobules, and bulbs are all linked by thin tubes called ducts. These ducts lead to the nipple in the center of a dark area of skin called the areola. Fat surrounds the lobules and ducts. There are no muscles in the breast, but muscles lie under each breast and cover the ribs. Each breast also contains blood vessels and lymph vessels. The lymph vessels carry colorless fluid called lymph, and lead to the lymph nodes. Clusters of lymph nodes are found near the breast in the axilla (under the arm), above the collarbone, and in the chest.


Breast tumors can be either benign or malignant. Benign tumors are not cancerous, they do not spread to other parts of the body, and are not a threat to life. They can usually be removed, and in most cases, do not come back. Malignant tumors are cancerous, and can invade and damage nearby tissues and organs. Malignant tumor cells may metastasize, entering the bloodstream or lymphatic system. When breast cancer cells metastasize outside the breast, they are often found in the lymph nodes under the arm (axillary lymph nodes). If the cancer has reached these nodes, it means that cancer cells may have spread to other lymph nodes or other organs, such as bones, liver, or lungs.


Major and intensive research has been focused on early detection, treatment and prevention. This has included an emphasis on determining the presence of precancerous or cancerous ductal epithelial cells. These cells are analyzed, for example, for cell morphology, for protein markers, for nucleic acid markers, for chromosomal abnormalities, for biochemical markers, and for other characteristic changes that would signal the presence of cancerous or precancerous cells. This has led to various molecular alterations that have been reported in breast cancer, few of which have been well characterized in human clinical breast specimens. Molecular alterations include presence/absence of estrogen and progesterone steroid receptors, HER-2 expression/amplification (Mark H F, et al. HER-2/neu gene amplification in stages I-IV breast cancer detected by fluorescent in situ hybridization. Genet Med; 1(3):98-103 1999), Ki-67 (an antigen that is present in all stages of the cell cycle except GO and used as a marker for tumor cell proliferation, and prognostic markers (including oncogenes, tumor suppressor genes, and angiogenesis markers) like p53, p27, Cathepsin D, pS2, multi-drug resistance (MDR) gene, and CD31.


van't Veer et al. (Nature 415:530-536, 2002) describe gene expression profiling of clinical outcome in breast cancer. They identified genes expressed in breast cancer tumors, the expression levels of which correlated either with patients afflicted with distant metastases within 5 years or with patients that remained metastasis-free after at least 5 years.


Ramaswamy et al. (Nature Genetics 33:49-54, 2003) describe the identification of a molecular signature of metastasis in primary solid tumors. The genes of the signature were identified based on gene expression profiles of 12 metastatic adenocarcinoma nodules of diverse origin (lung, breast, prostate, colorectal, uterus) compared to expression profiles of 64 primary adenocarcinomas representing the same spectrum of tumor types from different individuals. A 128 gene set was identified.


Both of the above described approaches, however, utilize heterogeneous populations of cells found in a tumor sample to obtain information on gene expression patterns. The use of such populations may result in the inclusion or exclusion of multiple genes that are differentially expressed in cancer cells. The gene expression patterns observed by the above described approaches may thus provide little confidence that the differences in gene expression are meaningfully associated with breast cancer recurrence or survival.


Citation of documents herein is not intended as an admission that any is pertinent prior art. All statements as to the date or representation as to the contents of documents is based on the information available to the applicant and does not constitute any admission as to the correctness of the dates or contents of the documents.


BRIEF SUMMARY OF THE INVENTION

The present invention relates to the identification and use of gene expression patterns (or profiles or “signatures”) which are clinically relevant to breast cancer. In particular, the identities of genes that are correlated with patient survival and breast cancer recurrence are provided. The gene expression profiles, whether embodied in nucleic acid expression, protein expression, or other expression formats, may be used to predict survival of subjects afflicted with breast cancer and the likelihood of breast cancer recurrence.


The invention thus provides for the identification and use of gene expression patterns (or profiles or “signatures”) which correlate with (and thus able to discriminate between) patients with good or poor survival outcomes. In one embodiment, the invention provides patterns that are able to distinguish patients with estrogen receptor (ER) positive breast tumors into those with a survival outcome poorer than that of patients with ER negative breast tumors and those with a better survival outcome than that of patients with ER positive breast tumors. These patterns are thus able to distinguish patients with ER positive breast tumors into at least two subtypes.


The invention also provides for the identification and use of gene expression patterns which correlate with the recurrence of breast cancer at the same location and/or in the form of metastases. The pattern is able to distinguish patients with breast cancer into at least those with good or poor survival outcomes.


In another aspect of the invention, the ability to identify the grade of invasive breast cancer by gene expression patterns of the invention is provided. In particular, gene expression patterns in a cell containing sample that distinguish “high-grade” (or “grade 3”) invasive breast tumors from “low-grade” (or grades “1” and “2”) invasive breast tumors are provided. The invention thus permits the distinguishing (or grading) of a subject's invasive tumors into two types which may be differentially treated based on the expected outcome associated with each type.


The present invention provides a non-subjective means for the identification of patients with breast cancer as likely to have a good or poor survival outcome by assaying for the expression patterns disclosed herein. Thus where subjective interpretation may have been previously used to determine the prognosis and/or treatment of breast cancer patients, the present invention provides objective gene expression patterns, which may used alone or in combination with subjective criteria to provide a more accurate assessment of breast cancer patient outcomes, including survival and the recurrence of cancer. The expression patterns of the invention thus provide a means to determine breast cancer prognosis. Furthermore, the expression patterns can also be used as a means to assay small, node negative tumors that are not readily assayed by other means.


The gene expression patterns comprise one or more than one gene capable of discriminating between breast cancer outcomes with significant accuracy. The gene(s) are identified as correlated with various breast cancer outcomes such that the levels of their expression are relevant to a determination of the preferred treatment protocols, of a breast cancer patient. Thus in one aspect, the invention provides a method to determine the outcome of a subject afflicted with, or suspected of having, breast cancer by assaying a cell containing sample from said subject for expression of one or more than one gene disclosed herein as correlated with breast cancer outcomes.


Gene expression patterns of the invention are identified as described below. Generally, a large sampling of the gene expression profile of a sample is obtained through quantifying the expression levels of mRNA corresponding to many genes. This profile is then analyzed to identify genes, the expression of which are positively, or negatively, correlated, with a breast cancer outcome. An expression profile of a subset of human genes may then be identified by the methods of the present invention as correlated with a particular breast cancer outcome. The use of multiple samples increases the confidence which a gene may be believed to be correlated with a particular survival outcome. Without sufficient confidence, it remains unpredictable whether a particular gene is actually correlated with a breast cancer outcome and also unpredictable whether a particular gene may be successfully used to identify the outcome for a breast cancer patient.


A profile of genes that are highly correlated with one outcome relative to another may be used to assay an sample from a subject afflicted with, or suspected of having, breast cancer to predict the outcome of the subject from whom the sample was obtained. Such an assay may be used as part of a method to determine the therapeutic treatment for said subject based upon the breast cancer outcome identified.


The correlated genes may be used singly with significant accuracy or in combination to increase the ability to accurately correlating a molecular expression phenotype with a breast cancer outcome. This correlation is a way to molecularly provide for the determination of survival outcomes as disclosed herein. Additional uses of the correlated gene(s) are in the classification of cells and tissues; determination of diagnosis and/or prognosis; and determination and/or alteration of therapy.


The ability to discriminate is conferred by the identification of expression of the individual genes as relevant and not by the form of the assay used to determine the actual level of expression. An assay may utilize any identifying feature of an identified individual gene as disclosed herein as long as the assay reflects, quantitatively or qualitatively, expression of the gene in the “transcriptome” (the transcribed fraction of genes in a genome) or the “proteome” (the translated fraction of expressed genes in a genome). Identifying features include, but are not limited to, unique nucleic acid sequences used to encode (DNA), or express (RNA), said gene or epitopes specific to, or activities of, a protein encoded by said gene. All that is required is the identity of the gene(s) necessary to discriminate between breast cancer outcomes and an appropriate cell containing sample for use in an expression assay.


In one embodiment, the invention provides for the identification of the gene expression patterns by analyzing global, or near global, gene expression from single cells or homogenous cell populations which have been dissected away from, or otherwise isolated or purified from, contaminating cells beyond that possible by a simple biopsy. Because the expression of numerous genes fluctuate between cells from different patients as well as between cells from the same patient sample, multiple data from expression of individual genes and gene expression patterns are used as reference data to generate models which in turn permit the identification of individual gene(s), the expression of which are most highly correlated with particular breast cancer outcomes.


In another aspect, the invention provides physical and methodological means for detecting the expression of gene(s) identified by the models generated by individual expression patterns. These means may be directed to assaying one or more aspect of the DNA template(s) underlying the expression of the gene(s), of the RNA used as an intermediate to express the gene(s), or of the proteinaceous product expressed by the gene(s).


In a further aspect, the gene(s) identified by a model as capable of discriminating between breast cancer outcomes may be used to identify the cellular state of an unknown sample of cell(s) from the breast. Preferably, the sample is isolated via non-invasive means. The expression of said gene(s) in said unknown sample may be determined and compared to the expression of said gene(s) in reference data of gene expression patterns correlated with breast cancer outcomes. Optionally, the comparison to reference samples may be by comparison to the model(s) constructed based on the reference samples.


One advantage provided by the present invention is that contaminating, non-breast cells (such as infiltrating lymphocytes or other immune system cells) are not present to possibly affect the genes identified or the subsequent analysis of gene expression to identify the survival outcomes of patients with breast cancer. Such contamination is present where a biopsy is used to generate gene expression profiles.


In another aspect, the invention provides the identification and use of four gene sequences the expression of which are significantly associated with tumor recurrence. Elevated expression of each one of the four gene sequences is correlated with increased likelihood of tumor recurrence and decreased patient survival. Therefore, the expression of each of these gene sequences may be used in the same manner as described herein for gene expression patterns.


The first set of sequences is that of mitotic spindle associated protein (also known as mitotic spindle coiled-coil related protein, ASTRIN or DEEPEST). Human DEEPEST protein has been characterized by Mack et al. (Proc Natl Acad Sci USA. 2001 98(25): 14434-9).


The second set of sequences is that of the “Rac GTPase activating protein 1” (RACGAP1).


The third set of sequences is that of the “zinc finger protein 145” or “PLZF” (Kruppel-like zinc finger protein, expressed in promyelocytic leukemia) which is also referred to as ZNF145.


The fourth set of sequences is that of “MS4A7” (membrane-spanning 4-domains, subfamily A, member 7).


While the present invention is described mainly in the context of human breast cancer, it may be practiced in the context of breast cancer of any animal known to be potentially afflicted by breast cancer. Preferred animals for the application of the present invention are mammals, particularly those important to agricultural applications (such as, but not limited to, cattle, sheep, horses, and other “farm animals”), animal models of breast cancer, and animals for human companionship (such as, but not limited to, dogs and cats).







DETAILED DESCRIPTION OF THE INVENTION

Definitions of terms as used herein:


A gene expression “pattern” or “profile” or “signature” refers to the relative expression of a gene between two or more breast cancer survival outcomes which is correlated with being able to distinguish between said outcomes.


A “gene” is a polynucleotide that encodes a discrete product, whether RNA or proteinaceous in nature. It is appreciated that more than one polynucleotide may be capable of encoding a discrete product. The term includes alleles and polymorphisms of a gene that encodes the same product, or a functionally associated (including gain, loss, or modulation of function) analog thereof, based upon chromosomal location and ability to recombine during normal mitosis.


The terms “correlate” or “correlation” or equivalents thereof refer to an association between expression of one or more genes and a physiologic state of a breast cell to the exclusion of one or more other state as identified by use of the methods as described herein. A gene may be expressed at higher or lower levels and still be correlated with one or more breast cancer state or outcome.


A “polynucleotide” is a polymeric form of nucleotides of any length, either ribonucleotides or deoxyribonucleotides. This term refers only to the primary structure of the molecule. Thus, this term includes double- and single-stranded DNA and RNA. It also includes known types of modifications including labels known in the art, methylation, “caps”, substitution of one or more of the naturally occurring nucleotides with an analog, and internucleotide modifications such as uncharged linkages (e.g., phosphorothioates, phosphorodithioates, etc.), as well as unmodified forms of the polynucleotide.


The term “amplify” is used in the broad sense to mean creating an amplification product can be made enzymatically with DNA or RNA polymerases. “Amplification,” as used herein, generally refers to the process of producing multiple copies of a desired sequence, particularly those of a sample. “Multiple copies” mean at least 2 copies. A “copy” does not necessarily mean perfect sequence complementarity or identity to the template sequence.


By corresponding is meant that a nucleic acid molecule shares a substantial amount of sequence identity with another nucleic acid molecule. Substantial amount means at least 95%, usually at least 98% and more usually at least 99%, and sequence identity is determined using the BLAST algorithm, as described in Altschul et al. (1990), J. Mol. Biol. 215:403-410 (using the published default setting, i.e. parameters w=4, t=17). Methods for amplifying mRNA are generally known in the art, and include reverse transcription PCR (RT-PCR) and those described in U.S. patent application Ser. No. 10/062,857 (filed on Oct. 25, 2001), as well as U.S. Provisional Patent Application 60/298,847 (filed Jun. 15, 2001) and 60/257,801 (filed Dec. 22, 2000), all of which are hereby incorporated by reference in their entireties as if fully set forth. Another method which may be used is quantitative PCR (or Q-PCR). Alternatively, RNA may be directly labeled as the corresponding cDNA by methods known in the art.


A “microarray” is a linear or two-dimensional array of preferably discrete regions, each having a defined area, formed on the surface of a solid support such as, but not limited to, glass, plastic, or synthetic membrane. The density of the discrete regions on a microarray is determined by the total numbers of immobilized polynucleotides to be detected on the surface of a single solid phase support, preferably at least about 50/cm2, more preferably at least about 100/cm2, even more preferably at least about 500/cm2, but preferably below about 1,000/cm2. Preferably, the arrays contain less than about 500, about 1000, about 1500, about 2000, about 2500, or about 3000 immobilized polynucleotides in total. As used herein, a DNA microarray is an array of oligonucleotides or polynucleotides placed on a chip or other surfaces used to hybridize to amplified or cloned polynucleotides from a sample. Since the position of each particular group of primers in the array is known, the identities of a sample polynucleotides can be determined based on their binding to a particular position in the microarray.


Because the invention relies upon the identification of genes that are over- or under-expressed, one embodiment of the invention involves determining expression by hybridization of mRNA, or an amplified or cloned version thereof, of a sample cell to a polynucleotide that is unique to a particular gene sequence. Preferred polynucleotides of this type contain at least about 20, at least about 22, at least about 24, at least about 26, at least about 28, at least about 30, or at least about 32 consecutive basepairs of a gene sequence that is not found in other gene sequences. The term “about” as used in the previous sentence refers to an increase or decrease of 1 from the stated numerical value. Even more preferred are polynucleotides of at least or about 50, at least or about 100, at least about or 150, at least or about 200, at least or about 250, at least or about 300, at least or about 350, or at least or about 400 basepairs of a gene sequence that is not found in other gene sequences. The term “about” as used in the preceding sentence refers to an increase or decrease of 10% from the stated numerical value. Such polynucleotides may also be referred to as polynucleotide probes that are capable of hybridizing to sequences of the genes, or unique portions thereof, described herein. Preferably, the sequences are those of mRNA encoded by the genes, the corresponding cDNA to such mRNAs, and/or amplified versions of such sequences. In preferred embodiments of the invention, the polynucleotide probes are immobilized on an array, other devices, or in individual spots that localize the probes.


Alternatively, and in another embodiment of the invention, gene expression may be determined by analysis of expressed protein in a cell sample of interest by use of one or more antibodies specific for one or more epitopes of individual gene products (proteins) in said cell sample. Such antibodies are preferably labeled to permit their easy detection after binding to the gene product.


The term “label” refers to a composition capable of producing a detectable signal indicative of the presence of the labeled molecule. Suitable labels include radioisotopes, nucleotide chromophores, enzymes, substrates, fluorescent molecules, chemiluminescent moieties, magnetic particles, bioluminescent moieties, and the like. As such, a label is any composition detectable by spectroscopic, photochemical, biochemical, immunochemical, electrical, optical or chemical means.


The term “support” refers to conventional supports such as beads, particles, dipsticks, fibers, filters, membranes and silane or silicate supports such as glass slides.


As used herein, a “breast tissue sample” or “breast cell sample” refers to a sample of breast tissue or fluid isolated from an individual suspected of being afflicted with, or at risk of developing, breast cancer. Such samples are primary isolates (in contrast to cultured cells) and may be collected by any non-invasive means, including, but not limited to, ductal lavage, fine needle aspiration, needle biopsy, the devices and methods described in U.S. Pat. No. 6,328,709, or any other suitable means recognized in the art. Alternatively, the “sample” may be collected by an invasive method, including, but not limited to, surgical biopsy.


“Expression” and “gene expression” include transcription and/or translation of nucleic acid material.


As used herein, the term “comprising” and its cognates are used in their inclusive sense; that is, equivalent to the term “including” and its corresponding cognates.


Conditions that “allow” an event to occur or conditions that are “suitable” for an event to occur, such as hybridization, strand extension, and the like, or “suitable” conditions are conditions that do not prevent such events from occurring. Thus, these conditions permit, enhance, facilitate, and/or are conducive to the event. Such conditions, known in the art and described herein, depend upon, for example, the nature of the nucleotide sequence, temperature, and buffer conditions. These conditions also depend on what event is desired, such as hybridization, cleavage, strand extension or transcription.


Sequence “mutation,” as used herein, refers to any sequence alteration in the sequence of a gene disclosed herein interest in comparison to a reference sequence. A sequence mutation includes single nucleotide changes, or alterations of more than one nucleotide in a sequence, due to mechanisms such as substitution, deletion or insertion. Single nucleotide polymorphism (SNP) is also a sequence mutation as used herein. Because the present invention is based on the relative level of gene expression, mutations in non-coding regions of genes as disclosed herein may also be assayed in the practice of the invention.


“Detection” includes any means of detecting, including direct and indirect detection of gene expression and changes therein. For example, “detectably less” products may be observed directly or indirectly, and the term indicates any reduction (including the absence of detectable signal). Similarly, “detectably more” product means any increase, whether observed directly or indirectly.


Unless defined otherwise all technical and scientific terms used herein have the same meaning as commonly understood to one of ordinary skill in the art to which this invention belongs.


Specific Embodiments

The present invention relates to the identification and use of gene expression patterns (or profiles or “signatures”) which discriminate between (or are correlated with) breast cancer survival and recurrence outcomes in a subject. Such patterns may be determined by the methods of the invention by use of a number of reference cell or tissue samples, such as those reviewed by a pathologist of ordinary skill in the pathology of breast cancer, which reflect breast cancer cells as opposed to normal or other non-cancerous cells. The outcomes experienced by the subjects from whom the samples may be correlated with expression data to identify patterns that correlate with the outcomes. Because the overall gene expression profile differs from person to person, cancer to cancer, and cancer cell to cancer cell, correlations between certain cells and genes expressed or underexpressed may be made as disclosed herein to identify genes that are capable of discriminating between breast cancer outcomes.


The present invention may be practiced with any number of the genes believed, or likely to be, differentially expressed with respect to breast cancer outcomes. The identification may be made by using expression profiles of various homogenous breast cancer cell populations, which were isolated by microdissection, such as, but not limited to, laser capture microdissection (LCM) of 100-1000 cells. The expression level of each gene of the expression profile may be correlated with a particular outcome. Alternatively, the expression levels of multiple genes may be clustered to identify correlations with particular outcomes.


Genes with significant correlations to breast cancer survival or recurrence outcomes may be used to generate models of gene expressions that would maximally discriminate between outcomes. Alternatively, genes with significant correlations may be used in combination with genes with lower correlations without significant loss of ability to discriminate between outcomes. Such models may be generated by any appropriate means recognized in the art, including, but not limited to, cluster analysis, supported vector machines, neural networks or other algorithm known in the art. The models are capable of predicting the classification of a unknown sample based upon the expression of the genes used for discrimination in the models. “Leave one out” cross-validation may be used to test the performance of various models and to help identify weights (genes) that are uninformative or detrimental to the predictive ability of the models. Cross-validation may also be used to identify genes that enhance the predictive ability of the models.


The gene(s) identified as correlated with particular breast cancer outcomes by the above models provide the ability to focus gene expression analysis to only those genes that contribute to the ability to identify a subject as likely to have a particular outcome relative to another. The expression of other genes in a breast cancer cell would be relatively unable to provide information concerning, and thus assist in the discrimination of, a breast cancer outcome.


As will be appreciated by those skilled in the art, the models are highly useful with even a small set of reference gene expression data and can become increasingly accurate with the inclusion of more reference data although the incremental increase in accuracy will likely diminish with each additional datum. The preparation of additional reference gene expression data using genes identified and disclosed herein for discriminating between different outcomes in breast cancer is routine and may be readily performed by the skilled artisan to permit the generation of models as described above to predict the status of an unknown sample based upon the expression levels of those genes.


To determine the (increased or decreased) expression levels of genes in the practice of the present invention, any method known in the art may be utilized. In one preferred embodiment of the invention, expression based on detection of RNA which hybridizes to the genes identified and disclosed herein is used. This is readily performed by any RNA detection or amplification+detection method known or recognized as equivalent in the art such as, but not limited to, reverse transcription-PCR, the methods disclosed in U.S. patent application Ser. No. 10/062,857 (filed on Oct. 25, 2001) as well as U.S. Provisional Patent Application 60/298,847 (filed Jun. 15, 2001) and 60/257,801 (filed Dec. 22, 2000), and methods to detect the presence, or absence, of RNA stabilizing or destabilizing sequences.


Alternatively, expression based on detection of DNA status may be used. Detection of the DNA of an identified gene as methylated or deleted may be used for genes that have decreased expression in correlation with a particular breast cancer outcome. This may be readily performed by PCR based methods known in the art, including, but not limited to, Q-PCR. Conversely, detection of the DNA of an identified gene as amplified may be used for genes that have increased expression in correlation with a particular breast cancer outcome. This may be readily performed by PCR based, fluorescent in situ hybridization (FISH) and chromosome in situ hybridization (CISH) methods known in the art.


Expression based on detection of a presence, increase, or decrease in protein levels or activity may also be used. Detection may be performed by any immunohistochemistry (IHC) based, blood based (especially for secreted proteins), antibody (including autoantibodies against the protein) based, exfoliate cell (from the cancer) based, mass spectroscopy based, and image (including used of labeled ligand) based method known in the art and recognized as appropriate for the detection of the protein. Antibody and image based methods are additionally useful for the localization of tumors after determination of cancer by use of cells obtained by a non-invasive procedure (such as ductal lavage or fine needle aspiration), where the source of the cancerous cells is not known. A labeled antibody or ligand may be used to localize the carcinoma(s) within a patient.


A preferred embodiment using a nucleic acid based assay to determine expression is by immobilization of one or more sequences of the genes identified herein on a solid support, including, but not limited to, a solid substrate as an array or to beads or bead based technology as known in the art. Alternatively, solution based expression assays known in the art may also be used. The immobilized gene(s) may be in the form of polynucleotides that are unique or otherwise specific to the gene(s) such that the polynucleotide would be capable of hybridizing to a DNA or RNA corresponding to the gene(s). These polynucleotides may be the full length of the gene(s) or be short sequences of the genes (up to one nucleotide shorter than the full length sequence known in the art by deletion from the 5′ or 3′ end of the sequence) that are optionally minimally interrupted (such as by mismatches or inserted non-complementary basepairs) such that hybridization with a DNA or RNA corresponding to the gene(s) is not affected. Preferably, the polynucleotides used are from the 3′ end of the gene. Polynucleotides containing mutations relative to the sequences of the disclosed genes may also be used so long as the presence of the mutations still allows hybridization to produce a detectable signal.


The immobilized gene(s) may be used to determine the state of nucleic acid samples prepared from sample breast cell(s) for which the outcome of the sample's subject (e.g. patient from whom the sample is obtained) is not known or for confirmation of an outcome that is already assigned to the sample's subject. Without limiting the invention, such a cell may be from a patient with breast cancer or alternatively suspected of being afflicted with, or at risk of developing, breast cancer. The immobilized polynucleotide(s) need only be sufficient to specifically hybridize to the corresponding nucleic acid molecules derived from the sample under suitable conditions. While even a single correlated gene sequence may to able to provide adequate accuracy in discriminating between two breast cancer outcomes, two or more, three or more, four or more, five or more, six or more, seven or more, eight or more, nine or more, ten or more, or eleven or more of the genes identified herein may be used as a subset capable of discriminating may be used in combination to increase the accuracy of the method. The invention specifically contemplates the selection of more than one, two or more, three or more, four or more, five or more, six or more, seven or more, eight or more, nine or more, ten or more, or eleven or more of the genes disclosed in the tables and figures herein for use as a subset in the identification of breast cancer survival outcome.


Of course 15 or more, 20 or more, 30 or more, 40 or more, 50 or more, 60 or more, 70 or more, 80 or more, 90 or more, 100 or more, 110 or more, 120 or more, 130 or more, 140 or more, or all the genes provided in Tables 2, 3, and/or 4 below may be used. “CloneID” as used in the context of Tables 2, 3, and 4 as well as the present invention refers to the IMAGE Consortium clone ID number of each gene, the sequences of which are hereby incorporated by reference in their entireties as they are available from the Consortium at http://image.llnl.gov/ as accessed on the filing date of the present application. Also provided in the tables are GenBank accession numbers which are comprised of letters, numbers and optionally underscores. P value refers to values assigned as described in the Examples below. The indications of “E-xx” where “xx” is a two digit number refers to alternative notation for exponential figures where “E-xx” is “10-xx”. Thus in combination with the numbers to the left of “E-xx”, the value being represented is the numbers to the left times 10-xx. Description provides a brief identifier of what the gene encodes.


Genes with a correlation identified by a p value below or about 0.02, below or about 0.01, below or about 0.005, or below or about 0.001 are preferred for use in the practice of the invention. The present invention includes the use of genes that identify different ER positive subtypes and breast cancer recurrence and invasive tumor grade to permit simultaneous identification of breast cancer survival outcome of a patient based upon assaying a breast cancer sample from said patient.


In embodiments where only one or a few genes are to be analyzed, the nucleic acid derived from the sample breast cancer cell(s) may be preferentially amplified by use of appropriate primers such that only the genes to be analyzed are amplified to reduce contaminating background signals from other genes expressed in the breast cell. Alternatively, and where multiple genes are to be analyzed or where very few cells (or one cell) is used, the nucleic acid from the sample may be globally amplified before hybridization to the immobilized polynucleotides. Of course RNA, or the cDNA counterpart thereof may be directly labeled and used, without amplification, by methods known in the art.


The above assay embodiments may be used in a number of different ways to identify or detect the invasive breast cancer grade, if any, of a breast cancer cell sample from a patient. In many cases, this would reflect a secondary screen for the patient, who may have already undergone mammography or physical exam as a primary screen. If positive, the subsequent needle biopsy, ductal lavage, fine needle aspiration, or other analogous methods may provide the sample for use in the above assay embodiments. The present invention may be used in combination with non-invasive protocols, such as ductal lavage or fine needle aspiration, to prepare a breast cell sample.


The present invention provides a more objective set of criteria, in the form of gene expression profiles of a discrete set of genes, to discriminate (or delineate) between breast cancer outcomes. In particularly preferred embodiments of the invention, the assays are used to discriminate between good and poor outcomes within 5, or about 5, years after surgical intervention to remove breast cancer tumors or within about 95 months after surgical intervention to remove breast cancer tumors. Comparisons that discriminate between outcomes after about 10, about 20, about 30, about 40, about 50, about 60, about 70, about 80, about 90, or about 100 months may also be performed.


While good and poor survival outcomes may be defined relatively in comparison to each other, a “good” outcome may be viewed as a better than 50% survival rate after about 60 months post surgical intervention to remove breast cancer tumor(s). A “good” outcome may also be a better than about 60%, about 70%, about 80% or about 90% survival rate after about 60 months post surgical intervention. A “poor” outcome may be viewed as a 50% or less survival rate after about 60 months post surgical intervention to remove breast cancer tumor(s). A “poor” outcome may also be about a 70% or less survival rate after about 40 months, or about a 80% or less survival rate after about 20 months, post surgical intervention.


In one embodiment of the invention, the isolation and analysis of a breast cancer cell sample may be performed as follows:


(1) Ductal lavage or other non-invasive procedure is performed on a patient to obtain a sample.


(2) Sample is prepared and coated onto a microscope slide. Note that ductal lavage results in clusters of cells that are cytologically examined as stated above.


(3) Pathologist or image analysis software scans the sample for the presence of non-normal and/or atypical breast cancer cells.


(4) If such cells are observed, those cells are harvested (e.g. by microdissection such as LCM).


(5) RNA is extracted from the harvested cells.


(6) RNA is purified, amplified, and labeled.


(7) Labeled nucleic acid is contacted with a microarray containing polynucleotides of the genes identified herein as correlated to discriminations between breast cancer outcomes under suitable hybridization conditions, then processed and scanned to obtain a pattern of intensities of each spot (relative to a control for general gene expression in cells) which determine the level of expression of the gene(s) in the cells.


(8) The pattern of intensities is analyzed by comparison to the expression patterns of the genes in known samples of breast cancer cells correlated with outcomes (relative to the same control).


A specific example of the above method would be performing ductal lavage following a primary screen, observing and collecting non-normal and/or atypical cells for analysis. The comparison to known expression patterns, such as that made possible by a model generated by an algorithm (such as, but not limited to nearest neighbor type analysis, SVM, or neural networks) with reference gene expression data for the different breast cancer survival outcomes, identifies the cells as being correlated with subjects with good or poor outcomes. Another example would be taking a breast tumor removed from a subject after surgical intervention, isolation and preparation of breast cancer cells from the tumor for determination/identification of atypical, non-normal, or cancer cells, and isolation of said cells followed by steps 5 through 8 above.


Alternatively, the sample may permit the collection of both normal as well as cancer cells for analysis. The gene expression patterns for each of these two samples will be compared to each other as well as the model and the normal versus individual comparisons therein based upon the reference data set. This approach can be significantly more powerful that the cancer cells only approach because it utilizes significantly more information from the normal cells and the differences between normal and cancer cells (in both the sample and reference data sets) to determine the breast cancer outcome of the patient based on gene expression in the cancer cells from the sample.


With use of the present invention, skilled physicians may prescribe treatments based on prognosis determined via non-invasive samples that they would have prescribed for a patient which had previously received a diagnosis via a solid tissue biopsy.


The above discussion is also applicable where a palpable lesion is detected followed by fine needle aspiration or needle biopsy of cells from the breast. The cells are plated and reviewed by a pathologist or automated imaging system which selects cells for analysis as described above.


The present invention may also be used, however, with solid tissue biopsies. For example, a solid biopsy may be collected and prepared for visualization followed by determination of expression of one or more genes identified herein to determine the breast cancer outcome. One preferred means is by use of in situ hybridization with polynucleotide or protein identifying probe(s) for assaying expression of said gene(s).


In an alternative method, the solid tissue biopsy may be used to extract molecules followed by analysis for expression of one or more gene(s). This provides the possibility of leaving out the need for visualization and collection of only cancer cells or cells suspected of being cancerous. This method may of course be modified such that only cells that have been positively selected are collected and used to extract molecules for analysis. This would require visualization and selection as a prerequisite to gene expression analysis.


In a further modification of the above, both normal cells and cancer cells are collected and used to extract molecules for analysis of gene expression. The approach, benefits and results are as described above using non-invasive sampling.


The genes identified herein may be used to generate a model capable of predicting the breast cancer survival and recurrence outcomes of an unknown breast cell sample based on the expression of the identified genes in the sample. Such a model may be generated by any of the algorithms described herein or otherwise known in the art as well as those recognized as equivalent in the art using gene(s) (and subsets thereof) disclosed herein for the identification of breast cancer outcomes. The model provides a means for comparing expression profiles of gene(s) of the subset from the sample against the profiles of reference data used to build the model. The model can compare the sample profile against each of the reference profiles or against a model defining delineations made based upon the reference profiles. Additionally, relative values from the sample profile may be used in comparison with the model or reference profiles.


In a preferred embodiment of the invention, breast cell samples identified as normal and cancerous from the same subject may be analyzed for their expression profiles of the genes used to generate the model. This provides an advantageous means of identifying survival and recurrence outcomes based on relative differences from the expression profile of the normal sample. These differences can then be used in comparison to differences between normal and individual cancerous reference data which was also used to generate the model.


The detection of gene expression from the samples may be by use of a single microarray able to assay gene expression from some or all genes disclosed herein for convenience and accuracy.


Other uses of the present invention include providing the ability to identify breast cancer cell samples as correlated with particular breast cancer survival or recurrence outcomes for further research or study. This provides a particular advantage in many contexts requiring the identification of cells based on objective genetic or molecular criteria.


The materials for use in the methods of the present invention are ideally suited for preparation of kits produced in accordance with well known procedures. The invention thus provides kits comprising agents for the detection of expression of the disclosed genes for identifying breast cancer outcomes. Such kits optionally comprising the agent with an identifying description or label or instructions relating to their use in the methods of the present invention, is provided. Such a kit may comprise containers, each with one or more of the various reagents (typically in concentrated form) utilized in the methods, including, for example, pre-fabricated microarrays, buffers, the appropriate nucleotide triphosphates (e.g., dATP, dCTP, dGTP and dTTP; or rATP, rCTP, rGTP and UTP), reverse transcriptase, DNA polymerase, RNA polymerase, and one or more primer complexes of the present invention (e.g., appropriate length poly(T) or random primers linked to a promoter reactive with the RNA polymerase). A set of instructions will also typically be included.


The methods provided by the present invention may also be automated in whole or in part. All aspects of the present invention may also be practiced such that they consist essentially of a subset of the disclosed genes to the exclusion of material irrelevant to the identification of breast cancer survival outcomes via a cell containing sample.


Having now generally described the invention, the same will be more readily understood through reference to the following examples which are provided by way of illustration, and are not intended to be limiting of the present invention, unless specified.


EXAMPLES
Example I
General

Clinical specimen collection and clinicopathological parameters. Laser capture microdissected invasive cancer cells from a total of 124 breast cancer biopsies were used to discover two sets of genes, the expression levels of which correlate with clinical breast cancer outcomes. These genes could thus be used either individually or in combination as prognostic factors for breast cancer management. The characteristics of the 124 patient profiles in the study are shown in Table 1.


Relative expression levels of ˜22000 genes were measured from the invasive cancer cells for each of the 124 patients. Genes varying by at least 3-fold from the median expression level across the 124 patients in at least 10 patients were selected, resulting in 7090 genes.


In particular, 4 genes (DEEPEST, RACGAP1, ZNF145, MS4A7) were shown to be strong prognostic factors individually for predicting tumor recurrence after surgery and adjuvant therapies.

TABLE 1GroupN%Age<=453024.245-552721.8>556754ERpositive6653.2negative5846.7Nodepositive5855.2negative4744.8Not avail.1915Grade189.622934.934655.4Not avail.4133.1


Example II
Identification of ER Positive Subtypes with Different Survival Outcomes

Hierarchical clustering, based on the 7090 genes described in Example 1, of the resulting gene expression matrix (7090×124) revealed a cluster of 67-genes (the Ki67 set) the expressions of which differentiates estrogen receptor positive patients into two subgroups with distinct clinical outcomes based on overall survival over time.


As shown in FIG. 1, left panel, a Kaplan-Meier curve on the left compares the disease-free survival of patients based on ER status, which shows slightly better survival for ER positive patients but with an insignificant p value (log-rank test). In contrast, and as shown in the right panel, when the ER positive patients are subdivided into two subgroups (A and B) based on the expression levels of the Ki67 signature genes, which are all expressed at levels above the median to define subgroup A and below the median to define subgroup B.


The three-group (ER+, subgroup A; ER+, subgroup B; and ER−) comparison shows significant differences in survival such that subgroup B subjects had significantly better survival outcomes than those of subgroup A. The ER− curve remains unchanged. This indicates that the Ki67 signature, and individual or groups of genes therein, can be used to subdivide ER positive patients into two clinically distinct subgroups based upon survival outcomes.


The identities of the genes in the Ki67 signature are shown in Table 2.

TABLE 2Genes, the expressions of which define two ER+ subgroupsCloneIDGeneDescription2967734BC007491EXO1 | exonuclease 1NM_000057BLM | Bloom syndrome2849551AW512559CDC25C | cell division cycle 25C3634656BC010044CDC20 | CDC20 cell division cycle 20 homolog (S. cerevisiae)2961114BC008718BIRC5 | baculoviral IAP repeat-containing 5 (survivin)NM_012112C20orf1 | chromosome 20 open reading frame 1AF399910DEEPEST | mitotic spindle coiled-coil related proteinNM_032997ZWINT | ZW10 interactorAF331796HCAP-G | chromosome condensation protein G2175265AI524385ANLN | “anillin, actin binding protein (scraps homolog, Drosophila)”3873367BC010658KIAA0008 | KIAA0008 gene product1338423AA810180FLJ10517 | hypothetical protein FLJ10517AL136794RACGAP1 | Rac GTPase activating protein 1AF334184FKSG42 | FKSG424420248BC017705KNSL5 | kinesin-like 5 (mitotic kinesin-like protein 1)AB035898KNSL7 | kinesin-like 71240937AA714213ESTs, Highly similar to T47163 hypothetical protein DKFZp762E1312.1 [H. sapiens]2820741BC001940DKFZp762E1312 | hypothetical protein DKFZp762E1312AF017790HEC | “highly expressed in cancer, rich in leucine heptad repeats”4048625BC013919TYMS | thymidylate synthetase1699365AI049877KIAA0186 | KIAA0186 gene product3139011BC001459RAD51 | “RAD51 homolog (RecA homolog, E. coli) (S. cerevisiae)”AF053306BUB1B | BUB1 budding uninhibited by benzimidazoles 1 homolog beta (yeast)3908972BC015050OIP5 | Opa-interacting protein 51911633AI268609ESPL1 | extra spindle poles like 1 (S. cerevisiae)NM_002417MKI67 | antigen identified by monoclonal antibody Ki-672988318BC013966FLJ10156 | hypothetical protein2964488BC013300STK12 | serine/threonine kinase 12NM_016343CENPF | “centromere protein F (350/400 kD, mitosin)”AF161499HSPC150 | HSPC150 protein similar to ubiquitin-conjugating enzymeAL136840MCM10 | MCM10 minichromosome maintenance deficient 10 (S. cerevisiae)3028566BC008947FLJ10540 | hypothetical protein FLJ105401986322AI273114ESTs, Weakly similar to I78885 serine/threonine-specific protein kinase [H. sapiens]3909951BC015706Homo sapiens, Similar to RIKEN cDNA 2810433K01 gene, clone MGC: 10200 IMAGE: 3909951,mRNA, complete cdsAF095289PTTG3 | pituitary tumor-transforming 32262695AI811894PTTG2 | pituitary tumor-transforming 21186167AA648922CDC25A | cell division cycle 25A3138951BC002551MGC2577 | hypothetical protein MGC2577U74612FOXM1 | forkhead box M13347875BC000703FLJ10468 | hypothetical protein FLJ104683347571BC008764KNSL6 | kinesin-like 6 (mitotic centromere-associated kinesin)2822981BC000404TRIP13 | thyroid hormone receptor interactor 131678629AI082049ESTsNM_003504CDC45L | CDC45 cell division cycle 45-like (S. cerevisiae)3345575BC007656UBE2C | ubiquitin-conjugating enzyme E2CD84212STK6 | serine/threonine kinase 6AF011468STK15 | serine/threonine kinase 153938081BC011000MGC16386 | similar to RIKEN cDNA2610036L13AF277375KIF4A | kinesin family member 4A3461992BC000881CENPA | centromere protein A (17 kD)AF108138PIF1 | DMA helicase homolog PIF1AF155827HELLS | “helicase, lymphoid-specific”NM_018492TOPK | T-LAK cell-originated protein kinase1686560AI088843ESTs1241465AA715810ESTs, Weakly similar to YK61_YEAST HYPOTHETICAL 39.6 KDA PROTEIN IN MTD1-NUP133INTERGENIC REGION [S. cerevisiae]2823731BC001068C20orf129 | chromosome 20 open reading frame 129AK026964FLJ23311 | hypothetical protein FLJ233113996265BC005389LOC51053 | geminin3901250BC010858EZH2 | enhancer of zeste homolog 2 (Drosophila)4547136BC014039KIAA0175 | likely ortholog of maternal embryonic leucine zipper kinase4091997BC017575CHEK1 | CHK1 checkpoint homolog (S. pombe)669114AA232651SUV39H2 | suppressor of variegation 3-9 (Drosophila) homolog 2; hypothetical protein FLJ23414NM_002497NEK2 | NIMA (never in mitosis gene a)-related kinase 24107592BC016330PIR51 | RAD51-interacting proteinAF025840POLE2 | “polymerase (DMA directed), epsilon 2”3510656BC007633EIF2C2 | “eukaryotic translation initiation factor 2C, 2”AL050151Homo sapiens mRNA; cDNA DKFZp586J0720 (from clone DKFZp586J0720)


Example III
Molecular Signature that Correlates with the Recurrence of Breast Cancer

A molecular signature that correlates with recurrence of breast cancer after removal of cancer by surgery was identified as follows. Each of the 7090 genes from Example 1 was used to fit a univariate Cox proportional hazard regression model using the survival information available for the patients in the study. A total of 143 genes with significant p values (p<0.01) in these univariate models were selected. Hierarchical clustering of patient samples by the 143 recurrence-associated genes identified them as having expression levels that correlated with the absence or presence of breast cancer recurrence.


These 143 genes are shown in Table 3. The sign of the coefficient values in Table 3 correspond to whether a gene is positively or negatively correlated with breast cancer recurrence. A positive coefficient means that the gene is positively correlated (overexpressed) in patients with a poor (shorter) survival outcome due to recurrence and negatively correlated (underexpressed) in patients with a good or better (longer) survival outcome due to the relative absence of recurrence. A negative coefficient means that the gene is positively correlated (overexpressed) in patients with a good or better (longer) survival outcome (due to the relative absence of cancer recurrence) and negatively correlated (underexpressed) in patients with a poor (shorter) survival outcome (due to cancer recurrence).


To validate this gene set, 22 of the top 27 genes from Table 3 (with the smallest p values) were mapped onto the microarray used by van't Veer et al. (Supra) via the Unigene database. The top 27 genes are provided in Table 4 while the mapping of genes are shown in Table 5 (showing identities of the genes via their GenBank ID, van't Veer et al. reference, and Unigene ID numbers). Thirteen of the 22 genes were filtered out due to low variance across the sample set, reducing the number of genes for cluster analysis to 9. The 27 gene set was used with the data from the patients of Example 1 to classify them as being in either the good prognosis or the poor prognosis group by hierarchical clustering based on disease-free survival. The results are shown in FIG. 2, left panel (Kaplan-Meier curves of patients stratified by the top 27 recurrence-associated genes).


The 9 genes not filtered out from the van't Veer et al. data were used to with the patient data therein to classify them as being in either the good prognosis or the poor prognosis group by hierarchical clustering based on disease-free survival. The results are shown in FIG. 2, right panel (Kaplan-Meier curves of patients stratified by 9 of the top 27 recurrence-associated genes).


Like FIG. 1, the horizontal axis of FIG. 2 is in time (months or years) and the vertical axis is in survival probability (where 1.0 is survival of 100% of the subjects and 0.5 is survival of 50% of the subjects). As shown in FIG. 2, differences in disease-free survival between the two groups in both datasets were highly significant.

TABLE 3Genes, the expressions of which correlate with breast cancer recurrenceClone IDgenepcoefdesc1184567AA6487777.58E−06−2.3582882MS4A7 | membrane-spanning 4-domains, subfamily A, member 72961112BC0058501.05E−04−1.845548CBFA2T1 | core-binding factor, runt domain, alpha subunit 2;translocated to, 1; cyclin D-related3565773BF4328131.65E−03−1.2898777KLRB1 | killer cell lectin-like receptor subfamily B, member 11352935AA8301317.60E−03−1.2502516ZNF80 | zinc finger protein 80 (pT17)3915193BC0170221.80E−03−0.9385143Homo sapiens, clone MGC: 8979 IMAGE: 3915193, mRNA, complete cds2630949AW1502672.87E−04−0.9172496C21orf9 | chromosome 21 open reading frame 92714519AW1379911.22E−03−0.9027559RELB | v-rel reticuloendotheliosis viral oncogene homolog B,nuclear factor of kappa light polypeptide gene enhancer in B-cells 3 (avian)2365891AI7417857.56E−03−0.8965429SLIT3 | slit homolog 3 (Drosophila)NM_0060068.59E−04−0.8880113ZNF145 | zinc finger protein 145 (Kruppel-like, expressed in promyelocytic leukemia)3645909BF4366563.08E−03−0.8506381MFAP4 | microfibrillar-associated protein 42349778AI8061093.89E−03−0.8414599KIAA1580 | KIAA1580 protein3504259BC0007237.95E−03−0.8413126CRAT | camitine acetyltransferase4342203BC0185382.92E−04−0.8358585ALOX5AP | arachidonate 5-lipoxygenase-activating proteinAL1220524.91E−04−0.8340079KIAA0793 | KIAA0793 gene productAK0250916.06E−03−0.8263796FLJ21438 | hypothetical protein FLJ214382612878AW1308882.14E−03−0.8099805PTK2B | protein tyrosine kinase 2 betaAF2441298.92E−03−0.807578LY9 | lymphocyte antigen 9AK0271208.04E−03−0.8071951FLJ23467 | hypothetical protein FLJ23467AF3674731.85E−04−0.8035463NYD-SP21 | testes development-related NYD-SP214214447BC0090324.52E−03−0.7976367PR48 | protein phosphatase 2A 48 kDa regulatory subunitAB0458325.30E−03−0.767161P53AIP1 | p53-regulated apoptosis-inducing protein 1NM_0005982.64E−03−0.7567474IGFBP3 | insulin-like growth factor binding protein 3AI9520555.06E−03−0.7474092ESTsNM_0037344.42E−03−0.7298537AOC3 | amine oxidase, copper containing 3 (vascular adhesion protein 1)4291158BC0083926.23E−03−0.7260119UCP3 | uncoupling protein 3 (mitochondrial, proton earner)3840457BC0129902.07E−03−0.7130147Homo sapiens, clone IMAGE: 3840457, mRNAAB0378862.20E−03−0.6839819NESH | NESH protein3622951BC0043006.00E−03−0.6820265VILL | villin-likeNM_0153856.90E−03−0.6555434SH3D5 | SH3-domain protein 5 (ponsin)289749N592844.93E−03−0.6497849ESTs3677098BC0048642.47E−03−0.6450665PPP3CC | protein phosphatase 3 (formerly 2B), catalytic subunit,gamma isoform (calcineurin A gamma)2254324AI6209651.08E−03−0.6415661ESTs1848897AI2479018.54E−03−0.6392025ESTs, Weakly similar to S23650 retrovirus-related hypothetical protein II [H. sapiens]4699374BC0178391.24E−03−0.6348455CASP4 | caspase 4, apoptosis-related cysteine proteaseU908783.45E−03−0.6303001PDLIM1 | PDZ and LIM domain 1 (elfin)2729801AW2938493.26E−03−0.6188886ESTs, Moderately similar to I54374 gene NF2 protein [H. sapiens]AL1376949.86E−03−0.6154861FLJ11286 | hypothetical protein FLJ112861884362AI2159028.92E−04−0.6139098ESTs, Highly similar to T50835 hypothetical protein [H. sapiens]3543310BC0016097.61E−03−0.6113849WBSCR5 | Williams-Beuren syndrome chromosome region 51869453AI2646445.93E−03−0.6105146KIAA0775 | KIAA0775 gene product3010091BC0061071.60E−03−0.6088258ARHGAP9 | Rho GTPase activating protein 9NM_0024055.53E−03−0.6041611MFNG | manic fringe homolog (Drosophila)AK0263431.86E−04−0.6018579FLJ22690 | hypothetical protein FLJ226902227051AI5831092.88E−03−0.5977938STAT5A | signal transducer and activator of transcription 5A1144648AA6135606.79E−03−0.5968032ALOX5 | arachidonate 5-lipoxygenase206683H595598.44E−03−0.5967191ESTsAK0216748.80E−03−0.5841158Homo sapiens cDNA FLJ11612 fis, clone HEMBA10040112364492AI7410864.47E−03−0.5804208ESTsBF7250078.40E−05−0.5750824ADRA2A | adrenergic, alpha-2A-, receptorAL0503913.94E−03−0.5724837Homo sapiens mRNA; cDNA DKFZp586A181 (from clone DKFZp586A181); partial cdsAF3674707.17E−03−0.5715338NYD-SP18 | testes development-related NYD-SP18293605AK0267475.13E−03−0.5606859LOC54103 | hypothetical proteinAK0257326.04E−03−0.5567897ASAH | N-acylsphingosine amidohydrolase (acid ceramidase)1645681AI0268384.77E−03−0.5536825ESTs, Weakly similar to NUCL_HUMAN NUCLEOLIN [H. sapiens]1670862AI0812358.40E−03−0.5523028CD53 | CD53 antigen3703127BF4336865.19E−03−0.5433968Homo sapiens cDNA FLJ32651 fis, clone SYNOV20015811837189AF3397814.94E−03−0.5432118GPR18 | G protein-coupled receptor 184263201BG2366451.63E−03−0.5426042ESTs4309471BC0099563.25E−03−0.5390891HLA-DPA1 | major histocompatibility complex, class II, DP alpha 131047R424632.11E−03−0.5357526ENTPD1 | ectonucleoside triphosphate diphosphohydrolase 1L027856.69E−03−0.5355299SLC26A3 | solute carrier family 26, member 3NM_0013372.41E−03−0.5279489CX3CR1 | chemokine (C-X3-C) receptor 1BC0167582.54E−03−0.5181311HCLS1 | hematopoietic cell-specific Lyn substrate 12214761AI5654895.26E−03−0.5077388PDE4A | phosphodiesterase 4A, cAMP-specific (phosphodiesteraseE2 dunce homolog, Drosophila)2483676BI4920734.82E−03−0.5068155ITM2A | integral membrane protein 2A128753R168389.27E−03−0.5059294ESTs4548935BC0141176.62E−03−0.5027082TBXAS1 | thromboxane A synthase 1 (platelet, cytochrome P450, subfamily V)1734062AI1916206.36E−03−0.5006683CDO1 | cysteine dioxygenase, type INM_0038208.71E−03−0.5005279TNFRSF14 | tumor necrosis factor receptor superfamily, member 14(herpesvirus entry mediator)AF3054283.48E−03−0.4933213APOL1 | apolipoprotein L, 13163446BC0087346.56E−03−0.4928143FCGRT | Fc fragment of IgG, receptor, transporter, alpha1964662AJ4205854.18E−04−0.4903945Homo sapiens mRNA full length insert cDNA clone EUROIMAGE 19646621351991AA8073465.64E−03−0.4893204Homo sapiens cDNA FLJ14296 fis, clone PLACE1008455NM_0052118.18E−03−0.485557CSF1R | colony stimulating factor 1 receptor, formerly McDonoughfeline sarcoma viral (v-fms) oncogene homolog2161081AI5802719.09E−03−0.4836746AFP | alpha-fetoprotein4862198BC0144567.46E−03−0.4809022CHRNA6 | cholinergic receptor, nicotinic, alpha polypeptide 62728733AW2951708.18E−03−0.4691647ESTs2163996AI4794615.00E−03−0.4666083CSR1 | CSR1 protein3086130BF5092354.52E−03−0.4567332KIAA1658 | KIAA1658 protein2364383AI7406712.10E−03−0.4565665Homo sapiens cDNA FLJ32430 fis, clone SKMUS2001129, weakly similarto NAD-DEPENDENT METHANOL DEHYDROGENASE (EC 1.1.1.244)1272059AA7432836.97E−03−0.4459972GZMK | granzyme K (serine protease, granzyme 3; tryptase II)3902651BC0168419.58E−03−0.440234RAB34 | RAB34, member RAS oncogene familyAB0587088.44E−03−0.4205688KIAA1805 | KIAA1805 protein40879R560533.88E−03−0.4205535ME3 | malic enzyme 3, NADP(+)-dependent, mitochondrial2222621AI5726057.12E−03−0.3955994HLA-DRA | major histocompatibility complex, class II, DR alpha2423726AI8603607.08E−03−0.3889526ESTs2586524AW0808312.88E−03−0.3552714SEC14L2 | SEC14-like 2 (S. cerevisiae)1056761AA5741748.03E−03−0.3505373CYP2A7 | cytochrome P450, subfamily IIA (phenobarbital-inducible), polypeptide 7NM_0333808.90E−03−0.3281259COL4A5 | collagen, type IV, alpha 5 (Alport syndrome)2148123AI4678468.94E−03−0.3273129IAN4L1 | immune associated nucleotide 4 like 1 (mouse)3026606BE0463256.30E−030.3235462IGFBP5 | insulin-like growth factor binding protein 53939513BC0138827.56E−030.3511008EYA2 | eyes absent homolog 2 (Drosophila)AK0573392.22E−030.359611LOC81569 | actin like proteinAF0071949.59E−030.3752148MUC3A | mucin 3A, intestinalAF2883952.13E−030.3941496C1orf15 | chromosome 1 open reading frame 153846346BC0170337.82E−030.4089227SQLE | squalene epoxidase3463613BC0036845.11E−030.42608CXADR | coxsackie virus and adenovirus receptor2190016AI5382268.41E−040.4327835GNG4 | guanine nucleotide binding protein 42138200AI5222158.62E−030.4352044KIAA1804 | KIAA1804 protein3087716BF5109797.81E−030.4409617DHDH | dihydrodiol dehydrogenase (dimeric)5677199BM1293937.66E−030.4526335GDF1 | growth differentiation factor 12144913AI4526344.23E−030.4561948GPR64 | G protein-coupled receptor 64M955853.62E−030.5164108HLF | hepatic leukemia factor3932186BC0053456.63E−030.5470575GTF2H2 | general transcription factor IIH, polypeptide 2 (44 kD subunit)2968940AW6138547.48E−030.5487125ESTs, Moderately similar to S02826 nonhistone chromosomal protein HMG-1 [H. sapiens]AF0177905.02E−030.5639359HEC | highly expressed in cancer, rich in leucine heptad repeatsAY0497376.54E−030.5703928NPM3 | nucleophosmin/nucleoplasmin, 3U877918.91E−030.5755378HBS1L | HBS1-like (S. cerevisiae)3915484BC0170539.97E−030.5772075ACOX3 | acyl-Coenzyme A oxidase 3, pristanoylAF1007519.93E−030.5901307LOC51661 | FK506-binding proteinAF2066739.98E−030.6026123BRF2 | BRF2, subunit of RNA polymerase III transcription initiation factor, BRF1-like1251833AA7312076.11E−030.6062738FLJ10858 | hypothetical protein FLJ108582507739AI9613692.77E−030.6190266INSIG1 | insulin induced gene 13504930BC0051415.35E−030.621928GALK2 | galactokinase 2AL1365706.98E−030.6225449LHX6 | LIM homeobox protein 62735278AW4507317.74E−030.6254488FLJ14642 | hypothetical protein FLJ14642AK0258207.09E−030.6282847FLJ22167 | hypothetical protein FLJ221672975886AW6291769.47E−040.6387345ESTs, Weakly similar to I38022 hypothetical protein [H. sapiens]AF1166704.34E−030.6509022NP | nucleoside phosphorylase1630968AI0186053.48E−030.6689059ESTs3996449BC0151077.48E−030.6910834FLJ13433 | hypothetical protein FLJ13433AB0491139.76E−030.7010916DUT | dUTP pyrophosphataseAK0255438.56E−030.707967KIAA1345 | KIAA1345 proteinAF0533066.31E−030.7227021BUB1B | BUB1 budding uninhibited by benzimidazoles 1 homolog beta (yeast)3010092BC0089549.78E−030.7272841SLC29A1 | solute carrier family 29 (nucleoside transporters), member 1NM_0016856.48E−030.7348505ATP5J | ATP synthase, H+ transporting, mitochondrial F0 complex, subunit F64838878BC0167515.11E−030.7473835PCDHB3 | protocadherin beta 32175265AI5243852.93E−030.7701449ANLN | anillin, actin binding protein (scraps homolog, Drosophila)3010727BE2060761.25E−040.7737357ALK | anaplastic lymphoma kinase (Ki-1)X625344.09E−030.7848284HMG2 | high-mobility group (nonhistone chromosomal) protein 2NM_0186693.59E−030.8006609WDR4 | WD repeat domain 4AB0358983.72E−030.8063743KNSL7 | kinesin-like 7NM_0067348.82E−030.8275957HIVEP2 | human immunodeficiency virus type I enhancer binding protein 2AF3999101.43E−030.8370408DEEPEST | mitotic spindle coiled-coil related proteinAF3317962.69E−030.8636081HCAP-G | chromosome condensation protein GAF0735186.93E−030.8902775SERF1A | small EDRK-rich factor 1A (telomeric)1870184AI2458071.20E−030.8964932MGC14798 | similar to RIKEN cDNA 5730421E18 gene4509200BC0129191.19E−030.9209921KLF7 | Kruppel-like factor 7 (ubiquitous)3926227BC0098551.51E−030.9417395FLJ14909 | hypothetical protein FLJ149091337864AA8113762.04E−031.0013265FLJ10545 | hypothetical protein FLJ105454109322BC0167822.22E−031.0106451KIAA0101 | KIAA0101 gene productAF3341841.66E−041.0332632FKSG42 | FKSG42AL1367941.98E−041.1631056RACGAP1 | Rac GTPase activating protein 1









TABLE 4










Top 27 genes










gene
p
coef
desc













BE206076
7.27446E−05
1.2270555
ALK | anaplastic lymphoma kinase (Ki-1)


AF040628
0.000116979
0.9344455
ED1 | ectodermal dysplasia 1, anhidrotic


BF725007
0.000331932
−0.5801744
ADRA2A | adrenergic, alpha-2A-, receptor


AF367473
0.00068377
−0.8560356
NYD-SP21 | testes development-related NYD-SP21


AI245807
0.000789859
1.0889127
MGC14798 | similar to RIKEN cDNA 5730421E18 gene


AI215902
0.00086535
−0.6485096
ESTs, Highly similar to T50835 hypothetical protein [H. sapiens]


AF334184
0.000883849
1.0194484
FKSG42 | FKSG42


AW137991
0.000969904
−1.1497247
RELB | v-rel reticuloendotheliosis viral oncogene homolog B,





nuclear factor of kappa light polypeptide gene enhancer in B-cells 3 (avian)


AA648777
0.001013976
−1.6577293
MS4A7 | membrane-spanning 4-domains, subfamily A, member 7


BC017053
0.001247656
0.8780177
ACOX3 | acyl-Coenzyme A oxidase 3, pristanoyl


AL136570
0.00133504
0.822904
LHX6 | LIM homeobox protein 6


AL136794
0.001396755
1.1506902
RACGAP1 | Rac GTPase activating protein 1


AF331796
0.001490968
1.2361318
HCAP-G | chromosome condensation protein G


BC005850
0.001567107
−1.5082099
CBFA2T1 | core-binding factor, runt domain, alpha subunit 2; translocated to, 1;





cyclin D-related


AF399910
0.00160607
0.8764045
DEEPEST | mitotic spindle coiled-coil related protein


AK057339
0.001660726
0.4461839
LOC81569 | actin like protein


NM_003265
0.001664793
−0.7178243
TLR3 | toll-like receptor 3


AK026343
0.001697437
−0.5688145
FLJ22690 | hypothetical protein FLJ22690


BC018538
0.001728603
−0.7262933
ALOX5AP | arachidonate 5-lipoxygenase-activating protein


AI806109
0.001789313
−1.0762434
KIAA1580 | KIAA1580 protein


AL122052
0.001810196
−0.9558644
KIAA0793 | KIAA0793 gene product


BC012919
0.002008228
0.9606709
KLF7 | Kruppel-like factor 7 (ubiquitous)


BC008392
0.002082612
−1.0314954
UCP3 | uncoupling protein 3 (mitochondrial, proton carrier)


BF432813
0.002110755
−1.2468785
KLRB1 | killer cell lectin-like receptor subfamily B, member 1


AI741086
0.002281532
−0.675948
ESTs


AK022729
0.002327927
−0.9643334
KIAA1681 | KIAA1681 protein


NM_006006
0.002390747
−1.0628839
ZNF145 | zinc finger protein 145 (Kruppel-like, expressed in promyelocytic leukemia)




















TABLE 5











GenBank
van't Veer et al.
UniGene









AA648777
AF201951
Hs.11090



AF367473
AL137391
Hs.28514



AK022729
Contig30485_RC
Hs.42656



AI741086
Contig39054_RC
Hs.115122



AK022729
Contig47136_RC
Hs.42656



AI215902
Contig52342_RC
Hs.88845



BF725007
Contig53357_RC
Hs.249159



BF725007
NM_000681
Hs.249159



AF040628
NM_001399
Hs.105407



BC018538
NM_001629
HS.100194



BF432813
NM_002258
Hs.169824



NM_003265
NM_003265
Hs.29499



BC008392
NM_003356
Hs.101337



BC017053
NM_003501
HS.12773



BC012919
NM_003709
Hs.21599



BE206076
NM_004304
Hs.278572



BC005850
NM_004349
Hs.31551



NM_006006
NM_006006
Hs.37096



AF399910
NM_006461
Hs.16244



AW137991
NM_006509
Hs.858



AL136794
NM_013277
Hs.23900



AL136570
NM_014368
Hs.103137



AL122052
NM_014808
Hs.301283



AK057339
U20582
Hs.2149










Example IV
Individual Genes that are Expressed at Higher than Median Levels and Correlated with the Recurrence of Breast Cancer

DEEPEST, RACGAP1, ZNF145 and MS4A7 were found to each be significantly associated with tumor recurrence. In both the datasets used in FIG. 2, patients were divided into high and low expression groups relative to the overall median for each gene across all patients, and their survival curves were compared (see FIG. 3, which shows Kaplan-Meier disease-free survival curves). The first six graphs in FIG. 3 display the results using the dataset from the 124 patients of Example 1; the X-axis is in months. The second six graphs in FIG. 3 display the results using the dataset from van't Veer et al. with the X-axis in years. The Y-axis for both are “survival probability” as described above. As control, MKI67 and CCNE1, two genes known to be associated with aggressive cancers were analyzed in the same manner.


Example V
Correlation with Tumor Grade

The expression pattern of the Ki67 genes was also found to be strongly correlated with tumor grade. All 67 genes were found to be relatively overexpressed in subjects with high-grade (grade 3) tumors and underexpressed in subjects with low-grade (grades 1 and 2) tumors.


Example VI
Cross-Validation Based on Recurrence Gene Signatures

As shown in FIG. 4, eighty-five selected EP positive (ERP) samples (training dataset) were evaluated for survival probability based upon 141 recurrence gene signatures. The horizontal axis of FIG. 4 is in time (months) and vertical axis is in survival probability.


Table 6 lists the 141 recurrence-associated genes. The sign of the coefficient values in Table 6 corresponds to whether a gene is positively or negatively correlated with breast cancer recurrence. A positive coefficient (score >0) means that the gene is positively correlated in patients with a poor (shorter) survival outcome and negatively correlated coefficients (score <0) mean that the gene is correlated in patients with better (longer) survival outcomes.


The 141 genes were identified from a starting gene pool of 180 genes, wherein the 141 genes had expression levels that correlated with the absence or presence of breast cancer recurrence.

TABLE 6Genes, the expression of which correlate with breast cancer recurrence.genepcoefdescAK0262160.0000124−0.5776051Homo sapiens cDNA: FLJ22563 fis, clone HSI01928BC0128890.00001570.4919877APLP1 | amyloid beta (A4) precursor-like protein 1N301580.0000207−0.5425882ESTsAL1097750.00005790.3522512SH3GL3 | SH3-domain GRB2-like 3AJ2759780.00007780.5278802CTAG1 | cancer/testis antigen 1AI3401910.000156−0.616326HSPC072 | HSPC072 proteinAJ4090650.000171−0.564874LEAP-2 | liver-expressed antimicrobial peptide 2AL1102790.000201−0.5886589H-L(3)MBT | lethal (3) malignant brain tumor l(3)mbt protein (Drosophila) homologAF0571640.000345−0.6014184SLC22A5 | solute carrier family 22 (organic cation transporter), member 5AA3987150.00035−0.550304Homo sapiens cDNA FLJ11529 fis, clone HEMBA1002629AI6090430.0003630.3483688ESTs, Highly similar to T50606 hypothetical protein DKFZp761J107.1 [H. sapiens]AI2770160.000367−0.4165628ESTsAK0247150.000376−0.5755967FLJ21062 | hypothetical protein FLJ21062AL1173960.000385−0.5432668DKFZP586M0622 | DKFZP586M0622 proteinAF2772900.0004650.4119716LOC81501 | DC-specific transmembrane proteinAW2948570.000543−0.5054436LOC51161 | g20 proteinAA8453380.000684−0.3511727FMO5 | flavin containing monooxygenase 5NM_0038420.00074−0.5120187TNFRSF10B | tumor necrosis factor receptor superfamily, member 10bAI4004020.0007540.4795643GBP2 | guanylate binding protein 2, interferon-inducibleAF0352810.00077−0.8496583Homo sapiens clone 23903 mRNA sequenceAK0262150.000794−0.4874715Homo sapiens cDNA: FLJ22562 fis, clone HSI01814BC0007150.0008510.3567116CLECSF9 | C-type (calcium dependent, carbohydrate-recognition domain) lectin, superfamily member 9AI7406710.00095−0.3584972Homo sapiens cDNA FLJ32430 fis, clone SKMUS2001129, weakly similar to NAD-DEPENDENTMETHANOL DEHYDROGENASE (EC 1.1.1.244)AK0249910.00099−0.5133878TRIP8 | thyroid hormone receptor interactor 8AL1367900.00101−0.5095425DKFZp434F1819 | hypothetical protein DKFZp434F1819BC0060000.001070.3974168MGC12536 | hypothetical protein MGC12536AI1394090.001110.5442699CDKN1C | cyclin-dependent kinase inhibitor 1C (p57, Kip2)AA7380430.00134−0.3124779SCAP1 | src family associated phosphoprotein 1AF0556340.001450.3174865UNC5C | unc-5 homolog B (C. elegans)AL1367210.00146−0.5008718DKFZP566K1946 | hypothetical protein DKFZp566K1946U626620.001540.3392979CHIT1 | chitinase 1 (chitotriosidase)BC0123810.001730.3432002FLJ10430 | hypothetical protein FLJ10430AW1955390.00178−0.4498633DDO | D-aspartate oxidaseAI2400720.00187−0.4635729Homo sapiens mRNA full length insert cDNA clone EUROIMAGE 2344436BC0005440.001960.2768966KCNJ8 | potassium inwardly-rectifying channel, subfamily J, member 8AA0294520.00196−0.4884253KCNQ1OT1 | KCNQ1 overlapping transcript 1BC0059390.002250.4954123PTGDS | prostaglandin D2 synthase (21 kD, brain)BC0007480.002260.3983599TUBB4 | tubulin, beta, 4BC0086320.002280.3394714Homo sapiens, clone IMAGE: 3464195, mRNAAI9222040.00246−0.2125933Homo sapiens cDNA FLJ30298 fis, clone BRACE2003172R557640.00247−0.5712032Homo sapiens cDNA FLJ33034 fis, clone THYMU2000236BF4779050.00254−0.4507438ESTsBC0058160.002550.3445679DTX1 | deltex homolog 1 (Drosophila)R560530.00267−0.303724ME3 | malic enzyme 3, NADP(+)-dependent, mitochondrialAK0271910.0027−0.4094706Homo sapiens cDNA: FLJ23538 fis, clone LNG08010, highly similar to BETA2 Human MEN1region clone epsilon/beta mRNAAB0286410.003090.241634SOX11 | SRY (sex determining region Y)-box 11AK0272510.0032−0.3436063FLJ23598 | hypothetical protein FLJ23598AF1533300.0032−0.3901168SLC19A2 | solute carrier family 19 (thiamine transporter), member 2AK0567200.00320.442301Homo sapiens cDNA FLJ32158 fis, clone PLACE6000231BE6714450.00323−0.5541909ESTsNM_0070230.00323−0.383125CAMP-GEFII | cAMP-regulated guanine nucleotide exchange factor IIBC0070450.00330.4541716MLF1 | myeloid leukemia factor 1AY0071140.00335−0.4164539Homo sapiens clone TCCCTA00151 mRNA sequenceAB0331100.00335−0.4816304KIAA1284 | KIAA1284 proteinAI1262710.00381−0.3325846Homo sapiens cDNA FLJ31235 fis, clone KIDNE2004681, moderately similar to Mus musculusperoxisomal long chain acyl-CoA thioesterase Ib (Pte1b) geneU660460.00396−0.3632168Homo sapiens clone 161455 breast expressed mRNA from chromosome XBC0174220.004−0.4540171Homo sapiens, clone MGC: 27375 IMAGE: 4688423, mRNA, complete cdsM765580.00411−0.3517135CACNA1D | calcium channel, voltage-dependent, L type, alpha 1D subunitAA6273580.00412−0.2554568ESTsAF2378130.00415−0.3188455NPD009 | NPD009 proteinBE5011030.00415−0.5697021Homo sapiens cDNA FLJ32173 fis, clone PLACE6000953BC0053620.0042−0.4013912ARHI | ras homolog gene family, member IBC0105630.00427−0.3934675Homo sapiens, clone MGC: 18111 IMAGE: 4152811, mRNA, complete cdsAB0281400.00446−0.5082397TMPRSS5 | transmembrane protease, serine 5 (spinesin)BC0026680.00446−0.3559138PECI | peroxisomal D3,D2-enoyl-CoA isomeraseBC0059480.004480.2903151SMPX | small muscle protein, X-linkedAI3829720.00458−0.5293407TPBG | trophoblast glycoproteinBE5511490.00459−0.4024775ESTsAK0248930.004590.3693678FLJ21240 | hypothetical protein FLJ21240AJ4204900.00461−0.2242605IL20RA | interleukin 20 receptor, alphaAF2616550.00479−0.3113893HMIC | 1,2-alpha-mannosidase ICAF2850890.00493−0.5135262LLT1 | lectin-like NK cell receptorAK0005200.00510.2643999FLJ20513 | hypothetical protein FLJ20513AF1321970.00517−0.4538855PRO1331 | hypothetical protein PRO1331AW1391560.005430.4096147CRMP5 | collapsin response mediator protein-5; CRMP3-associated moleculeU623250.00561−0.5801082APBB2 | amyloid beta (A4) precursor protein-binding, family B, member 2 (Fe65-like)AK0267400.005810.4834703Homo sapiens cDNA: FLJ23087 fis, clone LNG06994, highly similar to AF161368Homo sapiens HSPC105 mRNAAA9855200.00583−0.5217466ESTsAA6338450.00626−0.3464769ESTsNM_0142980.006280.3378903QPRT | quinolinate phosphoribosyltransferase (nicotinate-nucleotide pyrophosphorylase (carboxylating))AF2081110.00638−0.3116422IL17BR | interleukin 17B receptorBC0177330.006410.4047293MRAS | muscle RAS oncogene homologAF2223400.00649−0.4462726ARTS-1 | type 1 tumor necrosis factor receptor shedding aminopeptidase regulatorNM_0054300.006530.3978767WNT1 | wingless-type MMTV integration site family, member 1NM_0022810.006690.319678KRTHB1 | keratin, hair, basic, 1AF2851090.006710.315443-Sep | septin 3BC0185370.00683−0.4326731Homo sapiens, Similar to RIKEN cDNA 1300003P13 gene, clone MGC: 16810 IMAGE: 4340152, mRNA,complete cdsNM_0044300.00685−0.3053692EGR3 | early growth response 3AA5843060.00685−0.5508093GALNT5 | UDP-N-acetyl-alpha-D-galactosamine: polypeptide N-acetylgalactosaminyltransferase 5(GalNAc-T5)BC0121980.006860.394402Homo sapiens, clone MGC: 4408 IMAGE: 2906200, mRNA, complete cdsAI5727370.00702−0.4298417ESTsAI9509850.00726−0.3315239ESTsAI8297930.00734−0.5945521ESTX514200.007350.261491TYRP1 | tyrosinase-related protein 1AF0672230.007470.4465152PDE9A | phosphodiesterase 9AAF1300680.00766−0.4029209SERPINA1 | serine (or cysteine) proteinase inhibitor, clade A (alpha-1antiproteinase, antitrypsin), member 1NM_0035080.007690.2270396FZD9 | frizzled homolog 9 (Drosophila)AK0003970.007750.3705402FLJ10351 | likely ortholog of mouse piwi like homolog 1 (Drosophila)-likeAA5131330.00778−0.4115009ESTs, Weakly similar to 1209280A tropomyosin [H. sapiens]AA8119220.008190.4471232FLJ10140 | hypothetical protein FLJ10140AL3901290.008230.2753104ATP8A2 | ATPase, aminophospholipid transporter-like, Class I, type 8A, member 2BC0143360.00846−0.3339104HHEX | hematopoietically expressed homeoboxAJ0005340.008760.2307116SGCE | sarcoglycan, epsilonAW4448960.008910.2076823ESTsAL1367890.00893−0.4066268DKFZp434F1719 | hypothetical protein DKFZp434F1719AF2076640.008940.2878826ADAMTS1 | a disintegrin-like and metalloprotease (reprolysin type) with thrombospondin type 1 motif, 1AF0162670.00904−0.4665036TNFRSF10C | tumor necrosis factor receptor superfamily, member 10c, decoy without an intracellular domainAA8053170.00906−0.2905281HOXC4 | homeo box C4NM_0329610.009090.4222674PCDH10 | protocadherin 10BC0157060.009090.4310617Homo sapiens, Similar to RIKEN cDNA 2810433K01 gene, clone MGC: 10200 IMAGE: 3909951, mRNA,complete cdsAF1123450.00917−0.4117301ITGA10 | integrin, alpha 10AI6511020.00918−0.361676ESTsBC0014920.009540.2916312CNTFR | ciliary neurotrophic factor receptorBC0141890.00961−0.2332876MGC20702 | hypothetical protein MGC20702AA8901460.009640.486204ESTsAI4591340.009690.310778CD163|CD163 antigenD136260.009690.2913298GPR105 | G protein-coupled receptor 105AW0151710.0098−0.3731535KIAA0022 | KIAA0022 gene productBC0076310.009850.3549781MGC15827 | hypothetical protein MGC15827AF0705870.00988−0.3138695Homo sapiens clone 24741 mRNA sequenceBC0116280.00999−0.3693683EPHX2 | epoxide hydrolase 2, cytoplasmicBF1118830.01040.3647822bA430M15.1 | novel protein (ortholog of rat four repeat ion channel)AF0113330.0105−0.3938612LY75 | lymphocyte antigen 75AB0424100.01060.2578638GPR88 | G-protein coupled receptor 88AB0126430.01060.4935166ALPL | alkaline phosphatase, liver/bone/kidneyBC0005580.0109−0.2313315MAPT | microtubule-assotiated protein tauM818830.011−0.3791012GAD1 | glutamate decarboxylase 1 (brain, 67 kD)AA7652560.0112−0.3259152ESTs, Weakly similar to unnamed protein product [H. sapiens]AL3539480.01130.3592147Homo sapiens mRNA; cDNA DKFZp761P0114 (from clone DKFZp761P0114)AI6861600.0113−0.3480712ESTsBF1120170.01130.2891108KCNE1L | potassium voltage-gated channel, Isk-related family, member 1-likeAI3934760.01140.3087115RBP1 | retinol binding protein 1 , cellularR414290.0115−0.3424425ATP7B | ATPase, Cu++ transporting, beta polypeptide (Wilson disease)L204330.01170.4816472POU4F1 | POU domain, class 4, transcription factor 1BC0018470.01170.3254246MGC4504 | hypothetical protein MGC4504AK0256720.0118−0.4356503FLJ20761 | hypothetical protein FLJ20761BC0147940.0118−0.3526159FLJ20574 | hypothetical protein FLJ20574BC0024150.012−0.2664004GSTT2 | glutathione S-transferase theta 2AK0271790.012−0.4551677ATF7 | activating transcription factor 7AL3539440.0124−0.3365367Homo sapiens mRNA; cDNA DKFZp761J1112 (from clone DKFZp761J1112)NM_0121730.0125−0.4651839FBXO25 | F-box only protein 25
## NOTE.

P and coef are derived from Cox proportional hazard models. Genes with positive coefs are expressed higher in recurrence group, and those with negative coefs are expressed higher in no-recurrence group


Example VII
Prognosis of Recurrence Utilizing the 141 Recurrence Signature Genes

Sixty-six patients having ERP breast cancer (test dataset) were evaluated utilizing the identified 141 signature genes in order to predict survival outcomes, based upon recurrence of the breast cancer. The prognostic results are shown in FIG. 5.


Another group of patients were evaluated; this group contained both ERP and ER negative (ERN) patients, wherein the total number of patients evaluated was 162 (test dataset). The prognostic results for this second group of patients also are shown in FIG. 5.


All references cited herein, including patents, patent applications, and publications, are hereby incorporated by reference in their entireties, whether previously specifically incorporated or not.


Having now fully described this invention, it will be appreciated by those skilled in the art that the same can be performed within a wide range of equivalent parameters, concentrations, and conditions without departing from the spirit and scope of the invention and without undue experimentation.


While this invention has been described in connection with specific embodiments thereof, it will be understood that it is capable of further modifications. This application is intended to cover any variations, uses, or adaptations of the invention following, in general, the principles of the invention and including such departures from the present disclosure as come within known or customary practice within the art to which the invention pertains and as may be applied to the essential features hereinbefore set forth.

Claims
  • 1. (canceled)
  • 2. (canceled)
  • 3. (canceled)
  • 4. (canceled)
  • 5. A method to determine the survival outcome of a breast cancer afflicted subject or determine prognosis of a subject having breast cancer, said method comprising assaying a sample of breast cancer cells of said subject for the expression level(s) of one or more genes listed in Table 2, 3, 4, and/or 6.
  • 6. (canceled)
  • 7. A method of determining the prognosis of a subject having breast cancer correlated with the over or under expression of one or more genes in Table 2, 3, 4, and/or 6 said method comprising assaying for the expression level(s) of said one or more genes in a breast cancer cell from said subject.
  • 8. The method of claim 5 wherein said assaying comprises preparing RNA from said sample.
  • 9. The method of claim 8 wherein said RNA is used for quantitative PCR.
  • 10. The method of claim 5 wherein said assaying comprises using an array.
  • 11. The method of claim 5 wherein said sample is a ductal lavage or fine needle aspiration sample.
  • 12. The method of claim 11 wherein said sample is microdissected to isolate one or more cells suspected of being breast cancer cells.
  • 13. The method of claim 5 wherein said assaying comprises preparing RNA from said sample and optionally using said RNA for quantitative PCR.
  • 14. The method of claim 8 wherein said assaying comprises using an array.
  • 15. The method of claim 5 wherein said sample is a ductal lavage or fine needle aspiration sample, which sample is optionally microdissected to isolate one or more cells suspected of being breast cancer cells.
  • 16. The method of claim 7 wherein said assaying comprises preparing RNA from said cell and optionally using said RNA for quantitative PCR.
  • 17. The method of claim 7 wherein said assaying comprises using an array.
  • 18. The method of claim 7 wherein said cell is present in a ductal lavage or fine needle aspiration sample, which sample is optionally microdissected to isolate one or more cells suspected of being breast cancer cells.
  • 19. A method to determine the grade of breast cancer in a subject comprising assaying a sample of breast cancer cells of said subject for the expression level(s) of one or more genes listed in Table 2.
  • 20. A method to determine therapeutic treatment for a breast cancer patient based upon said patient's expected survival, said method comprising determining a survival outcome for said patient by assaying a sample of breast cancer cells from said patient for the expression level(s) of one or more one genes listed in Table 2, 3, 4, and/or 6; and selecting the appropriate treatment for a patient with such a survival outcome.
RELATED APPLICATIONS

This application claims benefit of priority from U.S. Provisional Patent Application 60/479,963, filed Jun. 18, 2003, and 60/545,810, filed Feb. 18, 2004, both of which are incorporated by reference as if fully set forth.

Provisional Applications (2)
Number Date Country
60545810 Feb 2004 US
60479963 Jun 2003 US