Search

This article focuses on screening the major secreted proteins by the ischemia-challenged cardiac stromal fibroblasts (CF), the assessment of their expression status and functional role in the post-ischemic left ventricle (LV) and in the ischemia-challenged CF culture and to phenotype CF at single cell resolution based on the positivity of the identified mediators. The expression level of CRSP2, HSP27, IL-8, Cofilin-1, and HSP90 in the LV tissues following coronary artery bypass graft (CABG) and myocardial infarction (MI) and CF cells followed the screening profile derived from the MS/MS findings. The histology data unveiled ECM disorganization, inflammation and fibrosis reflecting the ischemic pathology. CRSP2, HSP27, and HSP90 were significantly upregulated in the LV-CABG tissues with a concomitant reduction ion LV-MI whereas Cofilin-1, IL8, Nrf2, and Troponin I were downregulated in LV-CABG and increased in LV-MI. Similar trends were exhibited by ischemic CF. Single cell transcriptomics revealed multiple sub-phenotypes of CF based on their respective upregulation of CRSP2, HSP27, IL-8, Cofilin-1, HSP90, Troponin I and Nrf2 unveiling pathological and pro-healing phenotypes. Further investigations regarding the underlying signaling mechanisms and validation of sub-populations would offer novel translational avenues for the management of cardiac diseases.

Unsolved Mendelian cases often lack obvious pathogenic coding variants, suggesting potential non-coding etiologies. Here, we present a single cell multi-omic framework integrating embryonic mouse chromatin accessibility, histone modification, and gene expression assays to discover cranial motor neuron (cMN) cis-regulatory elements and subsequently nominate candidate non-coding variants in the congenital cranial dysinnervation disorders (CCDDs), a set of Mendelian disorders altering cMN development. We generate single cell epigenomic profiles for ~86,000 cMNs and related cell types, identifying ~250,000 accessible regulatory elements with cognate gene predictions for ~145,000 putative enhancers. We evaluate enhancer activity for 59 elements using an in vivo transgenic assay and validate 44 (75%), demonstrating that single cell accessibility can be a strong predictor of enhancer activity. Applying our cMN atlas to 899 whole genome sequences from 270 genetically unsolved CCDD pedigrees, we achieve significant reduction in our variant search space and nominate candidate variants predicted to regulate known CCDD disease genes MAFB, PHOX2A, CHN1, and EBF3 - as well as candidates in recurrently mutated enhancers through peak- and gene-centric allelic aggregation. This work delivers non-coding variant discoveries of relevance to CCDDs and a generalizable framework for nominating non-coding variants of potentially high functional impact in other Mendelian disorders.

In our cells, a limited number of RNA binding proteins (RBPs) are responsible for all aspects of RNA metabolism across the entire transcriptome. To accomplish this, RBPs form regulatory units that act on specific target regulons. However, the landscape of RBP combinatorial interactions remains poorly explored. Here, we perform a systematic annotation of RBP combinatorial interactions via multimodal data integration. We build a large-scale map of RBP protein neighborhoods by generating in vivo proximity-dependent biotinylation datasets of 50 human RBPs. In parallel, we use CRISPR interference with single-cell readout to capture transcriptomic changes upon RBP knockdowns. By combining these physical and functional interaction readouts, along with the atlas of RBP mRNA targets from eCLIP assays, we generate an integrated map of functional RBP interactions. We then use this map to match RBPs to their context-specific functions and validate the predicted functions biochemically for four RBPs. This study provides a detailed map of RBP interactions and deconvolves them into distinct regulatory modules with annotated functions and target regulons. This multimodal and integrative framework provides a principled approach for studying post-transcriptional regulatory processes and enriches our understanding of their underlying mechanisms.

Inferring gene regulatory networks from single-cell RNA-sequencing trajectories has been an active area of research yet methods are still needed to identify regulators governing cell transitions. We developed DREAMIT (Dynamic Regulation of Expression Across Modules in Inferred Trajectories) to annotate transcription-factor activity along single-cell trajectory branches, using ensembles of relations to target genes. Using a benchmark representing several different tissues, as well as external validation with ATAC-Seq and Perturb-Seq data on hematopoietic cells, the method was found to have higher tissue-specific sensitivity and specificity over competing approaches.

Cranial sutures separate neighboring skull bones and are sites of bone growth. A key question is how osteogenic activity is controlled to promote bone growth while preventing aberrant bone fusions during skull expansion. Using single-cell transcriptomics, lineage tracing, and mutant analysis in zebrafish, we uncover key developmental transitions regulating bone formation at sutures during skull expansion. In particular, we identify a subpopulation of mesenchyme cells in the mid-suture region that upregulate a suite of genes including BMP antagonists (e.g. grem1a) and pro-angiogenic factors. Lineage tracing with grem1a:nlsEOS reveals that this mid-suture subpopulation is largely non-osteogenic. Moreover, combinatorial mutation of BMP antagonists enriched in this mid-suture subpopulation results in increased BMP signaling in the suture, misregulated bone formation, and abnormal suture morphology. These data reveal establishment of a non-osteogenic mesenchyme population in the mid-suture region that restricts bone formation through local BMP antagonism, thus ensuring proper suture morphology.

The development of the retina is under tight temporal and spatial control. To gain insights into the molecular basis of this process, we generate a single-nuclei dual-omic atlas of the human developing retina with approximately 220,000 nuclei from 14 human embryos and fetuses aged between 8 and 23-weeks post-conception with matched macular and peripheral tissues. This atlas captures all major cell classes in the retina, along with a large proportion of progenitors and cell-type-specific precursors. Cell trajectory analysis reveals a transition from continuous progression in early progenitors to a hierarchical development during the later stages of cell type specification. Both known and unrecorded candidate transcription factors, along with gene regulatory networks that drive the transitions of various cell fates, are identified. Comparisons between the macular and peripheral retinae indicate a largely consistent yet distinct developmental pattern. This atlas offers unparalleled resolution into the transcriptional and chromatin accessibility landscapes during development, providing an invaluable resource for deeper insights into retinal development and associated diseases.

Merging diverse single-cell RNA sequencing (scRNA-seq) data from numerous experiments, laboratories and technologies can uncover important biological insights. Nonetheless, integrating scRNA-seq data encounters special challenges when the datasets are composed of diverse cell type compositions. Scanorama offers a robust solution for improving the quality and interpretation of heterogeneous scRNA-seq data by effectively merging information from diverse sources. Scanorama is designed to address the technical variation introduced by differences in sample preparation, sequencing depth and experimental batches that can confound the analysis of multiple scRNA-seq datasets. Here we provide a detailed protocol for using Scanorama within a Scanpy-based single-cell analysis workflow coupled with Google Colaboratory, a cloud-based free Jupyter notebook environment service. The protocol involves Scanorama integration, a process that typically spans 0.5-3 h. Scanorama integration requires a basic understanding of cellular biology, transcriptomic technologies and bioinformatics. Our protocol and new Scanorama-Colaboratory resource should make scRNA-seq integration more widely accessible to researchers.

