Single-Cell RNA Sequencing in Lymphatic Endothelial Cells: Comparison
Please note this is a comparison between Version 1 by Astrid Charlotte van Stigt and Version 2 by Vivi Li.

Lymphatic endothelial cells (LECs) line the lymphatic vasculature and play a central role in the immune response. LECs have abilities to regulate immune transport, to promote immune cell survival, and to cross present antigens to dendritic cells. Single-cell RNA sequencing (scRNA) technology has accelerated new discoveries in the field of lymphatic vascular biology. 

  • lymphatics
  • lymphatic endothelial cells
  • single cell RNA-sequencing

1. Introduction

The lymphatic vasculature plays an essential role in many physiological and pathological processes in the human body. Although the origin of the lymphatic vasculature is closely connected to the vascular system, it has unique tasks and properties. Amongst these are the regulation of interstitial fluid balance, clearance of tissue-derived debris, and lipids and immune cell transport. Thus far, tissue fluid homeostasis has been the best studied function of the lymphatic vasculature, but novel roles are being increasingly discovered [1][2][3][1,2,3]. Tissue inflammation or injury, featured by increased interstitial fluid and enhanced influx of immune cells, requires expansion of the lymphatic vasculature to facilitate transport of tissue debris, leukocytes, and antigens (Ag), either soluble or carried by Ag-presenting cells, to the lymph node (LN). These events lead to an efficient immune response, and, ultimately, resolution of inflammation. In the last 20 years, identification of lymphatic markers, new experimental in vivo models, and 2D and 3D culture assays with lymphatic endothelial cells (LECs) have advanced our knowledge on lymphatic vascular biology [4]. This has important implications for a wide spectrum of disease areas, such as wound healing, obesity, cancer, and chronic inflammation and autoimmunity [5][6][7][8][5,6,7,8]. Changes in lymphatic vasculature, such as dilated lymphatic vessels or decreased amount of lymphatic vessels, are also widely observed in chronic (rheumatic) inflammatory diseases, such as psoriasis, rheumatoid arthritis (RA), systemic sclerosis, and inflammatory bowel disease (IBD) [9][10][11][12][9,10,11,12].
The recently developed technique of single-cell RNA sequencing (scRNA-seq) is widely used to uncover the identity and heterogeneity of cell subpopulations, and it has, either alone or combined with other state-of-the-art techniques, increased our knowledge on new functions of LECs in health and disease. The scRNA-seq technique evolved from sequencing only a handful of cells to currently sequencing more than 100,000 cells per sample [13]. ScRNA-seq protocols usually consist of the following steps: single cell capture, single cell lysis, reverse transcription, preamplification, library preparation, and sequencing. Laser capture microdissection, fluorescence-activated cell sorting (FACS) in microtiter plates, microfluidics, and microdroplets can be used to capture single cells. Techniques can differ in the number of cells that can be included, sequencing depth, and costs [14]. High dimensional data analysis has been applied to obtain biologically meaningful results. Frequently used strategies include dimensionality reduction methods, such as principal component analysis (PCA), t-distributed stochastic neighbor embedding (t-SNE), and uniform manifold approximation and projection (UMAP), followed by clustering algorithms, such as hierarchical clustering or k-means clustering [15]. Other promising approaches, such as trajectory inference and cell–cell interaction plots, intend to predict the origin of cells and interactions with other cells [16][17][16,17].

2. Single-Cell Transcriptomic Insides in Embryonic Development

Studies on embryonic development can often explain complex pathological processes, as has been demonstrated by Burchill et al. using scRNA-seq of liver LECs in an inflammatory environment [18].

2.1. Peripheral Lymphatic Vessel Formation

