Multi-scale signaling and tumor evolution in high-grade gliomas.

Although genomic anomalies in glioblastoma (GBM) have been well studied for over a decade, its 5-year survival rate remains lower than 5%. We seek to expand the molecular landscape of high-grade glioma, composed of IDH-wildtype GBM and IDH-mutant grade 4 astrocytoma, by integrating proteomic, metabolomic, lipidomic, and post-translational modifications (PTMs) with genomic and transcriptomic measurements to uncover multi-scale regulatory interactions governing tumor development and evolution. Applying 14 proteogenomic and metabolomic platforms to 228 tumors (212 GBM and 16 grade 4 IDH-mutant astrocytoma), including 28 at recurrence, plus 18 normal brain samples and 14 brain metastases as comparators, reveals heterogeneous upstream alterations converging on common downstream events at the proteomic and metabolomic levels and changes in protein-protein interactions and glycosylation site occupancy at recurrence. Recurrent genetic alterations and phosphorylation events on PTPN11 map to important regulatory domains in three dimensions, suggesting a central role for PTPN11 signaling across high-grade gliomas.

Label-free single-vesicle based surface enhanced Raman spectroscopy: A robust approach for investigating the biomolecular composition of small extracellular vesicles.

Small extracellular vesicles (sEVs) are cell-released vesicles ranging from 30-150nm in size. They have garnered increasing attention because of their potential for both the diagnosis and treatment of disease. The diversity of sEVs derives from their biological composition and cargo content. Currently, the isolation of sEV subpopulations is primarily based on bio-physical and affinity-based approaches. Since a standardized definition for sEV subpopulations is yet to be fully established, it is important to further investigate the correlation between the biomolecular composition of sEVs and their physical properties. In this study, we employed a platform combining single-vesicle surface-enhanced Raman spectroscopy (SERS) and machine learning to examine individual sEVs isolated by size-exclusion chromatography (SEC). The biomolecular composition of each vesicle examined was reflected by its corresponding SERS spectral features (biomolecular "fingerprints"), with their roots in the composition of their collective Raman-active bonds. Origins of the SERS spectral features were validated through a comparative analysis between SERS and mass spectrometry (MS). SERS fingerprinting of individual vesicles was effective in overcoming the challenges posed by EV population averaging, allowing for the possibility of analyzing the variations in biomolecular composition between the vesicles of similar and/or different sizes. Using this approach, we uncovered that each of the size-based fractions of sEVs contained particles with predominantly similar SERS spectral features. Indeed, more than 84% of the vesicles residing within a particular group were clearly distinguishable from that of the other EV sub-populations, despite some spectral variations within each sub-population. Our results suggest the possibility that size-based EV fractionation methods produce samples where similarly eluted sEVs are correlated with their respective biochemical contents, as reflected by their SERS spectra. Our findings therefore highlight the possibility that the biogenesis and respective biological functionalities of the various sEV fractions may be inherently different.

The TLK-ASF1 histone chaperone pathway plays a critical role in IL-1ß-mediated AML progression.

ABSTRACT: Identifying and targeting microenvironment-driven pathways that are active across acute myeloid leukemia (AML) genetic subtypes should allow the development of more broadly effective therapies. The proinflammatory cytokine interleukin-1ß (IL-1ß) is abundant in the AML microenvironment and promotes leukemic growth. Through RNA-sequencing analysis, we identify that IL-1ß-upregulated ASF1B (antisilencing function-1B), a histone chaperone, in AML progenitors compared with healthy progenitors. ASF1B, along with its paralogous protein ASF1A, recruits H3-H4 histones onto the replication fork during S-phase, a process regulated by Tousled-like kinase 1 and 2 (TLKs). Although ASF1s and TLKs are known to be overexpressed in multiple solid tumors and associated with poor prognosis, their functional roles in hematopoiesis and inflammation-driven leukemia remain unexplored. In this study, we identify that ASF1s and TLKs are overexpressed in multiple genetic subtypes of AML. We demonstrate that depletion of ASF1s significantly reduces leukemic cell growth in both in vitro and in vivo models using human cells. Using a murine model, we show that overexpression of ASF1B accelerates leukemia progression. Moreover, Asf1b or Tlk2 deletion delayed leukemia progression, whereas these proteins are dispensable for normal hematopoiesis. Through proteomics and phosphoproteomics analyses, we uncover that the TLK-ASF1 pathway promotes leukemogenesis by affecting the cell cycle and DNA damage pathways. Collectively, our findings identify the TLK1-ASF1 pathway as a novel mediator of inflammatory signaling and a promising therapeutic target for AML treatment across diverse genetic subtypes. Selective inhibition of this pathway offers potential opportunities to intervene effectively, address intratumoral heterogeneity, and ultimately improve clinical outcomes in AML.

