Concordance between an FDA-approved companion diagnostic and an alternative assay kit for assessing homologous recombination deficiency in ovarian cancer.

OBJECTIVE: Authors evaluated the performance of a commercially available next-generation sequencing assay kit; this was based on genomic content from Illumina's TruSight™ Oncology 500 research assay that identifies BRCA variants and proprietary algorithms licensed from Myriad and, with additional genomic content, measures the homologous recombination deficiency (HRD) genomic instability score (GIS) in tumor tissue (TSO 500 HRD assay). METHODS: Data from the TSO 500 HRD assay were compared with data from the Myriad MyChoice®CDx PLUS assay (Myriad assay). Prevalence rates for overall HRD status and BRCA mutations (a deleterious or suspected deleterious BRCA1 or BRCA2 mutation or both) and assay agreement rates for HRD GIS and BRCA analysis were assessed in ovarian tumor samples. Pearson correlations of the continuous HRD GIS and analytic sensitivity and specificity were evaluated. RESULTS: The prevalence of overall HRD positivity was 51.2% (TSO 500 HRD assay) versus 49.2% (Myriad assay) and the prevalence of BRCA mutations was 27.6% (TSO 500 HRD assay) versus 25.5% (Myriad assay). After post-processing optimization, concordance of the HRD GIS was 0.980 in all samples and 0.976 in the non-BRCA mutation cohort; the area under the receiver operating characteristic curve was 0.995 and 0.992, respectively. CONCLUSIONS: Comparison between the Illumina and Myriad assays showed that overall HRD status, the individual components of BRCA analysis, and HRD GIS detection results were highly concordant (>93%), suggesting the TSO 500 HRD assay will approach the analytical accuracy of the FDA-approved Myriad assay.

The interplay between cancer type, panel size and tumor mutational burden threshold in patient selection for cancer immunotherapy.

The tumor mutational burden (TMB) is increasingly recognized as an emerging biomarker that predicts improved outcomes or response to immune checkpoint inhibitors in cancer. A multitude of technical and biological factors make it difficult to compare TMB values across platforms, histologies, and treatments. Here, we present a mechanistic model that explains the association between panel size, histology, and TMB threshold with panel performance and survival outcome and demonstrate the limitations of existing methods utilized to harmonize TMB across platforms.

Author Correction: Interaction between the microbiome and TP53 in human lung cancer.

Following publication of the original paper [1], the authors submitted a new Additional file 5 to replace the one containing formatting issues. The updated Additional file 5 is published in this correction.

Interaction between the microbiome and TP53 in human lung cancer.

BACKGROUND: Lung cancer is the leading cancer diagnosis worldwide and the number one cause of cancer deaths. Exposure to cigarette smoke, the primary risk factor in lung cancer, reduces epithelial barrier integrity and increases susceptibility to infections. Herein, we hypothesize that somatic mutations together with cigarette smoke generate a dysbiotic microbiota that is associated with lung carcinogenesis. Using lung tissue from 33 controls and 143 cancer cases, we conduct 16S ribosomal RNA (rRNA) bacterial gene sequencing, with RNA-sequencing data from lung cancer cases in The Cancer Genome Atlas serving as the validation cohort. RESULTS: Overall, we demonstrate a lower alpha diversity in normal lung as compared to non-tumor adjacent or tumor tissue. In squamous cell carcinoma specifically, a separate group of taxa are identified, in which Acidovorax is enriched in smokers. Acidovorax temporans is identified within tumor sections by fluorescent in situ hybridization and confirmed by two separate 16S rRNA strategies. Further, these taxa, including Acidovorax, exhibit higher abundance among the subset of squamous cell carcinoma cases with TP53 mutations, an association not seen in adenocarcinomas. CONCLUSIONS: The results of this comprehensive study show both microbiome-gene and microbiome-exposure interactions in squamous cell carcinoma lung cancer tissue. Specifically, tumors harboring TP53 mutations, which can impair epithelial function, have a unique bacterial consortium that is higher in relative abundance in smoking-associated tumors of this type. Given the significant need for clinical diagnostic tools in lung cancer, this study may provide novel biomarkers for early detection.

