Multiplexed 3D atlas of state transitions and immune interaction in colorectal cancer.

Advanced solid cancers are complex assemblies of tumor, immune, and stromal cells characterized by high intratumoral variation. We use highly multiplexed tissue imaging, 3D reconstruction, spatial statistics, and machine learning to identify cell types and states underlying morphological features of known diagnostic and prognostic significance in colorectal cancer. Quantitation of these features in high-plex marker space reveals recurrent transitions from one tumor morphology to the next, some of which are coincident with long-range gradients in the expression of oncogenes and epigenetic regulators. At the tumor invasive margin, where tumor, normal, and immune cells compete, T cell suppression involves multiple cell types and 3D imaging shows that seemingly localized 2D features such as tertiary lymphoid structures are commonly interconnected and have graded molecular properties. Thus, while cancer genetics emphasizes the importance of discrete changes in tumor state, whole-specimen imaging reveals large-scale morphological and molecular gradients analogous to those in developing tissues.

Targeting Pin1 renders pancreatic cancer eradicable by synergizing with immunochemotherapy.

Pancreatic ductal adenocarcinoma (PDAC) is characterized by notorious resistance to current therapies attributed to inherent tumor heterogeneity and highly desmoplastic and immunosuppressive tumor microenvironment (TME). Unique proline isomerase Pin1 regulates multiple cancer pathways, but its role in the TME and cancer immunotherapy is unknown. Here, we find that Pin1 is overexpressed both in cancer cells and cancer-associated fibroblasts (CAFs) and correlates with poor survival in PDAC patients. Targeting Pin1 using clinically available drugs induces complete elimination or sustained remissions of aggressive PDAC by synergizing with anti-PD-1 and gemcitabine in diverse model systems. Mechanistically, Pin1 drives the desmoplastic and immunosuppressive TME by acting on CAFs and induces lysosomal degradation of the PD-1 ligand PD-L1 and the gemcitabine transporter ENT1 in cancer cells, besides activating multiple cancer pathways. Thus, Pin1 inhibition simultaneously blocks multiple cancer pathways, disrupts the desmoplastic and immunosuppressive TME, and upregulates PD-L1 and ENT1, rendering PDAC eradicable by immunochemotherapy.

Antigen dominance hierarchies shape TCF1+ progenitor CD8 T cell phenotypes in tumors.

CD8 T cell responses against different tumor neoantigens occur simultaneously, yet little is known about the interplay between responses and its impact on T cell function and tumor control. In mouse lung adenocarcinoma, we found that immunodominance is established in tumors, wherein CD8 T cell expansion is predominantly driven by the antigen that most stably binds MHC. T cells responding to subdominant antigens were enriched for a TCF1+ progenitor phenotype correlated with response to immune checkpoint blockade (ICB) therapy. However, the subdominant T cell response did not preferentially benefit from ICB due to a dysfunctional subset of TCF1+ cells marked by CCR6 and Tc17 differentiation. Analysis of human samples and sequencing datasets revealed that CCR6+ TCF1+ cells exist across human cancers and are not correlated with ICB response. Vaccination eliminated CCR6+ TCF1+ cells and dramatically improved the subdominant response, highlighting a strategy to optimally engage concurrent neoantigen responses against tumors.

The Human Tumor Atlas Network: Charting Tumor Transitions across Space and Time at Single-Cell Resolution.

Crucial transitions in cancer-including tumor initiation, local expansion, metastasis, and therapeutic resistance-involve complex interactions between cells within the dynamic tumor ecosystem. Transformative single-cell genomics technologies and spatial multiplex in situ methods now provide an opportunity to interrogate this complexity at unprecedented resolution. The Human Tumor Atlas Network (HTAN), part of the National Cancer Institute (NCI) Cancer Moonshot Initiative, will establish a clinical, experimental, computational, and organizational framework to generate informative and accessible three-dimensional atlases of cancer transitions for a diverse set of tumor types. This effort complements both ongoing efforts to map healthy organs and previous large-scale cancer genomics approaches focused on bulk sequencing at a single point in time. Generating single-cell, multiparametric, longitudinal atlases and integrating them with clinical outcomes should help identify novel predictive biomarkers and features as well as therapeutically relevant cell types, cell states, and cellular interactions across transitions. The resulting tumor atlases should have a profound impact on our understanding of cancer biology and have the potential to improve cancer detection, prevention, and therapeutic discovery for better precision-medicine treatments of cancer patients and those at risk for cancer.

