No effect of 25-hydroxyvitamin D supplementation on the skeletal muscle transcriptome in vitamin D deficient frail older adults

Roland W J Hangelbroek, Anouk M M Vaes, Mark V Boekschoten, Lex B Verdijk, Guido J E J Hooiveld, Luc J C van Loon, Lisette C P G M de Groot, Sander Kersten, Roland W J Hangelbroek, Anouk M M Vaes, Mark V Boekschoten, Lex B Verdijk, Guido J E J Hooiveld, Luc J C van Loon, Lisette C P G M de Groot, Sander Kersten

Abstract

Objective: Vitamin D deficiency is common among older adults and has been linked to muscle weakness. Vitamin D supplementation has been proposed as a strategy to improve muscle function in older adults. The aim of this study was to investigate the effect of calcifediol (25-hydroxycholecalciferol) on whole genome gene expression in skeletal muscle of vitamin D deficient frail older adults.

Methods: A double-blind placebo-controlled trial was conducted in vitamin D deficient frail older adults (aged above 65), characterized by blood 25-hydroxycholecalciferol concentrations between 20 and 50 nmol/L. Subjects were randomized across the placebo group and the calcifediol group (10 μg per day). Muscle biopsies were obtained before and after 6 months of calcifediol (n = 10) or placebo (n = 12) supplementation and subjected to whole genome gene expression profiling using Affymetrix HuGene 2.1ST arrays.

Results: Expression of the vitamin D receptor gene was virtually undetectable in human skeletal muscle biopsies, with Ct values exceeding 30. Blood 25-hydroxycholecalciferol levels were significantly higher after calcifediol supplementation (87.3 ± 20.6 nmol/L) than after placebo (43.8 ± 14.1 nmol/L). No significant difference between treatment groups was observed on strength outcomes. The whole transcriptome effects of calcifediol and placebo were very weak, as indicated by the fact that correcting for multiple testing using false discovery rate did not yield any differentially expressed genes using any reasonable cut-offs (all q-values ~ 1). P-values were uniformly distributed across all genes, suggesting that low p-values are likely to be false positives. Partial least squares-discriminant analysis and principle component analysis was unable to separate treatment groups.

Conclusion: Calcifediol supplementation did not significantly affect the skeletal muscle transcriptome in frail older adults. Our findings indicate that vitamin D supplementation has no effects on skeletal muscle gene expression, suggesting that skeletal muscle may not be a direct target of vitamin D in older adults.

Trial registration: This study was registered at clinicaltrials.gov as NCT02349282 on January 28, 2015.

Keywords: 25-hydroxyvitamin D; Vitamin D; older adults; skeletal muscle; transcriptomics.

Conflict of interest statement

LdG declares to have filed a patent related to vitamin D and cognitive executive function. All other authors have nothing to declare.

Figures

