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.

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

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

The murine retinal pigment epithelium requires peroxisomal ß-oxidation to maintain lysosomal function and prevent dedifferentiation.

Retinal pigment epithelium (RPE) cells have to phagocytose shed photoreceptor outer segments (POS) on a daily basis over the lifetime of an organism, but the mechanisms involved in the digestion and recycling of POS lipids are poorly understood. Although it was frequently assumed that peroxisomes may play an essential role, this was never investigated. Here, we show that global as well as RPE-selective loss of peroxisomal ß-oxidation in multifunctional protein 2 (MFP2) knockout mice impairs the digestive function of lysosomes in the RPE at a very early age, followed by RPE degeneration. This was accompanied by prolonged mammalian target of rapamycin activation, lipid deregulation, and mitochondrial structural anomalies without, however, causing oxidative stress or energy shortage. The RPE degeneration caused secondary photoreceptor death. Notably, the deterioration of the RPE did not occur in an Mfp2/rd1 mutant mouse line, characterized by absent POS shedding. Our findings prove that peroxisomal ß-oxidation in the RPE is essential for handling the polyunsaturated fatty acids present in ingested POS and shed light on retinopathy in patients with peroxisomal disorders. Our data also have implications for gene therapy development as they highlight the importance of targeting the RPE in addition to the photoreceptor cells.

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.

PARP inhibition promotes endothelial-like traits in melanoma cells and modulates pericyte coverage dynamics during vasculogenic mimicry.

Vasculogenic mimicry (VM) describes the ability of highly aggressive tumor cells to develop pseudovascular structures without the participation of endothelial cells. PARP1 is implicated in the activation of hypoxia-inducible factors, which are crucial in tumor neovascularization. We have explored the role of hypoxia and PARP inhibition in VM. In uveal melanoma xenografts, the PARP inhibitor olaparib improved in vivo pericyte coverage specifically of VM channels. This was concomitant with reduced metastasis in olaparib-treated VM+ tumors. PARP inhibition and hypoxia modulated melanoma tube formation in vitro, inducing a more sparse and regular tubular architecture. Whole-transcriptome profiling revealed that olaparib treatment under hypoxic conditions modulated the expression of genes implicated in vasculogenesis during tube formation, enhancing the endothelial-like phenotype of VM+ uveal melanoma cells. PARP inhibition, especially during hypoxia, upregulated PDGFß, which is essential for pericyte recruitment. Our study indicates that PARP inhibitors may enhance the endothelial characteristics of VM+ cells, modulate pericyte coverage, and reduce metastatic spread in VM+ melanoma. © 2022 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

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.

Lipid droplet degradation by autophagy connects mitochondria metabolism to Prox1-driven expression of lymphatic genes and lymphangiogenesis.

Autophagy has vasculoprotective roles, but whether and how it regulates lymphatic endothelial cells (LEC) homeostasis and lymphangiogenesis is unknown. Here, we show that genetic deficiency of autophagy in LEC impairs responses to VEGF-C and injury-driven corneal lymphangiogenesis. Autophagy loss in LEC compromises the expression of main effectors of LEC identity, like VEGFR3, affects mitochondrial dynamics and causes an accumulation of lipid droplets (LDs) in vitro and in vivo. When lipophagy is impaired, mitochondrial ATP production, fatty acid oxidation, acetyl-CoA/CoA ratio and expression of lymphangiogenic PROX1 target genes are dwindled. Enforcing mitochondria fusion by silencing dynamin-related-protein 1 (DRP1) in autophagy-deficient LEC fails to restore LDs turnover and lymphatic gene expression, whereas supplementing the fatty acid precursor acetate rescues VEGFR3 levels and signaling, and lymphangiogenesis in LEC-Atg5-/- mice. Our findings reveal that lipophagy in LEC by supporting FAO, preserves a mitochondrial-PROX1 gene expression circuit that safeguards LEC responsiveness to lymphangiogenic mediators and lymphangiogenesis.

Cell Type-Selective Loss of Peroxisomal ß-Oxidation Impairs Bipolar Cell but Not Photoreceptor Survival in the Retina.

Retinal degeneration is a common feature in peroxisomal disorders leading to blindness. Peroxisomes are present in the different cell types of the retina; however, their precise contribution to retinal integrity is still unclear. We previously showed that mice lacking the central peroxisomal ß-oxidation enzyme, multifunctional protein 2 (MFP2), develop an early onset retinal decay including photoreceptor cell death. To decipher the function of peroxisomal ß-oxidation in photoreceptors, we generated cell type selective Mfp2 knockout mice, using the Crx promotor targeting photoreceptors and bipolar cells. Surprisingly, Crx-Mfp2-/- mice maintained photoreceptor length and number until the age of 1 year. A negative electroretinogram was indicative of preserved photoreceptor phototransduction, but impaired downstream bipolar cell signaling from the age of 6 months. The photoreceptor ribbon synapse was affected, containing free-floating ribbons and vesicles with altered size and density. The bipolar cell interneurons sprouted into the ONL and died. Whereas docosahexaenoic acid levels were normal in the neural retina, levels of lipids containing very long chain polyunsaturated fatty acids were highly increased. Crx-Pex5-/- mice, in which all peroxisomal functions are inactivated in photoreceptors and bipolar cells, developed the same phenotype as Crx-Mfp2-/- mice. In conclusion, the early photoreceptor death in global Mfp2-/- mice is not driven cell autonomously. However, peroxisomal ß-oxidation is essential for the integrity of photoreceptor ribbon synapses and of bipolar cells.

