Germline-mediated immunoediting sculpts breast cancer subtypes and metastatic proclivity.

Tumors with the same diagnosis can have different molecular profiles and response to treatment. It remains unclear when and why these differences arise. Somatic genomic aberrations occur within the context of a highly variable germline genome. Interrogating 5870 breast cancer lesions, we demonstrated that germline-derived epitopes in recurrently amplified genes influence somatic evolution by mediating immunoediting. Individuals with a high germline-epitope burden in human epidermal growth factor receptor 2 (HER2/ERBB2) are less likely to develop HER2-positive breast cancer compared with other subtypes. The same holds true for recurrent amplicons defining three aggressive estrogen receptor (ER)-positive subgroups. Tumors that overcome such immune-mediated negative selection are more aggressive and demonstrate an "immune cold" phenotype. These data show that the germline genome plays a role in dictating somatic evolution.

Spatially Segregated Macrophage Populations Predict Distinct Outcomes In Colon Cancer.

Tumor-associated macrophages are transcriptionally heterogeneous, but the spatial distribution and cell interactions that shape macrophage tissue roles remain poorly characterized. Here, we spatially resolve five distinct human macrophage populations in normal and malignant human breast and colon tissue and reveal their cellular associations. This spatial map reveals that distinct macrophage populations reside in spatially segregated micro-environmental niches with conserved cellular compositions that are repeated across healthy and diseased tissue. We show that IL4I1+ macrophages phagocytose dying cells in areas with high cell turnover and predict good outcome in colon cancer. In contrast, SPP1+ macrophages are enriched in hypoxic and necrotic tumor regions and portend worse outcome in colon cancer. A subset of FOLR2+ macrophages is embedded in plasma cell niches. NLRP3+ macrophages co-localize with neutrophils and activate an inflammasome in tumors. Our findings indicate that a limited number of unique human macrophage niches function as fundamental building blocks in tissue.

Spatial transcriptomics identifies candidate stromal drivers of benign prostatic hyperplasia.

Benign prostatic hyperplasia (BPH) is the nodular proliferation of the prostate transition zone in older men, leading to urinary storage and voiding problems that can be recalcitrant to therapy. Decades ago, John McNeal proposed that BPH originates with the "reawakening" of embryonic inductive activity by adult prostate stroma, which spurs new ductal proliferation and branching morphogenesis. Here, by laser microdissection and transcriptional profiling of the BPH stroma adjacent to hyperplastic branching ducts, we identified secreted factors likely mediating stromal induction of prostate glandular epithelium and coinciding processes. The top stromal factors were insulin-like growth factor 1 (IGF1) and CXC chemokine ligand 13 (CXCL13), which we verified by RNA in situ hybridization to be coexpressed in BPH fibroblasts, along with their cognate receptors (IGF1R and CXCR5) on adjacent epithelium. In contrast, IGF1 but not CXCL13 was expressed in human embryonic prostate stroma. Finally, we demonstrated that IGF1 is necessary for the generation of BPH-1 cell spheroids and patient-derived BPH cell organoids in 3D culture. Our findings partially support historic speculations on the etiology of BPH and provide what we believe to be new molecular targets for rational therapies directed against the underlying processes driving BPH.

Spatial Biology of Breast Cancer.

Spatial findings have shaped on our understanding of breast cancer. In this review, we discuss how spatial methods, including spatial transcriptomics and proteomics and the resultant understanding of spatial relationships, have contributed to concepts regarding cancer progression and treatment. In addition to discussing traditional approaches, we examine how emerging multiplex imaging technologies have contributed to the field and how they might influence future research.

Cell volume expansion and local contractility drive collective invasion of the basement membrane in breast cancer.

Breast cancer becomes invasive when carcinoma cells invade through the basement membrane (BM)-a nanoporous layer of matrix that physically separates the primary tumour from the stroma. Single cells can invade through nanoporous three-dimensional matrices due to protease-mediated degradation or force-mediated widening of pores via invadopodial protrusions. However, how multiple cells collectively invade through the physiological BM, as they do during breast cancer progression, remains unclear. Here we developed a three-dimensional in vitro model of collective invasion of the BM during breast cancer. We show that cells utilize both proteases and forces-but not invadopodia-to breach the BM. Forces are generated from a combination of global cell volume expansion, which stretches the BM, and local contractile forces that act in the plane of the BM to breach it, allowing invasion. These results uncover a mechanism by which cells collectively interact to overcome a critical barrier to metastasis.

