The present invention relates to a method for evaluating the probability of survival for an individual suffering from endometrial carcinoma. In another aspect, the present invention relates to the stratification of therapy regimen of endometrial tumor, ovarian cancer, breast cancer, non-small lung cancer or hormone refractory prostate cancer therapy in an individual or monitoring therapeutic efficacy in an individual suffering from the same based on the expression status of STMN1 gene or protein. Moreover, the present invention relates to a kit for use in any of the above referenced methods comprising a means for determining amplifications and deletions of chromosomal regions 3q26.32 and 12p12.1, determining alterations of the gene expression profile of the genes (gene signature): upregulation of the genes PLEKHK1, ATP10B, NMU, MMP1, ATAD2, NETO2, TNNI3, PHLDA2, OVOL1 and down-regulation of the genes: NDP, KIAA1434, MME, CFH, MOXD1, SLC47A1, RBP1, PDE8B, ASRGL1, ADAMTS19, EFHD1, ABCA5, NPAS3, SCML1, TNXB, ENTPD3, AMY1A, ENPP, RASL11B, PDZK3, or the expression status of the STMN1 gene or protein, respectively. Finally, the present invention provides a method for predicting the response to taxanes in an individual suffering from a disease treated with the taxanes based on the expression status of the STMN1 gene or protein.
With a 2-3% lifetime risk among women, endometrial cancer is the most common pelvic gynecologic malignancy in industrialized countries, and the incidence is increasing (Amant F et al. (2005), Lancet, 366:491-505). Approximately 75% of cases are diagnosed with the tumor confined to the uterine corpus, but 15%-20% of these recur after primary surgery with limited respond to systemic therapy. In light of these recurrences, patients with localized endometrial cancer have 2 major needs: (1) adjuvant therapies that will reduce the recurrence rate, and (2) the ability to target these therapies to the patients most likely to recur. In addition, women with metastatic disease require effective systemic therapy.
These needs, for effective systemic therapies and reliable prognostic markers, have been only partly addressed. The most common basis for determining risk of recurrent disease has been the categorization of endometrial cancer into two subtypes. The majority are type I, associated with good prognosis, low stage and grade, and endometrioid histology. In contrast, type II cancers are characterized by high stage and grade, non-endometrioid histology, and poor prognosis. However, the prognostic value of this distinction is limited as up to 20% of type I cancers recur, while half of type II cancers do not.
The molecular basis of the distinction between type I and II cancer is only partially understood. Type I cancer is associated with hyperestrogenic risk factors, is more often estrogen and progesterone receptor positive, diploid, microsatellite unstable, and KRAS or PTEN mutant. Type II cancer is more often aneuploid and harbors alterations in CDKN2A, TP53, and ERBB2. Such molecular alterations are of prognostic value but have not provided a basis for improved therapy Lax SF, 2004, Virchows Arch, 444:213-223). Hormone receptor status influences the choice of treatment in metastatic disease, but most aggressive tumors are receptor negative.
Recently, Saal et. al. PNAS, 2007, 104, 18, 7564 to 7569 report on observations that poor prognosis in carcinoma is associated with a gene expression signature of apparent PTEN tumor suppressant pathway activity. That is, expression of STMN1 has been shown previously to correlate with PI3K activity in breast cancer and can be measured by immunohistochemistry in paraffin-embedded tissue.
The present inventors hypothesized that tumors with an aggressive phenotype are likely to be distinguished by underlying genetic alterations reflected in distinct transcriptional signatures, and investigated whether tumors that recur share transcriptional signatures that suggest shared underlying genetic alterations.
Endometrial cancer is the most frequent gynaecological cancer in industrialised countries. Although the majority have a good prognosis, up to 20% recurs. To date there are few markers available to predict response to treatment of metastatic endometrial cancer. Patients with tumors expressing estrogen- and progestagen receptors have the best response to antihormonal treatment. Still, more markers are needed to predict response to other therapy modalities in patients with metastatic endometrium cancer.
Hence, the first object of the present invention is to provide methods allowing differentiation of endometrial carcinoma and other types of carcinoma in an individual in vie of treatment regimen, in particular, with respect to chemotherapy. Further, the present invention aims to provide a method of evaluating the probability of survival for an individual suffering from endometrial carcinoma or the clinical outcome thereof as well as providing a method for the stratification of endometrial tumor, ovarian cancer, breast cancer, non-small lung cancer or hormone refractory prostate cancer therapy in an individual or monitoring therapeutic efficacy in an individual suffering therefrom with respect to the usefulness of chemotherapy.
In a first aspect, the present invention relates to a method for differentiation of endometrial carcinoma in an individual for the responsiveness or susceptibility of whether said individual is responsive or susceptible to the treatment with chemotherapeutic drugs, in particular, chemotherapeutic drugs of disrupting microtubule function, comprising the steps of determining alterations, in particular, amplifications and deletions, of chromosomal regions 3q26.32 and 12p12.1, alterations of the gene expression profile of the genes (gene signature): upregulation of the genes PLEKHK1, ATP10B, NMU, MMP1, ATAD2, NETO2, TNNI3, PHLDA2, OVOL1 and down-regulation of the genes: NDP, KIAA1434, MME, CFH, MOXD1, SLC47A1, RBP1, PDE8B, ASRGL1, ADAMTS19, EFHD1, ABCA5, NPAS3, SCML1, TNXB, ENTPD3, AMY1A, ENPP, RASL11B, PDZK3, or the expression status of the STMN1 gene or protein, and determining the susceptibility or responsiveness of said individual to of a chemotherapeutic treatment, in particular, a chemotherapeutic treatment with a chemotherapy drug acting by disrupting microtubuli function, in particular, of taxanes.
In another aspect, the present invention relates to method for evaluating the probability of survival or the clinical outcome of an individual, intended to be or treated with chemotherapy drugs, in particular, taxanes whereby said individual suffering from endometrial carcinoma comprising the step of
Furthermore, the present invention provides a method for the stratification of a chemotherapeutic therapy of endometrial tumor, ovarian cancer, breast cancer, non-small lung cancer or hormone refractory prostate cancer in an individual or monitoring chemotherapeutic efficacy of said diseases in an individual comprising the steps of determining the expression status of the STMN1 gene or protein and stratifying the therapy or monitoring the efficacy of chemotherapy of the endometrial tumor, ovarian cancer, breast cancer, non-small lung cancer or hormone refractory prostate cancer in said individual.
In another aspect, the present invention relates to a kit for use in providing a differentiation of endometrial carcinoma in an individual, for the stratification of endometrial tumor therapy in an individual, monitoring therapeutic efficacy in an individual, or for evaluating the probability of survival for an individual suffering from endometrial carcinoma whereby said individual is treated or is intended to be treated with a chemotherapeutic drug comprising means for determining determining alterations, in particular, amplifications and deletions, of chromosomal regions 3q26.32 and 12p12.1, alterations of the gene expression profile of the genes (gene signature): upregulation of the genes PLEKHK1, ATP10B, NMU, MMP1, ATAD2, NETO2, TNNI3, PHLDA2, OVOL1 and down-regulation of the genes: NDP, KIAA1434, MME, CFH, MOXD1, SLC47A1, RBP1, PDE8B, ASRGL1, ADAMTS19, EFHD1, ABCA5, NPAS3, SCML1, TNXB, ENTPD3, AMY1A, ENPP, RASL11B, PDZK3, or the expression status of the STMN1 gene or protein.
Moreover, the present invention relates to a method for predicting the response or outcome of therapy with taxanes in an individual treated therewith based on the expression status of the STMN1 gene or protein.
The above methods are particularly useful for stratification of the therapy and for monitoring the therapy when treating metastatic cancer, in particular metastatic endometrial cancer.
Finally, the present invention relates to a method for the stratification of therapy or for monitoring the efficacy of the therapy based on chemotherapeutics, like PI3K inhibitors, Akt inhibitors, mTOR inhibitors or PTEN activators, in particular, chemotherapeutics disrupting microtubule function, like taxanes comprising the step of determining the expression status of the STMN1 gene or protein.
In a first aspect, the present invention relates to a method for the differentiation of endometrial carcinoma in an individual for the responsiveness or susceptibility of whether said individual is responsive or susceptible to the treatment with chemotherapeutic drugs, in particular, chemotherapeutic drugs of disrupting microtubule function, comprising the step of determining alterations, in particular, amplifications and deletions, of chromosomal regions 3q26.32 and 12p12.1, alterations of the gene expression profile of the genes (gene signature): upregulation of the genes PLEKHK1, ATP10B, NMU, MMP1, ATAD2, NETO2, TNNI3, PHLDA2, OVOL1 and down-regulation of the genes: NDP, KIAA1434, MME, CFH, MOXD1, SLC47A1, RBP1, PDE8B, ASRGL1, ADAMTS19, EFHD1, ABCA5, NPAS3, SCML1, TNXB, ENTPD3, AMY1A, ENPP, RASL11B, PDZK3, or the expression status of the STMN1 gene or protein, and determining the susceptibility or responsiveness of said individual to of a chemotherapeutic treatment, in particular, a chemotherapeutic treatment with a chemotherapy drug acting by disrupting microtubuli function, in particular, of taxanes.
That is, it is recognized that two major groups of tumors can be distinguished in patients suffering from endometrial carcinoma. Namely, two clusters allow to differentiate between two major groups of tumors whereby these clusters identify a two-fold or higher change for 138 significant genes of which 64 where upregulated and 74 downregulated in cluster 2. A set of 29 genes, validated by quantitative RT-PCR, predicted the clusters with 100% accuracy. Said clusters allow to differentiate the susceptibility or responsiveness of an individual in need of a treatment of endometrial cancer and other types of cancer as specified herein to chemotherapeutic drugs, in particular, drugs disrupting the microtubule function, like taxanes.
The expression clusters identified herein have strikingly different clinical and histopathologic characteristics. Cluster 2 contained more aggressive tumors containing almost all type II tumors. In addition, patients with tumors in Cluster 2 had significantly poorer recurrence-free survival. Segregation into Cluster 2 predicted recurrence better than known means in the art, like International Federation of Gynecology and Obstetrics (FIGO) stage, histologic grade, number of mitosis, presence of a non-endometrioid histologic subtype, tumor necrosis and vascular invasion.
Thus, the present inventors recognized that determining alterations, in particular, amplifications and deletions, of chromosomal regions 3q26.32 and 12p12.1, alterations of the gene expression profile of the genes (gene signature): upregulation of the genes PLEKHK1, ATP10B, NMU, MMP1, ATAD2, NETO2, TNNI3, PHLDA2, OVOL1 and down-regulation of the genes: NDP, KIAA1434, MME, CFH, MOXD1, SLC47A1, RBP1, PDE8B, ASRGL1, ADAMTS19, EFHD1, ABCA5, NPAS3, SCML1, TNXB, ENTPD3, AMY1A, ENPP, RASL11B, PDZK3, as well as determining the expression status of STMN 1 gene or protein in an individual in vivo or in vitro allows for the diagnosis or differentiation of endometrial carcinoma in said individual.
According to the present invention, the methods disclosed herein relates to in vitro and/or in vivo methods, respectively.
In a preferred embodiment, the method or differentiation of endometrial carcinoma in an individual comprise the steps of determining the PI3K activity in patients having aggressive endometrial carcinoma, in particular, based on the alterations in 3q26.32 or on the expression status of STMN1 gene or protein.
In another embodiment, it is preferred that the expression status of the STMN1 gene or protein is determined.
The method of the present invention allows to differentiate between high grade aggressive phenotype of endometrial cancer and low grade phenotype of endometrial cancer and thus, allow to differentiate or determine the susceptibility or responsiveness of an individual in need of a treatment of endometrial cancer and other types of cancer as specified herein to chemotherapeutic drugs, in particular, drugs disrupting the microtubule function, like taxanes.
As used herein, the term “taxanes” refers to diterpenes having cytostatic activity. Examples of suitable taxanes include paclitaxel and docetaxel. The skilled person is well aware of suitable forms of taxanes including salts and solvates thereof.
Hence, the present invention relates to methods allowing differentiation of endometrial carcinoma, in particular allowing to differentiate between low grade and high grade aggressive phenotype in endometrial carcinoma based on the STMN1 expression for determining the treatment regimen or the clinical outcome in an individual suffering therefrom. The present invention is directed to the prognosis as well as to the stratification of endometrial tumors and its therapy with respect to chemotherapeutic drugs. That is, in one aspect, the present invention relates to endometrial carcinoma and the importance of the PI3K pathway in patients having aggressive endometrial cancer. The STMN1 expression correlates with PI3K scores and, in addition, high STMN1 expression is associated with poor recurrence free survival and with poor recurrence free and overall survival in patients suffering from endometrial carcinomas. It is demonstrated herein that high STMN1 expression represents an independent prognostic indicator allowing to differentiate between high grade aggressive phenotype and low grade phenotype of endometrial cancer and the clinical outcome or the susceptibility or responsiveness of an individual in need of a treatment of endometrial cancer and other types of cancer as specified herein to chemotherapeutic drugs, in particular, drugs disrupting the microtubule function, like taxanes. In particular, high STMN1 expression is associated with poor prognosis and the otherwise low risk endometrioid subgroup.
The present inventors recognized that PI3K activity associates with poor prognosis, thus, indicating that measuring PI3K activity allows to improve prognostication of localized endometrial cancer.
The present invention covers the determination of STMN1 expression in methods allowing the differentiation of endometrial carcinoma as well as stratification of endometrial tumors and its therapy as well as monitoring the chemotherapy. Furthermore, the present invention provides a method for evaluating the probability of survival as well as methods for providing a prognosis of a subject afflicted with endometrial cancer based on PI3K activity and/or STMN1 expression and chemotherapy.
In further aspects, the present invention relates to methods including determining amplifications and deletions of specific chromosomal regions, like 3q and 12p, in particular 3q26.32 and 12p12.1 as detailed herein. In particular, the amplifications and deletions outlined in FIG. 2 allows to differentiate individuals afflicted with endometrial carcinoma in two clusters, namely cluster 1 and cluster 2 having significant differences in disease-free survival. Preferably, the methods according to the present invention includes determining expression of STMN1 in combination with determining at least one of the amplifications or deletions in the chromosomal regions identified herein or determining the gene signature of the genes PLEKHK1, ATP10B, NMU, MMP1, ATAD2, NETO2, TNNI3, PHLDA2, OVOL1 and NDP, KIAA1434, MME, CFH, MOXD1, SLC47A1, RBP1, PDE8B, ASRGL1, ADAMTS19, EFHD1, ABCA5, NPAS3, SCML1, TNXB, ENTPD3, AMY1A, ENPP, RASL11B, PDZK3.
To conclude, the present invention relates in another aspect to a method for evaluating the probability of survival for a patient with endometrial cancer, said method being characterized in that it comprises measuring the level or expression of STMN1 on nucleic acid or amino acid level in a sample obtained from said patient.
Moreover, in another preferred embodiment, the method according to the present invention comprises determining the expression status of STMN1. It has been recognized that high STMN1 expression is associated with poor recurrence-free survival and over survival in patients suffering from endometrial carcinoma. In particular, the STMN1 expression allows to differentiate between high grade aggressive phenotype and low grade phenotype of endometrial carcinoma whereby high STMN1 is associated with high grade aggressive phenotype of endometrial carcinoma and, in addition, allows to determine the susceptibility or responsiveness of an individual in need of a treatment of endometrial cancer and other types of cancer as specified herein to chemotherapeutic drugs, in particular, drugs disrupting the microtubule function, like taxanes.
In another preferred embodiment, the methods according to the present invention comprises the step of determining expression of the STMN1 gene in combination with determining alterations, in particular, the amplifications or deletions, in the chromosomal regions 3q 26.32 and 12p12.1, or altered expression of the gene signature of the genes: PLEKHK1, ATP10B, NMU, MMP1, ATAD2, NETO2, TNNI3, PHLDA2, OVOL1 (upregulation) and of the genes: NDP, KIAA1434, MME, CFH, MOXD1, SLC47A1, RBP1, PDE8B, ASRGL1, ADAMTS19, EFHD1, ABCA5, NPAS3, SCML1, TNXB, ENTPD3, AMY1A, ENPP, RASL11B, PDZK3 (downregulation).
Another aspect relates to a method for the stratification of the therapeutic regimen of a subject with endometrial carcinoma comprising
That is, the present inventors recognized that not only for the stratification of endometrial cancer and for monitoring therapeutic efficacy in the treatment of endometrial tumors and cancers but also in ovarian cancer, breast cancer, non-small lung cancer or hormone refractory prostate cancer therapy, STMN1 is a valuable biomarker.
In particular, the present inventors aimed in demonstrating that STMN1, also known as Stathmin, expression predicts the response to taxanes or chemotherapeutic drugs disrupting the microtubular function in metastatic endometrial cancer. Hence, Stathmin expression is useful as a marker for the treatment of metastatic endometrial cancer but also in endometrial cancer in general and ovarian cancer, breast cancer, non-small lung cancer or hormone refractory prostate cancer.
The method for the stratification of the therapeutic regimen or monitoring the therapeutic regimen or monitoring the therapeutic efficacy of an individual suffering from endometrial cancer, ovarian cancer, breast cancer, non-small lung cancer or hormone refractory prostate cancer comprises the step of determining the level or amount of STMN1 is a sample of said individual and determining the therapeutic regimen or strategy or monitoring the therapeutic efficacy based on the level or amount of STMN1, in particular, with respect to chemotherapeutic drugs, in particular, chemotherapeutic drugs of disrupting microtubular function, like taxanes.
Preferably, the STMN1 expression status is determined on nucleic acid or amino acid level in said individual.
The skilled person is well aware of suitable methods for determining the expression status of the gene STMN1 or the amplification and deletions in the chromosomal regions 3q 26.32 and 12p12.1, as well of determining alterations of the genes PLEKHK1, ATP 10B, NMU, MMP1, ATAD2, NETO2, TNNI3, PHLDA2, OVOL1 (upregulation) and of the genes: NDP, KIAA1434, MME, CFH, MOXD1, SLC47A1, RBP1, PDE8B, ASRGL1, ADAMTS19, EFHD1, ABCA5, NPAS3, SCML1, TNXB, ENTPD3, AMY1A, ENPP, RASL11B, PDZK3 (downregulation), respectively.
Preferred embodiments include the detection of nucleic acid level using PCR methods or hybridisation methods using suitable marker molecules.
On protein level, determining the expression status of the gene STMN1 may be effected by using appropriate antibodies and systems comprising the same. Suitable methods including ELISA, Western blot, immunohistochemical or immunofluorescence detection.
In another aspect, a kit for use in providing a differentiation of endometrial carcinoma in an individual, for the stratification of endometrial tumor therapy in an individual, monitoring therapeutic efficacy in an individual, or for evaluating the probability of survival for an individual suffering from endometrial carcinoma to allow to differentiate or determine the susceptibility or responsiveness of an individual in need of a treatment of endometrial cancer and other types of cancer as specified herein to chemotherapeutic drugs, in particular, drugs disrupting the microtubule function, like taxanes comprising means for determining amplifications and deletions of chromosomal regions 3q26.32 and 12p12.1, the expression status of the STMN1 gene or protein or means for determining amplification and deletions whereby said amplifications (upregulation) and deletions (downregulations) are amplifications of the genes: PLEKHK1, ATP10B, NMU, MMP1, ATAD2, NETO2, TNNI3, PHLDA2, OVOL1 and deletions of the genes: NDP, KIAA1434, MME, CFH, MOXD1, SLC47A1, RBP1, PDE8B, ASRGL1, ADAMTS19, EFHD1, ABCA5, NPAS3, SCML1, TNXB, ENTPD3, AMY1A, ENPP, RASL11B, PDZK3 is provided.
In another aspect, said kit comprises means for determining the PI3K activity in patients having aggressive endometrial carcinoma.
Particularly preferred, said kit according to the present invention is suitable for providing diagnosis or differentiation of endometrial carcinoma in an individual or for the stratification of the therapeutic regiment of monitoring the therapeutic efficacy comprising means for detecting STMN1 expression status.
Said kit is particularly useful for predicting the response to taxanes in an individual when treating the same considering a therapeutic regimen using taxanes in said individuals. In particular in case of the treatment of metastatic endometrial cancer, the method and kits according to the present invention are useful for stratifying the therapy thereof. For example, when taxanes are used for the treatment of metastatic cancer, like metastatic endometrial cancer, determining the STMN1 status allows to stratify and to diagnose therapeutic success of taxanes treatment.
That is, there are a few markers available to predict response to treatment of metastatic endometrial cancer. Patients with tumors expressing estrogen and progestagen receptors have the best response to antihormonal treatment. However, more markers are needed to predict the response to other therapy modalities in patients with metastatic endometrial cancer. It has been demonstrating herein that the level of stathmin expression (STMN1 expression) allows to predict response to tubuli stabilizing chemotherapy in cancer, like endometrial cancer. A typical example of tubuli stabilizing therapy includes Taxol, Taxotere, Eleutherobin, Sarcodicytin A, Sarcodicytin B, Epothilone A, Epothilone B, Discodermolide, Laulimalide, Isolaulimalide, Ixabepilone, Vinblastin, Vinkristin, Vinorelbin.
Finally, the present invention relates to a method for stratification of endometrial tumor or endometrial cancer, ovarian cancer, breast cancer, non-small lung cancer or hormone refractory prostate cancer therapy in an individual or monitoring therapeutic efficacy in an individual whereby the therapy, in particular, the endometrial tumor therapy based on PI3K inhibitors AKT inhibitors or mTOR inhibitors or PTEN activators comprising the step of determining the expression status of the STMN1 gene or protein and stratifying the therapy or monitoring the efficacy of the therapy accordingly.
For the primary investigation series, primary endometrial carcinomas were immediately frozen during hysterectomies conducted from 2001-2003. All samples were reviewed by a pathologist according to published criteria (Scully RE et al. (1994) Histological typing of female genital tract tumours. International histological classification of tumours. World Health Organization. Springer-Verlag, Berlin Heidelberg). Treatment included bilateral salpingo-oophorectomy and pelvic lymphadenectomy. Adjuvant therapy was recommended for patients with FIGO surgical stage 10B or higher disease or non-endometrioid histology. Patients were followed from primary surgery until June 2007 or death, with a median follow-up for survivors of 3.6 years (range 0.8-5.5). Deaths not attributable to endometrial cancer were censored. No patient was lost to follow-up.
RNA was extracted from biopsies with at least 50% (usually >80%) tumor content using the RNeasy kit (Qiagen). Quality and yield were assessed by agarose electrophoresis, the Agilent Bioanalyser 2100, and spectrophotometry. RNA was prepared in 2 batches and hybridized to Agilent 21K and 22K arrays respectively, according to manufacturer's instructions (www.agilent.com). Arrays were scanned using the Agilent Microarray Scanner Bundle.
Signal intensities were determined using J-Express (www.molmine.com) and filtered to remove genes with signal intensities below 2 standard deviations over background in either channel (Cy5, Cy3) in more than 30% of samples. Batch adjustment was performed as previously described (Engelsen IB et al. (2008) Br J Cancer 98:1662-1669). Genes were mean-centered across the tumor set.
Hierarchical clustering was performed using the 3500 genes with highest variance using weighted average linkage (WPGMA) and Pearson correlation as similarity measures. Clustering with more or fewer genes gave stable results (data not shown). A SAM analysis using these clusters as class labels identified 138 significantly changed genes, of which 29 were selected for their combined discriminatory power as described in SI Methods. Messenger RNA levels for these 29 genes and PTEN were validated by quantitative PCR using the TaqMan Low Density Array (Applied Biosystems) according to manufacturer's instructions (Engelsen IB et al. (2008) Br J Cancer 98:1662-1669).
For the external dataset (Affymetrix U133+2 arrays), individual probes were sequence-matched against Aceview (NCBI35) (Carter SL et al. (2006) Nat Genet 38:1043-1048) to construct transcript-level probesets. Summary expression levels were then derived by batch-normalization across samples via RMA (Irizarry RA et al. (2003) Nucleic Acids Res 31:el 5).
The PI3K score was obtained by comparing previously published expression data of 9 replicate transfections of activated PIK3CA to 5 GFP controls, and includes the 495 genes surpassing a Bonferroni-corrected 2-sided t-test p-value of 0.05. To evaluate this signature, expression data for each gene were normalized to a common mean and scaled to the same standard deviation. For each sample, the activation score is the sum of genes significantly upregulated in the cells with activated PIK3CA (relative to the cells with GFP control) minus genes significantly downregulated in those cells.
Genomic DNA was extracted from surgically dissected, fresh-frozen primary tumors and from nine cell lines: Ishicawa, Hec1A, KLE, AN3-CA, EFE184, MFE-280, MFE-296, MFE-319, RL-95-2. Tumors were needle dissected to ensure 80% purity.
PIK3CA, KRAS and PTEN were sequenced. Genomic DNA was analyzed by SNP arrays interrogating 116,204 SNP loci (Affymetrix) and the GISTIC algorithm, as previously described in Beroukhim R et al. (2007) Proc Natl Acad Sci USA 104:20007-20012. SNP, gene, and cytogenetic band locations are based on the hg16 (July 2003) genome build (genome.ucsc.edu).
For relations of molecular data to clinical phenotype, Pearson's chi-square- (χ2), Fisher's exact-, Mann-Whitney-, or Kruskal-Wallis tests were used as appropriate. P-values represent 2-sided tests except when testing the 1-sided hypothesis that 3 qamp correlates with measures of PI3K activation. Univariate survival analyses were performed by the Kaplan-Meier method. The log-rank (Mantel-Cox) test with Bonferroni correction was used to compare survival curves for different categories. Variables with significant impact on survival (p<0.05) were further examined by log-minus-log plot before incorporation in the Cox' proportional hazards regression model.
Genome-wide expression and clinical and histopathologic data from a random sampling of 57 endometrial carcinomas in a population-based tissue bank of gynaecologic cancer in Hordaland County, Norway, were collected. The characteristics of these patients were not significantly different from all patients diagnosed with endometrial carcinoma in a ten-year period from the same region, see Table 1.
An unsupervised analysis of these data distinguished two major groups of tumors (Clusters 1 and 2). SAM analysis (Tusher V G, Tibshirani, R & Chu G (2001) Proc Natl Acad Sci USA 98:5116-5121) between these clusters identified a two-fold or higher change for 138 significant genes, of which 64 were upregulated and 74 downregulated in Cluster 2. A set of 29 genes, validated by quantitative RT-PCR, predicted the clusters with 100% accuracy.
The two clusters had strikingly different clinical and histopathologic characteristics. Cluster 2 contained more aggressive tumors, with higher International Federation of Gynecology and Obstetrics (FIGO) stage, histologic grade, number of mitoses, presence of non-endometrioid histologic subtype, tumor necrosis and vascular invasion, (p<0.001 for presence of any of these; Table 2). Cluster 2 contains almost all the type II tumors (p<0.001) (Table 2), but it also contains almost one-third of the type I tumors, and these have more vascular invasion, necrosis, and frequent mitoses than the type I tumors in Cluster 1 (p=0.01). The 29-gene summary set was also significantly correlated with aggressive cancer (Table 2).
Most prominently, patients with tumors in Cluster 2 had significantly poorer recurrence-free survival (p=0.05). Segregation into Cluster 2 predicted recurrence better than FIGO stage, histologic subtype, or receptor status, and slightly poorer than grade, but did not exhibit independent prognostic impact, most likely due to the limited number of cases and events.
To identify the underlying somatic changes distinguishing aggressive tumors with the Cluster 2 signature, a genome-wide survey of copy-number changes and LOH among 84 tumors was performed. The majority exhibit a small number of amplifications (median of 4 in each tumor) and even fewer deletions (median of 1). Nevertheless, virtually every region of the genome is amplified or deleted in at least 1 tumor.
To distinguish copy-number changes associated with endometrial cancer from potentially random events, we applied the statistical method Genomic Identification of Significant Targets In Cancer (GISTIC) (Beroukhim R et al. (2007) Proc Natl Acad Sci USA 104:20007-20012). GISTIC assigns each region of the genome 2 G-scores, each representing the combined frequency and amplitude of either local amplifications or deletions. It then compares these to similar scores generated from random permutations of the data to determine False Discovery Rate q-values, representing the likelihood of obtaining the observed G-scores from chance events alone. The G-scores tend to be larger for amplifications than deletions due to the greater prevalence of amplifications. Conversely, deletions attain statistical significance (using a q-value threshold of 0.25 at lower prevalence due to their overall infrequency.
11 significantly amplified and 13 significantly deleted regions of the genome (Table 3) have been found. For each we selected the peak region, with the highest frequency and amplitude of events, as the region most likely to contain a cancer gene target was selected. Known oncogenes are located within these peaks for 8 amplified regions and known tumor suppressors are located within deletion peaks on chromosomes 1 and 3 (Table 3), but functional data tying any of these genes to endometrial carcinogenesis are lacking. Also, 14 regions contain no known cancer genes. These usually represent infrequent events (<17% of tumors), with the exception of lq amplification, where the gene target is unclear due to the large size of the amplicon. The consistent breadth of this amplicon, in fact, may suggest more than one target. LOH generally reflects deletions, with the exception of prevalent copy-neutral LOH on 10q containing the known endometrial tumor suppressor PTEN.
Amplifications of KRAS and PIK3CA Associate with Poor Prognosis
Among the 11 significant amplifications, only 2 (3q26.32 and 12p12.1) are associated with recurrence-free survival (in both cases poor survival) after correction for multiple hypotheses (Table 4). The amplifications due to the low prevalence of deletions were considered only. Amplification of 3q26.32 (3qamp) is also associated with non-endometrioid histology (44% vs 11% prevalence; p=0.02) and high grade (p<0.001). The association between 12p12.1 amplification and poor survival is surprising because mutations of KRAS, which is within the peak region, are known to associate with better survival. However, KRAS in 64 tumors were sequenced and found none of the 12p12.1 amplified samples had mutant KRAS, although mutations were seen in 4 unamplified samples. Amplification of 12p12.1 is also associated with high grade (p=0.02) and FIGO stage (p=0.04). Although 3q26.32 and 12p12.1 tended to be amplified in the same tumors (p=0.03), they usually did not coincide. We directed further analyses at 3qamp because all the samples with this amplification segregated into expression Cluster 2 (p=0.01), suggesting that the amplification could be associated with the Cluster 2 transcriptional profile (see below).
Integrated analyses associate markers of PI3 kinase activation with aggressive cancer
It should be investigated whether 3 qamp leads to an aggressive phenotype through activation of PIK3CA. Although PIK3CA has not been shown to be the 3qamp target, its suspected for four reasons: (1) PIK3CA is 1 of 36 genes within the peak region; (2) tumors with 3 qamp overexpress PIK3CA compared to unamplified tumors (p=0.003); (3) similar amplifications in ovarian cancer act through PIK3CA Shayesteh L et al. (1999) Nat Genet 21:99-102; and (4) the PI3 kinase (PI3K) pathway is frequently aberrant in endometrial cancer, including point mutations in PIK3CA.
Therefore it has been looked for wider effects of PIK3CA activation in the transcriptome of tumors with 3qamp. Published data (Potti A et al. (2006) Nat Med 12:1294-1300) from cell lines transfected with mutationally activated PIK3CA have been used to define a PI3K activation score (PI3K score), representing the expression levels of genes that correlate with activated PIK3CA (see Methods). Tumors with 3qamp scored higher than unamplified samples (p=0.05). However, the impact of this finding is limited by its borderline statistical significance and by the possibility that the PI3K score may not reflect PI3K activation generally, but only in the model systems in which it was measured.
To corroborate this finding it has been analyzed whether samples with 3qamp have an expression profile opposite that induced by PI3K pathway inhibition. To that end, the 50 most overexpressed and underexpressed genes in samples with 3qamp relative to unamplified samples have been queried using the Connectivity Map (Lamb J et al. (2006) Science 313:1929-1935). Among 164 small molecules represented in the Connectivity Map, the PI3K inhibitor LY-294002 (Vlahos CJ et al. (1994) J Biol Chem 269:5241-5248) had an expression signature most significantly anticorrelated with the 3qamp signature (SI FIG. 2B-C, p=0.003). LY-294002 is known to bind to additional kinases, raising the possibility that this anticorrelation is due to non-specific effects. The anticorrelation between the 3qamp signature and inhibitors of adenylate cyclase and Hsp90 also suggests potentially complex effects of the amplicon. Nevertheless, the findings that the 3qamp signature correlates with a PI3K activation signature and anticorrelates with the signature of a PI3K inhibitor support the hypothesis that one of the effects of 3qamp may be to increase PI3K activity.
Further, the correlation between PIK3CA amplification and the PI3K score in an independent expression dataset has been validated. First, amplification of 3q26-27 from local gene expression levels has been inferred, as reflected in a ‘functional amplification’ (FA) score. As expected, samples determined to have 3qamp by SNP array analysis also had high 3q26-27 FA scores (p<10−5), confirming the score as a meaningful assessment of amplification status. We then inferred 3q26-27 amplification levels in a publicly available expression dataset of 134 endometrial tumors (http://expo.intgen.org/geo/home.do). The correlations between 3qamp and both PIK3CA overexpression and the PI3K score validated (p=2×10−10 and 7×10−5, respectively).
In addition, the correlations between aggressive phenotype and both PIK3CA amplification and the Cluster 2 signature in this independent dataset has been validated. Although survival data were unavailable, both available markers of poor survival, high grade and non-endometrioid subtype, correlated with high 3q26-27 FA scores (p=0.001 and 0.005, respectively; and with high values of the 29 gene summary predictor for membership in Cluster 2 (p=3×10−4 and 0.004, respectively.
The finding that both PIK3CA amplification and the Cluster 2 expression profile indicate aggressive tumors, coupled with the association between PIK3CA amplification and the in vitro PI3K activation signature, suggested that the broader set of aggressive tumors in Cluster 2 might share the in vitro PI3K activation signature. This appears to be true: tumors in Cluster 2 without PIK3CA amplification have significantly higher PI3K scores than tumors in Cluster 1 (p<0.001) and equal to tumors with amplification of PIK3CA. Moreover, the Cluster 2 signature is highly anticorrelated with the signature of treatment with LY-294002 (p=0.02). Furthermore, tumors with high PI3K scores are associated with poor survival (p=0.03) and other markers of aggressive phenotype in both the test and validation datasets (p=0.01 and 0.001, respectively).
One possible cause of overexpression of the PI3K activation signature among tumors without PIK3CA amplification is decreased expression of the downstream PI3K pathway member PTEN. Decreased PTEN expression was associated with increased PI3K scores in both our test and validation datasets (p<0.001 and p=0.03 respectively), regardless of PIK3CA amplification status. Decreased PTEN expression was also associated with markers of aggressive disease (p=0.02).
Conversely, among the 45 tumors with expression data that we sequenced for PTEN, mutations did not associate with high PI3K scores (p=0.6). On the contrary, more mutations in the non-aggressive Cluster 1 than Cluster 2 (p=0.04) have been observed.
Overexpression and mutation of PIK3CA also appear to have different implications. Significantly higher PIK3CA expression in tumors with aggressive features, including those without PIK3CA amplification (p=0.05 and 0.0009 among test and validation data) have been found. However, among the 41 tumors with expression data that were sequenced for PIK3CA, mutations did not associate with high PI3K scores (p=0.8) or features of aggressive disease (p=0.5). Further, it cannot be confirm the finding that exon 20 mutations correspond to aggressive tumors (Catasus L et al. (2008) Mod Pathol 21:131-139). Although PIK3CA mutations have previously been noted primarily in endometrioid cancers (Ollikainen M et al. (2007) Int J Cancer 121:915-920), no correlation with histologic subtype (p=1) has been found. These results were surprising in light of evidence that overexpression of mutated, but not wild-type, PIK3CA leads to transformation, and suggest either of 2 possibilities: 1) PIK3CA suffers from prevalent cryptic mutations, or 2) the effects of wild-type PIK3CA overexpression in human tumors were not captured by the transformation assays.
The suggestion that PI3K activation associates with poor prognosis suggested that measuring PI3K activity might improve prognostication of localized endometrial cancer. Expression of STMN1 has previously been shown to correlate with PI3K activity in breast cancer (Saal LH et al. (2007) Proc Natl Acad Sci USA 104:7564-7569) and can be measured by immunohistochemistry in paraffin-embedded tissue. Herein, STMN1 expression by immunohistochemistry in 72 tumors, including 66 with SNP array and 53 with expression data has been determined. Although STMN1 is not a member of our PI3K activation signature, STMN1 expression correlated with PI3K scores (p=0.05). High STMN1 expression also correlated with PIK3CA amplification (p=0.04) and overexpression (p=0.04), and segregation in Cluster 2 (p=0.03), supporting our prior associations between these features and PI3K pathway activation.
High STMN1 expression was also associated with poor recurrence-free survival in our original tumor set (p=0.006) and with poor recurrence-free (p=0.01) and overall (p=0.01) survival in a validation set of 241 tumors from a population-based series of all endometrial carcinoma in Hordaland County from 1981-1990 (Salvesen HB, Iversen OE & Akslen L A (1999) J Clin Oncol 17:1382-1390; Salvesen HB et al. (2002) Cancer 94:2185-2191). In both tumor sets, STMN1 expression correlated with grade, mitotic rate, presence of necrosis or vascular invasion, and Type II status (Table 5). Nevertheless, across all 313 cases (minus 5 with missing clinical data), high STMN1 expression was an independent prognostic indicator to FIGO stage, histologic subtype, grade, and age (p=0.004; Table 6). In particular, high STMN1 expression was associated with poor prognosis in the otherwise low-risk endometrioid subgroup (p=0.007, data not shown).
Ultimately, the goals of integrated genomic analyses of localized tumors are to enable development of clinical assays to distinguish aggressive tumors requiring therapy beyond resection, and of effective therapeutics for such tumors. It is shown herein that both transcriptional and copy-number profiles of endometrial tumors contain prognostic information that is partly reflected in expression levels of PIK3CA, in vitro PI3K activation signatures, PTEN, and STMN1. Further, it is shown that PTEN and PIK3CA mutations appear to have different transcriptional and phenotypic correlates than changes in expression of these genes. These results suggest that further investigation of the specific consequences of mutation and altered expression is warranted. They also emphasize the potential utility of clinical assays for PI3K pathway activation to identify patients with aggressive disease, and the particular relevance of therapeutics that inhibit this pathway.
Aside from the PI3K pathway, the general survey of chromosomal changes in endometrial carcinoma also identified approximately twenty other regions of significant copy-number change. Most of these copy-number changes involve tens to hundreds of genes, so even in cases where known oncogenes or tumor suppressors are within the regions most affected by these copy-number changes, the genomic data are ambiguous as to the actual target. In many cases, including amplification of 3q26.32, the size of these events may suggest multiple targets. Moreover, functional data tying even known oncogenes and tumor suppressors to carcinogenesis in endometrial cancer model systems are for the most part lacking. The limited number of significant regions of copy-number change suggests that comprehensive, systematic experiments to identify these oncogenes and tumor suppressors in endometrial cancer are feasible. Such experiments point to therapeutic targets for women with all stages of endometrial carcinoma.
Selecting Gene Subsets with Good Combined Discriminatory Power
As our objective function to minimize for determination of maximal predictive power, we used the sum of squared residuals between the relative probability of the correct class label and one given by a diagonal linear discriminant analysis (DLDA) classifier. The relative probability given a DLDA classifier is the probability density for the correct class divided by the sum of probability densities over both classes. Ideally a classifier assigns relative probability 1 to the correct class label in all cases, but this will often not be the case in practice.
We then tested increasing numbers of genes using a forward feature subset selection method (Jonassen, B. T., 2002, Genome Biology, 3:1-0017.11) and found a 29 gene predictor gave the best results. These genes were therefore included in a gene set for validation with QRT-PCR.
5 μm tissue microarray sections of paraffin-embedded tissue were stained, using antigen retrieval for 10 min at 750 W and 15 min at 350 W in Citrate buffer (pH=6). Slides were incubated 1 hour at room temperature with polyoclonal STMN1 antibody (#3352, Cell Signaling) diluted 1:50. A staining index was calculated as the product of staining intensity (0-3) and area of positive tumor cells (1=<10%, 2=10%-50%, 3=>50%). Values in the upper quartile (which corresponded to indices of 6 and 9) were considered positive.
The association between 3q26.32 amplification and tumor recurrence suggests a causal relationship, with its functional effects leading to the aggressive phenotype. An alternative model would be that both 3q26.23 amplification and the aggressive phenotype are caused by a prior event, such as generalized aneuploidy in the cell, leading to an association but no direct causal link. Although this possibility cannot be ruled out, when aneuploidy in 59 of the tumors with SNP array data have been analysed, it has been found that 3q26.32 amplification remains significantly associated with recurrence-free survival after adjustment for the impact of ploidy (p=0.03). It therefore appears that amplification of 3q26.32 has an association with poor survival independent of the overall level of copy-number changes in the cell.
Ploidy was determined from DNA histograms based on measurement of 104-105 cells by flow cytometry, using fresh tumors and adjacent HE sections to confirm malignant histology.
Between 2001 and September 2010, 603 patients treated for endometrial cancer were recruited prospectively in a population based setting, Stathmin expression in primary tumors were measured by immunohistochemistry and linked to treatment response to taxanes in patients with metastatic disease. Response was evaluated by the RECIST criteria and analysed as partial-/complete response versus stable disease/progression.
Of the 603 patients a total of 116 either relapsed (n=79) or progressed (n=37) after their first line of treatment. Of these, 90 were treated with chemotherapy (n=33), radiation (n=38) or hormonal therapy (n=15). The remaining did not receive any further treatment or underwent surgery. Complete information regarding response to therapy according to the RECIST criteria was available in 57 patients. Stathmin expression in primary tumors predicted response to microtubule-stabilising chemotherapy (p=0.02, FE test): Amongst patients with low expression of stathmin 11 of 12 (92%) had partial-/complete response, whereas only 2 of 6 (33%) patients with a high level of stathmin had partial-/complete response (p=0.02, F.E.). Stathmin expression was not associated with response to other treatment modalities.
†Median value in millimeters
‡Median number of mitotic figures per 10 fields at original magnification × 40
§Defined as non-endometrioid, high-grade endometrioid, or lacking both estrogen receptor and progesterone receptor
0-3.5
0-4.2
†Frequency of amplification or deletion to any level. High-level amplifications were seen for LMYC, PIK3CA, EGFR, 6p21.2, EGFR, CCNE1 (1 case each), and MYC (2 cases).
†In years after primary surgery.
‡By log-rank test, after Bonferroni correction for 11 hypotheses.
†Using a Pearson chi-square test when otherwise not specified.
‡Median depth of myometrial infiltration in millimetres, number of mitotic figures per 10 fields at magnification 40×, Mann-Whitney U test.
§Defined as either non-endometrioid, high-grades endometrioid, or lacking both estrogen receptor and progesterone receptor.
¶N = 53 cases, 1-sided test.
†Using a log-ratio test.
‡Continuous variable with hazard ratio given per year.
Filing Document | Filing Date | Country | Kind | 371c Date |
---|---|---|---|---|
PCT/EP2011/072199 | 12/8/2011 | WO | 00 | 6/6/2013 |
Number | Date | Country | |
---|---|---|---|
61310109 | Mar 2010 | US |
Number | Date | Country | |
---|---|---|---|
Parent | 12962946 | Dec 2010 | US |
Child | 13991947 | US |