Search

In many species, mate choice continues after the mating via female- or egg-derived biochemical factors that induce selective changes in sperm pre-fertilization physiology and behaviour. Recent studies have indicated that gamete-mediated mate choice likely occurs also in humans, but the mechanistic basis of the process has remained virtually unexplored. Here, we investigated whether female-induced modifications in sperm protein SUMOylation (post-translational modification of the proteome) could serve as a novel mechanism for gamete-mediated mate choice in humans. We treated the sperm of ten males with the oocyte-surrounding bioactive liquid (follicular fluid) of five females and investigated motility, viability and global protein SUMOylation status of the sperm in all (n = 50) of these male-female combinations (full-factorial design). All the measured sperm traits were affected by male-female combinations, and sperm protein SUMOylation status was also negatively associated with sperm motility. Furthermore, our results indicate that female-induced sperm protein SUMOylation is selective, potentially allowing females to increase sperm motility in some males, whereas decreasing it in the others. Consequently, our findings suggest that follicular fluid may non-randomly modify the structure and function of sperm proteome and in this way facilitate gamete-mediated mate choice in humans and possibly many other species. However, due to the relatively low number of female subjects and their potential infertility problems, our results should be replicated with larger subset of fully fertile women.

Seminal plasma (SP) plays a crucial role in reproduction and contains a large number of proteins, many of which may potentially modify sperm functionality. To evaluate the effects of SP identity and its protein composition on human sperm function, we treated the sperm of several males with either their own or multiple foreign SPs in all possible sperm–SP combinations (full-factorial design). Then we recorded sperm motility and viability in these combinations and investigated whether the sperm performance is dependent on sperm and SP identity (or their interaction). Finally, we studied whether the above-mentioned sperm traits are affected by the abundance of three SP proteins, dipeptidyl peptidase-4 (DPP4), neutral endopeptidase (NEP), and aminopeptidase N (APN). The identity of the SP donor affected sperm swimming velocity, viability, and the proportion of hyperactivated sperm, but males’ own SP was not consistently more beneficial for sperm than foreign SPs. Furthermore, we show that sperm performance is also partly affected by the interaction between sperm and SP donor. Finally, we found that DPP4 and NEP levels in SP were positively associated with sperm swimming velocity and hyperactivation. Taken together, our results highlight the importance of seminal plasma as a potential source of biomarkers for diagnostics and therapeutic interventions for male-derived infertility.

Infertility is assumed to arise exclusively from male- and female-dependent pathological factors. However, recent studies have indicated that reproductive failure may also result from the reproductive incompatibility of the partners. Selection against such incompatibilities likely occurs via female-derived reproductive secretions, including follicular fluid (FF), that mediate gamete-level mate choice towards the sperm of specific males. To facilitate potential development of diagnostic tests for human reproductive incompatibility, we examined whether sperm physiological response to female serum indicate male–female compatibility in the presence of FF. We performed a full-factorial experiment, in which the sperm of 10 males were treated with the FF and serum of 6 healthy females. We found that sperm motility and viability in both biofluids were highly similar and that in 70% of the males, sperm serum treatment predicted male–female compatibility. We also identified male human leucocyte antigen (HLA) alleles and female (FF and serum) anti-HLA antibodies and tested whether the number of allele–antibody matches predict sperm physiological response to female fluids. However, no association was found between measured sperm traits and the number of allele–antibody matches. Overall, the present results may open novel possibilities for the future development of reproductive incompatibility tests and may pave the way towards more accurate infertility diagnostics and treatments.

Glucocorticoid (GC) receptor (GR) is a key transcription factor (TF) that regulates vital metabolic and antiinflammatory processes. We have identified BCL6 corepressor (BCOR) as a dexamethasone-stimulated interaction partner of GR. BCOR is a component of non-canonical polycomb repressor complex 1.1 (ncPCR1.1) and linked to different developmental disorders and cancers, but the role of BCOR in GC signaling is poorly characterized. Here, using ChIP-seq we show that, GC induces genome-wide redistribution of BCOR chromatin binding towards GR-occupied enhancers in HEK293 cells. As assessed by RNA-seq, depletion of BCOR altered the expression of hundreds of GC-regulated genes, especially the ones linked to TNF signaling, GR signaling and cell migration pathways. Biotinylation-based proximity mapping revealed that GR and BCOR share several interacting partners, including nuclear receptor corepressor NCOR1. ChIP-seq showed that the NCOR1 co-occurs with both BCOR and GR on a subset of enhancers upon GC treatment. Simultaneous depletion of BCOR and NCOR1 influenced GR target gene expression in a combinatorial and gene-specific manner. Finally, we show using live cell imaging that the depletion of BCOR together with NCOR1 markedly enhances cell migration. Collectively, our data suggest BCOR as an important gene and pathway selective coregulator of GR transcriptional activity.

