Steroid hormone bioavailability is controlled by the lymphatic system

Production of de‑novo steroid hormones from cholesterol in HLEC and dLEC. HLECs and dLECs were cultured in steroid-free medium and incubated with 3H-cholesterol to quantify the de-novo synthe- sis of progesterone, corticosterone, cortisol, testosterone and estradiol (HLEC), and of progesterone, 11-deoxy- cortisol, corticosterone, aldosterone and cortisol (dLEC), after an incubation time of 24 h. No significant de-novo steroid hormone production could be detected out of the precursor cholesterol by thin layer chromatography (TLC) in HLEC (Supplementary Fig. 1a online) and dLEC (Supplementary Fig. 1b online). Characterisation of progesterone metabolism in HLEC, dLEC and HTR‑8/SVneo by TLC. HLECs, dLECs and HTR-8/SVneo cells were cultured in steroid-free medium and incubated with 14C-progesterone for 4 h/8 h/24 h/48 h (HLEC/dLEC) or 1 h/4 h/8 h/24 h (HTR-8/SVneo). Steroids were extracted from superna- tants and the metabolites separated by TLC. Progesterone was time-dependently and significantly converted to one major steroid hormone, 6α-OH-pregnanolone, in all three cell lines with a production rate of 19.3% ± 3.6 (24 h) and 31.7% ± 6.3 (48 h) in HLEC, 13.6% ± 2.5 (24 h) and 23.9% ± 4.2 (48 h) in dLEC and 21.0% ± 11.5 (8 h) and 42.7% ± 17.1 (24  h) in HTR-8/SV neo (Fig.  1). One-way ANOV A, Dunnett’s multiple comparisons test. Mean ± SD, n = 5. Characterization of 6α‑OH‑pregnanolone in HLEC, dLEC and HTR‑8/SVneo by LC–MS. In order to validate the results obtained by TLC, steroidomic analysis of cell supernatants was performed using LC–MS. Samples were measured in untargeted mode and investigated for peaks that showed strongest signal increase over the time course experiment. Subsequently, the accurate mass of peaks of interest was compared with the theoretical value for 6α-OH-pregnanolone and both, fragmentation spectrum and retention time, were compared to an authentic standard (Figs. 2 and 3). The LC–MS data clearly confirms the production of 6α-OH- pregnanolone from progesterone as obtained by TLC. Gene expression analysis in HLEC, dLEC and HTR‑8/SVneo. HLECs, dLECs and HTR-8/SV neo cells were cultured as described in methods. RNA extraction and Real-time PCR was conducted using assay on demand primers or primers and probes from the Roche library (Tables  1, 2). Cyclophilin A served as endog- enous control. The mRNA expression of enzymes involved in progesterone metabolism (Table 3) was quantified first. HLEC and dLEC express high levels of SRD5A1, SRD5A3, AKR1C1, AKR1C2, AKR1C3, but they do not express SRD5A2, AKR1D1, and AKR1C4. No expression of HSD3B1/2 was found (Ct > 35). HTR-8/SV neo express high levels of SRD5A1, SRD5A3, AKR1C2 and AKR1C3, but rather low levels of SRD5A2 and AKR1C1. AKR1D1 and AKR1C4 as well as HSD3B1/2 are not expressed in HTR-8/SV neo (Ct > 35). Next, the mRNA expression of the progesterone regulated gene (PIBF1) and of steroid hormone receptors in HLEC, dLEC and HTR-8/SV neo was assessed (Table  4). All three cell lines significantly express PIBF1 as well as the glucocorticoid (NR3C1) and the mineralocorticoid receptor (NR3C2), the two membrane bound progesterone receptors type 1 and 2 (PGRMC1 and 2) and the estrogen receptor (GPER1). They show no or very low expression of the nuclear progesterone receptor (NR3C3). Steroidogenic enzymes involved in de-novo steroid hormone production were assessed thereafter (Table 5). HLEC and dLEC express no CYP17A1, CYP21A2, CYP11B1, CYP11B2, HSD11B1 and no HSD11B2, but they highly express HSD17B2. HTR-8/SV neo express no CYP17A1, CYP21A2, CYP11B1, CYP11B2, and no HSD11B1 but they express HSD11B2 and HSD17B2. Assessment of the mRNA expression of steroidogenic proteins involved in de-novo steroid hormone produc- tion (Table 6) revealed no significant expression of StAR in HLEC and dLEC, while FDXR, FDX1 and NR5A2 were strongly expressed. HTR-8/SV neo significantly expressed StAR, FDXR, FDX1 and NR5A2. 3 Vol.:(0123456789)Scientific Reports | (2021) 11:9666 | https://doi.org/10.1038/s41598-021-88508-w www.nature.com/scientificreports/ Derived cascade of progesterone metabolism in HLEC, dLEC and HTR‑8/SV neo. Based on our findings by Real-time PCR, LC–MS and Thin layer chromatography, progesterone is metabolized via 5α-DHP and allopregnanolone to 6α-OH-pregnanolone in HLECs, dLECs and HTR-8/SVneo (Fig.  4). Furthermore, HLEC and dLEC efficiently back-convert 20α-hydroxyprogesterone to progesterone due to high HSD17B2 expression. Inhibition of 6α‑OH‑pregnanolone formation in HLEC, dLEC and HTR‑8/SV neo using the SRD5A1‑3 inhibitor dutasteride. HLECs, dLECs and HTR-8/SV neo cells were cultured and incubated for 24 h with the 5α-reductase inhibitor dutasteride at 2 different concentrations, 10−5  M and 10−6  M for HLEC and HTR-8/SV neo and 10−6   M and 10−8   M for dLEC. Supernatants were analyzed by TLC to measure the conversion of 14C-progesterone to 6α-OH-pregnanolone (Fig.  5). Dutasteride significantly inhibited 6α-OH- pregnanolone formation at both concentrations in all three cell lines. In untreated HLECs (DMSO) 39.2% ± 6.3% 6α-OH-pregnanolone was produced from progesterone (Fig. 5a). This production was reduced to 6.3% ± 5.8%, respectively 6.2% ± 4.7% with dutasteride 10−5 M and 10−6 M. A similar inhibitory effect by dutasteride was observed in dLECs (Fig.  5b) and HTR-8/SV neo cells (Fig.  5c). In dLECs, production of 6α-OH-pregnanolone went from 14.5% ± 0.9% to 4.1% ± 0.6% with dutasteride 10−6 M and to 5.7% ± 0.2% with dutasteride 10−8 M. Similarly, HTR-8/SV neo cells produced 32.1% ± 9.1% 6α-OH-pregnanolone (DMSO), and this was reduced to 3.0% ± 2.7% with dutasteride 10−5 M and to 3.7% ± 2.9% with dutasteride 10−6 M. The involvement of the SRD5A1-3 for the production of 6α-OH-pregnanolone was confirmed by enzymatic activity and RT-PCR. One- way ANOV A, Dunnett’s multiple comparisons test. Mean ± SD, n = 3. Figure 1. Characterization of progesterone metabolism in HLEC, dLEC and HTR8/SV neo by TLC. Metabolism of progesterone in HLEC, dLEC and HTR-8/SV neo. (a) Characteristic phosphorimager pictures of time-dependent conversion of progesterone to a main metabolite, 6α-OH-pregnanolone. Cell-free controls were run for time point 0 h (HLEC, dLEC and HTR-8/SV neo), for time point 24 h (HTR-8/SVneo) and for time point 48 h (HLEC and dLEC). (b) Densitometric quantification of all performed experiments. Progesterone was time-dependently and significantly converted to 6α-OH-pregnanolone in all three cell lines. Data were normalized to condition 0 h = 100% progesterone. One-way ANOV A, Dunnett’s multiple comparisons test, n = 5. Production of 6α-OH-pregnanolone: HLEC: 8 h-24 h ***p < 0.0001, 24 h-48 h **p = 0.0002. dLEC: 8 h-24 h ** p = 0.0026, 24 h-48 h ** p = 0.0020. HTR-8/SV neo: 1 h-24 h ** p = 0.0002. Green rectangle/black column: substrate (progesterone). Red rectangle/white column: product (6α-OH-pregnanolone). * p < 0.05, ** p < 0.01, *** p < 0.0001, ns = not significant. 4 Vol:.(1234567890)Scientific Reports | (2021) 11:9666 | https://doi.org/10.1038/s41598-021-88508-w www.nature.com/scientificreports/ Figure 2. Identification of 6α-OH-pregnanolone by LC–MS/MS in HTR-8/SV neo cells: LC–MS trace of m/z 317.2469 corresponding to a sum formula of C21H34O3 in an authentic standard (a) and in HTR-8/SV neo cells (b). MS/MS fragment spectra recorded at LC–MS peak maximum for an authentic standard (c) and in HTR-8/ SVneo cells (d). Figure 3. Identification of 6α-OH-pregnanolone by LC–MS/MS in HLEC and dLEC cells. LC–MS trace of m/z 317.2469 corresponding to a sum formula of C21H34O3 in an authentic standard (a), in HLEC cells (b) and in dLEC cells (c). MS/MS fragment spectra recorded at LC–MS peak maximum in an authentic standard (d), in HLEC cells (e) and in dLEC cells (f). 5 Vol.:(0123456789)Scientific Reports | (2021) 11:9666 | https://doi.org/10.1038/s41598-021-88508-w www.nature.com/scientificreports/ Detection of HSD17B2 protein and its activity in HLEC, dLEC and HTR‑8/SV neo. Proteomics analysis revealed a strong expression of HSD17B2 protein in HLEC and dLEC. HTR-8/SV neo also expressed HSD17B2, but much weaker than the LECs. HSD17B2 activity was next quantified by LC–MS by measuring the conversion of testosterone to androstenedione, of 17β-estradiol to estrone, of androstenediol to dehydroepian- drosterone (DHEA) and DHEA-S and of 20α-OHP to progesterone. There was a significant HSD17B2 activity, both in HLEC (Fig.  6) and dLEC (data not shown). LECs, completely converted testosterone, 17β-estradiol, androstenediol and 20α-OHP to androstenedione, estrone, DHEA/DHEA-S and progesterone within 24  h.  Table 1. Human primers and probes used in Real-time PCR. Gene Accession number Primers Probe Amplicon size (nt) SRD5A1 NM_001324322.1 5′-ttg gag aaa tca tgg agt ggt-3’ 37 141 5′-act ctt caa att tcc gga ggt a-3’ SRD5A2 NM_000348.3 5′-cag cta cag gat tcc aca agg-3’ 50 72 5′-tca atg atc tca ccg agg aa-3’ SRD5A3 NM_024592.4 5′-ggc ttc atg gtt tgc tca g-3’ 24 96 5′-gca gcc aca gaa ata cta gca c-3’ AKR1D1 NM_001190906.1 5′-cca aaa tga aca cga agt tgg-3’ 7 108 5′-aca tga ttt gta gcc cat agc tt-3’ AKR1C1 NM_001353.5 5′-cat gcc tgt cct ggg att t-3’ 49 109 5′-aga atc aat atg gcg gaa gc-3’ AKR1C2 NM_001354.5 5′-cta tgc gcc tgc aga ggt-3’ 31 114 5′-acc tgc tcc tca tta ttg taa aca t-3’ AKR1C3 NM_003739.5 5′-cat tgg ggt gtc aaa ctt ca-3’ 27 112 5′-ccg gtt gaa ata cgg atg ac-3’ AKR1C4 NM_001818.3 5′-agg tga gac gcc act acc aa-3’ 53 96 5′-tcc tta cac ttc tcc atg acc tc-3’ HSD3B1 NM_000862.2 5′-atc atc cgc ctc ttg gtg-3’ 17 114 5′-cag ctt ggt ctt gtt ctg ga-3’ HSD3B2 NM_000198.3 5′-ctt gga caa ggc ctt cag ac-3’ 50 78 5′-tca agt aca gtc agc ttg gtc ct-3’ HSD17B2 NM_002153.2 5′-gcc aag aat tgt tac ctg tgg-3’ 8 88 5′-tcc aga tac ttg cac aaa gca-3’ HSD11B1 NM_001206741.1 5′-tct gtg ttc ttg gcc tca tag a-3’ 8 75 5′-gag ctg ctt gca tat gga cta tc-3’ FDXR NM_001258015.2 5′-tcc tac tga ccc cac ctg ag-3’ 8 77 5′-tcg ac tct gcc tca gta cac c-3’ NR5A2 NM_205860.2 5′-ccg aca agt ggt aca tgg aa-3’ 61 88 5′-tcc ggc ttg tga tgc tat ta-3’ Table 2. Human assay on demand primers (Applied Biosystems) used in Real-time PCR. Gene # Cat Nr PIBF1 Hs00197131_m1 NR3C1 Hs00353740_m1 NR3C2 Hs01031809_m1 NR3C3 Hs01556702_m1 PGRMC1 Hs009998344_m1 PGRMC2 Hs01128672_m1 GPER1 Hs00173506_m1 CYP17A1 Hs01124136_m1 CYP21A2 Hs00416901_g1 CYP11B1 Hs01596404_m1 CYP11B2 Hs01597732_m1 HSD11B2 Hs00388669_m1 StAR Hs00986559_g1 FDX1 Hs01070067_g1 6 Vol:.(1234567890)Scientific Reports | (2021) 11:9666 | https://doi.org/10.1038/s41598-021-88508-w www.nature.com/scientificreports/ Table 3. Gene expression analysis in HLEC, dLEC and HTR-8/SVneo. mRNA expression of enzymes involved in progesterone metabolism in HLEC, dLEC and HTR-8/SV neo. Expression levels are shown as ct values per 12.5 ng cDNA. Y ellow: significant expression (ct < 35). Blue: no or very low expression (ct ≥ 35). SRD5A1 (5α-reductase isoform 1), SRD5A2 (5α-reductase isoform 2), SRD5A3 (5α-reductase isoform 3), AKR1D1 (5β-reductase), AKR1C1 (20α-HSD), AKR1C2 (3α-HSD3), AKR1C3 (HSD17B5), AKR1C4 (3α-HSD1), HSD3B1, HSD3B2. Appropriate positive and negative controls were run. SRD5A1 SRD5A2 SRD5A3 AKR1D1 AKR1C1 AKR1C2 AKR1C3 AKR1C4 HSD3B1 HSD3B2 HLEC 28 38 26 > 40 31 26 25 > 40 > 40 > 40 dLEC 28 38 26 > 40 31 29 25 > 40 > 40 > 40 HTR-8/ SV neo 26 35 27 > 40 35 33 30 > 40 35 34 Table 4. Gene expression analysis in HLEC, dLEC and HTR-8/SVneo. mRNA expression of PIBF1, and of steroid hormone receptors in HLEC, dLEC and HTR-8/SV neo. Expression levels are shown as ct values per 12.5 ng cDNA. Y ellow: significant expression (ct < 35). Blue: no or very low expression (ct ≥ 35). PIBF1 (progesterone induced blocking factor 1), NR3C1 (glucocorticoid receptor), NR3C2 (mineralocorticoid receptor), NR3C3 (progesterone receptor), PGRMC1 (progesterone receptor membrane component 1), PGRMC2 (progesterone receptor membrane component 2), GPER1 (G protein-coupled estrogen receptor 1). Appropriate positive and negative controls were run. PIBF1 NR3C1 (GR) NR3C2 (MR) NR3C3 (PR) PGRMC1 PGRMC2 GPER1 HLEC 26 26 28 35 27 25 31 dLEC 26 26 30 > 40 27 26 34 HTR-8/ SV neo 25 24 33 36 26 24 32 Table 5. Gene expression analysis in HLEC, dLEC and HTR-8/SVneo. mRNA expression of steroidogenic enzymes involved in de-novo steroid hormone production in HLEC, dLEC and HTR-8/SV neo. Expression levels are shown as ct values per 12.5 ng cDNA. Y ellow: significant expression (ct < 35). Blue: no or very low expression (ct ≥ 35). CYP17A1 (steroid-17α-hydroxylase), CYP21A2 (steroid-21-hydroxylase), CYP11B1 (steroid-11β-hydroxylase), CYP11B2 (aldosterone synthase), HSD11B1 (11β-hydroxysteroid-dehydrogenase 1), HSD11B2 (11β-hydroxysteroid-dehydrogenase 2), HSD17B2 (17β-hydroxysteroid dehydrogenase 2). Appropriate positive and negative controls were run. CYP17A1 CYP21A2 CYP11B1 CYP11B2 HSD11B1 HSD11B2 HSD17B2 HLEC > 40 37 > 40 > 40 36 35 23 dLEC > 40 38 > 40 > 40 > 40 36 23 HTR-8/ SV neo 38 35 > 40 > 40 38 31 32 7 Vol.:(0123456789)Scientific Reports | (2021) 11:9666 | https://doi.org/10.1038/s41598-021-88508-w www.nature.com/scientificreports/ In HLECs, testosterone (0 h: 1000 ± 0; 4 h: 318.7 ± 160.6, 24 h: 10.3 ± 14.9) was time-dependently converted to androstenedione (0 h: 1.4 ± 1.0; 4 h: 776.6 ± 148.1; 24 h: 899.9 ± 114.4) (Fig. 6a); 17β-estradiol (0 h: 1000 ± 0; 4 h: 205.1 ± 159.7; 24 h: 7.7 ± 15.3) to estrone (0 h: 6.5 ± 4.4; 4 h: 659.4 ± 263.9; 24 h: 955.5 ± 62.0) (Fig. 6b); Andros- tenediol (0 h: 1000 ± 0; 4 h: 816.5 ± 217.2; 24 h: 14.3 ± 6.5) to DHEA (0 h: 2.2 ± 2.3; 4 h: 375.0 ± 277.6; 24 h: 946.3 ± 55.8) (Fig. 6c) and 20α-OHP (0 h: 1000 ± 0; 4 h: 618.9 ± 116.6; 24 h: 101.9 ± 31.9) to progesterone (0 h: 39.0 ± 51.2; 4 h: 542.3 ± 166.6; 24 h: 929.6 ± 104.6) (Fig. 6d). One-way ANOV A, Dunnett’s multiple comparisons test. Mean ± SD, n = 3. Progesterone metabolism in mouse lymph nodes. Mouse lymph nodes, isolated from male mice (n = 7), were incubated for 24  h and 48  h with 14C-progesterone. Progesterone metabolism was analysed in supernatants as described. Adrenal gland was used as positive control. Progesterone was converted into 7 differ- ent metabolites in all lymph nodes. The progesterone metabolism pattern was different as compared to the one found in the adrenal gland, confirming progesterone metabolism and not de-novo steroid hormone production (Fig. 7). Table 6. Gene expression analysis in HLEC, dLEC and HTR-8/SVneo. mRNA expression of steroidogenic proteins involved in de-novo steroid hormone production in HLEC, dLEC and HTR-8/SV neo. Expression levels are shown as ct values per 12.5 ng cDNA. Y ellow: significant expression (ct < 35). Blue: no or very low expression (ct ≥ 35). CYP17A1 (steroid-17α-hydroxylase), CYP21A2 (steroid-21-hydroxylase), CYP11B1 (steroid-11β- hydroxylase), CYP11B2 (aldosterone synthase), HSD11B1 (11β-hydroxysteroid-dehydrogenase 1), HSD11B2 (11β-hydroxysteroid-dehydrogenase 2), HSD17B2 (17β-hydroxysteroid dehydrogenase 2). Appropriate positive and negative controls were run. StAR FDXR FDX1 NR5A2 (LRH-1) HLEC 36 27 23 25 dLEC 37 28 23 25 HTR-8/ SV neo 33 26 24 30 Figure 4. Putative progesterone metabolism pathway in HLECs, dLECs and HTR-8/SV neo. Based on our mRNA, proteomics, TLC and LC–MS data, we propose the following pathway in green to take place. The pathway with the bulky arrows down to 6α-OH-pregnanolone is favored over the pathway to the 20α-hydroxyprogesterone in HLECs, dLECs and HTR-8/SV neo. HLEC and dLEC highly express HSD17B2 and efficiently back-convert 20α-hydroxyprogesterone to progesterone. 8 Vol:.(1234567890)Scientific Reports | (2021) 11:9666 | https://doi.org/10.1038/s41598-021-88508-w www.nature.com/scientificreports/ All assessed steroid hormones, enzymes and steroidogenic factors in this manuscript are summarized in Supplementary Fig. 2 online in a schematic view. Immune response of CD4+ and CD8+ T cells upon stimulation with progesterone metabo‑ lites. PBMCs of three different healthy female donors were isolated and pre-incubated with the steroid hor- mones progesterone, 5α-DHP , 6α-OH-pregnanolone, and with dexamethasone (Dexa) as positive control for 6 h before activation. Solvent control of all steroid hormones was EtOH. Activation of T cells was thereafter induced by PMA/Ionomycin. TNF-α and IFN-γ production was assessed 24 h later using FACS analysis. Progesterone and Dexa, strongly decreased TNF-α positive CD4+ (Fig. 8a) and CD8+ (Fig. 8b) T cell populations. 5α-DHP and 6α-OH-pregnanolone reduced the percentage of TNF-α positive CD4+ T cell subpopulations only in donor 3, while both progesterone metabolites reduced the percentage of TNF-α positive CD8+ T cells in all three donors (Fig. 8). Dexa strongly decreased IFN-γ positive CD4+ (Fig. 9a) and CD8+ (Fig. 9b) T cell populations. Progesterone reduced % of IFN-γ positive CD8+ T cells in all three donors, while it inhibited IFN-γ positive CD4+ T cells only in donor 1. The progesterone metabolites 5α-DHP and 6α-OH-pregnanolone showed a strong effect in CD8+ T cells of donor 3, where both steroids decreased IFN-γ positive CD8+ T cells by 20% (Fig. 9). Figure 5. Inhibitory effect of dutasteride on 5α-reductase (SRD5A1-3) in HLECs, dLECs and HTR-8/SV neo. Characteristic phorphorimager pictures of (a) HLEC, (b) dLEC, (c) HTR-8/SV neo and their densitometry (d) of three independent experiments. Cells were cultured without (lane 1) and with (lane 2 and 3) the indicated concentration of dutasteride and 14C-progesterone. Lane 4 is the cell-free control, where 14C-progesterone was incubated for 24 h without cells. HLEC and HTR-8/SV neo: (1) no dutasteride = DMSO control, (2) dutasteride 10−5 M, (3) dutasteride 10−6 M, (4) no cells; dLEC: (1) no dutasteride, (2) dutasteride 10−6 M, (3) dutasteride 10−8 M, (4) no cells. Progesterone was significantly metabolized to 6α-OH-pregnanolone in the DMSO control, while dutasteride significantly inhibited the conversion of progesterone to 6α-OH-pregnanolone in all three cell lines. Data were normalized to condition 0 h = 100% progesterone. One-way ANOV A, Dunnett’s multiple comparisons test, n = 3. HLEC: DMSO control *** p < 0.0001; dutasteride 10−5 M *** p < 0.0001; dutasteride 10−6 M *** p < 0.0001. dLEC: DMSO control *** p < 0.0001; dutasteride 10−6 M ** p = 0.0001; dutasteride 10−8 M ** p = 0.0006. HTR-8/SV neo: DMSO control *** p < 0.0001; dutasteride 10−5 M ** p = 0.0005; dutasteride 10−6 M ** p = 0.0006. Green rectangle/black column: substrate (progesterone). Red rectangle/white column: product (6α-OH-pregnanolone). * p < 0.05, ** p < 0.01, *** p < 0.0001, ns = not significant. 9 Vol.:(0123456789)Scientific Reports | (2021) 11:9666 | https://doi.org/10.1038/s41598-021-88508-w www.nature.com/scientificreports/

