Stromal Cells Covering Omental Fat-Associated Lymphoid Clusters Trigger Formation of Neutrophil Aggregates to Capture Peritoneal Contaminants.

The omentum is a visceral adipose tissue rich in fat-associated lymphoid clusters (FALCs) that collects peritoneal contaminants and provides a first layer of immunological defense within the abdomen. Here, we investigated the mechanisms that mediate the capture of peritoneal contaminants during peritonitis. Single-cell RNA sequencing and spatial analysis of omental stromal cells revealed that the surface of FALCs were covered by CXCL1+ mesothelial cells, which we termed FALC cover cells. Blockade of CXCL1 inhibited the recruitment and aggregation of neutrophils at FALCs during zymosan-induced peritonitis. Inhibition of protein arginine deiminase 4, an enzyme important for the release of neutrophil extracellular traps, abolished neutrophil aggregation and the capture of peritoneal contaminants by omental FALCs. Analysis of omental samples from patients with acute appendicitis confirmed neutrophil recruitment and bacterial capture at FALCs. Thus, specialized omental mesothelial cells coordinate the recruitment and aggregation of neutrophils to capture peritoneal contaminants.

A Macrophage-Pericyte Axis Directs Tissue Restoration via Amphiregulin-Induced Transforming Growth Factor Beta Activation.

The epidermal growth factor receptor ligand Amphiregulin has a well-documented role in the restoration of tissue homeostasis after injury; however, the mechanism by which Amphiregulin contributes to wound repair remains unknown. Here we show that Amphiregulin functioned by releasing bioactive transforming growth factor beta (TGF-ß) from latent complexes via integrin-αV activation. Using acute injury models in two different tissues, we found that by inducing TGF-ß activation on mesenchymal stromal cells (pericytes), Amphiregulin induced their differentiation into myofibroblasts, thereby selectively contributing to the restoration of vascular barrier function within injured tissue. Furthermore, we identified macrophages as a critical source of Amphiregulin, revealing a direct effector mechanism by which these cells contribute to tissue restoration after acute injury. Combined, these observations expose a so far under-appreciated mechanism of how cells of the immune system selectively control the differentiation of tissue progenitor cells during tissue repair and inflammation.

Opposing roles of hepatic stellate cell subpopulations in hepatocarcinogenesis.

Hepatocellular carcinoma (HCC), the fourth leading cause of cancer mortality worldwide, develops almost exclusively in patients with chronic liver disease and advanced fibrosis1,2. Here we interrogated functions of hepatic stellate cells (HSCs), the main source of liver fibroblasts3, during hepatocarcinogenesis. Genetic depletion, activation or inhibition of HSCs in mouse models of HCC revealed their overall tumour-promoting role. HSCs were enriched in the preneoplastic environment, where they closely interacted with hepatocytes and modulated hepatocarcinogenesis by regulating hepatocyte proliferation and death. Analyses of mouse and human HSC subpopulations by single-cell RNA sequencing together with genetic ablation of subpopulation-enriched mediators revealed dual functions of HSCs in hepatocarcinogenesis. Hepatocyte growth factor, enriched in quiescent and cytokine-producing HSCs, protected against hepatocyte death and HCC development. By contrast, type I collagen, enriched in activated myofibroblastic HSCs, promoted proliferation and tumour development through increased stiffness and TAZ activation in pretumoural hepatocytes and through activation of discoidin domain receptor 1 in established tumours. An increased HSC imbalance between cytokine-producing HSCs and myofibroblastic HSCs during liver disease progression was associated with increased HCC risk in patients. In summary, the dynamic shift in HSC subpopulations and their mediators during chronic liver disease is associated with a switch from HCC protection to HCC promotion.

Decoding myofibroblast origins in human kidney fibrosis.

