Spatially resolved multiomics of human cardiac niches.

The function of a cell is defined by its intrinsic characteristics and its niche: the tissue microenvironment in which it dwells. Here we combine single-cell and spatial transcriptomics data to discover cellular niches within eight regions of the human heart. We map cells to microanatomical locations and integrate knowledge-based and unsupervised structural annotations. We also profile the cells of the human cardiac conduction system1. The results revealed their distinctive repertoire of ion channels, G-protein-coupled receptors (GPCRs) and regulatory networks, and implicated FOXP2 in the pacemaker phenotype. We show that the sinoatrial node is compartmentalized, with a core of pacemaker cells, fibroblasts and glial cells supporting glutamatergic signalling. Using a custom CellPhoneDB.org module, we identify trans-synaptic pacemaker cell interactions with glia. We introduce a druggable target prediction tool, drug2cell, which leverages single-cell profiles and drug-target interactions to provide mechanistic insights into the chronotropic effects of drugs, including GLP-1 analogues. In the epicardium, we show enrichment of both IgG+ and IgA+ plasma cells forming immune niches that may contribute to infection defence. Overall, we provide new clarity to cardiac electro-anatomy and immunology, and our suite of computational approaches can be applied to other tissues and organs.

Cells of the human intestinal tract mapped across space and time.

The cellular landscape of the human intestinal tract is dynamic throughout life, developing in utero and changing in response to functional requirements and environmental exposures. Here, to comprehensively map cell lineages, we use single-cell RNA sequencing and antigen receptor analysis of almost half a million cells from up to 5 anatomical regions in the developing and up to 11 distinct anatomical regions in the healthy paediatric and adult human gut. This reveals the existence of transcriptionally distinct BEST4 epithelial cells throughout the human intestinal tract. Furthermore, we implicate IgG sensing as a function of intestinal tuft cells. We describe neural cell populations in the developing enteric nervous system, and predict cell-type-specific expression of genes associated with Hirschsprung's disease. Finally, using a systems approach, we identify key cell players that drive the formation of secondary lymphoid tissue in early human development. We show that these programs are adopted in inflammatory bowel disease to recruit and retain immune cells at the site of inflammation. This catalogue of intestinal cells will provide new insights into cellular programs in development, homeostasis and disease.

The mutational landscape of normal human endometrial epithelium.

All normal somatic cells are thought to acquire mutations, but understanding of the rates, patterns, causes and consequences of somatic mutations in normal cells is limited. The uterine endometrium adopts multiple physiological states over a lifetime and is lined by a gland-forming epithelium1,2. Here, using whole-genome sequencing, we show that normal human endometrial glands are clonal cell populations with total mutation burdens that increase at about 29 base substitutions per year and that are many-fold lower than those of endometrial cancers. Normal endometrial glands frequently carry 'driver' mutations in cancer genes, the burden of which increases with age and decreases with parity. Cell clones with drivers often originate during the first decades of life and subsequently progressively colonize the epithelial lining of the endometrium. Our results show that mutational landscapes differ markedly between normal tissues-perhaps shaped by differences in their structure and physiology-and indicate that the procession of neoplastic change that leads to endometrial cancer is initiated early in life.

Cells of the adult human heart.

Cardiovascular disease is the leading cause of death worldwide. Advanced insights into disease mechanisms and therapeutic strategies require a deeper understanding of the molecular processes involved in the healthy heart. Knowledge of the full repertoire of cardiac cells and their gene expression profiles is a fundamental first step in this endeavour. Here, using state-of-the-art analyses of large-scale single-cell and single-nucleus transcriptomes, we characterize six anatomical adult heart regions. Our results highlight the cellular heterogeneity of cardiomyocytes, pericytes and fibroblasts, and reveal distinct atrial and ventricular subsets of cells with diverse developmental origins and specialized properties. We define the complexity of the cardiac vasculature and its changes along the arterio-venous axis. In the immune compartment, we identify cardiac-resident macrophages with inflammatory and protective transcriptional signatures. Furthermore, analyses of cell-to-cell interactions highlight different networks of macrophages, fibroblasts and cardiomyocytes between atria and ventricles that are distinct from those of skeletal muscle. Our human cardiac cell atlas improves our understanding of the human heart and provides a valuable reference for future studies.

