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.

Atlas of clinically distinct cell states and ecosystems across human solid tumors.

Determining how cells vary with their local signaling environment and organize into distinct cellular communities is critical for understanding processes as diverse as development, aging, and cancer. Here we introduce EcoTyper, a machine learning framework for large-scale identification and validation of cell states and multicellular communities from bulk, single-cell, and spatially resolved gene expression data. When applied to 12 major cell lineages across 16 types of human carcinoma, EcoTyper identified 69 transcriptionally defined cell states. Most states were specific to neoplastic tissue, ubiquitous across tumor types, and significantly prognostic. By analyzing cell-state co-occurrence patterns, we discovered ten clinically distinct multicellular communities with unexpectedly strong conservation, including three with myeloid and stromal elements linked to adverse survival, one enriched in normal tissue, and two associated with early cancer development. This study elucidates fundamental units of cellular organization in human carcinoma and provides a framework for large-scale profiling of cellular ecosystems in any tissue.

A Structured Tumor-Immune Microenvironment in Triple Negative Breast Cancer Revealed by Multiplexed Ion Beam Imaging.

The immune system is critical in modulating cancer progression, but knowledge of immune composition, phenotype, and interactions with tumor is limited. We used multiplexed ion beam imaging by time-of-flight (MIBI-TOF) to simultaneously quantify in situ expression of 36 proteins covering identity, function, and immune regulation at sub-cellular resolution in 41 triple-negative breast cancer patients. Multi-step processing, including deep-learning-based segmentation, revealed variability in the composition of tumor-immune populations across individuals, reconciled by overall immune infiltration and enriched co-occurrence of immune subpopulations and checkpoint expression. Spatial enrichment analysis showed immune mixed and compartmentalized tumors, coinciding with expression of PD1, PD-L1, and IDO in a cell-type- and location-specific manner. Ordered immune structures along the tumor-immune border were associated with compartmentalization and linked to survival. These data demonstrate organization in the tumor-immune microenvironment that is structured in cellular composition, spatial arrangement, and regulatory-protein expression and provide a framework to apply multiplexed imaging to immune oncology.

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.

Diagnostic utility of FOXO1 immunohistochemistry for rhabdomyosarcoma classification.

AIMS: Rhabdomyosarcomas currently are classified into one of four subtypes (alveolar, embryonal, spindle cell/sclerosing, or pleomorphic) according to their morphological, immunohistochemical, and molecular genetic features. The alveolar subtype is characterised by a recurrent translocation involving PAX3 or PAX7 and FOXO1; identification of this translocation is important for appropriate classification and prognostication. In this study, we aimed to explore the diagnostic utility of FOXO1 immunohistochemistry for rhabdomyosarcoma classification. METHODS/RESULTS: A monoclonal antibody targeting a FOXO1 epitope retained in the fusion oncoprotein was used to study 105 rhabdomyosarcomas. FOXO1 was positive for expression by immunohistochemistry in all 25 alveolar rhabdomyosarcomas, with 84% showing diffuse expression in greater than 90% of neoplastic cells; the remainder of alveolar rhabdomyosarcomas displayed at least moderate staining in a minimum of 60% of lesional cells. Apart from three spindle cell rhabdomyosarcomas showing heterogeneous nuclear immunoreactivity in 40-80% of tumour cells, the 80 cases of embryonal, pleomorphic, and spindle cell/sclerosing rhabdomyosarcoma were negative for FOXO1 expression (96.3% specific) when using a threshold of nuclear staining in 20% of neoplastic cells to determine positivity. Variable cytoplasmic staining was present in a fraction of all rhabdomyosarcoma subtypes. Nonneoplastic lymphocytes, endothelial cells, and Schwann cells also showed variably intense nuclear anti-FOXO1 immunoreactivity. CONCLUSION: Taken together, our findings suggest that FOXO1 immunohistochemistry is a highly sensitive and relatively specific surrogate marker of the PAX3/7::FOXO1 fusion oncoprotein in rhabdomyosarcoma. Cytoplasmic immunoreactivity, expression in nonneoplastic tissues, and limited nuclear staining of nonalveolar rhabdomyosarcomas represent potential pitfalls in interpretation.

