# Unraveling the complexity of skeletal dysplasias in the national health system

Dorra Najjar, Aleš Maver, Ana Peterlin, Helena Jaklič, Borut Peterlin

► To cite this version:

Dorra Najjar, Aleš Maver, Ana Peterlin, Helena Jaklič, Borut Peterlin. Unraveling the complexity of skeletal dysplasias in the national health system. Frontiers in Endocrinology, 2025, 16, pp.1523737. 10.3389/fendo.2025.1523737. pasteur-05024270

HAL Id: pasteur-05024270

https://riip.hal.science/pasteur-05024270v1

Submitted on 8 Apr 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for updates

OPEN ACCESS

EDITED BY

Gretl Hendrickx,
KU Leuven, Belgium

REVIEWED BY

Emanuele Micaglio,
IRCCS San Donato Polyclinic, Italy
Noor Muhammad,
Kohat University of Science and
Technology, Pakistan

*CORRESPONDENCE

Borut Peterlin

<figure><img src="image_3.png" /></figure>

borut.peterlin@guest.arnes.si

†These authors have contributed
equally to this work and share
first authorship

RECEIVED 06 November 2024

ACCEPTED 19 February 2025

PUBLISHED 10 March 2025

CITATION

Najjar D, Maver A, Peterlin A, Jaklič H and
Peterlin B (2025) Unraveling the
complexity of skeletal dysplasias
in the national health system.
*Front. Endocrinol.* **16**:1523737.
doi: 10.3389/fendo.2025.1523737

COPYRIGHT

© 2025 Najjar, Maver, Peterlin, Jaklič and
Peterlin. This is an open-access article
distributed under the terms of the *Creative
Commons Attribution License (CC BY)*. The
use, distribution or reproduction in other
forums is permitted, provided the original
author(s) and the copyright owner(s) are
credited and that the original publication in
this journal is cited, in accordance with
accepted academic practice. No use,
distribution or reproduction is permitted
which does not comply with these terms.

Unraveling the complexity of
skeletal dysplasias in the
national health system

Dorra Najjar<fnref n="1" />,<fnref n="2" />†, Aleš Maver<fnref n="1" />†, Ana Peterlin<fnref n="1" />, Helena Jaklič<fnref n="1" /> and Borut Peterlin<fnref n="1" />*

**Introduction:** Skeletal dysplasia (SD) is a large and heterogeneous group of rare genetic disorders that affects bone and cartilage growth. These disorders are diagnosed based on radiographic, clinical, and molecular criteria. However, the diagnostics is challenging due to clinical and genetic heterogeneity. We present the experience of systematic use of comprehensive genetic testing in the national health system and the molecular epidemiology of SD in Slovenia.

**Methods:** We retrospectively reviewed 470 patients with clinical features of SD, including prenatal, childhood, and adult patients referred for diagnostic genetic evaluation to the national genetic reference center over ten years. In 262 patients, whole exome or whole genome sequencing was performed, while direct gene sequencing was performed in 208 patients with a specific clinical diagnosis.

**Results:** A definitive genetic diagnosis using NGS was achieved in 50% (n=131) of patients. Among the positive cases, 49.61% initially presented with a nonspecific diagnosis of SD, and genetic testing contributed to establishing the diagnosis. Moreover, we demonstrated high genetic heterogeneity in our SD cohort with 66 distinct causative genes, resulting in different types of SD. In detail, we detected 132 causative variants, of which 29 were novel, which expanded the mutational spectrum of SD. Furthermore, pathogenic copy number variants (CNVs) were identified in 4.55% of the total number of variants, highlighting the importance of CNV analysis in expanding the yield of molecular diagnosis of SD.

**Conclusion:** With the systematic use of WES and WGS, we have significantly improved the diagnostic yield of SD in the national health system and access to genetic testing. Moreover, we found significant genetic heterogeneity, and we report the genetic epidemiology of SD in the Slovenian population.

KEYWORDS

CNV, diagnostic yield, molecular pathology, NGS, prenatal diagnosis, rare genetic diseases, skeletal dysplasia. skeletal dysplasia

<footnote n="1">¹Clinical Institute of Genomic Medicine, University Medical Centre Ljubljana, Ljubljana, Slovenia,</footnote>

<footnote n="2">²Laboratory of Biomedical Genomics and Oncogenetics, LR16IPT05, Pasteur Institute of Tunisia, Tunis
El Manar University, El Manar I, Tunis, Tunisia</footnote>
Introduction

Bone tissue is a large compartment of the human body that
implicates a tightly regulated process for developing, growing, and
maintaining the human skeleton (cartilage and bone). SD is a
clinically and genetically heterogeneous group of rare disorders
that affects the process of bone and cartilage formation and
homeostasis (1). Despite the individual rarity of these diseases,
they present 5% of rare disorders (2), with a prevalence of at least 1/
5000 births (3). SD has significant consequences on patients' quality
of life due to its high morbidity, and the most severe forms can
reduce life expectancy (4, 5).

