AZGP1 deficiency promotes angiogenesis in prostate cancer.

BACKGROUND: Loss of AZGP1 expression is a biomarker associated with progression to castration resistance, development of metastasis, and poor disease-specific survival in prostate cancer. However, high expression of AZGP1 cells in prostate cancer has been reported to increase proliferation and invasion. The exact role of AZGP1 in prostate cancer progression remains elusive. METHOD: AZGP1 knockout and overexpressing prostate cancer cells were generated using a lentiviral system. The effects of AZGP1 under- or over-expression in prostate cancer cells were evaluated by in vitro cell proliferation, migration, and invasion assays. Heterozygous AZGP1± mice were obtained from European Mouse Mutant Archive (EMMA), and prostate tissues from homozygous knockout male mice were collected at 2, 6 and 10 months for histological analysis. In vivo xenografts generated from AZGP1 under- or over-expressing prostate cancer cells were used to determine the role of AZGP1 in prostate cancer tumor growth, and subsequent proteomics analysis was conducted to elucidate the mechanisms of AZGP1 action in prostate cancer progression. AZGP1 expression and microvessel density were measured in human prostate cancer samples on a tissue microarray of 215 independent patient samples. RESULT: Neither the knockout nor overexpression of AZGP1 exhibited significant effects on prostate cancer cell proliferation, clonal growth, migration, or invasion in vitro. The prostates of AZGP1-/- mice initially appeared to have grossly normal morphology; however, we observed fibrosis in the periglandular stroma and higher blood vessel density in the mouse prostate by 6 months. In PC3 and DU145 mouse xenografts, over-expression of AZGP1 did not affect tumor growth. Instead, these tumors displayed decreased microvessel density compared to xenografts derived from PC3 and DU145 control cells, suggesting that AZGP1 functions to inhibit angiogenesis in prostate cancer. Proteomics profiling further indicated that, compared to control xenografts, AZGP1 overexpressing PC3 xenografts are enriched with angiogenesis pathway proteins, including YWHAZ, EPHA2, SERPINE1, and PDCD6, MMP9, GPX1, HSPB1, COL18A1, RNH1, and ANXA1. In vitro functional studies show that AZGP1 inhibits human umbilical vein endothelial cell proliferation, migration, tubular formation and branching. Additionally, tumor microarray analysis shows that AZGP1 expression is negatively correlated with blood vessel density in human prostate cancer tissues. CONCLUSION: AZGP1 is a negative regulator of angiogenesis, such that loss of AZGP1 promotes angiogenesis in prostate cancer. AZGP1 likely exerts heterotypical effects on cells in the tumor microenvironment, such as stromal and endothelial cells. This study sheds light on the anti-angiogenic characteristics of AZGP1 in the prostate and provides a rationale to target AZGP1 to inhibit prostate cancer progression.

The retinoic acid receptor co-factor NRIP1 is uniquely upregulated and represents a therapeutic target in acute myeloid leukemia with chromosome 3q rearrangements.

Aberrant expression of Ecotropic Viral Integration Site 1 (EVI1) is a hallmark of acute myeloid leukemia (AML) with inv(3) or t(3;3), which is a disease subtype with especially poor outcome. In studying transcriptomes from AML patients with chromosome 3q rearrangements, we identified a significant upregulation of the Nuclear Receptor Interacting Protein 1 (NRIP1) as well as its adjacent non-coding RNA LOC101927745. Utilizing transcriptomic and epigenomic data from over 900 primary samples from patients as well as genetic and transcriptional engineering approaches, we have identified several mechanisms that can lead to upregulation of NRIP1 in AML. We hypothesize that the LOC101927745 transcription start site harbors a context-dependent enhancer that is bound by EVI1, causing upregulation of NRIP1 in AML with chromosome 3 abnormalities. Furthermore, we showed that NRIP1 knockdown negatively affects the proliferation and survival of 3qrearranged AML cells and increases their sensitivity to all-trans retinoic acid, suggesting that NRIP1 is relevant for the pathogenesis of inv(3)/t(3;3) AML and could serve as a novel therapeutic target in myeloid malignancies with 3q abnormalities.

XPNPEP3 is a novel transcriptional target of canonical Wnt/ß-catenin signaling.

