A Biallelic Variant of the RNA Exosome Gene, EXOSC4, Associated with Neurodevelopmental Defects Impairs RNA Exosome Function and Translation.

The RNA exosome is an evolutionarily conserved complex required for both precise RNA processing and decay. Pathogenic variants in EXOSC genes, which encode structural subunits of this complex, are linked to several autosomal recessive disorders. Here, we describe a missense allele of the EXOSC4 gene that causes a collection of clinical features in two affected siblings. This missense variant (NM_019037.3: exon3:c.560T>C), changes a leucine residue within a conserved region of EXOSC4 to proline (p.Leu187Pro). The two affected individuals show prenatal growth restriction, failure to thrive, global developmental delay, intracerebral and basal ganglia calcifications, and kidney failure. Homozygosity for the damaging variant was identified by exome sequencing with Sanger sequencing to confirm segregation. To explore the functional consequences of this amino acid change, we modeled EXOSC4-L187P in the corresponding budding yeast protein, Rrp41 (Rrp41-L187P). Cells that express Rrp41-L187P as the sole copy of the essential Rrp41 protein show growth defects. Steady-state levels of both Rrp41-L187P and EXOSC4-L187P are decreased compared to controls and EXOSC4-L187P shows decreased co-purification with other RNA exosome subunits. RNA exosome target transcripts accumulate in rrp41-L187P cells, including the 7S precursor of 5.8S rRNA. Polysome profiles show a decrease in actively translating ribosomes in rrp41-L187P cells as compared to control cells with incorporation of 7S pre-rRNA into polysomes. This work adds EXOSC4 to the structural subunits of the RNA exosome that have been linked to human disease and defines foundational molecular defects that could contribute to the adverse phenotypes caused by EXOSC pathogenic variants.

Recessive PRDM13 mutations cause fatal perinatal brainstem dysfunction with cerebellar hypoplasia and disrupt Purkinje cell differentiation.

Pontocerebellar hypoplasias (PCHs) are congenital disorders characterized by hypoplasia or early atrophy of the cerebellum and brainstem, leading to a very limited motor and cognitive development. Although over 20 genes have been shown to be mutated in PCHs, a large proportion of affected individuals remains undiagnosed. We describe four families with children presenting with severe neonatal brainstem dysfunction and pronounced deficits in cognitive and motor development associated with four different bi-allelic mutations in PRDM13, including homozygous truncating variants in the most severely affected individuals. Brain MRI and fetopathological examination revealed a PCH-like phenotype, associated with major hypoplasia of inferior olive nuclei and dysplasia of the dentate nucleus. Notably, histopathological examinations highlighted a sparse and disorganized Purkinje cell layer in the cerebellum. PRDM13 encodes a transcriptional repressor known to be critical for neuronal subtypes specification in the mouse retina and spinal cord but had not been implicated, so far, in hindbrain development. snRNA-seq data mining and in situ hybridization in humans show that PRDM13 is expressed at early stages in the progenitors of the cerebellar ventricular zone, which gives rise to cerebellar GABAergic neurons, including Purkinje cells. We also show that loss of function of prdm13 in zebrafish leads to a reduction in Purkinje cells numbers and a complete absence of the inferior olive nuclei. Altogether our data identified bi-allelic mutations in PRDM13 as causing a olivopontocerebellar hypoplasia syndrome and suggest that early deregulations of the transcriptional control of neuronal fate specification could contribute to a significant number of cases.

THUMPD1 bi-allelic variants cause loss of tRNA acetylation and a syndromic neurodevelopmental disorder.