PTPN2 Regulates Metabolic Flux to Affect ß-Cell Susceptibility to Inflammatory Stress.

Protein tyrosine phosphatase N2 (PTPN2) is a type 1 diabetes (T1D) candidate gene identified from human genome-wide association studies. PTPN2 is highly expressed in human and murine islets and becomes elevated upon inflammation and models of T1D, suggesting that PTPN2 may be important for ß-cell survival in the context of T1D. To test whether PTPN2 contributed to ß-cell dysfunction in an inflammatory environment, we generated a ß-cell-specific deletion of Ptpn2 in mice (PTPN2-ß knockout [ßKO]). Whereas unstressed animals exhibited normal metabolic profiles, low- and high-dose streptozotocin-treated PTPN2-ßKO mice displayed hyperglycemia and accelerated death, respectively. Furthermore, cytokine-treated Ptpn2-KO islets resulted in impaired glucose-stimulated insulin secretion, mitochondrial defects, and reduced glucose-induced metabolic flux, suggesting ß-cells lacking Ptpn2 are more susceptible to inflammatory stress associated with T1D due to maladaptive metabolic fitness. Consistent with the phenotype, proteomic analysis identified an important metabolic enzyme, ATP-citrate lyase, as a novel PTPN2 substrate.

Protection of ß cells against pro-inflammatory cytokine stress by the GDF15-ERBB2 signaling.

Aim/hypothesis: Growth/differentiation factor 15 (GDF15) is a therapeutic target for a variety of metabolic diseases, including type 1 diabetes (T1D). However, the nausea caused by GDF15 is a challenging point for therapeutic development. In addition, it is unknown why the endogenous GDF15 fails to protect from T1D development. Here, we investigate the GDF15 signaling in pancreatic islets towards opening possibilities for therapeutic targeting in ß cells and to understand why this protection fails to occur naturally. Methods: GDF15 signaling in islets was determined by proximity-ligation assay, untargeted proteomics, pathway analysis, and treatment of cells with specific inhibitors. To determine if GDF15 levels would increase prior to disease onset, plasma levels of GDF15 were measured in a longitudinal prospective study of children during T1D development (n=132 cases vs. n=40 controls) and in children with islet autoimmunity but normoglycemia (n=47 cases vs. n=40 controls) using targeted mass spectrometry. We also investigated the regulation of GDF15 production in islets by fluorescence microscopy and western blot analysis. Results: The proximity-ligation assay identified ERBB2 as the GDF15 receptor in islets, which was confirmed using its specific antagonist, tucatinib. The untargeted proteomics analysis and caspase assay showed that ERBB2 activation by GDF15 reduces ß cell apoptosis by downregulating caspase 8. In plasma, GDF15 levels were higher (p=0.0024) during T1D development compared to controls, but not in islet autoimmunity with normoglycemia. However, in the pancreatic islets GDF15 was depleted via sequestration of its mRNA into stress granules, resulting in translation halting. Conclusions/interpretation: GDF15 protects against T1D via ERBB2-mediated decrease of caspase 8 expression in pancreatic islets. Circulating levels of GDF15 increases pre-T1D onset, which is insufficient to promote protection due to its localized depletion in the islets. These findings open opportunities for targeting GDF15 downstream signaling for pancreatic ß cell protection in T1D and help to explain the lack of natural protection by the endogenous protein.