Reproducible, high-dimensional imaging in archival human tissue by multiplexed ion beam imaging by time-of-flight (MIBI-TOF).

Multiplexed ion beam imaging by time-of-flight (MIBI-TOF) is a form of mass spectrometry imaging that uses metal labeled antibodies and secondary ion mass spectrometry to image dozens of proteins simultaneously in the same tissue section. Working with the National Cancer Institute's (NCI) Cancer Immune Monitoring and Analysis Centers (CIMAC), we undertook a validation study, assessing concordance across a dozen serial sections of a tissue microarray of 21 samples that were independently processed and imaged by MIBI-TOF or single-plex immunohistochemistry (IHC) over 12 days. Pixel-level features were highly concordant across all 16 targets assessed in both staining intensity (R2 = 0.94 ± 0.04) and frequency (R2 = 0.95 ± 0.04). Comparison to digitized, single-plex IHC on adjacent serial sections revealed similar concordance (R2 = 0.85 ± 0.08) as well. Lastly, automated segmentation and clustering of eight cell populations found that cell frequencies between serial sections yielded an average correlation of R2 = 0.94 ± 0.05. Taken together, we demonstrate that MIBI-TOF, with well-vetted reagents and automated analysis, can generate consistent and quantitative annotations of clinically relevant cell states in archival human tissue, and more broadly, present a scalable framework for benchmarking multiplexed IHC approaches.

Transcriptome and genome evolution during HER2-amplified breast neoplasia.

BACKGROUND: The acquisition of oncogenic drivers is a critical feature of cancer progression. For some carcinomas, it is clear that certain genetic drivers occur early in neoplasia and others late. Why these drivers are selected and how these changes alter the neoplasia's fitness is less understood. METHODS: Here we use spatially oriented genomic approaches to identify transcriptomic and genetic changes at the single-duct level within precursor neoplasia associated with invasive breast cancer. We study HER2 amplification in ductal carcinoma in situ (DCIS) as an event that can be both quantified and spatially located via fluorescence in situ hybridization (FISH) and immunohistochemistry on fixed paraffin-embedded tissue. RESULTS: By combining the HER2-FISH with the laser capture microdissection (LCM) Smart-3SEQ method, we found that HER2 amplification in DCIS alters the transcriptomic profiles and increases diversity of copy number variations (CNVs). Particularly, interferon signaling pathway is activated by HER2 amplification in DCIS, which may provide a prolonged interferon signaling activation in HER2-positive breast cancer. Multiple subclones of HER2-amplified DCIS with distinct CNV profiles are observed, suggesting that multiple events occurred for the acquisition of HER2 amplification. Notably, DCIS acquires key transcriptomic changes and CNV events prior to HER2 amplification, suggesting that pre-amplified DCIS may create a cellular state primed to gain HER2 amplification for growth advantage. CONCLUSION: By using genomic methods that are spatially oriented, this study identifies several features that appear to generate insights into neoplastic progression in precancer lesions at a single-duct level.

Genomic landscape of ductal carcinoma in situ and association with progression.