Hypoxia-inducible factor (HIF), an αβ dimer, is the master regulator of oxygen homeostasis with hundreds of hypoxia-inducible target genes. Three HIF isoforms differing in the oxygen-sensitive α subunit exist in vertebrates. While HIF-1 and HIF-2 are known transcription activators, HIF-3 has been considered a negative regulator of the hypoxia response pathway. However, the humanHIF3AmRNA is subject to complex alternative splicing. It was recently shown that the long HIF-3α variants can form αβ dimers that possess transactivation capacity. Here, we show that overexpression of the long HIF-3α2 variant induces the expression of a subset of genes, including the erythropoietin (EPO) gene, while simultaneous downregulation of all HIF-3α variants by siRNA targeting a sharedHIF3Aregion leads to downregulation ofEPOand additional genes. EPO mRNA and protein levels correlated withHIF3Asilencing and HIF-3α2 overexpression. Chromatin immunoprecipitation analyses showed that HIF-3α2 binding associated with canonical hypoxia response elements in the promoter regions ofEPO. Luciferase reporter assays showed that the identified HIF-3α2 chromatin-binding regions were sufficient to promote transcription by all three HIF-α isoforms. Based on these data, HIF-3α2 is a transcription activator that directly regulatesEPOexpression.

Article, Contrary to cell cycle–associated cyclin-dependent kinases, CDK5 is best known for its regulation of signaling processes in differentiated cells and its destructive activation in Alzheimer's disease. Recently, CDK5 has been implicated in a number of different cancers, but how it is able to stimulate cancer-related signaling pathways remains enigmatic. Our goal was to study the cancer-promoting mechanisms of CDK5 in prostate cancer. We observed that CDK5 is necessary for proliferation of several prostate cancer cell lines. Correspondingly, there was considerable growth promotion when CDK5 was overexpressed. When examining the reasons for the altered proliferation effects, we observed that CDK5 phosphorylates S308 on the androgen receptor (AR), resulting in its stabilization and differential expression of AR target genes including several growth-priming transcription factors. However, the amplified cell growth was found to be separated from AR signaling, further corroborated by CDK5-depdent proliferation of AR null cells. Instead, we found that the key growth-promoting effect was due to specific CDK5-mediated AKT activation. Down-regulation of CDK5 repressed AKT phosphorylation by altering its intracellular localization, immediately followed by prominent cell cycle inhibition. Taken together, these results suggest that CDK5 acts as a crucial signaling hub in prostate cancer cells by controlling androgen responses through AR, maintaining and accelerating cell proliferation through AKT activation, and releasing cell cycle breaks., Publisher version, http://purl.org/eprint/status/PeerReviewed

Androgen receptor (AR) is a ligand-activated transcription factor that plays a central role in the development and growth of prostate carcinoma. PIAS1 is an AR- and SUMO-interacting protein and a putative transcriptional coregulator overexpressed in prostate cancer. To study the importance of PIAS1 for the androgen-regulated transcriptome of VCaP prostate cancer cells, we silenced its expression by RNAi. Transcriptome analyses revealed that a subset of the AR-regulated genes is significantly influenced, either activated or repressed, by PIAS1 depletion. Interestingly, PIAS1 depletion also exposed a new set of genes to androgen regulation, suggesting that PIAS1 can mask distinct genomic loci from AR access. In keeping with gene expression data, silencing of PIAS1 attenuated VCaP cell proliferation. ChIP-seq analyses showed that PIAS1 interacts with AR at chromatin sites harboring also SUMO2/3 and surrounded by H3K4me2; androgen exposure increased the number of PIAS1-occupying sites, resulting in nearly complete overlap with AR chromatin binding events. PIAS1 interacted also with the pioneer factor FOXA1. Of note, PIAS1 depletion affected AR chromatin occupancy at binding sites enriched for HOXD13 and GATA motifs. Taken together, PIAS1 is a genuine chromatin-bound AR coregulator that functions in a target gene selective fashion to regulate prostate cancer cell growth.

