Cross-domain information fusion for enhanced cell population delineation in single-cell spatial-omics data.

Cell population delineation and identification is an essential step in single-cell and spatial-omics studies. Spatial-omics technologies can simultaneously measure information from three complementary domains related to this task: expression levels of a panel of molecular biomarkers at single-cell resolution, relative positions of cells, and images of tissue sections, but existing computational methods for performing this task on single-cell spatial-omics datasets often relinquish information from one or more domains. The additional reliance on the availability of "atlas" training or reference datasets limits cell type discovery to well-defined but limited cell population labels, thus posing major challenges for using these methods in practice. Successful integration of all three domains presents an opportunity for uncovering cell populations that are functionally stratified by their spatial contexts at cellular and tissue levels: the key motivation for employing spatial-omics technologies in the first place. In this work, we introduce Cell S patio- and N eighborhood-informed A nnotation and P atterning (CellSNAP), a self-supervised computational method that learns a representation vector for each cell in tissue samples measured by spatial-omics technologies at the single-cell or finer resolution. The learned representation vector fuses information about the corresponding cell across all three aforementioned domains. By applying CellSNAP to datasets spanning both spatial proteomic and spatial transcriptomic modalities, and across different tissue types and disease settings, we show that CellSNAP markedly enhances de novo discovery of biologically relevant cell populations at fine granularity, beyond current approaches, by fully integrating cells' molecular profiles with cellular neighborhood and tissue image information.

Methionine and one carbon metabolism as a regulator of bone remodeling with fasting.

Caloric restriction improves metabolic health, but is often complicated by bone loss. We studied bone parameters in humans during a 10-day fast and identified candidate metabolic regulators of bone turnover. P1NP, a bone formation marker, decreased within 3 days of fasting. Whereas dual-energy X-ray absorptiometry measures of bone mineral density were unchanged after 10 days of fasting, high-resolution peripheral quantitative CT demonstrated remodeling of bone microarchitecture. Pathway analysis of longitudinal metabolomics data identified one-carbon metabolism as fasting-dependent. In cultured osteoblasts, we tested the functional significance of one-carbon metabolites modulated by fasting, finding that methionine - which surged after 3 days of fasting - impacted markers of osteoblast cell state in a concentration dependent manner, in some instances exhibiting a U-shaped response with both low and high concentrations driving putative anti-bone responses. Administration of methionine to mice for 5 days recapitulated some fasting effects on bone, including a reduction in serum P1NP. In conclusion, a 10-day fast in humans led to remodeling of bone microarchitecture, potentially mediated by a surge in circulating methionine. These data support an emerging model that points to a window of optimal methionine exposure for bone health.

Epigenetically active chromatin in neonatal iWAT reveals GABPα as a potential regulator of beige adipogenesis.

Background: Thermogenic beige adipocytes, which dissipate energy as heat, are found in neonates and adults. Recent studies show that neonatal beige adipocytes are highly plastic and contribute to >50% of beige adipocytes in adults. Neonatal beige adipocytes are distinct from recruited beige adipocytes in that they develop independently of temperature and sympathetic innervation through poorly defined mechanisms. Methods: We characterized the neonatal beige adipocytes in the inguinal white adipose tissue (iWAT) of C57BL6 postnatal day 3 and 20 mice (P3 and P20) by imaging, genome-wide RNA-seq analysis, ChIP-seq analysis, qRT-PCR validation, and biochemical assays. Results: We found an increase in acetylated histone 3 lysine 27 (H3K27ac) on the promoter and enhancer regions of beige-specific gene UCP1 in iWAT of P20 mice. Furthermore, H3K27ac ChIP-seq analysis in the iWAT of P3 and P20 mice revealed strong H3K27ac signals at beige adipocyte-associated genes in the iWAT of P20 mice. The integration of H3K27ac ChIP-seq and RNA-seq analysis in the iWAT of P20 mice reveal epigenetically active signatures of beige adipocytes, including oxidative phosphorylation and mitochondrial metabolism. We identify the enrichment of GA-binding protein alpha (GABPα) binding regions in the epigenetically active chromatin regions of the P20 iWAT, particularly on beige genes, and demonstrate that GABPα is required for beige adipocyte differentiation. Moreover, transcriptomic analysis and glucose oxidation assays revealed increased glycolytic activity in the neonatal iWAT from P20. Conclusions: Our findings demonstrate that epigenetic mechanisms regulate the development of peri-weaning beige adipocytes via GABPα. Further studies to better understand the upstream mechanisms that regulate epigenetic activation of GABPα and characterization of the metabolic identity of neonatal beige adipocytes will help us harness their therapeutic potential in metabolic diseases.