Cost-effective DNA methylation profiling by FML-seq.

Current methods for profiling DNA methylation require costly reagents, sequencing, and labor time. We introduce fragmentation at methylated loci and sequencing (FML-seq), a sequencing library protocol that greatly reduces all these costs. Relative to other techniques tested on the same human cell lines, FML-seq produces similar measurements of absolute and differential cytosine methylation at a fraction of the price. FML-seq enables inexpensive, high-throughput experimental designs for large-scale epigenetics research projects.

Synthetic whole-slide image tile generation with gene expression profile-infused deep generative models.

In this work, we propose an approach to generate whole-slide image (WSI) tiles by using deep generative models infused with matched gene expression profiles. First, we train a variational autoencoder (VAE) that learns a latent, lower-dimensional representation of multi-tissue gene expression profiles. Then, we use this representation to infuse generative adversarial networks (GANs) that generate lung and brain cortex tissue tiles, resulting in a new model that we call RNA-GAN. Tiles generated by RNA-GAN were preferred by expert pathologists compared with tiles generated using traditional GANs, and in addition, RNA-GAN needs fewer training epochs to generate high-quality tiles. Finally, RNA-GAN was able to generalize to gene expression profiles outside of the training set, showing imputation capabilities. A web-based quiz is available for users to play a game distinguishing real and synthetic tiles: https://rna-gan.stanford.edu/, and the code for RNA-GAN is available here: https://github.com/gevaertlab/RNA-GAN.

Germline-mediated immunoediting sculpts breast cancer subtypes and metastatic proclivity.

Cancer represents a broad spectrum of molecularly and morphologically diverse diseases. Individuals with the same clinical diagnosis can have tumors with drastically different molecular profiles and clinical response to treatment. It remains unclear when these differences arise during disease course and why some tumors are addicted to one oncogenic pathway over another. Somatic genomic aberrations occur within the context of an individual's germline genome, which can vary across millions of polymorphic sites. An open question is whether germline differences influence somatic tumor evolution. Interrogating 3,855 breast cancer lesions, spanning pre-invasive to metastatic disease, we demonstrate that germline variants in highly expressed and amplified genes influence somatic evolution by modulating immunoediting at early stages of tumor development. Specifically, we show that the burden of germline-derived epitopes in recurrently amplified genes selects against somatic gene amplification in breast cancer. For example, individuals with a high burden of germline-derived epitopes in ERBB2, encoding human epidermal growth factor receptor 2 (HER2), are significantly less likely to develop HER2-positive breast cancer compared to other subtypes. The same holds true for recurrent amplicons that define four subgroups of ER-positive breast cancers at high risk of distant relapse. High epitope burden in these recurrently amplified regions is associated with decreased likelihood of developing high risk ER-positive cancer. Tumors that overcome such immune-mediated negative selection are more aggressive and demonstrate an "immune cold" phenotype. These data show the germline genome plays a previously unappreciated role in dictating somatic evolution. Exploiting germline-mediated immunoediting may inform the development of biomarkers that refine risk stratification within breast cancer subtypes.

Cost-effective DNA methylation profiling by FML-seq.

Current methods for profiling DNA methylation require costly reagents, sequencing, or labor time. We introduce FML-seq, a sequencing library protocol that greatly reduces all these costs. Relative to other techniques tested on the same human cell lines, FML-seq produces similar measurements of absolute and differential cytosine methylation at a fraction of the price. FML-seq enables inexpensive, high-throughput experimental designs for large-scale epigenetics research projects.

A spatial map of human macrophage niches reveals context-dependent macrophage functions in colon and breast cancer.

Tumor-associated macrophages (TAMs) display heterogeneous phenotypes. Yet the exact tissue cues that shape macrophage functional diversity are incompletely understood. Here we discriminate, spatially resolve and reveal the function of five distinct macrophage niches within malignant and benign breast and colon tissue. We found that SPP1 TAMs reside in hypoxic and necrotic tumor regions, and a novel subset of FOLR2 tissue resident macrophages (TRMs) supports the plasma cell tissue niche. We discover that IL4I1 macrophages populate niches with high cell turnover where they phagocytose dying cells. Significantly, IL4I1 TAMs abundance correlates with anti-PD1 treatment response in breast cancer. Furthermore, NLRP3 inflammasome activation in NLRP3 TAMs correlates with neutrophil infiltration in the tumors and is associated with poor outcome in breast cancer patients. This suggests the NLRP3 inflammasome as a novel cancer immunetherapy target. Our work uncovers context-dependent roles of macrophage subsets, and suggests novel predictive markers and macrophage subset-specific therapy targets.