Androgen receptor (AR) plays an important regulatory role in prostate cancer. AR's transcriptional activity is regulated by androgenic ligands, but also by post-translational modifications, such as SUMOylation. To study the role of AR SUMOylation in genuine chromatin environment, we compared androgen-regulated gene expression and AR chromatin occupancy in PC-3 prostate cancer cell lines stably expressing wild-type (wt) or doubly SUMOylation site-mutated AR (AR-K386R,K520R). Our genome-wide gene expression analyses reveal that the SUMOylation modulates the AR function in a target gene and pathway selective manner. The transcripts that are differentially regulated by androgen and SUMOylation are linked to cellular movement, cell death, cellular proliferation, cellular development and cell cycle. Fittingly, SUMOylation mutant AR cells proliferate faster and are more sensitive to apoptosis. Moreover, ChIP-seq analyses show that the SUMOylation can modulate the chromatin occupancy of AR on many loci in a fashion that parallels their differential androgen-regulated expression. De novo motif analyses reveal that FOXA1, C/EBP and AP-1 motifs are differentially enriched at the wtAR- and the AR-K386R,K520R-preferred genomic binding positions. Taken together, our data indicate that SUMOylation does not simply repress the AR activity, but it regulates AR's interaction with the chromatin and the receptor's target gene selection.

Article, Androgen receptor (AR) is a male sex steroid-activated transcription factor (TF) that plays a critical role in prostate cancers, including castration-resistant prostate cancers (CRPC) that typically express amplified levels of the AR. CRPC-derived VCaP cells display an excessive number of chromatin AR-binding sites (ARBs) most of which localize to distal inter- or intragenic regions. Here, we analyzed direct transcription programs of the AR in VCaP cells using global nuclear run-on sequencing (GRO-seq) and integrated the GRO-seq data with the ARB and VCaP cell-specific TF-binding data. Androgen immediately activated transcription of hundreds of protein-coding genes, including IGF-1 receptor and EGF receptor. Androgen also simultaneously repressed transcription of a large number of genes, including MYC. As functional enhancers have been postulated to produce enhancer-templated non-coding RNAs (eRNAs), we also analyzed the eRNAs, which revealed that only a fraction of the ARBs reside at functional enhancers. Activation of these enhancers was most pronounced at the sites that also bound PIAS1, ERG and HDAC3, whereas binding of HDAC3 and PIAS1 decreased at androgen-repressed enhancers. In summary, our genome-wide data of androgen-regulated enhancers and primary target genes provide new insights how the AR can directly regulate cellular growth and control signaling pathways in CPRC cells, published version, peerReviewed

Article, Inflammatory processes and androgen signaling are critical for the growth of prostate cancer (PC), the most common cancer among males in Western countries. To understand the importance of potential interplay between pro-inflammatory and androgen signaling for gene regulation, we have interrogated the crosstalk between androgen receptor (AR) and NF-κB, a key transcriptional mediator of inflammatory responses, by utilizing genome-wide chromatin immunoprecipitation sequencing and global run-on sequencing in PC cells. Co-stimulation of LNCaP cells with androgen and pro-inflammatory cytokine TNFα invoked a transcriptome which was very distinct from that induced by either stimulation alone. The altered transcriptome that included gene programs linked to cell migration and invasiveness was orchestrated by significant remodeling of NF-κB and AR cistrome and enhancer landscape. Although androgen multiplied the NF-κB cistrome and TNFα restrained the AR cistrome, there was no general reciprocal tethering of the AR to the NF-κB on chromatin. Instead, redistribution of FOXA1, PIAS1 and PIAS2 contributed to the exposure of latent NF-κB chromatin-binding sites and masking of AR chromatin-binding sites. Taken together, concomitant androgen and pro-inflammatory signaling significantly remodels especially the NF-κB cistrome, reprogramming the PC cell transcriptome in fashion that may contribute to the progression of PC., published version, peerReviewed

