A deep intronic variant in MME causes autosomal recessive Charcot-Marie-Tooth neuropathy through aberrant splicing.

BACKGROUND: Loss-of-function variants in MME (membrane metalloendopeptidase) are a known cause of recessive Charcot-Marie-Tooth Neuropathy (CMT). A deep intronic variant, MME c.1188+428A>G (NM_000902.5), was identified through whole genome sequencing (WGS) of two Australian families with recessive inheritance of axonal CMT using the seqr platform. MME c.1188+428A>G was detected in a homozygous state in Family 1, and in a compound heterozygous state with a known pathogenic MME variant (c.467del; p.Pro156Leufs*14) in Family 2. AIMS: We aimed to determine the pathogenicity of the MME c.1188+428A>G variant through segregation and splicing analysis. METHODS: The splicing impact of the deep intronic MME variant c.1188+428A>G was assessed using an in vitro exon-trapping assay. RESULTS: The exon-trapping assay demonstrated that the MME c.1188+428A>G variant created a novel splice donor site resulting in the inclusion of an 83 bp pseudoexon between MME exons 12 and 13. The incorporation of the pseudoexon into MME transcript is predicted to lead to a coding frameshift and premature termination codon (PTC) in MME exon 14 (p.Ala397ProfsTer47). This PTC is likely to result in nonsense mediated decay (NMD) of MME transcript leading to a pathogenic loss-of-function. INTERPRETATION: To our knowledge, this is the first report of a pathogenic deep intronic MME variant causing CMT. This is of significance as deep intronic variants are missed using whole exome sequencing screening methods. Individuals with CMT should be reassessed for deep intronic variants, with splicing impacts being considered in relation to the potential pathogenicity of variants.

Advances and challenges in modeling inherited peripheral neuropathies using iPSCs.

Inherited peripheral neuropathies (IPNs) are a group of diseases associated with mutations in various genes with fundamental roles in the development and function of peripheral nerves. Over the past 10 years, significant advances in identifying molecular disease mechanisms underlying axonal and myelin degeneration, acquired from cellular biology studies and transgenic fly and rodent models, have facilitated the development of promising treatment strategies. However, no clinical treatment has emerged to date. This lack of treatment highlights the urgent need for more biologically and clinically relevant models recapitulating IPNs. For both neurodevelopmental and neurodegenerative diseases, patient-specific induced pluripotent stem cells (iPSCs) are a particularly powerful platform for disease modeling and preclinical studies. In this review, we provide an update on different in vitro human cellular IPN models, including traditional two-dimensional monoculture iPSC derivatives, and recent advances in more complex human iPSC-based systems using microfluidic chips, organoids, and assembloids.

Genome sequencing reanalysis increases the diagnostic yield in dystonia.

PURPOSE: We investigated the contribution of genomic data reanalysis to the diagnostic yield of dystonia patients who remained undiagnosed after prior genome sequencing. METHODS: Probands with heterogeneous dystonia phenotypes who underwent initial genome sequencing (GS) analysis in 2019 were included in the reanalysis, which was performed through gene-specific discovery collaborations and systematic genomic data reanalysis. RESULTS: Initial GS analysis in 2019 (n = 111) identified a molecular diagnosis in 11.7 % (13/111) of cases. Reanalysis between 2020 and 2023 increased the diagnostic yield by 7.2 % (8/111); 3.6 % (4/111) through focused gene-specific clinical correlation collaborative efforts [VPS16 (two probands), AOPEP and POLG], and 3.6 % (4/111) by systematic reanalysis completed in 2023 [NUS1 (two probands) and DDX3X variants, and a microdeletion encompassing VPS16]. Seven of these patients had a high phenotype-based dystonia score ≥3. Notable unverified findings in four additional cases included suspicious variants of uncertain significance in FBXL4 and EIF2AK2, and potential phenotypic expansion associated with SLC2A1 and TREX1 variants. CONCLUSION: GS data reanalysis increased the diagnostic yield from 11.7 % to 18.9 %, with potential extension up to 22.5 %. While optimal timing for diagnostic reanalysis remains to be determined, this study demonstrates that periodic re-interrogation of dystonia GS datasets can provide additional genetic diagnoses, which may have significant implications for patients and their families.

An Update on the Adult-Onset Hereditary Cerebellar Ataxias: Novel Genetic Causes and New Diagnostic Approaches.

