Macrophage-mediated myelin recycling fuels brain cancer malignancy.

Tumors growing in metabolically challenged environments, such as glioblastoma in the brain, are particularly reliant on crosstalk with their tumor microenvironment (TME) to satisfy their high energetic needs. To study the intricacies of this metabolic interplay, we interrogated the heterogeneity of the glioblastoma TME using single-cell and multi-omics analyses and identified metabolically rewired tumor-associated macrophage (TAM) subpopulations with pro-tumorigenic properties. These TAM subsets, termed lipid-laden macrophages (LLMs) to reflect their cholesterol accumulation, are epigenetically rewired, display immunosuppressive features, and are enriched in the aggressive mesenchymal glioblastoma subtype. Engulfment of cholesterol-rich myelin debris endows subsets of TAMs to acquire an LLM phenotype. Subsequently, LLMs directly transfer myelin-derived lipids to cancer cells in an LXR/Abca1-dependent manner, thereby fueling the heightened metabolic demands of mesenchymal glioblastoma. Our work provides an in-depth understanding of the immune-metabolic interplay during glioblastoma progression, thereby laying a framework to unveil targetable metabolic vulnerabilities in glioblastoma.

Distinct spatiotemporal dynamics of CD8+ T cell-derived cytokines in the tumor microenvironment.

Cells in the tumor microenvironment (TME) influence each other through secretion and sensing of soluble mediators, such as cytokines and chemokines. While signaling of interferon γ (IFNγ) and tumor necrosis factor α (TNFα) is integral to anti-tumor immune responses, our understanding of the spatiotemporal behavior of these cytokines is limited. Here, we describe a single cell transcriptome-based approach to infer which signal(s) an individual cell has received. We demonstrate that, contrary to expectations, CD8+ T cell-derived IFNγ is the dominant modifier of the TME relative to TNFα. Furthermore, we demonstrate that cell pools that show abundant IFNγ sensing are characterized by decreased expression of transforming growth factor ß (TGFß)-induced genes, consistent with IFNγ-mediated TME remodeling. Collectively, these data provide evidence that CD8+ T cell-secreted cytokines should be categorized into local and global tissue modifiers, and describe a broadly applicable approach to dissect cytokine and chemokine modulation of the TME.

CD4+ helper T cells endow cDC1 with cancer-impeding functions in the human tumor micro-environment.

Despite their low abundance in the tumor microenvironment (TME), classical type 1 dendritic cells (cDC1) play a pivotal role in anti-cancer immunity, and their abundance positively correlates with patient survival. However, their interaction with CD4+ T-cells to potentially enable the cytotoxic T lymphocyte (CTL) response has not been elucidated. Here we show that contact with activated CD4+ T-cells enables human ex vivo cDC1, but no other DC types, to induce a CTL response to cell-associated tumor antigens. Single cell transcriptomics reveals that CD4+ T-cell help uniquely optimizes cDC1 in many functions that support antigen cross-presentation and T-cell priming, while these changes don't apply to other DC types. We robustly identify "helped" cDC1 in the TME of a multitude of human cancer types by the overlap in their transcriptomic signature with that of recently defined, tumor-infiltrating DC states that prove to be positively prognostic. As predicted from the functional effects of CD4+ T-cell help, the transcriptomic signature of "helped" cDC1 correlates with tumor infiltration by CTLs and Thelper(h)-1 cells, overall survival and response to PD-1-targeting immunotherapy. These findings reveal a critical role for CD4+ T-cell help in enabling cDC1 function in the TME and may establish the helped cDC1 transcriptomic signature as diagnostic marker in cancer.

CD26-negative and CD26-positive tissue-resident fibroblasts contribute to functionally distinct CAF subpopulations in breast cancer.

