Transcriptomic signatures of progressive and regressive liver fibrosis and portal hypertension.

Persistent liver injury triggers a fibrogenic program that causes pathologic remodeling of the hepatic microenvironment (i.e., liver fibrosis) and portal hypertension. The dynamics of gene regulation during liver disease progression and early regression remain understudied. Here, we generated hepatic transcriptome profiles in two well-established liver disease models at peak fibrosis and during spontaneous regression after the removal of the inducing agents. We linked the dynamics of key disease readouts, such as portal pressure, collagen area, and transaminase levels, to differentially expressed genes, enabling the identification of transcriptomic signatures of progressive vs. regressive liver fibrosis and portal hypertension. These candidate biomarkers (e.g., Tcf4, Mmp7, Trem2, Spp1, Scube1, Islr) were validated in RNA sequencing datasets of patients with cirrhosis and portal hypertension, and those cured from hepatitis C infection. Finally, deconvolution identified major cell types and suggested an association of macrophage and portal hepatocyte signatures with portal hypertension and fibrosis area.

A high-resolution view of the heterogeneous aging endothelium.

Vascular endothelial cell (EC) aging has a strong impact on tissue perfusion and overall cardiovascular health. While studies confined to the investigation of aging-associated vascular readouts in one or a few tissues have already drastically expanded our understanding of EC aging, single-cell omics and other high-resolution profiling technologies have started to illuminate the intricate molecular changes underlying endothelial aging across diverse tissues and vascular beds at scale. In this review, we provide an overview of recent insights into the heterogeneous adaptations of the aging vascular endothelium. We address critical questions regarding tissue-specific and universal responses of the endothelium to the aging process, EC turnover dynamics throughout lifespan, and the differential susceptibility of ECs to acquiring aging-associated traits. In doing so, we underscore the transformative potential of single-cell approaches in advancing our comprehension of endothelial aging, essential to foster the development of future innovative therapeutic strategies for aging-associated vascular conditions.

The evolving functions of the vasculature in regulating adipose tissue biology in health and obesity.

Adipose tissue is an endocrine organ and a crucial regulator of energy storage and systemic metabolic homeostasis. Additionally, adipose tissue is a pivotal regulator of cardiovascular health and disease, mediated in part by the endocrine and paracrine secretion of several bioactive products, such as adipokines. Adipose vasculature has an instrumental role in the modulation of adipose tissue expansion, homeostasis and metabolism. The role of the adipose vasculature has been extensively explored in the context of obesity, which is recognized as a global health problem. Obesity-induced accumulation of fat, in combination with vascular rarefaction, promotes adipocyte dysfunction and induces oxidative stress, hypoxia and inflammation. It is now recognized that obesity-associated endothelial dysfunction often precedes the development of cardiovascular diseases. Investigations have revealed heterogeneity within the vascular niche and dynamic reciprocity between vascular and adipose cells, which can become dysregulated in obesity. Here we provide a comprehensive overview of the evolving functions of the vasculature in regulating adipose tissue biology in health and obesity.

Prioritization and functional validation of target genes from single-cell transcriptomics studies.

Translation of academic results into clinical practice is a formidable unmet medical need. Single-cell RNA-sequencing (scRNA-seq) studies generate long descriptive ranks of markers with predicted biological function, but without functional validation, it remains challenging to know which markers truly exert the putative function. Given the lengthy/costly nature of validation studies, gene prioritization is required to select candidates. We address these issues by studying tip endothelial cell (EC) marker genes because of their importance for angiogenesis. Here, by tailoring Guidelines On Target Assessment for Innovative Therapeutics, we in silico prioritize previously unreported/poorly described, high-ranking tip EC markers. Notably, functional validation reveals that four of six candidates behave as tip EC genes. We even discover a tip EC function for a gene lacking in-depth functional annotation. Thus, validating prioritized genes from scRNA-seq studies offers opportunities for identifying targets to be considered for possible translation, but not all top-ranked scRNA-seq markers exert the predicted function.

Deciphering endothelial heterogeneity in health and disease at single-cell resolution: progress and perspectives.