The hereditary cerebellar ataxias (HCAs) are rare, progressive neurologic disorders caused by variants in many different genes. Inheritance may follow autosomal dominant, autosomal recessive, X-linked or mitochondrial patterns. The list of genes associated with adult-onset cerebellar ataxia is continuously growing, with several new genes discovered in the last few years. This includes short-tandem repeat (STR) expansions in RFC1, causing cerebellar ataxia, neuropathy, vestibular areflexia syndrome (CANVAS), FGF14-GAA causing spinocerebellar ataxia type 27B (SCA27B), and THAP11. In addition, the genetic basis for SCA4, has recently been identified as a STR expansion in ZFHX3. Given the large and growing number of genes, and different gene variant types, the approach to diagnostic testing for adult-onset HCA can be complex. Testing methods include targeted evaluation of STR expansions (e.g. SCAs, Friedreich ataxia, fragile X-associated tremor/ataxia syndrome, dentatorubral-pallidoluysian atrophy), next generation sequencing for conventional variants, which may include targeted gene panels, whole exome, or whole genome sequencing, followed by various potential additional tests. This review proposes a diagnostic approach for clinical testing, highlights the challenges with current testing technologies, and discusses future advances which may overcome these limitations. Implementing long-read sequencing has the potential to transform the diagnostic approach in HCA, with the overall aim to improve the diagnostic yield.

Genetics of inherited peripheral neuropathies and the next frontier: looking backwards to progress forwards.

Inherited peripheral neuropathies (IPNs) encompass a clinically and genetically heterogeneous group of disorders causing length-dependent degeneration of peripheral autonomic, motor and/or sensory nerves. Despite gold-standard diagnostic testing for pathogenic variants in over 100 known associated genes, many patients with IPN remain genetically unsolved. Providing patients with a diagnosis is critical for reducing their 'diagnostic odyssey', improving clinical care, and for informed genetic counselling. The last decade of massively parallel sequencing technologies has seen a rapid increase in the number of newly described IPN-associated gene variants contributing to IPN pathogenesis. However, the scarcity of additional families and functional data supporting variants in potential novel genes is prolonging patient diagnostic uncertainty and contributing to the missing heritability of IPNs. We review the last decade of IPN disease gene discovery to highlight novel genes, structural variation and short tandem repeat expansions contributing to IPN pathogenesis. From the lessons learnt, we provide our vision for IPN research as we anticipate the future, providing examples of emerging technologies, resources and tools that we propose that will expedite the genetic diagnosis of unsolved IPN families.

A missense mutation in human INSC causes peripheral neuropathy.

PAR3/INSC/LGN form an evolutionarily conserved complex required for asymmetric cell division in the developing brain, but its post-developmental function and disease relevance in the peripheral nervous system (PNS) remains unknown. We mapped a new locus for axonal Charcot-Marie-Tooth disease (CMT2) and identified a missense mutation c.209 T > G (p.Met70Arg) in the INSC gene. Modeling the INSCM70R variant in Drosophila, we showed that it caused proprioceptive defects in adult flies, leading to gait defects resembling those in CMT2 patients. Cellularly, PAR3/INSC/LGN dysfunction caused tubulin aggregation and necrotic neurodegeneration, with microtubule-stabilizing agents rescuing both morphological and functional defects of the INSCM70R mutation in the PNS. Our findings underscore the critical role of the PAR3/INSC/LGN machinery in the adult PNS and highlight a potential therapeutic target for INSC-associated CMT2.

Role of the repeat expansion size in predicting age of onset and severity in RFC1 disease.

