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.

The molecular cytoarchitecture of the adult mouse brain.

The function of the mammalian brain relies upon the specification and spatial positioning of diversely specialized cell types. Yet, the molecular identities of the cell types and their positions within individual anatomical structures remain incompletely known. To construct a comprehensive atlas of cell types in each brain structure, we paired high-throughput single-nucleus RNA sequencing with Slide-seq1,2-a recently developed spatial transcriptomics method with near-cellular resolution-across the entire mouse brain. Integration of these datasets revealed the cell type composition of each neuroanatomical structure. Cell type diversity was found to be remarkably high in the midbrain, hindbrain and hypothalamus, with most clusters requiring a combination of at least three discrete gene expression markers to uniquely define them. Using these data, we developed a framework for genetically accessing each cell type, comprehensively characterized neuropeptide and neurotransmitter signalling, elucidated region-specific specializations in activity-regulated gene expression and ascertained the heritability enrichment of neurological and psychiatric phenotypes. These data, available as an online resource ( www.BrainCellData.org ), should find diverse applications across neuroscience, including the construction of new genetic tools and the prioritization of specific cell types and circuits in the study of brain diseases.

Spatiotemporal transcriptomic maps of whole mouse embryos at the onset of organogenesis.

Spatiotemporal orchestration of gene expression is required for proper embryonic development. The use of single-cell technologies has begun to provide improved resolution of early regulatory dynamics, including detailed molecular definitions of most cell states during mouse embryogenesis. Here we used Slide-seq to build spatial transcriptomic maps of complete embryonic day (E) 8.5 and E9.0, and partial E9.5 embryos. To support their utility, we developed sc3D, a tool for reconstructing and exploring three-dimensional 'virtual embryos', which enables the quantitative investigation of regionalized gene expression patterns. Our measurements along the main embryonic axes of the developing neural tube revealed several previously unannotated genes with distinct spatial patterns. We also characterized the conflicting transcriptional identity of 'ectopic' neural tubes that emerge in Tbx6 mutant embryos. Taken together, we present an experimental and computational framework for the spatiotemporal investigation of whole embryonic structures and mutant phenotypes.

An atlas of healthy and injured cell states and niches in the human kidney.

Understanding kidney disease relies on defining the complexity of cell types and states, their associated molecular profiles and interactions within tissue neighbourhoods1. Here we applied multiple single-cell and single-nucleus assays (>400,000 nuclei or cells) and spatial imaging technologies to a broad spectrum of healthy reference kidneys (45 donors) and diseased kidneys (48 patients). This has provided a high-resolution cellular atlas of 51 main cell types, which include rare and previously undescribed cell populations. The multi-omic approach provides detailed transcriptomic profiles, regulatory factors and spatial localizations spanning the entire kidney. We also define 28 cellular states across nephron segments and interstitium that were altered in kidney injury, encompassing cycling, adaptive (successful or maladaptive repair), transitioning and degenerative states. Molecular signatures permitted the localization of these states within injury neighbourhoods using spatial transcriptomics, while large-scale 3D imaging analysis (around 1.2 million neighbourhoods) provided corresponding linkages to active immune responses. These analyses defined biological pathways that are relevant to injury time-course and niches, including signatures underlying epithelial repair that predicted maladaptive states associated with a decline in kidney function. This integrated multimodal spatial cell atlas of healthy and diseased human kidneys represents a comprehensive benchmark of cellular states, neighbourhoods, outcome-associated signatures and publicly available interactive visualizations.

The cell type composition of the adult mouse brain revealed by single cell and spatial genomics.

The function of the mammalian brain relies upon the specification and spatial positioning of diversely specialized cell types. Yet, the molecular identities of the cell types, and their positions within individual anatomical structures, remain incompletely known. To construct a comprehensive atlas of cell types in each brain structure, we paired high-throughput single-nucleus RNA-seq with Slide-seq-a recently developed spatial transcriptomics method with near-cellular resolution-across the entire mouse brain. Integration of these datasets revealed the cell type composition of each neuroanatomical structure. Cell type diversity was found to be remarkably high in the midbrain, hindbrain, and hypothalamus, with most clusters requiring a combination of at least three discrete gene expression markers to uniquely define them. Using these data, we developed a framework for genetically accessing each cell type, comprehensively characterized neuropeptide and neurotransmitter signaling, elucidated region-specific specializations in activity-regulated gene expression, and ascertained the heritability enrichment of neurological and psychiatric phenotypes. These data, available as an online resource (BrainCellData.org) should find diverse applications across neuroscience, including the construction of new genetic tools, and the prioritization of specific cell types and circuits in the study of brain diseases.