Endothelial cells (ECs) constitute the inner lining of vascular beds in mammals and are crucial for homeostatic regulation of blood vessel physiology, but also play a key role in pathogenesis of many diseases, thereby representing realistic therapeutic targets. However, it has become evident that ECs are heterogeneous, encompassing several subtypes with distinct functions, which makes EC targeting and modulation in diseases challenging. The rise of the new single-cell era has led to an emergence of studies aimed at interrogating transcriptome diversity along the vascular tree, and has revolutionized our understanding of EC heterogeneity from both a physiological and pathophysiological context. Here, we discuss recent landmark studies aimed at teasing apart the heterogeneous nature of ECs. We cover driving (epi)genetic, transcriptomic, and metabolic forces underlying EC heterogeneity in health and disease, as well as current strategies used to combat disease-enriched EC phenotypes, and propose strategies to transcend largely descriptive heterogeneity towards prioritization and functional validation of therapeutically targetable drivers of EC diversity. Lastly, we provide an overview of the most recent advances and hurdles in single EC OMICs.

The pulmonary vasculature in lethal COVID-19 and idiopathic pulmonary fibrosis at single-cell resolution.

AIMS: Severe acute respiratory syndrome coronavirus-2 infection causes COVID-19, which in severe cases evokes life-threatening acute respiratory distress syndrome (ARDS). Transcriptome signatures and the functional relevance of non-vascular cell types (e.g. immune and epithelial cells) in COVID-19 are becoming increasingly evident. However, despite its known contribution to vascular inflammation, recruitment/invasion of immune cells, vascular leakage, and perturbed haemostasis in the lungs of severe COVID-19 patients, an in-depth interrogation of the endothelial cell (EC) compartment in lethal COVID-19 is lacking. Moreover, progressive fibrotic lung disease represents one of the complications of COVID-19 pneumonia and ARDS. Analogous features between idiopathic pulmonary fibrosis (IPF) and COVID-19 suggest partial similarities in their pathophysiology, yet, a head-to-head comparison of pulmonary cell transcriptomes between both conditions has not been implemented to date. METHODS AND RESULTS: We performed single-nucleus RNA-sequencing on frozen lungs from 7 deceased COVID-19 patients, 6 IPF explant lungs, and 12 controls. The vascular fraction, comprising 38 794 nuclei, could be subclustered into 14 distinct EC subtypes. Non-vascular cell types, comprising 137 746 nuclei, were subclustered and used for EC-interactome analyses. Pulmonary ECs of deceased COVID-19 patients showed an enrichment of genes involved in cellular stress, as well as signatures suggestive of dampened immunomodulation and impaired vessel wall integrity. In addition, increased abundance of a population of systemic capillary and venous ECs was identified in COVID-19 and IPF. COVID-19 systemic ECs closely resembled their IPF counterparts, and a set of 30 genes was found congruently enriched in systemic ECs across studies. Receptor-ligand interaction analysis of ECs with non-vascular cell types in the pulmonary micro-environment revealed numerous previously unknown interactions specifically enriched/depleted in COVID-19 and/or IPF. CONCLUSIONS: This study uncovered novel insights into the abundance, expression patterns, and interactomes of EC subtypes in COVID-19 and IPF, relevant for future investigations into the progression and treatment of both lethal conditions.

Neuroepithelial organoid patterning is mediated by a neighborhood watch mechanism.

During epithelial tissue patterning, morphogens operate across multiple length scales to instruct cell identities. However, how cell fate changes are coordinated over these scales to establish spatial organization remains poorly understood. Here, we use human neural tube organoids as models of epithelial patterning and develop an in silico approach to define conditions permissive to patterning. By systematically varying morphogen position, diffusivity, and fate-inducing concentration levels, we show that cells follow a "neighborhood watch" (NW) mechanism that is deterministically dictated by initial morphogen source positions, reflecting scale-invariant in vitro phenotypes. We define how the frequency and local bias of morphogen sources stabilize pattern orientation. The model predicts enhanced patterning through floor plate inhibition, and receptor-ligand interaction analysis of single-cell RNA sequencing (scRNA-seq) data identifies wingless-related integration site (WNT) and bone morphogenic protein (BMP) as inhibition modulators, which we validate in vitro. These results suggest that developing neuroepithelia employ NW-based mechanisms to organize morphogen sources, define cellular identity, and establish patterns.

Single cell atlas identifies lipid-processing and immunomodulatory endothelial cells in healthy and malignant breast.