A Genome-Wide Screen Reveals a Role for the HIR Histone Chaperone Complex in Preventing Mislocalization of Budding Yeast CENP-A.

Centromeric localization of the evolutionarily conserved centromere-specific histone H3 variant CENP-A (Cse4 in yeast) is essential for faithful chromosome segregation. Overexpression and mislocalization of CENP-A lead to chromosome segregation defects in yeast, flies, and human cells. Overexpression of CENP-A has been observed in human cancers; however, the molecular mechanisms preventing CENP-A mislocalization are not fully understood. Here, we used a genome-wide synthetic genetic array (SGA) to identify gene deletions that exhibit synthetic dosage lethality (SDL) when Cse4 is overexpressed. Deletion for genes encoding the replication-independent histone chaperone HIR complex (HIR1, HIR2, HIR3, HPC2) and a Cse4-specific E3 ubiquitin ligase, PSH1, showed highest SDL. We defined a role for Hir2 in proteolysis of Cse4 that prevents mislocalization of Cse4 to noncentromeric regions for genome stability. Hir2 interacts with Cse4 in vivo, and hir2∆ strains exhibit defects in Cse4 proteolysis and stabilization of chromatin-bound Cse4 Mislocalization of Cse4 to noncentromeric regions with a preferential enrichment at promoter regions was observed in hir2∆ strains. We determined that Hir2 facilitates the interaction of Cse4 with Psh1, and that defects in Psh1-mediated proteolysis contribute to increased Cse4 stability and mislocalization of Cse4 in the hir2∆ strain. In summary, our genome-wide screen provides insights into pathways that regulate proteolysis of Cse4 and defines a novel role for the HIR complex in preventing mislocalization of Cse4 by facilitating proteolysis of Cse4, thereby promoting genome stability.

Somatic mutations in murine models of leukemia and lymphoma: Disease specificity and clinical relevance.

Malignant transformation is a multistep process that is dictated by the acquisition of multiple genomic aberrations that provide growth and survival advantage. During the post genomic era, high throughput genomic sequencing has advanced exponentially, leading to identification of countless cancer associated mutations with potential for targeted therapy. Mouse models of cancer serve as excellent tools to examine the functionality of gene mutations and their contribution to the malignant process. However, it remains unclear whether the genetic events that occur during transformation are similar in mice and humans. To address that, we chose several transgenic mouse models of hematopoietic malignancies and identified acquired mutations in these mice by means of targeted re-sequencing of known cancer-associated genes as well as whole exome sequencing. We found that mutations that are typically found in acute myeloid leukemia or T cell acute lymphoblastic leukemia patients are also common in mouse models of the respective disease. Moreover, we found that the most frequent mutations found in a mouse model of lymphoma occur in a set of epigenetic modifier genes, implicating this pathway in the generation of lymphoma. These results demonstrate that genetically engineered mouse models (GEMM) mimic the genetic evolution of human cancer and serve as excellent platforms for target discovery and validation.

Progenitor B-1 B-cell acute lymphoblastic leukemia is associated with collaborative mutations in 3 critical pathways.

B-1 and B-2 lymphocytes are derived from distinct developmental pathways and represent layered arms of the innate and adaptive immune systems, respectively. In contrast to a majority of murine B-cell malignancies, which stain positive with the B220 antibody, we discovered a novel form of B-cell leukemia in NUP98-PHF23 (NP23) transgenic mice. The immunophenotype (Lin- B220- CD19+ AA4.1+) was identical to that of progenitor (pro) B-1 cells, and VH gene usage was skewed toward 3' V regions, similar to murine fetal liver B cells. Moreover, the gene expression profile of these leukemias was most similar to that of fetal liver pro-B fraction BC, a known source of B-1 B cells, further supporting a pro-B-1 origin of these leukemias. The NP23 pro-B-1 acute lymphoblastic leukemias (ALLs) acquired spontaneous mutations in both Bcor and Janus kinase (Jak) pathway (Jak1/2/3 and Stat5a) genes, supporting a hypothesis that mutations in 3 critical pathways (stem-cell self-renewal, B-cell differentiation, and cytokine signaling) collaborate to induce B-cell precursor (BCP) ALL. Finally, the thymic stromal lymphopoietin (Tslp) cytokine is required for murine B-1 development, and chromosomal rearrangements resulting in overexpression of the TSLP receptor (CRLF2) are present in some patients with high-risk BCP-ALL (referred to as CRLF2r ALL). Gene expression profiles of NP23 pro-B-1 ALL were more similar to that of CRLF2r ALL than non-CRLF2r ALL, and analysis of VH gene usage from patients with CRLF2r ALL demonstrated preferential usage of VH regions used by human B-1 B cells, leading to the suggestion that this subset of patients with BCP-ALL has a malignancy of B-1, rather than B-2, B-cell origin.