A Hitchhiker's guide to high-dimensional tissue imaging with multiplexed ion beam imaging.

Advancements in multiplexed tissue imaging technologies are vital in shaping our understanding of tissue microenvironmental influences in disease contexts. These technologies now allow us to relate the phenotype of individual cells to their higher-order roles in tissue organization and function. Multiplexed Ion Beam Imaging (MIBI) is one of such technologies, which uses metal isotope-labeled antibodies and secondary ion mass spectrometry (SIMS) to image more than 40 protein markers simultaneously within a single tissue section. Here, we describe an optimized MIBI workflow for high-plex analysis of Formalin-Fixed Paraffin-Embedded (FFPE) tissues following antigen retrieval, metal isotope-conjugated antibody staining, imaging using the MIBI instrument, and subsequent data processing and analysis. While this workflow is focused on imaging human FFPE samples using the MIBI, this workflow can be easily extended to model systems, biological questions, and multiplexed imaging modalities.

Nuclear speckle rejuvenation alleviates proteinopathies at the expense of YAP1.

Current treatments targeting individual protein quality control have limited efficacy in alleviating proteinopathies, highlighting the prerequisite for a common upstream druggable target capable of global proteostasis modulation. Building on our prior research establishing nuclear speckles as pivotal organelles responsible for global proteostasis transcriptional control, we aim to alleviate proteinopathies through nuclear speckle rejuvenation. We identified pyrvinium pamoate as a small-molecule nuclear speckle rejuvenator that enhances protein quality control while suppressing YAP1 signaling via decreasing the surface tension of nuclear speckle condensates through interaction with the intrinsically disordered region of nuclear speckle scaffold protein SON. In pre-clinical models, pyrvinium pamoate reduced tauopathy and alleviated retina degeneration by promoting autophagy and ubiquitin-proteasome system. Aberrant nuclear speckle morphology, reduced protein quality control and increased YAP1 activity were also observed in human tauopathies. Our study uncovers novel therapeutic targets for tackling protein misfolding disorders within an expanded proteostasis framework encompassing nuclear speckles and YAP1.

A Spatial Multi-Modal Dissection of Host-Microbiome Interactions within the Colitis Tissue Microenvironment.

The intricate and dynamic interactions between the host immune system and its microbiome constituents undergo dynamic shifts in response to perturbations to the intestinal tissue environment. Our ability to study these events on the systems level is significantly limited by in situ approaches capable of generating simultaneous insights from both host and microbial communities. Here, we introduce Microbiome Cartography (MicroCart), a framework for simultaneous in situ probing of host features and its microbiome across multiple spatial modalities. We demonstrate MicroCart by comprehensively investigating the alterations in both gut host and microbiome components in a murine model of colitis by coupling MicroCart with spatial proteomics, transcriptomics, and glycomics platforms. Our findings reveal a global but systematic transformation in tissue immune responses, encompassing tissue-level remodeling in response to host immune and epithelial cell state perturbations, and bacterial population shifts, localized inflammatory responses, and metabolic process alterations during colitis. MicroCart enables a deep investigation of the intricate interplay between the host tissue and its microbiome with spatial multiomics.

Epstein-Barr Virus Orchestrates Spatial Reorganization and Immunomodulation within the Classic Hodgkin Lymphoma Tumor Microenvironment.

