Kevin A. Kwei

Molecular Profiling of Breast Cancer Cell Lines Defines Relevant Tumor Models and Provides a Resource for Cancer Gene Discovery

Background Breast cancer cell lines have been used widely to investigate breast cancer pathobiology and new therapies. Breast cancer is a molecularly heterogeneous disease, and it is important to understand how well and which cell lines best model that diversity. In particular, microarray studies have identified molecular subtypes–luminal A, luminal B, ERBB2-associated, basal-like and normal-like–with characteristic gene-expression patterns and underlying DNA copy number alterations (CNAs). Here, we studied a collection of breast cancer cell lines to catalog molecular profiles and to assess their relation to breast cancer subtypes. Methods Whole-genome DNA microarrays were used to profile gene expression and CNAs in a collection of 52 widely-used breast cancer cell lines, and comparisons were made to existing profiles of primary breast tumors. Hierarchical clustering was used to identify gene-expression subtypes, and Gene Set Enrichment Analysis (GSEA) to discover biological features of those subtypes. Genomic and transcriptional profiles were integrated to discover within high-amplitude CNAs candidate cancer genes with coordinately altered gene copy number and expression. Findings Transcriptional profiling of breast cancer cell lines identified one luminal and two basal-like (A and B) subtypes. Luminal lines displayed an estrogen receptor (ER) signature and resembled luminal-A/B tumors, basal-A lines were associated with ETS-pathway and BRCA1 signatures and resembled basal-like tumors, and basal-B lines displayed mesenchymal and stem/progenitor-cell characteristics. Compared to tumors, cell lines exhibited similar patterns of CNA, but an overall higher complexity of CNA (genetically simple luminal-A tumors were not represented), and only partial conservation of subtype-specific CNAs. We identified 80 high-level DNA amplifications and 13 multi-copy deletions, and the resident genes with concomitantly altered gene-expression, highlighting known and novel candidate breast cancer genes. Conclusions Overall, breast cancer cell lines were genetically more complex than tumors, but retained expression patterns with relevance to the luminal-basal subtype distinction. The compendium of molecular profiles defines cell lines suitable for investigations of subtype-specific pathobiology, cancer stem cell biology, biomarkers and therapies, and provides a resource for discovery of new breast cancer genes.

Genomic profiling identifies TITF1 as a lineage-specific oncogene amplified in lung cancer

Lung cancer is a leading cause of cancer death, where the amplification of oncogenes contributes to tumorigenesis. Genomic profiling of 128 lung cancer cell lines and tumors revealed frequent focal DNA amplification at cytoband 14q13.3, a locus not amplified in other tumor types. The smallest region of recurrent amplification spanned the homeobox transcription factor TITF1 (thyroid transcription factor 1; also called NKX2-1), previously linked to normal lung development and function. When amplified, TITF1 exhibited increased expression at both the RNA and protein levels. Small interfering RNA (siRNA)-mediated knockdown of TITF1 in lung cancer cell lines with amplification led to reduced cell proliferation, manifested by both decreased cell-cycle progression and increased apoptosis. Our findings indicate that TITF1 amplification and overexpression contribute to lung cancer cell proliferation rates and survival and implicate TITF1 as a lineage-specific oncogene in lung cancer.

Genomic instability in breast cancer: Pathogenesis and clinical implications

Breast cancer is a heterogeneous disease, appreciable by molecular markers, gene‐expression profiles, and most recently, patterns of genomic alteration. In particular, genomic profiling has revealed three distinct patterns of DNA copy‐number alteration: a “simple” type with few gains or losses of whole chromosome arms, an “amplifier” type with focal high‐level DNA amplifications, and a “complex” type marked by numerous low‐amplitude changes and copy‐number transitions. The three patterns are associated with distinct gene‐expression subtypes, and preferentially target different loci in the genome (implicating distinct cancer genes). Moreover, the different patterns of alteration imply distinct underlying mechanisms of genomic instability. The amplifier pattern may arise from transient telomere dysfunction, although new data suggest ongoing “amplifier” instability. The complex pattern shows similarity to breast cancers with germline BRCA1 mutation, which also exhibit “basal‐like” expression profiles and complex‐pattern genomes, implicating a possible defect in BRCA1‐associated repair of DNA double‐strand breaks. As such, targeting presumptive DNA repair defects represents a promising area of clinical investigation. Future studies should clarify the pathogenesis of breast cancers with amplifier and complex‐pattern genomes, and will likely identify new therapeutic opportunities.

Genomic Profiling Identifies GATA6 as a Candidate Oncogene Amplified in Pancreatobiliary Cancer

