# EndoDB: a database of endothelial cell transcriptomics data

Shawez Khan, Federico Taverna, Katerina Rohlenova, Lucas Treps, Vincent Geldhof, Laura de Rooij, Liliana Sokol, Andreas Pircher, Lena-Christin Conradi, Joanna Kalucka, et al.

<figure><img src="image_1.png" /></figure>

To cite this version:

Shawez Khan, Federico Taverna, Katerina Rohlenova, Lucas Treps, Vincent Geldhof, et al. EndoDB: a database of endothelial cell transcriptomics data. Nucleic Acids Research, 2018, 47 (D1), pp.D736 - D744. 10.1093/nar/gky997. inserm-05169336

HAL Id: inserm-05169336

https://inserm.hal.science/inserm-05169336v1

Submitted on 18 Jul 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# EndoDB: a database of endothelial cell transcriptomics data

Shawez Khan<sup>1,2,3,†</sup>, Federico Taverna<sup>1,3,†</sup>, Katerina Rohlenova<sup>1,3,†</sup>, Lucas Treps<sup>1,3</sup>, Vincent Geldhof<sup>1,3</sup>, Laura de Rooij<sup>1,3</sup>, Liliana Sokol<sup>1,3</sup>, Andreas Pircher<sup>1,3</sup>, Lena-Christin Conradi<sup>1,3</sup>, Joanna Kalucka<sup>®1,3</sup>, Luc Schoonjans<sup>1,2,3</sup>, Guy Eelen<sup>1,3</sup>, Mieke Dewerchin<sup>1,3</sup>, Tobias Karakach<sup>1,3</sup>, Xuri Li<sup>2,*</sup>, Jermaine Goveia<sup>1,3,*</sup> and Peter Carmeliet<sup>1,2,3,*</sup>

<sup>1</sup>Department of Oncology and Leuven Cancer Institute (LKI), Laboratory of Angiogenesis and Vascular Metabolism, KU Leuven, 3000 Leuven, Belgium, <sup>2</sup>State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center, Sun Yat-Sen University, Guangzhou 510060, Guangdong, P.R. China and <sup>3</sup>Laboratory of Angiogenesis and Vascular Metabolism, Center for Cancer Biology, VIB, 3000 Leuven, Belgium

Received August 14, 2018; Revised September 20, 2018; Editorial Decision October 03, 2018; Accepted October 09, 2018

## ABSTRACT

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.

Endothelial cells (ECs) line the lumen of blood vessels, are metabolically active and orchestrate important processes such as vasomotor tone, coagulation, permeability, tissue vascularization and immune response (1). The inability of ECs to fulfill their physiological functions is a key feature of multiple prevalent diseases including hypertension, atherosclerosis, diabetes and cancer (2). In the past three decades, hypothesis-driven studies have started to unravel the molecular mechanisms that underlie EC function in health and disease, which resulted in clinically approved EC-targeting drugs such as bevacizumab and ranibizumab for cancer and age-related wet macular degeneration, respectively (3). Recent technological advances have made it possible to perform global profiling of transcript (transcriptomics), protein (proteomics) and metabolite levels (metabolomics), even in single cells. These technologies have provided unprecedented insight in EC biology and generated many novel hypotheses (4).

To facilitate the re-use of omics datasets, it is now a standard requirement of most scientific journals to make published data publicly available via repositories such as Ar-rayExpress (5) and Gene Expression Omnibus (GEO) (6).

<footnote>*To whom correspondence should be addressed. Tel: +32 16 373 204; Fax: +32 16 372 585; Email: peter.carmeliet@kuleuven.vib.be
Correspondence may also be addressed to Jermaine Goveia. Tel: +32 16 373 204; Fax: +32 16 372 585; Email: jermaine.goveia@kuleuven.vib.be
Correspondence may also be addressed to Xuri Li. Tel: +86 20 8733 1815; Fax: +86 20 8733 1815; Email: lixr6@mail.sysu.edu.cn
†The authors wish it to be known that, in their opinion, the first three authors should be regarded as Joint First Authors.</footnote>