RFC1 disease, caused by biallelic repeat expansion in RFC1, is clinically heterogeneous in terms of age of onset, disease progression and phenotype. We investigated the role of the repeat size in influencing clinical variables in RFC1 disease. We also assessed the presence and role of meiotic and somatic instability of the repeat. In this study, we identified 553 patients carrying biallelic RFC1 expansions and measured the repeat expansion size in 392 cases. Pearson's coefficient was calculated to assess the correlation between the repeat size and age at disease onset. A Cox model with robust cluster standard errors was adopted to describe the effect of repeat size on age at disease onset, on age at onset of each individual symptoms, and on disease progression. A quasi-Poisson regression model was used to analyse the relationship between phenotype and repeat size. We performed multivariate linear regression to assess the association of the repeat size with the degree of cerebellar atrophy. Meiotic stability was assessed by Southern blotting on first-degree relatives of 27 probands. Finally, somatic instability was investigated by optical genome mapping on cerebellar and frontal cortex and unaffected peripheral tissue from four post-mortem cases. A larger repeat size of both smaller and larger allele was associated with an earlier age at neurological onset [smaller allele hazard ratio (HR) = 2.06, P < 0.001; larger allele HR = 1.53, P < 0.001] and with a higher hazard of developing disabling symptoms, such as dysarthria or dysphagia (smaller allele HR = 3.40, P < 0.001; larger allele HR = 1.71, P = 0.002) or loss of independent walking (smaller allele HR = 2.78, P < 0.001; larger allele HR = 1.60; P < 0.001) earlier in disease course. Patients with more complex phenotypes carried larger expansions [smaller allele: complex neuropathy rate ratio (RR) = 1.30, P = 0.003; cerebellar ataxia, neuropathy and vestibular areflexia syndrome (CANVAS) RR = 1.34, P < 0.001; larger allele: complex neuropathy RR = 1.33, P = 0.008; CANVAS RR = 1.31, P = 0.009]. Furthermore, larger repeat expansions in the smaller allele were associated with more pronounced cerebellar vermis atrophy (lobules I-V ß = -1.06, P < 0.001; lobules VI-VII ß = -0.34, P = 0.005). The repeat did not show significant instability during vertical transmission and across different tissues and brain regions. RFC1 repeat size, particularly of the smaller allele, is one of the determinants of variability in RFC1 disease and represents a key prognostic factor to predict disease onset, phenotype and severity. Assessing the repeat size is warranted as part of the diagnostic test for RFC1 expansion.

Two New Families and a Literature Review of ELOVL4-Associated Spinocerebellar Ataxia Type 34.

Autosomal dominant variants in ELOVL4 cause spinocerebellar ataxia type 34 (SCA34; ATX-ELOVL4), classically associated with a skin condition known as erythrokeratoderma. Here, we report a large Italian-Maltese-Australian family with spinocerebellar ataxia. Notably, while there were dermatological manifestations (eczema), erythrokeratoderma was not present. Using a next-generation sequencing panel, we identified a previously reported ELOVL4 variant, NM_022726.4: c.698C > T p.(Thr233Met). The variant was initially classified as a variant of uncertain significance; however, through segregation studies, we reclassified the variant as likely pathogenic. We next identified an individual from another family (Algerian-Maltese-Australian) with the same ELOVL4 variant with spinocerebellar ataxia but without dermatological manifestations. We subsequently performed the first dedicated literature review of ELOVL4-associated ataxia to gain further insights into genotype-phenotype relationships. We identified a total of 60 reported cases of SCA34 to date. The majority had gait ataxia (88.3%), limb ataxia (76.7%), dysarthria (63.3%), and nystagmus (58.3%). Of note, skin lesions related to erythrokeratoderma were seen in a minority of cases (33.3%). Other extracerebellar manifestations included pyramidal tract signs, autonomic disturbances, retinitis pigmentosa, and cognitive impairment. For brain MRI data, cerebellar atrophy was seen in all cases (100%), whereas the hot cross bun sign (typically associated with multiple system atrophy type C) was seen in 32.4% of cases. Our family study and literature review highlight the variable phenotypic spectrum of SCA34. Importantly, it shows that erythrokeratoderma is not found in most cases and that, while a dermatological assessment may be helpful in these patients, SCA34 diagnosis should be considered irrespective of dermatological manifestations.

Upper and lower limb tremor in Charcot-Marie-Tooth neuropathy type 1A and the implications for standing balance.

BACKGROUND: Neuropathic tremor occurs in Charcot-Marie-Tooth neuropathy type 1A (CMT1A; hereditary motor and sensory neuropathy, HMSN), although the pathophysiological mechanisms remain to be elucidated. Separately, lower limb tremor has not been explored in CMT1A and could be associated with imbalance as in other neuropathies. The present study aimed to determine tremor characteristics in the upper and lower limbs in CMT1A and relate these findings to clinical disability, particularly imbalance. METHODS: Tremor and posturography studies were undertaken in phenotyped and genotyped CMT1A patients. Participants underwent detailed clinical assessment, tremor study recordings, and nerve conduction studies. Tremor stability index was calculated for upper limb tremor and compared to essential tremor. RESULTS: Seventeen patients were enrolled. Postural and kinetic upper limb tremors were evident in 65%, while postural and orthostatic lower limb tremors were seen in 35% of CMT1A patients. Peak upper limb frequencies were lower distally (~ 6 Hz) and higher proximally (~ 9 Hz), were unchanged by weight-loading, and not impacted by fatigue. The tremor stability index was significantly higher in CMT1A than in essential tremor. A 5-6 Hz lower limb tremor was recorded which did not vary along the limb and was unaffected by fatigue. Balance was impaired in patients with postural lower limb tremor. A high frequency peak on posturography was associated with 'good' balance. CONCLUSIONS: Tremor is a common clinical feature in CMT1A, distinct from essential tremor, mediated by a complex interaction between peripheral and central mechanisms. Postural lower limb tremor is associated with imbalance; strategies aimed at tremor modulation could be of therapeutic utility.

