Skip to main content

Genomic actions of 1,25-dihydroxyvitamin D3 on insulin receptor gene expression, insulin receptor number and insulin activity in the kidney, liver and adipose tissue of streptozotocin-induced diabetic rats

Abstract

Background

this study set out to examine the effects of the treatment with 1,25-dihydroxyvitamin D3 (1,25D3) [150 IU/Kg (3.75 μg/Kg) one a day, for 15 days] to non-diabetic rats and in rats rendered diabetic by a single injection of streptozotocin [65 mg/kg].

Results

treatment with 1,25D3 to non-diabetic rats did not affect the biochemical parameters measured in the plasma and urine of these animals. Likewise, insulin receptor expression in the kidney, liver, or adipose tissue and insulin-stimulated glucose transport in adipocytes from these animals were not affected either. Treatment with 1,25D3 to streptozotocin-induced diabetic rats did not correct the hyperglycemia, hypoinsulinemia, glycosuria or ketonemia induced by the diabetes, although it partially reversed the over-expression of the insulin receptor gene in the liver and adipose tissue, without altering the normal expression of this gene in the kidney. These effects were accompanied by a normalization of the number of insulin receptors without altering receptor affinity but improving the insulin response to glucose transport in adipocytes from these diabetic animals. Moreover, a computer search in the rat insulin receptor promoter revealed the existence of two candidate vitamin D response element (VDRE) sequences located at -256/-219 bp and -653/-620 bp, the first overlapped by three and the second by four AP-2-like sites.

Conclusion

these genomic actions of 1,25D3 could represent beneficial effects associated with the amelioration of diabetes via mechanisms that possibly involve direct transcriptional activation of the rat insulin receptor gene. The candidate VDREs identified may respond to 1,25D3 via activation of the vitamin D receptor, although this remains to be investigated.

Background

It is known that vitamin D and especially its activated metabolite 1,25-dihydroxyvitamin D3 (1,25D3), are involved in controlling the normal function of the endocrine pancreas, and particularly insulin secretion [1, 2]. Vitamin D deficiency inhibits rat pancreatic secretion and turnover of insulin, leading to impaired glucose tolerance, while replacement therapy with 1,25D3, is able to reverse these abnormalities [3, 4]. 1,25D3 also affects the insulin receptor (IR), the protein to which insulin must bind to carry out its multiple biological actions in the cells. In this respect, our group has reported the first demonstration that 1,25D3 increased human IR mRNA levels, the IR number, and the insulin response in U-937 human promonocytic cells through mechanisms that involve direct transcriptional activation of the human IR gene [57].

In non-obese diabetic mice, vitamin D deficiency accelerates type 1 diabetes [8], while chronic administration of 1,25D3 reduces the incidence of diabetes in these mice, principally by modulating immune mechanisms [911]. Streptozotocin-induced diabetes is another animal model of diabetes that comprises both toxic and inflammatory mechanisms [12]. Mononuclear infiltration and the altered morphology of islets coupled with the disappearance of beta cells are among the histological changes reported in the pancreas of these animals [12, 13]. Indeed, hyperglycemia and hypoinsulinemia, have also been described in this experimental model [14, 15]. This hypoinsulinemia was associated with increased insulin binding in the kidney [16, 17] and liver [14, 18] and with somewhat controversial insulin binding results in adipose tissue [15, 19, 20]. The increases in renal and hepatic IRs, were accompanied by elevated IR mRNA expression in both tissues, and it could be reversed by treatment with insulin [2123]. Despite the different alterations in IR mRNA levels and insulin binding, streptozotocin-induced diabetic rats characteristically display insulin resistance [15, 18, 19]. In this diabetic model, the administration of 1,25D3 for 8 weeks was reported to improve diabetes attenuating pancreatic islet damage and decreasing the insulin requirements [12].

1,25D3 acts in its genomic effects as a ligand for the vitamin D receptor (VDR, NR1I1) [24]. This receptor is a member of the superfamily of nuclear receptors that regulates gene expression as a vitamin D-dependent transcription factor. It exerts this action by binding, preferentially as a heterodimer with the retinoid × receptor (RXR), to vitamin D response elements (VDREs) in the promoter regions of target genes [25]. A VDRE generally consists of two direct imperfect repeats of six nucleotides separated by a three-nucleotide spacer. The VDR occupies the 3' half-site, while the RXR binds to the 5' half-site [26]. Sequence variations have been detected in the 3' half-element, the 5' half-element, the spacer, and in the sequences flanking the VDREs and these variations appear to be important in receptor-binding efficiency [27, 28]. The identification of VDREs has only been possible in a very limited number of vitamin D-regulated genes [26, 27]. Our group has reported the first identification of a functional VDRE overlapped by a downstream AP-2-like site that specifically binds VDR in the human IR gene promoter [29]. This VDRE accounted for the transcriptional induction of this gene by 1,25D3 in U-937 human cells [7, 29].

With these antecedents, the aim of the present investigation was to study the effects of the treatment with 1,25D3 [150 IU/Kg (3.75 μg/Kg) one a day, for 15 days] to non-diabetic and streptozotocin-induced diabetic rats. The results indicated that while treatment with 1,25D3 had practically no effect on non-diabetic rats, the same treatment in streptozotocin-induced diabetic rats corrected in part the over-expression of the IR gene in liver and adipose tissue, although it did not revert the hyperglycemia, hypoinsulinemia, glycosuria or ketonemia of these diabetic animals. At the same time, it produced normalization of the IR number without alterations in the receptor affinity and with an improvement of the insulin response in terms of glucose transport in isolated adipocytes of these diabetic animals. In addition, a computer search in the rat insulin receptor gene promoter revealed the existence of two DNA sequences: -256/-219 bp and -653/-620 bp, the first overlapped by three and the second by four AP-2-like sites. These sequences represent two candidate VDREs that could respond to 1,25D3 via activation of VDR, although this remains to be further investigated.

Results

Body weights and plasma values of glucose, insulin, 25-hidroxyvitamin D3, calcium, phosphorus and proteins

Treatment with 1,25D3 to non-diabetic rats increased the rat body weights but did not affect any parameter in the plasma of these animals (Table 1). The injection of streptozotocin reduced body weights and induced hyperglycemia and hypoinsulinemia (Table 1). Treatment with 1,25D3 to streptozotocin-induced diabetic rats increased body weights but did not revert the hyperglycemia and the hypoinsulinemia induced by the diabetes, while calcium and phosphorus plasma levels were increased and plasma levels of 25-hydroxyvitamin D3 were maintained in the normal range [30] (Table 1).

