Clinical exome sequencing by general pediatricians: high clinical utility and no evidence of inappropriate testing.

Background: Genetic disorders account for a large percentage of admissions and outpatient visits to children's hospitals around the world. Clinical exome sequencing (CES) is a valuable diagnostic tool in the workup of these disorders; however, it is not routinely requested by general pediatricians. This may represent a missed opportunity to increase patient access to this powerful diagnostic tool. In our institution, general pediatricians can directly order CES. In this context, this study aims to evaluate the appropriateness of CES and its clinical utility when ordered by general pediatricians. Methods: We retrospectively reviewed all CES tests ordered by general pediatricians in our institution between 2019 and 2023 and recorded their indications and results. General pediatricians were interviewed to evaluate how CES impacted the domains of clinical utility by assessing changes in management, communication, subsequent testing, and counseling. In addition, feedback was obtained, and barriers faced by general pediatricians to order CES were assessed. Results: The study cohort (n = 30) included children from the inpatient (60%) and outpatient (40%) departments. A positive finding (a pathogenic or likely pathogenic variant that explains the phenotype) was observed in 11 of 30 cases (37%), while 3 (10%) and 16 (53%) received ambiguous (variant of uncertain significance) and negative results, respectively. The indication was deemed appropriate in all 30 cases (100%). Clinical utility was reported in all 11 positive cases (100%). Reproductive counseling is a notable utility in this highly consanguineous population, as all variants identified, in the 11 positive cases, were autosomal recessive. Conclusion: We show that CES ordered by general pediatricians is appropriately indicated and provides a diagnostic yield comparable to that requested by specialists. In addition, we note the high clinical utility of positive results as judged by the ordering pediatricians. The findings of this study can empower general pediatricians to advocate for expanded CES adoption to improve patient access and shorten their diagnostic odyssey.

Biallelic NAA60 variants with impaired n-terminal acetylation capacity cause autosomal recessive primary familial brain calcifications.

Primary familial brain calcification (PFBC) is characterized by calcium deposition in the brain, causing progressive movement disorders, psychiatric symptoms, and cognitive decline. PFBC is a heterogeneous disorder currently linked to variants in six different genes, but most patients remain genetically undiagnosed. Here, we identify biallelic NAA60 variants in ten individuals from seven families with autosomal recessive PFBC. The NAA60 variants lead to loss-of-function with lack of protein N-terminal (Nt)-acetylation activity. We show that the phosphate importer SLC20A2 is a substrate of NAA60 in vitro. In cells, loss of NAA60 caused reduced surface levels of SLC20A2 and a reduction in extracellular phosphate uptake. This study establishes NAA60 as a causal gene for PFBC, provides a possible biochemical explanation of its disease-causing mechanisms and underscores NAA60-mediated Nt-acetylation of transmembrane proteins as a fundamental process for healthy neurobiological functioning.

Genomic analysis of presumed perinatal stroke in Saudi Arabia reveals a strong monogenic contribution.

Perinatal stroke is associated with significant short- and long-term morbidity and has been recognized as the most common cause of cerebral palsy in term infants. The diagnosis of presumed perinatal stroke (PPS) is made in children who present with neurological deficit and/or seizures attributable to focal chronic infarction on neuroimaging and have uneventful neonatal history. The underlying mechanism of presumed perinatal stroke remains unknown and thorough investigation of potential monogenic causes has not been conducted to date. Here, we describe the use of untargeted exome sequencing to investigate a cohort of eight patients from six families with PPS. A likely deleterious variant was identified in four families. These include the well-established risk genes COL4A2 and JAM3. In addition, we report the first independent confirmation of the recently described link between ESAM and perinatal stroke. Our data also highlight NID1 as a candidate gene for the condition. This study suggests that monogenic disorders are important contributors to the pathogenesis of PPS and should be investigated by untargeted sequencing especially when traditional risk factors are excluded.

Further delineation of the phenotypic and metabolomic profile of ALDH1L2-related neurodevelopmental disorder.

