The evolution of metastatic upper tract urothelial carcinoma through genomic-transcriptomic and single-cell protein markers analysis.

The molecular characteristics of metastatic upper tract urothelial carcinoma (UTUC) are not well understood, and there is a lack of knowledge regarding the genomic and transcriptomic differences between primary and metastatic UTUC. To address these gaps, we integrate whole-exome sequencing, RNA sequencing, and Imaging Mass Cytometry using lanthanide metal-conjugated antibodies of 44 tumor samples from 28 patients with high-grade primary and metastatic UTUC. We perform a spatially-resolved single-cell analysis of cancer, immune, and stromal cells to understand the evolution of primary to metastatic UTUC. We discover that actionable genomic alterations are frequently discordant between primary and metastatic UTUC tumors in the same patient. In contrast, molecular subtype membership and immune depletion signature are stable across primary and matched metastatic UTUC. Molecular and immune subtypes are consistent between bulk RNA-sequencing and mass cytometry of protein markers from 340,798 single cells. Molecular subtypes at the single-cell level are highly conserved between primary and metastatic UTUC tumors within the same patient.

Whole-Genome Sequencing Analysis of Male Breast Cancer Unveils Novel Structural Events and Potential Therapeutic Targets.

The molecular characterization of male breast cancer (MaBC) has received limited attention in research, mostly because of its low incidence rate, accounting for only 0.5% to 1% of all reported cases of breast cancer each year. Managing MaBC presents significant challenges, with most treatment protocols being adapted from those developed for female breast cancer. Utilizing whole-genome sequencing (WGS) and state-of-the-art analyses, the genomic features of 10 MaBC cases (n = 10) were delineated and correlated with clinical and histopathologic characteristics. Using fluorescence in situ hybridization, an additional cohort of 18 patients was interrogated to supplement WGS findings. The genomic landscape of MaBC uncovered significant genetic alterations that could influence diagnosis and treatment. We found common somatic mutations in key driver genes, such as FAT1, GATA3, SMARCA4, and ARID2. Our study also mapped out structural variants that impact cancer-associated genes, such as ARID1A, ESR1, GATA3, NTRK1, and NF1. Using a WGS-based classifier, homologous recombination deficiency (HRD) was identified in 2 cases, both presenting with deleterious variants in BRCA2. Noteworthy was the observation of FGFR1 amplification in 21% of cases. Altogether, we identified at least 1 potential therapeutic target in 8 of the 10 cases, including high tumor mutational burden, FGFR1 amplification, and HRD. Our study is the first WGS characterization of MaBC, which uncovered potentially relevant variants, including structural events in cancer genes, HRD signatures, and germline pathogenic mutations. Our results demonstrate unique genetic markers and potential treatment targets in MaBC, thereby underlining the necessity of tailoring treatment strategies for this understudied patient population. These WGS-based findings add to the growing knowledge of MaBC genomics and highlight the need to expand research on this type of cancer.

Cooperativity of c-MYC with Krüppel-Like Factor 6 Splice Variant 1 induces phenotypic plasticity and promotes prostate cancer progression and metastasis.

Metastasis remains a major cause of morbidity and mortality in men with prostate cancer, and the functional impact of the genetic alterations, alone or in combination, driving metastatic disease remains incompletely understood. The proto-oncogene c-MYC, commonly deregulated in prostate cancer. Transgenic expression of c-MYC is sufficient to drive the progression to prostatic intraepithelial neoplasia and ultimately to moderately differentiated localized primary tumors, however, c-MYC-driven tumors are unable to progress through the metastatic cascade, suggesting that a "second-hit" is necessary in the milieu of aberrant c-MYC-driven signaling. Here, we identified cooperativity between c-MYC and KLF6-SV1, an oncogenic splice variant of the KLF6 gene. Transgenic mice that co-expressed KLF6-SV1 and c-MYC developed progressive and metastatic prostate cancer with a histological and molecular phenotype like human prostate cancer. Silencing c-MYC expression significantly reduced tumor burden in these mice supporting the necessity for c-MYC in tumor maintenance. Unbiased global proteomic analysis of tumors from these mice revealed significantly enriched vimentin, a dedifferentiation and pro-metastatic marker, induced by KLF6-SV1. c-MYC-positive tumors were also significantly enriched for KLF6-SV1 in human prostate cancer specimens. Our findings provide evidence that KLF6-SV1 is an enhancer of c-MYC-driven prostate cancer progression and metastasis, and a correlated genetic event in human prostate cancer with potential translational significance.