AKT and EZH2 inhibitors kill TNBCs by hijacking mechanisms of involution.

Triple-negative breast cancer (TNBC) is the most aggressive breast cancer subtype and has the highest rate of recurrence1. The predominant standard of care for advanced TNBC is systemic chemotherapy with or without immunotherapy; however, responses are typically short lived1,2. Thus, there is an urgent need to develop more effective treatments. Components of the PI3K pathway represent plausible therapeutic targets; more than 70% of TNBCs have alterations in PIK3CA, AKT1 or PTEN3-6. However, in contrast to hormone-receptor-positive tumours, it is still unclear whether or how triple-negative disease will respond to PI3K pathway inhibitors7. Here we describe a promising AKT-inhibitor-based therapeutic combination for TNBC. Specifically, we show that AKT inhibitors synergize with agents that suppress the histone methyltransferase EZH2 and promote robust tumour regression in multiple TNBC models in vivo. AKT and EZH2 inhibitors exert these effects by first cooperatively driving basal-like TNBC cells into a more differentiated, luminal-like state, which cannot be effectively induced by either agent alone. Once TNBCs are differentiated, these agents kill them by hijacking signals that normally drive mammary gland involution. Using a machine learning approach, we developed a classifier that can be used to predict sensitivity. Together, these findings identify a promising therapeutic strategy for this highly aggressive tumour type and illustrate how deregulated epigenetic enzymes can insulate tumours from oncogenic vulnerabilities. These studies also reveal how developmental tissue-specific cell death pathways may be co-opted for therapeutic benefit.

The reprogramming of tumor stroma by HSF1 is a potent enabler of malignancy.

Stromal cells within the tumor microenvironment are essential for tumor progression and metastasis. Surprisingly little is known about the factors that drive the transcriptional reprogramming of stromal cells within tumors. We report that the transcriptional regulator heat shock factor 1 (HSF1) is frequently activated in cancer-associated fibroblasts (CAFs), where it is a potent enabler of malignancy. HSF1 drives a transcriptional program in CAFs that complements, yet is completely different from, the program it drives in adjacent cancer cells. This CAF program is uniquely structured to support malignancy in a non-cell-autonomous way. Two central stromal signaling molecules-TGF-ß and SDF1-play a critical role. In early-stage breast and lung cancer, high stromal HSF1 activation is strongly associated with poor patient outcome. Thus, tumors co-opt the ancient survival functions of HSF1 to orchestrate malignancy in both cell-autonomous and non-cell-autonomous ways, with far-reaching therapeutic implications.

Spatial genomics maps the structure, nature and evolution of cancer clones.

Genome sequencing of cancers often reveals mosaics of different subclones present in the same tumour1-3. Although these are believed to arise according to the principles of somatic evolution, the exact spatial growth patterns and underlying mechanisms remain elusive4,5. Here, to address this need, we developed a workflow that generates detailed quantitative maps of genetic subclone composition across whole-tumour sections. These provide the basis for studying clonal growth patterns, and the histological characteristics, microanatomy and microenvironmental composition of each clone. The approach rests on whole-genome sequencing, followed by highly multiplexed base-specific in situ sequencing, single-cell resolved transcriptomics and dedicated algorithms to link these layers. Applying the base-specific in situ sequencing workflow to eight tissue sections from two multifocal primary breast cancers revealed intricate subclonal growth patterns that were validated by microdissection. In a case of ductal carcinoma in situ, polyclonal neoplastic expansions occurred at the macroscopic scale but segregated within microanatomical structures. Across the stages of ductal carcinoma in situ, invasive cancer and lymph node metastasis, subclone territories are shown to exhibit distinct transcriptional and histological features and cellular microenvironments. These results provide examples of the benefits afforded by spatial genomics for deciphering the mechanisms underlying cancer evolution and microenvironmental ecology.