Cancer-associated fibroblasts (CAFs) are abundantly present in the microenvironment of virtually all tumors and strongly impact tumor progression. Despite increasing insight into their function and heterogeneity, little is known regarding the origin of CAFs. Understanding the origin of CAF heterogeneity is needed to develop successful CAF-based targeted therapies. Through various transplantation studies in mice, we show that CAFs in both invasive lobular breast cancer and triple-negative breast cancer originate from mammary tissue-resident normal fibroblasts (NFs). Single-cell transcriptomics, in vivo and in vitro studies reveal the transition of CD26+ and CD26- NF populations into inflammatory CAFs (iCAFs) and myofibroblastic CAFs (myCAFs), respectively. Functional co-culture experiments show that CD26+ NFs transition into pro-tumorigenic iCAFs which recruit myeloid cells in a CXCL12-dependent manner and enhance tumor cell invasion via matrix-metalloproteinase (MMP) activity. Together, our data suggest that CD26+ and CD26- NFs transform into distinct CAF subpopulations in mouse models of breast cancer.

E-Cadherin Expression Distinguishes Mouse from Human Hematopoiesis in the Basophil and Erythroid Lineages.

E-cadherin is a key regulator of epithelial cell-cell adhesion, the loss of which accelerates tumor growth and invasion. E-cadherin is also expressed in hematopoietic cells as well as epithelia. The function of hematopoietic E-cadherin is, however, mostly elusive. In this study, we explored the validity of mouse models to functionally investigate the role of hematopoietic E-cadherin in human hematopoiesis. We generated a hematopoietic-specific E-cadherin knockout mouse model. In mice, hematopoietic E-cadherin is predominantly expressed within the basophil lineage, the expression of which is dispensable for the generation of basophils. However, neither E-cadherin mRNA nor protein were detected in human basophils. In contrast, human hematopoietic E-cadherin marks the erythroid lineage. E-cadherin expression in hematopoiesis thereby revealed striking evolutionary differences between the basophil and erythroid cell lineage in humans and mice. This is remarkable as E-cadherin expression in epithelia is highly conserved among vertebrates including humans and mice. Our study therefore revealed that the mouse does not represent a suitable model to study the function of E-cadherin in human hematopoiesis and an alternative means to study the role of E-cadherin in human erythropoiesis needs to be developed.

Kinome capture sequencing of high-grade serous ovarian carcinoma reveals novel mutations in the JAK3 gene.

High-grade serous ovarian carcinoma (HGSOC) remains the deadliest form of epithelial ovarian cancer and despite major efforts little improvement in overall survival has been achieved. Identification of recurring "driver" genetic lesions has the potential to enable design of novel therapies for cancer. Here, we report on a study to find such new therapeutic targets for HGSOC using exome-capture sequencing approach targeting all kinase genes in 127 patient samples. Consistent with previous reports, the most frequently mutated gene was TP53 (97% mutation frequency) followed by BRCA1 (10% mutation frequency). The average mutation frequency of the kinase genes mutated from our panel was 1.5%. Intriguingly, after BRCA1, JAK3 was the most frequently mutated gene (4% mutation frequency). We tested the transforming properties of JAK3 mutants using the Ba/F3 cell-based in vitro functional assay and identified a novel gain-of-function mutation in the kinase domain of JAK3 (p.T1022I). Importantly, p.T1022I JAK3 mutants displayed higher sensitivity to the JAK3-selective inhibitor Tofacitinib compared to controls. For independent validation, we re-sequenced the entire JAK3 coding sequence using tagged amplicon sequencing (TAm-Seq) in 463 HGSOCs resulting in an overall somatic mutation frequency of 1%. TAm-Seq screening of CDK12 in the same population revealed a 7% mutation frequency. Our data confirms that the frequency of mutations in kinase genes in HGSOC is low and provides accurate estimates for the frequency of JAK3 and CDK12 mutations in a large well characterized cohort. Although p.T1022I JAK3 mutations are rare, our functional validation shows that if detected they should be considered as potentially actionable for therapy. The observation of CDK12 mutations in 7% of HGSOC cases provides a strong rationale for routine somatic testing, although more functional and clinical characterization is required to understand which nonsynonymous mutations alterations are associated with homologous recombination deficiency.

Endonuclease FEN1 Coregulates ERα Activity and Provides a Novel Drug Interface in Tamoxifen-Resistant Breast Cancer.