SLFN11 Is a Transcriptional Target of EWS-FLI1 and a Determinant of Drug Response in Ewing Sarcoma.

PURPOSE: SLFN11 was identified as a critical determinant of response to DNA-targeted therapies by analyzing gene expression and drug sensitivity of NCI-60 and CCLE datasets. However, how SLFN11 is regulated in cancer cells remained unknown. Ewing sarcoma, which is characterized by the chimeric transcription factor EWS-FLI1, has notably high SLFN11 expression, leading us to investigate whether EWS-FLI1 drives SLFN11 expression and the role of SLFN11 in the drug response of Ewing sarcoma cells. EXPERIMENTAL DESIGN: Binding sites of EWS-FLI1 on the SLFN11 promoter were analyzed by chromatin immunoprecipitation sequencing and promoter-luciferase reporter analyses. The relationship between SLFN11 and EWS-FLI1 were further examined in EWS-FLI1-knockdown or -overexpressing cells and in clinical tumor samples. RESULTS: EWS-FLI1 binds near the transcription start site of SLFN11 promoter and acts as a positive regulator of SLFN11 expression in Ewing sarcoma cells. EWS-FLI1-mediated SLFN11 expression is responsible for high sensitivity of Ewing sarcoma to camptothecin and combinations of PARP inhibitors with temozolomide. Importantly, Ewing sarcoma patients with higher SLFN11 expression showed better tumor-free survival rate. The correlated expression between SLFN11 and FLI1 extends to leukemia, pediatric, colon, breast, and prostate cancers. In addition, expression of other ETS members correlates with SLFN11 in NCI-60 and CCLE datasets, and molecular experiments demonstrate that ETS1 acts as a positive regulator for SLFN11 expression in breast cancer cells. CONCLUSIONS: Our results imply the emerging relevance of SLFN11 as an ETS transcription factor response gene and for therapeutic response to topoisomerase I inhibitors and temozolomide-PARP inhibitor combinations in ETS-activated cancers.

The CDX1-microRNA-215 axis regulates colorectal cancer stem cell differentiation.

The transcription factor caudal-type homeobox 1 (CDX1) is a key regulator of differentiation in the normal colon and in colorectal cancer (CRC). CDX1 activates the expression of enterocyte genes, but it is not clear how the concomitant silencing of stem cell genes is achieved. MicroRNAs (miRNAs) are important mediators of gene repression and have been implicated in tumor suppression and carcinogenesis, but the roles of miRNAs in differentiation, particularly in CRC, remain poorly understood. Here, we identified microRNA-215 (miR-215) as a direct transcriptional target of CDX1 by using high-throughput small RNA sequencing to profile miRNA expression in two pairs of CRC cell lines: CDX1-low HCT116 and HCT116 with stable CDX1 overexpression, and CDX1-high LS174T and LS174T with stable CDX1 knockdown. Validation of candidate miRNAs identified by RNA-seq in a larger cell-line panel revealed miR-215 to be most significantly correlated with CDX1 expression. Quantitative ChIP-PCR and promoter luciferase assays confirmed that CDX1 directly activates miR-215 transcription. miR-215 expression is depleted in FACS-enriched cancer stem cells compared with unsorted samples. Overexpression of miR-215 in poorly differentiated cell lines causes a decrease in clonogenicity, whereas miR-215 knockdown increases clonogenicity and impairs differentiation in CDX1-high cell lines. We identified the genome-wide targets of miR-215 and found that miR-215 mediates the repression of cell cycle and stemness genes downstream of CDX1. In particular, the miR-215 target gene BMI1 has been shown to promote stemness and self-renewal and to vary inversely with CDX1. Our work situates miR-215 as a link between CDX1 expression and BMI1 repression that governs differentiation in CRC.

