Comprehensive analysis of cutaneous and uveal melanoma liver metastases.

BACKGROUND: The profound disparity in response to immune checkpoint blockade (ICB) by cutaneous melanoma (CM) and uveal melanoma (UM) patients is not well understood. Therefore, we characterized metastases of CM and UM from the same metastatic site (liver), in order to dissect the potential underlying mechanism in differential response on ICB. METHODS: Tumor liver samples from CM (n=38) and UM (n=28) patients were analyzed at the genomic (whole exome sequencing), transcriptional (RNA sequencing) and protein (immunohistochemistry and GeoMx Digital Spatial Profiling) level. RESULTS: Comparison of CM and UM metastases from the same metastatic site revealed that, although originating from the same melanocyte lineage, CM and UM differed in somatic mutation profile, copy number profile, tumor mutational burden (TMB) and consequently predicted neoantigens. A higher melanin content and higher expression of the melanoma differentiation antigen MelanA was observed in liver metastases of UM patients. No difference in B2M and human leukocyte antigen-DR (HLA-DR) expression was observed. A higher expression of programmed cell death ligand 1 (PD-L1) was found in CM compared with UM liver metastases, although the majority of CM and UM liver metastases lacked PD-L1 expression. There was no difference in the extent of immune infiltration observed between CM and UM metastases, with the exception of a higher expression of CD163 (p<0.0001) in CM liver samples. While the extent of immune infiltration was similar for CM and UM metastases, the ratio of exhausted CD8 T cells to cytotoxic T cells, to total CD8 T cells and to Th1 cells, was significantly higher in UM metastases. CONCLUSIONS: While TMB was different between CM and UM metastases, tumor immune infiltration was similar. The greater dependency on PD-L1 as an immune checkpoint in CM and the identification of higher exhaustion ratios in UM may both serve as explanations for the difference in response to ICB. Consequently, in order to improve current treatment for metastatic UM, reversal of T cell exhaustion beyond programmed cell death 1 blockade should be considered.

The Ig heavy chain protein but not its message controls early B cell development.

Development of progenitor B cells (ProB cells) into precursor B cells (PreB cells) is dictated by immunoglobulin heavy chain checkpoint (IgHCC), where the IgHC encoded by a productively rearranged Igh allele assembles into a PreB cell receptor complex (PreBCR) to generate signals to initiate this transition and suppressing antigen receptor gene recombination, ensuring that only one productive Igh allele is expressed, a phenomenon known as Igh allelic exclusion. In contrast to a productively rearranged Igh allele, the Igh messenger RNA (mRNA) (IgHR) from a nonproductively rearranged Igh allele is degraded by nonsense-mediated decay (NMD). This fact prohibited firm conclusions regarding the contribution of stable IgHR to the molecular and developmental changes associated with the IgHCC. This point was addressed by generating the IghTer5H∆TM mouse model from IghTer5H mice having a premature termination codon at position +5 in leader exon of IghTer5H allele. This prohibited NMD, and the lack of a transmembrane region (∆TM) prevented the formation of any signaling-competent PreBCR complexes that may arise as a result of read-through translation across premature Ter5 stop codon. A highly sensitive sandwich Western blot revealed read-through translation of IghTer5H message, indicating that previous conclusions regarding a role of IgHR in establishing allelic exclusion requires further exploration. As determined by RNA sequencing (RNA-Seq), this low amount of IgHC sufficed to initiate PreB cell markers normally associated with PreBCR signaling. In contrast, the IghTer5H∆TM knock-in allele, which generated stable IgHR but no detectable IgHC, failed to induce PreB development. Our data indicate that the IgHCC is controlled at the level of IgHC and not IgHR expression.

CDK4/6 inhibitors target SMARCA4-determined cyclin D1 deficiency in hypercalcemic small cell carcinoma of the ovary.