Estrogen receptor α (ERα) is a key transcriptional regulator in the majority of breast cancers. ERα-positive patients are frequently treated with tamoxifen, but resistance is common. In this study, we refined a previously identified 111-gene outcome prediction-classifier, revealing FEN1 as the strongest determining factor in ERα-positive patient prognostication. FEN1 levels were predictive of outcome in tamoxifen-treated patients, and FEN1 played a causal role in ERα-driven cell growth. FEN1 impacted the transcriptional activity of ERα by facilitating coactivator recruitment to the ERα transcriptional complex. FEN1 blockade induced proteasome-mediated degradation of activated ERα, resulting in loss of ERα-driven gene expression and eradicated tumor cell proliferation. Finally, a high-throughput 465,195 compound screen identified a novel FEN1 inhibitor, which effectively blocked ERα function and inhibited proliferation of tamoxifen-resistant cell lines as well as ex vivo-cultured ERα-positive breast tumors. Collectively, these results provide therapeutic proof of principle for FEN1 blockade in tamoxifen-resistant breast cancer. SIGNIFICANCE: These findings show that pharmacologic inhibition of FEN1, which is predictive of outcome in tamoxifen-treated patients, effectively blocks ERα function and inhibits proliferation of tamoxifen-resistant tumor cells.

Fanconi anemia and homologous recombination gene variants are associated with functional DNA repair defects in vitro and poor outcome in patients with advanced head and neck squamous cell carcinoma.

Mutations in Fanconi Anemia or Homologous Recombination (FA/HR) genes can cause DNA repair defects and could therefore impact cancer treatment response and patient outcome. Their functional impact and clinical relevance in head and neck squamous cell carcinoma (HNSCC) is unknown. We therefore questioned whether functional FA/HR defects occurred in HNSCC and whether they are associated with FA/HR variants. We assayed a panel of 29 patient-derived HNSCC cell lines and found that a considerable fraction is hypersensitive to the crosslinker Mitomycin C and PARP inhibitors, a functional measure of FA/HR defects. DNA sequencing showed that these hypersensitivities are associated with the presence of bi-allelic rare germline and somatic FA/HR gene variants. We next questioned whether such variants are associated with prognosis and treatment response in HNSCC patients. DNA sequencing of 77 advanced stage HNSCC tumors revealed a 19% incidence of such variants. Importantly, these variants were associated with a poor prognosis (p = 0.027; HR = 2.6, 1.1-6.0) but favorable response to high cumulative cisplatin dose. We show how an integrated in vitro functional repair and genomic analysis can improve the prognostic value of genetic biomarkers. We conclude that repair defects are marked and frequent in HNSCC and are associated with clinical outcome.

BRCA1-like profile predicts benefit of tandem high dose epirubicin-cyclophospamide-thiotepa in high risk breast cancer patients randomized in the WSG-AM01 trial.

BRCA1 is an important protein in the repair of DNA double strand breaks (DSBs), which are induced by alkylating chemotherapy. A BRCA1-like DNA copy number signature derived from tumors with a BRCA1 mutation is indicative for impaired BRCA1 function and associated with good outcome after high dose (HD) and tandem HD DSB inducing chemotherapy. We investigated whether BRCA1-like status was a predictive biomarker in the WSG AM 01 trial. WSG AM 01 randomized high-risk breast cancer patients to induction (2× epirubicin-cyclophosphamide) followed by tandem HD chemotherapy with epirubicin, cyclophosphamide and thiotepa versus dose dense chemotherapy (4× epirubicin-cyclophospamide followed by 3× cyclophosphamide-methotrexate-5-fluorouracil). We generated copy number profiles for 143 tumors and classified them as being BRCA1-like or non-BRCA1-like. Twenty-six out of 143 patients were BRCA1-like. BRCA1-like status was associated with high grade and triple negative tumors. With regard to event-free-survival, the primary endpoint of the trial, patients with a BRCA1-like tumor had a hazard rate of 0.2, 95% confidence interval (CI): 0.07-0.63, p = 0.006. In the interaction analysis, the combination of BRCA1-like status and HD chemotherapy had a hazard rate of 0.19, 95% CI: 0.067-0.54, p = 0.003. Similar results were observed for overall survival. These findings suggest that BRCA1-like status is a predictor for benefit of tandem HD chemotherapy with epirubicin-thiotepa-cyclophosphamide.

ERBB2 Mutations Characterize a Subgroup of Muscle-invasive Bladder Cancers with Excellent Response to Neoadjuvant Chemotherapy.