With the increasing incidence of kidney diseases, there is an urgent need to develop therapeutic strategies to combat post-injury fibrosis. Immune cells, including platelets, play a pivotal role in this repair process, primarily through their released cytokines. However, the specific role of platelets in kidney injury and subsequent repair remains underexplored. Here, the detrimental role of platelets in renal recovery following ischemia/reperfusion injury and its contribution to acute kidney injury  to chronic kidney disease transition is aimed to investigated. In this study, it is shown that depleting platelets accelerates injury resolution and significantly reduces fibrosis. Employing advanced single-cell and spatial transcriptomic techniques, macrophages as the primary mediators modulated by platelet signals is identified. A novel subset of macrophages, termed cycling M2, which exhibit an M2 phenotype combined with enhanced proliferative activity is uncovered. This subset emerges in the injured kidney during the resolution phase and is modulated by platelet-derived thrombospondin 1 (THBS1) signaling, acquiring profibrotic characteristics. Conversely, targeted inhibition of THBS1 markedly downregulates the cycling M2 macrophage, thereby mitigating fibrotic progression. Overall, this findings highlight the adverse role of platelet THBS1-boosted cycling M2 macrophages in renal injury repair and suggest platelet THBS1 as a promising therapeutic target for alleviating inflammation and kidney fibrosis.

Standard single-cell (sc) proteomics of disease states inferred from multicellular organs or organoids cannot currently be related to single-cell physiology. Here, a scPatch-Clamp/Proteomics platform is developed on single neurons generated from hiPSCs bearing an Alzheimer's disease (AD) genetic mutation and compares them to isogenic wild-type controls. This approach provides both current and voltage electrophysiological data plus detailed proteomics information on single-cells. With this new method, the authors are able to observe hyperelectrical activity in the AD hiPSC-neurons, similar to that observed in the human AD brain, and correlate it to ≈1400 proteins detected at the single neuron level. Using linear regression and mediation analyses to explore the relationship between the abundance of individual proteins and the neuron's mutational and electrophysiological status, this approach yields new information on therapeutic targets in excitatory neurons not attainable by traditional methods. This combined patch-proteomics technique creates a new proteogenetic-therapeutic strategy to correlate genotypic alterations to physiology with protein expression in single-cells.

The contrast between the disruption of genome topology after cohesin loss and the lack of downstream gene expression changes instigates intense debates regarding the structure-function relationship between genome and gene regulation. Here, by analyzing transcriptome and chromatin accessibility at the single-cell level, we discover that, instead of dictating population-wide gene expression levels, cohesin supplies a general function to neutralize stochastic coexpression tendencies of cis-linked genes in single cells. Notably, cohesin loss induces widespread gene coactivation and chromatin co-opening tens of million bases apart in cis. Spatial genome and protein imaging reveals that cohesin prevents gene co-bursting along the chromosome and blocks spatial mixing of transcriptional hubs. Single-molecule imaging shows that cohesin confines the exploration of diverse enhancer and core promoter binding transcriptional regulators. Together, these results support that cohesin arranges nuclear topology to control gene coexpression in single cells.

Natural killer (NK) cells are innate lymphoid cells (ILCs) contributing to immune responses to microbes and tumors. Historically, their classification hinged on a limited array of surface protein markers. Here, we used single-cell RNA sequencing (scRNA-seq) and cellular indexing of transcriptomes and epitopes by sequencing (CITE-seq) to dissect the heterogeneity of NK cells. We identified three prominent NK cell subsets in healthy human blood: NK1, NK2 and NK3, further differentiated into six distinct subgroups. Our findings delineate the molecular characteristics, key transcription factors, biological functions, metabolic traits and cytokine responses of each subgroup. These data also suggest two separate ontogenetic origins for NK cells, leading to divergent transcriptional trajectories. Furthermore, we analyzed the distribution of NK cell subsets in the lung, tonsils and intraepithelial lymphocytes isolated from healthy individuals and in 22 tumor types. This standardized terminology aims at fostering clarity and consistency in future research, thereby improving cross-study comparisons.

Single-gene missense mutations remain challenging to interpret. Here, we deploy scalable functional screening by sequencing (SEUSS), a Perturb-seq method, to generate mutations at protein interfaces of RUNX1 and quantify their effect on activities of downstream cellular programs. We evaluate single-cell RNA profiles of 115 mutations in myelogenous leukemia cells and categorize them into three functionally distinct groups, wild-type (WT)-like, loss-of-function (LoF)-like, and hypomorphic, that we validate in orthogonal assays. LoF-like variants dominate the DNA-binding site and are recurrent in cancer; however, recurrence alone does not predict functional impact. Hypomorphic variants share characteristics with LoF-like but favor protein interactions, promoting gene expression indicative of nerve growth factor (NGF) response and cytokine recruitment of neutrophils. Accessible DNA near differentially expressed genes frequently contains RUNX1-binding motifs. Finally, we reclassify 16 variants of uncertain significance and train a classifier to predict 103 more. Our work demonstrates the potential of targeting protein interactions to better define the landscape of phenotypes reachable by missense mutations.

Archived patient-derived tissue specimens play a central role in understanding disease and developing therapies. To address specificity and sensitivity shortcomings of existing single-cell resolution proteoform analysis tools, we introduce a hybrid microfluidic platform (DropBlot) designed for proteoform analyses in chemically fixed single cells. DropBlot serially integrates droplet-based encapsulation and lysis of single fixed cells, with on-chip microwell-based antigen retrieval, with single-cell western blotting of target antigens. A water-in-oil droplet formulation withstands the harsh chemical (SDS, 6 M urea) and thermal conditions (98 °C, 1-2 hr) required for effective antigen retrieval, and supports analysis of retrieved protein targets by single-cell electrophoresis. We demonstrate protein-target retrieval from unfixed, paraformaldehyde-fixed (PFA), and methanol-fixed cells. Key protein targets (HER2, GAPDH, EpCAM, Vimentin) retrieved from PFA-fixed cells were resolved and immunoreactive. Relevant to biorepositories, DropBlot profiled targets retrieved from human-derived breast tumor specimens archived for six years, offering a workflow for single-cell protein-biomarker analysis of sparing biospecimens.

Melanoma clinical outcomes emerge from incompletely understood genetic mechanisms operating within the tumor and its microenvironment. Here, we used single-cell RNA-based spatial molecular imaging (RNA-SMI) in patient-derived archival tumors to reveal clinically relevant markers of malignancy progression and prognosis. We examined spatial gene expression of 203,472 cells inside benign and malignant melanocytic neoplasms, including melanocytic nevi and primary invasive and metastatic melanomas. Algorithmic cell clustering paired with intratumoral comparative two-dimensional analyses visualized synergistic, spatial gene signatures linking cellular proliferation, metabolism, and malignancy, validated by protein expression. Metastatic niches included up-regulation of CDK2 and FABP5, which independently predicted poor clinical outcome in 473 patients with melanoma via Cox regression analysis. More generally, our work demonstrates a framework for applying single-cell RNA-SMI technology toward identifying gene regulatory landscapes pertinent to cancer progression and patient survival.

