# Three-dimensional imaging of vascular development in the mouse epididymis

Christelle Damon-Soubeyrand, Antonino Bongiovanni, Areski Chorfa, Chantal Goubely, Nelly Pirot, Luc Pardanaud, Laurence Piboin-Fragner, Caroline Vachias, Stephanie Bravard, Rachel Guiton, et al.

► To cite this version:

Christelle Damon-Soubeyrand, Antonino Bongiovanni, Areski Chorfa, Chantal Goubely, Nelly Pirot, et al.. Three-dimensional imaging of vascular development in the mouse epididymis. eLife, 2023, 12, pp.e82748. 10.7554/eLife.82748 . inserm-04257565

HAL Id: inserm-04257565

https://inserm.hal.science/inserm-04257565

Submitted on 25 Oct 2023

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_1.png" /></figure>
<figure><img src="image_2.png" /></figure>

# Three-dimensional imaging of vascular development in the mouse epididymis

Christelle Damon-Soubeyrand¹†, Antonino Bongiovanni²†, Areski Chorfa¹†, Chantal Goubely¹, Nelly Pirot³, Luc Pardanaud⁴, Laurence Piboin-Fragner⁴, Caroline Vachias¹, Stephanie Bravard¹, Rachel Guiton¹, Jean-Leon Thomas⁵,⁶, Fabrice Saez¹, Ayhan Kocer¹, Meryem Tardivel²*, Joël R Drevet¹*, Joelle Henry-Berger¹

¹Institut « Génétique Reproduction & Développement », UMR CNRS 6293 - Inserm U1103 - Université Clermont Auvergne, Faculté de Médecine, Clermont-Ferrand, France; ²Université Lille, CNRS, Inserm, CHU Lille, US 41 - UAR 2014-PLBS, Lille, France; ³Réseau d'Histologie Expérimentale de Montpellier, Institut de Recherche en Cancérologie de Montpellier, Inserm U1194 - Université Montpellier - ICM, Campus Val d'Aurelle, Montpellier, France; ⁴Université de Paris, PARCC, INSERM, Paris, France; ⁵Department of Neurology, Yale University School of Medicine, New Haven, United States; ⁶Institut du Cerveau et de la Moelle, Inserm, Université Pierre et Marie Curie, Paris, France

**Abstract** Long considered an accessory tubule of the male reproductive system, the epididymis is proving to be a key determinant of male fertility. In addition to its secretory role in ensuring functional maturation and survival of spermatozoa, the epididymis has a complex immune function. Indeed, it must manage both peripheral tolerance to sperm antigens foreign to the immune system and the protection of spermatozoa as well as the organ itself against pathogens ascending the epididymal tubule. Although our knowledge of the immunobiology of this organ is beginning to accumulate at the molecular and cellular levels, the organization of blood and lymphatic networks of this tissue, important players in the immune response, remains largely unknown. In the present report, we have taken advantage of a VEGFR3:YFP transgenic mouse model. Using high-resolution three-dimensional (3D) imaging and organ clearing coupled with multiplex immunodetections of lymphatic (LYVE1, PDPN, PROX1) and/or blood (PLVAP/Meca32) markers, we provide a simultaneous deep 3D view of the lymphatic and blood epididymal vasculature in the mature adult mouse as well as during postnatal development.

## Editor's evaluation

This fundamental work substantially advances our understanding of vessel development in mouse epididymis. The evidence supporting the conclusions is compelling, with rigorous state-of-the-art microscopy involving high-resolution Three-dimensional imaging. The work will be of broad interest to cell biologists and male reproductive biologists.

## Introduction

Studies regarding the etiology of male infertility have revealed that at least 15% of cases have a clear immunological origin, a share that is most likely underestimated if idiopathic infertilities were to be included (Dohle et al., 2005; Jungwirth et al., 2012). In the male, they may be the result of chronic inflammatory situations, acute bacterial and/or viral infections, or autoimmune processes along the

*For correspondence: meryem.tardivel@univ-lille.fr (MT); joel.drevet@uca.fr (JRD)

†These authors contributed equally to this work

**Competing interest:** The authors declare that no competing interests exist.

**Funding:** See page 20

**Received:** 16 August 2022

**Preprinted:** 28 August 2022

**Accepted:** 18 May 2023

**Published:** 13 June 2023

**Reviewing Editor:** T Rajendra Kumar, University of Colorado, United States

© Copyright Damon-Soubeyrand, Bongiovanni, Chorfa et al. This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
genital tract leading to the production of anti-sperm antibodies (for a recent review, see *Shibahara et al., 2021*). These situations are difficult to diagnose and treat, in part because our knowledge of immune/inflammatory responses in male accessory organs is limited. Among the accessory organs, the epididymis plays an essential role in the maturation, storage, and protection of spermatozoa. A particularity of the mammalian epididymis is its highly segmented anatomical organization with three macroscopic territories (caput, corpus, and cauda as depicted in Figure 2), themselves sectorized into segments (10 for the mouse model, see Figure 2), some of which have distinct functional roles (*Turner et al., 2003; Johnston et al., 2005*). Unlike the testis, where mammalian evolution has chosen to tightly seal the germline from the immune system—conferring immune privilege to this tissue (*Meinhardt and Hedger, 2011*)—the epididymis tubule faces multiple immune challenges (*Guiton et al., 2013; Hedger, 2011*). On the one hand, the epididymis tubule is the testis gatekeeper, preventing ascending pathogens from reaching the immune-privileged seminiferous tubules. To do so, it should be endowed with all the power of a full and efficient immune response toward exogeneous antigens. On the other hand, the epididymis must maintain self-tolerance toward spermatic antigens that are unique to this cell lineage and appear at puberty long after the establishment of the self-immune repertoire (*Goodnow, 1996*). Understanding how this dual action of peripheral tolerance toward sperm antigens (*Mueller, 2010*) versus efficient immune survey toward non-self-antigens is performed is of particular interest for infertility issues, but also more largely in other settings where this situation occurs. This has prompted recent research aiming at identifying immune cells and molecules that could participate in the finely orchestrated epididymal immune physiology. This research has led to the characterization of a dense network of peritubular antigen-presenting cells (*Da Silva et al., 2011*), interstitial and intra-epithelial lymphocytes of distinct sub-families (*Voisin et al., 2018*), as well as of immunosuppressive players such as transforming growth factor beta (*Pierucci-Alves et al., 2018; Voisin et al., 2020*) and indoleamine/tryptophane dioxygenase activity (*Britan et al., 2006; Jrad-Lamine et al., 2011; Jrad-Lamine et al., 2013*).

Because the immune picture of the mammalian epididymis would not be complete without a clear view of the blood and lymphatic circulating networks—which are crucial in managing immune responses—we choose to explore these in the mouse. Reports describing mammalian epididymal blood and lymphatic vessels are rare and rather old (*Pérez-Clavier et al., 1982; Abe et al., 1984; McDonald and Scothorne, 1988; Itoh et al., 1998*). The technology used was not very efficient and consisted mainly of ink or contrast agent injection. It concerned only the peripheral epididymal lymphatic network with the organization of the collectors connecting the PP. More recently, *Hirai et al., 2010*, performed the first histological study of mouse epididymal blood and lymphatic vasculature using the markers CD31/PECAM1 and LYVE1, respectively. Although PECAM1 is used to visualize blood vessels, it is not an exclusive blood marker because it is also found expressed at button-like junctions of endothelial cells in the initial lymphatics (*Privratsky and Newman, 2014*). Furthermore, to characterize lymphatics, the use of LYVE1 alone limits the power of the analysis because this marker recognizes only a subpopulation of all lymphatics (*Kato et al., 2006*). In addition, conventional two-dimensional (2D) light microscopy combined with a horseradish peroxidase (HRP) reporter on paraffin-embedded tissue sections did not allow for a high level of accuracy in providing a holistic view of a structured network.

Characterization of the lymphatic vasculature in mammals is tricky because on the one hand, it shares common embryonic origins with the blood vessels (*Baluk and McDonald, 2008; Srinivasan et al., 2007*) and on the second hand, it concerns various different structures. There are blunted-end highly permeable initial lymphatics connected to pre-collector lymphatics and then to collector lymphatics linking to a lymph node (LN), each with specific characteristics (for a very recent review, see: *Lampejo et al., 2023*). A step forward was made with the identification of lymphangiogenic factors (vascular endothelial growth factors VEGF-C and VEGF-D) and their receptor (VEGFR3), which are now considered the major markers of lymphatics (*Kaipainen et al., 1995; Mäkinen et al., 2001*). To ensure confidence in the identification of lymphatic structures, one should use a combination of markers, including VEGFR3, the lymphatic vessel endothelial hyaluronan receptor 1 (LYVE1), the prospero-related homeobox 1 (PROX1) transcription factor, and podoplanin (PDPN) (*Banerji et al., 1999; Breiteneder-Geleff et al., 1999; Hamrah et al., 2003; Kivelä et al., 2016; Mouta Carreira et al., 2001; Partanen et al., 2000; Wigle and Oliver, 1999*). We have followed this rationale in the present study, using three lymphatic markers together with plasmalemma vesicle-associated protein
(PLVAP/MECA32) to specifically identify fenestrated blood vessels (Stan et al., 2004). In addition, we have taken advantage of a transgenic mouse model that expresses the fluorescent reporter YFP under the control of the VEGFR3 promoter (Calvo et al., 2011). Using the power of three-dimensional (3D) reconstruction of high-resolution imaging with light-sheet microscopy after epididymis 3DISCO clearing (Belle et al., 2014; Ertürk et al., 2012), we present here an in-depth analysis of both the blood and lymphatic networks of the epididymis in the adult mouse as well as during postnatal onto-genesis. We hope the data we bring will provide a better view of the systemic compartment that irrigates and drains this segmented organ, a prerequisite for a more precise understanding of its pathophysiology.

## Results

### An abundant lymphatic network drains the epididymis