PURPOSE: The detection rate of breast ductal carcinoma in situ (DCIS) has increased significantly, raising the concern that DCIS is overdiagnosed and overtreated. Therefore, there is an unmet clinical need to better predict the risk of progression among DCIS patients. Our hypothesis is that by combining molecular signatures with clinicopathologic features, we can elucidate the biology of breast cancer progression, and risk-stratify patients with DCIS. METHODS: Targeted exon sequencing with a custom panel of 223 genes/regions was performed for 125 DCIS cases. Among them, 60 were from cases having concurrent or subsequent invasive breast cancer (IBC) (DCIS + IBC group), and 65 from cases with no IBC development over a median follow-up of 13 years (DCIS-only group). Copy number alterations in chromosome 1q32, 8q24, and 11q13 were analyzed using fluorescence in situ hybridization (FISH). Multivariable logistic regression models were fit to the outcome of DCIS progression to IBC as functions of demographic and clinical features. RESULTS: We observed recurrent variants of known IBC-related mutations, and the most commonly mutated genes in DCIS were PIK3CA (34.4%) and TP53 (18.4%). There was an inverse association between PIK3CA kinase domain mutations and progression (Odds Ratio [OR] 10.2, p < 0.05). Copy number variations in 1q32 and 8q24 were associated with progression (OR 9.3 and 46, respectively; both p < 0.05). CONCLUSIONS: PIK3CA kinase domain mutations and the absence of copy number gains in DCIS are protective against progression to IBC. These results may guide efforts to distinguish low-risk from high-risk DCIS.

Discovery of recurrent structural variants in nasopharyngeal carcinoma.

We present the discovery of genes recurrently involved in structural variation in nasopharyngeal carcinoma (NPC) and the identification of a novel type of somatic structural variant. We identified the variants with high complexity mate-pair libraries and a novel computational algorithm specifically designed for tumor-normal comparisons, SMASH. SMASH combines signals from split reads and mate-pair discordance to detect somatic structural variants. We demonstrate a >90% validation rate and a breakpoint reconstruction accuracy of 3 bp by Sanger sequencing. Our approach identified three in-frame gene fusions (YAP1-MAML2, PTPLB-RSRC1, and SP3-PTK2) that had strong levels of expression in corresponding NPC tissues. We found two cases of a novel type of structural variant, which we call "coupled inversion," one of which produced the YAP1-MAML2 fusion. To investigate whether the identified fusion genes are recurrent, we performed fluorescent in situ hybridization (FISH) to screen 196 independent NPC cases. We observed recurrent rearrangements of MAML2 (three cases), PTK2 (six cases), and SP3 (two cases), corresponding to a combined rate of structural variation recurrence of 6% among tested NPC tissues.

Genome evolution during progression to breast cancer.

Cancer evolution involves cycles of genomic damage, epigenetic deregulation, and increased cellular proliferation that eventually culminate in the carcinoma phenotype. Early neoplasias, which are often found concurrently with carcinomas and are histologically distinguishable from normal breast tissue, are less advanced in phenotype than carcinomas and are thought to represent precursor stages. To elucidate their role in cancer evolution we performed comparative whole-genome sequencing of early neoplasias, matched normal tissue, and carcinomas from six patients, for a total of 31 samples. By using somatic mutations as lineage markers we built trees that relate the tissue samples within each patient. On the basis of these lineage trees we inferred the order, timing, and rates of genomic events. In four out of six cases, an early neoplasia and the carcinoma share a mutated common ancestor with recurring aneuploidies, and in all six cases evolution accelerated in the carcinoma lineage. Transition spectra of somatic mutations are stable and consistent across cases, suggesting that accumulation of somatic mutations is a result of increased ancestral cell division rather than specific mutational mechanisms. In contrast to highly advanced tumors that are the focus of much of the current cancer genome sequencing, neither the early neoplasia genomes nor the carcinomas are enriched with potentially functional somatic point mutations. Aneuploidies that occur in common ancestors of neoplastic and tumor cells are the earliest events that affect a large number of genes and may predispose breast tissue to eventual development of invasive carcinoma.

Breakpoint analysis of transcriptional and genomic profiles uncovers novel gene fusions spanning multiple human cancer types.