UNLABELLED: A pathologic complete response to neoadjuvant chemotherapy (NAC) containing platinum is a strong prognostic determinant for patients with muscle-invasive bladder cancer (MIBC). Despite comprehensive molecular characterization of bladder cancer, associations of molecular alterations with treatment response are still largely unknown. We selected pathologic complete responders (ypT0N0; n=38) and nonresponders (higher than ypT2; n=33) from a cohort of high-grade MIBC patients treated with NAC. DNA was isolated from prechemotherapy tumor tissue and used for next-generation sequencing of 178 cancer-associated genes (discovery cohort) or targeted sequencing (validation cohort). We found that 9 of 38 complete responders had erb-b2 receptor tyrosine kinase 2 (ERBB2) missense mutations, whereas none of 33 nonresponders had ERBB2 mutations (p=0.003). ERBB2 missense mutations in complete responders were mostly confirmed activating mutations. ERCC2 missense mutations, recently found associated with response to NAC, were more common in complete responders; however, this association did not reach statistical significance in our cohort. We conclude that ERBB2 missense mutations characterize a subgroup of MIBC patients with an excellent response to NAC. PATIENT SUMMARY: In this report we looked for genetic alterations that can predict the response to neoadjuvant chemotherapy (NAC) in bladder cancer. We found that mutations in the gene ERBB2 are exclusively present in patients responding to NAC.

Robust BRCA1-like classification of copy number profiles of samples repeated across different datasets and platforms.

Breast cancers with BRCA1 germline mutation have a characteristic DNA copy number (CN) pattern. We developed a test that assigns CN profiles to be 'BRCA1-like' or 'non-BRCA1-like', which refers to resembling a BRCA1-mutated tumor or resembling a tumor without a BRCA1 mutation, respectively. Approximately one third of the BRCA1-like breast cancers have a BRCA1 mutation, one third has hypermethylation of the BRCA1 promoter and one third has an unknown reason for being BRCA1-like. This classification is indicative of patients' response to high dose alkylating and platinum containing chemotherapy regimens, which targets the inability of BRCA1 deficient cells to repair DNA double strand breaks. We investigated whether this classification can be reliably obtained with next generation sequencing and copy number platforms other than the bacterial artificial chromosome (BAC) array Comparative Genomic Hybridization (aCGH) on which it was originally developed. We investigated samples from 230 breast cancer patients for which a CN profile had been generated on two to five platforms, comprising low coverage CN sequencing, CN extraction from targeted sequencing panels (CopywriteR), Affymetrix SNP6.0, 135K/720K oligonucleotide aCGH, Affymetrix Oncoscan FFPE (MIP) technology, 3K BAC and 32K BAC aCGH. Pairwise comparison of genomic position-mapped profiles from the original aCGH platform and other platforms revealed concordance. For most cases, biological differences between samples exceeded the differences between platforms within one sample. We observed the same classification across different platforms in over 80% of the patients and kappa values of at least 0.36. Differential classification could be attributed to CN profiles that were not strongly associated to one class. In conclusion, we have shown that the genomic regions that define our BRCA1-like classifier are robustly measured by different CN profiling technologies, providing the possibility to retro- and prospectively investigate BRCA1-like classification across a wide range of CN platforms.

REV7 counteracts DNA double-strand break resection and affects PARP inhibition.

Error-free repair of DNA double-strand breaks (DSBs) is achieved by homologous recombination (HR), and BRCA1 is an important factor for this repair pathway. In the absence of BRCA1-mediated HR, the administration of PARP inhibitors induces synthetic lethality of tumour cells of patients with breast or ovarian cancers. Despite the benefit of this tailored therapy, drug resistance can occur by HR restoration. Genetic reversion of BRCA1-inactivating mutations can be the underlying mechanism of drug resistance, but this does not explain resistance in all cases. In particular, little is known about BRCA1-independent restoration of HR. Here we show that loss of REV7 (also known as MAD2L2) in mouse and human cell lines re-establishes CTIP-dependent end resection of DSBs in BRCA1-deficient cells, leading to HR restoration and PARP inhibitor resistance, which is reversed by ATM kinase inhibition. REV7 is recruited to DSBs in a manner dependent on the H2AX-MDC1-RNF8-RNF168-53BP1 chromatin pathway, and seems to block HR and promote end joining in addition to its regulatory role in DNA damage tolerance. Finally, we establish that REV7 blocks DSB resection to promote non-homologous end-joining during immunoglobulin class switch recombination. Our results reveal an unexpected crucial function of REV7 downstream of 53BP1 in coordinating pathological DSB repair pathway choices in BRCA1-deficient cells.

