Advertisement
Research| Volume 102, ISSUE 5, P4072-4093, May 2019

Download started.

Ok

Phosphorus content of muscle tissue and muscle function in dairy cows fed a phosphorus-deficient diet during the transition period

Open AccessPublished:March 01, 2019DOI:https://doi.org/10.3168/jds.2018-15727

      ABSTRACT

      Phosphorus (P) deficiency and hypophosphatemia are believed to be associated with muscle function disturbances in dairy cows, particularly around parturition. The objective of this study was to determine the effect of dietary P deprivation during late gestation and early lactation on muscle P homeostasis and muscle function in periparturient dairy cows. Thirty-six multiparous dairy cows in late gestation were randomly assigned either to undergo dietary P depletion or to be offered a diet with adequate P content from 4 wk before to 4 wk after parturition. Phosphorus-deficient rations for dry and lactating cows contained 0.15 and 0.20% P on a dry matter basis, respectively. Blood and muscle tissue for biopsy were obtained and electromyographic examinations were conducted on biceps femoris and intercostal muscles in regular intervals throughout the study. Muscle tissue was analyzed for the total P, adenosine triphosphate, adenosine diphosphate, adenosine monophosphate, creatine phosphate, and tissue water content. Dietary P deprivation resulted in a pronounced and sustained decline of the plasma phosphate concentration, reaching a nadir at calving with mean values below 1.5 mg/dL and remaining below 2.0 mg/dL during the first 4 wk of lactation. Hypophosphatemia was not associated with signs of clinically apparent muscle weakness or disturbed muscle function and was not associated with a decline in the content of any of the studied P-containing compounds in muscle tissue. Accordingly, no association between plasma phosphate concentration and muscle tissue P content was found. Electromyographic examination identified subclinical effects on motor unit action potentials that are indicative of disturbed neuromuscular functionality. Increasing occurrence of pathologic spontaneous activity possibly resulting from membrane instability of nerve or muscle cells and suggestive of myopathy was also recorded as P deprivation progressed. These effects were predominantly observed in intercostal and to a lesser degree biceps femoris muscles. Electromyographic parameters affected by P deprivation were found to be associated primarily with the plasma phosphate and to a lesser extent with the amounts of energy storing P-containing compounds contained in muscle tissue. These results indicate that prolonged and pronounced dietary P deprivation in transition dairy cows leads to marked sustained hypophosphatemia without altering the muscle tissue P homeostasis or causing clinically apparent muscle function disturbances.

      Key words

      INTRODUCTION

      Hypophosphatemia is commonly encountered in dairy cows during the periparturient period and in early lactation (
      • Macrae A.I.
      • Whitaker D.A.
      • Burrough E.
      • Dowell A.
      • Kelly J.M.
      Use of metabolic profiles for the assessment of dietary adequacy in UK dairy herds.
      ,
      • Macrae A.I.
      • Burrough E.
      • Forrest J.
      Assessment of nutrition in dairy herds: Use of metabolic profiles.
      ). Although subnormal plasma inorganic phosphorus (P) concentrations ([Pi]) have been reported in at least 50% of clinically healthy dairy cows around parturition, hypophosphatemia is also diagnosed in the majority of recumbent fresh cows, frequently in association with subnormal plasma calcium concentrations ([Ca];
      • Grünberg W.
      Treatment of phosphorus balance disorders.
      ). Furthermore, hypophosphatemia was found to be more pronounced in recumbent dairy cattle not responding to intravenous calcium administration compared with recumbent cows making a full recovery after calcium supplementation (
      • Ménard L.
      • Thompson A.
      Milk fever and alert downer cows: Does hypophosphatemia affect the treatment response?.
      ). However, the proposed association between hypophosphatemia and recumbency in cattle is solely based on empirical evidence and has been discussed controversially for decades. Attempts to experimentally induce hypophosphatemic recumbency in cattle have been unsuccessful thus far, and reports about treatment efficacy of oral or parenteral P supplementation in cattle with presumed hypophosphatemic recumbency are inconsistent at best (
      • Rodehutscord M.
      • Pauen A.
      • Windhausen P.
      • Brintrup R.
      • Pfeffer E.
      Effects of drastic changes in P intake on P concentrations in blood and rumen fluid of lactating ruminants.
      ;
      • Grünberg W.
      Treatment of phosphorus balance disorders.
      ;
      • Puggaard L.
      • Lund P.
      • Liesegang A.
      • Sehested J.
      Long term effect of reduced dietary phosphorus on feed intake and milk yield in dry and lactating dairy cows.
      ;
      • Grünberg W.
      • Scherpenisse P.
      • Dobbelaar P.
      • Idink M.J.
      • Wijnberg I.D.
      The effect of transient, moderate dietary phosphorus deprivation on phosphorus metabolism, muscle content of different phosphorus-containing compounds, and muscle function in dairy cows.
      ). A potential mechanism through which P depletion or hypophosphatemia may result in disturbed muscle function has not been unraveled, but depletion of energy stores in muscle tissue in states of P deficiency has been proposed in combination with membrane instability (
      • Grünberg W.
      Treatment of phosphorus balance disorders.
      ;
      • Grünberg W.
      • Scherpenisse P.
      • Dobbelaar P.
      • Idink M.J.
      • Wijnberg I.D.
      The effect of transient, moderate dietary phosphorus deprivation on phosphorus metabolism, muscle content of different phosphorus-containing compounds, and muscle function in dairy cows.
      ). Muscle fibers depend on the availability of creatine phosphate (CPh) and adenosine phosphates as a source of energy, and it is thought that hypophosphatemia is indicative of P deficiency and thus is accompanied by depletion of intracellular P stores.
      Hypophosphatemia in early-lactating dairy cows is receiving increased attention because of recent incentives in many parts of the world aiming for the reduction of the P content in dairy cow rations, with the objective of containing environmental pollution with P of fecal origin (
      • Goff J.P.
      Invited review: Mineral absorption mechanisms, mineral interactions that affect acid-base and antioxidant status, and diet considerations to improve mineral status.
      ). Producers and veterinarians raised concerns that this practice may exacerbate disturbances of the P balance occurring in periparturient cows and thereby jeopardize health and productivity of the fresh dairy cow.
      The objective of this study was to determine the effects of pronounced dietary P deprivation during the transition period on muscle tissue P homeostasis and muscle function at a clinical and subclinical level. Based on the results of earlier studies conducted on nonperiparturient cattle, we hypothesized that dietary P deprivation would result in pronounced hypophosphatemia that is not associated with clinically apparent muscle function disturbances or alterations of the muscle tissue P balance.

      MATERIALS AND METHODS

      Ethics Statement

      The national and institutional guidelines for the care and use of experimental animals were followed and all experimental procedures were approved by the Utrecht University Institutional Animal Care and Use Committee (permit no. AVD108002016616). This study is part of a multi-institutional project evaluating the effect of dietary P deprivation in dairy cows. To further understand the molecular mechanisms through which P depletion may affect various organ systems of the dairy cow at different stages of the lactation cycle, other laboratories are currently conducting studies on liver and immune cell function from these transition cows. The present study focuses on the effect of dietary P deprivation on muscle function in the periparturient period.

      Animals and Housing

      A total of 36 healthy, multiparous, pregnant Holstein-Friesian or Holstein-Friesian crossbreed dairy cows in late gestation were used for this study. Cows were purchased from commercial dairy farms in the Netherlands at least 1 mo before enrollment. For logistical reasons, the study was conducted in 2 consecutive replicates consisting of 18 cows each. Cows of each replicate were required to have been artificially inseminated within a range of 5 d to narrow the calving period of each replicate as much as possible. Only cows found to be healthy based on physical, hematological, and blood biochemical examination were included in the study.
      Cows were housed in individual tiestalls with rubber bedding covered with sawdust in a temperature-controlled facility. Study animals were dried off at least 8 wk before the expected calving date (i.e., were dry for at least 2 wk before enrollment).

      Study Design and Experimental Rations

      The study was conceived as a randomized block design. Study animals of each replicate were blocked by lactation number and 305-d milk yield of the previous lactation and were then, within block, randomly assigned to a treatment with either deficient (LP) or adequate (AP) dietary P supply by drawing cards. The study period covered a minimum of 11 wk starting 6 wk before the expected calving date (Figure 1). The study was initiated with an acclimation period of 2 wk for both treatments that was followed by a period of deficient dietary P supply for LP cows and adequate dietary P supply for AP cows, extending from 4 wk before the expected week of calving to 4 wk after the expected week of calving (Figure 1). The LP cows of the first replicate underwent an additional 2-wk P repletion period immediately following the P-depletion period, whereas the AP cows of the first replicate were kept on the diet with adequate P supply. During acclimation both treatments received the same dry cow ration with adequate P content. During the P-repletion period of the first replicate, LP cows were switched to the AP diet and thus both groups were fed the same lactating cow ration. The objective of the repletion period was to determine whether effects that may have resulted from dietary P depletion in LP animals were reversed by dietary P supplementation over 2 wk.
      Figure thumbnail gr1
      Figure 1Experimental timeline for adequate-phosphorus (AP) and low-phosphorus (LP) treatments. Each shaded and open rectangle represents 1 wk of time. Open arrows represent the sampling times T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion), at which electromyograms were conducted and muscle biopsies were obtained. Small arrows represent the blood sampling time points. Only cows of replicate 1 underwent the repletion period following phosphorus depletion. During acclimation, cows in both treatments received the same dry cow ration with adequate P supply. During repletion and adequate P supply postpartum, both groups received the same lactating cow ration with adequate P supply; repletion and adequate P supply postpartum II applied only to the first replicate of the study.
      Feed was offered as a balanced TMR for dry or lactating cows depending on the stage of the study. The same base ration formulated to meet current recommendations for dairy cows, with exception of the dietary P content, was prepared for both groups (
      • NRC
      Nutrient Requirements of Dairy Cattle.
      ; Table 1). This P-deficient base ration was fed to LP cows, whereas AP cows were offered this ration supplemented with NaH2PO4 to obtain a dietary P content meeting or slightly exceeding current recommendations for dairy cows (
      • NRC
      Nutrient Requirements of Dairy Cattle.
      ). The P content of dry and lactating cow LP rations was 0.15 and 0.20% P on a DM basis, respectively. Dry cow and lactating cow rations of the AP treatment had a dietary P content of 0.28 and 0.44% P in DM, respectively (Table 1). The dietary P content of the LP rations fed during the depletion period was the lowest of an otherwise balanced ration we were able to achieve using standard ration ingredients. The dietary P content of the lactating cow AP ration slightly exceeded current recommendations for dietary P supply but met the P content commonly fed to high-yielding dairy cows in practice (
      • NRC
      Nutrient Requirements of Dairy Cattle.
      ). Orts were weighed back, and fresh feed was offered twice daily between 0600 and 0700 h and between 1800 and 1900 h. Access to feed was restricted to 12.5 kg of DM/d for dry cows and was ad libitum for lactating cows, and cows had ad libitum access to water. Lactating cows were milked twice daily between 0600 and 0700 h and between 1800 and 1900 h.
      Table 1Ingredients (% of DM) and composition (g/kg of DM unless noted) of experimental dry cow (antepartum) and lactating cow (postpartum) rations of low-phosphorus (LP) and adequate-phosphorus (AP) treatments
      Cows assigned to the LP treatment received the antepartum ration of AP during acclimation and the postpartum ration of AP during the repletion period.
      ItemAntepartumPostpartum
      APLPAPLP
      Ingredient
       Corn silage29.229.439.139.6
       Grass seed straw
      Grass seed straw from smooth meadow grass (Poa pratensis) known for its quality and palatability, similar to common grass hay.
      25.826.012.412.6
       Beet pulp35.235.536.136.6
       Soybean meal008.78.8
       Palatinose
      Palatinose was used as a substitute for molasses to stimulate rumen fermentation but without increasing the dietary K content.
      4.84.800
       Urea0.680.680.340.35
       Coated urea2.032.051.001.01
       Mineral mix1.56
      Mineral mix contained Ca (0 g/kg), Mg (180 g/kg), Na (100 g/kg), Cu (1.655 g/kg), Co (0.022 g/kg), Mn (0.9 g/kg), Zn (0.65 g/kg), I (0.01 g/kg), Se (0.015 g/kg), vitamin A (340,000 IU/kg), β-carotene (0.6 g/kg), vitamin D3 (140,000 IU/kg), and vitamin E (8,000 IU/kg).
      1.57
      Mineral mix contained Ca (0 g/kg), Mg (180 g/kg), Na (100 g/kg), Cu (1.655 g/kg), Co (0.022 g/kg), Mn (0.9 g/kg), Zn (0.65 g/kg), I (0.01 g/kg), Se (0.015 g/kg), vitamin A (340,000 IU/kg), β-carotene (0.6 g/kg), vitamin D3 (140,000 IU/kg), and vitamin E (8,000 IU/kg).
      1.10
      Mineral mix contained Ca (170 g/kg), Mg (38 g/kg), Na (120 g/kg), Cu (0.73 g/kg), Co (0.009 g/kg), Mn (1.2 g/kg), Zn (1.67 g/kg), I (0.042 g/kg), Se (0.017 g/kg), vitamin A (280,000 IU/kg), β-carotene (2.27 g/kg), vitamin D3 (90,000 IU/kg), and vitamin E (900 IU/kg).
      1.11
      Mineral mix contained Ca (170 g/kg), Mg (38 g/kg), Na (120 g/kg), Cu (0.73 g/kg), Co (0.009 g/kg), Mn (1.2 g/kg), Zn (1.67 g/kg), I (0.042 g/kg), Se (0.017 g/kg), vitamin A (280,000 IU/kg), β-carotene (2.27 g/kg), vitamin D3 (90,000 IU/kg), and vitamin E (900 IU/kg).
       NaH2PO40.6701.230
      Chemical analysis
       DM (% as fed)60.059.358.358.2
       NEL (MJ/kg of DM)5.905.886.486.54
       DVE
      Intestinal digestible protein (Tamminga et al., 1994).
      61.060.779.480.4
       OEB
      Degraded protein balance (Tamminga et al., 1994).
      (g/d)
      18.117.911.812.0
       Starch111110153154
       Ca4.354.385.85.7
       K13.513.414.114.5
       P2.831.524.412.03
       Mg3.473.511.831.8
       S1.241.251.361.41
       Na4.31.823.621.5
       DCAD (mEq/kg of DM)109.9105.0294.6290.8
      1 Cows assigned to the LP treatment received the antepartum ration of AP during acclimation and the postpartum ration of AP during the repletion period.
      2 Grass seed straw from smooth meadow grass (Poa pratensis) known for its quality and palatability, similar to common grass hay.
      3 Palatinose was used as a substitute for molasses to stimulate rumen fermentation but without increasing the dietary K content.
      4 Mineral mix contained Ca (0 g/kg), Mg (180 g/kg), Na (100 g/kg), Cu (1.655 g/kg), Co (0.022 g/kg), Mn (0.9 g/kg), Zn (0.65 g/kg), I (0.01 g/kg), Se (0.015 g/kg), vitamin A (340,000 IU/kg), β-carotene (0.6 g/kg), vitamin D3 (140,000 IU/kg), and vitamin E (8,000 IU/kg).
      5 Mineral mix contained Ca (170 g/kg), Mg (38 g/kg), Na (120 g/kg), Cu (0.73 g/kg), Co (0.009 g/kg), Mn (1.2 g/kg), Zn (1.67 g/kg), I (0.042 g/kg), Se (0.017 g/kg), vitamin A (280,000 IU/kg), β-carotene (2.27 g/kg), vitamin D3 (90,000 IU/kg), and vitamin E (900 IU/kg).
      6 Intestinal digestible protein (
      • Tamminga S.
      • Vanstraalen W.M.
      • Subnel A.P.J.
      • Meijer R.G.M.
      • Steg A.
      • Wever C.J.G.
      • Blok M.C.
      The Dutch protein evaluation system—The DVE/OEB-system.
      ).
      7 Degraded protein balance (
      • Tamminga S.
      • Vanstraalen W.M.
      • Subnel A.P.J.
      • Meijer R.G.M.
      • Steg A.
      • Wever C.J.G.
      • Blok M.C.
      The Dutch protein evaluation system—The DVE/OEB-system.
      ).

      Sample Collection and Experimental Procedures

      Animal Health, Feed Intake, Milk Production, and BW

      Attitude and behavior of the study animals were monitored daily, and a complete physical examination was conducted once a week. Cows were weighed on an electronic scale at the end of the acclimation period (T1), after 2 wk of dietary P depletion (T2), during the first week of lactation (T3), at the end of the depletion period (T4), and, for cows of replicate 1, at the end of the repletion phase (T5; Figure 1). Milk yield and DMI were recorded daily, and the average daily DMI and milk yield of the 7 d before each of the 5 sampling times were calculated for every cow.

      Blood and Feed Sample Collection

      Blood samples were collected by venipuncture of a jugular vein every Friday between 0700 and 0800 h. Blood was collected into blood tubes containing lithium–heparin as anticoagulant (Vacuette, Greiner Bio-One, Kremsmünster, Austria) and was centrifuged within 30 min of collection at 1,000 × g for 15 min at 6°C. Harvested plasma was stored at −21°C until analyzed as described below.
      Feed samples were collected once a week and were analyzed for their DM and P content to confirm the adequate P content of the experimental rations. Dry matter was determined by drying at 100°C to constant weight at atmospheric pressure; dietary P content was measured by inductively coupled plasma MS. The results of the P analysis of the experimental diets confirmed that the dietary P content of the dry cow LP ration was on average 0.15% P in DM (range: 0.14–0.18%). The mean dietary P content of the LP ration for lactating cows was 0.20% P in DM (range: 0.19–0.20%). The mean P content of AP diets was 0.28% (range: 0.26–0.32%) and 0.44% (range: 0.42–0.47%) for the dry cow and lactating cow rations, respectively.

      Electromyographic Examination

      Cows on study underwent repeated electromyographic examination to detect pathologic spontaneous activity as an indicator of muscular dysfunction and to conduct motor unit action potential (MUAP) and interference pattern analyses (
      • Finsterer J.
      • Mamoli B.
      • Fuglsang-Frederiksen A.
      Peak-ratio interference pattern analysis in the detection of neuromuscular disorders.
      ;
      • Rubin D.I.
      Needle electromyography: Basic concepts and patterns of abnormalities.
      ). Electromyographic examinations were conducted at the time points T1, T2, T3, T4, and, for cows of replicate 1, T5. Electromyographies were always conducted and analyzed by the same investigator as described in earlier studies (
      • Wijnberg I.D.
      • Back W.
      • De Jong M.
      • Zuidhof M.C.
      • Van Den Belt A.J.M.
      • Van Der Kolk J.H.
      The role of electromyography in clinical diagnosis of neuromuscular locomotor problems in the horse.
      ;
      • Grünberg W.
      • Scherpenisse P.
      • Dobbelaar P.
      • Idink M.J.
      • Wijnberg I.D.
      The effect of transient, moderate dietary phosphorus deprivation on phosphorus metabolism, muscle content of different phosphorus-containing compounds, and muscle function in dairy cows.
      ). Briefly, the electrical activity of intercostal (IC) and biceps femoris (BF) muscles was determined to evaluate pathologic spontaneous activity and to assess the properties of isolated MUAP and interference patterns (
      • Liguori R.
      • Fuglsang-Frederiksen A.
      • Nix W.
      • Fawcett P.R.W.
      • Andersen K.
      Electromyography in myopathy.
      ). Interference patterns of the IC muscle were recorded during breathing at rest and those of the BF muscle were recorded during weight bearing with an elevated contralateral limb (
      • Nirkko A.C.
      • Rosler K.M.
      • Hess C.W.
      Sensitivity and specificity of needle electromyography: A prospective study comparing automated interference pattern analysis with single motor unit potential analysis.
      ). Electromyographies were recorded using portable equipment (Viking Quest EMG system, Nicolet Biomedical Inc., Madison, WI;
      • Kimura J.
      Electronic systems and data analysis.
      ;
      • Daube J.R.
      • Rubin D.I.
      Needle electromyography.
      ) and were made using disposable concentric 26-gauge, 50-mm and 21-gauge, 100-mm electromyography needle electrodes (MEDcat, Emmen, the Netherlands). Cows were placed in stocks in a calm environment and remained without sedation for the examination. Motor unit action potentials of the IC muscle were recorded during normal breathing with the cow standing calmly in the stocks. No extra manipulation was required to induce MUAP in this muscle. Motor unit action potentials of the BF muscle were generated by inducing voluntary muscle contractions as described earlier by lifting the contralateral limb (
      • Grünberg W.
      • Scherpenisse P.
      • Dobbelaar P.
      • Idink M.J.
      • Wijnberg I.D.
      The effect of transient, moderate dietary phosphorus deprivation on phosphorus metabolism, muscle content of different phosphorus-containing compounds, and muscle function in dairy cows.
      ).

      Muscle Tissue Collection

      Muscle tissue for biopsy of the BF muscle was obtained from all animals on the days of electromyographic examination (T1–T5) immediately following electromyography and from the limb opposed to the one used for the foregoing electromyography. Tissue samples were obtained under local anesthesia (procaine hydrochloride 4%, V.M.D. n.v., Arendonk, Belgium) and following administration of a nonsteroidal anti-inflammatory drug (meloxicam 0.5 mg/kg s.c.; Metacam, Boehringer Ingelheim) using a suction-modified Bergstrom biopsy needle (6-mm diameter; Walter Veterinär-Instrumente, Baruth, Germany) as described earlier in cattle (
      • Grünberg W.
      • Scherpenisse P.
      • Dobbelaar P.
      • Idink M.J.
      • Wijnberg I.D.
      The effect of transient, moderate dietary phosphorus deprivation on phosphorus metabolism, muscle content of different phosphorus-containing compounds, and muscle function in dairy cows.
      ). Two specimens, each approximately 0.2 to 0.3 g, were obtained at each sampling time. Collected tissue was submerged in liquid nitrogen within 10 s of collection and stored in liquid nitrogen until processed for further analysis as described below.

      Analyses

      Plasma Biochemical Analysis

      Plasma was analyzed for [Pi] (ammonium molybdate method); total [Ca] (Arsenazo III method); the concentrations of creatinine ([Crea]; enzymatic method), nonesterified fatty acids ([NEFA]; acetyl-CoA-synthetase–acetyl-CoA-oxydase method), BHB ([BHB]; UV method), and total bilirubin (dichloraniline method); and the enzyme activity of aspartate aminotransferase (AST; UV method without pyridoxal phosphate) and creatine kinase (CK; N-acetylcysteine-creatine phosphate method). Plasma biochemical analyses were performed by a commercial veterinary diagnostic laboratory and conducted on an automated analyzer (ABX Pentra 400; Horiba, Europe GmbH, Langenhagen, Germany).

      Electromyographic Examination

      Insertional activity, pathologic spontaneous activity, MUAP, satellite potentials, and interference pattern analysis were recorded and studied as described in detail for an earlier study conducted in cattle (
      • Grünberg W.
      • Scherpenisse P.
      • Dobbelaar P.
      • Idink M.J.
      • Wijnberg I.D.
      The effect of transient, moderate dietary phosphorus deprivation on phosphorus metabolism, muscle content of different phosphorus-containing compounds, and muscle function in dairy cows.
      ). Briefly, the occurrence of pathologic spontaneous activity was determined outside the endplate region at the same site from where MUAP were obtained. Four types of pathologic spontaneous activity were differentiated: fibrillation potentials, positive sharp waves, complex repetitive discharges, and neuromyotonia. Pathologic spontaneous activity was counted as such when occurring in at least 2 locations within the same muscle; the incidence of pathologic spontaneous activity was categorized semiquantitatively into 4 groups from 0 (no pathologic spontaneous activity) to +++ (very frequent pathologic spontaneous activity;
      • Georgesco M.
      • Salerno A.
      Spontaneous activity in electromyography: Mechanisms and practical interest.
      ). For quantitative MUAP analysis, the amplitude, duration, number of phases, and number of turns, area under the action potential curve (referred to as “area”), and size index were calculated from at least 20 recorded MUAP per muscle (
      • Wijnberg I.D.
      • Franssen H.
      • van der Kolk J.H.
      • Breukink H.J.
      Quantitative analysis of motor unit action potentials in the subclavian muscle of healthy horses.
      ,
      • Wijnberg I.D.
      • Back W.
      • De Jong M.
      • Zuidhof M.C.
      • Van Den Belt A.J.M.
      • Van Der Kolk J.H.
      The role of electromyography in clinical diagnosis of neuromuscular locomotor problems in the horse.
      ). For interference pattern analysis of the IC muscle, a minimum of 20 contractions were evaluated by determining the maximum voluntary activity expressed as the ratio turns/second and amplitude/turns during spontaneous breathing patterns.

      Muscle Tissue Biochemical Analysis

      The content in muscle tissue wet weight (ww) of adenosine triphosphate (ATP), adenosine diphosphate (ADP), adenosine monophosphate (AMP), creatine, and CPh were determined as described earlier (
      • Teerlink T.
      • Hennekes M.
      • Bussemaker J.
      • Groeneveld J.
      Simultaneous determination of creatine compounds and adenine nucleotides in myocardial tissue by high-performance liquid chromatography.
      ;
      • Jiang Y.
      • Sun C.
      • Ding X.
      • Yuan D.
      • Chen K.
      • Gao B.
      • Chen Y.
      • Sun A.
      Simultaneous determination of adenine nucleotides, creatine phosphate and creatine in rat liver by high performance liquid chromatography-electrospray ionization-tandem mass spectrometry.
      ;
      • Grünberg W.
      • Scherpenisse P.
      • Dobbelaar P.
      • Idink M.J.
      • Wijnberg I.D.
      The effect of transient, moderate dietary phosphorus deprivation on phosphorus metabolism, muscle content of different phosphorus-containing compounds, and muscle function in dairy cows.
      ). Briefly, frozen muscle tissue was homogenized in 2 mL of ice-cold 0.66 M perchloric acid. After centrifugation at 3,500 × g and 4°C for 10 min, the supernatant was transferred into a 5-mL glass tube. The residue was extracted and centrifuged again as described above. Both supernatants were mixed and the volume was topped up with perchloric acid to 5 mL. Samples were stored at −20°C until analyzed as described below. Before analysis, 0.5 mL of extract was neutralized by mixing with 0.3 mL of 2 M KOH, 0.1 mL of 1 M KHCO3, and 0.1 mL of water. The sample was refrigerated for 20 min, centrifuged as described above, and finally diluted 5 times with water. The analysis of AMPww, creatineww, and CPhww was conducted on an HPLC-MS system consisting of 2 HPLC pumps, an autosampler with tray cooling set at 4°C (PE200 series, PESciex, Norwalk, CT), and a triple quadrupole MS detector (ABSciex, Framingham, MA) with a turbo ion spray source. The HPLC column was a hypercarb (5-µm porous graphitic carbon, 150 × 2.1 mm; Thermo Scientific, Breda, the Netherlands). Eluent 1 for pump A was 8 mM ammonium acetate, adjusted to pH 9.8 with 25% ammonia solution, and eluent 2 for pump B was 8 mM ammonium acetate, pH 9.8/acetonitrile (40/60; vol/vol). The flow rate was 200 µL/min, and a gradient run was made. The first 3 min after injection was discarded by a divert valve. The MS was operated in multiple reaction monitoring mode. Analysis of ADPww and ATPww was also done by gradient HPLC with UV detection at 259 nm (Shimadzu, ‘s-Hertogenbosch, the Netherlands). Eluent A was 50 mM KH2PO4 and 10 mM tetrabuyl ammonium hydrogen sulfate (pH 6.5), and eluent B was methanol. The column was a Luna C18 (3 µL, 150 × 2 mm; Phenomenex, Torrance, CA). The HPLC was controlled by Lab Solutions version 5.63 software (Shimadzu, Duisburg, Germany).
      The P and potassium (K) concentrations in the solution containing the tissue homogenate were determined by inductively coupled optical emission spectrometry. Potassium, like P, is a predominantly intracellular electrolyte and was analyzed to assess whether changes in muscle tissue P content over time are specific to P or whether they also occur with other intracellular electrolytes. The total amounts of P, K, ATP, ADP, AMP, and CPh per gram of wet weight muscle tissue were then calculated from the concentration of each compound in the homogenate and the mass of the wet weight of each tissue specimen. The DM content of each muscle tissue specimen was determined by drying to constant weight at 85°C at atmospheric pressure.

      Statistical Analysis

      Data are expressed as arithmetic mean ± standard deviation or as median and interquartile range. The significance level was set at P < 0.05. Values were log transformed when necessary to achieve normal distribution. Repeated-measures ANOVA with an autoregressive(1) covariance matrix with animal ID as subject was used to determine effects of replicate number, block, treatment, time (sampling time or, in the case of blood biochemical parameters, week relative to calving), and the interaction between treatment and time using PROC MIXED (SAS 9.4, SAS Institute Inc., Cary, NC). Repeated-measures ANOVA on ranks was conducted with the model as described above with severity of pathologic spontaneous activity as the dependent variable. For this purpose, categories of 0 to +++ were transformed to numerical values 0 to 3. If a significant replicate number effect was determined for one parameter, it was analyzed for each replicate separately. An autoregressive(1) covariance structure was chosen based on the lowest Akaike information criterion. Bonferroni-adjusted P-values were used to assess differences between treatments at specific sampling times whenever the F-test was significant.
      Forward stepwise regression analyses were conducted for each of the following parameters of the muscle tissue analysis as dependent variables for each of the 5 sampling times T1 to T5: Pww, Kww, ATPww, ADPww, AMPww, CPhww, and DM. The parameters included as independent variables with P = 0.2 as entry and P = 0.05 as exit, with exception of treatment and animal ID that were forced into the model, were chosen to identify possible factors affecting the intracellular P balance. These were lactation number, DMI, BW, age, and previous 305-d milk yield as cow-specific parameters; [Pi] as the parameter reflecting the extracellular P balance; [Ca] as the parameter associated with disturbed muscle function in periparturient cows (periparturient hypocalcemia); [Crea] as the parameter associated with muscle tissue metabolism; CK and AST as indicators for muscle cell injury; [NEFA] and [BHB] as indicators for disturbed energy balance; total bilirubin as an indicator for liver function; DM as an indicator for changes of muscle tissue water content; and Kww as an indicator for changes of intracellular electrolyte balance in muscle tissue that is not specific for P. The final model was checked for variance inflation by determining correlations between independent variables remaining in the model and screening tolerance and variance inflation factors of each variable in the final model. In case of significant associations between independent variables in the final model, the variable with the lowest coefficient of determination was removed and the analysis was rerun. Pearson correlation analyses were conducted to compare the associations of Pww, Kww, [Pi], and DM with other parameters of the intracellular P homeostasis (ATPww, ADPww, AMPww, CPhww) at the different sampling times.
      Forward multiple stepwise regression analyses were also conducted for each of the following parameters of the MUAP analysis of IC and BF muscles stratified by sampling time: MUAP duration, amplitude, number of phases, number of turns, area under the action potential curve, and size index. For the interference pattern analysis, dependent variables included in the analysis were turns per second (M), the amplitude per turn (T), and the turns per second:amplitude per turn ratio (M/T). Independent variables included in these analyses were treatment, animal ID, and, for IC muscles, the replicate number (all forced into the model) as well as lactation number, DMI, BW, previous 305-d milk yield, age, [Pi], [Ca], [Crea], [NEFA], [BHBA], total bilirubin, CK, AST, Pww, Kww, ATPww, ADPww, AMPww, creatine ww, and CPhww. Possible variance inflation was determined as described above. Pearson correlation analyses were again conducted to compare the associations of MUAP duration, amplitude, number of phases, number of turns, area under the action potential curve, and size index and IPA turns per second, amplitude per turn, and M/T with parameters of the intra- and extracellular P homeostasis such as ATPww, ADPww, AMPww, CPhww, and [Pi] as well as with [Ca] at the different sampling times.
      The required sample size for this study was estimated on the basis of results obtained from an earlier study investigating the effects of dietary P deprivation in mid-lactating dairy cows (
      • Grünberg W.
      • Scherpenisse P.
      • Dobbelaar P.
      • Idink M.J.
      • Wijnberg I.D.
      The effect of transient, moderate dietary phosphorus deprivation on phosphorus metabolism, muscle content of different phosphorus-containing compounds, and muscle function in dairy cows.
      ). All analyses were conducted with SAS software (SAS 9.4, SAS Institute Inc.).

      RESULTS

      Study Animals

      The LP and AP cows were 4.7 ± 1.2 and 4.5 ± 1.2 yr old, respectively, and in lactation 3 ± 1. The 305-d milk yield of the previous lactation was 8,935 ± 2,534 and 9,127 ± 1,731 kg for LP and AP cows, respectively. Because cows were blocked by lactation number and 305-d milk yield, differences of these parameters between groups were not significant. The mean BW for LP and AP cows at the time of enrollment was 692 ± 79 and 662 ± 69 kg, respectively, and did not differ between treatments. Cow-related parameters mentioned above were not different between first and second replicates.
      All 36 cows delivered healthy calves either spontaneously or with mild to moderate assistance. The cows calved within a range of 13 d and 24 calved within the predicted week of calving and thus were on the experimental dry cow ration and the experimental lactating cow ration for at least 4 wk each. Four AP cows calved between 1 and 3 d before the predicted week of calving, shortening the prepartum period on the experimental ration accordingly. Five LP and 3 AP cows calved between 1 and 3 d following the end of the predicted calving week.
      Clinical periparturient hypocalcemia, confirmed by blood biochemical analysis, occurred in 4 AP cows (2 in each replicate). In all cases cows were successfully treated with oral and parenteral administration of calcium salts. Four cows, all assigned to the LP treatment, were excluded from the study 2 to 4 wk after parturition, thus between 1 and 2 wk after T3. Three of these cows developed postparturient hemoglobinuria during the second week of lactation, which was characterized by severe intravascular hemolysis and ensuing profound anemia and hemoglobinuria. One cow that delivered twins had retained fetal membranes and later developed metritis and clinical ketosis and was excluded due to ongoing deterioration. Results from these cows obtained up to the first week of lactation (T3) were included in the data analysis. Two cows (1 from each group) with abomasal displacement in their second week of lactation underwent surgical correction of the condition immediately after diagnosis, recovered uneventfully, and completed the study. Other clinically apparent health events that occurred during the study period include 4 animals with febrile metritis (including the animal mentioned above; 3 LP and 1 AP) as well 3 cows with clinical ketosis (defined as depressed demeanor and marked ketonuria without other clinically apparent causes of disease; 2 LP and 1 AP). Treatment of metritis consisted of 3 consecutive days of parenteral administration of amoxicillin (10 mg/kg; Norobrittin, Norbrook Laboratories Ltd. Newry, Northern Ireland) and single or repeated administration of a nonsteroidal anti-inflammatory drug (meloxicam 0.5 mg/kg; Metacam, Boehringer Ingelheim, Ingelheim am Rhein, Germany). Treatment for ketosis consisted of a single or repeated (up to 3 times) intravenous administration of 500 mL of 30% dextrose solution (Glucose 30%, B. Braun, Melsungen, Germany) and oral administration of 250 mL of feed-grade propylene glycol twice daily until ketosis resolved. Dextrose treatments were never administered within 48 h before T3 or T4. Blood samples obtained within 48 h after dextrose administration were retrospectively excluded from the data analysis.

      Blood Biochemical Analysis

      The concentration–time curves for plasma [Pi] and [Ca] stratified by treatment are presented in Figure 2. Significant treatment, time, and treatment × time interaction effects (all P < 0.0001) for plasma [Pi] were identified. Mean concentrations determined in LP were below AP values throughout the entire depletion period (Figure 2), and mean plasma [Pi] of LP determined during the depletion period were below values measured during the acclimation period and below the reference for [Pi] in cattle (4–8 mg/dL; Figure 2;
      • Goff J.P.
      Treatment of calcium, phosphorus, and magnesium balance disorders.
      ). Values did not differ between treatments during the acclimation and repletion periods of the study (Figure 2). Significant treatment and time effects (P = 0.0154 and P < 0.0001, respectively) were identified for plasma [Ca], with higher concentrations in LP cows than AP cows (Figure 2).
      Figure thumbnail gr2
      Figure 2Mean ± SD plasma phosphorus concentrations ([Pi]; upper panel) and plasma calcium concentrations ([Ca]; lower panel) over time in cows assigned to the low-phosphorus (LP; closed circles/solid line) and adequate-phosphorus (AP; open squares/dashed line) diets. Vertical dashed lines mark the start and end of the P-depletion phase, and the vertical dotted line marks the time of parturition. The shaded area represents the reference range for plasma [Pi] and [Ca] in cattle. Time points labeled with an × differ significantly between groups (P < 0.05); time points labeled with an asterisk in the upper panel differ significantly from values measured during the acclimation period (i.e., wk −7 and wk −6) within each group (P < 0.05, Bonferroni corrected). Arrows mark the sampling times T1 to T5 (T1 = end of acclimation period, T2 = 2 wk after onset of P depletion, T3 = within 1 wk of parturition, T4 = end of P depletion, and T5 = end of P repletion).

      Muscle Tissue Biochemical Analysis

      The muscle tissue Pww, ATPww, ADPww, AMPww, and CPhww at the different sampling times are presented in Figure 3. A replicate number effect was not determined for any of these parameters so that results from both replicates of cows were pooled per treatment for the statistical analysis.
      Figure thumbnail gr3
      Figure 3Mean ± SD muscle tissue total phosphorus (P; upper left panel), total potassium (K; upper right panel), adenosine triphosphate (ATP; middle left panel), creatine phosphate (middle right panel), adenosine monophosphate (AMP; lower left panel), and adenosine diphosphate (ADP; lower right panel) in µmol/g of wet weight (ww) for cows on low-phosphorus (LP; closed circles/solid line) and adequate-phosphorus (AP; open squares/dashed line) diets for each of the 5 sampling times (T1 = end of acclimation period, T2 = 2 wk after onset of P depletion, T3 = within 1 wk of parturition, T4 = end of P depletion, and T5 = end of P repletion). Time points with different uppercase (LP) or lowercase (AP) letters (A–C and a,b, respectively) differ between sampling times within treatment (P < 0.05, Bonferroni corrected). Time points marked with an asterisk differ significantly between groups (P < 0.05).
      The Pww revealed neither a treatment nor time effect, whereas the treatment × time interaction effect was significant (P < 0.0001). The mean Pww of LP obtained at T5 was lower than values obtained at T1, T3, and T4 (Figure 3). At T5 the mean Pww of LP was below values determined in AP cows. The Pww showed a strong positive association with Kww for both treatments and at all sampling times. Accordingly, a treatment × time interaction effect (P = 0.0002) but no treatment or time effects were also identified for Kww. As for Pww values, the Kww measured at T5 were lower than at all other sampling times in group LP and were lower than AP values at T4 and T5 (Figure 3).
      The CPhww showed a significant time effect (P < 0.0001) but neither a treatment nor a treatment × time interaction effect. The mean CPhww of both groups combined declined significantly from 31.72 ± 7.16 µmol/g at T1 to 22.66 ± 5.04 µmol/g at T4 and increased again to 30.24 ± 6.94 µmol/g at T5 (Figure 3).
      For ATPww, a time effect (P = 0.0301) and a treatment × time interaction effect (P = 0.0089) were determined, but no treatment effect was determined. The ATPww content of both groups combined was higher at T3 (5.04 ± 0.71 µmol/g) and T4 (5.05 ± 0.72 µmol/g) compared with T1 (4.67 ± 0.58 µmol/g). In LP cows, the highest ATPww values were determined at T4 and were different from values measured at T1 and T5. Differences between treatments were not observed at any time.
      Significant time effects (P < 0.0001) and treatment × time interaction effects (P = 0.0097) were identified for ADPww. The most remarkable finding was a significant decline of ADPww in LP at T5 relative to all other sampling times (Figure 3). Muscle tissue AMPww showed a significant time effect (P = 0.0009) but neither a treatment nor a treatment × time interaction effect (Figure 3).
      Muscle tissue DM was subject to a significant time effect (P < 0.0001) but neither a treatment nor a treatment × time interaction effect. Values determined between T1 and T3 ranged between 23.4 ± 1.5% and 24.0 ± 1.3% and declined by more than 5% to 22.5 ± 1.1 and 22.1 ± 1.0% at the sampling times T4 and T5, respectively (Figure 4).
      Figure thumbnail gr4
      Figure 4Mean ± SD of muscle tissue DM content (% of wet weight) for cows on low-phosphorus (LP; closed circles/solid line) and adequate-phosphorus (AP; open squares/dashed line) diets for each of the 5 sampling times (T1 = end of acclimation period, T2 = 2 wk after onset of P depletion, T3 = within 1 wk of parturition, T4 = end of P depletion, and T5 = end of P repletion). Time points with different uppercase (LP) or lowercase (AP) letters (A–C and a,b, respectively) differ between sampling times within treatment (P < 0.05, Bonferroni corrected). Time points marked with an asterisk differ significantly between treatments (P < 0.05).
      The results of the multiple stepwise regression analysis conducted with Pww, CPhww, ATPww, ADPww, and AMPww as dependent variables for each sampling time are presented in Table 2. The most prominent finding was again a tight positive association between Pww and Kww at all sampling times. Positive associations with total bilirubin in plasma were identified for ATPww at T2, T3, and T5 as well as for ADPww at T5 (Table 2). Plasma [Pi] was not associated with any of the P-containing compounds at any of the sampling times.
      Table 2Results of multiple stepwise regression analyses with phosphorus (Pww), creatine phosphate (CPhww), adenosine triphosphate (ATPww), adenosine diphosphate (ADPww), and adenosine monophosphate (AMPww) in wet muscle tissue (ww) as well as muscle tissue DM as dependent variables
      Parameters included as independent variables were treatment and cow ID (both forced into the model) as well as age, lactation number, BW, 305-d milk yield of the previous lactation, and DMI; plasma concentrations of inorganic phosphate ([Pi]), calcium ([Ca]), potassium ([K]), total bilirubin ([TBil]), creatinine, nonesterified fatty acids ([NEFA]), and BHB ([BHB]); and the activity in plasma of creatinine kinase ([CK]) and aspartate aminotransferase. For the analysis of DM, the variables Pww, CPhww, ATPww, ADPww, and AMPww were also included as independent variables.
      Results are stratified by sampling times—T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion for the second replicate)—and given as partial coefficient of determination (pr2) and P-value. Associations are listed only when significant (P < 0.05) and (pr2) > = 0.1. The direction of the association is indicated by the sign preceding pr2(+ = positive; − = negative).
      ItemT1T2T3T4T5
      PwwKww: +pr
      Results are stratified by sampling times—T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion for the second replicate)—and given as partial coefficient of determination (pr2) and P-value. Associations are listed only when significant (P < 0.05) and (pr2) > = 0.1. The direction of the association is indicated by the sign preceding pr2(+ = positive; − = negative).
      = 0.968; P < 0.0001
      Kww: +pr
      Results are stratified by sampling times—T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion for the second replicate)—and given as partial coefficient of determination (pr2) and P-value. Associations are listed only when significant (P < 0.05) and (pr2) > = 0.1. The direction of the association is indicated by the sign preceding pr2(+ = positive; − = negative).
      = 0.896; P < 0.0001
      Kww: +pr
      Results are stratified by sampling times—T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion for the second replicate)—and given as partial coefficient of determination (pr2) and P-value. Associations are listed only when significant (P < 0.05) and (pr2) > = 0.1. The direction of the association is indicated by the sign preceding pr2(+ = positive; − = negative).
      = 0.816; P < 0.0001
      Kww: +pr
      Results are stratified by sampling times—T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion for the second replicate)—and given as partial coefficient of determination (pr2) and P-value. Associations are listed only when significant (P < 0.05) and (pr2) > = 0.1. The direction of the association is indicated by the sign preceding pr2(+ = positive; − = negative).
      = 0.264; P < 0.0001
      Kww: +pr
      Results are stratified by sampling times—T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion for the second replicate)—and given as partial coefficient of determination (pr2) and P-value. Associations are listed only when significant (P < 0.05) and (pr2) > = 0.1. The direction of the association is indicated by the sign preceding pr2(+ = positive; − = negative).
      = 0.719; P ≥ 0.0001
      CPhwwBW: −pr
      Results are stratified by sampling times—T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion for the second replicate)—and given as partial coefficient of determination (pr2) and P-value. Associations are listed only when significant (P < 0.05) and (pr2) > = 0.1. The direction of the association is indicated by the sign preceding pr2(+ = positive; − = negative).
      = 0.205; P = 0.008
      Kww: +pr
      Results are stratified by sampling times—T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion for the second replicate)—and given as partial coefficient of determination (pr2) and P-value. Associations are listed only when significant (P < 0.05) and (pr2) > = 0.1. The direction of the association is indicated by the sign preceding pr2(+ = positive; − = negative).
      = 0.304; P = 0.001
      Kww: +pr
      Results are stratified by sampling times—T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion for the second replicate)—and given as partial coefficient of determination (pr2) and P-value. Associations are listed only when significant (P < 0.05) and (pr2) > = 0.1. The direction of the association is indicated by the sign preceding pr2(+ = positive; − = negative).
      = 0.540; P < 0.0001
      NA
      No significant association.
      NA
      ATPww[Ca]: +pr
      Results are stratified by sampling times—T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion for the second replicate)—and given as partial coefficient of determination (pr2) and P-value. Associations are listed only when significant (P < 0.05) and (pr2) > = 0.1. The direction of the association is indicated by the sign preceding pr2(+ = positive; − = negative).
      = 0.134; P = 0.03
      Kww: +pr
      Results are stratified by sampling times—T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion for the second replicate)—and given as partial coefficient of determination (pr2) and P-value. Associations are listed only when significant (P < 0.05) and (pr2) > = 0.1. The direction of the association is indicated by the sign preceding pr2(+ = positive; − = negative).
      = 0.466; P < 0.0001 [TBil]: +pr
      Results are stratified by sampling times—T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion for the second replicate)—and given as partial coefficient of determination (pr2) and P-value. Associations are listed only when significant (P < 0.05) and (pr2) > = 0.1. The direction of the association is indicated by the sign preceding pr2(+ = positive; − = negative).
      = 0.112; P = 0.01
      Kww: +pr
      Results are stratified by sampling times—T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion for the second replicate)—and given as partial coefficient of determination (pr2) and P-value. Associations are listed only when significant (P < 0.05) and (pr2) > = 0.1. The direction of the association is indicated by the sign preceding pr2(+ = positive; − = negative).
      = 0.674; P < 0.0001 [TBil]: +pr
      Results are stratified by sampling times—T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion for the second replicate)—and given as partial coefficient of determination (pr2) and P-value. Associations are listed only when significant (P < 0.05) and (pr2) > = 0.1. The direction of the association is indicated by the sign preceding pr2(+ = positive; − = negative).
      = 0.154; P < 0.0001
      [CK]: −pr
      Results are stratified by sampling times—T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion for the second replicate)—and given as partial coefficient of determination (pr2) and P-value. Associations are listed only when significant (P < 0.05) and (pr2) > = 0.1. The direction of the association is indicated by the sign preceding pr2(+ = positive; − = negative).
      = 0.131; P = 0.02
      [TBil]: +pr
      Results are stratified by sampling times—T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion for the second replicate)—and given as partial coefficient of determination (pr2) and P-value. Associations are listed only when significant (P < 0.05) and (pr2) > = 0.1. The direction of the association is indicated by the sign preceding pr2(+ = positive; − = negative).
      = 0.454; P = 0.001
      ADPwwNA[Ca]: +pr
      Results are stratified by sampling times—T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion for the second replicate)—and given as partial coefficient of determination (pr2) and P-value. Associations are listed only when significant (P < 0.05) and (pr2) > = 0.1. The direction of the association is indicated by the sign preceding pr2(+ = positive; − = negative).
      = 0.1276; P = 0.03 Kww: −pr
      Results are stratified by sampling times—T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion for the second replicate)—and given as partial coefficient of determination (pr2) and P-value. Associations are listed only when significant (P < 0.05) and (pr2) > = 0.1. The direction of the association is indicated by the sign preceding pr2(+ = positive; − = negative).
      = 0.126; P = 0.02
      Kww: −pr
      Results are stratified by sampling times—T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion for the second replicate)—and given as partial coefficient of determination (pr2) and P-value. Associations are listed only when significant (P < 0.05) and (pr2) > = 0.1. The direction of the association is indicated by the sign preceding pr2(+ = positive; − = negative).
      = 0.354; P < 0.0001
      NA[TBil]: +pr
      Results are stratified by sampling times—T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion for the second replicate)—and given as partial coefficient of determination (pr2) and P-value. Associations are listed only when significant (P < 0.05) and (pr2) > = 0.1. The direction of the association is indicated by the sign preceding pr2(+ = positive; − = negative).
      = 0.325; P = 0.006
      AMPww[BHB]: −pr
      Results are stratified by sampling times—T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion for the second replicate)—and given as partial coefficient of determination (pr2) and P-value. Associations are listed only when significant (P < 0.05) and (pr2) > = 0.1. The direction of the association is indicated by the sign preceding pr2(+ = positive; − = negative).
      = 0.276; P = 0.003
      [BHB]: −pr
      Results are stratified by sampling times—T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion for the second replicate)—and given as partial coefficient of determination (pr2) and P-value. Associations are listed only when significant (P < 0.05) and (pr2) > = 0.1. The direction of the association is indicated by the sign preceding pr2(+ = positive; − = negative).
      = 0.273; P = 0.005 Kww: −pr
      Results are stratified by sampling times—T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion for the second replicate)—and given as partial coefficient of determination (pr2) and P-value. Associations are listed only when significant (P < 0.05) and (pr2) > = 0.1. The direction of the association is indicated by the sign preceding pr2(+ = positive; − = negative).
      = 0.203; P = 0.005
      [BHB]: −pr
      Results are stratified by sampling times—T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion for the second replicate)—and given as partial coefficient of determination (pr2) and P-value. Associations are listed only when significant (P < 0.05) and (pr2) > = 0.1. The direction of the association is indicated by the sign preceding pr2(+ = positive; − = negative).
      = 0.144; P = 0.0225
      DM: +pr
      Results are stratified by sampling times—T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion for the second replicate)—and given as partial coefficient of determination (pr2) and P-value. Associations are listed only when significant (P < 0.05) and (pr2) > = 0.1. The direction of the association is indicated by the sign preceding pr2(+ = positive; − = negative).
      = 0.141; P = 0.04
      [CK]: −pr
      Results are stratified by sampling times—T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion for the second replicate)—and given as partial coefficient of determination (pr2) and P-value. Associations are listed only when significant (P < 0.05) and (pr2) > = 0.1. The direction of the association is indicated by the sign preceding pr2(+ = positive; − = negative).
      = 0.353; P = 0.02
      DMNA[TBil]: +pr
      Results are stratified by sampling times—T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion for the second replicate)—and given as partial coefficient of determination (pr2) and P-value. Associations are listed only when significant (P < 0.05) and (pr2) > = 0.1. The direction of the association is indicated by the sign preceding pr2(+ = positive; − = negative).
      = 0.225; P = 0.002
      [NEFA]: +pr
      Results are stratified by sampling times—T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion for the second replicate)—and given as partial coefficient of determination (pr2) and P-value. Associations are listed only when significant (P < 0.05) and (pr2) > = 0.1. The direction of the association is indicated by the sign preceding pr2(+ = positive; − = negative).
      = 0.197; P = 0.01
      ATPww: +pr
      Results are stratified by sampling times—T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion for the second replicate)—and given as partial coefficient of determination (pr2) and P-value. Associations are listed only when significant (P < 0.05) and (pr2) > = 0.1. The direction of the association is indicated by the sign preceding pr2(+ = positive; − = negative).
      = 0.151; P = 0.02
      NA
      1 Parameters included as independent variables were treatment and cow ID (both forced into the model) as well as age, lactation number, BW, 305-d milk yield of the previous lactation, and DMI; plasma concentrations of inorganic phosphate ([Pi]), calcium ([Ca]), potassium ([K]), total bilirubin ([TBil]), creatinine, nonesterified fatty acids ([NEFA]), and BHB ([BHB]); and the activity in plasma of creatinine kinase ([CK]) and aspartate aminotransferase. For the analysis of DM, the variables Pww, CPhww, ATPww, ADPww, and AMPww were also included as independent variables.
      2 Results are stratified by sampling times—T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion for the second replicate)—and given as partial coefficient of determination (pr2) and P-value. Associations are listed only when significant (P < 0.05) and (pr2) > = 0.1. The direction of the association is indicated by the sign preceding pr2(+ = positive; − = negative).
      3 No significant association.
      The results of the Pearson correlation analysis are summarized in Table 3. The contents ATP, ADP, AMP, and CPh showed strong positive associations with Pww, particularly at T2 and T3, whereas no associations were identified during acclimation and repletion. Association of these concentrations with Kww were weaker than with Pww. Positive associations with muscle DM were identified for Pww and to a lesser extent with Kww as well as with ATPww, ADPww, and AMPww at T4 and for ATP at T2 (Table 3). The plasma [Pi] was strongly negatively associated with ATPww, ADPww, and muscle DM at T4 (Table 3).
      Table 3Results of the Pearson correlation analysis (r; P-values in parentheses) of muscle wet weight (ww) adenosine triphosphate (ATPww), adenosine diphosphate (ADPww), adenosine monophosphate (AMPww), and creatine phosphate (CPhww) with muscle phosphorus (Pww), muscle potassium (Kww), plasma inorganic phosphorus concentration [Pi], total bilirubin (TBil), and muscle tissue DM stratified by sampling time
      Sampling times were T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      Empty cells are either cells with the same parameter in row and column or duplicates of cells of the table.
      ItemTimePwwKww[Pi]TBilDMATPwwADPwwAMPww
      PwwT10.99 (P < 0.0001) NA
      NA = no significant association.
      NANA
      T20.97 (P < 0.0001)NA0.41 (P = 0.0205)NA
      T30.93 (P < 0.0001)NANANA
      T40.82 (P < 0.0001)−0.61 (P < 0.0001)NA0.42 (P = 0.02)
      T50.99 (P < 0.0001)NANANA
      ATPwwT1NANANANANA
      T20.68 (P < 0.0001)0.66 (P < 0.0001)NA0.59 (P = 0.0002)0.35 (P = 0.04)
      T30.88 (P < 0.0001)0.84 (P < 0.0001)NA0.56 (P = 0.0004)NA
      T4NANA−0.52 (P = 0.002)NA0.53 (P = 0.002)
      T5NANANA0.81 (P = 0.0001)NA
      ADPwwT10.37 (P = 0.03)0.35 (P = 0.0275)NANANANA
      T20.54 (P = 0.002)0.47 (P = 0.007)NANANA0.56 (P = 0.0004)
      T30.70 (P < 0.0001)0.62 (P < 0.0001)NA0.48 (P = 0.003)NA0.65 (P < 0.0001)
      T40.40 (P = 0.03)0.42 (P = 0.02)−0.59 (P = 0.0007)NA0.47 (P = 0.005)0.84 (P < 0.0001)
      T5NANANA0.76 (P = 0.0007)NA0.89 (P < 0.0001)
      AMPwwT1NANANA0.37 (P = 0.03)NANA0.48 (P = 0.004)
      T20.59 (P = 0.0005)0.50 (P = 0.004)NA0.37 (P = 0.03)NANANA
      T30.37 (P = 0.03)NANANANANA0.42 (P = 0.002)
      T4NANANANA0.35 (P = 0.05)0.53 (P < 0.0001)NA
      T5NANANANANANANA
      CPhwwT1NANANA−0.45 (P = 0.007)NANANANA
      T20.61 (P = 0.0002)0.57 (P = 0.0007)NA0.41 (P = 0.003)NA0.36 (P = 0.03)0.41 (P = 0.01)0.58 (P = 0.0003)
      T30.69 (P < 0.0001)0.74 (P < 0.0001)NANANA0.62 (P < 0.0001)0.34 (P = 0.04)NA
      T4NANANANANA0.45 (P = 0.0091)0.58 (P = 0.0005)0.45 (P = 0.01)
      T5NANANANANANANA
      DMT1NANANANA
      T2NANANA0.54 (P = 0.001)
      T3NANANA0.37 (P = 0.03)
      T40.42 (P = 0.02)0.36 (P = 0.04)−0.36 (P = 0.04)NA
      T5NANANANA
      1 Sampling times were T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      2 Empty cells are either cells with the same parameter in row and column or duplicates of cells of the table.
      3 NA = no significant association.

      Electromyographic Examination

      Pathologic Spontaneous Activity

      Occurrence of pathological spontaneous activity stratified by muscles and sampling times is summarized in Table 4. Minor pathological spontaneous activity of IC muscles was recorded in more than 20% of study animals of either treatment at T1 but occurred with lower frequency in BF muscles (Table 4). The frequency of minor pathologic spontaneous activity increased over time in IC muscles and to a lesser degree in BF muscles for both treatments. A marked increase in the occurrence of moderate (++) and pronounced (+++) pathologic spontaneous activity of IC muscles at T4 and T5 compared with earlier sampling times was observed in LP animals. Treatment and time effects were significant (P = 0.005 and P < 0.0001, respectively) for IC muscles with higher frequency in LP cow than in AP cows and higher frequency at T5 than at T1, T2, and T3 as well as higher frequency at T4 than at T1. Significant treatment and time effects were also identified in BF muscles (P = 0.005 and P < 0.0001, respectively) with higher frequency in LP cows compared with AP cows and higher frequency at T5 than at all other sampling times. Moderate or pronounced pathologic spontaneous activity was observed incidentally only in BF muscles of LP cows (Table 4). The most common types of pathologic spontaneous activity were fibrillation potentials followed by complex repetitive discharges and positive sharp waves.
      Table 4Frequency of pathological spontaneous activity in intercostal and biceps femoris muscles for treatment with low (LP) and adequate (AP) dietary P supply
      Sampling times were T1 (end of acclimation), T2 (after 2 wk of dietary P deprivation antepartum), T3 (within 1 wk postpartum), T4 (at the end of dietary P deprivation), and T5 (after 2 wk of dietary P repletion). At T4, 4 LP cows were removed from the experiment because of health issues; muscle tissue analyses during T4 excluded those cows. At T5, only cows in the first replicate of the experiment were included.
      MuscleGroupGrade
      Occurrence of pathologic spontaneous activity was graded as absent (0), mild (+), moderate (++), or pronounced (+++).
      T1T2T3T4T5
      No.Total%No.Total%No.Total%No.Total%No.Total%
      IntercostalLP0141878818449185021414070
      +318209185071838814571714
      ++118211861186314221714
      +++01800180118611475772
      AP0121866111861111861918503933
      +318176183361833918505956
      ++31817118611860180090
      +++01800180018001801911
      Biceps femorisLP018181001618891418781414783743
      +01802181141822414222729
      ++01800180018001401714
      +++01800180018001401714
      AP0181810017189414187818181006967
      +018011864182201803933
      ++0180018001800180090
      +++0180018001800180090
      1 Sampling times were T1 (end of acclimation), T2 (after 2 wk of dietary P deprivation antepartum), T3 (within 1 wk postpartum), T4 (at the end of dietary P deprivation), and T5 (after 2 wk of dietary P repletion). At T4, 4 LP cows were removed from the experiment because of health issues; muscle tissue analyses during T4 excluded those cows. At T5, only cows in the first replicate of the experiment were included.
      2 Occurrence of pathologic spontaneous activity was graded as absent (0), mild (+), moderate (++), or pronounced (+++).

      MUAP Analysis

      The MUAP parameters of IC muscles were analyzed for each replicate separately as a significant replicate number effect was determined on several parameters. Results are presented in Table 5, Table 6 for replicates 1 and 2, respectively. Time effects were identified on the MUAP amplitude (P = 0.0008 and P < 0.0001), the number of phases (P = 0.0004 and P = 0.0002), the number of turns (P = 0.0007 and P = 0.0492), and the size index (P = 0.0029 and P < 0.0001) for the first and second replicates, respectively (Table 5, Table 6). Treatment effects were only identified in the second replicate on MUAP amplitude (473 ± 153 µV in AP vs. 376 ± 145 µV in LP, P = 0.0002), the number of phases (3.13 ± 0.33 in AP vs. 2.92 ± 0.31 in LP, P = 0.0013), the number of turns (3.34 ± 0.62 in AP vs. 2.97 ± 0.50 in LP, P = 0.0029), and the size index (0.37 ± 0.23 in AP vs. 0.20 ± 0.21 in LP, P = 0.0003). Treatment × time interaction effects were also only apparent in the second replicate for the MUAP amplitude (P = 0.0010), area (P = 0.0002), and size index (P = 0.0084, Table 6). Specifically, the MUAP number of turns, MUAP area, and size index of IC muscles were higher at T1 compared with T2 to T4 in both replicates and tended to increase again between T4 and T5 in replicate 1 (Table 5, Table 6).
      Table 5Results of motor unit action potential analysis of intercostal muscles of cows of the first replicate on adequate (AP) and low (LP) dietary P supply stratified by sampling time
      Sampling times were T1 (end of acclimation), T2 (after 2 wk of dietary P deprivation antepartum), T3 (within 1 wk postpartum), T4 (at the end of dietary P deprivation), and T5 (after 2 wk of dietary P repletion).
      Results are presented as mean ± SD.
      ParameterT1T2T3T4T5
      AP (n = 9)LP (n = 9)AP (n = 9)LP (n = 9)AP (n = 9)LP (n = 9)AP (n = 9)LP (n = 7)AP (n = 9)LP (n = 7)
      Duration (ms)5.2 ± 0.96.2 ± 1.05.3 ± 1.24.9 ± 0.95.0 ± 1.05.1 ± 1.04.9 ± 1.15.0 ± 0.96.3 ± 2.45.7 ± 1.6
      Amplitude (μV)572 ± 98667 ± 203548 ± 161466 ± 142453 ± 138430 ± 145506 ± 139285 ± 80387 ± 103475 ± 303
      No. of phases3.4 ± 0.3
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      3.3 ± 0.5
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      3.1 ± 0.4
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      3.0 ± 0.2
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      3.2 ± 0.4
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      3.0 ± 0.2
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      3.1 ± 0.3
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      Values differ significantly between groups at 1 sampling time (P < 0.05).
      2.7 ± 0.3
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      2.9 ± 0.3
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      3.0 ± 0.4
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      No. of turns3.4 ± 0.53.7 ± 0.9
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      3.1 ± 0.63.0 ± 0.5
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      3.1 ± 0.72.8 ± 0.4
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      3.0 ± 0.42.4 ± 0.3
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      3.1 ± 0.53.0 ± 0.5
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      Size index0.48 ± 0.070.67 ± 0.27
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      0.49 ± 0.300.35 ± 0.21
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      0.39 ± 0.290.35 ± 0.27
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      0.26 ± 0.190.14 ± 0.08
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      0.27 ± 0.180.46 ± 0.45
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      Area (ms/mV)560 ± 151757 ± 300
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      565 ± 243456 ± 190
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      458 ± 223462 ± 220
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      459 ± 199265 ± 64
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      374 ± 156546 ± 438
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      A–C Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      1 Sampling times were T1 (end of acclimation), T2 (after 2 wk of dietary P deprivation antepartum), T3 (within 1 wk postpartum), T4 (at the end of dietary P deprivation), and T5 (after 2 wk of dietary P repletion).
      2 Results are presented as mean ± SD.
      * Values differ significantly between groups at 1 sampling time (P < 0.05).
      Table 6Results of motor unit action potential analysis of intercostal muscles of cows of the second replicate on adequate (AP) and low (LP) dietary P supply stratified by sampling time
      Sampling times were T1 (end of acclimation), T2 (after 2 wk of dietary P deprivation antepartum), T3 (within 1 wk postpartum), T4 (at the end of dietary P deprivation), and T5 (after 2 wk of dietary P repletion).
      Results are presented as mean ± SD.
      ParameterT1T2T3T4
      AP (n = 9)LP (n = 9)AP (n = 9)LP (n = 9)AP (n = 9)LP (n = 9)AP (n = 9)LP (n = 7)
      Duration (ms)5.2 ± 0.75.7 ± 0.95.5 ± 1.24.9 ± 0.76.0 ± 0.94.5 ± 0.47.3 ± 4.85.1 ± 1.4
      Amplitude (μV)476 ± 54
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      550 ± 85
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      429 ± 121
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      Values differ significantly between groups at 1 sampling time (P< 0.05).
      280 ± 62
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      623 ± 196
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      Values differ significantly between groups at 1 sampling time (P< 0.05).
      400 ± 121
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      368 ± 87
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      247 ± 48
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      No. of phases3.2 ± 0.33.1 ± 0.2
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      3.0 ± 0.22.8 ± 0.2
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      3.3 ± 0.33.1 ± 0.2
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      3.1 ± 0.5
      Values differ significantly between groups at 1 sampling time (P< 0.05).
      2.5 ± 0.3
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      No. of turns3.2 ± 0.33.2 ± 0.3
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      3.1 ± 0.32.8 ± 0.5
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      3.6 ± 0.43.2 ± 0.5
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      3.4 ± 1.12.6 ± 0.3
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      Size index0.36 ± 0.18
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      0.45 ± 0.24
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      0.26 ± 0.21
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      0.08 ± 0.07
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      0.55 ± 0.20
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      Values differ significantly between groups at 1 sampling time (P< 0.05).
      0.19 ± 0.15
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      0.29 ± 0.25
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      Values differ significantly between groups at 1 sampling time (P< 0.05).
      0.06 ± 0.05
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      Area (ms/mV)444 ± 101
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      552 ± 173
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      374 ± 140
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      240 ± 54
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      590 ± 134
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      335 ± 108
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      397 ± 133
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      209 ± 45
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      A–C Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      1 Sampling times were T1 (end of acclimation), T2 (after 2 wk of dietary P deprivation antepartum), T3 (within 1 wk postpartum), T4 (at the end of dietary P deprivation), and T5 (after 2 wk of dietary P repletion).
      2 Results are presented as mean ± SD.
      * Values differ significantly between groups at 1 sampling time (P< 0.05).
      The results of the MUAP analysis of BF muscles of both replicates of cows were pooled for statistical analysis as no replicate number effect was identified; the results are summarized in Table 7. A significant treatment effect was only observed for the MUAP size index (0.53 ± 0.20 in AP vs. 0.38 ± 0.21 in LP, P < 0.0001); the time effect was also significant (P = 0.0011).
      Table 7Results of motor unit action potential analysis of biceps femoris muscles of both replicates of cows on adequate (AP) and low (LP) dietary P supply stratified by sampling time
      Sampling times were T1 (end of acclimation), T2 (after 2 wk of dietary P deprivation antepartum), T3 (within 1 wk postpartum), T4 (at the end of dietary P deprivation), and T5 (after 2 wk of dietary P repletion).
      Results are presented as mean ± SD.
      ParameterT1T2T3T4T5
      AP (n = 18)LP (n = 18)AP (n = 18)LP (n = 18)AP (n = 18)LP (n = 18)AP (n = 18)LP (n = 14)AP (n = 9)LP (n = 7)
      Duration (ms)5.9 ± 1.16.1 ± 0.85.9 ± 1.06.0 ± 0.85.9 ± 1.15.9 ± 0.65.8 ± 1.06.2 ± 0.85.6 ± 096.2 ± 1.0
      Amplitude (μV)579 ± 147507 ± 126586 ± 152503 ± 123594 ± 160501 ± 118576 ± 154493 ± 111554 ± 147437 ± 56
      No. of phases3.2 ± 0.23.1 ± 0.23.2 ± 0.23.1 ± 0.43.2 ± 0.23.1 ± 0.23.2 ± 0.23.1 ± 0.23.2 ± 0.13.1 ± 0.1
      No. of turns3.3 ± 0.43.2 ± 0.43.3 ± 0.43.2 ± 0.43.3 ± 0.53.1 ± 0.43.3 ± 0.53.2 ± 0.33.3 ± 0.43.2 ± 0.3
      Size index0.51 ± 0.18
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      0.48 ± 0.23
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      0.48 ± 0.23
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      0.37 ± 0.16
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      0.59 ± 0.21
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      Values differ significantly between groups at 1 sampling time (P < 0.05).
      0.38 ± 0.18
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      0.52 ± 0.15
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      Values differ significantly between groups at 1 sampling time (P < 0.05).
      0.24 ± 0.18
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      0.66 ± 0.27
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      0.47 ± 0.24
      Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      Area (ms/mV)577 ± 142527 ± 146580 ± 129519 ± 138580 ± 147508 ± 135556 ± 130525 ± 141560 ± 142491 ± 85
      A,B Values of the same treatment with different superscripts at different sampling times are significantly different (P < 0.05, Bonferroni corrected).
      1 Sampling times were T1 (end of acclimation), T2 (after 2 wk of dietary P deprivation antepartum), T3 (within 1 wk postpartum), T4 (at the end of dietary P deprivation), and T5 (after 2 wk of dietary P repletion).
      2 Results are presented as mean ± SD.
      * Values differ significantly between groups at 1 sampling time (P < 0.05).
      The results of the multiple stepwise regression analysis with parameters of the MUAP analysis of IC and BF as dependent variables are summarized in Table 8. Plasma [Pi] was found to be positively associated with MUAP amplitude at T3 in BF muscles and with MUAP duration at T1 and T5 as well as with the number of turns at T5 in IC muscles (Table 8). Other positive associations were identified for plasma creatinine that was associated with MUAP duration, amplitude, and number of phases at various time points in BF muscles and plasma [Ca] at T5 that was associated with the MUAP area and size index of IC muscles (Table 8). In contrast, weak negative associations with plasma [Ca] were identified for MUAP duration and the number of turns in IC muscles at T2.
      Table 8Results of multiple stepwise regression analysis for biceps femoris and intercostal muscles with parameters of the motor unit action potential (MUAP) analysis as dependent variables stratified by sampling time
      Parameters included as independent variables were animal ID, treatment, and, for intercostal muscles, replicate number (all forced into the model) as well as age, lactation number (LN), BW, 305-d milk yield of the previous lactation, and DMI; plasma concentration of inorganic phosphate ([Pi]), calcium ([Ca]), creatinine ([Crea]), nonesterified fatty acids ([NEFA]), and BHB ([BHB]); total bilirubin concentration ([TBil]); the activity in plasma of creatinine kinase ([CK]) and aspartate aminotransferase ([AST]); and the total content of phosphorus (Pww), creatine phosphate (CPhww), creatine (Cww), ATP (ATPww), ADP (ADPww), and AMP (AMPww) in muscle tissue wet weight.
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      Results are given as partial coefficient of determination (pr2) and P-value. Associations are listed only when significant (P < 0.05) and partial coefficient of determination (pr2) ≥ 0.1. The direction of the association is indicated by the sign preceding pr2(+ = positive; − = negative).
      MUAP parameterT1T2T3T4T5
      Biceps femoris
       Duration (ms)[AST]: +pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.24; P = 0.004 [CK]: −pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.18; P = 0.005
      AMPww: −pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.22; P = 0.01
      [Crea]: +pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.14; P = 0.04
      [TBil]: +pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.30; P = 0.0009 [Crea]: +pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.12; P = 0.02
      NA
      No significant association.
       Amplitude (μV)NA[Crea]: −pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.17; P = 0.02
      [Pi]: +pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.29; P = 0.0002
      NANA
       No. of phasesNANAAMPww: −pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.12; P = 0.05
      [BHB]: +pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.16; P = 0.03
      BW: −pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.50; P = 0.004 [Crea]: −pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.15; P = 0.04
       No. of turnsNANANANAADPww: −pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.23; P = 0.03
       Size indexNANA[AST]: +pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.14; P = 0.02
      NA[AST]: −pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.27; P = 0.04 [CK]: +pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.37; P = 0.0006
       Area (ms/mV)NANADMI : +pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.15; P = 0.02 AMPww: −pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.11; P = 0.04
      NANA
      Intercostal
       Duration (ms)[Pi]: +pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.17; P = 0.006 LN: +pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.11; P = 0.02
      [Ca]: −pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.14; P = 0.002
      NALN: +pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.20; P = 0.01 [BHB]: +pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.12; P = 0.03
      CPhww: +pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.31; P = 0.04 [Pi]: +pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.25; P = 0.03
       Amplitude (μV)NANA[BHB]: +pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.16; P = 0.01
      NANA
       No. of phasesLN: +pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.29; P = 0.001
      DMI: −pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.21; P = 0.006
      NANA[NEFA]: −pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.39; P = 0.02 ATPww: +pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.21; P = 0.03
       No. of turnsLN: +pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.17; P = 0.02
      [Ca]: −pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.18; P = 0.02
      NA[NEFA]: +pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.14; P = 0.02
      [Pi]: −pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.36; P = 0.003 CPhww: +pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.29; P = 0.05 [AST]: −pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.16; P = 0.001
       Size indexNANA[BHB]: +pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.18; P = 0.007
      NA[Ca]: +pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.25; P = 0.04 ATPww: +pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.19; P = 0.02
       Area (ms/mV)NANA[BHB]: +pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.14; P = 0.02
      NA[Ca]: +pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.23; P = 0.03 [CK]: +pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.14; P = 0.05 [TBil]: −pr
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      = 0.14; P = 0.009
      1 Parameters included as independent variables were animal ID, treatment, and, for intercostal muscles, replicate number (all forced into the model) as well as age, lactation number (LN), BW, 305-d milk yield of the previous lactation, and DMI; plasma concentration of inorganic phosphate ([Pi]), calcium ([Ca]), creatinine ([Crea]), nonesterified fatty acids ([NEFA]), and BHB ([BHB]); total bilirubin concentration ([TBil]); the activity in plasma of creatinine kinase ([CK]) and aspartate aminotransferase ([AST]); and the total content of phosphorus (Pww), creatine phosphate (CPhww), creatine (Cww), ATP (ATPww), ADP (ADPww), and AMP (AMPww) in muscle tissue wet weight.
      2 Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion).
      3 Results are given as partial coefficient of determination (pr2) and P-value. Associations are listed only when significant (P < 0.05) and partial coefficient of determination (pr2) ≥ 0.1. The direction of the association is indicated by the sign preceding pr2(+ = positive; − = negative).
      4 No significant association.
      Results of the Pearson correlation analysis for IC muscles are summarized in Table 9. Strong positive associations with plasma [Pi] were identified for the MUAP amplitude, number of phases, number of turns, size index, and area at T4 and with MUAP duration and size index at T3 (Table 9). Positive associations were furthermore observed between CPhww and MUAP amplitude, number of phases, size index, and MUAP area at different sampling times (Table 9). Negative associations were identified between plasma [Ca] and the MUAP duration, the number of phases, and number of turns at T2 as well as with MUAP amplitude at T4 (Table 9). The ATPww was negatively associated with the number of phases at T2 and the MUAP area at T4 and T5 (Table 9). Neither Pww nor AMPww were found to be associated with any of the studied MUAP parameters.
      Table 9Results of the Pearson correlation analysis (r; P-values in parentheses) of parameters of the motor unit action potential (MUAP) analysis of intercostal muscles
      Parameters of the MUAP analysis were tested for associations with plasma inorganic phosphate ([Pi]), calcium ([Ca,]), and the content of adenosine triphosphate (ATPww), adenosine diphosphate (ADPww), and creatine phosphate (CPhww) in muscle tissue wet weight.
      Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion). The sample size at T1 to T3 was n = 36, at T4 was n = 32, and at T5 was n = 16.
      ParameterTime[Pi][Ca]ATPwwADPwwCPhww
      DurationT10.39 (P = 0.02) NA
      No association.
      NANANA
      T2NA−0.51 (P = 0.002)NANANA
      T30.39 (P = 0.02)NANANANA
      T4NANANANANA
      T5NANANANANA
      AmplitudeT1NANANANA0.36 (P = 0.03)
      T2NANANANANA
      T3NANANANANA
      T40.54 (P = 0.001)−0.48 (P = 0.005)NANA0.36 (P = 0.05)
      T5NANANANANA
      No. of phasesT1NANANANANA
      T2NA−0.41 (P = 0.01)−0.50 (P = 0.004)−0.57 (P = 0.0008)0.35 (P = 0.05)
      T3NANANANANA
      T40.48 (P = 0.006)NANANANA
      T5NANANANANA
      No. of turnsT1NANANANANA
      T2NA−0.44 (P = 0.007)NA0.38 (P = 0.03)NA
      T3NANANANANA
      T40.36 (P = 0.04)NANANANA
      T5NANANANANA
      Size indexT1NANANANANA
      T2NANANANA0.40 (P = 0.02)
      T30.36 (P = 0.03)NANANANA
      T40.39 (P = 0.03)NANA0.38 (P = 0.04)0.49 (P = 0.006)
      T5NANANANANA
      AreaT1NANANANA0.36 (P = 0.04)
      T2NANANANA0.40 (P = 0.02)
      T3NANANANANA
      T40.53 (P = 0.002)NA−0.39 (P = 0.03)−0.40 (P = 0.03)NA
      T5NANANANANA
      1 Parameters of the MUAP analysis were tested for associations with plasma inorganic phosphate ([Pi]), calcium ([Ca,]), and the content of adenosine triphosphate (ATPww), adenosine diphosphate (ADPww), and creatine phosphate (CPhww) in muscle tissue wet weight.
      2 Sampling times are T1 (end of acclimation period), T2 (2 wk after onset of P depletion), T3 (within 1 wk of parturition), T4 (end of P depletion), and T5 (end of P repletion). The sample size at T1 to T3 was n = 36, at T4 was n = 32, and at T5 was n = 16.
      3 No association.
      The Pearson correlation analysis between parameters of the MUAP analysis of BF muscles and parameters characterizing the extracellular and muscle tissue P homeostasis revealed strong positive associations of the plasma [Pi] with the MUAP size index at T4 (r = 0.68, P < 0.0001). Negative associations were identified between ATPww and MUAP size index at T4 (r = −0.37, P = 0.0471) as well as between AMPww and MUAP duration at T1 (r = −0.41, P = 0.0182) and T5 (r = −0.51, P = 0.0357).
      The results of the interference pattern analysis conducted on IC muscles of cows of the second replicate are summarized in Table 10. Time effects (P = 0.0031) and treatment effects were identified for the amplitude per turn of the action potentials [240 (209–279) for treatment AP vs. 234 (195–232) for treatment LP, P = 0.0055; interquartile range in parentheses] but not for the number of turns per seconds or the M/T. The stepwise regression analysis with number of turns per second as the dependent variable yielded significant associations only at T2, with the strongest associations with ADPww (partial R2 = 0.37, P = 0.04). For the interference pattern ratio (M/T), significant associations were identified for sampling time T1, with the strongest associations with ADPww (partial R2 = 0.38, P = 0.03) and plasma [NEFA] (partial R2 = 0.20, P = 0.05). The Pearson correlation analysis at T1 revealed significant associations of the interference pattern ratio (M/T) with plasma [Ca] (r = 0.60, P = 0.0403), ADPww (r = −0.58, P = 0.0484), and ATPww (r = 0.65, P = 0.0210). At T2, a significant correlation was determined between number of turns per second and ADPww (r = −0.60, P = 0.0237).
      Table 10Results of the interference pattern analysis of intercostal muscles of cows on adequate (AP) and low (LP) dietary P supply of the second replicate, stratified by sampling time
      Sampling times were T1 (acclimation), T2 (2 wk of P deprivation antepartum), T3 (first week of lactation), and T4 (end of P deprivation).
      ParameterT1T2T3T4
      AP (n = 9)LP (n = 9)AP (n = 9)LP (n = 9)AP (n = 9)LP (n = 9)AP (n = 9)LP (n = 7)
      No. of turns
      Values are mean ± SD.
      96 ± 3482 ± 2491 ± 3790 ± 2799 ± 3792 ± 3899 ± 3272 ± 13
      Amplitude
      Values are median (interquartile range).
      233 (205–293)
      Values for AP with different uppercase superscripts at different sampling times differ significantly (P < 0.05, Bonferroni corrected).
      233 (230–278)
      Values for LP with different lowercase superscripts at different sampling times differ significantly (P < 0.05, Bonferroni corrected).
      249 (206–255)
      Values for AP with different uppercase superscripts at different sampling times differ significantly (P < 0.05, Bonferroni corrected).
      Values at 1 sampling time differ significantly between treatments (P < 0.05).
      223 (203–275)
      Values for LP with different lowercase superscripts at different sampling times differ significantly (P < 0.05, Bonferroni corrected).
      222 (213–259)
      Values for AP with different uppercase superscripts at different sampling times differ significantly (P < 0.05, Bonferroni corrected).
      177 (162–217)
      Values for LP with different lowercase superscripts at different sampling times differ significantly (P < 0.05, Bonferroni corrected).
      238 (206–254)
      Values for AP with different uppercase superscripts at different sampling times differ significantly (P < 0.05, Bonferroni corrected).
      226 (209–243)
      Values for LP with different lowercase superscripts at different sampling times differ significantly (P < 0.05, Bonferroni corrected).
      Ratio
      Values are mean ± SD.
      0.37 ± 0.140.33 ± 0.130.41 ± 0.190.44 ± 0.120.40 ± 0.140.40 ± 0.140.46 ± 0.150.35 ± 0.07
      a,b Values for LP with different lowercase superscripts at different sampling times differ significantly (P < 0.05, Bonferroni corrected).
      A,B Values for AP with different uppercase superscripts at different sampling times differ significantly (P < 0.05, Bonferroni corrected).
      1 Sampling times were T1 (acclimation), T2 (2 wk of P deprivation antepartum), T3 (first week of lactation), and T4 (end of P deprivation).
      2 Values are mean ± SD.
      3 Values are median (interquartile range).
      * Values at 1 sampling time differ significantly between treatments (P < 0.05).

      DISCUSSION

      The most remarkable finding of this study was that dietary P deprivation that was associated with pronounced and sustained hypophosphatemia over several weeks was not associated with clinical signs of muscle weakness or with recumbency in any of the cows undergoing dietary P depletion. Clinical periparturient hypocalcemia was observed in 4 instances, all affecting AP cows, and significantly higher plasma [Ca] were measured in LP cows in the weeks before and after calving. The effects of P deprivation on calcium homeostasis have been studied in more detail on a subset of study animals, and the results have been published elsewhere (
      • Cohrs I.
      • Wilkens M.R.
      • Grunberg W.
      Short communication: Effect of dietary phosphorus deprivation in late gestation and early lactation on the calcium homeostasis of periparturient dairy cows.
      ). Although dietary P deprivation did not result in a clinically apparent disruption of muscle function, the results reported here revealed obvious subclinical effects of dietary P deprivation on muscle function as determined by electromyography. These are reflected in increased frequency of pathologic spontaneous activity and the decrease of several MUAP parameters that are indicative of muscle fiber functionality loss.
      Phosphorus deprivation in late gestation and early lactation was not associated with clinically apparent effects until after the first week of lactation, when feed intake depression and lower milk production compared with cows on a diet with adequate P supply were observed (data not shown). A similar outcome on feed intake and productivity has been reported in earlier studies (
      • Valk H.
      • Sebek L.B.J.
      Influence of long-term feeding of limited amounts of phosphorus on dry matter intake, milk production, and body weight of dairy cows.
      ;
      • Puggaard L.
      • Lund P.
      • Liesegang A.
      • Sehested J.
      Long term effect of reduced dietary phosphorus on feed intake and milk yield in dry and lactating dairy cows.
      ). In particular, feed intake depression can be considered as a well-accepted outcome of sustained P depletion in cattle and other species (
      • Knochel J.P.
      Pathophysiology and clinical characteristics of severe hypophosphatemia.
      ;
      • Valk H.
      • Sebek L.B.J.
      Influence of long-term feeding of limited amounts of phosphorus on dry matter intake, milk production, and body weight of dairy cows.
      ;
      • Puggaard L.
      • Lund P.
      • Liesegang A.
      • Sehested J.
      Long term effect of reduced dietary phosphorus on feed intake and milk yield in dry and lactating dairy cows.
      ). All AP cows completed the study, whereas 4 LP cows had to be removed prematurely between the second and fourth weeks of lactation. Three animals developed a potentially life-threatening episode of postparturient hemoglobinuria, requiring blood transfusion. Postparturient hemoglobinuria in dairy cows is characterized by severe intravascular hemolysis and hemoglobinuria and has been empirically associated with hypophosphatemia in fresh dairy cows (
      • Macwilliams P.S.
      • Searcy G.P.
      • Bellamy J.E.
      Bovine postparturient hemoglobinuria: A review of the literature.
      ).
      The results of the plasma biochemical analysis confirmed that the experimental rations used in this study were effective in inducing and sustaining pronounced hypophosphatemia throughout the entire depletion period. In LP cows, plasma [Pi] declined continuously until parturition, when a nadir with mean plasma [Pi] below 1.5 mg/dL and with values below the detection limit (0.6 mg/dL) in several cows was reached. No further decline was observed after calving despite ongoing dietary P depletion. Dietary P supplementation during the repletion period in LP cows of the first replicate was associated with a rapid increase in plasma [Pi], again confirming the efficacy of oral P supplementation to correct hypophosphatemia. Dietary P supplementation in LP cows was furthermore associated with an increase of MUAP parameters that apparently were affected by previous P deprivation. This observation suggests that the loss of contributing individual muscle fibers to the MUAP observed during P deprivation can be corrected at least partially with dietary P supplementation (
      • Kimura J.
      Electronic systems and data analysis.
      ;
      • Daube J.R.
      • Rubin D.I.
      Needle electromyography.
      ;
      • Wijnberg I.D.
      • Franssen H.
      The potential and limitations of quantitative electromyography in equine medicine.
      ).
      The biochemical analysis of muscle tissue revealed that P depletion and hypophosphatemia did not result in a remarkable disturbance of the muscle tissue P homeostasis. A treatment effect could not be identified for any of the biochemical parameters studied in muscle tissue. The Pww that comprises all P present in muscle tissue remained unchanged from acclimation through the entire depletion phase. Interestingly, a decline in Pww was observed at the end of the P supplementation phase in LP cows. This decline in Pww between T4 and T5 was associated with a decline in ATPww and ADPww, 2 compounds forming part of the muscle tissue Pww pool. Although there is no evident explanation for the changes observed during the P repletion phase, a similar decline in ATPww and ADPww during dietary P repletion following a depletion phase of 5 wk was observed in an earlier study (
      • Grünberg W.
      • Scherpenisse P.
      • Dobbelaar P.
      • Idink M.J.
      • Wijnberg I.D.
      The effect of transient, moderate dietary phosphorus deprivation on phosphorus metabolism, muscle content of different phosphorus-containing compounds, and muscle function in dairy cows.
      ). This observed decline in Pww during repletion was furthermore tightly associated with a concomitant decline of Kww of a similar extent. Parallel changes in the content of K and P in muscle tissue, both predominantly intracellular electrolytes, imply an underlying mechanism that is not specific to P. Such a mechanism could, for example, be an increase in cytosolic volume of muscle cells that would increase the distribution volume of compounds dissolved in the cytosol to a similar degree. In this study, the muscle tissue DM content was used as an indicator for changes in tissue water content, giving a crude approximation for possible changes in cytosol volume. Muscle tissue DM decreased and thus water content in muscle tissue increased between T3 and T5 in both groups. This effect contributed to the observed decrease in Pww, Kww, ATPww, and ADPww at T5 as implied by the results of the regression and correlation analyses. Changes in tissue water content in the form of an increase in empty body water (i.e., total body water excluding water contained in the digestive tract) in late gestation and early lactation have previously been reported in cattle and ewes (
      • Degen A.A.
      • Young B.A.
      Live-weight, total body-water and maternal body-solid changes in pregnant and lactating beef cows.
      ;
      • Robinson J.J.
      Changes in body composition during pregnancy and lactation.
      ). This effect has been attributed to increased water turnover at the onset of lactation as well as to the rapid loss of body fat around parturition that is tightly coupled with increased tissue water content.
      Creatine phosphate was the only P-containing compound in muscle tissue with a time curve crudely reflecting the different phases of dietary P supply in P-deprived animals. However, similar changes were observed in AP cows, again suggesting that these effects are not directly attributable to P deprivation.
      Associations between plasma [Pi] and intracellular P compounds were only identified at T4. At this time, the Pearson correlation analysis revealed correlations of plasma [Pi] with Pww, ATPww, and ADPww that were negative. These results are in agreement with earlier studies and indicate that plasma [Pi] is poorly suited to estimate muscle tissue P content. These findings specifically refute the concept that hypophosphatemia is indicative of intracellular P depletion of muscle tissue (
      • Grünberg W.
      • Scherpenisse P.
      • Dobbelaar P.
      • Idink M.J.
      • Wijnberg I.D.
      The effect of transient, moderate dietary phosphorus deprivation on phosphorus metabolism, muscle content of different phosphorus-containing compounds, and muscle function in dairy cows.
      ). Maintenance of constant P content in tissue despite variable P supply and plasma [Pi] has also been reported for other tissues such as liver and erythrocytes (
      • Grünberg W.
      • Staufenbiel R.
      • Constable P.D.
      • Dann H.M.
      • Morin D.E.
      • Drackley J.K.
      Liver phosphorus content in Holstein-Friesian cows during the transition period.
      ,
      • Grünberg W.
      • Mol J.A.
      • Teske E.
      Red blood cell phosphate concentration and osmotic resistance during dietary phosphate depletion in dairy cows.
      ).
      Plasma total bilirubin was found to be positively associated with all P compounds in muscle tissue at one or several sampling times. Animals with postparturient hemoglobinuria, a condition characterized by intravascular hemolysis and invariably associated with hyperbilirubinemia, are unlikely to have confounded these results because associations with total bilirubin were already observed antepartum. Furthermore, samples from animals with postparturient hemoglobinuria obtained at T3 and later were not included in this analysis. Plasma total bilirubin at T3 did not differ between groups but showed a transient mild increase in both groups during the second and third weeks of lactation. Mild hyperbilirubinemia may be indicative of impaired liver function frequently observed in early-lactating dairy cows. Because positive associations with total bilirubin were identified at T1 and T2, considerably before liver function is commonly hampered, hepatic hyperbilirubinemia also is deemed an improbable confounder. Muscle tissue contributes approximately 20% to plasma total bilirubin through its myoglobin turnover (
      • Kamga C.
      • Krishnamurthy S.
      • Shiva S.
      Myoglobin and mitochondria: A relationship bound by oxygen and nitric oxide.
      ). As the amount of muscle tissue myoglobin characterizes the capacity of aerobic glycolysis of myocytes, this may explain the positive association between compounds of adenosine phosphate metabolism and the metabolite of the myoglobin breakdown.
      The amounts of P, K, ATP, ADP, AMP, and CPh in muscle tissue in this study are in good agreement with values determined in muscle tissue from bovines and other species in earlier studies (
      • Fuller T.J.
      • Carter N.W.
      • Barcenas C.
      • Knochel J.P.
      Reversible changes of the muscle cell in experimental phosphorus deficiency.
      ;
      • Angerås U.
      • Hallangeras M.
      • Wagner K.R.
      • James H.
      • Hasselgren P.O.
      • Fischer J.E.
      Tissue metabolite levels in different types of skeletal muscle during sepsis.
      ;
      • Williams S.N.
      • McDowell L.R.
      • Warnick A.C.
      • Wilkinson N.S.
      • Lawrence L.A.
      Phosphorus concentrations in blood, milk, feces, bone and selected fluids and tissues of growing heifers as affected by dietary phosphorus.
      ;
      • Grünberg W.
      • Scherpenisse P.
      • Dobbelaar P.
      • Idink M.J.
      • Wijnberg I.D.
      The effect of transient, moderate dietary phosphorus deprivation on phosphorus metabolism, muscle content of different phosphorus-containing compounds, and muscle function in dairy cows.
      ;
      • Schneider S.
      • Muller A.
      • Wittek T.
      Concentration of potassium in plasma, erythrocytes, and muscle tissue in cows with decreased feed intake and gastrointestinal ileus.
      ).
      An important implication of the results reported here is that the degree and duration of P deprivation achieved in this study were not sufficient to negatively affect the muscle tissue P homeostasis or the adenosine phosphate metabolism in periparturient dairy cows. It is highly improbable that the degree of dietary P deprivation that was achieved in the present study would occur in dairy cattle under field conditions in Europe or North America. Studies conducted in other species and using P-deficient diets covering 10% of daily P requirements or less over several weeks reported a decline of the P content not only in blood but also in muscle tissue. Dietary P deprivation in these studies resulted in a reduction in muscle tissue P content in the range of 20% after 5 wk in dogs and of 45% after 4 wk in mice (
      • Fuller T.J.
      • Carter N.W.
      • Barcenas C.
      • Knochel J.P.
      Reversible changes of the muscle cell in experimental phosphorus deficiency.
      ;
      • Hettleman B.D.
      • Sabina R.L.
      • Drezner M.K.
      • Holmes E.W.
      • Swain J.L.
      Defective adenosine triphosphate synthesis—An explanation for skeletal muscle dysfunction in phosphate-deficient mice.
      ). It has been proposed that the intracellular P homeostasis is resistant to changes of the plasma [Pi] up to a certain degree due to the properties of the Na-linked P transport ensuring efficient cellular P uptake from the extracellular space despite an unfavorable chemical gradient (
      • Kemp G.J.
      Abnormalities of Pi concentration in plasma and in cells.
      ). Dietary P deprivation of comparable severity in dairy cattle was associated with rapid induction of counterregulation in the form of bone mobilization releasing Pi together with calcium from bone in mid-lactating dairy cows (
      • Grünberg W.
      • Scherpenisse P.
      • Dobbelaar P.
      • Idink M.J.
      • Wijnberg I.D.
      The effect of transient, moderate dietary phosphorus deprivation on phosphorus metabolism, muscle content of different phosphorus-containing compounds, and muscle function in dairy cows.
      ). In this previous study, plasma [Pi] was observed to increase after 10 d of dietary P depletion, which was attributed to the efficacy of counterregulation. In the present study, an increase in plasma [Pi] was not observed until P was supplemented orally. However, the lower incidence of clinical hypocalcemia and higher plasma [Ca] around parturition in LP cows compared with AP cows suggest that bone mobilization was indeed triggered by dietary P deprivation. It is conceivable that in the present study the combined challenge of P homeostasis through dietary P deprivation and the onset of lactation resulted in increased bone mobilization in combination with enhanced cellular P uptake to maintain intracellular P homeostasis. Such a scenario of increased P flux would preclude an increase of plasma [Pi] despite bone mobilization while intracellular P homeostasis is maintained.
      In the present study, experimental animals were fed dry cow rations with moderately elevated DCAD, which is more common in Europe than in the United States. Moderately alcalogenic diets are associated with an increased risk for clinical and subclinical periparturient hypocalcemia, and hypocalcemia is thought to contribute to and exacerbate periparturient hypophosphatemia through increased parathyroid (PTH) secretion in hypocalcemic animals (
      • Horst R.L.
      Regulation of calcium and phosphorus homeostasis in the dairy cow.
      ). Following this hypothesis, alcalogenic diets used in this study thus may have exacerbated periparturient hypophosphatemia. It should however be noted that results obtained in a subset of study animals and published elsewhere revealed that cows on LP diets had markedly lower plasma PTH concentrations and higher plasma [Ca] than AP cows around parturition, indicating that dietary P depletion in late gestation was associated with markedly hampered PTH release. This would at least have attenuated the hypothesized exacerbation of periparturient hypophosphatemia by high-DCAD diets in AP cows (
      • Cohrs I.
      • Wilkens M.R.
      • Grunberg W.
      Short communication: Effect of dietary phosphorus deprivation in late gestation and early lactation on the calcium homeostasis of periparturient dairy cows.
      ).
      Results of the electromyographic examination indicate that overall, functionality of muscle fibers was decreased in P-deprived animals, as can be concluded from the decrease in size index and lower MUAP amplitudes over time in P-deprived animals. Furthermore, pathologic spontaneous activity that is indicative of instable muscle fiber or nerve membranes became more frequent in LP cows as P deprivation progressed, especially in IC muscles. Increasing occurrence of pathologic spontaneous activity was previously reported in mid-lactating dairy cows that were fed a P-deficient diet over a course of 5 wk (
      • Grünberg W.
      • Scherpenisse P.
      • Dobbelaar P.
      • Idink M.J.
      • Wijnberg I.D.
      The effect of transient, moderate dietary phosphorus deprivation on phosphorus metabolism, muscle content of different phosphorus-containing compounds, and muscle function in dairy cows.
      ). Membrane instability developing in states of P deficiency may be attributable to disturbed synthesis of phospholipids that form an integral part of cell membranes that is well documented (
      • Subramanian R.
      • Khardori R.
      Severe hypophosphatemia—Pathophysiologic implications, clinical presentations, and treatment.
      ;
      • Chen Y.Y.
      • Kao T.W.
      • Chou C.W.
      • Wu C.J.
      • Yang H.F.
      • Lai C.H.
      • Wu L.W.
      • Chen W.L.
      Exploring the link between serum phosphate levels and low muscle strength, dynapenia, and sarcopenia.
      ). Dietary P supplementation over a course of 2 wk resulted in partial reversal of the observed effects of P deprivation on MUAP parameters, whereas the occurrence of pathologic spontaneous activity remained unchanged or even progressed during the repletion period of this study. These findings indicate that in contrast to disturbances of muscle fiber functionality that seem to respond to correction of P depletion, the effects of membrane instability reflecting in increased pathologic spontaneous activity are not easily reversed.
      Parameters of the MUAP analysis seemed to be associated with the plasma [Pi] rather than with the content of energy-carrying compounds of muscle tissue. This observation suggests that altered muscle fiber function is indeed associated with hypophosphatemia per se and thus extracellular P homeostasis rather than with depletion of energy stores of muscle tissue or disturbances of the intracellular P balance of muscle fibers. The negative association of ATP metabolites with several MUAP parameters identified in this study suggests that more and larger muscle fibers are activated with lower ATP availability to maintain comparable muscle strength (
      • Kimura J.
      Electronic systems and data analysis.
      ;
      • Wijnberg I.D.
      • Franssen H.
      The potential and limitations of quantitative electromyography in equine medicine.
      ). This interpretation is supported by the positive association of the number of turns per second and the M/T of the interference pattern analysis and ATPww, AMPww, CPhww, and plasma [Ca] determined at T1 and T2 but not at later time points. This proposed hypothesis implies a loss of relation between muscle strength and available ATP, AMP, and calcium in muscle tissue in states of P deprivation. Calcium, ATP, and AMP are all required for muscle contraction; the logically ensuing positive association of these parameters with interference pattern parameters characterizing muscle strength could be identified early in the study but was lost or even inversed as P depletion progressed. These results thus partly support the concept of disturbed muscle function in states of P deprivation.
      Several of the MUAP parameters were found to be positively associated with plasma [Pi], with the MUAP size index showing one of the most consistent responses to P depletion. It is this parameter that is considered one of the more robust in human electromyographic examination (
      • Sonoo M.
      • Stalberg E.
      The ability of MUP parameters to discriminate between normal and neurogenic MUPs in concentric EMG: Analysis of the MUP “thickness” and the proposal of “size index.”.
      ;
      • Wijnberg I.D.
      • Graubner C.
      • Auriemma E.
      • van de Belt A.J.
      • Gerber V.
      Quantitative motor unit action potential analysis in 2 paraspinal neck muscles in adult Royal Dutch Sport horses.
      ;
      • Wijnberg I.D.
      • Franssen H.
      The potential and limitations of quantitative electromyography in equine medicine.
      ). The observed correlations between plasma [Pi] and MUAP parameters at the end of the P-deprivation period indicate that the lower the plasma [Pi], the lower the number of functional muscle fibers. Indeed, several authors mention muscle cell membrane and structural integration disturbances as causes of myopathy noticed in hypophosphatemic patients (
      • Bugg N.C.
      • Jones J.A.
      Hypophosphataemia. Pathophysiology, effects and management on the intensive care unit.
      ;
      • Subramanian R.
      • Khardori R.
      Severe hypophosphatemia—Pathophysiologic implications, clinical presentations, and treatment.
      ;
      • Chen Y.Y.
      • Kao T.W.
      • Chou C.W.
      • Wu C.J.
      • Yang H.F.
      • Lai C.H.
      • Wu L.W.
      • Chen W.L.
      Exploring the link between serum phosphate levels and low muscle strength, dynapenia, and sarcopenia.
      ). If muscle fibers do not contribute to the MUAP, the summation of the individual action potential parameters decreases, which is reflected in lower values of MUAP parameters such as MUAP amplitude or size index. The positive association of MUAP parameters with plasma [Pi] at the time alterations of the MUAP were most prominent, whereas association related to intracellular P homeostasis such as ATPww or CPhww were absent, suggesting that disturbances of extracellular P homeostasis may etiologically contribute to the observed effects even without measurable changes of the intracellular P homeostasis. In this context, it should be noted that most prominent effects were identified on IC muscles while muscle tissue specimens for biopsy were obtained from BF muscle exclusively. In view of the obvious difference in sensitivity to P depletion of these 2 muscles, it is conceivable that the time–concentration curves for intracellular P compounds in BF muscles presented in this study may not be representative of IC muscles. The results of the correlation and regression analysis between parameters of the MUAP analysis in IC muscles and P-containing compounds of BF muscles presented in Table 8, Table 9 should thus be interpreted cautiously. Hypocalcemia was also found to have an effect on MUAP parameters as a result of increased temporal dispersion leading to increased MUAP duration, number of turns, and number of phases (
      • Wijnberg I.D.
      • van der Kolk J.H.
      • Franssen H.
      • Breukink H.J.
      Electromyographic changes of motor unit activity in horses with induced hypocalcemia and hypomagnesemia.
      ). To complicate the interpretation of the results of the regression and correlation analysis, it must be kept in mind that plasma [Pi] and muscle P content are not associated with each other, a finding corroborated by the results presented here (
      • Grünberg W.
      • Scherpenisse P.
      • Dobbelaar P.
      • Idink M.J.
      • Wijnberg I.D.
      The effect of transient, moderate dietary phosphorus deprivation on phosphorus metabolism, muscle content of different phosphorus-containing compounds, and muscle function in dairy cows.
      ). Furthermore, important diurnal variation in plasma [Pi] occurs in mammals and in cattle in particular (
      • Bugg N.C.
      • Jones J.A.
      Hypophosphataemia. Pathophysiology, effects and management on the intensive care unit.
      ;
      • Grünberg W.
      Treatment of phosphorus balance disorders.
      ).
      The differences between BF and IC muscles observed in this study might be explained by differences in fiber type composition of these muscles. Intercostal muscle that require aerobic activity at all times contain more oxidative type I fibers than BF muscles (
      • Greer J.J.
      • Martin T.P.
      Distribution of muscle fiber types and EMG activity in cat intercostal muscles.
      ), and hypophosphatemia has been associated with impaired oxygen release and tissue hypoxia in monogastric species (
      • Subramanian R.
      • Khardori R.
      Severe hypophosphatemia—Pathophysiologic implications, clinical presentations, and treatment.
      ). It is thus conceivable that aerobic muscle fibers are more susceptible to P depletion than more anaerobic muscle fibers. Also, IC muscles are activated consistently with constant strength during voluntary breathing, which is associated with a more consistent recruitment pattern of motor units than in the BF muscle. This more uniform recruitment pattern and the presumably smaller degree of variation associated with this constant recruitment pattern is likely to facilitate the disclosure of effects compared with BF muscles, where contraction had to be triggered by lifting the contralateral hind limb.
      The reported study has several limitations. Most important, housing study animals in tiestalls, which strongly limited physical exercise, may have contributed to a certain degree to the observed subclinical effects on muscle function. Part of the increase in frequency of pathologic spontaneous activity observed in both treatments could have been associated with changes occurring around parturition but may as well have been the result of prolonged exercise restriction. Also, the chosen study design does not allow even remotely estimating a threshold for dietary P depletion required to obtain observed subclinical effects.
      Marked and sustained P depletion in transition cows associated with severe hypophosphatemia did not cause clinically apparent signs of muscle weakness or recumbency despite electromyographic abnormalities suggestive of subclinical myopathic alterations. Furthermore, P depletion and hypophosphatemia did not measurably disturb the P balance of muscle tissue in this study, as determined in the example of the BF muscle. The increasing occurrence of pathologic spontaneous activity in P-deprived animals indicates disturbances of muscle fiber function as P deficiency progressed, possibly resulting from membrane instability of muscle or nerve cells. The association of MUAP parameters with the plasma [Pi] indicates that disturbed muscle function and altered strength of muscle fiber contractility are associated with disturbances of P homeostasis. The absence of clinically apparent muscle weakness despite severe hypophosphatemia does not support the concept of a role of P deprivation in periparturient recumbency of dairy cows. It remains to be determined whether the observed subclinical effects on neuromuscular function are relevant for the health and well-being of dairy cows.

      ACKNOWLEDGMENTS

      The authors acknowledge the dedicated technical assistance of the caretakers of the Department of Farm Animal Health at Utrecht University (Utrecht, the Netherlands). This study was funded by ZuivelNL (Zoetermeer, the Netherlands).

      REFERENCES

        • Angerås U.
        • Hallangeras M.
        • Wagner K.R.
        • James H.
        • Hasselgren P.O.
        • Fischer J.E.
        Tissue metabolite levels in different types of skeletal muscle during sepsis.
        Metabolism. 1991; 40 (1943743): 1147-1151
        • Bugg N.C.
        • Jones J.A.
        Hypophosphataemia. Pathophysiology, effects and management on the intensive care unit.
        Anaesthesia. 1998; 53 (9849285): 895-902
        • Chen Y.Y.
        • Kao T.W.
        • Chou C.W.
        • Wu C.J.
        • Yang H.F.
        • Lai C.H.
        • Wu L.W.
        • Chen W.L.
        Exploring the link between serum phosphate levels and low muscle strength, dynapenia, and sarcopenia.
        Sci. Rep. 2018; 8 (29476104): 3573
        • Cohrs I.
        • Wilkens M.R.
        • Grunberg W.
        Short communication: Effect of dietary phosphorus deprivation in late gestation and early lactation on the calcium homeostasis of periparturient dairy cows.
        J. Dairy Sci. 2018; 101 (30100496): 9591-9598
        • Daube J.R.
        • Rubin D.I.
        Needle electromyography.
        Muscle Nerve. 2009; 39 (19145648): 244-270
        • Degen A.A.
        • Young B.A.
        Live-weight, total body-water and maternal body-solid changes in pregnant and lactating beef cows.
        J. Agric. Sci. 1980; 95: 1-5
        • Finsterer J.
        • Mamoli B.
        • Fuglsang-Frederiksen A.
        Peak-ratio interference pattern analysis in the detection of neuromuscular disorders.
        Electroencephalogr. Clin. Neurophysiol. 1997; 105 (9363003): 379-384
        • Fuller T.J.
        • Carter N.W.
        • Barcenas C.
        • Knochel J.P.
        Reversible changes of the muscle cell in experimental phosphorus deficiency.
        J. Clin. Invest. 1976; 57 (947947): 1019-1024
        • Georgesco M.
        • Salerno A.
        Spontaneous activity in electromyography: Mechanisms and practical interest.
        Neurophysiol Clin. 2000; 30 (11126641): 289-305
        • Goff J.P.
        Treatment of calcium, phosphorus, and magnesium balance disorders.
        Vet. Clin. North Am. Food Anim. Pract. 1999; 15 (10573815): 619-639
        • Goff J.P.
        Invited review: Mineral absorption mechanisms, mineral interactions that affect acid-base and antioxidant status, and diet considerations to improve mineral status.
        J. Dairy Sci. 2018; 101 (29397180): 2763-2813
        • Greer J.J.
        • Martin T.P.
        Distribution of muscle fiber types and EMG activity in cat intercostal muscles.
        J. Appl. Physiol. 1990; 69: 1208-1211
        • Grünberg W.
        Treatment of phosphorus balance disorders.
        Vet. Clin. North Am. Food Anim. Pract. 2014; 30 (24980728): 383-408
        • Grünberg W.
        • Mol J.A.
        • Teske E.
        Red blood cell phosphate concentration and osmotic resistance during dietary phosphate depletion in dairy cows.
        J. Vet. Intern. Med. 2015; 29 (25407950): 395-399
        • Grünberg W.
        • Scherpenisse P.
        • Dobbelaar P.
        • Idink M.J.
        • Wijnberg I.D.
        The effect of transient, moderate dietary phosphorus deprivation on phosphorus metabolism, muscle content of different phosphorus-containing compounds, and muscle function in dairy cows.
        J. Dairy Sci. 2015; 98 (26026765): 5385-5400
        • Grünberg W.
        • Staufenbiel R.
        • Constable P.D.
        • Dann H.M.
        • Morin D.E.
        • Drackley J.K.
        Liver phosphorus content in Holstein-Friesian cows during the transition period.
        J. Dairy Sci. 2009; 92 (19389968): 2106-2117
        • Hettleman B.D.
        • Sabina R.L.
        • Drezner M.K.
        • Holmes E.W.
        • Swain J.L.
        Defective adenosine triphosphate synthesis—An explanation for skeletal muscle dysfunction in phosphate-deficient mice.
        J. Clin. Invest. 1983; 72 (6874957): 582-589
        • Horst R.L.
        Regulation of calcium and phosphorus homeostasis in the dairy cow.
        J. Dairy Sci. 1986; 69 (3517093): 604-616
        • Jiang Y.
        • Sun C.
        • Ding X.
        • Yuan D.
        • Chen K.
        • Gao B.
        • Chen Y.
        • Sun A.
        Simultaneous determination of adenine nucleotides, creatine phosphate and creatine in rat liver by high performance liquid chromatography-electrospray ionization-tandem mass spectrometry.
        J. Pharm. Biomed. Anal. 2012; 66 (22494519): 258-263
        • Kamga C.
        • Krishnamurthy S.
        • Shiva S.
        Myoglobin and mitochondria: A relationship bound by oxygen and nitric oxide.
        Nitric Oxide. 2012; 26: 251-258
        • Kemp G.J.
        Abnormalities of Pi concentration in plasma and in cells.
        Clin. Chem. 1993; 39 (8375098): 2028-2031
        • Kimura J.
        Electronic systems and data analysis.
        in: Kimura J. Electrodiagnosis in Diseases of Nerve and Muscle: Principles and Practice. Oxford University Press, Oxford, UK2001: 39-59
        • Knochel J.P.
        Pathophysiology and clinical characteristics of severe hypophosphatemia.
        Arch. Intern. Med. 1977; 137 (836118): 203-220
        • Liguori R.
        • Fuglsang-Frederiksen A.
        • Nix W.
        • Fawcett P.R.W.
        • Andersen K.
        Electromyography in myopathy.
        Neurophysiol. Clin. 1997; 27 (9260160): 200-203
        • Macrae A.I.
        • Burrough E.
        • Forrest J.
        Assessment of nutrition in dairy herds: Use of metabolic profiles.
        Cattle Pract. 2012; 20: 120-127
        • Macrae A.I.
        • Whitaker D.A.
        • Burrough E.
        • Dowell A.
        • Kelly J.M.
        Use of metabolic profiles for the assessment of dietary adequacy in UK dairy herds.
        Vet. Rec. 2006; 159 (17099173): 655-661
        • Macwilliams P.S.
        • Searcy G.P.
        • Bellamy J.E.
        Bovine postparturient hemoglobinuria: A review of the literature.
        Can. Vet. J. 1982; 23 (17422195): 309-312
        • Ménard L.
        • Thompson A.
        Milk fever and alert downer cows: Does hypophosphatemia affect the treatment response?.
        Can. Vet. J. 2007; 48 (17542366): 487-491
        • Nirkko A.C.
        • Rosler K.M.
        • Hess C.W.
        Sensitivity and specificity of needle electromyography: A prospective study comparing automated interference pattern analysis with single motor unit potential analysis.
        Electroencephalogr. Clin. Neurophysiol. 1995; 97 (7533715): 1-10
        • NRC
        Nutrient Requirements of Dairy Cattle.
        7th rev. ed. Natl. Acad. Sci., Washington, DC2001
        • Puggaard L.
        • Lund P.
        • Liesegang A.
        • Sehested J.
        Long term effect of reduced dietary phosphorus on feed intake and milk yield in dry and lactating dairy cows.
        Livest. Sci. 2014; 159: 18-28
        • Robinson J.J.
        Changes in body composition during pregnancy and lactation.
        Proc. Nutr. Soc. 1986; 45 (3703866): 71-80
        • Rodehutscord M.
        • Pauen A.
        • Windhausen P.
        • Brintrup R.
        • Pfeffer E.
        Effects of drastic changes in P intake on P concentrations in blood and rumen fluid of lactating ruminants.
        Zentralbl. Veterinarmed A. 1994; 41 (7732738): 611-619
        • Rubin D.I.
        Needle electromyography: Basic concepts and patterns of abnormalities.
        Neurol. Clin. 2012; 30 (22361369): 429-456
        • Schneider S.
        • Muller A.
        • Wittek T.
        Concentration of potassium in plasma, erythrocytes, and muscle tissue in cows with decreased feed intake and gastrointestinal ileus.
        J. Vet. Intern. Med. 2016; 30 (26875482): 679-685
        • Sonoo M.
        • Stalberg E.
        The ability of MUP parameters to discriminate between normal and neurogenic MUPs in concentric EMG: Analysis of the MUP “thickness” and the proposal of “size index.”.
        Electroencephalogr. Clin. Neurophysiol. 1993; 89 (7691568): 291-303
        • Subramanian R.
        • Khardori R.
        Severe hypophosphatemia—Pathophysiologic implications, clinical presentations, and treatment.
        Medicine (Baltimore). 2000; 79 (10670405): 1-8
        • Tamminga S.
        • Vanstraalen W.M.
        • Subnel A.P.J.
        • Meijer R.G.M.
        • Steg A.
        • Wever C.J.G.
        • Blok M.C.
        The Dutch protein evaluation system—The DVE/OEB-system.
        Livest. Prod. Sci. 1994; 40: 139-155
        • Teerlink T.
        • Hennekes M.
        • Bussemaker J.
        • Groeneveld J.
        Simultaneous determination of creatine compounds and adenine nucleotides in myocardial tissue by high-performance liquid chromatography.
        Anal. Biochem. 1993; 214 (8250235): 278-283
        • Valk H.
        • Sebek L.B.J.
        Influence of long-term feeding of limited amounts of phosphorus on dry matter intake, milk production, and body weight of dairy cows.
        J. Dairy Sci. 1999; 82 (10531601): 2157-2163
        • Wijnberg I.D.
        • Back W.
        • De Jong M.
        • Zuidhof M.C.
        • Van Den Belt A.J.M.
        • Van Der Kolk J.H.
        The role of electromyography in clinical diagnosis of neuromuscular locomotor problems in the horse.
        Equine Vet. J. 2004; 36 (15656503): 718-722
        • Wijnberg I.D.
        • Franssen H.
        The potential and limitations of quantitative electromyography in equine medicine.
        Vet. J. 2016; 209: 23-31
        • Wijnberg I.D.
        • Franssen H.
        • van der Kolk J.H.
        • Breukink H.J.
        Quantitative analysis of motor unit action potentials in the subclavian muscle of healthy horses.
        Am. J. Vet. Res. 2002; 63 (11843118): 198-203
        • Wijnberg I.D.
        • Graubner C.
        • Auriemma E.
        • van de Belt A.J.
        • Gerber V.
        Quantitative motor unit action potential analysis in 2 paraspinal neck muscles in adult Royal Dutch Sport horses.
        J. Vet. Intern. Med. 2011; 25 (21488965): 592-597
        • Wijnberg I.D.
        • van der Kolk J.H.
        • Franssen H.
        • Breukink H.J.
        Electromyographic changes of motor unit activity in horses with induced hypocalcemia and hypomagnesemia.
        Am. J. Vet. Res. 2002; 63 (12061532): 849-856
        • Williams S.N.
        • McDowell L.R.
        • Warnick A.C.
        • Wilkinson N.S.
        • Lawrence L.A.
        Phosphorus concentrations in blood, milk, feces, bone and selected fluids and tissues of growing heifers as affected by dietary phosphorus.
        Livest. Res. Rural Dev. 1991; 3: 1-14