Gene fusions, like BCR/ABL1 in chronic myelogenous leukemia, have long been recognized in hematologic and mesenchymal malignancies. The recent finding of gene fusions in prostate and lung cancers has motivated the search for pathogenic gene fusions in other malignancies. Here, we developed a "breakpoint analysis" pipeline to discover candidate gene fusions by tell-tale transcript level or genomic DNA copy number transitions occurring within genes. Mining data from 974 diverse cancer samples, we identified 198 candidate fusions involving annotated cancer genes. From these, we validated and further characterized novel gene fusions involving ROS1 tyrosine kinase in angiosarcoma (CEP85L/ROS1), SLC1A2 glutamate transporter in colon cancer (APIP/SLC1A2), RAF1 kinase in pancreatic cancer (ATG7/RAF1) and anaplastic astrocytoma (BCL6/RAF1), EWSR1 in melanoma (EWSR1/CREM), CDK6 kinase in T-cell acute lymphoblastic leukemia (FAM133B/CDK6), and CLTC in breast cancer (CLTC/VMP1). Notably, while these fusions involved known cancer genes, all occurred with novel fusion partners and in previously unreported cancer types. Moreover, several constituted druggable targets (including kinases), with therapeutic implications for their respective malignancies. Lastly, breakpoint analysis identified new cell line models for known rearrangements, including EGFRvIII and FIP1L1/PDGFRA. Taken together, we provide a robust approach for gene fusion discovery, and our results highlight a more widespread role of fusion genes in cancer pathogenesis.

Chromosomal copy number alterations for associations of ductal carcinoma in situ with invasive breast cancer.

INTRODUCTION: Screening mammography has contributed to a significant increase in the diagnosis of ductal carcinoma in situ (DCIS), raising concerns about overdiagnosis and overtreatment. Building on prior observations from lineage evolution analysis, we examined whether measuring genomic features of DCIS would predict association with invasive breast carcinoma (IBC). The long-term goal is to enhance standard clinicopathologic measures of low- versus high-risk DCIS and to enable risk-appropriate treatment. METHODS: We studied three common chromosomal copy number alterations (CNA) in IBC and designed fluorescence in situ hybridization-based assay to measure copy number at these loci in DCIS samples. Clinicopathologic data were extracted from the electronic medical records of Stanford Cancer Institute and linked to demographic data from the population-based California Cancer Registry; results were integrated with data from tissue microarrays of specimens containing DCIS that did not develop IBC versus DCIS with concurrent IBC. Multivariable logistic regression analysis was performed to describe associations of CNAs with these two groups of DCIS. RESULTS: We examined 271 patients with DCIS (120 that did not develop IBC and 151 with concurrent IBC) for the presence of 1q, 8q24 and 11q13 copy number gains. Compared to DCIS-only patients, patients with concurrent IBC had higher frequencies of CNAs in their DCIS samples. On multivariable analysis with conventional clinicopathologic features, the copy number gains were significantly associated with concurrent IBC. The state of two of the three copy number gains in DCIS was associated with a risk of IBC that was 9.07 times that of no copy number gains, and the presence of gains at all three genomic loci in DCIS was associated with a more than 17-fold risk (P = 0.0013). CONCLUSIONS: CNAs have the potential to improve the identification of high-risk DCIS, defined by presence of concurrent IBC. Expanding and validating this approach in both additional cross-sectional and longitudinal cohorts may enable improved risk stratification and risk-appropriate treatment in DCIS.

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.

CD142 Identifies Neoplastic Desmoid Tumor Cells, Uncovering Interactions Between Neoplastic and Stromal Cells That Drive Proliferation.

The interaction between neoplastic and stromal cells within a tumor mass plays an important role in cancer biology. However, it is challenging to distinguish between tumor and stromal cells in mesenchymal tumors because lineage-specific cell surface markers typically used in other cancers do not distinguish between the different cell subpopulations. Desmoid tumors consist of mesenchymal fibroblast-like cells driven by mutations stabilizing beta-catenin. Here we aimed to identify surface markers that can distinguish mutant cells from stromal cells to study tumor-stroma interactions. We analyzed colonies derived from single cells from human desmoid tumors using a high-throughput surface antigen screen, to characterize the mutant and nonmutant cells. We found that CD142 is highly expressed by the mutant cell populations and correlates with beta-catenin activity. CD142-based cell sorting isolated the mutant population from heterogeneous samples, including one where no mutation was previously detected by traditional Sanger sequencing. We then studied the secretome of mutant and nonmutant fibroblastic cells. PTX3 is one stroma-derived secreted factor that increases mutant cell proliferation via STAT6 activation. These data demonstrate a sensitive method to quantify and distinguish neoplastic from stromal cells in mesenchymal tumors. It identifies proteins secreted by nonmutant cells that regulate mutant cell proliferation that could be therapeutically. Significance: Distinguishing between neoplastic (tumor) and non-neoplastic (stromal) cells within mesenchymal tumors is particularly challenging, because lineage-specific cell surface markers typically used in other cancers do not differentiate between the different cell subpopulations. Here, we developed a strategy combining clonal expansion with surface proteome profiling to identify markers for quantifying and isolating mutant and nonmutant cell subpopulations in desmoid tumors, and to study their interactions via soluble factors.