Dissecting the immune suppressive human prostate tumor microenvironment via integrated single-cell and spatial transcriptomic analyses.

The treatment of low-risk primary prostate cancer entails active surveillance only, while high-risk disease requires multimodal treatment including surgery, radiation therapy, and hormonal therapy. Recurrence and development of metastatic disease remains a clinical problem, without a clear understanding of what drives immune escape and tumor progression. Here, we comprehensively describe the tumor microenvironment of localized prostate cancer in comparison with adjacent normal samples and healthy controls. Single-cell RNA sequencing and high-resolution spatial transcriptomic analyses reveal tumor context dependent changes in gene expression. Our data indicate that an immune suppressive tumor microenvironment associates with suppressive myeloid populations and exhausted T-cells, in addition to high stromal angiogenic activity. We infer cell-to-cell relationships from high throughput ligand-receptor interaction measurements within undissociated tissue sections. Our work thus provides a highly detailed and comprehensive resource of the prostate tumor microenvironment as well as tumor-stromal cell interactions.

TACCO unifies annotation transfer and decomposition of cell identities for single-cell and spatial omics.

Transferring annotations of single-cell-, spatial- and multi-omics data is often challenging owing both to technical limitations, such as low spatial resolution or high dropout fraction, and to biological variations, such as continuous spectra of cell states. Based on the concept that these data are often best described as continuous mixtures of cells or molecules, we present a computational framework for the transfer of annotations to cells and their combinations (TACCO), which consists of an optimal transport model extended with different wrappers to annotate a wide variety of data. We apply TACCO to identify cell types and states, decipher spatiomolecular tissue structure at the cell and molecular level and resolve differentiation trajectories using synthetic and biological datasets. While matching or exceeding the accuracy of specialized tools for the individual tasks, TACCO reduces the computational requirements by up to an order of magnitude and scales to larger datasets (for example, considering the runtime of annotation transfer for 1 M simulated dropout observations).

In-Hospital Virtual Peer-to-Peer Consultation to Increase Guideline-Directed Medical Therapy for Heart Failure: A Pilot Randomized Trial.

BACKGROUND: Guideline-directed medical therapy (GDMT) for heart failure with reduced ejection fraction (HFrEF) improves clinical outcomes and quality of life. Optimizing GDMT in the hospital is associated with greater long-term use in HFrEF. This study aimed to describe the efficacy of a multidisciplinary virtual HF intervention on GDMT optimization among patients with HFrEF admitted for any cause. METHODS: In this pilot randomized, controlled study, consecutive patients with HFrEF admitted to noncardiology medicine services for any cause were identified at a large academic tertiary care hospital between May to September 2021. Major exclusions were end-stage renal disease, hemodynamic instability, concurrent COVID-19 infection, and current enrollment in hospice care. Patients were randomized to a clinician-level virtual peer-to-peer consult intervention providing GDMT recommendations and information on medication costs versus usual care. Primary end points included (1) proportion of patients with new GDMT initiation or use and (2) changes to HF optimal medical therapy scores which included target dosing (range, 0-9). RESULTS: Of 242 patients identified, 91 (38%) were eligible and randomized to intervention (N=52) or usual care (N=39). Baseline characteristics were similar between intervention and usual care (mean age 63 versus 67 years, 23% versus 26% female, 46% versus 49% Black, mean ejection fraction 33% versus 31%). GDMT use on admission was also similar. There were greater proportions of patients with GDMT initiation or continuation with the intervention compared with usual care. After adjusting for optimal medical therapy score on admission, changes to optimal medical therapy score at discharge were higher for the intervention group compared with usual care (+0.44 versus -0.31, absolute difference +0.75, adjusted estimate 0.86±0.42; P=0.041). CONCLUSIONS: Among eligible patients with HFrEF hospitalized for any cause on noncardiology services, a multidisciplinary pilot virtual HF consultation increased new GDMT initiation and dose optimization at discharge.

Reactive organic carbon air emissions from mobile sources in the United States.

