Paget's disease of bone (PDB) is a chronic skeletal disorder that can affect one or several bones in individuals older than 55 y of age. PDB-like changes have been reported in archaeological remains as old as Roman, although accurate diagnosis and natural history of the disease is lacking. Six skeletons from a collection of 130 excavated at Norton Priory in the North West of England, which dates to medieval times, show atypical and extensive pathological changes resembling contemporary PDB affecting as many as 75% of individual skeletons. Disease prevalence in the remaining collection is high, at least 16% of adults, with age at death estimations as low as 35 y. Despite these atypical features, paleoproteomic analysis identified sequestosome 1 (SQSTM1) or p62, a protein central to the pathological milieu of PDB, as one of the few noncollagenous human sequences preserved in skeletal samples. Targeted proteomic analysis detected >60% of the ancient p62 primary sequence, with Western blotting indicating p62 abnormalities, including in dentition. Direct sequencing of ancient DNA excluded contemporary PDB-associated SQSTM1 mutations. Our observations indicate that the ancient p62 protein is likely modified within its C-terminal ubiquitin-associated domain. Ancient miRNAs were remarkably preserved in an osteosarcoma from a skeleton with extensive disease, with miR-16 expression consistent with that reported in contemporary PDB-associated bone tumors. Our work displays the use of proteomics to inform diagnosis of ancient diseases such as atypical PDB, which has unusual features presumably potentiated by yet-unidentified environmental or genetic factors.

Small for gestational age (SGA) infants are at increased risk of type 2 diabetes in adulthood. It is unknown whether any prenatal biomarkers are helpful for identifying SGA infants with altered metabolic health profile at birth or later life. In a nested study of 162 SGA (birth weight < 10th percentile) and 161 optimal birth weight (25th-75th percentiles) control infants in the 3D (design, develop and discover) birth cohort in Canada, we assessed whether maternal circulating placental growth factor (PlGF), a biomarker of placental function, is associated with metabolic health biomarkers in SGA infants. Main outcomes were cord plasma insulin, proinsulin, insulin-like growth factor-I (IGF-I), leptin, and high-molecular weight (HMW) adiponectin concentrations. Maternal PlGF concentrations at 32-35 weeks of gestation were substantially lower in SGA versus control infants (P < 0.001), so as were cord plasma proinsulin (P = 0.005), IGF-I (P < 0.001), leptin (P < 0.001), and HMW adiponectin (P = 0.002) concentrations. In SGA infants with both low (<25th percentile) and normal maternal PlGF concentrations, cord plasma IGF-I and leptin concentrations were lower than control infants, but the decreases were to a greater extent in SGA infants with low maternal PlGF. Cord blood leptin levels were lower comparing SGA infants with low vs. normal maternal PlGF levels (P = 0.01). SGA infants with low maternal circulating PlGF levels at late gestation were characterized by greater decreases in cord blood IGF-I and leptin concentrations. Maternal circulating PlGF appears to be associated with neonatal metabolic health profile in SGA infants.

Accelerated growth in postnatal life in low birth weight infants has been associated with insulin resistance and metabolic syndrome-related disorders in later life. Postnatal accelerated growth in also common in normal birth weight infants, but little is known about the impact on metabolic health. In a prospective cohort study of 203 term normal birth weight infants, we evaluated the impacts of accelerated (Δweight Z score > 0.5) or decelerated (Δweight ΔZ < -0.5) growth during early (0-3 months) and late (3-12 months) postnatal life on metabolic health indicators at age 1-year. The primary outcomes were homeostasis model assessment of insulin resistance (HOMA-IR), β-cell function [homeostasis model assessment of β-cell function (HOMA-β)], and fasting plasma lipids. Adjusting for maternal, paternal, and infant characteristics, accelerated growth during the first 3 months of life was associated with a 41.6% (95% confidence interval 8.9-84.2%) increase in HOMA-β, and a 8.3% (0.7-15.4%) decrease in fasting plasma total cholesterols, and was not associated with HOMA-IR in infants at age 1-year. Accelerated growth during 3-12 months was associated with a 30.9% (3.3-66.0%) increase in HOMA-IR and was not associated with HOMA-β. Neither accelerated nor decelerated growth was associated with fasting plasma triglycerides, high-density lipoprotein or low-density lipoprotein cholesterol concentrations in infants at age 1-year. Accelerated growth during early postnatal life may be beneficial for β-cell function, but during late postnatal life harmful for insulin sensitivity in normal birth weight infants.