The nuclear receptor vitamin D receptor (VDR) is known to associate with two vitamin D response element (VDRE) containing chromatin regions of the insulin-like growth factor binding protein 3 (IGFBP3) gene. In non-malignant MCF-10A human mammary cells, we show that the natural VDR ligand 1α,25-dihydroxyvitamin D(3) (1α,25(OH)(2)D(3)) causes cyclical IGFBP3 mRNA accumulation with a periodicity of 60 min, while in the presence of the potent VDR agonist Gemini the mRNA is continuously accumulated. Accordingly, VDR also showed cyclical ligand-dependent association with the chromatin regions of both VDREs. Histone deacetylases (HDACs) play an important role both in VDR signalling and in transcriptional cycling. From the 11 HDAC gene family members, only HDAC4 and HDAC6 are up-regulated in a cyclical fashion in response to 1α,25(OH)(2)D(3), while even these two genes do not respond to Gemini. Interestingly, HDAC4 and HDAC6 proteins show cyclical VDR ligand-induced association with both VDRE regions of the IGFBP3 gene, which coincides with histone H4 deacetylation on these regions. Moreover, combined silencing of HDAC4 and HDAC6 abolishes the cycling of the IGFBP3 gene. We assume that due to more efficient VDR interaction, Gemini induces longer lasting chromatin activation and therefore no transcriptional cycling but monotonically increasing IGFBP3 mRNA. In conclusion, 1α,25(OH)(2)D(3) regulates IGFBP3 transcription through short-term cyclical association of VDR, HDAC4 and HDAC6 to both VDRE-containing chromatin regions.

Primary 1alpha,25-dihydroxyvitamin D(3) (1alpha,25(OH)(2)D(3))-responding genes are controlled by the vitamin D receptor (VDR) binding to specific sites (VDREs) that are located within the regulatory regions of these genes. According to previous studies, the gene encoding 25-dihydroxyvitamin D(3) 24-hydroxylase, CYP24, which is the strongest known 1alpha,25(OH)(2)D(3)-responsive gene, has multiple VDREs that locate within the proximal and the distal promoter. However, it has remained unclear, what is the biological role of these regions and how they participate in the regulation of transcription. In this study, we found a different CYP24 expression profile in normal (MCF-10A) and malignant (MCF-7) human mammary cells. Moreover, CYP24 mRNA showed to be three times more stable in MCF-7 cells than in MCF-10A cells. We studied the mechanism of this difference using expression profiling, quantitative chromatin immunoprecipitation and chromosome conformation capture assays. Interestingly, the number of functional VDREs was higher in MCF-7 cells than in MCF-10A cells. Three functional VDREs in MCF-7 cells are connected to linear mRNA accumulation, whereas only one VDRE seems to lead to stepwise CYP24 mRNA accumulation in MCF-10A cells. The distal VDREs were involved in transcriptional regulation via ligand-dependent, dynamic chromatin looping, which brings cyclically the distal elements together either individually or simultaneously next to the transcription start site. In conclusion, our data suggest that in comparison to normal cells, clearing of 1alpha,25(OH)(2)D(3) is enhanced in malignant cells due to differences in transcriptional regulation of CYP24 and metabolism of CYP24 mRNA.

The anti-proliferative effects of histone deacetylase (HDAC) inhibitors and 1alpha,25-dihydroxyvitamin D3 [1alpha,25(OH)2D3] converge via the interaction of un-liganded vitamin D receptor (VDR) with co-repressors recruiting multiprotein complexes containing HDACs and via the induction of cyclin-dependent kinase inhibitor (CDKI) genes of the INK4 and Cip/Kip family. We investigated the effects of the HDAC inhibitor Trichostatin A (TSA) and 1alpha,25(OH)2D3 on the proliferation and CDKI gene expression in malignant and non-malignant mammary epithelial cell lines. TSA induced the INK4-family genes p18 and p19, whereas the Cip/Kip family gene p21 was stimulated by 1alpha,25(OH)2D3. Chromatin immunoprecipitation and RNA inhibition assays showed that the co-repressor NCoR1 and some HDAC family members complexed un-liganded VDR and repressed the basal level of CDKI genes, but their role in regulating CDKI gene expression by TSA and 1alpha,25(OH)2D3 were contrary. HDAC3 and HDAC7 attenuated 1alpha,25(OH)2D3-dependent induction of the p21 gene, for which NCoR1 is essential. In contrast, TSA-mediated induction of the p18 gene was dependent on HDAC3 and HDAC4, but was opposed by NCoR1 and un-liganded VDR. This suggests that the attenuation of the response to TSA by NCoR1 or that to 1alpha,25(OH)2D3 by HDACs can be overcome by their combined application achieving maximal induction of anti-proliferative target genes.

