Single-cell RNA sequencing of cystic fibrosis liver disease explants reveals endothelial complement activation.

BACKGROUND & AIMS: Cystic fibrosis (CF) is considered a multisystemic disorder in which CF-associated liver disease (CFLD) is the third most common cause of mortality. Currently, no effective treatment is available for CFLD because its pathophysiology is still unclear. Interestingly, CFLD exhibits identical vascular characteristics as non-cirrhotic portal hypertension, recently classified as porto-sinusoidal vascular disorders (PSVD). METHODS: Since endothelial cells (ECs) are an important component in PSVD, we performed single-cell RNA sequencing (scRNA-seq) on four explant livers from CFLD patients to identify differential endothelial characteristics which could contribute to the disease. We comprehensively characterized the endothelial compartment and compared it with publicly available scRNA-seq datasets from cirrhotic and healthy livers. Key gene signatures were validated ex vivo on patient tissues. RESULTS: We found that ECs from CF liver explants are more closely related to healthy than cirrhotic patients. In CF patients we also discovered a distinct population of liver sinusoidal ECs-coined CF LSECs-upregulating genes involved in the complement cascade and coagulation. Finally, our immunostainings further validated the predominant periportal location of CF LSECs. CONCLUSIONS: Our work showed novel aspects of human liver ECs at the single-cell level thereby supporting endothelial involvement in CFLD, and reinforcing the hypothesis that ECs could be a driver of PSVD. Therefore, considering the vascular compartment in CF and CFLD may help developing new therapeutic approaches for these diseases.

Detailed protocol for a corneal thermal cauterization-based (lymph-)angiogenesis assay in mice.

Angiogenesis and lymphangiogenesis, the formation of new blood or lymphatic vessels, respectively, from preexisting vasculature is essential during embryonic development, but also occurs during tissue repair and in pathological conditions (cancer; ocular disease; ischemic, infectious and inflammatory disorders), which are all characterized to a certain extent by inflammatory conditions. Hence, a rapid, inexpensive, feasible / technically easy, reliable assay of inflammation-induced (lymph-)angiogenesis is highly valuable. In this context, the corneal thermal cauterization assay in mice is a simple, low-cost, reproducible, insightful and labor-saving assay to gauge the role of inflammation in angiogenesis and lymphangiogenesis. However, to the best of our knowledge, there is no standardized protocol to perform this assay. Here, we provide a step-by-step description of the model's procedures, which include:•The thermal cauterization of the corneas,•Enucleation and dissection of the corneas,•Subsequent immunofluorescence staining of the neovasculature, and morphometric analysis. We also discuss ethical considerations and aspects related to animal welfare guidelines. Altogether, this paper will help to increase the reproducibility of the corneal thermal cauterization model and facilitate its use for angiogenesis and lymphangiogenesis research.

Intracellular BAPTA directly inhibits PFKFB3, thereby impeding mTORC1-driven Mcl-1 translation and killing MCL-1-addicted cancer cells.

Intracellular Ca2+ signals control several physiological and pathophysiological processes. The main tool to chelate intracellular Ca2+ is intracellular BAPTA (BAPTAi), usually introduced into cells as a membrane-permeant acetoxymethyl ester (BAPTA-AM). Previously, we demonstrated that BAPTAi enhanced apoptosis induced by venetoclax, a BCL-2 antagonist, in diffuse large B-cell lymphoma (DLBCL). This finding implied a novel interplay between intracellular Ca2+ signaling and anti-apoptotic BCL-2 function. Hence, we set out to identify the underlying mechanisms by which BAPTAi enhances cell death in B-cell cancers. In this study, we discovered that BAPTAi alone induced apoptosis in hematological cancer cell lines that were highly sensitive to S63845, an MCL-1 antagonist. BAPTAi provoked a rapid decline in MCL-1-protein levels by inhibiting mTORC1-driven Mcl-1 translation. These events were not a consequence of cell death, as BAX/BAK-deficient cancer cells exhibited similar downregulation of mTORC1 activity and MCL-1-protein levels. Next, we investigated how BAPTAi diminished mTORC1 activity and identified its ability to impair glycolysis by directly inhibiting 6-phosphofructo-2-kinase/fructose-2,6-bisphosphatase 3 (PFKFB3) activity, a previously unknown effect of BAPTAi. Notably, these effects were also induced by a BAPTAi analog with low affinity for Ca2+. Consequently, our findings uncover PFKFB3 inhibition as an Ca2+-independent mechanism through which BAPTAi impairs cellular metabolism and ultimately compromises the survival of MCL-1-dependent cancer cells. These findings hold two important implications. Firstly, the direct inhibition of PFKFB3 emerges as a key regulator of mTORC1 activity and a promising target in MCL-1-dependent cancers. Secondly, cellular effects caused by BAPTAi are not necessarily related to Ca2+ signaling. Our data support the need for a reassessment of the role of Ca2+ in cellular processes when findings were based on the use of BAPTAi.

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.