In mammals, both the origin of LECs lining the lymphatic vessels and the underlying molecular mechanisms and proteins involved in the embryonic development are not entirely clear. Historically, LECs were thought to exclusively originate from embryonic venous cells [19][20][21][22][19,20,21,22]. However, more recent studies also found that specific lymphatic vascular beds in the skin, the gut, and the heart can bear LECs from non-venous origins [23][24][25][26][23,24,25,26]. Molecular heterogeneity between LECs isolated from different organs also supports the possibility of multiple sources for LECs [27][28][29][27,28,29]. Unfortunately, the scRNA-seq technique has not yet been used to study the origin of LECs in the different lymphatic vascular beds.
The underlying molecular mechanisms and contributors of the venous-to-lymphatic transition have been best studied. This process starts with the transcription factor Nuclear Receptor Subfamily2 Group F Member 2 (NR2F2), also called COUP-TFII, cooperating with SRY-Box Transcription factor 18 (SOX18) to activate Prospero homeobox protein 1 (PROX1) expression [30]. PROX1 is the master transcription factor in LECs and drives expression of key markers, such as vascular endothelial growth factor receptor 3 (VEGFR-3), integrin-9alpha (adhesion molecule), podoplanin (cell migration), and chemokine (C-C motif) ligand 21 (CCL21) (chemokine) [20][31][32][33][20,31,32,33]. Downregulation of Prox1 and its target genes result in loss of LEC identity and function, as was observed in a scRNA-seq mouse liver study [18]. According to this study, LEC identity and function could be regained by exposing them to VEGF-C. VEGF-C is one of the five members of the VEGF family of angiogenic cytokines that signals through VEGFR-3, which, in turn, signals via the growth factor-regulated phosphoinositide 3-kinase (PI3K)–AKT signaling network or the Mitogen-activated protein kinase (MAPK)/extracellular signal-regulated kinases (ERK) pathways to regulate embryonic development of LECs and LEC identity, survival, proliferation, and migration in adulthood [34].
Based on data of scRNA-seq studies, new genes have been suggested to play a role in embryonic lymphatic development. A recent scRNA-seq study in a murine newborn lung model predicted, aside from Prox-1, Sox18, and Nr2f2, for example, homeobox d 8 (Hoxd8), Maf, T-box 1 (Tbx1), and ETS Transcription Factor 3 (Elk3) to be involved in LEC development [35]Hoxd8 is a homeodomain transcription factor and can induce and maintain Prox-1 expression [30]Maf is a transcription factor upregulated by VEGF-C/VEGFR3 signaling that plays a role in developmental lymphangiogenesis [36]Tbx1 belongs to the T-box transcription factors, which are mostly involved in developmental processes, and Tbx1 mouse mutants have severe problems with pharyngeal and cardiovascular development. Interestingly, conditional deletion of Tbx1 in mice also causes widespread lymphangiogenesis defects due to its role in regulating VEGFR3 expression [37]Elk3 is an Erythroblast Transformation Specific (ETS) binding domain transcription factor and has not yet been studied in the context of embryonic development. However, ELK3-expressing LECs promote breast cancer growth and migration of cancer cells [38]. As published earlier, NR2F2 can form heterodimers with PROX-1 to activate Notch signaling in LECs important for lymphangiogenesis [39]. Other unknown protein-coding genes found in the scRNA-seq in the newborn lung study include, for example, Profilin 1 (Pfn1) and AE Binding Protein 1(Aebp1) [35]. Future studies into these genes might confirm roles for these factors in embryonic development.
Strikingly, a recent study combining an scRNA-seq approach validated by experimental mouse models discovered that the protein Folliculin (FLCN), long known as a tumor suppressor protein, plays a major role in controlling venous-to-lymphatic cell transition. FLCN prevents expression of Prox-1 in venous endothelial cells (VECs) by binding of the Prox-1 regulatory element, the transcription factor E3 (Tfe3). FLCN deficiency results in disrupted lymphatic and venous development and dilatation of vessels [40]Figure 1A shows an overview of venous-to-lymphatic transdifferentiation. Many processes in embryonic development driven by these transcription factors have been implicated in post-natal lymphangiogenesis during inflammation. However, knowledge on aberrant expression of many of these transcription factors is lacking.
Figure 1. Lymphatic endothelial cell and lymph node embryonic development. (A) In embryonic development, lymphatic endothelial cells (LECs) mostly transdifferentiate from venous endothelial cells (VECs). Prospero homeobox protein 1 (Prox-1) is a master regulator of lymphatic fate. Homeobox d 8 (Hoxd8) can activate Prox-1 expression, and Prox-1 cooperates with Nuclear Receptor Subfamily2 Group F Member 2 (Nr2f2) to drive expression of other key peripheral LEC molecules, such as podoplanin(Pdpn), lymphatic vessel endothelial hyaluron receptor 1 (Lyve-1), chemokine (C-C motif) ligand 21 (Ccl21), and vascular endothelial growth factor receptor 3 (Vegfr3). Vegfr3 is also regulated by transcription factor T-box 1 (Tbx1) and can activate other transcription factors, such as Maf. Single cell RNA sequencing (scRNA-seq) data also suggested transcription factors, such as ETS Transcription Factor 3 (Elk3), or other protein coding RNAs, such as Profilin 1 (Pfn1) and AE binding protein 1 (Aebp1), to be involved in lymphatic development, but their function is not yet clear. (B) Lymph node (LN) development follows several steps: (1) mesenchymal cells secrete chemokine (C-X-C motif) ligand 13 (CXCL13) to attract lymphoid tissue inducer (LTi) cells. (2) LTi cells cluster together and activate receptor activator nuclear factor κ B (RANK) and RANK ligand (RANKL) signaling. (3) RANK-RANKL signaling increases lymphotoxin α1β2 (LT1β2), which activates lymphotoxin beta receptor (LTβR) on mesenchymal cells responding with more cytokine signaling and attracting more LTi cells. (4) Mesenchymal cells differentiate into stromal cell subsets of the LN, and LTi differentiate or attract LECs. Trajectory inference of scRNA-seq data showed that there are two developmental paths of the LN LECs.

