The gluconeogenesis enzyme PCK2 has a non-enzymatic role in proteostasis in endothelial cells.

Endothelial cells (ECs) are highly glycolytic, but whether they generate glycolytic intermediates via gluconeogenesis (GNG) in glucose-deprived conditions remains unknown. Here, we report that glucose-deprived ECs upregulate the GNG enzyme PCK2 and rely on a PCK2-dependent truncated GNG, whereby lactate and glutamine are used for the synthesis of lower glycolytic intermediates that enter the serine and glycerophospholipid biosynthesis pathways, which can play key roles in redox homeostasis and phospholipid synthesis, respectively. Unexpectedly, however, even in normal glucose conditions, and independent of its enzymatic activity, PCK2 silencing perturbs proteostasis, beyond its traditional GNG role. Indeed, PCK2-silenced ECs have an impaired unfolded protein response, leading to accumulation of misfolded proteins, which due to defective proteasomes and impaired autophagy, results in the accumulation of protein aggregates in lysosomes and EC demise. Ultimately, loss of PCK2 in ECs impaired vessel sprouting. This study identifies a role for PCK2 in proteostasis beyond GNG.

VEGF-B prevents excessive angiogenesis by inhibiting FGF2/FGFR1 pathway.

Although VEGF-B was discovered as a VEGF-A homolog a long time ago, the angiogenic effect of VEGF-B remains poorly understood with limited and diverse findings from different groups. Notwithstanding, drugs that inhibit VEGF-B together with other VEGF family members are being used to treat patients with various neovascular diseases. It is therefore critical to have a better understanding of the angiogenic effect of VEGF-B and the underlying mechanisms. Using comprehensive in vitro and in vivo methods and models, we reveal here for the first time an unexpected and surprising function of VEGF-B as an endogenous inhibitor of angiogenesis by inhibiting the FGF2/FGFR1 pathway when the latter is abundantly expressed. Mechanistically, we unveil that VEGF-B binds to FGFR1, induces FGFR1/VEGFR1 complex formation, and suppresses FGF2-induced Erk activation, and inhibits FGF2-driven angiogenesis and tumor growth. Our work uncovers a previously unrecognized novel function of VEGF-B in tethering the FGF2/FGFR1 pathway. Given the anti-angiogenic nature of VEGF-B under conditions of high FGF2/FGFR1 levels, caution is warranted when modulating VEGF-B activity to treat neovascular diseases.

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.

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.

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.

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).

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).

c-Rel orchestrates energy-dependent epithelial and macrophage reprogramming in fibrosis.

Fibrosis is a common pathological feature of chronic disease. Deletion of the NF-κB subunit c-Rel limits fibrosis in multiple organs, although the mechanistic nature of this protection is unresolved. Using cell-specific gene-targeting manipulations in mice undergoing liver damage, we elucidate a critical role for c-Rel in controlling metabolic changes required for inflammatory and fibrogenic activities of hepatocytes and macrophages and identify Pfkfb3 as the key downstream metabolic mediator of this response. Independent deletions of Rel in hepatocytes or macrophages suppressed liver fibrosis induced by carbon tetrachloride, while combined deletion had an additive anti-fibrogenic effect. In transforming growth factor-ß1-induced hepatocytes, c-Rel regulates expression of a pro-fibrogenic secretome comprising inflammatory molecules and connective tissue growth factor, the latter promoting collagen secretion from HMs. Macrophages lacking c-Rel fail to polarize to M1 or M2 states, explaining reduced fibrosis in RelΔLysM mice. Pharmacological inhibition of c-Rel attenuated multi-organ fibrosis in both murine and human fibrosis. In conclusion, activation of c-Rel/Pfkfb3 in damaged tissue instigates a paracrine signalling network among epithelial, myeloid and mesenchymal cells to stimulate fibrogenesis. Targeting the c-Rel-Pfkfb3 axis has potential for therapeutic applications in fibrotic disease.

DNA methylation repels binding of hypoxia-inducible transcription factors to maintain tumor immunotolerance.