Pancreatobiliary cancers have among the highest mortality rates of any cancer type. Discovering the full spectrum of molecular genetic alterations may suggest new avenues for therapy. To catalogue genomic alterations, we carried out array-based genomic profiling of 31 exocrine pancreatic cancers and 6 distal bile duct cancers, expanded as xenografts to enrich the tumor cell fraction. We identified numerous focal DNA amplifications and deletions, including in 19% of pancreatobiliary cases gain at cytoband 18q11.2, a locus uncommonly amplified in other tumor types. The smallest shared amplification at 18q11.2 included GATA6, a transcriptional regulator previously linked to normal pancreas development. When amplified, GATA6 was overexpressed at both the mRNA and protein levels, and strong immunostaining was observed in 25 of 54 (46%) primary pancreatic cancers compared to 0 of 33 normal pancreas specimens surveyed. GATA6 expression in xenografts was associated with specific microarray gene-expression patterns, enriched for GATA binding sites and mitochondrial oxidative phosphorylation activity. siRNA mediated knockdown of GATA6 in pancreatic cancer cell lines with amplification led to reduced cell proliferation, cell cycle progression, and colony formation. Our findings indicate that GATA6 amplification and overexpression contribute to the oncogenic phenotypes of pancreatic cancer cells, and identify GATA6 as a candidate lineage-specific oncogene in pancreatobiliary cancer, with implications for novel treatment strategies.

Genomic and functional analysis identifies CRKL as an oncogene amplified in lung cancer

DNA amplifications, leading to the overexpression of oncogenes, are a cardinal feature of lung cancer and directly contribute to its pathogenesis. To uncover such novel alterations, we performed an array-based comparative genomic hybridization survey of 128 non-small-cell lung cancer cell lines and tumors. Prominent among our findings, we identified recurrent high-level amplification at cytoband 22q11.21 in 3% of lung cancer specimens, with another 11% of specimens exhibiting low-level gain spanning that locus. The 22q11.21 amplicon core contained eight named genes, only four of which were overexpressed (by transcript profiling) when amplified. Among these, CRKL encodes an adapter protein functioning in signal transduction, best known as a substrate of the BCR-ABL kinase in chronic myelogenous leukemia. RNA-interference-mediated knockdown of CRKL in lung cancer cell lines with (but not without) amplification led to significantly decreased cell proliferation, cell-cycle progression, cell survival, and cell motility and invasion. In addition, overexpression of CRKL in immortalized human bronchial epithelial cells led to enhanced growth factor-independent cell growth. Our findings indicate that amplification and resultant overexpression of CRKL contribute to diverse oncogenic phenotypes in lung cancer, with implications for targeted therapy, and highlight a role of adapter proteins as primary genetic drivers of tumorigenesis.

SMURF1 amplification promotes invasiveness in pancreatic cancer.

Pancreatic cancer is a deadly disease, and new therapeutic targets are urgently needed. We previously identified DNA amplification at 7q21-q22 in pancreatic cancer cell lines. Now, by high-resolution genomic profiling of human pancreatic cancer cell lines and human tumors (engrafted in immunodeficient mice to enrich the cancer epithelial fraction), we define a 325 Kb minimal amplicon spanning SMURF1, an E3 ubiquitin ligase and known negative regulator of transforming growth factor β (TGFβ) growth inhibitory signaling. SMURF1 amplification was confirmed in primary human pancreatic cancers by fluorescence in situ hybridization (FISH), where 4 of 95 cases (4.2%) exhibited amplification. By RNA interference (RNAi), knockdown of SMURF1 in a human pancreatic cancer line with focal amplification (AsPC-1) did not alter cell growth, but led to reduced cell invasion and anchorage-independent growth. Interestingly, this effect was not mediated through altered TGFβ signaling, assayed by transcriptional reporter. Finally, overexpression of SMURF1 (but not a catalytic mutant) led to loss of contact inhibition in NIH-3T3 mouse embryo fibroblast cells. Together, these findings identify SMURF1 as an amplified oncogene driving multiple tumorigenic phenotypes in pancreatic cancer, and provide a new druggable target for molecularly directed therapy.

CAMK1D amplification implicated in epithelial–mesenchymal transition in basal-like breast cancer

Breast cancer exhibits clinical and molecular heterogeneity, where expression profiling studies have identified five major molecular subtypes. The basal‐like subtype, expressing basal epithelial markers and negative for estrogen receptor (ER), progesterone receptor (PR) and HER2, is associated with higher overall levels of DNA copy number alteration (CNA), specific CNAs (like gain on chromosome 10p), and poor prognosis. Discovering the molecular genetic basis of tumor subtypes may provide new opportunities for therapy. To identify the driver oncogene on 10p associated with basal‐like tumors, we analyzed genomic profiles of 172 breast carcinomas. The smallest shared region of gain spanned just seven genes at 10p13, including calcium/calmodulin‐dependent protein kinase ID (CAMK1D), functioning in intracellular signaling but not previously linked to cancer. By microarray, CAMK1D was overexpressed when amplified, and by immunohistochemistry exhibited elevated expression in invasive carcinomas compared to carcinoma in situ. Engineered overexpression of CAMK1D in non‐tumorigenic breast epithelial cells led to increased cell proliferation, and molecular and phenotypic alterations indicative of epithelial–mesenchymal transition (EMT), including loss of cell–cell adhesions and increased cell migration and invasion. Our findings identify CAMK1D as a novel amplified oncogene linked to EMT in breast cancer, and as a potential therapeutic target with particular relevance to clinically unfavorable basal‐like tumors.