Classic Hodgkin Lymphoma (cHL) is a tumor composed of rare malignant Hodgkin and Reed-Sternberg (HRS) cells nested within a T-cell rich inflammatory immune infiltrate. cHL is associated with Epstein-Barr Virus (EBV) in 25% of cases. The specific contributions of EBV to the pathogenesis of cHL remain largely unknown, in part due to technical barriers in dissecting the tumor microenvironment (TME) in high detail. Herein, we applied multiplexed ion beam imaging (MIBI) spatial pro-teomics on 6 EBV-positive and 14 EBV-negative cHL samples. We identify key TME features that distinguish between EBV-positive and EBV-negative cHL, including the relative predominance of memory CD8 T cells and increased T-cell dysfunction as a function of spatial proximity to HRS cells. Building upon a larger multi-institutional cohort of 22 EBV-positive and 24 EBV-negative cHL samples, we orthogonally validated our findings through a spatial multi-omics approach, coupling whole transcriptome capture with antibody-defined cell types for tu-mor and T-cell populations within the cHL TME. We delineate contrasting transcriptomic immunological signatures between EBV-positive and EBV-negative cases that differently impact HRS cell proliferation, tumor-immune interactions, and mecha-nisms of T-cell dysregulation and dysfunction. Our multi-modal framework enabled a comprehensive dissection of EBV-linked reorganization and immune evasion within the cHL TME, and highlighted the need to elucidate the cellular and molecular fac-tors of virus-associated tumors, with potential for targeted therapeutic strategies.

MAPS: pathologist-level cell type annotation from tissue images through machine learning.

Highly multiplexed protein imaging is emerging as a potent technique for analyzing protein distribution within cells and tissues in their native context. However, existing cell annotation methods utilizing high-plex spatial proteomics data are resource intensive and necessitate iterative expert input, thereby constraining their scalability and practicality for extensive datasets. We introduce MAPS (Machine learning for Analysis of Proteomics in Spatial biology), a machine learning approach facilitating rapid and precise cell type identification with human-level accuracy from spatial proteomics data. Validated on multiple in-house and publicly available MIBI and CODEX datasets, MAPS outperforms current annotation techniques in terms of speed and accuracy, achieving pathologist-level precision even for typically challenging cell types, including tumor cells of immune origin. By democratizing rapidly deployable and scalable machine learning annotation, MAPS holds significant potential to expedite advances in tissue biology and disease comprehension.

A spatial cell atlas of neuroblastoma reveals developmental, epigenetic and spatial axis of tumor heterogeneity.

Neuroblastoma is a pediatric cancer arising from the developing sympathoadrenal lineage with complex inter- and intra-tumoral heterogeneity. To chart this complexity, we generated a comprehensive cell atlas of 55 neuroblastoma patient tumors, collected from two pediatric cancer institutions, spanning a range of clinical, genetic, and histologic features. Our atlas combines single-cell/nucleus RNA-seq (sc/scRNA-seq), bulk RNA-seq, whole exome sequencing, DNA methylation profiling, spatial transcriptomics, and two spatial proteomic methods. Sc/snRNA-seq revealed three malignant cell states with features of sympathoadrenal lineage development. All of the neuroblastomas had malignant cells that resembled sympathoblasts and the more differentiated adrenergic cells. A subset of tumors had malignant cells in a mesenchymal cell state with molecular features of Schwann cell precursors. DNA methylation profiles defined four groupings of patients, which differ in the degree of malignant cell heterogeneity and clinical outcomes. Using spatial proteomics, we found that neuroblastomas are spatially compartmentalized, with malignant tumor cells sequestered away from immune cells. Finally, we identify spatially restricted signaling patterns in immune cells from spatial transcriptomics. To facilitate the visualization and analysis of our atlas as a resource for further research in neuroblastoma, single cell, and spatial-omics, all data are shared through the Human Tumor Atlas Network Data Commons at www.humantumoratlas.org.

T cell-mediated curation and restructuring of tumor tissue coordinates an effective immune response.