Bone marrow aspirate concentrate (BMAC) and adipose-derived stromal vascular fraction (ADSVF) are the most marketed stem cell therapies to treat a variety of conditions in the general population and elite athletes. Both tissues have been used interchangeably clinically even though their detailed composition, heterogeneity, and mechanisms of action have neither been rigorously inventoried nor compared. This lack of information has prevented investigations into ideal dosages and has facilitated anecdata and misinformation. Here, we analyzed single-cell transcriptomes, proteomes, and flow cytometry profiles from paired clinical-grade BMAC and ADSVF. This comparative transcriptional atlas challenges the prevalent notion that there is one therapeutic cell type present in both tissues. We also provide data of surface markers that may enable isolation and investigation of cell (sub)populations. Furthermore, the proteome atlas highlights intertissue and interpatient heterogeneity of injected proteins with potentially regenerative or immunomodulatory capacities. An interactive webtool is available online.

Glioblastoma (GBM) is the most common malignant primary brain tumor and remains incurable. Previous work has shown that systemic administration of Decitabine (DAC) induces sufficient expression of cancer-testis antigens (CTA) in GBM for targeting by adoptive T-cell therapy in vivo. However, the mechanisms by which DAC enhances immunogenicity in GBM remain to be elucidated. Using New York esophageal squamous cell carcinoma 1 (NY-ESO-1) as a representative inducible CTA, we demonstrate in patient tissue, immortalized glioma cells, and primary patient-derived gliomaspheres that basal CTA expression is restricted by promoter hypermethylation in gliomas. DAC treatment of glioma cells specifically inhibits DNA methylation silencing to render NY-ESO-1 and other CTA into inducible tumor antigens at single-cell resolution. Functionally, NY-ESO-1 T-cell receptor-engineered effector cell targeting of DAC-induced antigen in primary glioma cells promotes specific and polyfunctional T-cell cytokine profiles. In addition to induction of CTA, DAC concomitantly reactivates tumor-intrinsic human endogenous retroviruses, interferon response signatures, and MHC-I. Overall, we demonstrate that DAC induces targetable tumor antigen and enhances T-cell functionality against GBM, ultimately contributing to the improvement of targeted immune therapies in glioma.SignificanceThis study dissects the tumor-intrinsic epigenetic and transcriptional mechanisms underlying enhanced T-cell functionality targeting decitabine-induced cancer-testis antigens in glioma. Our findings demonstrate concomitant induction of tumor antigens, reactivation of human endogenous retroviruses, and stimulation of interferon signaling as a mechanistic rationale to epigenetically prime human gliomas to immunotherapeutic targeting.

The circadian clock is a critical regulator of immunity, and this circadian control of immune modulation has an essential function in host defense and tumor immunosurveillance. Here we use a single-cell RNA sequencing approach and a genetic model of colorectal cancer to identify clock-dependent changes to the immune landscape that control the abundance of immunosuppressive cells and consequent suppression of cytotoxic CD8+ T cells. Of these immunosuppressive cell types, PD-L1-expressing myeloid-derived suppressor cells (MDSCs) peak in abundance in a rhythmic manner. Disruption of the epithelial cell clock regulates the secretion of cytokines that promote heightened inflammation, recruitment of neutrophils and the subsequent development of MDSCs. We also show that time-of-day anti-PD-L1 delivery is most effective when synchronized with the abundance of immunosuppressive MDSCs. Collectively, these data indicate that circadian gating of tumor immunosuppression informs the timing and efficacy of immune checkpoint inhibitors.

Cardiomyocytes derived from human induced pluripotent stem cells (hiPSC-CMs) are powerful in vitro models to study the mechanisms underlying cardiomyopathies and cardiotoxicity. Quantification of the contractile function in single hiPSC-CMs at high-throughput and over time is essential to disentangle how cellular mechanisms affect heart function. Here, we present CONTRAX, an open-access, versatile, and streamlined pipeline for quantitative tracking of the contractile dynamics of single hiPSC-CMs over time. Three software modules enable: parameter-based identification of single hiPSC-CMs; automated video acquisition of >200 cells/hour; and contractility measurements via traction force microscopy. We analyze >4,500 hiPSC-CMs over time in the same cells under orthogonal conditions of culture media and substrate stiffnesses; +/- drug treatment; +/- cardiac mutations. Using undirected clustering, we reveal converging maturation patterns, quantifiable drug response to Mavacamten and significant deficiencies in hiPSC-CMs with disease mutations. CONTRAX empowers researchers with a potent quantitative approach to develop cardiac therapies.

Facioscapulohumeral muscular dystrophy (FSHD) is linked to abnormal derepression of the transcription activator DUX4. This effect is localized to a low percentage of cells, requiring single-cell analysis. However, single-cell/nucleus RNA-seq cannot fully capture the transcriptome of multinucleated large myotubes. To circumvent these issues, we use multiplexed error-robust fluorescent in situ hybridization (MERFISH) spatial transcriptomics that allows profiling of RNA transcripts at a subcellular resolution. We simultaneously examined spatial distributions of 140 genes, including 24 direct DUX4 targets, in in vitro differentiated myotubes and unfused mononuclear cells (MNCs) of control, isogenic D4Z4 contraction mutant and FSHD patient samples, as well as the individual nuclei within them. We find myocyte nuclei segregate into two clusters defined by the expression of DUX4 target genes, which is exclusively found in patient/mutant nuclei, whereas MNCs cluster based on developmental states. Patient/mutant myotubes are found in FSHD-hi and FSHD-lo states with the former signified by high DUX4 target expression and decreased muscle gene expression. Pseudotime analyses reveal a clear bifurcation of myoblast differentiation into control and FSHD-hi myotube branches, with variable numbers of DUX4 target-expressing nuclei found in multinucleated FSHD-hi myotubes. Gene coexpression modules related to extracellular matrix and stress gene ontologies are significantly altered in patient/mutant myotubes compared with the control. We also identify distinct subpathways within the DUX4 gene network that may differentially contribute to the disease transcriptomic phenotype. Taken together, our MERFISH-based study provides effective gene network profiling of multinucleated cells and identifies FSHD-induced transcriptomic alterations during myoblast differentiation.