Long-read RNA-seq atlas of novel microglia isoforms elucidates disease-associated genetic regulation of splicing.

Microglia, the innate immune cells of the central nervous system, have been genetically implicated in multiple neurodegenerative diseases. We previously mapped the genetic regulation of gene expression and mRNA splicing in human microglia, identifying several loci where common genetic variants in microglia-specific regulatory elements explain disease risk loci identified by GWAS. However, identifying genetic effects on splicing has been challenging due to the use of short sequencing reads to identify causal isoforms. Here we present the isoform-centric microglia genomic atlas (isoMiGA) which leverages the power of long-read RNA-seq to identify 35,879 novel microglia isoforms. We show that the novel microglia isoforms are involved in stimulation response and brain region specificity. We then quantified the expression of both known and novel isoforms in a multi-ethnic meta-analysis of 555 human microglia short-read RNA-seq samples from 391 donors, the largest to date, and found associations with genetic risk loci in Alzheimer's disease and Parkinson's disease. We nominate several loci that may act through complex changes in isoform and splice site usage.

Analysis of a macrophage carbamylated proteome reveals a function in post-translational modification crosstalk.

BACKGROUND: Lysine carbamylation is a biomarker of rheumatoid arthritis and kidney diseases. However, its cellular function is understudied due to the lack of tools for systematic analysis of this post-translational modification (PTM). METHODS: We adapted a method to analyze carbamylated peptides by co-affinity purification with acetylated peptides based on the cross-reactivity of anti-acetyllysine antibodies. We also performed immobilized-metal affinity chromatography to enrich for phosphopeptides, which allowed us to obtain multi-PTM information from the same samples. RESULTS: By testing the pipeline with RAW 264.7 macrophages treated with bacterial lipopolysaccharide, 7,299, 8,923 and 47,637 acetylated, carbamylated, and phosphorylated peptides were identified, respectively. Our analysis showed that carbamylation occurs on proteins from a variety of functions on sites with similar as well as distinct motifs compared to acetylation. To investigate possible PTM crosstalk, we integrated the carbamylation data with acetylation and phosphorylation data, leading to the identification 1,183 proteins that were modified by all 3 PTMs. Among these proteins, 54 had all 3 PTMs regulated by lipopolysaccharide and were enriched in immune signaling pathways, and in particular, the ubiquitin-proteasome pathway. We found that carbamylation of linear diubiquitin blocks the activity of the anti-inflammatory deubiquitinase OTULIN. CONCLUSIONS: Overall, our data show that anti-acetyllysine antibodies can be used for effective enrichment of carbamylated peptides. Moreover, carbamylation may play a role in PTM crosstalk with acetylation and phosphorylation, and that it is involved in regulating ubiquitination in vitro. Video Abstract.

Myalgic encephalomyelitis/chronic fatigue syndrome and fibromyalgia are indistinguishable by their cerebrospinal fluid proteomes.

BACKGROUND: Myalgic Encephalomyelitis/Chronic Fatigue Syndrome (ME/CFS) and fibromyalgia have overlapping neurologic symptoms particularly disabling fatigue. This has given rise to the question whether they are distinct central nervous system (CNS) entities or is one an extension of the other. MATERIAL AND METHODS: To investigate this, we used unbiased quantitative mass spectrometry-based proteomics to examine the most proximal fluid to the brain, cerebrospinal fluid (CSF). This was to ascertain if the proteome profile of one was the same or different from the other. We examined two separate groups of ME/CFS, one with (n = 15) and one without (n = 15) fibromyalgia. RESULTS: We quantified a total of 2083 proteins using immunoaffinity depletion, tandem mass tag isobaric labelling and offline two-dimensional liquid chromatography coupled to tandem mass spectrometry, including 1789 that were quantified in all the CSF samples. ANOVA analysis did not yield any proteins with an adjusted p value <.05. CONCLUSION: This supports the notion that ME/CFS and fibromyalgia as currently defined are not distinct entities.Key messageME/CFS and fibromyalgia as currently defined are not distinct entities.Unbiased quantitative mass spectrometry-based proteomics can be used to discover cerebrospinal fluid proteins that are biomarkers for a condition such as we are studying.

