Modulation of Fibroblast Activity via Vitamin D3 Is Dependent on Tumor Type—Studies on Mouse Mammary Gland Cancer

Simple Summary This study, which was conducted in healthy mice and mice bearing three mouse mammary gland cancers—4T1, 67NR, and E0771—showed that the divergent effects of vitamin D3 supplementation (5000 IU) or deficiency (100 IU of vitamin D3) observed in healthy mice led to the formation of various body microenvironments depending on the mouse strain. Developing tumors themselves modified the microenvironments by producing higher concentrations of osteopontin, SDF-1 (4T1), TGF-β (4T1 and E0771), CCL2, VEGF, FGF23 (E0771), and IL-6 (67NR), which influences the response to vitamin D3 supplementation/deficiency and calcitriol administration and leads to enhanced/decreased activation of lung fibroblasts and modulation of tumor tissue blood flow. Abstract Vitamin D3 and its analogs are known to modulate the activity of fibroblasts under various disease conditions. However, their impact on cancer-associated fibroblasts (CAFs) is yet to be fully investigated. The aim of this study was to characterize CAFs and normal fibroblasts (NFs) from the lung of mice bearing 4T1, 67NR, and E0771 cancers and healthy mice fed vitamin-D3-normal (1000 IU), -deficient (100 IU), and -supplemented (5000 IU) diets. The groups receiving control (1000 IU) and deficient diets (100 IU) were gavaged with calcitriol (+cal). In the 4T1-bearing mice from the 100 IU+cal group, increased NFs activation (increased α-smooth muscle actin, podoplanin, and tenascin C (TNC)) with a decreased blood flow in the tumor was observed, whereas the opposite effect was observed in the 5000 IU and 100 IU groups. CAFs from the 5000 IU group of E0771-bearing mice were activated with increased expression of podoplanin, platelet-derived growth factor receptor β, and TNC. In the 100 IU+cal group of E0771-bearing mice, a decreased blood flow was recorded with decreased expression of fibroblast growth factor 23 (FGF23) and C-C motif chemokine ligand 2 (CCL2) in tumors and increased expression of TNC on CAFs. In the 67NR model, the impact of vitamin D3 on blood flow or CAFs and lung NFs was not observed despite changes in plasma and/or tumor tissue concentrations of osteopontin (OPN), CCL2, transforming growth factor-β, vascular endothelial growth factor, and FGF23. In healthy mice, divergent effects of vitamin D3 supplementation/deficiency were observed, which lead to the creation of various body microenvironments depending on the mouse strain. Tumors developing in such microenvironments themselves modified the microenvironments by producing, for example, higher concentrations of OPN and stromal-cell-derived factor 1 (4T1), which influences the response to vitamin D3 supplementation/deficiency and calcitriol administration.


Introduction
The tumor microenvironment, which includes fibroblasts, vascular endothelial cells, and immune cells, is an important component of cancer and affects its progression, metastasis, and sensitivity to therapies. Vitamin D 3 -mainly its active form, calcitriol, 1,25-dihydroxycholecalciferol-can influence almost every cell in the body, including the tumor-building ones [1]. Thus, although cancer cells (e.g., murine cancer cell lines) are not sensitive to the antiproliferative activity of vitamin D 3 in vitro (despite the presence of the vitamin D receptor, VDR), they show significant sensitivity to vitamin D derivatives in vivo. One such example is 4T1 mouse breast cancer cells. Previous studies have shown that calcitriol and its analogs may increase the rate of tumor metastasis [2,3] or accelerate the progression of primary tumors in young mice [4]. However, in aged, ovariectomized animals (as a postmenopausal model), these compounds inhibited 4T1 tumor metastasis [5]. Besides the age of mice, other factors may also affect the outcome of the in vivo treatment of mammary gland tumors with vitamin D compounds. The abovementioned observations were from experiments in mice in which the therapy was started when tumors were already formed. On the contrary, an antimetastatic effect was observed in young mice bearing 4T1 tumors when calcitriol treatment was started prior to tumor inoculation [6]. The beneficial effect of vitamin D analogs in young mice was also evident when these substances were used in combination therapy with cyclophosphamide [7]. In the case of E0771, another mouse mammary gland cancer model that is also not sensitive to the antiproliferative effect of calcitriol [3], the antitumor effect of cholecalciferol was observed in normal young mice, whereas protumor activity was observed in obese mice. This effect was related to the impact of the treatment on tumor-infiltrating CD8+ T lymphocytes which are increased in number in normal mice and decreased in obese mice [8].
In recent clinical studies much attention has been paid on the effects of vitamin D in the development of various cancers, including breast cancer. These studies focused mainly on the relationship between tumor progression and plasma 25(OH)D concentrations (the most studied metabolite of vitamin D), vitamin D supplementation or the expression of various proteins related to vitamin D activity or metabolism (recently reviewed in [9]). However, based on these studies it is difficult to draw firm conclusions concerning beneficial or adverse effects of vitamin D on breast cancer development. For example, a systematic review and meta-analysis by Hossain et al. showed an inverse relationship between vitamin D intake and breast cancer occurrence along with a direct relationship between vitamin D deficiency and breast cancer risk [10]. Other studies reported the occurrence of vitamin D deficiency in newly diagnosed breast cancer patients [11]. They also indicate that vitamin D supplementation significantly reduced total cancer mortality but did not affect total cancer incidence [12]. On the other hand, e.g., O'Connor et al. indicated the lack of the relationship between vitamin D and breast cancer [13]. Moreover, Kanstrup et al. [14] found an inverse correlation between 25(OH)D plasma concentrations and breast cancer survival, however, this study also indicated poorer breast cancer survival among patients with high 25(OH)D concentrations [14]. Likewise, Ganji et al. reported an association of high 25(OH)D plasma concentrations (75-100 nmol/L) with greater risk of breast cancer in postmenopausal women [14,15]. Therefore, controlled trials are required to verify the data concerning vitamin D benefits for breast cancer patients.
Cancer-associated fibroblasts (CAFs) resembling myofibroblasts-activated spindleshaped fibroblasts-are the major cellular component in the tumor microenvironment [16]. The stimuli produced by tumor cells and stromal cells, such as fibroblasts, in the tumor microenvironment can activate fibroblasts and contain, among others, interleukin-1, interleukin-6 (IL-6), bone morphogenetic protein, sonic hedgehog, reactive oxygen species, transforming growth factor-β (TGF-β), platelet-derived growth factor, and tumor necrosis factor [17]. The activation of fibroblasts by these stimuli may result in quiescent, tumorrestraining, and tumor-promoting CAFs [18][19][20]. Orimo et al. showed that coimplantation of the tumor xenograft model with CAFs harvested from the breast carcinomas of patients promoted higher growth of breast cancer cells compared with coimplantation with normal mammary fibroblasts derived from the same patients [21]. CAFs promote tumor growth through the secretion of stromal-cell-derived factor 1 (SDF-1) and angiogenesis by recruiting endothelial progenitor cells into tumor tissue [21]. Thus, some studies have proposed targeting CAFs as a treatment option for triple-negative breast cancer [22]. Clinical retrospective studies on breast cancer patients focusing on the expression of caveolin-1, which acts as a tumor-suppressing molecule in CAFs, showed 72 months of cancer-specific survival in caveolin-1-positive patients, whereas only 29.5 months in caveolin-1-negative patients [23].
Osteopontin (OPN) is an important protein engaged in the crosstalk between cancer cells and stromal fibroblasts. OPN secreted by breast cancer cells induces differentiation of fibroblasts to myofibroblasts. Furthermore, OPN-driven CAFs secrete SDF-1, which in turn triggers epithelial-to-mesenchymal transition in tumor cells [24]. As reported in previous studies, calcitriol and other agonists of VDR stimulate the transcription of the OPN gene (Spp1) or secretion of OPN [25,26]. In our previous in vitro study, calcitriol and its analogs stimulated OPN secretion in 67NR mouse mammary gland cancer cells, which are sensitive to proliferation inhibition by calcitriol, but not in 4T1 cells [2]. It was also found that in vitro stimulation by calcitriol and its analogs led to increased OPN secretion by normal mouse fibroblast cells BALB/3T3, but not by the murine macrophage cell line RAW 264.7 [27]. Although much research has been conducted investigating the effects of vitamin D on fibroblasts in various types of diseases or physiological conditions [28][29][30][31], its effect on CAFs remains to be elucidated. Studies on fibroblasts derived from colon cancer [32,33], pancreatic cancer [34,35], and breast cancer [36,37] patients have indicated that calcitriol modulates gene expression and inhibits the protumoral properties of CAFs. Therefore, in the present study, we analyzed the effects of calcitriol treatment and vitamin D 3 (cholecalciferol) deprivation or supplementation on tumor and lung fibroblasts in three mouse mammary gland cancer models (4T1, 67NR, and E0771) with different metastatic potential.