Juvenile dermatomyositis (JDM) is one of several childhood-onset autoimmune disorders characterized by a type I IFN response and autoantibodies. Treatment options are limited due to an incomplete understanding of how the disease emerges from dysregulated cell states across the immune system. We therefore investigated the blood of patients with JDM at different stages of disease activity using single-cell transcriptomics paired with surface protein expression. By immunophenotyping peripheral blood mononuclear cells, we observed skewing of the B cell compartment toward an immature naive state as a hallmark of JDM at diagnosis. Furthermore, we find that these changes in B cells are paralleled by T cell signatures suggestive of Th2-mediated inflammation that persist despite disease quiescence. We applied network analysis to reveal that hyperactivation of the type I IFN response in all immune populations is coordinated with previously masked cell states including dysfunctional protein processing in CD4+ T cells and regulation of cell death programming in NK cells, CD8+ T cells, and γδ T cells. Together, these findings unveil the coordinated immune dysregulation underpinning JDM and provide insight into strategies for restoring balance in immune function.

The functional state of cells is dependent on their microenvironmental context. Prior studies described how polarizing cytokines alter macrophage transcriptomes and epigenomes. Here, we characterized the functional responses of 6 differentially polarized macrophage populations by measuring the dynamics of transcription factor nuclear factor κB (NF-κB) in response to 8 stimuli. The resulting dataset of single-cell NF-κB trajectories was analyzed by three approaches: (1) machine learning on time-series data revealed losses of stimulus distinguishability with polarization, reflecting canalized effector functions. (2) Informative trajectory features driving stimulus distinguishability (signaling codons) were identified and used for mapping a cell state landscape that could then locate macrophages conditioned by an unrelated condition. (3) Kinetic parameters, inferred using a mechanistic NF-κB network model, provided an alternative mapping of cell states and correctly predicted biochemical findings. Together, this work demonstrates that a single analytes dynamic trajectories may distinguish the functional states of single cells and molecular network states underlying them. A record of this papers transparent peer review process is included in the supplemental information.

The cellular components of tumors and their microenvironment play pivotal roles in tumor progression, patient survival, and the response to cancer treatments. Unveiling a comprehensive cellular profile within bulk tumors via single-cell RNA sequencing (scRNA-seq) data is crucial, as it unveils intrinsic tumor cellular traits that elude identification through conventional cancer subtyping methods. Our contribution, scBeacon, is a tool that derives cell-type signatures by integrating and clustering multiple scRNA-seq datasets to extract signatures for deconvolving unrelated tumor datasets on bulk samples. Through the employment of scBeacon on the The Cancer Genome Atlas (TCGA) cohort, we find cellular and molecular attributes within specific tumor categories, many with patient outcome relevance. We developed a tumor cell-type map to visually depict the relationships among TCGA samples based on the cell-type inferences.

In vivo molecular imaging tools are crucially important for elucidating how cells move through complex biological systems; however, achieving single-cell sensitivity over the entire body remains challenging. Here, we report a highly sensitive and multiplexed approach for tracking upward of 20 single cells simultaneously in the same subject using positron emission tomography (PET). The method relies on a statistical tracking algorithm (PEPT-EM) to achieve a sensitivity of 4 becquerel per cell and a streamlined workflow to reliably label single cells with over 50 becquerel per cell of 18F-fluorodeoxyglucose (FDG). To demonstrate the potential of the method, we tracked the fate of more than 70 melanoma cells after intracardiac injection and found they primarily arrested in the small capillaries of the pulmonary, musculoskeletal, and digestive organ systems. This study bolsters the evolving potential of PET in offering unmatched insights into the earliest phases of cell trafficking in physiological and pathological processes and in cell-based therapies.

MOTIVATION: Recent advances in spatial transcriptomics allow spatially resolved gene expression measurements with cellular or even sub-cellular resolution, directly characterizing the complex spatiotemporal gene expression landscape and cell-to-cell interactions in their native microenvironments. Due to technology limitations, most spatial transcriptomic technologies still yield incomplete expression measurements with excessive missing values. Therefore, gene imputation is critical to filling in missing data, enhancing resolution, and improving overall interpretability. However, existing methods either require additional matched single-cell RNA-seq data, which is rarely available, or ignore spatial proximity or expression similarity information. RESULTS: To address these issues, we introduce Impeller, a path-based heterogeneous graph learning method for spatial transcriptomic data imputation. Impeller has two unique characteristics distinct from existing approaches. First, it builds a heterogeneous graph with two types of edges representing spatial proximity and expression similarity. Therefore, Impeller can simultaneously model smooth gene expression changes across spatial dimensions and capture similar gene expression signatures of faraway cells from the same type. Moreover, Impeller incorporates both short- and long-range cell-to-cell interactions (e.g. via paracrine and endocrine) by stacking multiple GNN layers. We use a learnable path operator in Impeller to avoid the over-smoothing issue of the traditional Laplacian matrices. Extensive experiments on diverse datasets from three popular platforms and two species demonstrate the superiority of Impeller over various state-of-the-art imputation methods. AVAILABILITY AND IMPLEMENTATION: The code and preprocessed data used in this study are available at https://github.com/aicb-ZhangLabs/Impeller and https://zenodo.org/records/11212604.

The inability to scalably and precisely measure the activity of developmental cis-regulatory elements (CREs) in multicellular systems is a bottleneck in genomics. Here we develop a dual RNA cassette that decouples the detection and quantification tasks inherent to multiplex single-cell reporter assays. The resulting measurement of reporter expression is accurate over multiple orders of magnitude, with a precision approaching the limit set by Poisson counting noise. Together with RNA barcode stabilization via circularization, these scalable single-cell quantitative expression reporters provide high-contrast readouts, analogous to classic in situ assays but entirely from sequencing. Screening >200 regions of accessible chromatin in a multicellular in vitro model of early mammalian development, we identify 13 (8 previously uncharacterized) autonomous and cell-type-specific developmental CREs. We further demonstrate that chimeric CRE pairs generate cognate two-cell-type activity profiles and assess gain- and loss-of-function multicellular expression phenotypes from CRE variants with perturbed transcription factor binding sites. Single-cell quantitative expression reporters can be applied in developmental and multicellular systems to quantitatively characterize native, perturbed and synthetic CREs at scale, with high sensitivity and at single-cell resolution.

Intratumor heterogeneity underlies cancer evolution and treatment resistance, but targetable mechanisms driving intratumor heterogeneity are poorly understood. Meningiomas are the most common primary intracranial tumors and are resistant to all medical therapies, and high-grade meningiomas have significant intratumor heterogeneity. Here we use spatial approaches to identify genomic, biochemical and cellular mechanisms linking intratumor heterogeneity to the molecular, temporal and spatial evolution of high-grade meningiomas. We show that divergent intratumor gene and protein expression programs distinguish high-grade meningiomas that are otherwise grouped together by current classification systems. Analyses of matched pairs of primary and recurrent meningiomas reveal spatial expansion of subclonal copy number variants associated with treatment resistance. Multiplexed sequential immunofluorescence and deconvolution of meningioma spatial transcriptomes using cell types from single-cell RNA sequencing show decreased immune infiltration, decreased MAPK signaling, increased PI3K-AKT signaling and increased cell proliferation, which are associated with meningioma recurrence. To translate these findings to preclinical models, we use CRISPR interference and lineage tracing approaches to identify combination therapies that target intratumor heterogeneity in meningioma cell co-cultures.