High-throughput identification of antigen-specific TCRs by TCR gene capture.

The transfer of T cell receptor (TCR) genes into patient T cells is a promising approach for the treatment of both viral infections and cancer. Although efficient methods exist to identify antibodies for the treatment of these diseases, comparable strategies to identify TCRs have been lacking. We have developed a high-throughput DNA-based strategy to identify TCR sequences by the capture and sequencing of genomic DNA fragments encoding the TCR genes. We establish the value of this approach by assembling a large library of cancer germline tumor antigen-reactive TCRs. Furthermore, by exploiting the quantitative nature of TCR gene capture, we show the feasibility of identifying antigen-specific TCRs in oligoclonal T cell populations from either human material or TCR-humanized mice. Finally, we demonstrate the ability to identify tumor-reactive TCRs within intratumoral T cell subsets without knowledge of antigen specificities, which may be the first step toward the development of autologous TCR gene therapy to target patient-specific neoantigens in human cancer.

Drug-induced histone eviction from open chromatin contributes to the chemotherapeutic effects of doxorubicin.

DNA topoisomerase II inhibitors are a major class of cancer chemotherapeutics, which are thought to eliminate cancer cells by inducing DNA double-strand breaks. Here we identify a novel activity for the anthracycline class of DNA topoisomerase II inhibitors: histone eviction from open chromosomal areas. We show that anthracyclines promote histone eviction irrespective of their ability to induce DNA double-strand breaks. The histone variant H2AX, which is a key component of the DNA damage response, is also evicted by anthracyclines, and H2AX eviction is associated with attenuated DNA repair. Histone eviction deregulates the transcriptome in cancer cells and organs such as the heart, and can drive apoptosis of topoisomerase-negative acute myeloid leukaemia blasts in patients. We define a novel mechanism of action of anthracycline anticancer drugs doxorubicin and daunorubicin on chromatin biology, with important consequences for DNA damage responses, epigenetics, transcription, side effects and cancer therapy.

SERPINA6, BEX1, AGTR1, SLC26A3, and LAPTM4B are markers of resistance to neoadjuvant chemotherapy in HER2-negative breast cancer.

Response rates to chemotherapy remain highly variable in breast cancer patients. We set out to identify genes associated with chemotherapy resistance. We analyzed what is currently the largest single-institute set of gene expression profiles derived from breast cancers prior to a single neoadjuvant chemotherapy regimen (dose-dense doxorubicin and cyclophosphamide). We collected, gene expression-profiled, and analyzed 178 HER2-negative breast tumor biopsies ("NKI dataset"). We employed a recently developed approach for detecting imbalanced differential signal (DIDS) to identify markers of resistance to treatment. In contrast to traditional methods, DIDS is able to identify markers that show aberrant expression in only a small subgroup of the non-responder samples. We found a number of markers of resistance to anthracycline-based chemotherapy. We validated our findings in three external datasets, totaling 456 HER2-negative samples. Since these external sets included patients who received differing treatment regimens, the validated markers represent markers of general chemotherapy resistance. There was a highly significant overlap in the markers identified in the NKI dataset and the other three datasets. Five resistance markers, SERPINA6, BEX1, AGTR1, SLC26A3, and LAPTM4B, were identified in three of the four datasets (p value overlap < 1 × 10(-6)). These five genes identified resistant tumors that could not have been identified by merely taking ER status or proliferation into account. The identification of these genes might lead to a better understanding of the mechanisms involved in (clinically) observed chemotherapy resistance and could possibly assist in the recognition of breast cancers in which chemotherapy does not contribute to response or survival.

Glucosylated hydroxymethyluracil, DNA base J, prevents transcriptional readthrough in Leishmania.

