Search

Fast, nondestructive three-dimensional (3D) imaging of live suspension cells remains challenging without substrate treatment or fixation, precluding scalable single-cell morphometry with minimal alterations. While optical sectioning techniques achieve 3D live cell imaging, lateral versus depth resolution differences further complicate analysis. We present a scalable microfluidic method capable of 3D fluorescent isotropic imaging of live, nonadherent cells suspended inside picoliter droplets with high-speed single-cell volumetric readout (800 to 1,200 slices in 5 to 8 s) and near-diffraction limit resolution (~216 nm). The platform features a droplet trap array that leverages flow-induced droplet interfacial shear to generate intradroplet microvortices, which rotate single cells on their axis to enable optical projection tomography (OPT)-based imaging. This allows gentle (~1 mPa shear stress) observation of cells encapsulated inside nontoxic isotonic buffer droplets, facilitating scalable OPT acquisition by simultaneous spinning of hundreds of cells. We demonstrate 3D imaging of live myeloid and lymphoid cells in suspension, including K562 cells, as well as naive and activated T cells-small cells prone to movement in their suspended phenotype. Our fully suspended, orientation-independent cell morphometry, driven by isotropic imaging and spherical harmonic analysis, enabled the study of primary T cells across various immunological activation states. This approach unveiled six distinct nuclear content distributions, contrasting with conventional 2D images that typically portray spheroid and bean-like nuclear shapes associated with lymphocytes. Our arrayed-droplet OPT technology is capable of isotropic, single live-cell 3D imaging, with the potential to perform large-scale morphometry of immune cell effector function states while providing compatibility with microfluidic droplet operations.

Transcriptional control is fundamental to cellular function. However, despite knowing that transcription factors can repress or activate specific genes, how these functions are implemented at the molecular level has remained elusive, particularly in the endogenous context of developing animals. Here, we combine optogenetics, single-cell live-imaging, and mathematical modeling to study how a zinc-finger repressor, Knirps, induces switch-like transitions into long-lived quiescent states. Using optogenetics, we demonstrate that repression is rapidly reversible (~1 min) and memoryless. Furthermore, we show that the repressor acts by decreasing the frequency of transcriptional bursts in a manner consistent with an equilibrium binding model. Our results provide a quantitative framework for dissecting the in vivo biochemistry of eukaryotic transcriptional regulation.

Increasingly, scRNA-Seq studies explore cell populations across different samples and the effect of sample heterogeneity on organisms phenotype. However, relatively few bioinformatic methods have been developed which adequately address the variation between samples for such population-level analyses. We propose a framework for representing the entire single-cell profile of a sample, which we call a GloScope representation. We implement GloScope on scRNA-Seq datasets from study designs ranging from 12 to over 300 samples and demonstrate how GloScope allows researchers to perform essential bioinformatic tasks at the sample-level, in particular visualization and quality control assessment.

BACKGROUND: Genetic perturbation screens with single-cell readouts have enabled rich phenotyping of gene function and regulatory networks. These approaches have been challenging in vivo, especially in adult disease models such as cancer, which include mixtures of malignant and microenvironment cells. Glioblastoma (GBM) is a fatal cancer, and methods of systematically interrogating gene function and therapeutic targets in vivo, especially in combination with standard of care treatment such as radiotherapy, are lacking. RESULTS: Here, we iteratively develop a multiplex in vivo perturb-seq CRISPRi platform for single-cell genetic screens in cancer and tumor microenvironment cells that leverages intracranial convection enhanced delivery of sgRNA libraries into mouse models of GBM. Our platform enables potent silencing of drivers of in vivo growth and tumor maintenance as well as genes that sensitize GBM to radiotherapy. We find radiotherapy rewires transcriptional responses to genetic perturbations in an in vivo-dependent manner, revealing heterogenous patterns of treatment sensitization or resistance in GBM. Furthermore, we demonstrate targeting of genes that function in the tumor microenvironment, enabling alterations of ligand-receptor interactions between immune and stromal cells following in vivo CRISPRi perturbations that can affect tumor cell phagocytosis. CONCLUSION: In sum, we demonstrate the utility of multiplexed perturb-seq for in vivo single-cell dissection of adult cancer and normal tissue biology across multiple cell types in the context of therapeutic intervention, a platform with potential for broad application.

The drivers of immune evasion are not entirely clear, limiting the success of cancer immunotherapies. Here we applied single-cell spatial and perturbational transcriptomics to delineate immune evasion in high-grade serous tubo-ovarian cancer. To this end, we first mapped the spatial organization of high-grade serous tubo-ovarian cancer by profiling more than 2.5 million cells in situ in 130 tumors from 94 patients. This revealed a malignant cell state that reflects tumor genetics and is predictive of T cell and natural killer cell infiltration levels and response to immune checkpoint blockade. We then performed Perturb-seq screens and identified genetic perturbations-including knockout of PTPN1 and ACTR8-that trigger this malignant cell state. Finally, we show that these perturbations, as well as a PTPN1/PTPN2 inhibitor, sensitize ovarian cancer cells to T cell and natural killer cell cytotoxicity, as predicted. This study thus identifies ways to study and target immune evasion by linking genetic variation, cell-state regulators and spatial biology.