NIPBL::NACC1 Fusion Hepatic Carcinoma.

Several reports describing a rare primary liver tumor with histologic features reminiscent of follicular thyroid neoplasms have been published under a variety of descriptive terms including thyroid-like, solid tubulocystic, and cholangioblastic cholangiocarcinoma. Although these tumors are considered to represent histologic variants, they lack classic features of cholangiocarcinoma and have unique characteristics, namely immunoreactivity for inhibin and NIPBL::NACC1 fusions. The purpose of this study is to present clinicopathologic and molecular data for a large series of these tumors to better understand their pathogenesis. We identified 11 hepatic tumors with these features. Immunohistochemical and NACC1 and NIPBL fluorescence in situ hybridization assays were performed on all cases. Four cases had available material for whole-genome sequencing (WGS) analysis. Most patients were adult women (mean age: 42 y) who presented with abdominal pain and large hepatic masses (mean size: 14 cm). Ten patients had no known liver disease. Of the patients with follow-up information, 3/9 (33%) pursued aggressive behavior. All tumors were composed of bland cuboidal cells with follicular and solid/trabecular growth patterns in various combinations, were immunoreactive for inhibin, showed albumin mRNA by in situ hybridization, and harbored the NIPBL::NACC1 fusion by fluorescence in situ hybridization. WGS corroborated the presence of the fusion in all 4 tested cases, high tumor mutational burden in 2 cases, and over 30 structural variants per case in 3 sequenced tumors. The cases lacked mutations typical of conventional intrahepatic cholangiocarcinoma. In this report, we describe the largest series of primary inhibin-positive hepatic neoplasms harboring a NIPBL::NACC1 fusion and the first WGS analysis of these tumors. We propose to name this neoplasm NIPBL:NACC1 fusion hepatic carcinoma.

Evolution of structural rearrangements in prostate cancer intracranial metastases.

Intracranial metastases in prostate cancer are uncommon but clinically aggressive. A detailed molecular characterization of prostate cancer intracranial metastases would improve our understanding of their pathogenesis and the search for new treatment strategies. We evaluated the clinical and molecular characteristics of 36 patients with metastatic prostate cancer to either the dura or brain parenchyma. We performed whole genome sequencing (WGS) of 10 intracranial prostate cancer metastases, as well as WGS of primary prostate tumors from men who later developed metastatic disease (n = 6) and nonbrain prostate cancer metastases (n = 36). This first whole genome sequencing study of prostate intracranial metastases led to several new insights. First, there was a higher diversity of complex structural alterations in prostate cancer intracranial metastases compared to primary tumor tissues. Chromothripsis and chromoplexy events seemed to dominate, yet there were few enrichments of specific categories of structural variants compared with non-brain metastases. Second, aberrations involving the AR gene, including AR enhancer gain were observed in 7/10 (70%) of intracranial metastases, as well as recurrent loss of function aberrations involving TP53 in 8/10 (80%), RB1 in 2/10 (20%), BRCA2 in 2/10 (20%), and activation of the PI3K/AKT/PTEN pathway in 8/10 (80%). These alterations were frequently present in tumor tissues from other sites of disease obtained concurrently or sequentially from the same individuals. Third, clonality analysis points to genomic factors and evolutionary bottlenecks that contribute to metastatic spread in patients with prostate cancer. These results describe the aggressive molecular features underlying intracranial metastasis that may inform future diagnostic and treatment approaches.

Long-molecule scars of backup DNA repair in BRCA1- and BRCA2-deficient cancers.