Antigen-specific T cells traffic to, are influenced by, and create unique cellular microenvironments. Here we characterize these microenvironments over time with multiplexed imaging in a melanoma model of adoptive T cell therapy and human patients with melanoma treated with checkpoint inhibitor therapy. Multicellular neighborhood analysis reveals dynamic immune cell infiltration and inflamed tumor cell neighborhoods associated with CD8+ T cells. T cell-focused analysis indicates T cells are found along a continuum of neighborhoods that reflect the progressive steps coordinating the anti-tumor immune response. More effective anti-tumor immune responses are characterized by inflamed tumor-T cell neighborhoods, flanked by dense immune infiltration neighborhoods. Conversely, ineffective T cell therapies express anti-inflammatory cytokines, resulting in regulatory neighborhoods, spatially disrupting productive T cell-immune and -tumor interactions. Our study provides in situ mechanistic insights into temporal tumor microenvironment changes, cell interactions critical for response, and spatial correlates of immunotherapy outcomes, informing cellular therapy evaluation and engineering.

LEF1 isoforms regulate cellular senescence and aging.

The study of aging and its mechanisms, such as cellular senescence, has provided valuable insights into age-related pathologies, thus contributing to their prevention and treatment. The current abundance of high-throughput data combined with the surge of robust analysis algorithms has facilitated novel ways of identifying underlying pathways that may drive these pathologies. For the purpose of identifying key regulators of lung aging, we performed comparative analyses of transcriptional profiles of aged versus young human subjects and mice, focusing on the common age-related changes in the transcriptional regulation in lung macrophages, T cells, and B immune cells. Importantly, we validated our findings in cell culture assays and human lung samples. Our analysis identified lymphoid enhancer binding factor 1 (LEF1) as an important age-associated regulator of gene expression in all three cell types across different tissues and species. Follow-up experiments showed that the differential expression of long and short LEF1 isoforms is a key regulatory mechanism of cellular senescence. Further examination of lung tissue from patients with idiopathic pulmonary fibrosis, an age-related disease with strong ties to cellular senescence, revealed a stark dysregulation of LEF1. Collectively, our results suggest that LEF1 is a key factor of aging, and its differential regulation is associated with human and murine cellular senescence.

Integration of spatial and single-cell data across modalities with weakly linked features.

Although single-cell and spatial sequencing methods enable simultaneous measurement of more than one biological modality, no technology can capture all modalities within the same cell. For current data integration methods, the feasibility of cross-modal integration relies on the existence of highly correlated, a priori 'linked' features. We describe matching X-modality via fuzzy smoothed embedding (MaxFuse), a cross-modal data integration method that, through iterative coembedding, data smoothing and cell matching, uses all information in each modality to obtain high-quality integration even when features are weakly linked. MaxFuse is modality-agnostic and demonstrates high robustness and accuracy in the weak linkage scenario, achieving 20~70% relative improvement over existing methods under key evaluation metrics on benchmarking datasets. A prototypical example of weak linkage is the integration of spatial proteomic data with single-cell sequencing data. On two example analyses of this type, MaxFuse enabled the spatial consolidation of proteomic, transcriptomic and epigenomic information at single-cell resolution on the same tissue section.

LEF1 isoforms regulate cellular senescence and aging.

Background: The study of aging and its mechanisms, such as cellular senescence, has provided valuable insights into age-related pathologies, thus contributing to their prevention and treatment. The current abundance of high throughput data combined with the surge of robust analysis algorithms has facilitated novel ways of identifying underlying pathways that may drive these pathologies. Methods: With the focus on identifying key regulators of lung aging, we performed comparative analyses of transcriptional profiles of aged versus young human subjects and mice, focusing on the common age-related changes in the transcriptional regulation in lung macrophages, T cells, and B immune cells. Importantly, we validated our findings in cell culture assays and human lung samples. Results: We identified Lymphoid Enhancer Binding Factor 1 (LEF1) as an important age-associated regulator of gene expression in all three cell types across different tissues and species. Follow-up experiments showed that the differential expression of long and short LEF1 isoforms is a key regulatory mechanism of cellular senescence. Further examination of lung tissue from patients with Idiopathic Pulmonary Fibrosis (IPF), an age-related disease with strong ties to cellular senescence, we demonstrated a stark dysregulation of LEF1. Conclusions: Collectively, our results suggest that the LEF1 is a key factor of aging, and its differential regulation is associated with human and murine cellular senescence.