Interactions in CSF1-Driven Tenosynovial Giant Cell Tumors.

PURPOSE: A major component of cells in tenosynovial giant cell tumor (TGCT) consists of bystander macrophages responding to CSF1 that is overproduced by a small number of neoplastic cells with a chromosomal translocation involving the CSF1 gene. An autocrine loop was postulated where the neoplastic cells would be stimulated through CSF1R expressed on their surface. Here, we use single-cell RNA sequencing (scRNA-seq) to investigate cellular interactions in TGCT. EXPERIMENTAL DESIGN: A total of 18,788 single cells from three TGCT and two giant cell tumor of bone (GCTB) samples underwent scRNA-seq. The three TGCTs were additionally analyzed using long-read RNA sequencing. Immunofluorescence and IHC for a range of markers were used to validate and extend the scRNA-seq findings. RESULTS: Two recurrent neoplastic cell populations were identified in TGCT that are highly similar to nonneoplastic synoviocytes. We identified GFPT2 as a marker that highlights the neoplastic cells in TCGT. We show that the neoplastic cells themselves do not express CSF1R. We identified overlapping MAB features between the giant cells in TGCT and GCTB. CONCLUSIONS: The neoplastic cells in TGCT are highly similar to nonneoplastic synoviocytes. The lack of CSF1R on the neoplastic cells indicates they may be unaffected by current therapies. High expression of GFPT2 in the neoplastic cells is associated with activation of the YAP1/TAZ pathway. In addition, we identified expression of the platelet-derived growth factor receptor in the neoplastic cells. These findings suggest two additional pathways to target in this tumor.

Mesenchymal tumor cells drive adaptive resistance of Trp53-/- breast tumor cells to inactivated mutant Kras.

As precision medicine increases the response rate of treatment, tumors frequently bypass inhibition, and reoccur. In order for treatment to be effective long term, the mechanisms enabling treatment adaptation need to be understood. Here, we report a mouse model that, in the absence of p53 and the presence of oncogenic KrasG12D , develops breast tumors. Upon inactivation of KrasG12D , tumors initially regress and enter remission. Subsequently, the majority of tumors adapt to the withdrawal of KrasG12D expression and return. KrasG12D -independent tumor cells show a strong mesenchymal profile with active RAS-RAF-MEK-ERK (MAPK/ERK) signaling. Both KrasG12D -dependent and KrasG12D -independent tumors display a high level of genomic instability, and KrasG12D -independent tumors harbor numerous amplified genes that can activate the MAPK/ERK signaling pathway. Our study identifies both epithelial-mesenchymal transition (EMT) and active MAPK/ERK signaling in tumors that adapt to oncogenic KrasG12D withdrawal in a novel Trp53-/- breast cancer mouse model. To achieve long-lasting responses in the clinic to RAS-fueled cancer, treatment will need to focus in parallel on obstructing tumors from adapting to oncogene inhibition.

The microdissected gene expression landscape of nasopharyngeal cancer reveals vulnerabilities in FGF and noncanonical NF-κB signaling.

Nasopharyngeal cancer (NPC) is an Epstein-Barr virus (EBV)-positive epithelial malignancy with an extensive inflammatory infiltrate. Traditional RNA-sequencing techniques uncovered only microenvironment signatures, while the gene expression of the tumor epithelial compartment has remained a mystery. Here, we use Smart-3SEQ to prepare transcriptome-wide gene expression profiles from microdissected NPC tumors, dysplasia, and normal controls. We describe changes in biological pathways across the normal to tumor spectrum and show that fibroblast growth factor (FGF) ligands are overexpressed in NPC tumors, while negative regulators of FGF signaling, including SPRY1, SPRY2, and LGALS3, are down-regulated early in carcinogenesis. Within the NF-κB signaling pathway, the critical noncanonical transcription factors, RELB and NFKB2, are enriched in the majority of NPC tumors. We confirm the responsiveness of EBV-positive NPC cell lines to targeted inhibition of these pathways, reflecting the heterogeneity in NPC patient tumors. Our data comprehensively describe the gene expression landscape of NPC and unravel the mysteries of receptor tyrosine kinase and NF-κB pathways in NPC.