Homologous recombination (HR) deficiency is associated with DNA rearrangements and cytogenetic aberrations1. Paradoxically, the types of DNA rearrangements that are specifically associated with HR-deficient cancers only minimally affect chromosomal structure2. Here, to address this apparent contradiction, we combined genome-graph analysis of short-read whole-genome sequencing (WGS) profiles across thousands of tumours with deep linked-read WGS of 46 BRCA1- or BRCA2-mutant breast cancers. These data revealed a distinct class of HR-deficiency-enriched rearrangements called reciprocal pairs. Linked-read WGS showed that reciprocal pairs with identical rearrangement orientations gave rise to one of two distinct chromosomal outcomes, distinguishable only with long-molecule data. Whereas one (cis) outcome corresponded to the copying and pasting of a small segment to a distant site, a second (trans) outcome was a quasi-balanced translocation or multi-megabase inversion with substantial (10 kb) duplications at each junction. We propose an HR-independent replication-restart repair mechanism to explain the full spectrum of reciprocal pair outcomes. Linked-read WGS also identified single-strand annealing as a repair pathway that is specific to BRCA2 deficiency in human cancers. Integrating these features in a classifier improved discrimination between BRCA1- and BRCA2-deficient genomes. In conclusion, our data reveal classes of rearrangements that are specific to BRCA1 or BRCA2 deficiency as a source of cytogenetic aberrations in HR-deficient cells.

The importance of escalating molecular diagnostics in patients with low-grade pediatric brain cancer.

Pilocytic astrocytomas are the most common pediatric brain tumors, typically presenting as low-grade neoplasms. We report two cases of pilocytic astrocytoma with atypical tumor progression. Case 1 involves a 12-yr-old boy with an unresectable suprasellar tumor, negative for BRAF rearrangement but harboring a BRAF p.V600E mutation. He experienced tumor size reduction and stable disease following dabrafenib treatment. Case 2 describes a 6-yr-old boy with a thalamic tumor that underwent multiple resections, with no actionable driver detected using targeted next-generation sequencing. Whole-genome and RNA-seq analysis identified an internal tandem duplication in FGFR1 and RAS pathway activation. Future management options include FGFR1 inhibitors. These cases demonstrate the importance of escalating molecular diagnostics for pediatric brain cancer, advocating for early reflexing to integrative whole-genome sequencing and transcriptomic profiling when targeted panels are uninformative. Identifying molecular drivers can significantly impact treatment decisions and improve patient outcomes.

Immunogenomic Landscape of Neuroendocrine Prostate Cancer.

PURPOSE: Patients with neuroendocrine prostate cancer (NEPC) are often managed with immunotherapy regimens extrapolated from small-cell lung cancer (SCLC). We sought to evaluate the tumor immune landscape of NEPC compared with other prostate cancer types and SCLC. EXPERIMENTAL DESIGN: In this retrospective study, a cohort of 170 patients with 230 RNA-sequencing and 104 matched whole-exome sequencing data were analyzed. Differences in immune and stromal constituents, frequency of genomic alterations, and associations with outcomes were evaluated. RESULTS: In our cohort, 36% of the prostate tumors were identified as CD8+ T-cell inflamed, whereas the remaining 64% were T-cell depleted. T-cell-inflamed tumors were enriched in anti-inflammatory M2 macrophages and exhausted T cells and associated with shorter overall survival relative to T-cell-depleted tumors (HR, 2.62; P < 0.05). Among all prostate cancer types in the cohort, NEPC was identified to be the most immune depleted, wherein only 9 out of the 36 total NEPC tumors were classified as T-cell inflamed. These inflamed NEPC cases were enriched in IFN gamma signaling and PD-1 signaling compared with other NEPC tumors. Comparison of NEPC with SCLC revealed that NEPC had poor immune content and less mutations compared with SCLC, but expression of checkpoint genes PD-L1 and CTLA-4 was comparable between NEPC and SCLC. CONCLUSIONS: NEPC is characterized by a relatively immune-depleted tumor immune microenvironment compared with other primary and metastatic prostate adenocarcinoma except in a minority of cases. These findings may inform development of immunotherapy strategies for patients with advanced prostate cancer.

GLI1 -Rearranged Enteric Tumor : Expanding the Spectrum of Gastrointestinal Neoplasms With GLI1 Gene Fusions.

