Research Paper Volume 10, Issue 6 pp 1338—1355

Chronic vitamin D insufficiency impairs physical performance in C57BL/6J mice

Kenneth L. Seldeen1, , Manhui Pang1, , Merced M. Leiker1, , Jonathan E. Bard2, , Maria Rodríguez-Gonzalez1, , Mireya Hernandez1, , Zachary Sheridan1, , Norma Nowak2, , Bruce R. Troen1, ,

  • 1 Division of Geriatrics and Palliative Medicine, Jacobs School of Medicine and Biomedical Sciences, University at Buffalo and Research Service, Veterans Affairs Western New York Healthcare System, Buffalo, NY 14203, USA
  • 2 New York State Center of Excellence in Bioinformatics and Life Sciences and Department of Biochemistry, School of Medicine and Biomedical Sciences, University at Buffalo, Buffalo, NY 14214, USA

Received: April 27, 2018       Accepted: June 4, 2018       Published: June 14, 2018      

https://doi.org/10.18632/aging.101471
How to Cite

Copyright: © 2018 Seldeen et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 3.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Abstract

Vitamin D insufficiency (serum 25-OH vitamin D < 30 ng/ml) affects 70-80% of the general population, yet the long-term impacts on physical performance and the progression of sarcopenia are poorly understood. We therefore followed 6-month-old male C57BL/6J mice (n=6) consuming either sufficient (STD, 1000 IU) or insufficient (LOW, 125 IU) vitamin D3/kg chow for 12 months (equivalent to 20-30 human years). LOW supplemented mice exhibited a rapid decline of serum 25-OH vitamin D levels by two weeks that remained between 11-15 ng/mL for all time points thereafter. After 12 months LOW mice displayed worse grip endurance (34.6 ± 14.1 versus 147.5 ± 50.6 seconds, p=0.001), uphill sprint speed (16.0 ± 1.0 versus 21.8 ± 2.4 meters/min, p=0.0007), and stride length (4.4 ± 0.3 versus 5.1 ± 0.3, p=0.002). LOW mice also showed less lean body mass after 8 months (57.5% ± 5.1% versus 64.5% ± 4.0%, p=0.023), but not after 12 months of supplementation, as well as greater protein expression of atrophy pathway gene atrogin‑1. Additionally, microRNA sequencing revealed differential expression of mIR‑26a in muscle tissue of LOW mice. These data suggest chronic vitamin D insufficiency may be an important factor contributing to functional decline and sarcopenia.

Introduction

Sarcopenia affects between 10% and 35% of individuals over the age of 65 and underlies functional decline and increases fall risk in older individuals [1,2]. Although the underlying causes of sarcopenia remain unknown, sarcopenia is associated with vitamin D insufficiency (25-OH vitamin D >10 ng/ml and < 30 ng/ml) [3,4]. Vitamin D insufficiency is also prevalent worldwide, affecting as many as 70% of the population across all age demographics [5,6], raising the likelihood that some individuals may be vitamin D insufficient for decades or more.

Vitamin D insufficiency has been associated with a range of physical performance impairments across multiple studies, including: poor grip and leg strength, lower scores in the short physical performance battery (SPPB), low activities of daily living (ADL) scores, and reduced physical activity in community dwelling older individuals [3,7,8]. Low serum 25-OH vitamin D (< 12.5 ng/ml) was also associated with worse physical performance as a composite score of gait speed, chair-stand, and grip strength in the Women’s Health Initiative (WHI) study [9]. Low serum vitamin D (<12.0 ng/ml) has also been associated with poor physical performance, fracture risk, and fractures in the Longitudinal Aging Study Amsterdam (LASA) [10]. However, several studies involving older individuals have reported no association between vitamin D and physical performance [11-13], raising uncertainty regarding the role of vitamin D in physical performance.

Animal studies can provide insights free from genetic and lifestyle factors that might confound the interpretation of human studies, leading to conflicting reports. Vitamin D receptor (VDR) knockout mice exhibit declines in several aspects of performance including swimming endurance, rotarod, and open field activity [14-16]. Furthermore, VDR ablation in muscle stem cells has demonstrated its critical role in muscle metabolism including muscle anabolic signaling [17,18], satellite cell proliferation [19], and regulation of the uptake of 25-OH vitamin D into muscle cells [20]. Recently a study silencing VDR expression with siRNA in muscle stem cells (C2C12) demonstrated inhibition of myogenic differentiation [21]. Dietary elimination of vitamin D has also been reported to impair performance, including lesser grip strength in mice [22] and swim endurance [23], while 1α-OH vitamin D treatment in ovariectomized rats was associated with increases in muscle strength, but not with muscle fatigue [24].

Animal studies rarely include the physiological relevant condition of vitamin D insufficiency, opting for vitamin D receptor knockouts or complete dietary elimination. Furthermore, there are no rodent or human studies that follow the impacts of vitamin D insufficiency for long periods of time (≥1 year for mice or about 20-30 years for humans), particularly while investigating physical performance. Therefore, we established groups of mice at serum vitamin D sufficient and insufficient levels, and followed changes in body composition and physical performance over a 12‑month period. Our findings demonstrate that vitamin D insufficiency modulates muscle miRNA signaling, increases atrophy pathway proteins, and impairs specific aspects of physical performance.

Results

Depletion and repletion of serum 25-OH vitamin D occurs rapidly in response to changes in the amount of vitamin D3 supplementation, and levels remain consistent relative to the amount of supplementation

To understand the potential impacts of chronic vitamin D insufficiency we supplemented 6‑month old male C57BL6/J mice with standard (STD) facility levels of vitamin D3 in chow (1000 IU vitamin D3/kg) or reduced to a lower amount (LOW, 125 IU vitamin D3/kg) to induce vitamin D insufficiency. We observed that the STD level of supplementation results in stable serum 25‑OH vitamin D levels (30‑40 ng/ml) (Figure 1A, black line). We further observed that LOW supplementation leads to a rapid decline in serum 25‑OH vitamin D, reaching human equivalent levels of vitamin D insufficiency after just two weeks, and remaining consistently between 10-15 ng/ml for the remainder of the experiment (Figure 1A, dark gray line). After 4 months, we tested the rate of vitamin D repletion by switching a group of mice receiving 125 IU vitamin D3/kg to 1000 IU vitamin D3/kg chow (Figure 1A, light gray line), and found repletion also occurs within two weeks with resultant serum 25‑OH vitamin D levels in these mice similar to levels in the STD group. After 12 months, there were no statistically significant differences in serum calcium or levels of the active metabolite of vitamin D, serum 1,25‑(OH)2 vitamin D between the STD and LOW groups (Figure 1B & 1C, respectively). We did observe a trend towards higher serum intact parathyroid hormone in vitamin D insufficient mice (STD: 175.2 ± 18.4 pg/ml versus LOW: 279.7 ± 125.4 pg/ml, p=0.10, Figure 1D).

Longitudinal analysis of 25-OH vitamin D in mice given differential vitamin D3 supplementation. Six-month-old C57BL/6J male mice were given 1000 IU vitamin D3/kg chow (STD, black line, n=6), 125 IU (LOW, dark gray line, n=6), or were changed from 125 IU to 1000 IU after 2 months (light gray line, n=8) (A). After 12 months serum calcium (B), 1,25 (OH)2 vitamin D (C), and intact PTH (D), were assessed in STD and LOW mice using calorimetric assays or ELISA. Intact PTH was not significantly different (ns, p=0.10).

Figure 1. Longitudinal analysis of 25-OH vitamin D in mice given differential vitamin D3 supplementation. Six-month-old C57BL/6J male mice were given 1000 IU vitamin D3/kg chow (STD, black line, n=6), 125 IU (LOW, dark gray line, n=6), or were changed from 125 IU to 1000 IU after 2 months (light gray line, n=8) (A). After 12 months serum calcium (B), 1,25 (OH)2 vitamin D (C), and intact PTH (D), were assessed in STD and LOW mice using calorimetric assays or ELISA. Intact PTH was not significantly different (ns, p=0.10).

Chronic vitamin D insufficiency did not affect body weight, but accelerates the reduction in lean body mass and the increase in fat mass

We examined the impacts of chronic vitamin D insufficiency on body weight and composition. We found body weights to be similar between the two groups at all time points (Figure 2A), with equivalent overall weight gains (STD: 41.3% ± 10.6% versus LOW: 44.7% ± 15.6%, p=0.66), which is expected for mice of this age. Additionally, we examined lean body mass, fat mass, and bone mineral density using dual X-ray absorptiometry (DEXA). After 4 months of treatment (10 months of age), vitamin D insufficient mice trended towards a lower lean body mass and greater fat mass (p=0.08 for both, Figures 2B and 2C) and was significantly different from STD mice after 8 months (14 months of age, lean mass - STD: 64.5% ± 4.0% versus LOW: 57.5% ± 5.1%, p=0.0231 and fat mass - STD: 35.4% ± 4.0% versus LOW: 42.5% ± 5.2%, p=0.0243). However, there was no difference in either parameter after 12 months of treatment (18 months of age) (p=0.37). Bone mineral density of LOW mice was significantly lower than STD mice after 8 months (STD: 54.8 ± 0.7 mg/cm2 versus LOW: 52.8 ± 1.3 mg/cm2, p=0.0078, Figure 2D), but significant differences were not observed after 12 months.