Cells and Tissues Harvested from Mice Bearing Mammary Gland Tumors
Fibroblasts were isolated from lung and tumor tissues of mice bearing 4T1, 67NR, and E0771 tumors and of healthy BALB/c and C57BL/6 mice. Data regarding tumor growth, metastasis, and vitamin D metabolism in these mice have been previously published [3]. Origin of the cell lines used: 4T1 cells-the American Type Culture Collection (ATCC, Rockville, MD, USA); 67NR (nonmetastatic counterparts of 4T1) Barbara Ann Karmanos Cancer Institute (Detroit, MI, USA); E0771 cell line [38], gifted by Dr. Andreas Möller (School of Medicine, University of Queensland; Tumour Microenvironment Laboratory, QIMR Berghofer Medical Research Institute, Herston, Queensland, Australia).
The complete in vivo experimental design has been described earlier [3]. The animal study protocol was approved by the first Local Committee for Experiments with the Use of Laboratory Animals, Hirszfeld Institute of Immunology and Experimental Therapy, Wroclaw, Poland (permission number: 66/2018, 18 July 2018). Mice were divided into groups and fed ad libitum with an AIN67 synthetic diet (ZooLab, Sedziszow, Poland) for 6 weeks. The diet fed to groups had a varied content of vitamin D 3 as follows: control amount of vitamin D 3 (1000 IU/kg), supplemented vitamin D 3 (5000 IU/kg), and deficient vitamin D 3 (100 IU/kg). After 6 weeks (day 0), some mice were implanted orthotopically with tumor cells (1 × 104, 2 × 105, and 5 × 104 viable 4T1, 67NR, and E0771 cells, respectively), and the same diet was continued. Seven days after the implantation of tumor cells, calcitriol administration (1 µg/kg per os by gavage; p.o.) was started, which continued thrice a week in the groups receiving food with a normal, control amount of vitamin D 3 (1000 IU) and in the groups receiving a vitamin-D 3 -deficient diet (100 IU). On day 23 (C57BL/6) or 28 (BALB/c) after the inculcation of cancer cells, mice were subjected to isoflurane anesthesia and injected with buprenorphine (0.1 mg/kg of body weight) analgesia, subjected to blood collection, and then killed ( Figure 1). tively), and the same diet was continued. Seven days after the implantation of tumor cells, calcitriol administration (1 µg/kg per os by gavage; p.o.) was started, which continued thrice a week in the groups receiving food with a normal, control amount of vitamin D3 (1000 IU) and in the groups receiving a vitamin-D3-deficient diet (100 IU). On day 23 (C57BL/6) or 28 (BALB/c) after the inculcation of cancer cells, mice were subjected to isoflurane anesthesia and injected with buprenorphine (0.1 mg/kg of body weight) analgesia, subjected to blood collection, and then killed ( Figure 1). Figure 1. Scheme of the treatment of BALB/c and C57BL/6 mice. Mice of both strains were fed ad libitum with diets containing various amounts of vitamin D3 for 6 weeks. Next, on the day assigned as day 0, tumor cells were implanted orthotopically in selected mice (4T1 and 67NR to BABL/c mice, E0771 to C57BL/6 mice). The diet was continued for the next 7 days, and then calcitriol administration per os was started and continued thrice a week till day 28 (for BALB/c mice) or day 23 (for C57BL/6 mice). Calcitriol was administered to mice on control (1000 IU) and deficient (100 IU) diets at a dose of 1 µg/kg. During autopsy, blood, tumors, and lungs were harvested for further analyses.

Blood Flow Assessment
Tumor blood perfusion analysis was performed on day 21 (4T1 and 67NR models) or on day 19 (E0771 model) using the MicroMarker™ contrast agent (VisualSonics, Toronto, ON, Canada) and Vevo 2100 ultrasound imaging system (VisualSonics). Mice were anesthetized by continuous administration of 2-3% (v/v) isoflurane (Baxter, Deerfield, Germany) in synthetic air (600 mL/min) and immobilized. The tumor area was enclosed with air-bubble-free gel, and the central cross-section of the tumor was visualized in the transverse plane using an MS250 scanhead (VisualSonics). Then, 100 µL of the contrast agent was injected intravenously and the first imaging sequence (bolus) was recorded (ca. 15 fps, at least 1000 frames) after the contrast signal in the tumor tissue reached the steady state (ca. 50 s). Next, using the burst mode, contrast microbubbles were destroyed, and the second imaging sequence (replenishment) was recorded. Mice were maintained in a warm environment until fully awakened. The data were analyzed using the Vevo LAB 1.7.1 Software with VevoCQ modality (VisualSonics).

Tissue Microarrays
Tissue microarrays were prepared using TMA Grand Master (3DHistech, Budapest, Hungary), in accordance with the manufacturer's instructions. The representative spots for tissue microarrays were selected by a pathologist. Three core punches of 1.5 mm diameter from each tumor block were transferred into the recipient paraffin block.

Immunohistochemical Staining of Tissues
Using the automated immunostainer device Autostainer Link 48 (Dako, Glostrup, Denmark), immunohistochemistry analysis was performed on 4 µm formalin-fixed paraffin-embedded tissue microarray sections. To deparaffinize, rehydrate, and unmask the epitopes, the slides were boiled in EnVision™ FLEX Target Retrieval Solution of a high pH for 20 min at 97 °C using PT-Link (both from Dako). The endogenous peroxidase activity was blocked using the EnVision FLEX Peroxidase-Blocking Reagent (Dako). Then, the primary antibody against collagen type I alpha 1 chain (COL1A1) (dilution 1:1600, Figure 1. Scheme of the treatment of BALB/c and C57BL/6 mice. Mice of both strains were fed ad libitum with diets containing various amounts of vitamin D 3 for 6 weeks. Next, on the day assigned as day 0, tumor cells were implanted orthotopically in selected mice (4T1 and 67NR to BABL/c mice, E0771 to C57BL/6 mice). The diet was continued for the next 7 days, and then calcitriol administration per os was started and continued thrice a week till day 28 (for BALB/c mice) or day 23 (for C57BL/6 mice). Calcitriol was administered to mice on control (1000 IU) and deficient (100 IU) diets at a dose of 1 µg/kg. During autopsy, blood, tumors, and lungs were harvested for further analyses.

Blood Flow Assessment
Tumor blood perfusion analysis was performed on day 21 (4T1 and 67NR models) or on day 19 (E0771 model) using the MicroMarker™ contrast agent (VisualSonics, Toronto, ON, Canada) and Vevo 2100 ultrasound imaging system (VisualSonics). Mice were anesthetized by continuous administration of 2-3% (v/v) isoflurane (Baxter, Deerfield, Germany) in synthetic air (600 mL/min) and immobilized. The tumor area was enclosed with air-bubble-free gel, and the central cross-section of the tumor was visualized in the transverse plane using an MS250 scanhead (VisualSonics). Then, 100 µL of the contrast agent was injected intravenously and the first imaging sequence (bolus) was recorded (ca. 15 fps, at least 1000 frames) after the contrast signal in the tumor tissue reached the steady state (ca. 50 s). Next, using the burst mode, contrast microbubbles were destroyed, and the second imaging sequence (replenishment) was recorded. Mice were maintained in a warm environment until fully awakened. The data were analyzed using the Vevo LAB 1.7.1 Software with VevoCQ modality (VisualSonics).

Tissue Microarrays
Tissue microarrays were prepared using TMA Grand Master (3DHistech, Budapest, Hungary), in accordance with the manufacturer's instructions. The representative spots for tissue microarrays were selected by a pathologist. Three core punches of 1.5 mm diameter from each tumor block were transferred into the recipient paraffin block.

Immunohistochemical Staining of Tissues
Using the automated immunostainer device Autostainer Link 48 (Dako, Glostrup, Denmark), immunohistochemistry analysis was performed on 4 µm formalin-fixed paraffinembedded tissue microarray sections. To deparaffinize, rehydrate, and unmask the epitopes, the slides were boiled in EnVision™ (Santa Clara, CA, USA) FLEX Target Retrieval Solution of a high pH for 20 min at 97 • C using PT-Link (both from Dako). The endogenous peroxidase activity was blocked using the EnVision FLEX Peroxidase-Blocking Reagent (Dako). Then, the primary antibody against collagen type I alpha 1 chain (COL1A1) (dilution 1:1600, cat.no. PA5-86949; Invitrogen, Waltham, MA, USA) or against the smooth muscle actin (ready-to-use, cat. no. IR611; Dako, Glostrup, Denmark) were applied for 20 min at room temperature. Afterward, the slides were incubated with EnVision FLEX/HRP (Dako) for (20 min). Subsequently, the sections were incubated for 10 min at room temperature with the substrate for peroxidase, diaminobenzidine. Additionally, using the EnVision FLEX Hematoxylin (Dako), all slides were counterstained for 5 min. After dehydration in graded ethanol concentrations (70%, 96%, absolute) and in xylene, all slides were covered with coverslips in the SUB-X Mounting Medium. Two pathologists evaluated the expression of studied proteins with a routinely used immunoreactive scale of Remmele and Stegner [39], which is presented in Table 1.

Magnetic Separation
For magnetic separation, whole-cell suspensions (three mice from one group were pooled together for one separation) were used. CAFs were isolated using Anti-Fibroblast MicroBeads (Miltenyi Biotec, Auburn, CA, USA), according to the manufacturer's protocol. In brief, centrifuged pellets (7-10 min, 4 • C, 350× g) were resuspended in the separation buffer containing PBS of pH 7.2, 0.5% (v/v) bovine serum albumin, 2 mM ethylene diamine tetra-acetic acid (Sigma-Aldrich, Saint-Louis, MO, USA), and TruStain FcX (antimouse CD16/CD 3 2) antibody (BioLegend, San Diego, CA, USA) and then incubated for 10 min at 4 • C in order to block the Fc receptors (0.1 µg/100 µL volume). After blocking, magnetic beads were added (20 µL per 10 6 of total cells) to the cells and incubated for 30 min in the dark at room temperature. Subsequently, 1 mL of separation buffer was added to the cells and centrifuged (7-10 min, 4 • C, 350× g). The pellets were resuspended in 1 mL of separation buffer and applied onto activated MS columns (Miltenyi Biotec, Auburn, CA, USA) placed in the magnetic field of the MiniMACS Separator (Miltenyi Biotec, Auburn, CA, USA). After three washes with 500 µL of separation buffer, the columns were transferred into new sterile 10 mL tubes, and cells were flushed into collection tubes. The collected Cancers 2022, 14, 4585 6 of 25 cells were counted using a Bürker chamber in trypan blue solution (0.4% (w/v)) and used for further analysis.

Flow Cytometry Analysis of Isolated Fibroblasts
CAFs isolated from tumor tissues were analyzed immediately after the separation. Lung fibroblasts from tumor-bearing mice were processed directly after thawing, and those from healthy mice were subjected to cytometry analysis after stimulation with TGF-β1 or cancer-cell-conditioned medium. To assess cell viability, cell suspensions of 3-5 ×  ). The dishes were washed with PBS after the cells reached confluence, and cells were incubated for 24 h with an appropriate medium without serum. Then, CM were collected, centrifuged (7 min, 4 • C, 400× g), and applied in a 1:1 ratio (CM:NFs medium) for 72 h on NF cultures for stimulation.