Arachidonic acid (AA; C20∶4 n-6) and docosahexaenoic acid (DHA; C22∶6 n-3) are important long-chain polyunsaturated fatty acids (LC-PUFA) in maintaining pancreatic beta-cell structure and function. Newborns of gestational diabetic mothers are more susceptible to the development of type 2 diabetes in adulthood. It is not known whether low circulating AA or DHA is involved in perinatally "programming" this susceptibility. This study aimed to assess whether circulating concentrations of AA, DHA and other fatty acids are associated with fetal insulin sensitivity or beta-cell function, and whether low circulating concentrations of AA or DHA are involved in compromised fetal insulin sensitivity in gestational diabetic pregnancies.

Adults presenting with sporadic hypophosphatemia and elevations in circulating fibroblast growth factor-23 (FGF23) concentrations are usually investigated for an acquired disorder of FGF23 excess such as tumor induced osteomalacia (TIO). However, in some cases the underlying tumor is not detected, and such patients may harbor other causes of FGF23 excess. Indeed, coding-region and 3'UTR mutations of phosphate-regulating neutral endopeptidase (PHEX), which encodes a cell-surface protein that regulates circulating FGF23 concentrations, can lead to alterations in phosphate homeostasis, which are not detected until adulthood. Here, we report an adult female who presented with hypophosphatemic osteomalacia and raised serum FGF23 concentrations. The patient and her parents, who were her only first-degree relatives, had no history of rickets. The patient was thus suspected of having TIO. However, no tumor had been identified following extensive localization studies. Mutational analysis of the PHEX coding-region and 3'UTR was undertaken, and this revealed the patient to be heterozygous for a novel germline PHEX mutation (c.2158G>T; p.Ala720Ser). In vitro studies involving the expression of WT and mutant PHEX proteins in HEK293 cells demonstrated the Ala720Ser mutation to impair trafficking of PHEX, with ~20% of the mutant protein being expressed at the cell surface, compared to ~80% cell surface expression for WT PHEX (p<0.05). Thus, our studies have identified a pathogenic PHEX mutation in a sporadic case of adult-onset hypophosphatemic osteomalacia, and these findings highlight a role for PHEX gene analysis in some cases of suspected TIO, particularly when no tumor has been identified.

Osteoporosis is a recognised complication of axial spondyloarthritis (axSpA) and is thought to be due to functional impairment and the osteoclast-activating effects of proinflammatory cytokines. The development of autoantibodies to OPG (OPG-Ab) has been associated with severe osteoporosis and increased bone resorption in rheumatoid arthritis. In this study, we screened for the presence of OPG-Ab in axSpA and reviewed their clinical significance. We studied 134 patients, recruited from two centres in the United Kingdom. Their mean age was 47.5 years and 75% were male. Concentrations of OPG-Ab were related to bone mineral density (BMD) and fracture history using linear and logistic regression models adjusting for age, gender, disease duration and activity, body mass index and bisphosphonate use. We detected OPG-Ab in 11/134 patients (8.2%). Femoral neck and total hip BMD were significantly reduced in OPG-Ab positive patients (0.827 vs. 0.967 g/cm2, p = 0.008 and 0.868 vs. 1.028 g/cm2, p = 0.002, respectively). Regression analysis showed that the presence of OPG-Ab was independently associated with total hip osteopenia (ORadj 24.2; 95% CI 2.57, 228) and history of fractures (ORadj 10.5; 95% CI 2.07, 53.3). OPG-Ab concentration was associated with total hip BMD in g/cm2 (ß = -1.15; 95% CI -0.25, -0.04). There were no associations between OPG-Ab concentration and bone turnover markers, but free sRANKL concentrations were lower in OPG-Ab-positive patients (median 0.04 vs. 0.11 pmol/L, p = 0.050). We conclude that OPG-Ab are associated with hip BMD and fractures in axSpA suggesting that they may contribute to the pathogenesis of bone loss in some patients with this condition.