2.2. Lymph Node Development

The embryonic development of the lymph node (LN) is a step-by-step process and requires close interaction between mesenchymal cells and lymphoid tissue inducer (LTi) cells (Figure 1B). First, mesenchymal cells recruit LTi cells by expressing chemokine (C-X-C motif) ligand 13 (CXCL13) (Figure 1B, step 1). LTi cells cluster together and are identified by expression of chemokine (C-X-C motif) receptor 5 (CXCR5), the receptor for CXCL13, CD45, CD4, and CD127 (Figure 1B, step 2). LTi cells also express receptor activator nuclear factor κ B (RANK) and its ligand, RANKL, which induces high amounts of lymphotoxin α1β2 (LT1β2) [41][42][43][41,42,43]. LT1β2 activates the lymphtoxin beta receptor (LTβR) on mesenchymal cells. Subsequently, LTβR induces expression of cell adhesion molecules, such as intracellular adhesion molecule 1 (ICAM-1), and chemokines, such as CXCL13, CCL21, and CCL19, by mesenchymal cells. This initiates a positive feedback loop that attracts more LTi cells to the site (Figure 1B, step 3). The mesenchymal cells differentiate eventually in the stromal cells found in the lymph node. Whether LTi cells differentiate into LECs or, at their turn, attract LECs to the LN site is currently unknown [44]. Developed LNs consist of a subcapular sinus (SCS), which is connected to the afferent lymphatic vessels, and they collect migrating immune cells from draining tissues, including T- and B-lymphocytes and dendritic cells (DCs) [8]. Resident CD169+ macrophages and DCs in the SCS of the LN provide Ags, assembled from captured pathogens and molecules/proteins, to passing lymphocytes [45][46][47][45,46,47]. This is an essential interaction for educating lymphocytes towards proper immune responses. Lymphocytes are activated in LN parenchyma and travel through the cortical sinus (CS) and medullary sinus (MS) before they exit the LN into the peripheral blood [48].
Trajectory inference analysis of lymph node LECs predicted that LECs from the ceiling of the SCS and LECs from collecting lymphatic vessels have a shared origin, and they give rise to lymphatic valve cells (Figure 1B, step 4). LECs lining the floor of the SCS and in the MS also seem to have the same origin. No functional experiments have been performed yet to confirm these predicted developmental paths of LN LECs. Better insights in LN development are needed and will also help to further determine how tertiary lymphoid organs (TLOs) develop. TLOs, induced by LTi cells, are frequently found in chronically inflamed tissue, where they regulate inflammation and consist of many stromal cells, also found in lymph nodes [49]. Thus far, TLOs have not been investigated with scRNA-seq methodology yet.

3. The Heterogenous Family of LECs

3.1. Lymph Node LEC Subpopulations