To gain insight into the lymphatic network of the mouse epididymis, we used the VEGFR3:YFP mouse model. **Figure 1** shows representative in toto views of the lymphatic vascularization of the testis and epididymis of an adult (**Figure 1A**) and 10 days postnatal (DPN) mouse (**Figure 1B, C**). There is intense vascularization at the epididymis–testis interface regardless of the level of the epididymis (**caput**, **corpus**, or **cauda**). **Figure 1B** shows that the lymphatic drainage joins the main testicular lymphatic trunk at the **pampiniform plexus** (PP) and follows it to a highly fluorescent structure that corresponds to the most proximal LN. Inset 1C and **Figure 1D** show the connections of the collectors emanating from the **caput** (upper epididymal lymphatic collector [u-ELC] in **Figure 1C, D**) and the **corpus** (median epididymal lymphatic collector [m-ELC] in **Figure 1C**). A higher magnification micrograph (**Figure 1D** ①) shows the particularly prominent lymphatic vasculature of the efferent duct and its interconnection with that of the epididymal **caput**. It also points out the departure of the u-ELC to the PP surrounding the upper epididymal artery. It can be noted that the lymphatic collectors coincide well with the segmentation of the organ. A **corpus** lymphatic collector (Co-LC) follows the body of the epididymis on the side adjacent to the testis connecting the **caput** with the **cauda** (**Figure 1E**). The inset within **Figure 1E** presents a higher magnification of the Co-LC connections. This collector also connects to the main testicular lymphatic collector at the **corpus** level (**Figure 1E** ②) as well as at the **cauda** level (**Figure 1F** ③). Also noticeable in the enlargements provided in **Figure 1E** and in **Figure 1F** is the presence of fluorescence under the capsule, faint in the **caput** (**Figure 1E**) but clearly more visible in the most distal segment (S10) of the **cauda** epididymis (**Figure 1F and G**). Finally, **Figure 1G** presents a global complex organization of the **cauda** lymphatic collectors. Overall, our observations revealed a very dense lymphatic network in the mouse epididymis.

### Complexity of the blood and lymphatic vascularization of the epididymis

Because fenestrated blood vessels also express VEGFR3, a marker specific for these vessels (PLVAP/ Meca32) was used along with three additional markers for lymphatics, specifically LYVE1, PROX1, and PDPN. As the 3DISCO organ clearing procedures resulted in a loss of the endogenous transgene (VEGFR3:YFP) fluorescence, anti-GFP was used to strengthen the signal. The concordance of the labeling obtained with the anti-GFP and an anti-VEGFR3 antibodies was controlled (see **Figure 2—figure supplement 1A**). Similarly, **Figure 2—figure supplement 1B and C** show positive controls of the multiplex labeling as well as controls with the secondary antibody, respectively. **Figure 2A** and **Figure 2—videos 1–3** show the immunoreactive vasculature after 3DISCO clearing of an adult mouse epididymis obtained after light-sheet imaging and IMARIS 3D reconstruction. The **caput** and **cauda** territories show intense PLVAP reactivity revealing a dense blood network. The **corpus** of the epididymis is also well vascularized, as can be seen in **Figure 2—figure supplement 2A and Figure 2—video 2**. The video clearly shows that the lymphatic network is asymmetrical as it is very dense on the epididymal side next to the testis. In agreement with previous reports (Abe et al., 1984; Suzuki, 1982), the initial segment (IS) of the **caput** shows the highest reactivity (see **Figures 2B and 3A** as well as **Figure 2—figure supplement 2A and B**), with each IS tubule strongly reactive at their periphery (**Figure 2—figure supplement 3A and B**). Outside the IS, PLVAP reactivity is mainly found in the inter-tubular space as well as along the septum separating each epididymal segment. This is particularly true for the IS/S2–S3 and the S9/S10 septum junctions.
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Macroscopic view of the lymphatic vasculature of the epididymis and testis in the VEGFR3:YFP model. Representative image of the lymphatic vasculature of the adult epididymis and testis observed with a Leica binocular loupe (A). VEGFR3:YFP<sup>pos</sup> lymphatics are visible in the caput, corpus, and cauda regions. Lymphatics responsible for epididymal and testicular drainage follow the pampiniform plexus (PP) before reaching a lymph node (B). The superior and lateral epididymal lymphatic collectors (upperELC and medianELC, respectively) drain lymph from the caput and corpus, respectively,</figcaption></figure>

*Figure 1 continued on next page*
*Figure 1 continued*

before joining the main testicular lymphatic collector at the PP. The latter is connected to a main collector that surrounds the testis (arrowhead) and branches into a rich network (A–C). There are numerous lymphatic connections between the lymphatics of the epididymis and the main lymphatic collector of the testis, notably through a lymphatic network at the level of the efferent ducts (φ in D), corpus (κ in E), and cauda (λ in F). Through the capsule covering the epididymal duct, there is also fluorescence that outlines the tubules, weakly in the caput (D) and more intense in the cauda (F and G). The scale bar corresponds to 1 mm. AT = adipose tissue; Ca = caput; Cd = cauda; Co = corpus; DD = deferent duct; DDA = deferent duct artery; DLC = deferent lymphatic collector; DPN = days postnatal; ED = efferent duct; LELC = lower epididymal lymphatic collector; LN = lymph node; P=pampiniform plexus; SEA = superior epididymal artery; Cd-BV = caudal blood vessel; SELC = superior epididymal lymphatic collector; T = testis; TA = testis artery; TLC = testicular lymphatic collector. ①=caput–testis lymphatic connection; ②=corpus–testis lymphatic connection; ③=cauda–testis lymphatic connections between testis and epididymis.

With the four lymphatic markers (VEGFR3:YFP, LYVE1, PROX1, and PDPN), there is an intense and complex network (Figure 2B). This appears very similar in its organization to that of the blood vessels, especially in the IS (Figure 2B). However, outside the IS, lymphatic reactivity appears more extensive than that of the blood vessel marker (see Figure 2—figure supplement 3B). As was the case for the blood vessel marker and as expected because fenestrated vessels express VEGFR3 (Partanen et al., 2000), lymphatic marker reactivity is stronger in the IS (Figure 2B). As seen in Figure 1, large external lymphatic collectors follow the septa, particularly at the IS/S2-S3 and S9/S10 boundaries. These collectors can be seen to dip inside the organ and irrigate the intertubular compartments, especially in the cauda region (see Figure 2—video 3). The LYVE1, PROX1, and VEGFR3:YFP profiles are very similar for the caput and cauda (Figure 2B). However, when viewed at higher magnification (Figure 2—figure supplement 3A), they do not merge totally in either the caput or the cauda region, especially for the small lymphatic vessels separating the epididymal segments. The PDPN profile appears different from those of LYVE1, PROX1, and VEGFR3:YFP (Figure 2B). The PDPN reactivity is more homogeneous and is concentrated at the peritubular level (Figure 2B and Figure 2—figure supplement 3A). Figure 2B and Figure 2—figure supplement 3A also show that the lymphatic marker (PROX1) gives very similar results to LYVE1 in the caput and cauda regions. However, the pattern of PROX1 is slightly distinct from that of VEGFR3:YFP and very different from that of PDPN.

We performed a surface assessment to better evaluate the density of the blood and lymphatic networks in the most proximal and distal epididymal segments (IS and S10) compared with that observed in the rest of the *caput* or *cauda*. To measure only the specific density of the lymphatic network, voxels with signal in the 555 nm channel (PLVAP/MECA32) were set to 0 in all other channels. The densities presented in *Figure 3A and B* are the results of the vessel volume normalized by the total volume of the region of interest (IS, *caput*, S10, or *cauda*). As expected, PLVAP/MECA32<sup>pos</sup> blood vessel density in the IS was approximately 3.5 times greater than in the rest of the *caput* (p<0.001). Similarly, the VEGFR3<sup>pos</sup> lymphatic density was twice that of the rest of the organ (*Figure 3B*). A result of the same magnitude was obtained for LYVE1 and PROX1 (*Figure 3B*). Because of the expression profile of PDPN, surface rendering could not be used to calculate its density. In the *cauda*, there was no significant difference in the density of blood and lymphatic vessels. Comparison of vascular density between the *caput* (IS excluded) and the *cauda* showed a slight decrease in the *cauda* for both blood and lymphatic vessels (*Figure 2—figure supplement 2B*).

Dynamic development of epididymal vasculature related to organ maturation

To better understand the postnatal development of both networks, we followed their emergence during postnatal epididymal development from 10 to 40 DPN, when the organ is functionally mature (n=5). We present in **Figure 4A** (*caput*) and **Figure 4—figure supplement 1A** (*cauda*) representative high-resolution images of the blood and lymphatic systems using three of the five markers presented above (PLVAP/MECA32, VEGFR3:YFP/GFP, and LYVE1) at 10, 20, 30, and 40 DPN. We used IMARIS software to estimate the IS volume relative to the *caput* volume (**Figure 4B**) and the S10 volume relative to the *cauda* volume (**Figure 4—figure supplement 1B**) at each time point.

IS progresses four times faster than the *caput* between 10 DPN and adulthood (A), representing less than 10% of the *caput* at 10 DPN whereas it represents 40% of the adult *caput* (**Figure 4A and B**). Following the stronger postnatal development of the IS, we observed that the *caput* vascularization is rather homogeneous at 10 DPN, whereas from 20 DPN onward, the blood vessels seem to
<figure><img src="image_4.png" /><figcaption>**Figure 2.** High-resolution three-dimensional (3D) imaging of the blood and lymphatic vasculature of the mouse epididymis after organ clearing. (A) shows a schematic representation of the caput (red) and cauda (green) regionalization and segmentation commonly used to describe the mouse epididymis. (B) Representative multiplex immunostaining of the caput (left) and cauda (right) of the epididymis using five markers recognizing blood and lymphatic vessels is shown. An anti-MECA32 antibody (red) reveals PLVAP<sup>pos</sup> blood vessels, particularly fenestrated vessels. An anti-GFP</figcaption></figure>

*Figure 2 continued on next page*
Figure 2 continued

antibody revealing the VEGFR3:YFP transgene (green) and an anti-LYVE1 antibody (cyan). Anti-LYVE1 (cyan), anti-PROX1 (magenta), and anti-PDPN (white) antibodies were used as additional lymphatic markers. After organ clearance, high-resolution 3D imaging was performed with a light-sheet ultramicroscope (LaVision BioTec). The two extreme segments, IS (initial segment or S1) and S10, are delineated by a large dashed line, whereas the caput (left) and cauda (right) regions are delineated by a small dashed line. The scale bar is 500 µm. The PLVAP image corresponds to the contro-lateral epididymis used for the four lymphatic markers since it is not possible to use all five markers on the same organ.