A high-throughput screening campaign against PFKFB3 identified potential inhibitors with novel scaffolds.

The growth of solid tumors depends on tumor vascularization and the endothelial cells (ECs) that line the lumen of blood vessels. ECs generate a large fraction of ATP through glycolysis, and elevation of their glycolytic activity is associated with angiogenic behavior in solid tumors. 6-Phosphofructo-2-kinase/fructose-2,6-bisphosphatase 3 (PFKFB3) positively regulates glycolysis via fructose-2/6-bisphosphate, the product of its kinase activity. Partial inhibition of glycolysis in tumor ECs by targeting PFKFB3 normalizes the otherwise abnormal tumor vessels, thereby reducing metastasis and improving the outcome of chemotherapy. Although a limited number of tool compounds exist, orally available PFKFB3 inhibitors are unavailable. In this study we conducted a high-throughput screening campaign against the kinase activity of PFKFB3, involving 250,240 chemical compounds. A total of 507 initial hits showing >50% inhibition at 20 µM were identified, 66 of them plus 1 analog from a similarity search consistently displayed low IC50 values (<10 µM). In vitro experiments yielded 22 nontoxic hits that suppressed the tube formation of primary human umbilical vein ECs at 10 µM. Of them, 15 exhibited binding affinity to PFKFB3 in surface plasmon resonance assays, including 3 (WNN0403-E003, WNN1352-H007 and WNN1542-F004) that passed the pan-assay interference compounds screening without warning flags. This study provides potential leads to the development of new PFKFB3 inhibitors.

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.

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

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.

Serine metabolism remodeling after platinum-based chemotherapy identifies vulnerabilities in a subgroup of resistant ovarian cancers.

Resistance to platinum-based chemotherapy represents a major clinical challenge for many tumors, including epithelial ovarian cancer. Patients often experience several response-relapse events, until tumors become resistant and life expectancy drops to 12-15 months. Despite improved knowledge of the molecular determinants of platinum resistance, the lack of clinical applicability limits exploitation of many potential targets, leaving patients with limited options. Serine biosynthesis has been linked to cancer growth and poor prognosis in various cancer types, however its role in platinum-resistant ovarian cancer is not known. Here, we show that a subgroup of resistant tumors decreases phosphoglycerate dehydrogenase (PHGDH) expression at relapse after platinum-based chemotherapy. Mechanistically, we observe that this phenomenon is accompanied by a specific oxidized nicotinamide adenine dinucleotide (NAD+) regenerating phenotype, which helps tumor cells in sustaining Poly (ADP-ribose) polymerase (PARP) activity under platinum treatment. Our findings reveal metabolic vulnerabilities with clinical implications for a subset of platinum resistant ovarian cancers.

Skeletal progenitors preserve proliferation and self-renewal upon inhibition of mitochondrial respiration by rerouting the TCA cycle.

A functional electron transport chain (ETC) is crucial for supporting bioenergetics and biosynthesis. Accordingly, ETC inhibition decreases proliferation in cancer cells but does not seem to impair stem cell proliferation. However, it remains unclear how stem cells metabolically adapt. In this study, we show that pharmacological inhibition of complex III of the ETC in skeletal stem and progenitor cells induces glycolysis side pathways and reroutes the tricarboxylic acid (TCA) cycle to regenerate NAD+ and preserve cell proliferation. These metabolic changes also culminate in increased succinate and 2-hydroxyglutarate levels that inhibit Ten-eleven translocation (TET) DNA demethylase activity, thereby preserving self-renewal and multilineage potential. Mechanistically, mitochondrial malate dehydrogenase and reverse succinate dehydrogenase activity proved to be essential for the metabolic rewiring in response to ETC inhibition. Together, these data show that the metabolic plasticity of skeletal stem and progenitor cells allows them to bypass ETC blockade and preserve their self-renewal.