MAPS: Pathologist-level cell type annotation from tissue images through machine learning.

Highly multiplexed protein imaging is emerging as a potent technique for analyzing protein distribution within cells and tissues in their native context. However, existing cell annotation methods utilizing high-plex spatial proteomics data are resource intensive and necessitate iterative expert input, thereby constraining their scalability and practicality for extensive datasets. We introduce MAPS (Machine learning for Analysis of Proteomics in Spatial biology), a machine learning approach facilitating rapid and precise cell type identification with human-level accuracy from spatial proteomics data. Validated on multiple in-house and publicly available MIBI and CODEX datasets, MAPS outperforms current annotation techniques in terms of speed and accuracy, achieving pathologist-level precision even for challenging cell types, including tumor cells of immune origin. By democratizing rapidly deployable and scalable machine learning annotation, MAPS holds significant potential to expedite advances in tissue biology and disease comprehension.

Organization of the human intestine at single-cell resolution.

The intestine is a complex organ that promotes digestion, extracts nutrients, participates in immune surveillance, maintains critical symbiotic relationships with microbiota and affects overall health1. The intesting has a length of over nine metres, along which there are differences in structure and function2. The localization of individual cell types, cell type development trajectories and detailed cell transcriptional programs probably drive these differences in function. Here, to better understand these differences, we evaluated the organization of single cells using multiplexed imaging and single-nucleus RNA and open chromatin assays across eight different intestinal sites from nine donors. Through systematic analyses, we find cell compositions that differ substantially across regions of the intestine and demonstrate the complexity of epithelial subtypes, and find that the same cell types are organized into distinct neighbourhoods and communities, highlighting distinct immunological niches that are present in the intestine. We also map gene regulatory differences in these cells that are suggestive of a regulatory differentiation cascade, and associate intestinal disease heritability with specific cell types. These results describe the complexity of the cell composition, regulation and organization for this organ, and serve as an important reference map for understanding human biology and disease.

Expanded vacuum-stable gels for multiplexed high-resolution spatial histopathology.

Cellular organization and functions encompass multiple scales in vivo. Emerging high-plex imaging technologies are limited in resolving subcellular biomolecular features. Expansion Microscopy (ExM) and related techniques physically expand samples for enhanced spatial resolution, but are challenging to be combined with high-plex imaging technologies to enable integrative multiscaled tissue biology insights. Here, we introduce Expand and comPRESS hydrOgels (ExPRESSO), an ExM framework that allows high-plex protein staining, physical expansion, and removal of water, while retaining the lateral tissue expansion. We demonstrate ExPRESSO imaging of archival clinical tissue samples on Multiplexed Ion Beam Imaging and Imaging Mass Cytometry platforms, with detection capabilities of > 40 markers. Application of ExPRESSO on archival human lymphoid and brain tissues resolved tissue architecture at the subcellular level, particularly that of the blood-brain barrier. ExPRESSO hence provides a platform for extending the analysis compatibility of hydrogel-expanded biospecimens to mass spectrometry, with minimal modifications to protocols and instrumentation.

Preservation of ∼12-h ultradian rhythms of gene expression of mRNA and protein metabolism in the absence of canonical circadian clock.

Introduction: Besides the ∼24-h circadian rhythms, ∼12-h ultradian rhythms of gene expression, metabolism and behaviors exist in animals ranging from crustaceans to mammals. Three major hypotheses were proposed on the origin and mechanisms of regulation of ∼12-h rhythms, namely, that they are not cell-autonomous and controlled by a combination of the circadian clock and environmental cues, that they are regulated by two anti-phase circadian transcription factors in a cell autonomous manner, or that they are established by a cell-autonomous ∼12-h oscillator. Methods: To distinguish among these possibilities, we performed a post hoc analysis of two high temporal resolution transcriptome dataset in animals and cells lacking the canonical circadian clock. Results: In both the liver of BMAL1 knockout mice and Drosophila S2 cells, we observed robust and prevalent ∼12-h rhythms of gene expression enriched in fundamental processes of mRNA and protein metabolism that show large convergence with those identified in wild-type mice liver. Bioinformatics analysis further predicted ELF1 and ATF6B as putative transcription factors regulating the ∼12-h rhythms of gene expression independently of the circadian clock in both fly and mice. Discussion: These findings provide additional evidence to support the existence of an evolutionarily conserved 12-h oscillator that controls ∼12-h rhythms of gene expression of protein and mRNA metabolism in multiple species.

