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.

The utility of gene sequencing in identifying an underlying genetic disorder in prenatally suspected lower urinary tract obstruction.

OBJECTIVE: Fetal megacystis generally presents as suspected lower urinary tract obstruction (LUTO), which is associated with severe perinatal morbidity. Genetic etiologies underlying LUTO or a LUTO-like initial presentation are poorly understood. Our objectives are to describe single gene etiologies in fetuses initially ascertained to have suspected LUTO and to elucidate genotype-phenotype correlations. METHODS: A retrospective case series of suspected fetal LUTO positive for a molecular diagnosis was collected from five centers in the Fetal Sequencing Consortium. Demographics, sonograms, genetic testing including variant classification, and delivery outcomes were abstracted. RESULTS: Seven cases of initially prenatally suspected LUTO-positive for a molecular diagnosis were identified. In no case was the final diagnosis established as urethral obstruction that is, LUTO. All variants were classified as likely pathogenic or pathogenic. Smooth muscle deficiencies involving the bladder wall and interfering with bladder emptying were identified in five cases: MYOCD (2), ACTG2 (2), and MYH11 (1). Other genitourinary and/or non-genitourinary malformations were seen in two cases involving KMT2D (1) and BBS10 (1). CONCLUSION: Our series illustrates the value of molecular diagnostics in the workup of fetuses who present with prenatally suspected LUTO but who may have a non-LUTO explanation for their prenatal ultrasound findings.

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.

Brain monoamine vesicular transport disease caused by homozygous SLC18A2 variants: A study in 42 affected individuals.

PURPOSE: Brain monoamine vesicular transport disease is an infantile-onset movement disorder that mimics cerebral palsy. In 2013, the homozygous SLC18A2 variant, p.Pro387Leu, was first reported as a cause of this rare disorder, and dopamine agonists were efficient for treating affected individuals from a single large family. To date, only 6 variants have been reported. In this study, we evaluated genotype-phenotype correlations in individuals with biallelic SLC18A2 variants. METHODS: A total of 42 affected individuals with homozygous SLC18A2 variant alleles were identified. We evaluated genotype-phenotype correlations and the missense variants in the affected individuals based on the structural modeling of rat VMAT2 encoded by Slc18a2, with cytoplasm- and lumen-facing conformations. A Caenorhabditis elegans model was created for functional studies. RESULTS: A total of 19 homozygous SLC18A2 variants, including 3 recurrent variants, were identified using exome sequencing. The affected individuals typically showed global developmental delay, hypotonia, dystonia, oculogyric crisis, and autonomic nervous system involvement (temperature dysregulation/sweating, hypersalivation, and gastrointestinal dysmotility). Among the 58 affected individuals described to date, 16 (28%) died before the age of 13 years. Of the 17 patients with p.Pro237His, 9 died, whereas all 14 patients with p.Pro387Leu survived. Although a dopamine agonist mildly improved the disease symptoms in 18 of 21 patients (86%), some affected individuals with p.Ile43Phe and p.Pro387Leu showed milder phenotypes and presented prolonged survival even without treatment. The C. elegans model showed behavioral abnormalities. CONCLUSION: These data expand the phenotypic and genotypic spectra of SLC18A2-related disorders.

Biallelic Loss-of-Function NDUFA12 Variants Cause a Wide Phenotypic Spectrum from Leigh/Leigh-Like Syndrome to Isolated Optic Atrophy.

BACKGROUND: Biallelic loss-of-function NDUFA12 variants have hitherto been linked to mitochondrial complex I deficiency presenting with heterogeneous clinical and radiological features in nine cases only. OBJECTIVES: To fully characterize, both phenotypically and genotypically, NDUFA12-related mitochondrial disease. METHODS: We collected data from cases identified by screening genetic databases of several laboratories worldwide and systematically reviewed the literature. RESULTS: Nine unreported NDUFA12 cases from six pedigrees were identified, with presentation ranging from movement disorder phenotypes (dystonia and/or spasticity) to isolated optic atrophy. MRI showed basal ganglia abnormalities (n = 6), optic atrophy (n = 2), or was unremarkable (n = 1). All carried homozygous truncating NDUFA12 variants, three of which are novel. CONCLUSIONS: Our case series expands phenotype-genotype correlations in NDUFA12-associated mitochondrial disease, providing evidence of intra- and inter-familial clinical heterogeneity for the same variant. It confirms NDUFA12 variants should be included in the diagnostic workup of Leigh/Leigh-like syndromes - particularly with dystonia - as well as isolated optic atrophy.

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.

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.

Molecular autopsy by proxy in preconception counseling.

Monogenic diseases that result in early pregnancy loss or neonatal death are genetically and phenotypically highly variable. This often poses significant challenges in arriving at a molecular diagnosis for reproductive planning. Molecular autopsy by proxy (MABP) refers to the genetic testing of relatives of deceased individuals to deduce the cause of death. Here, we specifically tested couples who lost one or more children/pregnancies with no available DNA. We developed our testing strategy using whole exome sequencing data from 83 consanguineous Saudi couples. We detected the shared carrier state of 50 pathogenic variants/likely pathogenic variants in 43 families and of 28 variants of uncertain significance in 24 families. Negative results were seen in 16 couples after variant reclassification. In 10 families, the risk of more than one genetic disease was documented. Secondary findings were seen in 10 families: either genetic variants with potential clinical consequences for the tested individual or a female carrier for X-linked conditions. This couple-based approach has enabled molecularly informed genetic counseling for 52% (43/83 families). Given the predominance of autosomal recessive causes of pregnancy and child death in consanguineous populations, MABP can be a helpful approach to consanguineous couples who seek counseling but lack molecular data on their deceased offspring.

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.

Generation of Monogenic Candidate Genes for Human Nephrotic Syndrome Using 3 Independent Approaches.

INTRODUCTION: Steroid-resistant nephrotic syndrome (SRNS) is the second most common cause of chronic kidney disease during childhood. Identification of 63 monogenic human genes has delineated 12 distinct pathogenic pathways. METHODS: Here, we generated 2 independent sets of nephrotic syndrome (NS) candidate genes to augment the discovery of additional monogenic causes based on whole-exome sequencing (WES) data from 1382 families with NS. RESULTS: We first identified 63 known monogenic causes of NS in mice from public databases and scientific publications, and 12 of these genes overlapped with the 63 known human monogenic SRNS genes. Second, we used a set of 64 genes that are regulated by the transcription factor Wilms tumor 1 (WT1), which causes SRNS if mutated. Thirteen of these WT1-regulated genes overlapped with human or murine NS genes. Finally, we overlapped these lists of murine and WT1 candidate genes with our list of 120 candidate genes generated from WES in 1382 NS families, to identify novel candidate genes for monogenic human SRNS. Using this approach, we identified 7 overlapping genes, of which 3 genes were shared by all datasets, including SYNPO. We show that loss-of-function of SYNPO leads to decreased CDC42 activity and reduced podocyte migration rate, both of which are rescued by overexpression of wild-type complementary DNA (cDNA), but not by cDNA representing the patient mutation. CONCLUSION: Thus, we identified 3 novel candidate genes for human SRNS using 3 independent, nonoverlapping hypotheses, and generated functional evidence for SYNPO as a novel potential monogenic cause of NS.