HSF1 drives a transcriptional program distinct from heat shock to support highly malignant human cancers.

Heat-Shock Factor 1 (HSF1), master regulator of the heat-shock response, facilitates malignant transformation, cancer cell survival, and proliferation in model systems. The common assumption is that these effects are mediated through regulation of heat-shock protein (HSP) expression. However, the transcriptional network that HSF1 coordinates directly in malignancy and its relationship to the heat-shock response have never been defined. By comparing cells with high and low malignant potential alongside their nontransformed counterparts, we identify an HSF1-regulated transcriptional program specific to highly malignant cells and distinct from heat shock. Cancer-specific genes in this program support oncogenic processes: cell-cycle regulation, signaling, metabolism, adhesion and translation. HSP genes are integral to this program, however, many are uniquely regulated in malignancy. This HSF1 cancer program is active in breast, colon and lung tumors isolated directly from human patients and is strongly associated with metastasis and death. Thus, HSF1 rewires the transcriptome in tumorigenesis, with prognostic and therapeutic implications.

BRAF-MEK Inhibition in Newly Diagnosed Papillary Craniopharyngiomas.

BACKGROUND: Craniopharyngiomas, primary brain tumors of the pituitary-hypothalamic axis, can cause clinically significant sequelae. Treatment with the use of surgery, radiation, or both is often associated with substantial morbidity related to vision loss, neuroendocrine dysfunction, and memory loss. Genotyping has shown that more than 90% of papillary craniopharyngiomas carry BRAF V600E mutations, but data are lacking with regard to the safety and efficacy of BRAF-MEK inhibition in patients with papillary craniopharyngiomas who have not undergone previous radiation therapy. METHODS: Eligible patients who had papillary craniopharyngiomas that tested positive for BRAF mutations, had not undergone radiation therapy previously, and had measurable disease received the BRAF-MEK inhibitor combination vemurafenib-cobimetinib in 28-day cycles. The primary end point of this single-group, phase 2 study was objective response at 4 months as determined with the use of centrally determined volumetric data. RESULTS: Of the 16 patients in the study, 15 (94%; 95% confidence interval [CI], 70 to 100) had a durable objective partial response or better to therapy. The median reduction in the volume of the tumor was 91% (range, 68 to 99). The median follow-up was 22 months (95% CI, 19 to 30) and the median number of treatment cycles was 8. Progression-free survival was 87% (95% CI, 57 to 98) at 12 months and 58% (95% CI, 10 to 89) at 24 months. Three patients had disease progression during follow-up after therapy had been discontinued; none have died. The sole patient who did not have a response stopped treatment after 8 days owing to toxic effects. Grade 3 adverse events that were at least possibly related to treatment occurred in 12 patients, including rash in 6 patients. In 2 patients, grade 4 adverse events (hyperglycemia in 1 patient and increased creatine kinase levels in 1 patient) were reported; 3 patients discontinued treatment owing to adverse events. CONCLUSIONS: In this small, single-group study involving patients with papillary craniopharyngiomas, 15 of 16 patients had a partial response or better to the BRAF-MEK inhibitor combination vemurafenib-cobimetinib. (Funded by the National Cancer Institute and others; ClinicalTrials.gov number, NCT03224767.).

MCMICRO: a scalable, modular image-processing pipeline for multiplexed tissue imaging.

Highly multiplexed tissue imaging makes detailed molecular analysis of single cells possible in a preserved spatial context. However, reproducible analysis of large multichannel images poses a substantial computational challenge. Here, we describe a modular and open-source computational pipeline, MCMICRO, for performing the sequential steps needed to transform whole-slide images into single-cell data. We demonstrate the use of MCMICRO on tissue and tumor images acquired using multiple imaging platforms, thereby providing a solid foundation for the continued development of tissue imaging software.

Derivation of pre-X inactivation human embryonic stem cells under physiological oxygen concentrations.