Normal and pathogenic variation of RFC1 repeat expansions: implications for clinical diagnosis.

Cerebellar ataxia, neuropathy and vestibular areflexia syndrome (CANVAS) is an autosomal recessive neurodegenerative disease, usually caused by biallelic AAGGG repeat expansions in RFC1. In this study, we leveraged whole genome sequencing data from nearly 10 000 individuals recruited within the Genomics England sequencing project to investigate the normal and pathogenic variation of the RFC1 repeat. We identified three novel repeat motifs, AGGGC (n = 6 from five families), AAGGC (n = 2 from one family) and AGAGG (n = 1), associated with CANVAS in the homozygous or compound heterozygous state with the common pathogenic AAGGG expansion. While AAAAG, AAAGGG and AAGAG expansions appear to be benign, we revealed a pathogenic role for large AAAGG repeat configuration expansions (n = 5). Long-read sequencing was used to characterize the entire repeat sequence, and six patients exhibited a pure AGGGC expansion, while the other patients presented complex motifs with AAGGG or AAAGG interruptions. All pathogenic motifs appeared to have arisen from a common haplotype and were predicted to form highly stable G quadruplexes, which have previously been demonstrated to affect gene transcription in other conditions. The assessment of these novel configurations is warranted in CANVAS patients with negative or inconclusive genetic testing. Particular attention should be paid to carriers of compound AAGGG/AAAGG expansions when the AAAGG motif is very large (>500 repeats) or the AAGGG motif is interrupted. Accurate sizing and full sequencing of the satellite repeat with long-read sequencing is recommended in clinically selected cases to enable accurate molecular diagnosis and counsel patients and their families.

Transgenic mice overexpressing mutant TDP-43 show aberrant splicing of neurological disorders-associated gene Zmynd11 prior to onset of motor symptoms.

Mutations in TDP-43 are known to cause Amyotrophic Lateral Sclerosis (ALS) and Frontotemporal Dementia (FTD). TDP-43 binds to and regulates splicing of several RNA including Zmynd11 . Zmynd11 is a transcriptional repressor and a potential E3 ubiquitin ligase family member, known for its role in neuron and muscle differentiation. Mutations in Zmynd11 have been associated with autism with significant developmental motor delays, intellectual disability, and ataxia. Here, we show that Zmynd11 is aberrantly spliced in the brain and spinal cord of transgenic mice overexpressing a mutant human TDP-43 (A315T), and that these changes occur before the onset of motor symptoms.

A humanized yeast model reveals dominant-negative properties of neuropathy-associated alanyl-tRNA synthetase mutations.

Aminoacyl-tRNA synthetases (ARSs) are essential enzymes that ligate tRNA molecules to cognate amino acids. Heterozygosity for missense variants or small in-frame deletions in six ARS genes causes dominant axonal peripheral neuropathy. These pathogenic variants reduce enzyme activity without significantly decreasing protein levels and reside in genes encoding homo-dimeric enzymes. These observations raise the possibility that neuropathy-associated ARS variants exert a dominant-negative effect, reducing overall ARS activity below a threshold required for peripheral nerve function. To test such variants for dominant-negative properties, we developed a humanized yeast assay to co-express pathogenic human alanyl-tRNA synthetase (AARS1) mutations with wild-type human AARS1. We show that multiple loss-of-function AARS1 mutations impair yeast growth through an interaction with wild-type AARS1, but that reducing this interaction rescues yeast growth. This suggests that neuropathy-associated AARS1 variants exert a dominant-negative effect, which supports a common, loss-of-function mechanism for ARS-mediated dominant peripheral neuropathy.

Heterozygous Seryl-tRNA Synthetase 1 Variants Cause Charcot-Marie-Tooth Disease.