Covalent tRNA modifications play multi-faceted roles in tRNA stability, folding, and recognition, as well as the rate and fidelity of translation, and other cellular processes such as growth, development, and stress responses. Mutations in genes that are known to regulate tRNA modifications lead to a wide array of phenotypes and diseases including numerous cognitive and neurodevelopmental disorders, highlighting the critical role of tRNA modification in human disease. One such gene, THUMPD1, is involved in regulating tRNA N4-acetylcytidine modification (ac4C), and recently was proposed as a candidate gene for autosomal-recessive intellectual disability. Here, we present 13 individuals from 8 families who harbor rare loss-of-function variants in THUMPD1. Common phenotypic findings included global developmental delay, speech delay, moderate to severe intellectual deficiency, behavioral abnormalities such as angry outbursts, facial dysmorphism, and ophthalmological abnormalities. We demonstrate that the bi-allelic variants identified cause loss of function of THUMPD1 and that this defect results in a loss of ac4C modification in small RNAs, and of individually purified tRNA-Ser-CGA. We further corroborate this effect by showing a loss of tRNA acetylation in two CRISPR-Cas9-generated THUMPD1 KO cell lines. In addition, we also show the resultant amino acid substitution that occurs in a missense THUMPD1 allele identified in an individual with compound heterozygous variants results in a marked decrease in THUMPD1 stability and RNA-binding capacity. Taken together, these results suggest that the lack of tRNA acetylation due to THUMPD1 loss of function results in a syndromic form of intellectual disability associated with developmental delay, behavioral abnormalities, hearing loss, and facial dysmorphism.

Biallelic variants of the first Kunitz domain of SPINT2 cause a non-syndromic form of congenital diarrhea and tufting enteropathy.