Recent advances in bone and cartilage metabolism have clarified the molecular causes of approximately 771 rare skeletal disorders, which exhibit highly complex phenotypes (6, 7). However, diagnosing this large group of disorders is challenging (8–10). Firstly, clinicians have limited knowledge of these rare disorders. Secondly, SDs exhibit significant genetic heterogeneity and pleiotropy. Thirdly, distinctive skeletal manifestations either do not appear until skeletal maturity or tend to diminish over time. Finally, some disease categories are still not well characterized. Moreover, these conditions often exhibit overlapping phenotypes, necessitating a multidisciplinary approach involving comprehensive evaluations by specialists across various fields. This approach combines diverse radiographic imaging techniques, clinical assessments, and, in some cases, histological evaluation from bone biopsies, resulting in a significant socio-economic burden on the healthcare system (11, 12).

A molecular diagnosis enables identification of the specific
cause of the underlying SD; however, using hypothesis-based
diagnostic approaches makes it challenging to determine which
gene(s) should be investigated (7, 13). Additionally, the modes of
inheritance are variable (including autosomal dominant, autosomal
recessive, X-linked, and mitochondrial inheritance), and the
mutational spectrum may encompass both CNVs and SNVs. This
variability necessitates multiple tests, which are both time- and cost-
intensive (7, 14). Presently, genetic testing is considered an essential
element in the workup of patients with congenital SD because of the
impact of these malformations on the quality of life of the patient
and the risk of the lethality of some disorders. Genetic testing also
facilitates primary prevention, family planning, and access to
therapies or novel clinical trials for certain conditions with known
genetic etiologies.

Recently, next-generation sequencing (NGS) has significantly
advanced our understanding of the molecular pathology of SD,
identifying over 552 genes involved in ossification, chondrogenesis,
and bone metabolism (6). Nevertheless, translating novel genomic
technologies into clinical practice has been slow. In this study, we
present the genetic diagnostic strategy and evaluate the impact of
systematically applying NGS to diagnose and gain insights into the
molecular pathology and genetic epidemiology of SD in the
Slovenian population.

Materials and methods

Patients

We conducted a retrospective study of 470 patients with a
clinical diagnosis of SD, including prenatal, pediatric, and adult
cases, who were referred for genetic diagnostics to the Clinical
Institute of Medical Genetics at the University Medical Center
Ljubljana between January 2013 and December 2022. All patients
received genetic counseling and were referred by a clinical geneticist
for genetic testing. We utilized NGS for genetic diagnosis in 262
patients. All patients underwent exome sequencing (ES) as the
primary diagnostic approach, while for some families, trio design
involving the parents, was performed for urgent cases, such as
prenatal diagnoses, and Whole Genome Sequencing (WGS) was
selectively performed for families with a strong suspicion of a
genetic etiology, particularly in familial cases. For patients
referred with a well-established diagnosis or a family history for
specific diseases and variants, we used Sanger sequencing to identify
common variants in specific genes, such as *FGFR3* (OMIM:
\*134934) for thanatophoric dysplasia type 1 (OMIM: \*187600) or
hypo/achondroplasia (OMIM: \*100800 and OMIM: \*146000), and
*COL1A1* (OMIM: \*120150) and *COL1A2* (OMIM: \*120160) genes
for osteogenesis imperfecta.

Patients' data were obtained from the medical records,
including the age of detection of the first clinical signs, clinical
examination, standardized clinical diagnostic workup findings,
associated phenotypic features, non-genetic laboratory testing,
and family history. The clinical geneticist evaluated all these data,
and the phenotypic features were described according to the
Human Phenotype Ontology nomenclature (HPO) (15).

In this study, we analyzed retrospectively the results of genetic
testing previously performed as a part of routine clinical diagnostics
at our institution. Written informed consent was obtained from all
patients. No genetic testing was performed solely for the purpose of
this study.

Next generation sequencing

ES was performed using a standardized series of procedures.
Starting with an in-solution capture of exome sequences using
various capture kits, including Nextera Coding Exome capture kit
(manufactured by Illumina, San Diego, CA, USA) or Agilent
SureSelect Human All Exon (V2, V5, V6) capture kits
(manufactured by Agilent Technologies, Santa Clara, CA), Twist
Human-Core-Exome kit (manufactured by Twist Bioscience
Corporation, South San Francisco, CA, USA). As for the WGS, a
standardized sequence of procedures was performed following
PCR-free WGS library preparation protocol Illumina TrueSeq
DNA Nano (manufactured by Illumina, San Diego, CA, USA).
This was followed by sequencing on Illumina MiSeq, Illumina
NextSeq 550, Illumina HiSeq 2500 Illumina, or Illumina NovaSeq 6000 platforms.

The data obtained from ES or WGS was subsequently analyzed using an internally developed pipeline based on the combined disease and phenotype gene target definition approach, as we previously described (16-18). Basic analysis for the detection and annotation of Single Nucleotide Variant (SNV), canonical splice site variants, and small indels (insertion/deletion) was performed according to the GATK Best Practices workflow (19-21). Moreover, we employed additional data analysis methods to expand the spectrum of detected genetic variation and improve the identification of causative variants in cases of negative results. To achieve this, we utilized the following approaches: (a) detection of noncanonical splice site variants (22, 23). (b) Annotation and analysis of mitochondrial sequence (c) CNV analyses (d) Identification of breakpoints was analyzed by detecting clusters of soft-clipped reads in aligned sequence, and (e) finally, the identification of repeat expansions throughout the genome and the long runs of homozygosity was analyzed (24, 25).