EWS-FLI1 employs an E2F switch to drive target gene expression.

Cell cycle progression is orchestrated by E2F factors. We previously reported that in ETS-driven cancers of the bone and prostate, activating E2F3 cooperates with ETS on target promoters. The mechanism of target co-regulation remained unknown. Using RNAi and time-resolved chromatin-immunoprecipitation in Ewing sarcoma we report replacement of E2F3/pRB by constitutively expressed repressive E2F4/p130 complexes on target genes upon EWS-FLI1 modulation. Using mathematical modeling we interrogated four alternative explanatory models for the observed EWS-FLI1/E2F3 cooperation based on longitudinal E2F target and regulating transcription factor expression analysis. Bayesian model selection revealed the formation of a synergistic complex between EWS-FLI1 and E2F3 as the by far most likely mechanism explaining the observed kinetics of E2F target induction. Consequently we propose that aberrant cell cycle activation in Ewing sarcoma is due to the de-repression of E2F targets as a consequence of transcriptional induction and physical recruitment of E2F3 by EWS-FLI1 replacing E2F4 on their target promoters.

Identification of RECQ1-regulated transcriptome uncovers a role of RECQ1 in regulation of cancer cell migration and invasion.

The RECQ protein family of helicases has critical roles in protecting and stabilizing the genome. Three of the 5 known members of the human RecQ family are genetically linked with cancer susceptibility syndromes, but the association of the most abundant human RecQ homolog, RECQ1, with cellular transformation is yet unclear. RECQ1 is overexpressed in a variety of human cancers, indicating oncogenic functions. Here, we assessed genome-wide changes in gene expression upon knockdown of RECQ1 in HeLa and MDA-MB-231 cells. Pathway analysis suggested that RECQ1 enhances the expression of multiple genes that play key roles in cell migration, invasion, and metastasis, including EZR, ITGA2, ITGA3, ITGB4, SMAD3, and TGFBR2. Consistent with these results, silencing RECQ1 significantly reduced cell migration and invasion. In comparison to genome-wide annotated promoter regions, the promoters of genes downregulated upon RECQ1 silencing were significantly enriched for a potential G4 DNA forming sequence motif. Chromatin immunoprecipitation assays demonstrated binding of RECQ1 to the G4 motifs in the promoters of select genes downregulated upon RECQ1 silencing. In breast cancer patients, the expression of a subset of RECQ1-activated genes positively correlated with RECQ1 expression. Moreover, high RECQ1 expression was associated with poor prognosis in breast cancer. Collectively, our findings identify a novel function of RECQ1 in gene regulation and indicate that RECQ1 contributes to tumor development and progression, in part, by regulating the expression of key genes that promote cancer cell migration, invasion and metastasis.

Replicon: a software to accurately predict DNA replication timing in metazoan cells.

Eukaryotic DNA replication follows a strict temporal program where genomic loci are replicated at precise times during the S phase of the cell cycle. Yet, the mechanism in control of the timing program in metazoan cells is poorly understood. In a recent publication, the authors proposed an intuitive stochastic model of DNA replication and showed that it predicts replication timing with an accuracy approaching the level of experimental biological repeats. Here, we discuss an extended software implementation of the mechanistic model: Replicon. This package allows interested researchers to predict the global replication timing program in human cells from chromatin data.

A specific missense mutation in GTF2I occurs at high frequency in thymic epithelial tumors.

