Single-cell transcriptomic analysis of skeletal muscle regeneration across mouse lifespan identifies altered stem cell states associated with senescence.

Skeletal muscle regeneration is driven by the interaction of myogenic and non-myogenic cells. In aging, regeneration is impaired due to dysfunctions of myogenic and non-myogenic cells, but this is not understood comprehensively. We collected an integrated atlas of 273,923 single-cell transcriptomes from muscles of young, old, and geriatric mice (~5, 20, 26 months-old) at six time-points following myotoxin injury. We identified eight cell types, including T and NK cells and macrophage subtypes, that displayed accelerated or delayed response dynamics between ages. Through pseudotime analysis, we observed myogenic cell states and trajectories specific to old and geriatric ages. To explain these age differences, we assessed cellular senescence by scoring experimentally derived and curated gene-lists. This pointed to an elevation of senescent-like subsets specifically within the self-renewing muscle stem cells in aged muscles. This resource provides a holistic portrait of the altered cellular states underlying skeletal muscle regenerative decline across mouse lifespan.

Combinatorial treatment rescues tumour-microenvironment-mediated attenuation of MALT1 inhibitors in B-cell lymphomas.

Activated B-cell-like diffuse large B-cell lymphomas (ABC-DLBCLs) are characterized by constitutive activation of nuclear factor κB driven by the B-cell receptor (BCR) and Toll-like receptor (TLR) pathways. However, BCR-pathway-targeted therapies have limited impact on DLBCLs. Here we used >1,100 DLBCL patient samples to determine immune and extracellular matrix cues in the lymphoid tumour microenvironment (Ly-TME) and built representative synthetic-hydrogel-based B-cell-lymphoma organoids accordingly. We demonstrate that Ly-TME cellular and biophysical factors amplify the BCR-MYD88-TLR9 multiprotein supercomplex and induce cooperative signalling pathways in ABC-DLBCL cells, which reduce the efficacy of compounds targeting the BCR pathway members Bruton tyrosine kinase and mucosa-associated lymphoid tissue lymphoma translocation protein 1 (MALT1). Combinatorial inhibition of multiple aberrant signalling pathways induced higher antitumour efficacy in lymphoid organoids and implanted ABC-DLBCL patient tumours in vivo. Our studies define the complex crosstalk between malignant ABC-DLBCL cells and Ly-TME, and provide rational combinatorial therapies that rescue Ly-TME-mediated attenuation of treatment response to MALT1 inhibitors.

Large-scale integration of single-cell transcriptomic data captures transitional progenitor states in mouse skeletal muscle regeneration.

Skeletal muscle repair is driven by the coordinated self-renewal and fusion of myogenic stem and progenitor cells. Single-cell gene expression analyses of myogenesis have been hampered by the poor sampling of rare and transient cell states that are critical for muscle repair, and do not inform the spatial context that is important for myogenic differentiation. Here, we demonstrate how large-scale integration of single-cell and spatial transcriptomic data can overcome these limitations. We created a single-cell transcriptomic dataset of mouse skeletal muscle by integration, consensus annotation, and analysis of 23 newly collected scRNAseq datasets and 88 publicly available single-cell (scRNAseq) and single-nucleus (snRNAseq) RNA-sequencing datasets. The resulting dataset includes more than 365,000 cells and spans a wide range of ages, injury, and repair conditions. Together, these data enabled identification of the predominant cell types in skeletal muscle, and resolved cell subtypes, including endothelial subtypes distinguished by vessel-type of origin, fibro-adipogenic progenitors defined by functional roles, and many distinct immune populations. The representation of different experimental conditions and the depth of transcriptome coverage enabled robust profiling of sparsely expressed genes. We built a densely sampled transcriptomic model of myogenesis, from stem cell quiescence to myofiber maturation, and identified rare, transitional states of progenitor commitment and fusion that are poorly represented in individual datasets. We performed spatial RNA sequencing of mouse muscle at three time points after injury and used the integrated dataset as a reference to achieve a high-resolution, local deconvolution of cell subtypes. We also used the integrated dataset to explore ligand-receptor co-expression patterns and identify dynamic cell-cell interactions in muscle injury response. We provide a public web tool to enable interactive exploration and visualization of the data. Our work supports the utility of large-scale integration of single-cell transcriptomic data as a tool for biological discovery.

High-resolution single-cell transcriptomics reveals heterogeneity of self-renewing hair follicle stem cells.

Multipotent bulge stem cells (SCs) fuel the hair follicle (HF) cyclic growth during adult skin homeostasis, but their intrinsic molecular heterogeneity is not well understood. These hair follicle stem cells (HFSCs) engage in bouts of self-renewal, migration and differentiation during the hair cycle. Here, we perform high-resolution single-cell RNA sequencing (scRNA-seq) of HFSCs sorted as CD34+ /K14-H2BGFP+ from mouse skin at mid-anagen, the self-renewal stage. We dissect the transcriptomic profiles and unravel population-specific transcriptional heterogeneity. Unsupervised clustering reveals five major HFSC populations, which distinguished by known markers associated with both the bulge and the outer root sheath (ORS) underneath. These populations include quiescent bulge, ORS cellular states and proliferative cells. Lineage trajectory analysis predicted the prospective differentiation path of these cellular states and their corresponding self-renewing subpopulations. The bulge population itself can be further sub-divided into distinct subpopulations that can be mapped to the upper, mid and lower bulge regions, and present a decreasing quiescence score. Gene set enrichment analysis (GSEA) revealed new markers and suggested potentially distinct functions of the ORS and bulge subpopulations. This included communications between the upper bulge subpopulation and sensory nerves and between the upper ORS and skin vasculature, as well as enrichment of a bulge subset in cell migratory functions. The lower ORS enriched genes may potentially enable nutrients passing from the surrounding fat and vasculature cells towards the proliferating hair matrix cells. Thus, we provide a comprehensive account of HFSC molecular heterogeneity during their self-renewing stage, which enables future HF functional studies.