Mobile sources are responsible for a substantial controllable portion of the reactive organic carbon (ROC) emitted to the atmosphere, especially in urban environments of the United States. We update existing methods for calculating mobile source organic particle and vapor emissions in the United States with over a decade of laboratory data that parameterize the volatility and organic aerosol (OA) potential of emissions from on-road vehicles, nonroad engines, aircraft, marine vessels, and locomotives. We find that existing emission factor information from Teflon filters combined with quartz filters collapses into simple relationships and can be used to reconstruct the complete volatility distribution of ROC emissions. This new approach consists of source-specific filter artifact corrections and state-of-the-science speciation including explicit intermediate-volatility organic compounds (IVOCs), yielding the first bottom-up volatility-resolved inventory of US mobile source emissions. Using the Community Multiscale Air Quality model, we estimate mobile sources account for 20 %-25 % of the IVOC concentrations and 4.4 %-21.4 % of ambient OA. The updated emissions and air quality model reduce biases in predicting fine-particle organic carbon in winter, spring, and autumn throughout the United States (4.3 %-11.3 % reduction in normalized bias). We identify key uncertain parameters that align with current state-of-the-art research measurement challenges.

Spatial maps of T cell receptors and transcriptomes reveal distinct immune niches and interactions in the adaptive immune response.

T cells mediate antigen-specific immune responses to disease through the specificity and diversity of their clonotypic T cell receptors (TCRs). Determining the spatial distributions of T cell clonotypes in tissues is essential to understanding T cell behavior, but spatial sequencing methods remain unable to profile the TCR repertoire. Here, we developed Slide-TCR-seq, a 10-µm-resolution method, to sequence whole transcriptomes and TCRs within intact tissues. We confirmed the ability of Slide-TCR-seq to map the characteristic locations of T cells and their receptors in mouse spleen. In human lymphoid germinal centers, we identified spatially distinct TCR repertoires. Profiling T cells in renal cell carcinoma and melanoma specimens revealed heterogeneous immune responses: T cell states and infiltration differed intra- and inter-clonally, and adjacent tumor and immune cells exhibited distinct gene expression. Altogether, our method yields insights into the spatial relationships between clonality, neighboring cell types, and gene expression that drive T cell responses.

Cell type-specific inference of differential expression in spatial transcriptomics.

A central problem in spatial transcriptomics is detecting differentially expressed (DE) genes within cell types across tissue context. Challenges to learning DE include changing cell type composition across space and measurement pixels detecting transcripts from multiple cell types. Here, we introduce a statistical method, cell type-specific inference of differential expression (C-SIDE), that identifies cell type-specific DE in spatial transcriptomics, accounting for localization of other cell types. We model gene expression as an additive mixture across cell types of log-linear cell type-specific expression functions. C-SIDE's framework applies to many contexts: DE due to pathology, anatomical regions, cell-to-cell interactions and cellular microenvironment. Furthermore, C-SIDE enables statistical inference across multiple/replicates. Simulations and validation experiments on Slide-seq, MERFISH and Visium datasets demonstrate that C-SIDE accurately identifies DE with valid uncertainty quantification. Last, we apply C-SIDE to identify plaque-dependent immune activity in Alzheimer's disease and cellular interactions between tumor and immune cells. We distribute C-SIDE within the R package https://github.com/dmcable/spacexr .

Proper acquisition of cell class identity in organoids allows definition of fate specification programs of the human cerebral cortex.

Realizing the full utility of brain organoids to study human development requires understanding whether organoids precisely replicate endogenous cellular and molecular events, particularly since acquisition of cell identity in organoids can be impaired by abnormal metabolic states. We present a comprehensive single-cell transcriptomic, epigenetic, and spatial atlas of human cortical organoid development, comprising over 610,000 cells, from generation of neural progenitors through production of differentiated neuronal and glial subtypes. We show that processes of cellular diversification correlate closely to endogenous ones, irrespective of metabolic state, empowering the use of this atlas to study human fate specification. We define longitudinal molecular trajectories of cortical cell types during organoid development, identify genes with predicted human-specific roles in lineage establishment, and uncover early transcriptional diversity of human callosal neurons. The findings validate this comprehensive atlas of human corticogenesis in vitro as a resource to prime investigation into the mechanisms of human cortical development.

Machine learning to define phenotypes and outcomes of patients hospitalized for heart failure with preserved ejection fraction: Findings from ASCEND-HF.