OBJECTIVE: Despite the increasing number of genes associated with Charcot-Marie-Tooth (CMT) disease, many patients currently still lack appropriate genetic diagnosis for this disease. Autosomal dominant mutations in aminoacyl-tRNA synthetases (ARSs) have been implicated in CMT. Here, we describe causal missense mutations in the gene encoding seryl-tRNA synthetase 1 (SerRS) for 3 families affected with CMT. METHODS: Whole-exome sequencing was performed in 16 patients and 14 unaffected members of 3 unrelated families. The functional impact of the genetic variants identified was investigated using bioinformatic prediction tools and confirmed using cellular and biochemical assays. RESULTS: Combined linkage analysis for the 3 families revealed significant linkage (Zmax LOD = 6.9) between the genomic co-ordinates on chromosome 1: 108681600-110300504. Within the linkage region, heterozygous SerRS missense variants segregated with the clinical phenotype in the 3 families. The mutant SerRS proteins exhibited reduced aminoacylation activity and abnormal SerRS dimerization, which suggests the impairment of total protein synthesis and induction of eIF2α phosphorylation. INTERPRETATION: Our findings suggest the heterozygous SerRS variants identified represent a novel cause for autosomal dominant CMT. Mutant SerRS proteins are known to impact various molecular and cellular functions. Our findings provide significant advances on the current understanding of the molecular mechanisms associated with ARS-related CMT. ANN NEUROL 2023;93:244-256.

Novel gene-intergenic fusion involving ubiquitin E3 ligase UBE3C causes distal hereditary motor neuropathy.

Distal hereditary motor neuropathies (dHMNs) are a group of inherited diseases involving the progressive, length-dependent axonal degeneration of the lower motor neurons. There are currently 29 reported causative genes and four disease loci implicated in dHMN. Despite the high genetic heterogeneity, mutations in the known genes account for less than 20% of dHMN cases, with the mutations identified predominantly being point mutations or indels. We have expanded the spectrum of dHMN mutations with the identification of a 1.35 Mb complex structural variation (SV) causing a form of autosomal dominant dHMN (DHMN1 OMIM %182906). Given the complex nature of SV mutations and the importance of studying pathogenic mechanisms in a neuronal setting, we generated a patient-derived DHMN1 motor neuron model harbouring the 1.35 Mb complex insertion. The DHMN1 complex insertion creates a duplicated copy of the first 10 exons of the ubiquitin-protein E3 ligase gene (UBE3C) and forms a novel gene-intergenic fusion sense transcript by incorporating a terminal pseudo-exon from intergenic sequence within the DHMN1 locus. The UBE3C intergenic fusion (UBE3C-IF) transcript does not undergo nonsense-mediated decay and results in a significant reduction of wild-type full-length UBE3C (UBE3C-WT) protein levels in DHMN1 iPSC-derived motor neurons. An engineered transgenic Caenorhabditis elegans model expressing the UBE3C-IF transcript in GABA-ergic motor neurons shows neuronal synaptic transmission deficits. Furthermore, the transgenic animals are susceptible to heat stress, which may implicate defective protein homeostasis underlying DHMN1 pathogenesis. Identification of the novel UBE3C-IF gene-intergenic fusion transcript in motor neurons highlights a potential new disease mechanism underlying axonal and motor neuron degeneration. These complementary models serve as a powerful paradigm for studying the DHMN1 complex SV and an invaluable tool for defining therapeutic targets for DHMN1.

Mutations in MYO9B are associated with Charcot-Marie-Tooth disease type 2 neuropathies and isolated optic atrophy.

BACKGROUND AND PURPOSE: Charcot-Marie-Tooth disease (CMT) is a heterogeneous group of disorders caused by mutations in at least 100 genes. However, approximately 60% of cases with axonal neuropathies (CMT2) still remain without a genetic diagnosis. We aimed at identifying novel disease genes responsible for CMT2. METHODS: We performed whole exome sequencing and targeted next generation sequencing panel analyses on a cohort of CMT2 families with evidence for autosomal recessive inheritance. We also performed functional studies to explore the pathogenetic role of selected variants. RESULTS: We identified rare, recessive variants in the MYO9B (myosin IX) gene in two families with CMT2. MYO9B has not yet been associated with a human disease. MYO9B is an unconventional single-headed processive myosin motor protein with signaling properties, and, consistent with this, our results indicate that a variant occurring in the MYO9B motor domain impairs protein expression level and motor activity. Interestingly, a Myo9b-null mouse has degenerating axons in sciatic nerves and optic nerves, indicating that MYO9B plays an essential role in both peripheral nervous system and central nervous system axons, respectively. The degeneration observed in the optic nerve prompted us to screen for MYO9B mutations in a cohort of patients with optic atrophy (OA). Consistent with this, we found compound heterozygous variants in one case with isolated OA. CONCLUSIONS: Novel or very rare variants in MYO9B are associated with CMT2 and isolated OA.