Canonical Wnt/ß-catenin signaling plays important roles in embryonic development and adult tissue regeneration while aberrant Wnt activation is the major driver of sporadic colorectal cancer (CRC). Thus, it is important to characterize the complete ß-catenin target transcriptome. We previously performed microarray-based mRNA profiling of rectal cancer samples stratified for Wnt status. In addition to AXIN2 and EPHB2, XPNPEP3 transcripts were significantly elevated in tumors exhibiting activated Wnt/ß-catenin signaling, validated by Q-PCR. Three different cell lines supported elevated XPNPEP3 transcript levels upon activation of Wnt signaling, confirmed using promoter-luciferase assays. Ectopic expression of XPNPEP3 promoted tumorigenic properties in CRC cells. Immunohistochemistry on a CRC tissue microarray revealed significant correlation between ß-catenin nuclear localization and XPNPEP3 levels. More importantly, XPNPEP3 expression was upregulated compared to normal samples in published expression data sets from several cancers including CRC. Finally, XPNPEP3 expression correlated with poor survival in many cancers. Our results therefore suggest XPNPEP3 to be a transcriptional target of Wnt/ß-catenin pathway with particular significance for CRC.

Loss of Expression of AZGP1 Is Associated With Worse Clinical Outcomes in a Multi-Institutional Radical Prostatectomy Cohort.

BACKGROUND: Given the uncertainties inherent in clinical measures of prostate cancer aggressiveness, clinically validated tissue biomarkers are needed. We tested whether Alpha-2-Glycoprotein 1, Zinc-Binding (AZGP1) protein levels, measured by immunohistochemistry, and RNA expression, by RNA in situ hybridization (RISH), predict recurrence after radical prostatectomy independent of clinical and pathological parameters. METHODS: AZGP1 IHC and RISH were performed on a large multi-institutional tissue microarray resource including 1,275 men with 5 year median follow-up. The relationship between IHC and RISH expression levels was assessed using the Kappa analysis. Associations with clinical and pathological parameters were tested by the Chi-square test and the Wilcoxon rank sum test. Relationships with outcome were assessed with univariable and multivariable Cox proportional hazards models and the Log-rank test. RESULTS: Absent or weak expression of AZGP1 protein was associated with worse recurrence free survival (RFS), disease specific survival, and overall survival after radical prostatectomy in univariable analysis. AZGP1 protein expression, along with pre-operative serum PSA levels, surgical margin status, seminal vesicle invasion, extracapsular extension, and Gleason score predicted RFS on multivariable analysis. Similarly, absent or low AZGP1 RNA expression by RISH predicted worse RFS after prostatectomy in univariable and multivariable analysis. CONCLUSIONS: In our large, rigorously designed validation cohort, loss of AZGP1 expression predicts RFS after radical prostatectomy independent of clinical and pathological variables. Prostate 76:1409-1419, 2016. © 2016 Wiley Periodicals, Inc.

Breakpoint analysis of transcriptional and genomic profiles uncovers novel gene fusions spanning multiple human cancer types.

Gene fusions, like BCR/ABL1 in chronic myelogenous leukemia, have long been recognized in hematologic and mesenchymal malignancies. The recent finding of gene fusions in prostate and lung cancers has motivated the search for pathogenic gene fusions in other malignancies. Here, we developed a "breakpoint analysis" pipeline to discover candidate gene fusions by tell-tale transcript level or genomic DNA copy number transitions occurring within genes. Mining data from 974 diverse cancer samples, we identified 198 candidate fusions involving annotated cancer genes. From these, we validated and further characterized novel gene fusions involving ROS1 tyrosine kinase in angiosarcoma (CEP85L/ROS1), SLC1A2 glutamate transporter in colon cancer (APIP/SLC1A2), RAF1 kinase in pancreatic cancer (ATG7/RAF1) and anaplastic astrocytoma (BCL6/RAF1), EWSR1 in melanoma (EWSR1/CREM), CDK6 kinase in T-cell acute lymphoblastic leukemia (FAM133B/CDK6), and CLTC in breast cancer (CLTC/VMP1). Notably, while these fusions involved known cancer genes, all occurred with novel fusion partners and in previously unreported cancer types. Moreover, several constituted druggable targets (including kinases), with therapeutic implications for their respective malignancies. Lastly, breakpoint analysis identified new cell line models for known rearrangements, including EGFRvIII and FIP1L1/PDGFRA. Taken together, we provide a robust approach for gene fusion discovery, and our results highlight a more widespread role of fusion genes in cancer pathogenesis.

Ameloblastoma: a clinical review and trends in management.