Robust collection and processing for label-free single voxel proteomics.

With advanced mass spectrometry (MS)-based proteomics, genome-scale proteome coverage can be achieved from bulk tissues. However, such bulk measurement lacks spatial resolution and obscures important tissue heterogeneity, which make it impossible for proteome mapping of tissue microenvironment. Here we report an integrated wet collection of single tissue voxel and Surfactant-assisted One-Pot voxel processing method termed wcSOP for robust label-free single voxel proteomics. wcSOP capitalizes on buffer droplet-assisted wet collection of single tissue voxel dissected by LCM into the PCR tube cap and MS-compatible surfactant-assisted one-pot voxel processing in the collection cap. This convenient method allows reproducible label-free quantification of ∼900 and ∼4,600 proteins for single voxel from fresh frozen human spleen tissue at 20 µm × 20 µm × 10 µm (close to single cells) and 200 µm × 200 µm × 10 µm (∼100 cells), respectively. 100s-1000s of protein signatures with differential expression levels were identified to be spatially resolved between spleen red and white pulp regions depending on the voxel size. Region-specific signaling pathways were enriched from single voxel proteomics data. Antibody-based CODEX imaging was used to validate label-free MS quantitation for single voxel analysis. The wcSOP-MS method paves the way for routine robust single voxel proteomics and spatial proteomics.

Proteogenomic insights suggest druggable pathways in endometrial carcinoma.

We characterized a prospective endometrial carcinoma (EC) cohort containing 138 tumors and 20 enriched normal tissues using 10 different omics platforms. Targeted quantitation of two peptides can predict antigen processing and presentation machinery activity, and may inform patient selection for immunotherapy. Association analysis between MYC activity and metformin treatment in both patients and cell lines suggests a potential role for metformin treatment in non-diabetic patients with elevated MYC activity. PIK3R1 in-frame indels are associated with elevated AKT phosphorylation and increased sensitivity to AKT inhibitors. CTNNB1 hotspot mutations are concentrated near phosphorylation sites mediating pS45-induced degradation of ß-catenin, which may render Wnt-FZD antagonists ineffective. Deep learning accurately predicts EC subtypes and mutations from histopathology images, which may be useful for rapid diagnosis. Overall, this study identified molecular and imaging markers that can be further investigated to guide patient stratification for more precise treatment of EC.

Suspension Trapping-Based Sample Preparation Workflow for In-Depth Plant Phosphoproteomics.

Plant phosphoproteomics provides a global view of phosphorylation-mediated signaling in plants; however, it demands high-throughput methods with sensitive detection and accurate quantification. Despite the widespread use of protein precipitation for removing contaminants and improving sample purity, it limits the sensitivity and throughput of plant phosphoproteomic analysis. The multiple handling steps involved in protein precipitation lead to sample loss and process variability. Herein, we developed an approach based on suspension trapping (S-Trap), termed tandem S-Trap-IMAC (immobilized metal ion affinity chromatography), by integrating an S-Trap micro-column with a Fe-IMAC tip. Compared with a precipitation-based workflow, the tandem S-Trap-IMAC method deepened the coverage of the Arabidopsis (Arabidopsis thaliana) phosphoproteome by more than 30%, with improved number of multiply phosphorylated peptides, quantification accuracy, and short sample processing time. We applied the tandem S-Trap-IMAC method for studying abscisic acid (ABA) signaling in Arabidopsis seedlings. We thus discovered that a significant proportion of the phosphopeptides induced by ABA are multiply phosphorylated peptides, indicating their importance in early ABA signaling and quantified several key phosphorylation sites on core ABA signaling components across four time points. Our results show that the optimized workflow aids high-throughput phosphoproteome profiling of low-input plant samples.

Analysis of a macrophage carbamylated proteome reveals a function in post-translational modification crosstalk.