GLI1 encodes a transcription factor that targets cell cycle regulators affecting stem cell proliferation. GLI1 gene fusions were initially described in pericytomas with a t[7;12] translocation and more recently in gastric plexiform fibromyxomas and gastroblastomas. This study describes the clinicopathologic, immunohistochemical, and molecular features of three intestinal-based neoplasms harboring GLI1 gene fusions. We studied three unique mesenchymal small bowel tumors. Paraffin embedded tumor tissues from these cases and 62 additional tumor samples that included a plexiform fibromyxoma were sequenced using a targeted RNAseq method to detect fusion events. The study patients included two women and one man who were 52, 80, and 22 years of age at the time of diagnosis. The tumors involved the submucosa and muscularis propria of the duodenum, jejunum, and ileum. All 3 tumors contained a proliferation of monotonous oval or spindle cells with scattered, somewhat dilated vessels. Two cases showed epithelioid structures such as glands, tubules, or nests. Immunohistochemical analysis revealed cytokeratin expression in the epithelioid components of both tumors displaying these features, and variable numbers of mesenchymal cells. Diffuse CD56 positivity was seen in the mesenchymal component of 2 tumors and desmin and smooth muscle actin staining in the other tumor. Immunostains for S-100 protein, DOG-1, and CD117 were negative in all cases. GLI1 fusions with different partner genes were detected in all tumors, and in the plexiform fibromyxoma, used as a control. Validation by fluorescence in situ hybridization was performed. None of the tumors have recurred or metastasize after surgery. We describe novel GLI1 fusions in 3 mesenchymal neoplasms of the small intestine, including 2 with biphenotypic features. Thus far, all cases have pursued indolent clinical courses. We propose the term " GLI1 -rearranged enteric tumor" to encompass this group of unique neoplasms of the small intestine that harbor GLI1 gene fusions and expand the spectrum of gastrointestinal neoplasms with these alterations.

Deciphering the origin and therapeutic targets of cancer of unknown primary: a case report that illustrates the power of integrative whole-exome and transcriptome sequencing analysis.

Cancer of unknown primary (CUP) represents a significant diagnostic and therapeutic challenge, being the third to fourth leading cause of cancer death, despite advances in diagnostic tools. This article presents a successful approach using a novel genomic analysis in the evaluation and treatment of a CUP patient, leveraging whole-exome sequencing (WES) and RNA sequencing (RNA-seq). The patient, with a history of multiple primary tumors including urothelial cancer, exhibited a history of rapid progression on empirical chemotherapy. The application of our approach identified a molecular target, characterized the tumor expression profile and the tumor microenvironment, and analyzed the origin of the tumor, leading to a tailored treatment. This resulted in a substantial radiological response across all metastatic sites and the predicted primary site of the tumor. We argue that a comprehensive genomic and molecular profiling approach, like the BostonGene© Tumor Portrait, can provide a more definitive, personalized treatment strategy, overcoming the limitations of current predictive assays. This approach offers a potential solution to an unmet clinical need for a standardized approach in identifying the tumor origin for the effective management of CUP.

Whole-genome characterization of myoepithelial carcinomas of the soft tissue.

Myoepithelial carcinomas (MECs) of soft tissue are rare and aggressive tumors affecting young adults and children, but their molecular landscape has not been comprehensively explored through genome sequencing. Here, we present the whole-exome sequencing (WES), whole-genome sequencing (WGS), and RNA sequencing findings of two MECs. Patients 1 and 2 (P1, P2), both male, were diagnosed at 27 and 37 yr of age, respectively, with shoulder (P1) and inguinal (P2) soft tissue tumors. Both patients developed metastatic disease, and P2 died of disease. P1 tumor showed a rhabdoid cytomorphology and a complete loss of INI1 (SMARCB1) expression, associated with a homozygous SMARCB1 deletion. The tumor from P2 showed a clear cell/small cell morphology, retained INI1 expression and strong S100 positivity. By WES and WGS, tumors from both patients displayed low tumor mutation burdens, and no targetable alterations in cancer genes were detected. P2's tumor harbored an EWSR1::KLF15 rearrangement, whereas the tumor from P1 showed a novel ASCC2::GGNBP2 fusion. WGS evidenced a complex genomic event involving mainly Chromosomes 17 and 22 in the tumor from P1, which was consistent with chromoplexy. These findings are consistent with previous reports of EWSR1 rearrangements (50% of cases) in MECs and provide a genetic basis for the loss of SMARCB1 protein expression observed through immunohistochemistry in 10% of 40% of MEC cases. The lack of additional driver mutations in these tumors supports the hypothesis that these alterations are the key molecular events in MEC evolution. Furthermore, the presence of complex structural variant patterns, invisible to WES, highlights the novel biological insights that can be gained through the application of WGS to rare cancers.