BACKGROUND: Heart Failure with Preserved Ejection Fraction (HFpEF) is a heterogenous disease with few therapies proven to provide clinical benefit. Machine learning can characterize distinct phenotypes and compare outcomes among patients with HFpEF who are hospitalized for acute HF. METHODS: We applied hierarchical clustering using demographics, comorbidities, and clinical data on admission to identify distinct clusters in hospitalized HFpEF (ejection fraction >40%) in the ASCEND-HF trial. We separately applied a previously developed latent class analysis (LCA) clustering method and compared in-hospital and long-term outcomes across cluster groups. RESULTS: Of 7141 patients enrolled in the ASCEND-HF trial, 812 (11.4%) were hospitalized for HFpEF and met the criteria for complete case analysis. Hierarchical Cluster 1 included older women with atrial fibrillation (AF). Cluster 2 had elevated resting blood pressure. Cluster 3 had young men with obesity and diabetes. Cluster 4 had low resting blood pressure. Mortality at 180 days was lowest among Cluster 3 (KM event-rate 6.2 [95% CI: 3.5, 10.9]) and highest among Cluster 4 (18.8 [14.6, 24.0], P < .001). Twenty four-hour urine output was higher in Cluster 3 (2700 mL [1800, 3975]) than Cluster 4 (2100 mL [1400, 3055], P < .001). LCA also identified four clusters: A) older White or Asian women, B) younger men with few comorbidities, C) older individuals with AF and renal impairment, and D) patients with obesity and diabetes. Mortality at 180 days was lowest among LCA Cluster B (KM event-rate 5.5 [2.0, 10.3]) and highest among LCA Cluster C (26.3 [19.2, 35.4], P < .001). CONCLUSIONS: In patients hospitalized for HFpEF, cluster analysis demonstrated distinct phenotypes with differing clinical profiles and outcomes.

High-resolution Slide-seqV2 spatial transcriptomics enables discovery of disease-specific cell neighborhoods and pathways.

High-resolution spatial transcriptomics enables mapping of RNA expression directly from intact tissue sections; however, its utility for the elucidation of disease processes and therapeutically actionable pathways remains unexplored. We applied Slide-seqV2 to mouse and human kidneys, in healthy and distinct disease paradigms. First, we established the feasibility of Slide-seqV2 in tissue from nine distinct human kidneys, which revealed a cell neighborhood centered around a population of LYVE1+ macrophages. Second, in a mouse model of diabetic kidney disease, we detected changes in the cellular organization of the spatially restricted kidney filter and blood-flow-regulating apparatus. Third, in a mouse model of a toxic proteinopathy, we identified previously unknown, disease-specific cell neighborhoods centered around macrophages. In a spatially restricted subpopulation of epithelial cells, we discovered perturbations in 77 genes associated with the unfolded protein response. Our studies illustrate and experimentally validate the utility of Slide-seqV2 for the discovery of disease-specific cell neighborhoods.

Spatial transcriptomic reconstruction of the mouse olfactory glomerular map suggests principles of odor processing.

The olfactory system's ability to detect and discriminate between the vast array of chemicals present in the environment is critical for an animal's survival. In mammals, the first step of this odor processing is executed by olfactory sensory neurons, which project their axons to a stereotyped location in the olfactory bulb (OB) to form glomeruli. The stereotyped positioning of glomeruli in the OB suggests an importance for this organization in odor perception. However, because the location of only a limited subset of glomeruli has been determined, it has been challenging to determine the relationship between glomerular location and odor discrimination. Using a combination of single-cell RNA sequencing, spatial transcriptomics and machine learning, we have generated a map of most glomerular positions in the mouse OB. These observations significantly extend earlier studies and suggest an overall organizational principle in the OB that may be used by the brain to assist in odor decoding.

Physical Rehabilitation in Older Patients Hospitalized with Acute Heart Failure and Diabetes: Insights from REHAB-HF.

BACKGROUND: Prior studies showed an attenuated response to exercise training among patients with heart failure and type 2 diabetes mellitus. We explored the interaction between diabetes status and a novel, transitional, tailored, progressive rehabilitation intervention that improved physical function compared with usual care in the Rehabilitation Therapy in Older Acute Heart Failure Patients (REHAB-HF) trial. METHODS: The effect of the intervention on 3-month Short Physical Performance Battery (SPPB) (primary endpoint), 6-minute walk distance (6MWD), modified Fried frailty criteria, and quality-of-life scores (Kansas City Cardiomyopathy Questionnaire [KCCQ] and EuroQoL Visual Analogue Scale [VAS]) was compared between participants with and without diabetes. Differences in 6-month clinical outcomes were also explored. RESULTS: Of the 349 participants enrolled in REHAB-HF, 186 (53%) had diabetes. The prevalence of diabetes was higher in the intervention group (59% vs 48%). Participants with diabetes had worse baseline physical function by the SPPB and 6MWD, but similar frailty and quality-of-life scores. There was a consistent improvement with the intervention for 3-month SPPB, 6MWD, and VAS regardless of diabetes status (all interaction P value > .6), but participants with diabetes had significantly less improvement for frailty (P = .021) and a trend toward lower improvement in KCCQ (P = .11). There was no significant interaction by diabetes status for 6-month clinical event outcomes (all interaction P value > .3). CONCLUSIONS: Participants with diabetes had worse baseline physical function but showed similar clinically meaningful improvements from the intervention. There was less benefit for frailty with the intervention in participants with diabetes.