Identity-by-descent analysis of CMTX3 links three families through a common founder.

A large 78 kb insertion from chromosome 8q24.3 into Xq27.1 was identified as the cause of CMTX3 in three families of European descent from Australia (CMT193, CMT180) and New Zealand/United Kingdom (CMT623). Using the relatedness tool XIBD to perform genome-wide identity-by-descent (IBD) analysis on 16 affected individuals from the three families demonstrated they all share the CMTX3 disease locus identical-by-descent, confirming the mutation arose in a common ancestor. Relationship estimation from IBD segment data has genetically linked all three families through 6th and 7th degree relatives.

Whole exome sequencing identifies two novel variants in PHEX and DMP1 in Malaysian children with hypophosphatemic rickets.

BACKGROUND: Hypophosphatemic rickets (HR) is a genetic disease of phosphate wasting that is characterized by defective bone mineralization. The most common cause of the disease is mutations in the phosphate regulating gene with homologies to endopeptidases on the X chromosome (PHEX) gene. The aims of this study were to identify the gene variants responsible for HR in three cases of Malaysian origin from three independent families and to describe their clinical, biochemical, and radiological features. METHODS: Whole exome sequencing (WES) was performed on all patients and their parents, followed by Sanger sequencing validation. Bioinformatics tools were used to provide supporting evidence for pathogenicity of variants. To confirm that a mutation is de novo, paternity test was carried out. High resolution melting curve analysis was performed to assess the allele frequency in normal controls for mutations that were found in the patients. RESULTS: The patients showed typical characteristics of HR including lower limb deformity, hypophosphatemia, and elevated alkaline phosphatase. WES revealed two variants in the PHEX gene and one variant in the dentin matrix protein 1 (DMP1) gene. Two of the three variants were novel, including c.1946_1954del (p.Gly649_Arg651del) in PHEX and c.54 + 1G > A in DMP1. Our data suggests that the novel p.Gly649_Arg651del variant is likely pathogenic for HR disease. CONCLUSIONS: This study extends the variant spectrum of the PHEX and DMP1 genes. Our findings indicate that WES is an advantageous approach for diagnosis of genetic diseases which are heterogeneous.

Genetic analysis of failed male puberty using whole exome sequencing.

OBJECTIVES: Although at least 598 genes are involved in the development of the hypothalamo-pituitary-testicular (HPT) axis, mutations in only 75 genes have so far been shown to cause delayed puberty. METHODS: Six male patients with failed puberty, manifested as absence of pubertal changes by 18 years of age, underwent whole exome sequencing of genomic DNA with subsequent bioinformatics analysis and confirmation of selected variants by Sanger sequencing. Genes having plausibly pathogenic non-synonymous variants were characterized as group A (previously reported to cause delayed puberty), group B (expressed in the HPT-axis but no mutations therein were reported to cause delayed puberty) or group C (not reported previously to be connected with HPT-axis). RESULTS: We identified variants in genes involved in GnRH neuron differentiation (2 in group A, 1 in group C), GnRH neuron migration (2 each in groups A and C), development of GnRH neural connections with supra-hypothalamic and hypothalamic neurons (2 each in groups A and C), neuron homeostasis (1 in group C), molecules regulating GnRH neuron activity (2 each in groups B and C), receptors/proteins expressed on GnRH neurons (1 in group B), signaling molecules (3 in group C), GnRH synthesis (1 in group B), gonadotropins production and release (1 each in groups A, B, and C) and action of the steroid hormone (1 in group A). CONCLUSIONS: Non-synonymous variants were identified in 16 genes of the HPT-axis, which comprised 4 in group A that contains genes previously reported to cause delayed puberty, 4 in group B that are expressed along HPT-axis but no mutations therein were reported previously to cause delayed puberty and 8 in group C that contains novel candidate genes, suggesting wider genetic causes of failed male puberty.