A novel therapeutic strategy for skeletal disorders: Proof of concept of gene therapy for X-linked hypophosphatemia

Description


INTRODUCTION
X-linked hypophosphatemia [XLH; OMIM (Online Mendelian Inheritance in Man) # 307800; incidence, 3.9 of 100,000 births] (1) is the most common form of genetic rickets associated with increased levels of circulating fibroblast growth factor 23 (FGF23). Inactivating mutations in the phosphate-regulating endopeptidase homolog X-linked (PHEX) gene lead to increased FGF23 secretion from osteocytes and osteoblasts, inducing renal phosphate wasting, hypophosphatemia, and impairing the endogenous production of active vitamin D [1,25(OH) 2 vitamin D] via the inhibition of 1-hydroxylase and the activation of 24-hydroxylase (2). Chronic hypophosphatemia and insufficient synthesis of 1,25(OH) 2 vitamin D lead to altered bone mineralization that clinically manifests in children as rickets, i.e., severe skeletal deformities, bone pain, short stature, and severe dental abscesses. In adults, XLH is characterized by osteomalacia, pseudo-fractures, precocious and progressive osteoarthritis and enthesopathies, and poor dental health (1,(3)(4)(5)(6). This considerable disease burden has a considerable impact on the quality of life of affected individuals (7,8).
Standard treatment for XLH consists of multiple daily intakes of oral phosphate, to compensate the renal wasting, and active vitamin D analogs to counteract 1,25(OH) 2 vitamin D deficiency. Although commonly prescribed from infancy to the end of growth, the multiple daily doses of phosphate do not fully correct disease manifestations and are often associated with side effects and poor treatment compliance (1,3,4,9).
Recently, a new treatment based on FGF23 neutralization by a monoclonal antibody (burosumab, Crysvita, Ultragenyx) was approved for patients with XLH (9)(10)(11)(12)(13)(14)(15)(16). Bimonthly subcutaneous injection of burosumab in children with severe XLH resulted in amelioration of rickets with improvement of impaired biochemical parameters and physical ability and reduction in patient-reported pain and functional disability (9,12,13). A 6-to 12-month treatment with burosumab had an impact on the main features of XLH in adult patients by accelerating healing of active fractures and pseudo-fractures and by reducing stiffness and pain (11,15,16). However, the high costs (1,17,18) and the hurdles associated with the need for repeated administrations (1), together with the possible insurgence of antidrug antibodies (19), are potential limitations of the antibody treatment. The availability of predictive preclinical models (20)(21)(22), the extensive clinical data available, and the unmet medical need make XLH an ideal target for the development of innovative strategies to treat skeletal diseases.
Adeno-associated virus (AAV) vectors have been used for in vivo gene replacement in a number of clinical trials, which support the safety and efficacy of this technology (23)(24)(25)(26). Among the different tissues targeted with AAV vectors, the liver represents an ideal candidate due to the accessibility via the bloodstream, the efficient secretory machinery, and the protolerogenic hepatic environment (27). Conversely, AAV gene transfer directed to bone tissue represents an important limitation of this delivery platform. We and others previously demonstrated that liver gene therapy with AAV vectors can be exploited to secrete proteins and mediate cross-correction of enzyme deficiencies at a systemic level (28)(29)(30)(31)(32). On the basis of this, we hypothesized that hepatic gene transfer could be used to modulate signaling pathways in bone tissue and rescue a dominant disease.
Here, through engineering of the FGF23 sequence, we provided proof of concept of the use of the liver to secrete an interfering factor able to reduce the FGF23 pathway overactivation and rescue the bone phenotype in a murine model of XLH carrying a spontaneous deletion of Phex.

Design of AAV-cFGF23 vector and in vitro and in vivo efficacy evaluation
The C-terminal fragment of FGF23 (amino acids 180 to 251 of the human FGF23, cFGF23) competes with the native FGF23 for the binding to the FGF23 receptor/Klotho complex (33). In preclinical models, cFGF23 fusion with the Fc region of immunoglobulin G1 (IgG1) administered as recombinant protein was rapidly eliminated from the circulation and required multiple administrations per week (34). Given the low stability in the circulation of the cFGF23, an ideal method of delivery will require a stable infusion of the truncated factor. Hence, cFGF23 represents an ideal candidate for AAV liver gene therapy, which has a clinically proven capacity to deliver stable levels of transgenes in the circulation (23,26,35). To develop a stable, highly secreted form of cFGF23, we fused the human cFGF23 with the signal peptide of chymotrypsinogen B2 (sp7) (31) and with human albumin fused to cFGF23 either directly or through a cleavable linker derived from human coagulation factor IX (clFIX; Fig. 1A) (36). Those constructs were transfected in human hepatoma cells (Huh-7) in parallel with native human FGF23 and wild-type (WT) and codon-optimized versions of cFGF23 coding DNA fused with the FGF23 native signal peptide or sp7. Western blot analysis indicated that native FGF23 was efficiently produced and secreted in the conditioned medium of transfected Huh-7 cells (Fig. 1B and  fig. S1A). Low levels of the cFGF23 fragment were detected despite codon optimization (co) or fusion with the sp7 signal peptide, possibly due to its poor stability in the extracellular medium (Fig. 1C). Fusion of cFGF23co with sp7 and albumin with or without the clFIX linker greatly improved the secretion in the extracellular medium in vitro (Fig. 1, B and C, and fig. S1A).
cFGF23co fusions with sp7 and albumin with or without the clFIX linker were then tested in vivo and compared to cFGF23co fused with sp7. AAV8 vectors expressing those transgenes under the transcriptional control of the liver-specific apolipoprotein E enhancer, human 1-antitrypsin promoter (hAAT) (31) were injected in 6-week-old C57BL/6J mice at the dose of 1 × 10 12 vector genome (vg) per mouse ( Fig. 1, D to F, and fig. S1, B and C). One month after vector injection, mice were bled to measure the expression of the three cFGF23co variants by Western blot. No bands with a size compatible with cFGF23co were detectable in mice injected with the AAV expressing cFGF23co fused with sp7 ( Fig. 1D and fig. S1B). Bands corresponding to the molecular weight of cFGF23co-albumin chimeric proteins were revealed in mice injected with vectors expressing albumin fusion proteins ( Fig. 1D and fig. S1B). Mice receiving the AAV vector expressing the cFGF23co-albumin fusion containing the clFIX cleavable linker showed significantly higher levels of circulating fusion protein and cleaved FGF23 when compared to its counterpart lacking clFIX (Fig. 1, D and E, and fig. S1B) despite similar transduction as demonstrated by vector genome copy number (VGCN) per diploid genome in liver ( fig. S1C). Only mice injected with the AAV vector expressing the cFGF23co-albumin fusion with the clFIX cleavable linker showed increased phosphate levels compared to control animals, suggestive of decreased FGF23 signaling in the kidneys (Fig. 1F). These data support the hypothesis that the truncated FGF23 is unstable when expressed both in vitro and in vivo. They also suggest that cFGF23 fusion with albumin, despite its increased stability, has no effect on phosphatemia. In conclusion, we demonstrated that the release of cFGF23 from the albumin moiety through a cleavable linker is necessary to achieve stable circulating levels of active cFGF23. Although the exact mechanism behind the increased activity of cFGF23 after release from albumin is unknown, we can speculate that the fusion with the large protein may somehow hinder the interaction of cFGF23 peptide with its receptor.
Liver-targeting AAV-cFGF23 is an efficient delivery platform for the treatment of XLH Deletion of Phex leads to increased circulating FGF23, decreased sodium/phosphate cotransporter (Npt2a) expression in kidney, and the consequent phosphate wasting underlying bone disease in XLH (2,21). To verify whether liver-mediated secretion of the cFGF23albumin fusion was able to rescue phosphate wasting in the Hyp-Duk mouse, a murine model of XLH, an AAV8 vector expressing the cFGF23co-albumin fusion with the clFIX cleavable linker (from now on called AAV-cFGF23) was injected in 4-week-old Hyp-Duk mice at the dose of 1 × 10 12 vg per mouse (Hyp-Duk, AAV-cFGF23). Phosphate-buffered saline (PBS)-injected Hyp-Duk mice (Hyp-Duk, PBS) and WT littermates (WT, PBS) were used as controls ( Fig. 2A). Under these experimental conditions, a single injection of the liver-targeted AAV-cFGF23 resulted in efficient liver targeting ( fig. S1D) and restored the expression of Npt2a in kidney of AAV-injected animals, with levels of expression that were similar to those measured in WT littermates both at mRNA ( protein level (Fig. 2C). Unexpectedly, despite the complete restoration of Npt2a expression in kidney, we observed an only partial rescue of blood phosphate levels at 1 month after injection and a similar tendency at 3 months after injection of the AAV-cFGF23 in Hyp-Duk mice ( Fig. 2D and fig. S1E). We then evaluated whether other Na/Pi transporters were affected in Hyp-Duk mice and rescued by the AAV-cFGF23 treatment. The levels of expression of six Na/Pi transporters were evaluated in kidney. Only type II transporters were modified in Hyp-Duk mice and rescued by the treatment (fig. S1F). Similar to Npt2a, the expression of Npt2c, another type II Na/ Pi transporter, was increased in AAV-cFGF23-treated mice when compared to PBS-treated animals ( fig. S1F). Npt2b, a Na/Pi transporter only residually expressed in kidney (37,38), was increased in Hyp-Duk mice and rescued by the AAV-cFGF23 treatment ( fig. S1F). These results support the hypothesis of a different role for the three type II Na/Pi transporters in the pathogenesis of XLH. Regardless of their role, their levels in kidney are normalized by the treatment.
In XLH, inappropriately normal levels of circulating 1,25(OH) 2 vitamin D were reported. One potential secondary effect of the interference with FGF23, which inhibits the synthesis of 1,25(OH) 2 vitamin D in physiological condition, could be the overactivation of its synthesis. Proof of the potential risk of increasing the synthesis of calcitriol was obtained in Fgf23 knockout mice that have high blood levels of active vitamin D and ectopic calcifications. These calcifications disappeared when these animals were crossed with mice knockout for 1-hydroxylase (39), the main enzyme of 1,25(OH) 2 vitamin D synthesis, thus suggesting that high levels of calcitriol may be associated with the formation of ectopic calcifications. Importantly, the inhibition of the FGF23 signaling through the administration of cFGF23 as a peptide in XLH mice was associated with increased serum phosphate and normal levels of 1,25(OH) 2 vitamin D (34).
To assess the safety profile of our therapeutic strategy, we evaluated the levels of expression of 1-hydroxylase and 24-hydroxylase, key enzymes in 1,25(OH) 2 vitamin D synthesis and inactivation, respectively. Three months after AAV-cFGF23 vector injection, the levels of expression of 1-hydroxylase and 24-hydroxylase were measured in kidney, the primary site of production of these molecules. Both enzymes were increased in Hyp-Duk untreated mice and did not change after 3 months of AAV-cFGF23 treatment ( fig. S2, A and B). In line with these results, blood circulating levels of 1,25(OH) 2 vitamin D were unchanged across treatment groups ( fig. S2C). These findings suggest that the inappropriately normal levels of 1,25(OH) 2 vitamin D in hypophosphatemic Hyp-Duk mice are possibly the consequence of a balance between increased synthesis and degradation of active vitamin D. They also indicate that AAV-cFGF23 has no effect on the vitamin D metabolism in kidney.
Despite unchanged levels of active vitamin D ( fig. S2C), calcium levels in blood measured at 1 and 3 months after injection were significantly decreased in untreated Hyp-Duk mice and were restored in AAV-cFGF23-treated Hyp-Duk mice to levels similar to those AAV-cFGF23 treatment rescues bone microarchitecture, mineralization, and growth plate morphology Three months after vector injection, micro-computed tomography (micro-CT) analysis of the trabecular bone revealed complete restoration of the bone structure (Fig. 3A). In particular, AAV treatment in Hyp-Duk mice completely rescued three important parameters (i) the trabecular structure, namely, the bone volume-to-total volume ratio (Fig. 3B), (ii) the number of trabeculae per millimeter (Fig. 3C), and (iii) the trabecular separation ( fig. S3A), to values similar to those measured in WT animals. Significant rescue of the trabecular pattern in addition to substantial rescue of the trabecular thickness were also indicative of the improvement in bone microarchitecture ( fig. S3, B and C). AAV-cFGF23 treatment also improved the structure of the cortical bone (Fig. 3D) such as the ratio cortical to total cross-sectional area and the cortical thickness (Fig. 3, E and F). Improved bone structure was further confirmed by histological analyses. Von Kossa staining of proximal tibia showed increased bone mineralization in the ossification zone in AAV-treated Hyp-Duk mice (Fig. 3G). The enlarged osteoid observed in untreated Hyp-Duk mice was substituted by mineralized bone after treatment with AAV-cFGF23, with a level similar to that measured in WT animals (Fig. 3, G and H). Consequently, AAV-treated Hyp-Duk mice showed fainter alkaline phosphatase (ALP) staining, a marker of osteogenic cell activity, compared to untreated Hyp-Duk mice, which displayed areas of strong ALP activity at the periphery of the accumulated osteoid (fig. S3, D and E). In addition, tartrate-resistant acid phosphatase (TRAP) staining, a marker of osteoclasts, was observed along mineralized trabeculae in AAV-treated Hyp-Duk mice with a distribution similar to WT mice ( fig. S3F).

AAV-cFGF23 increases bone growth and rescues sacroiliac joint structure
Consistently with the amelioration of bone microarchitecture, we observed substantial improvement of the growth plate structure, responsible for bone elongation, as evidenced by Safranin O staining (Fig. 3I). In untreated Hyp-Duk mice, the columnar structure of chondrocytes in growth plates was severely disorganized, which resulted in the abnormal thickening of the growth plate structure. AAV-cFGF23 treatment restored the columnar organization of the chondrocytes and decreased the total thickness of the growth plate (Fig. 3I). Consistently, blood levels of total ALP, a good marker of skeletal remodeling (40,41), was significantly decreased in AAV-cFGF23-treated Hyp-Duk mice when compared with PBS-treated Hyp-Duk counterparts (Fig. 3J). Together, these results indicate that AAV-cFGF23 treatment in 1-month-old animals normalized bone mineralization and growth plate abnormalities.
We next evaluated the effect of liver gene therapy on the restoration of long bone growth. AAV-cFGF23 injection promoted the elongation of femur and tibia in treated mice together with a general amelioration of the bone structure and a reduction of the distorted epiphysis and diaphysis typically observed in Hyp-Duk mice (Fig. 4, A to D). Phenotypically, AAV-treated Hyp-Duk animals were larger (Fig. 4E) with significantly improved weight overtime (Fig. 4F). Consistent with the elongation observed in long bones (Fig. 4, A to D), we observed a partial rescue of body and tail length in Hyp-Duk mice after AAV-cFGF23 treatment, suggestive of a general positive effect on the longitudinal bone elongation (Fig. 4, G and H).
Osteoarticular rachitic manifestations are typical in adult patients with XLH (7) and were also reported in hypophosphatemic mice (6). We then evaluated the effect of the treatment on sacroiliac arthritis. As expected, untreated Hyp-Duk mice showed multiple erosions and an irregular aspect of the cortical board by micro-CT with higher sacroiliac score compared to WT animals (Fig. 4, I and J,

of 11
and table S1). In AAV-treated Hyp-Duk mice, the treatment rescued signs of sacroiliac arthritis both in terms of morphology and sacroiliac score (Fig. 4, I and J, and table S1). In addition, total body micro-CT scan showed that AAV-treated Hyp-Duk presented with skeleton features similar to WT mice with a restoration of vertebral body hypertrophy (fig. S4, A and B). Consistently with the general improvement in bone architecture, a significant reduction in the number of falls in four-limb hanging test was seen in AAV-treated Hyp-Duk mice 3 months after vector injection to levels similar to those measured in WT animals ( fig. S4C).

DISCUSSION
XLH is due to inactivating mutations in PHEX associated with increased secretion of FGF23 by osteocytes and osteoblasts. Correction of the genetic defect through PHEX replacement directed to those cells may prove difficult regarding their refractory to AAV vectors. Recently, several groups have been focused on the engineering of the AAV capsid to increase bone targeting, resulting also in increased risks to transduce nonbone tissues (42,43). Gene transfer strategies for skeletal diseases remain challenging nowadays.
Here, as an alternative to direct targeting of the bone cells, we propose an innovative strategy to bypass bone targeting by using liver as a biofactory to secrete a negative regulator of the FGF23 signaling pathway and reduce phosphate wasting, the main cause of poor bone mineralization in XLH. Through optimization of the delivery platform, we demonstrated that treatment with AAV-cFGF23 inhibited FGF23 signaling in kidney, resulting in increased expression of the Npt2a transporter and complete rescue of bone mineralization and osteomalacia in a mouse model of XLH. In addition, different from the approved anti-FGF23 monoclonal antibody treatment, in patients with XLH (10), AAV-cFGF23 did not affect blood levels of active vitamin D or the expression of key enzymes involved in its synthesis and degradation and did not trigger calcium deposition in the kidneys of Hyp-Duk mice. Moreover, AAV-cFGF23 completely rescued circulating calcium levels in animals. Future studies will help elucidate the molecular mechanism behind the selective effect of cFGF23 on circulating calcium levels, which, in turn, may constitute an additional safety advantage of gene therapy as a long-term treatment for XLH. Treatment with AAV-cFGF23 in 1-month-old mice resulted in complete rescue of bone structure but partially corrected bone length. This apparent discrepancy is possibly related to the fact that we treated mice at an age where the growth plate, the developmental center for endochondral ossification involved in linear growth, is almost completely mature (44). Although this age was selected because gene transfer to the liver results in stable expression of the transgene (45), we cannot exclude that AAV treatment at a younger age (i.e., 2 weeks) may still provide sufficient transgene expression and complete bone length restoration. Another option may be the repeated administration of AAV vectors starting from the postnatal period. The formation of anti-AAV neutralizing antibodies represents a limitation to this approach. We recently proposed the use of immunomodulation with rapamycin-containing nanoparticles (46) or by the IgG-degrading enzyme IdeS (47) as strategies to enable vector readministration. These findings are changing the paradigm of AAV gene therapy as a one-time treatment and, if successfully translated to the clinic, they may allow for the correction of XLH manifestations in young pediatric subjects.
One potential limitation of our study is the lack of data supporting efficacy in older mice. Future studies will help in evaluating whether the administration of AAV-cFGF23 can have beneficial effects at an advanced stage of the disease. This will provide crucial information on the optimal window of time in which gene therapy can be clinically relevant while generating insights into the disease pathophysiology and reversibility.
Although we cannot exclude that the partial rescue of phosphate levels in blood is due to phosphate reallocation to bone metabolism, one hypothesis is that it could be derived from the increased expression of Npt2b in Hyp-Duk mice kidney. Although the role of this Na/Pi cotransporter in the pathophysiology of XLH is still unknown, further investigations are ongoing to pinpoint its role in FGF23 signaling in kidney and in other organs where Npt2b transporter plays a central role in phosphate homeostasis, for example, the intestine (48).
In conclusion, our results provide proof of concept of efficacy of liver-targeted AAV gene therapy in XLH mice. They also support the use of the liver as a platform to secrete factors able to modulate pathways in tissues refractory to AAV gene therapy, such as bone or kidney, further expanding the therapeutic potential of the approach.

Cloning of the FGF23 constructs
The different FGF23 constructs were cloned into an expression cassette optimized for liver expression, including an enhancer derived from apolipoprotein E gene, hAAT, the human hemoglobin subunit  (HBB2)-derived synthetic intron, and the HBB2 polyadenylation signal (49). Different variants of FGF23 were generated on the basis of the coding sequence of human FGF23. Full-length, native FGF23 (NM_020638) was expressed including its natural signal peptide (nucleotides 1 to 72 of the FGF23 sequence). cFGF23 constructs corresponded to amino acids 180 to 251 of the human FGF23 protein (NP_065689) fused either with the native signal peptide (sp-nat, amino acids 1 to 24 of NP_065689) or with the chymotrypsinogen B2 signal peptide (sp7, amino acids 1 to 18 of NP_00102037) (31). Codon-optimized (co) cFGF23 fused with the sp7 signal peptide was expressed in fusion with human serum albumin without its signal peptide (albumin, amino acids 19 to 609 of NP_000468). A cleavable linker derived from human clFIX (amino acids 182 to 200 of CCA61112) was included in one of the sequences.

AAV vectors production
AAV vectors used in this study were produced using a slight modification of the adenovirus-free transient transfection methods and purified as described earlier. Briefly, adherent human embryonic kidney cells (HEK293) grown in roller bottles were transfected with the three plasmids containing the adenovirus helper proteins, the AAV Rep and Cap genes, and the inverted terminal repeat-flanked transgene expression cassette. After 72 hours of transfection, cells were harvested, lysed by sonication, and treated with Benzonase (Merck-Millipore, Darmstadt, Germany). Vectors were then purified using two successive ultracentrifugation rounds in cesium chloride density gradients. Full capsids were collected; the final product was formulated in sterile PBS containing 0.001% of Pluronic (Sigma-Aldrich, Saint Louis, MO) and stored at −80°C. Titers of AAV vector stocks were determined by SDS-polyacrylamide gel electrophoresis (SDS-PAGE) followed by SYPRO Ruby protein gel stain and band densitometry.

In vitro study
Human hepatocellular carcinoma cells (Huh-7) were transfected by Lipofectamine 2000 (Thermo Fisher Scientific, Waltham, MA) with 2 g of the FGF23 constructs. Two days after transfection, the medium was harvested and the levels of FGF23 were analyzed by Western blot. The experiment was performed three times.

Western blot
Huh-7 medium after transfection and serum from mice injected with the different AAV vectors were loaded onto a 4 to 12% gradient polyacrylamide gel to perform SDS-PAGE. After transfer on nitrocellulose, the membrane was blocked and incubated with an anti-FGF23 antibody (R&D System, Minneapolis, MN). The membrane was incubated with the appropriate secondary antibody (LI-COR Biosciences, Lincoln, NE) and then visualized by Odyssey imaging system and the Image Studio Lite (v5.2) software (LI-COR Biosciences).

In vivo study
Mouse studies were performed according to the French and European legislation on animal care and experimentation (2010/63/EU) and approved by the local Institutional Ethical Board (protocol number 2019-004). Six-week-old C57BL6/J mice were purchased from Charles River Laboratory (Ecully, France) and administered with AAV vectors at 1 × 10 12 vg per mouse intravenously via the tail vein. One month after vector injection, the levels of FGF23 isoforms were measured in circulation by Western blot.
All mice were housed in standard conditions of temperature (23° ± 2°C) in a light-controlled environment, with unlimited access to water and standard pelleted food. Given the small size of the animals, in some cases, mice had gel food supplementation (DietGel Recovery, Clear H 2 O, Portland, ME) from birth until 7 days after injection. DietGel Recovery composition indicated by the supplier is as follows: purified water, corn syrup, vegetable oil, hydrocolloids, wheat protein, electrolyte mix, food acid, and mineral mix. No phenotype differences were observed between groups with and without the supplementation; however, the number of pups surviving among the litter was markedly increased by adding the gel. Mouse genotype was determined by quantitative polymerase chain reaction (qPCR) using primers for Phex gene, on exon 14: forward, 5′-GATCAAATTTTCA GAATCGGACTAC-3′ and reverse, 5′-CCAGAAGAAATCAGACT-GTGC-3′. WT and HypDuk −/0 littermate males were used for this study.
The treatment consisted of intravenous injection of AAV-cFGF23 administered at the dose of 1 × 10 12 vg per mouse at 4 weeks of age. Mice were sacrificed after 3 months when they were at 4 months of age.

Growth parameters
Morphometric parameters, i.e., body weight, total body, and tail length, were measured on mice anesthetized with isoflurane before the treatment and thereafter, every month until sacrifice. Precision scales and a graduated ruler were used for weight and length measurements, respectively. Body length was considered from the muzzle to the beginning of the tail and tail length from the beginning of the tail to its end. Blood parameters Blood sample collection was performed from mice anesthetized with isoflurane before injection and 1 and 3 months later. At the latter times of follow-up, blood parameters were measured: phosphate, calcitriol [1,25(OH) 2 vitamin D], calcium (Ca 2+ ), and ALP.
Phosphatemia and ALP were measured in sera by a colorimetric reaction performed with the FUJI DRI-CHEM NX500 analyzer (FUJIFILM Corporation, Tokyo, Japan). 1,25(OH) 2 vitamin D blood levels were assessed on sera by radioimmunoassay (Immunodiagnostic Systems, Boldon, UK).

Functional evaluation
Four-limb hanging test was performed (after 3 months of treatment) by hanging the mouse on a grid that was then flipped over. The number of falls over a period of 3 min was recorded. The average number of falls per minute was reported for each animal.
Ex vivo study mRNA expression RNA was extracted from 50 to 100 mg of kidney tissue, using TRIzol (Thermo Fisher Scientific) according to the manufacturer's instructions. cDNA was obtained by real-time qPCR (RT-qPCR) performed using RevertAid H Minus First Strand cDNA Synthesis Kit (Thermo Fisher Scientific).

Micro-x-ray CT analysis
All the mice were scanned after sacrifice using a high-resolution x-ray micro-CT device (Quantum FX Caliper, Life Sciences, Perkin Elmer, Waltham, MA) hosted by the PIV (Plateforme Imageries du Vivant) Platform, URP2496, Montrouge, France. Acquisition and analysis were performed following published guidelines (50). Briefly, standard acquisition settings were applied (90 kV, 160 mA), and scans were performed with a field of view alternatively focused on the sacroiliac joint (120 s, 59-m 3 voxel size), the hip (180 s, 59-m 3 voxel size), and the spine (120-s scan time, 40-m 3 voxel size) and covering the full hip and tibia and the full body (26 s, 236-m 3 voxel size). Micro-CT datasets were analyzed using the inbuilt multiplanar reconstruction tool of OsiriX software 5.8 (Pixmeo, Switzerland) to achieve series of images aligned over time for each anatomical region of each animal.
Sagittal images of femur and tibia were used to take the measures of femur and tibia length. Images of full total body, of the dorsal spine, and of the sacroiliac joints were reconstructed. Sacroiliac osteoarthritis was defined as the presence of demineralization of iliac bone, erosions, and irregular board of sacroiliac joint. Erosion of the sacroiliac joints was assessed as described previously (6). The reader was blind to the status of the mouse. A semiquantitative score was established, ranging from 0 (normal) to 3 (most severe feature assessed) for sacroiliac erosions (table S2).
Trabecular and cortical bones were analyzed at distal metaphysis and mid-diaphysis of femur, respectively (180-s scan time, 10-m voxel size). For trabecular analysis, first, image stacks were oriented in the same direction using DataViewer (SkyScan, release 1.5.2.4; Kontich, Belgium) to determine the region of interest (ROI). Second, quantification of the distal metaphysis bone microarchitecture below the growth plate was performed using CT scan Analyser software (SkyScan, release 1.13.5.1). ROI was drawn in the transverse sections along the cortical zone. Batch processing was applied on our data including a median filter (with radius 1) for noise reduction, and adaptive thresholding was performed with a radius of 1. Following segmentation, batch processing calculated three-dimensional (3D) parameters of the binary image (50): bone volume/total volume, trabecular number, trabecular separation, trabecular thickness, and trabecular pattern factor (a marker of trabeculae connectedness). For cortical bone analysis, image stacks were only oriented in the same direction and batch processing was applied on our data, including a median filter (with radius 1) for noise reduction, and adaptive thresholding was performed with a radius of 1. Following segmentation, batch processing calculated 3D parameters of the binary image (50): ratio cortical/total cross-sectional area and cortical thickness. Since the femur length of Hyp-Duk mice is significantly shorter than in WT mice, the ROIs were adjusted according to the femur size (50). The average femur length of WT mouse in our study was equal to 15.3 mm. Two hundred slices in each ROI (distal metaphysis of femur for trabecular bone and mid-diaphysis of femur for cortical bone) taken in analyses correspond to 2 mm, which is equal to 13.1% of total femur length (2 mm/15.3 mm × 100 = 13.1%). Consequently, the ROIs for each femur were adjusted according to this calculation to obtain the 13.1% of total femur length. 3D visualization was performed with software CTVol (SkyScan, release 2.2.3.0).

Murine bone tissue preparation
After sacrifice, bones were fixed overnight at 4°C in 70% ethanol solution and dehydrated in a graded ethanol series. Undecalcified samples were embedded in methyl methacrylate (Merck, Rahway, NJ). Serial sections, 5 m thick, were cut on a microtome (Polycut E microtome, Leica, Wetzlar, Germany). Series of consecutive sections were stained with Von Kossa (5% silver nitrate solution; Sigma-Aldrich), counterstained with toluidine blue (pH 3.8), and stained with Safranin O Lillie's Trichrome (Sigma-Aldrich). Enzyme histochemistry TRAP was used to evaluate osteoclasts by using 2.5 mM naphthol AS-TR phosphate (Sigma-Aldrich), 0.36 M N-N-dimethyl-formamide (Sigma-Aldrich), and 4 mM salt in pH 5.2 acetate buffer. Nonosteoclastic acid phosphatase activity was inhibited with 100 mM l(+)-tartaric acid (Sigma-Aldrich) added to the substrate solution.
ALP was used to reveal the layer of osteogenic cells (pre-osteo/ cementoblasts and osteo/cementoblasts) by incubating the sections with naphthol-ASTR-phosphate (Sigma-Aldrich) and diazonium Fast Blue RR salt (Sigma-Aldrich) for 30 min at 37°C (pH 9) in the presence of MgCl 2 . Immunohistochemistry for Npt2a transporter expression Kidney sections embedded in paraffin were deparaffinized in toluene. After rehydration in a graded ethanol series to pure distilled water, the sections were incubated with antigen retrieval solution (sodium citrate buffer, pH 6) for 10 min at 90°C. Sections were washed in 1% PBS, treated with 1% glycine for 30 min, incubated with blocking solution (1% bovine serum albumin, 0.05% Tween 20, and 10% goat serum) at room temperature for 1 hour and then incubated with anti-Npt2a primary antibody (Slc34A1 antibody, Novus Biological, Centennial, CO) in a humidified dark chamber at 4°C overnight. Sections were washed and then incubated with polyclonal anti-rabbit biotinylated secondary antibody (VECTASTAIN Elite ABC Kit, Vector Laboratories, Burlingame, CA) diluted at 1:200 for 1.5 hours at room temperature in a dark chamber. Sections were washed and treated with 3% H 2 O 2 at 37°C for 30 min then with avidin biotin solution for 60 min (VECTASTAIN Elite ABC Kit), washed in 3% PBS and tris-HCL (pH 7.6), and reacted with diaminobenzidine substrate (Sigma-Aldrich). For nuclear counterstaining, sections were stained with hematoxylin. Control incubations to assess nonspecific staining consisted of the same procedure except that the primary antibody was replaced by nonimmune serum.

Histological quantification
For histological quantification, the images were processed using Fiji (Fiji Is Just ImageJ) software 2.1.0/1.53c, version for Windows. Osteoid quantification was performed on the images of left tibia's secondary spongiosa stained with Von Kossa. The ROI was focused on secondary spongiosa. All the ROIs were adjusted according the tibia size (described in the section of micro-CT analysis). Each ROI was split in square fields (from 7 to 12). The total bone area and the osteoid area were calculated in each square field, and the results were summarized. Then, the percentage of osteoid was calculated dividing the osteoid area to total bone area.

Statistical analysis
Statistical analysis and graphs were performed with GraphPad Prism for Windows, version 7.0. Results were expressed as means ± SD for continuous variables. The comparison of continuous variables was performed using Student's t test or one-way or two-way analysis of variance (ANOVA), as appropriate. P values of less than 0.05 were considered significant.

SUPPLEMENTARY MATERIALS
Supplementary material for this article is available at https://science.org/doi/10.1126/ sciadv.abj5018 View/request a protocol for this paper from Bio-protocol.