Transcriptional and Mutational Profiling of B-Other Acute Lymphoblastic Leukemia for Improved Diagnostics.

B-cell precursor acute lymphoblastic leukemia (BCP-ALL) is the most common cancer in children, and significant progress has been made in diagnostics and the treatment of this disease based on the subtypes of BCP-ALL. However, in a large proportion of cases (B-other), recurrent BCP-ALL-associated genomic alterations remain unidentifiable by current diagnostic procedures. In this study, we performed RNA sequencing and analyzed gene fusions, expression profiles, and mutations in diagnostic samples of 185 children with BCP-ALL. Gene expression clustering showed that a subset of B-other samples partially clusters with some of the known subgroups, particularly DUX4-positive. Mutation analysis coupled with gene expression profiling revealed the presence of distinctive BCP-ALL subgroups, characterized by the presence of mutations in known ALL driver genes, e.g., PAX5 and IKZF1. Moreover, we identified novel fusion partners of lymphoid lineage transcriptional factors ETV6, IKZF1 and PAX5. In addition, we report on low blast count detection thresholds and show that the use of EDTA tubes for sample collection does not have adverse effects on sequencing and downstream analysis. Taken together, our findings demonstrate the applicability of whole-transcriptome sequencing for personalized diagnostics in pediatric ALL, including tentative classification of the B-other cases that are difficult to diagnose using conventional methods.

Mutational mechanisms shaping the coding and noncoding genome of germinal center derived B-cell lymphomas.

B cells have the unique property to somatically alter their immunoglobulin (IG) genes by V(D)J recombination, somatic hypermutation (SHM) and class-switch recombination (CSR). Aberrant targeting of these mechanisms is implicated in lymphomagenesis, but the mutational processes are poorly understood. By performing whole genome and transcriptome sequencing of 181 germinal center derived B-cell lymphomas (gcBCL) we identified distinct mutational signatures linked to SHM and CSR. We show that not only SHM, but presumably also CSR causes off-target mutations in non-IG genes. Kataegis clusters with high mutational density mainly affected early replicating regions and were enriched for SHM- and CSR-mediated off-target mutations. Moreover, they often co-occurred in loci physically interacting in the nucleus, suggesting that mutation hotspots promote increased mutation targeting of spatially co-localized loci (termed hypermutation by proxy). Only around 1% of somatic small variants were in protein coding sequences, but in about half of the driver genes, a contribution of B-cell specific mutational processes to their mutations was found. The B-cell-specific mutational processes contribute to both lymphoma initiation and intratumoral heterogeneity. Overall, we demonstrate that mutational processes involved in the development of gcBCL are more complex than previously appreciated, and that B cell-specific mutational processes contribute via diverse mechanisms to lymphomagenesis.

Combined MEK and PI3K/p110ß Inhibition as a Novel Targeted Therapy for Malignant Mesothelioma Displaying Sarcomatoid Features.

Among malignant mesotheliomas (MM), the sarcomatoid subtype is associated with higher chemoresistance and worst survival. Due to its low incidence, there has been little progress in the knowledge of the molecular mechanisms associated with sarcomatoid MM, which might help to define novel therapeutic targets. In this work, we show that loss of PTEN expression is frequent in human sarcomatoid MM and PTEN expression levels are lower in sarcomatoid MM than in the biphasic and epithelioid subtypes. Combined Pten and Trp53 deletion in mouse mesothelium led to nonepithelioid MM development. In Pten;Trp53-null mice developing MM, the Gαi2-coupled receptor subunit activated MEK/ERK and PI3K, resulting in aggressive, immune-suppressed tumors. Combined inhibition of MEK and p110ß/PI3K reduced mouse tumor cell growth in vitro. Therapeutic inhibition of MEK and p110ß/PI3K using selumetinib (AZD6244, ARRY-142886) and AZD8186, two drugs that are currently in clinical trials, increased the survival of Pten;Trp53-null mice without major toxicity. This drug combination effectively reduced the proliferation of primary cultures of human pleural (Pl) MM, implicating nonepithelioid histology and high vimentin, AKT1/2, and Gαi2 expression levels as predictive markers of response to combined MEK and p110ß/PI3K inhibition. Our findings provide a rationale for the use of selumetinib and AZD8186 in patients with MM with sarcomatoid features. This constitutes a novel targeted therapy for a poor prognosis and frequently chemoresistant group of patients with MM, for whom therapeutic options are currently lacking. SIGNIFICANCE: Mesothelioma is highly aggressive; its sarcomatoid variants have worse prognosis. Building on a genetic mouse model, a novel combination therapy is uncovered that is relevant to human tumors.