Ameloblastoma is a rare odontogenic neoplasm of the mandible and maxilla, with multiple histologic variants, and high recurrence rates if improperly treated. The current mainstay of treatment is wide local excision with appropriate margins and immediate reconstruction. Here we review the ameloblastoma literature, using the available evidence to highlight the change in management over the past several decades. In addition, we explore the recent molecular characterization of these tumors which may point towards new potential avenues of personalized treatment.

Convergent structural alterations define SWItch/Sucrose NonFermentable (SWI/SNF) chromatin remodeler as a central tumor suppressive complex in pancreatic cancer.

Defining the molecular genetic alterations underlying pancreatic cancer may provide unique therapeutic insight for this deadly disease. Toward this goal, we report here an integrative DNA microarray and sequencing-based analysis of pancreatic cancer genomes. Notable among the alterations newly identified, genomic deletions, mutations, and rearrangements recurrently targeted genes encoding components of the SWItch/Sucrose NonFermentable (SWI/SNF) chromatin remodeling complex, including all three putative DNA binding subunits (ARID1A, ARID1B, and PBRM1) and both enzymatic subunits (SMARCA2 and SMARCA4). Whereas alterations of each individual SWI/SNF subunit occurred at modest-frequency, as mutational "hills" in the genomic landscape, together they affected at least one-third of all pancreatic cancers, defining SWI/SNF as a major mutational "mountain." Consistent with a tumor-suppressive role, re-expression of SMARCA4 in SMARCA4-deficient pancreatic cancer cell lines reduced cell growth and promoted senescence, whereas its overexpression in a SWI/SNF-intact line had no such effect. In addition, expression profiling analyses revealed that SWI/SNF likely antagonizes Polycomb repressive complex 2, implicating this as one possible mechanism of tumor suppression. Our findings reveal SWI/SNF to be a central tumor suppressive complex in pancreatic 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.

Spatial transcriptomics identifies candidate stromal drivers of benign prostatic hyperplasia.

Benign prostatic hyperplasia (BPH) is the nodular proliferation of the prostate transition zone in older men, leading to urinary storage and voiding problems that can be recalcitrant to therapy. Decades ago, John McNeal proposed that BPH originates with the "reawakening" of embryonic inductive activity by adult prostate stroma, which spurs new ductal proliferation and branching morphogenesis. Here, by laser microdissection and transcriptional profiling of the BPH stroma adjacent to hyperplastic branching ducts, we identified secreted factors likely mediating stromal induction of prostate glandular epithelium and coinciding processes. The top stromal factors were insulin-like growth factor 1 (IGF1) and CXC chemokine ligand 13 (CXCL13), which we verified by RNA in situ hybridization to be coexpressed in BPH fibroblasts, along with their cognate receptors (IGF1R and CXCR5) on adjacent epithelium. In contrast, IGF1 but not CXCL13 was expressed in human embryonic prostate stroma. Finally, we demonstrated that IGF1 is necessary for the generation of BPH-1 cell spheroids and patient-derived BPH cell organoids in 3D culture. Our findings partially support historic speculations on the etiology of BPH and provide what we believe to be new molecular targets for rational therapies directed against the underlying processes driving BPH.

Integrative genomic and functional profiling of the pancreatic cancer genome.

BACKGROUND: Pancreatic cancer is a deadly disease with a five-year survival of less than 5%. A better understanding of the underlying biology may suggest novel therapeutic targets. Recent surveys of the pancreatic cancer genome have uncovered numerous new alterations; yet systematic functional characterization of candidate cancer genes has lagged behind. To address this challenge, here we have devised a highly-parallel RNA interference-based functional screen to evaluate many genomically-nominated candidate pancreatic cancer genes simultaneously. RESULTS: For 185 candidate pancreatic cancer genes, selected from recurrently altered genomic loci, we performed a pooled shRNA library screen of cell growth/viability across 10 different cell lines. Knockdown-associated effects on cell growth were assessed by enrichment or depletion of shRNA hairpins, by hybridization to barcode microarrays. A novel analytical approach (COrrelated Phenotypes for On-Target Effects; COPOTE) was used to discern probable on-target knockdown, based on identifying different shRNAs targeting the same gene and displaying concordant phenotypes across cell lines. Knockdown data were integrated with genomic architecture and gene-expression profiles, and selected findings validated using individual shRNAs and/or independent siRNAs. The pooled shRNA library design delivered reproducible data. In all, COPOTE analysis identified 52 probable on-target gene-knockdowns. Knockdown of known oncogenes (KRAS, MYC, SMURF1 and CCNE1) and a tumor suppressor (CDKN2A) showed the expected contrasting effects on cell growth. In addition, the screen corroborated purported roles of PLEKHG2 and MED29 as 19q13 amplicon drivers. Most notably, the analysis also revealed novel possible oncogenic functions of nucleoporin NUP153 (ostensibly by modulating TGFß signaling) and Kruppel-like transcription factor KLF5 in pancreatic cancer. CONCLUSIONS: By integrating physical and functional genomic data, we were able to simultaneously evaluate many candidate pancreatic cancer genes. Our findings uncover new facets of pancreatic cancer biology, with possible therapeutic implications. More broadly, our study provides a general strategy for the efficient characterization of candidate genes emerging from cancer genome studies.