ALDH1L2, a mitochondrial enzyme in folate metabolism, converts 10-formyl-THF (10-formyltetrahydrofolate) to THF (tetrahydrofolate) and CO2. At the cellular level, deficiency of this NADP+-dependent reaction results in marked reduction in NADPH/NADP+ ratio and reduced mitochondrial ATP. Thus far, a single patient with biallelic ALDH1L2 variants and the phenotype of a neurodevelopmental disorder has been reported. Here, we describe another patient with a neurodevelopmental disorder associated with a novel homozygous missense variant in ALDH1L2, Pro133His. The variant caused marked reduction in the ALDH1L2 enzyme activity in skin fibroblasts derived from the patient as probed by 10-FDDF, a stable synthetic analog of 10-formyl-THF. Additional associated abnormalities in these fibroblasts include reduced NADPH/NADP+ ratio and pool of mitochondrial ATP, upregulated autophagy and dramatically altered metabolomic profile. Overall, our study further supports a link between ALDH1L2 deficiency and abnormal neurodevelopment in humans.

Beyond the exome: utility of long-read whole genome sequencing in exome-negative autosomal recessive diseases.

BACKGROUND: Long-read whole genome sequencing (lrWGS) has the potential to address the technical limitations of exome sequencing in ways not possible by short-read WGS. However, its utility in autosomal recessive Mendelian diseases is largely unknown. METHODS: In a cohort of 34 families in which the suspected autosomal recessive diseases remained undiagnosed by exome sequencing, lrWGS was performed on the Pacific Bioscience Sequel IIe platform. RESULTS: Likely causal variants were identified in 13 (38%) of the cohort. These include (1) a homozygous splicing SV in TYMS as a novel candidate gene for lethal neonatal lactic acidosis, (2) a homozygous non-coding SV that we propose impacts STK25 expression and causes a novel neurodevelopmental disorder, (3) a compound heterozygous SV in RP1L1 with complex inheritance pattern in a family with inherited retinal disease, (4) homozygous deep intronic variants in LEMD2 and SNAP91 as novel candidate genes for neurodevelopmental disorders in two families, and (5) a promoter SNV in SLC4A4 causing non-syndromic band keratopathy. Surprisingly, we also encountered causal variants that could have been identified by short-read exome sequencing in 7 families. The latter highlight scenarios that are especially challenging at the interpretation level. CONCLUSIONS: Our data highlight the continued need to address the interpretation challenges in parallel with efforts to improve the sequencing technology itself. We propose a path forward for the implementation of lrWGS sequencing in the setting of autosomal recessive diseases in a way that maximizes its utility.

A founder DBR1 variant causes a lethal form of congenital ichthyosis.

DBR1 encodes the only known human lariat debranching enzyme and its deficiency has been found to cause an autosomal recessive inborn error of immunity characterized by pediatric brainstem viral-induced encephalitis (MIM 619441). We describe a distinct allelic disorder caused by a founder recessive DBR1 variant in four families (DBR1(NM_016216.4):c.200A > G (p.Tyr67Cys)). Consistent features include prematurity, severe intrauterine growth deficiency, congenital ichthyosis-like presentation (collodion membrane, severe skin peeling and xerosis), and death before the first year of life. Patient-derived fibroblasts displayed the characteristic accumulation of intron lariats in their RNA as revealed by targeted and untargeted analysis, in addition to a marked reduction of DBR1 on immunoblot analysis. We propose a novel DBR1-related developmental disorder that is distinct from DBR1-related encephalitis susceptibility and highlight the apparent lack of correlation with the degree of DBR1 deficiency.

Diagnostic implications of pitfalls in causal variant identification based on 4577 molecularly characterized families.

