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.

Hematopoietic stem cells undergo a lymphoid to myeloid switch in early stages of emergency granulopoiesis.

Emergency granulopoiesis is the enhanced and accelerated production of granulocytes that occurs during acute infection. The contribution of hematopoietic stem cells (HSCs) to this process was reported; however, how HSCs participate in emergency granulopoiesis remains elusive. Here, using a mouse model of emergency granulopoiesis we observe transcriptional changes in HSCs as early as 4 h after lipopolysaccharide (LPS) administration. We observe that the HSC identity is changed towards a myeloid-biased HSC and show that CD201 is enriched in lymphoid-biased HSCs. While CD201 expression under steady-state conditions reveals a lymphoid bias, under emergency granulopoiesis loss of CD201 marks the lymphoid-to-myeloid transcriptional switch. Mechanistically, we determine that lymphoid-biased CD201+ HSCs act as a first response during emergency granulopoiesis due to direct sensing of LPS by TLR4 and downstream activation of NF-κΒ signaling. The myeloid-biased CD201- HSC population responds indirectly during an acute infection by sensing G-CSF, increasing STAT3 phosphorylation, and upregulating LAP/LAP* C/EBPß isoforms. In conclusion, HSC subpopulations support early phases of emergency granulopoiesis due to their transcriptional rewiring from a lymphoid-biased to myeloid-biased population and thus establishing alternative paths to supply elevated numbers of granulocytes.

Mitochondrial HER2 stimulates respiration and promotes tumorigenicity.

BACKGROUND: Amplification of HER2, a receptor tyrosine kinase and a breast cancer-linked oncogene, is associated with aggressive disease. HER2 protein is localised mostly at the cell membrane, but a fraction translocates to mitochondria. Whether and how mitochondrial HER2 contributes to tumorigenicity is currently unknown. METHODS: We enriched the mitochondrial (mt-)HER2 fraction in breast cancer cells using an N-terminal mitochondrial targeting sequence and analysed how this manipulation impacts bioenergetics and tumorigenic properties. The role of the tyrosine kinase activity of mt-HER2 was assessed in wild type, kinase-dead (K753M) and kinase-enhanced (V659E) mtHER2 constructs. RESULTS: We document that mt-HER2 associates with the oxidative phosphorylation system, stimulates bioenergetics and promotes larger respiratory supercomplexes. mt-HER2 enhances proliferation and invasiveness in vitro and tumour growth and metastatic potential in vivo, in a kinase activity-dependent manner. On the other hand, constitutively active mt-HER2 provokes excessive mitochondria ROS generation, sensitises to cell death, and restricts growth of primary tumours, suggesting that regulation of HER2 activity in mitochondria is required for the maximal pro-tumorigenic effect. CONCLUSIONS: mt-HER2 promotes tumorigenicity by supporting bioenergetics and optimal redox balance.

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.

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.

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.

The role of Her2 and other oncogenes of the PI3K/AKT pathway in mitochondria.

Altered metabolism and resistance to cell death are typical hallmarks of cancer phenotype. Mitochondria are organelles central to cellular metabolism as well as to cell death induction. Hyperactivation of pro-survival and pro-proliferative pathways such as PI3K/AKT leads to cancer initiation, which affects mitochondria. Growing body of evidence indicates that oncogenes such as HER2, EGFR and RAS, as well as the downstream members of the PI3K/AKT signaling pathway, directly regulate mitochondria by translocating to the organelle. Here we discuss evidence of this scenario and consider mechanisms for direct regulation of mitochondrial function. Being in close proximity to mitochondrial bioenergetics machinery as well as to the regulators/executors of programed cell death, oncogenes in mitochondria may be ideally placed to perform this task. This represents a thus far under-explored area, which may be relevant to better understanding of cancer initiation, progression and treatment.

Spatial Analysis of Nucleotide Metabolism: From CRISPR Knockout Cancer Cells to MALDI Imaging of Tumors.

