Citations

Abstract Purpose: The limited supply of oxygen and nutrients is thought to result in rigorous selection of cells that will eventually form the tumor. Experimental Design: Nupr1 expression pattern was analyzed in human tissue microarray (TMA) and correlated with survival time of the patient. Microarray analysis was conducted on MiaPaCa2 cells subjected to metabolic stress in Nupr1-silenced conditions. DNA repair and cell cycle–associated gene expression was confirmed by real-time quantitative PCR (qRT-PCR). Nupr1 and AURKA protective role were analyzed using RNA interference (RNAi) silencing or overexpression. DNA damage and autophagy were analyzed by Western blot analysis and immunofluorescence. Results: We showed that both Nupr1 and HIF1α are coexpressed in human pancreatic ductal adenocarcinoma (PDAC) samples and negatively correlate with survival time. PDAC-derived cells submitted to hypoxia and/or glucose starvation induce DNA damage–dependent cell death concomitantly to the overexpression of stress protein Nupr1. Affymetrix-based transcriptoma analysis reveals that Nupr1 knockdown enhances DNA damage and alters the expression of several genes involved in DNA repair and cell-cycle progression. Expression of some of these genes is common to hypoxia and glucose starvation, such as Aurka gene, suggesting that Nupr1 overexpression counteracts the transcriptional changes occurring under metabolic stress. The molecular mechanism by which hypoxia and glucose starvation induce cell death involves autophagy-associated, but not caspase-dependent, cell death. Finally, we have found that AURKA expression is partially regulated by Nupr1 and plays a major role in this response. Conclusions: Our data reveal that Nupr1 is involved in a defense mechanism that promotes pancreatic cancer cell survival when exposed to metabolic stress. Clin Cancer Res; 18(19); 5234–46. ©2012 AACR.

Abstract High throughput sequencing technologies open up a new dimension in cancer genomics by enabling the characterization of cancer genomes at base-pair resolution. To date, the majority of large-scale genomics projects rely on fresh-frozen biospecimens in their studies, which are difficult to obtain and often lack long-term clinical data. Unlike fresh-frozen samples, archived formalin-fixed paraffin-embedded (FFPE) tissues are more readily accessible, and are often associated with known clinical outcomes and more complete clinical annotations. However, sequencing library preparation methods need to be further optimized with regard to applicability to FFPE samples. To test whether next-generation sequencing technologies could overcome previously reported artifacts associated with formalin fixation and report accurate sequencing results, we compared whole-genome and targeted enrichment DNA sequencing data obtained from five FFPE tumor samples for which matching frozen tissues were available. We successfully generated sequencing libraries using 1 µg of genomic DNA as input material and a modified TruSeq™ sample preparation protocol. We designed a custom enrichment pool targeting exons from 215 cancer-related genes and utilized this pool to pull-down and sequence approximately 1.3 Mb region of the genome with an average 400x coverage depth to test if deep sequencing would improve validation of tumor-specific somatic mutations. We used the Illumina sequence analysis pipeline and the Windows-based second generation DNA sequencing software NextGENe® to analyze the data and identify sequence variations that are different from the human reference genome. CNV detection was overall higher among FFPE samples compared with fresh-frozen samples, demonstrating that tissue processing impacts sequencing data quality. We obtained good concordance in variant calls between matched FFPE and fresh-frozen samples. Discordant variant calls were mainly due to low depth of coverage in the regions where variant calls were made. Improvements to DNA sequencing methods for archived samples will significantly enhance cancer research and will result in more reliable prediction of individual cancer therapies. 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 3179. doi:1538-7445.AM2012-3179