Background. Lysine carbamylation is a biomarker of rheumatoid arthritis and kidney diseases. However, its cellular function is understudied due to the lack of tools for systematic analysis of this post-translational modification (PTM). Methods. We adapted a method to analyze carbamylated peptides by co-affinity purification with acetylated peptides based on the cross-reactivity of anti-acetyllysine antibodies. We integrated this method into a mass spectrometry-based multi-PTM pipeline to simultaneously analyze carbamylated and acetylated peptides in addition to phosphopeptides were enriched by sequential immobilized-metal affinity chromatography. Results. By testing the pipeline with RAW 264.7 macrophages treated with bacterial lipopolysaccharide, 7,299, 8,923 and 47,637 acetylated, carbamylated, and phosphorylated peptides were identified, respectively. Our analysis showed that carbamylation occurs on proteins from a variety of functions on sites with similar as well as distinct motifs compared to acetylation. To investigate possible PTM crosstalk, we integrated the carbamylation data with acetylation and phosphorylation data, leading to the identification 1,183 proteins that were modified by all 3 PTMs. Among these proteins, 54 had all 3 PTMs regulated by lipopolysaccharide and were enriched in immune signaling pathways, and in particular, the ubiquitin-proteasome pathway. We found that carbamylation of linear diubiquitin blocks the activity of the anti-inflammatory deubiquitinase OTULIN. Conclusions Overall, our data show that anti-acetyllysine antibodies can be used for effective enrichment of carbamylated peptides. Moreover, carbamylation may play a role in PTM crosstalk with acetylation and phosphorylation, and that it is involved in regulating ubiquitination in vitro .

exRNA-eCLIP intersection analysis reveals a map of extracellular RNA binding proteins and associated RNAs across major human biofluids and carriers.

Although the role of RNA binding proteins (RBPs) in extracellular RNA (exRNA) biology is well established, their exRNA cargo and distribution across biofluids are largely unknown. To address this gap, we extend the exRNA Atlas resource by mapping exRNAs carried by extracellular RBPs (exRBPs). This map was developed through an integrative analysis of ENCODE enhanced crosslinking and immunoprecipitation (eCLIP) data (150 RBPs) and human exRNA profiles (6,930 samples). Computational analysis and experimental validation identified exRBPs in plasma, serum, saliva, urine, cerebrospinal fluid, and cell-culture-conditioned medium. exRBPs carry exRNA transcripts from small non-coding RNA biotypes, including microRNA (miRNA), piRNA, tRNA, small nuclear RNA (snRNA), small nucleolar RNA (snoRNA), Y RNA, and lncRNA, as well as protein-coding mRNA fragments. Computational deconvolution of exRBP RNA cargo reveals associations of exRBPs with extracellular vesicles, lipoproteins, and ribonucleoproteins across human biofluids. Overall, we mapped the distribution of exRBPs across human biofluids, presenting a resource for the community.

Computational ranking-assisted identification of Plexin-B2 in homotypic and heterotypic clustering of circulating tumor cells in breast cancer metastasis.

Metastasis is the cause of over 90% of all deaths associated with breast cancer, yet the strategies to predict cancer spreading based on primary tumor profiles and therefore prevent metastasis are egregiously limited. As rare precursor cells to metastasis, circulating tumor cells (CTCs) in multicellular clusters in the blood are 20-50 times more likely to produce viable metastasis than single CTCs. However, the molecular mechanisms underlying various CTC clusters, such as homotypic tumor cell clusters and heterotypic tumor-immune cell clusters, are yet to be fully elucidated. Combining machine learning-assisted computational ranking with experimental demonstration to assess cell adhesion candidates, we identified a transmembrane protein Plexin- B2 (PB2) as a new therapeutic target that drives the formation of both homotypic and heterotypic CTC clusters. High PB2 expression in human primary tumors predicts an unfavorable distant metastasis-free survival and is enriched in CTC clusters compared to single CTCs in advanced breast cancers. Loss of PB2 reduces formation of homotypic tumor cell clusters as well as heterotypic tumor-myeloid cell clusters in triple-negative breast cancer. Interactions between PB2 and its ligand Sema4C on tumor cells promote homotypic cluster formation, and PB2 binding with Sema4A on myeloid cells (monocytes) drives heterotypic CTC cluster formation, suggesting that metastasizing tumor cells hijack the PB2/Sema family axis to promote lung metastasis in breast cancer. Additionally, using a global proteomic analysis, we identified novel downstream effectors of the PB2 pathway associated with cancer stemness, cell cycling, and tumor cell clustering in breast cancer. Thus, PB2 is a novel therapeutic target for preventing new metastasis.