Cancer cells depend on nucleotides for proliferation. Inhibition of nucleotide metabolism by antimetabolites is a well-established anticancer therapy. However, resistance and toxicity to antimetabolite treatments reduce their effectiveness. Here, we focus on the pyrimidine de novo synthesis pathway, which is crucial for cancer cell proliferation, yet its pharmacological targeting in cancer has been without much clinical success so far. Hence, it is important to understand how cancer cells cope with the insufficiency of this pathway. Here, we describe a procedure to prepare subcutaneous tumor model deficient in de novo pyrimidine synthesis. For examination of metabolic responses to de novo synthesis blockade in tumors, we propose application of MALDI imaging that allows spatially resolved examination of metabolic responses to de novo synthesis blockade in tumors.

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.

Mitochondrially targeted tamoxifen in patients with metastatic solid tumours: an open-label, phase I/Ib single-centre trial.

Background: Mitochondria present an emerging target for cancer treatment. We have investigated the effect of mitochondrially targeted tamoxifen (MitoTam), a first-in-class anti-cancer agent, in patients with solid metastatic tumours. Methods: MitoTam was tested in an open-label, single-centre (Department of Oncology, General Faculty Hospital, Charles University, Czech Republic), phase I/Ib trial in metastatic patients with various malignancies and terminated oncological therapies. In total, 75 patients were enrolled between May 23, 2018 and July 22, 2020. Phase I evaluated escalating doses of MitoTam in two therapeutic regimens using the 3 + 3 design to establish drug safety and maximum tolerated dose (MTD). In phase Ib, three dosing regimens were applied over 8 and 6 weeks to evaluate long-term toxicity of MitoTam as the primary objective and its anti-cancer effect as a secondary objective. This trial was registered with the European Medicines Agency under EudraCT 2017-004441-25. Findings: In total, 37 patients were enrolled into phase I and 38 into phase Ib. In phase I, the initial application of MitoTam via peripheral vein indicated high risk of thrombophlebitis, which was avoided by central vein administration. The highest dose with acceptable side effects was 5.0 mg/kg. The prevailing adverse effects (AEs) in phase I were neutropenia (30%), anaemia (30%) and fever/hyperthermia (30%), and in phase Ib fever/hyperthermia (58%) together with anaemia (26%) and neutropenia (16%). Serious AEs were mostly related to thromboembolic (TE) complications that affected 5% and 13% of patients in phase I and Ib, respectively. The only statistically significant AE related to MitoTam treatment was anaemia in phase Ib (p = 0.004). Of the tested regimens weekly dosing with 3.0 mg/kg for 6 weeks afforded the best safety profile with almost all being grade 1 (G1) AEs. Altogether, five fatalities occurred during the study, two of them meeting criteria for Suspected Unexpected Serious Adverse Events Reporting (SUSAR) (G4 thrombocytopenia and G5 stroke). MitoTam showed benefit evaluated as clinical benefit rate (CBR) in 37% patients with the largest effect in renal cell carcinoma (RCC) where four out of six patients reached disease stabilisation (SD), one reached partial response (PR) so that in total, five out of six (83%) patients showed CBR. Interpretation: In this study, the MTD was established as 5.0 mg/kg and the recommended dose of MitoTam as 3.0 mg/kg given once per week via central vein with recommended preventive anti-coagulation therapy. The prevailing toxicity included haematological AEs, hyperthermia/fever and TE complications. One fatal stroke and non-fatal G4 thrombocytopenia were recorded. MitoTam showed high efficacy against RCC. Funding: Smart Brain Ltd. Translation: For the Czech translation of the abstract see Supplementary Materials section.

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.

Effects of metabolic cancer therapy on tumor microenvironment.