Table 1 Body weights and plasma values of the four groups of rats under study

Values of specific gravity, pH, leukocytes, nitrite, protein, glucose, ketones, urobilinogen, bilirubin and blood in urine

Treatment with 1,25D3 to non-diabetic rats did not affect any of the above cited urine parameters under study. Streptozotocin injection induced glycosuria (> 500 mg/dl) and the appearance of leukocytes, ketones and blood in the urine of these animals. Treatment with 1,25D3to streptozotocin-induced diabetic rats did not modify the glycosuria and the ketonuria of the diabetes, but partially decreased the number of leukocytes in urine (from 10–25 to 5–10 leukocytes/μl).

Total DNA, RNA and protein content of the kidney, liver and epididymal adipose tissue

Treatment with 1,25D3 to non-diabetic rats increased both the protein content and the indicator of cell size (protein/DNA) in the kidney, liver and adipose tissue of these animals while the DNA and RNA content remained unaltered (Table 2). Streptozotocin injection incremented the protein content in these three tissues, accompanied in the case of kidney with an increment of the protein/DNA ratio and in the case of the liver with an increment of the DNA content. Treatment with 1,25D3 of streptozotocin-induced diabetic rats practically did not modify neither the high values of protein induced by the diabetes in the three tissues nor the increases in the DNA content in kidney and liver. In addition, this treatment did not alter the RNA content in any of these three tissues (Table 2).

Table 2 Total DNA, RNA and protein content in different tissues from the four groups of rats

Insulin receptor mRNA levels in the kidney, liver and epididymal adipose tissue

Northern blot assays of kidney (Figure 1, panels A and B), liver (Figure 1, panels C and D) and epididymal adipose tissue (Figure 1, panels E and F) from control animals revealed two major IR mRNA species of approximately 9.5 and 7.5 Kb in size. The relative amounts of the two IR mRNA species measured as 9.5 Kb/7.5 Kb ratio were: 2.2 ± 0.2 in kidney, 1.2 ± 0.2 in liver and 1.1 ± 0.03 in adipose tissue. Treatment with 1,25D3 to non-diabetic rats did not affect any of the two IR mRNA species in any of the three tissues studied (Figure 1, panels A to F). IR mRNA species were expressed per unit of RNA in view of the observation (Table 2) that treatment with 1,25D3 to both non-diabetic and diabetic rats did not alter RNA content per gram of tissue in any of the three tissues. The induction of diabetes by streptozotocin produced important increments in the levels of both IR mRNA species in liver and adipose tissue but not in kidney. In particular, 78% increase corresponding to the 9.5 Kb and 69% to the 7.5 Kb species in liver (Figure 1, panel D), and 58% increase corresponding to the 9.5 Kb and 92% increase to the 7.5 Kb species in the adipose tissue (Figure 1, panel F). Treatment with 1,25D3 to diabetic rats partially prevented the over-expression of the IR mRNA in the liver (Figure 1, panels C and D) and adipose tissue (Figure 1, panels E and F) of these animals. The 9.5 Kb species was particularly decreased in liver and the 7.5 kb species in the adipose tissue. However, treatment with 1,25D3 to diabetic rats did not affect IR expression in kidney (Figure 1, panels A and B).

Figure 1
figure 1

Renal, hepatic and adipose IR mRNA levels. Insulin receptor (IR) mRNA levels in the kidney (A, B), liver (C, D) and epididymal adipose tissue (E, F) of sham-treated rats (Control), rats treated with 1,25D3 [150 IU/Kg (3.75 μg/Kg) one a day, for 15 days](1,25D3), streptozotocin-induced diabetic rats (STZ) and streptozotocin-induced diabetic rats treated with 1,25D3 (STZ+1,25D3). On the right, the autoradiograph of a representative Northern blot experiment using the kidney (B), liver (D) and adipose tissue (F) of two rats. The sizes of the two major IR mRNA species are shown in the margin. On the left, densitometric analysis of three independent Northern blot experiments on kidney (A), liver (C), and adipose tissue (E). IR mRNA species were quantified separately, normalized to the respective 28S rRNA values and expressed per unit of RNA as a percentage of the values obtained in Control-rats (mean ± SEM). a p < 0.05 vs. Control-rats; b p < 0.05 vs. 1,25D3-rats; and c p < 0.05 vs. STZ-rats.

Insulin binding in isolated adipocytes

Insulin binding studies showed that treatment with 1,25D3 to non-diabetic rats did not alter IR number or IR affinity as reflected by the ED:50 value in adipocytes of these animals (Figure 2 and Table 3). Streptozotocin injection produced a 64% decrease in the number of IRs without altering IR affinity. Treatment with 1,25D3 to streptozotocin-induced diabetic rats inverted the decrease in the IR number induced by the diabetes to values not significantly different from those observed in control adipocytes and without altering IR affinity (Figure 2 and Table 3).

Table 3 Insulin binding parameters in epididymal adipocytes from the four groups of rats
Figure 2
figure 2

Insulin receptor number per cell. Insulin receptor number in epididymal adipocytes from sham-treated rats (Control), rats treated with 1,25D3 [150 IU/Kg (3.75 μg/Kg) one a day, for 15 days] (1,25D3), streptozotocin-induced diabetic rats (STZ) and streptozotocin-induced diabetic rats treated with 1,25D3 (STZ+1,25D3). Values are the mean ± SEM of 3–9 determinations in each group. a p < 0.05 vs. Control-rats; b p < 0.05 vs. 1,25D3-rats; and c p < 0.05 vs. STZ-rats.

Given that both the weight of the adipose tissue and the diameter of the adipocytes were markedly decreased by the diabetes (Table 3), we also calculated the IR number per cell surface area (μm2). In these conditions, although the effect of the diabetes decreasing IR number was somewhat minor (44%), the tendency of 1,25D3 to normalize the IR number was also evident (Table 3).

Glucose transport in isolated adipocytes