Transition to invasive breast cancer is associated with progressive changes in the structure and composition of tumor stroma.

Ductal carcinoma in situ (DCIS) is a pre-invasive lesion that is thought to be a precursor to invasive breast cancer (IBC). To understand the changes in the tumor microenvironment (TME) accompanying transition to IBC, we used multiplexed ion beam imaging by time of flight (MIBI-TOF) and a 37-plex antibody staining panel to interrogate 79 clinically annotated surgical resections using machine learning tools for cell segmentation, pixel-based clustering, and object morphometrics. Comparison of normal breast with patient-matched DCIS and IBC revealed coordinated transitions between four TME states that were delineated based on the location and function of myoepithelium, fibroblasts, and immune cells. Surprisingly, myoepithelial disruption was more advanced in DCIS patients that did not develop IBC, suggesting this process could be protective against recurrence. Taken together, this HTAN Breast PreCancer Atlas study offers insight into drivers of IBC relapse and emphasizes the importance of the TME in regulating these processes.

Immune cell topography predicts response to PD-1 blockade in cutaneous T cell lymphoma.

Cutaneous T cell lymphomas (CTCL) are rare but aggressive cancers without effective treatments. While a subset of patients derive benefit from PD-1 blockade, there is a critically unmet need for predictive biomarkers of response. Herein, we perform CODEX multiplexed tissue imaging and RNA sequencing on 70 tumor regions from 14 advanced CTCL patients enrolled in a pembrolizumab clinical trial (NCT02243579). We find no differences in the frequencies of immune or tumor cells between responders and non-responders. Instead, we identify topographical differences between effector PD-1+ CD4+ T cells, tumor cells, and immunosuppressive Tregs, from which we derive a spatial biomarker, termed the SpatialScore, that correlates strongly with pembrolizumab response in CTCL. The SpatialScore coincides with differences in the functional immune state of the tumor microenvironment, T cell function, and tumor cell-specific chemokine recruitment and is validated using a simplified, clinically accessible tissue imaging platform. Collectively, these results provide a paradigm for investigating the spatial balance of effector and suppressive T cell activity and broadly leveraging this biomarker approach to inform the clinical use of immunotherapies.

Acinar cell clonal expansion in pancreas homeostasis and carcinogenesis.

Pancreatic ductal adenocarcinoma (PDAC) is one of the leading causes of cancer deaths worldwide1. Studies in human tissues and in mouse models have suggested that for many cancers, stem cells sustain early mutations driving tumour development2,3. For the pancreas, however, mechanisms underlying cellular renewal and initiation of PDAC remain unresolved. Here, using lineage tracing from the endogenous telomerase reverse transcriptase (Tert) locus, we identify a rare TERT-positive subpopulation of pancreatic acinar cells dispersed throughout the exocrine compartment. During homeostasis, these TERThigh acinar cells renew the pancreas by forming expanding clones of acinar cells, whereas randomly marked acinar cells do not form these clones. Specific expression of mutant Kras in TERThigh acinar cells accelerates acinar clone formation and causes transdifferentiation to ductal pre-invasive pancreatic intraepithelial neoplasms by upregulating Ras-MAPK signalling and activating the downstream kinase ERK (phospho-ERK). In resected human pancreatic neoplasms, we find that foci of phospho-ERK-positive acinar cells are common and frequently contain activating KRAS mutations, suggesting that these acinar regions represent an early cancer precursor lesion. These data support a model in which rare TERThigh acinar cells may sustain KRAS mutations, driving acinar cell expansion and creating a field of aberrant cells initiating pancreatic tumorigenesis.

Unmasking the immune microecology of ductal carcinoma in situ with deep learning.