LECs mark the outline of SCS, CS, and MS in the LN, along with other stromal cells, such as fibroblastic reticular cells and perivascular cells [50]. ScRNA-seq studies have revealed transcriptomic differences between LEC populations lining these different parts of the LN, and they have identified subpopulations that delineate the floor of SCS (floor LECs), the ceiling of the SCS (ceiling LECs), the medullary sinus (medullary LECs), and the cortical sinus (cortical LECs) [51][52][53][51,52,53]Figure 2 gives an overview of the different subpopulations, their spatial location, markers, and proposed functions. Differences between mice and human LN subpopulations may be species-specific, but they are more likely the result of technical differences related to the scRNA-seq technique or sequencing depth.
Figure 2. Overview of the different lymph node lymphatic endothelial cells and peripheral lymphatic endothelial cell subpopulations. The lymph node (LN) consists of multiple sinuses to transport lymph and immune cells from the afferent lymphatic vessels to the efferent lymphatic vessel. Immune cells arrive in the subcapular sinus (SCS) and migrate over the floor of the SCS, lined by floor LECs, to the LN cortex. The SCS is a narrow sinus, and it is close to the LN capsule. The LN capsule side of the SCS is lined by ceiling LECs. In the LN cortex, immune cells also arrive via high endothelial venules (HEVs). Egress of immune cells from the cortex is via the cortical sinuses (CS) lined by cortical LECs. These cortical sinuses drain in the medulla. Medullary sinuses (MS), lined by medullary LECs, transport lymph and immune cells to the efferent lymphatic vessel. LECs lining the different sinuses were found to be transcriptionally different and to have a range of different functions.

3.2. Peripheral LEC Subpopulations

In the periphery, lymphatic capillaries resorb tissue fluid, antigens, and growth factors and enable transport of migration immune cells by means of their oak-leaf shape and button-like junctions and a discontinuous basal membrane by which leukocytes or macromolecules can pass (Figure 2) [47][54][55][47,62,63]. Smooth muscle cells control further propulsion of lymph to the collecting lymphatic vessels. Intraluminal valves ensure unidirectional flow of lymph in the collecting vessels [56][64].
ScRNA-seq data of peripheral LECs is limited. Takeda et al. sequenced a few LN surrounding collecting lymphatic vessels and lymphatic capillaries. The LECs from the afferent collecting vessels shared the expression of CD73 and CAV1 with ceiling LECs of the SCS. PDPN, LYVE-1, and CCL21 were highly expressed on lymphatic capillaries, and CLDN11 was found as a marker for valve LECs [53]. It is unclear if LEC subpopulations within capillaries or collecting vessels exist. Furthermore, if LECs populations differ between organs is unknown. However, in a large scRNA-seq study of endothelial cells in mouse organs, LECs from different organs clustered together and did not seem transcriptionally different [57][65]. Another large mouse endothelial scRNA-seq study also observed minimal differences between LECs from different tissues, except for LECs from the gut [58][66].
Multiple scRNA-seq studies sequencing whole organs or tissues have found an LEC population identified by the expression of LYVE-1, PROX-1, PDPN, VEGFR3, and CCL21 [58][59][60][61][62][66,67,68,69,70]. In these studies, a few additional molecules were found to be expressed by LECs. Murine LECs in multiple organs express other molecules, such as Thy1, Mmrn1, and Esam, mostly involved in extra cellular matrix remodeling and cell migration [58][66]. In the cornea, LECs were found to additionally express many major histocompatibility complex (MHC) I genes and angiogenic privilege factors [59][67]. Lung LECs highly express SEMD3D, a semaphorin that is known to guide entry of DCs into lymphatics [61][63][69,71]. In the mouse liver, LECs express Mmrn1 but also Tbx1, a molecule involved in lymphangiogenesis, and Ahnak2 [37]Ahnak2 is a large protein, and it has been found in endothelial cells in blood–tissue barriers [64][72]. The LECs in the mouse liver also highly expressed Cd36, indicating their role in LDL transport [18]. LECs were also found in a scRNA-seq study of human placentas. These LECs were specifically located in the maternal portion of the placenta, the decidua basalis. Pathways for cell migration and adhesion were highly enriched in these cells, indicating a role in regulating immune cell inflex in the decidua basalis [65][73]. The LEC clusters in most studies only consist of at most 1% of the sequenced cells. More in-depth scRNA-seq studies will hopefully provide more insights into LECs from different organs. Functional studies into some of the found expressed markers are also lacking and could potentially uncover new functions of LECs.