Despite large sequencing and data sharing efforts, previously characterized pathogenic variants only account for a fraction of Mendelian disease patients, which highlights the need for accurate identification and interpretation of novel variants. In a large Mendelian cohort of 4577 molecularly characterized families, numerous scenarios in which variant identification and interpretation can be challenging are encountered. We describe categories of challenges that cover the phenotype (e.g. novel allelic disorders), pedigree structure (e.g. imprinting disorders masquerading as autosomal recessive phenotypes), positional mapping (e.g. double recombination events abrogating candidate autozygous intervals), gene (e.g. novel gene-disease assertion) and variant (e.g. complex compound inheritance). Overall, we estimate a probability of 34.3% for encountering at least one of these challenges. Importantly, our data show that by only addressing non-sequencing-based challenges, around 71% increase in the diagnostic yield can be expected. Indeed, by applying these lessons to a cohort of 314 cases with negative clinical exome or genome reports, we could identify the likely causal variant in 54.5%. Our work highlights the need to have a thorough approach to undiagnosed diseases by considering a wide range of challenges rather than a narrow focus on sequencing technologies. It is hoped that by sharing this experience, the yield of undiagnosed disease programs globally can be improved.

SLC4A10 mutation causes a neurological disorder associated with impaired GABAergic transmission.

SLC4A10 is a plasma-membrane bound transporter that utilizes the Na+ gradient to drive cellular HCO3- uptake, thus mediating acid extrusion. In the mammalian brain, SLC4A10 is expressed in principal neurons and interneurons, as well as in epithelial cells of the choroid plexus, the organ regulating the production of CSF. Using next generation sequencing on samples from five unrelated families encompassing nine affected individuals, we show that biallelic SLC4A10 loss-of-function variants cause a clinically recognizable neurodevelopmental disorder in humans. The cardinal clinical features of the condition include hypotonia in infancy, delayed psychomotor development across all domains and intellectual impairment. Affected individuals commonly display traits associated with autistic spectrum disorder including anxiety, hyperactivity and stereotyped movements. In two cases isolated episodes of seizures were reported in the first few years of life, and a further affected child displayed bitemporal epileptogenic discharges on EEG without overt clinical seizures. While occipitofrontal circumference was reported to be normal at birth, progressive postnatal microcephaly evolved in 7 out of 10 affected individuals. Neuroradiological features included a relative preservation of brain volume compared to occipitofrontal circumference, characteristic narrow sometimes 'slit-like' lateral ventricles and corpus callosum abnormalities. Slc4a10 -/- mice, deficient for SLC4A10, also display small lateral brain ventricles and mild behavioural abnormalities including delayed habituation and alterations in the two-object novel object recognition task. Collapsed brain ventricles in both Slc4a10-/- mice and affected individuals suggest an important role of SLC4A10 in the production of the CSF. However, it is notable that despite diverse roles of the CSF in the developing and adult brain, the cortex of Slc4a10-/- mice appears grossly intact. Co-staining with synaptic markers revealed that in neurons, SLC4A10 localizes to inhibitory, but not excitatory, presynapses. These findings are supported by our functional studies, which show the release of the inhibitory neurotransmitter GABA is compromised in Slc4a10-/- mice, while the release of the excitatory neurotransmitter glutamate is preserved. Manipulation of intracellular pH partially rescues GABA release. Together our studies define a novel neurodevelopmental disorder associated with biallelic pathogenic variants in SLC4A10 and highlight the importance of further analyses of the consequences of SLC4A10 loss-of-function for brain development, synaptic transmission and network properties.

PRSS8, encoding prostasin, is mutated in patients with autosomal recessive ichthyosis.

Ichthyosis is a genetically heterogeneous genodermatosis characterized by severely rough, dry and scaly skin. We report two consanguineous families with congenital ichthyosis. Combined positional mapping and exome sequencing of the two families revealed novel homozygous likely deleterious variants in PRSS8 (encoding prostasin) within a linkage locus on chromosome 16. One variant involved a canonical splice site and was associated with reduced abundance of the normal transcript, while the other was a missense variant that altered a highly conserved residue. The phenotype of Prss8 knockout mouse bears a striking resemblance to the one we describe in human patients, including the skin histopathology. Our data suggest a novel PRSS8-related ichthyosis disorder.

Mitochondrial "dysmorphology" in variant classification.