Phosphorylation barcodes direct biased chemokine signaling at CXCR3.

G protein-coupled receptor (GPCR)-biased agonism, selective activation of certain signaling pathways relative to others, is thought to be directed by differential GPCR phosphorylation "barcodes." At chemokine receptors, endogenous chemokines can act as "biased agonists", which may contribute to the limited success when pharmacologically targeting these receptors. Here, mass spectrometry-based global phosphoproteomics revealed that CXCR3 chemokines generate different phosphorylation barcodes associated with differential transducer activation. Chemokine stimulation resulted in distinct changes throughout the kinome in global phosphoproteomics studies. Mutation of CXCR3 phosphosites altered ß-arrestin 2 conformation in cellular assays and was consistent with conformational changes observed in molecular dynamics simulations. T cells expressing phosphorylation-deficient CXCR3 mutants resulted in agonist- and receptor-specific chemotactic profiles. Our results demonstrate that CXCR3 chemokines are non-redundant and act as biased agonists through differential encoding of phosphorylation barcodes, leading to distinct physiological processes.

Phosphorylation barcodes direct biased chemokine signaling at CXCR3.

G protein-coupled receptor (GPCR) biased agonism, the activation of some signaling pathways over others, is thought to largely be due to differential receptor phosphorylation, or "phosphorylation barcodes." At chemokine receptors, ligands act as "biased agonists" with complex signaling profiles, which contributes to the limited success in pharmacologically targeting these receptors. Here, mass spectrometry-based global phosphoproteomics revealed that CXCR3 chemokines generate different phosphorylation barcodes associated with differential transducer activation. Chemokine stimulation resulted in distinct changes throughout the kinome in global phosphoproteomic studies. Mutation of CXCR3 phosphosites altered ß-arrestin conformation in cellular assays and was confirmed by molecular dynamics simulations. T cells expressing phosphorylation-deficient CXCR3 mutants resulted in agonist- and receptor-specific chemotactic profiles. Our results demonstrate that CXCR3 chemokines are non-redundant and act as biased agonists through differential encoding of phosphorylation barcodes and lead to distinct physiological processes.

Single-Cell Discovery and Multiomic Characterization of Therapeutic Targets in Multiple Myeloma.

Multiple myeloma (MM) is a highly refractory hematologic cancer. Targeted immunotherapy has shown promise in MM but remains hindered by the challenge of identifying specific yet broadly representative tumor markers. We analyzed 53 bone marrow (BM) aspirates from 41 MM patients using an unbiased, high-throughput pipeline for therapeutic target discovery via single-cell transcriptomic profiling, yielding 38 MM marker genes encoding cell-surface proteins and 15 encoding intracellular proteins. Of these, 20 candidate genes were highlighted that are not yet under clinical study, 11 of which were previously uncharacterized as therapeutic targets. The findings were cross-validated using bulk RNA sequencing, flow cytometry, and proteomic mass spectrometry of MM cell lines and patient BM, demonstrating high overall concordance across data types. Independent discovery using bulk RNA sequencing reiterated top candidates, further affirming the ability of single-cell transcriptomics to accurately capture marker expression despite limitations in sample size or sequencing depth. Target dynamics and heterogeneity were further examined using both transcriptomic and immuno-imaging methods. In summary, this study presents a robust and broadly applicable strategy for identifying tumor markers to better inform the development of targeted cancer therapy. SIGNIFICANCE: Single-cell transcriptomic profiling and multiomic cross-validation to uncover therapeutic targets identifies 38 myeloma marker genes, including 11 transcribing surface proteins with previously uncharacterized potential for targeted antitumor therapy.