Modulators of Sensitivity and Resistance to Inhibition of PI3K Identified in a Pharmacogenomic Screen of the NCI-60 Human Tumor Cell Line Collection

The phosphoinositide 3-kinase (PI3K) signaling pathway is significantly altered in a wide variety of human cancers, driving cancer cell growth and survival. Consequently, a large number of PI3K inhibitors are now in clinical development. To begin to improve the selection of patients for treatment with PI3K inhibitors and to identify de novo determinants of patient response, we sought to identify and characterize candidate genomic and phosphoproteomic biomarkers predictive of response to the selective PI3K inhibitor, GDC-0941, using the NCI-60 human tumor cell line collection. In this study, sixty diverse tumor cell lines were exposed to GDC-0941 and classified by GI50 value as sensitive or resistant. The most sensitive and resistant cell lines were analyzed for their baseline levels of gene expression and phosphorylation of key signaling nodes. Phosphorylation or activation status of both the PI3K-Akt signaling axis and PARP were correlated with in vitro response to GDC-0941. A gene expression signature associated with in vitro sensitivity to GDC-0941 was also identified. Furthermore, in vitro siRNA-mediated silencing of two genes in this signature, OGT and DDN, validated their role in modulating sensitivity to GDC-0941 in numerous cell lines and begins to provide biological insights into their role as chemosensitizers. These candidate biomarkers will offer useful tools to begin a more thorough understanding of determinants of patient response to PI3K inhibitors and merit exploration in human cancer patients treated with PI3K inhibitors.

Recurrent rearrangements of the Myb/SANT-like DNA-binding domain containing 3 gene (MSANTD3) in salivary gland acinic cell carcinoma

Pathogenic gene fusions have been identified in several histologic types of salivary gland neoplasia, but not previously in acinic cell carcinoma (AcCC). To discover novel gene fusions, we performed whole-transcriptome sequencing surveys of three AcCC archival cases. In one specimen we identified a novel HTN3-MSANTD3 gene fusion, and in another a novel PRB3-ZNF217 gene fusion. The structure of both fusions was consistent with the promoter of the 5’ partner (HTN3 or PRB3), both highly expressed salivary gland genes, driving overexpression of full-length MSANTD3 or ZNF217. By fluorescence in situ hybridization of an expanded AcCC case series, we observed MSANTD3 rearrangements altogether in 3 of 20 evaluable cases (15%), but found no additional ZNF217 rearrangements. MSANTD3 encodes a previously uncharacterized Myb/SANT domain-containing protein. Immunohistochemical staining demonstrated diffuse nuclear MSANTD3 expression in 8 of 27 AcCC cases (30%), including the three cases with MSANTD3 rearrangement. MSANTD3 displayed heterogeneous expression in normal salivary ductal epithelium, as well as among other histologic types of salivary gland cancer though without evidence of translocation. In a broader survey, MSANTD3 showed variable expression across a wide range of normal and neoplastic human tissue specimens. In preliminary functional studies, engineered MSANTD3 overexpression in rodent salivary gland epithelial cells did not enhance cell proliferation, but led to significant upregulation of gene sets involved in protein synthesis. Our findings newly identify MSANTD3 rearrangement as a recurrent event in salivary gland AcCC, providing new insight into disease pathogenesis, and identifying a putative novel human oncogene.

Integrative Genomics Implicates EGFR as a Downstream Mediator in NKX2-1 Amplified Non-Small Cell Lung Cancer.

NKX2-1, encoding a homeobox transcription factor, is amplified in approximately 15% of non-small cell lung cancers (NSCLC), where it is thought to drive cancer cell proliferation and survival. However, its mechanism of action remains largely unknown. To identify relevant downstream transcriptional targets, here we carried out a combined NKX2-1 transcriptome (NKX2-1 knockdown followed by RNAseq) and cistrome (NKX2-1 binding sites by ChIPseq) analysis in four NKX2-1-amplified human NSCLC cell lines. While NKX2-1 regulated genes differed among the four cell lines assayed, cell proliferation emerged as a common theme. Moreover, in 3 of the 4 cell lines, epidermal growth factor receptor (EGFR) was among the top NKX2-1 upregulated targets, which we confirmed at the protein level by western blot. Interestingly, EGFR knockdown led to upregulation of NKX2-1, suggesting a negative feedback loop. Consistent with this finding, combined knockdown of NKX2-1 and EGFR in NCI-H1819 lung cancer cells reduced cell proliferation (as well as MAP-kinase and PI3-kinase signaling) more than knockdown of either alone. Likewise, NKX2-1 knockdown enhanced the growth-inhibitory effect of the EGFR-inhibitor erlotinib. Taken together, our findings implicate EGFR as a downstream effector of NKX2-1 in NKX2-1 amplified NSCLC, with possible clinical implications, and provide a rich dataset for investigating additional mediators of NKX2-1 driven oncogenesis.