Single-cell assay for transposase-accessible chromatin by sequencing (scATAC-seq) has emerged as a powerful tool for dissecting regulatory landscapes and cellular heterogeneity. However, an exploration of systemic biases among scATAC-seq technologies has remained absent. In this study, we benchmark the performance of eight scATAC-seq methods across 47 experiments using human peripheral blood mononuclear cells (PBMCs) as a reference sample and develop PUMATAC, a universal preprocessing pipeline, to handle the various sequencing data formats. Our analyses reveal significant differences in sequencing library complexity and tagmentation specificity, which impact cell-type annotation, genotype demultiplexing, peak calling, differential region accessibility and transcription factor motif enrichment. Our findings underscore the importance of sample extraction, method selection, data processing and total cost of experiments, offering valuable guidance for future research. Finally, our data and analysis pipeline encompasses 169,000 PBMC scATAC-seq profiles and a best practices code repository for scATAC-seq data analysis, which are freely available to extend this benchmarking effort to future protocols.

Cell atlases serve as vital references for automating cell labeling in new samples, yet existing classification algorithms struggle with accuracy. Here we introduce SIMS (scalable, interpretable machine learning for single cell), a low-code data-efficient pipeline for single-cell RNA classification. We benchmark SIMS against datasets from different tissues and species. We demonstrate SIMS's efficacy in classifying cells in the brain, achieving high accuracy even with small training sets (

Single-cell genomics is a powerful tool for studying heterogeneous tissues such as the brain. Yet little is understood about how genetic variants influence cell-level gene expression. Addressing this, we uniformly processed single-nuclei, multiomics datasets into a resource comprising >2.8 million nuclei from the prefrontal cortex across 388 individuals. For 28 cell types, we assessed population-level variation in expression and chromatin across gene families and drug targets. We identified >550,000 cell type-specific regulatory elements and >1.4 million single-cell expression quantitative trait loci, which we used to build cell-type regulatory and cell-to-cell communication networks. These networks manifest cellular changes in aging and neuropsychiatric disorders. We further constructed an integrative model accurately imputing single-cell expression and simulating perturbations; the model prioritized ~250 disease-risk genes and drug targets with associated cell types.

BACKGROUND: Single-cell transcriptome sequencing (scRNA-Seq) has allowed new types of investigations at unprecedented levels of resolution. Among the primary goals of scRNA-Seq is the classification of cells into distinct types. Many approaches build on existing clustering literature to develop tools specific to single-cell. However, almost all of these methods rely on heuristics or user-supplied parameters to control the number of clusters. This affects both the resolution of the clusters within the original dataset as well as their replicability across datasets. While many recommendations exist, in general, there is little assurance that any given set of parameters will represent an optimal choice in the trade-off between cluster resolution and replicability. For instance, another set of parameters may result in more clusters that are also more replicable. RESULTS: Here, we propose Dune, a new method for optimizing the trade-off between the resolution of the clusters and their replicability. Our method takes as input a set of clustering results-or partitions-on a single dataset and iteratively merges clusters within each partitions in order to maximize their concordance between partitions. As demonstrated on multiple datasets from different platforms, Dune outperforms existing techniques, that rely on hierarchical merging for reducing the number of clusters, in terms of replicability of the resultant merged clusters as well as concordance with ground truth. Dune is available as an R package on Bioconductor: https://www.bioconductor.org/packages/release/bioc/html/Dune.html . CONCLUSIONS: Cluster refinement by Dune helps improve the robustness of any clustering analysis and reduces the reliance on tuning parameters. This method provides an objective approach for borrowing information across multiple clusterings to generate replicable clusters most likely to represent common biological features across multiple datasets.

Dorsal interneurons (dIs) in the spinal cord encode the perception of touch, pain, heat, itchiness and proprioception. Previous studies using genetic strategies in animal models have revealed important insights into dI development, but the molecular details of how dIs arise as distinct populations of neurons remain incomplete. We have developed a resource to investigate dI fate specification by combining a single-cell RNA-Seq atlas of mouse embryonic stem cell-derived dIs with pseudotime analyses. To validate this in silico resource as a useful tool, we used it to first identify genes that are candidates for directing the transition states that lead to distinct dI lineage trajectories, and then validated them using in situ hybridization analyses in the developing mouse spinal cord in vivo. We have also identified an endpoint of the dI5 lineage trajectory and found that dIs become more transcriptionally homogeneous during terminal differentiation. This study introduces a valuable tool for further discovery about the timing of gene expression during dI differentiation and demonstrates its utility in clarifying dI lineage relationships.

BACKGROUND: Alzheimers disease (AD) is a devastating neurodegenerative disorder affecting 44 million people worldwide, leading to cognitive decline, memory loss, and significant impairment in daily functioning. The recent single-cell sequencing technology has revolutionized genetic and genomic resolution by enabling scientists to explore the diversity of gene expression patterns at the finest resolution. Most existing studies have solely focused on molecular perturbations within each cell, but cells live in microenvironments rather than in isolated entities. Here, we leveraged the large-scale and publicly available single-nucleus RNA sequencing in the human prefrontal cortex to investigate cell-to-cell communication in healthy brains and their perturbations in AD. We uniformly processed the snRNA-seq with strict QCs and labeled canonical cell types consistent with the definitions from the BRAIN Initiative Cell Census Network. From ligand and receptor gene expression, we built a high-confidence cell-to-cell communication network to investigate signaling differences between AD and healthy brains. RESULTS: Specifically, we first performed broad communication pattern analyses to highlight that biologically related cell types in normal brains rely on largely overlapping signaling networks and that the AD brain exhibits the irregular inter-mixing of cell types and signaling pathways. Secondly, we performed a more focused cell-type-centric analysis and found that excitatory neurons in AD have significantly increased their communications to inhibitory neurons, while inhibitory neurons and other non-neuronal cells globally decreased theirs to all cells. Then, we delved deeper with a signaling-centric view, showing that canonical signaling pathways CSF, TGFβ, and CX3C are significantly dysregulated in their signaling to the cell type microglia/PVM and from endothelial to neuronal cells for the WNT pathway. Finally, after extracting 23 known AD risk genes, our intracellular communication analysis revealed a strong connection of extracellular ligand genes APP, APOE, and PSEN1 to intracellular AD risk genes TREM2, ABCA1, and APP in the communication from astrocytes and microglia to neurons. CONCLUSIONS: In summary, with the novel advances in single-cell sequencing technologies, we show that cellular signaling is regulated in a cell-type-specific manner and that improper regulation of extracellular signaling genes is linked to intracellular risk genes, giving the mechanistic intra- and inter-cellular picture of AD.

