Single-cell sequencing reveals the landscape of the human brain metastatic microenvironment.

Brain metastases is the most common intracranial tumor and account for approximately 20% of all systematic cancer cases. It is a leading cause of death in advanced-stage cancer, resulting in a five-year overall survival rate below 10%. Therefore, there is a critical need to identify effective biomarkers that can support frequent surveillance and promote efficient drug guidance in brain metastasis. Recently, the remarkable breakthroughs in single-cell RNA-sequencing (scRNA-seq) technology have advanced our insights into the tumor microenvironment (TME) at single-cell resolution, which offers the potential to unravel the metastasis-related cellular crosstalk and provides the potential for improving therapeutic effects mediated by multifaceted cellular interactions within TME. In this study, we have applied scRNA-seq and profiled 10,896 cells collected from five brain tumor tissue samples originating from breast and lung cancers. Our analysis reveals the presence of various intratumoral components, including tumor cells, fibroblasts, myeloid cells, stromal cells expressing neural stem cell markers, as well as minor populations of oligodendrocytes and T cells. Interestingly, distinct cellular compositions are observed across different samples, indicating the influence of diverse cellular interactions on the infiltration patterns within the TME. Importantly, we identify tumor-associated fibroblasts in both our in-house dataset and external scRNA-seq datasets. These fibroblasts exhibit high expression of type I collagen genes, dominate cell-cell interactions within the TME via the type I collagen signaling axis, and facilitate the remodeling of the TME to a collagen-I-rich extracellular matrix similar to the original TME at primary sites. Additionally, we observe M1 activation in native microglial cells and infiltrated macrophages, which may contribute to a proinflammatory TME and the upregulation of collagen type I expression in fibroblasts. Furthermore, tumor cell-specific receptors exhibit a significant association with patient survival in both brain metastasis and native glioblastoma cases. Taken together, our comprehensive analyses identify type I collagen-secreting tumor-associated fibroblasts as key mediators in metastatic brain tumors and uncover tumor receptors that are potentially associated with patient survival. These discoveries provide potential biomarkers for effective therapeutic targets and intervention strategies.

Bulk and Single-Cell Profiling of Breast Tumors Identifies TREM-1 as a Dominant Immune Suppressive Marker Associated With Poor Outcomes.

BACKGROUND: Triggering receptor expressed on myeloid cells (TREM)-1 is a key mediator of innate immunity previously associated with the severity of inflammatory disorders, and more recently, the inferior survival of lung and liver cancer patients. Here, we investigated the prognostic impact and immunological correlates of TREM1 expression in breast tumors. METHODS: Breast tumor microarray and RNAseq expression profiles (n=4,364 tumors) were analyzed for associations between gene expression, tumor immune subtypes, distant metastasis-free survival (DMFS) and clinical response to neoadjuvant chemotherapy (NAC). Single-cell (sc)RNAseq was performed using the 10X Genomics platform. Statistical associations were assessed by logistic regression, Cox regression, Kaplan-Meier analysis, Spearman correlation, Student's t-test and Chi-square test. RESULTS: In pre-treatment biopsies, TREM1 and known TREM-1 inducible cytokines (IL1B, IL8) were discovered by a statistical ranking procedure as top genes for which high expression was associated with reduced response to NAC, but only in the context of immunologically "hot" tumors otherwise associated with a high NAC response rate. In surgical specimens, TREM1 expression varied among tumor molecular subtypes, with highest expression in the more aggressive subtypes (Basal-like, HER2-E). High TREM1 significantly and reproducibly associated with inferior distant metastasis-free survival (DMFS), independent of conventional prognostic markers. Notably, the association between high TREM1 and inferior DMFS was most prominent in the subset of immunogenic tumors that exhibited the immunologically hot phenotype and otherwise associated with superior DMFS. Further observations from bulk and single-cell RNAseq analyses indicated that TREM1 expression was significantly enriched in polymorphonuclear myeloid-derived suppressor cells (PMN-MDSCs) and M2-like macrophages, and correlated with downstream transcriptional targets of TREM-1 (IL8, IL-1B, IL6, MCP-1, SPP1, IL1RN, INHBA) which have been previously associated with pro-tumorigenic and immunosuppressive functions. CONCLUSIONS: Together, these findings indicate that increased TREM1 expression is prognostic of inferior breast cancer outcomes and may contribute to myeloid-mediated breast cancer progression and immune suppression.