Kidney fibrosis is the hallmark of chronic kidney disease progression; however, at present no antifibrotic therapies exist1-3. The origin, functional heterogeneity and regulation of scar-forming cells that occur during human kidney fibrosis remain poorly understood1,2,4. Here, using single-cell RNA sequencing, we profiled the transcriptomes of cells from the proximal and non-proximal tubules of healthy and fibrotic human kidneys to map the entire human kidney. This analysis enabled us to map all matrix-producing cells at high resolution, and to identify distinct subpopulations of pericytes and fibroblasts as the main cellular sources of scar-forming myofibroblasts during human kidney fibrosis. We used genetic fate-tracing, time-course single-cell RNA sequencing and ATAC-seq (assay for transposase-accessible chromatin using sequencing) experiments in mice, and spatial transcriptomics in human kidney fibrosis, to shed light on the cellular origins and differentiation of human kidney myofibroblasts and their precursors at high resolution. Finally, we used this strategy to detect potential therapeutic targets, and identified NKD2 as a myofibroblast-specific target in human kidney fibrosis.

Pax6 limits the competence of developing cerebral cortical cells to respond to inductive intercellular signals.

The development of stable specialized cell types in multicellular organisms relies on mechanisms controlling inductive intercellular signals and the competence of cells to respond to such signals. In developing cerebral cortex, progenitors generate only glutamatergic excitatory neurons despite being exposed to signals with the potential to initiate the production of other neuronal types, suggesting that their competence is limited. Here, we tested the hypothesis that this limitation is due to their expression of transcription factor Pax6. We used bulk and single-cell RNAseq to show that conditional cortex-specific Pax6 deletion from the onset of cortical neurogenesis allowed some progenitors to generate abnormal lineages resembling those normally found outside the cortex. Analysis of selected gene expression showed that the changes occurred in specific spatiotemporal patterns. We then compared the responses of control and Pax6-deleted cortical cells to in vivo and in vitro manipulations of extracellular signals. We found that Pax6 loss increased cortical progenitors' competence to generate inappropriate lineages in response to extracellular factors normally present in developing cortex, including the morphogens Shh and Bmp4. Regional variation in the levels of these factors could explain spatiotemporal patterns of fate change following Pax6 deletion in vivo. We propose that Pax6's main role in developing cortical cells is to minimize the risk of their development being derailed by the potential side effects of morphogens engaged contemporaneously in other essential functions.

Cell atlas of the Atlantic salmon spleen reveals immune cell heterogeneity and cell-specific responses to bacterial infection.

The spleen is a conserved secondary lymphoid organ that emerged in parallel to adaptive immunity in early jawed vertebrates. Recent studies have applied single cell transcriptomics to reveal the cellular composition of spleen in several species, cataloguing diverse immune cell types and subpopulations. In this study, 51,119 spleen nuclei transcriptomes were comprehensively investigated in the commercially important teleost Atlantic salmon (Salmo salar L.), contrasting control animals with those challenged with the bacterial pathogen Aeromonas salmonicida. We identified clusters of nuclei representing the expected major cell types, namely T cells, B cells, natural killer-like cells, granulocytes, mononuclear phagocytes, endothelial cells, mesenchymal cells, erythrocytes and thrombocytes. We discovered heterogeneity within several immune lineages, providing evidence for resident macrophages and melanomacrophages, infiltrating monocytes, several candidate dendritic cell subpopulations, and B cells at distinct stages of differentiation, including plasma cells and an igt + subset. We provide evidence for twelve candidate T cell subsets, including cd4+ T helper and regulatory T cells, one cd8+ subset, three γδT subsets, and populations double negative for cd4 and cd8. The number of genes showing differential expression during the early stages of Aeromonas infection was highly variable across immune cell types, with the largest changes observed in macrophages and infiltrating monocytes, followed by resting mature B cells. Our analysis provides evidence for a local inflammatory response to infection alongside B cell maturation in the spleen, and upregulation of ccr9 genes in igt + B cells, T helper and cd8+ cells, and monocytes, consistent with the recruitment of immune cell populations to the gut to deal with Aeromonas infection. Overall, this study provides a new cell-resolved perspective of the immune actions of Atlantic salmon spleen, highlighting extensive heterogeneity hidden to bulk transcriptomics. We further provide a large catalogue of cell-specific marker genes that can be leveraged to further explore the function and structural organization of the salmonid immune system.

Mapping the cellular landscape of Atlantic salmon head kidney by single cell and single nucleus transcriptomics.

