Genetic evaluation using next-generation sequencing of children with short stature: a single tertiary-center experience

Article information

Ann Pediatr Endocrinol Metab. 2024;29(1):38-45
Publication date (electronic) : 2024 February 29
doi : https://doi.org/10.6065/apem.2346036.018
1Department of Pediatrics, Inha University Hospital, Inha University College of Medicine, Incheon, Korea
2Northwest Gyeonggi Regional Center for Rare Disease, Inha University Hospital, Incheon, Korea
3GC Genome, Yongin, Korea
4GC Labs, Yongin, Korea
Address for correspondence: Ji-Eun Lee Department of Pediatrics, Inha University Hospital, Inha University College of Medicine, Northwest Gyeonggi Regional Center for Rare Disease, 27 Inhang-ro, Jung-gu, Incheon 22332, Korea Email: anicca@inha.ac.kr
Received 2023 February 3; Revised 2023 February 24; Accepted 2023 February 28.

Abstract

Purpose

We used next-generation sequencing (NGS) to investigate the genetic causes of suspected genetic short stature in 37 patients, and we describe their phenotypes and various genetic spectra.

Methods

We reviewed the medical records of 50 patients who underwent genetic testing using NGS for suspected genetic short stature from June 2019 to December 2022. Patients with short stature caused by nongenetic factors or common chromosomal abnormalities were excluded. Thirty-seven patients from 35 families were enrolled in this study. We administered one of three genetic tests (2 targeted panel tests or whole exome sequencing) to patients according to their phenotypes.

Results

Clinical and molecular diagnoses were confirmed in 15 of the 37 patients, for an overall diagnostic yield of 40.5%. Fifteen pathogenic/likely pathogenic variants were identified in 13 genes (ACAN, ANKRD11, ARID1B, CEP152, COL10A1, COL1A2, EXT1, FGFR3, NIPBL, NRAS, PTPN11, SHOX, SLC16A2). The diagnostic rate was highest in patients who were small for their gestational age (7 of 11, 63.6%).

Conclusions

Genetic evaluation using NGS can be helpful in patients with suspected genetic short stature who have clinical and genetic heterogeneity. Further studies are needed to develop patient selection algorithms and panels containing growth-related genes.

Highlights

· Genetic evaluation can be helpful in patients with suspected genetic short stature. We performed next-generation sequencing (NGS) in 37 patients, and overall diagnostic yield was 40.5% (15 of 37). It is necessary to develop patient selection algorithms or guidelines for genetic testing.

Introduction

Short stature, defined as a height more than 2 standard deviations below the mean of the population adjusted for age and sex, is a common reason for referral to pediatric endocrinologists. Some children with severely short stature are vulnerable to various developmental, psychological, and social problems [1]. Growth is a complex process that is influenced by multiple factors, such as heredity, nutrition, hormones, and the environment. Nevertheless, genetic variations are important, accounting for 60%–80% of linear growth, and more than 700 genes that affect the growth process have been identified [2,3]. Genetic factors that influence height include chromosomal abnormalities, pathogenic variants of a single gene, and polygenetic predispositions [4,5]. With advances in nextgeneration sequencing (NGS) technologies, particularly targeted panel tests and whole exome sequencing (WES), genetic testing has become an important method in pediatric endocrinology clinics for identifying the causes of patients' short stature. The role of WES is expanding not only to increase the diagnostic rate but also to discover new causative genes in patients with short stature of unknown cause [6-8].

Short stature can also be caused by nongenetic causes, such as malnutrition, chronic systemic diseases, endocrine or metabolic disorders, and psychological deprivation. Some patients with short stature show other clinical abnormalities, including skeletal deformity, facial dysmorphism, developmental delay (DD), and other organ defects. In many cases, the genetic and nongenetic causes of short stature might not be clearly distinguishable. Because no clinical guidelines or diagnostic algorithms have yet been established to determine when and how to apply genetic testing to patients with short stature, the decision to perform genetic tests for such patients depends on their family history, accompanying clinical signs, and radiologic findings. The diagnostic rate of genetic testing using NGS in children with short stature has been reported to vary from 9% to 46.2%, depending on the characteristics of the patient group [6,9-13]. The choice of an appropriate genetic testing method, considering the diagnostic rate and cost, is also challenging for clinicians. To obtain a high diagnostic yield from genetic testing, it is important to select patients who will most likely benefit from it and implement a well-designed panel or WES [10,13]. In this study, we used NGS (targeted panel sequencing and WES) to investigate the genetic causes of short stature in 37 patients at a single tertiary center in Korea, and we describe their phenotypes and various genetic spectra.

Materials and methods

1. Selection and description of participants