Spatial and molecular profiling of the mononuclear phagocyte network in classic Hodgkin lymphoma.

Classic Hodgkin lymphoma (cHL) has a rich immune infiltrate, which is an intrinsic component of the neoplastic process. Malignant Hodgkin Reed-Sternberg cells (HRSCs) create an immunosuppressive microenvironment by the expression of regulatory molecules, preventing T-cell activation. It has also been demonstrated that mononuclear phagocytes (MNPs) in the vicinity of HRSCs express similar regulatory mechanisms in parallel, and their presence in tissue is associated with inferior patient outcomes. MNPs in cHL have hitherto been identified by a small number of canonical markers and are usually described as tumor-associated macrophages. The organization of MNP networks and interactions with HRSCs remains unexplored at high resolution. Here, we defined the global immune-cell composition of cHL and nonlymphoma lymph nodes, integrating data across single-cell RNA sequencing, spatial transcriptomics, and multiplexed immunofluorescence. We observed that MNPs comprise multiple subsets of monocytes, macrophages, and dendritic cells (DCs). Classical monocytes, macrophages and conventional DC2s were enriched in the vicinity of HRSCs, but plasmacytoid DCs and activated DCs were excluded. Unexpectedly, cDCs and monocytes expressed immunoregulatory checkpoints PD-L1, TIM-3, and the tryptophan-catabolizing protein IDO, at the same level as macrophages. Expression of these molecules increased with age. We also found that classical monocytes are important signaling hubs, potentially controlling the retention of cDC2 and ThExh via CCR1-, CCR4-, CCR5-, and CXCR3-dependent signaling. Enrichment of the cDC2-monocyte-macrophage network in diagnostic biopsies is associated with early treatment failure. These results reveal unanticipated complexity and spatial polarization within the MNP compartment, further demonstrating their potential roles in immune evasion by cHL.

Multi-omics and imaging mass cytometry characterization of human kidneys to identify pathways and phenotypes associated with impaired kidney function.

Despite the recent advances in our understanding of the role of lipids, metabolites, and related enzymes in mediating kidney injury, there is limited integrated multi-omics data identifying potential metabolic pathways driving impaired kidney function. The limited availability of kidney biopsies from living donors with acute kidney injury has remained a major constraint. Here, we validated the use of deceased transplant donor kidneys as a good model to study acute kidney injury in humans and characterized these kidneys using imaging and multi-omics approaches. We noted consistent changes in kidney injury and inflammatory markers in donors with reduced kidney function. Neighborhood and correlation analyses of imaging mass cytometry data showed that subsets of kidney cells (proximal tubular cells and fibroblasts) are associated with the expression profile of kidney immune cells, potentially linking these cells to kidney inflammation. Integrated transcriptomic and metabolomic analysis of human kidneys showed that kidney arachidonic acid metabolism and seven other metabolic pathways were upregulated following diminished kidney function. To validate the arachidonic acid pathway in impaired kidney function we demonstrated increased levels of cytosolic phospholipase A2 protein and related lipid mediators (prostaglandin E2) in the injured kidneys. Further, inhibition of cytosolic phospholipase A2 reduced injury and inflammation in human kidney proximal tubular epithelial cells in vitro. Thus, our study identified cell types and metabolic pathways that may be critical for controlling inflammation associated with impaired kidney function in humans.

Unique molecular and functional features of extramedullary hematopoietic stem and progenitor cell reservoirs in humans.

Rare hematopoietic stem and progenitor cell (HSPC) pools outside the bone marrow (BM) contribute to blood production in stress and disease but remain ill-defined. Although nonmobilized peripheral blood (PB) is routinely sampled for clinical management, the diagnosis and monitoring potential of PB HSPCs remain untapped, as no healthy PB HSPC baseline has been reported. Here we comprehensively delineate human extramedullary HSPC compartments comparing spleen, PB, and mobilized PB to BM using single-cell RNA-sequencing and/or functional assays. We uncovered HSPC features shared by extramedullary tissues and others unique to PB. First, in contrast to actively dividing BM HSPCs, we found no evidence of substantial ongoing hematopoiesis in extramedullary tissues at steady state but report increased splenic HSPC proliferative output during stress erythropoiesis. Second, extramedullary hematopoietic stem cells/multipotent progenitors (HSCs/MPPs) from spleen, PB, and mobilized PB share a common transcriptional signature and increased abundance of lineage-primed subsets compared with BM. Third, healthy PB HSPCs display a unique bias toward erythroid-megakaryocytic differentiation. At the HSC/MPP level, this is functionally imparted by a subset of phenotypic CD71+ HSCs/MPPs, exclusively producing erythrocytes and megakaryocytes, highly abundant in PB but rare in other adult tissues. Finally, the unique erythroid-megakaryocytic-skewing of PB is perturbed with age in essential thrombocythemia and ß-thalassemia. Collectively, we identify extramedullary lineage-primed HSPC reservoirs that are nonproliferative in situ and report involvement of splenic HSPCs during demand-adapted hematopoiesis. Our data also establish aberrant composition and function of circulating HSPCs as potential clinical indicators of BM dysfunction.

A proteomic approach towards understanding crypoprotective action of Me2SO on the CHO cell proteome.

Chinese hamster ovary (CHO) cell lines are the most widely used in vitro cells for research and production of recombinant proteins such as rhGH, tPA, and erythropoietin. We aimed to investigate changes in protein profiles after cryopreservation using 2D-DIGE MALDI-TOF MS and network pathway analysis. The proteome changes that occur in CHO cells between freshly prepared cells and cryopreserved cells with and without Me2SO were compared to determine the key proteins and pathways altered during recovery from cryopreservation. A total of 54 proteins were identified and successfully matched to 37 peptide mass fingerprints (PMF). 14 protein spots showed an increase while 23 showed decrease abundance in the Me2SO free group compared to the control. The proteins with increased abundance included vimentin, heat shock protein 60 kDa, mitochondrial, heat shock 70 kDa protein 9, protein disulfide-isomerase A3, voltage-dependent anion-selective channel protein 2. Those with a decrease in abundance were myotubularin, glutathione peroxidase, enolase, phospho glyceromutase, chloride intracellular channel protein 1. The main canonical functional pathway affected involved the unfolded protein response, aldosterone Signaling in Epithelial Cells, 14-3-3-mediated signaling. 2D-DIGE MALDI TOF mass spectrometry and network pathway analysis revealed the differential proteome expression of FreeStyle CHO cells after cryopreservation with and without 5% Me2SOto involve pathways related to post-translational modification, protein folding and cell death and survival (score = 56, 22 focus molecules). This study revealed, for the first time to our knowledge the proteins and their regulated pathways involved in the cryoprotective action of 5% Me2SO. The use of 5% Me2SO as a cryoprotectant maintained the CHO cell proteome in the cryopreserved cells, similar to that of fresh CHO cells.

FlowAtlas: an interactive tool for high-dimensional immunophenotyping analysis bridging FlowJo with computational tools in Julia.

As the dimensionality, throughput and complexity of cytometry data increases, so does the demand for user-friendly, interactive analysis tools that leverage high-performance machine learning frameworks. Here we introduce FlowAtlas: an interactive web application that enables dimensionality reduction of cytometry data without down-sampling and that is compatible with datasets stained with non-identical panels. FlowAtlas bridges the user-friendly environment of FlowJo and computational tools in Julia developed by the scientific machine learning community, eliminating the need for coding and bioinformatics expertise. New population discovery and detection of rare populations in FlowAtlas is intuitive and rapid. We demonstrate the capabilities of FlowAtlas using a human multi-tissue, multi-donor immune cell dataset, highlighting key immunological findings. FlowAtlas is available at https://github.com/gszep/FlowAtlas.jl.git.

Network-based prioritization and validation of regulators of vascular smooth muscle cell proliferation in disease.