A YAP/TAZ-TEAD signalling module links endothelial nutrient acquisition to angiogenic growth.

Angiogenesis, the process by which endothelial cells (ECs) form new blood vessels from existing ones, is intimately linked to the tissue's metabolic milieu and often occurs at nutrient-deficient sites. However, ECs rely on sufficient metabolic resources to support growth and proliferation. How endothelial nutrient acquisition and usage are regulated is unknown. Here we show that these processes are instructed by Yes-associated protein 1 (YAP)/WW domain-containing transcription regulator 1 (WWTR1/TAZ)-transcriptional enhanced associate domain (TEAD): a transcriptional module whose function is highly responsive to changes in the tissue environment. ECs lacking YAP/TAZ or their transcriptional partners, TEAD1, 2 and 4 fail to divide, resulting in stunted vascular growth in mice. Conversely, activation of TAZ, the more abundant paralogue in ECs, boosts proliferation, leading to vascular hyperplasia. We find that YAP/TAZ promote angiogenesis by fuelling nutrient-dependent mTORC1 signalling. By orchestrating the transcription of a repertoire of cell-surface transporters, including the large neutral amino acid transporter SLC7A5, YAP/TAZ-TEAD stimulate the import of amino acids and other essential nutrients, thereby enabling mTORC1 activation. Dissociating mTORC1 from these nutrient inputs-elicited by the loss of Rag GTPases-inhibits mTORC1 activity and prevents YAP/TAZ-dependent vascular growth. Together, these findings define a pivotal role for YAP/TAZ-TEAD in controlling endothelial mTORC1 and illustrate the essentiality of coordinated nutrient fluxes in the vasculature.

PHGDH heterogeneity potentiates cancer cell dissemination and metastasis.

Cancer metastasis requires the transient activation of cellular programs enabling dissemination and seeding in distant organs1. Genetic, transcriptional and translational heterogeneity contributes to this dynamic process2,3. Metabolic heterogeneity has also been observed4, yet its role in cancer progression is less explored. Here we find that the loss of phosphoglycerate dehydrogenase (PHGDH) potentiates metastatic dissemination. Specifically, we find that heterogeneous or low PHGDH expression in primary tumours of patients with breast cancer is associated with decreased metastasis-free survival time. In mice, circulating tumour cells and early metastatic lesions are enriched with Phgdhlow cancer cells, and silencing Phgdh in primary tumours increases metastasis formation. Mechanistically, Phgdh interacts with the glycolytic enzyme phosphofructokinase, and the loss of this interaction activates the hexosamine-sialic acid pathway, which provides precursors for protein glycosylation. As a consequence, aberrant protein glycosylation occurs, including increased sialylation of integrin αvß3, which potentiates cell migration and invasion. Inhibition of sialylation counteracts the metastatic ability of Phgdhlow cancer cells. In conclusion, although the catalytic activity of PHGDH supports cancer cell proliferation, low PHGDH protein expression non-catalytically potentiates cancer dissemination and metastasis formation. Thus, the presence of PHDGH heterogeneity in primary tumours could be considered a sign of tumour aggressiveness.

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.

Transcriptomic analysis of CFTR-impaired endothelial cells reveals a pro-inflammatory phenotype.