Inactivating mutations in SMARCA4 (BRG1), a key SWI/SNF chromatin remodelling gene, underlie small cell carcinoma of the ovary, hypercalcemic type (SCCOHT). To reveal its druggable vulnerabilities, we perform kinase-focused RNAi screens and uncover that SMARCA4-deficient SCCOHT cells are highly sensitive to the inhibition of cyclin-dependent kinase 4/6 (CDK4/6). SMARCA4 loss causes profound downregulation of cyclin D1, which limits CDK4/6 kinase activity in SCCOHT cells and leads to in vitro and in vivo susceptibility to CDK4/6 inhibitors. SCCOHT patient tumors are deficient in cyclin D1 yet retain the retinoblastoma-proficient/p16INK4a-deficient profile associated with positive responses to CDK4/6 inhibitors. Thus, our findings indicate that CDK4/6 inhibitors, approved for a breast cancer subtype addicted to CDK4/6 activation, could be repurposed to treat SCCOHT. Moreover, our study suggests a novel paradigm whereby critically low oncogene levels, caused by loss of a driver tumor suppressor, may also be exploited therapeutically.

Chromosomal Aberrations Associated with Sequential Steps of the Metastatic Cascade in Colorectal Cancer Patients.

BACKGROUND: Genomic information can help to identify colorectal tumors with high and low metastatic potential, thereby improving prediction of benefit of local and/or systemic treatment. Here we investigated chromosomal aberrations in relation to the different stages of the metastatic cascade: dissemination of tumor cells into the mesenteric vein, metastatic outgrowth in the liver, intravasation of the peripheral blood circulation, and development of further distant metastasis. METHODS: Peripheral and mesenteric blood from colorectal cancer patients (n = 72) were investigated for circulating tumor cells, and DNA extracted from their primary tumors was subjected to array comparative genomic hybridization profiling. The results were validated with an independent set of primary colorectal tumors (n = 53) by quantitative reverse transcription PCR. RESULTS: Mesenteric intravasation and liver metastasis were correlated with losses of chromosomes 16p (72%), 16q (27%), and 19 (54%), gain along 1q31 (45%) and 20q (60%), tumor cell infiltration into the peripheral blood circulation, and further distant metastasis with gain of chromosome 8q (59%) and 12 (47%, P < 0.01). Chromosome 12 gain was associated with poor overall survival in the initial (2.8 vs >7 years) and validation cohort (3.3 vs >6 years). The prospective study presented here is a hypothesis-generating study and confirmation with larger cohorts is required. CONCLUSIONS: This is the first study that investigated colorectal cancer in its different stages of metastasis in correlation with copy number changes of the primary tumor. This information might be helpful to identify patients with limited metastatic spread who may profit from liver metastasis resection and may lead to the discovery of new therapeutic targets.Microarray data have been deposited in NCBI's Gene Expression Omnibus and are accessible through GEO Series accession number GSE82228.

Defining chromosomal translocation risks in cancer.

Chromosomal translocations are a hallmark of cancer. Unraveling the molecular mechanism of these rare genetic events requires a clear distinction between correlative and causative risk-determinants, where technical and analytical issues can be excluded. To meet this goal, we performed in-depth analyses of publicly available genome-wide datasets. In contrast to several recent reports, we demonstrate that chromosomal translocation risk is causally unrelated to promoter stalling (Spt5), transcriptional activity, or off-targeting activity of the activation-induced cytidine deaminase. Rather, an open chromatin configuration, which is not promoter-specific, explained the elevated translocation risk of promoter regions. Furthermore, the fact that gene size directly correlates with the translocation risk in mice and human cancers further demonstrated the general irrelevance of promoter-specific activities. Interestingly, a subset of translocations observed in cancer patients likely initiates from double-strand breaks induced by an access-independent process. Together, these unexpected and novel insights are fundamental in understanding the origin of chromosome translocations and, consequently, cancer.

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.

BRCA1-like signature in triple negative breast cancer: Molecular and clinical characterization reveals subgroups with therapeutic potential.