Article, Background Cells have developed many ways to cope with external stress. One distinctive feature in acute proteotoxic stresses, such as heat shock (HS), is rapid post-translational modification of proteins by SUMOs (small ubiquitin-like modifier proteins; SUMOylation). While many of the SUMO targets are chromatin proteins, there is scarce information on chromatin binding of SUMOylated proteins in HS and the role of chromatin SUMOylation in the regulation of transcription. Results We mapped HS-induced genome-wide changes in chromatin occupancy of SUMO-2/3-modified proteins in K562 and VCaP cells using ChIP-seq. Chromatin SUMOylation was further correlated with HS-induced global changes in transcription using GRO-seq and RNA polymerase II (Pol2) ChIP-seq along with ENCODE data for K562 cells. HS induced a rapid and massive rearrangement of chromatin SUMOylation pattern: SUMOylation was gained at active promoters and enhancers associated with multiple transcription factors, including heat shock factor 1. Concomitant loss of SUMOylation occurred at inactive intergenic chromatin regions that were associated with CTCF-cohesin complex and SETDB1 methyltransferase complex. In addition, HS triggered a dynamic chromatin binding of SUMO ligase PIAS1, especially onto promoters. The HS-induced SUMOylation on chromatin was most notable at promoters of transcribed genes where it positively correlated with active transcription and Pol2 promoter-proximal pausing. Furthermore, silencing of SUMOylation machinery either by depletion of UBC9 or PIAS1 enhanced expression of HS-induced genes. Conclusions HS-triggered SUMOylation targets promoters and enhancers of actively transcribed genes where it restricts the transcriptional activity of the HS-induced genes. PIAS1-mediated promoter SUMOylation is likely to regulate Pol2-associated factors in HS., published version

Recently, insulin-like growth factor binding proteins (IGFBPs) have been found to be primary mediators of the anti-proliferative actions of the nuclear hormone 1alpha,25-dihydroxyvitamin D3 [1alpha,25(OH)2D3], but dependent on cellular context IGFBPs can also have a mitogenic effect. In this study, we performed expression profiling of all six human IGFBP genes in prostate and bone cancer cells and demonstrated that IGFBP1, 3 and 5 are primary 1alpha,25(OH)2D3 target genes. In silico screening of the 174 kb of genomic sequence surrounding all six IGFBP genes identified 15 candidate vitamin D response elements (VDREs) close to or in IGFBP1, 2, 3 and 5 but not in the IGFBP4 and 6 genes. The putative VDREs were evaluated in vitro by gelshift assays and in living cells by reporter gene and chromatin immuno-precipitation (ChIP) assays. Of these 10 VDREs appear to be functional. ChIP assays demonstrated for each of these an individual, stimulation time-dependent association profile not only with the vitamin D receptor, but also with first heterodimeric partner the retinoid X receptor, other regulatory complex components and phosphorylated RNA polymerase II. Some of the VDREs are located distantly from the transcription start sites of IGFBP1, 3 and 5, but all 10 VDREs seem to contribute to the regulation of the genes by 1alpha,25(OH)2D3. In conclusion, IGFBP1, 3 and 5 are primary 1alpha,25(OH)2D3 target genes that in intact cells are each under the control of multiple VDREs.

Aims G protein-coupled receptor 35 (GPR35) has been characterized to be one of the genes that are up-regulated in human heart failure. Since mechanisms controlling GPR35 expression are not known, we investigated the regulation of GPR35 gene and protein expression in cardiac myocytes and in the mouse models of cardiac failure. Methods and results In cardiac myocytes, GPR35 gene expression was found to be exceptionally sensitive to hypoxia and induced by hypoxia-inducible factor-1 (HIF-1) activation. HIF-1-dependent regulation was established by genetic (HIF-1/VP16, Inhibitory Per/Arnt/Sim domain protein) and chemical [desferrioxamine (DFO)] modulation of the HIF-1 pathway and further confirmed by mutation analysis of the GPR35 promoter and by demonstrating direct binding of endogenous HIF-1 to the gene promoter. Hypoxia increased the number and density of GPR35 receptors on the cardiomyocyte cell membranes. Chemical GPR35 agonist Zaprinast caused GPR35 activation and receptor internalization in cardiac myocytes. In addition, overexpressed GPR35 disrupted actin cytoskeleton arrangement and caused morphological changes in cultured cardiomyocytes. GPR35 gene and protein expressions were also induced in mouse models of cardiac failure; the acute phase of myocardial infarction and during the compensatory and decompensatory phase of pressure-load induced cardiac hypertrophy. Conclusions Cardiac expression of GPR35 is regulated by hypoxia through activation of HIF-1. The expression of GPR35 in mouse models of cardiac infarction and pressure load suggests that GPR35 could be used as an early marker of progressive cardiac failure.