The presence of two active X chromosomes (XaXa) is a hallmark of the ground state of pluripotency specific to murine embryonic stem cells (ESCs). Human ESCs (hESCs) invariably exhibit signs of X chromosome inactivation (XCI) and are considered developmentally more advanced than their murine counterparts. We describe the establishment of XaXa hESCs derived under physiological oxygen concentrations. Using these cell lines, we demonstrate that (1) differentiation of hESCs induces random XCI in a manner similar to murine ESCs, (2) chronic exposure to atmospheric oxygen is sufficient to induce irreversible XCI with minor changes of the transcriptome, (3) the Xa exhibits heavy methylation of the XIST promoter region, and (4) XCI is associated with demethylation and transcriptional activation of XIST along with H3K27-me3 deposition across the Xi. These findings indicate that the human blastocyst contains pre-X-inactivation cells and that this state is preserved in vitro through culture under physiological oxygen.

Integrating Systemic Therapies into the Multimodality Therapy of Patients with Craniopharyngioma.

OPINION STATEMENT: The integration of targeted therapy into the multimodal management of craniopharyngiomas represents a significant advancement in the field of neuro-oncology. Historically, the management of these tumors has been challenging due to their proximity to vital brain structures, necessitating a delicate balance between tumor control and the preservation of neurological function. Traditional treatment modalities, such as surgical resection and radiation, while effective, carry their own set of risks, including potential damage to surrounding healthy tissues and the potential for long-term side effects. Recent insights into the molecular biology of craniopharyngiomas, particularly the discovery of the BRAF V600E mutation in nearly all papillary craniopharyngiomas, have paved the way for a targeted systemic treatment approach. However, advances have been limited for adamantinomatous craniopharyngiomas. The success of BRAF/MEK inhibitors in clinical trials underscores the potential of these targeted therapies not only to control tumor growth but also to reduce the need for more invasive treatments, potentially minimizing treatment-related complications. However, the introduction of these novel therapies also brings forth new challenges, such as determining the optimal timing, sequencing, and duration of targeted treatments. Furthermore, there are open questions regarding which specific BRAF/MEK inhibitors to use, the potential need for combination therapy, and the strategies for managing intolerable adverse events. Finally, ensuring equitable access to these therapies, especially in healthcare systems with limited resources, is crucial to prevent widening healthcare disparities. In conclusion, targeted therapy with BRAF/MEK inhibitors holds great promise for improving outcomes and quality of life for patients with BRAF-mutated craniopharyngiomas. However, additional research is needed to address the questions that remain about its optimal use and integration into comprehensive treatment plans.

Chemical QuantArray: A Quantitative Tool for Mass Spectrometry Imaging.

Matrix-assisted laser desorption ionization mass spectrometry imaging (MALDI-MSI) is a powerful analytical technique that provides spatially preserved detection and quantification of analytes in tissue specimens. However, clinical translation still requires improved throughput, precision, and accuracy. To accomplish this, we created "Chemical QuantArray", a gelatin tissue microarray (TMA) mold filled with serial dilutions of isotopically labeled endogenous metabolite standards. The mold is then cryo-sectioned onto a tissue homogenate to produce calibration curves. To improve precision and accuracy, we automatically remove pixels outside of each TMA well and investigated several intensity normalizations, including the utilization of a second stable isotope internal standard (IS). Chemical QuantArray enables the quantification of several endogenous metabolites over a wide dynamic range and significantly improve over current approaches. The technique reduces the space needed on the MALDI slides for calibration standards by approximately 80%. Furthermore, removal of empty pixels and normalization to an internal standard or matrix peak provided precision (<20% RSD) and accuracy (<20% DEV). Finally, we demonstrate the applicability of Chemical QuantArray by quantifying multiple purine metabolites in 14 clinical tumor specimens using a single MALDI slide. Chemical QuantArray improves the analytical characteristics and practical feasibility of MALDI-MSI metabolite quantification in clinical and translational applications.

Stitching and registering highly multiplexed whole-slide images of tissues and tumors using ASHLAR.