Endothelial Zeb2 preserves the hepatic angioarchitecture and protects against liver fibrosis.

AIMS: Hepatic capillaries are lined with specialized liver sinusoidal endothelial cells (LSECs) which support macromolecule passage to hepatocytes and prevent fibrosis by keeping hepatic stellate cells (HSCs) quiescent. LSEC specialization is co-determined by transcription factors. The zinc-finger E-box-binding homeobox (Zeb)2 transcription factor is enriched in LSECs. Here, we aimed to elucidate the endothelium-specific role of Zeb2 during maintenance of the liver and in liver fibrosis. METHODS AND RESULTS: To study the role of Zeb2 in liver endothelium we generated EC-specific Zeb2 knock-out (ECKO) mice. Sequencing of liver EC RNA revealed that deficiency of Zeb2 results in prominent expression changes in angiogenesis-related genes. Accordingly, the vascular area was expanded and the presence of pillars inside ECKO liver vessels indicated that this was likely due to increased intussusceptive angiogenesis. LSEC marker expression was not profoundly affected and fenestrations were preserved upon Zeb2 deficiency. However, an increase in continuous EC markers suggested that Zeb2-deficient LSECs are more prone to dedifferentiation, a process called 'capillarization'. Changes in the endothelial expression of ligands that may be involved in HSC quiescence together with significant changes in the expression profile of HSCs showed that Zeb2 regulates LSEC-HSC communication and HSC activation. Accordingly, upon exposure to the hepatotoxin carbon tetrachloride (CCl4), livers of ECKO mice showed increased capillarization, HSC activation, and fibrosis compared to livers from wild-type littermates. The vascular maintenance and anti-fibrotic role of endothelial Zeb2 was confirmed in mice with EC-specific overexpression of Zeb2, as the latter resulted in reduced vascularity and attenuated CCl4-induced liver fibrosis. CONCLUSION: Endothelial Zeb2 preserves liver angioarchitecture and protects against liver fibrosis. Zeb2 and Zeb2-dependent genes in liver ECs may be exploited to design novel therapeutic strategies to attenuate hepatic fibrosis.

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.

AMPKα1 deletion in myofibroblasts exacerbates post-myocardial infarction fibrosis by a connexin 43 mechanism.

We have previously demonstrated that systemic AMP-activated protein kinase α1 (AMPKα1) invalidation enhanced adverse LV remodelling by increasing fibroblast proliferation, while myodifferentiation and scar maturation were impaired. We thus hypothesised that fibroblastic AMPKα1 was a key signalling element in regulating fibrosis in the infarcted myocardium and an attractive target for therapeutic intervention. The present study investigates the effects of myofibroblast (MF)-specific deletion of AMPKα1 on left ventricular (LV) adaptation following myocardial infarction (MI), and the underlying molecular mechanisms. MF-restricted AMPKα1 conditional knockout (cKO) mice were subjected to permanent ligation of the left anterior descending coronary artery. cKO hearts exhibit exacerbated post-MI adverse LV remodelling and are characterised by exaggerated fibrotic response, compared to wild-type (WT) hearts. Cardiac fibroblast proliferation and MF content significantly increase in cKO infarcted hearts, coincident with a significant reduction of connexin 43 (Cx43) expression in MFs. Mechanistically, AMPKα1 influences Cx43 expression by both a transcriptional and a post-transcriptional mechanism involving miR-125b-5p. Collectively, our data demonstrate that MF-AMPKα1 functions as a master regulator of cardiac fibrosis and remodelling and might constitute a novel potential target for pharmacological anti-fibrotic applications.

Peroxisomal Multifunctional Protein 2 Deficiency Perturbs Lipid Homeostasis in the Retina and Causes Visual Dysfunction in Mice.

Patients lacking multifunctional protein 2 (MFP2), the central enzyme of the peroxisomal ß-oxidation pathway, develop retinopathy. This pathway is involved in the metabolism of very long chain (VLCFAs) and polyunsaturated (PUFAs) fatty acids, which are enriched in the photoreceptor outer segments (POS). The molecular mechanisms underlying the retinopathy remain, however, elusive. Here, we report that mice with MFP2 inactivation display decreased retinal function already at the age of 3 weeks, which is accompanied by a profound shortening of the photoreceptor outer and inner segments, but with preserved photoreceptor ultrastructure. Furthermore, MFP2 deficient retinas exhibit severe changes in gene expression with downregulation of genes involved in the phototransduction pathway and upregulation of inflammation related genes. Lipid profiling of the mutant retinas revealed a profound reduction of DHA-containing phospholipids. This was likely due to a hampered systemic supply and retinal traffic of this PUFA, although we cannot exclude that the local defect of peroxisomal ß-oxidation contributes to this DHA decrease. Moreover, very long chain PUFAs were also reduced, with the exception of those containing ≥ 34 carbons that accumulated. The latter suggests that there is an uncontrollable elongation of retinal PUFAs. In conclusion, our data reveal that intact peroxisomal ß-oxidation is indispensable for retinal integrity, most likely by maintaining PUFA homeostasis.

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.