Hypercalciuria is a major cause of nephrolithiasis, and is a common and complex disorder involving genetic and environmental factors. Identification of genetic factors for monogenic forms of hypercalciuria is hampered by the limited availability of large families, and to facilitate such studies, we screened for hypercalciuria in mice from an N-ethyl-N-nitrosourea mutagenesis programme. We identified a mouse with autosomal dominant hypercalciuria (HCALC1). Linkage studies mapped the Hcalc1 locus to a 11.94 Mb region on chromosome 6 containing the transient receptor potential cation channel, subfamily V, members 5 (Trpv5) and 6 (Trpv6) genes. DNA sequence analysis of coding regions, intron-exon boundaries and promoters of Trpv5 and Trpv6 identified a novel T to C transition in codon 682 of TRPV5, mutating a conserved serine to a proline (S682P). Compared to wild-type littermates, heterozygous (Trpv5(682P/+)) and homozygous (Trpv5(682P/682P)) mutant mice had hypercalciuria, polyuria, hyperphosphaturia and a more acidic urine, and ∼10% of males developed tubulointerstitial nephritis. Trpv5(682P/682P) mice also had normal plasma parathyroid hormone but increased 1,25-dihydroxyvitamin D(3) concentrations without increased bone resorption, consistent with a renal defect for the hypercalciuria. Expression of the S682P mutation in human embryonic kidney cells revealed that TRPV5-S682P-expressing cells had a lower baseline intracellular calcium concentration than wild-type TRPV5-expressing cells, suggesting an altered calcium permeability. Immunohistological studies revealed a selective decrease in TRPV5-expression from the renal distal convoluted tubules of Trpv5(682P/+) and Trpv5(682P/682P) mice consistent with a trafficking defect. In addition, Trpv5(682P/682P) mice had a reduction in renal expression of the intracellular calcium-binding protein, calbindin-D(28K), consistent with a specific defect in TRPV5-mediated renal calcium reabsorption. Thus, our findings indicate that the TRPV5 S682P mutant is functionally significant and study of HCALC1, a novel model for autosomal dominant hypercalciuria, may help further our understanding of renal calcium reabsorption and hypercalciuria.

Multiple endocrine neoplasia type 1 (MEN1) is an autosomal dominant disorder characterized in man by parathyroid, pancreatic, pituitary and adrenal tumours. The MEN1 gene encodes a 610-amino acid protein (menin) which is a tumour suppressor. To investigate the in vivo role of menin, we developed a mouse model, by deleting Men1 exons 1 and 2 and investigated this for MEN1-associated tumours and serum abnormalities. Men1(+/-) mice were viable and fertile, and 220 Men1(+/-) and 94 Men1(+/+) mice were studied between the ages of 3 and 21 months. Survival in Men1(+/-) mice was significantly lower than in Men1(+/+) mice (<68% vs >85%, P<0.01). Men1(+/-) mice developed, by 9 months of age, parathyroid hyperplasia, pancreatic tumours which were mostly insulinomas, by 12 months of age, pituitary tumours which were mostly prolactinomas, and by 15 months parathyroid adenomas and adrenal cortical tumours. Loss of heterozygosity and menin expression was demonstrated in the tumours, consistent with a tumour suppressor role for the Men1 gene. Men1(+/-) mice with parathyroid neoplasms were hypercalcaemic and hypophosphataemic, with inappropriately normal serum parathyroid hormone concentrations. Pancreatic and pituitary tumours expressed chromogranin A (CgA), somatostatin receptor type 2 and vascular endothelial growth factor-A. Serum CgA concentrations in Men1(+/-) mice were not elevated. Adrenocortical tumours, which immunostained for 3-beta-hydroxysteroid dehydrogenase, developed in seven Men1(+/-) mice, but resulted in hypercorticosteronaemia in one out of the four mice that were investigated. Thus, these Men1(+/-) mice are representative of MEN1 in man, and will help in investigating molecular mechanisms and treatments for endocrine tumours.

Paget's disease of bone (PDB) is a common disorder with a strong genetic component characterized by focal increases in bone turnover, which in some cases is caused by mutations in SQSTM1. To identify additional susceptibility genes, we performed a genome-wide association study in 750 individuals with PDB (cases) without SQSTM1 mutations and 1,002 controls and identified three candidate disease loci, which were then replicated in an independent set of 500 cases and 535 controls. The strongest signal was with rs484959 on 1p13 near the CSF1 gene (P = 5.38 x 10(-24)). Significant associations were also observed with rs1561570 on 10p13 within the OPTN gene (P = 6.09 x 10(-13)) and with rs3018362 on 18q21 near the TNFRSF11A gene (P = 5.27 x 10(-13)). These studies provide new insights into the pathogenesis of PDB and identify OPTN, CSF1 and TNFRSF11A as candidate genes for disease susceptibility.

