A SEQUENCE LISTING in computer-readable form (.txt file) accompanies this application having SEQ ID NO:1 through SEQ ID NO:10. The computer-readable form (.txt file) of the SEQUENCE LISTING is incorporated by reference into this application.
This invention relates to a method of providing a treatment to a patient having non-small cell lung cancer comprising extracting total RNA from a formalin fixed and paraffin embedded tumor of non-small cell lung cancer of a patient after the surgical resection, generating complementary DNA (cDNA) of the extracted total RNA from said patient tumor, quantifying of mRNA expression of 7 genes of ABCC4 (SEQ ID NO:1), CCL19 (SEQ ID NO:2), SLC39A8 (SEQ ID NO:3), CD27 (SEQ ID NO:4), FUT7 (SEQ ID NO:5), ZNF71 (SEQ ID NO:6); and DAG1 (SEQ ID NO:7), normalizing of the quantification of said 7 genes with the quantification of a control gene UBC (SEQ ID NO:8) or a housekeeping gene, and utilizing said normalized 7 gene mRNA expression quantification to determine whether said patient will benefit from receiving adjuvant chemotherapy or not. This invention also relates to a method of providing a treatment to a patient having non-small cell lung cancer comprising providing protein expression of SEQ ID NO:9 (protein ZNF71) or SEQ ID NO:10 (protein CD27) and quantifying said protein expression with ELISA, or immunocytochemistry staining or immunohistochemistry staining correlated with either said CD27 mRNA expression or ZNF71 mRNA, respectively, in a patient tumor and a cancer-free tissue adjacent to said tumor, and determining a prognosis of said patient from said protein expression.
Lung cancer is the leading cause of cancer-related deaths in the world, and non-small cell lung cancer (NSCLC) accounts for almost 80% of lung cancer deaths [1]. Major histology of NSCLC includes lung adenocarcinoma and squamous cell lung carcinoma. Surgical resection is the major treatment for early stage NSCLC. However, about 22-38% of stage I NSCLC patients will develop tumor recurrence within five years following the surgery [2]. It is therefore important to select early stage NSCLC patients for more aggressive treatment. While adjuvant chemotherapy of stage II and stage III disease has resulted in 10-15% increased overall survival [3], the prognosis for early stage NSCLC remains poor [4]. Currently, there are no clinically available molecular assays to predict the risk for tumor recurrence and the clinical benefits of chemotherapy in NSCLC patients.
Immunotherapy has rapidly gained attention of oncologists as an effective and less toxic treatment than chemotherapy in patients with advanced lung cancers [5-7]. A recent study used paired single cell analysis to compare normal lung tissue and blood with tumor tissue in stage I NSCLC, and found that early-stage tumors had already begun to alter the immune cells in their microenvironment [8]. These results suggest that immunotherapy could potentially be used to treat early stage lung cancer patients. However, predictive biomarkers of immunotherapy are not well established except PD-1/PD-L1, and it is unlikely that a single marker is sufficient.
High-throughput technologies, such as microarray and RNA-seq, promise the discovery of novel biomarkers from genome-scale studies. The FDA conducted a systematic evaluation and suggested continued usefulness of legacy microarray data and established microarray biomarkers and predictive models in the forthcoming RNA-seq era [9]. However, several disadvantages have limited the application of high-throughput techniques in routine clinical tests, including costs, reproducibility, and data analyses [10]. Compared with microarray/RNA-seq, quantitative real-time RT-PCR (qRT-PCR) is more efficient, consistent, and able to measure gene expression over a greater dynamic range [11]. The combined use of real-time qRT-PCR with high-throughput analysis can overcome the inherent biases of the high-throughput techniques and is emerging as the optimal method of choice to translate genome research into clinical practice [12]. The protein expression validation of the identified MRNA biomarkers could substantiate their ultimate functional involvements in disease, and may lead to the discovery of potential proteomic biomarkers in abundant FFPE samples for broader applications in community hospitals.
DNA microarray-based studies identified gene expression-based NSCLC prognostic [13] and predictive biomarkers [14, 15]. A qRT-PCR based 14-gene assay by Kratz et al [16] is prognostic of non-squamous NSCLC outcome in FFPE tissues and is ready for wide-spread clinical applications. However, this 14-gene assay is limited to non-squamous NSCLC and is not shown to be predictive of the clinical benefits of chemotherapy.
The present invention provides a multi-gene assay predictive of the clinical benefits of chemotherapy in non-small cell lung cancer (NSCLC) patients, and provides for their protein expression as therapeutic targets.
This invention discloses a method using a 7-gene assay ABCC4 (SEQ ID NO:1), CCL19 (SEQ ID NO:2), SLC39A8 (SEQ ID NO:3), CD27 (SEQ ID NO:4), FUT7 (SEQ ID NO:5), ZNF71 (SEQ ID NO:6), and DAG1(SEQ ID NO:7) for selecting adjuvant chemotherapy treatment for a patient with non-small cell lung cancer after their surgery. This treatment method using the 7-gene assay can predict a patient's formalin fixed and paraffin embedded tumor as either with benefit from adjuvant chemotherapy or no benefit from adjuvant chemotherapy after receiving surgery. In the published data of the 7-gene assay, it is shown that for those patients who were predicted as with benefit from chemotherapy, their disease specific-survival was significantly (p<0.05) longer in those who actually received adjuvant chemotherapy compared with those who did not receive adjuvant chemotherapy. In the contrast, for those patients who were predicted with the 7-gene assay as no benefit from adjuvant chemotherapy, their disease-specific survival was actually shorter when they received adjuvant chemotherapy compared with those who did not receive adjuvant chemotherapy, due to unnecessary chemotherapeutic treatment and associated cytotoxicity side-effects (see
Within the 7-gene assay, 4 genes, ABCC4 (SEQ ID NO:1), FUT7 (SEQ ID NO:5), ZNF71 (SEQ ID NO:6), and SLC39A8 (SEQ ID NO:3), each individually predicted chemosensitivity or chemoresistance to specific adjuvant chemotherapy (see Table 2). Specifically, high expression of ABCC4 (SEQ ID NO:1) predicted chemoresistance to Carboplatin and Taxol, Taxol, Carboplatin and Taxotere, Cisplatin and Taxetere, and Cisplatin and Taxol. High expression of FUT7 (SEQ ID NO:5) predicted chemosensitivity to Carboplatin. High expression of ZNF71 (SEQ ID NO:6) predicted chemosentivity to Carboplatin and Taxol, Carboplatin and Taxotere, Cisplatin and Taxotere, and Cisplatin and Taxol. High expression of SLC39A8 (SEQ ID NO:3) predicted chemoresistance to Taxol and Alimta (pemetrexed).
The protein expression of ZNF71 (SEQ ID NO:9) quantified with automated quantitative analysis (AQUA) correlated with its mRNA expression in patient tumors. The protein expression of ZNF71 (SEQ ID NO:9) can independently classify patients into prognosis (longer survival) group or poor prognosis (shorter survival) group (see
The protein expression of CD27 (SEQ ID NO:10) quantified with ELISA had a significant correlation with its mRNA in patient tumors and adjacent normal lung tissues, and could be an independent protein biomarker for patient prognosis and treatment selection (see
An embodiment of this invention provides a method of providing a treatment to a patient having non-small cell lung cancer comprising extracting total RNA from a formalin fixed and paraffin embedded tumor of non-small cell lung cancer of a patient after the surgical resection; generating complementary DNA (cDNA) of the extracted total RNA from said patient tumor; quantifying of mRNA expression of 7 genes of ABCC4 (SEQ ID NO:1), CCL19 (SEQ ID NO:2), SLC39A8 (SEQ ID NO:3), CD27 (SEQ ID NO:4), FUT7 (SEQ ID NO:5), ZNF71 (SEQ ID NO:6); and DAG1 (SEQ ID NO:7); normalizing of the quantification of said 7 genes with the quantification of a control gene UBC (SEQ ID NO:8) or a housekeeping gene; and utilizing said normalized 7 gene mRNA expression quantification to determine whether said patient will benefit from receiving adjuvant chemotherapy or not. In a preferred embodiment of this method, the method further comprises administering to said patient a therapeutically effective amount of one of the following adjuvant chemotherapies (a) Cisplatin and Taxol, (b) Cisplatin and Taxotere, (c) Carboplatin, (d) Carboplatin and Taxol, (e) Carboplatin and Taxotere, (f) Taxol, and (g) Alimta (pemetrexed). In a more preferred embodiment of tis method, this method comprises the quantification of mRNA expression of three genes of ABCC4 (SEQ ID NO:1), CCL19 (SEQ ID NO:2), and SLC39A8 (SEQ ID NO:3). In another preferred embodiment of this invention, this method comprises the quantification of mRNA expression of four genes of CD27 (SEQ ID NO:4), FUT7 (SEQ ID NO:5), ZNF71 (SEQ ID NO:6), and DAG1 (SEQ ID NO:7).
Another embodiment of this invention provides the method, as described above, wherein said quantification of mRNA expression of ABCC4 (SEQ ID NO:1) and utilization of said normalized ABCC4 gene mRNA expression quantification to determine whether said patient will benefit from receiving adjuvant chemotherapy of one of (a) Cisplatin and Taxol, (b) Cisplatin and Taxotere, (d) Carboplatin and Taxol, (e) Carboplatin and Taxotere, and (f) Taxol.
Another embodiment of this invention provides the method, as described above, wherein said quantification of mRNA expression of FUT7 (SEQ ID NO:5) and utilization of said normalized FUT7 gene mRNA expression quantification to determine whether said patient will benefit from receiving adjuvant chemotherapy of Carboplatin.
Another embodiment of this invention provides the method, as described above, wherein said quantification of mRNA expression of ZNF71 (SEQ ID NO:6) and utilization of said normalized ZNF71 gene mRNA expression quantification to determine whether said patient will benefit from receiving adjuvant chemotherapy of (a) Carboplatin and Taxol, (b) Carboplatin and Taxotere, (c) Cisplatin and Taxotere, and (d) Cisplatin and Taxol.
Another embodiment of this invention provides the method, as described above, wherein said quantification of mRNA expression of SLC39A8 (SEQ ID NO:3) and utilization of said normalized SLC39A8 gene mRNA expression quantification to determine whether said patient will benefit from receiving adjuvant chemotherapy of one of (a) Taxol, and (b) Alimta (pemetrexed).
In another embodiment of this invention a method of providing a treatment to a patient having non-small cell lung cancer, is disclosed, comprising providing protein expression of ZNF71 (SEQ ID NO: 9); quantifying said protein expression of said ZNF71 (SEQ ID NO:9) with automated quantitative analysis (AQUA) correlated with said ZNF71 (SEQ ID NO:9) mRNA expression in a patient tumor; and determining a prognosis of said patient from said protein expression of said ZNF71 (SEQ ID NO:9). In another embodiment of this method, the method includes wherein said prognosis of said patient is either longer survival or shorter survival.
In another embodiment of this invention, a method of providing a treatment to a patient having non-small cell lung cancer, is disclosed, comprising providing protein expression of CD27 (SEQ ID NO:10); quantifying said protein expression of said CD27 (SEQ ID NO:10) with ELISA correlated with said CD27 (SEQ ID NO:10) mRNA expression in a patient tumor and a cancer-free tissue adjacent to said tumor; and determining a prognosis of said patient from said protein expression of said CD27 (SEQ ID NO:10). In another embodiment of this method, the method includes wherein said prognosis of said patient is either longer survival or shorter survival. This method, optionally, includes administering to said patient a therapeutically effective amount of an adjuvant chemotherapy.
PCT/US2019/036953 discloses a 7-gene assay using snap-frozen samples to identify patients at risk for tumor recurrence and metastasis and selection of optimal chemotherapeutic regimen in these patients. The following 7 genes were used in the gene assay based on qRT-PCR: ABCC4, CCL19, CD27, DAG1, FUT7, SLC39A8, and ZNF44. In this invention, a new protocol and algorithm for this new 7-gene assay in non-small lung cancer patient formalin fixed samples and paraffin embedded samples is provided. The present method comprises extracting total RNA from patient formalin fixed and patient paraffin embedded samples; quantifying mRNA expression profiles in qRT-PCR in formalin fixed and paraffin embedded samples and then compared with the matched snap-frozen tumor tissues from the same patient; based on the collected data a new algorithm and methodology is developed for prognosis and prediction of chemotherapy benefits in non-small cell lung cancer patients using the patient's formalin fixed and paraffin embedded sample. Protein expression of the identified biomarkers was also evaluated in patient formalin fixed or paraffin embedded tissue samples using immunohistochemistry (IHC). IHC is commonly used in pathology laboratories, and the IHC results of the method of the present invention provide prognosis of non-small lung cancer as independent companion tests. Immunohistochemistry is the most common application of immunostaining. IHC is a known technique used to determine the presence and level of specific cellular proteins. IHC involves the process of selectively identifying antigens (proteins) in cells of a tissue section by exploiting the principle of antibodies binding specifically to antigens in a biological tissue, here non-small cell lung cancer tissue. In general, immunohistochemistry comprises the the following steps: (1) fixation to keep the sample in place, (2) antigen retrieval to increase the availablility of proteinsfor detection, (3) bocking to minimize any background signals, and (4) antibody labeling and visualization. Generally, immunohistochemistry markers are monoclonal antibodies used to identify specific proteins in tissue samples. The antibody binds to the protein and a color reagent stains the protein, if in fact that protein is present in the tissue sample.
Specifically this invention provides ZNF71 protein expression in formalin fixed and paraffin embedded samples was a prognostic biobarker of non-small cell lung cancer using a technique known by those skilled in the art as AQUA. In this method, we use quantification results from IHC tests with new antibodies for ZNF71. This invention provides new protein expression assays for ZNF71 for non-small cell lung cancer prognosis. In addition CD27 was previously reported as a potential protein biomarker based upon snap-frozen samples (PCT/US2019/036953). This invention tests CD27 in formalin fixed and paraffin embedded samples with immunocytochemistry staining. This invention provides a prognostic model for non-small cell lung cancer using patient clinical, pathological, and demographic information to inform optimal treatment options. An online tool, PersonalizeRX (available at www.personalizedrx.org) is already used in clinics worldwide. This method integrates the mRNA 7-gene assay, protein based IHC tests, and the PersonalizeRX tool, into one algorithm to provide healthcare providers (i.e. physicians and clinicians) with an accurate estimate of a non-small cell lung cancer patient clinical outcomes. This method thus provides a tool for establishing precision therapy in non-small cell lung cancer patients.
PCT/US2019/036953 describes a 7-gene assay based on snap-frozen non-small cell lung cancer patient samples. The technology described in PCT/US2019/036953 is not applicable to formalin fixed or paraffin embedded samples that are abundant in the majority of community hospitals. PCT/US2019/036953 describes a protein biomarker ZNF71 based upon AQUA and a now discontinued antibody. The present invention provides a method to quantify ZNF71 with new antibodies using IHC in formalin fixed and paraffin embedded samples of non-small cell lung cancer tumors. Further, this 7-gene assay of the present invention and the IHC assay of ZNF71 is integrated with patient clinical, pathologies, and demographic information into one algorithm for selection of optimal treatment of the non-small cell lung cancer patient. The present invention provides a method that utilizes formalin fixed and paraffin embedded non-small cell lung cancer patient tissue samples for mRNA quantification.
This invention provides a (1) a mRNA based 7-gene assay and algorithm, (2) an IHC based ZNF71 and CD27 assays and algorithm, and (3) an integrated mRNA 7-gene assay, ZNF71 and CD27 IHC assays, and patient clinical information in one algorithm. Materials:
122 non-small cell lung cancer patient samples were obtained form Case Western Reserve University, 101 lung adenocarcinoma tumor specimens from University of Michigan Comprehensive Cancer Center, 65 non-small cell lung cancer tumor specimens from NorthShore University Health System Kellogg Center Cancer Center, and 49 specinens from West Virginia University Cancer Institute/Mary Babb Randolph Cancer Center.
Alexa 546-conjugated goat anti-mouse seconday antibody (Life Technologies), Cy5-Tyramide (Perkin Elmer).
Human CD27/TNFSRF7 DuoSet ELISA kit (contained antibodies).
ZNF71 antibodies from GeneTex and Sigma.
Data from JBR.10 clinical trial were obtained from the NCBI Gene Expression Omnibus website.
An embodiment of this invention provides a method of providing a treatment to a patient having non-small cell lung cancer comprising extracting total RNA from a formalin fixed or paraffin embedded tumor of non-small cell lung cancer of a patient after the surgical resection; generating complementary DNA (cDNA) of the extracted total RNA from said patient tumor; quantifying of mRNA expression of 7 genes of ABCC4 (SEQ ID NO:1), CCL19 (SEQ ID NO:2), SLC39A8 (SEQ ID NO:3), CD27 (SEQ ID NO:4), FUT7 (SEQ ID NO:5), ZNF71 (SEQ ID NO:6); and DAG1 (SEQ ID NO:7); normalizing of the quantification of said 7 genes with the quantification of a control gene UBC (SEQ ID NO:8) or a housekeeping gene; and utilizing said normalized 7 gene mRNA expression quantification to determine whether said patient will benefit from receiving adjuvant chemotherapy or not.
As used herein, a housekeeping gene is a typically constitutive genes that is required for the maintenance of basal cellular functions that are essential for the existence of a cell, regardless of its specific role in the tissue or organism. Thus, they are generally expressed in all cells of an organism under normal and patho-physiological conditions, irrespective of tissue type, developmental stage, cell cycle state, or external signal. For example, housekeeping genes are widely used as internal controls for experimental studies. The reliability of any relative RT-PCR experiment can be improved by including an invariant endogenous control (reference gene) in the assay to correct for sample to sample variations in RT-PCR efficiency and errors in sample quantification. A biologically meaningful reporting of target mRNA copy numbers requires accurate and relevant normalization to some standard and is recommended in quantitative RT-PCR. Many housekeeping genes are known to those persons skilled in the art, such as for example, but not limited to, 18S ribosomal RNA (RRN18S), RNA polymerase 2 subunit A (PolR2A), glyceraldehyde phosphate dehydrogenase (GAPDH), or β 2-microglobulin (B2M).
In a preferred embodiment of this method, the method further comprises administering to said patient a therapeutically effective amount of one of the following adjuvant chemotherapies (a) cisplatin and Taxol (paclitaxel), (b) cisplatin and Taxotere (docetaxel), (c) carboplatin, (d) carboplatin and Taxol, (e) carboplatin and Taxotere, (f) Taxol, and (g) Alimta (pemetrexed). Taxol is a registered trademark owned by Bristol-Myers Squibb Company, New York, N.Y., USA; Taxotere is a registered trademark owned by Aventis Pharma S.A., Cedex, France; and Alimnta is a registered trademark owned by Eli Lilly and Company, Indianapolis, Ind., USA. In a more preferred embodiment of this method, this method comprises the quantification of mRNA expression of three genes of ABCC4 (SEQ ID NO:1), CCL19 (SEQ ID NO:2), and SLC39A8 (SEQ ID NO:3). In another more preferred embodiment of this method, this method comprises the quantification of mRNA expression of four genes of CD27 (SEQ ID NO:4), FUT7 (SEQ ID NO:5), ZNF71 (SEQ ID NO:6), and DAG1 (SEQ ID NO:7).
Another embodiment of this invention provides the method, as described above, wherein said quantification of mRNA expression of ABCC4 (SEQ ID NO:1) and utilization of said normalized ABCC4 (SEQ ID NO:1) gene mRNA expression quantification to determine whether said patient will benefit from receiving adjuvant chemotherapy of one of (a) cisplatin and Taxol, (b) cisplatin and Taxotere, (d) carboplatin and Taxol, (e) carboplatin and Taxotere, and (f) Taxol.
Another embodiment of this invention provides the method, as described above, wherein said quantification of mRNA expression of FUT7 (SEQ ID NO:5) and utilization of said normalized FUT7 (SEQ ID NO:5) gene mRNA expression quantification to determine whether said patient will benefit from receiving adjuvant chemotherapy of Carboplatin.
Another embodiment of this invention provides the method, as described above, wherein said quantification of mRNA expression of ZNF71 (SEQ ID NO:6) and utilization of said normalized ZNF71 (SEQ ID NO:6) gene mRNA expression quantification to determine whether said patient will benefit from receiving adjuvant chemotherapy of (a) carboplatin and Taxol, (b) carboplatin and Taxotere, (c) cisplatin and Taxotere, and (d) cisplatin and Taxol.
Another embodiment of this invention provides the method, as described above, wherein said quantification of mRNA expression of SLC39A8 (SEQ ID NO:3) and utilization of said normalized SLC39A8 (SEQ ID NO:3) gene mRNA expression quantification to determine whether said patient will benefit from receiving adjuvant chemotherapy of one of (a) Taxol, and (b) Alimta (pemetrexed).
In another embodiment of this invention a method of providing a treatment to a patient having non-small cell lung cancer, is disclosed, comprising providing protein expression of ZNF71 (SEQ ID NO: 9); quantifying said protein expression of said ZNF71 (SEQ ID NO:9) with automated quantitative analysis (AQUA) correlated with said ZNF71 mRNA expression in a patient tumor; and determining a prognosis of said patient from said protein expression of said ZNF71 (SEQ ID NO:9). In another embodiment of this method, the method includes wherein said prognosis of said patient is either longer survival or shorter survival.
In another embodiment of this invention, a method of providing a treatment to a patient having non-small cell lung cancer, is disclosed, comprising providing protein expression of CD27 (SEQ ID NO:10); quantifying said protein expression of said CD27 (SEQ ID NO:10) with ELISA correlated with said CD27 mRNA expression in a patient tumor and a cancer-free tissue adjacent to said tumor; and determining a prognosis of said patient from said protein expression of said CD27 (SEQ ID NO:10). In another embodiment of this method, the method includes wherein said prognosis of said patient is either longer survival or shorter survival. This method further includes administering to said patient a therapeutically effective amount of an adjuvant chemotherapy.
Patients and methods: The mRNA expression of 160 genes identified from microarray was analyzed in qRT-PCR assays of independent formalin fixed and paraffin embedded (“FFPE”) non-small cell lung cancer (NSCLC) tumors to develop a predictive signature. A clinical trial JBR.10 was included in the validation. Hazard ratio was used to select genes, and decision-trees were used to construct the predictive model. Protein expression was quantified with AQUA in 500 FFPE NSCLC samples. Results: A 7-gene signature (of this invention) was identified from training cohort (n=83) with accurate patient stratification (P=0.0043) and was validated in independent patient cohorts (n=248, P<0.0001) in Kaplan-Meier analyses. In the predicted benefit group, there was a significantly better disease-specific survival in patients receiving adjuvant chemotherapy in both training (P=0.035) and validation (P=0.0049) sets. In the predicted non-benefit group, there was no survival benefit in patients receiving chemotherapy in either set. The protein expression of ZNF71 (SEQ ID NO:9) quantified with AQUA scores produced robust patient stratification in separate training (P=0.021) and validation (P=0.047)NSCLC cohorts. The protein expression of CD27 (SEQ ID NO:10) quantified with ELISA had a strong correlation with its mRNA expression in NSCLC tumors (Spearman coefficient-0.494, P<0.0088). Multiple signature genes had concordant DNA copy number variation, mRNA and protein expression in NSCLC progression.
Those persons skilled in the art will understand that this invention presents a predictive multi-gene assay and prognostic protein biomarkers clinically applicable for improving NSCLC treatment in patients suing formalin fixed or paraffin embedded tumor samples, with important implications in lung cancer chemotherapy/immunotherapy.
In this invention, a combined analysis of genome-wide transcriptional profiles and qRT-PCR was utilized to develop a multi-gene assay both prognostic of NSCLC outcome and predictive of the benefits of chemotherapy. Patient cohorts from multiple hospitals in the US and JBR.10 data [14] were used to validate this multi-gene assay. Protein expression of the identified biomarkers (using immunohistochemistrywas also evaluated in patient tissue samples and correlated with the mRNA expression and DNA copy number variation to substantiate their functional involvement and potential as therapeutic targets in chemotherapy and immunotherapy, in addition to companion tests.
Materials and Methods
Patient samples. Clinical characteristics of patient cohorts used in qRT-PCR assays are summarized in Table 1. All NSCLC patients were staged I, II, or IIIA at the time of diagnosis. Tumor tissues were collected in surgical resections and were then formalin fixed or paraffin embedded until used for RNA extraction. Tumor cell content was above 50% for qRT-PCR assays. Those with missing AJCC staging information, missing histology, death within 30 days of resection or from other disease conditions were excluded from further analysis. A total of 122 NSCLC patient samples were obtained from Case Western Reserve University (CWRU) Comprehensive Cancer Center. Total RNA of good quality was extracted from 89 tumor specimens. Good quality RNA from 101 lung adenocarcinoma tumor specimens was obtained from University of Michigan (UM) Comprehensive Cancer Center, with detailed description of patients, tissue specimens and mRNA quality check provided in [17]. A total of 65 NSCLC tumor specimens from NorthShore University HealthSystem Kellogg Cancer Center and 49 specimens from West Virginia University Cancer Institute [Mary Babb Randolph Cancer Center (MBRCC)] generated good quality mRNA. The tissue collection in this study was approved by an Institutional Review Board (IRB) at each institution.
RNA extraction, and quality and concentration assessments. Total RNA was extracted from formalin fixed or paraffin embedded tumor tissues using a RNeasy mini kit according the manufacturer's protocol (Qiagen, USA), followed by elution in 30 μl of RNase-free water and storage at −80° C. The quality and integrity of the RNA, the 28S to 18S ratio, and a visual image of the 28S and 18S bands were evaluated on the 2100 Bioanalyzer (Agilent Technologies, CA). RNA assessed as having good quality from 304 tumor samples was included for further analysis. The RNA concentration of each sample was assessed using a Nanodrop-1000 Spectrophotometer (NanoDrop Tech, Germany).
Generation of complementary DNA (cDNA). The reverse transcriptase polymerase chain reaction was used to convert the high-quality single-stranded RNA samples to double-stranded cDNA, using an Applied Biosystems GeneAmp® PCR 9700 machine (Foster City, Calif.). For standardization across all samples, one microgram of RNA was used to generate cDNA.
Real-time quantitative RT-PCR low-density arrays. Real-time qRT-PCR assays of independent patient cohorts of NSCLC tumor samples were used to further select biomarkers to form a multi-gene assay from prognostic genes identified from microarray data in our previous studies [18-21]. The identified prognostic genes were initially validated with multiple independent NSCLC microarray data publically available [18-21]. Based on the validation results, 160 prognostic genes and three housekeeping genes were included in the qRT-PCR experiments. The three housekeeping genes were 18S, UBC, and POLR2A due to their confirmed constant mRNA expressions across samples [18].
Three hundred thirty seven (337) tumor samples were analyzed with good RNA quality using TaqMan microfluidic low-density array (LDA) plates on an ABI 7900HT Fast RT-PCR instrument (Applied Biosystems). Total RNA samples were analyzed on an Agilent 2100 Bioanalyzer RNA 6000 Nano LabChip. The report was generated by the SDS2.3 software (Applied Biosystems). In the report, the number of cycles required to reach threshold fluorescence (Ct) and ΔCT for each sample relative to the control gene defines the expression pattern for a gene. The gene expression data were further analyzed using the 2ΔΔC
Statistical and computational analysis. Prognostic biomarkers were evaluated with Cox proportional hazard model. Hazard ratio was used in the evaluation of prognostic performance of biomarkers. If a biomarker gives a hazard ratio greater than 1, it means that patient samples predicted as high risk are more likely to have a poor outcome. In the evaluation of genes in qRT-PCR assays, ΔCT was used as a covariate in Cox model. If a gene as a hazard ratio greater than 1, it means that down-regulation of this gene is associated with a poor outcome and up-regulation of this gene is associated with a good outcome in NSCLC patients; otherwise, if a gene has a hazard ratio less than 1, it means that down-regulation of this gene is associated with a good outcome and up-regulation of this gene is associated with a poor outcome in NSCLC patients. During the evaluation, UBC (Hs00824723_ml) was chosen as the house keeping gene to normalize gene expression. The CWRU cohort was used as the training set, and seven genes were selected to form a prognostic classifier based on decision trees. These seven genes are ABCC4 (Hs00988717_ml) (SEQ ID NO:1), CCL19 (Hs00171149_ml) (SEQ ID NO:2), SLC39A8 (Hs00223357_ml) (SEQ IDNO:3), CD27 (Hs00154297_ml) (SEQ ID NO:4), FUT7 (Hs00237083_ml) (SEQ ID NO:5), ZNF71 (Hs00221893_ml) (SEQ ID NO:6), and DAG1 (Hs00189308_ml) (SEQ ID NO:7). The 7-gene prognostic method of this invention was validated with independent patient cohorts (UM, MBRCC, and NorthShore). In Kaplan-Meier analysis, log-rank tests or Wilcoxon tests were used to assess the difference in probability of survival of different prognostic groups. All the analyses were performed with packages in R or SAS unless otherwise specified.
Validation on clinical trial JBR.10. Data from JBR.10 was obtained from NCBI Gene Expression Omnibus with accession number GSE14814. A total of 133 non-small cell lung cancer samples were profiled for gene expression using Affymetrix 133A platform [14]. Patients were all in early stage (I or II). Patient samples assayed in the same batch with consecutive accession numbers ranging from GSM370913 to GSM371002 (n=90) were used in the validation of the 7-gene signature. Among these patient samples, those who died from other disease conditions were excluded from further analysis. ABCC4 (203196_at), CCL19 (210072_at), CD27 (206150_at), DAG1 (205417_s_at and 212128_s_at), FUT7 (210506_at and 217696_at), SLC39A8 (209266_s_at, 209267_s_at, 216504_s_at, and 219869_s_at), and ZNF444 (218707_at and 50376_at) were used in validating the qRT-PCR based multi-gene assay. For a gene with multiple probe sets, the one with the highest expression value (yielding the clearest signal) in each sample was chosen to represent the gene expression. ZNF71 was not available in the GSE14814 dataset. ZNF444 was chosen to replace ZNF71 to validate the qRT-PCR results, because both ZNF444 and ZNF71 are at locus NC 000019.10 in Chromosome 19 and belong to zinc finger protein family. To be compatible with the ΔCt values in qRT-PCR data, log2 transformed microarray data was used in the analysis, and the expression values of UBC minus those of selected probes were used in the normalization of the microarray data.
Tissue Microarrays (TMA). Samples from 2 retrospective collections of lung cancer were examined in TMA format from Yale University Pathology Archives; Cohort A (YTMA 250 [n=298]) and Cohort B (YTMA 79 [n=202]). TMAs consisted of 0.6 mm cores in 1 (Cohort A) and 2 fold (Cohort B) redundancy. TMAs were prepared according to standard methods. Cohort A comprises 314 serially collected NSCLC who underwent surgical resection of their primary tumor between 2004 and 2011. Cohort B comprises of 202 serially collected NSCLC patients who underwent surgical resection of their primary tumor between 1988 and 2003. All tissue was used after approval from the Yale Human Investigation Committee protocol #9505008219, which approved the patient consent forms or in some cases waiver of consent. The actual number of samples analyzed for each study is lower, due to unavoidable loss of tissue or the absence or limited tumor cells in some spots as is commonly seen in TMA studies. NSCLC patients in stage I, II, and IIIA were included in the analysis. Those who died with no evidence of disease were excluded from further analysis.
Quantitative immunofluorescence. FFPE whole-tissue sections, tissue microarrays (TMAs) and cell pellets were processed at Yale Cancer Center/Pathology Tissue Microarray Facility.
FFPE whole-tissue sections, tissue microarrays (TMAs) and cell pellets were processed as follows: briefly, sections were baked for 30 minutes at 60 degrees Centigrade and underwent two twenty minute wash cysles in xylenes. Slides were rehydrated in two 1-minute washes in 100% ethanol followed by one washing 70% ethanol and finally rinsed in streaming tap water for 5 minutes. Antigen retrieval was performed in sodium citrate buffer pH.6, for 20 minutes at 97 degrees Centigrade in a PT module (LabVision). Endogenous peroxidases were blocked by 30-minute incubation in 2.5% hydrogen peroxide in methanol. Nonspecific antigens were blocked by a 30 minute incubation in 0.3% BSA in TBST. Slides were then incubated with the target primary antibody (ZNF71 Abcam; ab87250), as well as pan cytokeratin (AE1/AE3) overnight at 4 degrees Centigrade diluted at 1:100 to define the tumor compartment.
Primary antibodies were followed by incubation with Alexa 546-conjugated goat anti-mouse secondary antibody (Life Technologies) diluted 1:100 in rabbit EnVision reagent (Dako) for 1 hour. ZNF71 signal was amplified with Cy5-Tyramide (Perkin Elmer) for 10 minutes, and then nuclei were stained with 0.05 mg DAPI in BSA-tween for 10 minutes. Slides were mounted with ProlongGold (Life Technologies). Two TBS-T and one TBS wash was performed between each step after the primary antibody.
Immunofluorescence was quantified using automated quantitative analysis (AQUA) Fluorescent images of DAPI, Cy3 (Alexa 546-cytokeratin), and Cy5 (ZNF71) for each TMA spot were collected. Image analysis was carried out using the AQUAnalysis software (Navigate Biopharma Inc.), which generated an AQUA score for each compartment by dividing the sum of target pixel intensities by the area of the compartment in which the target is measured. AQUA scores were normalized to the exposure time and bit depth at which the images were captured, allowing scores collected at different exposure times to be directly comparable. Specimens with less than 5% tumor area per region of interest were not included in AQUA analysis for not being representative of the corresponding tumor specimen.
Enzyme-Linked Immunosorbent Assay (ELISA). A total of 38 NSCLC patient tissue samples were selected for ELISA assays, including 29 tumor resections of lung adenocarcinoma and squamous cell lung cancer and 9 matched adjacent normal lung tissue samples. The DuoSet ELISA Development Systems from R&D Systems (Minneapolis, Minn.; catalog number: DY382-05) were used for quantifying protein expression of T-Cell Activation Antigen CD27 (CD27)/Tumor Necrosis Factor Receptor Superfamily, Member 7 (TNFRSF7) in NSCLC patient tissue samples, according to manufacturer's protocol. The ELISA assay results were quantified using the Synergy H1 Hybrid Multi-Mode Microplate Readers from BioTek Instruments, Inc. (Winooski, Vt.). Samples that yielded a positive OD values were included for further analysis. Statistical analysis was done using a two-sample t-test assuming unequal variances. The concordance between CD27 mRNA and protein expression was evaluated with Spearman correlation coefficient.
The NSCLC prognostic biomarkers identified with hybrid feature selection models [18, 19] and molecular network approach [20, 21] in our previous studies were validated with multiple independent microarray datasets. Based on the validation results in microarray data, 160 genes were selected for assays using low-density microfluidic qRT-PCR arrays. Among 160 genes analyzed in the qRT-PCR assays, a 7-gene signature of this invention was identified from training cohort obtained from Case Western Reserve University (CWRU; n=83). Details of the decision tree based 7-gene prognostic and predictive method of this invention are provided in
The chemoresponse prediction for specific therapeutic agents was examined in the identified 7 biomarkers. In particular, gene expression of ATP binding cassette subfamily C member 4 (ABCC4) was predictive of chemoresistance in patients receiving carboplatin, cisplatin, and Taxol, with under-expressed mRNA (higher ΔCt) value associated with significantly decreased hazard ratio of death from disease and tumor recurrence (see Table 2). In patients treated with carboplatin plus Taxol, using ΔCt value of ABCC4 in Cox model, the hazard ratio of death from disease of was 0.43 (95% CI: [0.208, 0.888], P=0.02) and the hazard ratio of recurrence was 0.343 (95% CI: [0.122, 0.968], P=0.04), both statistically significant. In patients treated with Taxol, the hazard ratio of death from disease of ABCC4 ΔCt value was 0.403 (95% CI: [0.194, 0.834], P=0.01, Cox model) and the hazard ratio of recurrence was 0.48 (95% CI: [0.253, 0.912], P =0.02, Cox model), both statistically significant. In patients treated with either carboplatin plus Taxol, carboplatin plus Taxotere, cisplatin plus Taxotere, or cisplatin plus Taxol, the hazard ratio of death from disease of ABCC4 ΔCt values was borderline significant (hazard ratio: 0.528 [0.271, 1.028], P=0.06, Cox model) and the hazard ratio of recurrence was significant at 0.545 (95% CI: [0.298, 0.998], P=0.049, Cox model; Table 2). The expression of fucosyltransferase 7 (FUT7) was predictive of chemosensitivity to carboplatin, with under-expressed mRNA (higher ΔCt value) associated with significantly increased hazard ratio of death from disease (hazard ratio: 1.605 [1.058, 2.435], P=0.026, Cox model; Table 2). The expression of zinc finger protein 71 (ZNF7 I)(SEQ ID NO:9) was also predictive of chemosensitivity in patients treated with either carboplatin plus Taxol, carboplatin plus Taxotere, cisplatin plus Taxotere, or cisplatin plus Taxol, with a significant hazard ratio of death from disease 1.986 (95% CI: [1.001, 3.938], P=0.049, Cox model; Table 2). Solute carrier family 39 member 8 (SLC39A8) was predictive of chemoresistance to Taxol, with a borderline significant hazard ratio of recurrence 0.584 (95% CI: [0.33, 1.03], P=0.06, Cox model; Table 2). The expression of SLC39A8 was also predictive of chemoresistance to Alimta (pemetrexed), with a borderline significant hazard ratio of recurrence 0.49 (95% CI: [0.219, 1.098], P=0.08, Cox model; Table 2).
The 7-gene NSCLC prognostic and predictive signature is involved in cell to cell signaling and interaction, inflammatory response, and cellular movement in Ingenuity Pathway Analysis (Qiagen, Redwood City, Calif.). Based on the molecular network of the 7 NSCLC biomarkers (
Protein expression of ZNF71 (SEQ ID NO:9) is prognostic of NSCLC outcome.
To substantiate the functional involvement of the identified 7 signature genes of the methods of this invention, protein expression of these biomarkers was evaluated with immunohistochemistry (IHC). Based on the IHC results, biomarkers with staining of good quality in FFPE NSCLC tumor tissues were further quantified with AQUA. Protein expression of ZNF71 (SEQ ID NO:9) was identified as prognostic of NSCLC outcome in two TMA cohorts (
Concordant mRNA and protein under-expression in NSCLC progression:
The protein expression level of CD27 (SEQ ID NO:10) was quantified with ELISA assays in FFPE NSCLC tumor tissues (n=29) and normal adjacent lung tissues (n=9). Spearman correlation coefficient between mRNA and protein expression of CD27 (SEQ ID NO:10) is 0.494 (P<0.0088;
Lung cancer is the second most common cancer in both men and women, and remains the highest cancer-related mortality with a death rate higher than colon, prostate, and breast cancer combined. Currently, there is no clinically available multi-gene assay to prognosticate and predict the benefits of chemotherapy in formalin fixed or paraffin embedded NSCLC tissues of NSCLC patients for improved personalized treatment. Immunotherapy is more effective and less toxic than chemotherapy in advanced lung cancers [5-8, 29, 30], and recent studies show promise of immunotherapy in early stage lung cancer patients [8]. Nevertheless, predictive biomarkers and therapeutic targets of immunotherapy are not well established.
There were abundant publically available microarray data generated in NSCLC patient tissues. Although microarray platforms are phasing out, the legacy data and biomarkers identified in microarray platforms are still useful in the RNA-seq era [9]. However, high-throughput platforms such as microarrays and RNA-seq are not suitable for routine clinical tests. Validation of biomarkers identified from high-throughput technologies with qRT-PCR emerges as the most promising experimental protocol for developing multi-gene assays for clinical applications.
NSCLC prognostic biomarkers were identified with hybrid feature selection models [18, 19, 31] and molecular network approach [20, 21] in our previous studies. The hybrid feature selection models [18, 19, 31] contain multiple layers of gene selection algorithms in the process of biomarker identification. This scheme takes advantage of different algorithms in different stages of gene shaving, in order to identify the gene signatures with the optimal performance. The molecular network approach [20, 21] constructs genome-scale co-expression networks in good-prognosis and poor-prognosis patient groups separately, and compares the network structures of these two patient groups to identify disease-specific network modules. Next, genes with concurrent co-expression with multiple major lung cancer signaling hallmarks were pinpointed from disease-specific network modules for further gene signature identification. This approach embedded biological relevance into biomarker identification. The signature genes identified with these sophisticated approaches were validated with multiple independent publically available microarray datasets. Genes with consistent expression patterns in multiple validation sets were included in qRT-PCR assays. The 7-gene signature of the methods of this invention identified in qRT-PCR assays was prognostic and predictive of chemoresponse in patient cohorts from multiple hospitals and JBR.10.
The identified 7 signature genes have interactions with major inflammatory and cancer signaling hallmarks including TNF, PI3K, NF-κB, and TGF-β (
The 7-gene signature identified in the methods of this invention does not overlap with the NSCLC gene signatures reported in recent studies [15, 16, 23-25]. However, several biomarker genes identified in this study belong to the same families or functional categories as the biomarkers identified in [14-16]. In particular, FUT7 from the current study and FUT3 from Kratz et al [16] are both fucosyltransferase and involved in metabolism. In the 12-gene prognostic and predictive signature from Tang et al [15], two genes belong to the same family or share similar functions as the 7-gene signature. Specifically, SLC35A5 from Tang et al [15] and SLC39A8 from this study both belong to solute carrier superfamily, and ATPase Phospholipid Transporting 8A1 (ATP8A1) from Tang et al [15] and ATP Binding Cassette Subfamily C Member 4 (ABCC4) from this study are both involved in energy metabolism. The 15-gene prognostic and predictive gene signature of JBR.10 [14] also contains two genes that share similar functions as the 7-gene signature. ATPase Na+/K+ Transporting Subunit Beta 1 (ATP8A1) from Zhu et al [14] and ABCC4 from this study are again involved in energy metabolism, and ZNF236 from Zhu et al [14] and ZNF71 identified in this study both belong to zinc finger protein family. Overall, the 7-gene signature presented in this invention and two previous gene signatures from Zhu et al [14] and Tang et al [15] are all prognostic of NSCLC outcome and predictive of the benefits of chemotherapy. These three gene signatures all contain a biomarker related to ATP activities and energy metabolism. Other shared gene families between the 7-gene signature of this invention and these two signatures include zinc finger protein and solute carrier superfamily. The 7-gene signature and the practical prognostic gene assay for non-squamous NSCLC by Kratz et al [16] both contain biomarkers from fucosyltransferase family. These common gene families shared by the NSCLC gene signatures with promise for clinical utility might be functionally involved in tumor metastasis with implications in lung cancer therapy.
The protein expression of the identified 7 signature genes was also validated in this study. In particular, ZNF71 protein expression quantified with AQUA was a prognostic biomarker in two NSCLC patient cohorts (n=191). Higher mRNA and protein expression of ZNF71 is both associated with good prognosis, and ZNF71 mRNA is predictive of chemosensitivity in Taxol (paclitaxel) plus platinum-based treatment in NSCLC patients, and docetaxel plus platinum-based treatment in NSCLC patients. These results demonstrate that ZNF71 mRNA and protein expression can both be used in prognostication of NSCLC in clinical applications and ZNF71 may be a therapeutic target. CD27 had highly correlated mRNA and protein expression, with significant under-expression in poor prognostic (high-risk) NSCLC patients. CD27 mRNA and protein expression could potentially be used as a biomarker and target in lung cancer immunotherapy. Protein expression of CCL19 was also confirmed with ELISA in NSCLC tumor and adjacent normal tissues. CCL19 protein was under-expressed in FFPE NSCLC tumor tissues compared with normal lung tissues, with no statistically significant difference (results not shown). CCL19 also had lower protein expression in poor-prognosis (high-risk) NSCLC patients compared with good-prognosis (low-risk) patients, with no statistically significant difference (results not shown). The trend of CCL19 protein expression was qualitatively concordant with its mRNA expression that higher expression of CCL19 is associated with good prognostic outcome of NSCLC. CCL19 had a 12.5% (34/271) of a loss of DNA copy number in the NSCLC patient cohort from Starczynowski et al [26] (n=271;
This invention provides a method of measuring the expression gene expression levels comprising determining the level of expression of the following multi-gene set consisting of ABCC4 (SEQ ID NO:1), CCL19 (SEQ ID NO:2), SLC39A8 (SEQ ID NO:3), CD27 (SEQ ID NO:4), FUT7 (SEQ ID NO:5), ZNF71 (SEQ ID NO:6), and DAG1 (SEQ ID NO:7). This method using this particular seven gene combination has never before been known to aid in the benefit of survival rates of patients afflicted with non-small cell lung cancer.
The method comprises the following steps: (1) extraction of total RNA from a formalin fixed or paraffin embedded tumor of non-small cell lung cancer after the surgical resection, (2) generation of complementary DNA (cDNA) of the extracted total RNA from a patient tumor, (3) quantification of mRNA expression of 7 genes: ABCC4 (SEQ ID NO:1), CCL19 (SEQ ID NO:2), SLC39A8 (SEQ ID NO:3) CD27 (SEQ ID NO:4), FUT7 (SEQ ID NO:5), ZNF71 (SEQ ID NO:6), and DAG1 (SEQ ID NO:7), (4) normalization of the quantification of the 7 genes with the quantification of a control gene UBC (SEQ ID NO:8), and (5) utilization of the normalized 7 gene mRNA expression quantification to predict whether a patient will benefit from receiving adjuvant chemotherapy or not. This method further comprises the step of predicting clinical benefit (i.e. prolonged disease-specific survival) of receiving adjuvant chemotherapy, including therapies selected from cisplatin and Taxol (paclitaxel), cisplatin and Taxotere (docetaxel), carboplatin, carboplatin and Taxol (paclitaxel), carboplatin and Taxotere (docetaxel), Taxol (paclitaxel), and Alimta (pemetrexed).
A preferred embodiment of this method includes use of a composition of only the following three: ABCC4 (SEQ ID NO:1), CCL19 (SEQ ID NO:2), and SLC39A8 (SEQ ID NO:3), within the 7-gene assays from this method, which also predicts the clinical benefit of receiving adjuvant chemotherapy. In another preferred embodiment of this method, the method includes use of a composition of only the following four genes: CD27 (SEQ ID NO:4), FUT7 (SEQ ID NO:5), ZNF71 (SEQ ID NO:6) and DAG1 (SEQ ID NO:7), within the 7-gene assays from the method, which also predicts the clinical benefit of receiving adjuvant chemotherapy.
Another method of this invention provides for the high expression of ABCC4 (SEQ ID NO:1) predicted chemoresistance to carboplatin and Taxol (paclitaxel), Taxol (paclitaxel), carboplatin and Taxotere (docetaxel), cisplatin and Taxotere (docetaxel), and cisplatin and Taxol (paclitaxel).
Another method of this invention provides for the high expression of FUT7 (SEQ ID NO:5) predicted chemosensitivity to carboplatin.
Another method of this invention provides for the high expression of ZNF71 (SEQ ID NO:6) predicted chemosentivity to carboplatin and Taxol (paclitaxel), carboplatin and Taxotere (docetaxel), cisplatin and Taxotere (docetaxol), and cisplatin and Taxol (paclitaxel).
Another method of this invention provides for the high expression of SLC39A8 (SEQ ID NO:3) predicted chemoresistance to Taxol (paclitaxel), and Alimta (pemetrexed).
Another method of this invention provides for the protein expression of ZNF71 (SEQ ID NO:9) quantified with automated quantitative analysis (AQUA) correlated with its mRNA expression in patient tumors. The protein expression of ZNF71 (SEQ ID NO:9) can independently classify patients into prognosis (longer survival) group or poor prognosis (shorter survival) group.
Another method of this invention provides for the protein expression of CD27 (SEQ ID NO:10) quantified with ELISA had a significant correlation with its mRNA in patient tumors and adjacent normal lung tissues, and could be an independent protein biomarker for patient prognosis and treatment selection.
#Hazard ratio borderline significant at p <0.08
This invention presents a method using a 7-gene predictive assay based on qRT-PCR to improve NSCLC treatment in clinics using formalin fixed or paraffin embedded samples. This method using a 7-gene assay provides accurate prognostication and prediction of the clinical benefits of chemotherapy in multiple patient cohorts from the US hospitals and the clinical trial JBR.10. The 7-gene assay is enriched in inflammatory response. The protein expression of ZNF71 (SEQ ID NO:9) is prognostic of NSCLC outcome in two independent patient cohorts, which is concordant with its mRNA expression. These results demonstrate that ZNF71 (SEQ ID NO:9) is a prognostic protein biomarker and a useful therapeutic target of NSCLC. The protein expression of CD27 (SEQ ID NO:10) was strongly correlated with its mRNA expression in NSCLC tumor tissues, and serves as a biomarker and target of immunotherapy in lung cancer. Multiple signature genes had concordant DNA copy number variation, mRNA and protein expression in NSCLC progression. The results presented in this invention are important for precision therapy in NSCLC patients, and further provides implications in developing new therapeutic strategies to combat this deadly disease.
This invention provides a method of treating a patient using a 7-gene assay that is predictive of clinical benefits of a patient receiving Alimta (pemetrexed for injection) and commercially available from Eli Lilly and Company, Indianapolis, Ind., USA. Alimta® product is a chemotherapy for the treatment of advanced nonsquamous non-small cell lung cancer (NSCLC). Alimta® is a registered trademark owned or licensed by Eli Lilly and Company.
This invention provides for the protein expression of ZNF71 (SEQ ID NO:9) that is a prognostic marker of non-small cell lung cancer. This invention provides a method of using the expression of ZNF71(SEQ ID NO:9) quantified with AQUA (i.e. Automated Quantitative Analysis ((AQUA)) of In Situ Protein Expression, to identify which patients having non-small cell lung cancer are likely to have good prognosis, and which patients are likely to be poor prognosis.
This invention provides an aid to help physicians determine which non-small cell lung cancer patients, who were initially treated with surgery, will benefit from chemotherapy or immunotherapy. The seven gene assay of the methods of this invention is an aid to predict which patients would benefit from chemotherapty and had significantly prolonged survival time compared to those patients who did not receive any chemotherapy, and which patients would not benefit from chemotherapy and whose long-term post surgical survival time was shorter compared to patients who also had surgery but did not receive any chemotherapy.
This invention provides a method for treating a patient having NSCLC comprising identifying two genes, CD27 (SEQ ID NO:4) and ZNF71 (SEQ ID NO:6), as useful in predicting patient outcomes and developing therapeutic targets in non-small cell lung cancer treatment.
It will be understood by those persons skilled in the art that this invention provides a multi-gene combination assay that provides guidance on the clinical benefits of providing chemotherapy to an individual having non-small cell lung cancer. This invention provides a method for providing precision medicine for lung cancer patients and provides therapeutic targets in both chemotherapy and immunotherapy.
This invention provides a method for improving personalized treatment of individuals having non-small cell lung cancer. Specifically, this invention provides a RT-PCR based method using a 7 gene assay for providing clinical benefits of chemotherapy to a patient having non-small cell lung cancer. This invention provides a prognostic protein biomarker ZNF71(SEQ ID NO:9) using AQUA technique. This invention provides a prognostic mRNA and protein biomarker CD27 (SEQ ID NO:10) with use in immunotherapy. This invention aids patients having non-small cell lung cancer who may benefit from chemotherapy. The protein biomarkers of this invention are new therapeutic targets in chemotherapy and immunotherapy.
IHC results on ZNF71 and CD27 in non-small cell lung cancer (NSCLC) formalin fixed paraffin embedded (FFPE) samples:
The immunohistochemistry assay was performed on ZNF71 and CD27 at Translational Pathology Research Laboratory at West Virginia University. A certified pathologist generated a score for the IHC staining in the following range: 0, 1, 2, 3, and 4, with 0 representing no staining and 4 the maximum staining.
The present invention provides a method of providing a treatment to a patient having non-small cell lung cancer comprising extracting total RNA from a formalin fixed and paraffin embedded tumor of a non-small cell lung cancer of a patient after the surgical resection; generating complementary DNA (cDNA) of the extracted total RNA from said formalin fixed or paraffin embedded patient tumor; quantifying of mRNA expression of 7 genes of ABCC4 (SEQ ID NO:1), CCL19 (SEQ ID NO:2), SLC39A8 (SEQ ID NO:3), CD27 (SEQ ID NO:4), FUT7 (SEQ ID NO:5), ZNF71 (SEQ ID NO:6); and DAG1 (SEQ ID NO:7) using qRT-PCR; normalizing of the quantification of said 7 genes with the quantification of a control gene UBC (SEQ ID NO:8) or a housekeeping gene; and utilizing said normalized 7 gene mRNA expression quantification to determine whether said patient will benefit from receiving adjuvant chemotherapy or not. This method further comprises administering to said patient a therapeutically effective amount of one of the following adjuvant chemotherapies (a) cisplatin and Taxol (paclitaxel), (b) cisplatin and Taxotere (docetaxel), (c) carboplatin, (d) carboplatin and Taxol (paclitaxel), (e) carboplatin and Taxotere (docetaxel), (f) Taxol (paclitaxel), and (g) Alimta (pemetrexed). In this method, the quantification of mRNA expression of three genes of ABCC4 (SEQ ID NO:1), CCL19 (SEQ ID NO:2), and SLC39A8 (SEQ ID NO:3) within said 7-genes. This method includes wherein said quantification of mRNA expression of four genes of CD27 (SEQ ID NO:4), FUT7 (SEQ ID NO:5), ZNF71 (SEQ ID NO:6), and DAG1 (SEQ ID NO:7) within said 7-genes. This method includes wherein said quantification of mRNA expression of ABCC4 (SEQ ID NO:1) and utilization of said normalized ABCC4 gene mRNA expression quantification to determine whether said patient will benefit from receiving adjuvant chemotherapy of one of (a) cisplatin and Taxol (paclitaxel), (b) cisplatin and Taxotere (docetaxel), (d) carboplatin and Taxol (paclitaxel), (e) carboplatin and Taxotere (docetaxel), and (f) Taxol (paclitaxel). This method includes wherein said quantification of mRNA expression of FUT7 (SEQ ID NO:5) and utilization of said normalized FUT7 gene mRNA expression quantification to determine whether said patient will benefit from receiving adjuvant chemotherapy of carboplatin. This method includes wherein said quantification of mRNA expression of ZNF71 (SEQ ID NO:6) and utilization of said normalized ZNF71 gene mRNA expression quantification to determine whether said patient will benefit from receiving adjuvant chemotherapy of (a) carboplatin and Taxol (paclitaxel), (b) carboplatin and Taxotere (docetaxel), (c) cisplatin and Taxotere (docetaxel), and (d) cisplatin and Taxol (paclitaxel). In another embodiment of this method, the method includes wherein said quantification of mRNA expression of SLC39A8 (SEQ ID NO:3) and utilization of said normalized SLC39A8 gene mRNA expression quantification to determine whether said patient will benefit from receiving adjuvant chemotherapy of one of (a) Taxol (paclitaxel), and (b) Alimta (pemetrexed).
In another embodiment of this invention, the method provides a treatment to a patient having non-small cell lung cancer comprising: providing protein expression of ZNF71 (SEQ ID NO: 9); quantifying said protein expression of said ZNF71 with automated quantitative analysis (AQUA) correlated with said ZNF71 mRNA expression in a formalin fixed and paraffin embedded patient tumor using immunohistochemistry staining; and determining a prognosis of said patient from said protein expression of said ZNF71. This method includes wherein said prognosis of said patient is either longer survival or shorter survival.
Another embodiment of this invention provides a treatment to a patient having non-small cell lung cancer comprising: providing protein expression of CD27 (SEQ ID NO:10); quantifying said protein expression of said CD27 with ELISA correlated with said CD27 mRNA expression in a formalin fixed and paraffin embedded patient tumor using immunocytochemistry staining and a cancer-free tissue adjacent to said tumor; and determining a prognosis of said patient from said protein expression of said CD27. This method includes wherein said prognosis of said patient is either longer survival or shorter survival. This method includes administering to said patient a therapeutically effective amount of an adjuvant chemotherapy.
Data on 34,203 lung adenocarcinoma and 26,967 SQCLC patients in linked SEER-Medicare databases were used to determine the contribution of COPD, cancer stage, age, gender, race, and tumor grade to prognostication in 30 treatment combinations. A Cox model including these variables was estimated on 1,000 bootstrap samples, with the resulting model assessed on ROC, Brier Score, Harrell's C, and Nagelkerke's R2 metrics. The comprehensive model was evaluated with two additional patient cohorts (n=1,994). The results show that combining patient information on COPD, cancer stage, age, gender, race, and tumor grade improves prognostication and prediction of treatment response in individual NSCLC patients receiving surgery, radiation, and chemotherapy, including Platinum-based, Platinum/Taxane, and Carboplatin/Paclitaxel/Avastin (Putila, Remick, and Guo 2011[51]; Putila and Guo 2014) [50]). An example of an online web-based model as a prognostic tool, is the model available at www.personalizedrx.org which has been used in the MBRCC clinic to advise patient treatment selection (see,
Integration of IHC of ZNF71, immunocytochemical staining of CD27, qRT-PCR of 7 gene assay in FFPE samples with PersonalizedRx:
The present invention provides a method that provides a comprehensive prognostic model combining COPD, age, gender, race, histology, AJCC staging edition, cancer stage and tumor grade using multivariate Cox model with SEER-Medicare data (Putila, Remick, and Guo 2011 [51]; Putila and Guo 2014 [50]). All the model covariates are available in our clinical cohorts. Since all NSCLC patients receive pulmonary function tests before surgery, we will capture FEV1 and DLCO parameters to refine the diagnosis of COPD and its coefficient in the comprehensive prognostic model. The molecular biomarkers will be integrated with this model as independent covariate(s) with coefficients of other clinical covariates adjusted in the Cox model using training clinical cohort; the new model parameters will be validated with multiple independent clinical cohorts. The IHC scores of ZNF71 will be used as a co-variate in the multivariate Cox proportional hazard model used to contruct the PersonalizedRx tool (web-based model). Similar, the immunocytochemical staining results of CD27 in T cells and B cells will be used as co-variates in the above model, so will be the output from the qRT-PCR of the 7-gene assay in FFPE patient samples. We have published such analysis in previous studies (Wan et al. 2012 [54]). Table 3 is an example of the results from the analysis.
An interactive web interface of the current comprehensive web-based model available at http://www.personalizedrx.org has been used in clinics to aid treatment intervention. Validated molecular biomarkers will be added into this tool for improved cancer care.
This represents an additional prognostication improvement over the use of cancer stage alone, which has already been validated with statistical rigor [51]. This comprehensive model enables refined prognosis and estimation of clinical outcome of treatment combinations in NSCLC patients, providing a useful tool in personalized clinical decision-making. The comprehensive web-based model online tool commercially available at www.personalizedRx.org is employed herein the method of this invention and is in use at clinics at Mary Babb Randolph Cancer Center (MBRCC), West Virginia University. With the advancement of clinical genomics research, this comprehensive prognostic model is integrated with the genomic biomarkers to predict NSCLC patient treatment response in the methods of this invention.
ψdenotes confidence interval.
Each co-morbid condition was assessed as an independent predictor of survival using Cox proportional hazards modeling in patients treated with surgery but without chemotherapy or radiation indicated in order to isolate the effects of comorbidity from those of disparate treatment benefit or treatment candidacy (Supplementary Table 1). Additionally, the presence of COPD as determined via the analysis of administrative records was assessed as an independent predictor of survival by testing for significant stratification of Kaplan-Meier survival curves. Patients were split into two outcome groups based on COPD status, and separate survival curves were estimated and plotted (
The distributions of COPD and other variables in the original model were assessed in patients with very long and very short survival, relative to other patients, to determine if certain characteristics were disparate between groups of patients with varied survival. This was accomplished by partitioning patients based on survival time and status, then using a t-test or test of proportions to compare the distributions of variables between each group. Again, only those patients who were treated with surgery without radiation or chemotherapy were included. Long survival was defined as greater than 60 months for the original 3rd Edition staging, and greater than 24 months for the 6th and recoded 7th Edition groups due to shortened follow-up. Short survival was defined as less than 24 months for the original 3rd Edition staging, and less than 12 months for the 6th and recoded 7th Edition groups (Supplementary Table 3).
COPD showed significant prognostic ability on multiple measures, both as an independent predictor and in the presence of other predictors. Other co-morbid conditions also showed promise as independent predictors in a Cox model (Supplementary Table 2). As an independent predictor, COPD status alone was able to significantly stratify patients into high and low-risk groups (p<0.05) in four of six groups (
The proportion of patients with COPD between those with relatively long and short survival was also assessed. Two survival cutoffs were used to split patients into three groups of short, intermediate, and long survival in order to test for differences in the distribution of each prognostic factor between groups of patients with relatively different survival. The difference in the prevalence of COPD between the short and long survival groups was assessed using a test of proportions. This test showed that COPD was much more prevalent in patients with relatively short survival when compared to those surviving relatively longer (Supplementary Table 4). This was true for each histology and coding scheme, despite differences in cutoffs and length of follow-up between the original and recoded staging systems. When the same test was performed for the other covariates similar results were seen, with factors previously seen to favor increased or diminished survival being disparate between the groups. These results are summarized in Supplementary Tables 4 through 9.
Patients were able to be split into high and low-risk groups with significantly different survival curves using COPD status alone in a variety of treatment and stage sub-groups. For adenocarcinoma patients staged using the original 3rd Edition staging, there was a significant difference in the survival of Stage I patients treated with surgery alone. There was also a significant difference seen in Stage 2 and 3a patients treated with surgery with a platinum and taxane or without systemic therapy (p<0.05). A fourth stratification in Stage 2 and 3a patients treated with radiation and systemic therapy was also present (p<0.05,
In the group of patients staged using the original 6th Edition, Stage 2 and 3a radiation patients treated both with and without systemic therapy trended toward a similar stratification but did not achieve significance (p>0.05). A similar trend was observed for Stage 2 and 3a surgical patients treated with systemic therapy, although the degree of separation was again not significant (p>0.05). The group of stage 1 surgical patients treated without systemic therapy did however achieve a significant stratification (p=0.0111,
In the group of adenocarcinoma patients recoded to AJCC 7th Edition staging, there was a significant difference in survival between patients with and without COPD in Stage I surgical patients treated without systemic therapy (p=0.0374). This same difference was also present in Stage 2 and 3a surgical patients treated without systemic therapy (p=0.0101). There was a similar trend in the group of Stage 2 and 3a patients treated with surgery and systemic therapy, but the difference in survival was not significant for any of the systemic therapy groups (p>0.05). The group of Stage 2 and 3a patients without COPD also experienced significantly better survival when treated with radiation without systemic therapy (p=0.0025,
In squamous cell patients staged using the original 3rd Edition the difference in survival when stratifying on COPD status for Stage I surgical patients treated without systemic therapy was significant, with patients having COPD again experiencing poorer survival (p=0.0002). There was also a significant stratification in the corresponding group of Stage I surgical patients treated with systemic therapy (p<0.05). COPD was able to produce a significant stratification in Stage 2 and 3a patients treated with systemic therapy in both the surgical and radiation groups (p<0.05,
In the recoded 7th Edition staging, there were two groups where COPD was able to stratify patients. The first was in Stage 2 and 3a patients treated with surgery without any systemic therapy, with COPD patients having poorer survival (p=0.0058,
It will be appreciated by those persons skilled in the art that changes could be made to embodiments of the present invention described herein without departing from the broad inventive concept thereof. It is understood, therefore, that this invention is not limited by any particular embodiments disclosed, but is intended to cover the modifications that are within the spirit and scope of the invention, as defined by the appended claims.
A SEQUENCE LISTING in computer-readable form (.txt file) accompanies this application having SEQ ID NO:1 through SEQ ID NO:10. The computer-readable form (.txt file) of the SEQUENCE LISTING is incorporated by reference into this application. The SEQUENCE LISTING in computer-readable form (.txt file) is electronically submitted along with the electronic submission of this application.
This patent application claims the benefit of co-pending PCT/US2019/036953, filed, Jun. 13, 2019. The entire contents of PCT/US2019/036953 is incorporated by reference into this patent application as if fully written herein.
This invention was made with government support under National Institute of Health Grants RO1 LM009500, R56 LM009500, RO1 ES021764, and P20 RR016440. The government has certain rights in the invention.
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/US2020/023597 | 3/19/2020 | WO |
Number | Date | Country | |
---|---|---|---|
Parent | PCT/US2019/036953 | Jun 2019 | US |
Child | 17906315 | US |