Tightening the retinal glia limitans attenuates neuroinflammation after optic nerve injury.

Increasing evidence suggests that functional impairments at the level of the neurovascular unit (NVU) underlie many neurodegenerative and neuroinflammatory diseases. While being part of the NVU, astrocytes have been largely overlooked in this context and only recently, tightening of the glia limitans has been put forward as an important neuroprotective response to limit these injurious processes. In this study, using the retina as a central nervous system (CNS) model organ, we investigated the structure and function of the glia limitans, and reveal that the blood-retina barrier and glia limitans function as a coordinated double barrier to limit infiltration of leukocytes and immune molecules. We provide in vitro and in vivo evidence for a protective response at the NVU upon CNS injury, which evokes inflammation-induced glia limitans tightening. Matrix metalloproteinase-3 (MMP-3) was found to be a crucial regulator of this process, thereby revealing its beneficial and immunomodulatory role in the CNS. in vivo experiments in which MMP-3 activity was deleted via genetic and pharmacological approaches, combined with a comprehensive study of tight junction molecules, glial end feet markers, myeloid cell infiltration, cytokine expression and neurodegeneration, show that MMP-3 attenuates neuroinflammation and neurodegeneration by tightening the glia limitans, thereby pointing to a prominent role of MMP-3 in preserving the integrity of the NVU upon injury. Finally, we gathered promising evidence to suggest that IL1b, which is also regulated by MMP-3, is at least one of the molecular messengers that induces glia limitans tightening in the injured CNS.

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.

Unraveling the transcriptional determinants of liver sinusoidal endothelial cell specialization.

Liver sinusoidal endothelial cells (LSECs) are the first liver cells to encounter waste macromolecules, pathogens, and toxins in blood. LSECs are highly specialized to mediate the clearance of these substances via endocytic scavenger receptors and are equipped with fenestrae that mediate the passage of macromolecules toward hepatocytes. Although some transcription factors (TFs) are known to play a role in LSEC specialization, information about the specialized LSEC signature and its transcriptional determinants remains incomplete.Based on a comparison of liver, heart, and brain endothelial cells (ECs), we established a 30-gene LSEC signature comprising both established and newly identified markers, including 7 genes encoding TFs. To evaluate the LSEC TF regulatory network, we artificially increased the expression of the 7 LSEC-specific TFs in human umbilical vein ECs. Although Zinc finger E-box-binding protein 2, homeobox B5, Cut-like homolog 2, and transcription factor EC (TCFEC) had limited contributions, musculoaponeurotic fibrosarcoma (C-MAF), GATA binding protein 4 (GATA4), and MEIS homeobox 2 (MEIS2) emerged as stronger inducers of LSEC marker expression. Furthermore, a combination of C-MAF, GATA4, and MEIS2 showed a synergistic effect on the increase of LSEC signature genes, including liver/lymph node-specific ICAM-3 grabbing non-integrin (L-SIGN) (or C-type lectin domain family member M (CLEC4M)), mannose receptor C-Type 1 (MRC1), legumain (LGMN), G protein-coupled receptor 182 (GPR182), Plexin C1 (PLXNC1), and solute carrier organic anion transporter family member 2A1 (SLCO2A1). Accordingly, L-SIGN, MRC1, pro-LGMN, GPR182, PLXNC1, and SLCO2A1 protein levels were elevated by this combined overexpression. Although receptor-mediated endocytosis was not significantly induced by the triple TF combination, it enhanced binding to E2, the hepatitis C virus host-binding protein. We conclude that C-MAF, GATA4, and MEIS2 are important transcriptional regulators of the unique LSEC fingerprint and LSEC interaction with viruses. Additional factors are however required to fully recapitulate the molecular, morphological, and functional LSEC fingerprint.NEW & NOTEWORTHY Liver sinusoidal endothelial cells (LSECs) are the first liver cells to encounter waste macromolecules, pathogens, and toxins in the blood and are highly specialized. Although some transcription factors are known to play a role in LSEC specialization, information about the specialized LSEC signature and its transcriptional determinants remains incomplete. Here, we show that Musculoaponeurotic Fibrosarcoma (C-MAF), GATA binding protein 4 (GATA4), and Meis homeobox 2 (MEIS2) are important transcriptional regulators of the unique LSEC signature and that they affect the interaction of LSECs with viruses.