This study was approved by the Institutional Review Board (2023-01-017-000) of Inha University Hospital. We reviewed the medical records of 50 patients who visited our pediatric endocrinology clinics from June 2019 to December 2022 and underwent genetic testing for suspected genetic short stature.

The exclusion criteria were as follows: short stature caused by (1) chromosomal aneuploidy or abnormalities (such as Down syndrome, Turner syndrome, chromosome 22q11.2 deletion); (2) multiple pituitary hormone deficiency (MPHD); (3) chronic systemic disorder such as chronic renal failure, chronic inflammatory disease, or inherited metabolic disorder; and (4) growth failure secondary to psychological factors.

We enrolled 37 patients from 35 families in this study. We obtained written informed consent from all parents or gardians before genetic testing using NGS. In all cases, conventional methods such as karyotyping and single gene sequencing failed to determine the genetic cause of short stature. All enrolled patients met one or more of the following criteria: (1) family history of short stature with similar phenotypes; (2) severe short stature, defined as a height-standard deviation score (Ht-SDS) of less than -3; (3) additional congenital anomalies or facial dysmorphisms; (4) evidence of skeletal deformities; and (5) associated intellectual disability (ID) or developmental delay (DD).

We assessed family history through three generations, including parental height for the calculation of midparental height (MPH), and evaluated the patient's detailed clinical phenotypes. If skeletal deformity was suspected, a radiological skeletal survey was performed.

Height and weight are expressed as standard deviation scores (Ht-SDS and Wt-SDS) and were calculated based on growth standards for Korean children and adolescents [14]. MPH was calculated as follows: for females, (maternal height + paternal height − 13 cm)/2, and for males, (maternal height + paternal height + 13 cm)/2. Small for gestational age (SGA) was defined as birth weight below the third percentile of a standard sexspecific birth weight curve for gestational age in South Korea [15]. Insulin-like growth factor-1 (IGF-1) SDS at the first visit/referral or before growth hormone (GH) treatment was calculated according to the Korean standard reference [16]. GH deficiency (GHD) was confirmed, via 2 different GH stimulation tests using insulin, arginine, or L-dopa, as a peak GH response of <10 ng/mL in both tests.

2. Genetic analysis

In this study, we administered 1 of 3 genetic tests (targeted panel sequencing using 1 of 2 different panels or WES) to each patient according to their phenotypes and accompanying symptoms. After obtaining written informed consent from parents or guardians, genomic DNA was extracted from the peripheral blood of the patients and, where possible, their parents and/or siblings. Among the 37 enrolled patients, 12 underwent panel sequencing using a short stature-targeted panel (129 genes, Supplementary Table 1), 7 patients with suspected skeletal deformity underwent panel sequencing using a skeletal dysplasia-targeted panel (433 genes, Supplementary Table 2), and 18 patients underwent WES.

For the targeted panel sequencing, all exons of the target genes were captured using a Celemics custom panel (Celemics, Seoul, Korea) or an Integrated DNA Technologies (IDT) custom panel (IDT, Coralville, IA, USA). Sequencing was performed on the Illumina NextSeq platform (Illumina, San Diego, CA, USA) and generated 2×150 bp paired-end reads. DNA sequence reads were aligned to the reference sequence based on the public human genome build UCSC GRCh37/hg19. Alignment was performed with BWA-MEM, and duplicate reads were marked with biobambam2. Base quality recalibration and variant calling were performed with the Genome Analysis Toolkit, and annotation was performed with Variant Effect Predictor and dbNSFP. Sequence variants were classified using the ACMG/ AMP guidelines [17].