Fluorescence Microscopy Analyses
A total of 2.5 × 10 3 cells/well were cultured for imaging on a Falcon ® 96-well Black/Clear Flat Bottom TC-treated Imaging Microplate (Corning, New York, NY, USA). Staining was carried out after stimulation with TGF-β1 or cancer-cell-conditioned medium. Before staining, cells were washed twice with PBS solution, fixed in freshly prepared 4% (v/v) paraformaldehyde (Avantor Performance Materials Poland, Gliwice, Poland) for 10-15 min, washed thrice with PBS, and permeabilized in 0.25% (v/v) Triton X-100 (Sigma-Aldrich, Saint-Louis, MO, USA) for 15 min at room temperature (only wells for VDR staining). After washing thrice in the PBS solution, cells were blocked for 30 min in 1% (w/v) bovine serum albumin (Sigma-Aldrich, Saint-Louis, MO, USA) solution in 0.1% (v/v) PBS/Tween 20 (Sigma-Aldrich, Saint-Louis, MO, USA) at room temperature. Then, the fixed cells were rinsed thrice with PBS (5 min each) and incubated with primary antibodies against VDR (bs-2987R, dilution 1:100; Bioss Antibodies, Woburn, MA, USA) and fibroblast activation protein (FAP) (ab28244, dilution 1:500; Abcam, Cambridge, UK) in blocking solution at 4 • C overnight. Next, after washing thrice with PBS, a secondary antibody (Anti-rabbit Antibody Alexa Fluor 488, ab150077; Abcam, Cambridge, UK) in the blocking solution was used for 1 h at room temperature. Then, the samples were rinsed thrice with PBS and stained with DAPI (dilution 1:1500; Cell-Signaling, Danvers, TX, USA) and DyLight™ 554 Phalloidin (dilution 1:100; Cell-Signaling, Danvers, MA, USA) in PBS solution for 15 min at room temperature. Cells were photographed using an Olympus IX81 fluorescence microscope (Olympus, Warsaw, Poland) with CellSense software (Olympus, Warsaw, Poland). Fluorescence was determined from the images using ImageJ according to the protocol of Luke Hammond (QBI, The University of Queensland, Australia; accessed on 18 October 2021; www.theolb.readthedocs.io/en/latest/imaging/measuringcell-fluorescence-using-imagej.html). The areas of interest (whole cell for FAP staining and cell nucleus for VDR staining) and areas of background were measured, and corrected total cell fluorescence (CTCF) was calculated using the following formula: CTCF = integrated density − (area of the selected cell × mean fluorescence of background readings)

Real-Time qPCR Analysis of Cultured NFs Isolated from Lungs of Healthy BALB/c or C57BL/6 Mice Fed Control Diet and Vitamin-D 3 -Deficient Diet and/or Treated with Calcitriol
From the stimulated ex vivo NFs, total RNA was extracted using 1 mL Tri-reagent (Sigma Aldrich, Saint-Louis, MO, USA) followed by RNA purification with Direct-zol™ RNA Miniprep (ZYMO RESEARCH, Tustin, CA, USA) according to the manufacturer's protocol. Using the iScript cDNA Synthesis Kit (Bio-Rad, Hercules, CA, USA), 1 µg of purified RNA was reverse transcribed into complementary DNA (cDNA). Then, the expression of the following genes was determined in cDNA with ready-to-use primers and probes (TaqMan ® Gene Expression Assays; Thermo Fisher Scientific, Waltham, MA, USA): Acta2 (Mm01546133_m1), Mmp9 (Mm00442991_m1), Spp1 (Mm00436767_m1), and Vdr (Mm00437297_m1). Real-time qPCR was carried out in a 10 µL reaction volume composed of 20× presented probes, 50 ng cDNA, and 2× TaqMan™ Gene Expression Master Mix (Thermo Fisher Scientific, Waltham, MA, USA) in a ViiA™ 7 Real-Time PCR System (Thermo Fisher Scientific, Waltham, MA, USA) using the following program: 10 min at 95 • C for initial denaturation and 40 cycles at 95 • C for 15 s and 60 • C for 1 min. Expression was calculated according to the comparative ∆∆Ct method in which Gadph (Mm99999915_g1) and Rps27a (Mm01180369_g1) were used as endogenous controls and normalized to each untreated control using QuantStudio™ Real-Time PCR Software and ExpressionSuite Software (Thermo Fisher Scientific, Waltham, MA, USA).

Tumor Tissue Preparation for Enzyme-Linked Immunosorbent Assays (ELISA)
Tumor specimens were frozen in liquid nitrogen, suspended in RIPA buffer (Sigma-Aldrich, Saint-Louis, MO, USA) that contained protease and phosphatase inhibitors (both Sigma-Aldrich, Saint-Louis, MO, USA), and then mechanically homogenized (MP Biomedi- cals, Santa Ana, CA, USA). The samples were centrifuged at 10,000× g for 10 min at 4 • C after homogenization, and supernatants were transferred to fresh Eppendorf tubes. Protein concentration in the homogenates was measured using the Quick Start™ Bio-Rad Protein Assay (Bio-Rad, Hercules, CA, USA). Tumor protein samples were used for ELISA.