Robust decomposition of cell type mixtures in spatial transcriptomics.

A limitation of spatial transcriptomics technologies is that individual measurements may contain contributions from multiple cells, hindering the discovery of cell-type-specific spatial patterns of localization and expression. Here, we develop robust cell type decomposition (RCTD), a computational method that leverages cell type profiles learned from single-cell RNA-seq to decompose cell type mixtures while correcting for differences across sequencing technologies. We demonstrate the ability of RCTD to detect mixtures and identify cell types on simulated datasets. Furthermore, RCTD accurately reproduces known cell type and subtype localization patterns in Slide-seq and Visium datasets of the mouse brain. Finally, we show how RCTD's recovery of cell type localization enables the discovery of genes within a cell type whose expression depends on spatial environment. Spatial mapping of cell types with RCTD enables the spatial components of cellular identity to be defined, uncovering new principles of cellular organization in biological tissue. RCTD is publicly available as an open-source R package at https://github.com/dmcable/RCTD .

Spatial genomics enables multi-modal study of clonal heterogeneity in tissues.

The state and behaviour of a cell can be influenced by both genetic and environmental factors. In particular, tumour progression is determined by underlying genetic aberrations1-4 as well as the makeup of the tumour microenvironment5,6. Quantifying the contributions of these factors requires new technologies that can accurately measure the spatial location of genomic sequence together with phenotypic readouts. Here we developed slide-DNA-seq, a method for capturing spatially resolved DNA sequences from intact tissue sections. We demonstrate that this method accurately preserves local tumour architecture and enables the de novo discovery of distinct tumour clones and their copy number alterations. We then apply slide-DNA-seq to a mouse model of metastasis and a primary human cancer, revealing that clonal populations are confined to distinct spatial regions. Moreover, through integration with spatial transcriptomics, we uncover distinct sets of genes that are associated with clone-specific genetic aberrations, the local tumour microenvironment, or both. Together, this multi-modal spatial genomics approach provides a versatile platform for quantifying how cell-intrinsic and cell-extrinsic factors contribute to gene expression, protein abundance and other cellular phenotypes.

In-Hospital Initiation of Sodium-Glucose Cotransporter-2 Inhibitors for Heart Failure With Reduced Ejection Fraction.

Sodium-glucose cotransporter-2 inhibitor therapy is well suited for initiation during the heart failure hospitalization, owing to clinical benefits that accrue rapidly within days to weeks, a strong safety and tolerability profile, minimal to no effects on blood pressure, and no excess risk of adverse kidney events. There is no evidence to suggest that deferring initiation to the outpatient setting accomplishes anything beneficial. Instead, there is compelling evidence that deferring in-hospital initiation exposes patients to excess risk of early postdischarge clinical worsening and death. Lessons from other heart failure with reduced ejection fraction therapies highlight that deferring initiation of guideline-recommended medications to the U.S. outpatient setting carries a >75% chance they will not be initiated within the next year. Recognizing that 1 in 4 patients hospitalized for worsening heart failure die or are readmitted within 30 days, clinicians should embrace the in-hospital period as an optimal time to initiate sodium-glucose cotransporter-2 inhibitor therapy and treat this population with the urgency it deserves.

Dissecting mammalian spermatogenesis using spatial transcriptomics.

Single-cell RNA sequencing has revealed extensive molecular diversity in gene programs governing mammalian spermatogenesis but fails to delineate their dynamics in the native context of seminiferous tubules, the spatially confined functional units of spermatogenesis. Here, we use Slide-seq, a spatial transcriptomics technology, to generate an atlas that captures the spatial gene expression patterns at near-single-cell resolution in the mouse and human testis. Using Slide-seq data, we devise a computational framework that accurately localizes testicular cell types in individual seminiferous tubules. Unbiased analysis systematically identifies spatially patterned genes and gene programs. Combining Slide-seq with targeted in situ RNA sequencing, we demonstrate significant differences in the cellular compositions of spermatogonial microenvironment between mouse and human testes. Finally, a comparison of the spatial atlas generated from the wild-type and diabetic mouse testis reveals a disruption in the spatial cellular organization of seminiferous tubules as a potential mechanism of diabetes-induced male infertility.