We analyzed 28 thymic epithelial tumors (TETs) using next-generation sequencing and identified a missense mutation (chromosome 7 c.74146970T>A) in GTF2I at high frequency in type A thymomas, a relatively indolent subtype. In a series of 274 TETs, we detected the GTF2I mutation in 82% of type A and 74% of type AB thymomas but rarely in the aggressive subtypes, where recurrent mutations of known cancer genes have been identified. Therefore, GTF2I mutation correlated with better survival. GTF2I ß and δ isoforms were expressed in TETs, and both mutant isoforms were able to stimulate cell proliferation in vitro. Thymic carcinomas carried a higher number of mutations than thymomas (average of 43.5 and 18.4, respectively). Notably, we identified recurrent mutations of known cancer genes, including TP53, CYLD, CDKN2A, BAP1 and PBRM1, in thymic carcinomas. These findings will complement the diagnostic assessment of these tumors and also facilitate development of a molecular classification and assessment of prognosis and treatment strategies.

A chromatin structure-based model accurately predicts DNA replication timing in human cells.

The metazoan genome is replicated in precise cell lineage-specific temporal order. However, the mechanism controlling this orchestrated process is poorly understood as no molecular mechanisms have been identified that actively regulate the firing sequence of genome replication. Here, we develop a mechanistic model of genome replication capable of predicting, with accuracy rivaling experimental repeats, observed empirical replication timing program in humans. In our model, replication is initiated in an uncoordinated (time-stochastic) manner at well-defined sites. The model contains, in addition to the choice of the genomic landmark that localizes initiation, only a single adjustable parameter of direct biological relevance: the number of replication forks. We find that DNase-hypersensitive sites are optimal and independent determinants of DNA replication initiation. We demonstrate that the DNA replication timing program in human cells is a robust emergent phenomenon that, by its very nature, does not require a regulatory mechanism determining a proper replication initiation firing sequence.

NUP98-PHF23 is a chromatin-modifying oncoprotein that causes a wide array of leukemias sensitive to inhibition of PHD histone reader function.

In this report, we show that expression of a NUP98-PHF23 (NP23) fusion, associated with acute myeloid leukemia (AML) in humans, leads to myeloid, erythroid, T-cell, and B-cell leukemia in mice. The leukemic and preleukemic tissues display a stem cell-like expression signature, including Hoxa, Hoxb, and Meis1 genes. The PHF23 plant homeodomain (PHD) motif is known to bind to H3K4me3 residues, and chromatin immunoprecipitation experiments demonstrated that the NP23 protein binds to chromatin at a specific subset of H3K4me3 sites, including at Hoxa, Hoxb, and Meis1. Treatment of NP23 cells with disulfiram, which inhibits the binding of PHD motifs to H3K4me3, rapidly and selectively killed NP23-expressing myeloblasts; cell death was preceded by decreased expression of Hoxa, Hoxb, and Meis1. Furthermore, AML driven by a related fusion gene, NUP98-JARID1A (NJL), was also sensitive to disulfiram. Thus, the NP23 mouse provides a platform to evaluate compounds that disrupt binding of oncogenic PHD proteins to H3K4me3.

A p21-ZEB1 complex inhibits epithelial-mesenchymal transition through the microRNA 183-96-182 cluster.

The tumor suppressor p21 acts as a cell cycle inhibitor and has also been shown to regulate gene expression by functioning as a transcription corepressor. Here, we identified p21-regulated microRNAs (miRNAs) by sequencing small RNAs from isogenic p21(+/+) and p21(-/-) cells. Three abundant miRNA clusters, miR-200b-200a-429, miR-200c-141, and miR-183-96-182, were downregulated in p21-deficient cells. Consistent with the known function of the miR-200 family and p21 in inhibition of the epithelial-mesenchymal transition (EMT), we observed EMT upon loss of p21 in multiple model systems. To explore a role of the miR-183-96-182 cluster in EMT, we identified its genome-wide targets and found that miR-183 and miR-96 repressed common targets, including SLUG, ZEB1, ITGB1, and KLF4. Reintroduction of miR-200, miR-183, or miR-96 in p21(-/-) cells inhibited EMT, cell migration, and invasion. Conversely, antagonizing miR-200 and miR-183-96-182 cluster miRNAs in p21(+/+) cells increased invasion and elevated the levels of VIM, ZEB1, and SLUG mRNAs. Furthermore, we found that p21 forms a complex with ZEB1 at the miR-183-96-182 cluster promoter to inhibit transcriptional repression of this cluster by ZEB1, suggesting a reciprocal feedback loop.