Triple negative (TN) breast cancers make up some 15% of all breast cancers. Approximately 10-15% are mutant for the tumor suppressor, BRCA1. BRCA1 is required for homologous recombination-mediated DNA repair and deficiency results in genomic instability. BRCA1-mutated tumors have a specific pattern of genomic copy number aberrations that can be used to classify tumors as BRCA1-like or non-BRCA1-like. BRCA1 mutation, promoter methylation, BRCA1-like status and genome-wide expression data was determined for 112 TN breast cancer samples with long-term follow-up. Mutation status for 21 known DNA repair genes and PIK3CA was assessed. Gene expression and mutation frequency in BRCA1-like and non-BRCA1-like tumors were compared. Multivariate survival analysis was performed using the Cox proportional hazards model. BRCA1 germline mutation was identified in 10% of patients and 15% of tumors were BRCA1 promoter methylated. Fifty-five percent of tumors classified as BRCA1-like. The functions of genes significantly up-regulated in BRCA1-like tumors included cell cycle and DNA recombination and repair. TP53 was found to be frequently mutated in BRCA1-like (P < 0.05), while PIK3CA was frequently mutated in non-BRCA1-like tumors (P < 0.05). A significant association with worse prognosis was evident for patients with BRCA1-like tumors (adjusted HR = 3.32, 95% CI = 1.30-8.48, P = 0.01). TN tumors can be further divided into two major subgroups, BRCA1-like and non-BRCA1-like with different mutation and expression patterns and prognoses. Based on these molecular patterns, subgroups may be more sensitive to specific targeted agents such as PI3K or PARP inhibitors.

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.

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.

Differential programming of B cells in AID deficient mice.

The Aicda locus encodes the activation induced cytidine deaminase (AID) and is highly expressed in germinal center (GC) B cells to initiate somatic hypermutation (SHM) and class switch recombination (CSR) of immunoglobulin (Ig) genes. Besides these Ig specific activities in B cells, AID has been implicated in active DNA demethylation in non-B cell systems. We here determined a potential role of AID as an epigenetic eraser and transcriptional regulator in B cells. RNA-Seq on different B cell subsets revealed that Aicda(-/-) B cells are developmentally affected. However as shown by RNA-Seq, MethylCap-Seq, and SNP analysis these transcriptome alterations may not relate to AID, but alternatively to a CBA mouse strain derived region around the targeted Aicda locus. These unexpected confounding parameters provide alternative, AID-independent interpretations on genotype-phenotype correlations previously reported in numerous studies on AID using the Aicda(-/-) mouse strain.

Dissecting the gray zone between follicular lymphoma and marginal zone lymphoma using morphological and genetic features.

Nodal marginal zone lymphoma is a poorly defined entity in the World Health Organization classification, based largely on criteria of exclusion and the diagnosis often remains subjective. Follicular lymphoma lacking t(14;18) has similar characteristics which results in a major potential diagnostic overlap which this study aims to dissect. Four subgroups of lymphoma samples (n=56) were analyzed with high-resolution array comparative genome hybridization: nodal marginal zone lymphoma, t(14;18)-negative follicular lymphoma, localized t(14:18)-positive follicular lymphoma and disseminated t(14;18)-positive follicular lymphoma. Gains on chromosomes 7, 8 and 12 were observed in all subgroups. The mean number of aberrations was higher in disseminated t(14;18)-positive follicular lymphoma than in localized t(14:18)-positive follicular lymphoma (P<0.01) and the majority of alterations in localized t(14:18)-positive follicular lymphoma were also found in disseminated t(14;18)-positive follicular lymphoma. Nodal marginal zone lymphoma was marked by 3q gains with amplifications of four genes. A different overall pattern of aberrations was seen in t(14;18)-negative follicular lymphoma compared to t(14;18)-positive follicular lymphoma. t(14;18)-negative follicular lymphoma is characterized by specific (focal) gains on chromosome 3, as observed in nodal marginal zone lymphoma. Our results support the notion that localized t(14:18)-positive follicular lymphoma represents an early phase of disseminated t(14;18)-positive follicular lymphoma. t(14;18)-negative follicular lymphoma bears aberrations that are more like those in nodal marginal zone lymphoma, suggesting a relation between these groups.

Deciphering the glycosylome of dystroglycanopathies using haploid screens for lassa virus entry.

Glycosylated α-dystroglycan (α-DG) serves as cellular entry receptor for multiple pathogens, and defects in its glycosylation cause hereditary Walker-Warburg syndrome (WWS). At least eight proteins are critical to glycosylate α-DG, but many genes mutated in WWS remain unknown. To identify modifiers of α-DG, we performed a haploid screen for Lassa virus entry, a hemorrhagic fever virus causing thousands of deaths annually that hijacks glycosylated α-DG to enter cells. In complementary screens, we profiled cells for absence of α-DG carbohydrate chains or biochemically related glycans. This revealed virus host factors and a suite of glycosylation units, including all known Walker-Warburg genes and five additional factors critical for the modification of α-DG. Our findings accentuate the complexity of this posttranslational feature and point out genes defective in dystroglycanopathies.