Cystic fibrosis (CF) is a life-threatening disorder characterised by decreased pulmonary mucociliary and pathogen clearance, and an exaggerated inflammatory response leading to progressive lung damage. CF is caused by bi-allelic pathogenic variants of the cystic fibrosis transmembrane conductance regulator (CFTR) gene, which encodes a chloride channel. CFTR is expressed in endothelial cells (ECs) and EC dysfunction has been reported in CF patients, but a role for this ion channel in ECs regarding CF disease progression is poorly described.We used an unbiased RNA sequencing approach in complementary models of CFTR silencing and blockade (by the CFTR inhibitor CFTRinh-172) in human ECs to characterise the changes upon CFTR impairment. Key findings were further validated in vitro and in vivo in CFTR-knockout mice and ex vivo in CF patient-derived ECs.Both models of CFTR impairment revealed that EC proliferation, migration and autophagy were downregulated. Remarkably though, defective CFTR function led to EC activation and a persisting pro-inflammatory state of the endothelium with increased leukocyte adhesion. Further validation in CFTR-knockout mice revealed enhanced leukocyte extravasation in lung and liver parenchyma associated with increased levels of EC activation markers. In addition, CF patient-derived ECs displayed increased EC activation markers and leukocyte adhesion, which was partially rescued by the CFTR modulators VX-770 and VX-809.Our integrated analysis thus suggests that ECs are no innocent bystanders in CF pathology, but rather may contribute to the exaggerated inflammatory phenotype, raising the question of whether normalisation of vascular inflammation might be a novel therapeutic strategy to ameliorate the disease severity of CF.

Macrophage-derived glutamine boosts satellite cells and muscle regeneration.

Muscle regeneration is sustained by infiltrating macrophages and the consequent activation of satellite cells1-4. Macrophages and satellite cells communicate in different ways1-5, but their metabolic interplay has not been investigated. Here we show, in a mouse model, that muscle injuries and ageing are characterized by intra-tissue restrictions of glutamine. Low levels of glutamine endow macrophages with the metabolic ability to secrete glutamine via enhanced glutamine synthetase (GS) activity, at the expense of glutamine oxidation mediated by glutamate dehydrogenase 1 (GLUD1). Glud1-knockout macrophages display constitutively high GS activity, which prevents glutamine shortages. The uptake of macrophage-derived glutamine by satellite cells through the glutamine transporter SLC1A5 activates mTOR and promotes the proliferation and differentiation of satellite cells. Consequently, macrophage-specific deletion or pharmacological inhibition of GLUD1 improves muscle regeneration and functional recovery in response to acute injury, ischaemia or ageing. Conversely, SLC1A5 blockade in satellite cells or GS inactivation in macrophages negatively affects satellite cell functions and muscle regeneration. These results highlight the metabolic crosstalk between satellite cells and macrophages, in which macrophage-derived glutamine sustains the functions of satellite cells. Thus, the targeting of GLUD1 may offer therapeutic opportunities for the regeneration of injured or aged muscles.

Basic and Therapeutic Aspects of Angiogenesis Updated.

All organisms growing beyond the oxygen diffusion limit critically depend on a functional vasculature for survival. Yet blood vessels are far more than passive, uniform conduits for oxygen and nutrient supply. A remarkable organotypic heterogeneity is brought about by tissue-specific differentiated endothelial cells (lining the blood vessels' lumen) and allows blood vessels to deal with organ-specific demands for homeostasis. On the flip side, when blood vessels go awry, they promote life-threatening diseases characterized by endothelial cells inappropriately adopting an angiogenic state (eg, tumor vascularization) or becoming dysfunctional (eg, diabetic microvasculopathies), calling respectively for antiangiogenic therapies and proangiogenic/vascular regenerative strategies. In solid tumors, despite initial enthusiasm, growth factor-based (mostly anti-VEGF [vascular endothelial growth factor]) antiangiogenic therapies do not sufficiently live up to the expectations in terms of efficiency and patient survival, in part, due to intrinsic and acquired therapy resistance. Tumors cunningly deploy alternative growth factors than the ones targeted by the antiangiogenic therapies to reinstigate angiogenesis or revert to other ways of securing blood flow, independently of the targeted growth factors. In trying to alleviate tissue ischemia and to repair dysfunctional or damaged endothelium, local in-tissue administration of (genes encoding) proangiogenic factors or endothelial (stem) cells harnessing regenerative potential have been explored. Notwithstanding evaluation in clinical trials, these approaches are often hampered by dosing issues and limited half-life or local retention of the administered agents. Here, without intending to provide an all-encompassing historical overview, we focus on some recent advances in understanding endothelial cell behavior in health and disease and identify novel molecular players and concepts that could eventually be considered for therapeutic targeting.

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.