Keyan Salari

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 Reveals Alternative Genetic Pathways of Prostate Tumorigenesis

Prostate cancer is clinically heterogeneous, ranging from indolent to lethal disease. Expression profiling previously defined three subtypes of prostate cancer, one (subtype-1) linked to clinically favorable behavior, and the others (subtypes-2 and -3) linked with a more aggressive form of the disease. To explore disease heterogeneity at the genomic level, we carried out array-based comparative genomic hybridization (array CGH) on 64 prostate tumor specimens, including 55 primary tumors and 9 pelvic lymph node metastases. Unsupervised cluster analysis of DNA copy number alterations (CNA) identified recurrent aberrations, including a 6q15-deletion group associated with subtype-1 gene expression patterns and decreased tumor recurrence. Supervised analysis further disclosed distinct patterns of CNA among gene-expression subtypes, where subtype-1 tumors exhibited characteristic deletions at 5q21 and 6q15, and subtype-2 cases harbored deletions at 8p21 (NKX3-1) and 21q22 (resulting in TMPRSS2-ERG fusion). Lymph node metastases, predominantly subtype-3, displayed overall higher frequencies of CNA, and in particular gains at 8q24 (MYC) and 16p13, and loss at 10q23 (PTEN) and 16q23. Our findings reveal that prostate cancers develop via a limited number of alternative preferred genetic pathways. The resultant molecular genetic subtypes provide a new framework for investigating prostate cancer biology and explain in part the clinical heterogeneity of the disease.

Population stratification confounds genetic association studies among Latinos

In the United States, asthma prevalence and mortality are the highest among Puerto Ricans and the lowest among Mexicans. Case-control association studies are a powerful strategy for identifying genes of modest effect in complex diseases. However, studies of complex disorders in admixed populations such as Latinos may be confounded by population stratification. We used ancestry informative markers (AIMs) to identify and correct for population stratification among Mexican and Puerto Rican subjects participating in case-control studies of asthma. Three hundred and sixty-two subjects with asthma (Mexican: 181, Puerto Rican: 181) and 359 ethnically matched controls (Mexican: 181, Puerto Rican: 178) were genotyped for 44 AIMs. We observed a greater than expected degree of association between pairs of AIMs on different chromosomes in Mexicans (P<0.00001) and Puerto Ricans (P<0.00002) providing evidence for population substructure and/or recent admixture. To assess the effect of population stratification on association studies of asthma, we measured differences in genetic background of cases and controls by comparing allele frequencies of the 44 AIMs. Among Puerto Ricans but not in Mexicans, we observed a significant overall difference in allele frequencies between cases and controls (P=0.0002); of 44 AIMs tested, 8 (18%) were significantly associated with asthma. However, after adjustment for individual ancestry, only two of these markers remained significantly associated with the disease. Our findings suggest that empirical assessment of the effects of stratification is critical to appropriately interpret the results of case-control studies in admixed populations.

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.

Recurrent deletion of CHD1 in prostate cancer with relevance to cell invasiveness

Though prostate cancer is often indolent, it is nonetheless a leading cause of cancer death. Defining the underlying molecular genetic alterations may lead to new strategies for prevention or treatment. Towards this goal, we performed array-based comparative genomic hybridization (CGH) on 86 primary prostate tumors. Among the most frequent alterations not associated with a known cancer gene, we identified focal deletions within 5q21 in 15 out of 86 (17%) cases. By high-resolution tiling array CGH, the smallest common deletion targeted just one gene, the chromatin remodeler chromodomain helicase DNA-binding protein 1 (CHD1). Expression of CHD1 was significantly reduced in tumors with deletion (P=0.03), and compared with normal prostate (P=0.04). Exon sequencing analysis also uncovered nonsynonymous mutations in 1 out of 7 (14%) cell lines (LAPC4) and in 1 out of 24 (4%) prostate tumors surveyed. RNA interference-mediated knockdown of CHD1 in two nontumorigenic prostate epithelial cell lines, OPCN2 and RWPE-1, did not alter cell growth, but promoted cell invasiveness, and in OPCN2-enhanced cell clonogenicity. Taken together, our findings suggest that CHD1 deletion may underlie cell invasiveness in a subset of prostate cancers, and indicate a possible novel role of altered chromatin remodeling in prostate tumorigenesis.

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.

Comparative Analyses of Chromosome Alterations in Soft-Tissue Metastases within and across Patients with Castration-Resistant Prostate Cancer

Androgen deprivation is the mainstay of therapy for progressive prostate cancer. Despite initial and dramatic tumor inhibition, most men eventually fail therapy and die of metastatic castration-resistant (CR) disease. Here, we characterize the profound degree of genomic alteration found in CR tumors using array comparative genomic hybridization (array CGH), gene expression arrays, and fluorescence in situ hybridization (FISH). Bycluster analysis, we show that the similarity of the genomic profiles from primary and metastatic tumors is driven by the patient. Using data adjusted for this similarity, we identify numerous high-frequency alterations in the CR tumors, such as 8p loss and chromosome 7 and 8q gain. By integrating array CGH and expression array data, we reveal genes whose correlated values suggest they are relevant to prostate cancer biology. We find alterations that are significantly associated with the metastases of specific organ sites, and others with CR tumors versus the tumors of patients with localized prostate cancer not treated with androgen deprivation. Within the high-frequency sites of loss in CR metastases, we find an overrepresentation of genes involved in cellular lipid metabolism, including PTEN. Finally, using FISH, we verify the presence of a gene fusion between TMPRSS2 and ERG suggested by chromosome 21 deletions detected by array CGH. We find the fusion in 54% of our CR tumors, and 81% of the fusion-positive tumors contain cells with multiple copies of the fusion. Our investigation lays the foundation for a better understanding of and possible therapeutic targets for CR disease, the poorly responsive and final stage of prostate cancer.

CDX2 is an amplified lineage-survival oncogene in colorectal cancer

The mutational activation of oncogenes drives cancer development and progression. Classic oncogenes, such as MYC and RAS, are active across many different cancer types. In contrast, “lineage-survival” oncogenes represent a distinct and emerging class typically comprising transcriptional regulators of a specific cell lineage that, when deregulated, support the proliferation and survival of cancers derived from that lineage. Here, in a large collection of colorectal cancer cell lines and tumors, we identify recurrent amplification of chromosome 13, an alteration highly restricted to colorectal-derived cancers. A minimal region of amplification on 13q12.2 pinpoints caudal type homeobox transcription factor 2 (CDX2), a regulator of normal intestinal lineage development and differentiation, as a target of the amplification. In contrast to its described role as a colorectal tumor suppressor, CDX2 when amplified is required for the proliferation and survival of colorectal cancer cells. Further, transcriptional profiling, binding-site analysis, and functional studies link CDX2 to Wnt/β-catenin signaling, itself a key oncogenic pathway in colorectal cancer. These data characterize CDX2 as a lineage-survival oncogene deregulated in colorectal cancer. Our findings challenge a prevailing view that CDX2 is a tumor suppressor in colorectal cancer and uncover an additional piece in the multistep model of colorectal tumorigenesis.