Some Ts in nuclear DNA of trypanosomes and Leishmania are hydroxylated and glucosylated to yield base J (ß-D-glucosyl-hydroxymethyluracil). In Leishmania, about 99% of J is located in telomeric repeats. We show here that most of the remaining J is located at chromosome-internal RNA polymerase II termination sites. This internal J and telomeric J can be reduced by a knockout of J-binding protein 2 (JBP2), an enzyme involved in the first step of J biosynthesis. J levels are further reduced by growing Leishmania JBP2 knockout cells in BrdU-containing medium, resulting in cell death. The loss of internal J in JBP2 knockout cells is accompanied by massive readthrough at RNA polymerase II termination sites. The readthrough varies between transcription units but may extend over 100 kb. We conclude that J is required for proper transcription termination and infer that the absence of internal J kills Leishmania by massive readthrough of transcriptional stops.

Interactions among Polycomb domains are guided by chromosome architecture.

Polycomb group (PcG) proteins bind and regulate hundreds of genes. Previous evidence has suggested that long-range chromatin interactions may contribute to the regulation of PcG target genes. Here, we adapted the Chromosome Conformation Capture on Chip (4C) assay to systematically map chromosomal interactions in Drosophila melanogaster larval brain tissue. Our results demonstrate that PcG target genes interact extensively with each other in nuclear space. These interactions are highly specific for PcG target genes, because non-target genes with either low or high expression show distinct interactions. Notably, interactions are mostly limited to genes on the same chromosome arm, and we demonstrate that a topological rather than a sequence-based mechanism is responsible for this constraint. Our results demonstrate that many interactions among PcG target genes exist and that these interactions are guided by overall chromosome architecture.

Gene expression profiles from formalin fixed paraffin embedded breast cancer tissue are largely comparable to fresh frozen matched tissue.

BACKGROUND AND METHODS: Formalin Fixed Paraffin Embedded (FFPE) samples represent a valuable resource for cancer research. However, the discovery and development of new cancer biomarkers often requires fresh frozen (FF) samples. Recently, the Whole Genome (WG) DASL (cDNA-mediated Annealing, Selection, extension and Ligation) assay was specifically developed to profile FFPE tissue. However, a thorough comparison of data generated from FFPE RNA and Fresh Frozen (FF) RNA using this platform is lacking. To this end we profiled, in duplicate, 20 FFPE tissues and 20 matched FF tissues and evaluated the concordance of the DASL results from FFPE and matched FF material. METHODOLOGY AND PRINCIPAL FINDINGS: We show that after proper normalization, all FFPE and FF pairs exhibit a high level of similarity (Pearson correlation >0.7), significantly larger than the similarity between non-paired samples. Interestingly, the probes showing the highest correlation had a higher percentage G/C content and were enriched for cell cycle genes. Predictions of gene expression signatures developed on frozen material (Intrinsic subtype, Genomic Grade Index, 70 gene signature) showed a high level of concordance between FFPE and FF matched pairs. Interestingly, predictions based on a 60 gene DASL list (best match with the 70 gene signature) showed very high concordance with the MammaPrint® results. CONCLUSIONS AND SIGNIFICANCE: We demonstrate that data generated from FFPE material with the DASL assay, if properly processed, are comparable to data extracted from the FF counterpart. Specifically, gene expression profiles for a known set of prognostic genes for a specific disease are highly comparable between two conditions. This opens up the possibility of using both FFPE and FF material in gene expressions analyses, leading to a vast increase in the potential resources available for cancer research.

Volatile anesthetics modulate gene expression in breast and brain tumor cells.

Gene expression is increasingly used for diagnostic, prognostic, and therapeutic purposes in clinical practice. We tested the hypothesis that volatile anesthetics (VA) affect gene expression of tumor cells. Cells from the neuronal cell line SH-SY5Y and from the breast cell line MCF-7 were exposed ex vivo to enflurane, isoflurane, desflurane, halothane, sevoflurane, or nitrous oxide. Microarray gene expression profiles were studied. We observed significant differences in gene expression levels of cell cultures and response in time when exposed to different VA. Some genes used for predictive genetic fingerprints for breast cancer were affected by VA. Our findings suggest that VA modulate gene expression in breast and brain tumor cell cultures in a unique and time-dependent manner.