Human BioMolecular Atlas Program (HuBMAP): 3D Human Reference Atlas Construction and Usage.

The Human BioMolecular Atlas Program (HuBMAP) aims to construct a reference 3D structural, cellular, and molecular atlas of the healthy adult human body. The HuBMAP Data Portal (https://portal.hubmapconsortium.org) serves experimental datasets and supports data processing, search, filtering, and visualization. The Human Reference Atlas (HRA) Portal (https://humanatlas.io) provides open access to atlas data, code, procedures, and instructional materials. Experts from more than 20 consortia are collaborating to construct the HRA's Common Coordinate Framework (CCF), knowledge graphs, and tools that describe the multiscale structure of the human body (from organs and tissues down to cells, genes, and biomarkers) and to use the HRA to understand changes that occur at each of these levels with aging, disease, and other perturbations. The 6th release of the HRA v2.0 covers 36 organs with 4,499 unique anatomical structures, 1,195 cell types, and 2,089 biomarkers (e.g., genes, proteins, lipids) linked to ontologies. In addition, three workflows were developed to map new experimental data into the HRA's CCF. This paper describes the HRA user stories, terminology, data formats, ontology validation, unified analysis workflows, user interfaces, instructional materials, application programming interface (APIs), flexible hybrid cloud infrastructure, and demonstrates first atlas usage applications and previews.

SEraster: a rasterization preprocessing framework for scalable spatial omics data analysis.

MOTIVATION: Spatial omics data demand computational analysis but many analysis tools have computational resource requirements that increase with the number of cells analyzed. This presents scalability challenges as researchers use spatial omics technologies to profile millions of cells. RESULTS: To enhance the scalability of spatial omics data analysis, we developed a rasterization preprocessing framework called SEraster that aggregates cellular information into spatial pixels. We apply SEraster to both real and simulated spatial omics data prior to spatial variable gene expression analysis to demonstrate that such preprocessing can reduce computational resource requirements while maintaining high performance, including as compared to other down-sampling approaches. We further integrate SEraster with existing analysis tools to characterize cell-type spatial co-enrichment across length scales. Finally, we apply SEraster to enable analysis of a mouse pup spatial omics dataset with over a million cells to identify tissue-level and cell-type-specific spatially variable genes as well as spatially co-enriched cell types that recapitulate expected organ structures. AVAILABILITY AND IMPLEMENTATION: SEraster is implemented as an R package on GitHub (https://github.com/JEFworks-Lab/SEraster) with additional tutorials at https://JEF.works/SEraster.

Gene count normalization in single-cell imaging-based spatially resolved transcriptomics.

BACKGROUND: Recent advances in imaging-based spatially resolved transcriptomics (im-SRT) technologies now enable high-throughput profiling of targeted genes and their locations in fixed tissues. Normalization of gene expression data is often needed to account for technical factors that may confound underlying biological signals. RESULTS: Here, we investigate the potential impact of different gene count normalization methods with different targeted gene panels in the analysis and interpretation of im-SRT data. Using different simulated gene panels that overrepresent genes expressed in specific tissue regions or cell types, we demonstrate how normalization methods based on detected gene counts per cell differentially impact normalized gene expression magnitudes in a region- or cell type-specific manner. We show that these normalization-induced effects may reduce the reliability of downstream analyses including differential gene expression, gene fold change, and spatially variable gene analysis, introducing false positive and false negative results when compared to results obtained from gene panels that are more representative of the gene expression of the tissue's component cell types. These effects are not observed with normalization approaches that do not use detected gene counts for gene expression magnitude adjustment, such as with cell volume or cell area normalization. CONCLUSIONS: We recommend using non-gene count-based normalization approaches when feasible and evaluating gene panel representativeness before using gene count-based normalization methods if necessary. Overall, we caution that the choice of normalization method and gene panel may impact the biological interpretation of the im-SRT data.

Single-cell morphology encodes functional subtypes of senescence in aging human dermal fibroblasts.

Cellular senescence is an established driver of aging, exhibiting context-dependent phenotypes across multiple biological length-scales. Despite its mechanistic importance, profiling senescence within cell populations is challenging. This is in part due to the limitations of current biomarkers to robustly identify senescent cells across biological settings, and the heterogeneous, non-binary phenotypes exhibited by senescent cells. Using a panel of primary dermal fibroblasts, we combined live single-cell imaging, machine learning, multiple senescence induction conditions, and multiple protein-based senescence biomarkers to show the emergence of functional subtypes of senescence. Leveraging single-cell morphologies, we defined eleven distinct morphology clusters, with the abundance of cells in each cluster being dependent on the mode of senescence induction, the time post-induction, and the age of the donor. Of these eleven clusters, we identified three bona-fide senescence subtypes (C7, C10, C11), with C10 showing the strongest age-dependence across a cohort of fifty aging individuals. To determine the functional significance of these senescence subtypes, we profiled their responses to senotherapies, specifically focusing on Dasatinib + Quercetin (D+Q). Results indicated subtype-dependent responses, with senescent cells in C7 being most responsive to D+Q. Altogether, we provide a robust single-cell framework to identify and classify functional senescence subtypes with applications for next-generation senotherapy screens, and the potential to explain heterogeneous senescence phenotypes across biological settings based on the presence and abundance of distinct senescence subtypes.

Cross-modality mapping using image varifolds to align tissue-scale atlases to molecular-scale measures with application to 2D brain sections.

This paper explicates a solution to building correspondences between molecular-scale transcriptomics and tissue-scale atlases. This problem arises in atlas construction and cross-specimen/technology alignment where specimens per emerging technology remain sparse and conventional image representations cannot efficiently model the high dimensions from subcellular detection of thousands of genes. We address these challenges by representing spatial transcriptomics data as generalized functions encoding position and high-dimensional feature (gene, cell type) identity. We map onto low-dimensional atlas ontologies by modeling regions as homogeneous random fields with unknown transcriptomic feature distribution. We solve simultaneously for the minimizing geodesic diffeomorphism of coordinates through LDDMM and for these latent feature densities. We map tissue-scale mouse brain atlases to gene-based and cell-based transcriptomics data from MERFISH and BARseq technologies and to histopathology and cross-species atlases to illustrate integration of diverse molecular and cellular datasets into a single coordinate system as a means of comparison and further atlas construction.

First-In-Human Phase I Study of Tinengotinib (TT-00420), a Multiple Kinase Inhibitor, as a Single Agent in Patients With Advanced Solid Tumors.

PURPOSE: This first-in-human phase I dose-escalation study evaluated the safety, pharmacokinetics, and efficacy of tinengotinib (TT-00420), a multi-kinase inhibitor targeting fibroblast growth factor receptors 1-3 (FGFRs 1-3), Janus kinase 1/2, vascular endothelial growth factor receptors, and Aurora A/B, in patients with advanced solid tumors. PATIENTS AND METHODS: Patients received tinengotinib orally daily in 28-day cycles. Dose escalation was guided by Bayesian modeling using escalation with overdose control. The primary objective was to assess dose-limiting toxicities (DLTs), maximum tolerated dose (MTD), and dose recommended for dose expansion (DRDE). Secondary objectives included pharmacokinetics and efficacy. RESULTS: Forty-eight patients were enrolled (dose escalation, n = 40; dose expansion, n = 8). MTD was not reached; DRDE was 12 mg daily. DLTs were palmar-plantar erythrodysesthesia syndrome (8 mg, n = 1) and hypertension (15 mg, n = 2). The most common treatment-related adverse event was hypertension (50.0%). In 43 response-evaluable patients, 13 (30.2%) achieved partial response (PR; n = 7) or stable disease (SD) ≥ 24 weeks (n = 6), including 4/11 (36.4%) with FGFR2 mutations/fusions and cholangiocarcinoma (PR n = 3; SD ≥ 24 weeks n = 1), 3/3 (100.0%) with hormone receptor (HR)-positive/HER2-negative breast cancer (PR n = 2; SD ≥ 24 weeks n = 1), 2/5 (40.0%) with triple-negative breast cancer (TNBC; PR n = 1; SD ≥ 24 weeks n = 1), and 1/1 (100.0%) with castrate-resistant prostate cancer (CRPC; PR). Four of 12 patients (33.3%; HR-positive/HER2-negative breast cancer, TNBC, prostate cancer, and cholangiocarcinoma) treated at DRDE had PRs. Tinengotinib's half-life was 28-34 hours. CONCLUSIONS: Tinengotinib was well tolerated with favorable pharmacokinetic characteristics. Preliminary findings indicated potential clinical benefit in FGFR inhibitor-refractory cholangiocarcinoma, HER2-negative breast cancer (including TNBC), and CRPC. Continued evaluation of tinengotinib is warranted in phase II trials.

Soluble PD-L1 reprograms blood monocytes to prevent cerebral edema and facilitate recovery after ischemic stroke.

Acute cerebral ischemia triggers a profound inflammatory response. While macrophages polarized to an M2-like phenotype clear debris and facilitate tissue repair, aberrant or prolonged macrophage activation is counterproductive to recovery. The inhibitory immune checkpoint Programmed Cell Death Protein 1 (PD-1) is upregulated on macrophage precursors (monocytes) in the blood after acute cerebrovascular injury. To investigate the therapeutic potential of PD-1 activation, we immunophenotyped circulating monocytes from patients and found that PD-1 expression was upregulated in the acute period after stroke. Murine studies using a temporary middle cerebral artery (MCA) occlusion (MCAO) model showed that intraperitoneal administration of soluble Programmed Death Ligand-1 (sPD-L1) significantly decreased brain edema and improved overall survival. Mice receiving sPD-L1 also had higher performance scores short-term, and more closely resembled sham animals on assessments of long-term functional recovery. These clinical and radiographic benefits were abrogated in global and myeloid-specific PD-1 knockout animals, confirming PD-1+ monocytes as the therapeutic target of sPD-L1. Single-cell RNA sequencing revealed that treatment skewed monocyte maturation to a non-classical Ly6Clo, CD43hi, PD-L1+ phenotype. These data support peripheral activation of PD-1 on inflammatory monocytes as a therapeutic strategy to treat neuroinflammation after acute ischemic stroke.

Gene count normalization in single-cell imaging-based spatially resolved transcriptomics.

Recent advances in imaging-based spatially resolved transcriptomics (im-SRT) technologies now enable high-throughput profiling of targeted genes and their locations in fixed tissues. Normalization of gene expression data is often needed to account for technical factors that may confound underlying biological signals. Here, we investigate the potential impact of different gene count normalization methods with different targeted gene panels in the analysis and interpretation of im-SRT data. Using different simulated gene panels that overrepresent genes expressed in specific tissue regions or cell types, we demonstrate how normalization methods based on detected gene counts per cell differentially impact normalized gene expression magnitudes in a region- or cell type-specific manner. We show that these normalization-induced effects may reduce the reliability of downstream analyses including differential gene expression, gene fold change, and spatially variable gene analysis, introducing false positive and false negative results when compared to results obtained from gene panels that are more representative of the gene expression of the tissue's component cell types. These effects are not observed with normalization approaches that do not use detected gene counts for gene expression magnitude adjustment, such as with cell volume or cell area normalization. We recommend using non-gene count-based normalization approaches when feasible and evaluating gene panel representativeness before using gene count-based normalization methods if necessary. Overall, we caution that the choice of normalization method and gene panel may impact the biological interpretation of the im-SRT data.

STalign: Alignment of spatial transcriptomics data using diffeomorphic metric mapping.

Spatial transcriptomics (ST) technologies enable high throughput gene expression characterization within thin tissue sections. However, comparing spatial observations across sections, samples, and technologies remains challenging. To address this challenge, we develop STalign to align ST datasets in a manner that accounts for partially matched tissue sections and other local non-linear distortions using diffeomorphic metric mapping. We apply STalign to align ST datasets within and across technologies as well as to align ST datasets to a 3D common coordinate framework. We show that STalign achieves high gene expression and cell-type correspondence across matched spatial locations that is significantly improved over landmark-based affine alignments. Applying STalign to align ST datasets of the mouse brain to the 3D common coordinate framework from the Allen Brain Atlas, we highlight how STalign can be used to lift over brain region annotations and enable the interrogation of compositional heterogeneity across anatomical structures. STalign is available as an open-source Python toolkit at https://github.com/JEFworks-Lab/STalign and as Supplementary Software with additional documentation and tutorials available at https://jef.works/STalign .

Single cell and spatial transcriptomics analysis of kidney double negative T lymphocytes in normal and ischemic mouse kidneys.

T cells are important in the pathogenesis of acute kidney injury (AKI), and TCR+CD4-CD8- (double negative-DN) are T cells that have regulatory properties. However, there is limited information on DN T cells compared to traditional CD4+ and CD8+ cells. To elucidate the molecular signature and spatial dynamics of DN T cells during AKI, we performed single-cell RNA sequencing (scRNA-seq) on sorted murine DN, CD4+, and CD8+ cells combined with spatial transcriptomic profiling of normal and post AKI mouse kidneys. scRNA-seq revealed distinct transcriptional profiles for DN, CD4+, and CD8+ T cells of mouse kidneys with enrichment of Kcnq5, Klrb1c, Fcer1g, and Klre1 expression in DN T cells compared to CD4+ and CD8+ T cells in normal kidney tissue. We validated the expression of these four genes in mouse kidney DN, CD4+ and CD8+ T cells using RT-PCR and Kcnq5, Klrb1, and Fcer1g genes with the NIH human kidney precision medicine project (KPMP). Spatial transcriptomics in normal and ischemic mouse kidney tissue showed a localized cluster of T cells in the outer medulla expressing DN T cell genes including Fcer1g. These results provide a template for future studies in DN T as well as CD4+ and CD8+ cells in normal and diseased kidneys.

Alignment of spatial transcriptomics data using diffeomorphic metric mapping.

Spatial transcriptomics (ST) technologies enable high throughput gene expression characterization within thin tissue sections. However, comparing spatial observations across sections, samples, and technologies remains challenging. To address this challenge, we developed STalign to align ST datasets in a manner that accounts for partially matched tissue sections and other local non-linear distortions using diffeomorphic metric mapping. We apply STalign to align ST datasets within and across technologies as well as to align ST datasets to a 3D common coordinate framework. We show that STalign achieves high gene expression and cell-type correspondence across matched spatial locations that is significantly improved over landmark-based affine alignments. Applying STalign to align ST datasets of the mouse brain to the 3D common coordinate framework from the Allen Brain Atlas, we highlight how STalign can be used to lift over brain region annotations and enable the interrogation of compositional heterogeneity across anatomical structures. STalign is available as an open-source Python toolkit at https://github.com/JEFworks-Lab/STalign and as supplementary software with additional documentation and tutorials available at https://jef.works/STalign.

A Universal Method for Crossing Molecular and Atlas Modalities using Simplex-Based Image Varifolds and Quadratic Programming.

This paper explicates a solution to the problem of building correspondences between molecular-scale transcriptomics and tissue-scale atlases. The central model represents spatial transcriptomics as generalized functions encoding molecular position and high-dimensional transcriptomic-based (gene, cell type) identity. We map onto low-dimensional atlas ontologies by modeling each atlas compartment as a homogeneous random field with unknown transcriptomic feature distribution. The algorithm presented solves simultaneously for the minimizing geodesic diffeomorphism of coordinates and latent atlas transcriptomic feature fractions by alternating LDDMM optimization for coordinate transformations and quadratic programming for the latent transcriptomic variables. We demonstrate the universality of the algorithm in mapping tissue atlases to gene-based and cell-based MERFISH datasets as well as to other tissue scale atlases. The joint estimation of diffeomorphisms and latent feature distributions allows integration of diverse molecular and cellular datasets into a single coordinate system and creates an avenue of comparison amongst atlas ontologies for continued future development.

Reference-free cell type deconvolution of multi-cellular pixel-resolution spatially resolved transcriptomics data.

Recent technological advancements have enabled spatially resolved transcriptomic profiling but at multi-cellular pixel resolution, thereby hindering the identification of cell-type-specific spatial patterns and gene expression variation. To address this challenge, we develop STdeconvolve as a reference-free approach to deconvolve underlying cell types comprising such multi-cellular pixel resolution spatial transcriptomics (ST) datasets. Using simulated as well as real ST datasets from diverse spatial transcriptomics technologies comprising a variety of spatial resolutions such as Spatial Transcriptomics, 10X Visium, DBiT-seq, and Slide-seq, we show that STdeconvolve can effectively recover cell-type transcriptional profiles and their proportional representation within pixels without reliance on external single-cell transcriptomics references. STdeconvolve provides comparable performance to existing reference-based methods when suitable single-cell references are available, as well as potentially superior performance when suitable single-cell references are not available. STdeconvolve is available as an open-source R software package with the source code available at https://github.com/JEFworks-Lab/STdeconvolve .

VeloViz: RNA velocity-informed embeddings for visualizing cellular trajectories.

MOTIVATION: Single-cell transcriptomics profiling technologies enable genome-wide gene expression measurements in individual cells but can currently only provide a static snapshot of cellular transcriptional states. RNA velocity analysis can help infer cell state changes using such single-cell transcriptomics data. To interpret these cell state changes inferred from RNA velocity analysis as part of underlying cellular trajectories, current approaches rely on visualization with principal components, t-distributed stochastic neighbor embedding and other 2D embeddings derived from the observed single-cell transcriptional states. However, these 2D embeddings can yield different representations of the underlying cellular trajectories, hindering the interpretation of cell state changes. RESULTS: We developed VeloViz to create RNA velocity-informed 2D and 3D embeddings from single-cell transcriptomics data. Using both real and simulated data, we demonstrate that VeloViz embeddings are able to capture underlying cellular trajectories across diverse trajectory topologies, even when intermediate cell states may be missing. By considering the predicted future transcriptional states from RNA velocity analysis, VeloViz can help visualize a more reliable representation of underlying cellular trajectories. AVAILABILITY AND IMPLEMENTATION: Source code is available on GitHub (https://github.com/JEFworks-Lab/veloviz) and Bioconductor (https://bioconductor.org/packages/veloviz) with additional tutorials at https://JEF.works/veloviz/. Datasets used can be found on Zenodo (https://doi.org/10.5281/zenodo.4632471). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Single cell biology-a Keystone Symposia report.

Single cell biology has the potential to elucidate many critical biological processes and diseases, from development and regeneration to cancer. Single cell analyses are uncovering the molecular diversity of cells, revealing a clearer picture of the variation among and between different cell types. New techniques are beginning to unravel how differences in cell state-transcriptional, epigenetic, and other characteristics-can lead to different cell fates among genetically identical cells, which underlies complex processes such as embryonic development, drug resistance, response to injury, and cellular reprogramming. Single cell technologies also pose significant challenges relating to processing and analyzing vast amounts of data collected. To realize the potential of single cell technologies, new computational approaches are needed. On March 17-19, 2021, experts in single cell biology met virtually for the Keystone eSymposium "Single Cell Biology" to discuss advances both in single cell applications and technologies.

Rewiring of human neurodevelopmental gene regulatory programs by human accelerated regions.

Human accelerated regions (HARs) are the fastest-evolving regions of the human genome, and many are hypothesized to function as regulatory elements that drive human-specific gene regulatory programs. We interrogate the in vitro enhancer activity and in vivo epigenetic landscape of more than 3,100 HARs during human neurodevelopment, demonstrating that many HARs appear to act as neurodevelopmental enhancers and that sequence divergence at HARs has largely augmented their neuronal enhancer activity. Furthermore, we demonstrate PPP1R17 to be a putative HAR-regulated gene that has undergone remarkable rewiring of its cell type and developmental expression patterns between non-primates and primates and between non-human primates and humans. Finally, we show that PPP1R17 slows neural progenitor cell cycle progression, paralleling the cell cycle length increase seen predominantly in primate and especially human neurodevelopment. Our findings establish HARs as key components in rewiring human-specific neurodevelopmental gene regulatory programs and provide an integrated resource to study enhancer activity of specific HARs.

Interactions between cancer cells and immune cells drive transitions to mesenchymal-like states in glioblastoma.

The mesenchymal subtype of glioblastoma is thought to be determined by both cancer cell-intrinsic alterations and extrinsic cellular interactions, but remains poorly understood. Here, we dissect glioblastoma-to-microenvironment interactions by single-cell RNA sequencing analysis of human tumors and model systems, combined with functional experiments. We demonstrate that macrophages induce a transition of glioblastoma cells into mesenchymal-like (MES-like) states. This effect is mediated, both in vitro and in vivo, by macrophage-derived oncostatin M (OSM) that interacts with its receptors (OSMR or LIFR) in complex with GP130 on glioblastoma cells and activates STAT3. We show that MES-like glioblastoma states are also associated with increased expression of a mesenchymal program in macrophages and with increased cytotoxicity of T cells, highlighting extensive alterations of the immune microenvironment with potential therapeutic implications.