Mitochondrial disorders are challenging to diagnose. Exome sequencing has greatly enhanced the diagnostic precision of these disorders although interpreting variants of uncertain significance (VUS) remains a formidable obstacle. Whether specific mitochondrial morphological changes can aid in the classification of these variants is unknown. Here, we describe two families (four patients), each with a VUS in a gene known to affect the morphology of mitochondria through a specific role in the fission-fusion balance. In the first, the missense variant in MFF, encoding a fission factor, was associated with impaired fission giving rise to a characteristically over-tubular appearance of mitochondria. In the second, the missense variant in DNAJA3, which has no listed OMIM phenotype, was associated with fragmented appearance of mitochondria consistent with its published deficiency states. In both instances, the highly specific phenotypes allowed us to upgrade the classification of the variants. Our results suggest that, in select cases, mitochondrial "dysmorphology" can be helpful in interpreting variants to reach a molecular diagnosis.

Lethal variants in humans: lessons learned from a large molecular autopsy cohort.

BACKGROUND: Molecular autopsy refers to DNA-based identification of the cause of death. Despite recent attempts to broaden its scope, the term remains typically reserved to sudden unexplained death in young adults. In this study, we aim to showcase the utility of molecular autopsy in defining lethal variants in humans. METHODS: We describe our experience with a cohort of 481 cases in whom the cause of premature death was investigated using DNA from the index or relatives (molecular autopsy by proxy). Molecular autopsy tool was typically exome sequencing although some were investigated using targeted approaches in the earlier stages of the study; these include positional mapping, targeted gene sequencing, chromosomal microarray, and gene panels. RESULTS: The study includes 449 cases from consanguineous families and 141 lacked family history (simplex). The age range was embryos to 18 years. A likely causal variant (pathogenic/likely pathogenic) was identified in 63.8% (307/481), a much higher yield compared to the general diagnostic yield (43%) from the same population. The predominance of recessive lethal alleles allowed us to implement molecular autopsy by proxy in 55 couples, and the yield was similarly high (63.6%). We also note the occurrence of biallelic lethal forms of typically non-lethal dominant disorders, sometimes representing a novel bona fide biallelic recessive disease trait. Forty-six disease genes with no OMIM phenotype were identified in the course of this study. The presented data support the candidacy of two other previously reported novel disease genes (FAAH2 and MSN). The focus on lethal phenotypes revealed many examples of interesting phenotypic expansion as well as remarkable variability in clinical presentation. Furthermore, important insights into population genetics and variant interpretation are highlighted based on the results. CONCLUSIONS: Molecular autopsy, broadly defined, proved to be a helpful clinical approach that provides unique insights into lethal variants and the clinical annotation of the human genome.

PLXNA2 as a candidate gene in patients with intellectual disability.

Intellectual disability (ID) is one of the most common disabilities in humans. In an effort to contribute to the expanding genetic landscape of ID, we describe a novel autosomal recessive ID candidate gene. Combined autozygome/exome analysis was performed in two unrelated consanguineous families with ID. Each of the two families had a novel homozygous likely deleterious variant in PLXNA2 and displayed the core phenotype of ID. PLXNA2 belongs to a family of transmembrane proteins that function as semaphorin receptors. Sema5A-PlexinA2 is known to regulate brain development in mouse, and Plxna2-/- mice display defective associative learning, sociability, and sensorimotor gating. We note the existence of variability in the phenotype among the three patients, including the existence of variable degree of ID, ranging from borderline intellectual functioning to moderate-severe ID, and the presence of cardiac anomalies in only one of the patients. We propose incomplete penetrance as a possible explanation of the observed difference in phenotypes. Future cases will be needed to support the proposed link between PLXNA2 and ID in humans.

Analysis of transcript-deleterious variants in Mendelian disorders: implications for RNA-based diagnostics.