Since a detailed inventory of endothelial cell (EC) heterogeneity in breast cancer (BC) is lacking, here we perform single cell RNA-sequencing of 26,515 cells (including 8433 ECs) from 9 BC patients and compare them to published EC taxonomies from lung tumors. Angiogenic ECs are phenotypically similar, while other EC subtypes are different. Predictive interactome analysis reveals known but also previously unreported receptor-ligand interactions between ECs and immune cells, suggesting an involvement of breast EC subtypes in immune responses. We also identify a capillary EC subtype (LIPEC (Lipid Processing EC)), which expresses genes involved in lipid processing that are regulated by PPAR-γ and is more abundant in peri-tumoral breast tissue. Retrospective analysis of 4648 BC patients reveals that treatment with metformin (an indirect PPAR-γ signaling activator) provides long-lasting clinical benefit and is positively associated with LIPEC abundance. Our findings warrant further exploration of this LIPEC/PPAR-γ link for BC treatment.

Generation of vessel co-option lung metastases mouse models for single-cell isolation of metastases-derived cells and endothelial cells.

Tumor vessel co-option, a process in which cancer cells "hijack" pre-existing blood vessels to grow and invade healthy tissue, is poorly understood but is a proposed resistance mechanism against anti-angiogenic therapy (AAT). Here, we describe protocols for establishing murine renal (RENCA) and breast (4T1) cancer lung vessel co-option metastases models. Moreover, we outline a reproducible protocol for single-cell isolation from murine lung metastases using magnetic-activated cell sorting as well as immunohistochemical stainings to distinguish vessel co-option from angiogenesis. For complete details on the use and execution of this protocol, please refer to Teuwen et al. (2021).

Endothelial cell heterogeneity and microglia regulons revealed by a pig cell landscape at single-cell level.

Pigs are valuable large animal models for biomedical and genetic research, but insights into the tissue- and cell-type-specific transcriptome and heterogeneity remain limited. By leveraging single-cell RNA sequencing, we generate a multiple-organ single-cell transcriptomic map containing over 200,000 pig cells from 20 tissues/organs. We comprehensively characterize the heterogeneity of cells in tissues and identify 234 cell clusters, representing 58 major cell types. In-depth integrative analysis of endothelial cells reveals a high degree of heterogeneity. We identify several functionally distinct endothelial cell phenotypes, including an endothelial to mesenchymal transition subtype in adipose tissues. Intercellular communication analysis predicts tissue- and cell type-specific crosstalk between endothelial cells and other cell types through the VEGF, PDGF, TGF-ß, and BMP pathways. Regulon analysis of single-cell transcriptome of microglia in pig and 12 other species further identifies MEF2C as an evolutionally conserved regulon in the microglia. Our work describes the landscape of single-cell transcriptomes within diverse pig organs and identifies the heterogeneity of endothelial cells and evolutionally conserved regulon in microglia.

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.

Protocols for endothelial cell isolation from mouse tissues: brain, choroid, lung, and muscle.

Endothelial cells (ECs) harbor distinct phenotypical and functional characteristics depending on their tissue localization and contribute to brain, eye, lung, and muscle diseases such as dementia, macular degeneration, pulmonary hypertension, and sarcopenia. To study their function, isolation of pure ECs in high quantities is crucial. Here, we describe protocols for rapid and reproducible blood vessel EC purification established for scRNA sequencing from murine tissues using mechanical and enzymatic digestion followed by magnetic and fluorescence-activated cell sorting. For complete details on the use and execution of these protocol, please refer to Kalucka et al. (2020), Rohlenova et al. (2020), and Goveia et al. (2020).

Protocols for endothelial cell isolation from mouse tissues: kidney, spleen, and testis.

Endothelial cells (ECs) exhibit phenotypic and functional tissue specificities, critical for studies in the vascular field and beyond. Thus, tissue-specific methods for isolation of highly purified ECs are necessary. Kidney, spleen, and testis ECs are relevant players in health and diseases such as chronic kidney disease, acute kidney injury, myelofibrosis, and cancer. Here, we provide tailored protocols for rapid and reproducible EC purification established for scRNA sequencing from these adult murine tissues using the combination of magnetic- and fluorescence-activated cell sorting. For complete details on the use and execution of these protocols, please refer to Kalucka et al. (2020) and Dumas et al. (2020).