<footnote>Present addresses:
Andreas Pircher, Department of Hematology and Oncology, Internal Medicine V, Medical University Innsbruck, Innsbruck, Austria.
Lena-Christin Conradi, Clinic of General, Visceral and Pediatric Surgery, University Medical Center Göttingen, Göttingen, Germany.</footnote>
However, the main purpose of these primary archives is to store raw data in a format that is accessible by bioinformaticians, not to enable bench scientists and experimentalists to directly explore and re-use data to answer biological questions (7). In the cancer field, added-value databases (databases that are specifically designed to allow bench scientists to explore pre-selected, highly curated and pre-analyzed data) have been transformative and facilitated several breakthrough discoveries (8–13). A few independent research groups have made the results of their EC profiling efforts web-accessible (14,15), but currently no comprehensive EC-specific added-value databases exist to facilitate the reuse of gene expression data by the vascular biology community. Here, we describe the development of the endothelial cell database (EndoDB), the first freely available added-value database of EC (single cell) transcriptomics studies, designed to allow vascular biologists and other bench scientists to unlock the untapped potential of publicly available data via an easy-to-use interactive web interface (https://vibcancer.be/software-tools/endodb). The user-interface and functionality of the EndoDB were developed in close collaboration between software developers, bioinformaticians and vascular biologists. The EndoDB is easy-to-use, allows to interactively explore data using powerful statistical and bioinformatics approaches, and has been field-tested and used in several recent publications (16,17).

# MATERIALS AND METHODS

## Retrieval of EC transcriptomics data and database content

We first aimed to collect and organize EC transcriptomics datasets available in the public domain. To do this, we used a two-step approach (Figure 1 and Supplementary Table S1). First, we used a broad and sensitive filter (‘endothelial OR endothelium’) to search ArrayExpress and GEO for EC transcriptomics datasets, which returned 1121 studies (as of 1 June 2018). We then screened the title, abstract and if necessary sample information to determine which of these studies performed gene expression profiling in ECs specifically. We excluded 663 studies based on the title abstract screen and 113 studies for which data was not made available in the public domain, resulting in a total of 345 studies comprising 357 bulk transcriptomics datasets (a single study can contain multiple independent datasets).

Second, we complemented our search with a PubMed screen to identify recently published single cell transcriptomics studies of ECs, which retrieved 10 single cell RNA sequencing (scRNA-seq) studies (as of 1 June 2018). To ensure high quality scRNA-seq data in the EndoDB, we only considered studies that aimed to perform scRNA-seq in ECs specifically, used isolation protocols optimized for this purpose and sequenced >1000 ECs because: (i) ECs are a minority cell type in all tissues, ranging from ~1–5% in tumor tissue to ~20% in highly vascularized organs such as the lung (18,19); (ii) generic dissociation protocols fail to adequately capture ECs, resulting in a further depletion of the EC fraction (20); (iii) ECs are sensitive to dissociation-induced artifacts (20), which biases data interpretation (21), together indicating that dissociation protocols have to be

specifically optimized to isolate a high number of high quality ECs; (iv) the reliability of single cell analyses depends on the number of cells analyzed (22,23); and (v) the number of ECs that can be sequenced with more advanced technology (e.g. droplet-based sequencing) will increase by orders of magnitude. Based on these criteria, we excluded eight studies (Supplementary Table S2), and at last included two studies comprising three datasets.

The final database contains 4741 bulk and 5847 single cell EC transcriptomes from six different model organisms (*Homo sapiens* [human], *Mus musculus* [mouse], *Rattus norvegicus* [rat], *Bos taurus* [cow], *Danio rerio* [zebrafish], *Sus scrofa* [pig]), generated using a variety of technologies (micro-arrays, (single cell) RNA sequencing) and platforms (Affymetrix, Agilent, Illumina).

## Pre-processing and quality control of bulk transcriptomics data

To maximize the number of datasets included in the EndoDB, we aimed to include data from all micro-array and RNA-sequencing platforms. For micro-arrays, we retrieved raw data files when possible, but included log2-transformed pre-processed data when raw data was not available. Raw data from the Affymetrix platform was downloaded as binary or text CEL files, normalized using the robust multi-array average using the R-packages *affy* and *oligo*, followed by quantile normalization, log2 transformation and probe-set summarization (24,25). Raw data files from Illumina BeadChips were normalized using the *normexp* algorithm using control probes (neqc) from the *limma* package (26). Raw data from Agilent platforms was normalized using the *normalizeWithinArrays* and *normalizeBetweenArrays* functions available from *limma* package (26). Probe set identifiers were mapped to Entrez identifiers, official gene symbols and gene names using annotation packages for the corresponding platform. We applied standard quality controls using the *arrayQualityMetrics* R package to identify low-quality arrays (27). In this approach, samples were scored on the basis of three metrics, the distances between the arrays, Bland–Altman plots (MA plots) and boxplots. Samples that were flagged as outliers based on more than one of these metrics were removed from the analysis in the EndoDB, but the raw data of all samples can be downloaded for downstream analyses in independent software pipelines. At last, we randomly selected 25% of studies and cross-checked our results for congruency with the original publication.

For RNA-sequencing, data were retrieved in fastq format, and alignment to the reference genome and quantification of transcript abundances was performed using the *Kallisto* software (28). We performed trimmed mean of M-values (TMM)-normalization using the R-package *EdgeR* (29) and VOOM-normalization using *limma* before downstream differential analysis.

## Pre-processing and quality control of single cell data

Recent technological breakthroughs have made it possible to RNA sequence transcriptomes of single cells (scRNA-seq) (30). These studies have great resource value and are
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Schematic overview of EndoDB construction and functionalities.</figcaption></figure>

beginning to elucidate EC heterogeneity at the single cell level. However, so far, no general searchable database of EC-specific scRNA-seq data exists, which limits the broad use of these large-scale sequencing efforts. To increase the resource value of scRNA-seq EC datasets, we have developed a dedicated scRNA-seq EndoDB module. For each dataset, we removed low quality and apoptotic cells based on the number of expressed genes (<200) and the percentage of reads assigned to mitochondrial genes (>5%) (31), and performed library size normalization followed by natural log transformation using log1p available via the *Seu-rat* package (32). To *in silico* select the ECs in a dataset, we used an unbiased clustering approach, in which we clustered the cells together based on their expression profile and then assessed which clusters express canonical EC marker genes (*PECAM1*, *CDH5*). We further excluded clusters that expressed markers for leukocytes (*PTPRC*), pericytes (*PDGFRB*) and fibroblasts (*COL1A1*). Using a clustering approach rather than marker gene gating on a cell per cell

basis is more robust, since it allows to overcome misclassification of cells as non-EC due to gene drop out.

## Metadata curation

Metadata was downloaded from GEO or ArrayExpress and manually curated by a team of expert vascular biologists. Repetitive and redundant information was removed, technical information such as the chip design and library preparation was standardized, and each array was manually annotated with biological information such as the model organism, vascular bed, EC type and the experimental conditions (Tables 1–3, Supplementary Table S4 and 5 and Figure 2). In a second step, we used the *OpenRefine* software (33) to standardize spelling across all arrays and datasets. In addition, a description of each study, the GEO or ArrayExpress identifiers and references to the original publication were included in the metadata.
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Most common experimental conditions in the EndoDB. Relative representation of the 15 most common sample types in the EndoDB grouped by experimental conditions. Abbreviations—EC: endothelial cell; VEGF: vascular endothelial growth factor.</figcaption></figure>

<table><caption>**Table 1.** List of species in the EndoDB</caption><thead><tr><th>Species</th><th>Colloquial species name</th><th>Number of datasets in EndoDB</th></tr></thead><tbody><tr><td>Homo sapiens</td><td>Human</td><td>270</td></tr><tr><td>Mus musculus</td><td>Mouse</td><td>76</td></tr><tr><td>Rattus norvegicus</td><td>Rat</td><td>5</td></tr><tr><td>Bos taurus</td><td>Cow</td><td>3</td></tr><tr><td>Danio rerio</td><td>Zebrafish</td><td>3</td></tr><tr><td>Sus scrofa</td><td>Pig</td><td>3</td></tr></tbody></table>

Complete list of species included in the EndoDB, the number of datasets per species is indicated.

<table><caption>**Table 2.** Top 10 most common tissue types in the EndoDB</caption><thead><tr><th>Organ</th><th>Number of datasets in EndoDB</th></tr></thead><tbody><tr><td>Umbilical cord</td><td>149</td></tr><tr><td>Skin</td><td>42</td></tr><tr><td>Aorta</td><td>28</td></tr><tr><td>Lung</td><td>28</td></tr><tr><td>Heart</td><td>25</td></tr><tr><td>Brain</td><td>22</td></tr><tr><td>Liver</td><td>16</td></tr><tr><td>Eye</td><td>15</td></tr><tr><td>Blood</td><td>13</td></tr><tr><td>Lymphatic system</td><td>13</td></tr></tbody></table>

Top 10 most common tissue types in the EndoDB, the number of datasets per tissue type is indicated.

## Data analysis and visualization

Data visualization greatly facilitates interpretation. For all bulk transcriptomics datasets independently, we performed: (i) principal component analysis (PCA) and visualized the results as 2D plots where individual samples are color-coded according to the experimental design; (ii) differential gene expression analysis (DGEA) (34) and visualized the results as volcano plots, bar plots and browsable tables; (iii) gene set enrichment analysis (GSEA) (35,36) and visualized the results as waterfall plots; and (iv) heatmap analysis to show gene expression changes of a panel of user-selected genes. Together, these analyses cover some of the most commonly used analyses and allow for detailed data interpretation. EC single cell data are provided as interactive t-SNE plots using the *Rtsne* and *plotly* packages (see Supplementary Table S3 for scRNA-seq pre-processing and visualization parameters). All analytical results are pre-calculated

<table><caption>**Table 3.** Top 10 most common EC types in the EndoDB</caption><thead><tr><th>Cell type</th><th>Number of datasets in EndoDB</th></tr></thead><tbody><tr><td>Umbilical vein ECs</td><td>142</td></tr><tr><td>Aorta ECs</td><td>29</td></tr><tr><td>Dermal ECs</td><td>23</td></tr><tr><td>Lymphatic ECs</td><td>21</td></tr><tr><td>Coronary artery ECs</td><td>17</td></tr><tr><td>Pulmonary ECs</td><td>16</td></tr><tr><td>Brain ECs</td><td>11</td></tr><tr><td>Retinal ECs</td><td>11</td></tr><tr><td>Dermal lymphatic ECs</td><td>10</td></tr><tr><td>Blood outgrowth ECs</td><td>9</td></tr></tbody></table>

Top 10 most common EC types in the EndoDB, the number of datasets per cell type is indicated.

and stored in the EndoDB to reduce computation time and to allow fast data retrieval.

## Website implementation

The EndoDB is implemented as an interactive web application using the R/Shiny web framework (37,38). We used the *plotly* package for data visualization and the *data.table* and *DT* packages to provide searchable tables (39).

## Data availability and downloads

All pre-processed data and curated metadata are available for download in comma-separated flat file format. Visualizations can be downloaded as PNG or HTML files that preserve the same interactivity as available from the web application, tables can be downloaded in comma-separated value format.

## Documentation and manual

A detailed user manual and video tutorials are available via the web application.

## EndoDB FUNCTIONALITY

The functionality of the EndoDB is centered on four analytical approaches to extract biologically meaningful information from (single cell) transcriptomics data. Below,
we briefly discuss each approach and provide illustrative examples focused on tumor ECs (TECs). TECs are targets of clinically approved anti-angiogenic therapy but remain poorly characterized. Isolating TECs from human or murine tumors is notoriously challenging and laborious, reuse of publicly available data to investigate gene expression reprogramming in TECs is therefore a potentially favorable and cost-effective complementary approach. The examples below are for illustrative purposes only, detailed transcriptome analysis and biological interpretation is beyond the scope of this manuscript.

## Study-centered data exploration

An important purpose of the EndoDB is to facilitate the re-use of publicly available datasets to address novel and unresolved questions, and to allow comparative analyses of similar studies to resolve study-specific biases. The study-centered data exploration options allow users to retrieve studies based on the study title, PubMed identifier, GEO and ArrayExpress identifiers, and manually annotated keywords. After selecting a study, users can download all data and metadata or explore the data using the results from PCA, DGEA, GSEA, heatmap analysis and bar plots.

**Example:** We used the study-centered data exploration functionality to re-analyze a previously published dataset on ECs isolated from normal murine hindbrain, Sonic hedgehog (Shh)-driven or Wnt-driven medulloblastoma (40). PCA showed that normal brain ECs (NECs) have a clearly distinct transcriptional profile compared to Shh- and Wnt-driven medulloblastoma TECs (Figure 3A). To determine which pathways are most upregulated in TECs, we assessed differences in gene expression signatures in NECs versus Shh-medulloblastoma TECs (Shh-TECs). This analysis revealed upregulation of cell cycle and DNA replication (Figure 3B). Consistently, DGEA showed that several proliferation markers (e.g. *Mki67*, *Top2a*) ranked in the top 25 most upregulated genes (Figure 3C-E). Together, these data suggest that Shh-TECs undergo an angiogenic switch and adopt a proliferating phenotype.

## Metabolism- and transcription factor-centered data exploration

Emerging evidence indicates that EC metabolism can overrule fundamental signaling cascades, and several metabolic enzymes and transporters have been identified as therapeutic targets to inhibit pathological angiogenesis in ocular disease, inflammation and cancer (41–43). While the importance of EC metabolism is becoming increasingly clear, still little is known about how metabolism supports EC functions in health and disease (44). Analyses of the subset of metabolic genes and gene sets can provide important insight in metabolic reprogramming (45,46). We therefore performed metabolic gene expression analysis and metabolic GSEA for all studies in the EndoDB. Similarly, to determine which transcription factors are deregulated in response to experimental manipulation, we performed analysis in the subset of genes encoding transcription factors.

**Example:** Performing metabolic GSEA in Shh-TECs versus NECs revealed that pyrimidine metabolism is the most

upregulated metabolic pathway. Consistently, metabolic GSEA and DGEA revealed upregulation of several pathways (pyrimidine metabolism) and genes (*Rrm2*, *Uck2*) involved in nucleotide biosynthesis (Figure 3F and G). Together, these results may indicate that proliferating Shh-TECs upregulate metabolic pathways to sustain nucleotide production required for DNA synthesis. Consistently, the cell cycle related transcription factor *E2f8* was among the most upregulated transcription factors (Figure 3H).

## Gene and pathway-centered data exploration

The EndoDB contains the results of pair-wise DGEA and GSEA results between all experimental conditions within all datasets. The gene- and pathway-centered data exploration functionality allows to unbiasedly determine under which (patho)-physiological and experimental conditions the expression of a particular gene or pathway of interest is most deregulated.

**Examples:** Having identified *Mki67* as a gene of interest in Shh-TECs, we interrogated the EndoDB to determine under which other conditions this gene is upregulated. To do this, we performed a gene-centered search, which showed that *Mki67* is highly upregulated in the brain endothelioma cell line (b.End5) compared with primary cultures of normal brain ECs (Figure 4A and B). DGEA and GSEA analysis further showed that *Mki67*, *Uck2* and DNA replication and pyrimidine metabolism are among the most upregulated genes and pathways in b.End5 endothelioma cells (Figure 4C-E). Together, these findings show that proliferation is a general phenotype of both TECs and malignant ECs.

## Single cell data exploration

An increasing number of studies investigates EC heterogeneity at the single cell level (14,47). scRNA-seq studies are expensive large-scale sequencing efforts that generate high-dimensional datasets with high resource value. We included a specific module in the EndoDB to enable exploration of transcriptional heterogeneity in scRNA-seq datasets via t-SNE plots color-coded for the expression of all detected genes.

**Example:** Based on the study- and gene-centered analyses described above, we hypothesized that ECs in tumors adopt a proliferative phenotype. From bulk transcriptomics analysis alone, it is not possible to determine whether all TECs, or only a subset, are actively proliferating. To determine whether the TEC-associated genes *Mki67*, *Top2a* and *Rrm2* are upregulated in a specific subset of TECs, we interrogated scRNA-seq data of TECs isolated from a colorectal cancer (COLO205) xenograft model (47). Interestingly, all three genes marked a specific subset of TECs suggesting that the observed upregulation of cell cycle-related genes in bulk transcriptomics data results from increased expression in a specific subset of TECs (Figure 5).

Together, these illustrative examples demonstrate the benefit and power of the EndoDB to capitalize on publicly available data to generate novel testable hypotheses in clinically relevant settings.
<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /></figure>

<figure><img src="image_9.png" /></figure>

<figure><img src="image_10.png" /></figure>

<figure><img src="image_11.png" /></figure>

<figure><img src="image_12.png" /></figure>

**Figure 3.** Study-centered data exploration. (A) PCA of normal hindbrain, Shh-medulloblastoma and Wnt-medulloblastoma ECs (study E-GEOD-73753). (B) Gene set enrichment analysis of normal hindbrain ECs versus Shh-medulloblastoma ECs. The upregulated gene sets are shown in red, the downregulated gene sets are shown in blue. (C) Differential analysis of normal hindbrain ECs versus Shh-medulloblastoma ECs shown in volcano plot; some highly deregulated proliferation-associated genes are indicated. (D) Gene expression heatmap of normal hindbrain ECs versus Shh-medulloblastoma and Wnt-medulloblastoma. The high-gene expression levels are shown in red, the low-expression levels in blue. (E) Expression of the indicated genes in normal hindbrain ECs and Shh-medulloblastoma ECs. (F) Metabolic gene set enrichment analysis of normal hindbrain ECs versus Shh-medulloblastoma ECs. The upregulated gene sets are shown in red, the downregulated gene sets are shown in blue. (G and H) Differential analysis of normal hindbrain ECs versus Shh-medulloblastoma ECs for the subset of metabolic genes (G) and transcription factors (H) shown as a volcano plot. Abbreviations—Ccna2: cyclin A2; Mki67: marker of proliferation ki67; NES: normalized enrichment score; Rrm2: ribonucleotide reductase regulatory subunit M2; Shh: Sonic Hedgehog; TEC: tumor endothelial cell; Top2a: topoisomerase 2 alpha; Uck2: uridine-cytidine kinase 2.

DISCUSSION

EndoDB facilitates bench scientists to unlock the untapped potential of publicly available transcriptomics data but can also be used by bioinformaticians as a resource of expert curated data. EndoDB has been constructed using field-tested methods, is tailored to extract biologically meaningful insight from gene expression profiling efforts, covers emerging

fields such as EC metabolism, and is a first-in-class resource of EC scRNA-seq data. EndoDB distinguishes itself from other databases by its flexible user interface, providing an exceptionally large variety of pre-calculated data and analyses in customizable plots and downloadable tables, allowing direct generation and implementation of working hypotheses.

The EndoDB was developed anticipating that an increasing number of EC profiling datasets will become available in
<figure><img src="image_13.png" /><figcaption>**Figure 4.** Gene and pathway-centered data exploration, example of *Mki67* gene search. (A) Dot plot showing log fold change of *Mki67* gene expression in all sample comparisons in EndoDB datasets. The dot indicated by the green arrow head in the enlarged circle represents the pair-wise comparison selected for further analysis (data from study E-GEOD-14375). (B) Expression of the *Mki67* gene in normal hindbrain ECs and endothelioma ECs. (C) Differential analysis of normal brain ECs versus endothelioma ECs for all genes shown in volcano plot. *Mki67* and *Top2a* genes are indicated. (D and E) Gene set enrichment analysis of normal brain ECs versus endothelioma ECs for all genes (D) and for the subset of detected metabolic genes (E). The upregulated gene sets are shown in red, the downregulated gene sets are shown in blue. Abbreviations—*Mki67*: marker of proliferation ki67; NES: normalized enrichment score; *Top2a*: topoisomerase 2 alpha.</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>**Figure 5.** Single cell data exploration. t-SNE plots of COLO205 TECs color-coded for the expression of the indicated proliferation-related genes (GSE110501). Abbreviations—*Mki67*: marker of proliferation ki67; *Rrm2*: ribonucleotide reductase regulatory subunit M2; *Top2a*: topoisomerase 2 alpha.</figcaption></figure>

the near future. The database can be further expanded by including datasets derived from other vascular cells, such as vascular smooth muscle cells and pericytes, but also to include other omics data such as metabolomics and proteomics.

Since ECs have key roles in health and disease, and (in)directly influence all cells in the body, we expect that EndoDB will be a useful tool, for both vascular biologists and the broader biomedical research community, to derive unprecedented insight in EC heterogeneity, identify new (metabolic) therapeutic targets and generate novel hypotheses. EndoDB is freely available at https://vibcancer.be/software-tools/endodb, and will be regularly updated.

## DATA AVAILABILITY

All pre-processed data and curated metadata are available for download in comma-separated flat file format.

## SUPPLEMENTARY DATA

Supplementary Data are available at NAR Online.

## ACKNOWLEDGEMENTS

We acknowledge the help of Mathias Declercq and Valentina Lykhopiy for their feedback on the database design and manuscript preparation.
**Author contributions:** J.G., K.R., S.K. searched primary archives and PubMed to identify EC transcriptomics studies; F.T., J.G., S.K. developed the R scripts for data retrieval and analysis; F.T. implemented the database; A.P., J.G., K.R., L.-C.C., L.dR., L.So., L.T., V.G. curated metadata; J.G., K.R., S.K., clustered single cell data; K.R., L.T., S.K. made the figures; J.G., P.C. conceptualized the EndoDB and supervised the study; J.G., K.R., S.K., P.C., wrote the manuscript; J.G., P.C., X.L., contributed to the execution, support and analysis of experiments, data interpretation and/or advice; G.E., J.K., M.D., L.Sc., T.K. provided advice and discussed results. All authors commented on the manuscript.

## FUNDING

Kom op Tegen Kanker (Stand up to Cancer, Flemish Cancer Society) (to S.K.); 'Fonds voor Wetenschappelijk Onderzoek' (FWO-Vlaanderen) (to K.R., J.K., J.G., P.C.); Marie Curie-IEF Fellowship (to L.T.); Bettencourt Schueller Foundation (to L.T.); Strategisch Basisonderzoek Fonds Wetenschappelijk Onderzoek-Vlaanderen (SB-FWO) (to V.G.); Austrian Science Fund (FWF) [J3730-B26 to A.P.]; Else Kröner-Fresenius-Stiftung Fellowship (to L.-C.C.); Fritz Thyssen Stiftung [10.16.2.017MN to L.-C.C.]; State Key Laboratory of Ophthalmology (to X.L.); Zhongshan Ophthalmic Center at the Sun Yat-Sen University (to X.L.); National Natural Science Foundation of China [81330021, 81670855 to X.L.]; VIB TechWatch program (to P.C.); Flemish Government (to P.C.); Foundation against Cancer [2012–175, 2016–078 to P.C.]; ERC Advanced Research Grant [EU-ERC743074 to P.C.]. Funding for open access charge: ERC Advanced Research Grant [EU-ERC743074 to P.C.].

*Conflict of interest statement.* None declared.

## REFERENCES

1. Aird, W.C. (2012) Endothelial cell heterogeneity. *Cold Spring Harb. Perspect. Med.*, **2**, a006429.

2. Goveia, J., Stapor, P. and Carmeliet, P. (2014) Principles of targeting endothelial cell metabolism to treat angiogenesis and endothelial cell dysfunction in disease. *EMBO Mol. Med.*, **6**, 1105–1120.

3. Carmeliet, P. and Jain, R.K. (2011) Molecular mechanisms and clinical applications of angiogenesis. *Nature*, **473**, 298–307.

4. Heath, J.R., Ribas, A. and Mischel, P.S. (2016) Single-cell analysis tools for drug discovery and development. *Nat. Rev. Drug Discov.*, **15**, 204–216.

5. Parkinson, H., Sarkans, U., Shojatalab, M., Abeygunawardena, N., Contrino, S., Coulson, R., Farne, A., Lara, G.G., Holloway, E., Kapushesky, M. et al. (2005) ArrayExpress—a public repository for microarray gene expression data at the EBI. *Nucleic Acids Res.*, **33**, D553–D555.

6. Barrett, T., Troup, D.B., Wilhite, S.E., Ledoux, P., Rudnev, D., Evangelista, C., Kim, I.F., Soboleva, A., Tomashevsky, M. and Edgar, R. (2007) NCBI GEO: mining tens of millions of expression profiles–database and tools update. *Nucleic Acids Res.*, **35**, D760–D765.

7. Rung, J. and Brazma, A. (2013) Reuse of public genome-wide gene expression data. *Nat. Rev. Genet.*, **14**, 89–99.

8. Rhodes, D.R., Yu, J.J., Shanker, K., Deshpande, N., Varambally, R., Ghosh, D., Barrette, T., Pandey, A. and Chinnaiyan, A.M. (2004) ONCOMINE: a cancer microarray database and integrated data-mining platform. *Neoplasia*, **6**, 1–6.

9. Tomlins, S.A., Rhodes, D.R., Perner, S., Dhanasekaran, S.M., Mehra, R., Sun, X.W., Varambally, S., Cao, X., Tchinda, J., Kuefer, R.

*et al.* (2005) Recurrent fusion of TMPRSS2 and ETS transcription factor genes in prostate cancer. *Science*, **310**, 644–648.

10. Bos, P.D., Zhang, X.H., Nadal, C., Shu, W., Gomis, R.R., Nguyen, D.X., Minn, A.J., van de Vijver, M.J., Gerald, W.L., Foekens, J.A. *et al.* (2009) Genes that mediate breast cancer metastasis to the brain. *Nature*, **459**, 1005–1009.

11. Qian, B.Z., Li, J., Zhang, H., Kitamura, T., Zhang, J., Campion, L.R., Kaiser, E.A., Snyder, L.A. and Pollard, J.W. (2011) CCL2 recruits inflammatory monocytes to facilitate breast-tumour metastasis. *Nature*, **475**, 222–225.

12. Tomlins, S.A., Mehra, R., Rhodes, D.R., Cao, X., Wang, L., Dhanasekaran, S.M., Kalyana-Sundaram, S., Wei, J.T., Rubin, M.A., Pienta, K.J. *et al.* (2007) Integrative molecular concept modeling of prostate cancer progression. *Nat. Genet.*, **39**, 41–51.

13. Erez, N., Truitt, M., Olson, P., Arron, S.T. and Hanahan, D. (2010) Cancer-associated fibroblasts are activated in incipient Neoplasia to orchestrate tumor-promoting inflammation in an NF-kappaB-Dependent manner. *Cancer Cell*, **17**, 135–147.

14. Vanlandewijck, M., He, L., Mae, M.A., Andrae, J., Ando, K., Del Gaudio, F., Nahar, K., Lebouvier, T., Lavina, B., Gouveia, L. *et al.* (2018) A molecular atlas of cell types and zonation in the brain vasculature. *Nature*, **554**, 475–480.

15. Salvucci, O., Ohnuki, H., Maric, D., Hou, X., Li, X., Yoon, S.O., Segarra, M., Eberhart, C.G., Acker-Palmer, A. and Tosato, G. (2015) EphrinB2 controls vessel pruning through STAT1-JNK3 signalling. *Nat. Commun.*, **6**, 6576.

16. Vandekeere, S., Dubois, C., Kalucka, J., Sullivan, M.R., Garcia-Caballero, M., Goveia, J., Chen, R., Diehl, F.F., Bar-Lev, L., Souffreaux, J. *et al.* (2018) Serine synthesis via PHGDH is essential for Heme production in endothelial cells. *Cell Metab.*, **28**, 573–587.

17. Bruning, U., Morales-Rodriguez, F., Kalucka, J., Goveia, J., Taverna, F., Queiroz, K.C.S., Dubois, C., Cantelmo, A.R., Chen, R., Loroch, S. *et al.* (2018) Impairment of angiogenesis by fatty acid synthase inhibition involves mTOR Malonylation. *Cell Metab.*, doi:10.1016/j.cmet.2018.07.019.

18. Singer, B.D., Mock, J.R., D'Alessio, F.R., Aggarwal, N.R., Mandke, P., Johnston, L. and Damarla, M. (2016) Flow-cytometric method for simultaneous analysis of mouse lung epithelial, endothelial, and hematopoietic lineage cells. *Am. J. Physiol. Lung Cell. Mol. Physiol.*, **310**, L796–L801.

19. Sedlacek, H.H. (2001) Pharmacological aspects of targeting cancer gene therapy to endothelial cells. *Crit. Rev. Oncol. Hematol.*, **37**, 169–215.

20. Lambrechts, D., Wauters, E., Boeckx, B., Aibar, S., Nittner, D., Burton, O., Bassez, A., Decaluwe, H., Pircher, A., Van den Eynde, K. *et al.* (2018) Phenotype molding of stromal cells in the lung tumor microenvironment. *Nat. Med.*, **24**, 1277–1289.

21. van den Brink, S.C., Sage, F., Vertese, A., Spanjaard, B., Peterson-Maduro, J., Baron, C.S., Robin, C. and van Oudenaarden, A. (2017) Single-cell sequencing reveals dissociation-induced gene expression in tissue subpopulations. *Nat. Methods*, **14**, 935–936.

22. Ziegenhain, C., Vieth, B., Parekh, S., Reinius, B., Guillaumet-Adkins, A., Smets, M., Leonhardt, H., Heyn, H., Hellmann, I. and Enard, W. (2017) Comparative analysis of Single-cell RNA sequencing methods. *Mol. Cell*, **65**, 631–643.

23. Bacher, R. and Kendziorski, C. (2016) Design and computational analysis of single-cell RNA-sequencing experiments. *Genome Biol.*, **17**, 63.

24. Bolstad, B.M., Irizarry, R.A., Astrand, M. and Speed, T.P. (2003) A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. *Bioinformatics*, **19**, 185–193.

25. Irizarry, R.A., Hobbs, B., Collin, F., Beazer-Barclay, Y.D., Antonellis, K.J., Scherf, U. and Speed, T.P. (2003) Exploration, normalization, and summaries of high density oligonucleotide array probe level data. *Biostatistics*, **4**, 249–264.

26. Ritchie, M.E., Phipson, B., Wu, D., Hu, Y., Law, C.W., Shi, W. and Smyth, G.K. (2015) limma powers differential expression analyses for RNA-sequencing and microarray studies. *Nucleic Acids Res.*, **43**, e47.

27. Kauffmann, A., Gentleman, R. and Huber, W. (2009) arrayQualityMetrics—a bioconductor package for quality assessment of microarray data. *Bioinformatics*, **25**, 415–416.
28. Bray,N.L., Pimentel,H., Melsted,P. and Pachter,L. (2016) Near-optimal probabilistic RNA-seq quantification. *Nat. Biotechnol.*, **34**, 525–527.

29. Robinson,M.D., McCarthy,D.J. and Smyth,G.K. (2010) edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. *Bioinformatics*, **26**, 139–140.

30. Linnarsson,S. and Teichmann,S.A. (2016) Single-cell genomics: coming of age. *Genome Biol.*, **17**, 97.

31. Lun,A.T., McCarthy,D.J. and Marioni,J.C. (2016) A step-by-step workflow for low-level analysis of single-cell RNA-seq data with Bioconductor [version 2; referees: 3 approved, 2 approved with reservations]. *F1000 Res.*, **5**, 2122.

32. Satija,R., Farrell,J.A., Gennert,D., Schier,A.F. and Regev,A. (2015) Spatial reconstruction of single-cell gene expression data. *Nat. Biotechnol.*, **33**, 495–502.

33. Ham,K. (2013) OpenRefine (version 2.5). Free, open-source tool for cleaning and transforming data. *J. Med. Libr. Assoc.*, **101**, 233–234.

34. Smyth,G.K. (2004) Linear models and empirical bayes methods for assessing differential expression in microarray experiments. *Stat. Appl. Genet. Mol. Biol.*, **3**, doi:10.2202/1544-6115.1027.

35. Subramanian,A., Tamayo,P., Mootha,V.K., Mukherjee,S., Ebert,B.L., Gillette,M.A., Paulovich,A., Pomeroy,S.L., Golub,T.R., Lander,E.S. et al. (2005) Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. *PNAS*, **102**, 15545–15550.

36. Yu,G., Wang,L.G., Han,Y. and He,Q.Y. (2012) clusterProfiler: an R package for comparing biological themes among gene clusters. *OMICS*, **16**, 284–287.

37. R Development Core Team (2018) R: A Language and Environment for Statistical Computing. R Foundation for Statistical Computing. Vienna.

38. Chang,W., Cheng,J., Allaire,J.J., Xie,Y. and McPherson,J. (2018) shiny: Web Application Framework for R.

39. P.T. Inc. (2015) Plotly Technologies Inc., Montréal.

40. Phoenix,T.N., Patmore,D.M., Boop,S., Boulos,N., Jacus,M.O., Patel,Y.T., Roussel,M.F., Finkelstein,D., Goumnerova,L., Perreault,S. et al. (2016) Medulloblastoma genotype dictates blood brain barrier phenotype. *Cancer Cell*, **29**, 508–522.

41. Cantelmo,A.R., Conradi,L.C., Brajic,A., Goveia,J., Kalucka,J., Pircher,A., Chaturvedi,P., Hol,J., Thienpont,B., Teuwen,L.A. et al. (2016) Inhibition of the glycolytic activator PFKFB3 in endothelium induces tumor vessel normalization, impairs metastasis, and improves chemotherapy. *Cancer Cell*, **30**, 968–985.

42. Schoors,S., Bruning,U., Missiaen,R., Queiroz,K.C.S., Borgers,G., Elia,I., Zecchin,A., Cantelmo,A.R., Christen,S., Goveia,J. et al. (2015) Fatty acid carbon is essential for dNTP synthesis in endothelial cells. *Nature*, **520**, 192–197.

43. Schoors,S., Cantelmo,A.R., Georgiadou,M., Stapor,P., Wang,X., Quaegebeur,A., Cauwenberghs,S., Wong,B.W., Bifari,F., Decimo,I. et al. (2014) Incomplete and transitory decrease of glycolysis: a new paradigm for anti-angiogenic therapy? *Cell Cycle*, **13**, 16–22.

44. Rohlenova,K., Veys,K., Miranda-Santos,I., De Bock,K. and Carmeliet,P. (2018) Endothelial cell metabolism in health and disease. *Trends Cell Biol.*, **28**, 224–236.

45. Nilsson,R., Jain,M., Madhusudhan,N., Sheppard,N.G., Strittmatter,L., Kampf,C., Huang,J., Asplund,A. and Mootha,V.K. (2014) Metabolic enzyme expression highlights a key role for MTHFD2 and the mitochondrial folate pathway in cancer. *Nat. Commun.*, **5**, 3128.

46. Hu,J., Locasale,J.W., Bielas,J.H., O'Sullivan,J., Sheahan,K., Cantley,L.C., Vander Heiden,M.G. and Vitkup,D. (2013) Heterogeneity of tumor-induced gene expression changes in the human metabolic network. *Nat. Biotechnol.*, **31**, 522–529.

47. Zhao,Q., Eichten,A., Parveen,A., Adler,C., Huang,Y., Wang,W., Ding,Y., Adler,A., Nevins,T., Ni,M. et al. (2018) Single-cell transcriptome analyses reveal endothelial cell heterogeneity in tumors and changes following antiangiogenic treatment. *Cancer Res.*, **78**, 2370–2382.