Oncogenic ETS fusions deregulate E2F3 target genes in Ewing sarcoma and prostate cancer.

Deregulated E2F transcription factor activity occurs in the vast majority of human tumors and has been solidly implicated in disturbances of cell cycle control, proliferation, and apoptosis. Aberrant E2F regulatory activity is often caused by impairment of control through pRB function, but little is known about the interplay of other oncoproteins with E2F. Here we show that ETS transcription factor fusions resulting from disease driving rearrangements in Ewing sarcoma (ES) and prostate cancer (PC) are one such class of oncoproteins. We performed an integrative study of genome-wide DNA-binding and transcription data in EWSR1/FLI1 expressing ES and TMPRSS2/ERG containing PC cells. Supported by promoter activity and mutation analyses, we demonstrate that a large fraction of E2F3 target genes are synergistically coregulated by these aberrant ETS proteins. We propose that the oncogenic effect of ETS fusion oncoproteins is in part mediated by the disruptive effect of the E2F-ETS interaction on cell cycle control. Additionally, a detailed analysis of the regulatory targets of the characteristic EWSR1/FLI1 fusion in ES identifies two functionally distinct gene sets. While synergistic regulation in concert with E2F in the promoter of target genes has a generally activating effect, EWSR1/FLI1 binding independent of E2F3 is predominantly associated with repressed differentiation genes. Thus, EWSR1/FLI1 appears to promote oncogenesis by simultaneously promoting cell proliferation and perturbing differentiation.

The exomes of the NCI-60 panel: a genomic resource for cancer biology and systems pharmacology.

The NCI-60 cell lines are the most frequently studied human tumor cell lines in cancer research. This panel has generated the most extensive cancer pharmacology database worldwide. In addition, these cell lines have been intensely investigated, providing a unique platform for hypothesis-driven research focused on enhancing our understanding of tumor biology. Here, we report a comprehensive analysis of coding variants in the NCI-60 panel of cell lines identified by whole exome sequencing, providing a list of possible cancer specific variants for the community. Furthermore, we identify pharmacogenomic correlations between specific variants in genes such as TP53, BRAF, ERBBs, and ATAD5 and anticancer agents such as nutlin, vemurafenib, erlotinib, and bleomycin showing one of many ways the data could be used to validate and generate novel hypotheses for further investigation. As new cancer genes are identified through large-scale sequencing studies, the data presented here for the NCI-60 will be an invaluable resource for identifying cell lines with mutations in such genes for hypothesis-driven research. To enhance the utility of the data for the greater research community, the genomic variants are freely available in different formats and from multiple sources including the CellMiner and Ingenuity websites.

Genome-wide profiles of CtBP link metabolism with genome stability and epithelial reprogramming in breast cancer.

The C-terminal binding protein (CtBP) is a NADH-dependent transcriptional repressor that links carbohydrate metabolism to epigenetic regulation by recruiting diverse histone-modifying complexes to chromatin. Here global profiling of CtBP in breast cancer cells reveals that it drives epithelial-to-mesenchymal transition, stem cell pathways and genome instability. CtBP expression induces mesenchymal and stem cell-like features, whereas CtBP depletion or caloric restriction reverses gene repression and increases DNA repair. Multiple members of the CtBP-targeted gene network are selectively downregulated in aggressive breast cancer subtypes. Differential expression of CtBP-targeted genes predicts poor clinical outcome in breast cancer patients, and elevated levels of CtBP in patient tumours predict shorter median survival. Finally, both CtBP promoter targeting and gene repression can be reversed by small molecule inhibition. These findings define broad roles for CtBP in breast cancer biology and suggest novel chromatin-based strategies for pharmacologic and metabolic intervention in cancer.