Patients
Patients affected by NF1 and followed at the Department of Translational Medicine, Federico II University of Naples, Pediatric Section, were enrolled in the study after the protocol was discussed with each patient (or legal tutor) and informed consent was obtained. Patients’ clinical data were obtained from medical records over the past 20 years.
A multidisciplinary team, including endocrinologists, geneticists and radiologists, evaluated all patients.
The enrollment was carried out according to the following inclusion criteria: (i) clinical diagnosis of NF1 based on recommendations of National Institute of Health [2]; (ii) informed consent expression. The exclusion criteria were: (i) concurrent disorders affecting bone metabolism, i.e., primary hyperparathyroidism, hyperthyroidism, renal or hepatic chronic disease, malabsorption, hypercortisolism, hypogonadism; (ii) previous or ongoing treatments with glucocorticoids, chemotherapies, bisphosphonates or other antiosteoporotic drugs; (iii) previous spine or hip fractures; (iv) surgical correction of scoliosis.
Collected clinical information included familiarity, presence or absence of CALMs, intertriginous SF, LN, cutaneous and subcutaneous NFs, plexiform NFs, spinal NFs, bone abnormalities, OPGs, cardiovascular malformations, endocrine system involvement, developmental delay/intellectual disability, cerebrovascular malformations, and occurrence of other neoplasms.
Patients were divided into three groups according to the severity of the phenotype using the classification proposed by Riccardi [20]. Patients presenting with CALMs, axillary freckling, LN, dermal and/or nodular NFs, and non-progressive scoliosis were classified as “mild”, those presenting with plexiform NFs, bone abnormalities, precocious or progressive scoliosis were classified as “moderate”, and patients with learning disability/mental retardation, OPG and/or other neoplasms, and/or cerebrovascular disease were classified as “severe”.
114 individuals diagnosed with NF1 according to National Institutes of Health criteria (2) were enrolled into the study. 69 patients were females and 44 males. The average age at time of diagnosis was 2.7 years (range 0.3–5 years), wherease the average age at observation was 11.9 years (range 1–23 years). 40 patients were children (aged between 1 and 11 years), 34 were in pubertal age (aged between 12 and 16) and 40 were adults (aged between 17 and 23).
Based on the phenotype, 31 patients were classified as “mild” (mean age 12.8 ± 6, range 2–22 years), 39 as “moderate” (mean age 12.3 ± 5.9, range 2–22 years), and 44 as “severe” (mean age 12.3 ± 5.5, range 1–23 years). Demographic and clinical characteristics are reported in Table 1.
All patients were screened for NF1 and SPRED1 mutations by parallel sequencing of the whole coding region and flanking splice sites (± 10 bp). Structural rearrangements were assessed by Multiplex Ligation Probe Amplification analysis using the MRC-Holland P295 probe set. A pathogenic or likely pathogenic NF1 variant was found in 90/114 (78.9%).
114 age- and sex-matched healthy controls were also enrolled. To minimize potential source of bias, none of the controls had disorders or treatments affecting bone metabolism or BMD.
To evaluate the efficacy of vitamin D therapy in patients with vitamin D insufficiency, 57 patients were followed for a mean period of 2,4 ± 0.9 years (range 2–5 years) after starting 25OH- D supplementation.
Methods
All methods were carried out in accordance with relevant guidelines and regulations. All experimental protocols were approved by “Comitato Etico Università Federico II”, protocol number 315/18
Clinical evaluation
Short stature was defined as height less than two standard deviations. Pubertal stage was assessed according to Tanner's pubertal stages; growth velocity was also evaluated.
Life style
In order to assess the potential contribution of “environmental” factors, physical activity, sun exposure, vitamin D and calcium intake were recordedduring the enrollment visit.
Physical activity was evaluated by administering the International Physical Activity Questionnaire (IPAQ) [23]. For each patient the time spent experiencing different kind of activities during the previous week was calculated: vigorous-intensity activity (hard physical effort and the patient breathe harder than normal), moderate-intensity activity (moderate physical effort, walking) and the time spent sitting.
Sun exposure was considered low with at least two of the following criteria: no arm and skin exposure during summer months; no sunbathing or holiday in sunny places; no working outdoors. Otherwise, it was deemed sufficient.
Diet calcium intake was defined low with at least two of the following criteria: milk assumption less than 100 mL a day; less than three yogurts in a week; eating cheese less than two times in a week. Otherwise, it was deemed sufficient.
Diet vitamin D intake was evaluated through a questionaire evaluating intake of vitamin D-rich-food as portion per week. These specific food included vitamin D enriched cereals, cheese, fish, lentils. Diet Vitamin D intake was defined low with less than 3 portions assumed in a week. Otherwise it was deemed sufficient.
Biochemical markers of bone metabolism
Bone metabolism was studied evaluating serum calcium, phosphorus, alkaline phosphatase, parathyroid hormone (PTH), calcitonin, 25OHD, C-terminal telopeptide of type I collagen (CTX) and osteocalcin (OC) and urinary calcium/creatinine ratio (UCa/UCr). CTX was measured using a serum cross-laps enzyme-linked immunosorbent assay (ELISA) kit (Immunodiagnosticsystems, Frankfurt, Germany) and OC was measured using a Microvue Osteocalcin ELISA kit (Quidel corporation, San Diego, CA, USA). ELISA tests were automated on a DSX system (Dynex technologies-Technogenetics-Milan-Italy) following the manufacturer’s instructions. The other markers of bone metabolism were measured using immunoassay with commercially available kits.
Considering the vitamin D levels, NF1 patients were classified according to the Endocrine Society’s latest guidelines for vitamin D levels [21]: deficiency was defined by 25OHD less than or equal to 20 ng/ml; insufficiency when it was between 21 and 29 mg/ml; normal if it was more than 30 ng/ml.
Whenever the presence of vitamin D deficiency was recorded, vitamin D supplementation was started and BMD data obtained during vitamin supplementation were compared with previous BMD data. Cholecalciferol was used at a dose of 2,000 UI /day in prepubertal patients and 5,000 UI/day in adolescent.
Hormonal studies
To address the impact of endocrine regulation on bone homeostasis, the different endocrine axes were evaluated. The somatotropic axis was evaluated by analyzing basal serum growth hormone (GH) and insulin-like growth factor 1 (IGF-1) via immune assay (CLIA); the thyrotropic axis function by analyzing serum thyroid-stimulating hormone (TSH), free triiodothyronine (fT3), free thyroxine (fT4), T3, T4; the adrenocorticotropic axis function by analyzing plasma adrenocorticotropic hormone (ACTH), serum cortisol, androstenedione, 17hydroxyprogesteron (17OHP), dehydroepiandrosterone sulphate (DHEA-S), renin, aldosterone levels and 24-hour urinary free cortisol (UFC); the gonadotropic axis function by analyzing serum follicle-stimulating hormone (FSH), luteinizing hormone (LH), 17β-estradiol, testosterone levels. The beta-cell function was analyzed by evaluating basal serum insulin levels, all measured by using immunoassay with commercially available kits.
Bone mineral density
BMD was studied using Dual-emission X-ray absorptiometry (DXA). DXA (Hologic QDR 1000, Hologic Inc., Waltham, USA). Since both children and adult participants were enrolled, BMD was measured at the L1-L4 vertebrae, considering that the hip is not a reliable site for measurement in growing children. Z-scores were calculated by comparing BMD with age and sex matched reference values according to the manufacturer’s internal reference database.
Osteopenia was defined as the presence of BMD between − 1.5 and − 2 SD while Osteoporosis was defined as BMD below − 2SD.
Statistical analysis
All data mentioned in the text or shown in the figure are expressed as mean ± standard deviation (SD). Sample size calculated using the following formula (with respect to the variable BMD): (r + 1/r) x [SD2 (Zβ-Zα/2)2]/d2 was 108 (r = 1, SD and d were retrieved from [26], Zβwas calculated for 80% power, Zα/2 was calculated at p < 0.05) Statistical analysis was performed using Statistical Package for Social Science (SPSS 18) for Windows Update; SPSS Inc., Chicago, Illinois, USA). The comparisons between numerical variables were performed by Student’s t test corrected for Fisher’s exact test. Pearson’s correlation test was performed to assess the relationship among variables with normal distribution (p < 0.05 was considered as statistically significant) whereas Spearman’s correlation test was performed to assess the relationship among variables with skewed distribution (a ρ < 0.005 was considered statistically significant). Logistic regression analysis was performed to determine independent predictors of patient outcome (between bone mineral marker and BMD). Univariate analysis of variance and covariance analysis were performed to eliminate the effects of age and gender.