24,25(OH)2D is the product of 25(OH)D catabolism by CYP24A1. The measurement of serum 24,25(OH)2D concentration may serve as an indicator of vitamin D catabolic status and the relative ratio with 25(OH)D can be used to identify patients with inactivating mutations in CYP24A1. We describe a LC-MS/MS method to determine: (1) the relationships between serum 24,25(OH)2D and 25(OH)D; (2) serum reference intervals in healthy individuals; (3) the diagnostic accuracy of 24,25(OH)2D measurement as an indicator for vitamin D status; 4) 24,25(OH)2D cut-off value for clinically significant change between inadequate and sufficient 25(OH)D status. Serum samples of healthy participants (n=1996) from Army recruits and patients (n=294) were analysed. The LC-MS/MS assay satisfied industry standards for method validation. We found a positive, concentration-dependent relationship between serum 24,25(OH)2D and 25(OH)2D concentrations. The 25(OH)D:24,25(OH)2D ratio was significantly higher (P<.001) at 25(OH)D<50 nmol/L. The reference interval for 25(OH)D:24,25(OH)2D ratio in healthy subjects was 7-23. Measurement of serum 24,25(OH)2D can be used as predictor of vitamin D status, a concentration of>4.2 nmol/L was identified as a diagnostic cut-off for 25(OH)D replete status. One patient sample with an elevated 25(OH)D:24,25(OH)2D ratio of 32 and hypercalcaemia who on genetic testing confirmed to have a biallelic mutation of CYP24A1. Our study demonstrated the feasibility of a combined 24,25(OH)2D and 25(OH)D assessment profile. Our established cut-off value for 24,25(OH)2D and ratio reference ranges can be useful to clinicians in the investigation of patients with an impaired calcium/phosphate metabolism and may point towards the existence of CYP24A1 gene abnormalities.

Short periods of basic military training increase the density and size of the tibia, but the adaptive response of bone microarchitecture, a key component of bone strength, is not fully understood.

HPV vaccination efficacy has been shown in clinical trials but it is important to verify population level vaccine effectiveness (VE). We aimed to explore VE and herd effect using HPV infection data from a cohort study of Canadian pregnant women. We analyzed the baseline data of the HERITAGE study, which includes pregnant women recruited in Montreal between 2010-2012 and 2015-2016. Cervicovaginal samples self-collected in the first trimester were tested for 36 HPV types. Vaccination status was self-reported. VE and 95% confidence intervals (CI) were estimated by comparing the prevalence of HPV between vaccinated and unvaccinated women. Herd effect was explored by comparing HPV prevalence in unvaccinated women between the 2 recruitment periods. Adjusted ORs (95%CI) were estimated using exact logistic regression. The proportion of vaccinated women with at least one dose of 4vHPV was 7.5%. Although most of them were vaccinated after the onset of sexual activity, a high VE was found for HPV-16/18 (86.1% (95%CI: 15.0-99.7)). For HPV-6/11/16/18 and for HPV-31/33/45, VE was 61.9% (-23.5-92.6) and 57.0% (-47.7-92.0%), respectively. We also observed a non-statistically significant reduction in the prevalence of HPV-6/11/16/18 and HPV-31/33/45 among unvaccinated women recruited during the second recruitment period (adjusted OR: 0.8 (0.4-1.8) and 0.8 (0.3-1.7), respectively).

Background: Preconception lifestyle interventions appear promising to reduce pregnancy complications, prevent adult cardiometabolic diseases, and prevent childhood obesity. These interventions have almost exclusively been studied in populations of obese infertile women. The development of preconception lifestyle interventions targeting a broader population of overweight and obese women without a history infertility and their partners is needed. Methods: This study is a multicenter open label parallel group randomized controlled trial. Sixty-eight non-infertile women with overweight or obesity in the preconception period and their partners will be recruited from the Sherbrooke and Quebec City regions. The couples will be randomized in a 1:1 ratio to receive the Healthy for my Baby intervention or standard care in the preconception period and pregnancy. Women and their partners will be invited to take part in this lifestyle intervention which includes motivational interviews and daily self-monitoring of lifestyle goals through a mobile phone application. The primary endpoint of this study is the diet quality of women during the preconception period, which will be evaluated using the C-HEI 2007 score at baseline, 2, 4- and 6-months following study enrolment. Women's dietary quality will also be evaluated through the measure of urinary biomarkers of habitual dietary intake at baseline and 2 months in preconception, and 24-26 weeks in pregnancy. Additional indicators of women's lifestyle as well as anthropometric measures will be documented in preconception and pregnancy. For the pregnancy period, the main secondary endpoint is the pattern of gestational weight gain. Pregnancy and neonatal complications will also be evaluated. For partners, diet quality, other lifestyle habits, and anthropometric measures will be documented in the preconception and pregnancy periods. Discussion: This study will evaluate the effectiveness of a low-cost intervention designed to improve diet and other lifestyle characteristics of women in the preconception period who are overweight or obese. If the Healthy for my Baby intervention is efficacious regarding dietary measures, larger trials will be needed to evaluate the impact of this intervention on the rates of pregnancy complications, childhood obesity, and adult cardiometabolic disease. Clinical Trial Registration:clinicaltrials.gov (NCT04242069).