Distinctive contact between CD34+ hematopoietic progenitors and CXCL12+ CD271+ mesenchymal stromal cells in benign and myelodysplastic bone marrow.

Mesenchymal stromal cells (MSCs) support hematopoiesis and are cytogenetically and functionally abnormal in myelodysplastic syndrome (MDS), implying a possible pathophysiologic role in MDS and potential utility as a diagnostic or risk-stratifying tool. We have analyzed putative MSC markers and their relationship to CD34+ hematopoietic stem/progenitor cells (HSPCs) within intact human bone marrow in paraffin-embedded bone marrow core biopsies of benign, MDS and leukemic (AML) marrows using tissue microarrays to facilitate scanning, image analysis and quantitation. We found that CD271+, ALP+ MSCs formed an extensive branching perivascular, periosteal and parenchymal network. Nestin was brightly positive in capillary/arteriolar endothelium and occasional subendothelial cells, whereas CD146 was most brightly expressed in SMA+ vascular smooth muscle/pericytes. CD271+ MSCs were distinct by double immunofluorescence from CD163+ macrophages and were in close contact with but distinct from brightly nestin+ and from brightly CD146+ vascular elements. Double immunofluorescence revealed an intimate spatial relationship between CD34+ HSPCs and CD271+ MSCs; remarkably, 86% of CD34+ HSPCs were in direct contact with CD271+ MSCs across benign, MDS and AML marrows, predominantly in a perivascular distribution. Expression of the intercrine chemokine CXCL12 was strong in the vasculature in both benign and neoplastic marrow, but was also present in extravascular parenchymal cells, particularly in MDS specimens. We identified these parenchymal cells as MSCs by ALP/CXCL12 and CD271/CXCL12 double immunofluorescence. The area covered by CXCL12+ ALP+ MSCs was significantly greater in MDS compared with benign and AML marrow (P=0.021, Kruskal-Wallis test). The preservation of direct CD271+ MSC/CD34+ HSPC contact across benign and neoplastic marrow suggests a physiologically important role for the CD271+ MSC/CD34+ HSPC relationship and possible abnormal exposure of CD34+ HSPCs to increased MSC CXCL12 expression in MDS.

CSF1 expression in nongynecological leiomyosarcoma is associated with increased tumor angiogenesis.

Leiomyosarcoma (LMS) is a malignant tumor of smooth muscle cells for which few effective therapies exist. A subset of LMS cases express macrophage colony-stimulating factor (CSF1) and the resultant tumor-associated macrophage (TAM) infiltration predicts poor clinical outcome. Further, TAMs have been shown to increase tumor angiogenesis. Here, we analyzed 149 LMS cases by immunohistochemistry for vascular marker CD34 and show that high microvessel density (MVD) in nongynecological LMS cases significantly predicts poor patient outcome. The majority of high MVD cases were also CSF1-positive, and when combining high MVD with CSF1 expression, an even stronger prognostic correlation with patient outcome was obtained. Gene expression profiling revealed that MVD has a stronger correlation with CSF1 expression than with expression of vascular endothelial growth factor isoforms, which have traditionally been used as markers of angiogenesis and as anti-angiogenic therapeutic targets. Finally, patterns of CSF1 expression and TAM recruitment remained consistent between primary tumors and their metastases, and between primary tumors and those grown as xenografts in mice, highlighting the stability of these features to the biology of LMS tumors. Together, these findings suggest an important role for CSF1 and the resulting TAM infiltration in the pathological neovascularization of LMS tumors and provide a rationale for CSF1-targeted therapies in LMS.