Transcriptomic Features of T Cell-Barren Tumors Are Conserved Across Diverse Tumor Types.

Background: Understanding how tumors subvert immune destruction is essential to the development of cancer immunotherapies. New evidence suggests that tumors limit anti-tumor immunity by exploiting transcriptional programs that regulate intratumoral trafficking and accumulation of effector cells. Here, we investigated the gene expression profiles that distinguish immunologically "cold" and "hot" tumors across diverse tumor types. Methods: RNAseq profiles of tumors (n = 8,920) representing 23 solid tumor types were analyzed using immune gene signatures that quantify CD8+ T cell abundance. Genes and pathways associated with a low CD8+ T cell infiltration profile (CD8-Low) were identified by correlation, differential expression, and statistical ranking methods. Gene subsets were evaluated in immunotherapy treatment cohorts and functionally characterized in cell lines and mouse tumor models. Results: Among different cancer types, we observed highly significant overlap of genes enriched in CD8-Low tumors, which included known immunomodulatory genes (e.g., BMP7, CMTM4, KDM5B, RCOR2) and exhibited significant associations with Wnt signaling, neurogenesis, cell-cell junctions, lipid biosynthesis, epidermal development, and cancer-testis antigens. Analysis of mutually exclusive gene clusters demonstrated that different transcriptional programs may converge on the T cell-cold phenotype as well as predict for response and survival of patients to Nivo treatment. Furthermore, we confirmed that a top-ranking candidate belonging to the TGF-ß superfamily, BMP7, negatively regulates CD8+ T cell abundance in immunocompetent murine tumor models, with and without anti-PD-L1 treatment. Conclusions: This study presents the first evidence that solid tumors of diverse anatomical origin acquire conserved transcriptional alterations that may be operative in the T cell-cold state. Our findings demonstrate the potential clinical utility of CD8-Low tumor-associated genes for predicting patient immunotherapy outcomes and point to novel mechanisms with potential for broad therapeutic exploitation.

Dissecting intratumoral myeloid cell plasticity by single cell RNA-seq.

Tumor-infiltrating myeloid cells are the most abundant leukocyte population within tumors. Molecular cues from the tumor microenvironment promote the differentiation of immature myeloid cells toward an immunosuppressive phenotype. However, the in situ dynamics of the transcriptional reprogramming underlying this process are poorly understood. Therefore, we applied single cell RNA-seq (scRNA-seq) to computationally investigate the cellular composition and transcriptional dynamics of tumor and adjacent normal tissues from 4 early-stage non-small cell lung cancer (NSCLC) patients. Our scRNA-seq analyses identified 11 485 cells that varied in identity and gene expression traits between normal and tumor tissues. Among these, myeloid cell populations exhibited the most diverse changes between tumor and normal tissues, consistent with tumor-mediated reprogramming. Through trajectory analysis, we identified a differentiation path from CD14+ monocytes to M2 macrophages (monocyte-to-M2). This differentiation path was reproducible across patients, accompanied by increased expression of genes (eg, MRC1/CD206, MSR1/CD204, PPARG, TREM2) with significantly enriched functions (Oxidative phosphorylation and P53 pathway) and decreased expression of genes (eg, CXCL2, IL1B) with significantly enriched functions (TNF-α signaling via NF-κB and inflammatory response). Our analysis further identified a co-regulatory network implicating upstream transcription factors (JUN, NFKBIA) in monocyte-to-M2 differentiation, and activated ligand-receptor interactions (eg, SFTPA1-TLR2, ICAM1-ITGAM) suggesting intratumoral mechanisms whereby epithelial cells stimulate monocyte-to-M2 differentiation. Overall, our study identified the prevalent monocyte-to-M2 differentiation in NSCLC, accompanied by an intricate transcriptional reprogramming mediated by specific transcriptional activators and intercellular crosstalk involving ligand-receptor interactions.

Systems biology approach to studying proliferation-dependent prognostic subnetworks in breast cancer.