Characterizing the physical genome.

The genome of an organism is a dynamic physical entity, comprising genomic DNA bound to many different proteins and organized into chromosomes. A thorough characterization of the physical genome is relevant to our understanding of processes such as the regulation of gene expression, DNA replication and repair, recombination, chromosome segregation, epigenetic inheritance and genomic instability. Methods based on microarrays are beginning to provide a detailed picture of this physical genome, and they complement the genome-wide studies of mRNA expression profiling that have previously been so successful.

A fused lasso latent feature model for analyzing multi-sample aCGH data.

Array-based comparative genomic hybridization (aCGH) enables the measurement of DNA copy number across thousands of locations in a genome. The main goals of analyzing aCGH data are to identify the regions of copy number variation (CNV) and to quantify the amount of CNV. Although there are many methods for analyzing single-sample aCGH data, the analysis of multi-sample aCGH data is a relatively new area of research. Further, many of the current approaches for analyzing multi-sample aCGH data do not appropriately utilize the additional information present in the multiple samples. We propose a procedure called the Fused Lasso Latent Feature Model (FLLat) that provides a statistical framework for modeling multi-sample aCGH data and identifying regions of CNV. The procedure involves modeling each sample of aCGH data as a weighted sum of a fixed number of features. Regions of CNV are then identified through an application of the fused lasso penalty to each feature. Some simulation analyses show that FLLat outperforms single-sample methods when the simulated samples share common information. We also propose a method for estimating the false discovery rate. An analysis of an aCGH data set obtained from human breast tumors, focusing on chromosomes 8 and 17, shows that FLLat and Significance Testing of Aberrant Copy number (an alternative, existing approach) identify similar regions of CNV that are consistent with previous findings. However, through the estimated features and their corresponding weights, FLLat is further able to discern specific relationships between the samples, for example, identifying 3 distinct groups of samples based on their patterns of CNV for chromosome 17.

DR-Integrator: a new analytic tool for integrating DNA copy number and gene expression data.

SUMMARY: DNA copy number alterations (CNA) frequently underlie gene expression changes by increasing or decreasing gene dosage. However, only a subset of genes with altered dosage exhibit concordant changes in gene expression. This subset is likely to be enriched for oncogenes and tumor suppressor genes, and can be identified by integrating these two layers of genome-scale data. We introduce DNA/RNA-Integrator (DR-Integrator), a statistical software tool to perform integrative analyses on paired DNA copy number and gene expression data. DR-Integrator identifies genes with significant correlations between DNA copy number and gene expression, and implements a supervised analysis that captures genes with significant alterations in both DNA copy number and gene expression between two sample classes. AVAILABILITY: DR-Integrator is freely available for non-commercial use from the Pollack Lab at http://pollacklab.stanford.edu/ and can be downloaded as a plug-in application to Microsoft Excel and as a package for the R statistical computing environment. The R package is available under the name 'DRI' at http://cran.r-project.org/. An example analysis using DR-Integrator is included as supplemental material. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

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.

KIT mutations confer a distinct gene expression signature in core binding factor leukaemia.

Core binding factor (CBF) leukaemias, characterized by either inv(16)(p13.1q22) or t(8;21)(q22;q22), constitute acute myeloid leukaemia (AML) subgroups with favourable prognosis. However, 40-50% of patients relapse, emphasizing the need for risk-adapted treatment approaches. In this regard, studying secondary genetic aberrations, such as mutations of the KIT gene, is of great interest, particularly as they can be targeted by receptor tyrosine kinase inhibitors (TKI). However, so far little is known about the biology underlying KIT-mutated CBF leukaemias. We analysed gene expression profiles of 83 CBF AML cases with known KIT mutation status in order to gain novel insights in KIT-mutated CBF pathogenesis. KIT-mutated cases were characterized by deregulation of genes belonging to the NFkB signalling complex suggesting impaired control of apoptosis. Notably, a subgroup of KIT wildtype cases was also characterized by the KIT mutation signature due to yet unknown aberrations. Our data suggest that this CBF leukaemia subgroup might profit from TKI therapy, however, the relevance of the KIT mutation-associated signature remains to be validated prior to clinical implementation. Nevertheless, the existence of such a signature supports the notion of relevant biological differences in CBF leukaemia and might serve as diagnostic tool in the future.