Analysis of body composition and bone mineral density. Body weight was measured every two weeks in mice given standard (STD) and low (LOW) vitamin D3 in chow (A). Four months following the initiation of treatment lean mass (B), fat mass (C), and bone mineral density (D) were assessed using dual X-ray absorptiometry, n=6 per group, “*” indicates p

Figure 2. Analysis of body composition and bone mineral density. Body weight was measured every two weeks in mice given standard (STD) and low (LOW) vitamin D3 in chow (A). Four months following the initiation of treatment lean mass (B), fat mass (C), and bone mineral density (D) were assessed using dual X-ray absorptiometry, n=6 per group, “*” indicates p < 0.05.

Vitamin D insufficiency impairs physical performance across multiple domains

To examine whether vitamin D insufficiency impacts physical performance, and in what ways, we performed a range of behavioral assessments (Figure 3). We observed that grip strength did not significantly change over the 12‑month period of our experiment in either cohort (Figure 3A). There was also no difference in rotarod fall latency between supplementation groups, an assessment of balance and coordination, although both groups exhibited declines with age (STD baseline time to fall: 272.0 ± 45.4 seconds versus endpoint: 183.6 ± 39.1 seconds, paired T-test: p=0.014; LOW baseline: 231.7 ± 36.6 seconds versus endpoint: 157.8 ± 54.7 seconds, p=0.005, Figure S1A). However, vitamin D insufficient mice showed a deficiency in grip endurance, as determined by both grip wire (STD 65.8 ± 18.6 seconds versus LOW: 35.4 ± 6.7 seconds, n = 6 and 5, respectively, p=0.0039, Figure 3B) and grip grid latency (STD 147.5 ± 50.6 seconds versus LOW: 34.6 ± 14.1 seconds, n = 6 and 5, respectively, p=0.001, Figure S1B). No difference was observed in treadmill endurance (Figure 3C), although a difference between the two groups was trending (p=0.11 at 8 months and 0.06 at 12 months). The treadmill assessment was performed at 0˚ inclination, and gradual increases in speed likely imposed greater utilization of anaerobic respiration. However, to more completely examine if vitamin D insufficiency impairs anaerobic performance, we devised an uphill treadmill assessment by setting the inclination to 25˚ and introducing periods of active recovery at low speed to maximize anaerobic response. We observed that after 48 weeks, vitamin D sufficient mice achieved greater time before exhaustion than did vitamin D insufficient mice (STD: 9.4 ± 1.7 mins versus LOW: 5.1 ± 0.8 m/min, p=0.0007, n = 6 and 5 respectively, Figure 3D).

Physical performance in vitamin D sufficient and insufficient mice. Vitamin D sufficient (STD) and insufficient (LOW) mice were assessed across a range of physical performance domains that include: grip strength assessed every 4 months as the best 3 of 5 trials on a grip strength meter, n=6 (A); grip endurance as the best of two trials timed for latency to fall from a wire, n=5 (B); aerobic endurance assessed as a single trial for time before exhaustion on a mouse treadmill, n=6 (C); anaerobic endurance assessed as a single trial of increasing intensity intervals on an inclined (25º) mouse treadmill, n=6, 5 respectively (D); exploratory behavior as a count of quadrant crossings and rearings over 5 minutes in an open field arena, n=6, 5 respectively (E); and gait as assessed by measurement of stride length normalized to femur length determined using dual X-ray absorptiometry, n=5, 6 respectively (F). Statistical significance indicated by “*” p “ns” indicating non‑significance.

Figure 3. Physical performance in vitamin D sufficient and insufficient mice. Vitamin D sufficient (STD) and insufficient (LOW) mice were assessed across a range of physical performance domains that include: grip strength assessed every 4 months as the best 3 of 5 trials on a grip strength meter, n=6 (A); grip endurance as the best of two trials timed for latency to fall from a wire, n=5 (B); aerobic endurance assessed as a single trial for time before exhaustion on a mouse treadmill, n=6 (C); anaerobic endurance assessed as a single trial of increasing intensity intervals on an inclined (25º) mouse treadmill, n=6, 5 respectively (D); exploratory behavior as a count of quadrant crossings and rearings over 5 minutes in an open field arena, n=6, 5 respectively (E); and gait as assessed by measurement of stride length normalized to femur length determined using dual X-ray absorptiometry, n=5, 6 respectively (F). Statistical significance indicated by “*” p < 0.05, “**” p< 0.01, “***” p<0.001, and “ns” indicating non‑significance.

We further analyzed functional capacity by assessing open field activity. Although we did not observe a difference in open field exploration between the two groups (open field quadrant crossings: STD: 23.5 ± 12.4 versus LOW: 16.5 ± 11.0, p=0.33, n = 6 and 5, respectively, Figure 3E), we observed that vitamin D insufficient mice reared less often (STD: 13.8 ± 6.6 versus LOW: 4.4 ± 2.1, p=0.0141, n = 6 and 5, respectively, Figure 3E). Additionally, vitamin D insufficient mice exhibited significantly shorter stride length at 8 months (p=0.0318) and at 12 months (p=0.0018) (Figure 3F) and these mice also exhibited further decline in stride length from 8 to 12 months (LOW: 4.9 ± 0.3 to 4.4 ± 0.3 stride : femur length, p=0.0026, n=6), which was not observed in vitamin D sufficient mice (STD: 5.3 ± 0.2 to 5.1 ± 0.2, p=0.1757, n=5).

Vitamin D insufficient mice exhibit trends for lower muscle fiber size and myofibrillar protein content, in addition to greater protein expression of atrophy associated atrogin-1

We next sought to identify factors that may be driving the observed differences in physical performance. NADH histological staining on quadriceps muscle (Figure 4A) showed a trend towards smaller cross- sectional area (CSA) in light stained fibers in vitamin D insufficient mice (LOW: 3,057.1 ± 708.3 versus STD: 4,046.7 ± 977.1 μm2 p=0.13, n=4 and 5, respectively,Figure 4B). We further examined if differences occurred in the content of specific muscle fractions by isolating myofibrillar, sarcoplasmic, and mitochondrial proteins using differential buffers and centrifugation (Figure 4C). We did not observe differences in sarcoplasmic protein content, yet there was a trend towards lower myofibrillar protein content in vitamin D insufficient mice (LOW: 17.4 ± 6.3 μg/mg tissue versus STD: 23.1 ± 4.9 μg/mg, p=0.11). Although these data are only suggestive, they are consistent with studies that report muscle decline in vitamin D receptor knockout models [25,26]. As aberrant vitamin D signaling is associated with greater atrophy pathway signaling [27], we investigated the expression of atrogin-1 and found greater expression in the quadriceps muscles of the vitamin D insufficient mice (STD: 2.33 ± 0.52 versus LOW: 2.88 ± 0.24, p=0.0393, Figure 4D). We did not observe differences in mitochondrial protein content (p=0.38, Figure 4C), the ratio of mitochondrial: nuclear DNA content in soleus muscle (STD: 0.99 ± 0.06 versus LOW: 0.89 ± 0.17, p=0.18, Figure 4E), or mitochondrial complex IV activity (p=0.82, Figure 4F).

Analysis of muscle quality in vitamin D sufficient and insufficient mice. To assess impacts of vitamin D insufficiency on muscle biology, tissues were harvested following 12 months of sufficient (STD) or insufficient (LOW) supplementation. Quadriceps muscle was then analyzed with NADH staining (A) allowing for quantification of the cross-sectional area (CSA) of light stain fibers that corresponds to fast twitch fibers, n=5, 4 respectively (B). Gastrocnemius muscle was analyzed by differential centrifugation to determine myofibrillar, sarcoplasmic, and mitochondrial protein content, n=6 (C). Atrogin-1 expression was also determined in gastrocnemius muscle with western blotting and relative expression (atrogin / tubulin) was quantified using ImageJ software, n=6 (D). Mitochondrial biomass (E) and activity (F) in soleus muscle, n=6, were determined by quantitative PCR and biochemical assays, respectively. Statistical significance indicated by “*” p “ns” indicating non‑significance.

Figure 4. Analysis of muscle quality in vitamin D sufficient and insufficient mice. To assess impacts of vitamin D insufficiency on muscle biology, tissues were harvested following 12 months of sufficient (STD) or insufficient (LOW) supplementation. Quadriceps muscle was then analyzed with NADH staining (A) allowing for quantification of the cross-sectional area (CSA) of light stain fibers that corresponds to fast twitch fibers, n=5, 4 respectively (B). Gastrocnemius muscle was analyzed by differential centrifugation to determine myofibrillar, sarcoplasmic, and mitochondrial protein content, n=6 (C). Atrogin-1 expression was also determined in gastrocnemius muscle with western blotting and relative expression (atrogin / tubulin) was quantified using ImageJ software, n=6 (D). Mitochondrial biomass (E) and activity (F) in soleus muscle, n=6, were determined by quantitative PCR and biochemical assays, respectively. Statistical significance indicated by “*” p < 0.05 and “ns” indicating non‑significance.