Abstract Little is known about the effects of many recognized breast cancer risk factors on the normal breast. We hypothesize that risk factors affect the number and morphometry of terminal ductal lobular units (TDLUs), the microanatomical structure from which most breast cancers arise. To assess this hypothesis, we evaluated relationships between breast cancer risk factors and TDLU number, TDLU diameter, and acini per TDLU in breast tissue donated by healthy volunteers. Subjects included 786 healthy women ages 18-84 years without a personal history of any cancer, who provided a risk factor questionnaire and donated a 10-guage needle core from the upper outer breast quadrant. A pathologist, masked to subject data, assessed TDLU number per section, TDLU diameter and acini per TDLU using digitized images of hematoxylin and eosin stained sections of tissue. Seventy-two randomly selected images were examined by two additional pathologists. We generated ordered categorical variables for the number of TDLUs per section (0-1, 2-8, 9+), median TDLU diameter (53-262, 263-377, 378-1375 microns) and median number of acini per TDLU (2-10, 11-20, 21+). We computed Spearman correlations between measures, and intraobserver and interobserver agreement. Associations between breast cancer risk factors and TDLU measures were assessed by estimating odds ratios (OR) and 95% confidence intervals (CI) using ordinal logistic regression models, adjusted for age. TDLU diameter and acini count were highly correlated (rho =0.71); TDLU number per section was less strongly related to the diameter (rho= 0.18) and acini measures (rho=0.20). Inter and intraobserver agreement for TDLU measures were highly correlated (rho range 0.68-0.98). All TDLU measurements were significantly lower among women age 50 years and older compared with younger women. Compared with Caucasians, African American women had more TDLUs per section (ORadj=2.03, 95% CI= 1.06-3.89). Greater numbers of TDLUs were also related to having a family history of two or more relatives with breast cancer compared with none (ORadj=2.01, 95% CI=0.91-4.46); and for parous compared with nulliparous women (ORadj=1.61, 95% CI=1.16-2.25). In contrast, parity was associated with reduced median number of acini per TDLU (ORadj=0.60, 95% CI = 0.40-0.90) and median diameter (ORadj=0.64, 95% CI = 0.42-0.96) compared with nulliparity. Our data suggest that TDLUs are more frequent among women with a positive family history of breast cancer or African Americans. Parity, an established protective factor for breast cancer risk, was related to increased numbers of TDLUs, but fewer acini per TDLU, suggesting that TDLUs may increase with pregnancy and then involute. Our data, along with published results linking TDLU involution with a reduction in breast cancer risk, suggest TDLU morphometry may reflect early biologic changes marking breast cancer susceptibility. 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 4465. doi:1538-7445.AM2012-4465

Abstract Background: Colon cancer is categorized by the acquisition of a number of molecular alterations in both oncogenes and in tumor suppressor genes. Although cMYC alterations have been reported in colon cancer cell lines, little other information exists in primary colon tumors. We sought to determine the frequency of cMYC amplification in primary colon cancer and its relationship to other well defined colon tumor alterations. Methods: An initial set of 392 tumors were evaluated. All samples were tested for the BRAF V600E mutation and all were classified as proficient mismatch repair (pMMR) based on both IHC and microsatellite analysis of tumor tissue. A second set of 522 tumors was selected to replicate associations identified in the first set of samples. Similar to the initial cohort, both the BRAF V600E and the mismatch repair status of the tumors were assessed. The second cohort, however, included tumors that were both pMMR and dMMR. In both sets of samples, fluorescent in situ hybridization was used to evaluate cMYC gene amplification. Results: In cohort 1, cMYC amplification was observed in 15 of the 362 (4%) pMMR tumors successfully analyzed (Table 1). cMYC amplification was detected in 10 of 323 (3%) tumors with normal BRAF status vs. 5 of 39 (13%) tumors with altered BRAF status (p=0.015). In the second cohort, cMYC amplification was found in 20 of the 522 (4%) tumors. Only 1 of the 152 (<1%) dMMR tumors had cMYC amplification vs. 19 of the 370 (5%) pMMR tumors (p=0.012). Of the pMMR tumors, cMYC amplification was detected in 13 of 333 (4%) tumors with normal BRAF status vs. 6 of 37 (16%) tumors with altered BRAF status (p=0.007), nearly identical to the first group. Conclusion: cMYC amplification in colon cancer is an infrequent event. However, when present, cMYC amplification occurs almost exclusively in tumors that are pMMR. Additionally, among pMMR tumors, cMYC amplification is associated with the presence of the BRAF V600E mutation. 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 1151. doi:1538-7445.AM2012-1151