Treatment with 1,25D3 to non-diabetic rats did not alter basal or insulin-stimulated glucose transport in adipocytes from these animals (Figure 3). The injection of streptozotocin clearly decreased both basal and insulin-stimulated glucose transport. The treatment with 1,25D3 to streptozotocin-induced diabetic rats, improved by 107% the decreased basal glucose transport and by 71% the insulin-stimulated glucose transport in adipocytes from these diabetic animals (Figure 3).

Figure 3
figure 3

Basal and insulin stimulated glucose transport. Basal and insulin stimulated glucose transport in epididymal adipocytes from sham-treated rats (Control), rats treated with 1,25D3 [150 IU/Kg (3.75 μg/Kg) one a day, for 15 days] (1,25D3), streptozotocin-induced diabetic rats (STZ) and streptozotocin-induced diabetic rats treated with 1,25D3 (STZ+1,25D3). Values are the mean ± SEM of 6–7 determinations within each group. a p < 0.05 vs. Control-rats; b p < 0.05 vs. 1,25D3-rats; and c p < 0.05 vs. STZ-rats.

Computer analysis of DNA sequences in the rat IR gene promoter

We performed a computer search of virtual VDRE sequences in the rat IR gene promoter by homology with two consensus VDRE sequences. The first consensus (5'GGGTCA NNGGGGGCA 3') was previously compiled and reported by us from a series of described functional VDRE sequences in various 1,25D3-stimulated promoters [29]. The search with this first consensus indicated no sequence identical to this VDRE sequence, and neither was any sequence found to show 1, 2, or 3 base variations. Nevertheless, we detected three virtual VDRE sequences showing a difference of four bases from this consensus: -249/-235, -637/-623 and -916/-902 (Figure 4, panel A). Comparison of the 3' half-element of each of these sequences with the 3' half-element of this first VDRE consensus, indicated complete homology only in the case of the -916/-902 sequence, and two base variations in the other two (Figure 4, panel A).

Figure 4
figure 4

Panel A: Computer search of vitamin D response elements (VDREs) in the rat insulin receptor (IR) gene promoter. Potential VDRE sequences in the rat IR gene promoter were identified by their homology with two consensus VDRE sequences using the SEQFIND programme developed in our laboratory. Nucleotides identical to the first VDRE consensus are labelled by filled circles () and nucleotides identical to the second VDRE consensus are labelled by empty circles (). Panel B: Computer search of AP-2-like sites in the rat insulin receptor (IR) gene promoter. Potential AP-2-like sites in the rat IR gene promoter were identified by their homology with the AP-2 consensus (5'GCCN3GGG3') [32] using the SEQFIND programme. We detected 17 AP-2-like sites showing a difference of only 1 base from this consensus. Of these, only 7 AP-2-like sites flanking or overlapping the potential VDRE sequences in panel A are shown. Nucleotides identical to the AP-2 consensus are labelled by stars (). Panel C: Candidate VDRE sequences overlapped by AP-2-like sites in the rat insulin receptor gene promoter. Nucleotides identical to the first VDRE consensus are labelled by filled circles () and nucleotides identical to the second VDRE consensus are labelled by empty circles ()

The second computer search of virtual VDREs, now by homology with the second consensus VDRE sequence (5'PuGGTCA NNPuPuGTTCA 3'), was also performed. Our consensus and the latter proposed by Colnot et al. [27], and Haussler et al. [26] have identical 5' half-elements, but ours contains two GG instead of two TT, in the 3' half-element. The spacers are identical, with a G or Pu at position 3, which appears to be important in VDR binding. The results of this computer search by homology with this second consensus indicated no sequence identical to this VDRE, and neither was any sequence found to show 1, 2, or 3 base variations. Nevertheless, we detected three virtual VDRE sequences showing a difference of four bases from this second consensus: -25/-11, -247/-233 and -881/-867 (Figure 4, panel A). Of these, the -247/-233 sequence was in part the same as that of the -249/-235 sequence detected with the first consensus (Figure 4, panel A). However, the 3' half-element of this second sequence has only one base variation instead of two that the first sequence had (Figure 4, panel A). Therefore, in principle there were six virtual VDRE sequences (Figure 4, panel A).

The next step was to check for the possibility of cis-elements related to these virtual VDRE sequences in the rat IR gene promoter. Thus, a computer search of AP-1 and AP-2 sites by homology with the AP-1 site consensus (5'TGAC/GTCA3') [31], and the AP-2 site consensus (5'GCCN3GGG3') [32], respectively, was performed.

Exploring the AP-1 sites, we detected two AP-1-like sites showing a difference of only one base from this consensus: the -374/-368 (5'TGTCTCA3') and the -963/-957 (5'TGAGCCA3'). However, these two sites are neither adjacent nor overlapped in any of our virtual VDRE sequences. Nevertheless, as far as we know, this is the first mention of these AP-1-like sites in the rat IR gene promoter. With respect to the AP-2 sites, we detected seventeen AP-2-like sites showing a difference of only one base from this consensus: -71/-63, -79/-71, -171/-163, -195/-187, -230/-222, -245/-237, -253/-245, -405/-397, -465/-457, -480/-472, -545/-537, -634/-626, -635/-627, -642/-634, -650/-642, -1087/-1079, -1088/-1080. Of these AP-2-like sites, only seven are included in panel B of the Figure 4. These were those flanking or overlapping some of the potential VDRE sequences indicated in panel A of the Figure 4. Of them, the -230/-222 (1) was located downstream, and the -245/-237 (2) and the -253/-245 (3) were overlapping the virtual VDRE sequences: -249/-235 and -247/-233, as shown in Figure 4, panel C. Moreover, the AP-2-like sites: -634/-626 (4) and -635/-627 (5) (Figure 4, panel B) were overlapping the virtual VDRE sequence: -637/-623 (Figure 4, panel C), and the AP-2-like sites: -642/-634 (6) and -650/-642 (7) (Figure 4, panel B) were located in tandem upstream of this last virtual VDRE sequence: -637/-623 (Figure 4, panel C).

Therefore, we postulate as candidate VDRE sequences overlapped by AP-2-like sites those represented in Figure 4, panel C. One, located between -256 and -219 bp with three AP-2-like sites, and the other extending from -653 and -620 bp of the rat IR gene promoter with four AP-2-like sites. Separately or together, these VDRE sequences could form a locus that may respond to 1,25D3 via activation of VDR.

Discussion

Effects of treatment with 1,25D3 to non-diabetic rats