BACKGROUND: Hypoxia is pervasive in cancer and other diseases. Cells sense and adapt to hypoxia by activating hypoxia-inducible transcription factors (HIFs), but it is still an outstanding question why cell types differ in their transcriptional response to hypoxia. RESULTS: We report that HIFs fail to bind CpG dinucleotides that are methylated in their consensus binding sequence, both in in vitro biochemical binding assays and in vivo studies of differentially methylated isogenic cell lines. Based on in silico structural modeling, we show that 5-methylcytosine indeed causes steric hindrance in the HIF binding pocket. A model wherein cell-type-specific methylation landscapes, as laid down by the differential expression and binding of other transcription factors under normoxia, control cell-type-specific hypoxia responses is observed. We also discover ectopic HIF binding sites in repeat regions which are normally methylated. Genetic and pharmacological DNA demethylation, but also cancer-associated DNA hypomethylation, expose these binding sites, inducing HIF-dependent expression of cryptic transcripts. In line with such cryptic transcripts being more prone to cause double-stranded RNA and viral mimicry, we observe low DNA methylation and high cryptic transcript expression in tumors with high immune checkpoint expression, but not in tumors with low immune checkpoint expression, where they would compromise tumor immunotolerance. In a low-immunogenic tumor model, DNA demethylation upregulates cryptic transcript expression in a HIF-dependent manner, causing immune activation and reducing tumor growth. CONCLUSIONS: Our data elucidate the mechanism underlying cell-type-specific responses to hypoxia and suggest DNA methylation and hypoxia to underlie tumor immunotolerance.

BIOMEX: an interactive workflow for (single cell) omics data interpretation and visualization.

The amount of biological data, generated with (single cell) omics technologies, is rapidly increasing, thereby exacerbating bottlenecks in the data analysis and interpretation of omics experiments. Data mining platforms that facilitate non-bioinformatician experimental scientists to analyze a wide range of experimental designs and data types can alleviate such bottlenecks, aiding in the exploration of (newly generated or publicly available) omics datasets. Here, we present BIOMEX, a browser-based software, designed to facilitate the Biological Interpretation Of Multi-omics EXperiments by bench scientists. BIOMEX integrates state-of-the-art statistical tools and field-tested algorithms into a flexible but well-defined workflow that accommodates metabolomics, transcriptomics, proteomics, mass cytometry and single cell data from different platforms and organisms. The BIOMEX workflow is accompanied by a manual and video tutorials that provide the necessary background to navigate the interface and get acquainted with the employed methods. BIOMEX guides the user through omics-tailored analyses, such as data pretreatment and normalization, dimensionality reduction, differential and enrichment analysis, pathway mapping, clustering, marker analysis, trajectory inference, meta-analysis and others. BIOMEX is fully interactive, allowing users to easily change parameters and generate customized plots exportable as high-quality publication-ready figures. BIOMEX is open source and freely available at https://www.vibcancer.be/software-tools/biomex.

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.

Single-Cell Transcriptome Atlas of Murine Endothelial Cells.

The heterogeneity of endothelial cells (ECs) across tissues remains incompletely inventoried. We constructed an atlas of >32,000 single-EC transcriptomes from 11 mouse tissues and identified 78 EC subclusters, including Aqp7+ intestinal capillaries and angiogenic ECs in healthy tissues. ECs from brain/testis, liver/spleen, small intestine/colon, and skeletal muscle/heart pairwise expressed partially overlapping marker genes. Arterial, venous, and lymphatic ECs shared more markers in more tissues than did heterogeneous capillary ECs. ECs from different vascular beds (arteries, capillaries, veins, lymphatics) exhibited transcriptome similarity across tissues, but the tissue (rather than the vessel) type contributed to the EC heterogeneity. Metabolic transcriptome analysis revealed a similar tissue-grouping phenomenon of ECs and heterogeneous metabolic gene signatures in ECs between tissues and between vascular beds within a single tissue in a tissue-type-dependent pattern. The EC atlas taxonomy enabled identification of EC subclusters in public scRNA-seq datasets and provides a powerful discovery tool and resource value.

An Integrated Gene Expression Landscape Profiling Approach to Identify Lung Tumor Endothelial Cell Heterogeneity and Angiogenic Candidates.

Heterogeneity of lung tumor endothelial cell (TEC) phenotypes across patients, species (human/mouse), and models (in vivo/in vitro) remains poorly inventoried at the single-cell level. We single-cell RNA (scRNA)-sequenced 56,771 endothelial cells from human/mouse (peri)-tumoral lung and cultured human lung TECs, and detected 17 known and 16 previously unrecognized phenotypes, including TECs putatively regulating immune surveillance. We resolved the canonical tip TECs into a known migratory tip and a putative basement-membrane remodeling breach phenotype. Tip TEC signatures correlated with patient survival, and tip/breach TECs were most sensitive to vascular endothelial growth factor blockade. Only tip TECs were congruent across species/models and shared conserved markers. Integrated analysis of the scRNA-sequenced data with orthogonal multi-omics and meta-analysis data across different human tumors, validated by functional analysis, identified collagen modification as a candidate angiogenic pathway.