Only the variants classified as pathogenic or likely pathogenic were considered in estimating positive yield in this study. All the variants reported were inspected visually in the IGV browser. Low-quality variants and variants with GATK quality score below 500 were confirmed with Sanger sequencing (26). For CNVs, the confirmation was performed using different approaches: CGH array, multiplex ligation-dependent probe amplification (MLPA), and targeted PCR amplification.

## Gene selection

Based on human phenotype ontology (HPO) (17), OMIM (27), Pubmed databases (28), PanelApp (29) and ClinGen (30) tools, we created a targeted gene panel composed of 615 genes to identify relevant genes for our genetic investigation of SD. This method

ensured the selection of relevant genes for our investigation, thus enhancing the accuracy of our genetic analysis by focusing on a small number of variants. For cases where no pathogenic variants were detected in the initial gene panels, we expanded our analysis to include all variants identified in the ES or WGS datasets.

## Classification

The pathogenicity of the variants was evaluated according to the standards and guidelines provided by ACMG/AMP (American College of Medical Genetics and Genomics and the Association for Molecular Pathology) and the Association for Clinical Genomic Science (ACGS) criteria (31, 32).

For the positive cases, the disease categories were classified according to the recent Orphanet classification of bone diseases (33).

## Results

### Genetic distribution and classification of the variants detected in our study

We retrospectively reviewed 470 patients with clinical features of SD, including fetal, pediatric (birth to 18 years), and adult (>18 years) patients who were referred for genetic investigation. We used NGS in 264 patients: ES for 244 patients, while trio-ES (11 cases) or WGS (7 cases) were performed on selected patients with negative results. The largest patient group consisted of pediatric patients (n=131), followed by adults (n=94), and 37 prenatal cases (37). Using NGS, we established a positive genetic diagnosis for 50% (n=131) of the referred patients. Of these, 38.93% were initially referred with an unspecified SD diagnosis, and 10.69% had been misdiagnosed. We identified 132 distinct causative variants across 66 different genes.

<figure><img src="image_4.png" /><figcaption>FIGURE 1
Diagnostic yield of diagnostic approaches and variant distribution. (A) High rate of positive diagnosis using NGS compared to Sanger sequencing. (B) Distribution of different types of variants detected in our cohort.</figcaption></figure>
Conversely, we used Sanger sequencing for patients referred with a well-established diagnosis linked to a specific gene. We identified causative variants in 37 out of 208 patients (17.79%), detecting 26 pathogenic or likely pathogenic variants across 11 different genes (Supplementary Table S1).

Comparing diagnostic approaches, we found that NGS identified 74% of the positive cases. (Figure 1A).

Missense variants were the most common in our cohort (64.22%), followed by null variants (30.86%); frameshift variants were the largest subgroup (14.45%), while nonsense variants accounted for 11.8% of positive cases, and splice site variants for 4.6%. Furthermore, in-frame variants (in-frame deletions or insertions) were detected in 2.65% of patients. In addition, CNV analysis identified 6 cases, increasing the total diagnostic yield by 2.27%. CNVs ranged from single-exon deletions or duplications to deletions or duplications of entire genes (Figure 1B).

Among all these variants, we identified 29 novel pathogenic variants associated with SD (Supplementary Table S2). Moreover, we identified ten variants classified as VUS in nine genes (Supplementary Table S3).

## Molecular pathology of SD diseases identified in our study