TRAF4-mediated nonproteolytic ubiquitination of androgen receptor promotes castration-resistant prostate cancer.

Castration-resistant prostate cancer (CRPC) poses a major clinical challenge with the androgen receptor (AR) remaining to be a critical oncogenic player. Several lines of evidence indicate that AR induces a distinct transcriptional program after androgen deprivation in CRPCs. However, the mechanism triggering AR binding to a distinct set of genomic loci in CRPC and how it promotes CRPC development remain unclear. We demonstrate here that atypical ubiquitination of AR mediated by an E3 ubiquitin ligase TRAF4 plays an important role in this process. TRAF4 is highly expressed in CRPCs and promotes CRPC development. It mediates K27-linked ubiquitination at the C-terminal tail of AR and increases its association with the pioneer factor FOXA1. Consequently, AR binds to a distinct set of genomic loci enriched with FOXA1- and HOXB13-binding motifs to drive different transcriptional programs including an olfactory transduction pathway. Through the surprising upregulation of olfactory receptor gene transcription, TRAF4 increases intracellular cAMP levels and boosts E2F transcription factor activity to promote cell proliferation under androgen deprivation conditions. Altogether, these findings reveal a posttranslational mechanism driving AR-regulated transcriptional reprogramming to provide survival advantages for prostate cancer cells under castration conditions.

Preservation of ∼12-hour ultradian rhythms of gene expression of mRNA and protein metabolism in the absence of canonical circadian clock.

Besides the ∼24-hour circadian rhythms, ∼12-hour ultradian rhythms of gene expression, metabolism and behaviors exist in animals ranging from crustaceans to mammals. Three major hypotheses were proposed on the origin and mechanisms of regulation of ∼12-hour rhythms, namely that they are not cell-autonomous and controlled by a combination of the circadian clock and environmental cues, that they are regulated by two anti-phase circadian transcriptional factors in a cell-autonomous manner, or that they are established by a cell-autonomous ∼12-hour oscillator. To distinguish among these possibilities, we performed a post-hoc analysis of two high temporal resolution transcriptome dataset in animals and cells lacking the canonical circadian clock. In both the liver of BMAL1 knockout mice and Drosophila S2 cells, we observed robust and prevalent ∼12-hour rhythms of gene expression enriched in fundamental processes of mRNA and protein metabolism that show large convergence with those identified in wild-type mice liver. Bioinformatics analysis further predicted ELF1 and ATF6B as putative transcription factors regulating the ∼12-hour rhythms of gene expression independently of the circadian clock in both fly and mice. These findings provide additional evidence to support the existence of an evolutionarily conserved 12-hour oscillator that controls ∼12-hour rhythms of gene expression of protein and mRNA metabolism in multiple species.

Evidence for conservation of primordial ~12-hour ultradian gene programs in humans under free-living conditions.

While circadian rhythms are entrained to the once daily light-dark cycle of the sun, many marine organisms exhibit ~12h ultradian rhythms corresponding to the twice daily movement of the tides. Although human ancestors emerged from circatidal environment millions of years ago, direct evidence of ~12h ultradian rhythms in humans is lacking. Here, we performed prospective, temporal transcriptome profiling of peripheral white blood cells and identified robust ~12h transcriptional rhythms from three healthy participants. Pathway analysis implicated ~12h rhythms in RNA and protein metabolism, with strong homology to the circatidal gene programs previously identified in Cnidarian marine species. We further observed ~12h rhythms of intron retention events of genes involved in MHC class I antigen presentation, synchronized to expression of mRNA splicing genes in all three participants. Gene regulatory network inference revealed XBP1, and GABP and KLF transcription factor family members as potential transcriptional regulators of human ~12h rhythms. These results suggest that human ~12h biological rhythms have a primordial evolutionary origin with important implications for human health and disease.