Vitamin D insufficient and sufficient mice exhibit a similar inflammatory profile

Both clinical [28] and cell culture studies [29,30] support a role for vitamin D in modulating inflammatory cytokines. We therefore set out to determine if 12 months of vitamin D insufficiency modulated serum cytokine levels (Figure 5A, Table 1). Surprisingly, we did not identify any serum cytokine concentrations as being significantly different due to treatment, which included IL-1α, IL‑1β, IL-6, IL‑10, IL‑15, IL-18, MCP, and TNFα. We note that the serum level of IL‑18 was trending lower in vitamin D insufficient mice (STD: 118.0 ± 9.9 pg/ml versus LOW: 109.6 ± 9.2 pg/ml, p=0.16). We further investigated tissue concentrations of IL‑6 in brain, heart, and epididymal adipose tissue, but did not find any elevation in these tissues (Figure 5B). However, we also note that adipose tissue IL‑6 trended higher in vitamin D insufficient mice (STD: 251.7 ± 75.5 ng/μg protein versus LOW: 425.5 ± 213.6 ng/μg protein, p=0.09).

Inflammatory profile of vitamin D sufficient and insufficient mice. Serum at endpoint and used for multiplex ELISA analysis of serum cytokines, n=6 (A), ns denotes non-significance between vitamin D sufficient (STD) and insufficient (LOW) mice (IL-18: p=0.16). Interleukin-6 content was further assessed in brain, heart, and adipose tissues, n=6 (B), ns denotes non-significance in adipose IL-6 content (p=0.09).

Figure 5. Inflammatory profile of vitamin D sufficient and insufficient mice. Serum at endpoint and used for multiplex ELISA analysis of serum cytokines, n=6 (A), ns denotes non-significance between vitamin D sufficient (STD) and insufficient (LOW) mice (IL-18: p=0.16). Interleukin-6 content was further assessed in brain, heart, and adipose tissues, n=6 (B), ns denotes non-significance in adipose IL-6 content (p=0.09).

Table 1. Inflammatory cytokines in vitamin D sufficient and insufficient mice. Serum cytokines were determined after 12 months of STD or LOW vitamin D supplementation using multiplex ELISA.

CytokineSTDLOWp‑value
IL-1α2.5 ± 1.02.0 ± 0.30.27
IL-1β25.1 ± 4.922.7 ± 5.30.44
IL-61.0 ± 0.30.9 ± 0.30.37
IL-100.8 ± 0.20.9 ± 0.10.44
MCP7.0 ± 2.46.6 ± 1.20.73
TNFα6.1 ± 1.05.7 ± 0.70.48
IL-1550.9 ± 14.851.0 ± 10.80.99
IL-18118.0 ± 9.9109.6 ± 9.20.16

Vitamin D insufficient mice exhibit a distinct muscle miRNA profile

MicroRNAs are powerful effectors underlying many physiological processes, including the pathophysiology of muscle mass decline [31]. To investigate whether vitamin D insufficiency modulates muscle miRNA signaling, we isolated RNA from tibialis anterior muscles and then prepared samples for Next-Gen RNA sequencing using an Illumina NextSeq 500 high-throughput sequencing system. RNA sequencing revealed an average of 11.7 ± 1.5 million and 11.9 ± 1.4 million total reads per mouse for STD and LOW, respectively. Of these 4.8 ± 1.0 x 105 and 4.9 ± 1.2 x 105 reads for each STD and LOW, respectively, were mapped to known miRNAs using mIRBase v21. To avoid using insufficient data due to low expression, we removed any miRNA that did not reach at least 10 reads across 66% or more of the samples (per Jung et al. [32]), yielding a total of 202 unique miRNAs for further analysis. Our data revealed 12 miRNAs with potential for differential expression between STD and LOW (Figure 6 and Table 2, p<0.05); however, after false discovery rate correction we identified only one differentially expressed miRNA, mIR26a-5p, with a q‑value = 0.0392.

Analysis of differential miRNA expression in vitamin D sufficient and insufficient mice. Total RNA was isolated from tibialis anterior muscle and analyzed using miRNA sequencing. A total of 202 unique miRNAs reached thresholds for analysis, and a heat map was generated from 15 miRNA with a pn=6). Further analysis for false positives (Q-val) revealed mir-26a as being differentially expressed between the two groups.

Figure 6. Analysis of differential miRNA expression in vitamin D sufficient and insufficient mice. Total RNA was isolated from tibialis anterior muscle and analyzed using miRNA sequencing. A total of 202 unique miRNAs reached thresholds for analysis, and a heat map was generated from 15 miRNA with a p<0.05 for differential expression between vitamin D sufficient (STD) and insufficient (LOW) mice (n=6). Further analysis for false positives (Q-val) revealed mir-26a as being differentially expressed between the two groups.

Table 2. Differential miRNA expression in vitamin D sufficient and insufficient mice.

MicroRNAmIRBASE ID#Log2FCP-valueQ-value
miR-26aMI00005730.070.0000.04
miR-128MI00001550.320.0030.15
miR-30cMI00005470.170.0050.15
Let-7cMI00005590.170.0050.15
Let-7aMI00005560.170.0160.15
miR-10bMI00002210.250.0180.36
miR-92MI00007190.110.0250.41
Let-7bMI00005580.550.0300.44
miR-204MI00002470.210.0320.44
miR-139MI00006930.210.0480.56
miR-145aMI0000169-0.140.0360.46
miR-146aMI0000170-0.560.0210.38

Discussion

Vitamin D insufficiency is a prevalent condition for which the long-term impacts are poorly understood. Here we demonstrate that mice kept in a serum 25‑OH vitamin D insufficient state for 12 months exhibit significant declines across multiple domains of physical performance. These domains include grip hang endurance, uphill sprint, and open field rearing, which together may be indicative of a decline in anaerobic capacity due to vitamin D insufficiency. The loss of anaerobic capacity is partially supported by a trend (p=0.13) towards lower smaller fast twitch fiber CSA in vitamin D insufficient mice and lower lean mass in vitamin D insufficient mice, although this latter difference was of limited duration. Belenchia et al. [33] also observed loss of lean mass in mice due to dietary vitamin D deficiency (serum 25‑OH vitamin D < 10 ng/ml) of approximately 40 weeks in female mice initiated at 8 weeks of age. We suspect our study was underpowered to identify such histological differences; however, the possibility that vitamin D insufficient mice exhibit smaller fast twitch fiber CSA is supported by our finding that vitamin D insufficient mice also exhibit greater expression of atrogin-1, for which this and other atrophy associated proteins have previously been linked to vitamin D signaling [27,34,35]. Thus these data support the notion that long-term vitamin D insufficiency may contribute to the progression of sarcopenia. Our finding that vitamin D insufficient mice exhibit shorter stride lengths further emphasizes the potential contribution of vitamin D status, as such gait disturbances are an integral component of functional capacity decline associated with sarcopenia [36,37].

Other areas of physical performance were not impacted in the time frame of our study. Surprisingly, we did not identify differences in grip strength in our mice in contrast to multiple human studies reporting an association between vitamin D and strength [38-41], but not all [42,43]. However, we did observe a trend towards lower myofibrillar protein content (p=0.11), which may be indicative of future grip strength decline had we continued our experiment beyond 12 months. We also did not observe significance differences in treadmill performance, although vitamin D insufficient mice were trending towards lower treadmill performance at both 8 months (p=0.11) and 12 months (p=0.06). We believe this may be due to a greater anaerobic requirement as opposed to a possibility of aerobic deficit, as supported by our finding of differences when assessing treadmill performance using uphill intervals. Our findings that neither mitochondrial biomass nor activity were affected by vitamin D insufficiency further supports the idea that both groups exhibit similar aerobic capacity.

Additionally, although the relationship between vitamin D and fall risk is suggestive, it remains inconclusive [44]. We anticipated vitamin D insufficient mice would exhibit worse rotarod performance, as rotarod is an indicator of balance and coordination [45]. Yet our data show both groups exhibit age-dependent declines. These findings are surprising in light of Sakai et al. [46] who reported vitamin D receptor ablation impairs rotarod performance, perhaps indicating that stark vitamin D deficiency or receptor ablation is necessary or that such differences would have appeared if we continued the experiment into advanced ages (> 24 months of age).