Unsupervised discovery of tissue architecture in multiplexed imaging.

Multiplexed imaging and spatial transcriptomics enable highly resolved spatial characterization of cellular phenotypes, but still largely depend on laborious manual annotation to understand higher-order patterns of tissue organization. As a result, higher-order patterns of tissue organization are poorly understood and not systematically connected to disease pathology or clinical outcomes. To address this gap, we developed an approach called UTAG to identify and quantify microanatomical tissue structures in multiplexed images without human intervention. Our method combines information on cellular phenotypes with the physical proximity of cells to accurately identify organ-specific microanatomical domains in healthy and diseased tissue. We apply our method to various types of images across healthy and disease states to show that it can consistently detect higher-level architectures in human tissues, quantify structural differences between healthy and diseased tissue, and reveal tissue organization patterns at the organ scale.

Case report: Metastatic urothelial cancer with an exceptional response to immunotherapy and comprehensive understanding of the tumor and the tumor microenvironment.

Although immune checkpoint inhibitors (ICIs) are increasingly used as second-line treatments for urothelial cancer (UC), only a small proportion of patients respond. Therefore, understanding the mechanisms of response to ICIs is critical to improve clinical outcomes for UC patients. The tumor microenvironment (TME) is recognized as a key player in tumor progression and the response to certain anti-cancer treatments. This study aims to investigate the mechanism of response using integrated genomic and transcriptomic profiling of a UC patient who was part of the KEYNOTE-045 trial and showed an exceptional response to pembrolizumab. Diagnosed in 2014 and receiving first-line chemotherapy without success, the patient took part in the KEYNOTE-045 trial for 2 years. She showed dramatic improvement and has now been free of disease for over 6 years. Recently described by Bagaev et al., the Molecular Functional (MF) Portrait was utilized to dissect genomic and transcriptomic features of the patient's tumor and TME. The patient's tumor was characterized as Immune Desert, which is suggestive of a non-inflamed microenvironment. Integrated whole-exome sequencing (WES) and RNA sequencing (RNA-seq) analysis identified an ATM mutation and high TMB level (33.9 mut/mb), which are both positive biomarkers for ICI response. Analysis further revealed the presence of the APOBEC complex, indicating the potential for use of APOBEC signatures as predictive biomarkers for immunotherapy response. Overall, comprehensive characterization of the patient's tumor and TME with the MF Portrait revealed important insights that could potentially be hypothesis generating to identify clinically useful biomarkers and improve treatment for UC patients.

Weakly-supervised tumor purity prediction from frozen H&E stained slides.

BACKGROUND: Estimating tumor purity is especially important in the age of precision medicine. Purity estimates have been shown to be critical for correction of tumor sequencing results, and higher purity samples allow for more accurate interpretations from next-generation sequencing results. Molecular-based purity estimates using computational approaches require sequencing of tumors, which is both time-consuming and expensive. METHODS: Here we propose an approach, weakly-supervised purity (wsPurity), which can accurately quantify tumor purity within a digitally captured hematoxylin and eosin (H&E) stained histological slide, using several types of cancer from The Cancer Genome Atlas (TCGA) as a proof-of-concept. FINDINGS: Our model predicts cancer type with high accuracy on unseen cancer slides from TCGA and shows promising generalizability to unseen data from an external cohort (F1-score of 0.83 for prostate adenocarcinoma). In addition we compare performance of our model on tumor purity prediction with a comparable fully-supervised approach on our TCGA held-out cohort and show our model has improved performance, as well as generalizability to unseen frozen slides (0.1543 MAE on an independent test cohort). In addition to tumor purity prediction, our approach identified high resolution tumor regions within a slide, and can also be used to stratify tumors into high and low tumor purity, using different cancer-dependent thresholds. INTERPRETATION: Overall, we demonstrate our deep learning model's different capabilities to analyze tumor H&E sections. We show our model is generalizable to unseen H&E stained slides from data from TCGA as well as data processed at Weill Cornell Medicine. FUNDING: Starr Cancer Consortium Grant (SCC I15-0027) to Iman Hajirasouliha.