Debate remains around the anatomical origins of specific brain cell subtypes and lineage relationships within the human forebrain1-7. Thus, direct observation in the mature human brain is critical for a complete understanding of its structural organization and cellular origins. Here we utilize brain mosaic variation within specific cell types as distinct indicators for clonal dynamics, denoted as cell-type-specific mosaic variant barcode analysis. From four hemispheres and two different human neurotypical donors, we identified 287 and 780 mosaic variants, respectively, that were used to deconvolve clonal dynamics. Clonal spread and allele fractions within the brain reveal that local hippocampal excitatory neurons are more lineage-restricted than resident neocortical excitatory neurons or resident basal ganglia GABAergic inhibitory neurons. Furthermore, simultaneous genome transcriptome analysis at both a cell-type-specific and a single-cell level suggests a dorsal neocortical origin for a subgroup of DLX1+ inhibitory neurons that disperse radially from an origin shared with excitatory neurons. Finally, the distribution of mosaic variants across 17 locations within one parietal lobe reveals that restriction of clonal spread in the anterior-posterior axis precedes restriction in the dorsal-ventral axis for both excitatory and inhibitory neurons. Thus, cell-type-resolved somatic mosaicism can uncover lineage relationships governing the development of the human forebrain.

Intermuscular adipose tissue (IMAT) is a relatively understudied adipose depot located between muscle fibers. IMAT content increases with age and BMI and is associated with metabolic and muscle degenerative diseases; however, an understanding of the biological properties of IMAT and its interplay with the surrounding muscle fibers is severely lacking. In recent years, single-cell and nuclei RNA sequencing have provided us with cell type-specific atlases of several human tissues. However, the cellular composition of human IMAT remains largely unexplored due to the inherent challenges of its accessibility from biopsy collection in humans. In addition to the limited amount of tissue collected, the processing of human IMAT is complicated due to its proximity to skeletal muscle tissue and fascia. The lipid-laden nature of the adipocytes makes it incompatible with single-cell isolation. Hence, single nuclei RNA sequencing is optimal for obtaining high-dimensional transcriptomics at single-cell resolution and provides the potential to uncover the biology of this depot, including the exact cellular composition of IMAT. Here, we present a detailed protocol for nuclei isolation and library preparation of frozen human IMAT for single nuclei RNA sequencing. This protocol allows for the profiling of thousands of nuclei using a droplet-based approach, thus providing the capacity to detect rare and low-abundant cell types.

In droplet-based single-cell and single-nucleus RNA-seq assays, systematic contamination of ambient RNA molecules biases the quantification of gene expression levels. Existing methods correct the contamination for all genes globally. However, there lacks specific evaluation of correction efficacy for varying contamination levels. Here, we show that DecontX and CellBender under-correct highly contaminating genes, while SoupX and scAR over-correct lowly/non-contaminating genes. Here, we develop scCDC as the first method to detect the contamination-causing genes and only correct expression levels of these genes, some of which are cell-type markers. Compared with existing decontamination methods, scCDC excels in decontaminating highly contaminating genes while avoiding over-correction of other genes.

CD4+ T cells recognize peptide antigens presented on class II major histocompatibility complex (MHC-II) molecules to carry out their function. The remarkable diversity of T cell receptor sequences and lack of antigen discovery approaches for MHC-II make profiling the specificities of CD4+ T cells challenging. We have expanded our platform of signaling and antigen-presenting bifunctional receptors to encode MHC-II molecules presenting covalently linked peptides (SABR-IIs) for CD4+ T cell antigen discovery. SABR-IIs can present epitopes to CD4+ T cells and induce signaling upon their recognition, allowing a readable output. Furthermore, the SABR-II design is modular in signaling and deployment to T cells and B cells. Here, we demonstrate that SABR-IIs libraries presenting endogenous and non-contiguous epitopes can be used for antigen discovery in the context of type 1 diabetes. SABR-II libraries provide a rapid, flexible, scalable and versatile approach for de novo identification of CD4+ T cell ligands from single-cell RNA sequencing data using experimental and computational approaches.

In recent years, data-driven inference of cell-cell communication has helped reveal coordinated biological processes across cell types. Here, we integrate two tools, LIANA and Tensor-cell2cell, which, when combined, can deploy multiple existing methods and resources to enable the robust and flexible identification of cell-cell communication programs across multiple samples. In this work, we show how the integration of our tools facilitates the choice of method to infer cell-cell communication and subsequently perform an unsupervised deconvolution to obtain and summarize biological insights. We explain how to perform the analysis step by step in both Python and R and provide online tutorials with detailed instructions available at https://ccc-protocols.readthedocs.io/. This workflow typically takes ∼1.5 h to complete from installation to downstream visualizations on a graphics processing unit-enabled computer for a dataset of ∼63,000 cells, 10 cell types, and 12 samples.

Dynamically organized chromatin complexes often involve multiplex chromatin interactions and sometimes chromatin-associated RNA1-3. Chromatin complex compositions change during cellular differentiation and ageing, and are expected to be highly heterogeneous among terminally differentiated single cells4-7. Here we introduce the multinucleic acid interaction mapping in single cells (MUSIC) technique for concurrent profiling of multiplex chromatin interactions, gene expression and RNA-chromatin associations within individual nuclei. When applied to 14 human frontal cortex samples from older donors, MUSIC delineated diverse cortical cell types and states. We observed that nuclei exhibiting fewer short-range chromatin interactions were correlated with both an older transcriptomic signature and Alzheimers disease pathology. Furthermore, the cell type exhibiting chromatin contacts between cis expression quantitative trait loci and a promoter tends to be that in which these cis expression quantitative trait loci specifically affect the expression of their target gene. In addition, female cortical cells exhibit highly heterogeneous interactions between XIST non-coding RNA and chromosome X, along with diverse spatial organizations of the X chromosomes. MUSIC presents a potent tool for exploration of chromatin architecture and transcription at cellular resolution in complex tissues.

Human induced pluripotent stem cell (hiPSC) to cardiomyocyte (CM) differentiation has reshaped approaches to studying cardiac development and disease. In this study, we employed a genome-wide CRISPR screen in a hiPSC to CM differentiation system and reveal here that BRD4, a member of the bromodomain and extraterminal (BET) family, regulates CM differentiation. Chemical inhibition of BET proteins in mouse embryonic stem cell (mESC)-derived or hiPSC-derived cardiac progenitor cells (CPCs) results in decreased CM differentiation and persistence of cells expressing progenitor markers. In vivo, BRD4 deletion in second heart field (SHF) CPCs results in embryonic or early postnatal lethality, with mutants demonstrating myocardial hypoplasia and an increase in CPCs. Single-cell transcriptomics identified a subpopulation of SHF CPCs that is sensitive to BRD4 loss and associated with attenuated CM lineage-specific gene programs. These results highlight a previously unrecognized role for BRD4 in CM fate determination during development and a heterogenous requirement for BRD4 among SHF CPCs.