Treatment with 1,25D3 to non-diabetic rats increased body weights but it did not affect any parameter measured in the plasma or urine of these animals, including the plasma levels of calcium, phosphorus and 25-hydroxyvitamin D3. Indeed, we did not detect significant differences in calcium and phosphorus, despite the increases of approximately 30%. In relation to 25-hydroxyvitamin D3 a number of studies [33] has shown that human serum values of 25-hydroxyvitamin D3 were maintained within a normal range across vitamin D supplies from 80 to 1000 IU daily, indicating the existence of an important homeostatic control system regulating human 25(OH)D3 serum concentration. Considering this, in the present study we hypothesize the existence of a similar counter regulatory mechanism in the rat. This mechanism could maintain rat plasma values of 25-hydroxyvitamin D3 within the observed normal range of 15–19 ng/ml, with vitamin D supplies of 150 IU daily administered by injection plus the vitamin D consumed in the diet.

Treatment with 1,25D3 increased the protein content and the indicator of cell size (protein/DNA) in the kidney, liver, and adipose tissue of non-diabetic rats, although the DNA and RNA content remained unaltered. These increments might indicate hypertrophy. However, although treatment with 1,25D3 increased rat body weights, it did not particularly augment the weight of these tissues. Another explanation is that 1,25D3 might cause hyperplasia. Regarding this, conflicting results have been reported [34, 35].

We also observed that treatment with 1,25D3 to non-diabetic rats did not affect the expression of any of the two 9.5 and 7.5 Kb IR mRNA species present in kidney, liver, and adipose tissue. Our group and others have previously reported the presence of these species and their relative proportions in rat tissues [21, 36]. Since there is only one rat IR gene, variation in transcript length may reflect alternative RNA processing events, such as polyadenylation at different 3' end sites in the final untranslated exon [37].

In addition, we detected that treatment with 1,25D3 to non-diabetic rats affected neither the IR number nor the insulin response in terms of glucose transport in isolated adipocytes from these animals.

The lack of effects of 1,25D3 on IR expression and insulin activity in tissues of non-diabetic rats could be related to the absence of VDR in these tissues. However, both the VDR protein and VDR mRNA have been detected in the rat kidney [38, 39] and liver [40, 41], and while a VDR-like protein was identified in 3T3-L1 adipose cells [42] VDR mRNA has been found in rat perirenal adipose tissue and 3T3-L1 adipose cells [43]. Another possibility could be the absence of regulation of VDR by its ligand in these tissues. Such lack of regulation has been described in the kidney [39]. In the liver there are no previous data, while in adipose cells treatment with 1,25D3 led to an increase in VDR mRNA levels [43] and a decrease in VDR protein [35]. Regarding this and given our previous experience in determining both VDR protein and VDR mRNA expression and their regulation by 1,25D3 in human cells [6, 7], we determined the expression of this receptor and its regulation by 1,25D3 at the protein level in the kidney, liver and adipose tissue of non-diabetic rats. Our results showed a VDR protein only in the kidney (data not shown). More experiments are necessary to confirm this.

Effects of treatment with 1,25D3 to streptozotocin-induced diabetic rats

Treatment with 1,25D3 to streptozotocin-induced diabetic rats increased their body weight, but it did not revert the hyperglycemia and the hypoinsulinemia provoked by the streptozotocin. We also detected normal levels of 25-hydroxyvitamin D3 in streptozotocin-induced diabetic rats in accordance with other authors [30]. These unchanged levels of 25-hydroxyvitamin D3 in streptozotocin-induced diabetic rats neither were altered by the treatment with 1,25D3, possibly due to the postulated plasma homeostatic control system that could regulate 25-hydroxyvitamin D3 plasma levels also in the absence of insulin in both groups of diabetic rats. In addition, we observed increased plasma levels of calcium and phosphorus after the treatment with 1,25D3. With regard to these last parameters, it is known that calcium per se is important for insulin secretion, as well as for correction of glucose intolerance [44, 45]. Moreover, calcium and phosphorus have been described as regulators of VDR in renal and intestinal tissues [39, 46].

Treatment with 1,25D3 to streptozotocin-diabetic rats did not modify the glycosuria and the ketonuria induced by the diabetes, but partially decreased the number of leukocytes in urine. This latter effect could be due to the broadly reported anti-inflammatory effects of 1,25D3 and/or to a possible anti-proliferative effect of this hormone, as observed by our group and others in cell culture [6, 33].

With respect to tissue effects of 1,25D3 in streptozotocin-induced diabetic rats, we observed the permanence of high protein and DNA content values induced by the diabetes, accompanied by slight increases of the protein/DNA ratio values in the kidney and liver, but not in the adipose tissue. This suggests a certain hypertrophy produced by the action of 1,25D3 in these two tissues. However, although treatment with 1,25D3 to streptozotocin-diabetic rats increased rat body weights, it did not increase the weight of the kidney, liver or adipose tissue. The other explanation, hyperplasia, has conflicting results [33, 34].

The induction of diabetes by streptozotocin did not affect IR mRNA species in the kidney. This finding was not in accord with previous results of increased renal IR mRNA levels, quantified by slot blot, in a similar rat model [22]. This may reflect the failure to alter VDR number in the kidney of streptozotocin-diabetic rats [34]. The treatment with 1,25D3 to streptozotocin-induced diabetic rats neither affected IR mRNA species in the kidney.

Contrarily, streptozotocin produced an important increment in the levels of both IR mRNA species in liver. These findings agree with previous results of increased IR mRNA levels quantified by slot blot [22], Northern blot analysis [21, 47] and increased IR gene transcription [21] in the liver of streptozotocin-induced diabetic rats. Streptozotocin also increased the levels of both IR mRNA species in the adipose tissue. This increment represents the first demonstration of an in vivo regulation of IR mRNA levels in the adipose tissue of streptozotocin-induced diabetic rats.

In addition, we observed that the treatment with 1,25D3 to streptozotocin-induced diabetic rats partially prevented the over-expression of IR mRNA induced by the diabetes in the liver and adipose tissue of these animals. The 9.5 Kb IR mRNA species was particularly decreased by 1,25D3 in the liver, and the 7.5 Kb species in the adipose tissue. Therefore, these results provide the first demonstration of an in vivo tissue-specific regulation of rat IR mRNA levels by 1,25D3 in streptozotocin-induced diabetic rats.