Dynamics of Transcription Regulation in Human Bone Marrow Myeloid Differentiation to Mature Blood Neutrophils.

Neutrophils are short-lived blood cells that play a critical role in host defense against infections. To better comprehend neutrophil functions and their regulation, we provide a complete epigenetic overview, assessing important functional features of their differentiation stages from bone marrow-residing progenitors to mature circulating cells. Integration of chromatin modifications, methylation, and transcriptome dynamics reveals an enforced regulation of differentiation, for cellular functions such as release of proteases, respiratory burst, cell cycle regulation, and apoptosis. We observe an early establishment of the cytotoxic capability, while the signaling components that activate these antimicrobial mechanisms are transcribed at later stages, outside the bone marrow, thus preventing toxic effects in the bone marrow niche. Altogether, these data reveal how the developmental dynamics of the chromatin landscape orchestrate the daily production of a large number of neutrophils required for innate host defense and provide a comprehensive overview of differentiating human neutrophils.

Transcriptional regulation by NR5A2 links differentiation and inflammation in the pancreas.

Chronic inflammation increases the risk of developing one of several types of cancer. Inflammatory responses are currently thought to be controlled by mechanisms that rely on transcriptional networks that are distinct from those involved in cell differentiation. The orphan nuclear receptor NR5A2 participates in a wide variety of processes, including cholesterol and glucose metabolism in the liver, resolution of endoplasmic reticulum stress, intestinal glucocorticoid production, pancreatic development and acinar differentiation. In genome-wide association studies, single nucleotide polymorphisms in the vicinity of NR5A2 have previously been associated with the risk of pancreatic adenocarcinoma. In mice, Nr5a2 heterozygosity sensitizes the pancreas to damage, impairs regeneration and cooperates with mutant Kras in tumour progression. Here, using a global transcriptomic analysis, we describe an epithelial-cell-autonomous basal pre-inflammatory state in the pancreas of Nr5a2+/- mice that is reminiscent of the early stages of pancreatitis-induced inflammation and is conserved in histologically normal human pancreases with reduced expression of NR5A2 mRNA. In Nr5a2+/-mice, NR5A2 undergoes a marked transcriptional switch, relocating from differentiation-specific to inflammatory genes and thereby promoting gene transcription that is dependent on the AP-1 transcription factor. Pancreatic deletion of Jun rescues the pre-inflammatory phenotype, as well as binding of NR5A2 to inflammatory gene promoters and the defective regenerative response to damage. These findings support the notion that, in the pancreas, the transcriptional networks involved in differentiation-specific functions also suppress inflammatory programmes. Under conditions of genetic or environmental constraint, these networks can be subverted to foster inflammation.

c-Myc downregulation is required for preacinar to acinar maturation and pancreatic homeostasis.

BACKGROUND AND AIMS: c-Myc is highly expressed in pancreatic multipotent progenitor cells (MPC) and in pancreatic cancer. The transition from MPC to unipotent acinar progenitors is associated with c-Myc downregulation; a role for c-Myc in this process, and its possible relationship to a role in cancer, has not been established. DESIGN: Using coimmunoprecipitation assays, we demonstrate that c-Myc and Ptf1a interact. Using reverse transcriptase qPCR, western blot and immunofluorescence, we show the erosion of the acinar programme. To analyse the genomic distribution of c-Myc and Ptf1a and the global transcriptomic profile, we used ChIP-seq and RNA-seq, respectively; validation was performed with ChIP-qPCR and RT-qPCR. Lineage-tracing experiments were used to follow the effect of c-Myc overexpression in preacinar cells on acinar differentiation. RESULTS: c-Myc binds and represses the transcriptional activity of Ptf1a. c-Myc overexpression in preacinar cells leads to a massive erosion of differentiation. In adult Ela1-Myc mice: (1) c-Myc binds to Ptf1a, and Tcf3 is downregulated; (2) Ptf1a and c-Myc display partially overlapping chromatin occupancy but do not bind the same E-boxes; (3) at the proximal promoter of genes coding for digestive enzymes, we find reduced PTF1 binding and increased levels of repressive chromatin marks and PRC2 complex components. Lineage tracing of committed acinar precursors reveals that c-Myc overexpression does not restore multipotency but allows the persistence of a preacinar-like cell population. In addition, mutant KRas can lead to c-Myc overexpression and acinar dysregulation. CONCLUSIONS: c-Myc repression during development is crucial for the maturation of preacinar cells, and c-Myc overexpression can contribute to pancreatic carcinogenesis through the induction of a dedifferentiated state.