We found long term impacts of serum 25‑OH vitamin D on body weight and body composition, which is consistent with our 6-month study examining the impacts of alterations of serum 25‑OH vitamin D levels in lean and obese mice [47]. However, Belenchia et al. reported declines in body weight in vitamin D deficient mice after 6 months until the endpoint of the study after 10 months, at which point these mice also exhibited decreased fat mass and lean body mass [33]. In contrast, our study identified a decrease in lean mass, and an increase in fat mass, after 8 months of insufficiency that was not observed after 12 months. However, the similarities and differences in our body composition findings compared to Belenchia et al. may be explained by the use of males versus females, insufficiency versus deficiency, and that our study continued two months longer. It was also surprising that our study failed to show differences in the inflammatory milieu in light of the reported relationships between serum vitamin D and chronic inflammation [28]. In particular, lower serum vitamin D was previously shown to be associated with increased IL‑6 expression [28,48,49], yet we did not observe any increases in serum, brain, or cardiac tissues. We did observe a trend of greater IL‑6 expression in adipose tissue of vitamin D insufficient mice (0.09), which would be consistent with mechanistic reports of the impacts of vitamin D in adipocytes [50].

Our analysis of miRNA sequencing revealed a single miRNA, miR-26, was differentially expressed in vitamin D insufficient mice. miR-26 has previously been shown to be important for skeletal muscle differentiation [51,52], and was found in two human studies to be differentially expressed in response to exercise [53,54]. Additionally, miR-26a was differentially expressed in two separate pilot studies that included 5 subjects supplemented with high dose vitamin D for a 12‑month period [55]. However, the authors further reported that mIR-26a was not differentially expressed in a larger scale study [55], indicating the need for additional studies to better elucidate the impact of vitamin D status on miRNA profiles and expression. Our study also identified other miRNAs with reported roles in skeletal muscle biology including miR-204 [56], miR-139 [57,58], miR-146 [58-60], and miR‑30 [58,61], although none of these met the threshold for false discovery rate (q < 0.05). We think this may be in part due to insufficient power, which may have also prevented us from identifying differences in other parameters (i.e. fiber size, inflammation, physical performance). Additionally, our low sequencing depth (of just under 500,000 mapped miRNA reads per sample) may have also restricted our ability to discern differences in only those miRNAs with high expression [62].

Our study also confirms the findings that altering vitamin D supplementation results in a rapid shift (both depletion and repletion) in serum 25-OH vitamin D levels (within 2 weeks) that is sustained relative to the amount of supplementation [33,47,63]. Our 12-month study was approximately 6 weeks longer than Belenchia et al. [33], which was performed in female mice, and together these studies demonstrate little to no impact of gender on the relationship between vitamin D supplementation and serum 25-OH vitamin D concentration. Interestingly, we did not observe differences in serum 1,25‑(OH)2 vitamin D between 25‑OH vitamin D sufficient and insufficient mice. Although this finding is consistent with our previous study when we induced vitamin D insufficiency for 6 months in mice [47], it remains surprising in light of the functional role of 1,25‑(OH)2 vitamin D in skeletal muscle regulation [64-66]. In particular, Hassan‑Smith et al. have reported that serum 1,25‑(OH)2 vitamin D concentrations correlated with leg power, velocity, and jump height in 20-74 year old males and females, parameters that were not found to be correlated with serum 25-OH vitamin D [67]. Serum 1,25‑(OH)2 vitamin D was also found to be correlated with low muscle mass in cross-sectional analyses of men and women aged 21-97 years, as well as knee extension force in women [11]. Additionally, both low serum 25‑OH and 1,25‑(OH)2 vitamin D were independently found to be associated with the incidence of sarcopenia at a 5 year follow-up in men >70 years of age [68]. Yet, Boonen et al. reported no association between 1,25‑(OH)2 vitamin D and knee extension strength in women aged 70 to 90 [69], and Gielen et al. also reported no association between 1,25‑(OH)2 vitamin D and physical performance, specifically grip strength and gait speed in men aged 70 and older [70]. These studies were observational and did not involve long-term abatement of dietary vitamin D3. It is therefore possible that the duration and the degree of serum 25-OH vitamin D reduction may affect outcomes. Another possibility is that serum measures of 1,25‑(OH)2 vitamin D in our year long study did not reflect the actual 1,25‑(OH)2 vitamin D within muscle since skeletal muscle expresses the 25‑hydroxyvitamin D3 1-alpha-hydroxylase [20,71].

With regards to other serum markers, Belenchia et al. reported changes in serum calcium levels in vitamin D deficient mice, which was not observed in this study, our previous study [47], or by Mallya et al. [63]. Likewise, Belenchia et al. was the only study to report significant differences in PTH concentrations; however, intact PTH was trending towards elevation in our vitamin D insufficient mice (p=0.10), and our study may have been insufficiently powered to observe this. Despite these differences, the behavior of serum 25-OH vitamin D in response to altered vitamin D supplementation was consistent between these studies and supports the use of the dietary supplementation and deprivation to examine serum 25-OH vitamin D related phenomena.

Conclusion

Serum 25-OH vitamin D declines rapidly and remains consistently depressed in response to low supplementation. Prolonged vitamin D insufficiency induces characteristics of sarcopenia that include poor anaerobic capacity, lower lean mass, and a trend towards smaller fast twitch fiber CSA, as well as gait disturbance. Vitamin D insufficient mice also exhibited increased expression of atrophy-associated Atrogin-1 and differential expression of muscle regulation associated miR-26a. These data suggest a role for chronic vitamin D insufficiency in the development of sarcopenia, highlighting the need for further animal and human studies to investigate the impacts of vitamin D during aging.

Materials and Methods

Animals

Twelve C57BL6/J mice (5 months old) were purchased from Jackson labs (Bar Harbor, ME). After 1 month the mice were randomly sorted into groups (n=6) that received AIN‑93G chow (Dyets Inc., Bethleham, PA) supplemented with either the standard facility amount of 1000 IU vitamin D3 / kg chow (STD) to maintain serum 25-OH vitamin D sufficiency or 125 IU vitamin D3 / kg chow (LOW) to induce serum vitamin D insufficiency over a period of 12 months (Table 3). Additionally, eight mice were initially supplemented with 125 IU, but then switched after two months to 1000 IU to examine the rate of serum 25-OH vitamin D repletion. Food and water were provided ad libitum, and mice were housed in large shoebox animal cages containing 6 or 8 mice per cage. Lighting was on a 12 hour on / 12 hour off cycle, and cages were shielded to reduce exposure to facility lighting. Body weight was measured every two weeks. All studies and experimental protocols were approved by and in compliance with guidelines of the Miami VA and VA Western New York Animal Care and Use Committees.

Table 3. Diet compositions.

ComponentContent (per kg chow)
Vitamin Free Casein200 g
L-Cystine3 g
Dyetrose35 g
Sucrose349 g
Cornstarch3.6 g
Soybean oil9 g
Lard20 g
t-Butylhydroquinone (TBHQ)0.005 g
Cellulose50 g
Dicalcium Phosphate13 g
Calcium Carbonate5.5 g
Potassium Citrate16.5 g
Choline Bitartrate2 g
Vitamin/Mineral Mix (Vit D free)20 g
Vitamin D3 (cholecalciferol)STD: 1000 IU (25 μg)
LOW: 125 IU (3.1μg)

ELISA and Colorimetric assays

Blood was collected through the sub-mandibular vein using a mouse lancet (MEDIpoint, inc., Mineola, NY.) into microcentrifuge tubes. Samples were held at room temperature for 10 minutes to allow coagulation and then centrifuged at 16,0000 x g for 10 minutes at 4 °C to allow separation of serum. Analysis of serum was performed using ELISA kits for 25-OH vitamin D (ImmunoDiagnostic Systems, Inc., Scottsdale, AZ), 1,25-(OH)2 vitamin D (MyBioSource, San Diego, CA), and intact PTH (MyBioSource). Colorimetric assays were performed to assess serum calcium concentration (Biovision, San Francisco, CA.) according to manufacturer protocols. Multiplex ELISA was performed using a multi-analyte ELISA plate (Biorad, Hercules, CA) that includes IL-1α, IL‑1β, IL-6, IL-10, IL-15, IL-18, MCP, and TNFα, which was then analyzed using a Bio-plex Magpix (Biorad, Hercules, CA).

Dual-energy X-Ray Absorptiometry (DEXA)

Analysis of bone mineral density, body fat % and lean mass were performed using a Lunar PIXImus II (GE Healthcare, United Kingdom). Animals were anesthetized using a ketamine/xylazine cocktail and then analyzed with a single scan after 4 months of treatment and every 4 months thereafter.

Physical performance assessments

A single investigator, blinded to the group designations of the mice, performed all animal assessments. Additionally, all experiments were performed during lighted hours and at the same time of day at each assessment time point. Protocols for each assessment were as follows:

Grip strength

Maximal grip strength data was generated using a Columbus instruments grip force meter (Columbus, OH) as the average of the best 3 of 5 trials for each mouse. For each trial the mouse was held firmly near the base of the tail and placed with all four paws upon a metal grid attached to a force meter. The mouse was then pulled such that the body of the mouse was parallel to the ground until the mouse lost grip. Mice were given 10 seconds of rest between trials.

Treadmill endurance

