key: cord-0982275-qr8bzp4r
authors: Kutuzova, Galina D.; DeLuca, Hector F.
title: Gene expression profiles in rat intestine identify pathways for 1,25-dihydroxyvitamin D3 stimulated calcium absorption and clarify its immunomodulatory properties
date: 2004-12-15
journal: Archives of Biochemistry and Biophysics
DOI: 10.1016/j.abb.2004.09.004
sha: eb8fd1687558675ae897364349cfdb7cff603ea9
doc_id: 982275
cord_uid: qr8bzp4r

Abstract Microarray technology has been used to discover 1,25-dihydroxyvitamin D3 (1,25-(OH)2D3) induced gene expression changes in rat small intestine in vivo. Here, we report gene expression changes related to intestinal absorption or transport, the immune system and angiogenesis in response to 1,25-(OH)2D3. Vitamin D deficient rats were intrajugularly given vehicle or vehicle containing 730ng of 1,25-(OH)2D3/kg of body weight. Intestinal mRNA was harvested from duodenal mucosa at 15min, 1, 3, and 6h post-injection and studied by Affymetrix microarrays. Genes significantly affected by 1,25-(OH)2D3 were confirmed by quantitative RT-PCR with remarkable agreement. The most strongly affected gene in intestine was CYP24 with 97-fold increase at 6h post-1,25-(OH)2D3 treatment. Intestinal calcium absorption genes: TRPV5, TRPV6, calbindin D9k, and Ca2+ dependent ATPase all were up-regulated in response to 1,25-(OH)2D3, supporting the currently accepted mechanism of 1,25-(OH)2D3 induced transcellular calcium transport. However, a 1,25-(OH)2D3 suppression of several intra-/intercellular matrix modeling proteins such as sodium/potassium ATPase, claudin 3, aquaporin 8, cadherin 17, and RhoA suggests a vitamin D regulation of tight junction permeability and paracellular calcium transport. Several other genes related to the immune system and angiogenesis whose expression was changed in response to 1,25-(OH)2D3 provided evidence for an immunomodulatory and anti-angiogenic role of 1,25-(OH)2D3.

The active form of vitamin D 3 , 1,25-dihydroxyvitamin D 3 or calcitriol (1,25-(OH) 2 D 3 ), 1 is a seco-steroid hormone that in association with high affinity vitamin D 3 receptor (VDR), a ligand-activated transcription factor, transactivates or transrepresses a variety of genes [1] [2] [3] [4] . For decades the major physiologic function of vitamin D 3 was believed to be a stimulation of intestinal calcium and phosphorous absorption [1] . However, the detailed mechanism of vitamin D-induced intestinal cal- cium absorption is still not fully delineated. The current model of vitamin D-induced transcellular calcium absorption in intestine is the induction of (1) calcium channel proteins, TRPV5 and TRPV6 (transient receptor potential channels vaniloid subfamily, former names ECaC1 and ECaC2 [5] [6] [7] [8] ); (2) cellular calcium transfer protein, calbindin-D 9K [1, 4, 9] ; and (3) a calcium extrusion protein, calcium ATPase (PMCA 1b ) [10] . PMCA 1b has been implicated in the regulation of calcium absorption by vitamin D, but its expression is not as tightly regulated by 1,25-(OH) 2 D 3 as that of TRPV5 and TRPV6 and calbindin-D 9K . Evidence for this regulation being transcriptional is only partly substantiated [11] .

The data on regulation of epithelial calcium channels (TRPV5 and TRPV6), calbindin-D 9K , and PMCA 1b expression by 1,25-(OH) 2 D 3 seem clear and these genes have VDRE in their promoter region [12] [13] [14] [15] . However, administration of 1,25-(OH) 2 D 3 did not change the expression of calbindin D 9K and PMCA 1b mRNA levels in the ovariectomized rats [16] . Also, there is no correlation between calbindin D 9K expression and plasma level of 1,25-(OH) 2 D 3 in dairy cattle [17] . A report has also appeared on the absence of classical calbindin D 9k -mediated mechanisms for active Ca 2+ transport in sheep rumen [18] . Vitamin D receptor (VDR) knockout mice revealed VDR-independent regulation of intestinal calcium absorption and levels of TRPV6 and calbindin D 9K mRNA, raising doubts about the role of calbindin D 9K as calcium ferry protein [19, 20] . High dietary Ca 2+ intake in VDR knockout mice resulted in a decreased expression of both TRPV5 and TRPV6 and in a reduction in calbindin D 9K and PMCA1b expression [21] . It is currently recognized that other genes/proteins could be involved in calcium uptake process [10] .

The present study represent a comprehensive gene expression study in vivo of the effect of 1,25-(OH) 2 D 3 on rat intestine. To specify the role of vitamin D in intestinal calcium absorption and to identify new genes that may be involved in these processes, we have used highdensity oligonucleotide arrays and studied the dynamics of gene expression in rat duodenum within 6 h postintrajugular injection of 1,25-(OH) 2 D 3 .

Animals were maintained and research was conducted in accordance with guidelines set forth by the Animal Care and Research Committee (University of Wisconsin, Madison, WI). Holtzman male weanling rats were obtained from Sprague-Dawley (Madison, WI) and maintained on a highly purified vitamin D-deficient diet, containing 0.47% calcium and 0.3% phosphorus (Pi) supplemented three times a week with 500 lg DL DL-a-tocopherol, 60 lg menadione, and 40 lg b-carotene in 0.1 ml soybean oil (AEK). Rats were housed in hanging wire cages and maintained on a 12 h light/dark cycle. Rats fed the vitamin D-deficient diet were maintained in a room with incandescent lighting, and all potential sources of ultraviolet light and vitamin D were excluded. At 14 weeks of age, blood was taken from the tail for measurement of serum calcium concentrations to assess vitamin D depletion.

Blood samples were obtained from the tail artery. Whole blood was centrifuged at 1100g for 15 min at 25°C to yield serum. Serum calcium concentration was determined using a 3110 atomic absorption spectrometer (Perkin-Elmer, Norwalk, CT) on serum diluted 1:40 with 1 g/L LaCl 3 [22] .

Vitamin D-deficient rats were given intrajugularly one dose of 730 ng of 1,25-(OH) 2 D 3 /kg of body weight in ethanol or vehicle (for control group) and a sample of blood was taken immediately before the injection for serum calcium concentration measurement. Groups of three rats per time point were deeply anesthetized with isoflurane and decapitated at 15 min, 1, 3, and 6 h after injection. Blood was collected at the same time for determination of changes in serum calcium concentration. The first 15 cm of intestine (duodenum) was removed, slit open longitudinally and scraped with the glass slide. The mucosa was homogenized with Power-Gen 700 (Fisher Scientific, Pittsburgh, PA) in guanidine thiocyanate (GTC) extraction buffer, supplemented with 2% b-mercaptoethanol (PolyATtract System 1000, Promega, Madison, WI), flash frozen in liquid N 2 , and stored at À80°C. Experiments were done in duplicate.