Single-cell transcriptomics is the current gold standard for global gene expression profiling, not only in mammals and model species, but also in non-model fish species. This is a rapidly expanding field, creating a deeper understanding of tissue heterogeneity and the distinct functions of individual cells, making it possible to explore the complexities of immunology and gene expression on a highly resolved level. In this study, we compared two single cell transcriptomic approaches to investigate cellular heterogeneity within the head kidney of healthy farmed Atlantic salmon (Salmo salar). We compared 14,149 cell transcriptomes assayed by single cell RNA-seq (scRNA-seq) with 18,067 nuclei transcriptomes captured by single nucleus RNA-Seq (snRNA-seq). Both approaches detected eight major cell populations in common: granulocytes, heamatopoietic stem cells, erythrocytes, mononuclear phagocytes, thrombocytes, B cells, NK-like cells, and T cells. Four additional cell types, endothelial, epithelial, interrenal, and mesenchymal cells, were detected in the snRNA-seq dataset, but appeared to be lost during preparation of the single cell suspension submitted for scRNA-seq library generation. We identified additional heterogeneity and subpopulations within the B cells, T cells, and endothelial cells, and revealed developmental trajectories of heamatopoietic stem cells into differentiated granulocyte and mononuclear phagocyte populations. Gene expression profiles of B cell subtypes revealed distinct IgM and IgT-skewed resting B cell lineages and provided insights into the regulation of B cell lymphopoiesis. The analysis revealed eleven T cell sub-populations, displaying a level of T cell heterogeneity in salmon head kidney comparable to that observed in mammals, including distinct subsets of cd4/cd8-negative T cells, such as tcrγ positive, progenitor-like, and cytotoxic cells. Although snRNA-seq and scRNA-seq were both useful to resolve cell type-specific expression in the Atlantic salmon head kidney, the snRNA-seq pipeline was overall more robust in identifying several cell types and subpopulations. While scRNA-seq displayed higher levels of ribosomal and mitochondrial genes, snRNA-seq captured more transcription factor genes. However, only scRNA-seq-generated data was useful for cell trajectory inference within the myeloid lineage. In conclusion, this study systematically outlines the relative merits of scRNA-seq and snRNA-seq in Atlantic salmon, enhances understanding of teleost immune cell lineages, and provides a comprehensive list of markers for identifying major cell populations in the head kidney with significant immune relevance.

Single-cell RNA sequencing redefines the mesenchymal cell landscape of mouse endometrium.

The endometrium is a dynamic tissue that exhibits remarkable resilience to repeated episodes of differentiation, breakdown, regeneration, and remodeling. Endometrial physiology relies on a complex interplay between the stromal and epithelial compartments with the former containing a mixture of fibroblasts, vascular, and immune cells. There is evidence for rare populations of putative mesenchymal progenitor cells located in the perivascular niche of human endometrium, but the existence of an equivalent cell population in mouse is unclear. We used the Pdgfrb-BAC-eGFP transgenic reporter mouse in combination with bulk and single-cell RNA sequencing to redefine the endometrial mesenchyme. In contrast to previous reports we show that CD146 is expressed in both PDGFRß + perivascular cells and CD31 + endothelial cells. Bulk RNAseq revealed cells in the perivascular niche which express the high levels of Pdgfrb as well as genes previously identified in pericytes and/or vascular smooth muscle cells (Acta2, Myh11, Olfr78, Cspg4, Rgs4, Rgs5, Kcnj8, and Abcc9). scRNA-seq identified five subpopulations of cells including closely related pericytes/vascular smooth muscle cells and three subpopulations of fibroblasts. All three fibroblast populations were PDGFRα+/CD34 + but were distinct in their expression of Ngfr/Spon2/Angptl7 (F1), Cxcl14/Smoc2/Rgs2 (F2), and Clec3b/Col14a1/Mmp3 (F3), with potential functions in the regulation of immune responses, response to wounding, and organization of extracellular matrix, respectively. Immunohistochemistry was used to investigate the spatial distribution of these populations revealing F1/NGFR + cells in most abundance beside epithelial cells. We provide the first definitive analysis of mesenchymal cells in the adult mouse endometrium identifying five subpopulations providing a platform for comparisons between mesenchymal cells in endometrium and other adult tissues which are prone to fibrosis.

Transcriptional dynamics of pluripotent stem cell-derived endothelial cell differentiation revealed by single-cell RNA sequencing.