Background The histone demethylase GASC1 (JMJD2C) is an epigenetic factor suspected of involvement in development of different cancers, including breast cancer. It is thought to be overexpressed in the more aggressive breast cancer types based on mRNA expression studies on cell lines and meta analysis of human breast cancer sets. This study aimed to evaluate the prognostic and predictive value of GASC1 for women with invasive breast cancer. Methods All the 355 cases were selected from a cohort enrolled in the Kuopio Breast Cancer Project between April 1990 and December 1995. The expression of GASC1 was studied by immunohistochemistry (IHC) on tissue microarrays. Additionally relative GASC1 mRNA expression was measured from available 57 cases. Results In our material, 56% of the cases were GASC1 negative and 44% positive in IHC staining. Women with GASC1 negative tumors had two years shorter breast cancer specific survival and time to relapse than the women with GASC1 positive tumors (p=0.017 and p=0.034 respectively). The majority of GASC1 negative tumors were ductal cases (72%) of higher histological grade (84% of grade II and III altogether). When we evaluated estrogen receptor negative and progesterone receptor negative cases separately, there was 2 times more GASC1 negative than GASC1 positive tumors in each group (chi2, p= 0.033 and 0.001 respectively). In the HER2 positive cases, there was 3 times more GASC1 negative cases than GASC1 positives (chi2, p= 0.029). Patients treated with radiotherapy (n=206) and hormonal treatment (n=62) had better breast cancer specific survival, when they were GASC1 positive (Cox regression: HR=0.49, p=0.007 and HR=0.33, p=0.015, respectively). The expression of GASC1 mRNA was in agreement with the protein analysis. Conclusions This study indicates that the GASC1 is both a prognostic and a predictive factor for women with invasive breast cancer. GASC1 negativity is associated with tumors of more aggressive histopathological types (ductal type, grade II and III, ER negative, PR negative). Patients with GASC1 positive tumors have better breast cancer specific survival and respond better to radiotherapy and hormonal treatment.

The candidate human tumor suppressor gene cyclin C is a primary target of the anti-proliferative hormone 1alpha,25-dihydroxyvitamin D3 [1alpha,25(OH)2D3], but binding sites for the 1alpha,25(OH)2D3 receptor (VDR), so-called 1alpha,25(OH)2D3 response elements (VDREs), have not yet been identified in the promoter of this gene. We screened various cancer cell lines by quantitative PCR and found that the 1alpha,25(OH)2D3 inducibility of cyclin C mRNA expression, in relationship with the 24-hydroxylase (CYP24) gene, was best in MCF-7 human breast cancer cells. To characterize the molecular mechanisms, we analyzed 8.4 kb of the cyclin C promoter by using chromatin immunoprecipitation assays (ChIP) with antibodies against acetylated histone 4, VDR and its partner receptor, retinoid X receptor (RXR). The histone 4 acetylation status of all 23 investigated regions of the cyclin C promoter did not change significantly in response to 1alpha,25(OH)2D3, but four independent promoter regions showed a consistent, 1alpha,25(OH)2D3-dependent association with VDR and RXR over a time period of 240 min. Combined in silico/in vitro screening identified in each of these promoter regions a VDRE and reporter gene assays confirmed their functionality. Moreover, re-ChIP assays monitored simultaneous association of VDR with RXR, coactivator, mediator and RNA polymerase II proteins on these regions. Since cyclin C protein is associated with those mediator complexes that display transcriptional repressive properties, this study contributes to the understanding of the downregulation of a number of secondary 1alpha,25(OH)2D3-responding genes.