Poly(A) + RNA was isolated from pooled homogenized mucosa from three rats at each time point. The mRNA was isolated using the PolyATtract System 1000 (Promega, Madison, WI) and purified using an RNeasy kit (Qiagen, Chatsworth, CA). The quality, integrity, and quantity of the poly(A) + RNA was determined by agarose gel electrophoresis, UV absorption spectrophotometry, and Agilent Bioanalyser 2100 (Agilent Technologies, Palo Alto, CA).

Double stranded cDNA was synthesized from 3 lg of polyadenylated poly(A) + RNA using the Superscript Choice system (Invitrogen Life Technologies, Carlsbad, CA), all according to the Affymetrix Gene Expression manual (Affymetrix, Santa Clara, CA). Following phenol/chloroform extraction and ethanol precipitation, a biotin-labeled in vitro transcription reaction was performed using the cDNA template and BioArray High Yield In Vitro Transcription kit (Enzo Life Sciences, Farmingdale, NY). The cRNA was fragmented at 0.7 lg/ll final concentration in 1· fragmentation buffer (40 mM Tris-acetate, pH 8.1, 100 mM potassium acetate, and 30 mM magnesium acetate). The size range of cRNA before (0.5 kb and longer) and after base fragments) fragmentation was checked by agarose gel electrophoresis.