Dosage-dependent tumor suppression by histone deacetylases 1 and 2 through regulation of c-Myc collaborating genes and p53 function.

Histone deacetylases (HDACs) are epigenetic erasers of lysine-acetyl marks. Inhibition of HDACs using small molecule inhibitors (HDACi) is a potential strategy in the treatment of various diseases and is approved for treating hematological malignancies. Harnessing the therapeutic potential of HDACi requires knowledge of HDAC-function in vivo. Here, we generated a thymocyte-specific gradient of HDAC-activity using compound conditional knockout mice for Hdac1 and Hdac2. Unexpectedly, gradual loss of HDAC-activity engendered a dosage-dependent accumulation of immature thymocytes and correlated with the incidence and latency of monoclonal lymphoblastic thymic lymphomas. Strikingly, complete ablation of Hdac1 and Hdac2 abrogated lymphomagenesis due to a block in early thymic development. Genomic, biochemical and functional analyses of pre-leukemic thymocytes and tumors revealed a critical role for Hdac1/Hdac2-governed HDAC-activity in regulating a p53-dependent barrier to constrain Myc-overexpressing thymocytes from progressing into lymphomas by regulating Myc-collaborating genes. One Myc-collaborating and p53-suppressing gene, Jdp2, was derepressed in an Hdac1/2-dependent manner and critical for the survival of Jdp2-overexpressing lymphoma cells. Although reduced HDAC-activity facilitates oncogenic transformation in normal cells, resulting tumor cells remain highly dependent on HDAC-activity, indicating that a critical level of Hdac1 and Hdac2 governed HDAC-activity is required for tumor maintenance.

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.

Search for a gene expression signature of breast cancer local recurrence in young women.

PURPOSE: A gene expression signature, predictive for local recurrence after breast-conserving treatment, has previously been identified from a series of 165 young patients with breast cancer. We evaluated this signature on both another platform and an independent series, compared its performance with other published gene-sets, and investigated the gene expression profile of a larger data set. EXPERIMENTAL DESIGN: Gene expression tumor profiles were obtained on 148 of the initial 165 Dutch patients and on an independent validation series of 195 French patients. Both unsupervised and supervised classifications were used to study the gene expression profile of the 343 breast cancers and to identify subgroups that differ for their risk of local recurrence. RESULTS: The previous local recurrence signature was validated across platforms. However, when applied to the French patients, the signature did not reproduce its reported performance and did not better classify the patients than other published gene sets. Hierarchical clustering of all 343 breast cancers did not show any grouping reflecting local recurrence status. Genes related to proliferation were found differentially expressed between patients with or without local recurrence only in triple-negative tumors. Supervised classification revealed no significant gene set predictive for local recurrence or able to outperform classification based on clinical variables. CONCLUSIONS: Although the previously identified local recurrence signature was robust on another platform, we were neither able to validate it on an independent data set, nor able to define a strong gene expression classifier for local recurrence using a larger data set. We conclude that there are no significant differences in gene expression pattern in tumors from patients with and without local recurrence after breast-conserving treatment.

A barcode screen for epigenetic regulators reveals a role for the NuB4/HAT-B histone acetyltransferase complex in histone turnover.

Dynamic modification of histone proteins plays a key role in regulating gene expression. However, histones themselves can also be dynamic, which potentially affects the stability of histone modifications. To determine the molecular mechanisms of histone turnover, we developed a parallel screening method for epigenetic regulators by analyzing chromatin states on DNA barcodes. Histone turnover was quantified by employing a genetic pulse-chase technique called RITE, which was combined with chromatin immunoprecipitation and high-throughput sequencing. In this screen, the NuB4/HAT-B complex, containing the conserved type B histone acetyltransferase Hat1, was found to promote histone turnover. Unexpectedly, the three members of this complex could be functionally separated from each other as well as from the known interacting factor and histone chaperone Asf1. Thus, systematic and direct interrogation of chromatin structure on DNA barcodes can lead to the discovery of genes and pathways involved in chromatin modification and dynamics.

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.