Single-Cell RNA Sequencing Reveals Renal Endothelium Heterogeneity and Metabolic Adaptation to Water Deprivation.

BACKGROUND: Renal endothelial cells from glomerular, cortical, and medullary kidney compartments are exposed to different microenvironmental conditions and support specific kidney processes. However, the heterogeneous phenotypes of these cells remain incompletely inventoried. Osmotic homeostasis is vitally important for regulating cell volume and function, and in mammals, osmotic equilibrium is regulated through the countercurrent system in the renal medulla, where water exchange through endothelium occurs against an osmotic pressure gradient. Dehydration exposes medullary renal endothelial cells to extreme hyperosmolarity, and how these cells adapt to and survive in this hypertonic milieu is unknown. METHODS: We inventoried renal endothelial cell heterogeneity by single-cell RNA sequencing >40,000 mouse renal endothelial cells, and studied transcriptome changes during osmotic adaptation upon water deprivation. We validated our findings by immunostaining and functionally by targeting oxidative phosphorylation in a hyperosmolarity model in vitro and in dehydrated mice in vivo. RESULTS: We identified 24 renal endothelial cell phenotypes (of which eight were novel), highlighting extensive heterogeneity of these cells between and within the cortex, glomeruli, and medulla. In response to dehydration and hypertonicity, medullary renal endothelial cells upregulated the expression of genes involved in the hypoxia response, glycolysis, and-surprisingly-oxidative phosphorylation. Endothelial cells increased oxygen consumption when exposed to hyperosmolarity, whereas blocking oxidative phosphorylation compromised endothelial cell viability during hyperosmotic stress and impaired urine concentration during dehydration. CONCLUSIONS: This study provides a high-resolution atlas of the renal endothelium and highlights extensive renal endothelial cell phenotypic heterogeneity, as well as a previously unrecognized role of oxidative phosphorylation in the metabolic adaptation of medullary renal endothelial cells to water deprivation.

Leigh Syndrome Mouse Model Can Be Rescued by Interventions that Normalize Brain Hyperoxia, but Not HIF Activation.

Leigh syndrome is a devastating mitochondrial disease for which there are no proven therapies. We previously showed that breathing chronic, continuous hypoxia can prevent and even reverse neurological disease in the Ndufs4 knockout (KO) mouse model of complex I (CI) deficiency and Leigh syndrome. Here, we show that genetic activation of the hypoxia-inducible factor transcriptional program via any of four different strategies is insufficient to rescue disease. Rather, we observe an age-dependent decline in whole-body oxygen consumption. These mice exhibit brain tissue hyperoxia, which is normalized by hypoxic breathing. Alternative experimental strategies to reduce oxygen delivery, including breathing carbon monoxide (600 ppm in air) or severe anemia, can reverse neurological disease. Therefore, unused oxygen is the most likely culprit in the pathology of this disease. While pharmacologic activation of the hypoxia response is unlikely to alleviate disease in vivo, interventions that safely normalize brain tissue hyperoxia may hold therapeutic potential.

EndoDB: a database of endothelial cell transcriptomics data.

Endothelial cells (ECs) line blood vessels, regulate homeostatic processes (blood flow, immune cell trafficking), but are also involved in many prevalent diseases. The increasing use of high-throughput technologies such as gene expression microarrays and (single cell) RNA sequencing generated a wealth of data on the molecular basis of EC (dys-)function. Extracting biological insight from these datasets is challenging for scientists who are not proficient in bioinformatics. To facilitate the re-use of publicly available EC transcriptomics data, we developed the endothelial database EndoDB, a web-accessible collection of expert curated, quality assured and pre-analyzed data collected from 360 datasets comprising a total of 4741 bulk and 5847 single cell endothelial transcriptomes from six different organisms. Unlike other added-value databases, EndoDB allows to easily retrieve and explore data of specific studies, determine under which conditions genes and pathways of interest are deregulated and assess reprogramming of metabolism via principal component analysis, differential gene expression analysis, gene set enrichment analysis, heatmaps and metabolic and transcription factor analysis, while single cell data are visualized as gene expression color-coded t-SNE plots. Plots and tables in EndoDB are customizable, downloadable and interactive. EndoDB is freely available at https://vibcancer.be/software-tools/endodb, and will be updated to include new studies.