Automatic identification of informative regions with epigenomic changes associated to hematopoiesis.

Hematopoiesis is one of the best characterized biological systems but the connection between chromatin changes and lineage differentiation is not yet well understood. We have developed a bioinformatic workflow to generate a chromatin space that allows to classify 42 human healthy blood epigenomes from the BLUEPRINT, NIH ROADMAP and ENCODE consortia by their cell type. This approach let us to distinguish different cells types based on their epigenomic profiles, thus recapitulating important aspects of human hematopoiesis. The analysis of the orthogonal dimension of the chromatin space identify 32,662 chromatin determinant regions (CDRs), genomic regions with different epigenetic characteristics between the cell types. Functional analysis revealed that these regions are linked with cell identities. The inclusion of leukemia epigenomes in the healthy hematological chromatin sample space gives us insights on the healthy cell types that are more epigenetically similar to the disease samples. Further analysis of tumoral epigenetic alterations in hematopoietic CDRs points to sets of genes that are tightly regulated in leukemic transformations and commonly mutated in other tumors. Our method provides an analytical approach to study the relationship between epigenomic changes and cell lineage differentiation. Method availability: https://github.com/david-juan/ChromDet.

Wolf-Hirschhorn Syndrome Candidate 1 Is Necessary for Correct Hematopoietic and B Cell Development.

Immunodeficiency is one of the most important causes of mortality associated with Wolf-Hirschhorn syndrome (WHS), a severe rare disease originated by a deletion in chromosome 4p. The WHS candidate 1 (WHSC1) gene has been proposed as one of the main genes responsible for many of the alterations in WHS, but its mechanism of action is still unknown. Here, we present in vivo genetic evidence showing that Whsc1 plays an important role at several points of hematopoietic development. Particularly, our results demonstrate that both differentiation and function of Whsc1-deficient B cells are impaired at several key developmental stages due to profound molecular defects affecting B cell lineage specification, commitment, fitness, and proliferation, demonstrating a causal role for WHSC1 in the immunodeficiency of WHS patients.

GATA6 regulates EMT and tumour dissemination, and is a marker of response to adjuvant chemotherapy in pancreatic cancer.

BACKGROUND AND AIMS: The role of GATA factors in cancer has gained increasing attention recently, but the function of GATA6 in pancreatic ductal adenocarcinoma (PDAC) is controversial. GATA6 is amplified in a subset of tumours and was proposed to be oncogenic, but high GATA6 levels are found in well-differentiated tumours and are associated with better patient outcome. By contrast, a tumour-suppressive function of GATA6 was demonstrated using genetic mouse models. We aimed at clarifying GATA6 function in PDAC. DESIGN: We combined GATA6 silencing and overexpression in PDAC cell lines with GATA6 ChIP-Seq and RNA-Seq data, in order to understand the mechanism of GATA6 functions. We then confirmed some of our observations in primary patient samples, some of which were included in the ESPAC-3 randomised clinical trial for adjuvant therapy. RESULTS: GATA6 inhibits the epithelial-mesenchymal transition (EMT) in vitro and cell dissemination in vivo. GATA6 has a unique proepithelial and antimesenchymal function, and its transcriptional regulation is direct and implies, indirectly, the regulation of other transcription factors involved in EMT. GATA6 is lost in tumours, in association with altered differentiation and the acquisition of a basal-like molecular phenotype, consistent with an epithelial-to-epithelial (ET2) transition. Patients with basal-like GATA6low tumours have a shorter survival and have a distinctly poor response to adjuvant 5-fluorouracil (5-FU)/leucovorin. However, modulation of GATA6 expression in cultured cells does not directly regulate response to 5-FU. CONCLUSIONS: We provide mechanistic insight into GATA6 tumour-suppressive function, its role as a regulator of canonical epithelial differentiation, and propose that loss of GATA6 expression is both prognostic and predictive of response to adjuvant therapy.