Aberrant vascular smooth muscle cell (VSMC) homeostasis and proliferation characterize vascular diseases causing heart attack and stroke. Here we elucidate molecular determinants governing VSMC proliferation by reconstructing gene regulatory networks from single-cell transcriptomics and epigenetic profiling. We detect widespread activation of enhancers at disease-relevant loci in proliferation-predisposed VSMCs. We compared gene regulatory network rewiring between injury-responsive and nonresponsive VSMCs, which suggested shared transcription factors but differing target loci between VSMC states. Through in silico perturbation analysis, we identified and prioritized previously unrecognized regulators of proliferation, including RUNX1 and TIMP1. Moreover, we showed that the pioneer transcription factor RUNX1 increased VSMC responsiveness and that TIMP1 feeds back to promote VSMC proliferation through CD74-mediated STAT3 signaling. Both RUNX1 and the TIMP1-CD74 axis were expressed in human VSMCs, showing low levels in normal arteries and increased expression in disease, suggesting clinical relevance and potential as vascular disease targets.

Network-based prioritization and validation of regulators of vascular smooth muscle cell proliferation in disease.

Aberrant vascular smooth muscle cell (VSMC) homeostasis and proliferation characterize vascular diseases causing heart attack and stroke. Here we elucidate molecular determinants governing VSMC proliferation by reconstructing gene regulatory networks from single-cell transcriptomics and epigenetic profiling. We detect widespread activation of enhancers at disease-relevant loci in proliferation-predisposed VSMCs. We compared gene regulatory network rewiring between injury-responsive and nonresponsive VSMCs, which suggested shared transcription factors but differing target loci between VSMC states. Through in silico perturbation analysis, we identified and prioritized previously unrecognized regulators of proliferation, including RUNX1 and TIMP1. Moreover, we showed that the pioneer transcription factor RUNX1 increased VSMC responsiveness and that TIMP1 feeds back to promote VSMC proliferation through CD74-mediated STAT3 signaling. Both RUNX1 and the TIMP1-CD74 axis were expressed in human VSMCs, showing low levels in normal arteries and increased expression in disease, suggesting clinical relevance and potential as vascular disease targets.

An integrated single-cell reference atlas of the human endometrium.

The complex and dynamic cellular composition of the human endometrium remains poorly understood. Previous endometrial single-cell atlases profiled few donors and lacked consensus in defining cell types. We introduce the Human Endometrial Cell Atlas (HECA), a high-resolution single-cell reference atlas (313,527 cells) combining published and new endometrial single-cell transcriptomics datasets of 63 women with and without endometriosis. HECA assigns consensus and identifies previously unreported cell types, mapped in situ using spatial transcriptomics and validated using a new independent single-nuclei dataset (312,246 nuclei, 63 donors). In the functionalis, we identify intricate stromal-epithelial cell coordination via transforming growth factor beta (TGFß) signaling. In the basalis, we define signaling between fibroblasts and an epithelial population expressing progenitor markers. Integration of HECA with large-scale endometriosis genome-wide association study data pinpoints decidualized stromal cells and macrophages as most likely dysregulated in endometriosis. The HECA is a valuable resource for studying endometrial physiology and disorders, and for guiding microphysiological in vitro systems development.

Human skeletal muscle aging atlas.

