Get 20M+ Full-Text Papers For Less Than $1.50/day. Start a 14-Day Trial for You or Your Team.

Learn More →

SNX10 gene mutation leading to osteopetrosis with dysfunctional osteoclasts

SNX10 gene mutation leading to osteopetrosis with dysfunctional osteoclasts www.nature.com/scientificreports OPEN SNX10 gene mutation leading to osteopetrosis with dysfunctional osteoclasts Received: 8 November 2016 1,2 3 2 4 Eva-Lena Stattin , Petra Henning , Joakim Klar , Emma McDermott , Christina Stecksen- 5 6 7 7 8 Accepted: 12 April 2017 Blicks , Per-Erik Sandström , Therese G. Kellgren , Patrik Rydén , Göran Hallmans , Torsten 9 2 4 4 2 Published: xx xx xxxx Lönnerholm , Adam Ameur , Miep H. Helfrich , Fraser P. Coxon , Niklas Dahl , Johan 9 3,10 Wikström & Ulf H. Lerner Autosomal recessive osteopetrosis (ARO) is a heterogeneous disorder, characterized by defective osteoclastic resorption of bone that results in increased bone density. We have studied nine individuals with an intermediate form of ARO, from the county of Västerbotten in Northern Sweden. All afflicted individuals had an onset in early infancy with optic atrophy, and in four patients anemia was present at diagnosis. Tonsillar herniation, foramen magnum stenosis, and severe osteomyelitis of the jaw were common clinical features. Whole exome sequencing, verified by Sanger sequencing, identified a splice site mutation c.212 + 1 G > T in the SNX10 gene encoding sorting nexin 10. Sequence analysis of the SNX10 transcript in patients revealed activation of a cryptic splice site in intron 4 resulting in a frame shift and a premature stop (p.S66Nfs * 15). Haplotype analysis showed that all cases originated from a single mutational event, and the age of the mutation was estimated to be approximately 950 years. Functional analysis of osteoclast progenitors isolated from peripheral blood of patients revealed that stimulation with receptor activator of nuclear factor kappa-B ligand (RANKL) resulted in a robust formation of large, multinucleated osteoclasts which generated sealing zones; however these osteoclasts exhibited defective ruffled borders and were unable to resorb bone in vitro. In healthy individuals, normal bone homeostasis is maintained by balancing bone resorption by osteoclasts with new bone formation by osteoblasts. Osteopetrosis constitutes a heterogeneous group of rare disorders characterized by increased bone mass caused either by an inability to form osteoclasts or by loss of osteoclast 1–4 function . Bone resorbing multinucleated osteoclasts are formed by the differentiation and fusion of mono- nuclear hematopoietic precursor cells in the monocyte lineage. Their differentiation is mediated by macrophage colony-stimulating factor (M-CSF), which is required for proliferation and survival of the mononucleated cells, and receptor activator of nuclear factor kB ligand (RANKL), which triggers differentiation along the osteoclastic lineage. The main clinical features of osteopetrosis are an increased bone density with brittle bones and multiple fractures, narrowing of foramina in the skull with neurological symptoms, and reduced bone marrow cavities 4, 5 leading to impaired hematopoiesis, anemia and recurrent infections . Department of Medical Biosciences, Medical and Clinical Genetics, Umeå University, 901 87, Umeå, Sweden. Department of Immunology, Genetics and Pathology, Science for Life Laboratory, Uppsala University, 751 85, Uppsala, Sweden. Centre for Bone and Arthritis Research, Department of internal medicine and clinical nutrition, Institute of Medicine, Sahlgrenska Academy, University of Gothenburg, 405 30, Gothenburg, Sweden. Arthritis and Musculoskeletal Medicine Programme, Institute of Medical Sciences, University of Aberdeen, Foresterhill, Aberdeen, AB25 2ZD, UK. Pediatric Dentistry, Department of Odontology, Faculty of Medicine, Umeå University, 901 87, 6 7 Umeå, Sweden. Department of Pediatrics, Umeå University, 901 87, Umeå, Sweden. Department of Mathematics and Mathematical Statistics, Computational Life science Cluster (CLiC), Umeå University, 901 87, Umeå, Sweden. 8 9 Department of Biobank Research, Umeå University, 901 87, Umeå, Sweden. Department of Surgical Sciences, Radiology, Uppsala University, 751 85, Uppsala, Sweden. Molecular Periodontology, Department of Odontology, Faculty of Medicine, Umeå University, 901 87, Umeå, Sweden. Eva-Lena Stattin and Petra Henning contributed equally to this work. Correspondence and requests for materials should be addressed to E.-L.S. (email: eva-lena. stattin@igp.uu.se) or P.H. (email: petra.henning@gu.se) Scientific Repo R ts | 7: 3012 | DOI:10.1038/s41598-017-02533-2 1 www.nature.com/scientificreports/ Autosomal recessive osteopetrosis (ARO) comprises the infantile, or “malignant” form, and a less severe inter- mediate form (IARO). In malignant ARO, the skeletal defects have an onset before birth and the resulting nerve damage due to entrapment is irreversible. Treatment with hematopoietic stem cell transplantation (HSCT) is 4, 6 critical in order to prevent disease progression and early death . To date, at least 10 genes coding for osteo- clast function or osteoclast differentiation have been demonstrated to be causally involved in the pathogenesis of 1–4 osteopetrosis . The incidence of ARO is estimated to be 1 in 250,000 births, but can be much higher in isolated 7–9 consanguineous populations . Recently, a missense mutation in the gene encoding sorting nexin 10 (SNX10 [MIM614780]) was shown to be the cause of ARO (OPTB8 [MIM615085]) in consanguineous families of Palestinian origin . Since then, dif- 11, 12 ferent mutations in the SNX10 gene have been reported in osteopetrosis patients of diverse ethnicities . At the initiation of the present study, we were aware of several cases with osteopetrosis of unknown origin in the county of Västerbotten, Sweden. A genealogy study revealed that the patients belonged to the same family. We therefore performed exome sequencing with the aim to find the pathogenic sequence variant causing the disease and, in addition, performed cellular and functional studies to assess the role of the mutation in bone remodeling. During the course of our studies, three patients with osteopetrosis from the Västerbotten County were reported to have a mutation in the SNX10 gene . We confirmed that all osteopetrosis patients in the Västerbotten County carry the same SNX10 mutation. Here, we describe the clinical features and natural course of the disease in the largest cohort of patients with IARO reported thus far. We identified the disease-causing mutation in SNX10, as well as its origin and age and report detailed functional analyses on osteoclasts carrying this pathogenic sequence variant. Results Subjects. We studied nine patients, v fi e males and four females, with IARO from the county of Västerbotten in Northern Sweden. The clinical data are summarized in Table  1. Manifestation of the disease was seen in early infancy in all patients, but the age at correct diagnosis varied from 7 weeks (Pt1) to 3 years (Pt9). The patient (Pt9) diagnosed at the age of 3 years was first thought to have a hereditary optic atrophy. Delayed tooth eruption prompted radiology and the diagnosis of osteopetrosis became apparent (Fig. 1A–D). Skeletal features, radiographs and fractures. Birth length was within ±2 SD (range 48–52 cm), and head circumference ±1 SD (range 34–35 cm). All patients had short final height (males: 132 cm, (−6.7 SD), 160 cm, (−2.5 SD), females: 149 cm, (−2.5 SD), 159 cm, (−1.1 SD) (Table 1)). In all patients, skeletal radiographs show a generalized increase in bone density with metaphyseal modelling defects, transverse bands of greater and lesser density in tubular bones, short wide femoral neck, and”bone-within-bone” pattern of the phalanges (Fig. 1A–D). There were no radiological signs of rickets in five evaluated patients (Pt1, Pt5, and Pt7–9). Seven of nine patients suffered from recurrent fractures up to about 30 times (Pt1), even aer min ft or trauma. During child- hood, almost all fractures were transverse diaphyseal fractures in the extremities without pronounced dislocation and they healed normally with periosteal callus. e Th adult patients (Pt1–5) also had femoral shaft fractures with prolonged healing or pseudoarthrosis despite external or internal fixation. One boy (Pt7), who received hemato- poietic stem cell transplantation (HSCT), did not fracture and his growth normalized (Table 1). Neurological symptoms. All patients developed visual impairment before the age of one year, and five patients became blind (Pt1, Pt4, Pt6, Pt8, and Pt9). Decompression of the optic nerve was not performed in any of the patients. Retinal atrophy was not observed. All but one patient suffers or suffered from hearing impairment, and one patient (Pt6) developed a profound deafness (Table 1). Aer b ft irth, the head circumference increased rapidly in six patients. One girl (Pt9) was found to have hydro- cephalus and cerebellar tonsillar herniation, requiring a ventriculo-peritoneal shunt. She was operated at the age of eight years with decompression of the brainstem, and she died unexpectedly, two years later, possibly due to compression of a postsurgical cephalocele. At that time, she was hospitalized for evaluation of severe back pain, and headache episodes. Two patients (Pt8, and Pt5) developed stenosis of the foramen magnum; one (Pt8) died at the age of 12 years due to spinal cord compression while the other (Pt5) is still alive at the age of 24 years. One boy (Pt6) showed intellectual disability in addition to early blindness and profound deafness. Facial paralysis was present in four patients (Pt1–4). Evaluation of six magnetic resonance imaging (MRI), and three computed tomography (CT) scans showed a thickened and sclerotic bone of the skull in all but one patient (Table 2). Scans did not show brain malforma- tions (Fig. 1E,F). Hematological symptoms. Anemia was present early in life in four patients (Pt1, Pt2, Pt7, and Pt9) due to bone marrow encroachment. Consumption and hemolysis of erythrocytes, followed by a low leucocyte count and throm- bocytopenia, was seen in two patients (Pt8, and Pt9); haptoglobin and carboxyhemoglobin (CoHb) were below the detection limit. Splenectomies were performed on these patients and the hematopoiesis almost normalized aer ft wards. Six of nine patients require or required repeated blood transfusions due to normocytic anemia (Table 1). All patients were oer ff ed stem cell transplantation, most of them during adolescence, but HSCT was accepted only by the parents of the youngest one (Pt7). He received HSCT at the age of 33 months from an HLA-matched unrelated donor; he is still alive and healthy at the age of 15 years. However, he suffers from irreversible tissue damage that occurred prior to transplantation: he has visual deterioration as well as missing and malformed teeth. Two patients (Pt4, and Pt6) died due to a combination of pancytopenia and septicemia; one (Pt6) died at the age of 19, due to septicemia originating from osteomyelitis of the jaw, and the other (Pt4) died at the age of 35, due to chronic osteomyelitis of one femur aer a ft n unhealed fracture. Autopsy was not performed in any of the patients (Table  1). Scientific Repo R ts | 7: 3012 | DOI:10.1038/s41598-017-02533-2 2 www.nature.com/scientificreports/ Patient number 1 2 3 4 5 6 7 8 9 Gender Male Female Male Female Female Male Male Male Female age (year) 47 43 46 35 24 19 15 12 10 51 cm 50 cm 52 cm NA 48 cm 52 cm NA (±0SD) NA 48 cm Birth length, 4000 g 3620 g 3910 g 3110 g 2930 g 3500 g 3094 g NA 3305 g weight, OFC 34 cm NA NA NA NA 34.5 cm NA NA 35 cm Normal growth, 132 cm 149 cm 160 cm 159 cm NA NA — — Final height aer ft HSCT (age year) (SD) (17)(−6.7) (20)(−2.5) (20)(−2.9) (19)(−1.1) Age at onset Early infancy Early infancy Early infancy Early infancy Early infancy Early infancy Early infancy Early infancy Early infancy 36 months. 5 months. 28 months. Retained 6 months. 19 months. Optic atrophy. teeth. Nasal 7 weeks. Optic atrophy. Optic atrophy. Age at 2 years. Optic 12 months. Days, disease Macrocephaly. Macrocephaly Optic atrophy. Optic atrophy. Macrocephaly congestion. diagnosis, atrophy. Optic atrophy, known in the 46 cm. Macrocephaly Macrocephaly 55.5 cm. Optic atrophy. Anaemia ( Hb symptoms Macrocephaly blindness family 6.64 mmol/l) Anaemia (Hb 44.5 cm 50 cm Anaemia (Hb Macrocephaly. 6.09 mmol/l) 6.33 mmol/l) Anaemia (Hb 6.6 mmol/l) Stem cell HSCT at the age — — — — — — — — transplantation of 33 month Deceased (age) — — — 35 years — 19 years — 12 years 10 years Blindness. Hydrocephalus, Blindness. Impaired vision VP-shunt. Impaired vision, Impaired Impaired vision, and hearing. Blindness. Decompression Blindness. hearing. Blindness. Impaired (blind on one (blind on one Foramen Profound op of the Neurological hearing, unilateral eye). Impaired Impaired Died due eye). Impaired magnum deafness. Impaired vision brainstem at hearing. Facial to foramen symptoms profound deafness. hearing, tuba hearing. Facial stenosis, op Intellectual the age of 8, Facial paralysis aperta. Facial paralysis magnum paralysis decompression disability died due to paralysis stenosis at the at the age of 14 compression of age of 12 a postsurgical cephalocele Anaemia, hepato- Anaemia, Anaemia, splenomegaly. Anaemia, hepato- Anaemia, Anaemia, hepato- hepato- Blood splenomegaly, hepato- hepato- Anaemia, splenomegaly, splenomegaly, Haematological transfusions Almost normal Normal splenectomy at splenomegaly. splenomegaly. hepato- splenectomy splenectomy symptoms Died due to haematopoiesis. haematopoiesis. the age of 5. Blood Blood Blood splenomegaly. at the age at the age septicaemia transfusions. transfusions. transfusions. of 9. Blood of 5. Blood originating transfusions. transfusions. from osteomyelitis. At the age of 28 months, two At the age of Eruption of At the age of At the age of At the age of At the age of primary teeth 31 months, primary teeth NA   4-6 months, ten 6 months, five NA   6 months, 12 NA 13 months, six Malformed and eight primary (month), Malformed and primary teeth primary teeth Malformed and primary teeth Malformed primary teeth retained but a teeth and teeth retained Malformed and Malformed and retained Malformed and and retained Malformed lot better than Malformed and development retained retained retained and retained patients not retained transplanted At the age of 6. At the age of Osteonecrosis of the At the age of 13. At the age of 30. 12. Died due ONJ (year) ear canal, and the Osteonecrosis Osteonecrosis ONJ (age NA) At the age of 18. to septicaemia None None None skull with associated of the ear canal of the ear canal originating epidural abscess from ONJ Narrow nasal, Nasal Early nasal Narrow nasal, Narrow nasal, and and pharynx obstruction, obstruction, and pharynx pharynx cavities Narrow nasal, Narrow nasal, Narrow nasal, cavities before the age Airway before the age cavities CPAP Tracheostomy at the and pharynx and pharynx and pharynx None Apnoea during of 12 month, symptoms of 5 months. during night age of 22. CPAP during cavities cavities cavities the night adenoidectomy CPAP during since the age night Tracheostomy Tracheostomy night of 35. at the age of 9. at the age of 8. Multiple Multiple Multiple Multiple fractures, Multiple Multiple Multiple Bone fractures fractures, fractures, fractures, None Few fractures pseudoarthrosis fractures fractures fractures pseudoarthrosis pseudoarthrosis pseudoarthrosis Exophthalmos keratitis, Exophthalmos endophthalmitis, Exophthalmos, Exophthalmos, Exophthalmos Drooling, Exophthalmos enucleation of one Bilat pes Mitral valve Inguinal Exophthalmos, Exophthalmos, treated Media otitis, Other features eye. Ear canal stenosis, calcaneovalgus, insufficiency, hernia, Exophthalmos Drooling, Severe back with Botox, Wheelchair Kidney tumour, Hypertension, Wheelchair Wheelchair Anosmia pain Wheelchair bound Urinary stones, Weight Walking frame bound bound bound loss, Wheelchair bound Table 1. Clinical features of nine individuals with a splice site mutation in the SNX10 gene, and intermediate A B autosomal recessive osteopetrosis (IARO). OFC = occipital-frontal circumference. NA = Not available. C D E SD = Standard deviations. HSCT = hematopoietic stem cell transplantation. Symptoms that lead to F G H the diagnosis. Hb = hemoglobin. VP = ventriculo-peritoneal shunt. ONJ = osteonecrosis of the jaw. CPAP = continuous positive airway pressure. Scientific Repo R ts | 7: 3012 | DOI:10.1038/s41598-017-02533-2 3 www.nature.com/scientificreports/ Figure 1. A Radiographs of subjects with the Västerbotten form of intermediate osteopetrosis, and sequencing of the SNX10 gene. (A–D) Skeletal radiographs of a girl at the age of three years (Pt9), showing characteristic features of intermediate osteopetrosis. (A) Skull: Increased bone density of the skull base, frontal bossing and VP-shunt (B) AP view of the legs: “Erlenmeyer flask” shape of the femoral diaphysis. Longitudinal and transverse bands of lesser density in the metaphyses. Bone-within-bone pattern in the epiphyses in the knee. (C) Lateral view of the spine: Very dens endplates of the vertebral bodies, so called “sandwich vertebrae”. The ribs (seen in part) are broad and dens. (D) AP view of the left arm: No fractures or sign of rickets, “bone- within-bone” pattern of the phalanges. (E,F) Magnetic resonance imaging of the brain, and spine in an 11-year-old boy. MRI shows cerebellar tonsillar descent through the foramen magnum, a foramen magnum narrowing, and signs of brain stem compression with resulting syringo-hydromyelia. (G,H) Consequence of the c.212 + 1 G > T variant on (G) genomic level and (H) transcript level. (G) Sequencing of DNA from a control (top), a heterozygous carrier (middle) and a homozygous patient (bottom) for the c.212 + 1 G > T variant. The variant changes the donor splice site at the 5′ end of intron 4 from GT to TT. This leads to a predicted use of an alternative donor splice site 16 nucleotides upstream of intron 4. (H) Sequencing of cDNA from a control (top) and a patient (bottom) confirms the use of an alternative splice site and the skipping of 16 bp leading to a frame shift and a premature stop codon (p.S66Nfs * 15) at the protein level. Symptoms of the airways, jaw and teeth. All patients had narrow nasal, and pharyngeal cavities, and two patients (Pt2, and Pt9) had documented nasal congestion before the age of one year. Three patients (Pt1, Pt8, and Pt9) needed a tracheostomy at the age of 8, 9 and 22 years, respectively. Three adult patients (Pt1–3) are treated with continuous positive airway pressure (CPAP) during nighttime, due to sleep apnea syndrome. Tooth Scientific Repo R ts | 7: 3012 | DOI:10.1038/s41598-017-02533-2 4 www.nature.com/scientificreports/ Patient number 1 1 1 4 5 7 8 8 9 Gender (male/female) Male Male Male Female Female Male Male Male Female Age at the time of imaging (year) 29 34 33 35 14 2 11 11 8 MRI/CT-scan MRI MRI CT CT MRI MRI MRI CT MRI Corpus callosum aplasia/hypoplasia 0 0 X X 0 0 0 X 0 Sclerotic bone 1 1 1 1 1 0 1 1 1 Thickened bone 1 1 1 1 1 0 1 1 1 Ventriculomegaly 0 0 0 0 0 0 1 1 0 Remodelled inner table 1 1 1 0 1 0 0 1 0 Tonsillar herniation 1 1 X 0 1 0 1 X 1 Optic nerve sheath dilatation 0 0 0 0 1 0 1 1 0 Proptosis 1 1 1 1 0 0 1 X 1 Brain atrophy 0 0 0 0 0 0 0 0 0 Optic nerve atrophy 0 0 0 X X 1 1 X 1 Optic canal stenosis 0 1 1 1 0 0 1 0 0 Ear fluid 1 1 X X 0 1 0 X 0 Dural venous sinus stenosis X X X 0 1 X 1 X 1 Foramen magnum stenosis 0 0 X X 1 0 1 X 1 Posteriorly angulated dens with stenosis X X X X 1 0 1 X 1 Cervical spinal stenosis in included upper 1 1 0 0 0 0 0 0 0 part Subdural hematoma or intra-cerebral 0 0 0 0 0 0 0 0 0 bleeding Prominent collateral veins between sinus 0 0 0 0 0 0 1 0 0 transversus and extra-cranial veins Lemon shape of the scull (similar to what is seen in Chari II malformation), syrinx 0 0 0 0 0 0 1 0 0 of the cervical cord, large occipital horns Calcification in basal ganglia and 0 0 0 0 0 0 0 1 0 subcortical frontal lobes Table 2. Magnetic resonance imaging (MRI) and computed tomography scan (CT) of the skull and brain, in six individuals with SNX10 related intermediate autosomal recessive osteopetrosis (IARO). 1 = present, 0 = not present, X = feature not possible to evaluate. eruption was delayed in three patients (Pt7–9), and teeth are congenitally missing, malformed or retained in all patients. Severe osteomyelitis and osteonecrosis of the jaw is seen in all adult patients (Pt1–6), with abscess for- mations and fistulas to the orbits and skin. One patient (Pt1) presented with osteonecrosis of the jaw at the age of six years, other patients (Pt2, Pt3, Pt5, and Pt6) at the age of 12, 13, 18 and 30 years, respectively. Osteonecrosis of the ear canal is seen in three of the eldest patients (Pt1–3) (Table 1). Exome sequencing, data filtering, DNA and cDNA sequencing, carrier frequency. Exome sequencing of three patients (Pt1, Pt2, and Pt5) revealed 115,283 unique SNPs. Among those, 25,301 were com- mon, 24,867 remained aer ex ft cluding SNPs on chromosomes X and Y, and 9,401 were homozygous. Nineteen of those were not present in the control population and three were non-synonymous, and located in three different genes. One of the three genes was the SNX10 gene. Homozygous sequence variants at the donor splice site of exon 4 (c.212 + 1 G > T) were confirmed in all ae ff cted patients using Sanger sequencing, and the parents were found to be heterozygous for the splice site mutation (Fig. 1G–H). The splice site variant is not present in any of >1000 13, 14 exomes present in our local Canvas DB, or in controls in Exome Seq. project or 1000 Genomes . Sequence analysis of cDNA from cultured patient peripheral blood cells revealed that an alternative splice site, located 16 bases upstream from intron 4, was used, resulting in a frame shift and a stop codon (p.S66Nfs * 15). The carrier frequency of the sequence variant in the county of Västerbotten is 1:93 or 1.1% of the population. Genealogy and haplotype analysis. Genealogic data revealed a common ancestor 8–12 generations th back, in the early 19 century. He was born 1820 in Bureå, Skellefteå, County of Västerbotten, northern Sweden. Haplotype analysis was performed to evaluate the age of the sequence variant. The size of the homozygous regions in Pts1–5 is 4044451 bp, 3494094 bp, 3515023 bp, 5632740 bp, 4628862 bp, and 4628862 bp, respectively. Pt9 showed homozygosity between marker rs477644 and rs6462058, corresponding to 11307603 bp, which was included in the age calculations (Fig. S1). e Th linkage disequilibrium in six ARO probands extends over a mean homozygous region of 5.3 cM, which equals 38 meioses/generations. With an assumed generation time of 25 years, this corresponds to a shared ancestral origin of the SNX10 gene variant 950 years ago . Monocytes from patients form large osteoclasts in vitro. We isolated CD14 monocytes from the peripheral blood of three patients, (Pt1, Pt2 and Pt5) and seven healthy controls, designated C1–C7. We assessed Scientific Repo R ts | 7: 3012 | DOI:10.1038/s41598-017-02533-2 5 www.nature.com/scientificreports/ Figure 2. Osteoclast differentiation of CD14 monocytes from peripheral blood. Photographs of TRAP staining (A) and quantification of surface area (B ) of CD14 cells cultured in M-CSF (M) or M-CSF and RANKL (M/RL) for four days (C1, C2, Pt2, Pt5) or five days (C1, lower photos, and Pt1). Surfaces are presented as mean ± SEM. (C) Actin ring staining of cells cultured in M/RL for four days. Scale bars: 100 μm. *** P ≤ 0.001. osteoclastogenesis on both plastic and bone slices, by incubating CD14 monocytes with RANKL. Cells from patients were tested alongside control cells in six different independent experiments. The data are presented sep- arately for these individual experiments instead of being pooled because of the large inter-individual variability that is commonly seen in primary osteoclast cultures from humans. Multinucleated cells that stained positive for tartrate resistant acid phosphatase (TRAP ) were formed on + + plastic dishes in all cultures from patients and controls. The TRAP multinucleated cells (TRAP MuOCL) from patients were larger and stained less intensely for TRAP than those from controls (Fig. 2A,B). Accurate counting of TRAP MuOCL in patient cultures was difficult as osteoclasts in patient cultures were large and pale. The patient TRAP MuOCL formed peripheral belts of F-actin-rich podosomes similar to controls. The cir - cumference of these actin belts was greater in the patient osteoclasts than the control osteoclasts, reflecting the greater size of the patient cells (Fig. 2C). Expression of osteoclastic and osteoclastogenic genes in cells from patients. We next assessed if the expression of genes associated with osteoclast differentiation and function was altered in the patient oste- oclasts. Gene expression analyses were performed twice with cells from two patients (Pt2 and Pt5) at one time point during the culture, and with two different controls (C1–C4) in each analysis (Figs  3 and S2). The reason data are provided separately for different controls and different patients is that the logistics of the cell collections and experiments did not allow us to compare the same control with the same patient in each experiment. We Scientific Repo R ts | 7: 3012 | DOI:10.1038/s41598-017-02533-2 6 www.nature.com/scientificreports/ Figure 3. Gene expression in osteoclasts from patients (Pt2, Pt5), and controls (C1, C2, C3, C4). Fold induction of gene expression between cells treated with M-CSF (M) with and without RANKL (RL) for three days. All values are given as mean ± SEM (n = 4). *P ≤ 0.05, **P ≤ 0.005, ***P ≤ 0.001 between indicated groups. performed separate experiments comparing cells from different controls and from different patients. We repeated the experiments at different times to ensure reproducibility of our observations. e mRN Th A expression of ACP5 (encoding the enzyme TRAP) and CTSK (encoding cathepsin K), both robust markers of the osteoclastic phenotype, was strongly upregulated by RANKL in cells from both patients and con- trols (Fig. 3A,B). There was no difference between patients and controls in the degree of upregulation of NFATc1 (encoding nuclear factor of activated T-cells, the master osteoclastogenic transcription factor), (Fig. 3C). No die ff rences between patients and controls were found in RANKL-induced upregulation of TCIRG1 (encoding the proton pump subunit, Atp6i, essential in acidification process) (Fig.  3D), or the mRNA expression of ITGAV and ITGB3, (encoding the α and β subunits of the vitronectin receptor, essential for osteoclast attachment to v 3 bone) (Fig. 3E,F). The expression of CLCN7 (encoding chloride channel 7, with an essential role in acidification) was consistent throughout osteoclastogenesis in cells from both patients and controls (Fig. 3G). Although some genes showed a significant difference between patient and controls in one of the two independent experiments, none of the analyzed genes were significantly and reproducibly changed in both patients compared to controls. es Th e data indicate that the SNX10 mutation does not ae ff ct the expression of key genes associated with osteoclast differentiation and function. e Th ability of osteoclasts with the SNX10 mutation to resorb bone was analysed by culturing CD14 cells from patients on slices of devitalized bovine bone in six different experiments, each comparing cells from one of 3 patients with several of 7 controls. In control cultures, RANKL stimulated the formation of many TRAP MuOCL. In patient cultures we observed TRAP MuOCL, but, as in cultures on plastic, we also noticed several much larger cells with weak TRAP staining, making accurate counting of the number of osteoclasts impossible (Fig. 4A). We, therefore, also used TRAP5b released into the culture media as a parameter of mature osteoclast formation , but found no consistent difference between patients and controls (data not shown). Phalloidin staining demonstrated that control osteoclasts frequently exhibited crescent-shaped actin rings, which were associated with extensive resorption pits and resorption trails in the bone slices (Fig. 4B,C) and release of CTX (a degradation fragment from collagen type I) (Fig. 4D). Patient osteoclasts polarized their actin Scientific Repo R ts | 7: 3012 | DOI:10.1038/s41598-017-02533-2 7 www.nature.com/scientificreports/ Figure 4. Patient-derived osteoclasts form actin rings but are unable to resorb bone. (A) TRAP staining of CD14 cells cultured in M/RL on bone for eight days. Scale bar 40× panel: 500 μm, scale bar 100× light microscope and reflective light panel: 200 μm, scale bar 100× toluidine blue panel: 100 μm. (B,C) Osteoclasts were analyzed by confocal microscopy; (B) stained with TRITC-phalloidin (red), DAPI (blue) and the bone surface using a fluorescent bisphosphonate (AF-ALN; green). Resorption pits stain brightly with AF-ALN and can be seen in the control only. Scale bar = 100 μm; (C) stained with TRITC-phalloidin (red), DAPI (blue) AF- ALN (magenta) and the cell membranes with wheat germ agglutinin (green). AF-ALN omitted from upper panel of Pt2 for clarity. Black patches in C2 xy image indicates resorption pits. Bottom panel of Pt2 is a zoom of the zx section illustrating lack of resorption. Scale bar = 20 μm; (D) CTX-1 release in media during day 6 to 8 of culture in M/RL. C1, C6 = controls, Pt2, Pt5 = patients (n = 4–5). Scientific Repo R ts | 7: 3012 | DOI:10.1038/s41598-017-02533-2 8 www.nature.com/scientificreports/ cytoskeleton in rings at the periphery of the cells, similar to the actin belts seen in cells on plastic. Patient cultures, however, showed no evidence of resorption. er Th e were no resorption pits (Fig.  4A–C) and there was no release of CTX (Fig. 4D). There was no evidence of vacuolation in any of the patient osteoclasts (Fig.  4B). Increased size of patient osteoclasts is independent of precursor fusion. Individual patient oste- oclasts occupied larger areas of the bone surface than control osteoclasts (Fig. 4A–C), similar to the observations on plastic (Fig. 2A–C). However, osteoclasts on bone from patients contained similar numbers of nuclei as those from healthy controls (mean nuclei/osteoclast = 7.3 for C1, 7.8 for C2, 7.9 for Pt2). Patient osteoclasts were much flatter than the control osteoclasts; in cross-sectional views most osteoclasts were no deeper than a single nucleus (Fig. 4C). This indicates that the larger size of the patient osteoclasts is the result of increased spreading, rather than the result of increased fusion. In agreement with these findings, we were unable to observe any significant differences replicated in two separate experiments between expressions of several genes suggested to be involved in osteoclast fusion (Fig. S2). 17–19 RANKL-induced upregulation of the osteoclast fusion genes DC-STAMP and OC-STAMP was similar in cells from patients and controls (Fig. S2). Syncytin-1 and its receptor Asct2 have also been suggested to be involved in osteoclast fusion but the mRNA expression of the genes encoding these proteins (ERVW-1 and SLC1A5, respec- tively) was not different between patients and controls (Fig.  S2). Osteoclasts exhibiting SNX10 mutation have an impaired ruffled border. Since osteoclast forma- tion and F-actin ring formation was unae ff cted, we speculated that the lack of resorption by osteoclasts derived from the patients would be due to defective ruffled border formation, similar to well-characterized cases of osteoclast-rich ARO . This would be consistent with the likely role of SNX10 in vesicular trafficking processes 22, 23 in osteoclasts. While control osteoclasts on bone accumulated PNA-lectin (a marker of the ruffled border) within the actin ring, patient osteoclasts showed no PNA-lectin staining (Fig. 5A), suggesting impaired ruffled border formation in the patient osteoclasts. Transmission electron microscopy was used to further examine the sealing zone (F-actin rings) and the ruf- fled border in cultured osteoclasts. This ultrastructural analysis confirmed that patient osteoclasts retained the capacity to form sealing zones (Fig. 5E), but that 65% of patient osteoclasts lacked ruffled borders (Grade 0) and showed no evidence of bone resorption. In contrast, control osteoclasts showed evidence of extensive bone resorption and only 18% of these cells lacked ruffled border (Fig.  5B,D and F). Fifteen percent of patient osteo- clasts showed moderate/extensive ruffled borders, compared to 50% of control osteoclasts. The ultrastructural analysis confirmed that patient osteoclasts were very large and flat compared to control cells, with sealing zones located at the cell periphery. Gene expression of wild type and mutated SNX10. SNX10 mRNA levels were measured in RANKL-stimulated CD14 cells using two different TaqMan assays, spanning the junction between exon 1–2 and 4–5, respectively. SNX10 mRNA, analyzed with the exon 1–2 assay, was present in patients and controls and significantly upregulated by RANKL in controls and in one of the two patients (Fig.  6A). However, the mRNA levels were lower in patient cells. The assay spanning the exon 4–5 junction also showed upregulation of SNX10 by RANKL in control cells (Fig. 6B). However, as the mutation is located at the splice donor site of exon 4, this domain was barely detected in cells from patients. Discussion We herein report the clinical features and natural course of the Västerbotten form of IARO [MIM 615085], the disease-causing mutation, its origin and age. We also present functional analyses of the effects of this mutation on osteoclast differentiation and function. The Västerbotten form of IARO is caused by a splice site mutation in the SNX10 gene [MIM 614780] c.212 + 1 G > T resulting in a frame shift and a stop codon (p.S66Nfs * 15). The carrier frequency was found to be 1/93 in the population of Västerbotten. According to haplotype analysis, the shared ancestral origin of the SNX10 gene variant is approximately 950 years old. 1–3, 21, 24 e Th phenotype of IARO, previously associated with mutations in the PLEKHM1 , CLCN7, and CA2 genes , is defined as asymptomatic at birth, with patients exhibiting fractures by the end of their first decade of life, 24–27 while bone marrow failure and hepato-splenomegaly are rare . In contrast, our data show that patients with SNX10-related IARO of Västerbotten type described here have a more severe phenotype, with age at onset in early infancy, presenting with optic atrophy, macrocephaly, frontal bossing, proptosis, and anemia. Splenomegaly, due to compensatory extra-medullary hematopoiesis, was present before the age of five, and splenectomy was performed in three patients (Pt1, Pt8, and Pt9), the youngest at the age of five years (Pt1). Although the patients displayed many of the symptoms seen in malignant ARO, they all survived early childhood without being trans- planted. We, therefore, placed the patients in the group of IARO, rather than in the group of malignant ARO in which patients do not survive childhood if not transplanted. Osteopetrosis presenting with hydrocephalus, hind-brain posterior fossa crowding, foramen magnum nar- 28–31 rowing, and syringohydromyelia is rare, but has been described in the literature . Cure et al. studied cranial MRI of osteopetrosis, but did not identify ventriculomegaly, tonsillar herniation, or remodelled inner table in IARO . In the Västerbotten intermediate form of osteopetrosis these features seem to be common; in our study, four patients (Pt1, Pt5, Pt8, and Pt9) had cerebellar tonsillar descent through the foramen magnum, and three (Pt1, Pt5, and Pt8) of them showed remodelled inner table. Three subjects (Pt5, Pt8, and Pt9) had foramen mag- num stenosis, posteriorly angulated dens with stenosis and signs of brain stem compression. Structural brain malformation is rarely reported as part of the clinical spectrum of osteopetrosis, but cerebral atrophy, corpus callosum agenesis, hypoplasia of the hippocampus, and Dandy-Walker malformation has been Scientific Repo R ts | 7: 3012 | DOI:10.1038/s41598-017-02533-2 9 www.nature.com/scientificreports/ Figure 5. Patient-derived osteoclasts formed defective ruffled borders. (A) Phalloidin staining for acting rings (left) and lectin staining for ruffled border (middle) of CD14 cells cultured in M/RL on bone for 11 days. (B–E)-Representative TEM micrographs of control C2 (B,C) and patient Pt2 (D,E) -derived osteoclasts show extensive ruffled border (grade 3) and resorption in the control, compared to lack of ruffled border (grade 0) and resorption in the patient. Patient osteoclasts formed sealing zones (E). N = nucleus, SZ = sealing zone, RB = ruffled border. Scale bars in left panels = 2 μm. Scale bars in right panels = 500 nm. (F) The percentage of osteoclasts categorized into each ruffled border grade (control, n = 38; patient, n = 34). 33–36 reported . However, none of our SNX10 related IARO patients had intra-cerebral malformations. MRI of the brain showed calcification of the basal ganglia, and subcortical frontal lobes in one male patient (Pt8), who died due to spinal cord compression caused by a combination of prominent calvarial and skull base thickening, Scientific Repo R ts | 7: 3012 | DOI:10.1038/s41598-017-02533-2 10 www.nature.com/scientificreports/ Figure 6. Expression of SNX10 mRNA in cells treated with M or M/RL for three days, analyzed with an assay covering the exon 1–2 (A) and exon 4–5 (B) junction. Values are given as mean ± SD (n = 2–4). ***P ≤ 0.001 between M and M/RL treated cells. toncillar herniation, foramen magnum narrowing and posteriorly angulated dens with stenosis. He also had ven- triculomegaly, dural venous sinus stenosis, and remodelled inner table as a sign of elevated intra-cranial pressure. MRI also showed collateral veins between sinus transversus and extra-cranial veins, possibly caused by a jugular foramen stenosis, occlusion of the jugular vein and thrombosis of the transverse and sagittal sinuses, which is not previously described in subjects with IARO. However, no brain phenotype was reported in Snx10 knock out mice . Calcification of basal ganglia, dura mater, tentorium cerebelli and periventricular calcification extending into the white matter (Marble brain disorder) has previously been described with or without association to renal 38–41 tubular acidosis due to carbonic anhydrase II deficiency . The reason for calcification of the brain is unclear. Results of studies performed in Snx10 knock out mice suggest an additional role for Snx10 in the gastric epithelium. Snx10 knock out mice have impaired gastric acidification and are severely hypocalcaemic com- pared to wild type littermates, resulting in rickets as well as osteopetrosis (osteopetrorickets). In contrast, osteoclast-specic fi knockout of Snx10 resulted in osteopetrosis without rickets or alterations in calcium balance . Interestingly, the patients in the present study have no signs of rickets. Also, SNX10 osteopetrotic patients have been shown to be almost completely rescued from the sclerotic phenotype by hematopoietic stem cell transplan- tation . Considering the results of these studies, there seems to be a difference in the importance of SNX10 in different organs in mice and humans. Using RANKL-stimulated peripheral blood monocytes we found that osteoclasts could be generated from cells isolated from the patients bearing the SNX10 splice site mutation. This is consistent with previous stud- ies, which were also able to generate osteoclasts from peripheral blood of patients with different SNX10 muta - 10 37, 43 tions , and from spleen or bone marrow cells of Snx10 KO mice . In contrast, silencing of Snx10 in the mouse pre-osteoclastic cell line Raw 264.7 was shown to inhibit RANKL-induced osteoclast differentiation . Moreover, SNX10 mRNA levels were increased during RANKL-stimulated osteoclastogenesis, consistent with a likely role in osteoclast function rather than formation. Although we have not been able to obtain bone biopsies from the patients, these findings suggest that the patients have an osteoclast-rich form of osteopetrosis. e os Th teoclasts generated from patients had normal expression of key osteoclast genes including those asso- ciated with osteoclast fusion and bone attachment, however, the osteoclasts were larger than in healthy controls. e Th re was not die ff rence in the number of nuclei per osteoclasts between patients and controls, indicating that the larger size is the result of increased spreading, similar to what is seen in other cases of IARO . By contrast, Aker et Scientific Repo R ts | 7: 3012 | DOI:10.1038/s41598-017-02533-2 11 www.nature.com/scientificreports/ al. reported that osteoclasts generated from patients with a missense mutation (Arg51Gln) in SNX10 were fewer in number and smaller than in controls and, similar to osteoclasts from the Snx10 null mouse, had a markedly 10, 37, 43 decreased capacity to dissolve calcium phosphate . Our data show that osteoclasts cultured from patients were unable to resorb bone, evidenced by both an absence of resorption pits and lack of CTX in the culture super- natant. This finding is identical to earlier findings with osteoclasts cultured from IARO patients bearing muta- 21 10 tions in PLEKHM1 . Interestingly, human osteoclasts with a missense SNX10 mutation , and mouse osteoclasts 37, 43 from Snx10 null mice although having reduced activity, still had some capacity to dissolve calcium phosphate film, whereas human osteoclasts from our patients did not show any evidence of bone resorbing activity. This discrepancy most likely is due to the use of crystalline calcium phosphate coating instead of bone as substrate for the assessments of the bone resorbing activity. Electron microscopy revealed that the inability of osteoclasts from osteopetrotic patients in the Västerbotten County to resorb bone is associated with defective ruffled border formation, very similar to patients with ARO 4, 21, 44 caused by mutations in other genes . Using phalloidin staining, we found that the SNX10 patient osteoclasts were still able to form F-actin rich podosome belts in plastic dishes and normal F-actin rings when cultured on bone. In agreement with these observations, electron microscopy demonstrated these osteoclasts formed a seal- ing zones with normal appearance. These observations indicate that the defect lies in the endosomal/lysosomal trafficking pathways associated with ruffled border formation. This is perhaps not surprising, since SNX10 is a member of the sorting nexin family of proteins that play crucial roles in cargo sorting in the endosomal path- 45, 46 4 way , moreover; other genetic mutations that cause osteoclast-rich ARO also ae ff ct this pathway . In contrast to these findings, osteoclasts generated from bone marrow macrophages with deletion of the Snx10 gene have been reported to form an abnormal F-actin positive podosome belt in plastic dishes . We have no explanation for these differences between mouse and human osteoclasts. As discussed below, we do not believe that the difference might be due to expression of a truncated, functional form of SNX10 in our patients. er Th e seem to be clear differences in the functional consequences for osteoclasts in patients with the Arg51Gln mutated SNX10 and those in the present study. u Th s, the Arg51Gln mutation results in decreased osteoclast for - mation and decreased osteoclast activity, whereas the SNX10 mutation found in the Västerbotten form of IARO impairs resorption by disrupting ruffled border formation in mature osteoclasts. Furthermore, osteoclasts with the Arg51Gln mutation exhibit extensive cytoplasmic vacuolation, which was not seen in osteoclasts with the Västerbotten mutation. Increased vacuolation has also been observed by overexpressing SNX10 in several human cell lines , which raises the possibility that the Arg51Gln mutation results in gain of function; however, overex- pression of the Arg51Gln form of SNX10 does not cause vacuolation , therefore, the functional consequences of this mutation remain unclear. No full length SNX10 protein can be expressed in osteoclasts with the Västerbotten mutation, as the mutation leads to a stop codon upstream from intron 4. We do not know whether a truncated form of SNX10 is expressed but we have not been able to conclusively show the presence of a truncated SNX10 protein using four different antibodies in Western blots, including antibodies recognizing the amino-terminal part of SNX10. SNX proteins contain a conserved phox-homology (PX) domain, which can bind phosphorylated phosphatidylinositols (PI), thereby targeting SNXs to PI-enriched membranes. The PX domain extends from c.30–381 in exons 3–6. It seems as if the PX domain in SNX10 is even further extended (PXe) in the C-terminal 48, 49 part . Since PX or PXe is the functional domain in SNXs, it is unlikely, even if a truncated SNX10 protein is expressed in osteoclasts carrying the Västerbotten mutation, that this truncated protein has any biological activity. e diff Th erence between the Arg51Gln mutation and the Västerbotten mutation is most likely responsible for the differences seen in osteoclast formation and function. In conclusion, the Västerbotten intermediate form of osteopetrosis is caused by a splice site sequence variant in the SNX10 gene and differs in clinical characteristics from other intermediate types of ARO. The age at onset is early in infancy, optic atrophy and anemia is present at diagnosis and tonsillar herniation, foramen magnum stenosis, and severe osteomyelitis/osteonecrosis of the jaw are common. Osteoclasts generated from peripheral blood monocytes of the patients showed defects in ruffled border formation and failed to resorb bone. These data extend the clinical spectrum of IARO and demonstrate a critical role for SNX10 in ruffled border formation. Materials and Methods Subjects. We collected clinical data, medical records, x-ray, computed tomography scan (CT) and magnetic resonance imaging (MRI) of the brain, and blood samples from five male and four female patients with IARO. We also collected blood samples from their parents and siblings. Patients descended from the county of Västerbotten, northern Sweden. Exome sequencing. DNA from three patients and 15 non-related individuals (controls) from the county of Västerbotten were sequenced. Exome enrichment was performed using 3 μg genomic DNA. DNA samples were sheared by sonication with the Covaris S2 instrument (Covaris, Inc.). Fragment libraries were constructed from the sheared samples using the AB Library Builder System (Life Technologies) and target enrichment was performed according to the manufacturer’s protocols (Agilent SureSelect Human All Exon v3 kit). Exome capture was conducted by hybridizing the DNA libraries with biotinylated RNA baits for 24 h followed by extraction using streptavidin coated magnetic beads. Captured DNA was then amplified followed by emulsion PCR using the EZ Bead System (Life Technologies) and sequenced on the SOLiD5500xl system. Individual libraries were labelled by a post-hybridization barcoding procedure (Agilent, Santa Clara, CA, USA; barcodes compatible with SOLiD sequencing technology). Alignment and variant calling. Alignment of color space reads to the human reference genome (GRCh37/ hg19) was performed using v2.1 of the Lifescope Software (Life Technologies). Single nucleotide variants (SNVs), Scientific Repo R ts | 7: 3012 | DOI:10.1038/s41598-017-02533-2 12 www.nature.com/scientificreports/ and small insertions and deletions (indels) were subsequently identified by the diBayes algorithm available within the Lifescope software. All identified SNVs and indels were imported into a local installation of the CanvasDB database system , for annotation and further analysis of the variants. DNA and cDNA sequencing, haplotype analysis, carrier frequency. DNA was extracted from whole blood using a standard salting-out method. SNX10 gene amplification and bidirectional sequencing of exons and intron-exon boundaries were performed in samples from all patients and their parents. RNA was extracted from cultured patient peripheral blood monocytes and converted to cDNA as described below for gene expression analysis. Bidirectional Sanger sequencing (Applied Biosystems BigDye Terminator v3.1 Cycle Sequencing Kit, Applied Biosystems) on a 3730xl DNA Analyzer (Applied Biosystems) was used for analysis of SNX10 cDNA (NM_001199837.1) using primer in exon 3 to 5 to confirm the mutation and to identify the effect of the splicing. Sequence analysis was performed using the Sequencher software (Gene Codes Corporation, Ann Arbor, MI, USA). Primers sequences where designed using Primer 3 Plus software (primer3plus.com/) and are available upon request. Five patients (Pt1–5) were analyzed using the GeneChip Mapping 250 K array (Affymetrix, Santa Clara, CA), and one patient (Pt9) was analyzed using the HumanOmniExpressExome-8v1array (Illumina), for haplotype analysis, and estimation of the age of the sequence variant. To calculate the carrier frequency, DNA from 1000 randomly chosen individuals from the county of Västerbotten was used (The medical bio-bank, Umeå university hospital, Umeå, Sweden) . In vitro culture of CD14 monocytes from patients and controls. Peripheral blood mononuclear cells were isolated by Ficoll-Paque PLUS separation and CD14 cells were labelled with CD14 MicroBeads and extracted using a MACS column according to the manufacturer’s instructions (Miltenyi). Cells were then seeded 5 2 in 96-well plates (3 × 10 cells/cm ) in complete α-MEM medium (Gibco cat no 22561-021) supplemented with 10% heat inactivated fetal bovine serum (FBS, Sigma cat no F7524), 2 mM GlutaMAX (Gibco cat no 35050-038), 50 µ g/ml gentamicin (Gibco cat no 15750-037), 100 U/ml penicillin and 100 µ g/ml streptomycin (Gibco cat no 15140-148) and 30 ng/ml human M-CSF (M, R&D Systems cat no 216-MC-025/CF). For osteoclast generation, 2 ng/ml recombinant mouse RANKL (RL, R&D Systems cat no 462-TEC-010) was added from the start of the experiments (day 0). Media was replenished every three days and the cells were stained for TRAP at different days of culture using the Acid Phosphatase, Leukocyte (TRAP) Kit from Sigma (cat no 387A). TRAP cells containing three or more nuclei were counted as TRAP multinucleated osteoclast (MuOCL). Analysis of osteoclast differentiation, gene expression and bone resorption. RANKL-stimulated osteoclast differentiation of CD14 cells was analyzed on plastic and on devitalized bovine bone discs (IDS Immunodiagnostics Systems cat no TDT-1BON1000-96). The surface area of TRAP positive cells on plastic TM were analyzed using the Osteomeasure software (v3.2.1.0, Osteometrics Inc.). Total RNA was purified from cells at day 3 using RNeasy Micro Kit (Qiagen). cDNA synthesis was performed using the High Capacity cDNA Reverse Transcription Kit from Applied Biosystems. Gene expression was analyzed in two different experiments with in total four healthy controls and two patients using predesigned TaqMan Assays (Life Technologies) and the StepOnePlus Real-Time PCR system (Applied Biosystems). We analyzed ACP5 (Hs00356261_m1), CTSK (Hs01080388_m1), NFATc1 (Hs00542678_m1), TCIRG1 (ATP6i, HS00990751_m1), CLC7 (Hs01126462_ m1), ITGAV (Hs00233808_m1), ITGB3 (Hs01001469_m1), DC-STAMP (Hs00229255_m1), OC-STAMP (Hs00875776_m1), ERVW-1 (Syncytin-1, Hs02341206_g1), SLC1A5 (Asct2, Hs00194540_m1), SNX10 exon 2–3 (Hs01007224_m1) and SNX10 exon 4–5 (HS01007226_m1). The mRNA abundance of each gene was adjusted for the expression of 18S (Life Technologies 4310893E) used as an internal control. In most figures, gene expres- sion is presented as the fold increase between M/RL and M treated cells. Statistical analysis was performed using One-Way ANOVA with Tukey’s post hoc test. The bone resorptive activity of osteoclasts on bone slices was ana - lyzed by measuring the release of C-terminal telopeptides of type I collagen (CTX) in culture supernatants using a commercial ELISA (IDS Immunodiagnostics Systems cat no AC-07F1) and by Toluidine blue staining of positive excavations. TRAP5b was analyzed in culture media using commercial ELISA (IDS Immunodiagnostics Systems cat no SB-TR201A). Staining for the actin ring and ruffled border. Actin organization was studied at day 4 on plastic and day 11 on bone. Cells were fixed in 4% paraformaldehyde, washed with PBS and then permeabilised using 0.1% Triton X-100 in PBS for 10 min at 4 °C. After washing again with PBS, the cells were stained with rhodamine conjugated phalloidin (Life Technologies cat no R415), 5 U/ml in 2% BSA/PBS, for 20 min at 4 °C. Following washing in PBS the cells were mounted in Prolong Gold Mountant (Life Technologies cat no P-36931). Cells cultured on bone were stained for phalloidin as above but including 100 µ g/ml FITC conjugated peanut aggluti- 22, 23 nin (PNA)-lectin (Sigma cat no L7381) to stain the ruffled border . Photographs were taken using widefield uo fl rescence microscopy. Confocal microscopy. Osteoclasts cultured on bovine bone discs were prepared as above, and, in addition to phalloidin, were stained with 5 μg/ml Alexa Fluor 633-conjugated wheat germ agglutinin to stain cell membranes and 10 μM Alexa Fluor 488-conjugated alendronate (AF-ALN, ref. 51) to stain bone surfaces. Bone discs were mounted on to glass slides in VectaShield with DAPI (to stain nuclei) mounting medium before analyzing on a Zeiss LSM700 confocal microscope. To assess potential differences in fusion, the number of nuclei was counted in 80 randomly selected osteoclasts from two control cultures and one patient culture. Scientific Repo R ts | 7: 3012 | DOI:10.1038/s41598-017-02533-2 13 www.nature.com/scientificreports/ Transmission Electron Microscopy investigation of ruffled border morphology. To investigate whether patient osteoclasts were capable of forming mature ruffled borders, cells cultured on bovine bone were analyzed by transmission electron microscopy (TEM). Briey fl , cell cultured with M-CSF and RANKL for eight days were fixed in 2.5% glutaraldehyde in 0.1 M cacodylate buffer. The bone discs were then demineralized in 0.15 M EDTA before being postfixed in 1% osmium tetroxide and embedded in Epon resin. Bone discs were sectioned perpendicular to the bone surface to obtain cross sections of osteoclasts and their ruffled border areas. Ultrathin sections were imaged using a JEOL Jem-1400Plus TEM equipped with an AMT UltraVUE camera. Following image acquisition, osteoclasts were graded based on their ruffled border morphological characteristics, similarly to a previous study, ref. 52: grade 0 - no ruffled border apparent; grade 1 - rudimentary ruffled border, with few poorly-developed membrane folds apparent; grade 2 - moderate ruffled border, with folded membrane apparent at some areas; grade 3 - extensive, highly-folded ruffled border. Statistics. For experiments with two patients and two controls, statistical analysis was performed using One-Way ANOVA with Tukey’s post hoc test. When one patient was compared to one control, a two-sided Student’s t test was used. A die ff rence was considered statistically significant with P < 0.05. Due to the limited number of patients and the large inter-patient variation, it was not possible to pool the results and make statistical comparisons based on the patients as one group. Instead, the results are shown for each patient separately. For gene expression analyses (Figs 3 and S1) each independent experiment is shown. Study approval. e R Th egional Ethical Review Board at Umeå University, Umeå, Sweden, approved the study. Patients, parents or first-degree relatives gave their written informed consent. Methods using human tissue were in accordance with the Helsinki declaration. References 1. Frattini, A. et al. Defects in TCIRG1 subunit of the vacuolar proton pump are responsible for a subset of human autosomal recessive osteopetrosis. Nature genetics 25, 343–346, doi:10.1038/77131 (2000). 2. Kornak, U. et al. Loss of the ClC-7 chloride channel leads to osteopetrosis in mice and man. Cell 104, 205–215, doi:10.1016/S0092- 8674(01)00206-9 (2001). 3. Sobacchi, C. et al. The mutational spectrum of human malignant autosomal recessive osteopetrosis. Human molecular genetics 10, 1767–1773, doi:10.1093/hmg/10.17.1767 (2001). 4. Sobacchi, C., Schulz, A., Coxon, F. P., Villa, A. & Helfrich, M. H. Osteopetrosis: genetics, treatment and new insights into osteoclast function. Nature reviews. Endocrinology 9, 522–536, doi:10.1038/nrendo.2013.137 (2013). 5. Villa, A., Guerrini, M. M., Cassani, B., Pangrazio, A. & Sobacchi, C. Infantile malignant, autosomal recessive osteopetrosis: the rich and the poor. Calcie fi d tissue international 84, 1–12, doi:10.1007/s00223-008-9196-4 (2009). 6. Askmyr, M. K., Fasth, A. & Richter, J. Towards a better understanding and new therapeutics of osteopetrosis. British journal of haematology 140, 597–609, doi:10.1111/j.1365-2141.2008.06983.x (2008). 7. Fasth, A. Osteopetrosis–more than only a disease of the bone. American journal of hematology 84, 469–470, doi:10.1002/ajh.21454 (2009). 8. Bliznetz, E. A. et al. Genetic analysis of autosomal recessive osteopetrosis in Chuvashiya: the unique splice site mutation in TCIRG1 gene spread by the founder effect. European journal of human genetics: EJHG 17, 664–672, doi:10.1038/ejhg.2008.234 (2009). 9. Phadke, S. R. et al. Malignant recessive osteopetrosis. Indian pediatrics 36, 69–74 (1999). 10. Aker, M. et al. An SNX10 mutation causes malignant osteopetrosis of infancy. Journal of medical genetics 49, 221–226, doi:10.1136/ jmedgenet-2011-100520 (2012). 11. Megarbane, A. et al. Homozygous stop mutation in the SNX10 gene in a consanguineous Iraqi boy with osteopetrosis and corpus callosum hypoplasia. European journal of medical genetics 56, 32–35, doi:10.1016/j.ejmg.2012.10.010 (2013). 12. Pangrazio, A. et al. SNX10 mutations define a subgroup of human autosomal recessive osteopetrosis with variable clinical severity. Journal of bone and mineral research: the official journal of the American Society for Bone and Mineral Research 28, 1041–1049, doi:10.1002/jbmr.1849 (2013). 13. Ameur, A., Bunikis, I., Enroth, S. & Gyllensten, U. CanvasDB: a local database infrastructure for analysis of targeted- and whole genome re-sequencing projects. Database: the journal of biological databases and curation 2014, bau098–bau098, doi:10.1093/ database/bau098 (2014). 14. Genomes Project, C. et al. A global reference for human genetic variation. Nature 526, 68–74, doi:10.1038/nature15393 (2015). 15. Hastbacka, J. et al. Linkage disequilibrium mapping in isolated founder populations: diastrophic dysplasia in Finland. Nature genetics 2, 204–211, doi:10.1038/ng1192-204 (1992). 16. Alatalo, S. L., Halleen, J. M., Hentunen, T. A., Monkkonen, J. & Vaananen, H. K. Rapid screening method for osteoclast differentiation in vitro that measures tartrate-resistant acid phosphatase 5b activity secreted into the culture medium. Clinical chemistry 46, 1751–1754 (2000). 17. Kukita, T. et al. RANKL-induced DC-STAMP is essential for osteoclastogenesis. e J Th ournal of experimental medicine 200, 941–946, doi:10.1084/jem.20040518 (2004). 18. Yagi, M. et al. DC-STAMP is essential for cell-cell fusion in osteoclasts and foreign body giant cells. The Journal of experimental medicine 202, 345–351, doi:10.1084/jem.20050645 (2005). 19. Yang, M., Birnbaum, M. J., Mackay, C. A., Mason-Savas, A., Thompson, B. & Odgren, P. R. Osteoclast stimulatory transmembrane protein (OC-STAMP), a novel protein induced by RANKL that promotes osteoclast differentiation. Journal of Cellular Physiology 215, 497–505, doi:10.1002/jcp.21331 (2007).  20. Soe, K. et al. Involvement of human endogenous retroviral syncytin-1 in human osteoclast fusion. Bone 48, 837–846, doi:10.1016/j. bone.2010.11.011 (2011). 21. Van Wesenbeeck, L. et al. Involvement of PLEKHM1 in osteoclastic vesicular transport and osteopetrosis in incisors absent rats and humans. e J Th ournal of clinical investigation 117, 919–930, doi:10.1172/JCI30328 (2007). 22. Takagi, M., Yagasaki, H., Baba, T. & Baba, H. Ultrastructural visualization of selective peanut agglutinin binding sites in rat osteoclasts. The journal of histochemistry and cytochemistry: official journal of the Histochemistry Society 36, 95–101, doi:10.1177/36.1.2447153 (1988). 23. Salo, J., Metsikko, K., Palokangas, H., Lehenkari, P. & Vaananen, H. K. Bone-resorbing osteoclasts reveal a dynamic division of basal plasma membrane into two different domains. Journal of cell science 109(Pt 2), 301–307 (1996). 24. Palagano, E. et al. Buried in the Middle, But Guilty: Intronic Mutations in the TCIRG1 Gene Cause Human Autosomal Recessive Osteopetrosis. Journal of bone and mineral research: the official journal of the American Society for Bone and Mineral Research, doi:10.1002/jbmr.2517 (2015). Scientific Repo R ts | 7: 3012 | DOI:10.1038/s41598-017-02533-2 14 www.nature.com/scientificreports/ 25. Kahler, S. G., Burns, J. A. & Aylsworth, A. S. A mild autosomal recessive form of osteopetrosis. American journal of medical genetics 17, 451–464, doi:10.1002/ajmg.1320170208 (1984). 26. Whyte, M. P. Primer on the Metabolic Bone Diseases and Disorders od Mineral Metabolism. 4th edn, 367–383 (Lippincott Williams & Wilkins, 1999). 27. Campos-Xavier, A. B., Saraiva, J. M., Ribeiro, L. M., Munnich, A. & Cormier-Daire, V. Chloride channel 7 (CLCN7) gene mutations in intermediate autosomal recessive osteopetrosis. Human genetics 112, 186–189, doi:10.1007/s00439-002-0861-9 (2003). 28. McCleary, L., Rovit, R. L. & Murali, R. Case report: myelopathy secondary to congenital osteopetrosis of the cervical spine. Neurosurgery 20, 487–489, doi:10.1227/00006123-198703000-00025 (1987). 29. Elster, A. D., Theros, E. G., Key, L. L. & Chen, M. Y. Cranial imaging in autosomal recessive osteopetrosis. Part I. Facial bones and calvarium. Radiology 183, 129–135, doi:10.1148/radiology.183.1.1549658 (1992). 30. Sari, A. & Demirci, A. Radiographic type I autosomal dominant osteopetrosis with syringohydromyelia. Neuroradiology 38, 532–533, doi:10.1007/BF00626090 (1996). 31. Al-Tamimi, Y. Z., Tyagi, A. K., Chumas, P. D. & Crimmins, D. W. Patients with autosomal-recessive osteopetrosis presenting with hydrocephalus and hindbrain posterior fossa crowding. Journal of neurosurgery. Pediatrics 1, 103–106, doi:10.3171/PED-08/01/103 (2008). 32. Cure, J. K., Key, L. L., Goltra, D. D. & VanTassel, P. Cranial MR imaging of osteopetrosis. AJNR. American journal of neuroradiology 21, 1110–1115 (2000). 33. Stark, Z., Pangrazio, A., McGillivray, G. & Fink, A. M. Association of severe autosomal recessive osteopetrosis and structural brain abnormalities: a case report and review of the literature. European journal of medical genetics 56, 36–38, doi:10.1016/j. ejmg.2012.10.001 (2013). 34. Rees, H., Ang, L. C., Casey, R. & George, D. H. Association of infantile neuroaxonal dystrophy and osteopetrosis: a rare autosomal recessive disorder. Pediatric neurosurgery 22, 321–327, doi:10.1159/000120923 (1995). 35. Ben Hamouda, H. et al. Association of severe autosomal recessive osteopetrosis and Dandy-Walker syndrome with agenesis of the corpus callosum. Acta orthopaedica Belgica 67, 528–532 (2001). 36. Abinun, M. & Pieniazek, P. Successful haematopoietic stem cell transplantation for osteopetrosis due to TCRIG1 mutation. Archives of disease in childhood 95, 984–984, doi:10.1136/adc.2010.194928 (2010). 37. Ye, L. et al. Osteopetrorickets due to Snx10 deficiency in mice results from both failed osteoclast activity and loss of gastric acid- dependent calcium absorption. PLoS genetics 11, e1005057, doi:10.1371/journal.pgen.1005057 (2015). 38. Whyte, M. P. et al. Osteopetrosis, renal tubular acidosis and basal ganglia calcification in three sisters. The American journal of medicine 69, 64–74, doi:10.1016/0002-9343(80)90501-X (1980). 39. Bourke, E., Delaney, V. B., Mosawi, M., Reavey, P. & Weston, M. Renal tubular acidosis and osteopetrosis in siblings. Nephron 28, 268–272, doi:10.1159/000182216 (1981). 40. Ohlsson, A., Stark, G. & Sakati, N. Marble brain disease: recessive osteopetrosis, renal tubular acidosis and cerebral calcification in three Saudi Arabian families. Developmental medicine and child neurology 22, 72–84, doi:10.1111/j.1469-8749.1980.tb04307.x (1980). 41. Patel, P. J., Kolawole, T. M., al-Mofada, S., Malabarey, T. M. & Hulailah, A. Osteopetrosis: brain ultrasound and computed tomography findings. European journal of pediatrics 151, 827–828, doi:10.1007/BF01957934 (1992). 42. Zhu, C. H., Morse, L. R. & Battaglino, R. A. SNX10 is required for osteoclast formation and resorption activity. Journal of cellular biochemistry 113, 1608–1615, doi:10.1002/jcb.24029 (2012). 43. Zhou, C. et al. Deficiency of sorting nexin 10 prevents bone erosion in collagen-induced mouse arthritis through promoting NFATc1 degradation. Ann Rheum Dis 75, 1211–1218, doi:10.1136/annrheumdis-2014-207134 (2016). 44. Bruder, E., Stallmach, T., Peier, K., Superti-Furga, A. & Vezzoni, P. Osteoclast morphology in autosomal recessive malignant osteopetrosis due to a TCIRG1 gene mutation. Pediatric pathology & molecular medicine 22, 3–9 (2003). 45. Teasdale, R. D. & Collins, B. M. Insights into the PX (phox-homology) domain and SNX (sorting nexin) protein families: structures, functions and roles in disease. Biochem J 441, 39–59, doi:10.1042/BJ20111226 (2012). 46. Carlton, J., Bujny, M., Rutherford, A. & Cullen, P. Sorting nexins–unifying trends and new perspectives. Traffic 6, 75–82, doi:10.1111/j.1600-0854.2005.00260.x (2005). 47. Qin, B., He, M., Chen, X. & Pei, D. Sorting nexin 10 induces giant vacuoles in mammalian cells. e J Th ournal of biological chemistry 281, 36891–36896, doi:10.1074/jbc.M608884200 (2006). 48. Xu, T. et al. Structure of human SNX10 reveals insights into its role in human autosomal recessive osteopetrosis. Proteins 82, 3483–3489, doi:10.1002/prot.24689 (2014). 49. Xu, J. et al. Structure of sorting nexin 11 (SNX11) reveals a novel extended phox homology (PX) domain critical for inhibition of SNX10-induced vacuolation. e J Th ournal of biological chemistry 288, 16598–16605, doi:10.1074/jbc.M112.449306 (2013). 50. Hallmans, G. et al. Cardiovascular disease and diabetes in the Northern Sweden Health and Disease Study Cohort - evaluation of risk factors and their interactions. Scandinavian journal of public health. Supplement 61, 18–24, doi:10.1080/14034950310001432 (2003). 51. o Th mpson, K., Rogers, M. J., Coxon, F. P. & Crockett, J. C. Cytosolic entry of bisphosphonate drugs requires acidification of vesicles aer fl ft uid-phase endocytosis. Molecular pharmacology 69, 1624–1632, doi:10.1124/mol.105.020776 (2006). 52. Taylor, A. et al. Impaired prenylation of Rab GTPases in the gunmetal mouse causes defects in bone cell function. Small GTPases 2, 131–142, doi:10.4161/sgtp.2.3.16488 (2011). Acknowledgements We sincerely thank the patients and family members who participated in this study. We would also like to thank Stefan Esher, Umeå University, for help with genealogy, and Anna Westerlund for excellent technical assistance. This work was supported by grants from the FOU, at the Umeå university hospital, and the Medical Faculty at Umeå University. The work at University of Gothenburg was supported by grants from The Swedish Research Council, the Swedish Rheumatism Association, the Royal 80-Year Fund of King Gustav V, ALF/LUA research grant from Sahlgrenska University Hospital in Gothenburg and the Lundberg Foundation. The work at the University of Gothenburg and the University of Aberdeen was supported by Euroclast, a Marie Curie FP7-People- 2013-ITN: # 607446. Author Contributions E.L.S., P.H., J.K., N.D. and U.H.L. designed the study. E.L.S., P.H., J.K., E.M., C.S.B., P.E.S., T.L., M.H.H., F.C., G.H., J.W., N.D. and U.H.L. acquired the data. E.L.S., T.L. and J.W. re-evaluated radiological examinations. J.K. and A.A. performed exome sequencing, Sanger sequencing and haplotype analysis. A.A., T.G.K. and P.R. filtered exome data. P.H., E.M., M.H.H., F.C. and U.H.L. conducted experiments. E.L.S., P.H., J.K., E.M., C.S.B., T.L., A.A., M.H.H., F.C., N.D., J.W. and U.H.L. wrote the manuscript. Scientific Repo R ts | 7: 3012 | DOI:10.1038/s41598-017-02533-2 15 www.nature.com/scientificreports/ Additional Information Supplementary information accompanies this paper at doi:10.1038/s41598-017-02533-2 Competing Interests: The authors declare that they have no competing interests. Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Cre- ative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not per- mitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/. © The Author(s) 2017 Scientific Repo R ts | 7: 3012 | DOI:10.1038/s41598-017-02533-2 16 http://www.deepdyve.com/assets/images/DeepDyve-Logo-lg.png Scientific Reports Springer Journals