Tumor proliferative capacity is a major biological correlate of breast tumor metastatic potential. In this paper, we developed a systems approach to investigate associations among gene expression patterns, representative protein-protein interactions, and the potential for clinical metastases, to uncover novel survival-related subnetwork signatures as a function of tumor proliferative potential. Based on the statistical associations between gene expression patterns and patient outcomes, we identified three groups of survival prognostic subnetwork signatures (SPNs) corresponding to three proliferation levels. We discovered 8 SPNs in the high proliferation group, 8 SPNs in the intermediate proliferation group, and 6 SPNs in the low proliferation group. We observed little overlap of SPNs between the three proliferation groups. The enrichment analysis revealed that most SPNs were enriched in distinct signaling pathways and biological processes. The SPNs were validated on other cohorts of patients, and delivered high accuracy in the classification of metastatic vs non-metastatic breast tumors. Our findings indicate that certain biological networks underlying breast cancer metastasis differ in a proliferation-dependent manner. These networks, in combination, may form the basis of highly accurate prognostic classification models and may have clinical utility in guiding therapeutic options for patients.

Trk activation of the ERK1/2 kinase pathway stimulates intermediate chain phosphorylation and recruits cytoplasmic dynein to signaling endosomes for retrograde axonal transport.

The retrograde transport of Trk-containing endosomes from the axon to the cell body by cytoplasmic dynein is necessary for axonal and neuronal survival. We investigated the recruitment of dynein to signaling endosomes in rat embryonic neurons and PC12 cells. We identified a novel phosphoserine on the dynein intermediate chains (ICs), and we observed a time-dependent neurotrophin-stimulated increase in intermediate chain phosphorylation on this site in both cell types. Pharmacological studies, overexpression of constitutively active MAP kinase kinase, and an in vitro assay with recombinant proteins demonstrated that the intermediate chains are phosphorylated by the MAP kinase ERK1/2, extracellular signal-regulated kinase, a major downstream effector of Trk. Live cell imaging with fluorescently tagged IC mutants demonstrated that the dephosphomimic mutants had significantly reduced colocalization with Trk and Rab7, but not a mitochondrial marker. The phosphorylated intermediate chains were enriched on immunoaffinity-purified Trk-containing organelles. Inhibition of ERK reduced the amount of phospho-IC and the total amount of dynein that copurified with the signaling endosomes. In addition, inhibition of ERK1/2 reduced the motility of Rab7- and TrkB-containing endosomes and the extent of their colocalization with dynein in axons. NGF-dependent survival of sympathetic neurons was significantly reduced by the overexpression of the dephosphomimic mutant IC-1B-S80A, but not WT IC-1B, further demonstrating the functional significance of phosphorylation on this site. These results demonstrate that neurotrophin binding to Trk initiates the recruitment of cytoplasmic dynein to signaling endosomes through ERK1/2 phosphorylation of intermediate chains for their subsequent retrograde transport in axons.

ISG15 disrupts cytoskeletal architecture and promotes motility in human breast cancer cells.

The interferon-stimulated gene 15 (ISG15) pathway is highly elevated in breast cancer; however, very little is known about how the ISG15 pathway contributes to breast tumorigenesis. In the current study, using the gene disruption approach, we demonstrate that both ISG15 and UbcH8 (ISG15-specific conjugating enzyme) disrupt F-actin architecture and formation of focal adhesions in ZR-75-1 breast cancer cells. In addition, ISG15 and UbcH8 promote breast cancer cell migration. We also demonstrate that ISG15 inhibits ubiquitin/26S proteasome-mediated turnover of proteins implicated in tumor cell motility, invasion and metastasis. Together, our results suggest that the aberrant activation of the ISG15 pathway confers a motile phenotype to breast cancer cells by disrupting cell architecture and stabilizing proteins involved in cell motility, invasion and metastasis. Because the cellular architecture is conserved and the ISG15 pathway is constitutively activated in tumor cells of different lineages, it is reasonable to assume that our observations in breast cancer must hold true for many other tumors.

Extracellular signal-regulated kinase promotes Rho-dependent focal adhesion formation by suppressing p190A RhoGAP.

Cell migration is critical for normal development and for pathological processes including cancer cell metastasis. Dynamic remodeling of focal adhesions and the actin cytoskeleton are crucial determinants of cell motility. The Rho family and the mitogen-activated protein kinase (MAPK) module consisting of MEK-extracellular signal-regulated kinase (ERK) are important regulators of these processes, but mechanisms for the integration of these signals during spreading and motility are incompletely understood. Here we show that ERK activity is required for fibronectin-stimulated Rho-GTP loading, Rho-kinase function, and the maturation of focal adhesions in spreading cells. We identify p190A RhoGAP as a major target for ERK signaling in adhesion assembly and identify roles for ERK phosphorylation of the C terminus in p190A localization and activity. These observations reveal a novel role for ERK signaling in adhesion assembly in addition to its established role in adhesion disassembly.