For WES, genomic DNA was extracted from proband blood. All exon regions of all human genes (~22,000) were captured by a Twist Human Core Exome kit (Twist Bioscience, South San Francisco, CA, United States). The captured regions of the genome were sequenced using a Novaseq 6,000 sequencing machine (Illumina). The raw genome sequencing data analysis, including alignment to the GRCh37/hg19 human reference genome, variant calling, and annotation, was conducted using open-source bioinformatics tools and in house software. We extracted evidence data on the pathogenicity of variants from previous studies and disease databases, particularly ClinVar (https://www.ncbi.nlm.nih.gov/clinvar/) and UniProt (https://www.uniprot.org/). For the screening and reading of genetic mutations, EVIDENCE v.3.218)), an automatic mutation analysis system developed by 3billion, was used. Sequence variants were classified using the ACMG/AMP guidelines [17,19,20].

All suspected pathogenic/likely pathogenic (P/LP) variants were confirmed by Sanger sequencing, and copy number variants were identified by chromosomal microarray (CMA). According to the inheritance patterns of the disease, results were considered positive if one or 2 P/LP variants were identified in the disease-related genes. We report cases without a variant of uncertain significance or a P/LP variant as negative results.

3. Statistical analysis

Fisher exact test and the Wilcoxon rank-sum test were performed between the groups. Results with P<0.05 were considered statistically significant. All analyses were performed using IBM SPSS Statistics ver. 26.0 (IBM Co., Armonk, NY, USA).

4. Ethical statement

This study was performed according to the Helsinki Declaration (http://www.wma.net/en/30publications/10policies/b3/) and was approved by the Institutional Review Board of Inha University Hospital (2023-01-017-000). Written informed consent was obtained from all subjects.

Results

1. Patient characteristics

We included 37 Korean patients (17 females and 20 males) from 35 nonconsanguineous families in this study. Their median age at diagnosis was 10.5 years (range, 6 months to 17.5 years), and their mean Ht-SDS at diagnosis was -3.00±0.85. Patient demographics, auxological profiles, and clinical characteristics are listed in Table 1. Severe short stature (18 of 37, 48.6%) and facial dysmorphism (17 of 37, 45.9%) were the most common clinical manifestations. Other clinical manifestations included skeletal deformities (13 of 37, 35.1%), ID/DD (13 of 37, 35.1%), SGA (11 of 37, 29.7%), family history (10 of 37, 27%), and congenital anomalies (6 of 37, 16.2%). The 37 patients with various clinical features were divided into 2 groups, one containing those with congenital anomalies, facial dysmorphism, skeletal deformities, or ID/DD (syndromic group), and one containing those with only a family history or severe short stature (nonsyndromic group). Their diagnostic rates and clinical characteristics are compared in Table 1. The basal characteristics did not differ significantly between the groups. The diagnostic rate was higher in the syndromic group than the nonsyndromic group (46.1% vs. 27.2%, P<0.01) (Table 1).

Clinical characteristics of all patients

2. Molecular diagnosis of patients

Genetic analyses using NGS and further validation studies (Sanger sequencing or CMA) identified 13 causative genes in 15 patients (7 females, 8 males). The overall diagnostic yield was 40.5%, with the highest diagnostic yield of 57.1% for targeted skeletal dysplasia panel sequencing, 41.6% for targeted short stature panel sequencing, and 27.7% for WES. The genetic variants of the 15 patients with confirmed clinical and genetic diagnoses are summarized in Table 2. The identified genes and final diagnoses of those patients were as follows: 2 Noonan syndrome, (MIM#164790 and MIM#163950) caused by an NRAS variant and PTPN11 variant, respectively; 2 hypochondroplasia (MIM#146000) by FGFR3 variants; 2 Seckel syndrome 5 (MIM#3613823) in siblings by CEP152 variants; short stature and advanced bone age, with or without early-onset osteoarthritis and/ or osteochondritis dissecans (MIM#165800) by an ACAN variant; KBG syndrome (MIM#148050) by an ANKRD11 variant; Coffin-Siris syndrome 1 (MIM#135900) by an ARID1B variant; osteogenesis imperfecta type 3 (MIM#259420) by a COL1A2 variant; metaphyseal chondrodysplasia, Schmid type (MIM#120110) by a COL10A1 variant; exostoses, multiple, type 1 (MIM#133700) by an EXT1 variant; Cornelia De Lange syndrome 1 (MIM#122470) by an NIPBL variant; Leri-Weill dyschondrosteosis (MIM# 127300) by a SHOX variant; and X-linked Allan-Herndon-Dudley syndrome (MIM#300523) by an SLC16A2 variant. One gene (CEP152) was inherited in an autosomally recessive manner, 1 gene (SLC16A2) was inherited in an X-linked manner, 1 gene (SHOX) was inherited in a pseudo-autosomally dominant manner, and the remaining 10 genes were inherited in an autosomally dominant manner.

Genetic variants of 15 patients with confirmed diagnoses

Table 3 shows the clinical characteristics and phenotypes of the positive and negative groups. The basal characteristics and auxological profiles did not differ significantly between them. In the positive group, only 2 patients (13.3%) received GH stimulation tests, and neither of them was diagnosed with GHD. On the other hand, most of the negative group received GH stimulation tests, and 64.7% of them were diagnosed with GHD. Among the major clinical features, SGA (46.6% vs. 18.1%, P=0.031) and skeletal deformity (53.3% vs. 22.7%, P=0.027) occurred more frequently in the positive group. Other clinical characteristics did not differ significantly between the groups (Table 3).

Comparison of clinical characteristics between groups

Discussion

This study used NGS to explore diagnoses of genetic short stature in 37 patients with short stature and various clinical features at a single tertiary center in South Korea. In this study, we established a molecular diagnosis in 15 patients, who together had variations in 13 different genes. These results suggest that growth is affected by various genetic factors and that patients with genetic short stature form a highly heterogeneous group. This genetic heterogeneity makes it difficult to perform genetic testing in patients with short stature and is a cause of low diagnostic yield.

Several studies have used NGS to identify genetic causes in patients with short stature [6,8-13,21,22]. Because the number of causative genes is large and patient phenotypes are diverse, the diagnostic yield differed depending on the genes included in the panel and the patient group selected. Hauer et al. [8] reported an additional diagnostic yield of 16.5% when WES was performed on 200 patients (and their families) whose disease cause was not identified by systemic phenotyping. In a Chinese study, targeted panel sequencing of 166 growth-related genes was applied to 91 patients with idiopathic short stature (ISS), and the genetic diagnostic rate was low (9%) [11]. In a recently published domestic study of 114 Korean patients with ISS and isolated GHD, targeted panel sequencing of 96 genes was performed, and the diagnostic rate was 10%, similar to the previous study [9].

However, Kim et al. [13] conducted targeted exome sequencing on 15 Korean patients with suspected syndromic growth disorder and reported a diagnostic rate as high as 46.2%. Recently, Li et al. [6] published a study in which targeted NGS and WES were performed in 814 patients with short stature and other clinical symptoms, the largest cohort to date. Their overall diagnostic yield was 44.3%, with 111 causative genes in 361 patients. Among the various clinical features associated with short stature, ID/DD had the highest diagnostic rate at 70%, followed by skeletal deformity at 64.7%, microcephaly at 56.3%, congenital anomalies or dysmorphic features at 56.2%, MPHD at 36.4%, severe isolated GHD at 25%, and SGA without catch-up growth at 20.5%. Patients with severe short stature without additional phenotypes had the lowest diagnostic rate (11.2%). In this study, the overall diagnostic rate was 40.5%, which is similar to that reported in other studies that targeted patients with suspected syndromic genetic stature. Although we included a small number of patients, the diagnostic rate of SGA (7 of 11, 63.6%) was the highest in this study, followed by skeletal deformities (8 of 13, 61.5%), ID/DD (5 of 13, 38.4%), and congenital anomalies (2 of 6, 33.3%) or facial dysmorphism (5 of 17, 29.4%).

Genes affect short stature by multiple complex mechanisms that can generally be classified as follows [23,24]: (1) Defects in the GH–IGF-1 axis; (2) defects in paracrine signaling of the growth plate; (3) defects in cartilage extracellular matrix (ECM) maintenance; (4) defects in fundamental cellular (intracellular) pathways, including RASopathies, DNA repair mechanisms, and transcriptional factors; and (5) chromosomal abnormalities, copy number variations, and imprinting disorders. In this study, patients with isolated severe GHD, MPHD, and GH insensitivity (GHI) were excluded, and genes related to the GH– IGF-1 axis were not identified. Defects in fundamental cellular (intracellular) pathways were the most important genetic cause of short stature in this study. This classification contains a large number of genes, and the associated pathophysiologic mechanisms and clinical diseases are heterogeneous. In previous studies, this pathway accounted for the largest portion of genetic causes in patients with short stature [6,8,10,11]. The phenotypes of these patient groups were very diverse, but they were often syndromic short stature accompanied by clinical features other than those of ISS [10]. Among them, abnormalities in the Ras-MAPK pathway, including atypical Noonan syndrome, is an important determinant of genetic short stature, and 2 of the patients in this study belonged to that category.

Another important genetic cause of short stature is defects in cartilage ECM construction and maintenance. Five patients with that cause were included in this study, and most of them had skeletal deformity. The ACAN gene, which encodes the proteoglycan core protein aggrecan, is an important causative gene for short stature and has recently received attention. Aggrecan is mainly expressed in cartilage growth plates, and the pathogenic variant of ACAN can accelerate bone mutations and cause early growth cessation [25]. A heterozygous variant of the ACAN gene was found at a high frequency (2.5%) in both nonsyndromic and syndromic short stature [8,9,11] and was similar to the SHOX gene (2.4%), which is known to be the most common single gene defect in short stature [26]. These patients might have severe short stature with accelerated bone age, and other symptoms, such as early-onset osteoarthritis, should be noted.

Short children with SGA also form a clinically heterogeneous group, and various factors are known to be involved in the failure of catch-up growth [27]. As mentioned above, disorders of the GH–IGF axis, paracrine signaling, cartilage ECM, and fundamental cellular (intracellular) pathways can also be affected in children with SGA without catch-up growth. The highest diagnostic rate associated with SGA in this study is thought to reflect our choice to perform genetic testing on SGA patients only when they also had other clinical signs, such as skeletal deformity, facial dysplasia, and ID/DD. Therefore, genetic testing can be considered for short children with SGA accompanied by other symptoms after ruling out maternal/placental factors and common disorders related to SGA through careful examination.27,28

Genetic diagnosis is important in short stature, especially in children, because it can influence treatment decisions and improve prognosis by preventing possible accompanying symptoms. GH treatment (GHT) is the most important treatment for children with short stature. However, many cases of genetic short stature do not have GHD; in such cases, GHT is often not helpful in improving the patient's adult height. GHT can be burdensome to families because of its costs and need for daily injections for children. GHT might not be an appropriate treatment option for patients with GHI or skeletal deformities such as severe scoliosis. In our cohort, 51.4% of the patients were tested with GH stimulation, and the prevalence of GHD was as low as 29.7%. Only 13.3% of the positive group was tested with GH stimulation, and no patient was diagnosed with GHD. On the other hand, in the negative group, most of the patients received the GH stimulation test, and 64.7% of them were diagnosed with GHD and treated with GH. These findings suggest that appropriate genetic evaluations can help to determine whether GH stimulation tests or GHT should be performed. In addition, some of the patients with definite molecular diagnoses received multidisciplinary care and additional tests for possible accompanying symptoms, including developmental disorders, other endocrine disorders, and skeletal abnormalities. Accordingly, prognosis and followup management were planned, and genetic counseling was performed.

When applying a genetic test using NGS, the decision to perform a targeted panel test or WES is important for maximizing the diagnostic rate. In this study, 2 panels were used based on the accompanying clinical features of the patient, and WES was performed when a targeted panel could not be specified due to the variety of clinical features. In previous studies, the diagnostic rates of panel tests were similar to or slightly higher than those of WES [6,8,12]; however, in our study, the diagnostic rate of the panel tests was considerably higher. This finding indicates that performing genetic evaluations according to patient phenotype is an important way to optimize the diagnostic rate.

For patients with multiple clinical features, WES can be a suitable alternative. However, WES is not available except for research purposes in Korea, so it is important to develop panels that include a wide range of genes that can accompany short stature. Furthermore, genes related to growth are still being discovered; therefore, panel designs will need to be updated regularly to increase the diagnostic rate.

Our study has some limitations. First, this study was conducted in a single center for a relatively short period of time, and the number of enrolled patients was not sufficient to demonstrate the various genetic spectra of short stature. In addition, because this study was conducted retrospectively, the clinical characteristics of the patients were heterogeneous, and selection bias was a possibility. Finally, the non-uniformity of NGS testing methods used might have affected the diagnostic rate.

In conclusion, we performed genetic evaluations using NGS in 37 patients with suspected genetic short stature. We established a molecular diagnosis in 15 (P/LP variants in 13 causative genes) of those 37 patients, for an overall diagnostic yield of 40.5%. It is necessary to develop patient selection algorithms or guidelines for genetic testing, develop a panel that includes growth-related genes with high diagnostic rates, and conduct large-scale studies to understand the characteristics of patients with genetic short stature.

Supplementary Material

Supplementary Tables 1 and 2 can be found via https://doi.org/10.6065/apem.2346036.018.

Supplementary Table 1.

Lists of including genes in a short stature targeted panel

apem-2346036-018-Supplementary-Table-1.pdf
Supplementary Table 2.

Lists of including genes in a skeletal dysplasia targeted panel

apem-2346036-018-Supplementary-Table-2.pdf

Notes

Conflicts of interest

No potential conflict of interest relevant to this article was reported.

Funding

This work was supported by Inha University/Inha University Hospital Research Grant.

Data availability

The data that support the findings of this study can be provided by the corresponding author upon reasonable request.

Author contribution

Conceptualization: SJK; Data curation: SJK, EJ; Formal analysis: SJK; Funding acquisition: SJK; Methodology: SJK, CAS; Project administration: SJK, JEL; Visualization: SJK; Writing - original draft: SJK; Writing - review & editing: SJK, JEL

Acknowledgements

The authors thank all the patients and their caregivers for their involvement in this study.

References

1. Quitmann JH, Bullinger M, Sommer R, Rohenkohl AC, Bernardino Da Silva NM. Associations between psychological problems and quality of life in pediatric short stature from patients' and parents' perspectives. PLoS One 2016;11e0153953.
2. McEvoy BP, Visscher PM. Genetics of human height. Econ Hum Biol 2009;7:294–306.
3. Dauber A, Rosenfeld RG, Hirschhorn JN. Genetic evaluation of short stature. J Clin Endocrinol Metab 2014;99:3080–92.
4. Wood AR, Esko T, Yang J, Vedantam S, Pers TH, Gustafsson S, et al. Defining the role of common variation in the genomic and biological architecture of adult human height. Nat Genet 2014;46:1173–86.
5. Hirschhorn JN, Lettre G. Progress in genome-wide association studies of human height. Horm Res 2009;71 Suppl 2:5–13.
6. Li X, Yao R, Chang G, Li Q, Song C, Li N, et al. Clinical profiles and genetic spectra of 814 Chinese children with short stature. J Clin Endocrinol Metab 2022;107:972–85.
7. Costantini A, Muurinen MH, Makitie O. New gene discoveries in skeletal diseases with short stature. Endocr Connect 2021;10:R160–74.
8. Hauer NN, Popp B, Schoeller E, Schuhmann S, Heath KE, Hisado-Oliva A, et al. Clinical relevance of systematic phenotyping and exome sequencing in patients with short stature. Genet Med 2018;20:630–8.
9. Ahn J, Oh J, Suh J, Song K, Kwon A, Chae HW, et al. Nextgeneration sequencing-based mutational analysis of idiopathic short stature and isolated growth hormone deficiency in Korean pediatric patients. Mol Cell Endocrinol 2022;544:111489.
10. Fan X, Zhao S, Yu C, Wu D, Yan Z, Fan L, et al. Exome sequencing reveals genetic architecture in patients with isolated or syndromic short stature. J Genet Genomics 2021;48:396–402.
11. Yang L, Zhang C, Wang W, Wang J, Xiao Y, Lu W, et al. Pathogenic gene screening in 91 Chinese patients with short stature of unknown etiology with a targeted nextgeneration sequencing panel. BMC Med Genet 2018;19:212.
12. Huang Z, Sun Y, Fan Y, Wang L, Liu H, Gong Z, et al. Genetic evaluation of 114 Chinese short stature children in the next generation era: a single center study. Cell Physiol Biochem 2018;49:295–305.
13. Kim YM, Lee YJ, Park JH, Lee HD, Cheon CK, Kim SY, et al. High diagnostic yield of clinically unidentifiable syndromic growth disorders by targeted exome sequencing. Clin Genet 2017;92:594–605.
14. Kim JH, Yun S, Hwang SS, Shim JO, Chae HW, Lee YJ, et al. The 2017 Korean National Growth Charts for children and adolescents: development, improvement, and prospects. Korean J Pediatr 2018;61:135–49.
15. Lee JK, Jang HL, Kang BH, Lee KS, Choi YS, Shim KS, et al. Percentile Distributions of Birth Weight according to Gestational Ages in Korea (2010-2012). J Korean Med Sci 2016;31:939–49.
16. Hyun SE, Lee BC, Suh BK, Chung SC, Ko CW, Kim HS, et al. Reference values for serum levels of insulin-like growth factor-I and insulin-like growth factor binding protein-3 in Korean children and adolescents. Clin Biochem 2012;45:16–21.
17. 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–24.
18. Seo GH, Kim T, Choi IH, Park JY, Lee J, Kim S, et al. Diagnostic yield and clinical utility of whole exome sequencing using an automated variant prioritization system, EVIDENCE. Clin Genet 2020;98:562–70.
19. Riggs ER, Andersen EF, Cherry AM, Kantarci S, Kearney H, Patel A, et al. Technical standards for the interpretation and reporting of constitutional copy-number variants: a joint consensus recommendation of the American College of Medical Genetics and Genomics (ACMG) and the Clinical Genome Resource (ClinGen). Genet Med 2020;22:245–57.
20. McCormick EM, Lott MT, Dulik MC, Shen L, Attimonelli M, Vitale O, et al. Specifications of the ACMG/AMP standards and guidelines for mitochondrial DNA variant interpretation. Hum Mutat 2020;41:2028–57.
21. Sentchordi-Montane L, Benito-Sanz S, Aza-Carmona M, Diaz-Gonzalez F, Modamio-Hoybjor S, de la Torre C, et al. High prevalence of variants in skeletal dysplasia associated genes in individuals with short stature and minor skeletal anomalies. Eur J Endocrinol 2021;185:691–705.
22. Guo MH, Shen Y, Walvoord EC, Miller TC, Moon JE, Hirschhorn JN, et al. Whole exome sequencing to identify genetic causes of short stature. Horm Res Paediatr 2014;82:44–52.
23. Wit JM, Oostdijk W, Losekoot M, van Duyvenvoorde HA, Ruivenkamp CA, Kant SG. Mechanisms in endocrinology: novel genetic causes of short stature. Eur J Endocrinol 2016;174:R145–73.
24. Andrade AC, Jee YH, Nilsson O. New genetic diagnoses of short stature provide insights into local regulation of childhood growth. Horm Res Paediatr 2017;88:22–37.
25. Nilsson O, Guo MH, Dunbar N, Popovic J, Flynn D, Jacobsen C, et al. Short stature, accelerated bone maturation, and early growth cessation due to heterozygous aggrecan mutations. J Clin Endocrinol Metab 2014;99:E1510–8.
26. Rappold GA, Fukami M, Niesler B, Schiller S, Zumkeller W, Bettendorf M, et al. Deletions of the homeobox gene SHOX (short stature homeobox) are an important cause of growth failure in children with short stature. J Clin Endocrinol Metab 2002;87:1402–6.
27. Finken MJJ, van der Steen M, Smeets CCJ, Walenkamp MJE, de Bruin C, Hokken-Koelega ACS, et al. Children born small for gestational age: differential diagnosis, molecular genetic evaluation, and implications. Endocr Rev 2018;39:851–94.
28. Netchine I, van der Steen M, Lopez-Bermejo A, Koledova E, Maghnie M. New horizons in short children born small for gestational age. Front Pediatr 2021;9:655931.

Article information Continued

Table 1.

Clinical characteristics of all patients

Characteristic All (n=37) Nonsyndromic (n=11) Syndromic (n=26) P-value
Sex, male:female 20:17 5:6 15:11 0.248
Age at diagnosis (yr) 9.39±4.41 7.90±4.28 10.03±4.31 0.201
Ht-SDS -3.00±0.85 -2.81±0.57 -3.07±0.93 0.325
Wt-SDS -2.18±1.34 -2.33±1.02 -2.11±1.44 0.604
MPH-SDS -0.82±1.00 -1.38±1.07 -0.56±0.84 0.048
IGF-1 SDS -0.11±1.40 -0.17±0.79 -0.08±1.55 0.838
GH stimulation test 19 (51.4) 7 (63.6) 12 (46.1) 0.165
GH peak (ng/mL) 11.00±8.72 11.25±6.80 10.84±9.66 0.921
GHD 11 (29.7) 3 (27.2) 8 (30.7) 0.584
FHx of short stature 10 (27.0) 4 (36.3) 6 (23.0) 0.202
SGA 11 (29.7) 2 (18.1) 9 (34.6) 0.157
Ht-SDS<-3 18 (48.6) 5 (45.4) 13 (50.0) 0.399
Congenital anomalies 6 (16.2) 0 (0) 6 (23.0) <0.01
Facial dysmorphism 17 (45.9) 0 (0) 17 (65.3) <0.01
Skeletal deformities 13 (35.1) 0 (0) 13 (50.0) <0.01
ID/DD 13 (35.1) 0 (0) 13 (50.0) <0.01
Diagnostic rate 15 (40.5) 3 (27.2) 12 (46.1) <0.01

Values are presented as mean±standard deviation or number (%).

Ht-SDS, height-standard deviation score; Wt-SDS, weight-standard deviation score; MPH, midparental height; IGF-1, insulin-like growth factor-1; GH, growth hormone; GHD, growth hormone deficiency; FHx, family history; SGA, small for gestational age; ID/DD, intellectual disability/developmental delay.

Table 2.

Genetic variants of 15 patients with confirmed diagnoses

Case No. Gene Nucleotide change Amino acid change Zygosity/inheritance ACMG Genetic test method Clinical characteristics Final diagnosis
1 ACAN NM_001374828.1:c.1429+1G>T Aberrant splicing (SpliceAI score 1, ADA score 0.99, RF score 0.94) Het/AD LP NGS_SS Severe SS Short stature and advanced bone age, with/without early-onset osteoarthritis and/or osteochondritis dissecans
2 ANKRD11 NM_013275.6:c.2197C>T p.Arg733Ter Het/AD P WES Facial dysmorphism LD/DD KBG syndrome
3 ARID1B NM_001374828.1:c.5161del p.Gln1721AsnfsTer16 Het/AD P WES Severe SS, facial dysmorphism, SGA, congenital anomaly (congenital heart disease), ID/DD Coffin-Siris syndrome 1
4&5* CEP152 NM_001194998.2:c.314G>A p.Tyr678Ter Het/AR P WES FHx, facial dysmorphism, SGA, skeletal deformities Seckel syndrome 5
NM_001194998.2:c.2034T>G p.Trp105Ter Het/AR P
6 COL10A1 NM_000493.3:c.1765T>G p.Phe589Val Het/AD LP NGS_SD SGA, skeletal deformities Metaphyseal chondrodysplasia, Schmid type
7 COL1A2 NM_000089.4:c.3250G>T p.Gly1084Cys Het/AD P NGS_SD FHx, facial dysmorphism, skeletal deformities Osteogenesis imperfecta
8 EXT1 NM_000127.2:c.637C>T p.Gln213Ter Het/AD P NGS_SD Severe SS, skeletal deformities Exostoses, multiple, type 1
9 FGFR3 NM_000142.5c.1620C>G p.Asn540Lys Het/AD P NGS_SS Severe SS Hypochondroplasia
10 FGFR3 NM_000142.5:c.1619A>C p.Asn540Thr Het/AD P WES Severe SS Hypochondroplasia
11 NIPBL NM_133433.4:c.8377C>T p.Arg2793Ter Het/AD P NGS_SS Severe SS, facial dysmorphism, skeletal deformities Cornellia De Lange syndrome
12 NRAS NM_002524.5:c.68G>C p.Gly13Ala Het/AD LP NGS_SS Facial dysmorphism, congenital anomaly (congenital heart disease), ID/DD Noonan syndrome
13 PTPN11 NM_002834.4:c.794G>A p.Arg265Gln Het/AD P NGS_SS Severe SS, FHx, SGA Noonan syndrome
14 SHOX NM_000451.4:c.785_796dup p.Val262_Ala265dup Het/PD LP NGS_SD FHx, SGA, skeletal deformities Leri-Weill dyschondrosteosis
15 SLC16A2 del(X)(q13.2) NC_000023.11:g.74420675-74423535 delins AG Hemi/XL P WES, CMA Severe SS, ID/DD X-linked Allan-Herndon-Dudley syndrome

ACMG, American College of Medical Genetics and Genomics; Het, heterozygous; AD, autosomally dominant inherited; LP, likely pathogenic; NGS_SS, next-generation sequencing panel_short stature; SS, short stature; P, pathogenic; WES, whole exome sequencing; LD, learning difficulty; DD, developmental delay; Hemi, hemizygous; AR, autosomally recessive inherited; ID/DD, intellectual disability/developmental delay; FHx, family history; NGS_SD, next-generation sequencing panel_skeletal dysplasia; SGA, small for gestational age; PD, pseudo-autosomally dominant inherited; CMA, chromosomal microarray.

*

Case numbers 4 and 5 are siblings.

Table 3.

Comparison of clinical characteristics between groups

Clinical characteristics All (n=37) Positive (n=15) Negative (n=22) P-value
Sex, male:female 20:17 8:7 12:10 0.528
Age at diagnosis (yr) 9.39±4.41 8.09±4.82 10.28±3.86 0.166
Ht-SDS -3.00±0.85 -3.00±0.96 -3.00±0.77 0.991
Wt-SDS -2.18±1.34 -1.87±1.23 -2.39±1.37 0.245
MPH-SDS -0.82±1.00 -0.71±1.12 -0.89±0.89 0.615
IGF-1 SDS -0.11±1.40 -0.04±0.87 -0.16±1.65 0.755
GH stimulation test 19 (51.4) 2 (13.3) 17 (77.2) <0.01
GH peak 11.00±8.72 19.40±1.80 10.01±8.68 0.002
GHD 11 (29.7) 0/2 (0) 11/17 (64.7) <0.01
Age at short stature onset
 Neonate 11 (29.7) 5 (33.3) 6 (27.2) 0.345
 Infancy 19 (51.4) 7 (46.6) 12 (54.5) 0.681
 Childhood 7 (18.9) 3 (20.0) 4 (18.1) 0.442
FHx of short stature 10 (27.0) 5 (33.3) 5 (22.7) 0.237
SGA 11 (29.7) 7 (46.6) 4 (18.1) 0.031
Ht-SDS<-3 18 (48.6) 7 (46.6) 11 (50.0) 0.580
Congenital anomalies 6 (16.2) 2 (13.3) 4 (18.1) 0.301
Facial dysmorphism 17 (45.9) 5 (33.3) 12 (54.5) 0.101
Skeletal deformity 13 (35.1) 8 (53.3) 5 (22.7) 0.027
ID/DD 13 (35.1) 5 (33.3) 8 (36.3) 0.574

Values are presented as mean±standard deviation or number (%).

Ht-SDS, height-standard deviation score; Wt-SDS, weight-standard deviation score; MPH, midparental height; IGF-1, insulin-like growth factor-1; GH, growth hormone; GHD, growth hormone deficiency; FHx, family history; SGA, small for gestational age; ID/DD, intellectual disability/developmental delay.