From single-cell RNA-sequencing (scRNA-seq) and spatial transcriptomics (ST), one can extract high-dimensional gene expression patterns that can be described by intercellular communication networks or decoupled gene modules. These two descriptions of information flow are often assumed to occur independently. However, intercellular communication drives directed flows of information that are mediated by intracellular gene modules, in turn triggering outflows of other signals. Methodologies to describe such intercellular flows are lacking. We present FlowSig, a method that infers communication-driven intercellular flows from scRNA-seq or ST data using graphical causal modeling and conditional independence. We benchmark FlowSig using newly generated experimental cortical organoid data and synthetic data generated from mathematical modeling. We demonstrate FlowSig's utility by applying it to various studies, showing that FlowSig can capture stimulation-induced changes to paracrine signaling in pancreatic islets, demonstrate shifts in intercellular flows due to increasing COVID-19 severity and reconstruct morphogen-driven activator-inhibitor patterns in mouse embryogenesis.

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.

Despite recent work linking mixed phenotype acute leukemia (MPAL) to certain genetic lesions, specific driver mutations remain undefined for a significant proportion of patients and no genetic subtype is predictive of clinical outcomes. Moreover, therapeutic strategy for MPAL remains unclear, and prognosis is overall poor. We performed multiomic single cell profiling of 14 newly diagnosed adult MPAL patients to characterize the inter- and intra-tumoral transcriptional, immunophenotypic, and genetic landscapes of MPAL. We show that neither genetic profile nor transcriptome reliably correlate with specific MPAL immunophenotypes. Despite this, we find that MPAL blasts express a shared stem cell-like transcriptional profile indicative of high differentiation potential. Patients with the highest differentiation potential demonstrate inferior survival in our dataset. A gene set score, MPAL95, derived from genes highly enriched in the most stem-like MPAL cells, is applicable to bulk RNA sequencing data and is predictive of survival in an independent patient cohort, suggesting a potential strategy for clinical risk stratification.

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.

Understanding cell fate decision-making during complex biological processes is an open challenge that is now aided by high-resolution single-cell sequencing technologies. Specifically, it remains challenging to identify and characterize transition states corresponding to "tipping points" whereby cells commit to new cell states. Here, we present a computational method that takes advantage of single-cell transcriptomics data to infer the stability and gene regulatory networks (GRNs) along cell lineages. Our method uses the unspliced and spliced counts from single-cell RNA sequencing data and cell ordering along lineage trajectories to train an RNA splicing multivariate model, from which cell-state stability along the lineage is inferred based on spectral analysis of the model's Jacobian matrix. Moreover, the model infers the RNA cross-species interactions resulting in GRNs and their variation along the cell lineage. When applied to epithelial-mesenchymal transition in ovarian and lung cancer-derived cell lines, our model predicts a saddle-node transition between the epithelial and mesenchymal states passing through an unstable, intermediate cell state. Furthermore, we show that the underlying GRN controlling epithelial-mesenchymal transition rearranges during the transition, resulting in denser and less modular networks in the intermediate state. Overall, our method represents a flexible tool to study cell lineages with a combination of theory-driven modeling and single-cell transcriptomics data.

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 (

The molecular organization of the human neocortex historically has been studied in the context of its histological layers. However, emerging spatial transcriptomic technologies have enabled unbiased identification of transcriptionally defined spatial domains that move beyond classic cytoarchitecture. We used the Visium spatial gene expression platform to generate a data-driven molecular neuroanatomical atlas across the anterior-posterior axis of the human dorsolateral prefrontal cortex. Integration with paired single-nucleus RNA-sequencing data revealed distinct cell type compositions and cell-cell interactions across spatial domains. Using PsychENCODE and publicly available data, we mapped the enrichment of cell types and genes associated with neuropsychiatric disorders to discrete spatial domains.

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.

Neural crest cells exemplify cellular diversification from a multipotent progenitor population. However, the full sequence of early molecular choices orchestrating the emergence of neural crest heterogeneity from the embryonic ectoderm remains elusive. Gene-regulatory-networks (GRN) govern early development and cell specification toward definitive neural crest. Here, we combine ultradense single-cell transcriptomes with machine-learning and large-scale transcriptomic and epigenomic experimental validation of selected trajectories, to provide the general principles and highlight specific features of the GRN underlying neural crest fate diversification from induction to early migration stages using Xenopus frog embryos as a model. During gastrulation, a transient neural border zone state precedes the choice between neural crest and placodes which includes multiple converging gene programs. During neurulation, transcription factor connectome, and bifurcation analyses demonstrate the early emergence of neural crest fates at the neural plate stage, alongside an unbiased multipotent-like lineage persisting until epithelial-mesenchymal transition stage. We also decipher circuits driving cranial and vagal neural crest formation and provide a broadly applicable high-throughput validation strategy for investigating single-cell transcriptomes in vertebrate GRNs in development, evolution, and disease.

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