The online version of this article includes the following video and figure supplement(s) for figure 2:

**Figure supplement 1.** Positive and negative immunohistochemical controls.

**Figure supplement 2.** Three-dimensional (3D) view of blood and lymphatic vascularization of whole epididymis after 3DISCO clarification.

**Figure supplement 3.** Cross-sectional views of the clarified caput and cauda epididymides.

**Figure 2—video 1.** Three-dimensional (3D) imaging of multiplex labeling of blood and lymphatic networks in the adult mouse *caput* epididymis.

https://elifesciences.org/articles/82748/figures#fig2video1

**Figure 2—video 2.** Three-dimensional (3D) imaging of multiplex labeling of blood and lymphatic networks in the adult mouse *corpus* epididymis.

https://elifesciences.org/articles/82748/figures#fig2video2

**Figure 2—video 3.** Three-dimensional (3D) imaging of multiplex labeling of blood and lymphatic networks in the adult mouse *cauda* epididymis.

https://elifesciences.org/articles/82748/figures#fig2video3

concentrate in the IS segment. This is also well illustrated in the densitometric quantification shown in **Figure 4C**. This figure also shows that the density of PLVAP<sup>pos</sup> blood vessels in the IS segment is rather stable during postnatal development despite the increase in the segment volume. In contrast, there is a strong decrease between 10 and 20 DPN in the rest of the *caput*. In the *cauda* region (**Figure 4—figure supplement 1A–C**), PLVAP<sup>pos</sup> blood vessel expansion follows organ growth without a specific segmental pattern. However, as was the case in the *caput*, between 10 and 20 DPN, there is a decrease in PLVAP<sup>pos</sup> blood vessel density in the caudal region.

Looking at the two lymphatic markers (VEGFR3:YFP and LYVE1), we observed very similar kinetics of postnatal lymphatic development in the *caput* (**Figure 4A**) and *cauda* (**Figure 4—figure supplement 1A**). There is a progressive increase in VEGFR3:YFP/LYVE1 lymphatic vessels that follows postnatal organ growth in a fairly linear fashion at least between 10 and 40 DPN (**Figure 4D** and **Figure 4—figure supplement 1D**). Beyond 40 DPN, there appears to be a greater expansion of lymphatics in the *caput* and *cauda*, in the IS and S10, respectively, compared with the rest of the organ (**Figure 4D** and **Figure 4—figure supplement 1D**). This is particularly true in the IS because the size of this territory is increased approximately four times between the 10 and 40 DPN stages.

Because blood and lymphatic vascular density data suggest neovascularization during postnatal epididymal development, we performed localization and quantification of angiogenic (VEGF-A) and lymphangiogenic (VEGF-C and VEGF-D) factors, whose receptors are VEGFR2 (also known as KDR and FLK1) and VEGFR3, respectively (**Figures 5 and 6**). Quantification of the VEGFR ligands (VEGF-A, VEGF-C, and VEGF-D) was performed by western blotting as shown in **Figure 5A**. For the angiogenic ligand (VEGF-A), there was a decrease between 10 and 20 DPN followed by a plateau, and then an increase at 40 DPN; by adulthood, VEGF-A had returned to the initial level (**Figure 5B**). For the lymphangiogenic ligands (VEGF-C and VEGF-D), we observed a linear increase from 30 DPN for both ligands, consistent with the accumulation of their corresponding receptor VEGFR3 (**Figure 5B**). An analysis of covariance of the different ligand isoforms of VEGF (see **Figure 5—figure supplement 1**) with respect to the VEGFR3 receptor (see **Figure 5C**) shows strong statistical significance for the active ligand isoforms (i.e., the 21 kDa form for VEGF-C and the 24 kDa form for VEGF-D). In addition, the tissue localization of VEGFR2 is very similar to that of VEGFR3 (**Figure 6A**), being present mainly in a peritubular position at the level of the IS whereas in the more distal part of the organ it seems to be restricted to some interstitial vessels (used with permission from Dr. A. Medvinsky; VEGFR2:GFP mouse strain; Xu *et al.*, 2010). In accordance with previous work (**Korpelainen *et al.*, 1998**), the tissue localization of its angiogenic ligand (VEGF-A) seems to correspond correctly (**Figure 6A**), supporting the intense blood vascularization observed in this very proximal segment of the epididymis. Both of the lymphangiogenic ligands—VEGF-C and VEGF-D—show a similar tissue localization, being weakly present in the IS and increasingly present when moving toward the tail of the epididymis, mainly in an interstitial position (**Figure 6B**). A co-localization (Pearson coefficient = 0.684 ± 0.038) and
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Densitometry of the blood and lymphatic vasculature of the mouse epididymis. Panel A shows the surface rendering performed with IMARIS software of blood vessels (red, left) and lymphatic vessels (green, right) in the caput. The surface rendering of vessels in S1 is in white in both cases. Vessel densities shown in the graphs below correspond to the ratio of the volume occupied by vessels in the S1 or caput (minus S1) normalized to the Figure 3 continued on next page</figcaption></figure>
Figure 3 continued

total volume of the S1 or caput (minus S1), respectively. Panel B shows the surface rendering of blood (red, left) and lymphatic (green, right) vessels in the cauda. Blood or lymphatic vessel densities are measured as described in A for the caput. The Mann-Whitney test was used to determine statistical significance (***p<0.0001, **p<0.001, *p<0.01, NS = not significant).