GRM1 Immunohistochemistry Distinguishes Chondromyxoid Fibroma From its Histologic Mimics.

Chondromyxoid fibroma (CMF) is a rare benign bone neoplasm that manifests histologically as a lobular proliferation of stellate to spindle-shaped cells in a myxoid background, exhibiting morphologic overlap with other cartilaginous and myxoid tumors of bone. CMF is characterized by recurrent genetic rearrangements that place the glutamate receptor gene GRM1 under the regulatory control of a constitutively active promoter, leading to increased gene expression. Here, we explore the diagnostic utility of GRM1 immunohistochemistry as a surrogate marker for GRM1 rearrangement using a commercially available monoclonal antibody in a study of 230 tumors, including 30 CMF cases represented by 35 specimens. GRM1 was positive by immunohistochemistry in 97% of CMF specimens (34/35), exhibiting moderate to strong staining in more than 50% of neoplastic cells; staining was diffuse (>95% of cells) in 25 specimens (71%). Among the 9 CMF specimens with documented exposure to acid decalcification, 4 (44%) exhibited diffuse immunoreactivity (>95%) for GRM1, whereas all 15 CMF specimens (100%) with lack of exposure to decalcification reagents were diffusely immunoreactive ( P =0.003). High GRM1 expression at the RNA level was previously observed by quantitative reverse transcription polymerase chain reaction in 9 CMF cases that were also positive by immunohistochemistry; low GRM1 expression was observed by quantitative reverse transcription polymerase chain reaction in the single case of CMF that was negative by immunohistochemistry. GRM1 immunohistochemistry was negative (<5%) in histologic mimics of CMF, including conventional chondrosarcoma, enchondroma, chondroblastoma, clear cell chondrosarcoma, giant cell tumor of the bone, fibrous dysplasia, chondroblastic osteosarcoma, myoepithelial tumor, primary aneurysmal bone cyst, brown tumor, phosphaturic mesenchymal tumor, CMF-like osteosarcoma, and extraskeletal myxoid chondrosarcoma. These results indicate that GRM1 immunohistochemistry may have utility in distinguishing CMF from its histologic mimics.

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.

Multiomics Analysis of Spatially Distinct Stromal Cells Reveals Tumor-Induced O-Glycosylation of the CDK4-pRB Axis in Fibroblasts at the Invasive Tumor Edge.

The invasive leading edge represents a potential gateway for tumor metastasis. The role of fibroblasts from the tumor edge in promoting cancer invasion and metastasis has not been comprehensively elucidated. We hypothesize that cross-talk between tumor and stromal cells within the tumor microenvironment results in activation of key biological pathways depending on their position in the tumor (edge vs. core). Here we highlight phenotypic differences between tumor-adjacent-fibroblasts (TAF) from the invasive edge and tumor core fibroblasts from the tumor core, established from human lung adenocarcinomas. A multiomics approach that includes genomics, proteomics, and O-glycoproteomics was used to characterize cross-talk between TAFs and cancer cells. These analyses showed that O-glycosylation, an essential posttranslational modification resulting from sugar metabolism, alters key biological pathways including the cyclin-dependent kinase 4 (CDK4) and phosphorylated retinoblastoma protein axis in the stroma and indirectly modulates proinvasive features of cancer cells. In summary, the O-glycoproteome represents a new consideration for important biological processes involved in tumor-stroma cross-talk and a potential avenue to improve the anticancer efficacy of CDK4 inhibitors. SIGNIFICANCE: A multiomics analysis of spatially distinct fibroblasts establishes the importance of the stromal O-glycoproteome in tumor-stroma interactions at the leading edge and provides potential strategies to improve cancer treatment. See related commentary by De Wever, p. 537.

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.