To our knowledge, this is the first report showing detailed steroid hormone metabolism in the lymphatic sys- tem. The lymphatic endothelium could be identified as an important regulator of steroid hormone bioavail- ability, yet also as a target of steroid hormones by expressing the respective steroid hormone receptors. De-novo steroid hormone synthesis from cholesterol is absent in lymphatic tissue, which is consistent with the lack of steroidogenic acute regulatory protein (StAR) in LECs. However, FDXR, FDX1 and LRH-1, genes involved in steroidogenesis40–42, are expressed in LECs. Progesterone, a steroid hormone critical in regulating immune responses and pregnancy, is intensely metabolized. In primary human LECs derived either from lymph nodes or from dermal lymphatic vessels, progesterone is mainly converted to a single metabolite, which was identified as 6α-OH-pregnanolone (5α-pregnan-3α, 6α-diol-20-one). Enzymatic steps involved are the two isoforms of the 5α-reductase, SRD5A1 and SRD5A3, the isoform 1 of 3α-HSD (AKR1C1) and a 6α-hydroxylase. Besides converting progesterone to the down-stream metabolites 5α-dihydroprogesterone (5α-DHP), allopregnanolone and 6α-OH-pregnanolone, HLECs and dLECs were able to reactivate progesterone from the less potent metabo- lite 20α-OHP . Further steroidogenic activity of HLEC and dLEC included deactivation of testosterone and 17β-estradiol to androstenedione and estrone as well as androstenediol to DHEA and DHEA-S. These findings are pivotal in understanding the progesterone-dependent immunomodulation observed earlier by various groups7,10,11,43. The data show that progesterone is not only metabolized to the primary and active metabolite 5α-DHP which plays a critical role in pregnancy44, but is also restored from 20α-OHP . Figure 6. Activity of HSD17B2 in HLECs. Time-dependent conversion of testosterone to androstenedione, 17β-estradiol to estrone, androstenediol to DHEA-S and DHEA as well as 20α-hydroxyprogesterone to progesterone was assessed by LC–MS. Data were normalized to condition at 0 h = 1000 nM = 10−6 M for each compound. Y-axis shows steroid hormone concentrations in nM. Testosterone was significantly converted to androstenedione (testosterone 0 h/4 h/24 h: *** p < 0.0001/** p = 0.0019; androstenedione 0 h/4 h/24 h: *** p < 0.0001/ns). 17β-estradiol was significantly converted to estrone (17β-estradiol 0 h/4 h/24 h: *** p < 0.0001/ ns; estrone 0 h/4 h/24 h: *** p < 0.0001/* p = 0.01). Androstenediol was significantly converted to DHEA (androstenediol 0 h/4 h/24 h: ns/*** p < 0.0001; DHEA 0 h/4 h/24 h: ** p = 0.0087/** p = 0.0002). 20α-OHP was significantly converted to progesterone (20α-OHP 0 h/4 h/24 h: *** p < 0.0001/*** p < 0.0001; progesterone 0 h/4 h/24 h: *** p < 0.0001/*** p < 0.0001). One-way ANOV A, Dunnett’s multiple comparisons test, n = 3. Black rectangle/column: substrate. White rectangle/column: product. * p < 0.05, ** p < 0.01, *** p < 0.0001, ns = not significant. 10 Vol:.(1234567890)Scientific Reports | (2021) 11:9666 | https://doi.org/10.1038/s41598-021-88508-w www.nature.com/scientificreports/ Progesterone metabolism to 6α-OH-pregnanolone and the reactivation of progesterone from 20α-OHP by LECs serves to adapt immune tolerance and to control local and/or systemic progesterone bioavailability in high progesterone conditions as found during pregnancy. As the exact mechanism for this phenomenon is unknown and given the huge surface of the lymphatic endothelium throughout the body, it is supposable that local pro - gesterone metabolism in LECs could be the key regulator in preventing cancer cells from being destroyed. The direct immunosuppressive role of progesterone in both reproduction and tumor progression was impres- sively shown by Szekeres-Bartho and Polgar45. Immunosuppression by progesterone is mediated via the proges- terone regulated gene PIBF18,11,13,46, which has now also been found to be highly expressed in LECs. A reduced expression of PIBF in threatened pregnancies47–51 and an increased production of PIBF in cancer52,53 are clearly linked to a failed immune modulation. Production of TNF-α and IFN-γ in activated CD4+ and CD8+ T cells from 3 different donors was modulated upon stimulation with progesterone and its metabolites. It needs to be investigated if the interaction of LECs and the immune cells is PIBF1 regulated. Steroid hormone effector mechanisms such as the glucocorticoid receptor (GR), the mineralocorticoid recep- tor (MR), the G protein-coupled estrogen receptor 1 (GPER1), and the progesterone receptor membrane com- ponents PGRMC1 and PGRMC2 are present in HLECs and dLECs, while no expression of nuclear progesterone receptors (PRA and PRB) was found. Most of the known anti-inflammatory effects of progesterone are transmitted through the GR, with proges- terone binding to the GR even though progesterone can also bind to PGRMCs54,55. PGRMCs were overexpressed in the maternal–fetal interface and in the embryonic/fetal trophectoderm in pregnancy56,57 as well as in T cells during pregnancy58,59. Since many steroid hormone receptors are expressed by immune cells12,60–62 it is most likely that the LECs not only have an autocrine, but also a paracrine role in controlling immune responses through the action of steroid hormones via binding to a receptor. The serum concentrations of progesterone are according to the literature much too low to support the concept of a generalized immunosuppression63,64 and are only sufficient to inhibit peripheral natural killer cell activity in normal pregnant women. The presence of lymphatic vessels throughout the body, together with the immu - nosuppressive role of the placenta, might be a much better explanation for the systemic immunosuppression. The mononuclear phagocyte system stimulates VEGF-C release and regulates thereby lymphangiogenesis65. High progesterone levels seem to coincide with high VEGF-C expression66. Activated macrophages and mono- cytes releasing VEGF-C thus enhance LEC growth and could well contribute to immune regulatory mechanisms via progesterone. Figure 7. Progesterone metabolism in male mice lymph nodes and adrenal gland. Phosphorimager picture of TLC showing time-dependent conversion of progesterone to down-stream metabolites in lymph nodes. For positive control, adrenal glands are used under similar experimental conditions. Tissue-free controls were run for both time points. Line 1 at timepoint 24 h and line 2 at 48 h were entire lymph nodes, while all other lymph nodes were cut in half. The figure shows results of 7 isolated lymph nodes. Each lymph node was from a separate mouse. * p < 0.05, ** p < 0.01, *** p < 0.0001, ns = not significant. 11 Vol.:(0123456789)Scientific Reports | (2021) 11:9666 | https://doi.org/10.1038/s41598-021-88508-w www.nature.com/scientificreports/ In 2009, Zhang et  al. found, that in isolated mature adipocytes progesterone was converted to 20α-hydroxyprogesterone as the main metabolite, most likely through the activity of aldo–keto reductases AKR1C1, AKR1C2 and AKR1C3 (20α-HSD, 3α-HSD type 3 and HSD17B5), respectively67. Even though HLECs and dLECs express the same enzymes, 6α-OH-pregnanolone was found instead as major metabolite of progester- one, and not 20α-hydroxyprogesterone. That LECs favour the 5α-reductase/3α-HSD pathway over the 20α-HSD pathway might be explained by the different tissues investigated and eventually the limitations of some cofactors. The strength of this study is the thorough analysis of different steroid hormone metabolites using state-of-the- art technology. It shows enzyme expression on mRNA and on protein level, as well as functional assays including inhibitory experiments. All results are referred to our self-designed and validated model of a positive control, the placental cell line HTR-8/SV neo. Since fast, significant and similar progesterone metabolism was found in HTR-8/SV neo cells, they were used as a model of reference for all experiments. Taking a placental cell line as model is obvious, as the placenta is the primary organ exposed to high progesterone levels. A direct immunologic effect of the lymphatic progesterone metabolites on CD4+ and CD8+ T cells could be shown as well. Investigation of further immunomodulatory effects of the lymphatic progesterone metabolism and proges - terone reactivation described herein on immune and cancer cells is the goal of future studies. Verification of our observations in an in vivo model is difficult due to the systemic presence of steroid hormones, further complicated by the delicate structure of lymphatic vessels, therefore precluding differential concentration measurements. But, the capacity of steroid hormone metabolism and reactivation in lymph nodes isolated from both sexes, including pregnancy, should be explored. Additionally, future detailed assessment in patients treated with the 5α-reductase inhibitors dutasteride or finasteride could be investigated with more accuracy. In conclusion, we identified and characterized the lymphatic system as major steroid hormone metabolizing tissue, which is additionally able to reactivate progesterone from 20α-OHP and therefore adapts progesterone availability. This is important in high progesterone conditions such as in pregnancy in order to adjust local immune tolerance. The impact of these findings on cancer and pregnancy research might be meaningful and could provide new targets for the treatment of cancer and pregnancy related diseases. Furthermore, these results could also help to elucidate disease mechanisms in autoimmunity and allergy. Figure 8. TNF-α production in CD4+ and CD8+ T cells upon stimulation with progesterone metabolites. Progesterone (Prog), 5α-dihydroprogesterone (5α-DHP), 6α-hydroxypregnanolone (6α-OH-Pregn) and dexamethasone (Dexa), were added at a concentration of 10−3  M to the PBMCs obtained from 3 different female donors for 24 h. 6 h after first steroid hormone contact, PBMCs were activated with PMA/Ionomycin. TNF-α positive CD4+ and CD8+ T cells were counted by FACS analysis. The y-axis shows the % of CD4+ (a) and CD8+ (b) T cells positively staining for TNF-α. * p < 0.05, ** p < 0.01, *** p < 0.0001, ns = not significant. 12 Vol:.(1234567890)Scientific Reports | (2021) 11:9666 | https://doi.org/10.1038/s41598-021-88508-w www.nature.com/scientificreports/

Material and cell lines. Collagen I coated petridishes were from Corning (Milian, Nesselnbach, Swit- zerland). The primary HLEC (Cat. No. 2500, Lots. 19,399, 19,394 and 19,415) were from Sciencell (Chemie Brunschwig, Basel, Switzerland) and the primary human dLEC (CC-2810, Lots. 4F3029 and 4F3037) were from Lonza (Ruwag, Bettlach, Switzerland) or from PromoCell Heidelberg, Germany (C-12217, Lot 431Z021.3). Pri- mary cells were obtained from female donors. The HTR-8/SVneo cells were a gift from Charles H. Graham (Queen’s University, Kingston, Ontario, Canada). The vascular cell basal medium (PCS-100-030) and the sup- plementary factors (PCS-100-041) were from ATCC (LGC Standards GmbH, Wesel, Germany). The Endothelial Cell Growth Medium EGM-2 (CC-3162), including EBM-2 Basal Medium (CC-3156) and EGM-2 Supplements (CC-4176), was from Lonza (Ruwag, Bettlach, Switzerland). The PromoCell MV-2 medium with supplements (C-22121) was from PromoCell (Heidelberg, Germany). Medium RPMI 1640 (Cat.Nr. 21,875) was from Gibco (ThermoFisher Scientific, Reinach, Switzerland). Penicillin/streptomycin (P/S, Cat. No.  0503) and FBS (Cat. No. 0025) were from Sciencell (Chemie Brunschwig, Basel, Switzerland). 3H-cholesterol (ART-0255), 3H-pro- gesterone (ART-0795) and 14C-progesterone (ARC-1398) were from American radiolabeled chemicals (ARC) Inc. (Anawa Trading SA, Wangen ZH, Switzerland). Scintillation solution Optiphase Supermix (Cat. No. 1200- 439) was from Perkin Elmer Schweiz AG (Schwerzenbach, Schweiz). The SRD5A-inhibitor dutasteride (SML 1221) was purchased from Sigma-Aldrich (Buchs, Switzerland). Cholesterol, progesterone, and all steroid hormone standards for LC–MS and TLC analysis were purchased from Steraloids, Inc. (Newport, RI, USA). Pre-coated TLC plates SIL G25 UV254 (REF 809,023) and pre-coated TLC sheets Polygram SIL G/UV254 (REF 805,023) were from Macherey–Nagel GmbH, Düren, Germany. Dichlo- romethane (1.06050), methanol (1.06009) and choloroform (1.02445) were from Merck Millipore, Zug, Swit- zerland. Microbeta 2 Microplate counter 2450 was from Perkin Elmer Schweiz AG (Schwerzenbach, Schweiz). Phosphorimager Typhoon FLA 7000 and the Storage Phosphor Screen (BAS-IP TR 2040 E Tritium Screen) were from GE Healthcare Life Sciences, Glattbrugg, Switzerland. Prime Script RT Reagent Kit Cat. RR037A was from Takara Bio Europe (Saint-Germain-en-Laye, France). Figure 9. IFN-γ production in CD4+ and CD8+ T cells upon stimulation with progesterone metabolites. Progesterone (Prog), 5α-dihydroprogesterone (5α-DHP), 6α-hydroxypregnanolone (6α-OH-Pregn) and dexamethasone (Dexa), were added at a concentration of 10−3  M to the PBMCs of 3 different female donors for 24 h. 6 h after first steroid hormone contact, PBMCs were activated with PMA/Ionomycin. IFN-γ positive CD4+ and CD8+ T cells were counted by FACS analysis. The y-axis shows the % of CD4+ (a) and CD8+ (b) T cells positively staining for IFN-γ. * p < 0.05, ** p < 0.01, *** p < 0.0001, ns = not significant. 13 Vol.:(0123456789)Scientific Reports | (2021) 11:9666 | https://doi.org/10.1038/s41598-021-88508-w www.nature.com/scientificreports/ For LC–MS analysis, a Vanquish UHPLC (equipped with an ACQUITY UPLC HSS T3 Column, 100 Å, 1.8 µm, 1 mm X 100 mm column; Waters, Switzerland) was coupled to a Q Exactive Plus Orbitrap (both Thermo Fisher Scientific, Reinach, Switzerland). Separation was achieved using gradient elution over 11 min using water and methanol both supplemented with 0.1% formic acid (all Sigma-Aldrich, Buchs, Switzerland) as mobile phases. Data analysis was performed using TraceFinder 4.1 (Thermo Fisher Scientific, Reinach, Switzerland). All steroids analyzed by LC–MS and their systematic names are listed in Table 7. Primer pairs for SRD5A1, SRD5A2, SRD5A3, AKR1D1, AKR1C1, AKR1C2, AKR1C3, AKR1C4, HSD3B1, HSD3B2, HSD17B2, HSD11B1, FDXR, and NR5A2 (Table  1) were designed by means of the Universal Pro- beLibrary Assay Design Center software from Roche. All primers were intron-spanning and were synthesized by Microsynth AG (Balgach, Switzerland). Universal ProbeLibrary hydrolysis probes were purchased from Sigma-Aldrich (Buchs, Switzerland). Assay on demand primers for Cyclophilin A (4326316E), PIBF1, NR3C1, NR3C2, NR3C3, PGRMC1, PGRMC2, GPER1, CYP17A1, CYP21A2, CYP11B1, CYP11B2, HSD11B2, StAR, FDX1 (Table 2), and the 7500 Fast Real-time PCR system machine were from Applied Biosystems (Thermo Fisher Scientific, Reinach, Switzerland). The MS part was done by the Proteomic Mass Spectrometry Core Facility PMSCF at DBMR in Bern, Switzerland. For FACS analysis, GolgiPlug, anti-CD3 V500 (clone SK7), anti-CD4 PerCP (clone SK3) and anti-CD8 PE-Cy7 (clone SK1) were purchased from BD Biosciences. The fixable Live/Dead dye was obtained from Invit- rogen. Anti-TNF-α AF647 (clone Mab11) and anti-IFN-γ BV421 (clone 4sB3) and their matched isotype controls were bought from Biolegend. PMA/Ionomycin was obtained from Sigma Aldrich. All methods used in this manuscript were carried out in accordance with relevant guidelines and regulations. Cell culture. Primary HLECs and dLECs were cultured in collagen I coated cell ware. Medium for HLECs was the vascular cell basal medium from ATCC with all supplementary factors but without cortisol (hydro- cortisone). FBS concentration was 5%. Medium for dLECs was the EBM-2 Basal Medium from Lonza with all EGM-2 supplements but without cortisol. FBS concentration was 2%. Primary HLECs and dLECs were used up to passage 6. HTR-8/SV neo cells were cultured in RPMI 1640 5% FBS. In order to minimize steroid hormone contamination in our medium or FBS, all experiments were performed in cortisol-free medium containing 5% or 2% of charcoal treated (ct) FBS. Table 7. List of all steroids analyzed or mentioned in this study with their short and systematic name. Class Short name Systematic name Progesterones Pregnenolone 5-Pregnen-3β-ol-20-one Progesterone 4-Pregnen-3, 20-dione 17OH-Progesterone 4-Pregnen-17-ol-3, 20-dione 20α-Hydroxyprogesterone 4-Pregnen-20α-ol-3-one 5α-Dihydroprogesterone 5α-Pregnan-3, 20-dione 5β-Dihydroprogesterone 5β-Pregnan-3, 20-dione Pregnanolone 5β-Pregnan-3α-ol-20-one Epipregnanolone 5β-Pregnan-3β-ol-20-one Allopregnanolone 5α-Pregnan-3α-ol-20-one Isopregnanolone 5α-Pregnan-3β-ol-20-one 6α-OH-Pregnanolone 5α-Pregnan-3α, 6α-diol-20-one 6α-OH-Epipregnanolone 5α-Pregnan-3β,6α-diol-20-one Allopregnane-3,20-diol 5α-pregnan-3β,20α-diol Mineralocorticoids 11-deoxycorticosterone 4-Pregnen-21-ol-3, 20-dione Corticosterone 4-Pregnen-11β, 21-diol-3, 20-dione Aldosterone 4-Pregnen-11β, 21-diol-3, 18, 20-trione Estrogens 17β-Estradiol 1,3,5(10)-estratriene-3,17β-diol Estriol 1,3,5(10)-estratriene-3,16α,17β-triol Glucocorticoids 11-deoxycortisol 4-Pregnen-17α,21-diol-3,20-dione Cortisone 17,21-dihydroxy-4-pregnene-3,11,20-trione Cortisol 11β,17,21-trihydroxy-4-pregnene-3,20-dione Androgens Androsterone 3α-hydroxy-5α-androstan-17-one Etiocholanolone 3α-Hydroxy-5β-androstan-17-one Dehydroepiandrosterone 5-androsten-3β-ol-17-one Dehydroepiandrosterone-sulphate 5-androsten-3β-ol-17-one sulphate 5α-Dihydrotestosterone 17β-hydroxy-5α-androstan-3-one Testosterone 17β-hydroxy-4-androsten-3-one Androstenedione 4-androsten-3, 17-dione 14 Vol:.(1234567890)Scientific Reports | (2021) 11:9666 | https://doi.org/10.1038/s41598-021-88508-w www.nature.com/scientificreports/ PBMCs were cultured in RPMI1640 containing 10% FBS after isolation. During experiments, PBMCs were maintained in steroid-free RPMI1640. Thin layer chromatography. 14C-progesterone and 3H-cholesterol as substrate for HLEC and dLEC. HLECs and dLECs were cultured in collagen I coated 6-well plates with medium containing 5% or 2% ct FBS for 24 h. After 24 h cells were washed twice with DPBS and 0.05 μCi 14C-progesterone or 1 μCi 3H-cholesterol were added in steroid free medium. Cell culture supernatant was collected into glass vials at the indicated time points. 5 ml ethyl acetate was added and glass vials were vortexed for exactly 1 min. Thereafter, supernatants were stored at − 20 °C until phase separation. The upper phase was transferred into new glass vials and evaporated under a nitrogen stream at 56 °C. The pellet was suspended in 30ul of cold progesterone dissolved in EtOH (10 mg/ml) for the 14C-experiments and in 30ul of the cold mixture cholesterol, progesterone, corticosterone, cortisol, tes- tosterone and estradiol (all 10 mg/ml EtOH) for the 3H-cholesterol experiments in HLEC. For the 3H-cholesterol experiments in dLEC, the pellet was dissolved in 30ul of the cold mixture of cholesterol, progesterone, 11-de- oxycortisol, corticosterone, aldosterone and cortisol (all 10 mg/ml EtOH). A 10ul aliquot of each sample was loaded on TLC glass plates or sheets. Running medium was dichloromethane, methanol, H2O (150:10:1). For Microbeta 2 measurements, spots were visualized and marked under the UV lamp. Marked spots were scratched into scintillation vials. 3.5 ml of scintillation fluid was added and vials were counted in a Microbeta 2. Individual steroid hormones were run on the same plate to allow for localization of the steroid hormones. A control without cells (0 h/48 h for HLEC and dLEC; 0 h/24 h for HTR-8/SV neo) was loaded as baseline/background. Counts are displayed in counts per minute (CPM = CCPM1). Conversion of 14C-progesterone to 6α-OH-pregnanolone or of 3H-cholesterol to progesterone, corticosterone, cortisol, testosterone and estradiol (HLEC), or of 3H-cholesterol to progesterone, 11-deoxycortisol, corticosterone, aldosterone and cortisol (dLEC), was calculated. For visuali- zation by a phorphorimager, TLC sheets were exposed to a Storage Phosphor Screen. Spots were visualized by the Phosphorimager Typhoon FLA 7000. Individual steroid hormones were run on the same sheet to allow for localization of the steroid hormones. A control without cells (0 h/48 h for HLEC and dLEC; 0 h/24 h for HTR- 8/SV neo) was loaded as baseline/background. Quantification was done using ImageJ software. All conversion rates were compared to the controls without cells. Progesterone substrate availability at timepoint 0 h (control without cells) was taken as 100% for all experiments. SRD5A-inhibition with dutasteride. HLECs and dLECs were cultured as described above and incubated with dutasteride dissolved in DMSO and used at a final concentration of 10−5  M and 10−6  M for HLECs and HTR-8/ SV neo and 10−6   M and 10−8   M for dLECs. Total incubation time was 24  h. Conversion of progesterone to 6a-OH-pregnanolone was quantified by TLC. 14C-progesterone as substrate for HTR-8/SV neo. HTR-8/SV neo were cultured in RPMI 1640 medium contain- ing 5% ct FBS for 24 h. After 24 h cells were washed twice with DPBS and 0.05 μCi 14C-progesterone was added in steroid free medium. Cell culture supernatant was collected into glass vials at timepoints 1 h, 4 h, 8 h and 24 h. 5 ml ethyl acetate was added and glass vials were vortexed for exactly 1 h. The rest of the procedure was done exactly as described above for the HLEC and dLEC. Real‑time PCR. mRNA expression of the progesterone metabolizing enzymes (SRD5A1, SRD5A2, SRD5A3, AKR1D1, AKR1C1, AKR1C2, AKR1C3, AKR1C4, HSD3B1, HSD3B2), of PIBF1, of the steroid hormone recep- tors (NR3C1, NR3C2, NR3C3, PGRMC1, PGRMC2, GPER1) of the steroidogenic enzymes (CYP17A1, CYP21A2, CYP11B1, CYP11B2, HSD11B1, HSD11B2, HSD17B2) and of the steroidogenic proteins StAR, FDXR, FDX1, NR5A2 in HLEC, dLEC and HTR-8/SV neo. HLECs and dLECs were cultured in collagen I coated 6-well plates with medium containing 5% or 2% FBS for 24 h. HTR-8/SV neo were cultured in RPMI 1640 medium containing 5% FBS for 24 h. After the incubation time total extraction of RNA was performed using the Trizol method. RNA was reverse transcribed using the Prime Script RT Reagent kit from Takara. Real-time PCR was performed using assay on demand primers or primers and probes designed with the Roche Library in order to detect the following genes: SRD5A1, SRD5A2, SRD5A3, AKR1D1, AKR1C1, AKR1C2, AKR1C3, AKR1C4, HSD3B1, HSD3B2, PIBF1, NR3C1, NR3C2, NR3C3, PGRMC1, PGRMC2, GPER1, CYP17A1, CYP21A2, CYP11B1, CYP11B2, HSD11B1, HSD11B2, HSD17B2, StAR, FDXR, FDX1, and NR5A2 (Tables  1, 2). Cyclophilin A served as endogenous control. Assays were performed in duplicates. Proteomics. HLEC and dLEC were cultured in their medium containing 5% or 2% FBS. Upon confluency, cells were washed 2 × with DPBS, detached with trypsine/EDTA (1 × concentrated) and centrifuged at 800 rpm. Cell pellet was washed 2 × with DPBS and the dried pellet was lysed in 8 M UREA buffer containing a protease inhibitor cocktail. Protein amount was measured by the Qubit Protein Assay. The MS part was done by the Pro- teomic Mass Spectrometry Core Facility PMSCF using standard procedure. LC–MS. HLECs and dLECs were cultured in collagen I coated 6-well plates with medium containing 5% ct FBS for 24  h. After 24  h cells were washed twice with DPBS and the cold testosterone, 17β-estradiol, 5α-dihydroprogesterone and androstenediol were added separately in steroid-free medium at a concentration of 10−6  M. A time-course experiment was performed with incubation times 0 h as baseline, 4 h and 24 h for HLEC and 0 h and 24 h for dLEC. Steroids were extracted from the supernatant by solid phase extraction using Oasis HLB SPE plates (Waters, Switzerland). Remaining substrate concentrations and metabolites were assessed using 15 Vol.:(0123456789)Scientific Reports | (2021) 11:9666 | https://doi.org/10.1038/s41598-021-88508-w www.nature.com/scientificreports/ LC–MS. 20 uL of the extract were injected and separation was achieved using gradient elution over 11 min using water and methanol both supplemented with 0.1% formic acid as mobile phases. 6α-OH-pregnanolone was identified based on accurate mass measurement, comparison of MS/MS spectra and retention times between the analyte in the samples and an authentic standard. Data was normalized to condition at 0 h = 1000 nM = 10−6  M for each compound. Mouse lymph nodes and adrenals. Animal experimentation was approved by the Ethics Committee for Animal Experiments of the Veterinary Administration of the Canton of Berne, Switzerland and conformed to the rules of the Swiss Federal Act on Animal Protection (BE58/19). C57/Bl6 mice were bred according to the rules in the central animal facility of the university of Bern. They were maintained under 12-h dark–light cycles with unrestricted access to food and water. Males used in experiment were fed for 6–8 weeks with regular chow diet. On the day of sacrifice, lymph nodes and adrenal glands were isolated and put on ice until arrival at the lab. Thereafter, organs were washed in PBS and put into warm, steroid-free medium (PCS-100-030) containing 5% ct FBS. The fat was dissected, and the adrenal gland, and half of the lymph nodes were cut in half. The rest of the lymph nodes were uncut and contained a capsule. Every lymph node (cut or entirely), and the adrenal gland were placed in 48-well plates and incubated with 500ul of steroid-free medium containing 5% ct FBS and 0.05 µCi 14C-progesterone for 24 h or 48 h at 37 °C. For controls, medium from wells incubated under the same conditions and incubation time but without tissue was used. After 24 h and 48 h supernatants were collected into glass vials. 5 ml ethyl acetate was added and vials were vortexed exactly for 1 min. The next steps were performed as described in the section thin layer chromatography. Products of progesterone metabolism were visualized with a Phosphorimager. FACS. Peripheral blood mononuclear cells (PBMCs) were isolated from three healthy, non-pregnant female donors using Ficoll gradient centrifugation. After isolation, PBMCs were kept for 30 min in RPMI1640 medium containing 10  mM HEPES, 10% heat-inactivated FBS and 100  μg/ml Streptomycin, 100U/ml penicillin (all from Gibco, Life Technologies). After 30  min, PBMCs were stimulated with progesterone, 5α-DHP , 6α-OH- pregnanolone and dexamethasone for 6 h (priming). Control medium was RPMI1640 containing EtOH, the solvent of the steroid hormones. Final concentration of all steroid hormones on PBMCs was 10−3  M. After the 6 h pre-incubation period with steroid hormones, lymphocytes were activated with 15 ng/ml PMA and 1 μg/ml Ionomycin. Simultaneously, GolgiPlug was added in order to inhibit cytokine release. After a total incubation period of 24 h, PBMCs were stained for CD3, CD4, CD8, IFN-γ and TNF-α and analyzed by FACS as follows: Cells were stained with fixable Live/Dead dye for 30 min one ice. Subsequently cells were washed using PBS and fixed by adding Medium A from the Fix & Perm kit (Invitrogen) for 15 min at RT. Afterwards cells were washed and the permeabilization Medium B and the antibody cocktail containing anti-CD3, anti-CD4, anti-CD8, anti- IFN-γ and anti-TNF-α was added. Cells were incubated for 20 min, washed and resuspended in PBS contain- ing 2% heat-inactivated FBS, 0.05% sodium azide. Data was acquired using a FACSCanto II (BD Biosciences). Matched isotype controls were used for anti-IFN-γ and anti-TNF-α Percentage of IFN-γ and TNFα positive cells was determined in alive cells, which were CD3 positive. The % of CD4+ and CD8+ T cells owning the intracellular cytokines IFN-γ and TNF-α is shown in Fig.  8 (TNF-α) and Fig. 9 (IFN-γ). Statistics. All experiments were performed at least 3 times. Pictures in figures are unprocessed. Representa- tive blots are shown. Densitometry was performed for all phosphorimager pictures using Image J software. Conversion and production rates are displayed as mean ± SD using one-way ANOV A with Dunnetts multiple comparison test. Significance was assigned at p < 0.05. TLC experiments were additionally controlled by CCPM1 measurements in a microbeta2 instrument. For TaqMan results ct values were calculated as mean of all experi- ments. All statistical analyses were performed using GraphPad PRISM version 8 (PRISM, USA). Equipment and settings. For the phosphorimager pictures in Figs. 1, 5 and 7 the image acquisition tool Image J was used. No manipulations have been made to the pictures. Data availability All data generated or analysed during this study are included in this published article (and its supplementary information files online). Received: 28 July 2020; Accepted: 13 April 2021

1. Swartz, M. A. & Lund, A. W . Lymphatic and interstitial flow in the tumour microenvironment: linking mechanobiology with immunity. Nat. Rev. Cancer 12(3), 210–219 (2012). 2. Lund, A. W . et al. VEGF-C promotes immune tolerance in B16 melanomas and cross-presentation of tumor antigen by lymph node lymphatics. Cell Rep. 1(3), 191–199 (2012). 3. Cohen, J. N. et al. Lymph node-resident lymphatic endothelial cells mediate peripheral tolerance via Aire-independent direct antigen presentation. J. Exp. Med. 207(4), 681–688 (2010). 4. Tewalt, E. F . et al. Lymphatic endothelial cells induce tolerance via PD-L1 and lack of costimulation leading to high-level PD-1 expression on CD8 T cells. Blood 120(24), 4772–4782 (2012). 5. Hirosue, S. et al. Steady-state antigen scavenging, cross-presentation, and CD8+ T cell priming: a new role for lymphatic endothelial cells. J. Immunol. 192(11), 5002–5011 (2014). 16 Vol:.(1234567890)Scientific Reports | (2021) 11:9666 | https://doi.org/10.1038/s41598-021-88508-w www.nature.com/scientificreports/ 6. Fournel, S. et al. Cutting edge: soluble HLA-G1 triggers CD95/CD95 ligand-mediated apoptosis in activated CD8+ cells by inter- acting with CD8. J. Immunol. 164(12), 6100–6104 (2000). 7. Kyurkchiev, D., Ivanova-Todorova, E. & Kyurkchiev, S. D. New target cells of the immunomodulatory effects of progesterone. Reprod. Biomed. Online 21(3), 304–311 (2010). 8. Szekeres-Bartho, J. & Wegmann, T. G. A progesterone-dependent immunomodulatory protein alters the Th1/Th2 balance. J. Reprod. Immunol. 31(1–2), 81–95 (1996). 9. Szekeres-Bartho, J., Par, G., Szereday, L., Smart, C. Y . & Achatz, I. Progesterone and non-specific immunologic mechanisms in pregnancy. Am. J. Reprod. Immunol. 38(3), 176–182 (1997). 10. Szekeres-Bartho, J. et al. Progesterone-dependent immunomodulation. Chem. Immunol. Allergy 89, 118–125 (2005). 11. Cohen, R. A., Check, J. H. & Dougherty, M. P . Evidence that exposure to progesterone alone is a sufficient stimulus to cause a precipitous rise in the immunomodulatory protein the progesterone induced blocking factor (PIBF). J. Assist. Reprod. Genet. 33(2), 221–229 (2016). 12. Shah, N. M., Lai, P . F ., Imami, N. & Johnson, M. R. Progesterone-Related Immune Modulation of Pregnancy and Labor. Front. Endocrinol. (Lausanne) 10, 198 (2019). 13. Srivastava, M. D., Thomas, A., Srivastava, B. I. & Check, J. H. Expression and modulation of progesterone induced blocking fac - tor (PIBF) and innate immune factors in human leukemia cell lines by progesterone and mifepristone. Leuk. Lymphoma. 48(8), 1610–1617 (2007). 14. Pluchino, N., Ansaldi, Y . & Genazzani, A. R. Brain intracrinology of allopregnanolone during pregnancy and hormonal contracep- tion. Horm. Mol. Biol. Clin. Investig. 37(1) (2019). 15. Pelegrina, L. T. et al. Effect of progesterone and first evidence about allopregnanolone action on the progression of epithelial human ovarian cancer cell lines. J. Steroid Biochem. Mol. Biol. 196, 105492 (2020). 16. Chantilis, S., Dombroski, R., Shackleton, C. H., Casey, M. L. & MacDonald, P . C. Metabolism of 5 alpha-dihydroprogesterone in women and men: 3 beta- and 3 alpha-,6 alpha-dihydroxy-5 alpha-pregnan-20-ones are major urinary metabolites. J. Clin. Endo - crinol. Metab. 81(10), 3644–3649 (1996). 17. Horwitz, K. B., Pike, A. W ., Gonzalez-Aller, C. & Fennessey, P . V . Progesterone metabolism in T47Dco human breast cancer cells–II Intracellular metabolic path of progesterone and synthetic progestins. J. Steroid Biochem. 25(6), 911–916 (1986). 18. Norberg, L., Wahlstrom, G. & Backstrom, T. The anaesthetic potency of 3 alpha-hydroxy-5 alpha-pregnan-20-one and 3 alpha- hydroxy-5 beta-pregnan-20-one determined with an intravenous EEG-threshold method in male rats. Pharmacol. Toxicol. 61(1), 42–47 (1987). 19. Parizek, A. et al. Neuroactive pregnanolone isomers during pregnancy. J. Clin. Endocrinol. Metab. 90(1), 395–403 (2005). 20. Hill, M. et al. Circulating levels of pregnanolone isomers during the third trimester of human pregnancy. J. Steroid Biochem. Mol. Biol. 105(1–5), 166–175 (2007). 21. Brunton, P . J., Russell, J. A. & Hirst, J. J. Allopregnanolone in the brain: protecting pregnancy and birth outcomes. Prog. Neurobiol. 113, 106–136 (2014). 22. Y awno, T., Y an, E. B., Walker, D. W . & Hirst, J. J. Inhibition of neurosteroid synthesis increases asphyxia-induced brain injury in the late gestation fetal sheep. Neuroscience 146(4), 1726–1733 (2007). 23. Zamora-Sanchez, C. J. et al. Allopregnanolone alters the gene expression profile of human glioblastoma cells. Int. J. Mol. Sci. 19(3), 684 (2018). 24. Dombroski, R., Casey, M. L. & Macdonald, P . C. Human saturated steroid 6alpha-hydroxylase. J. Clin. Endocrinol. Metab. 82(5), 1338–1344 (1997). 25. O’Hare, M. J., Nice, E. C., McIlhinney, R. A. & Capp, M. Progesterone synthesis, secretion and metabolism by human teratoma- derived cell-lines. Steroids 38(6), 719–737 (1981). 26. Antila, E., Fotsis, T., Wartiovaara, J. & Adlercreutz, H. Steroid metabolism in human teratocarcinoma cell line PA 1. J. Steroid Biochem. 19(5), 1583–1590 (1983). 27. Fennessey, P . V ., Pike, A. W ., Gonzalez-Aller, C. & Horwitz, K. B. Progesterone metabolism in T47Dco human breast cancer cells–I. 5 alpha-pregnan-3 beta,6 alpha-diol-20-one is the secreted product. J. Steroid Biochem. 25(5A), 641–648 (1986). 28. Hilborn, E., Stal, O., Alexeyenko, A. & Jansson, A. The regulation of hydroxysteroid 17beta-dehydrogenase type 1 and 2 gene expression in breast cancer cell lines by estradiol, dihydrotestosterone, microRNAs, and genes related to breast cancer. Oncotarget 8(37), 62183–62194 (2017). 29. Wang, C. T. et al. High expression of 17beta-hydroxysteroid dehydrogenase type 2 is associated with a better prognosis in urothelial carcinoma of the urinary tract. J. Cancer 7(15), 2221–2230 (2016). 30. Lamb, J. C., Levy, M. A., Johnson, R. K. & Isaacs, J. T. Response of rat and human prostatic cancers to the novel 5 alpha-reductase inhibitor, SK&F 105657. Prostate 21(1), 15–34 (1992). 31. Khan, M. A. & Partin, A. W . Finasteride and prostate cancer. Rev. Urol. 6(2), 97–98 (2004). 32. Liss, M. A. & Thompson, I. M. Prostate cancer prevention with 5-alpha reductase inhibitors: concepts and controversies. Curr. Opin. Urol. 28(1), 42–45 (2018). 33. Pinacho-Garcia, L. M. et al. The effect of finasteride and dutasteride on the synthesis of neurosteroids by glioblastoma cells. Steroids 155, 108556 (2020). 34. Byrns, M. C. Regulation of progesterone signaling during pregnancy: implications for the use of progestins for the prevention of preterm birth. J. Steroid Biochem. Mol. Biol. 139, 173–181 (2014). 35. Nowak, F . V . Distribution and metabolism of 20 alpha-hydroxylated progestins in the female rat. J. Steroid Biochem. Mol. Biol. 80(4–5), 469–479 (2002). 36. Wu, L. et al. Expression cloning and characterization of human 17 beta-hydroxysteroid dehydrogenase type 2, a microsomal enzyme possessing 20 alpha-hydroxysteroid dehydrogenase activity. J. Biol. Chem. 268(17), 12964–12969 (1993). 37. Moeller, G. & Adamski, J. Multifunctionality of human 17beta-hydroxysteroid dehydrogenases. Mol. Cell Endocrinol. 248(1–2), 47–55 (2006). 38. Mindnich, R., Moller, G. & Adamski, J. The role of 17 beta-hydroxysteroid dehydrogenases. Mol. Cell Endocrinol. 218(1–2), 7–20 (2004). 39. Hilborn, E., Stal, O. & Jansson, A. Estrogen and androgen-converting enzymes 17beta-hydroxysteroid dehydrogenase and their involvement in cancer: with a special focus on 17beta-hydroxysteroid dehydrogenase type 1, 2, and breast cancer. Oncotarget 8(18), 30552–30562 (2017). 40. Hanukoglu, I. Steroidogenic enzymes: structure, function, and role in regulation of steroid hormone biosynthesis. J. Steroid Bio - chem. Mol. Biol. 43(8), 779–804 (1992). 41. Voutilainen, R., Picado-Leonard, J., DiBlasio, A. M. & Miller, W . L. Hormonal and developmental regulation of adrenodoxin mes- senger ribonucleic acid in steroidogenic tissues. J. Clin. Endocrinol. Metab. 66(2), 383–388 (1988). 42. Fayard, E., Auwerx, J. & Schoonjans, K. LRH-1: an orphan nuclear receptor involved in development, metabolism and steroido- genesis. Trends Cell Biol. 14(5), 250–260 (2004). 43. Szekeres-Bartho, J. et al. The role of gamma/delta T cells in progesterone-mediated immunomodulation during pregnancy: a review. Am. J. Reprod. Immunol. 42(1), 44–48 (1999). 44. Dombroski, R. A., Casey, M. L. & MacDonald, P . C. 5-Alpha-dihydroprogesterone formation in human placenta from 5alpha- pregnan-3beta/alpha-ol-20-ones and 5-pregnan-3beta-yl-20-one sulfate. J. Steroid Biochem. Mol. Biol. 63(1–3), 155–163 (1997). 17 Vol.:(0123456789)Scientific Reports | (2021) 11:9666 | https://doi.org/10.1038/s41598-021-88508-w www.nature.com/scientificreports/ 45. Szekeres-Bartho, J. & Polgar, B. PIBF: the double edged sword. Pregnancy and tumor. Am. J. Reprod. Immunol. 64(2), 77–86 (2010). 46. Szekeres-Bartho, J. et al. Immunoregulatory effects of a suppressor factor from healthy pregnant women’s lymphocytes after pro- gesterone induction. Cell Immunol. 122(2), 281–294 (1989). 47. Polgar, B., Nagy, E., Miko, E., Varga, P . & Szekeres-Bartho, J. Urinary progesterone-induced blocking factor concentration is related to pregnancy outcome. Biol. Reprod. 71(5), 1699–1705 (2004). 48. Hudic, I. et al. Lower urinary and serum progesterone-induced blocking factor in women with preterm birth. J. Reprod. Immunol. 117, 66–69 (2016). 49. Szekeres-Bartho, J. et al. The mechanism of the inhibitory effect of progesterone on lymphocyte cytotoxicity: I. Progesterone- treated lymphocytes release a substance inhibiting cytotoxicity and prostaglandin synthesis. Am. J. Reprod. Immunol. Microbiol. 9(1), 15–18 (1985). 50. Hudic, I. et al. Progesterone-induced blocking factor and cytokine profile in women with threatened pre-term delivery. Am. J. Reprod. Immunol. 61(5), 330–337 (2009). 51. Hudic, I. & Fatusic, Z. Progesterone - induced blocking factor (PIBF) and Th(1)/Th(2) cytokine in women with threatened spon- taneous abortion. J. Perinat. Med. 37(4), 338–342 (2009). 52. Gutierrez-Rodriguez, A., Hansberg-Pastor, V . & Camacho-Arroyo, I. Proliferative and invasive effects of progesterone-induced blocking factor in human glioblastoma cells. Biomed. Res. Int. 2017, 1295087 (2017). 53. Anderle, C. et al. Human trophoblast cells express the immunomodulator progesterone-induced blocking factor. J. Reprod. Immu- nol. 79(1), 26–36 (2008). 54. Meyer, C., Schmid, R., Scriba, P . C. & Wehling, M. Purification and partial sequencing of high-affinity progesterone-binding site(s) from porcine liver membranes. Eur. J. Biochem. 239(3), 726–731 (1996). 55. Kaluka, D., Batabyal, D., Chiang, B. Y ., Poulos, T. L. & Y eh, S. R. Spectroscopic and mutagenesis studies of human PGRMC1. Biochemistry 54(8), 1638–1647 (2015). 56. Zhang, L. et al. Expression of progesterone receptor membrane component 1 and its partner serpine 1 mRNA binding protein in uterine and placental tissues of the mouse and human. Mol. Cell Endocrinol. 287(1–2), 81–89 (2008). 57. Keator, C. S., Mah, K. & Slayden, O. D. Alterations in progesterone receptor membrane component 2 (PGRMC2) in the endome- trium of macaques afflicted with advanced endometriosis. Mol. Hum. Reprod. 18(6), 308–319 (2012). 58. Lissauer, D. et al. Progesterone promotes maternal-fetal tolerance by reducing human maternal T-cell polyfunctionality and induc- ing a specific cytokine profile. Eur. J. Immunol. 45(10), 2858–2872 (2015). 59. Areia, A. et al. Membrane progesterone receptors in human regulatory T cells: a reality in pregnancy. BJOG 122(11), 1544–1550 (2015). 60. Lu, J., Reese, J., Zhou, Y . & Hirsch, E. Progesterone-induced activation of membrane-bound progesterone receptors in murine macrophage cells. J. Endocrinol. 224(2), 183–194 (2015). 61. Rubinow, K. B. An intracrine view of sex steroids, immunity, and metabolic regulation. Mol. Metab. 15, 92–103 (2018). 62. Buttgereit, F . & Scheffold, A. Rapid glucocorticoid effects on immune cells. Steroids 67(6), 529–534 (2002). 63. Csapo, A. I. & Wiest, W . G. An examination of the quantitative relationship between progesterone and the maintenance of preg- nancy. Endocrinology 85(4), 735–746 (1969). 64. Schiff, R. I., Mercier, D. & Buckley, R. H. Inability of gestational hormones to account for the inhibitory effects of pregnancy plasmas on lymphocyte responses in vitro. Cell Immunol. 20(1), 69–80 (1975). 65. Machnik, A. et al. Macrophages regulate salt-dependent volume and blood pressure by a vascular endothelial growth factor-C- dependent buffering mechanism. Nat. Med. 15(5), 545–552 (2009). 66. Chen, X. et al. Differential expression of vascular endothelial growth factor angiogenic factors in different endometrial compart- ments in women who have an elevated progesterone level before oocyte retrieval, during in vitro fertilization-embryo transfer treatment. Fertil. Steril. 104(4), 1030–1036 (2015). 67. Zhang, Y . et al. Progesterone metabolism in adipose cells. Mol. Cell Endocrinol. 298(1–2), 76–83 (2009).

We thank Sofia Bachmann and the Centre of Laboratory Medicine, Inselspital Bern for their help with the immunologic experiments, and Sophie Braga Lagache and Natasha Buchs for the proteomics analysis. CGM has drawn Figure 4 and Supplementary Figure 2. Author contributions R.K. performed experiments and reviewed the manuscript. M.G. designed the steroidomic analysis, analyzed the LC–MS results and edited the manuscript, N.S. performed the LC–MS experiments, M.F . conducted the FACS experiments and contributed to scientific discussions, G.E. and S.V . provided the mouse lymph nodes and adrenal glands and edited the manuscript. H.J. performed experiments. B.V . reviewed the manuscript and contributed to scientific discussions. M.G.M. supported with scientific discussions, edited and reviewed the manuscript. C.G.M. designed and analyzed the experiments, performed some and wrote the manuscript. Funding This research was supported by the Swiss National Science Foundation (personal grant 32-135596 to MGM) and B. Vogt is supported by «Fonds pour la recherche thérapeutique, Lausanne, Switzerland». Competing interests The authors declare no competing interests. Additional information Supplementary Information The online version contains supplementary material available at https:// doi. org/ 10. 1038/ s41598- 021- 88508-w. Correspondence and requests for materials should be addressed to C.G.-M. Reprints and permissions information is available at www.nature.com/reprints. Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. 18 Vol:.(1234567890)Scientific Reports | (2021) 11:9666 | https://doi.org/10.1038/s41598-021-88508-w www.nature.com/scientificreports/ Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/. © The Author(s) 2021