We demonstrated considerable molecular heterogeneity in the Slovenian population, identifying 66 altered genes associated with various SD disorders (Figure 2). De novo variants in FGFR3 (OMIM: *134934) were the most commonly identified, found in 30 patients with achondroplasia (OMIM: #100800) and hypochondroplasia (OMIM:

#146000). Additionally, we detected FGFR3 variants in five patients with Thanatophoric dysplasia (types 1 (OMIM: #187600) and type 2 (OMIM: #187601) and three with Muenke syndrome (OMIM: #602849). In addition, 20 patients were diagnosed with Osteogenesis Imperfecta (OI) (OMIM: #166210, #259420, #166220), with ten patients showing variants in the COL1A2 gene (OMIM: *120160) and eight patients with variants in the COL1A1 gene (OMIM: *120150). In the FGFR2 gene (OMIM: *176943), five patients were diagnosed with Apert syndrome (OMIM: 101200), two with Crouzon syndrome (OMIM: #123500), and one with Pfeiffer syndrome (OMIM: 101600). Among patients with variants in the COL2A1 gene (OMIM: *120140) (n=13), ten were diagnosed with spondyloepiphyseal dysplasia congenita (OMIM: #183900) and spondyloepimetaphyseal dysplasia (OMIM: #184250), two with Stickler syndrome (OMIM: #108300), and one with achondrogenesis type II (OMIM: #200610). Six patients had pathogenic variants in the EXT1 gene (OMIM: *608177) causing hereditary multiple osteochondromas (OMIM: #608177).

## Novel genes and genotype/phenotype associations

We identified a null frameshift duplication c.295_296dup, p.(Pro100Leu fs*3) in the CBFB gene (Core-Binding Factor, Beta Subunit) (OMIM: *121360) associated with a new skeletal disorder resembling cleidocranial dysplasia syndrome (34).

The second variant c.2768T>G, p.(Leu923*), was detected in the MIA3 gene (Melanoma Inhibitory Activity Family, Member 3) (OMIM: *613455), also known as TANGO gene (Transport And Golgi Organization Gene 1) in a fetus who presented with short

<figure><img src="image_5.png" /><figcaption>FIGURE 2
Genetic heterogeneity and distribution of the mutated genes in SD patients in our cohort.</figcaption></figure>
bones of extremities (7 percentile), fibular aplasia, bilateral radial aplasia, tibial aplasia, hypoplastic nasal bone, delayed ossification, and congenital contractures. Variants in this gene were previously reported in a single case of Odontochondrodysplasia 2 with hearing loss and diabetes disease (OMIM # 619269).

Diagnosis and classification of skeletal dysplasias detected in our cohort

We found significant clinical and genetic heterogeneity in our cohort, identifying different SD syndromes. To provide a clear overview of the prevalence and distribution of these various diseases, we used the Orphanet classification for rare bone diseases (ORPHA:93419). Each identified disorder was assigned to a specific group and subgroup, resulting in a comprehensive categorization of our findings into four groups. Primary bone dysplasias were the most common group of disorders among our patients (71.6%), followed by dysostoses (21.89%), rare bone tumors represented by hereditary multiple osteochondromas (4.73%), and lysosomal storage diseases (1.78%) (Figure 3).

The primary bone dysplasias represented the most diverse group, with 17 different subgroups identified. Primary bone dysplasia with micromelia was the most commonly diagnosed subgroup (n=28), followed closely by spondyloepiphyseal and spondyloepimetaphyseal dysplasias, which also formed a significant subgroup (n=23). We identified 22 patients diagnosed with primary bone dysplasia with decreased bone density. In the dysostosis group, craniosynostosis was the largest subgroup with the highest number of diagnosed cases (n=22), followed by

dysostosis with limb anomalies as the second largest subgroup, with 13 diagnosed patients.

Among the diagnosed SD cases, autosomal dominant (AD) inheritance was the most common mode (82.99%), followed by autosomal recessive (AR) inheritance (14.72%) and X-linked dominant (XLD) inheritance (2.29%).

Discussion

In this study, we present a systematic strategy for SD diagnostics in the Slovenian health care system (population of 2 million), which provides insight into molecular pathology and genetic epidemiology of SD.

SDs are clinically and genetically heterogeneous. The recent update of the Nosology of genetic skeletal disorders identified 552 genes associated with SD, which covers 95.6% of these disorders (6). Using NGS, we achieved a genetic diagnostic yield of 50% (131/262). This diagnostic yield is comparable to the findings of other studies in the literature, with an average diagnostic yield of 44.76% (10, 35, 36). The genetic spectrum of SDs is continuously expanding, as the genetic etiology of many disorders remains partially understood. New candidate genes and novel phenotypes are being discovered regularly (13). According to a recent recommendation by the American College of Medical Genetics and Genomics, ES or WGS should be utilized as a first- or second-tier test for patients with congenital abnormalities (37). Moreover, with the advancement of NGS, 67% of previously undiagnosed disease cases were resolved using ES (34), making it the most efficient diagnostic tool in terms of time and cost (35).

<figure><img src="image_6.png" /><figcaption>FIGURE 3
Distribution and classification of the subgroups of SD detected for positive cases in our cohort. The number associated with each subgroup represents the number of patients identified with positive genetic results.</figcaption></figure>
While our diagnostic approach primarily utilized next-generation sequencing, other genomic alterations may contribute to genetic disorders. Integrating additional methodologies, such as optical genome mapping, RNA sequencing, and methylation analysis, could enhance the diagnostic accuracy for rare diseases (38).

In this study, we identified 149 different pathogenic or likely pathogenic variants in 66 different genes associated with SD, underscoring the genetic heterogeneity present within our patient population. Using NGS, we identified 132 different variants across 66 genes, including 29 novel variants, which contribute to expanding genotype-phenotype correlations. Given that CNVs and mitochondrial variants have been reported to contribute to SD (5, 13), we systematically analyzed our SD cohort for these genetic variations. Using NGS, we achieved a diagnostic yield of 4.55% (n=6) for CNVs but did not detect any mitochondrial variants. Similar results were reported by scocchia et al. with 5.7% of positive results presenting CNV variants (35).

Throughout the course of our research, using ES, we identified two potential genes, *CBFB* and *MIA3*, associated with novel phenotypes of SDs. We recently published the *CBFB* gene alteration associated with a new skeletal disorder resembling cleidocranial dysplasia (34). Moreover, we detected a novel biallelic loss of function variant c.2768T>G, p.(Leu923*) in the *MIA3* gene in a fetus with a severe skeletal dysplasia phenotype. Homozygous variants in the *MIA3* gene were associated with Odontochondrodysplasia 2 with hearing loss and diabetes disease. This was discovered in one family presenting a synonymous variant in the *MIA3* gene leading to exon skipping associated with a skeletal dysplasia phenotype characterized by short stature, various skeletal abnormalities, severe dentinogenesis imperfecta, and insulin-dependent diabetes mellitus in four children.

These two genes are still not listed in the Nosology of genetic skeletal disorders despite these studies. Studies on a large cohort of patients using ES and WGS have, so far, notably, contributed to the molecular pathology of SD (7, 39–41) and have always contributed to the discovery of new genes responsible for SD, which lead to an increasing number of genes. This was presented by the increased number of identified genes in the updated revisions of the nosology of genetic skeletal disorders from 226 in 2010 (1) to 552 disease-causing genes that cover 737 skeletal disorders in the last revision in 2023 (6).

In addition to the genetic heterogeneity, there's a large diversity of clinical manifestations of skeletal dysplasia. The recent update of the Nosology of genetic skeletal disorders reported 771 different disorders classified into 42 different groups (6). In our study, we categorized the detected disorders into four major groups: primary bone dysplasias, dysostosis, rare bone tumors, and lysosomal storage diseases. Primary bone dysplasias emerged as the most common and heterogeneous group in our cohort.

This clinical heterogeneity spanned diverse age categories within our cohort, with a particular emphasis on prenatal cases and carrier screening, both of which have significant implications for the healthcare system. This category showed a high diagnostic rate of 38.51% (n=181), with a positive yield of 23.2% (n=42) using both diagnostic approaches. Consistent with the findings of South America, fetal cases constitute 40% of all cases with genetic SDs (42). Using NGS, we achieved a diagnostic yield of prenatal diagnosis of 54.05% (n=20). Similar results using NGS for prenatal diagnosis

were observed in the literature, summarized in the literature review, with an average positive diagnostic yield of 69% (43).

Fetal congenital anomalies not only increase infant morbidity and mortality but also cause intangible suffering for families. Ultrasound detection rates of fetal abnormalities are observed in 2–3% of pregnancies (44). Timely, accurate diagnoses and appropriate interventions for congenital anomalies are therefore crucial. Given that SDs are among the most common fetal malformations (45), prenatal diagnosis is essential. From a public health perspective, experts and policymakers have opted that newborn and carrier screening should be conducted as a state-run program, especially for rare diseases (46, 47). This will allow early detection of genetic conditions and congenital disorders, reducing health disparities and mortality at a young age and early parental decision-making (48).

Clinical and genetic diagnosis has an important impact on the patient's clinical care and participates in initiating the appropriate patient treatment management. Recently, precise therapies such as gene therapies and authorized orphan drugs have been developed, particularly for rare diseases that are targeted to some specific altered genes. A definitive molecular diagnosis is essential for these specialized treatments, underscoring the role of genetic diagnostics in advancing precision medicine. To support potential therapeutic strategies for our patients, we examined the availability of specific medical interventions and orphan drugs associated with the identified diseases and genes. This analysis included resources such as the NHGRI Clinical Genomics Database (CGD) clinical categorization (49), Orphadata for orphan drugs (50), and RxGenes (51). We found registered drugs with orphan designation for 11 genes associated with different disorders, of which six are already in use as they have marketing authorization, with 2 of them using enzyme replacement therapy. In addition, gene therapy has been developed for five disorders. We also identified mechanism-based therapeutics, therapies targeting pathological mechanisms identified in specific diseases for five genes associated with different disorders. Finally, 31 disorders were identified with potentially beneficial interventions targeting certain clinical features, which could be advantageous for our patients.

In detail, we cite here examples of treatments for diseases that are highly presented in our cohort. The vosoritide is an approved drug for disorders related to *FGFR3* pathogenic variants (43, 52). Moreover, for this treatment, phase 2 and 3 clinical trials are currently ongoing for children between 0 to < 60 months and 5 to 18 years (2, 44). In our study, 38 patients present pathogenic variants in this gene, which could benefit from this treatment. The Asfotase alfa, which has been recently approved by the European Medicines Agency (EMA), should be started before the age of 5 years old for the most severe forms of hypophosphatasia (perinatal/infancy) (53). This could be an enzyme replacement therapeutic target for four patients in our cohort positively diagnosed. We cite also, the clinical trial started to evaluate the Denosumab treatment with moderate to severe OI use in children between the age of 5–10 years (ClinicalTrials.gov identifier: NCT01799798). It is also used as a potent and effective treatment for bone resorption pathologies such as osteoporosis and different types of bone tumors (54). OI is the most common disease detected in our patients (20 patients); in addition, infantile osteoporosis and osteochondromas are detected in 10 patients in our cohort. Hence, the discovery of the genetic etiology of our patients could help for better
clinical management and treatment guidance for the referral clinicians
and the patients.

The strength of our approach lies in the systematic
implementation of next-generation sequencing within the
national healthcare system, facilitating broad access to genetic
testing. However, the generalizability of our findings regarding
diagnostic yield and molecular pathology may be influenced by
the unique genetic characteristics of the Slovenian population.In
conclusion, this study outlines a national strategy for the genetic
diagnosis of SDs in Slovenia, presenting data on the genetic
epidemiology of these disorders within the Slovenian population
and revealing substantial clinical and genetic variability. Our
findings contribute to expanding the mutational spectrum by
identifying novel variants and genes associated with extremely
rare SD entities, supporting enhanced genotype-phenotype
correlations. Furthermore, our experience underscores the utility
of NGS for achieving molecular diagnoses and improving
diagnostic yield. By integrating both SNV and CNV analyses
from NGS data, our approach demonstrates a robust diagnostic
yield, advocating for the inclusion of CNV analysis to aid in
diagnosing unresolved cases. These genetic investigations refine
clinical diagnostics and enable the application of emerging, targeted
treatments including enzyme replacement therapy and gene therapy
that show promising results for specific conditions, based on patient
genetic profiles. Additionally, many clinical trials show promising
treatments for skeletal dysplasia; however, further research is
needed to refine and adapt these therapies, integrating genetic
findings with therapeutic developments for broader patient benefit.

Data availability statement

The datasets presented in this study can be found in online
repositories. The names of the repository/repositories and accession
number(s) can be found in the article/Supplementary Material.

Ethics statement

The studies involving humans were approved by the Slovenian
National Medical Ethics Committee, Ministry of Health of the
Republic of Slovenia. The studies were conducted in accordance
with the local legislation and institutional requirements. Written
informed consent for participation in this study was provided by the
participants' legal guardians/next of kin.

Author contributions

DN: Data curation, Formal analysis, Investigation, Methodology, Visualization, Writing – original draft, Writing – review & editing.

AM: Data curation, Investigation, Software, Supervision, Validation, Writing - review & editing. AP: Data curation, Investigation, Writing - review & editing. HJ: Data curation, Investigation, Writing - review & editing. BP: Conceptualization, Funding acquisition, Methodology, Project administration, Resources, Supervision, Validation, Visualization, Writing - review & editing.

Funding

The author(s) declare financial support was received for the
research, authorship, and/or publication of this article. Slovenian
Research Agency (ARRS), grant number P3-0326 and J3-4517.

Acknowledgments

We would like to thank the participating patients and family
members for their collaboration, as well as the clinicians who
referred patients and provided genetic counseling.

Conflict of interest

The authors declare that the research was conducted in the
absence of any commercial or financial relationships that could be
construed as a potential conflict of interest.

Generative AI statement

The author(s) declare that no Generative AI was used in the
creation of this manuscript.

Publisher’s note

All claims expressed in this article are solely those of the
authors and do not necessarily represent those of their affiliated
organizations, or those of the publisher, the editors and the
reviewers. Any product that may be evaluated in this article, or
claim that may be made by its manufacturer, is not guaranteed or
endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online
at: https://www.frontiersin.org/articles/10.3389/fendo.2025.
1523737/full#supplementary-material
## References

1. Warman ML, Cormier-Daire V, Hall C, Krakow D, Lachman R, LeMerrer M, et al. Nosology and classification of genetic skeletal disorders: 2010 revision. *Am J Med Genet A.* (2011) 155:943–68. doi: 10.1002/ajmg.a.155.5

2. Sabir AH, Cole T. The evolving therapeutic landscape of genetic skeletal disorders. *Orphanet J Rare Diseases.* (2019) 14:300. doi: 10.1186/s13023-019-1222-2

3. Kolambage YD, Walpita YN, Liyanage UA, Dayaratne BM, Dissanayake VH. The burden of hospital admissions for skeletal dysplasias in Sri Lanka: a population-based study. *Orphanet J Rare Diseases.* (2023) 18:279. doi: 10.1186/s13023-023-02884-2

4. Krakow D. Skeletal dysplasias. *Clinics perinatology.* (2015) 42:301–19. doi: 10.1016/j.clp.2015.03.003

5. Geister KA, Camper SA. Advances in skeletal dysplasia genetics. *Annu Rev Genomics Hum Genet.* (2015) 16:199–227. doi: 10.1146/annurev-genom-090314-045904

6. Unger S, Ferreira CR, Mortier GR, Ali H, Bertola DR, Calder A, et al. Nosology of genetic skeletal disorders: 2023 revision. *Am J Med Genet A.* (2023) 191:1164–209. doi: 10.1002/ajmg.a.191.5

7. Mortier GR, Cohn DH, Cormier-Daire V, Hall C, Krakow D, Mundlos S, et al. Nosology and classification of genetic skeletal disorders: 2019 revision. *Am J Med Genet Part A.* (2019) 179:2393–419. doi: 10.1002/ajmg.a.61366

8. Cui Y, Zhao H, Liu Z, Liu C, Luan J, Zhou X, et al. A systematic review of genetic skeletal disorders reported in Chinese biomedical journals between 1978 and 2012. *Orphanet J Rare Diseases.* (2012) 7:55. doi: 10.1186/1750-1172-7-55

9. Faruqi T, Dhawan N, Bahl J, Gupta V, Vohra S, Tu K, et al. Molecular, phenotypic aspects and therapeutic horizons of rare genetic bone disorders. *BioMed Res Int.* (2014) 2014:670842. doi: 10.1155/2014/670842

10. Bae J-S, Kim NKD, Lee C, Kim SC, Lee HR, Song H-R, et al. Comprehensive genetic exploration of skeletal dysplasia using targeted exome sequencing. *Genet Med.* (2016) 18:563–9. doi: 10.1038/gim.2015.129

11. Mody KS, Henstenburg J, Herman MJ. The health & Economic disparities of congenital musculoskeletal disease worldwide: an analysis of 25 years (1992-2017). *Glob Pediatr Health.* (2021) 8:2333794x2199498. doi: 10.1177/2333794X2199498

12. Seefried L, Smyth M, Keen R, Harvengt P. Burden of disease associated with X-linked hypophosphataemia in adults: a systematic literature review. *Osteoporosis Int.* (2021) 32:7–22. doi: 10.1007/s00198-020-05548-0

13. Calder AD. The changing world of skeletal dysplasia. *Lancet Child Adolesc Health.* (2020) 4:253–4. doi: 10.1016/S2352-4642(20)30056-0

14. Falardeau F, Camurri MV, Campeau PM. Genomic approaches to diagnose rare bone disorders. *Bone.* (2017) 102:5–14. doi: 10.1016/j.bone.2016.07.020

15. Briggs MD, Bell PA, Wright MJ, Pirog KA. New therapeutic targets in rare genetic skeletal diseases. *Expert Opin orphan Drugs.* (2015) 3:1137–54. doi: 10.1517/21678707.2015.1083853

16. Bergant G, Maver A, Lovrecic L, Čuturilo G, Hodzic A, Peterlin B. Comprehensive use of extended exome analysis improves diagnostic yield in rare disease: a retrospective survey in 1,059 cases. *Genet Med.* (2018) 20:303–12. doi: 10.1038/gim.2017.142

17. Ales M, Luca L, Marija V, Gorazd R, Karin W, Ana B, et al. Phenotype-driven gene target definition in clinical genome-wide sequencing data interpretation. *Genet Med.* (2016) 18:1102–10. doi: 10.1038/gim.2016.22

18. Bergant G, Maver A, Peterlin B. Whole-genome sequencing in diagnostics of selected Slovenian undiagnosed patients with rare disorders. *Life.* (2021) 11:205. doi: 10.3390/life11030205

19. Wang K, Li M, Hakonarson H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. *Nucleic Acids Res.* (2010) 38:e164–e. doi: 10.1093/nar/gkq603

20. McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The Genome Analysis Toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data. *Genome Res.* (2010) 20:1297–303. doi: 10.1101/gr.107524.110

21. Molina-Ramirez LP, Kyle C, Ellingford JM, Wright R, Taylor A, Bhaskar SS, et al. Personalised virtual gene panels reduce interpretation workload and maintain diagnostic rates of proband-only clinical exome sequencing for rare disorders. *J Med Genet.* (2022) 59:393–8. doi: 10.1136/jmedgenet-2020-107303

22. Jian X, Boerwinkle E, Liu X. In silico prediction of splice-altering single nucleotide variants in the human genome. *Nucleic Acids Res.* (2014) 42:13534–44. doi: 10.1093/nar/gku1206

23. Jaganathan K, Panagiotopoulou SK, McRae JF, Darbandi SF, Knowles D, Li YI, et al. Predicting splicing from primary sequence with deep learning. *Cell.* (2019) 176:535–48. doi: 10.1016/j.cell.2018.12.015

24. Dolzhenko E, Deshpande V, Schlesinger F, Krusche P, Petrovski R, Chen S, et al. ExpansionHunter: a sequence-graph-based tool to analyze variation in short tandem repeat regions. *Bioinformatics.* (2019) 35:4754–6. doi: 10.1093/bioinformatics/btz431

25. Van der Auwera GA, Carneiro MO, Hartl C, Poplin R, del Angel G, Levy-Moonshine A, et al. From fastQ data to high-confidence variant calls: the genome analysis toolkit best practices pipeline. *Curr Protoc Bioinf.* (2013) 43(1):11.10.11-11.10.33. doi: 10.1002/0471250953.2013.43.issue-1

26. Strom SP, Lee H, Das K, Vilain E, Nelson SF, Grody WW, et al. Assessing the necessity of confirmatory testing for exome-sequencing results in a clinical molecular diagnostic laboratory. *Genet Med.* (2014) 16:510–5. doi: 10.1038/gim.2013.183

27. Amberger JS, Bocchini CA, Schiettecatte F, Scott AF, Hamosh A. OMIM.org: Online Mendelian Inheritance in Man (OMIM®), an online catalog of human genes and genetic disorders. *Nucleic Acids Res.* (2015) 43:D789–98. doi: 10.1093/nar/gku1205

28. White J. PubMed 2.0. *Medical Reference Services Quarterly.* (2020) 39(4):382–7. doi: 10.1080/02763869.2020.1826228

29. Martin AR, Williams E, Foulger RE, Leigh S, Daugherty LC, Niblock O, et al. PanelApp crowdsources expert knowledge to establish consensus diagnostic gene panels. *Nat Genet.* (2019) 51:1560–5. doi: 10.1038/s41588-019-0528-2

30. ClinGen the clinical genome resource (2015). Available online at: https://clinicalgenome.org/ (Accessed December 2022).

31. Durkie M, Cassidy E-J, Berry I, Owens M, Turnbull C, Scott RH, et al. ACGS best practice guidelines for variant classification in rare. United Kingdom: Association for Clinical Genomic Science (ACGS). (2023).

32. Richards S, Aziz N, Bale S, Bick D, Das S, Gastier-Foster J, et al. Standards and guidelines for the interpretation of sequence variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics and the Association for Molecular Pathology. *Genet Med.* (2015) 17:405–23. doi: 10.1038/gim.2015.30

33. Pavan S, Rommel K, Mateo Marquina ME, Höhn S, Lanneau V, Rath A. Clinical practice guidelines for rare diseases: the orphanet database. *PloS One.* (2017) 12:e0170365. doi: 10.1371/journal.pone.0170365

34. Beyltjens T, Boudin E, Revencu N, Boeckx N, Bertrand M, Schütz L, et al. Heterozygous pathogenic variants involving CBFB cause a new skeletal disorder resembling cleidocranial dysplasia. *J Med Genet.* (2023) 60:498–504. doi: 10.1136/jmg-2022-108739

35. Scocchia A, Kangas-Kontio T, Irving M, Hero M, Saarinen I, Pelttari L, et al. Diagnostic utility of next-generation sequencing-based panel testing in 543 patients with suspected skeletal dysplasia. *Orphanet J Rare Diseases.* (2021) 16:412. doi: 10.1186/s13023-021-02025-7

36. Retterer K, Juusola J, Cho MT, Vitazka P, Millan F, Gibellini F, et al. Clinical application of whole-exome sequencing across clinical indications. *Genet Med.* (2016) 18:696–704. doi: 10.1038/gim.2015.148

37. Manickam K, McClain MR, Demmer LA, Biswas S, Kearney HM, Malinowski J, et al. Exome and genome sequencing for pediatric patients with congenital anomalies or intellectual disability: an evidence-based clinical guideline of the American College of Medical Genetics and Genomics (ACMG). *Genet Med.* (2021) 23:2029–37. doi: 10.1038/s41436-021-01242-6

38. Vinkšel M, Writl K, Maver A, Peterlin B. Improving diagnostics of rare genetic diseases with NGS approaches. *J Community Genet.* (2021) 12:247–56. doi: 10.1007/s12687-020-00500-5

39. Shaheen R, Patel N, Shamseldin H, Alzahrani F, Al-Yamany R, ALMoisheer A, et al. Accelerating matchmaking of novel dysmorphology syndromes through clinical and genomic characterization of a large cohort. *Genet Med.* (2016) 18:686–95. doi: 10.1038/gim.2015.147

40. Maddirevula S, Alsahli S, Alhabeeb L, Patel N, Alzahrani F, Shamseldin HE, et al. Expanding the phenome and variome of skeletal dysplasia. *Genet Med.* (2018) 20:1609–16. doi: 10.1038/gim.2018.50

41. Uttarilli A, Shah H, Bhavani GS, Upadhyai P, Shukla A, Girisha KM. Phenotyping and genotyping of skeletal dysplasias: Evolution of a center and a decade of experience in India. *Bone.* (2019) 120:204–11. doi: 10.1016/j.bone.2018.10.026

42. Barbosa-Buck CO, Orioli IM, da Graça Dutra M, Lopez-Camelo J, Castilla EE, Cavalcanti DP. Clinical epidemiology of skeletal dysplasias in South America. *Am J Med Genet Part A.* (2012) 158:1038–45. doi: 10.1002/ajmg.a.v158a.5

43. Tse KY, Surya IU, Irwinda R, Leung KY, Ting YH, Cao Y, et al. Diagnostic yield of exome sequencing in fetuses with sonographic features of skeletal dysplasias but normal karyotype or chromosomal microarray analysis: A systematic review. *Genes.* (2023) 14:1203. doi: 10.37766/inplasy2023.5.0048

44. Calzolari E, Barisic I, Loane M, Morris J, Wellesley D, Dolk H, et al. Epidemiology of multiple congenital anomalies in europe: a EUROCAT population-based registry study. *Birth Defects Res Part A: Clin Mol Teratology.* (2014) 100(4):270–6. doi: 10.1002/bdra.23240

45. Liu W, Cao J, Shi X, Li Y, Qiao F, Wu Y. Genetic testing and diagnostic strategies of fetal skeletal dysplasia: a preliminary study in wuhan, china. *Orphanet J Rare Dis* (2023) 18(1):336. doi: 10.1186/s13023-023-02955-4

46. Loeber JG, Platis D, Zetterström RH, Almashanu S, Boemer F, Bonham JR, et al. Neonatal screening in europe revisited: an ISNS perspective on the current state and developments since 2010. *Int J Neonatal Screen.* (2021) 7:15. doi: 10.3390/ijns7010015

47. Brosco JP, Grosse SD, Ross LF. Universal state newborn screening programs can reduce health disparities. *JAMA pediatrics.* (2015) 169:7–8. doi: 10.1001/jamapediatrics.2014.2465
48. Urv T, Parisi M. Newborn screening: beyond the spot. In: *Rare diseases epidemiology: update overview*. Springer, Berlin/Heidelberg, Germany (2017).

49. Solomon BD, Nguyen A-D, Bear KA, Wolfsberg TG. Clinical genomic database. *Proc Natl Acad Sci* (2013) 110(24):9851-5. doi: 10.1073/pnas.1302575110

50. Orphanet AS. Orpha data. (2011). Available at: http://wwworphanadataorg/php (Accessed October 23, 2012).

51. Bick D, Bick SL, Dimmock DP, Fowler TA, Caulfield MJ, Scott RH. An online compendium of treatable genetic disorders. *Am J Med Genet Part C* (2021) 187C:48-54. doi: 10.1002/ajmg.c.31874

52. Barbosa-Buck CO, Orioli IM, da Graça Dutra M, Lopez-Camelo J, Castilla EE, Cavalcanti DP. Clinical epidemiology of skeletal dysplasias in south america†. *Am J Med Genet.* (2012) 158A:1038-45. doi: 10.1002/ajmg.a.35246

53. Whyte MP, Rockman-Greenberg C, Ozono K, Riese R, Moseley S, Melian A, et al. Asfotase alfa treatment improves survival for perinatal and infantile hypophosphatasia. *J Clin Endocrinology*. (2016) 101:334-42. doi: 10.1210/jc.2015-3462

54. Argente J, Martos Moreno G. Displasias esqueléticas: nuevos tratamientos médicos. *An. pediatr.* (2003. *Ed. impr.*). (2016) 85:1-3. doi: 10.1016/j.anpedi.2016.03.015