A role for VDR in such genomic actions of 1,25D3 in liver and adipose tissue is possible, although to our knowledge the regulation of VDR by 1,25D3 has not yet been studied in these tissues of streptozotocin-diabetic rats.

The correction by 1,25D3 of the over-expression of IR expression at the RNA level in the liver and adipose tissue of streptozotocin-diabetic animals could represent a beneficial fact of amelioration of the diabetes. In this sense, other facts of amelioration of diabetes were the almost normalization of the number of IRs in adipocytes of streptozotocin-diabetic animals treated with 1,25D3, without alterations in receptor affinity, and the improvement of both basal and insulin-stimulated glucose transport in adipocytes of these diabetic animals treated with 1,25D3.

The present in vivo results do not allow to determine whether the prevention by 1,25D3 of the over-expression of IR mRNA induced by the diabetes in liver and adipose tissue of streptozotocin-diabetic animals is due to transcriptional and/or post-transcriptional regulation. Nevertheless, our previous in vitro results demonstrated that 1,25D3 increased IR mRNA levels via mechanisms involving direct transcriptional activation of the human IR gene [57]. Moreover, we identified an active VDRE in the human IR gene promoter that accounted for the transcriptional induction of this gene by 1,25D3 in U-937 cells [29].

Therefore, in support of the possible participation of a transcriptional regulation by 1,25D3 in the present in vivo processes, we have detected by computer search in the rat IR gene promoter two candidate VDREs overlapped by various AP-2-like sites. One VDRE is located between -256 and -219 bp, while the other extends from -653 and -620 bp of the rat IR gene promoter (Figure 4, panel C). Individually or in conjunction, these potential VDREs could form a locus that may respond to 1,25D3 via activation of VDR. This locus may mediate the cross-talk between vitamin D and insulin signalling, although this remains to be determined.

Conclusion

The present results indicate that while treatment with 1,25D3 had practically no effect on non-diabetic rats, the same treatment in streptozotocin-induced diabetic rats corrected in part the over-expression of the IR gene in liver and adipose tissue, although it did not revert the hyperglycemia, hypoinsulinemia, glycosuria or ketonemia of these diabetic animals. At the same time, it produced normalization of the IR number without alterations in the receptor affinity, but with an improvement of the insulin response in terms of glucose transport in adipocytes of these diabetic animals. These genomic actions of 1,25D3 could represent beneficial facts of amelioration of diabetes via mechanisms, possibly involving direct transcriptional activation of the rat IR gene. A computer search in the rat IR promoter revealed the existence of two candidate vitamin D response element (VDRE) sequences, located at -256/-219 bp and -653/-620 bp and overlapped by three and four AP-2-like sites respectively. These candidate VDREs may respond to 1,25D3 via activation of the vitamin D receptor, although this remains to be investigated.

Methods

Animals and sampling

Four groups of male Wistar rats weighing 200–220 g at the onset were used in this study. They were kept under standard conditions of light and temperature and given free access to tap water and standard laboratory chow (Panlab, A04) containing 8.8 mg/g of calcium, 5.9 mg/g of phosphorus and 1.5 IU/g of vitamin D. The first group comprised of non-diabetic rats receiving sham-treatments during the 15 days of the experimental period (Control-rats). The second group included non-diabetic rats i.p. injected with 1,25D3 (Calcijex, Abbot) [150 IU/Kg (3.75 μg/Kg) one a day, for 15 days] (1,25D3-rats). The third group consisted of rats rendered diabetic by a single injection of streptozotocin (Sigma) [65 mg/kg] on the 8th day of the 15-day period (STZ-rats). Finally, the fourth group comprised of rats rendered diabetic by streptozotocin on the 8th day, and i.p. injected with 1,25D3 [150 IU/Kg (3.75 μg/Kg) one a day, for 15 days] (STZ+1,25D3-rats). The National Guide for the Care and Use of Laboratory Animals was strictly followed during this study.

Urine samples were collected before decapitation without anaesthesia and trunk blood samples were recovered for plasma measurements. On sacrifice, the kidney, liver and epididymal adipose tissue were immediately removed, individually packed and frozen in liquid nitrogen for nucleic acid and protein determinations. Isolated adipocytes were obtained from fresh epididymal adipose tissue for insulin binding and insulin activity assays.

Analytical methods

Different parameters in urine were measured using the Combur Test (Roche). Plasma insulin levels were determined by radioimmunoassay using rat insulin as standard. Plasma concentrations of glucose, 25-hydroxyvitamin D3, calcium, phosphorus and proteins were estimated using commercially available techniques. The total DNA and RNA contents of the homogenates from individual tissues were estimated by spectrofluorometry. The protein content was determined by the Bradford method.

RNA blot assays

For Northern blot assays, RNA samples (40 μg) of kidney, liver and epididymal adipose tissue, were denatured, electrophoresed in 1.1% agarose-formaldehyde gels and blotted onto nylon membranes [36]. Ethidium bromide staining of the 28S and 18S ribosomal RNAs was routinely checked before blotting as a control of sample loading and after blotting as a control of RNA transfer. RNA blots were prehybridized, hybridized with excess [32P]-labelled probe (the 0.98-Kb rat IR specific EcoR1 fragment of the p16 clone, gift from Prof. Goldstein), washed under stringent conditions and finally autoradiographed, as described previously [36]. The autoradiographs were scanned with a laser densitometer and the readings normalized with the respective amounts of 28S rRNA, as revealed by ethidium bromide.

Insulin binding assays

Binding assays were carried out as described previously [48]. In brief, isolated adipocytes (0.25 × 106 cells/ml) were incubated at 30°C for 30 min in 400 μl Krebs-Hepes buffer pH 7.6 with mono-[125I]-insulin (0.2 × 10-9 M) (NEN Life Sciences), in the presence or absence of increasing concentrations (0.2 × 10-10 to 0.5 × 10-7 M) of unlabelled insulin (Sigma). Non-specific binding was determined in the presence of 1 × 10-7 M unlabelled insulin and subtracted from each point. The adipocyte number was determined with a Neubauer-type hemocytometer and cell viability (assessed by 0.2% trypan blue exclusion) was greater than 90%. Measurements of the cell diameter were also performed. The data were analysed by the Scatchard method using the LIGAND program of Munson and Rodbard [49].