We examined the effects of collagen peptides (CP) supplementation on exercise-induced gastrointestinal (GI) stress.

In the initial stages of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) COVID-19 pandemic, a plethora of new serology tests were developed and introduced to the global market. Many were not evaluated rigorously, and there is a significant lack of concordance in results across methods. To enable meaningful clinical decisions to be made, robustly evaluated, quantitative serology methods are needed. These should be harmonized to a primary reference material, allowing for the comparison of trial data and improved clinical decision making. A comprehensive evaluation of the new Abbott IgG II anti-SARS-CoV-2 IgG method was undertaken using CLSI-based protocols. Two different candidate primary reference materials and verification panels were assessed with a goal to move toward harmonization. The Abbott IgG II method performed well across a wide range of parameters with excellent imprecision (<3.5%) and was linear throughout the positive range (tested to 38,365 AU/ml). The sensitivity (based on ≥14-day post-positive reverse transcription-PCR [RT-PCR] samples) and specificity were 98.3% (90.6% to 100.0%) and 99.5% (97.1% to 100%), respectively. The candidate reference materials showed poor correlation across methods, with mixed responses noted in methods that use the spike protein versus the nucleocapsid proteins as their binding antigen. The Abbott IgG II anti-SARS-CoV-2 measurement appears to be the first linear method potentially capable of monitoring the immune response to natural infection, including from new emerging variants. The candidate reference materials assessed did not generate uniform results across several methods, and further steps are needed to enable the harmonization process.

Military training increases tibial density and size, but it is unknown if men and women adapt similarly to the same arduous training. Seventy-seven men and 57 women not using hormonal contraceptives completed this study. Tibial volumetric bone mineral density (vBMD) and geometry were measured by peripheral quantitative computed tomography (4%, 14%, 38%, and 66% sites) at the start (week 1) and end (week 14) of British Army basic training. Training increased trabecular vBMD (4% site in men; 4% and 14% sites in women), cortical vBMD (38% site), total area (14% and 38% sites), trabecular area (14% site), cortical area and thickness (14%, 38%, and 66% sites), periosteal perimeter (14%, 38%, and 66% sites), and all indices of estimated strength (14%, 38%, and 66% sites); and, decreased endosteal perimeter (66% site) in men and women (all p ≤ 0.045). The increase in trabecular vBMD (4% and 14% sites) was greater in women and the increases in cortical area and strength (38% site) were greater in men (sex × time interactions, all p ≤ 0.047). P1NP increased and βCTX and sclerostin decreased during training in men and women, consistent with adaptive bone formation. PTH decreased in men but increased in women. Arduous weight-bearing activity increased the density and size of the tibia after 14 weeks. Women experienced similar tibial adaptations as men, however, a greater increase in trabecular vBMD in women compared with men could be due to higher loading at this skeletal site in women, whereas the small increase in cortical area could be due to inhibitory effects of oestradiol.