Affymetrix high-density rat oligonucleotide arrays (GeneChips RG-U34A) were synthesized photolithographically by the manufacturer using the UniGene 34 set of sequence clusters. Two sets of 16-25 base oligonucleotides each were used to probe each target sequenceperfect match (PM) and mismatch (MM) probe sets. Perfect match (PM) probe set and mismatch (MM) probe set were the same except MM contained a mismatched base in the center of the oligonucleotide. The MM probe set was used to control for non-specific hybridization of related sequences. The chip contained around 8800 probe sets, with about 10% of the (longer) sequences represented by more than one probe set. Target and probe set sequences were obtained from the Netaffx Analysis Center (http://www.affymetrix.com/ analysis/index.affx; Affymetrix).

The hybridization reaction and the automated hybridization procedure were performed by the Gene Expression Center at the Biotechnology Center of the University of Wisconsin-Madison. Each probe sample was tested on an Affymetrix Test3 Array and the quality of the cDNA and cRNA syntheses was determined by the 3 0 /5 0 ratio of housekeeping genes within the array (ubiquitin, rat glyceraldehyde 3-phosphate dehydrogenase, b-actin, and hexokinase). If the sample passed the quality control on the Affymetrix Test3 Array, it was hybridized to an Affymetrix high-density rat oligonucleotide array GeneChip U34A per protocol recommendation in the Affymetrix GeneChip Expression Analysis Technical Manual [see: http://www.affymetrix. com/support/technical/manual/expressionmanual.affx]. GeneChips were hybridized with 15 lg of biotin labeled cRNA in 300 ll, in 1· hybridization buffer [100 mM Mes, 1 M NaCl (Ambion), 20 mM EDTA (Ambion), and 0.01% Tween 20 (Pierce Chemical), 50 pM control oligonucleotide B2 (Affymetrix), 0.1 mg/ml Herring Sperm DNA (Promega), 0.5 mg/ml acetylated BSA (Invitrogen Life Technologies), and 1· Eukaryotic Hybridization Controls (BioB, BioC, BioD, and Cre at 1.5, 5, 25, and 100 pM, respectively) (Affymetrix)] for 16 h at 45°C on a rotisserie at 60 rpm. Prior to application to the GeneChip, samples were heated at 95°C for 5 min, followed by incubation at 45°C for 5 min and spun at 14,000g for 5 min. Following hybridization, the labeled samples were removed from the GeneChip, stored in the appropriate vial at À20°C, and immediately filled with non-stringent buffer A which contains 6· SSPE [0.9 M sodium chloride, 60 mM sodium phosphate, and 6 mM EDTA (Ambion)] and 0.01% Tween 20. All GeneChips were post-processed using the automated Affymetrix GeneChip Fluidics Station 400.

The post-processing protocol for the RG_U34 Genome GeneChip is as follows: Wash#1: 10 

To quantify the fluorescent signal from each feature on the GeneChip, all GeneChips were scanned at a wavelength of 570 nm using the Hewlett-Packard Gene Array Scanner. Fluorescent signals corresponding to hybridization intensities were analyzed with the Affymetrix Microarray Suite (MAS) 5.0 software using the following settings: Detection Algorithm defaults, a1, 0.04; a2, 0.06; sf, 0.015 and Comparison Algorithm defaults, c1H, 0.0025; c1L, 0.0025; c2H, 0.003; c2L, 0.003; Perturbation, 1.1. In all the analyses, all probe sets were scaled to a target signal of 1000 using the ''Scale'' function in the GeneChip software. To estimate the range of linearity, four different bacterial mRNAs were added to the hybridization cocktail at the following concentrations: BioB (1.5 pM), BioC (5 pM), BioD (25 pM), and Cre (100 pM) [23] .

Expression data were analyzed using the Affymetrix Microarrray Suite software version 5.0 (MAS 5.0). Comparison tables for each time point for 1,25-(OH) 2 D 3 vs. vehicle-treated rats were generated in EX-CEL (Microsoft). For each comparison, e.g., 1,25-(OH) 2 D 3 treated relative to control (vehicle treated), and for each cDNA represented in the array, a ratio (e.g., 1,25-(OH) 2 D 3 /control) and an absolute difference of intensities for 1,25-(OH) 2 D 3 and vehicle treated were calculated.

For real-time quantitative PCR (Q-PCR), reverse transcription reactions were carried out in parallel on DNaseI digested pooled mRNA from 1,25-(OH) 2 D 3 or vehicle treated rats at different time points. Prior to reverse transcription, total RNA was confirmed to be free of contaminating DNA sequences by PCR using rat b-actin specific primer pairs designed to differentiate between cDNA, genomic DNA, and pseudogene genomic DNA. Reverse transcription reactions (five replicas) were done for each time point for 1,25-(OH) 2 D 3 -or vehicle-treated rats. Reverse transcription reactions were as follows: 50 ng mRNA, 2.5 ll oligo(dT) 16 (500 ng/ll), 5 ll dNTPmix (10 mM each dNTP), 1 pg artificial transcript of human b-actin (IDT, Coralville, IA), and H 2 O to 29.75 ll. Samples were incubated at 70°C for 10 min, put briefly on ice, and then incubated at 42°C for 2-5 min. Mix#2 (10 ll 5· Superscript II first strand buffer, 5 ll of 0.1 M DTT, 4 ll of 25 mM MgCl 2 , and 1.25 ll Superscript II RNAse H $ Reverse Transcriptase 200 U/ll) was added, mixed, and samples were immediately returned to incubate at 42°C for 1 h. Reactions were stopped by incubation at 70°C for 15 min.

To confirm uniformity of reverse transcription within sample and control groups, Q-PCRs were done in the presence of artificial transcript (human b-actin). Then, reverse transcription reactions with similar quantities of human b-actin within a sample group were pooled. Quantitative PCRs were performed using a GeneAmp 5700 Sequence Detection System (Perkin-Elmer, Norwalk, CT), using the ''standard-curvequantitation'' method [24] . Each reaction contained target-specific forward and reverse primers (200-750 nM final concentration, Table 1 ), 2· SYBR Green Master Mix (Applied Biosystems, Foster City, CA), 5ll of a 1:10 dilution of pooled reverse transcription product and H 2 O to a total volume of 25 ll. A two-step PCR profile was used: 10 min at 95°C denaturation and Amplitaq Gold activation, followed by 40 cycles alternating between 95°C for 15 s and 60°C for 60 s. Dilution series (1:2; 1:10; 1:50; 1:250; and 1:1250) standard curves were performed in quadruplicates for each primer pair using reverse transcription products described above. PCR was done in five replicas for each sample and relative quantities were determined based Table 1 Sequences of primers and fold change in expression of selected genes chosen for confirmation study by real-time quantitative PCR GenBank accession number Description Real-time PCR primers and probe sets were selected for each cDNA by using PRIMER EXPRESS software (Ver. 1.0; Applied Biosystems, Foster City, CA) and are presented in Table 1 .

We studied differential gene expression profiles in rat intestine after a single intrajugular injection of 1,25-(OH) 2 D 3 with the purpose of identifying novel genes involved in intestinal Ca 2+ and other nutrient absorption. It was shown previously [25] that serum concentration of Ca 2+ in the plasma begins to increase 3 h after treatment with 1,25-(OH) 2 D 3 , peaks at about 6 h, and declines at 12 h. We, therefore, examined gene expression in rat intestine at: 15 min, 1, 3, and 6 h.

We used Affymetrix Rat GeneChips U-34A array that contains 8799 known rat transcripts (77%) and ESTs (23%). In comparison, tables (sample vs. control) of gene expression (MAS 5.0), only genes considered (P) with a statistically valid signal increase (change ''I'') were considered genes upregulated by 1,25-(OH) 2 D 3 . Only genes present (P) in control with a statistically valid signal decrease in the sample (change ''D'') were considered as down-regulated. To identify genes that were differentially expressed between 1,25-(OH) 2 D 3 (sample) and vehicle (control) treated animals for each time point, we arbitrarily set up cut-off values to 1.5 for the fold change in ratio. In some cases, it was hard to assign the reliable fold change for the genes that were absent (A) in control and become present (P) in the sample or vise versa.

We used RT-PCR to confirm the effect of 1,25-(OH) 2 D 3 on regulated genes. The list of genes confirmed, maximum fold change in their expression after the stimulation with 1,25-(OH) 2 D 3 , and primers used are presented in Table 1 . The absolute magnitude of fold regulation detected with RT-PCR technique was always similar or higher than the fold change detected by microarray analysis, except for plasma membrane CaATPase 2b (Tables 1 and 2 ). Microarray technology can reliably detect changes in gene expression as subtle as 1.3-to 2-fold [26] . We were anxious not to overlook genes that could be very important in understanding of vitamin D mechanism of action that may only change by a factor of 2 (cut-off value currently accepted by users). An example is the well-established vitamin D responsive calbindin D 9k gene. Our GeneChip data showed its maximal up-regulation only 1.6-fold at 3 h after 1,25-(OH) 2 D 3 treatment ( Table 2) .

All 1,25-(OH) 2 D 3 regulated genes that passed the selection criteria (see above) were classified in terms of their function by referring to the literature and Affymetrix Analysis Center Web site and links (see Materials and methods). The data on 1,25-(OH) 2 D 3 stimulated gene expression are presented on separate tables and are given at the time of expression maximum fold change. We did not provide the fold change at other time points, which we observed, to avoid the complexity in presentation and to conserve space.

In this paper, we have restricted our presentation to 1,25-(OH) 2 D 3 -stimulated differential expression of genes dealing with digestion, absorption, and the immune system.

In our experiment, 1,25-(OH) 2 D 3 stimulated the highest level of expression of 25-hydroxyvitamin D 3 24-hydroxylase (CYP24), the major enzyme of 1,25-(OH) 2 D 3 degradation pathway, compared to all other transcripts which is consistent with the previous findings on strong up-regulation of this enzyme by 1,25-(OH) 2 D 3 both in vivo and in vitro [1] . As we observed, CYP24 mRNA was undetectable in the intestine of vehicle treated rats but after 1,25-(OH) 2 D 3 injection its level increased 84-fold at 3 h and 97-fold at 6 h.

We observed the increased expression of genes considered to be directly involved in the intestinal Ca 2+ absorption. The maximum fold change of the expression level of calbindin D 9k -the vitamin D-dependent cytosolic calcium binding protein within 6 h after the treatment, was 1.6-fold at 3 h (at 1 h after injection there was 1.4-fold increase) ( Table 2 ). Plasma membrane Ca 2+ ATPase transcript (EST AI103671) was not detectable at all time points in the control (vehicle treated) rats, or at 15 min and 1 h in 1,25-(OH) 2 D 3 -treated animals and had 8.6-fold expression increase at 3 h (2-fold by Q-PCR) followed by a further increase in transcript level at 6 h. The activity of this Ca 2+ ATPase is regulated by calmodulin, which also showed a maximal 1.6-fold increase at 3 h (Table 2) . Calmodulin, as a major intracellular Ca 2+ sensor and modulator, is involved in numerous calcium signaling pathways by interaction with diverse group of cellular proteins [27] . Calmodulin antagonists significantly reduced 1,25-(OH) 2 D 3 stimulated calcium uptake in Caco-2 colon derived cells, implying its involvement in Ca 2+ influx [28] . It is well known that blood calcium is regulated by several calcitropic hormones, e.g., calcitonin, parathy-roid hormone (PTH), and 1,25-(OH) 2 D 3 . Caldecrin, a serum calcium-decreasing factor, is a chymotrypsin-type serine protease, which belongs to the elastase family and inhibits parathyroid hormone or parathyroid hormonerelated, peptide-induced bone resorption. Caldecrin is synthesized as preprocaldecrin and is secreted from the cell. Preprocaldecrin does not possess serum calcium-decreasing activity but acquires it as well as protease activity, upon trypsin treatment [29] . In our experiments, 1,25-(OH) 2 D 3 stimulated preprocaldecrin expression 1.5-fold at 3 h.

The Affymetrix Rat Genome U34A Array used in our study did not have probe sets for the epithelial calcium channels TRPV5 and TRPV6, which are considered to be the major channels for calcium entry in intestine [5] [6] [7] . We analyzed the regulation of expression of TRPV5 and TRPV6 channels by 1,25-(OH) 2 D 3 in rat intestine within the time frame of our study (6 h) using real-time quantitative PCR (Table 1) . A slight increase in expression of TRPV5 was detected at 3 h, but at 6 h, it increased more than 10-fold ( Fig. 1) . Expression of TRPV6 channel started to increase at 1 h (3-fold increase) and continued to increase up to 9.8-fold at 6 h ( Fig. 1 ).

In Table 3 , we present the list of genes involved in intestinal transport of different compounds that were differentially expressed in intestine of rats within 6 h after administration of 1,25-(OH) 2 D 3 . Prepro-uroguanylin À1.6 a These genes also showed up-or down-regulation with other probe sets derived from different GenBank Accession numbers of the same protein.

At 1 h, expression of very long-chain acyl-CoA synthetase (VLACS) mRNA was decreased 2.7-fold (Table  3) . VLACS catalyzes the activation of very long-chain fatty acids (VLCFAs) to their CoA thioesters in peroxisomes. VLACS is identical with fatty acid transporter protein 2 (FATP2) and is involved in cellular uptake and metabolism of very long-chain fatty acids [30] .

At 3 h, up-regulation by 1,25-(OH) 2 D 3 was observed for genes that are currently known to be involved in vesicle transport (endocytosis and exocytosis) between secretory compartments like Cis-Golgi p28 protein (1.7-fold increase) and Sip9, a zymogen granule protein (syncollin, 1.5-fold increase), though their other possible functions remain to be elucidated (Table 3 ). p28 is the protein involved in a docking and fusion stage of protein transport from endoplasmic reticulum to Golgi [31] . Sip9 may function as a calcium-sensitive regulator of exocytosis in exocrine tissues triggered by the release of Ca 2+ from intracellular stores. Membrane-associated syncollin is a component of lipid rafts, where it binds to lipid bilayers in a cholesterol-dependent manner and may form pores [32] .

At 3 h, 1,25-(OH) 2 D 3 down-regulated expression the number of transporters and channels. They include a-1 and b-1 subunits of sodium/potassium ATPase (Na,K-ATPase), voltage dependent anion channel 1, transmembrane protein rnp21.4, fatty acid translocase/CD36 (FAT/36), claudin-3, urate transporter/channel, hypothetical ATPase, and aquaporin-8 ( Table 3) .

The Na,K-ATPase, or the sodium pump, consists of a non-covalently linked a-and b-subunits and is a key enzyme that regulates the intracellular Na + and K + homeostasis in animal cells [33] . We observed a 3-fold decrease in expression of Na,K-ATPase a-1 subunit and a 2.0-fold decrease in expression of b-1 subunit (Table 3) . It was shown previously that inhibition of Na,K-ATPase in tight monolayers of epithelial cells with ouabain results in an increase in tight junction permeability to ions and non-ionic molecules [34] .

1,25-(OH) 2 D 3 caused a 2.7-fold down-regulation of voltage-dependent anion-selective channels (VDACs). Voltage-dependent anion-selective channel (VDAC), or mitochondrial porin, forms a channel through the mitochondrial outer membrane and also through the plasma membrane, which allows the diffusion of a small anionic molecules [35] .

At 3 h, intestinal expression of transmembrane protein rnp21.4, involved in vesicular trafficking and protein transport [36] , decreased 2.5-fold (Table 3 ). Galectin-9 that showed a 2-fold expression decrease in the response to 1,25-(OH) 2 D 3 is a 36 kDa b-galactoside binding lectin and was identified as a specific urate transporter/channel. Galectin-9 is also involved in a variety of biological functions such as cell activation, proliferation, adhesion, and apoptosis [37] . The expression of hypothetical putative P-type ATPase that possi-bly catalyzes the transport of heavy metal ions was suppressed by 1,25-(OH) 2 D 3 (Table 3) .

Fatty acid translocase/CD36 (FAT/36), which mediates the long-chain fatty acid (LCFA) uptake, was down-regulated by 1,25-(OH) 2 D 3 1.9-fold at 3 h (Table  3 ). In rodent models of obesity and Type 1 diabetes, increased expression of FAT/CD36 (Type 1 diabetes) has increased the LCFA uptake into the heart and muscle [38] .

We also observed a 1.5-fold decrease in the expression of aquaporin-8 (Table 3) . Aquaporins are a conserved group of small hydrophobic integral membrane proteins that form a barrel-like channels for water transport [39] .

At 6 h after 1,25-(OH) 2 D 3 treatment, we observed a significant increase in expression of two transporters ( Table 3) . Organic cation transporter (OCT1) expression was up-regulated more than 2-fold with different probe sets. OCT1 is the member of the polyspecific potentialsensitive organic cation transporter gene family and belongs to the class of proteins responsible for transport of sugar and organic cations that include endogenous compounds such as monoamine neurotransmitters, choline, and coenzymes, but also numerous drugs and xenobiotics [40] .

1,25-(OH) 2 D 3 caused more than a 2-fold increase at 6 h in the expression of the multidrug transporter 1 or ATP-binding cassette (ABCB) transporter or the multidrug-resistance/transporter associated with the antigen processing (MDR/TAP1), which is involved in the transport of peptide antigens from the cytoplasm into a membrane-bound compartment of endoplasmic reticulum (ER) for association with MHC class I molecules. MDR/TAP1 is also the part of ER chaperone complex and functions in association with calnexin and calreticulin-the two lectins (carbohydrate binding chaperones), which interact with and assist the folding of proteins that carry monoglucosylated N-linked glycans. Calnexin expression was also increased 2.2-fold, 6 h in response to 1,25-(OH) 2 D 3 (Table 3 ). In some carcinomas and malignant tumors, the transcription of TAP1 is drastically reduced [41] .

Natural resistance-associated macrophage protein 2 (Nramp2) or solute carrier family 11 member 2 expression was reduced 2.2-fold by 1,25-(OH) 2 D 3 at 6 h ( Table 3 ) and this was seen with multiple probe sets. Nramp2 is a broad specificity divalent-metal transporter and is expressed at the duodenum brush border where it is responsible for transferrin-independent uptake of dietary iron from the intestinal lumen [42] .

1,25-(OH) 2 D 3 stimulated differential expression of genes involved in intra-/intercellular matrix modeling 1,25-(OH) 2 D 3 regulated the expression of several genes involved in intracellular and intercellular structure formation (Table 4) . (Table 4 ). It caused a 1.7-fold increase in expression of the lysosomal-associated membrane protein 2 (LAMP-2). LAMP-2 functions as the receptor for the selective uptake and degradation of cytosolic proteins by lysosomes and is involved in chaperonemediated autophagy and lysosomal biogenesis. Remarkably, LAMP-2 deficiency in humans leads to Danon disease, a lysosomal glycogen storage disease, cardiomyopathy, and myopathy [43] .

At 3 h, 1,25-(OH) 2 D 3 treatment increased the expression of other cytoskeleton proteins. These include a tropomyosin non-muscle isoform NM3 (TPM-c), a component of the actin microfilament complex; cytoplasmic dynein, involved in numerous intracellular motility events, including the transport and positioning of vesicles and organelles, spindle assembly and morphogenesis, and ankyrin, the adaptor molecule that links ion channels and cell adhesion molecules to the spectrin-based skeleton in specialized membrane domains (Table 4) .

At the same time (3 h), 1,25-(OH) 2 D 3 induced significant down-regulation of several genes involved in cellcell interactions and cell motility ( Table 4 ).

The expression of cadherin-17 was decreased 2.6-fold by 1,25-(OH) 2 D 3 ( Table 4 ). The cadherinÕs superfamily of transmembrane glycoproteins is calcium-dependent cell adhesion molecules that play an active role in tissue morphogenesis, patterning and cell adhesion [44] .

1,25-(OH) 2 D 3 treatment inhibited the expression of the brush border myosin-I (BBMI) (2.5-fold at 3 h) (Table 4). BBMI is a major component of the actin assembly in the microvilli of intestinal cells and is involved in endocytosis. Non-functional truncated BBMI affected the rate of transferrin recycling and the rate of transepi-thelial transport of dipeptidyl-peptidase IV from the basolateral plasma membrane to the apical plasma membrane [45] .

RhoA GTPase, a member of the family GTP-binding proteins (small GTPases), is the key regulator of actin cytoskeletal dynamics in cells. Its expression was significantly inhibited (2.5-fold, 3 h) by 1,25-(OH) 2 D 3 ( Table  4 ). The Rho proteins, Cdc42, Rac1, and RhoA, regulate the actin cytoskeleton. Cdc42 and Rac1 are primarily involved in the formation of protrusive structures, while RhoA generates myosin-based contractility. RhoA GTPase is involved in the formation of stress fibers and in the regulation of tight junctions structure and function as well as in generation of myosin-based contractility. All mutations in RhoA induced time-dependent disruptions in epithelial gate function [46] . At 3 h, claudin-3 expression decreased 2.2-fold after 1,25-(OH) 2 D 3 treatment (Table 4 ). Claudins are membrane proteins involved in the formation of tight junctions strands. Elevated expression of claudin-3 was detected in prostate adenocarcinoma and in the primary ovarian carcinoma cells [47] . It has been suggested that claudin-3 is a central component determining the integrity of blood-brain barrier tight junctions in vivo [48] .

At 6 h, 1,25-(OH) 2 D 3 increased expression of the following proteins (Table 4 ): (1) dynactin subunit p62, an accessory factor that mediates targeting of cytoplasmic dynein (dynein expression was increased at 3 h as we showed above) to a wide array of intracellular substrates; (2) mucin 3-the protein, that interacts with actin cytoskeleton; and (3) cell division cycle 42 homolog (Cdc42), a member of Rho family GTPases, which regulate the actin cytoskeleton and are primarily involved in the formation of protrusive structures. In mammalian Highly similar to myosin light chain alkali, smooth-muscle isoform À1.7 a These genes also showed up-or down-regulation with other probe sets derived from different GenBank Accession numbers of the same protein.

cells, Cdc42 has been implicated in the regulation of an actin rearrangement, receptor-mediated signal transduction pathways, cell cycle progression, apoptosis and multiple membrane trafficking events, including phagocytosis, exocytosis, and endocytosis. Cdc42 may also govern pathways required for establishment and maintenance of cellular polarity [49] . Expression of skeletal muscle a-actin (a-SMA) and thymosin b-10 (TB-10)-the intracellular matrix related genes-was strongly inhibited at 6 h by 1,25-(OH) 2 D 3 (Table 4 ). Our data are similar with published data on 1,25-(OH) 2 D 3 suppression of the expression of a-SMA, that was accompanied by a significant inhibition of mesangial cell proliferation [50] . Elevated expression of thymosin b-10 was shown to be an indicator of a malignant phenotype of human tumor cells, including mammary carcinomas. Suppression of thymosin b-10 protein synthesis reduced the growth of human thyroid carcinoma cells [51] . Strong inhibition of expression of thymosin b-10 seen in our experiment (Table 4 ) may in part explain the antiproliferative activity of 1,25-(OH) 2 D 3 .

At 15 min, 1,25-(OH) 2 D 3 stimulated a strong up-regulation (more than 2-fold) of variant regions of multiple immunoglobulin chains (data not shown).

At 3 h, 1,25-(OH) 2 D 3 caused a strong increase in the expression of interferon-c inducing factor isoform a precursor (IGIF) or IL-18, increased expression of CD59 glycoprotein precursor and chemokine CX3C (Table 5) .

Both IL-18 and its precursor were strongly induced by 1,25-(OH) 2 D 3 (2.7-fold at 3 h and 2.1-fold at 6 h, respectively) ( Table 5 ). IL-18 is a pleiotropic factor that shares structural features with IL-1 and functional activities with IL-12 and it was identified due to ability to induce interferon-gamma (IFN-c) production by T cells. The regulation of IL-18 by vitamin D receptor (VDR) or by its ligand, 1,25-(OH) 2 D 3 , was proposed recently [52] . IL-18 was shown to inhibit the osteoclastogenesis in vitro via granulocyte-macrophage colony-stimulating factor (GM-CSF) [53] and to possess the antiangiogenic and antitumor effects, particularly in combination with IL-2 or IL-12 [54] .

The expression of CD59 (protectin) precursor, a phosphatidylinositol-anchored glycoprotein, was increased by 1,25-(OH) 2 D 3 (1.5-fold, 3 h) ( Table 5 ). CD59 is a member of cell membrane-bound complement regulatory proteins that inhibit the formation of the terminal membrane attack complex of complement. Loss of CD59 offers a selective advantage for breast cancers, resulting in more aggressive tumors and conferring a poor prognosis for patients [55] . Increased expression of CD59 precursor after 1,25-(OH) 2 D 3 treatment, shown in our experiment, may contribute to the known anticancer activity of 1,25-(OH) 2 D 3 ( Table 5 ) [1] .

At 3 h, 1,25-(OH) 2 D 3 strongly inhibited the expression of interferon-c receptor (2.5-fold) and decreased the expression of MHC class II-associated invariant chain and proteasome subunit RCX (Table 5) . Interferon gamma receptor (IFN-cR) is the prototypical Th-1 cytokine (Table 5) . IFNc receptor knockout non-obese diabetic (NOD) mice showed a marked inhibition of insulitis and no signs of diabetes, in a sharp contrast to IFN-c knockout NOD mice, that had no change in insulitis and only a relatively minor delay in diabetes [56] . IFN-c receptor null mice also displayed a reduced pain-related behavior after nerve injury compared to wild-type mice. The decreased expression of MHC class II-associated invariant chain we observed (Table 5 ) is in agreement with previous findings on sig- Pancreatitis-associated protein precursor (PAP) À2.0 nificant reductions in class II MHC expression in monocytes exposed to 1,25-(OH) 2 D 3 [53] . Expression of small inducible cytokine a20 precursor (CCL20) and cytokine IL-15 was strongly increased by 1,25-(OH) 2 D 3 at 6 h (Table 5 ). CCL20 or macrophage inflammatory protein-3a is a CC-type chemokine that highly specifically binds to and activates CC chemokine receptor-6 (CCR6) and acts as a chemoattractant for memory/differentiated T-cells, B-cells, and immature dendritic cells, and possesses the antibacterial activity of a greater potency than b-defensins -1 and -2 [57] . IL-15 is expressed in multiple tissues and was able to induce the proliferation of activated T cells. It plays an important role in the development of memory CD8+ T cells and natural killer (NK) cells, for which IL-15 serves as the ''fuel.'' NK cells spontaneously kill tumor cell lines in vitro. IL-15 null mice displayed reduced numbers of T cells and lack of NK and NK T cells proving that this cytokine is crucial for murine NK cell development and suggesting the potential use of IL-15 therapy for expansion of NK cells in patients [58] .

Interestingly at 6 h, 1,25-(OH) 2 D 3 strongly inhibited the expression of both pancreatitis-associated protein (PAP) and its precursor (2.9-and 2.0-fold, respectively) ( Table 5 ). PAPs are the products of the C-type lectin supergene family of sugar-binding proteins with ability to agglutinate cells, to form antigen-antibody like precipitates with glycoconjugates and to induce mitosis in cells, which are normally not dividing. PAP is considered to be a stress protein and is constitutively expressed in the epithelial cells of the small intestine. PAP has been described as a marker protein for pancreatitis and for cystic fibrosis in neonates. Overexpression of PAP in human pancreatic ductal adenocarcinoma indicates tumor aggressiveness [59] .

The maximum change in the expression of these enzymes was observed at 3 h (Table 6 ). 1,25-(OH) 2 D 3 stimulated the increased expression of all three types of trypsin precursors (I and II-both anionic forms and III-cationic form) (Table 6) . Trypsinogens, the precursors to the serine protease trypsin, are found in the pancreas and mediate the digestive proteolysis in the small intestine. Anionic and cationic trypsinogens are approximately 90% identical in their primary structure. Expression of both pancreatic trypsin inhibitor type I and II (PSTI-I and PSTI-II) was increased 1.7-and 2.5-fold, respectively (Table 6 ). Both PSTI-I and II are expressed in pancreas, liver, and small intestine.

Expression of caspase 3 was strongly inhibited by 1,25-(OH) 2 D 3 (2.2-fold) and this was seen with different probe sets (Table 6) . Caspase 3 cleaves a variety of important cellular proteins and is considered to be a primary executioner of apoptosis or programmed cell death that can be initiated by a number of stimuli. Studies in caspase-3 null mice showed that this protease is essential for brain development [60] .

1,25-(OH) 2 D 3 strongly suppressed the expression of angiotensin-converting enzymes: CD13/aminopeptidase N (3.6-fold, Table 6 ) and kininase II or angiotensin Iconverting enzyme (ACE) (3.5-fold, Table 6 ).

CD13/aminopeptidase N (CD13/APN) is a type II membrane-bound metalloprotease that is expressed on the endothelial cells of angiogenic, but not normal, vasculature. It is essential for later stages of neovascular formation and is an important angiogenic activator, indicating that CD13/APN plays a functional role in tumorigenesis [61] . The cell surface aminopeptidase N is overexpressed in tumor cells. It is now generally agreed that conversion (degradation) of ANG III that causes high blood pressure to the hexapeptide ANG IV is aminopeptidase N dependent [62] . Intestine brush-border cells present a high concentration of aminopeptidase N that plays a role in the final digestion of peptides generated from hydrolysis of proteins by gastric and pancreatic proteases. Human CD13/APN is the receptor for coronaviruses; thus, its inhibitors might protect again SARS [63] . Our data are in concert with previous finding on reduction of cell surface CD13/ APN expression in the phagocytic cells by 1,25-(OH) 2 D 3 [64] and suggest 1,25-(OH) 2 D 3 as the potential inhibitor Table 6 1,25-(OH) 2 D 3 stimulated differential expression at 3 h of proteases, their inhibitors, and peptidases genes Interleukin-1b-converting enzyme-related protease CPP32 (caspase 3) À2.3 U49930 a ICE-like cysteine protease (Lice) or caspase 3 À2.2 a These genes also showed up-or down-regulation with other probe sets derived from different GenBank Accession numbers of the same protein.

of CD13/APN expression. Interestingly, in our experiment 1,25-(OH) 2 D 3 simultaneously increased the expression of transcription factor c-Maf (Table 5) , which was shown to suppress the CD13/APN expression (85-95% reduction) in human immature myeloblastic cells [65] . This could be the explanation for 1,25-(OH) 2 D 3 stimulated down-regulation of CD13/APN expression observed in our case. Angiotensin I-converting enzyme (ACE) plays a central role in the renin-angiotensin system. ACE is a carboxypeptidase that hydrolyzes the amino acid peptide angiotensin I into the potent vasoconstrictor angiotensin II. It was reported that angiotensin II stimulates angiogenesis in vivo, and angiotensin-converting enzyme (ACE) inhibitors block angiogenesis [66] . In addition to inducing vasoconstriction, angiotensin II increases blood pressure by a variety of physiological actions, including renal salt and water retention. ACE may affect blood pressure through the production of the vasoconstrictor angiotensin II and the inactivation of the vasodilator bradykinin. ACE inhibitors block the formation of angiotensin II and have been used to treat hypertension and heart failure [67] . ACE null mice have low blood pressure and the inability to concentrate urine [68] . Further, it has been reported that vitamin D 3 supplementation reduces blood pressure in patients with essential hypertension [69] , which may be in part due to its ability to down-regulate ACE.

Vitamin D is known not only as the principal regulator of bone development and calcium homeostasis but also for its widespread effects in many other tissues in vivo. It is a potent inhibitor of the variety of cancer cells [1] . 1,25-(OH) 2 D 3 exerts its effects through binding to a high affinity receptor (VDR) that following heterodimerization with the retinoid X receptor acts as a ligand-activated transcription factor and selectively binds to recognition sequences in the promoter region of target genes (vitamin D response elements-VDREs), transactivating or transrepressing a large variety of genes [2] [3] [4] . It is known that among the vast majority of genes that are regulated by vitamin D or its analogs only relatively few contain VDRE in their promoter region and thus are under direct transcriptional control of 1,25-(OH) 2 D 3 . For the majority of genes that do not have the VDREs in their promoter sequences, their regulation is indirect but important response to the cascade of reactions induced by 1,25-(OH) 2 D 3 . In this study, the microarray technique allowed us to investigate the role of vitamin D in intestinal absorption of calcium and other nutrients and provide the comprehensive knowledge of genes regulated by vitamin D in intestine, thus clarifying its widespread effects in vivo. In vitro micro-array technique was applied to study the 1,25-(OH) 2 D 3 stimulated gene expression in several cell lines: in mouse osteoblasts [70] , in squamous carcinoma cells [71] , and human colon carcinoma cells [72] . Although there is some similarity in regulation of expression of some genes by 1,25-(OH) 2 D 3 in our system and the squamous carcinoma and human colon carcinoma cells [71, 72] (in strong up-regulation of CYP24, in up-regulation of calmodulin, and in some other genes not presented in this paper), our studies were done in vivo in highly differentiated tissue that is responsible for nutrient absorption. We do not expect the same pattern of gene expression in immortal cell lines treated with high and unphysiological concentrations of 1,25-(OH) 2 D 3 as we see in vivo in a functional tissue carrying out intestinal absorption.

The most interesting for us was to identify 1,25-(OH) 2 D 3 regulated genes involved in Ca 2+ homeostasis and also genes involved in nutrient absorption in general. Our microarray and Q-PCR data showed the increase in the expression level of calcium homeostasis genes, and the differential expression of transporters and channels starting at 1 h after 1,25-(OH) 2 D 3 treatment with the expression maximum fold increase at 3 and 6 h (Tables 2 and 3 ). Our data confirm previously published data that 1,25-(OH) 2 D 3 up-regulates expression of transcellular calcium transport genes such as calbindin D 9k , plasma membrane Ca 2+ ATPase, epithelial calcium channels, TRPV5, and TRPV6 (Table 2 and Fig. 1 ) [1, 4, 7, 8, [12] [13] [14] [15] .

Molecules cross the intestinal epithelium into the systemic circulation primarily by three pathways: passive diffusion across the cell membranes (transcellular pathway), passive diffusion between adjacent cells (paracellular pathway), or carrier-mediated transport (carrier-mediated transcellular pathway). Lipophilic molecules easily cross the cell membrane via transcellular diffusion. Hydrophilic molecules, if not recognized by a carrier, traverse the epithelial barrier via the paracellular pathway, which is severely restricted by the presence of tight junctions. Historically, a simplified view of this absorptive process was that transcellular movement of nutrients and water via specific pumps, transporters, and channels would account for absorption, while an impermeable tight junction seal adjoining epithelial cells for the requisite barrier function. It has now become clear that transjunctional solute movement occurs in a regulated fashion, and that its regulation may be coupled to transcellular absorptive events. Thus epithelial solute transport and tight junction barrier function have to be viewed as related coordinated events [73] .

Tight junctions (TJ) are the contact points between the apical and basolateral membranes that limit paracel-lular flux of ions, proteins, and other macromolecules (gate function) and serve as a fence between the apical and basolateral plasma membranes in epithelial cells (fence function). The restricted movement of ions across the tight junctions gives rise to transepithelial electrical resistance (TER), which is often used as an index of the integrity of the tight junctions in an epithelial or endothelial tissue. Tight junctions are crucial in maintaining the polarized phenotype and the vectorial transport functions of epithelial cells. They are also a specialized membrane microdomains that might function as a molecular platform involved in cell signaling, vesicle protein docking, actin organization, and cell polarity in epithelial cells. Each particle of the tight junction strand is composed of transmembrane proteins and cytoplasmic plaque proteins connected to the actin cytoskeleton. Occludins, claudins, and the junctional adhesion molecule (JAM) are the three classes of transmembrane proteins localized to the tight junctions. TJs also concentrate tumor suppressor proteins and cell polarity proteins [34] .

Tight junctions have a highly dynamic structure so their permeability, assembly or disassembly can be regulated by a variety of cellular and metabolic mediators including cytokines [74] . The controlled and reversible opening of the tight junctions could be an attractive approach to increase the absorption of hydrophilic drugs across the intestinal epithelium because it could be applied to many different hydrophilic drugs and thus would allow them to escape the degradation by intracellular enzymes.

We identified the number of 1,25-(OH) 2 D 3 target genes with the maximum fold change in the expression at 3 h that are (co)localized in the proximal vicinity of intestinal tight junctions and thus could influence their integrity and permeability. These include not only transporters and channels, but also several intra-/intercellular matrix related genes and G-proteins (Tables 3 and 4 ). The following genes down-regulated by 1,25-(OH) 2 D 3 are involved in the regulation of tight junctions permeability (Tables 3 and 4) .

Sodium/potassium ATPase. Both a-1 and b-1 subunits of sodium/potassium ATPase (Na,K-ATPase) were significantly down-regulated by 1,25-(OH) 2 D 3 (Table 3) . Na,K-ATPase catalyzes an ATP-dependent transport of three Na + ions out and two K + ions into the cell per pump cycle. In epithelial cells, Na,K-ATPase was also involved in the formation of tight junctions through RhoA GTPase and stress fibers. The inhibition of Na,K-ATPase in tight monolayers of epithelial cells resulted in an increased permeability of tight junctions to ions and non-ionic molecules [34] .

Claudin-3. We observed a decreased expression of claudin-3 (2.2-fold, Table 4 ). The claudins, tight junction-specific adhesion molecules, create paracellular channels and their first extracellular domain is sufficient to determine both paracellular charge selectivity and transepithelial electrical resistance (TER). The tight junctions charge selectivity towards cations or anions in epithelial cells could be reversed by expression of different claudins [75] . The selective loss of claudin-3 was the cause for ''leakiness'' of blood-brain barrierÕs tight junctions at experimental autoimmune encephalomyelitis [48] .

Water channel aquaporin 8. Aquaporin 8 is the tight junction channel; was down-regulated by 1,25-(OH) 2 D 3 (Table 3) . Thus, the inhibition of its expression might change the TJ selectivity towards cations.

Cadherin-17. It plays a very important role in cell-tocell adhesion and was down-regulated 2.6-fold by 1,25-(OH) 2 D 3 (Table 4 ) [44] .

RhoA. It is the small GTP-ase that regulates remodeling of the actin cytoskeleton during cell morphogenesis and motility. It was shown that RhoA GTPase is an essential component downstream of Na,K-ATPasemediated regulation of tight junctions [34] .

Thus, 1,25-(OH) 2 D 3 may increase intestinal epithelial tight junction permeability or modulate their selectivity towards Ca 2+ and other cations by regulation of expression of proteins structurally involved in tight junction formation. The increased tight junction permeability and/or selectivity, regulated by 1,25-(OH) 2 D 3 , could route Ca 2+ absorption through the tight junction-regulated paracellular pathway in the intestinal epithelia. Our proposal is in agreement with published data on the 1,25-(OH) 2 D 3 stimulated increase of tight junction conductance and increased paracellular Ca 2+ , Na + , Rb + , and mannitol transport in enterocyte-like cell line Caco-2, while no significant contribution of the Ca 2+ -ATPase-mediated transcellular pathway to overall transepithelial Ca 2+ transport was detected [76] . The evidence has accumulated since late 80th for in vitro intestinal model for Ca 2+ and Pi transport [77] and recently for Ca 2+ transport in Caco-2 cells [78] that 1,25-(OH) 2 D 3 enhanced both cell-mediated active and passive paracellular ion movement.

So based on our microarray data we propose that 1,25-(OH) 2 D 3 regulates the intestinal absorption of Ca 2+ in vivo through both transcellular and paracellular pathways by the stimulation or suppression of the expression the group of genes and also identified 1,25-(OH) 2 D 3 target genes possibly involved in regulation of tight junctionÕs permeability and/or selectivity.

At 3 h, 1,25-(OH) 2 D 3 caused more down-regulation of transporters mRNA and channels genes than up-regulation (Table 3) . It was the time of a maximal decrease in the expression for both a-1 and b-1 subunits of sodium/potassium ATPase (Na,K-ATPase), voltage dependent anion channel 1, transmembrane protein rnp21.4, claudin-3, urate transporter/channel, hypothetical ATPase, and aquaporin-8 after the 1,25-(OH) 2 D 3 treatment (Table 3) . At the same time, 1,25-(OH) 2 D 3 increased the expression of p28 and Sip9 (syncollin) involved in vesicle transport of proteins (endocytosis and exocytosis) though some of their functions remain unknown ( Table 3) .

The expression of both organic cation transporter (OCT1) responsible for transport of sugar and organic cations like drugs and xenobiotics and expression of multidrug transporter 1, ATP-binding cassette (ABCB) transporter, which is involved in the transport of peptide antigens for association with MHC class I molecules was significantly increased at 6 h post-1,25-(OH) 2 D 3 treatment (Table 3 ). At the same time, the expression of Nramp2-transporter for other than Ca 2+ divalent cations, like iron, was significantly decreased (Table 3) .

The actions of 1,25-(OH) 2 D 3 on the immune system are thought to be targeted primarily at T helper (Th) cells. Our data showed that transcription factor c-maf which is involved in the regulation of the Th1/Th2 balance was up-regulated by 1,25-(OH) 2 D 3 (see Table 5 ). Our data are in agreement with published data on 1,25-(OH) 2 D 3 acting directly on Th cells, enhancing the development of a Th2 phenotype and augmenting the expression of the transcription factor c-maf [79] .

Some of the 1,25-(OH) 2 D 3 target genes that we identified have also antitumor effect. These include IL-18 (and its precursor), IL-15, chemokines CX3C and CC chemokine ST38 precursor and CD59 precursor ( Table  5) . Murine models have shown that IL-18 has both antiangiogenic and antitumor effects particularly in combination with IL-2 or IL-12 and also induces the high levels of IFN-c secretion [54] . In a new gene therapy approach for treatment of osteosarcoma, recombinant IL-18 was proposed for use in combination with plasmids encoding the IL-12 [80] . IL-15 is involved in the development of natural killer (NK) cells that can spontaneously kill tumor cell lines in vitro [81] . The loss of CD59 offered a selective advantage for development of breast cancers, resulting in more aggressive tumors [55] . At the same time overexpression of PAP (its expression was suppressed by 1,25-(OH) 2 D 3 in our experiment) in human pancreatic ductal adenocarcinoma indicated tumor aggressiveness [59] .

The agents that interfere with blood vessel formation can be used to block tumor progression. Our data showed that 1,25-(OH) 2 D 3 strongly down-regulated enzymes involved in angiogenesis which is mandatory for tumor progression: CD13/aminopeptidase N and angiotensin I-converting enzyme (ACE) or kininase II. The cell surface aminopeptidase N (APN/CD13) was overex-pressed in tumor cells. However, potent, selective, and, particularly, non-cytotoxic inhibitors of this protein were lacking. Angiotensin I-converting enzyme (ACE) hydrolyzed the amino acid peptide angiotensin I into the potent vasoconstrictor angiotensin II. It was reported that angiotensin II stimulates angiogenesis in vivo, and ACE inhibitors inhibit angiogenesis [66] . Our data suggest 1,25-(OH) 2 D 3 as the potential inhibitor of both CD13/APN and ACE expression which also can contribute to 1,25-(OH) 2 D 3 anticancer activity. Further, this could account for the well-known effect of vitamin D on the reduction of hypertension [69] .

Microarray data have provided major insight into gene transcription profiles in rat intestine in response to 1,25-(OH) 2 D 3 thus creating a snapshot of molecular events following secosteroid intervention. We proposed that 1,25-(OH) 2 D 3 regulates not only established transcellular calcium absorption but also paracellular calcium transport as well. We showed that 1,25-(OH) 2 D 3 modulated the expression of different classes of genes in rat intestine, not only those directly involved in the absorption of nutrients in small intestine but also genes involved in immune response and angiogenesis. Since many genes may not have a VDRE in the promoter region, their regulation by 1,25-(OH) 2 D 3 could be indirect through other proteins/factors expressed early in response to 1,25-(OH) 2 D 3 or through increased intracellular Ca 2+ concentration.

In addition to its central role in the maintenance of extracellular calcium level and bone mineralization, 1,25-(OH) 2 D 3 also acts as a modulator of cell growth and differentiation in a number of cell types, including breast cancer cells. Particularly important to us was to discover possible biochemical grounds for anti-proliferative and anticancer effects of 1,25-(OH) 2 D 3 by induction of expression IL-15, IL-18, CD59 (protectin), CX3C chemokine, and inhibition of the expression of thymosin-b-10 and both angiogenesis promoting enzymes CD13/APN and ACE. The down-regulation of ACE may also account in part for the anti-hypertensive actions of vitamin D. These data may help to extend the potential use of 1,25-(OH) 2 D 3 and its analogs in the treatment or prevention of a number of diseases.

Smith, Wendy Hellwig, Maggie Highland, and Margaret Clagett-Dame from the Biochemistry Department, UW-Madison for their assistance and valuable advice with this project and Pat Mings from the Biochemistry Department, UW-Madison for her help with manuscript preparation.

Methods in Enzymology

We cordially thank Wayne Davis and Sandra Splinter BonDurant from the Gene Expression Center at the Biotechnology Center of UW-Madison, Christina Gutierrez and Chiara Cirelli from the Psychiatry Institute at UW-Madison, Stan Trask from Affymetrix, Connie