Chromatin profiles classify castration-resistant prostate cancers suggesting therapeutic targets.

In castration-resistant prostate cancer (CRPC), the loss of androgen receptor (AR) dependence leads to clinically aggressive tumors with few therapeutic options. We used ATAC-seq (assay for transposase-accessible chromatin sequencing), RNA-seq, and DNA sequencing to investigate 22 organoids, six patient-derived xenografts, and 12 cell lines. We identified the well-characterized AR-dependent and neuroendocrine subtypes, as well as two AR-negative/low groups: a Wnt-dependent subtype, and a stem cell-like (SCL) subtype driven by activator protein-1 (AP-1) transcription factors. We used transcriptomic signatures to classify 366 patients, which showed that SCL is the second most common subtype of CRPC after AR-dependent. Our data suggest that AP-1 interacts with the YAP/TAZ and TEAD proteins to maintain subtype-specific chromatin accessibility and transcriptomic landscapes in this group. Together, this molecular classification reveals drug targets and can potentially guide therapeutic decisions.

Identifying synergistic high-order 3D chromatin conformations from genome-scale nanopore concatemer sequencing.

High-order three-dimensional (3D) interactions between more than two genomic loci are common in human chromatin, but their role in gene regulation is unclear. Previous high-order 3D chromatin assays either measure distant interactions across the genome or proximal interactions at selected targets. To address this gap, we developed Pore-C, which combines chromatin conformation capture with nanopore sequencing of concatemers to profile proximal high-order chromatin contacts at the genome scale. We also developed the statistical method Chromunity to identify sets of genomic loci with frequencies of high-order contacts significantly higher than background ('synergies'). Applying these methods to human cell lines, we found that synergies were enriched in enhancers and promoters in active chromatin and in highly transcribed and lineage-defining genes. In prostate cancer cells, these included binding sites of androgen-driven transcription factors and the promoters of androgen-regulated genes. Concatemers of high-order contacts in highly expressed genes were demethylated relative to pairwise contacts at the same loci. Synergies in breast cancer cells were associated with tyfonas, a class of complex DNA amplicons. These results rigorously link genome-wide high-order 3D interactions to lineage-defining transcriptional programs and establish Pore-C and Chromunity as scalable approaches to assess high-order genome structure.

Inhibition of FGF receptor blocks adaptive resistance to RET inhibition in CCDC6-RET-rearranged thyroid cancer.

Genetic alterations in RET lead to activation of ERK and AKT signaling and are associated with hereditary and sporadic thyroid cancer and lung cancer. Highly selective RET inhibitors have recently entered clinical use after demonstrating efficacy in treating patients with diverse tumor types harboring RET gene rearrangements or activating mutations. In order to understand resistance mechanisms arising after treatment with RET inhibitors, we performed a comprehensive molecular and genomic analysis of a patient with RET-rearranged thyroid cancer. Using a combination of drug screening and proteomic and biochemical profiling, we identified an adaptive resistance to RET inhibitors that reactivates ERK signaling within hours of drug exposure. We found that activation of FGFR signaling is a mechanism of adaptive resistance to RET inhibitors that activates ERK signaling. Combined inhibition of FGFR and RET prevented the development of adaptive resistance to RET inhibitors, reduced cell viability, and decreased tumor growth in cellular and animal models of CCDC6-RET-rearranged thyroid cancer.