Type 2 diabetes (T2D) is a heterogeneous disease that develops through diverse pathophysiological processes1,2 and molecular mechanisms that are often specific to cell type3,4. Here, to characterize the genetic contribution to these processes across ancestry groups, we aggregate genome-wide association study data from 2,535,601 individuals (39.7% not of European ancestry), including 428,452 cases of T2D. We identify 1,289 independent association signals at genome-wide significance (P

In single-cell RNA sequencing (scRNA-Seq), gene expression is assessed individually for each cell, allowing the investigation of developmental processes, such as embryogenesis and cellular differentiation and regeneration, at unprecedented resolution. In such dynamic biological systems, cellular states form a continuum, e.g., for the differentiation of stem cells into mature cell types. This process is often represented via a trajectory in a reduced-dimensional representation of the scRNA-Seq dataset. While many methods have been suggested for trajectory inference, it is often unclear how to handle multiple biological groups or conditions, e.g., inferring and comparing the differentiation trajectories of wild-type and knock-out stem cell populations. In this manuscript, we present condiments, a method for the inference and downstream interpretation of cell trajectories across multiple conditions. Our framework allows the interpretation of differences between conditions at the trajectory, cell population, and gene expression levels. We start by integrating datasets from multiple conditions into a single trajectory. By comparing the cells conditions along the trajectorys path, we can detect large-scale changes, indicative of differential progression or fate selection. We also demonstrate how to detect subtler changes by finding genes that exhibit different behaviors between these conditions along a differentiation path.

Cerebrospinal fluid (CSF) in the subarachnoid space around the brain has long been known to drain through the lymphatics to cervical lymph nodes1-17, but the connections and regulation have been challenging to identify. Here, using fluorescent CSF tracers in Prox1-GFP lymphatic reporter mice18, we found that the nasopharyngeal lymphatic plexus is a major hub for CSF outflow to deep cervical lymph nodes. This plexus had unusual valves and short lymphangions but no smooth-muscle coverage, whereas downstream deep cervical lymphatics had typical semilunar valves, long lymphangions and smooth muscle coverage that transported CSF to the deep cervical lymph nodes. α-Adrenergic and nitric oxide signalling in the smooth muscle cells regulated CSF drainage through the transport properties of deep cervical lymphatics. During ageing, the nasopharyngeal lymphatic plexus atrophied, but deep cervical lymphatics were not similarly altered, and CSF outflow could still be increased by adrenergic or nitric oxide signalling. Single-cell analysis of gene expression in lymphatic endothelial cells of the nasopharyngeal plexus of aged mice revealed increased type I interferon signalling and other inflammatory cytokines. The importance of evidence for the nasopharyngeal lymphatic plexus functioning as a CSF outflow hub is highlighted by its regression during ageing. Yet, the ageing-resistant pharmacological activation of deep cervical lymphatic transport towards lymph nodes can still increase CSF outflow, offering an approach for augmenting CSF clearance in age-related neurological conditions in which greater efflux would be beneficial.

CRISPR-based gene activation (CRISPRa) is a strategy for upregulating gene expression by targeting promoters or enhancers in a tissue/cell-type specific manner. Here, we describe an experimental framework that combines highly multiplexed perturbations with single-cell RNA sequencing (sc-RNA-seq) to identify cell-type-specific, CRISPRa-responsive cis-regulatory elements and the gene(s) they regulate. Random combinations of many gRNAs are introduced to each of many cells, which are then profiled and partitioned into test and control groups to test for effect(s) of CRISPRa perturbations of both enhancers and promoters on the expression of neighboring genes. Applying this method to a library of 493 gRNAs targeting candidate cis-regulatory elements in both K562 cells and iPSC-derived excitatory neurons, we identify gRNAs capable of specifically upregulating intended target genes and no other neighboring genes within 1 Mb, including gRNAs yielding upregulation of six autism spectrum disorder (ASD) and neurodevelopmental disorder (NDD) risk genes in neurons. A consistent pattern is that the responsiveness of individual enhancers to CRISPRa is restricted by cell type, implying a dependency on either chromatin landscape and/or additional trans-acting factors for successful gene activation. The approach outlined here may facilitate large-scale screens for gRNAs that activate genes in a cell type-specific manner.

BackgroundMany approaches have been developed to overcome technical noise in single cell RNA-sequencing (scRNAseq). As researchers dig deeper into data-looking for rare cell types, subtleties of cell states, and details of gene regulatory networks-there is a growing need for algorithms with controllable accuracy and fewer ad hoc parameters and thresholds. Impeding this goal is the fact that an appropriate null distribution for scRNAseq cannot simply be extracted from data in which ground truth about biological variation is unknown (i.e., usually).ResultsWe approach this problem analytically, assuming that scRNAseq data reflect only cell heterogeneity (what we seek to characterize), transcriptional noise (temporal fluctuations randomly distributed across cells), and sampling error (i.e., Poisson noise). We analyze scRNAseq data without normalization-a step that skews distributions, particularly for sparse data-and calculate p values associated with key statistics. We develop an improved method for selecting features for cell clustering and identifying gene-gene correlations, both positive and negative. Using simulated data, we show that this method, which we call BigSur (Basic Informatics and Gene Statistics from Unnormalized Reads), captures even weak yet significant correlation structures in scRNAseq data. Applying BigSur to data from a clonal human melanoma cell line, we identify thousands of correlations that, when clustered without supervision into gene communities, align with known cellular components and biological processes, and highlight potentially novel cell biological relationships.ConclusionsNew insights into functionally relevant gene regulatory networks can be obtained using a statistically grounded approach to the identification of gene-gene correlations.

Rheumatoid arthritis (RA) management leans toward achieving remission or low disease activity. In this study, we conducted single-cell RNA sequencing (scRNA-Seq) of peripheral blood mononuclear cells (PBMCs) from 36 individuals (18 patients with RA and 18 matched controls, accounting for age, sex, race, and ethnicity), to identify disease-relevant cell subsets and cell type-specific signatures associated with disease activity. Our analysis revealed 18 distinct PBMC subsets, including an IFN-induced transmembrane 3-overexpressing (IFITM3-overexpressing) IFN-activated monocyte subset. We observed an increase in CD4+ T effector memory cells in patients with moderate-high disease activity (DAS28-CRP ≥ 3.2) and a decrease in nonclassical monocytes in patients with low disease activity or remission (DAS28-CRP < 3.2). Pseudobulk analysis by cell type identified 168 differentially expressed genes between RA and matched controls, with a downregulation of proinflammatory genes in the γδ T cell subset, alteration of genes associated with RA predisposition in the IFN-activated subset, and nonclassical monocytes. Additionally, we identified a gene signature associated with moderate-high disease activity, characterized by upregulation of proinflammatory genes such as TNF, JUN, EGR1, IFIT2, MAFB, and G0S2 and downregulation of genes including HLA-DQB1, HLA-DRB5, and TNFSF13B. Notably, cell-cell communication analysis revealed an upregulation of signaling pathways, including VISTA, in both moderate-high and remission-low disease activity contexts. Our findings provide valuable insights into the systemic cellular and molecular mechanisms underlying RA disease activity.

The spatial organization of molecules in a cell is essential for their functions. While current methods focus on discerning tissue architecture, cell-cell interactions, and spatial expression patterns, they are limited to the multicellular scale. We present Bento, a Python toolkit that takes advantage of single-molecule information to enable spatial analysis at the subcellular scale. Bento ingests molecular coordinates and segmentation boundaries to perform three analyses: defining subcellular domains, annotating localization patterns, and quantifying gene-gene colocalization. We demonstrate MERFISH, seqFISH + , Molecular Cartography, and Xenium datasets. Bento is part of the open-source Scverse ecosystem, enabling integration with other single-cell analysis tools.

OBJECTIVES: To determine the dysregulated signaling pathways of head and neck squamous cell carcinoma associated with circulating tumor cells (CTCs) via single-cell molecular characterization. INTRODUCTION: Head and neck squamous cell carcinoma (HNSCC) has a significant global burden and is a disease with poor survival. Despite trials exploring new treatment modalities to improve disease control rates, the 5 year survival rate remains low at only 60%. Most cancer malignancies are reported to progress to a fatal phase due to the metastatic activity derived from treatment-resistant cancer cells, regarded as one of the most significant obstacles to develope effective cancer treatment options. However, the molecular profiles of cancer cells have not been thoroughly studied. METHODS: Here, we examined in-situ HNSCC tumors and pairwisely followed up with the downstream circulating tumor cells (CTCs)-based on the surrogate biomarkers to detect metastasis that is established in other cancers - not yet being fully adopted in HNSCC treatment algorithms. RESULTS: Specifically, we revealed metastatic HNSCC patients have complex CTCs that could be defined through gene expression and mutational gene profiling derived from completed single-cell RNASeq (scRNASeq) that served to confirm molecular pathways inherent in these CTCs. To enhance the reliability of our findings, we cross-validated those molecular profiles with results from previously published studies. CONCLUSION: Thus, we identified 5 dysregulated signaling pathways in CTCs to derive HNSCC biomarker panels for screening HNSCC in situ tumors.

This study focuses on understanding the transcriptional heterogeneity of activated platelets and its impact on diseases such as sepsis, COVID-19, and systemic lupus erythematosus (SLE). Recognizing the limited knowledge in this area, our research aims to dissect the complex transcriptional profiles of activated platelets to aid in developing targeted therapies for abnormal and pathogenic platelet subtypes. We analyzed single-cell transcriptional profiles from 47,977 platelets derived from 413 samples of patients with these diseases, utilizing Deep Neural Network (DNN) and eXtreme Gradient Boosting (XGB) to distinguish transcriptomic signatures predictive of fatal or survival outcomes. Our approach included source data annotations and platelet markers, along with SingleR and Seurat for comprehensive profiling. Additionally, we employed Uniform Manifold Approximation and Projection (UMAP) for effective dimensionality reduction and visualization, aiding in the identification of various platelet subtypes and their relation to disease severity and patient outcomes. Our results highlighted distinct platelet subpopulations that correlate with disease severity, revealing that changes in platelet transcription patterns can intensify endotheliopathy, increasing the risk of coagulation in fatal cases. Moreover, these changes may impact lymphocyte function, indicating a more extensive role for platelets in inflammatory and immune responses. This study identifies crucial biomarkers of platelet heterogeneity in serious health conditions, paving the way for innovative therapeutic approaches targeting platelet activation, which could improve patient outcomes in diseases characterized by altered platelet function.

INTRODUCTION: Hypoxia is a common pathological driver contributing to various forms of pulmonary vascular diseases leading to pulmonary hypertension (PH). Pulmonary interstitial macrophages (IMs) play pivotal roles in immune and vascular dysfunction, leading to inflammation, abnormal remodeling, and fibrosis in PH. However, IMs response to hypoxia and their role in PH progression remain largely unknown. We utilized a murine model of hypoxia-induced PH to investigate the repertoire and functional profiles of IMs in response to acute and prolonged hypoxia, aiming to elucidate their contributions to PH development. METHODS: We conducted single-cell transcriptomic analyses to characterize the repertoire and functional profiles of murine pulmonary IMs following exposure to hypobaric hypoxia for varying durations (0, 1, 3, 7, and 21 days). Hallmark pathways from the mouse Molecular Signatures Database were utilized to characterize the molecular function of the IM subpopulation in response to hypoxia. RESULTS: Our analysis revealed an early acute inflammatory phase during acute hypoxia exposure (Days 1-3), which was resolved by Day 7, followed by a pro-remodeling phase during prolonged hypoxia (Days 7-21). These phases were marked by distinct subpopulations of IMs: MHCIIhiCCR2+EAR2+ cells characterized the acute inflammatory phase, while TLF+VCAM1hi cells dominated the pro-remodeling phase. The acute inflammatory phase exhibited enrichment in interferon-gamma, IL-2, and IL-6 pathways, while the pro-remodeling phase showed dysregulated chemokine production, hemoglobin clearance, and tissue repair profiles, along with activation of distinct complement pathways. DISCUSSION: Our findings demonstrate the existence of distinct populations of pulmonary interstitial macrophages corresponding to acute and prolonged hypoxia exposure, pivotal in regulating the inflammatory and remodeling phases of PH pathogenesis. This understanding offers potential avenues for targeted interventions, tailored to specific populations and distinct phases of the disease. Moreover, further identification of triggers for pro-remodeling IMs holds promise in unveiling novel therapeutic strategies for pulmonary hypertension.

As the circadian clock regulates fundamental biological processes, disrupted clocks are often observed in patients and diseased tissues. Determining the circadian time of the patient or the tissue of focus is essential in circadian medicine and research. Here we present tauFisher, a computational pipeline that accurately predicts circadian time from a single transcriptomic sample by finding correlations between rhythmic genes within the sample. We demonstrate tauFisher's performance in adding timestamps to both bulk and single-cell transcriptomic samples collected from multiple tissue types and experimental settings. Application of tauFisher at a cell-type level in a single-cell RNAseq dataset collected from mouse dermal skin implies that greater circadian phase heterogeneity may explain the dampened rhythm of collective core clock gene expression in dermal immune cells compared to dermal fibroblasts. Given its robustness and generalizability across assay platforms, experimental setups, and tissue types, as well as its potential application in single-cell RNAseq data analysis, tauFisher is a promising tool that facilitates circadian medicine and research.