The mice were given a single assessment at each time point on a Columbus instruments treadmill set with no inclination (flat - 0°). In the trial, the belt slowly accelerates from 5 to 25 m/min over 60 minutes and the mouse is timed until exhaustion, defined as 10 visits to a shock pad (54V, 0.72mA), 20 total shocks, or having remained on the treadmill belt for 60 minutes. Prior to the initial assessment, mice were given 3 similar trials (separated two weeks apart) to acclimate the animals to the device.

Uphill sprint assessment

To assess uphill sprint endurance the treadmill was inclined at 25° and mice were given a warm up period of 1 minute at 5 m/min. This was then followed by intervals that started at 10 m/min for 20 seconds and increased by 1 m/min increments after each 20 second active recovery period at 5 m/min. The mouse continued until exhaustion, defined as visiting the shock pad 5 times or receiving 10 total shocks.

Grip wire endurance

Mice were placed on a 5 mm thick neoprene wire and timed until fall. Dividers were placed on either side of the wire to prevent mice from leaving the apparatus. Mice were given three attempts to attain a minimum 15 seconds per trial, and the better score of two trials was used for each mouse.

Stride length

To measure stride length, the front paws of the mouse were coated in dye (Bradford reagent) and the mouse was then placed on one end of an apparatus (75 cm x 30 cm wide) that was lined with paper. A mouse shelter from the home cage was inserted on the other end of the apparatus to motivate the mouse to walk. Stride length was measured as the average of 5-6 steps per mouse as the distance between the centers of the paws. The stride length was further normalized by femur length as determined by measuring images of the femur generated by DEXA.

Open field activity

To assess spontaneous activity the mouse was placed into a 90 cm x 90 cm apparatus that was divided equally into 4 quadrants. An investigator was positioned approximately 3 feet away and manually counted crossings into new quadrants and rearings (standing on hind legs) over a 5‑minute period.

Muscle histology