A streamlined tandem tip-based workflow for sensitive nanoscale phosphoproteomics.

Effective phosphoproteome of nanoscale sample analysis remains a daunting task, primarily due to significant sample loss associated with non-specific surface adsorption during enrichment of low stoichiometric phosphopeptide. We develop a tandem tip phosphoproteomics sample preparation method that is capable of sample cleanup and enrichment without additional sample transfer, and its integration with our recently developed SOP (Surfactant-assisted One-Pot sample preparation) and iBASIL (improved Boosting to Amplify Signal with Isobaric Labeling) approaches provides a streamlined workflow enabling sensitive, high-throughput nanoscale phosphoproteome measurements. This approach significantly reduces both sample loss and processing time, allowing the identification of >3000 (>9500) phosphopeptides from 1 (10) µg of cell lysate using the label-free method without a spectral library. It also enables precise quantification of ~600 phosphopeptides from 100 sorted cells (single-cell level input for the enriched phosphopeptides) and ~700 phosphopeptides from human spleen tissue voxels with a spatial resolution of 200 µm (equivalent to ~100 cells) in a high-throughput manner. The new workflow opens avenues for phosphoproteome profiling of mass-limited samples at the low nanogram level.

Machine learning-assisted elucidation of CD81-CD44 interactions in promoting cancer stemness and extracellular vesicle integrity.

Tumor-initiating cells with reprogramming plasticity or stem-progenitor cell properties (stemness) are thought to be essential for cancer development and metastatic regeneration in many cancers; however, elucidation of the underlying molecular network and pathways remains demanding. Combining machine learning and experimental investigation, here we report CD81, a tetraspanin transmembrane protein known to be enriched in extracellular vesicles (EVs), as a newly identified driver of breast cancer stemness and metastasis. Using protein structure modeling and interface prediction-guided mutagenesis, we demonstrate that membrane CD81 interacts with CD44 through their extracellular regions in promoting tumor cell cluster formation and lung metastasis of triple negative breast cancer (TNBC) in human and mouse models. In-depth global and phosphoproteomic analyses of tumor cells deficient with CD81 or CD44 unveils endocytosis-related pathway alterations, leading to further identification of a quality-keeping role of CD44 and CD81 in EV secretion as well as in EV-associated stemness-promoting function. CD81 is coexpressed along with CD44 in human circulating tumor cells (CTCs) and enriched in clustered CTCs that promote cancer stemness and metastasis, supporting the clinical significance of CD81 in association with patient outcomes. Our study highlights machine learning as a powerful tool in facilitating the molecular understanding of new molecular targets in regulating stemness and metastasis of TNBC.

Ablating Lgr5-expressing prostatic stromal cells activates the ERK-mediated mechanosensory signaling and disrupts prostate tissue homeostasis.

Functional implication of stromal heterogeneity in the prostate remains incompletely understood. Using lineage tracing and light-sheet imaging, we show that some fibroblast cells at the mouse proximal prostatic ducts and prostatic urethra highly express Lgr5. Genetic ablation of these anatomically restricted stromal cells, but not nonselective ablation of prostatic stromal cells, rapidly induces prostate epithelial turnover and dedifferentiation that are reversed following spontaneous restoration of the Lgr5+ stromal cells. RNA sequencing (RNA-seq) analysis indicates that ablating the Lgr5+ stromal cells activates a mechanosensory response. Ablating the Lgr5+ stromal cells impairs the control of prostatic ductal outlet, increases prostate tissue stiffness, and activates the mitogen-activated protein kinase (MAPK). Suppressing MAPK overrides the elevated epithelial proliferation. In summary, the Lgr5+ stromal cells regulate prostate tissue homeostasis and maintain its functional integrity in a long-distance manner. Our study implies that the cells at organ junctions most likely control organ homeostasis by sustaining a balanced mechanoforce.