An FLT3 gene-expression signature predicts clinical outcome in normal karyotype AML.

Acute myeloid leukemia with normal karyotype (NK-AML) represents a cytogenetic grouping with intermediate prognosis but substantial molecular and clinical heterogeneity. Within this subgroup, presence of FLT3 (FMS-like tyrosine kinase 3) internal tandem duplication (ITD) mutation predicts less favorable outcome. The goal of our study was to discover gene-expression patterns correlated with FLT3-ITD mutation and to evaluate the utility of a FLT3 signature for prognostication. DNA microarrays were used to profile gene expression in a training set of 65 NK-AML cases, and supervised analysis, using the Prediction Analysis of Microarrays method, was applied to build a gene expression-based predictor of FLT3-ITD mutation status. The optimal predictor, composed of 20 genes, was then evaluated by classifying expression profiles from an independent test set of 72 NK-AML cases. The predictor exhibited modest performance (73% sensitivity; 85% specificity) in classifying FLT3-ITD status. Remarkably, however, the signature outperformed FLT3-ITD mutation status in predicting clinical outcome. The signature may better define clinically relevant FLT3 signaling and/or alternative changes that phenocopy FLT3-ITD, whereas the signature genes provide a starting point to dissect these pathways. Our findings support the potential clinical utility of a gene expression-based measure of FLT3 pathway activation in AML.

MYC stimulates EZH2 expression by repression of its negative regulator miR-26a.

The MYC oncogene, which is commonly mutated/amplified in tumors, represents an important regulator of cell growth because of its ability to induce both proliferation and apoptosis. Recent evidence links MYC to altered miRNA expression, thereby suggesting that MYC-regulated miRNAs might contribute to tumorigenesis. To further analyze the impact of MYC-regulated miRNAs, we investigated a murine lymphoma model harboring the MYC transgene in a Tet-off system to control its expression. Microarray-based miRNA expression profiling revealed both known and novel MYC targets. Among the miRNAs repressed by MYC, we identified the potential tumor suppressor miR-26a, which possessed the ability to attenuate proliferation in MYC-dependent cells. Interestingly, miR-26a was also found to be deregulated in primary human Burkitt lymphoma samples, thereby probably being of clinical relevance. Although today only few miRNA targets have been identified in human disease, we could show that ectopic expression of miR-26a influenced cell cycle progression by targeting the bona fide oncogene EZH2, a Polycomb protein and global regulator of gene expression yet unknown to be regulated by miRNAs. Thus, in addition to directly targeting protein-coding genes, MYC modulates genes important to oncogenesis via deregulation of miRNAs, thereby vitally contributing to MYC-induced lymphomagenesis.

Gene Expression Profiling of Head and Neck Tumors Identifies FOXP1 and SOX10 Expression as Useful for Distinguishing Ameloblastoma From Basaloid Salivary Gland Tumors.

Odontogenic tumors show considerable morphologic heterogeneity and at times the diagnosis can be challenging. Ameloblastoma, the most common odontogenic tumor, can have morphologic similarity to some salivary gland tumors and therefore we sought to identify biomarkers that might aid in the diagnosis by performing transcriptome wide gene expression profiling of 80 odontogenic and salivary gland neoplasms. These data identified the FOXP1/SOX10 expression profile as characteristic of many odontogenic tumors including ameloblastoma but largely absent in salivary gland tumors. We then assessed 173 salivary gland tumors and 108 odontogenic tumors by immunohistochemistry for FOXP1 and SOX10 expression and found that 34/35 (97%) cases of ameloblastomas were diffusely positive for FOXP1 but completely negative for SOX10. None of the basaloid salivary neoplasms (basal cell adenoma, adenoid cystic carcinoma, polymorphous adenocarcinoma, and myoepitheloma) demonstrated FOXP1/SOX10 expression pattern. Taken together, the results of this study suggest that the FOXP1/SOX10 immunophenotype is common in odontogenic tumors including ameloblastoma and might be useful distinguishing these from similar appearing basaloid salivary gland tumors.