NADH histological analysis of quadriceps muscle (rectus femoris) was assessed as described previously [72]. Briefly, 10 µm frozen muscle sections were submerged in a solution containing 1 mg/ml NADH (Sigma, St. Louis, MO), 1 mg/ml Nitro Blue Tetrazolium (VWR #TCD0844), and 0.2 M Tris‑HCl buffer at pH 7.4 for 45 minutes at 37° C. Sections were then immersed in a series of acetone baths, rinsed in distilled water, and dehydrated by immersing in ethanol and xylene, before finally mounting on a cover slip with Cytoseal (Fisher #23-244257). An investigator, blinded to the identity of the mice, identified and tallied fiber types and also measured cross sectional area (CSA) of the fibers using Motic software (Motic, Hong Kong).

Mitochondrial biomass

Total DNA was isolated from soleus and anterior tibialis muscle using a Qiagen Tissue Quick mini-prep kit (Qiagen, Germantown, MD). Primers were designed to amplify mitochondrial DNA (forward: 5’-CCGCAAGGGAAAGATGAAAGA-3’, reverse: 5’-TCGTTTGGTTTCGGGGTTTC-3’) and nuclear DNA (hexokinase gene, forward: 5’-CCCTGTCATGTCCCTTTGTT-3’, reverse: 5’‑GCCACCAGCTCAGTTAAAGG-3’) and then amplified using quantitative PCR (LightCycler 2.0, Roche).

Mitochondrial activity

Mitochondrial activity was assessed as described previously [72]. Briefly, soleus muscle was homogenized in Chappel-Perry isolation medium (100 mM KCl, 50 mM Tris‑HCl, 5 mM MgCl2, 1 mM ATP, 1 mM EGTA, pH 7.5) and centrifuged 10 minutes at 900 x g. The supernatant was collected and centrifuged for 10 minutes at 9,000 x g at 4 °C, and the mitochondria rich pellet was then washed once and resuspended in ice‑cold SHE buffer (250 mM sucrose, 10 mM Hepes, 1 mM EGTA, pH 7.2). Protein concentration was determined using a Bradford assay. Complex IV activity was determined as the rate of absorbance (412 nm) decline induced by adding 1 μg of mitochondria into a cuvette containing 1 mL of ice‑cold reaction buffer (10 mM KH2PO4, 250 mM sucrose, 1 mg/ml BSA, 10 µM reduced cytochrome C (reduced using sodium hydrosulfite), 2.5 mM lauryl maltoside, pH 6.5). Citrate activity was determined as the rate of absorbance (412 nm) increase by adding 1 μg of mitochondria extract into a cuvette containing 500 μL of ice‑cold reaction buffer (100 mM Tris, 2 mM DTNB (5,5-dithio-bis-2-nitrobenzoic acid), 4 mM oxaloacetic acid, 1 mM acetyl Co-A, pH 8.0).

Analysis of muscle myofibrillar, sarcoplasmic, and mitochondrial fractions

Our methodology to assess relative percentages of protein fractions was adapted from previous studies [73-75]. Approximately 30-40 mg of gastrocnemius muscle was homogenized using a pestle tissue homogenizer in 3.0 mL of ice-cold analysis buffer (20 mM tris-HCl, 250 mM sucrose, 100 mM KCl, 5 mM EDTA, pH 6.8). The homogenate was centrifuged at 1,000 x g for 15 minutes. To isolate the myofibrillar fraction, the pellet was washed twice with 5 mL of ice cold wash buffer (20 mM tris‑HCl, 175 mM KCl, 5 mM EDTA, 0.5% (v/v) Triton-X100 pH 6.8) and centrifuged at 1,000 x g for 10 minutes at 4°C. The pellet was suspended in 1.0 mL of ice‑cold analysis buffer for subsequent protein determination. To isolate the mitochondrial and sarcoplasmic fractions, the initial supernatant was centrifuged at 9,000 x g for 20 minutes at 4°C. The supernatant was then collected as the sarcoplasmic fraction for subsequent protein determination. The pellet was washed in 5 mL of ice cold SHE buffer (250 mM sucrose, 10 mM Hepes, 1 mM EGTA, pH 7.2) and centrifuged at 9,000 x g for 20 minutes at 4°C. The pellet was then suspended in 100 μL of ice‑cold SHE buffer. Determination of protein concentration was performed using a Bradford assay and used to determine the protein content of each fraction normalized to mg of wet tissue weight.

Western blot analysis

Gastrocnemius muscle was homogenized in ice‑cold extraction buffer (20 mM HEPES, 1 mM EDTA, 5 mM EGTA, 1.5 mM NaVO4, 10 mM MgCl2, 50 mM glycerolphosphate, 2 mM DTT, 10 mM NaF, 1% Triton X-100, 100 µM PMSF, pH 7.4) and protein concentration was determined using Bradford assay. Proteins were resolved by protein electrophoresis by applying 30 μg of each sample onto a 10% Biorad gel (Ready Gel Tris-HCl -101-1394). Assessment of atrogin-1 expression was performed using overnight incubation of sample with antibodies against atrogin-1 (ThermoFisher Scientific, Cat# PA5-19056, Waltham, MA) and tubulin (ThermoFisher Scientific, Cat #62204) for normalization. Quantitation was performed using ImageJ software (NIH).

MicroRNA analysis

Total RNA was harvested from tibialis anterior muscle using Qiagen miRNAeasy Purification Kit (Qiagen, Germantown, MD) according to the manufacturer’s instructions. RNA libraries for sequencing were established using the NEBNext Multiplex Small RNA Library preparation kit, and the MiRNA libraries were sequenced on the Illumina NextSeq 500 generating 76-cycle single reads. Demultiplexing was performed with Illumina’s bcl2fastq version 2.17.1.14. General sequence quality was evaluated with FastQC, and reads were trimmed of adapters using trim galore v0.4.4. Subsequently, reads were aligned to the Ensembl GRCm38 genome build using bowtie2 v2.2.8 with the very-sensitive-local parameter set [76]. Aligned reads were quantified using featureCounts [77] against miRBase v21 miRNA database, and the resulting counts were tested in R using the Bioconductor package DESeq2 [78]. MiRNA with counts of less than 10 reads in 66% of the samples, post normalization, were removed from the analysis. Statistical analysis was performed using DESeq2, which includes a Benjamini‑Hochberg correction for false positives [79].

Statistics

Statistical analysis was performed using XLStat statistical software (Addinsoft, New York, NY). A Student’s t-test was used for all comparisons of standard (STD) supplementation versus insufficient (LOW) supplementation. All data were screened for outliers using a Grubbs outlier test with alpha equal to 0.05. The cut-off for significant comparisons was p < 0.05. All data are presented as mean ± standard deviation.

Supplementary Materials

Supplementary File

Abbreviations

CSA: Cross sectional area; iPTH: intact Parathyroid hormone; IL-6, -10: Interleukin-6, -10; LOW: Low supplementation; miRNA: micro ribonucleic acid; STD: Standard supplementation.

Author Contributions

KLS, MP, MML, JB, MRG, MH, and ZS performed experiments. KLS, NN, and BRT contributed to experimental design. KLS, JB, and BRT analyzed data. KLS and BRT prepared and wrote manuscript.

Acknowledgements

The authors thank Guy Howard of the University of Miami, Jessy Alexander at the University at Buffalo, and Brandon Marzullo at the New York State Center of Excellence in Bioinformatics at the University at Buffalo for their suggestions and insights during the study and preparation of this manuscript. The authors also wish to thank the University at Buffalo and the research service of the VA Western New York Healthcare System.

Conflicts of Interest

The authors of this manuscript declare no conflicts of interest.

Funding

This work was funded by Veteran Affairs Rehabilitation Research and Development Grant RX001066 and the Indian Trail Foundation.

References

  • 1. Cruz-Jentoft AJ, Landi F, Schneider SM, Zúñiga C, Arai H, Boirie Y, Chen LK, Fielding RA, Martin FC, Michel JP, Sieber C, Stout JR, Studenski SA, et al. Prevalence of and interventions for sarcopenia in ageing adults: a systematic review. Report of the International Sarcopenia Initiative (EWGSOP and IWGS). Age Ageing. 2014; 43:748–59. https://doi.org/10.1093/ageing/afu115 [PubMed]
  • 2. Volpato S, Bianchi L, Cherubini A, Landi F, Maggio M, Savino E, Bandinelli S, Ceda GP, Guralnik JM, Zuliani G, Ferrucci L. Prevalence and clinical correlates of sarcopenia in community-dwelling older people: application of the EWGSOP definition and diagnostic algorithm. J Gerontol A Biol Sci Med Sci. 2014; 69:438–46. https://doi.org/10.1093/gerona/glt149 [PubMed]
  • 3. Scott D, Blizzard L, Fell J, Ding C, Winzenberg T, Jones G. A prospective study of the associations between 25-hydroxy-vitamin D, sarcopenia progression and physical activity in older adults. Clin Endocrinol (Oxf). 2010; 73:581–87. https://doi.org/10.1111/j.1365-2265.2010.03858.x [PubMed]
  • 4. Cherniack EP, Levis S, Troen BR. Hypovitaminosis D: a widespread epidemic. Geriatrics. 2008; 63:24–30. [PubMed]
  • 5. Wei MY, Giovannucci EL. Vitamin D and multiple health outcomes in the Harvard cohorts. Mol Nutr Food Res. 2010; 54:1114–26. [PubMed]
  • 6. Maloney S, Buckler J, Gil-Rivas V, Goolkasian P. A Retrospective Study of Serum Vitamin D Levels in Veterans. Fed Pract. 2011; 28:31–37.
  • 7. Houston DK, Tooze JA, Davis CC, Chaves PH, Hirsch CH, Robbins JA, Arnold AM, Newman AB, Kritchevsky SB. Serum 25-hydroxyvitamin D and physical function in older adults: the Cardiovascular Health Study All Stars. J Am Geriatr Soc. 2011; 59:1793–801. https://doi.org/10.1111/j.1532-5415.2011.03601.x [PubMed]
  • 8. Houston DK, Tooze JA, Hausman DB, Johnson MA, Nicklas BJ, Miller ME, Neiberg RH, Marsh AP, Newman AB, Blair SN, Kritchevsky SB. Change in 25-hydroxyvitamin D and physical performance in older adults. J Gerontol A Biol Sci Med Sci. 2011; 66:430–36. https://doi.org/10.1093/gerona/glq235 [PubMed]
  • 9. Michael YL, Smit E, Seguin R, Curb JD, Phillips LS, Manson JE. Serum 25-hydroxyvitamin D and physical performance in postmenopausal women. J Womens Health (Larchmt). 2011; 20:1603–08. https://doi.org/10.1089/jwh.2010.2606 [PubMed]
  • 10. Van Schoor NM, Heymans MW, Lips P. Vitamin D status in relation to physical performance, falls and fractures in the Longitudinal Aging Study Amsterdam: A reanalysis of previous findings using standardized serum 25-hydroxyvitamin D values. J Steroid Biochem Mol Biol. 2018; 177:255–60. https://doi.org/10.1016/j.jsbmb.2017.09.012 [PubMed]
  • 11. Marantes I, Achenbach SJ, Atkinson EJ, Khosla S, Melton LJ3rd, Amin S. Is vitamin D a determinant of muscle mass and strength? J Bone Miner Res. 2011; 26:2860–71. https://doi.org/10.1002/jbmr.510 [PubMed]
  • 12. Matheï C, Van Pottelbergh G, Vaes B, Adriaensen W, Gruson D, Degryse JM. No relation between vitamin D status and physical performance in the oldest old: results from the Belfrail study. Age Ageing. 2013; 42:186–90. https://doi.org/10.1093/ageing/afs186 [PubMed]
  • 13. McDermott MM, Liu K, Ferrucci L, Tian L, Guralnik J, Kopp P, Tao H, Van Horn L, Liao Y, Green D, Kibbe M, Criqui MH. Vitamin D status and functional performance in peripheral artery disease. Vasc Med. 2012; 17:294–302. https://doi.org/10.1177/1358863X12448457 [PubMed]
  • 14. Burne TH, McGrath JJ, Eyles DW, Mackay-Sim A. Behavioural characterization of vitamin D receptor knockout mice. Behav Brain Res. 2005; 157:299–308. https://doi.org/10.1016/j.bbr.2004.07.008 [PubMed]
  • 15. Keisala T, Minasyan A, Lou YR, Zou J, Kalueff AV, Pyykkö I, Tuohimaa P. Premature aging in vitamin D receptor mutant mice. J Steroid Biochem Mol Biol. 2009; 115:91–97. https://doi.org/10.1016/j.jsbmb.2009.03.007 [PubMed]
  • 16. Minasyan A, Keisala T, Zou J, Zhang Y, Toppila E, Syvälä H, Lou YR, Kalueff AV, Pyykkö I, Tuohimaa P. Vestibular dysfunction in vitamin D receptor mutant mice. J Steroid Biochem Mol Biol. 2009; 114:161–66. https://doi.org/10.1016/j.jsbmb.2009.01.020 [PubMed]
  • 17. Buitrago CG, Ronda AC, de Boland AR, Boland R. MAP kinases p38 and JNK are activated by the steroid hormone 1alpha,25(OH)2-vitamin D3 in the C2C12 muscle cell line. J Cell Biochem. 2006; 97:698–708. https://doi.org/10.1002/jcb.20639 [PubMed]
  • 18. Salles J, Chanet A, Giraudet C, Patrac V, Pierre P, Jourdan M, Luiking YC, Verlaan S, Migné C, Boirie Y, Walrand S. 1,25(OH)2-vitamin D3 enhances the stimulating effect of leucine and insulin on protein synthesis rate through Akt/PKB and mTOR mediated pathways in murine C2C12 skeletal myotubes. Mol Nutr Food Res. 2013; 57:2137–46. https://doi.org/10.1002/mnfr.201300074 [PubMed]
  • 19. Domingues-Faria C, Chanet A, Salles J, Berry A, Giraudet C, Patrac V, Denis P, Bouton K, Goncalves-Mendes N, Vasson MP, Boirie Y, Walrand S. Vitamin D deficiency down-regulates Notch pathway contributing to skeletal muscle atrophy in old wistar rats. Nutr Metab (Lond). 2014; 11:47. https://doi.org/10.1186/1743-7075-11-47 [PubMed]
  • 20. Girgis CM, Mokbel N, Cha KM, Houweling PJ, Abboud M, Fraser DR, Mason RS, Clifton-Bligh RJ, Gunton JE. 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:3227–37. https://doi.org/10.1210/en.2014-1016 [PubMed]
  • 21. Tanaka M, Kishimoto KN, Okuno H, Saito H, Itoi E. Vitamin D receptor gene silencing effects on differentiation of myogenic cell lines. Muscle Nerve. 2014; 49:700–08. https://doi.org/10.1002/mus.23950 [PubMed]
  • 22. Girgis CM, Cha KM, Houweling PJ, Rao R, Mokbel N, Lin M, Clifton-Bligh RJ, Gunton JE. Vitamin D Receptor Ablation and Vitamin D Deficiency Result in Reduced Grip Strength, Altered Muscle Fibers, and Increased Myostatin in Mice. Calcif Tissue Int. 2015; 97:602–10. https://doi.org/10.1007/s00223-015-0054-x [PubMed]
  • 23. Groves NJ, Bradford D, Sullivan RK, Conn KA, Aljelaify RF, McGrath JJ, Burne TH. Behavioural Effects of Adult Vitamin D Deficiency in BALB/c Mice Are not Associated with Proliferation or Survival of Neurons in the Adult Hippocampus. PLoS One. 2016; 11:e0152328. https://doi.org/10.1371/journal.pone.0152328 [PubMed]
  • 24. Kasukawa Y, Miyakoshi N, Maekawa S, Nozaka K, Noguchi H, Shimada Y. Effects of alfacalcidol on muscle strength, muscle fatigue, and bone mineral density in normal and ovariectomized rats. Biomed Res. 2010; 31:273–79. https://doi.org/10.2220/biomedres.31.273 [PubMed]
  • 25. Camperi A, Pin F, Costamagna D, Penna F, Menduina ML, Aversa Z, Zimmers T, Verzaro R, Fittipaldi R, Caretti G, Baccino FM, Muscaritoli M, Costelli P. Vitamin D and VDR in cancer cachexia and muscle regeneration. Oncotarget. 2017; 8:21778–93. https://doi.org/10.18632/oncotarget.15583 [PubMed]
  • 26. Endo I, Inoue D, Mitsui T, Umaki Y, Akaike M, Yoshizawa T, Kato S, Matsumoto T. Deletion of vitamin D receptor gene in mice results in abnormal skeletal muscle development with deregulated expression of myoregulatory transcription factors. Endocrinology. 2003; 144:5138–44. https://doi.org/10.1210/en.2003-0502 [PubMed]
  • 27. Bhat M, Kalam R, Qadri SS, Madabushi S, Ismail A. Vitamin D deficiency-induced muscle wasting occurs through the ubiquitin proteasome pathway and is partially corrected by calcium in male rats. Endocrinology. 2013; 154:4018–29. https://doi.org/10.1210/en.2013-1369 [PubMed]
  • 28. Yin K, Agrawal DK. Vitamin D and inflammatory diseases. J Inflamm Res. 2014; 7:69–87. https://doi.org/10.2147/JIR.S63898 [PubMed]
  • 29. Calton EK, Keane KN, Newsholme P, Soares MJ. The Impact of Vitamin D Levels on Inflammatory Status: A Systematic Review of Immune Cell Studies. PLoS One. 2015; 10:e0141770. https://doi.org/10.1371/journal.pone.0141770 [PubMed]
  • 30. Choi H, Myung K. Calcitriol enhances fat synthesis factors and calpain activity in co-cultured cells. Cell Biol Int. 2014; 38:910–17. https://doi.org/10.1002/cbin.10281 [PubMed]
  • 31. Brown DM, Goljanek-Whysall K. microRNAs: modulators of the underlying pathophysiology of sarcopenia? Ageing Res Rev. 2015; 24:263–73. https://doi.org/10.1016/j.arr.2015.08.007 [PubMed]
  • 32. Jung HJ, Lee KP, Milholland B, Shin YJ, Kang JS, Kwon KS, Suh Y. Comprehensive miRNA Profiling of Skeletal Muscle and Serum in Induced and Normal Mouse Muscle Atrophy During Aging. J Gerontol A Biol Sci Med Sci. 2017; 72:1483–91. https://doi.org/10.1093/gerona/glx025 [PubMed]
  • 33. Belenchia AM, Johnson SA, Kieschnick AC, Rosenfeld CS, Peterson CA. Time Course of Vitamin D Depletion and Repletion in Reproductive-age Female C57BL/6 Mice. Comp Med. 2017; 67:483–90. [PubMed]
  • 34. Testerink J, Jaspers RT, Rittweger J, de Haan A, Degens H. Effects of alfacalcidol on circulating cytokines and growth factors in rat skeletal muscle. J Physiol Sci. 2011; 61:525–35. https://doi.org/10.1007/s12576-011-0174-7 [PubMed]
  • 35. Yakabe M, Ogawa S, Ota H, Iijima K, Eto M, Ouchi Y, Akishita M. Inhibition of interleukin-6 decreases atrogene expression and ameliorates tail suspension-induced skeletal muscle atrophy. PLoS One. 2018; 13:e0191318. https://doi.org/10.1371/journal.pone.0191318 [PubMed]
  • 36. Dodds RM, Roberts HC, Cooper C, Sayer AA. The Epidemiology of Sarcopenia. J Clin Densitom. 2015; 18:461–66. https://doi.org/10.1016/j.jocd.2015.04.012 [PubMed]
  • 37. Stoever K, Heber A, Eichberg S, Brixius K. Sarcopenia and Predictors of Skeletal Muscle Mass in Elderly Men With and Without Obesity. Gerontol Geriatr Med. 2017; 3:2333721417713637. https://doi.org/10.1177/2333721417713637 [PubMed]
  • 38. Grimaldi AS, Parker BA, Capizzi JA, Clarkson PM, Pescatello LS, White MC, Thompson PD. 25(OH) vitamin D is associated with greater muscle strength in healthy men and women. Med Sci Sports Exerc. 2013; 45:157–62. https://doi.org/10.1249/MSS.0b013e31826c9a78 [PubMed]
  • 39. Haslam A, Johnson MA, Hausman DB, Cress ME, Houston DK, Davey A, Poon LW, and Georgia Centenarian Study. Vitamin D status is associated with grip strength in centenarians. J Nutr Gerontol Geriatr. 2014; 33:35–46. https://doi.org/10.1080/21551197.2013.867825 [PubMed]
  • 40. Muir SW, Montero-Odasso M. Effect of vitamin D supplementation on muscle strength, gait and balance in older adults: a systematic review and meta-analysis. J Am Geriatr Soc. 2011; 59:2291–300. https://doi.org/10.1111/j.1532-5415.2011.03733.x [PubMed]
  • 41. Tomlinson PB, Joseph C, Angioi M. Effects of vitamin D supplementation on upper and lower body muscle strength levels in healthy individuals. A systematic review with meta-analysis. J Sci Med Sport. 2015; 18:575–80. https://doi.org/10.1016/j.jsams.2014.07.022 [PubMed]
  • 42. Glendenning P, Zhu K, Inderjeeth C, Howat P, Lewis JR, Prince RL. Effects of three-monthly oral 150,000 IU cholecalciferol supplementation on falls, mobility, and muscle strength in older postmenopausal women: a randomized controlled trial. J Bone Miner Res. 2012; 27:170–76. https://doi.org/10.1002/jbmr.524 [PubMed]
  • 43. Stockton KA, Mengersen K, Paratz JD, Kandiah D, Bennell KL. Effect of vitamin D supplementation on muscle strength: a systematic review and meta-analysis. Osteoporos Int. 2011; 22:859–71. [PubMed]
  • 44. Duque G, Daly RM, Sanders K, Kiel DP. Vitamin D, bones and muscle: myth versus reality. Australas J Ageing. 2017 (Suppl 1); 36:8–13. https://doi.org/10.1111/ajag.12408 [PubMed]
  • 45. Carter RJ, Morton J, Dunnett SB. Motor coordination and balance in rodents.. Curr Protoc Neurosci. 2001; Chapter 8:Unit 8.12. https://doi.org/10.1002/0471142301.ns0812s15 [PubMed]
  • 46. Sakai S, Suzuki M, Tashiro Y, Tanaka K, Takeda S, Aizawa K, Hirata M, Yogo K, Endo K. Vitamin D receptor signaling enhances locomotive ability in mice. J Bone Miner Res. 2015; 30:128–36. https://doi.org/10.1002/jbmr.2317 [PubMed]
  • 47. Seldeen KL, Pang M, Rodríguez-Gonzalez M, Hernandez M, Sheridan Z, Yu P, Troen BR. A mouse model of vitamin D insufficiency: is there a relationship between 25(OH) vitamin D levels and obesity? Nutr Metab (Lond). 2017; 14:26. https://doi.org/10.1186/s12986-017-0174-6 [PubMed]
  • 48. Chen W, Jiao X, Zhang J, Wang L, Yu X. Vitamin D deficiency and high serum IL-6 concentration as risk factors for tubal factor infertility in Chinese women. Nutrition. 2018; 49:24–31. https://doi.org/10.1016/j.nut.2017.11.016 [PubMed]
  • 49. Gubatan J, Mitsuhashi S, Longhi MS, Zenlea T, Rosenberg L, Robson S, Moss AC. Higher serum vitamin D levels are associated with protective serum cytokine profiles in patients with ulcerative colitis. Cytokine. 2018; 103:38–45. https://doi.org/10.1016/j.cyto.2017.12.023 [PubMed]
  • 50. Lira FS, Rosa JC, Cunha CA, Ribeiro EB, do Nascimento CO, Oyama LM, Mota JF. Supplementing alpha-tocopherol (vitamin E) and vitamin D3 in high fat diet decrease IL-6 production in murine epididymal adipose tissue and 3T3-L1 adipocytes following LPS stimulation. Lipids Health Dis. 2011; 10:37. https://doi.org/10.1186/1476-511X-10-37 [PubMed]
  • 51. Dey BK, Gagan J, Yan Z, Dutta A. miR-26a is required for skeletal muscle differentiation and regeneration in mice. Genes Dev. 2012; 26:2180–91. https://doi.org/10.1101/gad.198085.112 [PubMed]
  • 52. Lee SW, Yang J, Kim SY, Jeong HK, Lee J, Kim WJ, Lee EJ, Kim HS. MicroRNA-26a induced by hypoxia targets HDAC6 in myogenic differentiation of embryonic stem cells. Nucleic Acids Res. 2015; 43:2057–73. https://doi.org/10.1093/nar/gkv088 [PubMed]
  • 53. Davidsen PK, Gallagher IJ, Hartman JW, Tarnopolsky MA, Dela F, Helge JW, Timmons JA, Phillips SM. High responders to resistance exercise training demonstrate differential regulation of skeletal muscle microRNA expression. J Appl Physiol (1985). 2011; 110:309–17. https://doi.org/10.1152/japplphysiol.00901.2010 [PubMed]
  • 54. Margolis LM, Lessard SJ, Ezzyat Y, Fielding RA, Rivas DA. Circulating MicroRNA Are Predictive of Aging and Acute Adaptive Response to Resistance Exercise in Men. J Gerontol A Biol Sci Med Sci. 2017; 72:1319–26. https://doi.org/10.1093/gerona/glw243 [PubMed]
  • 55. Jorde R, Svartberg J, Joakimsen RM, Coucheron DH. Plasma profile of microRNA after supplementation with high doses of vitamin D3 for 12 months. BMC Res Notes. 2012; 5:245. https://doi.org/10.1186/1756-0500-5-245 [PubMed]
  • 56. Cheng X, Du J, Shen L, Tan Z, Jiang D, Jiang A, Li Q, Tang G, Jiang Y, Wang J, Li X, Zhang S, Zhu L. MiR-204-5p regulates C2C12 myoblast differentiation by targeting MEF2C and ERRγ. Biomed Pharmacother. 2018; 101:528–35. https://doi.org/10.1016/j.biopha.2018.02.096 [PubMed]
  • 57. Mi L, Li Y, Zhang Q, Zhao C, Peng Y, Yang G, Zheng X. MicroRNA-139-5p regulates C2C12 cell myogenesis through blocking Wnt/β-catenin signaling pathway. Biochem Cell Biol. 2015; 93:8–15. https://doi.org/10.1139/bcb-2014-0079 [PubMed]
  • 58. Nielsen S, Åkerström T, Rinnov A, Yfanti C, Scheele C, Pedersen BK, Laye MJ. The miRNA plasma signature in response to acute aerobic exercise and endurance training. PLoS One. 2014; 9:e87308. https://doi.org/10.1371/journal.pone.0087308 [PubMed]
  • 59. Kuang W, Tan J, Duan Y, Duan J, Wang W, Jin F, Jin Z, Yuan X, Liu Y. Cyclic stretch induced miR-146a upregulation delays C2C12 myogenic differentiation through inhibition of Numb. Biochem Biophys Res Commun. 2009; 378:259–63. https://doi.org/10.1016/j.bbrc.2008.11.041 [PubMed]
  • 60. Sun Y, Li Y, Wang H, Li H, Liu S, Chen J, Ying H. miR-146a-5p acts as a negative regulator of TGF-β signaling in skeletal muscle after acute contusion. Acta Biochim Biophys Sin (Shanghai). 2017; 49:628–34. https://doi.org/10.1093/abbs/gmx052 [PubMed]
  • 61. Guess MG, Barthel KK, Harrison BC, Leinwand LA. miR-30 family microRNAs regulate myogenic differentiation and provide negative feedback on the microRNA pathway. PLoS One. 2015; 10:e0118229. https://doi.org/10.1371/journal.pone.0118229 [PubMed]
  • 62. Sun Z, Evans J, Bhagwate A, Middha S, Bockol M, Yan H, Kocher JP. CAP-miRSeq: a comprehensive analysis pipeline for microRNA sequencing data. BMC Genomics. 2014; 15:423. https://doi.org/10.1186/1471-2164-15-423 [PubMed]
  • 63. Mallya SM, Corrado KR, Saria EA, Yuan FF, Tran HQ, Saucier K, Atti E, Tetradis S, Arnold A. Modeling vitamin D insufficiency and moderate deficiency in adult mice via dietary cholecalciferol restriction. Endocr Res. 2016; 41:290–99. https://doi.org/10.3109/07435800.2016.1141937 [PubMed]
  • 64. Buitrago C, Pardo VG, Boland R. Role of VDR in 1α,25-dihydroxyvitamin D3-dependent non-genomic activation of MAPKs, Src and Akt in skeletal muscle cells. J Steroid Biochem Mol Biol. 2013; 136:125–30. https://doi.org/10.1016/j.jsbmb.2013.02.013 [PubMed]
  • 65. Garcia LA, King KK, Ferrini MG, Norris KC, Artaza JN. 1,25(OH)2vitamin D3 stimulates myogenic differentiation by inhibiting cell proliferation and modulating the expression of promyogenic growth factors and myostatin in C2C12 skeletal muscle cells. Endocrinology. 2011; 152:2976–86. https://doi.org/10.1210/en.2011-0159 [PubMed]
  • 66. Ryan ZC, Craig TA, Folmes CD, Wang X, Lanza IR, Schaible NS, Salisbury JL, Nair KS, Terzic A, Sieck GC, Kumar R. 1α,25-Dihydroxyvitamin D3 Regulates Mitochondrial Oxygen Consumption and Dynamics in Human Skeletal Muscle Cells. J Biol Chem. 2016; 291:1514–28. https://doi.org/10.1074/jbc.M115.684399 [PubMed]
  • 67. Hassan-Smith ZK, Jenkinson C, Smith DJ, Hernandez I, Morgan SA, Crabtree NJ, Gittoes NJ, Keevil BG, Stewart PM, Hewison M. 25-hydroxyvitamin D3 and 1,25-dihydroxyvitamin D3 exert distinct effects on human skeletal muscle function and gene expression. PLoS One. 2017; 12:e0170665. https://doi.org/10.1371/journal.pone.0170665 [PubMed]
  • 68. Hirani V, Cumming RG, Naganathan V, Blyth F, Le Couteur DG, Hsu B, Handelsman DJ, Waite LM, Seibel MJ. Longitudinal Associations Between Vitamin D Metabolites and Sarcopenia in Older Australian men: The Concord Health and Aging in Men Project. J Gerontol A Biol Sci Med Sci. 2017; 73:131–38. https://doi.org/10.1093/gerona/glx086 [PubMed]
  • 69. Boonen S, Lysens R, Verbeke G, Joosten E, Dejaeger E, Pelemans W, Flamaing J, Bouillon R. Relationship between age-associated endocrine deficiencies and muscle function in elderly women: a cross-sectional study. Age Ageing. 1998; 27:449–54. https://doi.org/10.1093/ageing/27.4.449 [PubMed]
  • 70. Gielen E, O’Neill TW, Pye SR, Adams JE, Wu FC, Laurent MR, Claessens F, Ward KA, Boonen S, Bouillon R, Vanderschueren D, Verschueren S. Endocrine determinants of incident sarcopenia in middle-aged and elderly European men. J Cachexia Sarcopenia Muscle. 2015; 6:242–52. https://doi.org/10.1002/jcsm.12030 [PubMed]
  • 71. Srikuea R, Zhang X, Park-Sarge OK, Esser KA. VDR and CYP27B1 are expressed in C2C12 cells and regenerating skeletal muscle: potential role in suppression of myoblast proliferation. Am J Physiol Cell Physiol. 2012; 303:C396–405. https://doi.org/10.1152/ajpcell.00014.2012 [PubMed]
  • 72. Seldeen KL, Lasky G, Leiker MM, Pang M, Personius KE, Troen BR. High Intensity Interval Training Improves Physical Performance and Frailty in Aged Mice. J Gerontol A Biol Sci Med Sci. 2018; 73:429–37. https://doi.org/10.1093/gerona/glx120 [PubMed]
  • 73. Donges CE, Burd NA, Duffield R, Smith GC, West DW, Short MJ, Mackenzie R, Plank LD, Shepherd PR, Phillips SM, Edge JA. Concurrent resistance and aerobic exercise stimulates both myofibrillar and mitochondrial protein synthesis in sedentary middle-aged men. J Appl Physiol (1985). 2012; 112:1992–2001. https://doi.org/10.1152/japplphysiol.00166.2012 [PubMed]
  • 74. Haus JM, Carrithers JA, Carroll CC, Tesch PA, Trappe TA. Contractile and connective tissue protein content of human skeletal muscle: effects of 35 and 90 days of simulated microgravity and exercise countermeasures. Am J Physiol Regul Integr Comp Physiol. 2007; 293:R1722–27. https://doi.org/10.1152/ajpregu.00292.2007 [PubMed]
  • 75. Nakanishi T, Jarmakani JM. Developmental changes in myocardial mechanical function and subcellular organelles. Am J Physiol. 1984; 246:H615–25. https://doi.org/10.1152/ajpheart.1984.246.4.H615 [PubMed]
  • 76. Ziemann M, Kaspi A, El-Osta A. Evaluation of microRNA alignment techniques. RNA. 2016; 22:1120–38. https://doi.org/10.1261/rna.055509.115 [PubMed]
  • 77. Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics. 2014; 30:923–30. https://doi.org/10.1093/bioinformatics/btt656 [PubMed]
  • 78. Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 2014; 15:550. https://doi.org/10.1186/s13059-014-0550-8 [PubMed]
  • 79. Benjamini Y, Hochberg Y. Controlling the False Discovery Rate - a Practical and Powerful Approach to Multiple Testing. J R Stat Soc B. 1995; 57:289–300.