Glucose transport determinations

Glucose transport was measured in isolated adipocytes by determining the uptake of D-[14C(U)]-glucose (NEN Life Sciences) at trace concentrations using our modification of the method of Kashiwagi et al. [50]. This method is based on the premise that glucose uptake provides a measure of glucose transport when studies are carried out at very low glucose concentrations. In brief, 0.25 × 106 adipocytes/ml were incubated with 5 × 10-7 M D-[14C(U)]-glucose for 1 h at 37°C in the presence and absence of a concentration of insulin (10-8 M) that gave maximal glucose transport in earlier experiments [6]. Assays were stopped by washing the cells with ice-cold PBS. The cells were solubilized in 0.5% SDS and 0.1 N NaOH and the associated radioactivity was counted.

Computer analysis of DNA sequences in the rat IR gene promoter

The identification of virtual VDRE sequences in the Rattus norvegicus partial IR promoter [GenBank:AJ006071] was carried out by homology with two consensus VDRE sequences using the SEQFIND programme developed in our laboratory and previously employed in the computer analysis of various DNA sequences including VDREs [29]. The first consensus (5'GGGTCA NNGGGGGCA 3') had been compiled and reported by us from a series of described functional VDRE sequences in various 1,25D3-stimulated promoters [29]. The second consensus VDRE sequence (5'PuGGTCA NNPuPuGTTCA 3') had been proposed by Colnot et al. [27], and Haussler et al. [26]. The presence of AP-1 and AP-2 sites flanking or overlapping the potential VDREs were identified by their homology with the AP-1 site consensus (5'TGAC/GTCA3') [31] and the AP-2 site consensus (5'GCCN3GGG3') [32] respectively, also using the SEQFIND programme.

Statistical analysis

Unless otherwise indicated, the data are expressed as the mean ± SEM. A comparison between the groups was carried out using the two-tailed, unpaired Student t-test and/or ANOVA comparison, as appropriate. Differences were considered statistically significant when p < 0.05.

Abbreviations

1,25D3:

1,25-dihydroxyvitamin D3

IR:

insulin receptor

VDR:

vitamin D receptor

RXR:

retinoid × receptor

VDRE:

vitamin D response element

STZ:

streptozotocin.