BACKGROUND: At least 50% of patients with suspected Mendelian disorders remain undiagnosed after whole-exome sequencing (WES), and the extent to which non-coding variants that are not captured by WES contribute to this fraction is unclear. Whole transcriptome sequencing is a promising supplement to WES, although empirical data on the contribution of RNA analysis to the diagnosis of Mendelian diseases on a large scale are scarce. RESULTS: Here, we describe our experience with transcript-deleterious variants (TDVs) based on a cohort of 5647 families with suspected Mendelian diseases. We first interrogate all families for which the respective Mendelian phenotype could be mapped to a single locus to obtain an unbiased estimate of the contribution of TDVs at 18.9%. We examine the entire cohort and find that TDVs account for 15% of all "solved" cases. We compare the results of RT-PCR to in silico prediction. Definitive results from RT-PCR are obtained from blood-derived RNA for the overwhelming majority of variants (84.1%), and only a small minority (2.6%) fail analysis on all available RNA sources (blood-, skin fibroblast-, and urine renal epithelial cells-derived), which has important implications for the clinical application of RNA-seq. We also show that RNA analysis can establish the diagnosis in 13.5% of 155 patients who had received "negative" clinical WES reports. Finally, our data suggest a role for TDVs in modulating penetrance even in otherwise highly penetrant Mendelian disorders. CONCLUSIONS: Our results provide much needed empirical data for the impending implementation of diagnostic RNA-seq in conjunction with genome sequencing.

An exome-first approach to aid in the diagnosis of primary ciliary dyskinesia.

Unlike disorders of primary cilium, primary ciliary dyskinesia (PCD) has a much narrower clinical spectrum consistent with the limited tissue distribution of motile cilia. Nonetheless, PCD diagnosis can be challenging due to the overlapping features with other disorders and the requirement for sophisticated tests that are only available in specialized centers. We performed exome sequencing on all patients with a clinical suspicion of PCD but for whom no nasal nitric oxide test or ciliary functional assessment could be ordered. Among 81 patients (56 families), in whom PCD was suspected, 68% had pathogenic or likely pathogenic variants in established PCD-related genes that fully explain the phenotype (20 variants in 11 genes). The major clinical presentations were sinopulmonary infections (SPI) (n = 58), neonatal respiratory distress (NRD) (n = 2), laterality defect (LD) (n = 6), and combined LD/SPI (n = 15). Biallelic likely deleterious variants were also encountered in AKNA and GOLGA3, which we propose as novel candidates in a lung phenotype that overlaps clinically with PCD. We also encountered a PCD phenocopy caused by a pathogenic variant in ITCH, and a pathogenic variant in CEP164 causing Bardet-Biedl syndrome and PCD presentation as a very rare example of the dual presentation of these two disorders of the primary and motile cilia. Exome sequencing is a powerful tool that can help "democratize" the diagnosis of PCD, which is currently limited to highly specialized centers.

Multiple Family Members With Delayed Cord Separtion and Combined Immunodeficiency With Novel Mutation in IKBKB.

Background: Inhibitor of kappa kinase 2 (IKK2) deficiency is a recently described combined immunodeficiency. It undermines the nuclear factor-kappa B (NF-κB) activation pathway. Methods: The clinical and immunological data of four patients diagnosed with combined immunodeficiency (CID) from two related Saudi families were collected. Autozygosity mapping of all available members and whole exome sequencing of the index case were performed to define the genetic etiology. Results: The patients had early onset (2-4 months of age) severe infections caused by viruses, bacteria, mycobacteria, and fungi. They all had hypogammaglobulinemia and low absolute lymphocyte count. Their lymphocytes failed to respond to PHA mitogen stimulation. A novel homozygous non-sense mutation in the IKBKB gene, c.850C>T (p. Arg284*) was identified in the index patient and segregated with the disease in the rest of the family. He underwent hematopoietic stem cell transplantation (HSCT) from a fully matched sibling with no conditioning. The other three patients succumbed to their disease. Interestingly, all patients had delayed umbilical cord separation. Conclusion: IKK2 deficiency causes CID with high mortality. Immune reconstitution with HSCT should be considered as early as possible. Delayed umbilical cord separation in CID patients may be a clue to IKK2 deficiency.

The morbid genome of ciliopathies: an update.