Biallelic SPINT2 pathogenic variants cause a syndromic form of congenital diarrhea and enteropathy (OMIM 270420). To date, 35 patients have been reported and all presented with additional extra-intestinal features, apart from one case. We report on a 5-year-old girl who presented early in life with diarrhea and was found to have a novel homozygous variant in SPINT2. Pathological studies confirmed tufting enteropathy, and during her 5 years of life, she has not developed any extra-intestinal features. Molecular analysis detected a homozygous variant (NM_021102.4: c.203A>G (p. [Tyr68Cys]) in SPINT2. This is the first missense variant reported in the first Kunitz domain (KD1) of SPINT2 in humans. In vitro functional studies of this variant confirmed the deleterious effect leading to the loss of inhibitory activity of the intestinal serine proteases. This is the first description of SPINT2-related diarrhea in a patient who lived without long-term total parenteral nutrition. This study expands the clinical and molecular characteristics of SPINT2-related conditions.

A homozygous loss-of-function C1S mutation is associated with Kikuchi-Fujimoto disease.

BACKGROUND: Kikuchi-Fujimoto disease (KFD) is a self-limited inflammatory disease of unknown pathogenesis. Familial cases have been described and defects in classical complement components C1q and C4 have been identified in some patients. MATERIAL AND METHODS: We describe genetic and immune investigations of a 16 years old Omani male, a product of consanguineous marriage, who presented with typical clinical and histological features of KFD. RESULTS: We identified a novel homozygous single base deletion in C1S (c.330del; p. Phe110LeufsTer23) resulting in a defect in the classical complement pathway. The patient was negative for all serological markers of SLE. In contrast, two female siblings (also homozygous for the C1S mutation), one has autoimmune thyroid disease (Hashimoto thyroiditis) and a positive ANA and the other sibling has serology consistent with SLE. CONCLUSION: We report the first association between C1s deficiency and KFD.

Homozygous loss-of-function variants in FILIP1 cause autosomal recessive arthrogryposis multiplex congenita with microcephaly.

Arthrogryposis multiplex congenita forms a broad group of clinically and etiologically heterogeneous disorders characterized by congenital joint contractures that involve at least two different parts of the body. Neurological and muscular disorders are commonly underlying arthrogryposis. Here, we report five affected individuals from three independent families sharing an overlapping phenotype with congenital contractures affecting shoulder, elbow, hand, hip, knee and foot as well as scoliosis, reduced palmar and plantar skin folds, microcephaly and facial dysmorphism. Using exome sequencing, we identified homozygous truncating variants in FILIP1 in all patients. FILIP1 is a regulator of filamin homeostasis required for the initiation of cortical cell migration in the developing neocortex and essential for the differentiation process of cross-striated muscle cells during myogenesis. In summary, our data indicate that bi-allelic truncating variants in FILIP1 are causative of a novel autosomal recessive disorder and expand the spectrum of genetic factors causative of arthrogryposis multiplex congenita.

Recessive, Deleterious Variants in SMG8 Expand the Role of Nonsense-Mediated Decay in Developmental Disorders in Humans.

We have previously described a heart-, eye-, and brain-malformation syndrome caused by homozygous loss-of-function variants in SMG9, which encodes a critical component of the nonsense-mediated decay (NMD) machinery. Here, we describe four consanguineous families with four different likely deleterious homozygous variants in SMG8, encoding a binding partner of SMG9. The observed phenotype greatly resembles that linked to SMG9 and comprises severe global developmental delay, microcephaly, facial dysmorphism, and variable congenital heart and eye malformations. RNA-seq analysis revealed a general increase in mRNA expression levels with significant overrepresentation of core NMD substrates. We also identified increased phosphorylation of UPF1, a key SMG1-dependent step in NMD, which most likely represents the loss of SMG8--mediated inhibition of SMG1 kinase activity. Our data show that SMG8 and SMG9 deficiency results in overlapping developmental disorders that most likely converge mechanistically on impaired NMD.

Genetic Causes, Clinical Features, and Survival of Underlying Inborn Errors of Immunity in Omani Patients: a Single-Center Study.

PURPOSE: Early identification of inborn errors of immunity (IEIs) is crucial due to the significant risk of morbidity and mortality. This study aimed to describe the genetic causes, clinical features, and survival rate of IEIs in Omani patients. METHODS: A prospective study of all Omani patients evaluated for immunodeficiency was conducted over a 17-year period. Clinical features and diagnostic immunological findings were recorded. Targeted gene testing was performed in cases of obvious immunodeficiency. For cases with less conclusive phenotypes, a gene panel was performed, followed by whole-exome sequencing if necessary. RESULTS: A total of 185 patients were diagnosed with IEIs during the study period; of these, 60.5% were male. Mean ages at symptom onset and diagnosis were 30.0 and 50.5 months, respectively. Consanguinity and a family history of IEIs were present in 86.9% and 50.8%, respectively. Most patients presented with lower respiratory infections (65.9%), followed by growth and development manifestations (43.2%). Phagocytic defects were the most common cause of IEIs (31.9%), followed by combined immunodeficiency (21.1%). Overall, 109 of 132 patients (82.6%) who underwent genetic testing received a genetic diagnosis, while testing was inconclusive for the remaining 23 patients (17.4%). Among patients with established diagnoses, 37 genes and 44 variants were identified. Autosomal recessive inheritance was present in 81.7% of patients with gene defects. Several variants were novel. Intravenous immunoglobulin therapy was administered to 39.4% of patients and 21.6% received hematopoietic stem cell transplantation. The overall survival rate was 75.1%. CONCLUSION: This study highlights the genetic causes of IEIs in Omani patients. This information may help in the early identification and management of the disease, thereby improving survival and quality of life.

Inactivation of DRG1, encoding a translation factor GTPase, causes a recessive neurodevelopmental disorder.

PURPOSE: Developmentally regulated Guanosine-5'-triphosphate-binding protein 1 (DRG1) is a highly conserved member of a class of GTPases implicated in translation. Although the expression of mammalian DRG1 is elevated in the central nervous system during development, and its function has been implicated in fundamental cellular processes, no pathogenic germline variants have yet been identified. Here, we characterize the clinical and biochemical consequences of DRG1 variants. METHODS: We collate clinical information of 4 individuals with germline DRG1 variants and use in silico, in vitro, and cell-based studies to study the pathogenicity of these alleles. RESULTS: We identified private germline DRG1 variants, including 3 stop-gained p.Gly54∗, p.Arg140∗, p.Lys263∗, and a p.Asn248Phe missense variant. These alleles are recessively inherited in 4 affected individuals from 3 distinct families and cause a neurodevelopmental disorder with global developmental delay, primary microcephaly, short stature, and craniofacial anomalies. We show that these loss-of-function variants (1) severely disrupt DRG1 messenger RNA/protein stability in patient-derived fibroblasts, (2) impair its GTPase activity, and (3) compromise its binding to partner protein ZC3H15. Consistent with the importance of DRG1 in humans, targeted inactivation of mouse Drg1 resulted in preweaning lethality. CONCLUSION: Our work defines a new Mendelian disorder of DRG1 deficiency. This study highlights DRG1's importance for normal mammalian development and underscores the significance of translation factor GTPases in human physiology and homeostasis.

Clinical, neuroradiological, and molecular characterization of mitochondrial threonyl-tRNA-synthetase (TARS2)-related disorder.

PURPOSE: Biallelic variants in TARS2, encoding the mitochondrial threonyl-tRNA-synthetase, have been reported in a small group of individuals displaying a neurodevelopmental phenotype but with limited neuroradiological data and insufficient evidence for causality of the variants. METHODS: Exome or genome sequencing was carried out in 15 families. Clinical and neuroradiological evaluation was performed for all affected individuals, including review of 10 previously reported individuals. The pathogenicity of TARS2 variants was evaluated using in vitro assays and a zebrafish model. RESULTS: We report 18 new individuals harboring biallelic TARS2 variants. Phenotypically, these individuals show developmental delay/intellectual disability, regression, cerebellar and cerebral atrophy, basal ganglia signal alterations, hypotonia, cerebellar signs, and increased blood lactate. In vitro studies showed that variants within the TARS2301-381 region had decreased binding to Rag GTPases, likely impairing mTORC1 activity. The zebrafish model recapitulated key features of the human phenotype and unraveled dysregulation of downstream targets of mTORC1 signaling. Functional testing of the variants confirmed the pathogenicity in a zebrafish model. CONCLUSION: We define the clinico-radiological spectrum of TARS2-related mitochondrial disease, unveil the likely involvement of the mTORC1 signaling pathway as a distinct molecular mechanism, and establish a TARS2 zebrafish model as an important tool to study variant pathogenicity.

Expanding the spectrum of EEF1D neurodevelopmental disorders: Biallelic variants in the guanine exchange domain.

Protein translation is an essential cellular process and dysfunctional protein translation causes various neurodevelopmental disorders. The eukaryotic translation elongation factor 1A (eEF1A) delivers aminoacyl-tRNA to the ribosome, while the eEF1B complex acts as a guanine exchange factor (GEF) of GTP for GDP indirectly catalyzing the release of eEF1A from the ribosome. The gene EEF1D encodes the eEF1Bδ subunit of the eEF1B complex. EEF1D is alternatively spliced giving rise to one long and three short isoforms. Two different homozygous, truncating variants in EEF1D had been associated with severe intellectual disability and microcephaly in two families. The published variants only affect the long isoform of EEF1D that acts as a transcription factor of heat shock element proteins. By exome sequencing, we identified two different homozygous variants in EEF1D in two families with severe developmental delay, severe microcephaly, spasticity, and failure to thrive with optic atrophy, poor feeding, and recurrent aspiration pneumonia. The EEF1D variants reported in this study are localized in the C-terminal GEF domain, suggesting that a disturbed protein translation machinery might contribute to the neurodevelopmental phenotype. Pathogenic variants localized in both the alternatively spliced domain or the GEF domain of EEF1D cause a severe neurodevelopmental disorder with microcephaly and spasticity.

TMEM63C mutations cause mitochondrial morphology defects and underlie hereditary spastic paraplegia.

The hereditary spastic paraplegias (HSP) are among the most genetically diverse of all Mendelian disorders. They comprise a large group of neurodegenerative diseases that may be divided into 'pure HSP' in forms of the disease primarily entailing progressive lower-limb weakness and spasticity, and 'complex HSP' when these features are accompanied by other neurological (or non-neurological) clinical signs. Here, we identified biallelic variants in the transmembrane protein 63C (TMEM63C) gene, encoding a predicted osmosensitive calcium-permeable cation channel, in individuals with hereditary spastic paraplegias associated with mild intellectual disability in some, but not all cases. Biochemical and microscopy analyses revealed that TMEM63C is an endoplasmic reticulum-localized protein, which is particularly enriched at mitochondria-endoplasmic reticulum contact sites. Functional in cellula studies indicate a role for TMEM63C in regulating both endoplasmic reticulum and mitochondrial morphologies. Together, these findings identify autosomal recessive TMEM63C variants as a cause of pure and complex HSP and add to the growing evidence of a fundamental pathomolecular role of perturbed mitochondrial-endoplasmic reticulum dynamics in motor neurone degenerative diseases.

BCL10 loss-of-function novel mutation leading to atypical severe combined immunodeficiency.

BACKGROUND: Severe combined immunodeficiency (SCID) is characterized by severe, early-onset infection in infants. B-cell lymphoma/leukemia (BCL) 10 defects causing SCID have been reported previously in two patients. MATERIAL & METHODS: A seven-month-old female infant was admitted with bilateral pneumonia requiring ventilatory support. She had a history of recurrent infections starting from four months of age. The patient was investigated for primary immunodeficiency. RESULTS: Immunological investigations revealed hypogammaglobulinemia with normal CD4 and CD8 lymphocyte counts, while a lymphocyte proliferation assay showed absent response to phytohemagglutinin stimulation, thereby establishing the diagnosis of an atypical form of SCID. Genetic testing revealed a homozygous mutation in the BCL10 gene, with both parents demonstrating a heterozygous state (NM_003921.5:c.271A > C:p.[Thr91Pro]). The patient died before bone marrow transplantation due to severe disseminated adenovirus disease. CONCLUSION: We report the first patient from the Middle East with a novel homozygous mutation in the BCL10 gene causing SCID.

Further phenotypic delineation of Alazami syndrome.

Alazami syndrome (AS) is an autosomal recessive condition characterized by the cardinal features of severe growth restriction, moderate to severe intellectual disability, and distinctive facial features. Biallelic pathogenic variants of the LARP7, encoding a chaperone of 7SK noncoding RNA, is implicated in this disease. There are <35 reported cases in the literature. All reported cases share the same three cardinal features of the syndrome. Herein, we report on 12 patients with a confirmed diagnosis of AS from eight unrelated families. The cohort shares the same key feature of the syndrome. Moreover, we report additional phenotypic features, including genito-renal anomalies, ophthalmological abnormalities, and congenital heart disease. Whole-exome sequencing was used in all reported cases, implicating a clinical under-recognition of the syndrome. This report further expands the clinical and molecular characteristics of Alazami syndrome.

Deficiency of acyl-CoA synthetase 5 is associated with a severe and treatable failure to thrive of neonatal onset.

Failure to thrive (FTT) causes significant morbidity, often without clear etiologies. Six individuals of a large consanguineous family presented in the neonatal period with recurrent vomiting and diarrhea, leading to severe FTT. Standard diagnostic work up did not ascertain an etiology. Autozygosity mapping and whole exome sequencing identified homozygosity for a novel genetic variant of the long chain fatty acyl-CoA synthetase 5 (ACSL5) shared among the affected individuals (NM_203379.1:c.1358C>A:p.(Thr453Lys)). Autosomal recessive genotype-phenotype segregation was confirmed by Sanger sequencing. Functional in vitro analysis of the ACSL5 variant by immunofluorescence, western blotting and enzyme assay suggested that Thr453Lys is a loss-of-function mutation without any remaining activity. ACSL5 belongs to an essential enzyme family required for lipid metabolism and is known to contribute the major activity in the mouse intestine. Based on the function of ACSL5 in intestinal long chain fatty acid metabolism and the gastroenterological symptoms, affected individuals were treated with total parenteral nutrition or medium-chain triglyceride-based formula restricted in long-chain triglycerides. The patients responded well and follow up suggests that treatment is only required during early life.

Recessive mutations in SCYL2 cause a novel syndromic form of arthrogryposis in humans.

Arthrogryposis multiplex congenita (AMC) is an important birth defect with a significant genetic contribution. Many syndromic forms of AMC have been described, but remain unsolved at the molecular level. In this report, we describe a novel syndromic form of AMC in two multiplex consanguineous families from Saudi Arabia and Oman. The phenotype is highly consistent, and comprises neurogenic arthrogryposis, microcephaly, brain malformation (absent corpus callosum), optic atrophy, limb fractures, profound global developmental delay, and early lethality. Whole-exome sequencing revealed a different homozygous truncating variant in SCYL2 in each of the two families. SCYL2 is a component of clathrin-coated vesicles, and deficiency of its mouse ortholog results in a severe neurological phenotype that largely recapitulates the phenotype observed in our patients. Our results suggest that severe neurogenic arthrogryposis with brain malformation is the human phenotypic consequence of SCYL2 loss of function mutations.

Phenotypic spectrum of ALPK3-related cardiomyopathy.

Cardiomyopathies are clinically heterogeneous disorders and are the leading cause of cardiovascular morbidity and mortality. Different etiologies have a significant impact on prognosis. Recently, novel biallelic loss-of-function pathogenic variants in alpha-kinase 3 (ALPK3) were implicated in causing early-onset pediatric cardiomyopathy (cardiomyopathy, familial hypertrophic 27; OMIM 618052). To date, eight patients, all presented during early childhood, were reported with biallelic ALPK3 pathogenic variants. We describe the molecular and clinical phenotype characterization of familial cardiomyopathy on one family with six affected individuals. We identified homozygosity for an ALPK3 deleterious sequence variant (NM_020778.4:c.639G>A:p.Trp213*) in all the affected individuals. They presented with either dilated cardiomyopathy that progressed to hypertrophic cardiomyopathy (HCM) or HCM with left ventricular noncompaction. The age of presentation in our cohort extends between infancy to the fourth decade. The phenotypic severity decreases with the progression of age.

Homozygosity for FARSB mutation leads to Phe-tRNA synthetase-related disease of growth restriction, brain calcification, and interstitial lung disease.

Aminoacyl-tRNA synthetases (ARSs) canonical function is to conjugate specific amino acids to cognate tRNA that are required for the first step of protein synthesis. Genetic mutations that cause dysfunction or absence of ARSs result in various neurodevelopmental disorders. The human phenylalanine-tRNA synthetase (PheRS) is a tetrameric protein made of two subunits coded by FARSA gene and two subunits coded by FARSB gene. We describe eight affected individuals from an extended family with a multisystemic recessive disease manifest as a significant growth restriction, brain calcifications, and interstitial lung disease. Genome-wide linkage analysis and whole exome sequencing identified homozygosity for a FARSB mutation (NM_005687.4:c.853G > A:p.Glu285Lys) that co-segregate with the disease and likely cause loss-of-function. This study further implicates FARSB mutations in a multisystem, recessive, neurodevelopmental phenotype that share clinical features with the previously known aminoacyl-tRNA synthetase-related diseases.

Mutations in PYCR2, Encoding Pyrroline-5-Carboxylate Reductase 2, Cause Microcephaly and Hypomyelination.

Despite recent advances in understanding the genetic bases of microcephaly, a large number of cases of microcephaly remain unexplained, suggesting that many microcephaly syndromes and associated genes have yet to be identified. Here, we report mutations in PYCR2, which encodes an enzyme in the proline biosynthesis pathway, as the cause of a unique syndrome characterized by postnatal microcephaly, hypomyelination, and reduced cerebral white-matter volume. Linkage mapping and whole-exome sequencing identified homozygous mutations (c.355C>T [p.Arg119Cys] and c.751C>T [p.Arg251Cys]) in PYCR2 in the affected individuals of two consanguineous families. A lymphoblastoid cell line from one affected individual showed a strong reduction in the amount of PYCR2. When mutant cDNAs were transfected into HEK293FT cells, both variant proteins retained normal mitochondrial localization but had lower amounts than the wild-type protein, suggesting that the variant proteins were less stable. A PYCR2-deficient HEK293FT cell line generated by genome editing with the clustered regularly interspaced short palindromic repeat (CRISPR)-Cas9 system showed that PYCR2 loss of function led to decreased mitochondrial membrane potential and increased susceptibility to apoptosis under oxidative stress. Morpholino-based knockdown of a zebrafish PYCR2 ortholog, pycr1b, recapitulated the human microcephaly phenotype, which was rescued by wild-type human PYCR2 mRNA, but not by mutant mRNAs, further supporting the pathogenicity of the identified variants. Hypomyelination and the absence of lax, wrinkly skin distinguishes this condition from that caused by previously reported mutations in the gene encoding PYCR2's isozyme, PYCR1, suggesting a unique and indispensable role for PYCR2 in the human CNS during development.

Reanalysis of exome sequencing data of intellectual disability samples: Yields and benefits.

Recently, with the advancement in next generation sequencing (NGS) along with the improvement of bioinformatics tools, whole exome sequencing (WES) has become the most efficient diagnostic test for patients with intellectual disability (ID). This study aims to estimate the yield of a reanalysis of ID negative exome cases after data reannotation. Total of 50 data files of exome sequencing, representing 50 samples were collected. The inclusion criteria include ID phenotype, and previous analysis indicated a negative result (no abnormality detected). These files were pre-processed and reannotated using ANNOVAR tool. Prioritized variants in the 50 cases studied were classified into three groups, (1) disease-causative variants (2) possible disease-causing variants and (3) variants in novel genes. Reanalysis resulted in the identification of pathogenic/likely pathogenic variants in six cases (12%). Thirteen cases (26%) were classified as having possible disease-causing variants. Candidate genes requiring future functional studies were detected in seven cases (14%). Improvement in bioinformatics tools, update in the genetic databases and literature, and patients' clinical phenotype update were the main reasons for identification of these variants in this study.