Primary bone cancer (PBC) comprises several subtypes each underpinned by distinctive genetic drivers. This driver diversity produces novel morphological features and clinical behaviour that serendipitously makes PBC an excellent metastasis model. Here, we report that some transfer RNA-derived small RNAs termed tRNA fragments (tRFs) perform as a constitutive tumour suppressor mechanism by blunting a potential pro-metastatic protein-RNA interaction. This mechanism is reduced in PBC progression with a gradual loss of tRNAGlyTCC cleavage into 5' end tRF-GlyTCC when comparing low-grade, intermediate-grade and high-grade patient tumours. We detected recurrent activation of miR-140 leading to upregulated RUNX2 expression in high-grade patient tumours. Both tRF-GlyTCC and RUNX2 share a sequence motif in their 3' ends that matches the YBX1 recognition site known to stabilise pro-metastatic mRNAs. Investigating some aspects of this interaction network, gain- and loss-of-function experiments using small RNA mimics and antisense LNAs, respectively, showed that ectopic tRF-GlyTCC reduced RUNX2 expression and dispersed 3D micromass architecture in vitro. iCLIP sequencing revealed YBX1 physical binding to the 3' UTR of RUNX2. The interaction between YBX1, tRF-GlyTCC and RUNX2 led to the development of the RUNX2 inhibitor CADD522 as a PBC treatment. CADD522 assessment in vitro revealed significant effects on PBC cell behaviour. In xenograft mouse models, CADD522 as a single agent without surgery significantly reduced tumour volume, increased overall and metastasis-free survival and reduced cancer-induced bone disease. Our results provide insight into PBC molecular abnormalities that have led to the identification of new targets and a new therapeutic.

Autosomal dominant hypocalcemia type 1 (ADH1) is caused by germline gain-of-function mutations of the calcium-sensing receptor (CaSR) and may lead to symptomatic hypocalcemia, inappropriately low serum PTH concentrations and hypercalciuria. Negative allosteric CaSR modulators, known as calcilytics, have been shown to normalize the gain-of-function associated with ADH-causing CaSR mutations in vitro and represent a potential targeted therapy for ADH1. However, the effectiveness of calcilytic drugs for the treatment of ADH1-associated hypocalcemia remains to be established. We have investigated NPS 2143, a calcilytic compound, for the treatment of ADH1 by in vitro and in vivo studies involving a mouse model, known as Nuf, which harbors a gain-of-function CaSR mutation, Leu723Gln. Wild-type (Leu723) and Nuf mutant (Gln723) CaSRs were expressed in HEK293 cells, and the effect of NPS 2143 on their intracellular calcium responses was determined by flow cytometry. NPS 2143 was also administered as a single ip bolus to wild-type and Nuf mice and plasma concentrations of calcium and PTH, and urinary calcium excretion measured. In vitro administration of NPS 2143 decreased the intracellular calcium responses of HEK293 cells expressing the mutant Gln723 CaSR in a dose-dependent manner, thereby rectifying the gain-of-function associated with the Nuf mouse CaSR mutation. Intraperitoneal injection of NPS 2143 in Nuf mice led to significant increases in plasma calcium and PTH without elevating urinary calcium excretion. These studies of a mouse model with an activating CaSR mutation demonstrate NPS 2143 to normalize the gain-of-function causing ADH1 and improve the hypocalcemia associated with this disorder.

Pregnant women are an especially important population to monitor for environmental exposures given the vulnerability of the developing fetus. During pregnancy and lactation chemical body burdens may change due to the significant physiological changes that occur. Developmental exposures to some persistent organic pollutants (POPs) have been linked with adverse health outcomes.

Vitamin D status, assessed by serum concentration of 25(OH)D, is the prime candidate marker for many disease-association studies, but the interplay between the subsequent 1,25-dihydroxyvitamin D (1,25(OH)2D) and 24,25-dihydroxyvitamin D (24,25(OH)2D) metabolites is unclear. In this study, we conducted an analysis from a large cohort of healthy, physically fit, young army recruits (n = 940). We found a significant, inverse relationship between serum 25(OH)D and 1,25(OH)2D:24,25(OH)2D vitamin D metabolite ratio (VMR) (r2Exp = 0.582, p < 0.0001), and demonstrated a significant association with increasing PTH concentration (p < 0.001). Circannual rhythms were evident for all vitamin D metabolites and VMRs except for 1,25(OH)2D when fitted to Cosinor curves. We estimated 1,25(OH)2D:24,25(OH)2D VMR of ≥35 to be the threshold value for vitamin D insufficiency, and ≥51 to be predictive of vitamin D deficiency. Our three-dimensional model provides mechanistic insight into the vitamin D-PTH endocrine system, and further substantiates the role of 24,25(OH)2D in human physiology. The model sets a new paradigm for vitamin D treatment strategy, and may help the establishment of vitamin D-adjusted PTH reference intervals. The study was approved by the UK Ministry of Defence research ethics committee (MODREC 165/Gen/10 and 692/MoDREC/15). ClinicalTrials.gov Identifier NCT02416895.