Differential requirement for MEK Partner 1 in DU145 prostate cancer cell migration.

ERK signaling regulates focal adhesion disassembly during cell movement, and increased ERK signaling frequently contributes to enhanced motility of human tumor cells. We previously found that the ERK scaffold MEK Partner 1 (MP1) is required for focal adhesion disassembly in fibroblasts. Here we test the hypothesis that MP1-dependent ERK signaling regulates motility of DU145 prostate cancer cells. We find that MP1 is required for motility on fibronectin, but not for motility stimulated by serum or EGF. Surprisingly, MP1 appears not to function through its known binding partners MEK1 or PAK1, suggesting the existence of a novel pathway by which MP1 can regulate motility on fibronectin. MP1 may function by regulating the stability or expression of paxillin, a key regulator of motility.

Scaffold mediated regulation of MAPK signaling and cytoskeletal dynamics: a perspective.

Cell migration is critical for many physiological processes and is often misregulated in developmental disorders and pathological conditions including cancer and neurodegeneration. MAPK signaling and the Rho family of proteins are known regulators of cell migration that exert their influence on cellular cytoskeleton during cell adhesion and migration. Here we review data supporting the view that localized ERK signaling mediated through recently identified scaffold proteins may regulate cell migration.

In vivo membrane trafficking role for an insect N-ethylmaleimide-sensitive factor which is developmentally regulated in endocrine cells.

The hexameric ATPase, N-ethylmaleimide-sensitive factor (NSF) is implicated in the release of neurotransmitters and in mediating fusion between intracellular membranes. Due to the conservation of proteins in constitutive and regulated membrane fusion reactions, NSF and its downstream targets have been predicted also to participate in fusion reactions underlying endocrine function, but there is little experimental evidence to support such a role for NSF in insect neuroendocrine secretion. Here we have characterized the NSF orthologue (MsNSF) from the endocrine model for development Manduca sexta. MsNSF is developmentally regulated in endocrine organs of the protocerebral complex. Enrichment of MsNSF in corpora cardiaca (CC) and not in corpora allata (CA) indicates that it might play a preferential role in releasing hormones produced in CC. Endocrine/paracrine cells of the enteric system in M. sexta exhibit selective MsNSF enrichment. Together the data point to a more selective participation of MsNSF in development of M. sexta by its involvement in a subset of factors, whereas other as-yet-unidentified homolog(s) might regulate secretion from CA and a large set of endocrine/paracrine cells. We further characterized the in vivo role of MsNSF by heterologous expression. In contrast to vertebrate NSF, MsNSF is functional in yeast membrane fusion in vivo. MsNSF rectifies defects in SEC18 (yeast NSF homologue) at nearly all discernible steps where Sec18p has been implicated in the biosynthetic route. This underscores the utility of our approach to delineate functional roles for proteins from systems that are not currently amenable to in vitro reconstitution.

Heliothis virescens and Manduca sexta lipid rafts are involved in Cry1A toxin binding to the midgut epithelium and subsequent pore formation.

Lipid rafts are characterized by their insolubility in nonionic detergents such as Triton X-100 at 4 degrees C. They have been studied in mammals, where they play critical roles in protein sorting and signal transduction. To understand the potential role of lipid rafts in lepidopteran insects, we isolated and analyzed the protein and lipid components of these lipid raft microdomains from the midgut epithelial membrane of Heliothis virescens and Manduca sexta. Like their mammalian counterparts, H. virescens and M. sexta lipid rafts are enriched in cholesterol, sphingolipids, and glycosylphosphatidylinositol-anchored proteins. In H. virescens and M. sexta, pretreatment of membranes with the cholesterol-depleting reagent saponin and methyl-beta-cyclodextrin differentially disrupted the formation of lipid rafts, indicating an important role for cholesterol in lepidopteran lipid rafts structure. We showed that several putative Bacillus thuringiensis Cry1A receptors, including the 120- and 170-kDa aminopeptidases from H. virescens and the 120-kDa aminopeptidase from M. sexta, were preferentially partitioned into lipid rafts. Additionally, the leucine aminopeptidase activity was enriched approximately 2-3-fold in these rafts compared with brush border membrane vesicles. We also demonstrated that Cry1A toxins were associated with lipid rafts, and that lipid raft integrity was essential for in vitro Cry1Ab pore forming activity. Our study strongly suggests that these microdomains might be involved in Cry1A toxin aggregation and pore formation.