Abstract Background: Estrogen receptors are over-expressed in around 70% of breast cancer cases. The genetic changes that occur during aromatase inhibitor (AI) treatment are not well understood and may differ depending upon the patient's response phenotype. Methods: We performed whole genome sequencing (WGS) of matched blood, pre-treatment, and post-treatment biopsy samples from 22 estrogen receptor positive breast cancer patients treated with neoadjuvant aromatase inhibitors. For 5 cases, we performed the whole genome sequencing (WGS) on patients’ matched normal, two pre AI-treatment, and two post AI-treatment DNA isolates from biopsy samples. We validated all putative coding and non-coding somatic mutations using deep sequencing. By comparing the validated somatic mutations from pre- and post- AI treatment biopsy samples, we were able to determine the alterations in the tumor genomes. In every case we defined the clonal architecture of each pair of pre-treatment and post-treatment biopsy samples by comparing the variant allele frequencies from thousands of validated somatic mutations. Results: Comparisons of the two pre AI-treatment biopsy samples from the same patient indicates that the variant allele frequencies of mutations showed high concordances in all 5 cases, 0.74 to 0.95 range of correlation coefficient. Only a small percentage of somatic mutations were detected in one pre-treatment sample and not the other (4.65% overall). In comparing the somatic variations between pre-treatment and matched post-treatment biopsy samples in 22 cases, we found that patients with good clinical response to AI treatment retained known driver mutations only in their pre-treatment tumors. Conversely, those patients with poor clinical response presented new driver mutations in their post-treatment samples. Furthermore, the variant allele frequency for most mutated genes decreased in post AI treatment samples for patients with good AI treatment response; on the contrary, the variant allele frequency increased for patients with poor clinical response. Conclusions: From WGS of matched normal, pre-treatment, and post-treatment biopsy samples, we identified new driver genes mutated in patients with poor clinical response, while patients with good clinical response had lost mutated driver genes in their post-treatment biopsy samples. The genetic landscape revealed by WGS of pre-treatment and post-treatment biopsy samples reveals mutational repertoires are remodeled by AI therapy. This finding suggests deep sequencing of AI treated samples will be necessary to reveal the complete complement of mutations present in a patient's tumor. 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 LB-423. doi:1538-7445.AM2012-LB-423

Abstract A substantial proportion of familial colorectal cancer (CRC) is not a consequence of known susceptibility loci, such as mismatch repair (MMR) genes, supporting the existence of additional loci. To identify novel CRC loci, we conducted a genome-wide linkage scan in 356 white families with no evidence of defective MMR (i.e., no tumors with loss of expression of MMR proteins; no microsatellite instability (MSI)-high tumors, and no evidence of linkage to MMR genes). This represents the largest linkage analysis of proficient MMR (pMMR) familial CRC to date. Families were ascertained via the NCI-supported Colon Cancer Family Registry multi-site (Colon CFR) consortium, the City of Hope, and Memorial University of Newfoundland. A total of 1,612 individuals (average 5.0 per family, range 2-10 with an average of 2.2 affected and 2.8 unaffected individuals) were genotyped using genome-wide single nucleotide polymorphism linkage arrays; parametric and non-parametric linkage analysis used MERLIN in a priori-defined family groups. Five lod scores greater than 3.0 were observed, accounting for heterogeneity, in four chromosomal regions. The greatest lod scores were at 4q21.1 among families with mean age of diagnosis less than 50 years (dominant HLOD=4.51, α=0.84, 145.40 cM, rs10518142) and at 12q24.32 among all families (dominant HLOD=3.60, α=0.48, 285.15 cM, rs952093). Among families with four or more affected individuals and among clinic-based families, a common peak was observed at 15q22.31 (101.40 cM, rs1477798; dominant HLOD=3.07, α=0.29; dominant HLOD=3.03, α=0.32, respectively). Analysis of families with only two affected individuals yielded a peak at 8q13.2 (recessive HLOD=3.02, α=0.51, 132.52 cM, rs1319036). These linkage regions comprise previously unreported loci, supporting the hypothesis that novel loci contribute to CRC in familial pMMR cases. These loci are also largely distinct from those identified in case-control GWAS. To follow up on these findings, a custom Agilent Sure Select Target Enrichment panel has been designed for genes in three of these regions, in order to sequence probands using an Illumina HiSeq 2000. This study demonstrates the utility of family-based data and provides evidence for additional alleles for familial CRC. 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 2930. doi:1538-7445.AM2012-2930

