1. Introduction
The lymphatic vascular system was already known in ancient times, although its functions were misinterpreted.
Venae albae or
Ductus lactei were already known to Hippocrates of Kos (460-370 BC), Aristoteles (384 – 322 BC), and the doctors of the Alexandrian school (approx. 300 BC – 600 AC) (cited from [
1]). And when Gaspare Aselli (re)discovered the lymphatics [
2], he knew very well what he had to look for. He realized that the
Venae albae were easy to find in the mesenteries e.g. of dogs, but he was probably the first to make a connection of these milky vessels with food intake. The anatomical representation of the lymphatic vascular system was especially advanced at Italian universities [
3], although the function of the system was still very mysterious. One of the first to recognize the importance of directed fluid transport in the lymphatics was [
4]. Contrary to the prevailing doctrine, which postulated a lymph flow into the liver, he described the connection to the central venous system. He also recognized a valve system in lymphatics and the coagulability of the lymph.
We should be very humble in evaluating ancient knowledge, because one thing is certain: we are still far from fully understanding the functions of the lymphatic vascular system, although more and more functions are being discussed [
5,
6,
7,
8]. The complexity is immediately apparent when we visualize the heterogeneity of lymphatic endothelial cells (LECs) [
9], and the long lists of molecules expressed in LECs [
10,
11], lymph collectors [
12] or lymph nodes [
13,
14]. Thereby, development and behavior of LECs is critically regulated by the transcription factor PROX1 [
15,
16,
17]; and the histological characterization of human LECs is typically performed with antibodies against PROX1 and CD31 [
18], the Vascular endothelial growth factor receptor-3 (VEGFR3, FLT4) [
19], the Lymphatic vessel endothelial hyaluronic acid receptor 1 (LYVE1) [
20], and the type-I integral membrane glycoprotein Podoplanin (PDPN) [
21].
Congenital malformations of the blood vascular system are very common, those of the lymphatic system are much rarer, mostly located in the head-neck region and can be life threatening. Most of the lymphatic malformations are caused by somatic mutations in genes involved in the VEGFR3 signaling pathway [
22,
23,
24,
25]. However, most genetic causes of lymphatic malformations and primary lymphedema have not yet been discovered. This illustrates that there are still a very large number of molecules whose significance for the development and function of the lymphatic vessels is not yet known. We have previously used RNASeq to study human foreskin-derived LECs under normoxic and hypoxic conditions, and defined 162 genes that are significantly regulated by hypoxia [
10], as well as highly expressed genes that influence the composition of the extracellular matrix and may be involved in lymphedema-induced fibrosis [
26]. Here, we went through the RNASeq list and studied the large number of highly expressed LEC genes. We compared RNA expression with protein expression by systematically studying the Human Protein Atlas [
27] (
https://www.proteinatlas.org/). In this way, we sought to define additional molecules that can be used for deeper characterization of lymphatic vessels and to define further functions of LECs.
2. Results and Discussion
We recently performed expression analyses of three well characterized human foreskin-derived LEC lines under normoxic and hypoxic conditions [
10,
26], and received a list of approx. 16,000 LEC-expressed genes. Only a very small number of these molecules have been studied in LECs so far. We have tried to systematically match RNA expression with protein expression using the Human Protein Atlas [
27]. We have concentrated on molecules whose function has not yet been investigated in LECs. However, the number of matches between RNA and protein expression was not very high, but there can be many technical reasons (fixation time of tissue, paraffine permeability of the antibodies, etc.) for this.
One of the main functions of the lymphatics is regulation of fluid homeostasis. The fluid is mainly absorbed via the flexible microvalves of the initial lymphatic vessels. However, part of the fluid uptake also occurs by means of active transcytosis, which has also been observed for the uptake of chylomicrons by lacteals [
28,
29,
30]. Very recently the uptake of myofibroblast-derived microvesicles (MVs) was studied in human dermal BECs and LECs [
31]. The authors describe that MVs cross an LEC layer but not a BEC layer
in vitro. Internalization into BECs is mainly via the caveolin-dependent pathway, while LECs mainly use the macropinocytosis pathway. High expression of caveolin 1 (CAV1) and CAV2 (
Table 1), for endo- and transcytosis, seems to be in line with the free passage of MVs through LECs; and antibodies against CAV1 and CAV2 clearly stain lymphatics, and selected blood vessels (
Figure 1 and
Figure 2).
Very active transcytosis by LECs is also reflected by the RNA expression the ferlin family members: myoferlin (
MYOF; very high), dysferlin (
DYSF; high), and otoferlin (
OTOF; moderate) (
Table 1). Ferlins regulate membrane fusion, and fusion of vesicles to cell membranes e.g. for exocytosis or membrane regeneration [
32]. Of the three ferlins, we observed clear immunostaining against DYSF (dystrophy-associated fer-1-like protein) in various organs (
Figure 3) and performed immunodouble staining with CD31 in human foreskin (
Figure 4). We observed DYSF in dermal lymphatics. Blood vessels were mostly negative, except for the subepithelial capillary plexus. We tested various OTOF antibodies, but did not receive a positive result; and, in contrast to the high RNA expression, we found no immunopositivity for MYOF in the Human Protein Atlas. However, DYSF is clearly present in LECs, which is in line with high transcytotic activity.
Both the initial lymphatic vessels and the lymph collectors exhibit strong functional fluctuations in their diameter, and they have valves, whose elasticity is of great importance for reliable valve closure. Spectrins have originally been identified as the major elastic component of erythrocytes, linking the actin cytoskeleton to the cell membrane. However, nonerythroid spectrins have also been identified [
33,
34,
35]. Spectrins form tetrameric proteins of alpha and beta subunits. SPTAN1 (spectrin, alpha, nonerythrocytic 1) and SPTBN1 (spectrin, beta, nonerythrocytic 1) are highly expressed in LECs (
Table 1), and they are well detectable at protein level in lymphatics of various organs (
Figure 5 and
Figure 6). Neuropathies related to malfunction of the two molecules have been observed, however, lymphedema has not been described yet.
Among the genes most highly expressed in LECs is Multimerin 1 (MMRN1) (
Table 1). Lymph contains fibrinogen and can coagulate [
36], and we have previously pointed out that LECs are an important source for factor VIII (F8, antihemophilic globulin A) and its carrier protein von-Willebrand factor (VWF) [
8]. MMRN1 is a specific coagulation factor V binding platelet protein with a role in hemostasis and coagulation, and accordingly it is highly expressed in megakaryocytes [
37]. Besides that, it is found with extremely high specificity in endothelial cells, however, a clear distinction between BECs and LECs has not been made [
38]. However, it appears that LECs are clearly stronger positive than BECs (
https://www.proteinatlas.org/ENSG00000138722-MMRN1/single+cell+type). We found immunopositivity for MMRN1 in lymphatics of various organs (
Figure 7). During lymphostasis, the accumulation of pro-coagulatory factors (F8, VWF, MMRN1) and the decrease in the anti-coagulatory and anti-inflammatory 5'-nucleotidase (CD73) of the LECs [
39] can be important reasons for the increased tendency to thrombosis and inflammation.
The MARCKS-like protein-1 (MARCKSL1 = MLP) has great similarity with the myristoylated, alanine-rich protein MARCKS [
40,
41], a substrate for protein kinase C. MARCKSL1 has mainly been studied in neural development, and in cancer cell migration [
42]. High MARCKSL1 expression has a strong prognostic value in lymph node-negative breast cancer patients [
43]. Upon phosphorylation, MARCKSL1 induces actin bundling and inhibits cell migration [
42]. We observed high RNA expression of MARCKSL1 in LECs (
Table 1), and immune-positive lymphatics in various organs (
Figure 8). The influence of MARCKSL1 on LEC stability and function has not been studied yet.
Another molecule associated with the actin cytoskeletal system is calponin 3 (CNN3). It consists of an acidic C terminus and a basic N terminus [
44,
45]. CNN3 regulates contractility of actomyosin-containing stress fibers of non-muscle cells [
46]. Thereby, control of stress fiber contractility by CNN3 was found to be associated with mechanosensitive Yap/Taz (Yes-associated protein/transcriptional coactivator with PDZ binding motif) transcriptional activation [
47]. We observed high RNA expression of CNN3 in LECs (
Table 1) and immune-positive lymphatics in various organs (
Figure 9), making LECs an attractive model for further studies on the function of CNN3. The importance of the Yap/Taz signaling in the Hippo pathway for lymphangiogenesis has been well recognized [
48,
49].
The cytoskeleton of LECs appears to be very stable, and accordingly, the invasiveness of LECs is usually extremely low; the only exception being the vanishing bone disease (Gorham-Stout disease, GHS) with approx. 350 cases being reported [
50,
51,
52]. The ankyrin repeat domain protein 37 (ANKRD37) is not only associated with preeclampsia during pregnancy. Its knock-down enhances trophoblast invasiveness, migration, and regulation of key invasion proteins [
53]. ANKRD37 is expressed at very low levels in LECs and significantly upregulated by hypoxia (
Table 1) [
26]. It is a hypoxia-inducible factor-1 (HIF1) target gene [
54]. Its expression in GHS has, to the best of our knowledge, not been studied yet. We observed immune-positivity for ANKRD37 in LECs (Figure 10), however, due to the large number of ANKRD family members expressed in LECs (data not shown), specificity of antibody staining might be problematic.
Another molecule involved in actin stress fiber formation, and containing an ankyrin repeat domain, is KANK3 (KN motif- and ankyrin repeat domain-containing protein 3) [
55]. It is highly expressed in LECs (
Table 1), and immune-positivity is seen in lymphatics of various organs (
Figure 11). The four members of the KANK family regulate integrin-mediated adhesion, actomyosin contractility, and link focal adhesions to the cortical microtubule stabilization complex [
56]. Specific expression of KANK3 in endothelial cells has been noted [
56], but functional studies are lacking.
Lymph contains high amounts of immunoglobulin G (IgG) heavy chain [
36]. Antigen-IgG complexes represent a starting point for the classical complement system, which can stimulate phagocytes, inflammation and the cell-killing membrane attack complex (MAC) [
57]. LECs express high amounts of CD59, aka protectin (
Table 1), a glycoprotein functioning as a membrane-bound inhibitor of MAC [
58]. Malfunctioning of CD59 causes hemolytic anemia with immune-mediated polyneuropathy [
59]. The specific function in LECs has not been investigated yet. Immune-positivity is seen in lymphatics and blood vessels (
Figure 12), suggesting an important immune-suppressive function in vessels. Strong expression of CD59 in colorectal cancer is associated with higher incidence of lymph node metastasis [
55]. In ulcerative colitis, a chronic inflammatory disease of the colon, CD59 belongs to the complement components that have shown to safeguard the intestinal barrier and reduce intestinal inflammation [
60].
It was previously shown that the non-canonical WNT (Wingless-type MMTV integration site) signaling pathway is of great importance for the elongation of lymphatics during embryonic lymphangiogenesis [
61,
62]. Nucleoredoxin (NXN) is a 48kDa protein and a redox-dependent negative regulator of the Wnt signaling pathway [
63]. NXN mRNA is highly expressed in LECs (
Table 1), and by immunohistology it can be detected in lymphatics of various organs (
Figure 13). Similar to mutations in
WNT5A [
64] and other WNT signaling members, mutations in
NXN cause the Robinow syndrome [
65], which is mainly characterized by dysmorphic facial features and short-limbed dwarfism. To our knowledge, morphology and function of the lymphatics have not yet been investigated in Robinow patients.
Figure 1.
Immunostaining of CAV1 in lymphatics (arrows) of human A) oral mucosa, Antibody CAB003791 B) rectum, Antibody CAB003791 C) Breast, Antibody CAB003791 and D) parietal layer of lymph node marginal sinus, Antibody HPA049326. From: The Human Protein Atlas.
Figure 1.
Immunostaining of CAV1 in lymphatics (arrows) of human A) oral mucosa, Antibody CAB003791 B) rectum, Antibody CAB003791 C) Breast, Antibody CAB003791 and D) parietal layer of lymph node marginal sinus, Antibody HPA049326. From: The Human Protein Atlas.
Figure 2.
Immunostaining of CAV2 in lymphatics (arrows) of human A) Esophagus, Antibody HPA044810 B) Colon, Antibody CAB013488, C) Oral mucosa, Antibody HPA044810 and D) Skin, Antibody HPA044810. From: The Human Protein Atlas.
Figure 2.
Immunostaining of CAV2 in lymphatics (arrows) of human A) Esophagus, Antibody HPA044810 B) Colon, Antibody CAB013488, C) Oral mucosa, Antibody HPA044810 and D) Skin, Antibody HPA044810. From: The Human Protein Atlas.
Figure 3.
Immunostaining of DYSF in lymphatics (arrows) of human A) Breast, Antibody CAB002510, B) Duodenum, Antibody CAB002510, C) Colon, Antibody HPA017071, D) Skin anal, Antibody CAB002510. From: The Human Protein Atlas.
Figure 3.
Immunostaining of DYSF in lymphatics (arrows) of human A) Breast, Antibody CAB002510, B) Duodenum, Antibody CAB002510, C) Colon, Antibody HPA017071, D) Skin anal, Antibody CAB002510. From: The Human Protein Atlas.
Figure 4.
Immunostaining of DYSF (red) and CD31 (green) in lymphatics of human foreskin. Blood vessels strongly express CD31; lymphatics have a punctate weak staining. Nuclei are stained blue with Dapi. In lymphatics the red DYSF staining is dominant, but subepithelial capillaries also express DYSF. A-C) 10x objective; Bar = 100 µm, D) 40x objective; Bar = 25µm.
Figure 4.
Immunostaining of DYSF (red) and CD31 (green) in lymphatics of human foreskin. Blood vessels strongly express CD31; lymphatics have a punctate weak staining. Nuclei are stained blue with Dapi. In lymphatics the red DYSF staining is dominant, but subepithelial capillaries also express DYSF. A-C) 10x objective; Bar = 100 µm, D) 40x objective; Bar = 25µm.
Figure 5.
Immunostaining of SPTAN1 in lymphatics (arrows) of human A) Oral mucosa, Antibody HPA007927, B) Small intestine. Antibody CAB004581, C) Breast Antibody HPA007927, D) Epididymis, Antibody HPA007927. From: The Human Protein Atlas.
Figure 5.
Immunostaining of SPTAN1 in lymphatics (arrows) of human A) Oral mucosa, Antibody HPA007927, B) Small intestine. Antibody CAB004581, C) Breast Antibody HPA007927, D) Epididymis, Antibody HPA007927. From: The Human Protein Atlas.
Figure 6.
Immunostaining of SPTBN1 in lymphatics (arrows) of human A) Esophagus, Antibody HPA013149, B) Colon, Antibody HPA013149, C) Breast, Antibody HPA013149. From: The Human Protein Atlas.
Figure 6.
Immunostaining of SPTBN1 in lymphatics (arrows) of human A) Esophagus, Antibody HPA013149, B) Colon, Antibody HPA013149, C) Breast, Antibody HPA013149. From: The Human Protein Atlas.
Figure 7.
Immunostaining of MMRN1 in lymphatics (arrows) of human A) Skin, Antibody HPA035769, B) Oral mucosa, Antibody HPA035769, C) Urinary bladder, Antibody HPA035769. From: The Human Protein Atlas.
Figure 7.
Immunostaining of MMRN1 in lymphatics (arrows) of human A) Skin, Antibody HPA035769, B) Oral mucosa, Antibody HPA035769, C) Urinary bladder, Antibody HPA035769. From: The Human Protein Atlas.
Figure 8.
Immunostaining of MARCKSL1 in lymphatics (arrows) of human A) Oral mucosa, Antibody HPA030528, B) Esophagus, Antibody HPA030528, C) Duodenum, Antibody HPA030528, D) Cervix uteri, Antibody HPA030528. From: The Human Protein Atlas.
Figure 8.
Immunostaining of MARCKSL1 in lymphatics (arrows) of human A) Oral mucosa, Antibody HPA030528, B) Esophagus, Antibody HPA030528, C) Duodenum, Antibody HPA030528, D) Cervix uteri, Antibody HPA030528. From: The Human Protein Atlas.
Figure 9.
Immunostaining of CNN3 in lymphatics (arrows) of human A) Oral mucosa, Antibody HPA051237, B) Esophagus, Antibody CAB009849, C) Skin, Antibody CAB009849, D) Colon, Antibody CAB009849. From: The Human Protein Atlas.
Figure 9.
Immunostaining of CNN3 in lymphatics (arrows) of human A) Oral mucosa, Antibody HPA051237, B) Esophagus, Antibody CAB009849, C) Skin, Antibody CAB009849, D) Colon, Antibody CAB009849. From: The Human Protein Atlas.
Figure 10.
Immunostaining of ANKRD37 in lymphatics (arrows) of human A) Breast, Antibody HPA036626, B) Esophagus, Antibody HPA036626, C) Colon, Antibody HPA036626, D) Lymph node, Antibody HPA036626. From: The Human Protein Atlas.
Figure 10.
Immunostaining of ANKRD37 in lymphatics (arrows) of human A) Breast, Antibody HPA036626, B) Esophagus, Antibody HPA036626, C) Colon, Antibody HPA036626, D) Lymph node, Antibody HPA036626. From: The Human Protein Atlas.
Figure 11.
Immunostaining of KANK3 in lymphatics (arrows) of human A) Skin, (obviously edematous), Antibody HPA051153, B) Colon, Antibody HPA051153, C) Rectum, Antibody HPA051153, D) Lymph node, Antibody HPA051153. From: The Human Protein Atlas.
Figure 11.
Immunostaining of KANK3 in lymphatics (arrows) of human A) Skin, (obviously edematous), Antibody HPA051153, B) Colon, Antibody HPA051153, C) Rectum, Antibody HPA051153, D) Lymph node, Antibody HPA051153. From: The Human Protein Atlas.
Figure 12.
Immunostaining of CD59 in lymphatics (arrows) of human A) Skin, Antibody HPA026494, B) Breast, Antibody HPA0264949, C) Oral mucosa, Antibody HPA026494, D) Colon, Antibody HPA026494. From: The Human Protein Atlas.
Figure 12.
Immunostaining of CD59 in lymphatics (arrows) of human A) Skin, Antibody HPA026494, B) Breast, Antibody HPA0264949, C) Oral mucosa, Antibody HPA026494, D) Colon, Antibody HPA026494. From: The Human Protein Atlas.
Figure 13.
Immunostaining of NXN in lymphatics (arrows) of human A) Skin, Antibody HPA023566, B) Breast, Antibody HPA023566, C) Duodenum, Antibody HPA023566, D) Colon, Antibody HPA023566. From: The Human Protein Atlas.
Figure 13.
Immunostaining of NXN in lymphatics (arrows) of human A) Skin, Antibody HPA023566, B) Breast, Antibody HPA023566, C) Duodenum, Antibody HPA023566, D) Colon, Antibody HPA023566. From: The Human Protein Atlas.
Table 1.
RNASeq analysis of 3 defined human dermal LEC lines under normoxia and hypoxia.
Table 1.
RNASeq analysis of 3 defined human dermal LEC lines under normoxia and hypoxia.
gene_id |
gene_name |
chromos. |
start |
end |
width |
strand |
gene_biotype |
HDLEC-5 |
HDLEC-6 |
HDLEC-7 |
Hypox-5 |
Hypox-6 |
Hypox-7 |
ENSG00000186352 |
ANKRD37 |
4 |
1.85E+08 |
1.85E+08 |
4608 |
+ |
protein_cod. |
92.68 |
117.12 |
124.35 |
312.56 |
549.57 |
650.6 |
ENSG00000105974 |
CAV1 |
7 |
1.17E+08 |
1.17E+08 |
36184 |
+ |
protein_cod. |
28254.02 |
22816.55 |
40041.39 |
41886.36 |
29070.25 |
44229.96 |
ENSG00000105971 |
CAV2 |
7 |
1.16E+08 |
1.17E+08 |
221162 |
+ |
protein_cod. |
5923.02 |
6362.42 |
8338.77 |
8258.02 |
7425.45 |
10627.94 |
ENSG00000085063 |
CD59 |
11 |
33703010 |
33736491 |
33482 |
- |
protein_cod. |
25996.8 |
30542.27 |
33419.64 |
25877.99 |
36425.42 |
34222.24 |
ENSG00000117519 |
CNN3 |
4 |
94896949 |
94927223 |
30275 |
- |
protein_cod. |
16690.62 |
24868.18 |
11453.55 |
14356.97 |
20710.53 |
8937.79 |
ENSG00000135636 |
DYSF |
2 |
71453722 |
71686768 |
233047 |
+ |
protein_cod. |
2655.52 |
1698.22 |
4384.61 |
3353.58 |
2204.57 |
5193.08 |
ENSG00000186994 |
KANK3 |
19 |
8322584 |
8343262 |
20679 |
- |
protein_cod. |
2257.23 |
4007.72 |
2948.58 |
3945.85 |
4201.93 |
5177.89 |
ENSG00000175130 |
MARCKSL1 |
1 |
32333839 |
32336233 |
2395 |
- |
protein_cod. |
8786.11 |
9698.24 |
8003.97 |
8408.33 |
7819.15 |
4811.13 |
ENSG00000138722 |
MMRN1 |
4 |
89879532 |
89954629 |
75098 |
+ |
protein_cod. |
19949464 |
204322.7 |
16868.85 |
263839.6 |
303443.4 |
38908.39 |
ENSG00000138119 |
MYOF |
10 |
93306429 |
93482334 |
175906 |
- |
protein_cod. |
12506.28 |
4986.79 |
5959.34 |
11865.43 |
2877.56 |
11681.51 |
ENSG00000167693 |
NXN |
17 |
799310 |
979776 |
180467 |
- |
protein_cod. |
5646.52 |
7277.79 |
2263.45 |
4812.87 |
5743.42 |
2401.48 |
ENSG00000115155 |
OTOF |
2 |
26457203 |
26558698 |
101496 |
- |
protein_cod. |
1187.21 |
1087.97 |
463.93 |
498.71 |
109.01 |
146 |
ENSG00000197694 |
SPTAN1 |
9 |
1.29E+08 |
1.29E+08 |
81105 |
+ |
protein_cod. |
10798.23 |
8180.84 |
6609.8 |
8411.32 |
9523.71 |
8408.68 |
ENSG00000115306 |
SPTBN1 |
2 |
54456317 |
54671446 |
215130 |
+ |
protein_cod. |
23119.16 |
28282.09 |
11418.88 |
25384.26 |
39159.74 |
14476.62 |