Fig. 1
Fig. 1
Very low expression of VDR in human skeletal muscle. qPCR Amplification curves of VDR in muscle biopsies. cDNA obtained from human muscle biopsies was PCR amplified using primers against human VDR. The different lines represent different subjects
Fig. 2
Fig. 2
Minimal effect of calcifediol and placebo on skeletal muscle gene expression in vitamin D deficient frail older adults. Volcano plot showing the relation between signal log ratio (log2[fold-change], x-axis) and the -log10 of the IBMT P-value (y-axis) for the effect of calcifediol and placebo on the skeletal muscle transcriptome. The dotted lines show the threshold for fold-change of 2, and P-value of 0.001
Fig. 3
Fig. 3
Minimal effect of calcifediol and placebo on skeletal muscle gene expression in vitamin D deficient frail older adults. Top left (a): P-values for the change in gene expression for all genes after filtering; before and after calcifediol supplementation. Uniform distribution (i.e. no increased frequency of genes for lower p-values) indicates an absence of an effect. Top right (b): values for the change in gene expression; before and after in the placebo group. Bottom left (c): P-values for the interaction effect (change in the calcifediol group vs. change in placebo group). Bottom right (d): receiver operator characteristic curve for the PLS-DA model
Fig. 4
Fig. 4
Multilevel principle component analysis (PCA) of gene expression data could not distinguish the four experimental groups. Sample plot from an unsupervised multilevel PCA. In multilevel PCA the paired structure of the data, i.e. measurements performed before and after intervention on the same subject, is taken into account to eliminate the inherent between-subject variation. Results revealed no separation between the experimental groups
Fig. 5
Fig. 5
Lack of effect of calcifediol supplementation on expression of putative VDR-dependent genes. Heatmap of gene expression changes (signal log ratio, SLR) of putative genomic targets of VDR
Fig. 6
Fig. 6
Partially overlapping effect of calcifediol and placebo on muscle transcriptome at the pathway level. Gene set enrichment analysis was performed separately for the effect of calcifediol treatment and placebo treatment. Venn diagram showing overlap in upregulated gene sets (FDR q-value

References

    1. Nair KS. Aging muscle. Am J Clin Nutr. 2005;81(5):953–963. doi: 10.1093/ajcn/81.5.953.
    1. Boucher BJ. The problems of vitamin d insufficiency in older people. Aging Dis. 2012;3(4):313–329.
    1. Casey C, Woodside JV, McGinty A, Young IS, McPeake J, Chakravarthy U, et al. Factors associated with serum 25-hydroxyvitamin D concentrations in older people in Europe: the EUREYE study. Eur J Clin Nutr. 2019;73(2):319–328. doi: 10.1038/s41430-018-0353-1.
    1. Tajar A, Lee DM, Pye SR, O'Connell MDL, Ravindrarajah R, Gielen E, et al. The association of frailty with serum 25-hydroxyvitamin D and parathyroid hormone levels in older European men. Age Ageing. 2013;42(3):352–359. doi: 10.1093/ageing/afs162.
    1. Bischoff-Ferrari HA, Dawson-Hughes B, Willett WC, Staehelin HB, Bazemore MG, Zee RY, et al. Effect of vitamin D on falls: a meta-analysis. J Am Med Assoc. 2004;291(16):1999–2006. doi: 10.1001/jama.291.16.1999.
    1. Visser M, Deeg DJH, Lips P. Low vitamin D and high parathyroid hormone levels as determinants of loss of muscle strength and muscle mass (sarcopenia): the Longitudinal Aging Study Amsterdam. J Clin Endocrinol Metab. 2003;88(12):5766–5772. doi: 10.1210/jc.2003-030604.
    1. Bischoff-Ferrari HA, Dietrich T, Orav EJ, Hu FB, Zhang Y, Karlson EW, et al. Higher 25-hydroxyvitamin D concentrations are associated with better lower-extremity function in both active and inactive persons aged ≥60 y. Am J Clin Nutr. 2004;80(3):752–758. doi: 10.1093/ajcn/80.3.752.
    1. Houston DK, Cesari M, Ferrucci L, Cherubini A, Maggio D, Bartali B, et al. Association between vitamin D status and physical performance: the InCHIANTI study. J Gerontol A Biol Sci Med Sci. 2007;62(4):440–446. doi: 10.1093/gerona/62.4.440.
    1. Tanner SB, Harwell SA. More than healthy bones: a review of vitamin D in muscle health. Ther Adv Musculoskelet Dis. 2015;7(4):152–159. doi: 10.1177/1759720X15588521.
    1. Boland R. Role of vitamin d in skeletal muscle function. Endocr Rev. 1986;7(4):434–448. doi: 10.1210/edrv-7-4-434.
    1. Boland R, de Boland AR, Marinissen MJ, Santillan G, Vazquez G, Zanello S. Avian muscle cells as targets for the secosteroid hormone 1,25-dihydroxy-vitamin D3. Mol Cell Endocrinol. 1995;114(1–2):1–8. doi: 10.1016/0303-7207(95)03650-V.
    1. Girgis CM, Clifton-Bligh RJ, Mokbel N, Cheng K, Gunton JE. Vitamin D signaling regulates proliferation, differentiation, and Myotube size in C2C12 skeletal muscle cells. Endocrinology. 2014;155(2):347–357. doi: 10.1210/en.2013-1205.
    1. Costa EM, Blau HM, Feldman D. 1,25-dihydroxyvitamin d3 receptors and hormonal responses in cloned human skeletal muscle cells. Endocrinology. 1986;119(5):2214–2220. doi: 10.1210/endo-119-5-2214.
    1. Simpson RU, Thomas GA, Arnold AJ. Identification of 1,25-dihydroxyvitamin D3 receptors and activities in muscle. J Biol Chem. 1985;260(15):8882–8891.
    1. Srikuea R, Zhang X, Park-Sarge O-K, Esser KA. VDR and CYP27B1 are expressed in C2C12 cells and regenerating skeletal muscle: potential role in suppression of myoblast proliferation. Am J Phys Cell Phys. 2012;303(4):C396–C405. doi: 10.1152/ajpcell.00014.2012.
    1. Endo I, Inoue D, Mitsui T, Umaki Y, Akaike M, Yoshizawa T, et al. Deletion of vitamin D receptor gene in mice results in abnormal skeletal muscle Development with deregulated expression of Myoregulatory transcription factors. Endocrinology. 2003;144(12):5138–5144. doi: 10.1210/en.2003-0502.
    1. Olsson K, Saini A, Strömberg A, Alam S, Lilja M, Rullman E, et al. Evidence for vitamin D receptor expression and direct effects of 1α,25(OH)2D3 in human skeletal muscle precursor cells. Endocrinology. 2016;157(1):98–111. doi: 10.1210/en.2015-1685.
    1. Ceglia L, Harris SS. Vitamin D and its role in skeletal muscle. Calcif Tissue Int. 2013;92(2):151–162. doi: 10.1007/s00223-012-9645-y.
    1. Pike JW. Closing in on vitamin D action in skeletal muscle: early activity in muscle stem cells? Endocrinology. 2016;157(1):48–51. doi: 10.1210/en.2015-2009.
    1. Wang Y, DeLuca HF. Is the vitamin D receptor found in muscle? Endocrinology. 2011;152(2):354–363. doi: 10.1210/en.2010-1109.
    1. Vaes AMM, Tieland M, Toussaint N, Nilwik R, Verdijk LB, van Loon LJC, et al. Cholecalciferol or 25-hydroxycholecalciferol supplementation does not affect muscle strength and physical performance in Prefrail and frail older adults. J Nutr. 2018;148(5):712–720. doi: 10.1093/jn/nxy024.
    1. Okazaki R, Ozono K, Fukumoto S, Inoue D, Yamauchi M, Minagawa M, et al. Assessment criteria for vitamin D deficiency/insufficiency in Japan: proposal by an expert panel supported by the Research Program of Intractable Diseases, Ministry of Health, Labour and Welfare, Japan, the Japanese Society for Bone and Mineral Research and the Japan Endocrine Society [Opinion] J Bone Miner Metab. 2017;35(1):1–5. doi: 10.1007/s00774-016-0805-4.
    1. Fried LP, Tangen CM, Walston J, Newman AB, Hirsch C, Gottdiener J, et al. Frailty in older adults: evidence for a phenotype. J Gerontol Ser A Biol Med Sci. 2001;56(3):M146–MM57. doi: 10.1093/gerona/56.3.M146.
    1. Vaes AMM, Tieland M, de Regt MF, Wittwer J, van Loon LJC, de Groot L. Dose-response effects of supplementation with calcifediol on serum 25-hydroxyvitamin D status and its metabolites: a randomized controlled trial in older adults. Clin Nutr. 2018;37(3):808–814. doi: 10.1016/j.clnu.2017.03.029.
    1. Bergstrom J. Percutaneous needle biopsy of skeletal muscle in physiological and clinical research. Scand J Clin Lab Invest. 1975;35(7):609–616. doi: 10.3109/00365517509095787.
    1. R Development Core Team. R . A Language and Environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2015.
    1. Bates D, Mächler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015;67(1):48. doi: 10.18637/jss.v067.i01.
    1. Lin K, Kools H, de Groot PJ, Gavai AK, Basnet Ram K, Cheng F, et al. MADMAX – Management and analysis database for multiple ~omics experiments. Journal of Integrative Bioinformatics. 2011;8(2):160.
    1. Bolstad BM, Irizarry RA, Astrand M, Speed TP. A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003;19(2):185–193. doi: 10.1093/bioinformatics/19.2.185.
    1. Dai M, Wang P, Boyd AD, Kostov G, Athey B, Jones EG, et al. Evolving gene/transcript definitions significantly alter the interpretation of GeneChip data. Nucleic Acids Res. 2005;33(20):e175.
    1. Piccolo SR, Sun Y, Campbell JD, Lenburg ME, Bild AH, Johnson WE. A single-sample microarray normalization method to facilitate personalized-medicine workflows. Genomics. 2012;100(6):337–344. doi: 10.1016/j.ygeno.2012.08.003.
    1. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015;43(7):e47.
    1. Sartor MA, Tomlinson CR, Wesselkamper SC, Sivaganesan S, Leikauf GD, Medvedovic M. Intensity-based hierarchical Bayes method improves testing for differentially expressed genes in microarray experiments. BMC Bioinformatics. 2006;7(1):1–17. doi: 10.1186/1471-2105-7-538.
    1. Chen EY, Tan CM, Kou Y, Duan Q, Wang Z, Meirelles GV, et al. Enrichr: interactive and collaborative HTML5 gene list enrichment analysis tool. BMC Bioinformatics. 2013;14(1):128. doi: 10.1186/1471-2105-14-128.
    1. Kuhn M. Building predictive models in R using the caret package. J Stat Softw. 2008;28(5):1–26. doi: 10.18637/jss.v028.i05.
    1. Sing T, Sander O, Beerenwinkel N, Lengauer T. ROCR: visualizing classifier performance in R. Bioinformatics. 2005;21(20):3940–3941. doi: 10.1093/bioinformatics/bti623.
    1. Gu Z, Eils R, Schlesner M. Complex heatmaps reveal patterns and correlations in multidimensional genomic data. Bioinformatics. 2016;32(18):2847–2849. doi: 10.1093/bioinformatics/btw313.
    1. Girgis CM, Mokbel N, Cha KM, Houweling PJ, Abboud M, Fraser DR, et al. The vitamin D receptor (VDR) is expressed in skeletal muscle of male mice and modulates 25-Hydroxyvitamin D (25OHD) uptake in Myofibers. Endocrinology. 2014;155(9):3227–3237. doi: 10.1210/en.2014-1016.
    1. Saksa N, Neme A, Ryynänen J, Uusitupa M, de Mello VDF, Voutilainen S, et al. Dissecting high from low responders in a vitamin D3 intervention study. J Steroid Biochem Mol Biol. 2015;148:275–282. doi: 10.1016/j.jsbmb.2014.11.012.
    1. Hassan-Smith ZK, Jenkinson C, Smith DJ, Hernandez I, Morgan SA, Crabtree NJ, et al. 25-hydroxyvitamin D3 and 1,25-dihydroxyvitamin D3 exert distinct effects on human skeletal muscle function and gene expression. PLoS One. 2017;12(2):e0170665. doi: 10.1371/journal.pone.0170665.
    1. Vaes AMM, Tieland M, Toussaint N, Nilwik R, Verdijk LB, van Loon LJC, et al. The effect of calcifediol or vitamin D3 supplementation on muscle strength and physical performance in pre-frail and frail older adults: a randomized placebo-controlled trial. In preparation. 2017.
    1. Tuoresmäki P, Väisänen S, Neme A, Heikkinen S, Carlberg C. Patterns of genome-wide VDR locations. PLoS One. 2014;9(4):e96105. doi: 10.1371/journal.pone.0096105.
    1. Zanello SB, Collins ED, Marinissen MJ, Norman AW, Boland RL. Vitamin D receptor expression in chicken muscle tissue and cultured myoblasts. Horm Metab Res. 1997;29(5):231–236. doi: 10.1055/s-2007-979027.
    1. Bischoff HA, Borchers M, Gudat F, Duermueller U, Theiler R, Stähelin HB, et al. In situ detection of 1,25-dihydroxyvitamin D receptor in human skeletal muscle tissue. Histochem J. 2001;33(1):19–24. doi: 10.1023/A:1017535728844.
    1. Ceglia L, da Silva Morais M, Park LK, Morris E, Harris SS, Bischoff-Ferrari HA, et al. Multi-step immunofluorescent analysis of vitamin D receptor loci and myosin heavy chain isoforms in human skeletal muscle. J Mol Histol. 2010;41(2–3):137–142. doi: 10.1007/s10735-010-9270-x.
    1. Girgis CM, Clifton-Bligh RJ, Hamrick MW, Holick MF, Gunton JE. The roles of vitamin D in skeletal muscle: form, function, and metabolism. Endocr Rev. 2013;34(1):33–83. doi: 10.1210/er.2012-1012.
    1. Schubert L, DeLuca HF. Hypophosphatemia is responsible for skeletal muscle weakness of vitamin D deficiency. Arch Biochem Biophys. 2010;500(2):157–161. doi: 10.1016/j.abb.2010.05.029.

Source: PubMed

3
Prenumerera