AIMS: Pluripotent stem cell-derived endothelial cell products possess therapeutic potential in ischaemic vascular disease. However, the factors that drive endothelial differentiation from pluripotency and cellular specification are largely unknown. The aims of this study were to use single-cell RNA sequencing (scRNA-seq) to map the transcriptional landscape and cellular dynamics of directed differentiation of human embryonic stem cell-derived endothelial cells (hESC-EC) and to compare these cells to mature endothelial cells from diverse vascular beds. METHODS AND RESULTS: A highly efficient directed 8-day differentiation protocol was used to generate a hESC-derived endothelial cell product (hESC-ECP), in which 66% of cells co-expressed CD31 and CD144. We observed largely homogeneous hESC and mesodermal populations at Days 0 and 4, respectively, followed by a rapid emergence of distinct endothelial and mesenchymal populations. Pseudotime trajectory identified transcriptional signatures of endothelial commitment and maturation during the differentiation process. Concordance in transcriptional signatures was verified by scRNA-seq analysis using both a second hESC line RC11, and an alternative hESC-EC differentiation protocol. In total, 105 727 cells were subjected to scRNA-seq analysis. Global transcriptional comparison revealed a transcriptional architecture of hESC-EC that differs from freshly isolated and cultured human endothelial cells and from organ-specific endothelial cells. CONCLUSION: A transcriptional bifurcation into endothelial and mesenchymal lineages was identified, as well as novel transcriptional signatures underpinning commitment and maturation. The transcriptional architecture of hESC-ECP was distinct from mature and foetal human EC.

Single-cell transcriptome analyses reveal novel targets modulating cardiac neovascularization by resident endothelial cells following myocardial infarction.

AIMS: A better understanding of the pathways that regulate regeneration of the coronary vasculature is of fundamental importance for the advancement of strategies to treat patients with heart disease. Here, we aimed to investigate the origin and clonal dynamics of endothelial cells (ECs) associated with neovascularization in the adult mouse heart following myocardial infarction (MI). Furthermore, we sought to define murine cardiac endothelial heterogeneity and to characterize the transcriptional profiles of pro-angiogenic resident ECs in the adult mouse heart, at single-cell resolution. METHODS AND RESULTS: An EC-specific multispectral lineage-tracing mouse (Pdgfb-iCreERT2-R26R-Brainbow2.1) was used to demonstrate that structural integrity of adult cardiac endothelium following MI was maintained through clonal proliferation by resident ECs in the infarct border region, without significant contributions from bone marrow cells or endothelial-to-mesenchymal transition. Ten transcriptionally discrete heterogeneous EC states, as well as the pathways through which each endothelial state is likely to enhance neovasculogenesis and tissue regeneration following ischaemic injury were defined. Plasmalemma vesicle-associated protein (Plvap) was selected for further study, which showed an endothelial-specific and increased expression in both the ischaemic mouse and human heart, and played a direct role in regulating human endothelial proliferation in vitro. CONCLUSION: We present a single-cell gene expression atlas of cardiac specific resident ECs, and the transcriptional hierarchy underpinning endogenous vascular repair following MI. These data provide a rich resource that could assist in the development of new therapeutic interventions to augment endogenous myocardial perfusion and enhance regeneration in the injured heart.

Increased linear bone growth by GH in the absence of SOCS2 is independent of IGF-1.

Growth hormone (GH) signaling is essential for postnatal linear bone growth, but the relative importance of GHs actions on the liver and/or growth plate cartilage remains unclear. The importance of liver derived insulin like-growth factor-1 (IGF-1) for endochondral growth has recently been challenged. Here, we investigate linear growth in Suppressor of Cytokine Signaling-2 (SOCS2) knockout mice, which have enhanced growth despite normal systemic GH/IGF-1 levels. Wild-type embryonic ex vivo metatarsals failed to exhibit increased linear growth in response to GH, but displayed increased Socs2 transcript levels (P < 0.01). In the absence of SOCS2, GH treatment enhanced metatarsal linear growth over a 12 day period. Despite this increase, IGF-1 transcript and protein levels were not increased in response to GH. In accordance with these data, IGF-1 levels were unchanged in GH-challenged postnatal Socs2(-/-) conditioned medium despite metatarsals showing enhanced linear growth. Growth-plate Igf1 mRNA levels were not elevated in juvenile Socs2(-/-) mice. GH did however elevate IGF-binding protein 3 levels in conditioned medium from GH challenged metatarsals and this was more apparent in Socs2(-/-) metatarsals. GH did not enhance the growth of Socs2(-/-) metatarsals when the IGF receptor was inhibited, suggesting that IGF receptor mediated mechanisms are required. IGF-2 may be responsible as IGF-2 promoted metatarsal growth and Igf2 expression was elevated in Socs2(-/-) (but not WT) metatarsals in response to GH. These studies emphasise the critical importance of SOCS2 in regulating GHs ability to promote bone growth. Also, GH appears to act directly on the metatarsals of Socs2(-/-) mice, promoting growth via a mechanism that is independent of IGF-1.