References

  1. Kadowaki S, Noman AW: Dietary vitamin D is essential for normal insulin secretion from the perfused rat pancreas. J Clin Invest 1984, 73: 759-766. 10.1172/JCI111269.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  2. Lips P: Vitamin D physiology. Prog Biophys Mol Bio 2006, 92: 4-8. 10.1016/j.pbiomolbio.2006.02.016.

    Article  CAS  Google Scholar 

  3. Norman AW, Frankel JB, Heldt AM, Grodsky GM: Vitamin D deficiency inhibits pancreatic secretion of insulin. Science 1980, 209: 823-825. 10.1126/science.6250216.

    Article  CAS  PubMed  Google Scholar 

  4. Bourlon PM, Billaudel B, Faure-Dussert A: Influence of vitamin D 3 deficiency and 1, 25 dihydroxyvitamin D 3 on the novo insulin biosynthesis in the islets of the rat endocrine pancreas. J Endocrinol 1999, 160: 87-95. 10.1677/joe.0.1600087 10.1677/joe.0.1600087

    Article  CAS  PubMed  Google Scholar 

  5. Leal MA, Aller P, Más A, Calle C: The effect of 1, 25-dihydroxyvitamin D 3 on insulin binding, insulin receptor mRNA levels, and isotype RNA pattern in U-937 human promonocytic cells. Exp Cell Res 1995, 217: 189-194. 10.1006/excr.1995.1078 10.1006/excr.1995.1078

    Article  CAS  PubMed  Google Scholar 

  6. Maestro B, Campión J, Dávila N, Calle C: Stimulation by 1, 25 dihydroxyvitamin D 3 on insulin receptor expression and insulin responsiveness for glucose transport in U-937 human promonocytic cells. Endocr J 2000, 47: 383-391. 10.1507/endocrj.47.383 10.1507/endocrj.47.383

    Article  CAS  PubMed  Google Scholar 

  7. Maestro B, Molero S, Bajo S, Dávila N, Calle C: Transcriptional activation of the human insulin receptor gene by 1, 25-dihydroxyvitamin D 3 . Cell Biochem Funct 2002, 20: 227-232. 10.1002/cbf.951 10.1002/cbf.951

    Article  CAS  PubMed  Google Scholar 

  8. Giulietti A, Gysemans C, Stoffels K, Van Etten E, Decallonne B, Overbergh L: Vitamin D deficiency in early life accelerates Type 1 diabetes in non-obese diabetic mice. Diabetologia 2004, 47: 451-462. 10.1007/s00125-004-1329-3.

    Article  CAS  PubMed  Google Scholar 

  9. Mathieu C, Laureys J, Sobis H, Vandeputte M, Waer M, Bouillon R: 1, 25-dihydroxyvitamin D 3 prevents insulitis in NOD mice. Diabetes 1992, 41: 1491-1495. 10.2337/diabetes.41.11.1491 10.2337/diabetes.41.11.1491

    Article  CAS  PubMed  Google Scholar 

  10. Mathieu C, Waer M, Laureys J, Rutgeerts O, Bouillon R: Prevention of autoimmune diabetes in NOD mice by 1, 25-dihydroxyvitamin D 3 . Diabetologia 1994, 37: 552-558. 10.1007/BF00403372 10.1007/BF00403372

    Article  CAS  PubMed  Google Scholar 

  11. Mathieu C, Gysemans C, Giulietti A, Bouillon R: Vitamin D and diabetes. Diabetologia 2005, 48: 1247-1257. 10.1007/s00125-005-1802-7.

    Article  CAS  PubMed  Google Scholar 

  12. del Pino-Montes J, Benito GE, Fernández-Salazar MP, Coveñas R, Calvo JJ, Bouillon R, Quesada JM: Calcitriol improves streptozotocin-induced diabetes and recovers bone mineral density in diabetic rats. Diabetes 1997, 46: 1346-1353. 10.2337/diabetes.46.8.1346.

    Article  Google Scholar 

  13. Steiner H, Oelz O, Zahnd G, Froesch ER: Studies on islet cell regeneration, hyperplasia and intrainsular cellular interrelation in long lasting streptozotocin diabetes in rats. Diabetologia 1970, 6: 558-564. 10.1007/BF00418221.

    Article  CAS  PubMed  Google Scholar 

  14. Davidson MB, Kaplan SA: Increased insulin binding by hepatic plasma membranes from diabetic rats: normalization by insulin therapy. J Clin Invest 1977, 59: 22-30. 10.1172/JCI108618.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  15. Kasuga MY, Iwamoto Y, KosaKa K: Insulin binding and glucose metabolism in adipocytes of streptozotocin-diabetic rats. Am J Physiol 1978, 235: E175-E182.

    CAS  PubMed  Google Scholar 

  16. Papachistodoulou DK, Bass PS, Davey PG, Thomas JH: Insulin binding and degradation by kidney cell membranes of streptozotocin-diabetic rats. Horm Metab Res 1982, 14: 345-350.

    Article  Google Scholar 

  17. Rabkin R, Hirayama P, Roth RA, Frank BH: Effect of experimental diabetes on insulin binding by renal basolateral membranes. Kidney Int 1986, 30: 348-354. 10.1038/ki.1986.191.

    Article  CAS  PubMed  Google Scholar 

  18. Kadowaki T, Kasuga M, Akanuma Y, Ezaki O, Takaku F: Decreased autophosphorylation of the insulin receptor-kinase in streptozotocin-diabetic rats. J Biol Chem 1984, 259(22):14208-14216.

    CAS  PubMed  Google Scholar 

  19. Kobayashi M, Olefsky JM: Effects of streptozotocin-induced diabetes on insulin binding, glucose transport, and intracellular glucose metabolism in isolated rat adipocytes. Diabetes 1979, 28: 87-95.

    Article  CAS  PubMed  Google Scholar 

  20. Whittaker J, Cuthbert C, Hammond V, Alberti KG: Impaired insulin binding in isolated adipocytes in experimental diabetic ketoacidosis. Diabetologia 1981, 21(6):563-568.

    CAS  PubMed  Google Scholar 

  21. Tozzo E, Desbuquois B: Effect of STZ-induced diabetes and fasting on insulin receptor mRNA expression and insulin receptor gene transcription in rat liver. Diabetes 1992, 41: 1609-1616. 10.2337/diabetes.41.12.1609.

    Article  CAS  PubMed  Google Scholar 

  22. Sechi IA, Griffin CA, Grady E, Grunfeld C, Kalinyak JE, Schambelam M: Tissue-specific regulation of insulin receptor mRNA in rats with STZ-induced diabetes mellitus. Diabetes 1992, 41: 1113-1118. 10.2337/diabetes.41.9.1113.

    Article  CAS  PubMed  Google Scholar 

  23. Vidal H, Auboeuf D, Beylot M, Riou JP: Regulation of insulin receptor mRNA splicing in rat tissues. Effect of fasting, aging and diabetes. Diabetes 1995, 44: 1196-1201. 10.2337/diabetes.44.10.1196.

    Article  CAS  PubMed  Google Scholar 

  24. Nuclear Receptors Nomenclature Committee: A unified nomenclature system for the nuclear receptor superfamily. Cell 1999, 97: 161-163. 10.1016/S0092-8674(00)80726-6.

    Article  Google Scholar 

  25. Norman AW: Vitamin D receptor: New assignments for an already busy receptor. Endocrinology 2006, 147: 5542-5548. 10.1210/en.2006-0946.

    Article  CAS  PubMed  Google Scholar 

  26. Haussler MR, Whitfield GK, Haussler CA, Hsieh J-C, Thompson PD, Selznick SH, Encinas Dominguez C, Jurutka PW: The nuclear vitamin D receptor: biological and molecular regulatory properties revealed. J. Bone Miner Res 1998, 13: 325-349. 10.1359/jbmr.1998.13.3.325.

    Article  CAS  PubMed  Google Scholar 

  27. Colnot S, Lambert M, Blin C, Tomasset C, Perret C: Identification of DNA sequences that bind retinoid × receptor-1, 25(OH) 2 D 3 -receptor heterodimers with high affinity. Mol Cell Endocrinol 1995, 113: 89-98. 10.1016/0303-7207(95)03618-H 10.1016/0303-7207(95)03618-H

    Article  CAS  PubMed  Google Scholar 

  28. Koszewski NJ, Malluche HH, Russell J: Vitamin D receptor interactions with positive and negative DNA response elements: an interference footprint comparison. J Steroid Biochem Mol Biol 2000, 72: 125-132. 10.1016/S0960-0760(00)00022-4.

    Article  CAS  PubMed  Google Scholar 

  29. Maestro B, Dávila N, Carranza MC, Calle C: Identification of a vitamin D response element in the human insulin receptor gene promoter. J Steroid Biochem Mol Biol 2003, 84: 223-230. 10.1016/S0960-0760(03)00032-3.

    Article  CAS  PubMed  Google Scholar 

  30. Nyomba BL, Bouillon R, Lissens W, Van Baelen H, De Moor P: 1, 25-Dihydroxyvitamin D and vitamin D-binding protein are both decreased in streptozotocin-diabetic rats. Endocrinology 1985, 116(6):2483-2488.

    Article  CAS  PubMed  Google Scholar 

  31. Lee W, Mitchell R, Tjian R: Purified transcription factor AP-1 interacts with TPA-inducible enhancer elements. Cell 1987, 49: 741-752. 10.1016/0092-8674(87)90612-X.

    Article  CAS  PubMed  Google Scholar 

  32. Hilger-Eversheim K, Moser M, Schorle H, Buettner R: Regulatory roles of AP-2 transcription factors in vertebrate development, apoptosis and cell-cycle control. Gene 2000, 260: 1-12. 10.1016/S0378-1119(00)00454-6.

    Article  CAS  PubMed  Google Scholar 

  33. Chatterjee M: Vitamin D and genomic stability. Mutat Res Fund Mol M 2001, 475: 69-88. 10.1016/S0027-5107(01)00080-X.

    Article  CAS  Google Scholar 

  34. Stone LA, Weaver VM, Bruns ME, Welsh J: Vitamin D receptors in intestine, kidney and thymus of streptozotocin diabetic rats. Diabetes Res 1990, 15(4):165-172.

    CAS  PubMed  Google Scholar 

  35. Vu D, Ong JM, Clemens TL, Kern PA: 1, 25-Dihydroxyvitamin D induces lipoprotein lipase expression in 3T3-L1 cells in association with adipocyte differentiation. Endocrinology 1996, 137: 1540-1544. 10.1210/en.137.5.1540.

    CAS  PubMed  Google Scholar 

  36. Campión J, Aller P, Dávila N, Carranza MC, de Miguel R, Calle C: Tissue-specific modulation of insulin receptor mRNA levels in adrenaline-treated rats. Mol Cell Biochem 1997, 169: 165-169. 10.1023/A:1006872309179.

    Article  PubMed  Google Scholar 

  37. Tewari M, Tewari DS, Taub R: Posttranscriptional mechanism account for differences in steady state levels of insulin receptor messenger RNA in different cells. Mol Endocrinol 1991, 5: 653-660.

    Article  CAS  PubMed  Google Scholar 

  38. Costa EM, Feldman D: Homologous up-regulation of the 1, 25(OH) 2 vitamin D 3 receptor in rats. Biochem Bioph Res Co 1986, 2: 742-747. 10.1016/0006-291X(86)91141-1 10.1016/0006-291X(86)91141-1

    Article  Google Scholar 

  39. Zineb R, Zhor B, Odile W, Marthe R-R: Distinct, tissue specific regulation of vitamin D receptor in the intestine, kidney, and skin by dietary calcium and vitamin D. Endocrinology 1998, 139: 1844-1852. 10.1210/en.139.4.1844.

    Article  CAS  PubMed  Google Scholar 

  40. Sandgren ME, Bronnegard M, DeLuca HF: Tissue distribution of the 1, 25-dihydroxyvitamin D 3 receptor in the male rat. Biochem Bioph Res Co 1991, 181: 611-616. 10.1016/0006-291X(91)91234-4 10.1016/0006-291X(91)91234-4

    Article  CAS  Google Scholar 

  41. Mee AP, Davenport LK, Hoyland JA, Davies M, Mawer EB: Novel and sensitive detection systems for the vitamin D receptor-in situ-reverse transcriptase-polymerase chain reaction and immunogold cytochemistry. J Mol Endocrinol 1996, 16(2):183-195.

    Article  CAS  PubMed  Google Scholar 

  42. Sato M, Hiragun A: Demonstration of 1 alpha, 25-dihydroxyvitamin D3 receptor-like molecule in ST13 and 3T3L1 preadipocytes and its inhibitory effects on preadipocyte differentiation. J Cell Physiol 1988, 135: 545-550. 10.1002/jcp.1041350326.

    Article  CAS  PubMed  Google Scholar 

  43. Kamei Y, Kawada T, Kazuki R, Ono T, Kato S, Sugimoto E: Vitamin D receptor gene expression is up-regulated by 1, 25-dihydroxyvitamin D3 in 3T3-L1 preadipocytes. Biochem Bioph Res Co 1993, 193: 948-955. 10.1006/bbrc.1993.1717.

    Article  CAS  Google Scholar 

  44. Hough S, Russell JE, Teitelbaum SL, Avioli LV: Calcium homeostasis in chronic streptozotocin-induced diabetes mellitus in the rat. Am J Physiol 1982, 242(6):E451-6.

    CAS  PubMed  Google Scholar 

  45. Beaulieu C, Kestekian R, Havrankova J, Gascon-Barre M: Calcium is essential in normalizing intolerance to glucose that accompanies vitamin D depletion in vivo . Diabetes 1993, 42: 35-43. 10.2337/diabetes.42.1.35.

    Article  CAS  PubMed  Google Scholar 

  46. Sriussadaporn S, Wong M-S, Pike JW, Favus M: Tissue specificity and mechanism of vitamin D receptor up-regulation during dietary phosphorus restriction in the rat. J Bone Miner Res 1995, 10(2):271-280.

    Article  CAS  PubMed  Google Scholar 

  47. Amessou M, Bortoli S, Liemans V, Collinet M, Desbuquois B, Brichard S, Girard J: Treatment of streptozotocin-induced diabetic rats with vanadate and phlorizin prevents the over-expression of the liver insulin receptor gene. Eur J Endocrinol 1999, 140: 79-86. 10.1530/eje.0.1400079.

    Article  CAS  PubMed  Google Scholar 

  48. Cabrera R, Mayor P, Fernandez-Ruiz J, Calle C: Insulin binding and action on adipocytes from female rats with experimentally induced chronic hyperprolactinemia. Mol Cell Endocrinol 1988, 58: 167-173. 10.1016/0303-7207(88)90151-7.

    Article  CAS  PubMed  Google Scholar 

  49. Munson PJ, Rodbard D: LIGAND: a versatile computerized approach for characterization of ligand-binding systems. Anal Biochem 1980, 107: 220-239. 10.1016/0003-2697(80)90515-1.

    Article  CAS  PubMed  Google Scholar 

  50. Kashiwagi A, Verso MA, Andrews J, Vazquez B, Reaven G, Foley JE: In vitro insulin resistance of human adipocytes isolated from subjects with noninsulin-dependent diabetes mellitus. J Clin Invest 1983, 72: 1246-1254. 10.1172/JCI111080.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We are most grateful to Prof. B. J. Goldstein, for providing the rat IR cDNA probe. This work was supported by Funds from the Comunidad de Madrid/Universidad Complutense (14134/2005) and the Ministerio de Sanidad y Consumo (P1061744/2006).

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Consuelo Calle, Begoña Maestro or Moisés García-Arencibia.

Additional information

Authors' contributions

CC conceived the concept and design of the study, carried out the binding and activity studies, analyzed the computer data of DNA sequences and provided drafting of the article. BM carried out the rat models and the determinations in plasma and urine of these animals and performed the RNA blot assays. She also supplied statistical expertise, collected and assembled the data. MG-A participated in the identification by computer analysis of potential VDRE sequences and provided computer support. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Calle, C., Maestro, B. & García-Arencibia, M. Genomic actions of 1,25-dihydroxyvitamin D3 on insulin receptor gene expression, insulin receptor number and insulin activity in the kidney, liver and adipose tissue of streptozotocin-induced diabetic rats. BMC Molecular Biol 9, 65 (2008). https://doi.org/10.1186/1471-2199-9-65

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2199-9-65

Keywords