co-occurrence analysis (Mander's coefficient VEGF-C=0.934 ± 0.028 and VEGF-D=0.569 ± 0.031) supports the hypothesis of heterodimerization of the two ligands.

Lymphatics expressing all four markers (VEGFR3:YFP, LYVE1, PROX1, and PDPN) are particularly visible at the septa (Figure 7A). We show in Figure 7B (upper panels) a dynamic 3D view of the cauda epididymis lymphangiogenesis as well as our interpretation of what is happening, presented in the lower panels of Figure 7B. LYVE1<sup>pos</sup> cubblestone cells are observed at 10 DPN while very few VEGFR3-YFP<sup>pos</sup> cells are visible. Later, cubblestone cells appear to aggregate to form tubular structures. Concomitantly, we record an increase in VEGFR3-YFP<sup>pos</sup> structures associated with a decrease in PLVAP<sup>pos</sup> vessels. This kinetics of events is consistent with active lympho-vasculogenesis during postnatal epididymal development. The schematic shown in Figure 8 presents the detailed organization of the epididymis lymphatics with respect to the tubules, septa, and orientation of the organ, with the anterior side of the epididymis tubule in close contact with the testis.

Discussion

Using whole-mount multiplex immunolabeling, organ clearing, ultramicroscopy imaging, and 3D reconstruction with a highly adapted lymphatic-YFP transgenic reporter mouse model, the sensitivity and power of our study has allowed us to present for the first time an in-depth analysis of blood and lymphatic networks both in the adult and during postnatal epididymal ontogeny. With the combination of lymphatic markers used in this work, we have shown that a dense network of conventional lymphatics can be observed in the epididymis, with initial lymphatics present in the interstitial compartment and collecting lymphatics present at the septa. We have shown that the lymphatic network is denser on the anterior face of the epididymis (face adjacent to the testis) and progresses in the organ toward the posterior face. We found that the collecting lymphatics are mainly observed on the anterior surface and radiate inward toward the posterior surface following the septa. The initial lymphatics drain the interstitial compartment and join the collectors at the septa, resulting in an asymmetric distribution (anterior/posterior) of lymphatic vascularization. This centrifuge development of the lymphatic network is easily visible (Figure 4A and Figure 2—video 2). This tree structure is logical if we consider that lymphatic vascularization progresses during mouse development following the PP before reaching the rete testis and sprouting laterally, as had been suggested earlier (Svingen et al., 2012). Each epididymal septum has lymphatic drainage, but three septa appear to be more involved, including the S1/S2–S3, S6/S7, and S8–9/S10 junctions, with the latter having the largest lymphatic collector. We expected this finding because it corresponds to the position of the arteries and veins that irrigate the epididymis (Suzuki, 1982), and it is known that the large lymphatic collectors usually follow the path of the arteries (Gruffaz, 1984). Previous researchers have hypothesized that this strict epididymal segmentation is involved in the control of ascending pathogens, although they have not provided clues as to how this might be accomplished (Stammler et al., 2015; Turner et al., 2003). Supporting this hypothesis is the observation that large lymphatic collectors that are in direct and rapid connection with LNs and the blood supply that allows for the arrival of leukocytes are present at these sites. These three septa could then be considered important checkpoints (Figure 8). This view corroborates recent data showing that immune responses in autoimmune epididymitis are concentrated in the epididymal cauda and corpus, where we do see the major lymphatic collectors (Wijayarathna et al., 2020). Although we have arbitrarily chosen to focus our discussion on the specific involvement of the epididymal vascular compartment in inflammation and immune responses, it should be mentioned that it could also serve other processes such as tissue metabolism and the stem/progenitor cell niche microenvironment associated with different regions of this organ.

Looking at the postnatal development of the epididymal lymphatic network, our data show a progressive increase that follows the postnatal growth of the organ between 10 and 40 DPN. Beyond this time point, the lymphatic network progresses more in the S1/IS and S10 compartments. This
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Evolution of the blood and lymphatic vasculature during postnatal epididymal development. Panel A shows representative three-dimensional (3D) images of blood and lymphatic networks at different postnatal stages during caput ontogenesis. Immunostaining was done with the blood vessel marker MECA32/PLVAP antibody (red) and the lymphatic marker LYVE1 antibody (cyan). The VEGFR3:YFP transgene is revealed by an anti-GFP antibody.</figcaption></figure>

Figure 4 continued on next page
Figure 4 continued

Only LYVE1 and VEGFR3 are presented here for better resolution. Immunostaining was also performed with PDPN (not shown in the figure) but PDPN data can be seen in the data sources (figure 4 data sources 1,2 and **Figure 4—figure supplement 1** data source 1&2 at https://www.ebi.ac.uk/biostudies/bioimages/studies/S-BIAD618). The dotted line indicates the contours of the initial segment (IS = S1). DPN: days postnatal. The scale bar is 200 µm for 10 DPN, 300 µm for 20 and 30 DPN, and 400 µm for 40 DPN. Panel B shows the evolution of the volume (in log10) of the S1 segment (light blue curve) and the caput region (dark blue) during postnatal development. The superimposed histogram gives the proportion of volume occupied by the S1 relative to the caput region at different stages of postnatal development. Surface rendering of blood vessels and lymphatics was performed using IMARIS software. Blood and lymphatic vessel densities were calculated as described previously (Figures 2 and 3). Quantification was performed on five replicates for the different postnatal developmental stages. Panels C and D are graphs presenting densities (mean and standard error of the mean) of blood and lymphatic vasculature, respectively. The Kruskal–Wallis test with Dunn's posttest correction was used to determine statistical significance (***p<0.0001, **p<0.001, *p<0.01).

The online version of this article includes the following figure supplement(s) for figure 4:

**Figure supplement 1.** Evolution of the blood and lymphatic vascularization of the epididymal cauda during postnatal development.

argues for active lymphangiogenesis or/and lymphangio-vasculogenesis, the two modalities by which lymphatics develop, in the caput and cauda territories. Lymphangiogenesis corresponds to the extension of existing vessels, also known as 'sprouting,' while lymphangio-vasculogenesis involves the participation of precursor cells that create clusters later integrating already formed lymphatic vessels (Gutierrez-Miranda and Yaniv, 2020; Jafree et al., 2021). Lymphangio-vasculogenesis has been reported in the development of the mesenteric lymphatic vessels (Stanczuk et al., 2015) and for the heart and kidney lymphatic networks (Jafree et al., 2019; Lioux et al., 2020), two organs with a septal organization. The lymphatic neovascularization we observed during postnatal epididymal development is supported by the observation that lymphangiogenic ligands (VEGF-C and VEGF-D) and their corresponding receptor (VEGFR3) are concordantly present in the developing tissue. Active lymphangio-vasculogenesis is also supported by our observation that foci expressing LYVE1 can be detected between 10 and 20 DPN, and eventually extend and join the initial lymphatics (see Figure 7B).

Suzuki, 1982, reported that a testis-like blood network is present in most of the epididymis (segments 2–9) with intertubular vessels connected by perpendicular micro-vessels in a rope ladder organization. Distinctly, in the most proximal (IS/S1) and distal (S10) segments of the epididymis, we found micro-vessels encircling each tubule in a spider web organization with capillaries that penetrate the epithelial sublayer. Because the IS and S10 compartments of the epididymis have distinct embryonic origins and are functionally different (Abe et al., 1984; Johnston et al., 2007), it is difficult to explain this common vascular organization. We hypothesize that the common set-up of these two territories where capillaries penetrate the sub-epithelial layer could be associated with their immune function. The S10 compartment must monitor ascending pathogens and maintain self-tolerance to sperm antigens that accumulate in this storage part of the organ. Similarly, the S1 compartment is the final gatekeeper preventing ascending pathogens from reaching the immune-privileged seminiferous tubules, while at the same time it must be tolerant to sperm-specific antigens entering the epididymal tubule and considered non-self. Therefore, in both compartments, a spider web blood network can be expected to allow immune cells easy access to the epididymal tubule. This should be particularly true for the S1 segment because, independently of infectious situations, it is constantly solicited in the adult animal by the arrival of sperm antigens and by its function of reabsorption of Sertolian fluid requiring very permeable vessels. This hypervascularization of the S1 epididymal compartment is further supported by the fact that it has been shown to be PLVAP<sup>pos</sup>, a marker of fenestrated vessels allowing high trans-endothelial transport (Guo et al., 2016). In other contexts, PLVAP has been shown to be required for fenestron biogenesis and organization, controlling fenestron permeability, angiogenesis, as well as leukocyte, antigen migration, and immomodulatory processes (Rantakari et al., 2015; Amersfoort et al., 2022). In addition, the presence of VEGFR2, VEGFR3, and VEGF-A in the microvascularization of the epididymis S1 compartment argues for active angiogenesis (Bussolati et al., 2003). This has been corroborated by the recent single-cell data (Rinaldi et al., 2020) showing that expression of angiogenic markers such as those mentioned above occurs in endothelial cells of
<figure><img src="image_7.png" /><figcaption>**Figure 5.** VEGF-A, VEGF-C, and VEGF-D levels vary during postnatal epididymal development. Panel A shows the expression profile obtained in total extracts of epididymal proteins at different stages of development. The profiles obtained are presented in the following order: VEGFR3, VEGF-A, VEGF-C, and VEGF-D. GAPDH was used for normalization. Panel B shows the quantification of epididymal proteins extracted from six mice. Quantification of VEGFR3 is shown in black, and VEGF-A, VEGF-C, and VEGF-D are shown in red, blue, and green, respectively. The Kruskal-Wallis test and Dunn's posttest correction were used to determine statistical significance (***p<0.0001, **p<0.001, *p<0.01). Panel C shows the comparison of different isoform profiles. Covariance is assessed by a two-way analysis of variance.</figcaption></figure>

The online version of this article includes the following figure supplement(s) for figure 5:

**Figure supplement 1.** Quantification of all forms of the hemangiogenic and lymphangiogenic ligands VEGF-A, VEGF-C, and VEGF-D.

the caput mouse epididymis. While active angiogenesis during postnatal epididymal development makes sense because it accompanies postnatal growth of the organ, it is more surprising that angiogenesis is still active in the adult IS epididymal compartment. In adult tissues, angiogenesis is most commonly associated with vascular remodeling in tissue repair processes and tumor progression (Komi et al., 2020). To date, self-renewal of the epididymal epithelium is collectively assumed to be rather slow, a situation that should not be associated with high angiogenic activity. However, very recent data could modify this general assumption because it has been reported that the basal cells of the epididymis share common properties with adult stem cells (Dufresne et al., 2022). This could
<figure><img src="image_8.png" /><figcaption>**Figure 6.** VEGF receptors and ligands in the mouse epididymis. Panel A (upper row) shows images (Zeiss Axio-Imager) of the expression of the VEGFR2:GFP transgene (green) (used with permission from Prof. A. Medvinsky) and its ligand VEGF-A (red) obtained from paraffin sections of adult epididymis. The lower row shows confocal views (SP8, Leica) of the same VEGF-A labeling. Panel B shows in (1) a mosaic view (AxioVision scanner) of an epididymis after immunolabeling with VEGF-C (green) and VEGF-D (cyan). Photographs 2–9 were taken with a confocal (SP8, Leica) at the level of *Figure 6 continued on next page*</figcaption></figure>
Figure 6 continued

the initial segment (IS)/S1 (2), caput (3), corpus (4), and cauda (5). Notable differences in expression of these two ligands are shown in photographs 6–8, which respectively concerns the efferent duct and the IS/S1 boundary, the capsule, and a septum. Photograph 9 is a negative control. Pearson correlation and Mander's co-occurrence were used to analyze the relation between the two ligands with the IMARIS co-localization module. The values shown represent means with standard error of the mean.

support the idea that the epithelium of the adult epididymis may be in a perpetual process of self-renewal that could be promoted by the quasi-inflammatory situation inherent to immune stimulation mediated by sperm antigens.

It is interesting to note that during postnatal ontogeny of the epididymis, PLVAP expression increases in the S1 territory but decreases in the other regions of the organ as early as 15 DPN.

<figure><img src="image_9.png" /><figcaption>**Figure 7.** The septa between epididymal segments are intimately linked to lymphatic vasculature. Panel A shows the intimate connection between the septa and the lymphatic vasculature in a region of interest of the caput shown in **Figure 2**. Surface rendering of the lymphatics was performed as described previously using IMARIS software. Panel B (top panel) shows a representative three-dimensional (3D) view of the cauda region at the S8–S9/S10 septa during postnatal epididymal development after surface rendering. Images were obtained using the light-sheet ultramicroscope with a 20× objective lens. Putative LYVE1+ lymphatic precursors are in cyan and VEGFR3:YFP lymphatics are in green, whereas PLVAP+ blood vessels are in red. The clipping plane is positioned at the septum level and the yellow arrow-cursor is directed to the S10. The lower panel shows our interpretation of the events observed concomitantly with a progression via sprouting lymphangiogenesis (green to yellow arrows) of peripheral lymphatics that progress and radiate into the organ at the level of the septum, from the anterior side (adjacent to the testis) to the posterior side of the epididymis and via lymphangio-vasculogenesis. It develops in four steps: (1) swarming of precursors mainly LYVE+ (cyan) and VEGFR3+ (green); (2) a stage of LECP expansion; (3) grouping or clustering step; (4) then assembly and fusion with the lymphatics located at the septum from lymphangiogenesis (green to yellow arrow). For clarity, PLVAP+ blood vessels have not been represented here (A=anterior side adjacent to the testis, P=posterior side, S=septum, LN = lymph node). Scale bar = 30 µm.</figcaption></figure>
<figure><img src="image_10.png" /><figcaption>**Figure 8.** The proposed model of the lympho-septa of the epididymis. The left scheme presents our view of the lymphatic vasculature at the level of a 'checkpoint' septum. The large collectors enter at the septum into the epididymis and radiate within the adjacent segments at the interstitial level, where initial lymphatics are found. The septum where the tubule crosses from one segment to another one would be the most 'monitored' site and therefore the richest in lymphatics. The close association of lymphatics with the septa creates both a physical and immunogenic barrier to preserve the organ from ascending infections and thus limited uncontrolled progression of pathogens protecting the epididymis and ultimately the testis from orchitis deleterious to male fertility. e-CLV=external collector lymphatic vessel; i-CLV=internal collector lymphatic vessel; ILV = initial lymphatic vessel; LN = lymph node; S=segment.</figcaption></figure>

Whether this is due to blood vessels loss or to PLVAP<sup>pos</sup> permeable fenestrons becoming PLVAP<sup>neg</sup>, less permeable vessels will require further study. In the mouse model, it has been shown that the development of the S1 territory occurs between 15 and 20 DPN concomitantly with the arrival of testicular fluid (Jegou et al., 1982), which has been shown to stimulate the Raf/Mek/Erk pathway known to be activated by VEGFs (Xu et al., 2016). Therefore, it is possible to consider that VEGF-A could be one of the lumicrine factors triggering angiogenesis and differentiation of S1/IS. In agreement with this, it has been reported that mice overexpressing VEGF-A showed dilation of the caput and corpus epididymis, with enlarged and more permeable blood vessels (Korpelainen et al., 1998). More recently, Pawlak and Caron, 2020, in their review have suggested that such vessels could be hybrid lymphatics also recently called 'mosaic vessels' (Lampejo et al., 2023) expressing both blood and lymphatic markers. Such hybrid vessels have been found in distinct vascular beds, including high endothelial venules in secondary lymphoid organs, liver sinusoidal endothelial cells, Schlemm's canal endothelial cells of the eye, the ascending vasa recta (AVR) of the kidney inner medulla region, and the remodeled spiral arteries of the placenta decidua (Gola et al., 2021; Kenig-Kozlovsky et al., 2018; Kim et al., 2017; Pawlak et al., 2019; Russell et al., 2019; Stamataki and Swadling, 2020) where they perform specialized exchange functions. As shown in Figure 9, these structures share many markers, including those tested in the present work, suggesting that some of the assumed blood vessels found in the S1 epididymal territory might be hybrid vessels. Given the close mesonephric origin and similar fluid reabsorption function of the IS/S1 epididymal territory and the AVR, it is possible that the two structures share such hybrid vessels. Tie2 expression, a marker of all currently described hybrid lymphatics (Kenig-Kozlovsky et al., 2018, see Figure 9), supports this hypothesis as it localizes to the S1 peritubular microvascularization together with VEGFR3 and PVLAP (not shown).

In conclusion, based on our extensive study of the blood and lymphatic vascularization of the epididymis in mice, we propose that the epididymis presents two types of lymphatic vessels. On the one hand, there are peritubular PLVAP<sup>pos</sup> hybrid vessels, very permeable, particularly well represented
<figure><img src="image_11.png" /><figcaption>**Figure 9.** Expression of hybrid lymphatic markers in the mouse epididymis. This table compares markers associated with hybrid lymphatic vessels described in the literature (for a review, see *Pawlak and Caron, 2020*) with (column 2) their expression at different segments of the epididymis (source: Mammalian Reproductive Genetic; https://www.mrgd.org). Multiple results for the same gene correspond to the use of different Oligosets in the microarrays used by the MRG. Column 3 summarizes our present results at the initial segment level. The asterisk refers to the publication by *Korpelainen et al., 1998*.</figcaption></figure>
<figure><img src="image_12.png" /><figcaption>**Figure 10.** Schematic summary of the expansion of the conventional and hybrid lymphatic vasculature during postnatal development of the murine epididymis.</figcaption></figure>

in the S1 compartment that, according to our hypothesis, participate both in the reabsorption function and in the immune surveillance of this territory. On the other hand, we have more conventional lymphatics mainly in the interstitial compartment. The various segments of the epididymis are drained by initial and collecting lymphatic vessels, the latter being closely associated with the septa and then connecting via the PP to the most proximal LN. Three checkpoints are represented by the IS/S2-S3, S6/S7, and S8-9/S10 septa. It is not yet clear to us that the lymphatic drainage of the S8-9/S10 compartment goes to the PP because it could also be drained by the lymphatic circuit monitoring the vas deferens. This detailed knowledge of the blood and lymphatic circuits of the epididymis highlights how open the epididymis is to the systemic compartment and how different its immune/inflammatory regulation is expected to be from that of the testis, where lymphatic collectors are restricted to the conjunctiva and never reach the interstitial compartment (Hirai et al., 2012; Svingen et al., 2012). A detailed schematic of epididymal vascular ontogeny in the mouse is proposed in **Figure 10**. If these observations in the mouse model are somehow translatable to the human epididymis, they could help to understand the complex inflammatory and immune contexts of posttesticular sperm maturation and their impact on male fertility.

# Materials and methods

## Mice

In this study, we have used 3- to 6-month-old male C57BL/6 mice (Janvier Labs, France); the male VEGFR3-YFP mice came from the J.L. Thomas laboratory (UPMC-Inserm, Paris, France). The
Vegfr3TYFP construct was generated from a BAC clone (RP23-65D23, Chori BACPAC resources) containing a 238kb geneomic fragment spanning the mouse Vegfr3 locus (Calvo et al., 2011). All animals were housed according to institutional guideline with a 12 hr photoperiod and food and water available ad libitum. Measures were taken to keep animal suffering to a minimum. The Auvergne Animal Experiment Ethics Committee (C2E2A) and the French Ministry approved all the following procedures for research (APAFIS authorization #12376-2017112913113962V3).

## Antibodies

The primary and secondary antibodies used in clearing methods are indicated in **Figure 11B**. Additional antibodies used are: rat anti-mouse VEGFR3 (#NB110-61018, Novus Biologicals, LLC, Bio-Techne SAS, France); goat anti-mouse VEGF-A (#AF-493-SP, Novus Biologicals, LLC, Bio-Techne SAS, France); rabbit anti-mouse VEGF-C (#NB110-61022, Novus Biologicals, LLC, Bio-Techne SAS, France); goat anti-VEGF-D (#AF469, R&D Systems, Bio-Techne SAS, France); and rabbit anti-GAPDH (SAB2108668, Sigma-Aldrich Chimie Sarl, France).

## Western blot analysis

Soluble epididymal proteins were prepared, as described previously (Chorfa et al., 2021), separated with sodium dodecyl sulfate-polyacrylamide gel electrophoresis (12% gel), and transferred to polyvinylidene fluoride membranes (Hybond ECL, Amersham Biosciences, Germany). Primary antibodies against the following protein were used: VEGFR3, VEGF-A, VEGF-C, VEGF-D, and GAPDH (served as a loading control). The appropriate HRP-conjugated secondary antibodies, goat anti-rabbit IgG or goat anti-mouse IgG (Abliance, France), were used to visualize the protein bands. Immunoreactive bands were detected by chemiluminescence (Clarity Western ECL Substrate Bio-Rad, France) using the ChemiDoc MP imaging system (Bio-Rad). Protein quantification was performed with ImageJ software. Protein amounts are expressed as relative values to the GAPDH.

## Whole-mount multiplex immunolabeling and clearing procedure

Epididymides were collected from mice after cervical dislocation and fixed in 4% paraformaldehyde for a period of time determined by their size (see **Figure 11**). Then, they were subjected to saturation/permeabilization for 1–5 days depending on their size in PBSGTS solution (1x saline buffer phosphate containing 2% gelatin, 0.5% Triton X-100, and saponin 1 µg/ml) on a rotary shaker (100 rpm) at room temperature. Incubation with primary and secondary antibodies was performed under rotation at 37°C in the same buffer according to the schedule described in **Figure 11A**. Concentrations of primary and secondary antibodies used are reported in **Figure 11B**. The immunolabeled epididymides were then embedded in a 1.5% agarose cube in 0.5x Tris-Acetate-EDTA buffer and cleared according to the 3DISCO method (Belle et al., 2014). The samples were stored in a dark place at room temperature until observation.

## Immunofluorescence

Paraffin-embedded 5 µm sections of epididymides were subjected to heat-induced antigen recovery and then permeabilized for 30 min at room temperature with PBS supplemented with 0.3% Tween 20 and saponin (1 mg/mL). The sections were incubated for 1 hr at room temperature in the same solution with 10% serum before overnight incubation with the primary antibody at 4°C. Their detection was then carried out with appropriate secondary antibodies conjugated to Alexa Fluor 488, Alexa Fluor 555, or Alexa Fluor 647 (Invitrogen, Thermo Fisher Scientific, France). Nuclei were stained with Hoechst solution (1 µg/µL). The sections were embedded in Mowiol 4-88 (Sigma-Aldrich Chimie Sarl, France) and stored at 4°C until observation.

## Image acquisition

### Macroscopic microscopy

Macroscopic observations of endogenous VEGFR3:YFP transgene expression of the epididymis and testis were performed using a Leica binocular magnifier (with a 1x objective).

### Confocal microscopy

Multiplex immunofluorescence acquisition was performed using a SP8 confocal laser scanning microscope equipped with a Plan Apo λ 40× Oil objective (Leica, Germany). The following parameters
A

1 Blood + 3 Lymphatic markers

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

4 Lymphatic markers

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

B

3DISCO antibodies

<table><thead><tr><th>Primary antibodies</th><th>Host species</th><th>Working concentration</th><th>working vomume</th><th>Distributors</th><th>Catalog #</th></tr></thead><tbody><tr><td>Anti-GFP</td><td>RABBIT</td><td>0.6 to 4µg/ml</td><td rowspan="6">1.5 to 2ml</td><td>ThermoFisher</td><td>A-11122</td></tr><tr><td>Anti-mouse LYVE-1</td><td>RAT</td><td>0.5 to 1µg/ml</td><td>NovusBio</td><td>NBP1-43411</td></tr><tr><td>Anti-mouse MECA32</td><td>RAT</td><td>0.03 to 0.06µg/ml</td><td>BD Biosciences</td><td>550563</td></tr><tr><td>Anti-mouse,PDPN</td><td>HAMSTER</td><td>0.2 to 0.5µg/ml</td><td>ReliaTech</td><td>103-M40</td></tr><tr><td>Anti-mouse LYVE-1 biotinylated</td><td>GOAT</td><td>0.2 to 0.4µg/ml</td><td>R&amp;D System</td><td>BAF2125</td></tr><tr><td>Anti-Human PROX1 biotinylated</td><td>GOAT</td><td>0.2 to 0.4µg/ml</td><td>R&amp;D System</td><td>BAF2727</td></tr><tr><th>Secondary antibodies</th><th>Working concentration</th><th>working vomume</th><th>Distributors</th><th>Catalog #</th></tr><tr><td>Alexa Fluor 488 F(ab') goat anti-rabbit IgG (H+L)</td><td>0.2ng/ml</td><td rowspan="4">1.5 to 2ml</td><td>Invitrogen</td><td>A11070</td></tr><tr><td>Alexa Fluor 555 F(ab') goat anti-rat IgG (H+L)</td><td>0.2ng/ml</td><td>Invitrogen</td><td>A21434</td></tr><tr><td>Alexa Fluor 647 F(ab') goat anti-hamster IgG (H+L)</td><td>0.2ng/ml</td><td>Invitrogen</td><td>A21451</td></tr><tr><td>Chromeo™ 494 conjugated to streptavidin</td><td>0.2ng/ml</td><td>Active Motif</td><td>15113</td></tr></tbody></table>

**Figure 11.** Multiplex labeling workflow used to visualize the blood and lymphatic vasculature of the mouse epididymis clarified by the 3DISCO method. Panel A presents the multiplex labeling schedule for blood and lymphatic immunodetection. Each square corresponds to 24 hr except for the fixation and streptavidin/biotin treatment, where the time is indicated. Panel B provides detailed information regarding the various antibodies used in the course of the study.
were used: pinhole size of 1 airy unit (AU); Z-step: 0.5 µm. Four acquisition sequences were executed: (1) for Alexa 647 detection (638 nm laser 659–698 nm window, using a photomultiplier tubes (PMT) detector at 700 V gain); (2) for Alexa 488 detection (488 nm laser; 493–541 nm window, using a PMT detector at 565 V gain) and for Chromeo 494 detection (488 nm laser; 660–700 nm window, using a PMT detector with a gain of 800 V); (3) for Alexa 555 detection (laser 552 nm; 560–570 nm window, using a PMT detector with a gain of 700 V); and (4) for Hoechst detection (laser 405 nm; 423–518 nm window, using a Hybrid detector with 10% gain).

## Light-sheet ultramicroscopy

3D imaging was achieved with an ultramicroscope (LaVision BioTec Miltenyi, Germany) using Imspec- torPro software. The images were obtained with either an MI PLAN 2x/NA 0.5, an MI PLAN 4x/0.35 objectives (MVPLAPO, Olympus), or a 20x/0.95 objective (Leica). Each sample was placed in a 100% quartz imaging reservoir filled with dibenzyl ether and illuminated from the side by the laser sheet. Images were acquired with a Andor Neo SCMOS CCD camera (2160x2,160 pixels, LaVision BioTec). All 3D acquisitions were performed with a step size between each image was fixed at 1 µm (1.6 µm/ pixel). Lasers at 488, 561, and 635 nm were used to obtain images using two light sources for the caput, corpus, or cauda epididymis.

## Image processing

Ultramicroscopy datasets were uploaded to IMARIS version 9.7 (Bitplane, Oxford Instruments, England). The stacks were converted to IMARIS files (.ims) and 3D visualizations of z-stack images were generated using the volume rendering function. The vessels were transformed into surface rendering and the volume calculated automatically by IMARIS. Vessel densities were estimated as the ratio of the vessel volume to the total volume of the region. The videos were made IMARIS without deconvolution of the image.

## Statistical analysis

All experiments were repeated at least three times and representative images are shown. All statistical analyses were performed with Prism software (GraphPad, USA). The Mann-Whitney test was used to compare two groups. The Kruskal-Wallis test followed by Dunn's posttest was used to compare more than two groups. A two-way analysis of variance test with a Bonferroni posttest was used to analyze the covariance of ligands (VEGF-A, VEGF-C, and VEGF-D) with the VEGFR3 receptor. The co-localiza- tion of the ligands VEGF-C and VEGF-D was evaluated by Pearson correlation and Mander's overlap coefficient with IMARIS co-localization modules (Bitplane). The values for co-localization analysis represent the mean ± standard error of the mean.

## Acknowledgements

The authors would like to thank the CLIC confocal microscopy platform (GReD Institute) for their expert assistance during this study.

## Additional information

### Funding

No external funding was received for this work.

### Author contributions

Christelle Damon-Soubeyrand, Areski Chorfa, Chantal Goubely, Nelly Pirot, Laurence Piboin-Fragner, Caroline Vachias, Stephanie Bravard, Ayhan Kocer, Investigation; Antonino Bongiovanni, Software, Investigation; Luc Pardanaud, Rachel Guiton, Jean-Leon Thomas, Fabrice Saez, Validation; Meryem Tardivel, Validation, Investigation; Joël R Drevet, Conceptualization, Supervision, Validation, Writing – original draft, Writing – review and editing; Joelle Henry-Berger, Data curation, Formal analysis, Investigation, Writing – original draft, Writing – review and editing
## Author ORCIDs

Antonino Bongiovanni http://orcid.org/0000-0001-9876-2940

Joël R Drevet http://orcid.org/0000-0003-3077-6558

Joelle Henry-Berger http://orcid.org/0000-0002-5107-8663

## Ethics

The Auvergne Animal Experiment Ethics Committee (C2E2A) and the French Ministry approved all the following procedures for research (APAFIS authorization # 12376-2017112913113962 V3).

## Decision letter and Author response

Decision letter https://doi.org/10.7554/eLife.82748.sa1

Author response https://doi.org/10.7554/eLife.82748.sa2

## Additional files

### Supplementary files

* MDAR checklist

* Supplementary file 1. Describes all the data source legends in (*.ims) (Imaris viewer) and (*.czi) (Zeiss AxioScan viewer) format that can be found in the Biolmage Archive Biostudies EBI repository.

* Supplementary file 2. Is an Excel file listing all the data sources for the figures provided.

### Data availability

Data sources are available at https://www.ebi.ac.uk/biostudies/bioimages/studies/S-BIAD618. Data source legends are provided as Supplementary files 1 & 2. These files list the source data files available via the Biolmage Archive repository BioStudies accession number S-BIAD618.

The following dataset was generated:

<table><thead><tr><td>Author(s)</td><td>Year</td><td>Dataset title</td><td>Dataset URL</td><td>Database and Identifier</td></tr></thead><tbody><tr><td>Bongiovanni A, Damon-Soubeyrand C, Chorfa A, Goubely C, Pirot N, Pardanaud L, Pibouin-Fragner L, Vachias C, Bravard S, Guiton R, Thomas J, Saez F, Kocer A, Tardivel M, Drevet JR, Henry-Berger J</td><td>2023</td><td>Three-dimensional imaging of vascular development in the mouse epididymis: a prerequisite to better understand the post-testicular immune context of spermatozoa</td><td>https://www.ebi.ac.uk/biostudies/bioimages/studies/S-BIAD618</td><td>EBI Biolmage Archive, SBIAD618</td></tr></tbody></table>

## References

* Abe K, Takano H, Ito T. 1984. Microvasculature of the mouse Epididymis, with special reference to Fenestrated Capillaries localized in the initial segment. *The Anatomical Record* **209**:209–218. DOI: https://doi.org/10.1002/ar.1092090208, PMID: 6465531

* Amersfoort J, Eelen G, Carmeliet P. 2022. Immunomodulation by endothelial cells - Partnering up with the immune system. *Nature Reviews. Immunology* **22**:576–588. DOI: https://doi.org/10.1038/s41577-022-00694-4, PMID: 35288707

* Baluk P, McDonald DM. 2008. Markers for microscopic imaging of Lymphangiogenesis and angiogenesis. *Annals of the New York Academy of Sciences* **1131**:1–12. DOI: https://doi.org/10.1196/annals.1413.001, PMID: 18519955

* Banerji S, Ni J, Wang SX, Clasper S, Su J, Tammi R, Jones M, Jackson DG. 1999. Lyve1, a new Homologue of the Cd44 Lycoprotein, is a lymph-specific receptor for Hyaluronan. *The Journal of Cell Biology* **144**:789–801. DOI: https://doi.org/10.1083/jcb.144.4.789, PMID: 10037799

* Belle M, Godefroy D, Dominici C, Heitz-Marchaland C, Zelina P, Hellal F, Bradke F, Chédotal A. 2014. A simple method for 3D analysis of Immunolabeled axonal tracts in a transparent nervous system. *Cell Reports* **9**:1191–1201. DOI: https://doi.org/10.1016/j.celrep.2014.10.037, PMID: 25456121

* Breiteneder-Geleff S, Soleiman A, Kowalski H, Horvat R, Amann G, Kriehuber E, Diem K, Weninger W, Tschachler E, Alitalo K, Kerjaschki D. 1999. Angiosarcomas Express mixed endothelial phenotypes of blood and Lymphatic Capillaries: Podoplanin as a specific marker for Lymphatic Endothelium. *The American Journal of Pathology* **154**:385–394. DOI: https://doi.org/10.1016/S0002-9440(10)65285-6, PMID: 10027397
Britan A, Maffre V, Tone S, Drevet JR. 2006. Quantitative and spatial differences in the expression of Tryptophan-Metabolizing enzymes in mouse Epididymis. *Cell and Tissue Research* **324**:301–310. DOI: https://doi.org/10.1007/s00441-005-0151-7, PMID: 16450123

Bussolati B, Deambrosia I, Russo S, Deregibus MC, Camussi G. 2003. Altered angiogenesis and survival in human tumor-derived endothelial cells. *FASEB Journal* **17**:1159–1161. DOI: https://doi.org/10.1096/fj.02-0557fje, PMID: 12709414

Calvo C-F, Fontaine RH, Soueid J, Tammela T, Makinen T, Alfaro-Cervello C, Bonnaud F, Miguez A, Benhaim L, Xu Y, Barallobre M-J, Moutkine I, Lyttikä J, Tatlisumak T, Pytowski B, Zalc B, Richardson W, Kessaris N, Garcia-Verdugo JM, Alitalo K, et al. 2011. Vascular endothelial growth factor receptor 3 directly regulates murine Neurogenesis. *Genes & Development* **25**:831–844. DOI: https://doi.org/10.1101/gad.615311, PMID: 21498572

Chorfa A, Goubely C, Henry-Berger J, Guiton R, Drevet JR, Saez F. 2021. Identification of Arvicola Terrestris Scherman Sperm antigens for immune contraceptive purposes. *International Journal of Molecular Sciences* **22**:9965. DOI: https://doi.org/10.3390/ijms22189965, PMID: 34576131

Da Silva N, Cortez-Retamozo V, Reinecker H-C, Wildgruber M, Hill E, Brown D, Swirski FK, Pittet MJ, Breton S. 2011. A dense network of Dendritic cells Populates the murine Epididymis. *REPRODUCTION* **141**:653–663. DOI: https://doi.org/10.1530/REP-10-0493

Dohle GR, Colpi GM, Hargreave TB, Papp GK, Jungwirth A, Weidner W, EAU Working Group on Male Infertility. 2005. EAU guidelines on male infertility. *European Urology* **48**:703–711. DOI: https://doi.org/10.1016/j.eururo.2005.06.002, PMID: 16005562

Dufresne J, Gregory M, Pinel L, Cyr DG. 2022. Differential gene expression and hallmarks of Stemness in epithelial cells of the developing rat Epididymis. *Cell and Tissue Research* **389**:327–349. DOI: https://doi.org/10.1007/s00441-022-03634-9, PMID: 35590013

Ertürk A, Becker K, Jährling N, Mauch CP, Hojer CD, Egen JG, Hellal F, Bradke F, Sheng M, Dodt HU. 2012. Three-dimensional imaging of solvent-cleared organs using 3Disco. *Nature Protocols* **7**:1983–1995. DOI: https://doi.org/10.1038/nprot.2012.119, PMID: 23060243

Gola A, Dorrington MG, Speranza E, Sala C, Shih RM, Radtke AJ, Wong HS, Baptista AP, Hernandez JM, Castellani G, Fraser IDC, Germain RN. 2021. Commensal-driven immune Zonation of the liver promotes host defense. *Nature* **589**:131–136. DOI: https://doi.org/10.1038/s41586-020-2977-2

Goodnow CC. 1996. Balancing immunity and tolerance: deleting and tuning lymphocyte Repertoires. *PNAS* **93**:2264–2271. DOI: https://doi.org/10.1073/pnas.93.6.2264

Gruffaz J. 1984. Anatomy and physiology of the Lymphatic system. *Phlebologie* **37**:189–193 PMID: 6473526.

Guiton R, Henry-Berger J, Drevet JR. 2013. The Immunobiology of the mammalian Epididymis: the black box is now open! basic Clin Basic and Clinical Andrology **23**:1–10. DOI: https://doi.org/10.1186/2051-4190-23-8

Guo L, Zhang H, Hou Y, Wei T, Liu J. 2016. Plasmalemma Vesicle-associated protein: a crucial component of vascular homeostasis. *Experimental and Therapeutic Medicine* **12**:1639–1644. DOI: https://doi.org/10.3892/etm.2016.3557

Gutierrez-Miranda L, Yaniv K. 2020. Cellular origins of the Lymphatic Endothelium: implications for cancer Lymphangiogenesis. *Frontiers in Physiology* **11**:577584. DOI: https://doi.org/10.3389/fphys.2020.577584, PMID: 33071831

Hamrah P, Chen L, Zhang Q, Dana MR. 2003. Novel expression of vascular endothelial growth factor receptor (VEGFR)-3 and VEGF-C on corneal Dendritic cells. *The American Journal of Pathology* **163**:57–68. DOI: https://doi.org/10.1016/S0002-9440(10)63630-9, PMID: 12819011

Hedger MP. 2011. Immunophysiology and pathology of inflammation in the Testis and Epididymis. *Journal of Andrology* **32**:625–640. DOI: https://doi.org/10.2164/jandrol.111.012989, PMID: 21764900

Hirai S, Naito M, Terayama H, Ning Q, Miura M, Shirakami G, Itoh M. 2010. Difference in abundance of blood and Lymphatic Capillaries in the murine Epididymis. *Medical Molecular Morphology* **43**:37–42. DOI: https://doi.org/10.1007/s00795-009-0473-8, PMID: 20340004

Hirai S., Naito M, Terayama H, Qu N, Kuerban M, Musha M, Ikeda A, Miura M, Itoh M. 2012. The origin of Lymphatic Capillaries in murine Testes. *Journal of Andrology* **33**:745–751. DOI: https://doi.org/10.2164/jandrol.111.015156

Itoh M, Li XQ, Yano A, Xie Q, Takeuchi Y. 1998. Patterns of Efferent Lymphatics of the mouse Testis. *Journal of Andrology* **19**:466–472 PMID: 9733149.

Jafree DJ, Moulding D, Kolatsi-Joannou M, Perretta Tejedor N, Price KL, Milmoe NJ, Walsh CL, Correra RM, Winyard PJ, Harris PC, Ruhrberg C, Walker-Samuel S, Riley PR, Woolf AS, Scambler PJ, Long DA. 2019. Spatiotemporal Dynamics and heterogeneity of renal Lymphatics in mammalian development and cystic kidney disease. *eLife* **8**:e48183. DOI: https://doi.org/10.7554/eLife.48183, PMID: 31808745

Jafree DJ, Long DA, Scambler PJ, Ruhrberg C. 2021. Mechanisms and cell lineages in Lymphatic vascular development. *Angiogenesis* **24**:271–288. DOI: https://doi.org/10.1007/s10456-021-09784-8, PMID: 33825109

Jegou B, Le Gac F, de Kretser DM. 1982. Seminiferous tubule fluid and interstitial fluid production I effects of age and hormonal regulation in immature rats. *Biology of Reproduction* **27**:590–595. DOI: https://doi.org/10.1095/biolreprod27.3.590, PMID: 6814545

Johnston DS, Jelinsky SA, Bang HJ, DiCandeloro P, Wilson E, Kopf GS, Turner TT. 2005. The mouse Epididymal Transcriptome: transcriptional profiling of segmental gene expression in the Epididymis. *Biology of Reproduction* **73**:404–413. DOI: https://doi.org/10.1095/biolreprod.105.039719, PMID: 15878890
Johnston DS, Turner TT, Finger JN, Owtscharuk TL, Kopf GS, Jelinsky SA. 2007. Identification of Epididymis-specific transcripts in the mouse and rat by transcriptional profiling. *Asian Journal of Andrology* **9**:522–527.
DOI: https://doi.org/10.1111/j.1745-7262.2007.00317.x, PMID: 17589790

Jrad-Lamine A, Henry-Berger J, Gourbeyre P, Damon-Soubeyrand C, Lenoir A, Combaret L, Saez F, Kocer A, Tone S, Fuchs D, Zhu W, Oefner PJ, Munn DH. Mellor AL, Gharbi N, Cadet R, Aitken RJ, Drevet JR. 2011. Deficient Tryptophan catabolism along the Kynurenine pathway reveals that the Epididymis is in a unique Tolerogenic state. *Journal of Biological Chemistry* **286**:8030–8042. DOI: https://doi.org/10.1074/jbc.M110.172114

Jrad-Lamine A, Henry-Berger J, Damon-Soubeyrand C, Saez F, Kocer A, Janny L, Pons-Rejraji H, Munn DH, Mellor AL, Gharbi N, Cadet R, Guiton R, Aitken RJ, Drevet JR. 2013. Indoleamine 2,3-Dioxygenase 1 (Ido1) is involved in the control of mouse Caput Epididymis immune environment. *PLOS ONE* **8**:e66494. DOI: https://doi.org/10.1371/journal.pone.0066494, PMID: 23840489

Jungwirth A, Giwercman A, Tournaye H, Diemer T, Kopa Z, Dohle G, Krausz C, European Association of Urology Working Group on Male Infertility. 2012. European Association of Urology guidelines on male infertility: the 2012 update. *European Urology* **62**:324–332. DOI: https://doi.org/10.1016/j.eururo.2012.04.048, PMID: 22591628

Kaipainen A, Korhonen J, Mustonen T, van Hinsbergh VW, Fang GH, Dumont D, Breitman M, Alitalo K. 1995. Expression of the FMS-like tyrosine kinase 4 Gene becomes restricted to Lymphatic Endothelium during development. *PNAS* **92**:3566–3570. DOI: https://doi.org/10.1073/pnas.92.8.3566, PMID: 7724599

Kato S, Shimoda H, Ji R-C, Miura M. 2006. Lymphangiogenesis and expression of specific molecules as Lymphatic endothelial cell markers. *Anatomical Science International* **81**:71–83. DOI: https://doi.org/10.1111/j.1447-073X.2006.00142.x

Kenig-Kozlovsky Y, Scott RP, Onay T, Carota IA, Thomson BR, Gil HJ, Ramirez V, Yamaguchi S, Tanna CE, Heinen S, Wu C, Stan RV, Klein JD, Sands JM, Oliver G, Quaggin SE. 2018. Ascending vasa recta are Angiopoietin/Tie2-dependent lymphatic-like vessels. *Journal of the American Society of Nephrology* **29**:1097–1107. DOI: https://doi.org/10.1681/ASN.2017090962, PMID: 29237738

Kim J, Park D-Y, Bae H, Park DY, Kim D, Lee C, Song S, Chung T-Y, Lim DH, Kubota Y, Hong Y-K, He Y, Augustin HG, Oliver G, Koh GY. 2017. Impaired Angiopoietin/Tie2 signaling compromises Schlemm's Canal integrity and induces glaucoma. *Journal of Clinical Investigation* **127**:3877–3896. DOI: https://doi.org/10.1172/JCI94668

Kivelä R, Salmela I, Nguyen YH, Petrova TV, Koistinen HA, Wiener Z, Alitalo K. 2016. The transcription factor Prox1 is essential for satellite cell differentiation and muscle fibre-type regulation. *Nature Communications* **7**:13124. DOI: https://doi.org/10.1038/ncomms13124

Komi DEA, Khomtchouk K, Santa Maria PL. 2020. A review of the contribution of mast cells in wound healing: involved molecular and cellular mechanisms. *Clinical Reviews in Allergy & Immunology* **58**:298–312. DOI: https://doi.org/10.1007/s12016-019-08729-w, PMID: 30729428

Korpelainen EI, Karkkainen MJ, Tenhunen A, Lakso M, Rauvala H, Vierula M, Parvinen M, Alitalo K. 1998. Overexpression of VEGF in Testis and Epididymis causes infertility in transgenic mice: evidence for Nonendothelial targets for VEGF. *The Journal of Cell Biology* **143**:1705–1712. DOI: https://doi.org/10.1083/jcb.143.6.1705, PMID: 9852161

Lampejo AO, Ghavimi SAA, Hägerling R, Agarwal S, Murfee WL. 2023. Lymphatic/blood vessel plasticity: motivation for a future research area based on present and past observations. *American Journal of Physiology. Heart and Circulatory Physiology* **324**:H109–H121. DOI: https://doi.org/10.1152/ajpheart.00612.2022, PMID: 36459445

Lioux G, Liu X, Temiño S, Oxendine M, Ayala E, Ortega S, Kelly RG, Oliver G, Torres M. 2020. A second heart field-derived Vasculogenic niche contributes to cardiac Lymphatics. *Developmental Cell* **52**:350–363. DOI: https://doi.org/10.1016/j.devcel.2019.12.006, PMID: 31928974

Mäkinen T, Veikkola T, Mustjoki S, Karpanen T, Catimel B, Nice EC, Wise L, Mercer A, Kowalski H, Kerjaschki D, Stacker SA, Achen MG, Alitalo K. 2001. Isolated Lymphatic endothelial cells Transduce growth, survival and migratory signals via the VEGF-C/D receptor VEGFR-3. *The EMBO Journal* **20**:4762–4773. DOI: https://doi.org/10.1093/emboj/20.17.4762, PMID: 11532940

McDonald SW, Scothorne RJ. 1988. The Lymphatic drainage of the Epididymis and of the Ductus Deferens of the rat, with reference to the immune response to Vasectomy. *Journal of Anatomy* **158**:57–64. DOI: 3225225.

Meinhardt A, Hedger MP. 2011. Immunological, Paracrine and endocrine aspects of Testicular immune privilege. *Molecular and Cellular Endocrinology* **335**:60–68. DOI: https://doi.org/10.1016/j.mce.2010.03.022

Mouta Carreira C, Nasser SM, di Tomaso E, Padera TP, Boucher Y, Tomarev SI, Jain RK. 2001. Lyve1 is not restricted to the lymph vessels: expression in normal liver blood sinusoids and down-regulation in human liver cancer and cirrhosis. *Cancer Research* **61**:8079–8084. DOI: 11719431.

Mueller DL. 2010. Mechanisms maintaining peripheral tolerance. *Nature Immunology* **11**:21–27. DOI: https://doi.org/10.1038/ni.1817, PMID: 20016506

Partanen TA, Arola J, Saaristo A, Jussila L, Ora A, Miettinen M, Stacker SA, Achen MG, Alitalo K. 2000. VEGF-C and VEGF-D expression in Neuroendocrine cells and their receptor, VEGFR-3, in Fenestrated blood vessels in human tissues. *FASEB Journal* **14**:2087–2096. DOI: https://doi.org/10.1096/fj.99-1049com, PMID: 11023993

Pawlak JB, Bálint L, Lim L, Ma W, Davis RB, Benyó Z, Soares MJ, Oliver G, Kahn ML, Jakus Z, Caron KM. 2019. Lymphatic Mimicry in maternal endothelial cells promotes Placental spiral artery remodeling. *Journal of Clinical Investigation* **129**:4912–4921. DOI: https://doi.org/10.1172/JCI120446
Pawlak JB, Caron KM. 2020. Lymphatic programming and specialization in hybrid vessels. *Frontiers in Physiology* **11**:114. DOI: https://doi.org/10.3389/fphys.2020.00114, PMID: 32153423

Pérez-Clavier R, Harrison RG, Macmillian EW. 1982. The pattern of the Lymphatic drainage of the rat Epididymis. *Journal of Anatomy* **134**:667-675 PMID: 7138628.

Pierucci-Alves F, Midura-Kiela MT, Fleming SD, Schultz BD, Kiela PR. 2018. Transforming growth factor beta signaling in Dendritic cells is required for Immunotolerance to Sperm in the Epididymis. *Frontiers in Immunology* **9**:1882. DOI: https://doi.org/10.3389/fimmu.2018.01882, PMID: 30166986

Privratsky JR, Newman PJ. 2014. PECAM-1: regulator of endothelial junctional integrity. *Cell and Tissue Research* **355**:607-619. DOI: https://doi.org/10.1007/s00441-013-1779-3

Rantakari P, Auvinen K, Jäppinen N, Kapraali M, Valtonen J, Karikoski M, Gerke H, Iftakhar-E-Khuda I, Keuschnigg J, Umemoto E, Tohya K, Miyasaka M, Elima K, Jalkanen S, Salmi M. 2015. The endothelial protein PLVAP in Lymphatics controls the entry of lymphocytes and antigens into lymph nodes. *Nature Immunology* **16**:386-396. DOI: https://doi.org/10.1038/ni.3101

Rinaldi VD, Donnard E, Gellatly K, Rasmussen M, Kucukural A, Yukselen O, Garber M, Sharma U, Rando OJ. 2020. An Atlas of cell types in the mouse Epididymis and Vas Deferens. *eLife* **9**:e55474. DOI: https://doi.org/10.7554/eLife.55474, PMID: 32729827

Russell PS, Hong J, Windsor JA, Itkin M, Phillips ARJ. 2019. Renal Lymphatics, anatomy, physiology, and clinical implications. *Frontiers in Physiology* **10**:251. DOI: https://doi.org/10.3389/fphys.2019.00251, PMID: 30923503

Shibahara H, Wakimoto Y, Fukui A, Hasegawa A. 2021. Anti-Sperm antibodies and reproductive failures. *American Journal of Reproductive Immunology* **85**:e13337. DOI: https://doi.org/10.1111/aji.13337, PMID: 32885505

Srinivasan RS, Dillard ME, Lagutin OV, Lin FJ, Tsai S, Tsai MJ, Samokhvalov IM, Oliver G. 2007. Lineage tracing demonstrates the venous origin of the mammalian Lymphatic vasculature. *Genes & Development* **21**:2422-2432. DOI: https://doi.org/10.1101/gad.1588407, PMID: 17908929

Stamataki Z, Swadling L. 2020. The liver as an immunological barrier redefined by single-cell analysis. *Immunology* **160**:157-170. DOI: https://doi.org/10.1111/imm.13193, PMID: 32176810

Stammler A, Hau T, Bhushan S, Meinhardt A, Jonigk D, Lippmann T, Pilatz A, Schneider-Hüther I, Middendorff R. 2015. Epididymitis: ascending infection restricted by segmental boundaries. *Human Reproduction* **30**:1557-1565. DOI: https://doi.org/10.1093/humrep/dev112, PMID: 25994666

Stan RV, Tkachenko E, Niesman IR. 2004. Pv1 is a key structural component for the formation of the Stomatal and Fenestral diaphragms. *Molecular Biology of the Cell* **15**:3615-3630. DOI: https://doi.org/10.1091/mbc. e03-08-0593, PMID: 15155804

Stanczuk L, Martinez-Corral I, Ulvmar MH, Zhang Y, Laviña B, Fruttiger M, Adams RH, Saur D, Betsholtz C, Ortega S, Alitalo K, Graupera M, Mäkinen T. 2015. cKit lineage Hemogenic Endothelium-derived cells contribute to mesenteric Lymphatic vessels. *Cell Reports* **10**:1708-1721. DOI: https://doi.org/10.1016/j.celrep. 2015.02.026, PMID: 25772358

Suzuki F. 1982. Microvasculature of the mouse Testis and Excurrent duct system. *The American Journal of Anatomy* **163**:309-325. DOI: https://doi.org/10.1002/aja.1001630404, PMID: 7091016

Svingen T, François M, Wilhelm D, Koopman P. 2012. Three-dimensional imaging of Prox1-EGFP transgenic mouse Gonads reveals divergent modes of Lymphangiogenesis in the Testis and ovary. *PLOS ONE* **7**:e52620. DOI: https://doi.org/10.1371/journal.pone.0052620, PMID: 23285114

Turner TT, Bomgardner D, Jacobs JP, Nguyen QAT. 2003. Association of Segmentation of the Epididymal Interstitium with Segmented tubule function in rats and mice. *Reproduction* **125**:871-878. DOI: https://doi.org/10.1530/rep.0.1250871, PMID: 12773110

Ulvmar MH, Mäkinen T. 2016. Heterogeneity in the Lymphatic vascular system and its origin. *Cardiovascular Research* **111**:310-321. DOI: https://doi.org/10.1093/cvr/cvw175, PMID: 27357637

Voisin A, Whitfield M, Damon-Soubeyrand C, Goubely C, Henry-Berger J, Saez F, Kocer A, Drevet JR, Guiton R. 2018. Comprehensive overview of murine Epididymal mononuclear phagocytes and lymphocytes: unexpected populations arise. *Journal of Reproductive Immunology* **126**:11-17. DOI: https://doi.org/10.1016/j.jri.2018.01.003, PMID: 29421624

Voisin A, Damon-Soubeyrand C, Bravard S, Saez F, Drevet JR, Guiton R. 2020. Differential expression and Localisation of TGF-B Isoforms and receptors in the murine Epididymis. *Scientific Reports* **10**:995. DOI: https://doi.org/10.1038/s41598-020-57839-5, PMID: 31969637

Wigle JT, Oliver G. 1999. Prox1 function is required for the development of the murine Lymphatic system. *Cell* **98**:769-778. DOI: https://doi.org/10.1016/s0092-8674(00)81511-1, PMID: 10499794

Wijayarathna R, Pasalic A, Nicolas N, Biniwale S, Ravinthiran R, Genovese R, Muir JA, Loveland KL, Meinhardt A, Fijak M, Hedger MP. 2020. Region-specific immune responses to autoimmune Epididymitis in the murine reproductive tract. *Cell and Tissue Research* **381**:351-360. DOI: https://doi.org/10.1007/s00441-020-03215-8, PMID: 32383098

Xu Y, Yuan L, Mak J, Pardanaud L, Caunt M, Kasman I, Larrivée B, del Toro R, Suchting S, Medvinsky A, Silva J, Yang J, Thomas J-L, Koch AW, Alitalo K, Eichmann A, Bagri A. 2010. Neuropilin-2 mediates VEGF-C-induced Lymphatic sprouting together with Vegfr3. *Journal of Cell Biology* **188**:115-130. DOI: https://doi.org/10.1083/jcb.200903137

Xu B, Washington AM, Hinton BT. 2016. Initial segment differentiation begins during a critical window and is dependent upon Lumicrine factors and SRC proto-Oncogene (SRC) in the mouse. *Biology of Reproduction* **95**:15. DOI: https://doi.org/10.1095/biolreprod.116.138388, PMID: 27281706