Targeting tumor metabolism for cancer therapy is an old strategy. In fact, historically the first effective cancer therapeutics were directed at nucleotide metabolism. The spectrum of metabolic drugs considered in cancer increases rapidly - clinical trials are in progress for agents directed at glycolysis, oxidative phosphorylation, glutaminolysis and several others. These pathways are essential for cancer cell proliferation and redox homeostasis, but are also required, to various degrees, in other cell types present in the tumor microenvironment, including immune cells, endothelial cells and fibroblasts. How metabolism-targeted treatments impact these tumor-associated cell types is not fully understood, even though their response may co-determine the overall effectivity of therapy. Indeed, the metabolic dependencies of stromal cells have been overlooked for a long time. Therefore, it is important that metabolic therapy is considered in the context of tumor microenvironment, as understanding the metabolic vulnerabilities of both cancer and stromal cells can guide new treatment concepts and help better understand treatment resistance. In this review we discuss recent findings covering the impact of metabolic interventions on cellular components of the tumor microenvironment and their implications for metabolic cancer therapy.

Mitochondrial respiration supports autophagy to provide stress resistance during quiescence.

Mitochondrial oxidative phosphorylation (OXPHOS) generates ATP, but OXPHOS also supports biosynthesis during proliferation. In contrast, the role of OXPHOS during quiescence, beyond ATP production, is not well understood. Using mouse models of inducible OXPHOS deficiency in all cell types or specifically in the vascular endothelium that negligibly relies on OXPHOS-derived ATP, we show that selectively during quiescence OXPHOS provides oxidative stress resistance by supporting macroautophagy/autophagy. Mechanistically, OXPHOS constitutively generates low levels of endogenous ROS that induce autophagy via attenuation of ATG4B activity, which provides protection from ROS insult. Physiologically, the OXPHOS-autophagy system (i) protects healthy tissue from toxicity of ROS-based anticancer therapy, and (ii) provides ROS resistance in the endothelium, ameliorating systemic LPS-induced inflammation as well as inflammatory bowel disease. Hence, cells acquired mitochondria during evolution to profit from oxidative metabolism, but also built in an autophagy-based ROS-induced protective mechanism to guard against oxidative stress associated with OXPHOS function during quiescence.Abbreviations: AMPK: AMP-activated protein kinase; AOX: alternative oxidase; Baf A: bafilomycin A1; CI, respiratory complexes I; DCF-DA: 2',7'-dichlordihydrofluorescein diacetate; DHE: dihydroethidium; DSS: dextran sodium sulfate; ΔΨmi: mitochondrial inner membrane potential; EdU: 5-ethynyl-2'-deoxyuridine; ETC: electron transport chain; FA: formaldehyde; HUVEC; human umbilical cord endothelial cells; IBD: inflammatory bowel disease; LC3B: microtubule associated protein 1 light chain 3 beta; LPS: lipopolysaccharide; MEFs: mouse embryonic fibroblasts; MTORC1: mechanistic target of rapamycin kinase complex 1; mtDNA: mitochondrial DNA; NAC: N-acetyl cysteine; OXPHOS: oxidative phosphorylation; PCs: proliferating cells; PE: phosphatidylethanolamine; PEITC: phenethyl isothiocyanate; QCs: quiescent cells; ROS: reactive oxygen species; PLA2: phospholipase A2, WB: western blot.

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

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.

Heterogeneous Effects of Calorie Content and Nutritional Components Underlie Dietary Influence on Pancreatic Cancer Susceptibility.

Pancreatic cancer is a rare but fatal form of cancer, the fourth highest in absolute mortality. Known risk factors include obesity, diet, and type 2 diabetes; however, the low incidence rate and interconnection of these factors confound the isolation of individual effects. Here, we use epidemiological analysis of prospective human cohorts and parallel tracking of pancreatic cancer in mice to dissect the effects of obesity, diet, and diabetes on pancreatic cancer. Through longitudinal monitoring and multi-omics analysis in mice, we found distinct effects of protein, sugar, and fat dietary components, with dietary sugars increasing Mad2l1 expression and tumor proliferation. Using epidemiological approaches in humans, we find that dietary sugars give a MAD2L1 genotype-dependent increased susceptibility to pancreatic cancer. The translation of these results to a clinical setting could aid in the identification of the at-risk population for screening and potentially harness dietary modification as a therapeutic measure.

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.