Arhgef2 regulates mitotic spindle orientation in hematopoietic stem cells and is essential for productive hematopoiesis.

How hematopoietic stem cells (HSCs) coordinate their divisional axis and whether this orientation is important for stem cell-driven hematopoiesis is poorly understood. Single-cell RNA sequencing data from patients with Shwachman-Diamond syndrome (SDS), an inherited bone marrow failure syndrome, show that ARHGEF2, a RhoA-specific guanine nucleotide exchange factor and determinant of mitotic spindle orientation, is specifically downregulated in SDS hematopoietic stem and progenitor cells (HSPCs). We demonstrate that transplanted Arhgef2-/- fetal liver and bone marrow cells yield impaired hematopoietic recovery and a production deficit from long-term HSCs, phenotypes that are not the result of differences in numbers of transplanted HSCs, their cell cycle status, level of apoptosis, progenitor output, or homing ability. Notably, these defects are functionally restored in vivo by overexpression of ARHGEF2 or its downstream activated RHOA GTPase. By using live imaging of dividing HSPCs, we show an increased frequency of misoriented divisions in the absence of Arhgef2. ARHGEF2 knockdown in human HSCs also impairs their ability to regenerate hematopoiesis, culminating in significantly smaller xenografts. Together, these data demonstrate a conserved role for Arhgef2 in orienting HSPC division and suggest that HSCs may divide in certain orientations to establish hematopoiesis, the loss of which could contribute to HSC dysfunction in bone marrow failure.

Tumor vessel co-option probed by single-cell analysis.

Tumor vessel co-option is poorly understood, yet it is a resistance mechanism against anti-angiogenic therapy (AAT). The heterogeneity of co-opted endothelial cells (ECs) and pericytes, co-opting cancer and myeloid cells in tumors growing via vessel co-option, has not been investigated at the single-cell level. Here, we use a murine AAT-resistant lung tumor model, in which VEGF-targeting induces vessel co-option for continued growth. Single-cell RNA sequencing (scRNA-seq) of 31,964 cells reveals, unexpectedly, a largely similar transcriptome of co-opted tumor ECs (TECs) and pericytes as their healthy counterparts. Notably, we identify cell types that might contribute to vessel co-option, i.e., an invasive cancer-cell subtype, possibly assisted by a matrix-remodeling macrophage population, and another M1-like macrophage subtype, possibly involved in keeping or rendering vascular cells quiescent.

Protocols for endothelial cell isolation from mouse tissues: small intestine, colon, heart, and liver.

Endothelial cells (ECs) from the small intestine, colon, liver, and heart have distinct phenotypes and functional adaptations that are dependent on their physiological environment. Gut ECs adapt to low oxygen, heart ECs to contractile forces, and liver ECs to low flow rates. Isolating high-purity ECs in sufficient quantities is crucial to study their functions. Here, we describe protocols combining magnetic and fluorescent activated cell sorting for rapid and reproducible EC purification from four adult murine tissues. For complete details on the use and execution of these protocols, please refer to Kalucka et al. (2020).

Single-Cell RNA Sequencing Maps Endothelial Metabolic Plasticity in Pathological Angiogenesis.

Endothelial cell (EC) metabolism is an emerging target for anti-angiogenic therapy in tumor angiogenesis and choroidal neovascularization (CNV), but little is known about individual EC metabolic transcriptomes. By single-cell RNA sequencing 28,337 murine choroidal ECs (CECs) and sprouting CNV-ECs, we constructed a taxonomy to characterize their heterogeneity. Comparison with murine lung tumor ECs (TECs) revealed congruent marker gene expression by distinct EC phenotypes across tissues and diseases, suggesting similar angiogenic mechanisms. Trajectory inference predicted that differentiation of venous to angiogenic ECs was accompanied by metabolic transcriptome plasticity. ECs displayed metabolic transcriptome heterogeneity during cell-cycle progression and in quiescence. Hypothesizing that conserved genes are important, we used an integrated analysis, based on congruent transcriptome analysis, CEC-tailored genome-scale metabolic modeling, and gene expression meta-analysis in cross-species datasets, followed by in vitro and in vivo validation, to identify SQLE and ALDH18A1 as previously unknown metabolic angiogenic targets.