SNX10 gene mutation leading to osteopetrosis with dysfunctional osteoclasts

Download PDF
16 pages

Loading next page...
 
/lp/springer_journal/snx10-gene-mutation-leading-to-osteopetrosis-with-dysfunctional-dkAJTB302P

References (73)

Publisher
Springer Journals
Copyright
Copyright © 2017 by The Author(s)
Subject
Science, Humanities and Social Sciences, multidisciplinary; Science, Humanities and Social Sciences, multidisciplinary; Science, multidisciplinary
eISSN
2045-2322
DOI
10.1038/s41598-017-02533-2
Publisher site
See Article on Publisher Site

Abstract

www.nature.com/scientificreports OPEN SNX10 gene mutation leading to osteopetrosis with dysfunctional osteoclasts Received: 8 November 2016 1,2 3 2 4 Eva-Lena Stattin , Petra Henning , Joakim Klar , Emma McDermott , Christina Stecksen- 5 6 7 7 8 Accepted: 12 April 2017 Blicks , Per-Erik Sandström , Therese G. Kellgren , Patrik Rydén , Göran Hallmans , Torsten 9 2 4 4 2 Published: xx xx xxxx Lönnerholm , Adam Ameur , Miep H. Helfrich , Fraser P. Coxon , Niklas Dahl , Johan 9 3,10 Wikström & Ulf H. Lerner Autosomal recessive osteopetrosis (ARO) is a heterogeneous disorder, characterized by defective osteoclastic resorption of bone that results in increased bone density. We have studied nine individuals with an intermediate form of ARO, from the county of Västerbotten in Northern Sweden. All afflicted individuals had an onset in early infancy with optic atrophy, and in four patients anemia was present at diagnosis. Tonsillar herniation, foramen magnum stenosis, and severe osteomyelitis of the jaw were common clinical features. Whole exome sequencing, verified by Sanger sequencing, identified a splice site mutation c.212 + 1 G > T in the SNX10 gene encoding sorting nexin 10. Sequence analysis of the SNX10 transcript in patients revealed activation of a cryptic splice site in intron 4 resulting in a frame shift and a premature stop (p.S66Nfs * 15). Haplotype analysis showed that all cases originated from a single mutational event, and the age of the mutation was estimated to be approximately 950 years. Functional analysis of osteoclast progenitors isolated from peripheral blood of patients revealed that stimulation with receptor activator of nuclear factor kappa-B ligand (RANKL) resulted in a robust formation of large, multinucleated osteoclasts which generated sealing zones; however these osteoclasts exhibited defective ruffled borders and were unable to resorb bone in vitro. In healthy individuals, normal bone homeostasis is maintained by balancing bone resorption by osteoclasts with new bone formation by osteoblasts. Osteopetrosis constitutes a heterogeneous group of rare disorders characterized by increased bone mass caused either by an inability to form osteoclasts or by loss of osteoclast 1–4 function . Bone resorbing multinucleated osteoclasts are formed by the differentiation and fusion of mono- nuclear hematopoietic precursor cells in the monocyte lineage. Their differentiation is mediated by macrophage colony-stimulating factor (M-CSF), which is required for proliferation and survival of the mononucleated cells, and receptor activator of nuclear factor kB ligand (RANKL), which triggers differentiation along the osteoclastic lineage. The main clinical features of osteopetrosis are an increased bone density with brittle bones and multiple fractures, narrowing of foramina in the skull with neurological symptoms, and reduced bone marrow cavities 4, 5 leading to impaired hematopoiesis, anemia and recurrent infections . Department of Medical Biosciences, Medical and Clinical Genetics, Umeå University, 901 87, Umeå, Sweden. Department of Immunology, Genetics and Pathology, Science for Life Laboratory, Uppsala University, 751 85, Uppsala, Sweden. Centre for Bone and Arthritis Research, Department of internal medicine and clinical nutrition, Institute of Medicine, Sahlgrenska Academy, University of Gothenburg, 405 30, Gothenburg, Sweden. Arthritis and Musculoskeletal Medicine Programme, Institute of Medical Sciences, University of Aberdeen, Foresterhill, Aberdeen, AB25 2ZD, UK. Pediatric Dentistry, Department of Odontology, Faculty of Medicine, Umeå University, 901 87, 6 7 Umeå, Sweden. Department of Pediatrics, Umeå University, 901 87, Umeå, Sweden. Department of Mathematics and Mathematical Statistics, Computational Life science Cluster (CLiC), Umeå University, 901 87, Umeå, Sweden. 8 9 Department of Biobank Research, Umeå University, 901 87, Umeå, Sweden. Department of Surgical Sciences, Radiology, Uppsala University, 751 85, Uppsala, Sweden. Molecular Periodontology, Department of Odontology, Faculty of Medicine, Umeå University, 901 87, Umeå, Sweden. Eva-Lena Stattin and Petra Henning contributed equally to this work. Correspondence and requests for materials should be addressed to E.-L.S. (email: eva-lena. stattin@igp.uu.se) or P.H. (email: petra.henning@gu.se) Scientific Repo R ts | 7: 3012 | DOI:10.1038/s41598-017-02533-2 1 www.nature.com/scientificreports/ Autosomal recessive osteopetrosis (ARO) comprises the infantile, or “malignant” form, and a less severe inter- mediate form (IARO). In malignant ARO, the skeletal defects have an onset before birth and the resulting nerve damage due to entrapment is irreversible. Treatment with hematopoietic stem cell transplantation (HSCT) is 4, 6 critical in order to prevent disease progression and early death . To date, at least 10 genes coding for osteo- clast function or osteoclast differentiation have been demonstrated to be causally involved in the pathogenesis of 1–4 osteopetrosis . The incidence of ARO is estimated to be 1 in 250,000 births, but can be much higher in isolated 7–9 consanguineous populations . Recently, a missense mutation in the gene encoding sorting nexin 10 (SNX10 [MIM614780]) was shown to be the cause of ARO (OPTB8 [MIM615085]) in consanguineous families of Palestinian origin . Since then, dif- 11, 12 ferent mutations in the SNX10 gene have been reported in osteopetrosis patients of diverse ethnicities . At the initiation of the present study, we were aware of several cases with osteopetrosis of unknown origin in the county of Västerbotten, Sweden. A genealogy study revealed that the patients belonged to the same family. We therefore performed exome sequencing with the aim to find the pathogenic sequence variant causing the disease and, in addition, performed cellular and functional studies to assess the role of the mutation in bone remodeling. During the course of our studies, three patients with osteopetrosis from the Västerbotten County were reported to have a mutation in the SNX10 gene . We confirmed that all osteopetrosis patients in the Västerbotten County carry the same SNX10 mutation. Here, we describe the clinical features and natural course of the disease in the largest cohort of patients with IARO reported thus far. We identified the disease-causing mutation in SNX10, as well as its origin and age and report detailed functional analyses on osteoclasts carrying this pathogenic sequence variant. Results Subjects. We studied nine patients, v fi e males and four females, with IARO from the county of Västerbotten in Northern Sweden. The clinical data are summarized in Table  1. Manifestation of the disease was seen in early infancy in all patients, but the age at correct diagnosis varied from 7 weeks (Pt1) to 3 years (Pt9). The patient (Pt9) diagnosed at the age of 3 years was first thought to have a hereditary optic atrophy. Delayed tooth eruption prompted radiology and the diagnosis of osteopetrosis became apparent (Fig. 1A–D). Skeletal features, radiographs and fractures. Birth length was within ±2 SD (range 48–52 cm), and head circumference ±1 SD (range 34–35 cm). All patients had short final height (males: 132 cm, (−6.7 SD), 160 cm, (−2.5 SD), females: 149 cm, (−2.5 SD), 159 cm, (−1.1 SD) (Table 1)). In all patients, skeletal radiographs show a generalized increase in bone density with metaphyseal modelling defects, transverse bands of greater and lesser density in tubular bones, short wide femoral neck, and”bone-within-bone” pattern of the phalanges (Fig. 1A–D). There were no radiological signs of rickets in five evaluated patients (Pt1, Pt5, and Pt7–9). Seven of nine patients suffered from recurrent fractures up to about 30 times (Pt1), even aer min ft or trauma. During child- hood, almost all fractures were transverse diaphyseal fractures in the extremities without pronounced dislocation and they healed normally with periosteal callus. e Th adult patients (Pt1–5) also had femoral shaft fractures with prolonged healing or pseudoarthrosis despite external or internal fixation. One boy (Pt7), who received hemato- poietic stem cell transplantation (HSCT), did not fracture and his growth normalized (Table 1). Neurological symptoms. All patients developed visual impairment before the age of one year, and five patients became blind (Pt1, Pt4, Pt6, Pt8, and Pt9). Decompression of the optic nerve was not performed in any of the patients. Retinal atrophy was not observed. All but one patient suffers or suffered from hearing impairment, and one patient (Pt6) developed a profound deafness (Table 1). Aer b ft irth, the head circumference increased rapidly in six patients. One girl (Pt9) was found to have hydro- cephalus and cerebellar tonsillar herniation, requiring a ventriculo-peritoneal shunt. She was operated at the age of eight years with decompression of the brainstem, and she died unexpectedly, two years later, possibly due to compression of a postsurgical cephalocele. At that time, she was hospitalized for evaluation of severe back pain, and headache episodes. Two patients (Pt8, and Pt5) developed stenosis of the foramen magnum; one (Pt8) died at the age of 12 years due to spinal cord compression while the other (Pt5) is still alive at the age of 24 years. One boy (Pt6) showed intellectual disability in addition to early blindness and profound deafness. Facial paralysis was present in four patients (Pt1–4). Evaluation of six magnetic resonance imaging (MRI), and three computed tomography (CT) scans showed a thickened and sclerotic bone of the skull in all but one patient (Table 2). Scans did not show brain malforma- tions (Fig. 1E,F). Hematological symptoms. Anemia was present early in life in four patients (Pt1, Pt2, Pt7, and Pt9) due to bone marrow encroachment. Consumption and hemolysis of erythrocytes, followed by a low leucocyte count and throm- bocytopenia, was seen in two patients (Pt8, and Pt9); haptoglobin and carboxyhemoglobin (CoHb) were below the detection limit. Splenectomies were performed on these patients and the hematopoiesis almost normalized aer ft wards. Six of nine patients require or required repeated blood transfusions due to normocytic anemia (Table 1). All patients were oer ff ed stem cell transplantation, most of them during adolescence, but HSCT was accepted only by the parents of the youngest one (Pt7). He received HSCT at the age of 33 months from an HLA-matched unrelated donor; he is still alive and healthy at the age of 15 years. However, he suffers from irreversible tissue damage that occurred prior to transplantation: he has visual deterioration as well as missing and malformed teeth. Two patients (Pt4, and Pt6) died due to a combination of pancytopenia and septicemia; one (Pt6) died at the age of 19, due to septicemia originating from osteomyelitis of the jaw, and the other (Pt4) died at the age of 35, due to chronic osteomyelitis of one femur aer a ft n unhealed fracture. Autopsy was not performed in any of the patients (Table  1). Scientific Repo R ts | 7: 3012 | DOI:10.1038/s41598-017-02533-2 2 www.nature.com/scientificreports/ Patient number 1 2 3 4 5 6 7 8 9 Gender Male Female Male Female Female Male Male Male Female age (year) 47 43 46 35 24 19 15 12 10 51 cm 50 cm 52 cm NA 48 cm 52 cm NA (±0SD) NA 48 cm Birth length, 4000 g 3620 g 3910 g 3110 g 2930 g 3500 g 3094 g NA 3305 g weight, OFC 34 cm NA NA NA NA 34.5 cm NA NA 35 cm Normal growth, 132 cm 149 cm 160 cm 159 cm NA NA — — Final height aer ft HSCT (age year) (SD) (17)(−6.7) (20)(−2.5) (20)(−2.9) (19)(−1.1) Age at onset Early infancy Early infancy Early infancy Early infancy Early infancy Early infancy Early infancy Early infancy Early infancy 36 months. 5 months. 28 months. Retained 6 months. 19 months. Optic atrophy. teeth. Nasal 7 weeks. Optic atrophy. Optic atrophy. Age at 2 years. Optic 12 months. Days, disease Macrocephaly. Macrocephaly Optic atrophy. Optic atrophy. Macrocephaly congestion. diagnosis, atrophy. Optic atrophy, known in the 46 cm. Macrocephaly Macrocephaly 55.5 cm. Optic atrophy. Anaemia ( Hb symptoms Macrocephaly blindness family 6.64 mmol/l) Anaemia (Hb 44.5 cm 50 cm Anaemia (Hb Macrocephaly. 6.09 mmol/l) 6.33 mmol/l) Anaemia (Hb 6.6 mmol/l) Stem cell HSCT at the age — — — — — — — — transplantation of 33 month Deceased (age) — — — 35 years — 19 years — 12 years 10 years Blindness. Hydrocephalus, Blindness. Impaired vision VP-shunt. Impaired vision, Impaired Impaired vision, and hearing. Blindness. Decompression Blindness. hearing. Blindness. Impaired (blind on one (blind on one Foramen Profound op of the Neurological hearing, unilateral eye). Impaired Impaired Died due eye). Impaired magnum deafness. Impaired vision brainstem at hearing. Facial to foramen symptoms profound deafness. hearing, tuba hearing. Facial stenosis, op Intellectual the age of 8, Facial paralysis aperta. Facial paralysis magnum paralysis decompression disability died due to paralysis stenosis at the at the age of 14 compression of age of 12 a postsurgical cephalocele Anaemia, hepato- Anaemia, Anaemia, splenomegaly. Anaemia, hepato- Anaemia, Anaemia, hepato- hepato- Blood splenomegaly, hepato- hepato- Anaemia, splenomegaly, splenomegaly, Haematological transfusions Almost normal Normal splenectomy at splenomegaly. splenomegaly. hepato- splenectomy splenectomy symptoms Died due to haematopoiesis. haematopoiesis. the age of 5. Blood Blood Blood splenomegaly. at the age at the age septicaemia transfusions. transfusions. transfusions. of 9. Blood of 5. Blood originating transfusions. transfusions. from osteomyelitis. At the age of 28 months, two At the age of Eruption of At the age of At the age of At the age of At the age of primary teeth 31 months, primary teeth NA   4-6 months, ten 6 months, five NA   6 months, 12 NA 13 months, six Malformed and eight primary (month), Malformed and primary teeth primary teeth Malformed and primary teeth Malformed primary teeth retained but a teeth and teeth retained Malformed and Malformed and retained Malformed and and retained Malformed lot better than Malformed and development retained retained retained and retained patients not retained transplanted At the age of 6. At the age of Osteonecrosis of the At the age of 13. At the age of 30. 12. Died due ONJ (year) ear canal, and the Osteonecrosis Osteonecrosis ONJ (age NA) At the age of 18. to septicaemia None None None skull with associated of the ear canal of the ear canal originating epidural abscess from ONJ Narrow nasal, Nasal Early nasal Narrow nasal, Narrow nasal, and and pharynx obstruction, obstruction, and pharynx pharynx cavities Narrow nasal, Narrow nasal, Narrow nasal, cavities before the age Airway before the age cavities CPAP Tracheostomy at the and pharynx and pharynx and pharynx None Apnoea during of 12 month, symptoms of 5 months. during night age of 22. CPAP during cavities cavities cavities the night adenoidectomy CPAP during since the age night Tracheostomy Tracheostomy night of 35. at the age of 9. at the age of 8. Multiple Multiple Multiple Multiple fractures, Multiple Multiple Multiple Bone fractures fractures, fractures, fractures, None Few fractures pseudoarthrosis fractures fractures fractures pseudoarthrosis pseudoarthrosis pseudoarthrosis Exophthalmos keratitis, Exophthalmos endophthalmitis, Exophthalmos, Exophthalmos, Exophthalmos Drooling, Exophthalmos enucleation of one Bilat pes Mitral valve Inguinal Exophthalmos, Exophthalmos, treated Media otitis, Other features eye. Ear canal stenosis, calcaneovalgus, insufficiency, hernia, Exophthalmos Drooling, Severe back with Botox, Wheelchair Kidney tumour, Hypertension, Wheelchair Wheelchair Anosmia pain Wheelchair bound Urinary stones, Weight Walking frame bound bound bound loss, Wheelchair bound Table 1. Clinical features of nine individuals with a splice site mutation in the SNX10 gene, and intermediate A B autosomal recessive osteopetrosis (IARO). OFC = occipital-frontal circumference. NA = Not available. C D E SD = Standard deviations. HSCT = hematopoietic stem cell transplantation. Symptoms that lead to F G H the diagnosis. Hb = hemoglobin. VP = ventriculo-peritoneal shunt. ONJ = osteonecrosis of the jaw. CPAP = continuous positive airway pressure. Scientific Repo R ts | 7: 3012 | DOI:10.1038/s41598-017-02533-2 3 www.nature.com/scientificreports/ Figure 1. A Radiographs of subjects with the Västerbotten form of intermediate osteopetrosis, and sequencing of the SNX10 gene. (A–D) Skeletal radiographs of a girl at the age of three years (Pt9), showing characteristic features of intermediate osteopetrosis. (A) Skull: Increased bone density of the skull base, frontal bossing and VP-shunt (B) AP view of the legs: “Erlenmeyer flask” shape of the femoral diaphysis. Longitudinal and transverse bands of lesser density in the metaphyses. Bone-within-bone pattern in the epiphyses in the knee. (C) Lateral view of the spine: Very dens endplates of the vertebral bodies, so called “sandwich vertebrae”. The ribs (seen in part) are broad and dens. (D) AP view of the left arm: No fractures or sign of rickets, “bone- within-bone” pattern of the phalanges. (E,F) Magnetic resonance imaging of the brain, and spine in an 11-year-old boy. MRI shows cerebellar tonsillar descent through the foramen magnum, a foramen magnum narrowing, and signs of brain stem compression with resulting syringo-hydromyelia. (G,H) Consequence of the c.212 + 1 G > T variant on (G) genomic level and (H) transcript level. (G) Sequencing of DNA from a control (top), a heterozygous carrier (middle) and a homozygous patient (bottom) for the c.212 + 1 G > T variant. The variant changes the donor splice site at the 5′ end of intron 4 from GT to TT. This leads to a predicted use of an alternative donor splice site 16 nucleotides upstream of intron 4. (H) Sequencing of cDNA from a control (top) and a patient (bottom) confirms the use of an alternative splice site and the skipping of 16 bp leading to a frame shift and a premature stop codon (p.S66Nfs * 15) at the protein level. Symptoms of the airways, jaw and teeth. All patients had narrow nasal, and pharyngeal cavities, and two patients (Pt2, and Pt9) had documented nasal congestion before the age of one year. Three patients (Pt1, Pt8, and Pt9) needed a tracheostomy at the age of 8, 9 and 22 years, respectively. Three adult patients (Pt1–3) are treated with continuous positive airway pressure (CPAP) during nighttime, due to sleep apnea syndrome. Tooth Scientific Repo R ts | 7: 3012 | DOI:10.1038/s41598-017-02533-2 4 www.nature.com/scientificreports/ Patient number 1 1 1 4 5 7 8 8 9 Gender (male/female) Male Male Male Female Female Male Male Male Female Age at the time of imaging (year) 29 34 33 35 14 2 11 11 8 MRI/CT-scan MRI MRI CT CT MRI MRI MRI CT MRI Corpus callosum aplasia/hypoplasia 0 0 X X 0 0 0 X 0 Sclerotic bone 1 1 1 1 1 0 1 1 1 Thickened bone 1 1 1 1 1 0 1 1 1 Ventriculomegaly 0 0 0 0 0 0 1 1 0 Remodelled inner table 1 1 1 0 1 0 0 1 0 Tonsillar herniation 1 1 X 0 1 0 1 X 1 Optic nerve sheath dilatation 0 0 0 0 1 0 1 1 0 Proptosis 1 1 1 1 0 0 1 X 1 Brain atrophy 0 0 0 0 0 0 0 0 0 Optic nerve atrophy 0 0 0 X X 1 1 X 1 Optic canal stenosis 0 1 1 1 0 0 1 0 0 Ear fluid 1 1 X X 0 1 0 X 0 Dural venous sinus stenosis X X X 0 1 X 1 X 1 Foramen magnum stenosis 0 0 X X 1 0 1 X 1 Posteriorly angulated dens with stenosis X X X X 1 0 1 X 1 Cervical spinal stenosis in included upper 1 1 0 0 0 0 0 0 0 part Subdural hematoma or intra-cerebral 0 0 0 0 0 0 0 0 0 bleeding Prominent collateral veins between sinus 0 0 0 0 0 0 1 0 0 transversus and extra-cranial veins Lemon shape of the scull (similar to what is seen in Chari II malformation), syrinx 0 0 0 0 0 0 1 0 0 of the cervical cord, large occipital horns Calcification in basal ganglia and 0 0 0 0 0 0 0 1 0 subcortical frontal lobes Table 2. Magnetic resonance imaging (MRI) and computed tomography scan (CT) of the skull and brain, in six individuals with SNX10 related intermediate autosomal recessive osteopetrosis (IARO). 1 = present, 0 = not present, X = feature not possible to evaluate. eruption was delayed in three patients (Pt7–9), and teeth are congenitally missing, malformed or retained in all patients. Severe osteomyelitis and osteonecrosis of the jaw is seen in all adult patients (Pt1–6), with abscess for- mations and fistulas to the orbits and skin. One patient (Pt1) presented with osteonecrosis of the jaw at the age of six years, other patients (Pt2, Pt3, Pt5, and Pt6) at the age of 12, 13, 18 and 30 years, respectively. Osteonecrosis of the ear canal is seen in three of the eldest patients (Pt1–3) (Table 1). Exome sequencing, data filtering, DNA and cDNA sequencing, carrier frequency. Exome sequencing of three patients (Pt1, Pt2, and Pt5) revealed 115,283 unique SNPs. Among those, 25,301 were com- mon, 24,867 remained aer ex ft cluding SNPs on chromosomes X and Y, and 9,401 were homozygous. Nineteen of those were not present in the control population and three were non-synonymous, and located in three different genes. One of the three genes was the SNX10 gene. Homozygous sequence variants at the donor splice site of exon 4 (c.212 + 1 G > T) were confirmed in all ae ff cted patients using Sanger sequencing, and the parents were found to be heterozygous for the splice site mutation (Fig. 1G–H). The splice site variant is not present in any of >1000 13, 14 exomes present in our local Canvas DB, or in controls in Exome Seq. project or 1000 Genomes . Sequence analysis of cDNA from cultured patient peripheral blood cells revealed that an alternative splice site, located 16 bases upstream from intron 4, was used, resulting in a frame shift and a stop codon (p.S66Nfs * 15). The carrier frequency of the sequence variant in the county of Västerbotten is 1:93 or 1.1% of the population. Genealogy and haplotype analysis. Genealogic data revealed a common ancestor 8–12 generations th back, in the early 19 century. He was born 1820 in Bureå, Skellefteå, County of Västerbotten, northern Sweden. Haplotype analysis was performed to evaluate the age of the sequence variant. The size of the homozygous regions in Pts1–5 is 4044451 bp, 3494094 bp, 3515023 bp, 5632740 bp, 4628862 bp, and 4628862 bp, respectively. Pt9 showed homozygosity between marker rs477644 and rs6462058, corresponding to 11307603 bp, which was included in the age calculations (Fig. S1). e Th linkage disequilibrium in six ARO probands extends over a mean homozygous region of 5.3 cM, which equals 38 meioses/generations. With an assumed generation time of 25 years, this corresponds to a shared ancestral origin of the SNX10 gene variant 950 years ago . Monocytes from patients form large osteoclasts in vitro. We isolated CD14 monocytes from the peripheral blood of three patients, (Pt1, Pt2 and Pt5) and seven healthy controls, designated C1–C7. We assessed Scientific Repo R ts | 7: 3012 | DOI:10.1038/s41598-017-02533-2 5 www.nature.com/scientificreports/ Figure 2. Osteoclast differentiation of CD14 monocytes from peripheral blood. Photographs of TRAP staining (A) and quantification of surface area (B ) of CD14 cells cultured in M-CSF (M) or M-CSF and RANKL (M/RL) for four days (C1, C2, Pt2, Pt5) or five days (C1, lower photos, and Pt1). Surfaces are presented as mean ± SEM. (C) Actin ring staining of cells cultured in M/RL for four days. Scale bars: 100 μm. *** P ≤ 0.001. osteoclastogenesis on both plastic and bone slices, by incubating CD14 monocytes with RANKL. Cells from patients were tested alongside control cells in six different independent experiments. The data are presented sep- arately for these individual experiments instead of being pooled because of the large inter-individual variability that is commonly seen in primary osteoclast cultures from humans. Multinucleated cells that stained positive for tartrate resistant acid phosphatase (TRAP ) were formed on + + plastic dishes in all cultures from patients and controls. The TRAP multinucleated cells (TRAP MuOCL) from patients were larger and stained less intensely for TRAP than those from controls (Fig. 2A,B). Accurate counting of TRAP MuOCL in patient cultures was difficult as osteoclasts in patient cultures were large and pale. The patient TRAP MuOCL formed peripheral belts of F-actin-rich podosomes similar to controls. The cir - cumference of these actin belts was greater in the patient osteoclasts than the control osteoclasts, reflecting the greater size of the patient cells (Fig. 2C). Expression of osteoclastic and osteoclastogenic genes in cells from patients. We next assessed if the expression of genes associated with osteoclast differentiation and function was altered in the patient oste- oclasts. Gene expression analyses were performed twice with cells from two patients (Pt2 and Pt5) at one time point during the culture, and with two different controls (C1–C4) in each analysis (Figs  3 and S2). The reason data are provided separately for different controls and different patients is that the logistics of the cell collections and experiments did not allow us to compare the same control with the same patient in each experiment. We Scientific Repo R ts | 7: 3012 | DOI:10.1038/s41598-017-02533-2 6 www.nature.com/scientificreports/ Figure 3. Gene expression in osteoclasts from patients (Pt2, Pt5), and controls (C1, C2, C3, C4). Fold induction of gene expression between cells treated with M-CSF (M) with and without RANKL (RL) for three days. All values are given as mean ± SEM (n = 4). *P ≤ 0.05, **P ≤ 0.005, ***P ≤ 0.001 between indicated groups. performed separate experiments comparing cells from different controls and from different patients. We repeated the experiments at different times to ensure reproducibility of our observations. e mRN Th A expression of ACP5 (encoding the enzyme TRAP) and CTSK (encoding cathepsin K), both robust markers of the osteoclastic phenotype, was strongly upregulated by RANKL in cells from both patients and con- trols (Fig. 3A,B). There was no difference between patients and controls in the degree of upregulation of NFATc1 (encoding nuclear factor of activated T-cells, the master osteoclastogenic transcription factor), (Fig. 3C). No die ff rences between patients and controls were found in RANKL-induced upregulation of TCIRG1 (encoding the proton pump subunit, Atp6i, essential in acidification process) (Fig.  3D), or the mRNA expression of ITGAV and ITGB3, (encoding the α and β subunits of the vitronectin receptor, essential for osteoclast attachment to v 3 bone) (Fig. 3E,F). The expression of CLCN7 (encoding chloride channel 7, with an essential role in acidification) was consistent throughout osteoclastogenesis in cells from both patients and controls (Fig. 3G). Although some genes showed a significant difference between patient and controls in one of the two independent experiments, none of the analyzed genes were significantly and reproducibly changed in both patients compared to controls. es Th e data indicate that the SNX10 mutation does not ae ff ct the expression of key genes associated with osteoclast differentiation and function. e Th ability of osteoclasts with the SNX10 mutation to resorb bone was analysed by culturing CD14 cells from patients on slices of devitalized bovine bone in six different experiments, each comparing cells from one of 3 patients with several of 7 controls. In control cultures, RANKL stimulated the formation of many TRAP MuOCL. In patient cultures we observed TRAP MuOCL, but, as in cultures on plastic, we also noticed several much larger cells with weak TRAP staining, making accurate counting of the number of osteoclasts impossible (Fig. 4A). We, therefore, also used TRAP5b released into the culture media as a parameter of mature osteoclast formation , but found no consistent difference between patients and controls (data not shown). Phalloidin staining demonstrated that control osteoclasts frequently exhibited crescent-shaped actin rings, which were associated with extensive resorption pits and resorption trails in the bone slices (Fig. 4B,C) and release of CTX (a degradation fragment from collagen type I) (Fig. 4D). Patient osteoclasts polarized their actin Scientific Repo R ts | 7: 3012 | DOI:10.1038/s41598-017-02533-2 7 www.nature.com/scientificreports/ Figure 4. Patient-derived osteoclasts form actin rings but are unable to resorb bone. (A) TRAP staining of CD14 cells cultured in M/RL on bone for eight days. Scale bar 40× panel: 500 μm, scale bar 100× light microscope and reflective light panel: 200 μm, scale bar 100× toluidine blue panel: 100 μm. (B,C) Osteoclasts were analyzed by confocal microscopy; (B) stained with TRITC-phalloidin (red), DAPI (blue) and the bone surface using a fluorescent bisphosphonate (AF-ALN; green). Resorption pits stain brightly with AF-ALN and can be seen in the control only. Scale bar = 100 μm; (C) stained with TRITC-phalloidin (red), DAPI (blue) AF- ALN (magenta) and the cell membranes with wheat germ agglutinin (green). AF-ALN omitted from upper panel of Pt2 for clarity. Black patches in C2 xy image indicates resorption pits. Bottom panel of Pt2 is a zoom of the zx section illustrating lack of resorption. Scale bar = 20 μm; (D) CTX-1 release in media during day 6 to 8 of culture in M/RL. C1, C6 = controls, Pt2, Pt5 = patients (n = 4–5). Scientific Repo R ts | 7: 3012 | DOI:10.1038/s41598-017-02533-2 8 www.nature.com/scientificreports/ cytoskeleton in rings at the periphery of the cells, similar to the actin belts seen in cells on plastic. Patient cultures, however, showed no evidence of resorption. er Th e were no resorption pits (Fig.  4A–C) and there was no release of CTX (Fig. 4D). There was no evidence of vacuolation in any of the patient osteoclasts (Fig.  4B). Increased size of patient osteoclasts is independent of precursor fusion. Individual patient oste- oclasts occupied larger areas of the bone surface than control osteoclasts (Fig. 4A–C), similar to the observations on plastic (Fig. 2A–C). However, osteoclasts on bone from patients contained similar numbers of nuclei as those from healthy controls (mean nuclei/osteoclast = 7.3 for C1, 7.8 for C2, 7.9 for Pt2). Patient osteoclasts were much flatter than the control osteoclasts; in cross-sectional views most osteoclasts were no deeper than a single nucleus (Fig. 4C). This indicates that the larger size of the patient osteoclasts is the result of increased spreading, rather than the result of increased fusion. In agreement with these findings, we were unable to observe any significant differences replicated in two separate experiments between expressions of several genes suggested to be involved in osteoclast fusion (Fig. S2). 17–19 RANKL-induced upregulation of the osteoclast fusion genes DC-STAMP and OC-STAMP was similar in cells from patients and controls (Fig. S2). Syncytin-1 and its receptor Asct2 have also been suggested to be involved in osteoclast fusion but the mRNA expression of the genes encoding these proteins (ERVW-1 and SLC1A5, respec- tively) was not different between patients and controls (Fig.  S2). Osteoclasts exhibiting SNX10 mutation have an impaired ruffled border. Since osteoclast forma- tion and F-actin ring formation was unae ff cted, we speculated that the lack of resorption by osteoclasts derived from the patients would be due to defective ruffled border formation, similar to well-characterized cases of osteoclast-rich ARO . This would be consistent with the likely role of SNX10 in vesicular trafficking processes 22, 23 in osteoclasts. While control osteoclasts on bone accumulated PNA-lectin (a marker of the ruffled border) within the actin ring, patient osteoclasts showed no PNA-lectin staining (Fig. 5A), suggesting impaired ruffled border formation in the patient osteoclasts. Transmission electron microscopy was used to further examine the sealing zone (F-actin rings) and the ruf- fled border in cultured osteoclasts. This ultrastructural analysis confirmed that patient osteoclasts retained the capacity to form sealing zones (Fig. 5E), but that 65% of patient osteoclasts lacked ruffled borders (Grade 0) and showed no evidence of bone resorption. In contrast, control osteoclasts showed evidence of extensive bone resorption and only 18% of these cells lacked ruffled border (Fig.  5B,D and F). Fifteen percent of patient osteo- clasts showed moderate/extensive ruffled borders, compared to 50% of control osteoclasts. The ultrastructural analysis confirmed that patient osteoclasts were very large and flat compared to control cells, with sealing zones located at the cell periphery. Gene expression of wild type and mutated SNX10. SNX10 mRNA levels were measured in RANKL-stimulated CD14 cells using two different TaqMan assays, spanning the junction between exon 1–2 and 4–5, respectively. SNX10 mRNA, analyzed with the exon 1–2 assay, was present in patients and controls and significantly upregulated by RANKL in controls and in one of the two patients (Fig.  6A). However, the mRNA levels were lower in patient cells. The assay spanning the exon 4–5 junction also showed upregulation of SNX10 by RANKL in control cells (Fig. 6B). However, as the mutation is located at the splice donor site of exon 4, this domain was barely detected in cells from patients. Discussion We herein report the clinical features and natural course of the Västerbotten form of IARO [MIM 615085], the disease-causing mutation, its origin and age. We also present functional analyses of the effects of this mutation on osteoclast differentiation and function. The Västerbotten form of IARO is caused by a splice site mutation in the SNX10 gene [MIM 614780] c.212 + 1 G > T resulting in a frame shift and a stop codon (p.S66Nfs * 15). The carrier frequency was found to be 1/93 in the population of Västerbotten. According to haplotype analysis, the shared ancestral origin of the SNX10 gene variant is approximately 950 years old. 1–3, 21, 24 e Th phenotype of IARO, previously associated with mutations in the PLEKHM1 , CLCN7, and CA2 genes , is defined as asymptomatic at birth, with patients exhibiting fractures by the end of their first decade of life, 24–27 while bone marrow failure and hepato-splenomegaly are rare . In contrast, our data show that patients with SNX10-related IARO of Västerbotten type described here have a more severe phenotype, with age at onset in early infancy, presenting with optic atrophy, macrocephaly, frontal bossing, proptosis, and anemia. Splenomegaly, due to compensatory extra-medullary hematopoiesis, was present before the age of five, and splenectomy was performed in three patients (Pt1, Pt8, and Pt9), the youngest at the age of five years (Pt1). Although the patients displayed many of the symptoms seen in malignant ARO, they all survived early childhood without being trans- planted. We, therefore, placed the patients in the group of IARO, rather than in the group of malignant ARO in which patients do not survive childhood if not transplanted. Osteopetrosis presenting with hydrocephalus, hind-brain posterior fossa crowding, foramen magnum nar- 28–31 rowing, and syringohydromyelia is rare, but has been described in the literature . Cure et al. studied cranial MRI of osteopetrosis, but did not identify ventriculomegaly, tonsillar herniation, or remodelled inner table in IARO . In the Västerbotten intermediate form of osteopetrosis these features seem to be common; in our study, four patients (Pt1, Pt5, Pt8, and Pt9) had cerebellar tonsillar descent through the foramen magnum, and three (Pt1, Pt5, and Pt8) of them showed remodelled inner table. Three subjects (Pt5, Pt8, and Pt9) had foramen mag- num stenosis, posteriorly angulated dens with stenosis and signs of brain stem compression. Structural brain malformation is rarely reported as part of the clinical spectrum of osteopetrosis, but cerebral atrophy, corpus callosum agenesis, hypoplasia of the hippocampus, and Dandy-Walker malformation has been Scientific Repo R ts | 7: 3012 | DOI:10.1038/s41598-017-02533-2 9 www.nature.com/scientificreports/ Figure 5. Patient-derived osteoclasts formed defective ruffled borders. (A) Phalloidin staining for acting rings (left) and lectin staining for ruffled border (middle) of CD14 cells cultured in M/RL on bone for 11 days. (B–E)-Representative TEM micrographs of control C2 (B,C) and patient Pt2 (D,E) -derived osteoclasts show extensive ruffled border (grade 3) and resorption in the control, compared to lack of ruffled border (grade 0) and resorption in the patient. Patient osteoclasts formed sealing zones (E). N = nucleus, SZ = sealing zone, RB = ruffled border. Scale bars in left panels = 2 μm. Scale bars in right panels = 500 nm. (F) The percentage of osteoclasts categorized into each ruffled border grade (control, n = 38; patient, n = 34). 33–36 reported . However, none of our SNX10 related IARO patients had intra-cerebral malformations. MRI of the brain showed calcification of the basal ganglia, and subcortical frontal lobes in one male patient (Pt8), who died due to spinal cord compression caused by a combination of prominent calvarial and skull base thickening, Scientific Repo R ts | 7: 3012 | DOI:10.1038/s41598-017-02533-2 10 www.nature.com/scientificreports/ Figure 6. Expression of SNX10 mRNA in cells treated with M or M/RL for three days, analyzed with an assay covering the exon 1–2 (A) and exon 4–5 (B) junction. Values are given as mean ± SD (n = 2–4). ***P ≤ 0.001 between M and M/RL treated cells. toncillar herniation, foramen magnum narrowing and posteriorly angulated dens with stenosis. He also had ven- triculomegaly, dural venous sinus stenosis, and remodelled inner table as a sign of elevated intra-cranial pressure. MRI also showed collateral veins between sinus transversus and extra-cranial veins, possibly caused by a jugular foramen stenosis, occlusion of the jugular vein and thrombosis of the transverse and sagittal sinuses, which is not previously described in subjects with IARO. However, no brain phenotype was reported in Snx10 knock out mice . Calcification of basal ganglia, dura mater, tentorium cerebelli and periventricular calcification extending into the white matter (Marble brain disorder) has previously been described with or without association to renal 38–41 tubular acidosis due to carbonic anhydrase II deficiency . The reason for calcification of the brain is unclear. Results of studies performed in Snx10 knock out mice suggest an additional role for Snx10 in the gastric epithelium. Snx10 knock out mice have impaired gastric acidification and are severely hypocalcaemic com- pared to wild type littermates, resulting in rickets as well as osteopetrosis (osteopetrorickets). In contrast, osteoclast-specic fi knockout of Snx10 resulted in osteopetrosis without rickets or alterations in calcium balance . Interestingly, the patients in the present study have no signs of rickets. Also, SNX10 osteopetrotic patients have been shown to be almost completely rescued from the sclerotic phenotype by hematopoietic stem cell transplan- tation . Considering the results of these studies, there seems to be a difference in the importance of SNX10 in different organs in mice and humans. Using RANKL-stimulated peripheral blood monocytes we found that osteoclasts could be generated from cells isolated from the patients bearing the SNX10 splice site mutation. This is consistent with previous stud- ies, which were also able to generate osteoclasts from peripheral blood of patients with different SNX10 muta - 10 37, 43 tions , and from spleen or bone marrow cells of Snx10 KO mice . In contrast, silencing of Snx10 in the mouse pre-osteoclastic cell line Raw 264.7 was shown to inhibit RANKL-induced osteoclast differentiation . Moreover, SNX10 mRNA levels were increased during RANKL-stimulated osteoclastogenesis, consistent with a likely role in osteoclast function rather than formation. Although we have not been able to obtain bone biopsies from the patients, these findings suggest that the patients have an osteoclast-rich form of osteopetrosis. e os Th teoclasts generated from patients had normal expression of key osteoclast genes including those asso- ciated with osteoclast fusion and bone attachment, however, the osteoclasts were larger than in healthy controls. e Th re was not die ff rence in the number of nuclei per osteoclasts between patients and controls, indicating that the larger size is the result of increased spreading, similar to what is seen in other cases of IARO . By contrast, Aker et Scientific Repo R ts | 7: 3012 | DOI:10.1038/s41598-017-02533-2 11 www.nature.com/scientificreports/ al. reported that osteoclasts generated from patients with a missense mutation (Arg51Gln) in SNX10 were fewer in number and smaller than in controls and, similar to osteoclasts from the Snx10 null mouse, had a markedly 10, 37, 43 decreased capacity to dissolve calcium phosphate . Our data show that osteoclasts cultured from patients were unable to resorb bone, evidenced by both an absence of resorption pits and lack of CTX in the culture super- natant. This finding is identical to earlier findings with osteoclasts cultured from IARO patients bearing muta- 21 10 tions in PLEKHM1 . Interestingly, human osteoclasts with a missense SNX10 mutation , and mouse osteoclasts 37, 43 from Snx10 null mice although having reduced activity, still had some capacity to dissolve calcium phosphate film, whereas human osteoclasts from our patients did not show any evidence of bone resorbing activity. This discrepancy most likely is due to the use of crystalline calcium phosphate coating instead of bone as substrate for the assessments of the bone resorbing activity. Electron microscopy revealed that the inability of osteoclasts from osteopetrotic patients in the Västerbotten County to resorb bone is associated with defective ruffled border formation, very similar to patients with ARO 4, 21, 44 caused by mutations in other genes . Using phalloidin staining, we found that the SNX10 patient osteoclasts were still able to form F-actin rich podosome belts in plastic dishes and normal F-actin rings when cultured on bone. In agreement with these observations, electron microscopy demonstrated these osteoclasts formed a seal- ing zones with normal appearance. These observations indicate that the defect lies in the endosomal/lysosomal trafficking pathways associated with ruffled border formation. This is perhaps not surprising, since SNX10 is a member of the sorting nexin family of proteins that play crucial roles in cargo sorting in the endosomal path- 45, 46 4 way , moreover; other genetic mutations that cause osteoclast-rich ARO also ae ff ct this pathway . In contrast to these findings, osteoclasts generated from bone marrow macrophages with deletion of the Snx10 gene have been reported to form an abnormal F-actin positive podosome belt in plastic dishes . We have no explanation for these differences between mouse and human osteoclasts. As discussed below, we do not believe that the difference might be due to expression of a truncated, functional form of SNX10 in our patients. er Th e seem to be clear differences in the functional consequences for osteoclasts in patients with the Arg51Gln mutated SNX10 and those in the present study. u Th s, the Arg51Gln mutation results in decreased osteoclast for - mation and decreased osteoclast activity, whereas the SNX10 mutation found in the Västerbotten form of IARO impairs resorption by disrupting ruffled border formation in mature osteoclasts. Furthermore, osteoclasts with the Arg51Gln mutation exhibit extensive cytoplasmic vacuolation, which was not seen in osteoclasts with the Västerbotten mutation. Increased vacuolation has also been observed by overexpressing SNX10 in several human cell lines , which raises the possibility that the Arg51Gln mutation results in gain of function; however, overex- pression of the Arg51Gln form of SNX10 does not cause vacuolation , therefore, the functional consequences of this mutation remain unclear. No full length SNX10 protein can be expressed in osteoclasts with the Västerbotten mutation, as the mutation leads to a stop codon upstream from intron 4. We do not know whether a truncated form of SNX10 is expressed but we have not been able to conclusively show the presence of a truncated SNX10 protein using four different antibodies in Western blots, including antibodies recognizing the amino-terminal part of SNX10. SNX proteins contain a conserved phox-homology (PX) domain, which can bind phosphorylated phosphatidylinositols (PI), thereby targeting SNXs to PI-enriched membranes. The PX domain extends from c.30–381 in exons 3–6. It seems as if the PX domain in SNX10 is even further extended (PXe) in the C-terminal 48, 49 part . Since PX or PXe is the functional domain in SNXs, it is unlikely, even if a truncated SNX10 protein is expressed in osteoclasts carrying the Västerbotten mutation, that this truncated protein has any biological activity. e diff Th erence between the Arg51Gln mutation and the Västerbotten mutation is most likely responsible for the differences seen in osteoclast formation and function. In conclusion, the Västerbotten intermediate form of osteopetrosis is caused by a splice site sequence variant in the SNX10 gene and differs in clinical characteristics from other intermediate types of ARO. The age at onset is early in infancy, optic atrophy and anemia is present at diagnosis and tonsillar herniation, foramen magnum stenosis, and severe osteomyelitis/osteonecrosis of the jaw are common. Osteoclasts generated from peripheral blood monocytes of the patients showed defects in ruffled border formation and failed to resorb bone. These data extend the clinical spectrum of IARO and demonstrate a critical role for SNX10 in ruffled border formation. Materials and Methods Subjects. We collected clinical data, medical records, x-ray, computed tomography scan (CT) and magnetic resonance imaging (MRI) of the brain, and blood samples from five male and four female patients with IARO. We also collected blood samples from their parents and siblings. Patients descended from the county of Västerbotten, northern Sweden. Exome sequencing. DNA from three patients and 15 non-related individuals (controls) from the county of Västerbotten were sequenced. Exome enrichment was performed using 3 μg genomic DNA. DNA samples were sheared by sonication with the Covaris S2 instrument (Covaris, Inc.). Fragment libraries were constructed from the sheared samples using the AB Library Builder System (Life Technologies) and target enrichment was performed according to the manufacturer’s protocols (Agilent SureSelect Human All Exon v3 kit). Exome capture was conducted by hybridizing the DNA libraries with biotinylated RNA baits for 24 h followed by extraction using streptavidin coated magnetic beads. Captured DNA was then amplified followed by emulsion PCR using the EZ Bead System (Life Technologies) and sequenced on the SOLiD5500xl system. Individual libraries were labelled by a post-hybridization barcoding procedure (Agilent, Santa Clara, CA, USA; barcodes compatible with SOLiD sequencing technology). Alignment and variant calling. Alignment of color space reads to the human reference genome (GRCh37/ hg19) was performed using v2.1 of the Lifescope Software (Life Technologies). Single nucleotide variants (SNVs), Scientific Repo R ts | 7: 3012 | DOI:10.1038/s41598-017-02533-2 12 www.nature.com/scientificreports/ and small insertions and deletions (indels) were subsequently identified by the diBayes algorithm available within the Lifescope software. All identified SNVs and indels were imported into a local installation of the CanvasDB database system , for annotation and further analysis of the variants. DNA and cDNA sequencing, haplotype analysis, carrier frequency. DNA was extracted from whole blood using a standard salting-out method. SNX10 gene amplification and bidirectional sequencing of exons and intron-exon boundaries were performed in samples from all patients and their parents. RNA was extracted from cultured patient peripheral blood monocytes and converted to cDNA as described below for gene expression analysis. Bidirectional Sanger sequencing (Applied Biosystems BigDye Terminator v3.1 Cycle Sequencing Kit, Applied Biosystems) on a 3730xl DNA Analyzer (Applied Biosystems) was used for analysis of SNX10 cDNA (NM_001199837.1) using primer in exon 3 to 5 to confirm the mutation and to identify the effect of the splicing. Sequence analysis was performed using the Sequencher software (Gene Codes Corporation, Ann Arbor, MI, USA). Primers sequences where designed using Primer 3 Plus software (primer3plus.com/) and are available upon request. Five patients (Pt1–5) were analyzed using the GeneChip Mapping 250 K array (Affymetrix, Santa Clara, CA), and one patient (Pt9) was analyzed using the HumanOmniExpressExome-8v1array (Illumina), for haplotype analysis, and estimation of the age of the sequence variant. To calculate the carrier frequency, DNA from 1000 randomly chosen individuals from the county of Västerbotten was used (The medical bio-bank, Umeå university hospital, Umeå, Sweden) . In vitro culture of CD14 monocytes from patients and controls. Peripheral blood mononuclear cells were isolated by Ficoll-Paque PLUS separation and CD14 cells were labelled with CD14 MicroBeads and extracted using a MACS column according to the manufacturer’s instructions (Miltenyi). Cells were then seeded 5 2 in 96-well plates (3 × 10 cells/cm ) in complete α-MEM medium (Gibco cat no 22561-021) supplemented with 10% heat inactivated fetal bovine serum (FBS, Sigma cat no F7524), 2 mM GlutaMAX (Gibco cat no 35050-038), 50 µ g/ml gentamicin (Gibco cat no 15750-037), 100 U/ml penicillin and 100 µ g/ml streptomycin (Gibco cat no 15140-148) and 30 ng/ml human M-CSF (M, R&D Systems cat no 216-MC-025/CF). For osteoclast generation, 2 ng/ml recombinant mouse RANKL (RL, R&D Systems cat no 462-TEC-010) was added from the start of the experiments (day 0). Media was replenished every three days and the cells were stained for TRAP at different days of culture using the Acid Phosphatase, Leukocyte (TRAP) Kit from Sigma (cat no 387A). TRAP cells containing three or more nuclei were counted as TRAP multinucleated osteoclast (MuOCL). Analysis of osteoclast differentiation, gene expression and bone resorption. RANKL-stimulated osteoclast differentiation of CD14 cells was analyzed on plastic and on devitalized bovine bone discs (IDS Immunodiagnostics Systems cat no TDT-1BON1000-96). The surface area of TRAP positive cells on plastic TM were analyzed using the Osteomeasure software (v3.2.1.0, Osteometrics Inc.). Total RNA was purified from cells at day 3 using RNeasy Micro Kit (Qiagen). cDNA synthesis was performed using the High Capacity cDNA Reverse Transcription Kit from Applied Biosystems. Gene expression was analyzed in two different experiments with in total four healthy controls and two patients using predesigned TaqMan Assays (Life Technologies) and the StepOnePlus Real-Time PCR system (Applied Biosystems). We analyzed ACP5 (Hs00356261_m1), CTSK (Hs01080388_m1), NFATc1 (Hs00542678_m1), TCIRG1 (ATP6i, HS00990751_m1), CLC7 (Hs01126462_ m1), ITGAV (Hs00233808_m1), ITGB3 (Hs01001469_m1), DC-STAMP (Hs00229255_m1), OC-STAMP (Hs00875776_m1), ERVW-1 (Syncytin-1, Hs02341206_g1), SLC1A5 (Asct2, Hs00194540_m1), SNX10 exon 2–3 (Hs01007224_m1) and SNX10 exon 4–5 (HS01007226_m1). The mRNA abundance of each gene was adjusted for the expression of 18S (Life Technologies 4310893E) used as an internal control. In most figures, gene expres- sion is presented as the fold increase between M/RL and M treated cells. Statistical analysis was performed using One-Way ANOVA with Tukey’s post hoc test. The bone resorptive activity of osteoclasts on bone slices was ana - lyzed by measuring the release of C-terminal telopeptides of type I collagen (CTX) in culture supernatants using a commercial ELISA (IDS Immunodiagnostics Systems cat no AC-07F1) and by Toluidine blue staining of positive excavations. TRAP5b was analyzed in culture media using commercial ELISA (IDS Immunodiagnostics Systems cat no SB-TR201A). Staining for the actin ring and ruffled border. Actin organization was studied at day 4 on plastic and day 11 on bone. Cells were fixed in 4% paraformaldehyde, washed with PBS and then permeabilised using 0.1% Triton X-100 in PBS for 10 min at 4 °C. After washing again with PBS, the cells were stained with rhodamine conjugated phalloidin (Life Technologies cat no R415), 5 U/ml in 2% BSA/PBS, for 20 min at 4 °C. Following washing in PBS the cells were mounted in Prolong Gold Mountant (Life Technologies cat no P-36931). Cells cultured on bone were stained for phalloidin as above but including 100 µ g/ml FITC conjugated peanut aggluti- 22, 23 nin (PNA)-lectin (Sigma cat no L7381) to stain the ruffled border . Photographs were taken using widefield uo fl rescence microscopy. Confocal microscopy. Osteoclasts cultured on bovine bone discs were prepared as above, and, in addition to phalloidin, were stained with 5 μg/ml Alexa Fluor 633-conjugated wheat germ agglutinin to stain cell membranes and 10 μM Alexa Fluor 488-conjugated alendronate (AF-ALN, ref. 51) to stain bone surfaces. Bone discs were mounted on to glass slides in VectaShield with DAPI (to stain nuclei) mounting medium before analyzing on a Zeiss LSM700 confocal microscope. To assess potential differences in fusion, the number of nuclei was counted in 80 randomly selected osteoclasts from two control cultures and one patient culture. Scientific Repo R ts | 7: 3012 | DOI:10.1038/s41598-017-02533-2 13 www.nature.com/scientificreports/ Transmission Electron Microscopy investigation of ruffled border morphology. To investigate whether patient osteoclasts were capable of forming mature ruffled borders, cells cultured on bovine bone were analyzed by transmission electron microscopy (TEM). Briey fl , cell cultured with M-CSF and RANKL for eight days were fixed in 2.5% glutaraldehyde in 0.1 M cacodylate buffer. The bone discs were then demineralized in 0.15 M EDTA before being postfixed in 1% osmium tetroxide and embedded in Epon resin. Bone discs were sectioned perpendicular to the bone surface to obtain cross sections of osteoclasts and their ruffled border areas. Ultrathin sections were imaged using a JEOL Jem-1400Plus TEM equipped with an AMT UltraVUE camera. Following image acquisition, osteoclasts were graded based on their ruffled border morphological characteristics, similarly to a previous study, ref. 52: grade 0 - no ruffled border apparent; grade 1 - rudimentary ruffled border, with few poorly-developed membrane folds apparent; grade 2 - moderate ruffled border, with folded membrane apparent at some areas; grade 3 - extensive, highly-folded ruffled border. Statistics. For experiments with two patients and two controls, statistical analysis was performed using One-Way ANOVA with Tukey’s post hoc test. When one patient was compared to one control, a two-sided Student’s t test was used. A die ff rence was considered statistically significant with P < 0.05. Due to the limited number of patients and the large inter-patient variation, it was not possible to pool the results and make statistical comparisons based on the patients as one group. Instead, the results are shown for each patient separately. For gene expression analyses (Figs 3 and S1) each independent experiment is shown. Study approval. e R Th egional Ethical Review Board at Umeå University, Umeå, Sweden, approved the study. Patients, parents or first-degree relatives gave their written informed consent. Methods using human tissue were in accordance with the Helsinki declaration. References 1. Frattini, A. et al. Defects in TCIRG1 subunit of the vacuolar proton pump are responsible for a subset of human autosomal recessive osteopetrosis. Nature genetics 25, 343–346, doi:10.1038/77131 (2000). 2. Kornak, U. et al. Loss of the ClC-7 chloride channel leads to osteopetrosis in mice and man. Cell 104, 205–215, doi:10.1016/S0092- 8674(01)00206-9 (2001). 3. Sobacchi, C. et al. The mutational spectrum of human malignant autosomal recessive osteopetrosis. Human molecular genetics 10, 1767–1773, doi:10.1093/hmg/10.17.1767 (2001). 4. Sobacchi, C., Schulz, A., Coxon, F. P., Villa, A. & Helfrich, M. H. Osteopetrosis: genetics, treatment and new insights into osteoclast function. Nature reviews. Endocrinology 9, 522–536, doi:10.1038/nrendo.2013.137 (2013). 5. Villa, A., Guerrini, M. M., Cassani, B., Pangrazio, A. & Sobacchi, C. Infantile malignant, autosomal recessive osteopetrosis: the rich and the poor. Calcie fi d tissue international 84, 1–12, doi:10.1007/s00223-008-9196-4 (2009). 6. Askmyr, M. K., Fasth, A. & Richter, J. Towards a better understanding and new therapeutics of osteopetrosis. British journal of haematology 140, 597–609, doi:10.1111/j.1365-2141.2008.06983.x (2008). 7. Fasth, A. Osteopetrosis–more than only a disease of the bone. American journal of hematology 84, 469–470, doi:10.1002/ajh.21454 (2009). 8. Bliznetz, E. A. et al. Genetic analysis of autosomal recessive osteopetrosis in Chuvashiya: the unique splice site mutation in TCIRG1 gene spread by the founder effect. European journal of human genetics: EJHG 17, 664–672, doi:10.1038/ejhg.2008.234 (2009). 9. Phadke, S. R. et al. Malignant recessive osteopetrosis. Indian pediatrics 36, 69–74 (1999). 10. Aker, M. et al. An SNX10 mutation causes malignant osteopetrosis of infancy. Journal of medical genetics 49, 221–226, doi:10.1136/ jmedgenet-2011-100520 (2012). 11. Megarbane, A. et al. Homozygous stop mutation in the SNX10 gene in a consanguineous Iraqi boy with osteopetrosis and corpus callosum hypoplasia. European journal of medical genetics 56, 32–35, doi:10.1016/j.ejmg.2012.10.010 (2013). 12. Pangrazio, A. et al. SNX10 mutations define a subgroup of human autosomal recessive osteopetrosis with variable clinical severity. Journal of bone and mineral research: the official journal of the American Society for Bone and Mineral Research 28, 1041–1049, doi:10.1002/jbmr.1849 (2013). 13. Ameur, A., Bunikis, I., Enroth, S. & Gyllensten, U. CanvasDB: a local database infrastructure for analysis of targeted- and whole genome re-sequencing projects. Database: the journal of biological databases and curation 2014, bau098–bau098, doi:10.1093/ database/bau098 (2014). 14. Genomes Project, C. et al. A global reference for human genetic variation. Nature 526, 68–74, doi:10.1038/nature15393 (2015). 15. Hastbacka, J. et al. Linkage disequilibrium mapping in isolated founder populations: diastrophic dysplasia in Finland. Nature genetics 2, 204–211, doi:10.1038/ng1192-204 (1992). 16. Alatalo, S. L., Halleen, J. M., Hentunen, T. A., Monkkonen, J. & Vaananen, H. K. Rapid screening method for osteoclast differentiation in vitro that measures tartrate-resistant acid phosphatase 5b activity secreted into the culture medium. Clinical chemistry 46, 1751–1754 (2000). 17. Kukita, T. et al. RANKL-induced DC-STAMP is essential for osteoclastogenesis. e J Th ournal of experimental medicine 200, 941–946, doi:10.1084/jem.20040518 (2004). 18. Yagi, M. et al. DC-STAMP is essential for cell-cell fusion in osteoclasts and foreign body giant cells. The Journal of experimental medicine 202, 345–351, doi:10.1084/jem.20050645 (2005). 19. Yang, M., Birnbaum, M. J., Mackay, C. A., Mason-Savas, A., Thompson, B. & Odgren, P. R. Osteoclast stimulatory transmembrane protein (OC-STAMP), a novel protein induced by RANKL that promotes osteoclast differentiation. Journal of Cellular Physiology 215, 497–505, doi:10.1002/jcp.21331 (2007).  20. Soe, K. et al. Involvement of human endogenous retroviral syncytin-1 in human osteoclast fusion. Bone 48, 837–846, doi:10.1016/j. bone.2010.11.011 (2011). 21. Van Wesenbeeck, L. et al. Involvement of PLEKHM1 in osteoclastic vesicular transport and osteopetrosis in incisors absent rats and humans. e J Th ournal of clinical investigation 117, 919–930, doi:10.1172/JCI30328 (2007). 22. Takagi, M., Yagasaki, H., Baba, T. & Baba, H. Ultrastructural visualization of selective peanut agglutinin binding sites in rat osteoclasts. The journal of histochemistry and cytochemistry: official journal of the Histochemistry Society 36, 95–101, doi:10.1177/36.1.2447153 (1988). 23. Salo, J., Metsikko, K., Palokangas, H., Lehenkari, P. & Vaananen, H. K. Bone-resorbing osteoclasts reveal a dynamic division of basal plasma membrane into two different domains. Journal of cell science 109(Pt 2), 301–307 (1996). 24. Palagano, E. et al. Buried in the Middle, But Guilty: Intronic Mutations in the TCIRG1 Gene Cause Human Autosomal Recessive Osteopetrosis. Journal of bone and mineral research: the official journal of the American Society for Bone and Mineral Research, doi:10.1002/jbmr.2517 (2015). Scientific Repo R ts | 7: 3012 | DOI:10.1038/s41598-017-02533-2 14 www.nature.com/scientificreports/ 25. Kahler, S. G., Burns, J. A. & Aylsworth, A. S. A mild autosomal recessive form of osteopetrosis. American journal of medical genetics 17, 451–464, doi:10.1002/ajmg.1320170208 (1984). 26. Whyte, M. P. Primer on the Metabolic Bone Diseases and Disorders od Mineral Metabolism. 4th edn, 367–383 (Lippincott Williams & Wilkins, 1999). 27. Campos-Xavier, A. B., Saraiva, J. M., Ribeiro, L. M., Munnich, A. & Cormier-Daire, V. Chloride channel 7 (CLCN7) gene mutations in intermediate autosomal recessive osteopetrosis. Human genetics 112, 186–189, doi:10.1007/s00439-002-0861-9 (2003). 28. McCleary, L., Rovit, R. L. & Murali, R. Case report: myelopathy secondary to congenital osteopetrosis of the cervical spine. Neurosurgery 20, 487–489, doi:10.1227/00006123-198703000-00025 (1987). 29. Elster, A. D., Theros, E. G., Key, L. L. & Chen, M. Y. Cranial imaging in autosomal recessive osteopetrosis. Part I. Facial bones and calvarium. Radiology 183, 129–135, doi:10.1148/radiology.183.1.1549658 (1992). 30. Sari, A. & Demirci, A. Radiographic type I autosomal dominant osteopetrosis with syringohydromyelia. Neuroradiology 38, 532–533, doi:10.1007/BF00626090 (1996). 31. Al-Tamimi, Y. Z., Tyagi, A. K., Chumas, P. D. & Crimmins, D. W. Patients with autosomal-recessive osteopetrosis presenting with hydrocephalus and hindbrain posterior fossa crowding. Journal of neurosurgery. Pediatrics 1, 103–106, doi:10.3171/PED-08/01/103 (2008). 32. Cure, J. K., Key, L. L., Goltra, D. D. & VanTassel, P. Cranial MR imaging of osteopetrosis. AJNR. American journal of neuroradiology 21, 1110–1115 (2000). 33. Stark, Z., Pangrazio, A., McGillivray, G. & Fink, A. M. Association of severe autosomal recessive osteopetrosis and structural brain abnormalities: a case report and review of the literature. European journal of medical genetics 56, 36–38, doi:10.1016/j. ejmg.2012.10.001 (2013). 34. Rees, H., Ang, L. C., Casey, R. & George, D. H. Association of infantile neuroaxonal dystrophy and osteopetrosis: a rare autosomal recessive disorder. Pediatric neurosurgery 22, 321–327, doi:10.1159/000120923 (1995). 35. Ben Hamouda, H. et al. Association of severe autosomal recessive osteopetrosis and Dandy-Walker syndrome with agenesis of the corpus callosum. Acta orthopaedica Belgica 67, 528–532 (2001). 36. Abinun, M. & Pieniazek, P. Successful haematopoietic stem cell transplantation for osteopetrosis due to TCRIG1 mutation. Archives of disease in childhood 95, 984–984, doi:10.1136/adc.2010.194928 (2010). 37. Ye, L. et al. Osteopetrorickets due to Snx10 deficiency in mice results from both failed osteoclast activity and loss of gastric acid- dependent calcium absorption. PLoS genetics 11, e1005057, doi:10.1371/journal.pgen.1005057 (2015). 38. Whyte, M. P. et al. Osteopetrosis, renal tubular acidosis and basal ganglia calcification in three sisters. The American journal of medicine 69, 64–74, doi:10.1016/0002-9343(80)90501-X (1980). 39. Bourke, E., Delaney, V. B., Mosawi, M., Reavey, P. & Weston, M. Renal tubular acidosis and osteopetrosis in siblings. Nephron 28, 268–272, doi:10.1159/000182216 (1981). 40. Ohlsson, A., Stark, G. & Sakati, N. Marble brain disease: recessive osteopetrosis, renal tubular acidosis and cerebral calcification in three Saudi Arabian families. Developmental medicine and child neurology 22, 72–84, doi:10.1111/j.1469-8749.1980.tb04307.x (1980). 41. Patel, P. J., Kolawole, T. M., al-Mofada, S., Malabarey, T. M. & Hulailah, A. Osteopetrosis: brain ultrasound and computed tomography findings. European journal of pediatrics 151, 827–828, doi:10.1007/BF01957934 (1992). 42. Zhu, C. H., Morse, L. R. & Battaglino, R. A. SNX10 is required for osteoclast formation and resorption activity. Journal of cellular biochemistry 113, 1608–1615, doi:10.1002/jcb.24029 (2012). 43. Zhou, C. et al. Deficiency of sorting nexin 10 prevents bone erosion in collagen-induced mouse arthritis through promoting NFATc1 degradation. Ann Rheum Dis 75, 1211–1218, doi:10.1136/annrheumdis-2014-207134 (2016). 44. Bruder, E., Stallmach, T., Peier, K., Superti-Furga, A. & Vezzoni, P. Osteoclast morphology in autosomal recessive malignant osteopetrosis due to a TCIRG1 gene mutation. Pediatric pathology & molecular medicine 22, 3–9 (2003). 45. Teasdale, R. D. & Collins, B. M. Insights into the PX (phox-homology) domain and SNX (sorting nexin) protein families: structures, functions and roles in disease. Biochem J 441, 39–59, doi:10.1042/BJ20111226 (2012). 46. Carlton, J., Bujny, M., Rutherford, A. & Cullen, P. Sorting nexins–unifying trends and new perspectives. Traffic 6, 75–82, doi:10.1111/j.1600-0854.2005.00260.x (2005). 47. Qin, B., He, M., Chen, X. & Pei, D. Sorting nexin 10 induces giant vacuoles in mammalian cells. e J Th ournal of biological chemistry 281, 36891–36896, doi:10.1074/jbc.M608884200 (2006). 48. Xu, T. et al. Structure of human SNX10 reveals insights into its role in human autosomal recessive osteopetrosis. Proteins 82, 3483–3489, doi:10.1002/prot.24689 (2014). 49. Xu, J. et al. Structure of sorting nexin 11 (SNX11) reveals a novel extended phox homology (PX) domain critical for inhibition of SNX10-induced vacuolation. e J Th ournal of biological chemistry 288, 16598–16605, doi:10.1074/jbc.M112.449306 (2013). 50. Hallmans, G. et al. Cardiovascular disease and diabetes in the Northern Sweden Health and Disease Study Cohort - evaluation of risk factors and their interactions. Scandinavian journal of public health. Supplement 61, 18–24, doi:10.1080/14034950310001432 (2003). 51. o Th mpson, K., Rogers, M. J., Coxon, F. P. & Crockett, J. C. Cytosolic entry of bisphosphonate drugs requires acidification of vesicles aer fl ft uid-phase endocytosis. Molecular pharmacology 69, 1624–1632, doi:10.1124/mol.105.020776 (2006). 52. Taylor, A. et al. Impaired prenylation of Rab GTPases in the gunmetal mouse causes defects in bone cell function. Small GTPases 2, 131–142, doi:10.4161/sgtp.2.3.16488 (2011). Acknowledgements We sincerely thank the patients and family members who participated in this study. We would also like to thank Stefan Esher, Umeå University, for help with genealogy, and Anna Westerlund for excellent technical assistance. This work was supported by grants from the FOU, at the Umeå university hospital, and the Medical Faculty at Umeå University. The work at University of Gothenburg was supported by grants from The Swedish Research Council, the Swedish Rheumatism Association, the Royal 80-Year Fund of King Gustav V, ALF/LUA research grant from Sahlgrenska University Hospital in Gothenburg and the Lundberg Foundation. The work at the University of Gothenburg and the University of Aberdeen was supported by Euroclast, a Marie Curie FP7-People- 2013-ITN: # 607446. Author Contributions E.L.S., P.H., J.K., N.D. and U.H.L. designed the study. E.L.S., P.H., J.K., E.M., C.S.B., P.E.S., T.L., M.H.H., F.C., G.H., J.W., N.D. and U.H.L. acquired the data. E.L.S., T.L. and J.W. re-evaluated radiological examinations. J.K. and A.A. performed exome sequencing, Sanger sequencing and haplotype analysis. A.A., T.G.K. and P.R. filtered exome data. P.H., E.M., M.H.H., F.C. and U.H.L. conducted experiments. E.L.S., P.H., J.K., E.M., C.S.B., T.L., A.A., M.H.H., F.C., N.D., J.W. and U.H.L. wrote the manuscript. Scientific Repo R ts | 7: 3012 | DOI:10.1038/s41598-017-02533-2 15 www.nature.com/scientificreports/ Additional Information Supplementary information accompanies this paper at doi:10.1038/s41598-017-02533-2 Competing Interests: The authors declare that they have no competing interests. Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Cre- ative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not per- mitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/. © The Author(s) 2017 Scientific Repo R ts | 7: 3012 | DOI:10.1038/s41598-017-02533-2 16

Journal

Scientific Reports – Springer Journals

Published: Jun 7, 2017

There are no references for this article.