Despite increasing evidence supporting the clinical relevance of tumour infiltrating lymphocytes (TILs) in invasive breast cancer, TIL spatial variability within ductal carcinoma in situ (DCIS) samples and its association with progression are not well understood. To characterise tissue spatial architecture and the microenvironment of DCIS, we designed and validated a new deep learning pipeline, UNMaSk. Following automated detection of individual DCIS ducts using a new method IM-Net, we applied spatial tessellation to create virtual boundaries for each duct. To study local TIL infiltration for each duct, DRDIN was developed for mapping the distribution of TILs. In a dataset comprising grade 2-3 pure DCIS and DCIS adjacent to invasive cancer (adjacent DCIS), we found that pure DCIS cases had more TILs compared to adjacent DCIS. However, the colocalisation of TILs with DCIS ducts was significantly lower in pure DCIS compared to adjacent DCIS, which may suggest a more inflamed tissue ecology local to DCIS ducts in adjacent DCIS cases. Our study demonstrates that technological developments in deep convolutional neural networks and digital pathology can enable an automated morphological and microenvironmental analysis of DCIS, providing a new way to study differential immune ecology for individual ducts and identify new markers of progression.

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.

Integrating genomic features for non-invasive early lung cancer detection.

Radiologic screening of high-risk adults reduces lung-cancer-related mortality1,2; however, a small minority of eligible individuals undergo such screening in the United States3,4. The availability of blood-based tests could increase screening uptake. Here we introduce improvements to cancer personalized profiling by deep sequencing (CAPP-Seq)5, a method for the analysis of circulating tumour DNA (ctDNA), to better facilitate screening applications. We show that, although levels are very low in early-stage lung cancers, ctDNA is present prior to treatment in most patients and its presence is strongly prognostic. We also find that the majority of somatic mutations in the cell-free DNA (cfDNA) of patients with lung cancer and of risk-matched controls reflect clonal haematopoiesis and are non-recurrent. Compared with tumour-derived mutations, clonal haematopoiesis mutations occur on longer cfDNA fragments and lack mutational signatures that are associated with tobacco smoking. Integrating these findings with other molecular features, we develop and prospectively validate a machine-learning method termed 'lung cancer likelihood in plasma' (Lung-CLiP), which can robustly discriminate early-stage lung cancer patients from risk-matched controls. This approach achieves performance similar to that of tumour-informed ctDNA detection and enables tuning of assay specificity in order to facilitate distinct clinical applications. Our findings establish the potential of cfDNA for lung cancer screening and highlight the importance of risk-matching cases and controls in cfDNA-based screening studies.

Primary mammary angiosarcomas harbor frequent mutations in KDR and PIK3CA and show evidence of distinct pathogenesis.

Angiosarcoma (AS) is the most frequent primary sarcoma of the breast but nevertheless remains uncommon, accounting for <0.05% of breast malignancies. Secondary mammary AS arise following radiation therapy for breast cancer, in contrast to primary AS which occur sporadically. Essentially all show aggressive clinical behavior independent of histologic grade and most are treated by mastectomy. MYC amplification is frequently identified in radiation-induced AS but only rarely in primary mammary AS (PMAS). As a heterogeneous group, AS from various anatomic sites have been shown to harbor recurrent alterations in TP53, MAP kinase pathway genes, and genes involved in angiogenic signaling including KDR (VEGFR2) and PTPRB. In part due to its rarity, the pathogenesis of PMAS has not been fully characterized. In this study, we examined the clinical, pathologic, and genomic features of ten cases of PMAS, including one patient with bilateral disease. Recurrent genomic alterations were identified in KDR (70%), PIK3CA/PIK3R1 (70%), and PTPRB (30%), each at higher frequencies than reported in AS across all sites. Six tumors harbored a KDR p.T771R hotspot mutation, and all seven KDR-mutant cases showed evidence suggestive of biallelism (four with loss of heterozygosity and three with two aberrations). Of the seven tumors with PI3K alterations, six harbored pathogenic mutations other than in the canonical PIK3CA residues which are most frequent in breast cancer. Three AS were hypermutated (≥10 mutations/megabase (Mb)); hypermutation was seen concurrent with KDR or PIK3CA mutations. The patient with bilateral disease demonstrated shared alterations, indicative of contralateral metastasis. No MYC or TP53 aberrations were detected in this series. Immunohistochemistry for VEGFR2 was unable to discriminate between KDR-mutant tumors and benign vascular lesions of the breast. These findings highlight the underrecognized frequency of KDR and PIK3CA mutation in PMAS, and a significant subset with hypermutation, suggesting a pathogenesis distinct from other AS.