CD8+ T Cells from Human Neonates Are Biased toward an Innate Immune Response.

To better understand why human neonates show a poor response to intracellular pathogens, we compared gene expression and histone modification profiles of neonatal naive CD8+ T cells with that of their adult counterparts. We found that neonatal lymphocytes have a distinct epigenomic landscape associated with a lower expression of genes involved in T cell receptor (TCR) signaling and cytotoxicity and a higher expression of genes involved in the cell cycle and innate immunity. Functional studies corroborated that neonatal CD8+ T cells are less cytotoxic, transcribe antimicrobial peptides, and produce reactive oxygen species. Altogether, our results show that neonatal CD8+ T cells have a specific genetic program biased toward the innate immune response. These findings will contribute to better diagnosis and management of the neonatal immune response.

The BLUEPRINT Data Analysis Portal.

The impact of large and complex epigenomic datasets on biological insights or clinical applications is limited by the lack of accessibility by easy, intuitive, and fast tools. Here, we describe an epigenomics comparative cyber-infrastructure (EPICO), an open-access reference set of libraries to develop comparative epigenomic data portals. Using EPICO, large epigenome projects can make available their rich datasets to the community without requiring specific technical skills. As a first instance of EPICO, we implemented the BLUEPRINT Data Analysis Portal (BDAP). BDAP provides a desktop for the comparative analysis of epigenomes of hematopoietic cell types based on results, such as the position of epigenetic features, from basic analysis pipelines. The BDAP interface facilitates interactive exploration of genomic regions, genes, and pathways in the context of differentiation of hematopoietic lineages. This work represents initial steps toward broadly accessible integrative analysis of epigenomic data across international consortia. EPICO can be accessed at https://github.com/inab, and BDAP can be accessed at http://blueprint-data.bsc.es.

Recruitment of the Mammalian Histone-modifying EMSY Complex to Target Genes Is Regulated by ZNF131.

Recent work from others and us revealed interactions between the Sin3/HDAC complex, the H3K4me3 demethylase KDM5A, GATAD1, and EMSY. Here, we characterize the EMSY/KDM5A/SIN3B complex in detail by quantitative interaction proteomics and ChIP-sequencing. We identify a novel substoichiometric interactor of the complex, transcription factor ZNF131, which recruits EMSY to a large number of active, H3K4me3 marked promoters. Interestingly, using an EMSY knock-out line and subsequent rescue experiments, we show that EMSY is in most cases positively correlated with transcriptional activity of its target genes and stimulates cell proliferation. Finally, by immunohistochemical staining of primary breast tissue microarrays we find that EMSY/KDM5A/SIN3B complex subunits are frequently overexpressed in primary breast cancer cases in a correlative manner. Taken together, these data open venues for exploring the possibility that sporadic breast cancer patients with EMSY amplification might benefit from epigenetic combination therapy targeting both the KDM5A demethylase and histone deacetylases.

BPTF is required for c-MYC transcriptional activity and in vivo tumorigenesis.

c-MYC oncogene is deregulated in most human tumours. Histone marks associated with transcriptionally active genes define high-affinity c-MYC targets. The mechanisms involved in their recognition by c-MYC are unknown. Here we report that c-MYC interacts with BPTF, a core subunit of the NURF chromatin-remodelling complex. BPTF is required for the activation of the full c-MYC transcriptional programme in fibroblasts. BPTF knockdown leads to decreased c-MYC recruitment to DNA and changes in chromatin accessibility. In Bptf-null MEFs, BPTF is necessary for c-MYC-driven proliferation, G1-S progression and replication stress, but not for c-MYC-driven apoptosis. Bioinformatics analyses unveil that BPTF levels correlate positively with c-MYC-driven transcriptional signatures. In vivo, Bptf inactivation in pre-neoplastic pancreatic acinar cells significantly delays tumour development and extends survival. Our findings uncover BPTF as a crucial c-MYC co-factor required for its biological activity and suggest that the BPTF-c-MYC axis is a potential therapeutic target in cancer.

Choline Kinase Alpha (CHKα) as a Therapeutic Target in Pancreatic Ductal Adenocarcinoma: Expression, Predictive Value, and Sensitivity to Inhibitors.