A versatile nuclei extraction protocol for single nucleus sequencing in non-model species-Optimization in various Atlantic salmon tissues.

The use of single cell sequencing technologies has exploded over recent years, and is now commonly used in many non-model species. Sequencing nuclei instead of whole cells has become increasingly popular, as it does not require the processing of samples immediately after collection. Here we present a highly effective nucleus isolation protocol that outperforms previously available method in challenging samples in a non-model specie. This protocol can be successfully applied to extract nuclei from a variety of tissues and species.

An autocrine signaling circuit in hepatic stellate cells underlies advanced fibrosis in nonalcoholic steatohepatitis.

Advanced hepatic fibrosis, driven by the activation of hepatic stellate cells (HSCs), affects millions worldwide and is the strongest predictor of mortality in nonalcoholic steatohepatitis (NASH); however, there are no approved antifibrotic therapies. To identify antifibrotic drug targets, we integrated progressive transcriptomic and morphological responses that accompany HSC activation in advanced disease using single-nucleus RNA sequencing and tissue clearing in a robust murine NASH model. In advanced fibrosis, we found that an autocrine HSC signaling circuit emerged that was composed of 68 receptor-ligand interactions conserved between murine and human NASH. These predicted interactions were supported by the parallel appearance of markedly increased direct stellate cell-cell contacts in murine NASH. As proof of principle, pharmacological inhibition of one such autocrine interaction, neurotrophic receptor tyrosine kinase 3-neurotrophin 3, inhibited human HSC activation in culture and reversed advanced murine NASH fibrosis. In summary, we uncovered a repertoire of antifibrotic drug targets underlying advanced fibrosis in vivo. The findings suggest a therapeutic paradigm in which stage-specific therapies could yield enhanced antifibrotic efficacy in patients with advanced hepatic fibrosis.

Human and murine fibroblast single-cell transcriptomics reveals fibroblast clusters are differentially affected by ageing and serum cholesterol.

AIMS: Specific fibroblast markers and in-depth heterogeneity analysis are currently lacking, hindering functional studies in cardiovascular diseases (CVDs). Here, we established cell-type markers and heterogeneity in murine and human arteries and studied the adventitial fibroblast response to CVD and its risk factors hypercholesterolaemia and ageing. METHODS AND RESULTS: Murine aorta single-cell RNA-sequencing analysis of adventitial mesenchymal cells identified fibroblast-specific markers. Immunohistochemistry and flow cytometry validated platelet-derived growth factor receptor alpha (PDGFRA) and dipeptidase 1 (DPEP1) across human and murine aorta, carotid, and femoral arteries, whereas traditional markers such as the cluster of differentiation (CD)90 and vimentin also marked transgelin+ vascular smooth muscle cells. Next, pseudotime analysis showed multiple fibroblast clusters differentiating along trajectories. Three trajectories, marked by CD55 (Cd55+), Cxcl chemokine 14 (Cxcl14+), and lysyl oxidase (Lox+), were reproduced in an independent RNA-seq dataset. Gene ontology (GO) analysis showed divergent functional profiles of the three trajectories, related to vascular development, antigen presentation, and/or collagen fibril organization, respectively. Trajectory-specific genes included significantly more genes with known genome-wide associations (GWAS) to CVD than expected by chance, implying a role in CVD. Indeed, differential regulation of fibroblast clusters by CVD risk factors was shown in the adventitia of aged C57BL/6J mice, and mildly hypercholesterolaemic LDLR KO mice on chow by flow cytometry. The expansion of collagen-related CXCL14+ and LOX+ fibroblasts in aged and hypercholesterolaemic aortic adventitia, respectively, coincided with increased adventitial collagen. Immunohistochemistry, bulk, and single-cell transcriptomics of human carotid and aorta specimens emphasized translational value as CD55+, CXCL14+ and LOX+ fibroblasts were observed in healthy and atherosclerotic specimens. Also, trajectory-specific gene sets are differentially correlated with human atherosclerotic plaque traits. CONCLUSION: We provide two adventitial fibroblast-specific markers, PDGFRA and DPEP1, and demonstrate fibroblast heterogeneity in health and CVD in humans and mice. Biological relevance is evident from the regulation of fibroblast clusters by age and hypercholesterolaemia in vivo, associations with human atherosclerotic plaque traits, and enrichment of genes with a GWAS for CVD.