Statistical Analysis
Statistical analysis was performed using the GraphPad Prism 7.1 software. Using the Shapiro-Wilk data normality test, the normality of the data distribution was analyzed (assumption of significance of the test for p < 0.05). The datasets that did not meet normality requirements (did not pass the Shapiro-Wilk test with α = 0.05) were further tested using the Kruskal-Wallis test using Dunn's post-test for multiple comparisons (if not otherwise stated in the dataset description). Datasets that met the normality requirement were further analyzed using a one-way ANOVA followed by Sidak's post-hoc test for multiple comparisons (since only selected, biologically/experimentally meaningful comparisons were challenged Sidak's post-test provides better correction for multiple comparisons, thus, provides a better test power). The statistical analysis of individual data, depending on their distribution, is presented in table and figure legends. Differences between groups for which p < 0.05 were considered statistically significant.

Impact of Vitamin D on Surface Markers of Lung NFs and CAFs
NFs isolated from the lungs of 4T1 tumor-bearing BALB/c mice (both on normal and on vitamin-D 3 -deficient diet) treated with calcitriol showed significantly increased expression of α-SMA. Compared with the deficiency group, the expression of podoplanin and TNC was increased on fibroblasts from mice fed the vitamin-D 3 -deficient diet and treated with calcitriol. Furthermore, the expression of these two proteins, as well as PDGFRβ, was found to be decreased in mice fed with 5000 IU vitamin-D 3 -supplemented and 100 IU vitamin-D 3 -deficient diets ( Table 2). Administration of different types of diets or calcitriol did not change the expression of these proteins on lung fibroblasts from 67NRbearing BALB/c mice ( Table 2). In C57BL/6 mice bearing E0771 tumors, only podoplanin was decreased in mice fed with diets containing 5000 IU vitamin D 3 and in the 100 IU+cal group (Table 2). Similar characteristics were observed for NFs isolated from healthy BALB/c and C57BL/6 mice, with only a decrease in α-SMA in all treatment groups of BALB/c mice (Supplementary Figure S1).
Differences in vitamin D status did not change the expression of surface antigens on CAFs derived from 4T1 and 67NR tumors (Table 3). However, an exception was the decreased expression of podoplanin on CAFs isolated from 4T1 tumor-bearing mice fed with 5000 IU of vitamin D 3 and mice fed with a control 1000 IU diet and administered with calcitriol by gavage (Table 3). A decrease in α-SMA was observed on CAFs from E0771 tumors isolated from the 1000 IU+cal group, similar to that observed for the 5000 IU group. A reverse tendency was observed for podoplanin and PDGFRβ as the expression of these proteins was increased as compared to the control group. Increased expression was also observed for TNC in all groups supplemented with vitamin D or calcitriol (Table 3).
α-SMA-α-smooth muscle actin; EpCAM-epithelial cell adhesion molecule; PDGFRβ-platelet-derived growth factor receptor beta; TNC-tenascin C. Data are shown as normalized median fluorescence intensity (MFI) calculated by dividing the MFI of the stained sample by the MFI of isotype control. Steps of CD 3 1 − EpCAM − CD45 − fibroblast gating and example dot-plots after staining of CD 3 1 − EpCAM − CD45 − fibroblasts with isotype controls and appropriate antibodies are shown in Supplementary Figure S2. Mice were fed with diets containing various amounts of vitamin D for 6 weeks. Next, on the day assigned as day 0, tumor cells were implanted orthotopically. Diets were continued for the next 7 days, and then calcitriol gavage was started and continued thrice a week till day 23 (for C57BL/6 mice) or day 28 (for BALB/c mice). Calcitriol was administered to mice on control (1000 IU) and deficient (100 IU) diets at a dose of 1 µg/kg. During autopsy, lungs were harvested and NFs were isolated. Flow cytometry analyses were performed on NFs after thawing to determine the expression of α-SMA, podoplanin, PDGFRβ, and TNC. N = 5-7. Statistical analysis: Kruskal-Wallis test followed by Dunn's test for multiple comparisons. Data are shown as mean ± SD. a ≤ 0.05 as compared to 1000 IU, b ≤ 0.05 as compared to 5000 IU, c ≤ 0.05 as compared to 100 IU.

Analysis of Blood Flow in Tumor Tissue from 4T1, 67NR, and E0771 Tumor-Bearing Mice
To evaluate the blood flow in tumor tissue, the following time-intensity curve (TIC) parameters were estimated: peak enhancement (PE)-the maximum intensity in the TIC (blood volume); rise time (RT)-calculated from the beginning of enhancement to PE; mean transit time (mTT)-corresponding to the center of gravity of the best-fit function of echopower (or fitted signal); wash-in area under the TIC (WiAUC); wash-in rate, maximum slope between the time of onset of contrast inflow and the time of PE on the TIC (WiR); wash-in perfusion index (WiPI = WiAUC/RT)-representative of blood flow, relative blood volume (rBV = amplitude of the plateau + offset amplitude), and relative blood flow (rBF = rBV/mTT) ( Table 4, Supplementary Figure S4).
Both vitamin D 3 supplementation and calcitriol treatment influenced the perfusion of 4T1 tumors. Significantly reduced WiR was observed in calcitriol-receiving groups. This indicates the substantially reduced blood inflow into tumor tissue, which is confirmed by reduced rBF in those groups. Overall tumor perfusion was also significantly reduced, as indicated by reduced WiAUC and PE. At the same time, supplementation with a high amount of vitamin D 3 (5000 IU) caused a significant increase in the blood inflow rate (WiR), but a significantly lower WiAUC compared with the 1000 IU group. The high level of WiAUC in the 100 IU group, along with the long RT, indicates the slower blood supply in tumors in the mice fed a vitamin-D 3 -deficient (100 IU) diet, but equally good supply as observed in the 1000 IU reference group. Thus, it can be concluded that in the 4T1 model, vitamin D 3 supply modulates the kinetic, but not overall, tumor perfusion, whereas calcitriol addition reduces tumor perfusion. This finding is supported by WiPI, which was significantly reduced only in calcitriol-receiving groups (Table 4). Interestingly, a similar analysis conducted with E0771-bearing mice showed distinct observations. Vitamin D 3 deficiency (100 IU group) increased the WiR (along with PE and WiAUC), which resulted in an increase in WiPI, indicating high blood perfusion (significantly higher than in the 1000 IU and 5000 IU groups). Surprisingly, calcitriol addition further increased the WiR when combined with 100 IU vitamin D 3 (with negligible influence when combined with 1000 IU). At the same time, a significant reduction in RT along with a substantial decrease in PE, rBV, and WiAUC, was observed, resulting in the WiPI level comparable with other groups, indicating no changes in tumor perfusion. Thus, it can be concluded that in the E0771 model, vitamin D 3 deficiency increases tumor perfusion, and calcitriol addition counteracts it (Table 4).
In mice bearing 67NR tumors, no statistically significant changes in blood flow were observed (Supplementary Figure S4). Altogether, both vitamin D 3 supply and calcitriol addition affect tumor perfusion differently, depending on cancer cell lines used in the study.

COL1A1 and α-SMA Expression in Tumor Tissue
COL1A1 staining of tumor cells did not change significantly in 4T1 and E0771 tumors, but in 67NR tumors, the deficiency diet led to a significant increase in COL1A1 expression as compared to the control and 100 IU+cal groups (Supplementary Figure S6A-D). The expression of α-SMA did not change significantly in all tumor models (Supplementary Figure S6E-H).

Plasma and Tumor Tissue Expression of Selected Cytokines and Growth Factors
The plasma concentrations of OPN were diminished in the 100 IU vitamin-D 3 -deficient group in 4T1-bearing mice and in 100 IU+cal in 67NR tumor-bearing mice, but in healthy BALB/c and C57BL6 mice, the concentration was the highest in the 100 IU+cal group (Table 5). A higher OPN plasma concentration was observed in 4T1 tumor-bearing mice as compared to healthy BALB/c mice and mice bearing 67NR or E0771 tumors. The OPN plasma concentration was also significantly higher in E0771 tumor-bearing mice as compared to C57BL/6 healthy mice (Supplementary Figure S7A). ELISA tests performed in tumor tissue lysates showed that 4T1 tumors were characterized by higher OPN expression as compared to 67NR and E0771, irrespective of the treatment group. No significant changes in OPN expression were observed between the treatment groups with the exception of E0771, where OPN expression in tumor tissue was diminished in the 5000 IU group (Supplementary Figure S8A). The expression of TGF-β was diminished in the plasma of mice bearing 4T1 tumors from 1000 IU+cal, 5000 IU, and 100 IU groups. A similar tendency of TGF-β expression was observed in the plasma of 67NR-bearing mice from the 100 IU group, whereas an opposite tendency was observed in the 1000 IU+cal and 5000 IU groups (Table 5). In C57BL/6 mice bearing E0771 tumors, as well as in healthy mice, the plasma concentration of TGF-β did not change (Table 5). Furthermore, the TGF-β plasma concentration did not differ between tumor models or healthy mice fed with vitamin-D 3 -sufficient diet (1000 IU), whereas in mice bearing E0771 tumors, the expression of TGF-β was significantly higher as compared to 4T1 tumor-bearing mice. Among the 5000 IU groups of BALB/c mice, the concentration of TGF-β was the lowest in 4T1 tumor-bearing mice. However, in deficiency groups (with and without calcitriol), the TGF-β plasma concentration was higher in healthy BALB/c mice as compared to 4T1 tumor-or 67NR tumor-bearing mice (Supplementary Figure S7B). The TGF-β concentration measured in tumor tissue was significantly higher in 67NR tumors from the 100 IU deficiency group as compared to all supplemented groups. A similar tendency was observed in calcitriol-treated 4T1 tumors as compared to the 100 IU-treated group, but the difference was not statistically significant. 4T1 and E0771 tumors showed a similar tumor tissue concentration of TGF-β, but significantly higher than that in 67NR tumors (Supplementary Figure S8B). The plasma concentration of CCL2 in 4T1 tumor-bearing mice was diminished in the 100 IU+cal group as compared to 100 IU and 1000 IU+cal groups, while no change was observed in 67NR tumor-bearing mice and an increased concentration in healthy BALB/c mice in calcitriol-treated groups ( Table 5). In C57BL/6 healthy mice, the plasma concentration of CCL2 was diminished in the deficiency group as compared to the control group, and a similar tendency was preserved in E0771 tumor-bearing mice (p = 0.0779) ( Table 5). The highest plasma concentration of CCL2 was observed in C57BL/6 mice bearing E0771 tumors irrespective of the treatment (as compared to C57BL/6 healthy mice and all BALB/c mice). Among BALB/c mice, the highest CCL2 concentrations were observed in 67NR tumor-bearing (1000 IU and 1000 IU+cal) and healthy (100 IU+cal) mice (Supplementary Figure S7C). An increased plasma concentration of CCL2 was noticed in the tumor tissues of 67NR tumors in the 100 IU group, and in 4T1 tumors from the same group, a tendency to decrease was observed. In E0771 tumors, a significantly decreased CCL2 plasma concentration was observed in the 100 IU+cal group. The CCL2 plasma concentration was higher in E0771 tumors as compared to 4T1 and 67NR (Supplementary Figure S8C). The plasma concentration of VEGF and IL-6 was also analyzed. The VEGF plasma concentration was decreased in 67NR mice in all treatment groups and was almost undetectable with ELISA in BALB/c healthy and C57BL/6 healthy and tumor-bearing mice (Supplementary Figure S7D,E). The highest tumor tissue concentration of VEGF was observed in E0771-bearing mice (Supplementary Figure S8D). In turn, the tumor tissue concentration of VEGF was not changed significantly in 4T1 or E0771 tumor-bearing mice, while it was found to be increased in 67NR tumors in the 100 IU+cal group as compared to the 1000 IU+cal group. In the 5000 IU group, the VEGF plasma concentration was higher in 4T1 than in 67NR tumors (Supplementary Figure S7D). The IL-6 plasma concentration was changed significantly upon treatments only in healthy BALB/c mice; a significant increase in the IL-6 plasma concentration was observed as compared to the control in the 100 IU+cal group (Supplementary Figure S7F). In mice bearing 4T1 and 67NR tumors, the IL-6 plasma concentration was higher as compared to healthy BALB/c mice in control groups; in the 1000 IU+cal, 100 IU, and 100 IU+cal groups, it was higher only in 4T1 tumors. In the 100 IU+cal group, the IL-6 plasma concentration was higher in E0771 tumor-bearing mice as compared to C57BL/6 healthy mice (Supplementary Figure S7G). SDF-1, the tumor concentration of which was the highest in 4T1 tissue (except for the 100 IU group), did not change significantly with respect to treatments applied (Supplementary Figure S8F). The tumor tissue concentration of FGF23 was significantly increased in 100 IU+cal 4T1 tissue; in 67NR tumors, it was increased in the 5000 IU and 100 IU groups as compared to the control 1000 IU group; and it was diminished in the 100 IU+cal as compared to the 100 IU group. In E0771 tumors, the FGF23 concentration was the lowest in the 100 IU+cal group ( Table 5). The FGF23 tumor tissue concentration was higher in E0771 tumors as compared to 4T1 and 67NR and was also higher in 4T1 tissue as compared to 67NR in the 5000 IU and 100 IU+cal groups (Supplementary Figure S8G). Plasma concentration of osteopontin (OPN), transforming growth factor β (TGF-β), and C-C motif chemokine ligand 2 (CCL2) in BALB/c and C57BL/6 mice. Tumor tissue concentration of fibroblast growth factor 23 (FGF23). Mice were fed with diets containing various amounts of vitamin D 3 for 6 weeks. Next, on the day assigned as day 0, tumor cells were implanted orthotopically. Diets were continued for the next 7 days, and then calcitriol gavage was started and continued thrice a week till day 23 (for C57BL/6 mice) or day 28 (for BALB/c mice). Calcitriol was administered by gavage to mice on control (1000 IU) and deficient (100 IU) diets at a dose of 1 µg/kg. During autopsy, blood and tumors were harvested. N = 5-7. Statistical analysis: Kruskal-Wallis test followed by Dunn's test for multiple comparisons. Data are shown as mean ± SD. a ≤ 0.05 as compared to 1000 IU, b ≤ 0.05 as compared to 1000 IU+cal, c ≤ 0.05 as compared to 5000 IU, d ≤ 0.05 as compared to 100 IU.

Effect of TGF-β and Cancer Cell CM on Lung Fibroblasts Harvested from Healthy Mice Fed with Different Vitamin D 3 Diets
Lung fibroblasts obtained from healthy mice fed with different vitamin D 3 diets and treated or not treated with calcitriol were incubated with TGF-β and CM from 4T1, 67NR, and E0771 cell cultures. Then, the fibroblasts were subjected to further analysis. A comparison of the phenotype analysis of untreated cells is presented in Supplementary  Figure S1. 3.5.1. Flow Cytometry Analyses for α-SMA, Podoplanin, PDGFRβ, and TNC Expression TGF-β, but not CM from 4T1 or 67NR cells, decreased the expression of α-SMA in fibroblasts harvested from the control group of BALB/c mice. The MFI analysis showed a similar tendency in the remaining groups ( Table 6). The effect was more pronounced and statistically significant when considering the percentage of positive staining (Supplementary Figure S9A). Podoplanin expression was downregulated by TGF-β (in all groups of mice) but upregulated by 4T1 CM (in all groups except for the control). On the other hand, it was not influenced by 67NR CM (Table 6). PDGFRβ and TNC were downregulated by both TGF-β and 4T1 CM, but only on fibroblasts obtained from mice fed with the control diet (1000 IU) ( Table 6). 3.0 ± 0.6 2.7 ± 0.5 3.1 ± 0.5 3.0 ± 0.5 100 IU 3.0 ± 0.5 2.5 ± 0.3 2.9 ± 0.5 2.9 ± 0.5 100 IU+cal 3.0 ± 0.6 2.5 ± 0.4 2.9 ± 0.7 3.2 ± 0.8 α-SMA-α-smooth muscle actin; PDGFRβ-platelet-derived growth factor receptor beta; TNC-tenascin C. Data are shown as normalized median fluorescence intensity (MFI) calculated by dividing the MFI of the stained sample by the MFI of the isotype control. Healthy mice were fed with diets containing various amounts of vitamin D for 7 weeks. Then, calcitriol gavage was started and continued thrice a week till day 28. Calcitriol was administered to mice receiving control (1000 IU) or deficient (100 IU) diets at a dose of 1 µg/kg. During autopsy, lungs were harvested and NFs were isolated. Flow cytometry analyses were performed after the culture of NFs with TGF-β, 4T1, and 67NR culture media to determine the expression of smooth muscle α-actin (α-SMA), podoplanin, platelet-derived growth factor receptor β (PDGFRβ), and tenascin C (TNC). Representative dot-plots are shown in Supplementary Figure S1. N = 6-7. Statistical analysis: one-way analysis of variance test followed by Sidak's test for multiple comparisons. Data are shown as mean ± SD a ≤ 0.05 as compared to untreated control in the same mice group (1000 IU, 1000 IU+cal, 5000 IU, 100 IU, or 100 IU+cal), x ≤ 0.05 as compared to the same treatment in the 1000 IU group.
TGF-β or E0771 CM did not significantly affect the expression of α-SMA and podoplanin on lung fibroblasts harvested from C57BL/6 mice (Supplementary Figures S9B and S10A,B). In addition, TGF-β or E0771 CM did not influence the expression of TNC (Supplementary Figures S9B and S10C), with the exception of fibroblasts from mice fed with the 5000 IU vitamin D diet, where the percentage of TNC-positive cells significantly increased after TGF-β stimulation (Supplementary Figure S9B). In the case of PDGFRβ (Supplementary Figures S9B and S10C), especially when considering the percentage of positive stained cells in all groups (Supplementary Figure S9B), TGF-β stimulation led to an increased expression of this receptor. E0771 CM had no significant influence on PDGFRβ expression ( Supplementary Figures S9B and S10C).

Fluorescence Microscopy Analysis of VDR and FAP Expression in Lung Fibroblasts
The expression of VDR and FAP was also assessed in cultured lung NFs from healthy BALB/c (Figure 2A-L) and C57BL/6 ( Figure 3A-J) mice. Lung NFs from BALB/c mice administered with calcitriol showed higher expression of VDR than that of appropriate controls (1000 IU and 100 IU), whereas lung NFs from mice in both vitamin D deficiency and supplementation groups showed lower VDR expression than that of control (1000 IU) mice (Figure 2A). TGF-β stimulation of lung fibroblasts did not change the effect of calcitriol but increased the expression of VDR in fibroblasts from the 5000 IU group ( Figure 2B,K). Stimulation with 4T1 CM led to a similar profile of VDR expression as stimulation with TGFβ, but the effect was weaker ( Figure 2C,K). Stimulation with 67NR CM did not significantly alter the relationships observed in nonstimulated cells (Figure 2A vs. Figure 2D). However, an analysis of the effects of stimulations used in the control group (1000 IU) revealed that none of them affected VDR expression ( Figure 2I). In lung fibroblasts obtained from mice on a control diet and administered with calcitriol (1000 IU+cal), the expression of VDR significantly decreased upon all stimulations used. However, in fibroblasts from the 5000 IU, 100 IU, and 100 IU+cal groups, stimulation with TGF-β and 4T1 CM caused a significant increase in VDR expression. 67NR CM had no effect on VDR expression in these groups ( Figure 2I).
FAP expression was the highest in lung fibroblasts from mice on a deficient diet administered with calcitriol (100 IU+cal), and in general, all treatments used resulted in increased expression of FAP in lung fibroblasts derived from BALB/c mice ( Figure 2E). Stimulation with TGF-β led to a reduction in or elimination of differences observed in nontreated cells derived from various groups of BALB/c mice. At the same time, TGF-β led to a significant decrease in FAP staining in cells derived from the 1000 IU+cal group as compared to the control (1000 IU) group ( Figure 2F). A similar effect was observed in cells incubated with 4T1 CM: a significantly lower expression of FAP in both groups treated with calcitriol and the deficient group ( Figure 2G,L). In the case of 67NR CM stimulation, lung fibroblasts from mice administered with calcitriol and vitamin-D-deficient mice showed higher FAP staining than that of control mice. Fibroblasts from mice on the vitamin-D-supplemented diet showed lower FAP staining than control mice upon 67NR CM stimulation ( Figure 2H). While comparing the effects of stimulation within each group of mice, it was observed that TGF-β stimulation increased FAP expression in the control (mice fed with the 1000 IU vitamin D diet) and deficiency groups. However, in the deficiency group treated with calcitriol, the opposite effect was observed. Stimulation with 4T1 CM led to increased FAP fluorescence in fibroblasts from control mice and decreased FAP fluorescence in fibroblasts from the 1000 IU+cal, 100 IU, and 100 IU+cal groups. 67NR CM decreased FAP staining in the 5000 IU and 100 IU+cal groups and increased it in the deficient group ( Figure 2J,L).
For control mice, the highest VDR expression was observed in cultured ex vivo lung NFs from healthy C57BL/6 mice. Calcitriol administration to the vitamin-D-deficient group led to a significant increase in VDR expression ( Figure 3A). Stimulation with TGF-β or E0771 CM did not significantly alter these relationships ( Figure 3B,C,I). However, both stimulation with TGF-β and E0771 CM of lung fibroblasts from control mice decreased VDR expression. VDR expression in lung fibroblasts from mice fed with a control diet and treated with calcitriol did not change upon stimulations, but it was increased by TGF-β in the 5000 IU, 100 IU, and 100 IU+cal groups ( Figure 3G). FAP expression was the highest in lung fibroblasts from mice fed with the 5000 IU diet and the lowest in the deficiency group administered with calcitriol (100 IU+cal) ( Figure 3D,J). Stimulation with TGF-β or E0771 CM exacerbated the effect of calcitriol in lowering the FAP expression in both groups ( Figure 3E,F). However, the effect of TGF-β stimulation on FAP expression was opposite to that of E0771 CM stimulation, causing an increase and a decrease in FAP expression, respectively. Therefore, the lowest FAP expression was noticed in fibroblasts harvested from mice fed a deficient diet and administered with calcitriol (100 IU+cal) ( Figure 3H,J). NFs were treated ex vivo with transforming growth factor-β (TGF-β) or culture media from the culture of 4T1 and 67NR cells. Cells were reviewed and photographed using an Olympus IX81 fluorescence microscope. Fluorescence was determined from images using ImageJ according to the protocol of Luke Hammond (QBI, The University of Queensland, Australia; accessed on 18 October 2021; https://theolb.readthedocs.io/en/latest/imaging/measuring-cell-fluorescence-using-imagej.html). The areas of interest (whole cell for FAP staining and cell nucleus for VDR staining) and areas of background were measured, and Corrected Total Cell Fluorescence (CTCF) was calculated using the following formula: CTCF = integrated density-(area of selected cell × mean fluorescence of background readings). Data are shown as (A-H) mean with individual data; (I-J) mean ± standard error of mean. Statistical analysis: Kruskal-Wallis test followed by Dunn's test. * p < 0.05; ** p < 0.01; *** p < 0.001; **** p < 0.0001.
For control mice, the highest VDR expression was observed in cultured ex vivo lung NFs from healthy C57BL/6 mice. Calcitriol administration to the vitamin-D-deficient group led to a significant increase in VDR expression ( Figure 3A). Stimulation with TGFβ or E0771 CM did not significantly alter these relationships ( Figure 3B,C and I). However, NFs were treated ex vivo with transforming growth factor-β (TGF-β) or culture media from the culture of 4T1 and 67NR cells. Cells were reviewed and photographed using an Olympus IX81 fluorescence microscope. Fluorescence was determined from images using ImageJ according to the protocol of Luke Hammond (QBI, The University of Queensland, Australia; accessed on 18 October 2021; https: //theolb.readthedocs.io/en/latest/imaging/measuring-cell-fluorescence-using-imagej.html). The areas of interest (whole cell for FAP staining and cell nucleus for VDR staining) and areas of background were measured, and Corrected Total Cell Fluorescence (CTCF) was calculated using the following formula: CTCF = integrated density-(area of selected cell × mean fluorescence of background readings). Data are shown as (A-H) mean with individual data; (I-J) mean ± standard error of mean. Statistical analysis: Kruskal-Wallis test followed by Dunn's test. * p < 0.05; ** p < 0.01; *** p < 0.001; **** p < 0.0001. ciency group administered with calcitriol (100 IU+cal) ( Figure 3D and J). Stimulation with TGF-β or E0771 CM exacerbated the effect of calcitriol in lowering the FAP expression in both groups ( Figure 3E,F). However, the effect of TGF-β stimulation on FAP expression was opposite to that of E0771 CM stimulation, causing an increase and a decrease in FAP expression, respectively. Therefore, the lowest FAP expression was noticed in fibroblasts harvested from mice fed a deficient diet and administered with calcitriol (100 IU+cal) (Figure 3H and J).  NFs were treated ex vivo with transforming growth factor-β (TGF-β) or culture media from the culture of E0771 cells. Cells were reviewed and photographed using an Olympus IX81 fluorescence microscope. Fluorescence was determined from images using ImageJ according to the protocol of Luke Hammond (QBI, The University of Queensland, Australia; accessed on 18 October 2021; https://theolb.readthedocs.io/en/latest/imaging/measuring-cell-fluorescence-using-imagej. htmL). The areas of interest (whole cell for FAP staining and cell nucleus for VDR staining) and areas of background were measured and Corrected Total Cell Fluorescence (CTCF) was calculated using the following formula: CTCF = integrated density − (area of selected cell × mean fluorescence of background readings). Data are shown as (A-H) mean with individual data; (I-J) mean ± standard error of mean. Statistical analysis: Kruskal-Wallis test followed by Dunn's test. * p < 0.05; ** p < 0.01; *** p < 0.001; **** p < 0.0001.

Acta2, Mmp9, Spp1, and Vdr Expression in Cultured NFs Isolated from Lungs of Healthy BALB/c or C57BL/6 Mice Fed with Control and Vitamin-D 3 -Deficient Diet and/or Treated with Calcitriol
The expression of selected genes was measured in the control (1000 IU) and deficiency groups (100 IU and 100 IU+cal), in which the highest effects on fibroblasts were noticed ( Figure 4). In NFs from BALB/c mice, the expression of Acta2 (encoding α-SMA) increased upon stimulation with TGF-β, but this effect was not significant in the 100 IU group ( Figure 4A). In C57BL/6 NFs, TGF-β stimulation increased Acta2 expression in all groups; however, this effect was the lowest in the 100 IU+cal group. E0771 CM also stimulated Acta2 expression, but only in NFs from mice treated with 100 IU+cal ( Figure 4B). In NFs from the BALB/c deficiency group, the expression of Mmp9 was significantly stimulated by TGF-β and 4T1 CM. In the 100 IU+cal group, this effect was preserved only after stimulation with 4T1 CM ( Figure 4A). 4T1 CM also stimulated the expression of Spp1 (encoding OPN), but only in NFs from control (1000 IU) BALB/c mice ( Figure 4A). Interestingly, in NFs from C57BL/6 control mice (1000 IU), decreased Spp1 expression was observed after stimulation with E0771 CM and TGF-β. In NFs from deficiency group stimulation with E0771 CM led to a significant increase in Spp1 expression in comparison to the 1000 IU group ( Figure 4B). TGF-β stimulation decreased the expression of Vdr in NFs from all groups of BALB/c and C57BL/6 mice ( Figure 4A,B). In the 1000 IU groups, 4T1 CM increased, while E0771 decreased, the Vdr expression. However, in deficiency groups, an increase in Vdr expression was observed only after stimulation with 67NR CM. Moreover, Vdr expression was the highest after 67NR stimulation in the 100 IU+cal group, similar to that observed after E0771 CM stimulation ( Figure 4A,B).

Discussion
To further understand the effect of varying vitamin D3 status on the growth and metastasis of breast cancer, we implanted two metastatic mouse breast cancer cells (4T1 and E0771) and nonmetastatic cells (67NR) into mice fed with vitamin-D3-normal, vitamin-D3deficient, and vitamin-D3-supplemented diets. Implantation was performed when differ- . Gene expression in cultured NFs isolated from lungs of healthy BALB/c or C57BL/6 mice fed a control and vitamin-D 3 -deficient diet and/or treated with calcitriol. NFs were treated ex vivo with transforming growth factor-β (TGF-β) or culture media from the culture of 4T1, 67NR (A), or E0771 (B) cells. Expression of Acta2 (α-smooth muscle actin), Mmp9 (matrix metalloproteinase 9), Spp1 (osteopontin) and Vdr (vitamin D receptor) was calculated according to the comparative ∆∆Ct method with Gadph and Rps27a as endogenous controls and normalized to each untreated control using QuantStudio™ Real-Time PCR Software and ExpressionSuite Software. Data are shown as mean ± standard error of mean. Statistical analysis: Kruskal-Wallis test followed by Dunn's test. * p < 0.05; ** p < 0.01; *** p < 0.001; **** p < 0.0001.

Discussion
To further understand the effect of varying vitamin D 3 status on the growth and metastasis of breast cancer, we implanted two metastatic mouse breast cancer cells (4T1 and E0771) and nonmetastatic cells (67NR) into mice fed with vitamin-D 3 -normal, vitamin-D 3 -deficient, and vitamin-D 3 -supplemented diets. Implantation was performed when differences in the plasma 25(OH)D 3 concentrations statistically differed between normal and experimental diets, which was achieved after 6 weeks of feeding [3]. On day 7 after tumor transplantation, calcitriol gavage (1 µg/kg/day) was started in mice fed with normal and deficient diets and the same diets were continued [3]. As discussed earlier, the effect of the above treatments on tumor tissue blood flow was the highest in deficiency groups in both 4T1 and E0771 tumors: vitamin D 3 deficiency significantly increased selected bloodflow parameters, whereas a deficient diet followed by calcitriol gavage led to a significant decrease in most of them. However, in 4T1 tumors (not in E0771), some parameters describing the blood flow were also decreased, but other parameters were increased in mice fed with a normal diet and administered with calcitriol. The modulation of tumor vascularization is not only important in the context of metastasis but may also have an impact on the response to treatment. Decreased blood flow in the tumor may cause a decrease in the response to chemotherapy [40]. A previous study showed that calcitriol treatment led to decreased vascularization of MCF-7 human breast cancer tumors with overexpression of VEGF [41]. However, other studies have shown that calcitriol stimulated in vitro expression of VEGF and downregulated the expression of thrombospondin 1. In T47D human breast cancer, calcitriol in vivo inhibited tumor growth, but had no effect on vascular density. Moreover, calcitriol showed a stimulating effect on the plasma VEGF concentration in the tested mice [42]. In the present study, no effect of the used treatment on VEGF expression was observed in 4T1 or E0771 tumor tissue or in mice plasma; however, in 67NR tumor-bearing mice, in which no effect of treatments on blood flow was noted, the VEGF plasma concentration was decreased irrespective of the treatment as compared to control mice. Therefore, it can be concluded that the effect of vitamin D 3 on tumor vascularization or angiogenic factors may depend on the type of tumor. Our studies have indicated that the basal expression of molecules, which can be implicated in the angiogenesis process, differ significantly between mouse mammary gland tumors and the plasma of tumor-bearing mice. For instance, in this study, the concentration of CCL2 was significantly higher in the tumor-bearing mice and in the plasma of E0771-bearing mice; however, OPN was high in concentration in 4T1 plasma and tumor tissue. Both cytokines can modify the tumor microenvironment, either by directly influencing angiogenesis [43,44] or via recruiting macrophages (CCL2) [45] and activating fibroblasts (OPN) [24]. This effect may also depend on the dose and schedule of treatment with calcitriol or vitamin D 3 (as can be seen from different effects of calcitriol gavage on, for example, blood flow in tumor tissue in mice fed with vitamin-D 3 -normal or vitamin-D 3 -deficient, or vitamin-D 3supplemented diet) as well as the observation time. In our previous studies, a decrease in the 4T1 tumor tissue VEGF concentration was found on day 14 of tumor progression in mice injected subcutaneously (0.5 µg/kg/day) with calcitriol or its analogs, with increased PE and TTP observed in young mice on day 24 [2]. However, the effect on 4T1 cancer metastasis was the same, namely in mice from which tissue samples were collected in the present study [3], as well as in our previous work [2]: in all groups of mice in which calcitriol was administered or vitamin D 3 content in the diet was increased, an enhanced metastatic process was observed [2,3].
As evidenced from previous studies, CAFs play an important role in invasion and metastasis as well as in angiogenesis of tumors [20,46]. Tumor cells can stimulate stromal fibroblasts to differentiate into CAFs, and the crosstalk between cancer cells and stromal fibroblasts leads to tumor progression [24]. In the present study, the phenotype of CAFs as well as lung NFs from tumor-bearing and healthy mice was analyzed. Upon various vitamin D 3 or calcitriol treatments, significant changes in the phenotype of CAFs were observed only in E0771 tumors. Interestingly, in lung NFs from healthy C57BL/6 mice, calcitriol did not influence the expression of studied molecules, but in NFs from BALB/c mice, α-SMA expression was decreased in all treatment groups. However, in tumor-bearing mice, the effect of various treatments on NFs was dependent on the tumors implanted. No significant effects were observed on NFs from 67NR nonmetastatic and E0771 metastatic tumors. Whereas in highly metastatic 4T1 tumor-bearing mice, signs of fibroblast activation were observed in both groups of mice administered with calcitriol but were more pronounced in mice on a deficient diet (100 IU+cal). NFs from this group of mice showed increased expression of α-SMA and podoplanin, as well as TNC, which suggests increased activation of NFs [47]. Interestingly, the expression of these proteins was significantly higher in the 100 IU+cal group also as compared to the supplemented group. In this tumor model, metastases were located in the lungs and were significantly increased in both groups treated with calcitriol as well as in the 5000 IU group but by a lower number (count for macrometastases presented earlier [3] as mean number ± standard deviation is as follows-control 40 ± 4; 1000 IU+cal 67 ± 27; 5000 IU 52 ± 11; 100 IU 44.5 ± 15; 100 IU+cal 77 ± 27) [3]. This proves the impact of the mammary gland tumor type on the biological effects of vitamin D 3 . For example, the expression of OPN and SDF-1, which are some of the important molecules secreted from cancer cells stimulating the differentiation of fibroblasts to myofibroblasts [24,48], was the highest in plasma and/or tumor tissue of 4T1 as compared to 67NR and E0771 tumor-bearing mice. In the present study, no significant effects of vitamin D 3 on plasma OPN were observed, but previously, in later stages of 4T1 tumor progression (days 28 and 33), calcitriol injections led to an increased plasma concentration of OPN [2]. Moreover, increased features of NF activation noted in 4T1 tumor-bearing mice, which were the most visible in the 100 IU+cal group, may also be an effect of the increased tumor tissue concentration of IL-6, a cytokine that may activate fibroblasts through the induction of STAT3 phosphorylation [49]. As a consequence of myofibroblast activation, the increased expression of α-SMA on NFs may result in macrophage recruitment and M2 polarization as well as myeloid-derived suppressor cells (MDSC) and regulatory T lymphocytes (Tregs) recruitment and their differentiation with Th2 polarization [50,51]. Similarly, the expression of podoplanin on fibroblasts may be responsible for M2 polarization and immunosuppression of the tumor microenvironment [52]. In addition, TNC-expressing fibroblasts modulate the recruitment of monocytes/macrophages [53]. Our previous study has shown that treatment of 4T1 tumor-bearing young mice with calcitriol or its analogs (mice on a diet containing a normal concentration of vitamin D 3 ) led to an increase in the percentage of monocytes in the blood [54] and an increase in the ratio of inflammatory (Ly6C high CXCR1 low CCR2 + ) spleen monocytes to monocytes with repair and anti-inflammatory properties (Ly6C low CXCR1 high ) [27]. Moreover, previous ex vivo studies confirmed that cultured media from 4T1 cells enhanced the M2 polarization stimulated by calcitriol; moreover, M2 macrophages differentiated in the presence of calcitriol stimulated the migration of 4T1 and 67NR cells in vitro [55]. Treg prevalence with Th2 polarization was also observed in calcitriol-treated mice bearing 4T1 tumors [4,54]. Together, these factors lead to the modulation of blood flow and/or increased metastasis as an effect of calcitriol treatment of mice bearing 4T1 tumors observed in our studies [3].
Further analysis of the effects of calcitriol on fibroblast activation was performed on lung NFs from healthy BALB/c and C57Bl/6 mice treated with various diets and/or calcitriol. It is worth emphasizing that NFs from these two strains of mice react differently to various vitamin D 3 or calcitriol treatments. The differences between these two strains of mice, for example, in the concentration of vitamin D metabolites were already reported. For example, higher plasma concentrations of 25(OH)D 3 , 24,25(OH) 2 D 3 , and 3-epi-25(OH)D 3 were observed in C57BL/6 than in BALB/c mice in the control (1000 IU), 1000 IU+cal, and 5000 IU groups [3]. In the present study, decreased α-SMA expression on NFs from BALB/c mice was observed in all treatment groups, but no effect was observed in C57BL/6 mice. In both strains of healthy mice, vitamin D 3 deficiency combined with calcitriol gavage (100 IU+cal) led to an increased plasma concentration of OPN, but CCL2 and IL-6 were increased in these groups only in BALB/c mice; in C57BL/6 mice, the opposite tendency was observed. NFs growing in such microenvironments modified by varying vitamin D 3 treatments were stimulated ex vivo with TGF-β or cancer cell CM. Because the most significant changes in fibroblasts were observed in deficiency groups, the gene expression in NFs of these mice was also analyzed. However, these ex vivo studies did not fully confirm the observations of lung NF activation in 4T1 tumor-bearing mice. In the ex vivo experimental conditions of the present study, there were limited components of the full microenvironment in which CAFs interact with tumor cells and with TAMs, and their interactions are important for their final phenotype [56]. In NFs from deficiency groups, 4T1 CM stimulated Mmp9 expression as well. Circulating MMP9 was previously described as inversely correlated with the vitamin D plasma concentration, with the highest concentrations in vitamin-D-deficient healthy adults [57]. Furthermore, E0771 CM stimulated the expression of Spp1. Our previous study demonstrated that calcitriol stimulated OPN secretion in vitro in a mouse BALB/3T3 NF cell line [27], and others reported a similar effect of calcitriol on other mouse and human fibroblast cell lines [32]. This effect is dependent on the interaction of VDR with the Spp1 gene [58]. Moreover, the authors of these studies reported that CAFs and NFs expressed VDR, and the expression of selected genes such as CD82 and S100A4 correlated with stromal VDR expression and clinical outcome in patients with colorectal patients [32]. VDR expression was also analyzed in stimulated lung NFs in our study. In E0771-stimulated NFs, the expression of Vdr mRNA was the highest in the 100 IU+cal group, in which the expression of Spp1 was also the highest, but on the protein level, the expression of VDR was lower in this group of mice as compared to the control. Thus, in this study, VDR expression did not correspond to the observed effects of vitamin D deficiency/supplementation, also when measured in whole tumor tissue, as we previously described [3].
The differences between two metastatic mouse mammary gland cancer models observed in response to vitamin D 3 can be clearly visualized. Of particular note are the differences in pulmonary NFs responses to different modes of vitamin D 3 delivery in the 4T1 tumor model. Mice fed with a vitamin-D 3 -supplemented diet (5000 IU) showed decreased expression of proteins responsible for protumoral activity of myofibroblasts (α-SMA, PDGFRβ, and TNC). This effect was accompanied by high plasma concentrations of vitamin D metabolites 25(OH)D 3 and 24,25(OH) 2 D 3 at the beginning of tumor progression, but significantly decreased during continued supplementation in the late stages of tumor growth. High concentrations of 3-epi-25(OH)D 3 were observed during the whole period of observation [3]. In the E0771 tumor model, no change in lung NF activity was observed, but CAFs in mice fed with a vitamin-D-supplemented diet (5000 IU) were activated as podoplanin, TNC, and PDGFRβ increased. In these mice, however, all vitamin D metabolites measured were high during the whole experiment [3]. Thus, it can be noticed that vitamin D 3 supplementation may lead to decreased or increased fibroblast activation, which is dependent on the tumor model/microenvironment and may also be dependent on the mice strain studied [3]. Interestingly, vitamin D deficiency (which in general does not lead to increased metastatic process, like vitamin D 3 supplementation or calcitriol administration in the 4T1 model, but shows a tendency for an increased number of micrometastases [3]) also led to decreased NFs activation in the 4T1 model (decreased podoplanin and TNC) and to some signs of increased activation of CAFs in the E0771 model (increased podoplanin). In these deficient mice, a low concentration of vitamin D metabolites was observed during the whole experiment [3]. The most interesting observation is the comparison of calcitriol administration in mice fed with a normal vs. deficient diet, especially in the 4T1 model. Only α-SMA induction on NFs was observed in the former, but in the latter, podoplanin and TNC also increased significantly. Mice fed with a deficient diet and treated with calcitriol showed the lowest concentrations of all vitamin D metabolites, also as compared to mice on a deficient diet [3].
In summary, the impact of vitamin D 3 supplementation or calcitriol administration on blood flow or CAFs and lung NFs is not visible in less invasive tumor models, such as 67NR, despite induced changes in plasma and/or tumor tissue concentrations of OPN, CCL2, TGF-β, VEGF, and FGF23. Moreover, differences between invasive tumor models were observed in response to the above treatments. In E0771-tumor-bearing mice, the influence on the metastatic process was not visible, but a tendency of tumor growth inhibition was observed in the 5000 IU diet as well as in 100 IU+cal groups [3]. In the latter, a decrease in the blood flow was also observed, which was accompanied by a decrease in FGF23 and CCL2 in tumor tissue, but an increase in TNC on CAFs. Interestingly, similar effects on blood flow (decrease) were also observed in 4T1 tumors from the 100 IU+cal group, but it did not affect the CAF phenotype. At the same time, increasing α-SMA, podoplanin, and TNC on lung NFs and increasing metastasis.

Conclusions
Divergent effects of vitamin D 3 supplementation or deficiency in healthy mice lead to the formation of various microenvironments in the body, which depends on the mouse strain. Tumors develop in such microenvironments and themselves modify the microenvironments, producing higher concentrations of OPN, SDF-1 (4T1), TGF-β (4T1 and E0771), CCL2, VEGF, FGF23 (E0771), and IL-6 (67NR), which influence the response to vitamin D 3 supplementation/deficiency and calcitriol administration. Based on the results presented here, it can be concluded that the safest approach may be to constantly maintain the appropriate concentration of vitamin D because additional supplementation, especially in the case of primary deficiency, may lead to undesirable effects.
Supplementary Materials: The following supporting information can be downloaded at: https: //www.mdpi.com/article/10.3390/cancers14194585/s1, Figure S1. Flow cytometry analysis of lung NFs from healthy mice fed with different amounts of vitamin D in diets and gavaged with calcitriol. Figure S2. Flow cytometry analysis of lung NFs from mice bearing 4T1, 67NR, and E0771 tumors. Figure S3. Flow cytometry analysis of CAFs from mice bearing 4T1, 67NR, and E0771 tumors. Figure  S4. Blood flow in tumor tissue from 67NR-tumor-bearing mice. Figure S5. Blood flow in tumor tissue from 4T1-and E0771-tumor-bearing mice. Figure S6. Tumor tissue expression of COL1A1 and α-SMA. Figure S7. Plasma level of cytokines. Figure S8. Tumor tissue level of cytokines. Figure  S9. Phenotypic characteristics of NFs from healthy BALB/C and C57BL/6 mice treated ex vivo with TGF-β or CM from the culture of 4T1, 67NR (BALB/c), and E0771 (C57BL/6) cells. Figure S10. Phenotypic characteristics of NFs from healthy C57BL/6 mice treated ex vivo with TGF-β or CM from the culture of E0771 cells.

Conflicts of Interest:
The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.