PURPOSE: Ciliopathies are highly heterogeneous clinical disorders of the primary cilium. We aim to characterize a large cohort of ciliopathies phenotypically and molecularly. METHODS: Detailed phenotypic and genomic analysis of patients with ciliopathies, and functional characterization of novel candidate genes. RESULTS: In this study, we describe 125 families with ciliopathies and show that deleterious variants in previously reported genes, including cryptic splicing variants, account for 87% of cases. Additionally, we further support a number of previously reported candidate genes (BBIP1, MAPKBP1, PDE6D, and WDPCP), and propose nine novel candidate genes (CCDC67, CCDC96, CCDC172, CEP295, FAM166B, LRRC34, TMEM17, TTC6, and TTC23), three of which (LRRC34, TTC6, and TTC23) are supported by functional assays that we performed on available patient-derived fibroblasts. From a phenotypic perspective, we expand the phenomenon of allelism that characterizes ciliopathies by describing novel associations including WDR19-related Stargardt disease and SCLT1- and CEP164-related Bardet-Biedl syndrome. CONCLUSION: In this cohort of phenotypically and molecularly characterized ciliopathies, we draw important lessons that inform the clinical management and the diagnostics of this class of disorders as well as their basic biology.

Bi-allelic Variants in RALGAPA1 Cause Profound Neurodevelopmental Disability, Muscular Hypotonia, Infantile Spasms, and Feeding Abnormalities.

Ral (Ras-like) GTPases play an important role in the control of cell migration and have been implicated in Ras-mediated tumorigenicity. Recently, variants in RALA were also described as a cause of intellectual disability and developmental delay, indicating the relevance of this pathway to neuropediatric diseases. Here, we report the identification of bi-allelic variants in RALGAPA1 (encoding Ral GTPase activating protein catalytic alpha subunit 1) in four unrelated individuals with profound neurodevelopmental disability, muscular hypotonia, feeding abnormalities, recurrent fever episodes, and infantile spasms . Dysplasia of corpus callosum with focal thinning of the posterior part and characteristic facial features appeared to be unifying findings. RalGAPA1 was absent in the fibroblasts derived from two affected individuals suggesting a loss-of-function effect of the RALGAPA1 variants. Consequently, RalA activity was increased in these cell lines, which is in keeping with the idea that RalGAPA1 deficiency causes a constitutive activation of RalA. Additionally, levels of RalGAPB, a scaffolding subunit of the RalGAP complex, were dramatically reduced, indicating a dysfunctional RalGAP complex. Moreover, RalGAPA1 deficiency clearly increased cell-surface levels of lipid raft components in detached fibroblasts, which might indicate that anchorage-dependence of cell growth signaling is disturbed. Our findings indicate that the dysregulation of the RalA pathway has an important impact on neuronal function and brain development. In light of the partially overlapping phenotype between RALA- and RALGAPA1-associated diseases, it appears likely that dysregulation of the RalA signaling pathway leads to a distinct group of genetic syndromes that we suggest could be named RALopathies.

Confirming the recessive inheritance of PERP-related erythrokeratoderma.

Erythrokeratoderma (EK) is heterogeneous clinical entity characterized by excessive scaling with resulting erythrokeratotic plaques. Several genes have been linked to EK and they encode a number of proteins that are important for the integrity of the keratinocyte layer of the epidermis. PERP is a transcription factor that is activated by both p53 and p63. However, its deficiency in a mouse model appears to only recapitulate p63-mediated role in skin development and organization. We report an extended multiplex consanguineous family in which an EK phenotype with a striking similarity to that observed in Perp-/- mice, is mapped to an autozygous region on chromosome 6 that spans PERP. Whole-exome sequencing revealed a novel variant in PERP that fully segregated with the phenotype. Functional analysis of patient- and control-derived keratinocytes revealed a deleterious effect of the identified variant on the intracellular localization of PERP. A previous report showed that PERP mutation causes a dominant form of keratoderma but a single patient in that report with a homozygous variant in PERP suggests that recessive inheritance is also possible. Our results, therefore, support the establishment of an autosomal recessive PERP-related EK phenotype in humans.