Single-cell RNA sequencing and lineage tracing confirm mesenchyme to epithelial transformation (MET) contributes to repair of the endometrium at menstruation.

The human endometrium experiences repetitive cycles of tissue wounding characterised by piecemeal shedding of the surface epithelium and rapid restoration of tissue homeostasis. In this study, we used a mouse model of endometrial repair and three transgenic lines of mice to investigate whether epithelial cells that become incorporated into the newly formed luminal epithelium have their origins in one or more of the mesenchymal cell types present in the stromal compartment of the endometrium. Using scRNAseq, we identified a novel population of PDGFRb + mesenchymal stromal cells that developed a unique transcriptomic signature in response to endometrial breakdown/repair. These cells expressed genes usually considered specific to epithelial cells and in silico trajectory analysis suggested they were stromal fibroblasts in transition to becoming epithelial cells. To confirm our hypothesis we used a lineage tracing strategy to compare the fate of stromal fibroblasts (PDGFRa+) and stromal perivascular cells (NG2/CSPG4+). We demonstrated that stromal fibroblasts can undergo a mesenchyme to epithelial transformation and become incorporated into the re-epithelialised luminal surface of the repaired tissue. This study is the first to discover a novel population of wound-responsive, plastic endometrial stromal fibroblasts that contribute to the rapid restoration of an intact luminal epithelium during endometrial repair. These findings form a platform for comparisons both to endometrial pathologies which involve a fibrotic response (Asherman's syndrome, endometriosis) as well as other mucosal tissues which have a variable response to wounding.

The human uterus is a formidable organ. From puberty to menopause, it completely sheds off its internal lining every 28 days or so, creating what is in effect a large open wound. Unlike the skin or other parts of the body, however, this tissue can quickly repair itself without scarring. This fascinating process remains poorly understood, partly because human samples and animal models that mimic human menstruation are still lacking. This makes it difficult to grasp how various types of uterine cells get mobilised for healing. To fill this gap, Kirkwood et al. focused on fibroblasts, a heterogenous cell population which helps to support the epithelial cells lining the inside of the uterus. How these cells responded to the advent of menstruation was examined in female mice genetically manipulated to have human-like periods. A method known as single-cell RNAseq was used to track which genes were active in each of these cells before, one day and two days after period onset. This revealed the existence of a subpopulation of cells which only appeared when wound healing was most needed. These 'repair-specific' fibroblasts expressed a mixture of genes; those typical of fibroblasts but also some known to be active in the epithelial cells lining the uterus. This suggests that the cells were in the process of changing their identity so they could remake the uterine layer lost during a period. And indeed, labelling these fibroblasts with a fluorescent tag showed that, during healing, they had migrated from within the uterine tissue to become part of its newly restored internal surface. These results represent the first evidence that fibroblasts play a direct role in repairing the uterus during menstruation. From endometriosis to infertility, the lives of millions of people around the world are impacted by disorders which affect the uterine lining. A better understanding of how the uterus can fix itself month after month could help to find new treatments for these conditions. This knowledge could also be useful for to address abnormal wound healing in the skin and other tissues, as this process often involves fibroblasts.

Mapping the developing human cardiac endothelium at single-cell resolution identifies MECOM as a regulator of arteriovenous gene expression.

AIMS: Coronary vasculature formation is a critical event during cardiac development, essential for heart function throughout perinatal and adult life. However, current understanding of coronary vascular development has largely been derived from transgenic mouse models. The aim of this study was to characterize the transcriptome of the human foetal cardiac endothelium using single-cell RNA sequencing (scRNA-seq) to provide critical new insights into the cellular heterogeneity and transcriptional dynamics that underpin endothelial specification within the vasculature of the developing heart. METHODS AND RESULTS: We acquired scRNA-seq data of over 10 000 foetal cardiac endothelial cells (ECs), revealing divergent EC subtypes including endocardial, capillary, venous, arterial, and lymphatic populations. Gene regulatory network analyses predicted roles for SMAD1 and MECOM in determining the identity of capillary and arterial populations, respectively. Trajectory inference analysis suggested an endocardial contribution to the coronary vasculature and subsequent arterialization of capillary endothelium accompanied by increasing MECOM expression. Comparative analysis of equivalent data from murine cardiac development demonstrated that transcriptional signatures defining endothelial subpopulations are largely conserved between human and mouse. Comprehensive characterization of the transcriptional response to MECOM knockdown in human embryonic stem cell-derived EC (hESC-EC) demonstrated an increase in the expression of non-arterial markers, including those enriched in venous EC. CONCLUSIONS: scRNA-seq of the human foetal cardiac endothelium identified distinct EC populations. A predicted endocardial contribution to the developing coronary vasculature was identified, as well as subsequent arterial specification of capillary EC. Loss of MECOM in hESC-EC increased expression of non-arterial markers, suggesting a role in maintaining arterial EC identity.

Single cell transcriptomics of Atlantic salmon (Salmo salar L.) liver reveals cellular heterogeneity and immunological responses to challenge by Aeromonas salmonicida.

The liver is a multitasking organ with essential functions for vertebrate health spanning metabolism and immunity. In contrast to mammals, our understanding of liver cellular heterogeneity and its role in regulating immunological status remains poorly defined in fishes. Addressing this knowledge gap, we generated a transcriptomic atlas of 47,432 nuclei isolated from the liver of Atlantic salmon (Salmo salar L.) contrasting control fish with those challenged with a pathogenic strain of Aeromonas salmonicida, a problematic bacterial pathogen in global aquaculture. We identified the major liver cell types and their sub-populations, revealing poor conservation of many hepatic cell marker genes utilized in mammals, while identifying novel heterogeneity within the hepatocyte, lymphoid, and myeloid lineages. This included polyploid hepatocytes, multiple T cell populations including γδ T cells, and candidate populations of monocytes/macrophages and dendritic cells. A dominant hepatocyte population radically remodeled its transcriptome following infection to activate the acute phase response and other defense functions, while repressing routine functions such as metabolism. These defense-specialized hepatocytes showed strong activation of genes controlling protein synthesis and secretion, presumably to support the release of acute phase proteins into circulation. The infection response further involved up-regulation of numerous genes in an immune-cell specific manner, reflecting functions in pathogen recognition and killing, antigen presentation, phagocytosis, regulation of inflammation, B cell differentiation and T cell activation. Overall, this study greatly enhances our understanding of the multifaceted role played by liver immune and non-immune cells in host defense and metabolic remodeling following infection and provides many novel cell-specific marker genes to empower future studies of this organ in fishes.

Deciphering Mesenchymal Drivers of Human Dupuytren's Disease at Single-Cell Level.

Dupuytren's disease (DD) is a common, progressive fibroproliferative disease affecting the palmar fascia of the hands, causing fingers to irreversibly flex toward the palm with significant loss of function. Surgical treatments are limited; therefore, effective new therapies for DD are urgently required. To identify the key cellular and molecular pathways driving DD, we employed single-cell RNA sequencing, profiling the transcriptomes of 35,250 human single cells from DD, nonpathogenic fascia, and healthy dermis. We identify a DD-specific population of pathogenic PDPN+/FAP+ mesenchymal cells displaying an elevated expression of fibrillar collagens and profibrogenic genes. In silico trajectory analysis reveals resident fibroblasts to be the source of this pathogenic population. To resolve the processes governing DD progression, genes differentially expressed during fibroblast differentiation were identified, including upregulated TNFRSF12A and transcription factor SCX. Knockdown of SCX and blockade of TNFRSF12A inhibited the proliferation and altered the profibrotic gene expression of cultured human FAP+ mesenchymal cells, demonstrating a functional role for these genes in DD. The power of single-cell RNA sequencing is utilized to identify the major pathogenic mesenchymal subpopulations driving DD and the key molecular pathways regulating the DD-specific myofibroblast phenotype. Using this precision medicine approach, inhibition of TNFRSF12A has shown potential clinical utility in the treatment of DD.

Single-cell RNA sequencing profiling of mouse endothelial cells in response to pulmonary arterial hypertension.

AIMS: Endothelial cell (EC) dysfunction drives the initiation and pathogenesis of pulmonary arterial hypertension (PAH). We aimed to characterize EC dynamics in PAH at single-cell resolution. METHODS AND RESULTS: We carried out single-cell RNA sequencing (scRNA-seq) of lung ECs isolated from an EC lineage-tracing mouse model in Control and SU5416/hypoxia-induced PAH conditions. EC populations corresponding to distinct lung vessel types, including two discrete capillary populations, were identified in both Control and PAH mice. Differential gene expression analysis revealed global PAH-induced EC changes that were confirmed by bulk RNA-seq. This included upregulation of the major histocompatibility complex class II pathway, supporting a role for ECs in the inflammatory response in PAH. We also identified a PAH response specific to the second capillary EC population including upregulation of genes involved in cell death, cell motility, and angiogenesis. Interestingly, four genes with genetic variants associated with PAH were dysregulated in mouse ECs in PAH. To compare relevance across PAH models and species, we performed a detailed analysis of EC heterogeneity and response to PAH in rats and humans through whole-lung PAH scRNA-seq datasets, revealing that 51% of up-regulated mouse genes were also up-regulated in rat or human PAH. We identified promising new candidates to target endothelial dysfunction including CD74, the knockdown of which regulates EC proliferation and barrier integrity in vitro. Finally, with an in silico cell ordering approach, we identified zonation-dependent changes across the arteriovenous axis in mouse PAH and showed upregulation of the Serine/threonine-protein kinase Sgk1 at the junction between the macro- and microvasculature. CONCLUSION: This study uncovers PAH-induced EC transcriptomic changes at a high resolution, revealing novel targets for potential therapeutic candidate development.

Selective Myeloid Depletion of Galectin-3 Offers Protection Against Acute and Chronic Lung Injury.

Rationale: Galectin-3 (Gal-3) is an immune regulator and an important driver of fibrosis in chronic lung injury, however, its role in acute lung injury (ALI) remains unknown. Previous work has shown that global deletion of galectin-3 reduces collagen deposition in a bleomycin-induced pulmonary fibrosis model (MacKinnon et al., Am. J. Respir. Crit. Care Med., 2012, 185, 537-46). An inhaled Gal-3 inhibitor, GB0139, is undergoing Phase II clinical development for idiopathic pulmonary fibrosis (IPF). This work aims to elucidate the role of Gal-3 in the myeloid and mesenchymal compartment on the development of acute and chronic lung injury. Methods: LgalS3 fl/fl mice were generated and crossed with mice expressing the myeloid (LysM) and mesenchymal (Pdgfrb) cre drivers to yield LysM-cre +/- /LgalS3 fl/fl and Pdgfrb-cre +/- /LgalS3 fl/fl mice. The response to acute (bleomycin or LPS) or chronic (bleomycin) lung injury was compared to globally deficient Gal-3 -/- mice. Results: Myeloid depletion of Gal-3 led to a significant reduction in Gal-3 expression in alveolar macrophages and neutrophils and a reduction in neutrophil recruitment into the interstitium but not into the alveolar space. The reduction in interstitial neutrophils corelated with decreased levels of pulmonary inflammation following acute bleomycin and LPS administration. In addition, myeloid deletion decreased Gal-3 levels in bronchoalveolar lavage (BAL) and reduced lung fibrosis induced by chronic bleomycin. In contrast, no differences in BAL Gal-3 levels or fibrosis were observed in Pdgfrb-cre +/- /LgalS3 fl/fl mice. Conclusions: Myeloid cell derived Galectin-3 drives acute and chronic lung inflammation and supports direct targeting of galectin-3 as an attractive new therapy for lung inflammation.