Abstract Mi(cro)RNAs are short non-coding RNA molecules that play a key role in carcinogenesis by regulating tumor suppressors and oncogenes. Published data has implicated miRNAs in ovarian cancer (OC) development and progression, and we hypothesize that single nucleotide polymorphisms (SNPs) in sites of miRNA: messenger RNA (mRNA) binding may influence OC risk by altering expression levels of targeted mRNAs. To evaluate associations between SNPs in miRNA binding sites and OC risk, we used bioinformatics tools and public databases to identify miRNA binding SNPs relevant to ovarian cancer. We then evaluated the frequency of approximately 1,100 SNPs in ∼20,000 epithelial OC cases and 20,000 controls represented in the international Collaborative Oncological Gene-environment Study. Unconditional logistic regression was used to estimate odds ratios (OR) and 95% confidence intervals (CI) between genotypes and case status, with adjustment for the first five principal components representing European ancestry. Log-additive genetic models were applied to each SNP, modeling the number of copies of the minor allele. Subgroup analysis was conducted for serous adenocarcinomas, the most predominant histologic subtype of epithelial OC. Preliminary analysis of 11,215 invasive epithelial OCs and 15,910 controls included four SNPs with P<10−6 and 1 SNP with P<10−5. Analysis of 6,215 invasive serous adenocarcinomas revealed two SNPs with P<10−7 and three SNPs with P<10−6. Noteworthy is rs6104808, a C>T SNP predicted to reside in a highly conserved miRNA binding site of BTBD3 (BTB POZ domain containing 3), a transcription factor suggested to be involved in cell proliferation, cancer progression, and response to platinum-based chemotherapy. rs6104808 was associated with a decreased risk of OC overall (OR (95% CI): 0.42 (0.28-0.62), P=1.26 x 10−5)), and the association was stronger among cases with serous adenocarcinoma (OR (95% CI): 0.30 (0.17-0.50), P=9.89 x 10−6)). This represents the largest, most comprehensive epidemiologic study to date to evaluate associations between miRNA binding site SNPs and OC susceptibility. Although based on a preliminary analysis, these data suggest a significant role and putative biological mechanism for common germline variants in miRNA binding sites leading to cancer development. 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 2927. doi:1538-7445.AM2012-2927

Abstract Purpose: Pancreatic adenocarcinoma is among the most lethal forms of all malignancies. While significant research has implicated multiple genes in disease pathogenesis, identification of therapeutic leads has been difficult. The majority of currently available therapies provide only marginal benefit. To address this issue, our goal was to assess global somatic events across multiple pancreatic tumors by using next-generation sequencing technologies (NGS). Experimental Design: We performed whole genome sequencing (WGS) on paired tumor and normal samples from 3 pancreatic adenocarcinoma patients. Additionally, we performed whole transcriptome sequencing (WTS) on 2 of these patients’ tumors for which tumor RNA was available. Results: We generated on average 132 billion mappable bases across all patients using WGS, and identified 140 somatic coding events including point mutations, insertion/deletions, and chromosomal copy number variants. Data from WGS and WTS were integrated and analyzed to identify the most heavily affected pancreatic cancer pathways across all patients. The KRAS signaling pathway was identified as the most significant map. Unbiased global pathway analysis was also performed to identify processes that may be impacted. Consistent with the previous analysis, the most significant pathway from this analysis is also KRAS signaling. This integrated analysis specifically identifies heavily affected pathways in pancreatic cancer and pinpoints genes and processes that may be considered as targets for more effective therapies. Conclusions: While sequencing of additional patients is necessary to strengthen selection of key pathways, we demonstrate the utility of using NGS to understand the processes driving pancreatic adenocarcinoma so that standardized treatments for patients may be improved. 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 5104. doi:1538-7445.AM2012-5104

Abstract Background. Autophagy maintains energy metabolism and enables tumor cell survival during cellular stress. Beclin 1 is involved in autophagosome formation during which membrane-bound LC3 binds to p62/sequestosome 1. We reported that ectopic p62 can promote 5-fluorouracil-induced cell death (JBC 286:40002-12, 2011). It is unknown whether these autophagy markers can predict clinical outcome. Materials & Methods. Beclin 1, LC3B, and p62 proteins were analyzed by immunohistochemistry in tissue microarrays containing stage II and III colon carcinomas (N=151) from patients receiving 5-FU or surgery alone. The mean weighted score (intensity x extent) was dichotomized for survival analysis using Cox models. Results. Increased levels of cytoplasmic LC3, Beclin 1, and p62 expression were detected in 21%, 14%, and 68% of cancers, respectively, and were not associated with tumor site, stage or grade. By multivariable analysis, increased LC3 was significantly associated with worse disease-free survival (DFS) (p=0.022), but not overall survival (OS) [p=0.069](Table). Beclin 1 was not prognostic, but increased p62 was significantly associated with improved DFS (p=0.038), but not OS (p= 0.176)[Table]. When restricted to 5-FU-treated patients, increased LC3 was associated with worse DFS [HR= 2.38 (1.19, 4.75, p=0.020), elevated Beclin 1 predicted worse OS [HR=2.83 (1.27, 6.34), p = 0.019], and p62 was not prognostic. Conclusion. LC3 and Beclin 1 overexpression were independently associated with adverse outcome. Increased p62 was associated with better outcome consistent with its ability to promote cell death. These data support strategies to inhibit autophagy for therapeutic advantage. 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 4112. doi:1538-7445.AM2012-4112