MOTIVATION: Stitching microscope images into a mosaic is an essential step in the analysis and visualization of large biological specimens, particularly human and animal tissues. Recent approaches to highly multiplexed imaging generate high-plex data from sequential rounds of lower-plex imaging. These multiplexed imaging methods promise to yield precise molecular single-cell data and information on cellular neighborhoods and tissue architecture. However, attaining mosaic images with single-cell accuracy requires robust image stitching and image registration capabilities that are not met by existing methods. RESULTS: We describe the development and testing of ASHLAR, a Python tool for coordinated stitching and registration of 103 or more individual multiplexed images to generate accurate whole-slide mosaics. ASHLAR reads image formats from most commercial microscopes and slide scanners, and we show that it performs better than existing open-source and commercial software. ASHLAR outputs standard OME-TIFF images that are ready for analysis by other open-source tools and recently developed image analysis pipelines. AVAILABILITY AND IMPLEMENTATION: ASHLAR is written in Python and is available under the MIT license at https://github.com/labsyspharm/ashlar. The newly published data underlying this article are available in Sage Synapse at https://dx.doi.org/10.7303/syn25826362; the availability of other previously published data re-analyzed in this article is described in Supplementary Table S4. An informational website with user guides and test data is available at https://labsyspharm.github.io/ashlar/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A comprehensive genomic study of 390 H3F3A-mutant pediatric and adult diffuse high-grade gliomas, CNS WHO grade 4.

Malignant brain tumors, known as H3K27-altered diffuse midline glioma (DMG) and H3G34-mutant diffuse hemispheric glioma (DHG), can affect individuals of all ages and are classified as CNS WHO grade 4. We comprehensively characterized 390 H3F3A-mutant diffuse gliomas (201 females, 189 males) arising in pediatric patients (under 20 years old) and adults (20 years and older) evaluated by the CGP program at Foundation Medicine between 2013 and 2020. We assessed information from pathology reports, histopathology review, and clinical data. The cohort included 304 H3K27M-mutant DMG (156 females, 148 males) and 86 H3G34-mutant DHG (45 females, 41 males). Median patient age was 20 years (1-74 years). The frequency of H3K27M-mutant DMG was similar in both pediatric and adult patients in our cohort-48.6% of the patients were over 20 years old, 31.5% over 30, and 18% over 40 at initial diagnosis. FGFR1 hotspot point mutations (N546K and K656E) were exclusively identified in H3K27M-mutant DMG tumors (64/304, 21%; p = 0.0001); these tend to occur in older patients (median age: 32.5 years) and mainly arose in the diencephalon. H3K27M-mutant DMG had higher rates of mutations in NF1 (31.0 vs 8.1%; p = 0.0001) and PIK3CA/PIK3R1 (27.9% vs 15.1%; p = 0.016) compared to H3G34-mutant DHG. However, H3G34-mutant DHG had higher rates of targetable alterations in cell-cycle pathway genes (CDK4 and CDK6 amplification; CDKN2A/B deletion) (27.0 vs 9.0%). Potentially targetable PDGFRA alterations were identified in ~ 20% of both H3G34-mutant DHG and H3K27M-mutant DMG. Overall, in the present study H3K27M-mutant DMG occurred at similar rates in both adult and patient patients. Through our analysis, we were able to identify molecular features characteristic of DMG and DHG. By identifying the recurrent co-mutations including actionable FGFR1 point mutations found in nearly one-third of H3K27M-mutant DMG in young adults, our findings can inform clinical translational studies, patient diagnosis, and clinical trial design.

PD-L1 and PD-1 expression in pediatric central nervous system germ cell tumors.

