- Open Access
1,25-Dihydroxyvitamin D3 modulates calcium transport in goat mammary epithelial cells in a dose- and energy-dependent manner
Journal of Animal Science and Biotechnology volume 7, Article number: 41 (2016)
Calcium is a vital mineral and an indispensable component of milk for ruminants. The regulation of transcellular calcium transport by 1,25-dihydroxyvitamin D3 (1,25-(OH)2D3, the active form of vitamin D) has been confirmed in humans and rodents, and regulators, including vitamin D receptor (VDR), calcium binding protein D9k (calbindin-D9k), plasma membrane Ca2+-ATPase 1b (PMCA1b), PMAC2b and Orai1, are involved in this process. However, it is still unclear whether 1,25-(OH)2D3 could stimulate calcium transport in the ruminant mammary gland. The present trials were conducted to study the effect of 1,25-(OH)2D3 supplementation and energy availability on the expression of genes and proteins related to calcium secretion in goat mammary epithelial cells.
An in vitro culture method for goat secreting mammary epithelial cells was successfully established. The cells were treated with different doses of 1,25-(OH)2D3 (0, 0.1, 1.0, 10.0 and 100.0 nmol/L) for calcium transport research, followed by a 3-bromopyruvate (3-BrPA, an inhibitor of glucose metabolism) treatment to determine its dependence on glucose availability. Cell proliferation ratios, glucose consumption and enzyme activities were measured with commercial kits, and real-time quantitative polymerase chain reaction (RT-qPCR), and western blots were used to determine the expression of genes and proteins associated with mammary calcium transport in dairy goats, respectively.
1,25-(OH)2D3 promoted cell proliferation and the expression of genes involved in calcium transport in a dose-dependent manner when the concentration did not exceed 10.0 nmol/L. In addition, 100.0 nmol/L 1,25-(OH)2D3 inhibited cell proliferation and the expression of associated genes compared with the 10.0 nmol/L treatment. The inhibition of hexokinase 2 (HK2), a rate-limiting enzyme in glucose metabolism, decreased the expression of PMCA1b and PMCA2b at the mRNA and protein levels as well as the transcription of Orai1, indicating that glucose availability was required for goat mammary calcium transport. The optimal concentration of 1,25-(OH)2D3 that facilitated calcium transport in this study was 10.0 nmol/L.
Supplementation with 1,25-(OH)2D3 influenced cell proliferation and regulated the expression of calcium transport modulators in a dose- and energy-dependent manner, thereby highlighting the role of 1,25-(OH)2D3 as an efficacious regulatory agent that produces calcium-enriched milk in ruminants when a suitable energy status was guaranteed.
As a crucial macro-mineral for animals, calcium has functions in many physiological processes, including skeletal formation, nerve pulse transmission, muscle contraction, blood clotting, stimulus secretion coupling, and is an indispensable component of milk [1–3]. Milk is a naturally calcium-rich fluid produced by animals and humans. Actually, the total calcium concentration in ruminant milk is approximately 30 mmol/L . It was reported that a substantial calcium flux was generated from blood to milk during the lactation period [5–8]. Accordingly, there must be a precise regulatory mechanism involved in the modulation of calcium transport in the mammary glands of dairy animals.
It is not entirely understood how mammary epithelial cells (MECs) extract large quantities of ionized calcium from plasma and produce a calcium-rich secretion, particularly for ruminants. The blood total calcium levels of dairy cows have a narrow range (approximately 2.0 to 2.5 mmol/L) ; thus, the process of calcium transport in the mammary gland occurs against a tremendous concentration gradient. Moreover, VanHouten and Wysolmerski  reported the existence of transcellular calcium transport and summarized this process in human MECs. Consequently, it can be extrapolated that the transcellular process is involved in calcium transport during milk secretion in ruminants.
Calcium-transport proteins, such as calcium binding protein-D9k (calbindin-D9k), plasma membrane Ca2+-ATPase 1b (PMCA1b) and 2b (PMCA2b), have been confirmed as essential elements for transcellular calcium transport [5, 7, 10, 11]. According to recent research, Orai1, a pore subunit of the Ca2+ release-activated Ca2+ (CRAC) channels, is essential for calcium entry into cells and calcium homeostasis [12–14], but no trial has been conducted in mammary epithelial cells from dairy goats. Evidence circumstantiated that 1,25-dihydroxyvitamin D3 (1,25-(OH)2D3), the active form of vitamin D, was the most critical regulator of transcellular calcium transport and body calcium homeostasis [1, 15, 16]. 1,25-(OH)2D3 stimulated mammary calcium transport to elevate the milk calcium content by upregulating calbindin-D9k and PMCA2b in lactating mice; knockout mice were used in this study . Furthermore, 1,25-(OH)2D3 has been reported to facilitate the synthesis of epithelial calcium channels, increase the expression of plasma membrane calcium pumps, and induce the formation of calbindin in humans, rats and other species [18–20]. In addition, Kohler et al.  measured the blood concentrations of 1,25-(OH)2D3 in lactating goats at different altitudes, but the potential regulatory effects of 1,25-(OH)2D3 on mammary calcium transport and milk secretion, such as the expression of key regulators, were not studied. In summary, few research studies called attention to goat mammary calcium transport, and it has not been fully elucidated whether 1,25-(OH)2D3 regulates calcium transport in goat MECs.
Therefore, we hypothesized that 1,25-(OH)2D3 supplementation could modulate the expression of genes involved in calcium transport in goat MECs in a dose-dependent manner. Meanwhile, as an active transport process, calcium transport might be influenced by the cellular energy status.
In the present research, all the procedures and operation were approved by the Animal Welfare Committee of Institute of Animal Nutrition and Feed Science, College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, P.R. China.
In vitro culture of goat mammary epithelial cells
Dulbecco’s Modified Eagle Medium F12 (DMEM/F-12), fetal bovine serum (FBS), epidermal growth factor (EGF) and 0.25 % trypsin were purchased from Life Technologies (Carlsbad, California, USA). Penicillin, streptomycin, insulin and hydrocortisone were obtained from Sigma-Aldrich (Shanghai, China). The other materials used for cell culture were provided by Dr. Xiaofei Wang from the Institute of Animal Nutrition and Feed Science, Northwest A&F University, China.
Three healthy China Guanzhong dairy goats that had been raised in the livestock farm of Northwest A&F University since birth were selected for this study and used during the second parity and at peak lactation (day in milk (DIM) = 60 d). In detail, a 1 cm3 sample of the parenchymal tissue of the mammary gland was collected and placed in sterilized tubes containing ice-cold D-Hanks’ balanced salt solution (D-HBSS; pH = 7.4) after official approval for scientific sampling, and the tubes were immediately and aseptically transported to the laboratory immediately and aseptically. The tissue samples were washed with D-HBSS several times until the washing buffer was transparent, then sheared into 0.5 to 1.0 mm3 cubic fragments with a sterilized surgical scissor, and washed until clean. These fragments were placed in empty 60 mm cell culture dishes (Corning, New York, USA), maintaining an approximate distance of 0.5 cm between pieces, and the dishes were incubated in a cell incubator (Thermo Scientific, Massachusetts, USA) at 37 °C in 5 % CO2 and 95 % air for 30 min. Then, 1 mL of basal medium was added and incubated for 2 h, followed by the addition of another 1 mL of basal medium and incubation for an additional 48 h. The basal media contained 90 % DMEM/F-12 and 10 % FBS, and the concentrations of penicillin, streptomycin, insulin, hydrocortisone and EGF were 100.0 U/mL, 100.0 μg/mL, 5.0 μg/L, 1.0 μg/L and 1.0 μg/L, respectively. The medium was substituted for fresh basal medium every 48 h. When 90-95 % of the dish was occupied by visible cells under an inverted microscope (Nikon, Tokyo, Japan), the cells could be passaged. The cells were digested with 0.25 % trypsin for 5 min and passaged to new dishes. Subsequently, the medium was transferred to separate new culture plates 40 min later and incubated for 48 h to remove the fibroblasts. The adhesion time for fibroblasts (30 to 40 min) was shorter than that of MECs; hence, purified MECs were procured after the last procedure was repeated 5 times. The MECs were previously characterized by Wang et al.  in our college.
Purified MECs passaged to 7–12 generations were used in this study. The cells were seeded in 24-well flat-bottom culture plates (Corning, New York, USA) at a density of 2.0 × 104 cells per well. Afterward, 700 μL of basal medium was added to each well and incubated for 24 h. The medium was removed, the cells were washed with sterilized phosphate-buffered saline (PBS; pH = 7.4) 3 times, and then 700 μL/well of treatment medium containing 1,25-(OH)2D3 (Sigma-Aldrich, Shanghai, China) was added. The final concentrations of 1,25-(OH)2D3 in the medium were 0, 0.1, 1.0, 10.0 and 100.0 nmol/L, respectively. Each treatment was conducted on 6 replicates with 1 replicate per passage to avoid the potential effects of different passages. Culture dishes were incubated under the same conditions described above for 24 h, and then the subsequent steps and analyses were implemented.
A specific inhibitor of hexokinase 2 (HK2), 3-bromopyruvate (3-BrPA; Sigma-Aldrich, Shanghai, China), was added to the medium to investigate the potential effects of the cellular energy status on calcium transport. HK2 phosphorylates glucose to generate glucose-6-phosphate (G6P), the first step in the cellular glucose catabolism, and HK2 inhibition is usually used to study the effect of energy status on metabolic processes . The concentrations of 1,25-(OH)2D3 were 0 or 10.0 nmol/L, and the 3-BrPA concentrations were 0 or 50.0 μmol/L, respectively. The other procedures were consistent with the 1,25-(OH)2D3 treatment.
Cell proliferation measurement
A commercial 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) Kit was obtained from Jiancheng Bioengineering Institute (Nanjing, China) to measure cell proliferation. Briefly, the MECs were seeded in a 96-well plate (2.0 × 104 cells/well; Corning, New York, USA) and were incubated with basal medium (200 μL/well) at 37 °C in 5 % CO2 and 95 % air for 24 h. Subsequently, the basal medium was replaced with treatment medium (150 μL/well) and incubated under standard conditions for 24 h. Then, 1× MTT (50 μL/well) was added and incubated under the same conditions for 4 h. The supernatant was removed carefully and 150 μL dimethyl sulfoxide (DMSO; Amresco, OH 44139, USA) was added to each well, followed by an 8 min mixing process using a Tablet Shaker (Kylin-Bell Lab Instruments Co., Ltd., Jiangsu, China). The absorbance at 570 nm was determined using a Microplate Reader (Power Wave XS2, Bio Tek, USA).
The glucose content in the medium was determined via a Glucose Assay Reagent Kit (Jiancheng, Nanjing, China) based on the glucose oxidase/peroxidase colorimetric method. Medium samples were collected in each well of culture dishes. The reaction reagent (1,000 μL) and liquid sample (10 μL) were mixed in a pure plastic tube, incubated at 37 °C for 15 min, and then the optical density (OD) at 505 nm was read on a Microplate Reader (Power Wave XS2, Bio Tek, USA). The OD of a tube with a standard glucose (Sigma-Aldrich, Shanghai, China) solution was determined using the same method as the test wells. The glucose concentration is presented in millimoles per liter (mmol/L).
Total protein assay of MECs
The total protein content of the treated MECs was determined using a Coomassie Protein Assay Reagent (Jiancheng, Nanjing, China). The cells were lysed using a repeated freeze-thaw fragmentation method. Accordingly, the MECs were frozen at −80 °C for 60 min and transferred to a 37 °C water bath for 15 min to thaw the cells, which was repeated 3 times. Samples of the cell debris and contents were collected by adding 300 μL of a 0.9 % sodium chloride (NaCl) solution to each well. Double distilled water (blank control), a standard protein solution and sample liquid with an equal volume (50 μL) were mixed with 3.0 mL of reagent and incubated at room temperature for 10 min. Finally, the OD was recorded at a specific wavelength (595 nm) and optical path (1 cm) using a U3900 Spectrophotometer (Hitachi, Tokyo, Japan).
Real Time Quantitative Polymerase Chain Reaction (RT-qPCR)
Total RNA was extracted from the MECs using an RNAprep Pure Cell/Bacteria Kit (TIANGEN, Beijing, China). The purity and concentration of the total RNA was determined using a NanoDrop 2000 UV–vis Spectrophotometer (Thermo Scientific, Massachusetts, USA). Reverse transcription was performed with a PrimeScript® RT reagent Kit (Takara Biotechnology, Dalian, China), and the cDNA samples were stored at −20 °C until further analysis. The mRNA expression levels of the facilitative Na+-independent glucose transporters (GLUT1 and GLUT12), vitamin D receptor (VDR), calbindin-D 9k , PMCA1b, PMCA2b and Orai1 were measured using a SYBR® Premix Ex Taq™ II (Takara Biotechnology, Dalian, China). Briefly, a 20 μL reaction system was used that consisted of 10 μL of SYBR Premix Ex Taq II (2×), 0.8 μL of forward primer (10.0 μmol/L), 0.8 μL of reverse primer (10.0 μmol/L), 1 μL (500 ng) of cDNA, and 7.4 μL of RNase-free water. The reaction procedure was performed using an iCycler iQ5 multicolor real-time PCR detection system (Bio-Rad Laboratories, Hercules, CA) with the following program: 95 °C for 5 min; 35 cycles of 95 °C for 10 s, 60 °C for 30 s, and 72 °C for 30 s; and 72 °C for 5 min. All samples were run in triplicate, and the 2-△△Ct method, which was previously established by Livak , was adopted to analyze the gene expression data. The primers are presented in Table 1, and β-actin was used as a reference gene in this study.
After treatments, the supernatant fluid was removed and the cells were washed three times. Total protein was extracted using a High Performance RIPA buffer (Solarbio Science & Technology Co., Ltd., Beijing, China) in which the final concentration of phenylmethylsulfonyl fluoride (PMSF; Roche, Shanghai, China) was 1.0 mmol/L. The cells were collected in a 4 °C-precooled Eppendorf tube using a cell scraper, and the cells were lysed for 30 min at 4 °C. Afterward, the turbid liquid was centrifuged at a speed of 13,000 r/min for 10 min at 4 °C. The supernatant contained the total protein and was collected for further analysis. The western blot analysis was conducted according to the protocols reported by Xu et al. . Briefly, the protein content was determined using a Pierce™ bicinchoninic acid (BCA) Protein Assay Kit (Thermo Scientific, Rockford, USA), according to the manufacturer’s instructions. The total proteins were separated by SDS-PAGE, transferred to nitrocellulose membranes (Millipore, Billerica, USA), and then probed with the primary antibodies anti-PMCA1b, anti-PMCA2b and anti-β-actin, which were all purchased from Abcam (Cambridge, UK). Goat anti-rabbit IgG (Abcam, Cambridge, UK) was used as a secondary antibody. The chemiluminescent ECL western blot assay system (Thermo, Rockford, USA) was used to detect the signals.
Enzyme activity assay
A Hexokinase Test Kit (Jiancheng, Nanjing, China) was used to detect the HK activity of the solutions containing cell debris, and the samples were collected according to the user’s manual. The prepared reagent was pre-warmed at 37 °C for 10 min, and then 50 μL of liquid sample and 960 μL of reagent were immediately mixed in a tube to start the reaction. The absorbance at 340 nm (optical path: 0.5 cm) was recorded after 30 s (OD1) using a U3900 Spectrophotometer (Hitachi, Tokyo, Japan). Subsequently, the liquid was transferred back to the previous tube and warmed in a 37 °C water bath for 2 min. The absorbance was measured again under the same conditions and denoted as OD2. The HK activity was calculated using the following formula:
where “6.22” represents the millimolar extinction coefficient, “2” represents the reaction time (min), “0.5” represents the optical path (cm), and “1.01/0.05” refers to the dilution factor.
The Na+K+-ATPase and Ca2+Mg2+-ATPase activities were detected with a Trace ATPase Test Kit (Jiancheng, Nanjing, China). Protein samples were mixed with the appropriate reagents (different reagents for these two enzymes) and heated in a 37 °C water bath for 10 min; then, another reagent was added to the reaction system and centrifuged at 3,500 r/min for 15 min. The supernatants were collected to determine the inorganic phosphate (Pi) concentration. The Pi samples were treated with the appropriate reagents at room temperature for 2 min. Afterward, a final reagent was added and incubated at room temperature for 5 min. The OD values at 636 nm (optical path: 1 cm), including blank control (ODblank), control (ODcontrol), standard product (ODstandard) and sample (ODsample), were read using a Microplate Reader (Power Wave XS2, Bio Tek, USA). The formula to determine the protein concentration is as follows:
where “0.02” represents the concentration of the standard Pi solution (μmol/mL), “6” represents the reaction time (min), and “7.8” represents the dilution factor.
The data were subjected to one-way analysis of variance (ANOVA) using Statistical Product and Service Solutions 21.0 (SPSS 21.0; IBM SPSS Statistics, USA), and multiple comparisons were performed using Duncan’s method . The values were presented as the means ± SE (standard error). The results were declared significantly different if P < 0.05.
Supplementation with 1,25-(OH)2D3 significantly promoted MEC proliferation as the concentration increased from 0.1 to 10.0 nmol/L (P < 0.05, Fig. 1a), and no difference was observed between the control and the 0.1 nmol/L group (P > 0.05). Compared with the control, the rates of cell proliferation at the concentration of 0.1, 1.0, 10.0 and 100.0 nmol/L were increased by 3.79 %, 9.16 %, 15.99 % and 8.09 %, respectively. The cell proliferation rate in the 100.0 nmol/L group (P < 0.05) was lower than the 10.0 nmol/L group. In addition, the proliferation rate in the 100.0 nmol/L group was statistically equal to the 1.0 nmol/L group (P > 0.05).
Cell proliferation was inhibited in the 3-BrPA-supplemented group and the 3-BrPA plus 1,25-(OH)2D3 group (P < 0.05, Fig. 1b), and proliferation decreased by 37.85 % and 31.64 %, respectively. Increased cell proliferation was observed in the 1,25-(OH)2D3 group without 3-BrPA supplementation (P < 0.05). Whether or not the 1,25-(OH)2D3 was supplemented, no difference was observed in the MECs treated with 3-BrPA (P > 0.05).
The 0.1 nmol/L 1,25-(OH)2D3 treatment did not affect the glucose consumption by the goat MECs (P > 0.05, Fig. 2). The glucose uptake was significantly promoted when the 1,25-(OH)2D3 concentration increased from 0.1 to 10.0 nmol/L (P < 0.05). In accordance with cell proliferation, 100.0 nmol/L 1,25-(OH)2D3 decreased glucose consumption compared with the 10.0 nmol/L treatment (P < 0.05), and no differences were observed between 1.0 and 100.0 nmol/L (P > 0.05).
The expression of genes related to calcium transport in goat MECs were presented in Fig. 3. An increase in VDR expression was observed as the 1,25-(OH)2D3 levels increased from 0 to 10.0 nmol/L (P < 0.05), whereas no effect was observed between 10.0 and 100.0 nmol/L (P > 0.05). The same trend was observed for calbindin-D 9k , with the exception of an insignificant difference at 0.1 nmol/L compared with the control. In addition, supplementation with 10.0 and 100.0 nmol/L 1,25-(OH)2D3 increased PMCA1b expression (P < 0.05), and the peak PMCA1b expression level appeared at 10.0 nmol/L (P < 0.05). However, 1,25-(OH)2D3 had no influence on PMCA1b expression at concentrations of 0 and 1.0 nmol/L (P > 0.05).
The 1,25-(OH)2D3 supplementation altered the GLUT1 and GLUT12 gene expression levels as well (Fig. 4). There was an increase in GLUT1 mRNA abundance as the 1,25-(OH)2D3 levels increased from 0.1 to 10.0 nmol/L (P < 0.05, Fig. 4a). No difference was observed between the control and 0.1 nmol/L. However, compared with 10.0 nmol/L 1,25-(OH)2D3, the 100 nmol/L treatment did not increase GLUT1 expression (P > 0.05). Inconsistently, supplementation with 1,25-(OH)2D3 had no influence on GLUT12 expression when the concentration was less than 1.0 nmol/L (P > 0.05, Fig. 4b). The 10.0 nmol/L treatment promoted GLUT12 expression compared to the 1.0 nmol/L treatment (P < 0.05), and there was no difference between the 10.0 and 100.0 nmol/L treatments (P > 0.05).
Supplementation with 3-BrPA down-regulated PMCA1b and PMCA2b expression (P < 0.05, Fig. 5a and b), regardless of whether 1,25-(OH)2D3 was added. The expression levels of PMCA1b and PMCA2b in group D (10.0 nmol/L 1,25-(OH)2D3) were higher than those of the control. Specifically, the 1,25-(OH)2D3 treatment up-regulated PMCA1b expression in the 3-BrPA-supplemented groups (P < 0.05, Fig. 5a), but no difference in PMCA2b expression was observed (P > 0.05, Fig. 5b). As we could see from the immunoblots (Fig. 5c), the changes in the levels of the PMCA1b and PMCA2b proteins in the supplemented groups were similar to the changes in the transcripts.
As shown in Fig. 6, the expression levels of GLUT1 and Orai1 were increased by 1,25-(OH)2D3 supplementation (P < 0.05) and reduced by the addition of 3-BrPA (P < 0.05). No difference was observed between the 3-BrPA-supplemented group and 1,25-(OH)2D3 plus 3-BrPA group (P > 0.05).
Cell metabolic enzymes
As a whole, the enzyme activities, including HK, Ca2+Mg2+-ATPase and Na+K+-ATPase, were increased when the 1,25-(OH)2D3 levels increased from 0 to 10.0 nmol/L (Table 2). Compared with the 10.0 nmol/L treatment, decreased activities were detected in the 100.0 nmol/L group (P < 0.05). The HK activity in the 100.0 nmol/L group was statistically equal to the 0.1 nmol/L and control groups (P > 0.05). Supplementation with 0.1 nmol/L 1,25-(OH)2D3 did not affect the Ca2+Mg2+-ATPase and Na+K+-ATPase activities (P > 0.05), and no difference in Ca2+Mg2+-ATPase activity was observed between the 0.1 and 1.0 nmol/L groups (P > 0.05). The Na+K+-ATPase activity in the 100.0 nmol/L group was equivalent to the control (P > 0.05). Moreover, the Ca2+Mg2+-ATPase activity presented a sudden decrease at the highest 1,25-(OH)2D3 concentration, which was even lower than the control (P < 0.05).
1,25-Dihydroxyvitamin D3 is a natural ligand of the vitamin D receptor (VDR), and plays an important role in anti-inflammatory processes and calcium transport [31, 32]. It has been reported that 1,25-(OH)2D3 could activate the VDR to modulate gene transcription and mineral ion homeostasis [33, 34]. Vitamin D-facilitated calcium transport is a complicated process, including the up-regulation and down-regulation of associated genes. Calbindin-D9k, PMCAs and Orai were considered essential elements for transcellular calcium transport following stimulation with 1,25-(OH)2D3 [10, 11, 35, 36]. Our data showed that 1,25-(OH)2D3 influenced the expression of the VDR, calbindin-D 9k , PMCA1b, PMCA2b and Orai1 genes in goat MECs in a dose-dependent manner, which indicated enhanced calcium transport. Furthermore, we could infer that this process was closely related to cellular energy availability, based on the changes in GLUT1 and GLUT12 expression and the responses after the inhibition of HK2.
Supplementation with 1,25-(OH)2D3 improved cell proliferation in a concentration-dependent manner, with the exception of a relative decrease at 100.0 nmol/L. Our results were inconsistent with the results reported by Rayalam et al. , who found that 1,25-(OH)2D3 enhanced preadipocyte viability generated from 3 T3-L1 mouse embryo fibroblasts in a dose-dependent manner from 0.1 to 10.0 nmol/L, but no significant difference existed between the 10.0 and 100.0 nmol/L treatments. However, the proliferation of human airway smooth muscle cells (HASMCs) was gradually inhibited by increasing levels of 1,25-(OH)2D3 in another experiment . These variant effects might result from different cell types and functions as well as from the tolerated doses. Due to the high calcium content of milk, MECs assimilate large amounts of calcium from plasma. In addition, calcium is an essential element for cell growth, differentiation and maintenance. Consequently, it is plausible that the 1,25-(OH)2D3-induced promotion of calcium uptake can enhance MECs proliferation. To our knowledge, this was the first study in which 1,25-(OH)2D3-stimulated cell proliferation of secreting MECs was investigated.
Mammary lactation is a complicated biological process that is sustained by a variety of nutrients, among which glucose acts as the supreme precursor for lactose synthesis as well as an energy resource of metabolic activities . Hence, glucose plays an essential role in mammary milk secretion. It has been testified that glucose transporters (GLUTs) are the main tools for glucose uptake by mammary epithelial cells, and GLUT1 was the major transporter, although GLUT12 is involved as well [23, 39]. Previous studies rarely called attention to the effects of 1,25-(OH)2D3 on glucose uptake and metabolism. In our present study, 1,25-(OH)2D3 increased cell glucose consumption and up-regulated GLUT1 and GLUT12 expression, indicating that more glucose was utilized for cell metabolism or component synthesis. In addition, intracellular glucose phosphorylation catalyzed by HK is the first step in energy metabolism and is a rate-limiting process. Consequently, the increased HK activity was another persuasive indicator of glucose utilization [23, 40]. The main reason for the enhanced glucose consumption might be that 1,25-(OH)2D3-induced calcium transport led to the promotion of milk secretion in goat MECs. In addition, several studies have shown that 1,25-(OH)2D3 regulated the immune response in ruminants [41–43], which also required energy to sustain the process.
1,25-(OH)2D3 is a flexible secosteroid and exerts its regulatory functions by binding to VDR, a specific nuclear receptor and DNA-binding transcription factor . A series of biological processes, such as maintaining calcium homeostasis and mediating inflammation responses, are triggered by the binding between ligand and receptor . We found that 0 to 10.0 nmol/L 1,25-(OH)2D3 promoted VDR expression, with no difference between the 10.0 and 100.0 nmol/L treatments. This finding indicated that 1,25-(OH)2D3 could increase the number of VDRs in a dose-dependent manner, with an optimal concentration of 10.0 nmol/L. Haussler et al.  noted that the activation and function of VDR were induced by 1,25-(OH)2D3, but saturation was not mentioned. From the authors’ point of view, the cell metabolic capacity was limited and could not be induced in an unlimited manner. This hypothesis was supported by the results from a previous study by Rayalam et al. , who discovered that 1,25-(OH)2D3 could no longer promote adipocyte growth when the concentration exceeded 10.0 nmol/L.
The diffusion of intracellular calcium from the apical side to basolateral side depends on its binding to calbindin-D9k, and calcium passes through the basolateral side via PMCA1b [1, 5, 9, 10]. An overall increase in the calbindin-D9k and PMCA1b transcripts was detected when the 1,25-(OH)2D3 concentrations ranged from 0 to 10.0 nmol/L, which was a marker to distinguish the enhanced calcium transport. According to previous findings, both calbindin-D9k and PMCA1b had a vitamin D response element (VDRE) in their promoter region, and the VDRE was the direct binding site of VDR [37, 46, 47], which may be why 1,25-(OH)2D3 could regulate transcellular calcium transport. Moreover, there are other proteins that regulate cellular calcium transport. Using a null mutation mouse model, Reinhardt et al.  showed that the activity of PMCA2b, another isoform of PMCA, was required for the secretion of milk calcium, and Ji et al.  showed that 1,25-(OH)2D3 could stimulate PMCA2b expression to regulate mammary calcium transport. Davis et al.  suggested that Orai1, a novel channel, was important for mammary calcium transport during lactation. Orai1 is a key component of the CRAC channels and plays an extremely important role in the transmembrane influx of calcium [13, 14, 36]. The biology and molecular mechanism of Orai1 have been reviewed by Cahalan et al.  and Hogan et al. . The 1,25-(OH)2D3-stimulated up-regulation of PMCA2b and Orai1, together with their down-regulation by the inhibition of glucose metabolism, indicated that calcium transport in goat MECs could be regulated by 1,25-(OH)2D3 availability and the cellular energy status.
Plasma membrane Ca2+-ATPase is a transcellular Ca2+ transporter encoded by the PMCA gene family that plays a vital role in regulating cellular calcium metabolism and maintaining intracellular Ca2+ homeostasis [28, 50]. Ca2+Mg2+-ATPase activity showed a similar trend as the expression of PMCA1b and PMCA2b, indirectly indicating that calcium secretion was promoted when the 1,25-(OH)2D3 concentration did not exceed 10.0 nmol/L. Additionally, there was recent evidence showing that Na+/Ca2+ exchangers (NCX) on the mammalian plasma membrane co-modulated calcium transport with PMCA [50, 51]. Moreover, Zanatta et al.  found that 1,25-(OH)2D3 mediated transcellular calcium transport by stimulating NCX activation in rat Sertoli cells. Our data also showed an increase in Na+K+-ATPase activity as the 1,25-(OH)2D3 levels increased from 0 to 10.0 nmol/L. However, NCX expression was not examined in this study; therefore, we could not verify its regulatory role in the Ca2+ transcellular transport process.
Previous studies showed that 3-BrPA inhibited glycolysis in a dose-dependent manner by decreasing HK activity, particularly HK2; thus it has been widely used to investigate the impact of cellular energy status on biological processes [53, 54]. In our trials, the effect of energy availability on calcium transport in goat MECs was studied by supplementing the cells with 3-BrPA. Accordingly, cell proliferation and GLUT1 expression decreased, which was most likely due to the inhibition of glucose metabolism. In support of our findings, Yun et al.  described that glycolysis inhibitors, such as 3-BrPA, could inhibit cell and tumor growth at proper dosages. The decrease in PMCA1b and PMCA2b expression at the mRNA and protein levels, as well as down-regulated Orai1 transcription, attested that calcium transport was inhibited in goat MECs. Hence, 1,25-(OH)2D3 promoted calcium transport in goat MECs, and this process depended on the intracellular availability of glucose. It is well known that glucose is the main energy source of many metabolic activities, and active nutrient transport is a process that expends energy. Therefore, the inhibition of glycolysis reduced PMCA and Orai1 expression.
Compared with the 3-BrPA group, the 3-BrPA plus 1,25-(OH)2D3 group exhibited higher PMCA1b expression, whereas GLUT1 expression showed no difference, indicating that 1,25-(OH)2D3 could still enhance calcium transport when glucose uptake was suppressed in goat MECs. To our knowledge, this was a novel discovery. Many substances, such as clenbuterol  and conjugated linoleic acids (CLAs) , have been proven to induce nutrient repartition. We speculated that the stimulation of 1,25-(OH)2D3 repartitioned cellular energy for calcium secretion, but this assumption required convincing support. More trials are required to explore the roles of PMCAs, Orai1, NCX and other potential proteins. From the authors’ point of view, mammary calcium secretion is a complicated system, and multiple, cross-linked networks should be established via transcriptomics and proteomics technologies to better understand milk calcium synthesis. In addition, the isotope tracer technology should be used to directly reflect mammary calcium transport in dairy goats.
Suitable concentrations of 1,25-(OH)2D3 promoted proliferation and glucose utilization in goat MECs in a dose-dependent manner. Supplementation with 1,25-(OH)2D3 could modulate calcium transport by altering the expression of VDR, calbindin-D 9k , PMCA1b, PMCA2b and Orai1 in a dose- and energy-dependent manner. In the present study, the optimal concentration of 1,25-(OH)2D3 that stimulated the expression of calcium transport indicators in goat MECs was 10.0 nmol/L. Our findings highlighted the role of 1,25-(OH)2D3 as a potential regulatory agent to produce calcium-enriched milk in ruminants when sufficient intracellular energy was available.
Christakos S. Recent advances in our understanding of 1,25-dihydroxyvitamin D3 regulation of intestinal calcium absorption. Arch Biochem Biophys. 2012;523:73–6.
Christakos S. Mechanism of action of 1,25-dihydroxyvitamin D3 on intestinal calcium absorption. Rev Endocr Metab Disord. 2012;13:39–44.
Lee WJ, Monteith GR, Roberts-Thomson SJ. Calcium transport and signaling in the mammary gland: targets for breast cancer. BBA-Rev Cancer. 2006;1765:235–55.
Neville MC. Calcium secretion into milk. J Mammary Gland Biol. 2005;10:119–28.
Van Houten JN, Wysolmerski JJ. The calcium-sensing receptor in the breast. Best Pract Res Clin Endocr Metab. 2013;27:403–14.
Ardeshirpour L, Dumitru C, Dann P, Sterpka J, VanHouten J, Kim W, et al. OPG treatment prevents bone loss during lactation but does not affect milk production or maternal calcium metabolism. Endocrinology. 2015;156:2762–73.
Mamillapalli R, VanHouten J, Dann P, Bikle D, Chang WH, Brown E, et al. Mammary-specific ablation of the calcium-sensing receptor during lactation alters maternal calcium metabolism, milk calcium transport, and neonatal calcium accrual. Endocrinology. 2013;154:3031–42.
DeGaris PJ, Lean IJ. Milk fever in dairy cows: a review of pathophysiology and control principles. Vet J. 2008;176:58–69.
VanHouten JN, Wysolmerski JJ. Transcellular calcium transport in mammary epithelial cells. J Mammary Gland Biol. 2007;12:223–35.
VanHouten J, Sullivan C, Bazinet C, Ryoo T, Camp R, Rimm DL, et al. PMCA2 regulates apoptosis during mammary gland involution and predicts outcome in breast cancer. Proc Natl Acad Sci U S A. 2010;107:11405–10.
Jeong J, VanHouten JN, Dann P, Kim W, Sullivan C, Yu H, et al. PMCA2 regulates HER2 protein kinase localization and signaling and promotes HER2-mediated breast cancer. Proc Natl Acad Sci U S A. 2016;113:E282–90.
Cahalan MD, Zhang SL, Yeromin AV, Ohlsen K, Roos J, Stauderman KA. Molecular basis of the CRAC channel. Cell Calcium. 2007;42:133–44.
Hendron E, Wang X, Zhou Y, Cai X, Goto JI, Mikoshiba K, et al. Potent functional uncoupling between STIM1 and Orai1 by dimeric 2-aminodiphenyl borinate analogs. Cell Calcium. 2014;56:482–92.
Feske S, Gwack Y, Prakriya M, Srikanth S, Puppel SH, Tanasa B, et al. A mutation in Orai1 causes immune deficiency by abrogating CRAC channel function. Nature. 2006;441:179–85.
Lippolis JD, Reinhardt TA, Sacco RA, Nonnecke BJ, Nelson CD. Treatment of an intramammary bacterial infection with 25-hydroxyvitamin D3. PLoS One. 2011;6:e25479.
Chow EC, Quach HP, Vieth R, Pang KS. Temporal changes in tissue 1α,25-dihydroxyvitamin D3, vitamin D receptor target genes, and calcium and PTH levels after 1,25(OH)2D3 treatment in mice. Am J Physiol-Endocr Metab. 2013;304:E977–89.
Ji J, Lu RN, Zhou XJ, Xue YB, Shi CM, Goltzman D, et al. 1,25-Dihydroxyvitamin D3 contributes to regulating mammary calcium transport and modulates neonatal skeletal growth and turnover cooperatively with calcium. Am J Physiol-Endoc M. 2011;301:E889–900.
Leyssens C, Verlinden L, Verstuyf A. Antineoplastic effects of 1, 25 (OH) 2D3 and its analogs in breast, prostate and colorectal cancer. Endocr-relat Cancer. 2013;20:R31–47.
Mehta RG, Peng X, Alimirah F, Murillo G, Mehta R. Vitamin D and breast cancer: Emerging concepts. Cancer Lett. 2013;334:95–100.
Christakos S, Gabrielides C, Rhoten WB. Vitamin D-dependent calcium binding proteins: chemistry, distribution, functional considerations, and molecular biology. Endocr Rev. 1989;10:3–26.
Kohler M, Leiber F, Willems H, Merbold L, Liesegang A. Influence of altitude on vitamin D and bone metabolism of lactating sheep and goats. J Anim Sci. 2013;91:5259–68.
Wang Z, Luo J, Wang W, Zhao WS, Lin XZ. Characterization and culture of isolated primary dairy goat mammary gland epithelial cells. Chinese J Biotechnol. 2010;26:1123–7.
Zhao K, Liu HY, Wang HF, Zhou MM, Liu JX. Effect of glucose availability on glucose transport in bovine mammary epithelial cells. Animal. 2012;6:488–93.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2− ΔΔCT method. Methods. 2001;25:402–8.
Herm G, Muscher-Banse A, Breves A, Schröder B, Wilkens M. Renal mechanisms of calcium homeostasis in sheep and goats. J Anim Sci. 2015;93:1608–21.
Wilkens MR, Richter J, Fraser DR, Liesegang A, Breves G, Schröder B. In contrast to sheep, goats adapt to dietary calcium restriction by increasing intestinal absorption of calcium. Comp Biochem Phys A. 2012;163:396–406.
Sun YS. Regulation of calcium metabolism by calcium-sensing receptor (CaSR) in mammary gland. Mater Thesis. Yangling: Northwest A & F University; 2014.
Davis FM, Janoshazi A, Janardhan KS, Steinckwich N, D’Agostin DM, Petranka JG, et al. Essential role of Orai1 store-operated calcium channels in lactation. Proc Natl Acad Sci U S A. 2015;112:5827–32.
Xu HF, Luo J, Zhao WS, Yang YC, Tian HB, Shi HB, et al. Overexpression of SREBP1 (sterol regulatory element binding protein 1) promotes de novo fatty acid synthesis and triacylglycerol accumulation in goat mammary epithelial cells. J Dairy Sci. 2016;99:783–95.
IBM SPSS. SPSS Base 21.0 User’s Guide. Chicago: IBM SPSS Statistics; 2012.
Rosli S, Shintani T, Hayashido Y, Toratani S, Usui E, Okamoto T. 1α,25(OH)2D3 down-regulates HBp17/FGFBP-1 expression via NF-kB pathway. J Steroid Biochem Mol Biol. 2012;136:98–101.
Nelson C, Beitz DC, Reinhardt TA, Lippolis J. 1,25-Dihydroxyvitamin D3 enhances bovine mammary epithelial innate immune responses. Anim Ind Rep. 2010;656:9.
Tarroni P, Villa I, Mrak E, Zolezzi F, Mattioli M, Gattuso C, et al. Microarray analysis of 1,25(OH)2D3 regulated gene expression in human primary osteoblasts. J Cell Biochem. 2012;113:640–9.
Bouillon R, Carmeliet G, Verlinden L, Van Etten E, Verstuyf A, Luderer HF, et al. Vitamin D and human health: lessons from vitamin D receptor null mice. Endocr Rev. 2008;29:726–76.
Kutuzova GD, DeLuca HF. Gene expression profiles in rat intestine identify pathways for 1,25-dihydroxyvitamin D3 stimulated calcium absorption and clarify its immunomodulatory properties. Arch Biochem Biophys. 2004;432:152–66.
Niemeyer BA. Changing calcium: CRAC channel (STIM and Orai) expression, splicing, and posttranslational modifiers. Am J Physiol-Cell Physiol. 2016;310:C701–9.
Rayalam S, Della-Fera MA, Ambati S, Yang JY, Park HJ, Baile CA. Enhanced effects of 1,25 (OH) 2D3 plus genistein on adipogenesis and apoptosis in 3T3-L1 adipocytes. Obesity. 2012;16:539–46.
Song YF, Qi HW, Wu CG. Effect of 1,25-(OH)2D3 (a vitamin D analogue) on passively sensitized human airway smooth muscle cells. Respirology. 2007;12:486–94.
Yu QH, Zhu LQ, Lin J, Zhang Q, Tian Q, Hu WW, et al. Functional analyse of GLUT1 and GLUT12 in glucose uptake in goat mammary gland epithelial cells. PLoS One. 2013;8:e65013.
Liu HY, Zhao K, Liu JX. Effects of glucose availability on expression of the key genes involved in synthesis of milk fat, lactose and glucose metabolism in bovine mammary epithelial cells. PLoS One. 2013;8:e66092.
Waters W, Nonnecke B, Rahner T, Palmer M, Whipple D, Horst R. Modulation of mycobacterium bovis-specific responses of bovine peripheral blood mononuclear cells by 1, 25-dihydroxyvitamin D3. Clin Diagn Lab Immun. 2001;8:1204–12.
Nelson CD, Reinhardt TA, Lippolis JD, Sacco RE, Nonnecke BJ. Vitamin D signaling in the bovine immune system: a model for understanding human vitamin D requirements. Nutrients. 2012;4:181–96.
Nelson C, Reinhardt T, Thacker T, Beitz D, Lippolis J. Modulation of the bovine innate immune response by production of 1α,25-dihydroxyvitamin D3 in bovine monocytes. J Dairy Sci. 2010;93:1041–9.
Haussler MR, Jurutka PW, Mizwicki M, Norman AW. Vitamin D receptor (VDR)-mediated actions of 1α,25(OH)2vitamin D3: genomic and non-genomic mechanisms. Best Pract Res Clin Endocr Metab. 2011;25:543–59.
Wu SP, Liao AP, Xia YL, Li YC, Li JD, Sartor RB, et al. Vitamin D receptor negatively regulates bacterial-stimulated NF-kB activity in intestine. Am J Pathol. 2010;177:686–97.
Nemazannikova N, Antonas K, Dass CR. Vitamin D: metabolism, molecular mechanisms, and mutations to malignancies. Mol Carcinogen. 2013;53:421–31.
Ma YY, Yu WD, Hidalgo AA, Luo W, Delansorne R, Johnson CS, et al. Inecalcitol, an analog of 1,25D3, displays enhanced antitumor activity through the induction of apoptosis in a squamous cell carcinoma model system. Cell Cycle. 2013;12:743–52.
Reinhardt TA, Lippolis JD, Shull GE, Horst RL. Null mutation in the gene encoding plasma membrane Ca2+-ATPase isoform 2 impairs calcium transport into milk. J Biol Chem. 2004;279:42369–73.
Hogan PG, Lewis RS, Rao A. Molecular basis of calcium signaling in lymphocytes: STIM and ORAI. Annu Rev Immunol. 2010;28:491–533.
Okumura R, Shibukawa Y, Muramatsu T, Hashimoto S, Nakagawa KI, Tazaki M, et al. Sodium-calcium exchangers in rat ameloblasts. J Pharmacol Sci. 2010;112:223–30.
Yang H, Kim TH, Lee HH, Choi KC, Jeung EB. Distinct expression of the calcium exchangers, NCKX3 and NCX1, and their regulation by steroid in the human endometrium during the menstrual cycle. Reprod Sci. 2011;18:577–85.
Zanatta L, Zamoner A, Gonçalves R, Zanatta AP, Bouraïma-Lelong H, Carreau S, et al. 1α,25-Dihydroxyvitamin D3 signaling pathways on calcium uptake in 30-day-old rat Sertoli cells. Biochemistry. 2011;50:10284–92.
Yun J, Rago C, Cheong I, Pagliarini R, Angenendt P, Rajagopalan H, et al. Glucose deprivation contributes to the development of KRAS pathway mutations in tumor cells. Science. 2009;325:1555–9.
Cardaci S, Rizza S, Filomeni G, Bernardini R, Bertocchi F, Mattei M, et al. Glutamine deprivation enhances antitumor activity of 3-bromopyruvate through the stabilization of monocarboxylate transporter-1. Cancer Res. 2012;72:4526–36.
Mokhtari B, Pourabdollah K. Preparation and characterization of bonded-phases of nano-baskets bearing sulfonyl-carboxamide. J Liq Chromatogr Relat Technol. 2012;35:2910–22.
Stangl GI. Conjugated linoleic acids exhibit a strong fat-to-lean partitioning effect, reduce serum VLDL lipids and redistribute tissue lipids in food-restricted rats. J Nutr. 2000;130:1140–6.
We really appreciated Dr. Xiaofei Wang from the Institute of Animal Nutrition and Feed Science, Northwest A&F University, China, for providing materials of cell culture. We also expressed the heartfelt gratitude to Dr. Kang Yu from the Faculty of Medicine and Dentistry, University of Alberta, Canada, for the assistance in the isolation of goat MECs and determinations of gene and protein expressions.
The research was supported by the National Key Technologies R&D Program of China (2012BAD12B02 and 2012BAD39B05-2), the National Funds for Natural Science of China (31472122), and Northwest A&F University Ph.D. Research Start-up funds (Z111021309).
Availability of data and materials
All the datasets were presented in the main manuscript and available to readers.
FFS conceived and designed the experiments. FFS and YCC conducted the experiments. CY and XSW assisted with the analysis of cell proliferation and enzyme activities. YCC performed the statistical analysis of the experimental data. Finally, the paper was written by FFS and modified by JHY. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate