Abstract Motivation: Overexpression of HER2 (the product of the ERBB2 gene) occurs in about 15% of all breast tumors. We have undertaken to use next generation transcriptome sequencing technology to identify genomic features that are unique to HER2+ tumors. Interactome mapping was then used to integrate the genes associated with these features into a transcriptome landscape model, with a view towards identifying nodes of interaction that might be targeted in HER2+ tumors. Methods: We performed 50nt paired-end RNA-sequencing for 24 breast tumors: 8 each HER2+, ER+, triple negatives (TN). In addition to breast adenocarcinomas, we also sequenced 8 early passage non-transformed HMEC cell lines as normal controls. Reads from RNA sequencing were aligned to the genome and exon junctions using TopHat software. Gene counts were summarized and annotations were performed using our in-house programs. Tukey's test was used to obtain genes or transcripts that are specific to HER2 tumor group compared to other tumors/normal. A combination of bioinformatics softwares and algorithms were used to identify SNVs. Results: Only 13527 genes with median gene count greater than 16 reads in at least one of the 4 groups were considered for differential gene expression and splicing analysis. Some 696 genes were differentially expressed in HER2+ tumors compared to ER+, TN tumors and HMECs. We identified 272 alternately spliced transcripts for which the HER2+ tumors exhibited a mean transcript expression ratio significantly different from the means of other tumor groups. We also identified 4735 expressed single nucleotide variants that are uniquely associated with HER2+ tumors compared to other tumors/normal groups. Among these 3579/4735 sequence polymorphisms were not present in the 1000 genome germline sequence database or in the dbSNP132 database of naturally occurring germline polymorphisms. Integration of all the genes obtained from genomic feature analyses (differential expression, alternative splicing, single nucleotide variance) has been carried out to indentify biological processes that are specific to the HER2+ tumor subtype. Key nodes and pathways that are specific to HER2+ tumors will be evaluated for association with clinical outcome in a large series of patients who have received HER2-targeted therapy. Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 103rd Annual Meeting of the American Association for Cancer Research; 2012 Mar 31-Apr 4; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2012;72(8 Suppl):Abstract nr 4926. doi:1538-7445.AM2012-4926
Six Cycles of Doxorubicin and Cyclophosphamide or Paclitaxel Are Not Superior to Four Cycles As Adjuvant Chemotherapy for Breast Cancer in Women With Zero to Three Positive Axillary Nodes: Cancer and Leukemia Group B 40101
Journal of Clinical Oncology | July 24, 2012 | Volume: 30 | Issue: 33 | Pages: 4071-4076
PurposeThe ideal duration of adjuvant chemotherapy for patients with lower risk primary breast cancer is not known. Cancer and Leukemia Group B trial 40101 was conducted using a phase III factorial design to define whether six cycles of a chemotherapy regimen are superior to four cycles. We also sought to determine whether paclitaxel (T) is as efficacious as doxorubicin/cyclophosphamide (AC), but with reduced toxicity.Patients and MethodsBetween 2002 and 2008, the study enrolled women with operable breast cancer and zero to three positive nodes. Patients were randomly assigned to either four or six cycles of either AC or T. Study stratifiers were estrogen receptor/progesterone receptor (ER/PgR), human epidermal growth factor receptor 2 (HER2), and menopausal status. After 2003, all treatment was administered in dose-dense fashion. The primary efficacy end point was relapse-free survival (RFS).ResultsA total of 3,171 patients were enrolled; 94% were node-negative and 6% had one to three positive nodes. At a median follow-up of 5.3 years, the 4-year RFS was 90.9% and 91.8% for six and four cycles, respectively. The adjusted hazard ratio (HR) of six to four cycles regarding RFS was 1.03 (95% CI, 0.84 to 1.28; P = .77). The 4-year OS was 95.3% and 96.3% for six and four cycles, respectively, with an HR of six to four cycles of 1.12 (95% CI, 0.84 to 1.49; P = .44). There was no interaction between treatment duration and chemotherapy regimen, ER/PgR, or HER2 status on RFS or OS.ConclusionFor women with resected primary breast cancer and zero to three positive nodes, we found no evidence that extending chemotherapy regimens of AC or single-agent T from four to six cycles improves clinical outcome.
Reexpression of Tumor Suppressor, sFRP1, Leads to Antitumor Synergy of Combined HDAC and Methyltransferase Inhibitors in Chemoresistant Cancers
Molecular Cancer Therapeutics | July 24, 2012 | Volume: 11 | Issue: 10 | Pages: 2105-2115
Abstract Motivation: KRAS is commonly mutated in a variety of cancers including lung cancer. The KRAS gene is frequently mutated at codons 12 and 13 in lung adenocarcinomas in patients with a history of smoking. Tumors harboring an activating KRAS mutation are aggressive and are often resistant to available therapies. In the present study, we set out to identify pathways that are specifically altered in lung cancer patients whose tumors harbor a KRAS mutation using next generation sequencing technology. Methods: We performed 50nt paired-end RNA-sequencing in 15 lung adenocarcinoma samples (8 and 7 samples with and without KRAS mutation). Reads were aligned to genome and exon junctions using Illumina's alignment tool ELAND_RNA. CASAVA and Genome Studio data analysis software were used to obtain read counts for genes, exons and exon junctions. We applied a variety of computational methods (Casper R package, Bowtie, TopHat) and softwares (Partek, R statistical software, JMP, Ingenuity Pathway Anlaysis, DAVID) to carry out our analyses. Results: Differential gene expression and splicing analysis were performed between the two groups. We identified 115 genes that were differentially expressed and 112 genes that consist of splicing variants with at least 2 fold changes and p-value < 0.01. We randomly selected 6/15 samples and performed real-time qRT-PCR for 6 differentially expressed genes. Significant correlations were observed ranging from 0.44 to 0.83, when qRT-PCR results were compared with RNA-sequencing expression data. Pathway analysis with 115 differentially expressed genes and 112 splicing variants revealed that the most significant pathway is composed entirely of NFκB focus genes, suggesting that there is a direct connection between oncogenic KRAS and activation of the NFκB signaling pathway. Conclusions: NFκB has been implicated in KRas-mediated formation of early stage lung adenomas in the LSL-Kas mouse model. Our data indicate that this connection also occurs in human lung cancer and establish NFκB activation as a key manifestation of KRAS mutation in human lung cancer. Hence drugs targeting NFκB and NFκB related genes may potentially be helpful for the treatment of patients with oncogenic KRAS mutations. (Supported in part by grants from the 26.2 with Donna Foundation and the Florida Department of Health Bankhead/Coley program.) Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 102nd Annual Meeting of the American Association for Cancer Research; 2011 Apr 2-6; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2011;71(8 Suppl):Abstract nr 4975. doi:10.1158/1538-7445.AM2011-4975
Abstract 4908: Deep sequence analysis of the relationship between gene expression, CpG island methylation, and gene copy number in breast cancer cells
Cancer Research | June 27, 2012 | Volume: 71 | Issue: 8_Supplement | Pages: 4908-4908
Abstract Background: Estrogen receptor (ER) expression in breast caner is an important biomarker for targeted therapy and outcome prediction. Nevertheless, the molecular mechanism that controls the phenotypic difference is little known. We used deep sequence technology to profile the transcriptome, gene copy number, and CpG island methylation simultaneously in eight commonly used breast cell lines to develop a model for how these genomic features are integrated in ER+ and ER- breast cancer. Materials and methods: Total mRNA sequence (mRNA-seq), gene copy number (DNA-seq), and genomic CpG island methylation (Methyl-seq) were carried out using the Illumina Genome Analyzer. Sequences were mapped to the human genome (hg18) to obtain digitized gene expression data, DNA copy number in reference to the non-tumor cell line (MCF10A), and methylation status of 21,570 CpG islands; differentially expressed genes between ER+ and ER- cell lines were selected and then correlated with methylation status of these genes’ CpG islands within 5kb of transcript start or copy number changes. The genes that appeared controlled by methylation or copy number changes in cell lines were further evaluated for their expression in the dataset of 129 primary breast tumors. Results: ER+ breast cancer cells had very different gene expression pattern from ER- cancer cells. The two cell line types formed two distinct clusters in unsupervised cluster and 1,873 genes were differentially expressed by moderated t statistics. We identified 149 differentially expressed genes that exhibited differential methylation of one or more CpG islands within 5kb of the 5’ end of the gene and for which mRNA abundance was inversely correlated with CpG island methylation status. Eighty nine of these 149 genes were also differentially expressed in the primary tumor samples and 84 of them were consistent with cell line gene expression and methylation data. The set of 149 genes were significantly enriched in the ER+ tumors by Gene Set Enrichment Analysis. Part of chromosome 8 was deleted in all ER- cells and part of chromosome 17 amplified in all ER+ cells. These loci encoded 30 genes that were overexpressed in ER+ cells; however, only 9 of these genes were overexpressed in ER+ tumors. Many of the methylation or CNA affected genes were functionally significant and some are known to be associated with breast cancer outcome such as GATA3 and LYN. Conclusions: A global pattern of differential CpG island methylation influences expression of a cohort of genes that contribute to the transcriptome landscape of ER+ and ER- breast cancer cells and tumors. The role of gene amplification/deletion in defining the transcriptome landscape of ER+ and ER- cells appears to more modest, although several potentially significant genes appear to be globally regulated by copy number aberrations. Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 102nd Annual Meeting of the American Association for Cancer Research; 2011 Apr 2-6; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2011;71(8 Suppl):Abstract nr 4908. doi:10.1158/1538-7445.AM2011-4908
Abstract LB-276: Fusion transcripts in breast cancer cell lines and tumors
Cancer Research | June 27, 2012 | Volume: 71 | Issue: 8_Supplement | Pages: LB-276-LB-276
Abstract The role of fusion gene transcripts such as BRC-ABL1 has been long appreciated in hematopoietic malignancies. Recent evidence from next generation sequencing projects has suggested that gene fusion events and the resultant chimaeric transcripts may be expressed in solid tumors. We have developed a novel bioinformatics analytical pipeline to detect fusion gene transcripts from paired-end mRNA-seq data. The bioinformatics tool used for fusion transcript discovery employs multiple steps of false positive filtering and nominates the fusion transcript candidates with high confidence (approaching 100%). The screening and validation of the fusion candidates were quickly carried out using the recommended primer design regions as one of the outputs of the SnowShoes-FD pipeline. This pipeline was used to analyze the transcriptome of 22 breast cancer cell lines and 9 non-transformed breast cell populations. Fifty-four fusion candidates were nominated, all of which were validated using reverse transcription PCR and Sanger sequencing. Five fusion products were predicted to have a second fusion junction point between two fusion partners. These 54 fusion transcripts consist of 103 partner genes that are involved in cell cycle and/or nuclear receptor signaling. No fusion transcripts were identified from the non-transformed breast cell lines. We subsequently analyzed a panel of primary breast tumors, 8 each HER2+, triple negative, and ER+. Fifteen novel fusion transcript candidates have been nominated. Validation and functional analysis of these candidates is in progress. Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 102nd Annual Meeting of the American Association for Cancer Research; 2011 Apr 2-6; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2011;71(8 Suppl):Abstract nr LB-276. doi:10.1158/1538-7445.AM2011-LB-276
Abstract C114: Androstenedione levels in postmenopausal women with resected early-stage breast cancer are associated with SNPs in CYP11B1 and CYP11B2 identified by a genome-wide association study (GWAS).
Molecular Cancer Therapeutics | June 27, 2012 | Volume: 10 | Issue: 11_Supplement | Pages: C114-C114
Abstract Introduction: Androstenedione is a precursor for estrogen biosynthesis catalyzed by aromatase. Plasma androstenedione concentrations in postmenopausal women with early stage breast cancer display large interindividual variation (Ingle et al., Cancer Res. 70:3278–86, 2010). We performed a GWAS to identify single nucleotide polymorphisms (SNPs) associated with plasma levels of androstenedione in 776 postmenopausal women with resected early stage breast cancer. Methods: Plasma androstenedione concentrations were measured by GC-MS/MS. GWAS was performed using Illumina Human610-Quad BeadChips. The GWAS results were controlled for population stratification (PS) by Eigenstrat analyses. Plasma androstenedione concentrations were associated with genotype by using a quasi-likelihood model with adjustment for PS, BMI, age and other clinical variables associated with androstenedione concentrations (p<0.01). Imputation using HapMap2 data was performed within 200 kb on either side of regions containing SNPs with p<E-05. Results: 563,945 SNPs were analyzed, followed by imputation. The genotyped SNP with the lowest p value, rs4736317 (p=4.61E-06), was on chromosome 8 near the CYP11B1 and CYP11B2 genes. This SNP had a multiplicative effect of 1.128 (95% Cl: 1.07, 1.19) and a MAF of 0.40. After imputation, 20 SNPs with apparent p<E-06 were observed in this same region, with the lowest p value for rs9773022 (p=8.1E-08). We then performed functional genomic studies, and observed that two of the SNPs with low p values (rs4736312, p=1.49E-05 and rs4736350, p=6.34E-06) created apparent Oct-1 transcription factor binding sites (TFBS). Differential nuclear protein binding of Oct-1 to probes for these SNPs was verified with EMS and supershift assays using nuclear extract from A549 and NCI-H295R cells. Treatment of these two cell lines with 10 nM androstenedione led to the disruption of Oct-1 binding with the A (wild type) allele on rs4736312. rs4736312 is located in the 3'-UTR of CYP11B1 15 bp distant from response elements for the progesterone receptor and the glucocorticoid receptor. We also demonstrated that CYP11B1 and CYP11B2 were induced by estradiol in U2OS cell lines that express ER receptor and by androstenedione in H295R cells that express the androgen receptor. Transfection of a reporter gene construct with the rs4736312 SNP sequences into COS-1 and H295R cells suppressed reporter gene activity for the variant SNP sequence as compared to wild type, an observation that may correlate with the higher plasma levels of androstenedione in women (549× 328 vs 410×208 mean×SD, M/ml) with the variant allele. Conclusion: We have identified a series of SNPs in a region of chromosome 8 in and around CYP11B1 and CYP11B2 that are associated with elevated plasma androstenedione concentrations in postmenopausal women with resected early stage breast cancer. Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the AACR-NCI-EORTC International Conference: Molecular Targets and Cancer Therapeutics; 2011 Nov 12-16; San Francisco, CA. Philadelphia (PA): AACR; Mol Cancer Ther 2011;10(11 Suppl):Abstract nr C114.
Choosing the Best Trastuzumab-Based Adjuvant Chemotherapy Regimen: Should We Abandon Anthracyclines?
Journal of Clinical Oncology | May 22, 2012 | Volume: 30 | Issue: 18 | Pages: 2179-2182
Purpose NCCTG (North Central Cancer Treatment Group) N9831 is the only randomized phase III trial evaluating trastuzumab added sequentially or used concurrently with chemotherapy in resected stages I to III invasive human epidermal growth factor receptor 2–positive breast cancer. Patients and Methods Patients received doxorubicin and cyclophosphamide every 3 weeks for four cycles, followed by paclitaxel weekly for 12 weeks (arm A), paclitaxel plus sequential trastuzumab weekly for 52 weeks (arm B), or paclitaxel plus concurrent trastuzumab for 12 weeks followed by trastuzumab for 40 weeks (arm C). The primary end point was disease-free survival (DFS). Results Comparison of arm A (n = 1,087) and arm B (n = 1,097), with 6-year median follow-up and 390 events, revealed 5-year DFS rates of 71.8% and 80.1%, respectively. DFS was significantly increased with trastuzumab added sequentially to paclitaxel (log-rank P < .001; arm B/arm A hazard ratio [HR], 0.69; 95% CI, 0.57 to 0.85). Comparison of arm B (n = 954) and arm C (n = 949), with 6-year median follow-up and 313 events, revealed 5-year DFS rates of 80.1% and 84.4%, respectively. There was an increase in DFS with concurrent trastuzumab and paclitaxel relative to sequential administration (arm C/arm B HR, 0.77; 99.9% CI, 0.53 to 1.11), but the P value (.02) did not cross the prespecified O'Brien-Fleming boundary (.00116) for the interim analysis. Conclusion DFS was significantly improved with 52 weeks of trastuzumab added to adjuvant chemotherapy. On the basis of a positive risk-benefit ratio, we recommend that trastuzumab be incorporated into a concurrent regimen with taxane chemotherapy as an important standard-of-care treatment alternative to a sequential regimen.