Central nervous system (CNS) germ cell tumors (GCTs) represent 2-3% of all primary CNS tumors. The majority are germinomas, which are radiosensitive and have an excellent prognosis. Contrarily, CNS non-germinomatous GCTs (NGGCTs) have less favorable prognosis and require more aggressive treatment. The expression of checkpoint/immune markers in CNS GCTs, particularly NGGCTs, is unknown. We previously reported a case of a patient whose intracranial NGGCT (predominantly choriocarcinoma) responded to immune checkpoint inhibition therapy. This case led us to evaluate our archive of intracranial GCTs for expression of PD-L1 and PD-1. With IRB approval, we searched the pathology archives at our institution for CNS GCTs. Demographic, radiologic, clinical, and histologic information was extracted from the medical records. Immunohistochemistry for lymphocytic markers (CD4, CD8, CD20), PD-1, and PD-L1 was performed. PD-L1 was considered positive if greater than 1% of tumor cells were positive and PD-1 was reported as a percentage of positive inflammatory cells. Fifty cases were identified, including 28 germinomas (mean age at diagnosis: 15.5 years; 17 males, 11 females), and 22 NGGCTs (mean age at diagnosis: 12.0 years, 21 males, 1 female). Germinomas were mostly suprasellar (17/28) and NGGCTs were predominantly pineal (17/22). Twenty-two germinomas (79%) were positive for PD-L1 expression, and 13 NGGCTs (57%) were positive for PD-L1. Cases of choriocarcinoma showed the most diffuse PD-L1 expression. PD-1 expression was seen in lymphocytes among 27/28 of the germinomas and 20/23 of the NGGCTs (ranging from 1-40% of lymphocytes). As expected, larger quantities of inflammatory cells were present in cases of germinoma. We demonstrate immune activity in CNS GCTs, and our results suggest that immune checkpoint inhibitors may be efficacious in the treatment of intracranial GCTs. Among NGGCTs, cases of choriocarcinoma showed the highest expression of PD-L1 in tumor cells, suggesting that this subtype may have the greatest benefit from checkpoint blockade.

Author Correction: Mitochondrial metabolism promotes adaptation to proteotoxic stress.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Mitochondrial metabolism promotes adaptation to proteotoxic stress.

The mechanisms by which cells adapt to proteotoxic stress are largely unknown, but are key to understanding how tumor cells, particularly in vivo, are largely resistant to proteasome inhibitors. Analysis of cancer cell lines, mouse xenografts and patient-derived tumor samples all showed an association between mitochondrial metabolism and proteasome inhibitor sensitivity. When cells were forced to use oxidative phosphorylation rather than glycolysis, they became proteasome-inhibitor resistant. This mitochondrial state, however, creates a unique vulnerability: sensitivity to the small molecule compound elesclomol. Genome-wide CRISPR-Cas9 screening showed that a single gene, encoding the mitochondrial reductase FDX1, could rescue elesclomol-induced cell death. Enzymatic function and nuclear-magnetic-resonance-based analyses further showed that FDX1 is the direct target of elesclomol, which promotes a unique form of copper-dependent cell death. These studies explain a fundamental mechanism by which cells adapt to proteotoxic stress and suggest strategies to mitigate proteasome inhibitor resistance.

Telomere length alterations and ATRX/DAXX loss in pituitary adenomas.

Telomeres are nucleoprotein complexes located at the termini of eukaryotic chromosomes that prevent exonucleolytic degradation and end-to-end chromosomal fusions. Cancers often have critically shortened, dysfunctional telomeres contributing to genomic instability. Telomere shortening has been reported in a wide range of precancerous lesions and invasive carcinomas. However, the role of telomere alterations, including the presence of alternative lengthening of telomeres (ALT), has not been studied in pituitary adenomas. Telomere length and the presence of ALT were assessed directly at the single cell level using a telomere-specific fluorescence in situ hybridization assay in tissue microarrays. Tumors were characterized as either ALT-positive or having short, normal, or long telomere lengths and then these categories were compared with clinicopathological characteristics. ATRX and DAXX expression was studied through immunohistochemistry. We characterized a discovery set of 106 pituitary adenomas including both functional and nonfunctional subsets (88 primary, 18 recurrent). Telomere lengths were estimated and we observed 64 (59.4%) cases with short, 39 (36.8%) cases with normal, and 0 (0%) cases with long telomeres. We did not observe significant differences in the clinicopathological characteristics of the group with abnormally shortened telomeres compared to the group with normal telomeres. However, three pituitary adenomas were identified as ALT-positive of which two were recurrent tumors. Two of these three ALT-positive cases had alterations in either of the chromatin remodeling proteins, ATRX and DAXX, which are routinely altered in other ALT-positive tumor subtypes. In a second cohort of 32 recurrent pituitary adenomas from 22 patients, we found that the tumors from 36% of patients (n = 8) were ALT-positive. This study demonstrates that short telomere lengths are prevalent in pituitary adenomas and that ALT-positive pituitary adenomas are enriched in recurrent disease.