Skeletal muscle aging is a key contributor to age-related frailty and sarcopenia with substantial implications for global health. Here we profiled 90,902 single cells and 92,259 single nuclei from 17 donors to map the aging process in the adult human intercostal muscle, identifying cellular changes in each muscle compartment. We found that distinct subsets of muscle stem cells exhibit decreased ribosome biogenesis genes and increased CCL2 expression, causing different aging phenotypes. Our atlas also highlights an expansion of nuclei associated with the neuromuscular junction, which may reflect re-innervation, and outlines how the loss of fast-twitch myofibers is mitigated through regeneration and upregulation of fast-type markers in slow-twitch myofibers with age. Furthermore, we document the function of aging muscle microenvironment in immune cell attraction. Overall, we present a comprehensive human skeletal muscle aging resource ( https://www.muscleageingcellatlas.org/ ) together with an in-house mouse muscle atlas to study common features of muscle aging across species.

Chromatin accessibility profiling of targeted cell populations with laser capture microdissection coupled to ATAC-seq.

Spatially resolved omics technologies reveal context-dependent cellular regulatory networks in tissues of interest. Beyond transcriptome analysis, information on epigenetic traits and chromatin accessibility can provide further insights on gene regulation in health and disease. Nevertheless, compared to the enormous advancements in spatial transcriptomics technologies, the field of spatial epigenomics is much younger and still underexplored. In this study, we report laser capture microdissection coupled to ATAC-seq (LCM-ATAC-seq) applied to fresh frozen samples for the spatial characterization of chromatin accessibility. We first demonstrate the efficient use of LCM coupled to in situ tagmentation and evaluate its performance as a function of cell number, microdissected areas, and tissue type. Further, we demonstrate its use for the targeted chromatin accessibility analysis of discrete contiguous or scattered cell populations in tissues via single-nuclei capture based on immunostaining for specific cellular markers.

A spatially resolved atlas of the human lung characterizes a gland-associated immune niche.

Single-cell transcriptomics has allowed unprecedented resolution of cell types/states in the human lung, but their spatial context is less well defined. To (re)define tissue architecture of lung and airways, we profiled five proximal-to-distal locations of healthy human lungs in depth using multi-omic single cell/nuclei and spatial transcriptomics (queryable at lungcellatlas.org ). Using computational data integration and analysis, we extend beyond the suspension cell paradigm and discover macro and micro-anatomical tissue compartments including previously unannotated cell types in the epithelial, vascular, stromal and nerve bundle micro-environments. We identify and implicate peribronchial fibroblasts in lung disease. Importantly, we discover and validate a survival niche for IgA plasma cells in the airway submucosal glands (SMG). We show that gland epithelial cells recruit B cells and IgA plasma cells, and promote longevity and antibody secretion locally through expression of CCL28, APRIL and IL-6. This new 'gland-associated immune niche' has implications for respiratory health.

Optimal location for fibular osteotomy to provide maximal compression to the tibia in the management of delayed union and hypertrophic non-union of the tibia.

BACKGROUND: Tibial shaft fractures are the commonest long bone fracture, with early weight-bearing improving the rate of bony union. However, an intact fibula can act as a strut that splints the tibial segments and holds them apart. A fibular osteotomy, in which a 2.5 cm length of fibula is removed, has been used to treat delayed and hypertrophic non-union by increasing axial tibial loading. However, there is no consensus on the optimal site for the partial fibulectomy. METHODS: Nine leg specimens were obtained from formalin-embalmed cadavers. Transverse mid-shaft tibial fractures were created using an oscillating saw. A rig was designed to compress the legs with an adjustable axial load and measure the force within the fracture site in order to ascertain load transmission through the tibia over a range of weights. After 2.5cm-long fibulectomies were performed at one of three levels on each specimen, load transmission through the tibia was re-assessed. A beam structure model of the intact leg was designed to explain the findings. RESULTS: With an intact fibula, mean tibial loading at 34 kg was 15.52 ± 3.26 kg, increasing to 17.42 ± 4.13 kg after fibular osteotomy. This increase was only significant where the osteotomy was performed proximal to or at the level of the tibial fracture. Modelling midshaft tibial loading using the Euler-Bernoulli beam theory showed that 80.5% of the original force was transmitted through the tibia with an intact fibula, rising to 81.1% after a distal fibulectomy, and 100% with a proximal fibulectomy. CONCLUSION: This study describes a novel method of measuring axial tibial forces. We demonstrated that a fibular osteotomy increases axial tibial loading regardless of location, with the greatest increase after proximal fibular osteotomy. A contributing factor for this can be explained by a simple beam model. We therefore recommend a proximal fibular osteotomy when it is performed in the treatment of delayed and non-union of tibial midshaft fractures.

Management of post appendicectomy intra-abdominal collections: A volumetric cut off for drainage?

BACKGROUND: Intra-abdominal collections (IAC) are a common complication following appendicectomy, one of the most commonly performed emergency abdominal procedures in childhood. The option to drain a collection is frequently available but not always required. AIM: The aim of this study was to compare the outcomes of medically and procedurally-managed post appendicectomy IACs and suggest a method of standardising the need for intervention. METHODS: A single centre, retrospective review of children aged ≤ 16 years presenting between 2014 and 2019 was performed. Patient demographics, management, and outcome data were collected. IAC volume and surface area were calculated assuming a prolate spheroid or true ellipsoid depending on the number of dimensions reported. RESULTS: 60 patients (18%) of 334 patients developed an IAC post appendicectomy. Medical management was undertaken in 44 (73%), drainage in 12 (20%), and surgical washout in 4 (7%). Collection size was associated with failure of medical management: maximum diameter (p = 0.028), volume (p = 0.002), and surface area (p = 0.001). Collections with a volume of 2 ml/kg were significantly less likely to fail medical management than larger collections (0/33 vs 6/11; p < 0.0001). DISCUSSION: Not all post appendicectomy IACs require drainage. The relationship between collection volume and need for drainage is more closely assessed using a volume calculation rather than a single dimension measurement, particularly when adjusted for weight of the child. A cut off of 2 ml/kg appears to be a good objective measure for intervention and provides a communication tool for discussion amongst the multidisciplinary team. Prospectively collected multicentre data on this subject would be timely. LEVEL OF EVIDENCE: III.

Two subsets of human marginal zone B cells resolved by global analysis of lymphoid tissues and blood.

B cells generate antibodies that are essential for immune protection, but their subgroups are poorly defined. Here, we perform undirected deep profiling of B cells in matched human lymphoid tissues from deceased transplant organ donors and blood. In addition to identifying unanticipated features of tissue-based B cell differentiation, we resolve two subsets of marginal zone B (MZB) cells differing in cell surface and transcriptomic profiles, clonal relationships to other subsets, enrichment of genes in the NOTCH pathway, distribution bias within splenic marginal zone microenvironment, and immunoglobulin repertoire diversity and hypermutation frequency. Each subset is present in spleen, gut-associated lymphoid tissue, mesenteric lymph nodes, and blood. MZB cells and the lineage from which they are derived are depleted in lupus nephritis. Here, we show that this depletion is of only one MZB subset. The other remains unchanged as a proportion of total B cells compared with health. Thus, it is important to factor MZB cell heterogeneity into studies of human B cell responses and pathology.

Group 3 innate lymphocytes make a distinct contribution to type 17 immunity in bladder defence.

Bladder infection affects a hundred million people annually, but our understanding of bladder immunity is incomplete. We found type 17 immune response genes among the most up-regulated networks in mouse bladder following uropathogenic Escherichia coli (UPEC) challenge. Intravital imaging revealed submucosal Rorc+ cells responsive to UPEC challenge, and we found increased Il17 and IL22 transcripts in wild-type and Rag2 -/- mice, implicating group 3 innate lymphoid cells (ILC3s) as a source of these cytokines. NCR-positive and negative ILC3 subsets were identified in murine and human bladders, with local proliferation increasing IL17-producing ILC3s post infection. ILC3s made a more limited contribution to bladder IL22, with prominent early induction of IL22 evident in Th17 cells. Single-cell RNA sequencing revealed bladder NCR-negative ILC3s as the source of IL17 and identified putative ILC3-myeloid cell interactions, including via lymphotoxin-ß-LTBR. Altogether, our data provide important insights into the orchestration and execution of type 17 immunity in bladder defense.

Single-cell atlas of human liver development reveals pathways directing hepatic cell fates.

The liver has been studied extensively due to the broad number of diseases affecting its vital functions. However, therapeutic advances have been hampered by the lack of knowledge concerning human hepatic development. Here, we addressed this limitation by describing the developmental trajectories of different cell types that make up the human liver at single-cell resolution. These transcriptomic analyses revealed that sequential cell-to-cell interactions direct functional maturation of hepatocytes, with non-parenchymal cells playing essential roles during organogenesis. We utilized this information to derive bipotential hepatoblast organoids and then exploited this model system to validate the importance of signalling pathways in hepatocyte and cholangiocyte specification. Further insights into hepatic maturation also enabled the identification of stage-specific transcription factors to improve the functionality of hepatocyte-like cells generated from human pluripotent stem cells. Thus, our study establishes a platform to investigate the basic mechanisms directing human liver development and to produce cell types for clinical applications.