Choline kinase α (CHKα) plays a crucial role in the regulation of membrane phospholipid synthesis and has oncogenic properties in vitro. We have analyzed the expression of CHKα in cell lines derived from pancreatic ductal adenocarcinoma (PDAC) and have found increased CHKα expression, associated with differentiation. CHKα protein expression was directly correlated with sensitivity to MN58b, a CHKα inhibitor that reduced cell growth through the induction of apoptosis. Accordingly, CHKα knockdown led to reduced drug sensitivity. In addition, we found that gemcitabine-resistant PDAC cells displayed enhanced sensitivity to CHKα inhibition and, in vitro, MN58b had additive or synergistic effects with gemcitabine, 5-fluorouracil, and oxaliplatin, three active drugs in the treatment of PDAC. Using tissue microarrays, CHKα was found to be overexpressed in 90% of pancreatic tumors. While cytoplasmic CHKα did not relate to survival, nuclear CHKα distribution was observed in 43% of samples and was associated with longer survival, especially among patients with well/moderately differentiated tumors. To identify the mechanisms involved in resistance to CHKα inhibitors, we cultured IMIM-PC-2 cells with increasingly higher concentrations of MN58b and isolated a subline with a 30-fold higher IC50. RNA-Seq analysis identified upregulation of ABCB1 and ABCB4 multidrug resistance transporters, and functional studies confirmed that their upregulation is the main mechanism involved in resistance. Overall, our findings support the notion that CHKα inhibition merits further attention as a therapeutic option in patients with PDAC and that expression levels may predict response.

DNA methylome analysis in Burkitt and follicular lymphomas identifies differentially methylated regions linked to somatic mutation and transcriptional control.

Although Burkitt lymphomas and follicular lymphomas both have features of germinal center B cells, they are biologically and clinically quite distinct. Here we performed whole-genome bisulfite, genome and transcriptome sequencing in 13 IG-MYC translocation-positive Burkitt lymphoma, nine BCL2 translocation-positive follicular lymphoma and four normal germinal center B cell samples. Comparison of Burkitt and follicular lymphoma samples showed differential methylation of intragenic regions that strongly correlated with expression of associated genes, for example, genes active in germinal center dark-zone and light-zone B cells. Integrative pathway analyses of regions differentially methylated in Burkitt and follicular lymphomas implicated DNA methylation as cooperating with somatic mutation of sphingosine phosphate signaling, as well as the TCF3-ID3 and SWI/SNF complexes, in a large fraction of Burkitt lymphomas. Taken together, our results demonstrate a tight connection between somatic mutation, DNA methylation and transcriptional control in key B cell pathways deregulated differentially in Burkitt lymphoma and other germinal center B cell lymphomas.

The UBC-40 Urothelial Bladder Cancer cell line index: a genomic resource for functional studies.

BACKGROUND: Urothelial bladder cancer is a highly heterogeneous disease. Cancer cell lines are useful tools for its study. This is a comprehensive genomic characterization of 40 urothelial bladder carcinoma (UBC) cell lines including information on origin, mutation status of genes implicated in bladder cancer (FGFR3, PIK3CA, TP53, and RAS), copy number alterations assessed using high density SNP arrays, uniparental disomy (UPD) events, and gene expression. RESULTS: Based on gene mutation patterns and genomic changes we identify lines representative of the FGFR3-driven tumor pathway and of the TP53/RB tumor suppressor-driven pathway. High-density array copy number analysis identified significant focal gains (1q32, 5p13.1-12, 7q11, and 7q33) and losses (i.e. 6p22.1) in regions altered in tumors but not previously described as affected in bladder cell lines. We also identify new evidence for frequent regions of UPD, often coinciding with regions reported to be lost in tumors. Previously undescribed chromosome X losses found in UBC lines also point to potential tumor suppressor genes. Cell lines representative of the FGFR3-driven pathway showed a lower number of UPD events. CONCLUSIONS: Overall, there is a predominance of more aggressive tumor subtypes among the cell lines. We provide a cell line classification that establishes their relatedness to the major molecularly-defined bladder tumor subtypes. The compiled information should serve as a useful reference to the bladder cancer research community and should help to select cell lines appropriate for the functional analysis of bladder cancer genes, for example those being identified through massive parallel sequencing.