Novel dominant and recessive variants in human ROBO1 cause distinct neurodevelopmental defects through different mechanisms.

The Roundabout (Robo) receptors, located on growth cones of neurons, induce axon repulsion in response to the extracellular ligand Slit. The Robo family of proteins controls midline crossing of commissural neurons during development in flies. Mono- and bi-allelic variants in human ROBO1 (HGNC: 10249) have been associated with incomplete penetrance and variable expressivity for a breath of phenotypes, including neurodevelopmental defects such as strabismus, pituitary defects, intellectual impairment, as well as defects in heart and kidney. Here, we report two novel ROBO1 variants associated with very distinct phenotypes. A homozygous missense p.S1522L variant in three affected siblings with nystagmus; and a monoallelic de novo p.D422G variant in a proband who presented with early-onset epileptic encephalopathy. We modeled these variants in Drosophila and first generated a null allele by inserting a CRIMIC T2A-GAL4 in an intron. Flies that lack robo1 exhibit reduced viability but have very severe midline crossing defects in the central nervous system. The fly wild-type cDNA driven by T2A-Gal4 partially rescues both defects. Overexpression of the human reference ROBO1 with T2A-GAL4 is toxic and reduces viability, whereas the recessive p.S1522L variant is less toxic, suggesting that it is a partial loss-of-function allele. In contrast, the dominant variant in fly robo1 (p.D413G) affects protein localization, impairs axonal guidance activity and induces mild phototransduction defects, suggesting that it is a neomorphic allele. In summary, our studies expand the phenotypic spectrum associated with ROBO1 variant alleles.

High prevalence of multilocus pathogenic variation in neurodevelopmental disorders in the Turkish population.

Neurodevelopmental disorders (NDDs) are clinically and genetically heterogenous; many such disorders are secondary to perturbation in brain development and/or function. The prevalence of NDDs is > 3%, resulting in significant sociocultural and economic challenges to society. With recent advances in family-based genomics, rare-variant analyses, and further exploration of the Clan Genomics hypothesis, there has been a logarithmic explosion in neurogenetic "disease-associated genes" molecular etiology and biology of NDDs; however, the majority of NDDs remain molecularly undiagnosed. We applied genome-wide screening technologies, including exome sequencing (ES) and whole-genome sequencing (WGS), to identify the molecular etiology of 234 newly enrolled subjects and 20 previously unsolved Turkish NDD families. In 176 of the 234 studied families (75.2%), a plausible and genetically parsimonious molecular etiology was identified. Out of 176 solved families, deleterious variants were identified in 218 distinct genes, further documenting the enormous genetic heterogeneity and diverse perturbations in human biology underlying NDDs. We propose 86 candidate disease-trait-associated genes for an NDD phenotype. Importantly, on the basis of objective and internally established variant prioritization criteria, we identified 51 families (51/176 = 28.9%) with multilocus pathogenic variation (MPV), mostly driven by runs of homozygosity (ROHs) - reflecting genomic segments/haplotypes that are identical-by-descent. Furthermore, with the use of additional bioinformatic tools and expansion of ES to additional family members, we established a molecular diagnosis in 5 out of 20 families (25%) who remained undiagnosed in our previously studied NDD cohort emanating from Turkey.

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.

A biallelic frameshift indel in PPP1R35 as a cause of primary microcephaly.

Protein phosphatase 1 regulatory subunit 35 (PPP1R35) encodes a centrosomal protein required for recruiting microtubule-binding elongation machinery. Several proteins in this centriole biogenesis pathway correspond to established primary microcephaly (MCPH) genes, and multiple model organism studies hypothesize PPP1R35 as a candidate MCPH gene. Here, using exome sequencing (ES) and family-based rare variant analyses, we report a homozygous, frameshifting indel deleting the canonical stop codon in the last exon of PPP1R35 [Chr7: c.753_*3delGGAAGCGTAGACCinsCG (p.Trp251Cysfs*22)]; the variant allele maps in a 3.7 Mb block of absence of heterozygosity (AOH) in a proband with severe MCPH (-4.3 SD at birth, -6.1 SD by 42 months), pachygyria, and global developmental delay from a consanguineous Turkish kindred. Droplet digital PCR (ddPCR) confirmed mutant mRNA expression in fibroblasts. In silico prediction of the translation of mutant PPP1R35 is expected to be elongated by 18 amino acids before encountering a downstream stop codon. This complex indel allele is absent in public databases (ClinVar, gnomAD, ARIC, 1000 genomes) and our in-house database of 14,000+ exomes including 1800+ Turkish exomes supporting predicted pathogenicity. Comprehensive literature searches for PPP1R35 variants yielded two probands affected with severe microcephaly (-15 SD and -12 SD) with the same homozygous indel from a single, consanguineous, Iranian family from a cohort of 404 predominantly Iranian families. The lack of heterozygous cases in two large cohorts representative of the genetic background of these two families decreased our suspicion of a founder allele and supports the contention of a recurrent mutation. We propose two potential secondary structure mutagenesis models for the origin of this variant allele mediated by hairpin formation between complementary GC rich segments flanking the stop codon via secondary structure mutagenesis.

Biallelic pathogenic variants in roundabout guidance receptor 1 associate with syndromic congenital anomalies of the kidney and urinary tract.

Congenital anomalies of the kidney and urinary tract (CAKUT) represent the most common cause of chronic kidney failure in children. Despite growing knowledge of the genetic causes of CAKUT, the majority of cases remain etiologically unsolved. Genetic alterations in roundabout guidance receptor 1 (ROBO1) have been associated with neuronal and cardiac developmental defects in living individuals. Although Slit-Robo signaling is pivotal for kidney development, diagnostic ROBO1 variants have not been reported in viable CAKUT to date. By next-generation-sequencing methods, we identified six unrelated individuals and two non-viable fetuses with biallelic truncating or combined missense and truncating variants in ROBO1. Kidney and genitourinary manifestation included unilateral or bilateral kidney agenesis, vesicoureteral junction obstruction, vesicoureteral reflux, posterior urethral valve, genital malformation, and increased kidney echogenicity. Further clinical characteristics were remarkably heterogeneous, including neurodevelopmental defects, intellectual impairment, cerebral malformations, eye anomalies, and cardiac defects. By in silico analysis, we determined the functional significance of identified missense variants and observed absence of kidney ROBO1 expression in both human and murine mutant tissues. While its expression in multiple tissues may explain heterogeneous organ involvement, variability of the kidney disease suggests gene dosage effects due to a combination of null alleles with mild hypomorphic alleles. Thus, comprehensive genetic analysis in CAKUT should include ROBO1 as a new cause of recessively inherited disease. Hence, in patients with already established ROBO1-associated cardiac or neuronal disorders, screening for kidney involvement is indicated.

Novel RETREG1 (FAM134B) founder allele is linked to HSAN2B and renal disease in a Turkish family.

Hereditary sensory and autonomic neuropathy type 2B (HSAN2B) is a rare autosomal recessive peripheral neuropathy caused by biallelic variants in RETREG1 (formerly FAM134B). HSAN2B is characterized by sensory impairment resulting in skin ulcerations, amputations, and osteomyelitis as well as variable weakness, spasticity, and autonomic dysfunction. Here, we report four affected individuals with recurrent osteomyelitis, ulceration, and amputation of hands and feet, sensory neuropathy, hyperhidrosis, urinary incontinence, and renal failure from a family without any known shared parental ancestry. Due to the history of chronic recurrent multifocal osteomyelitis and microcytic anemia, a diagnosis of Majeed syndrome was considered; however, sequencing of LPIN2 was negative. Family-based exome sequencing (ES) revealed a novel homozygous ultrarare RETREG1 variant NM_001034850.2:c.321G>A;p.Trp107Ter. Electrophysiological studies of the proband demonstrated axonal sensorimotor neuropathy predominantly in the lower extremities. Consistent with the lack of shared ancestry, the coefficient of inbreeding calculated from ES data was low (F = 0.002), but absence of heterozygosity (AOH) analysis demonstrated a 7.2 Mb AOH block surrounding the variant consistent with a founder allele. Two of the four affected individuals had unexplained renal failure which has not been reported in HSAN2B cases to date. Therefore, this report describes a novel RETREG1 founder allele and suggests renal failure may be an unrecognized feature of the RETREG1-disease spectrum.

Risk of sudden cardiac death in EXOSC5-related disease.

The RNA exosome is a multi-subunit complex involved in the processing, degradation, and regulated turnover of RNA. Several subunits are linked to Mendelian disorders, including pontocerebellar hypoplasia (EXOSC3, MIM #614678; EXOSC8, MIM #616081: and EXOSC9, MIM #618065) and short stature, hearing loss, retinitis pigmentosa, and distinctive facies (EXOSC2, MIM #617763). More recently, EXOSC5 (MIM *606492) was found to underlie an autosomal recessive neurodevelopmental disorder characterized by developmental delay, hypotonia, cerebellar abnormalities, and dysmorphic facies. An unusual feature of EXOSC5-related disease is the occurrence of complete heart block requiring a pacemaker in a subset of affected individuals. Here, we provide a detailed clinical and molecular characterization of two siblings with microcephaly, developmental delay, cerebellar volume loss, hypomyelination, with cardiac conduction and rhythm abnormalities including sinus node dysfunction, intraventricular conduction delay, atrioventricular block, and ventricular tachycardia (VT) due to compound heterozygous variants in EXOSC5: (1) NM_020158.4:c.341C > T (p.Thr114Ile; pathogenic, previously reported) and (2) NM_020158.4:c.302C > A (p.Thr101Lys; novel variant). A review of the literature revealed an additional family with biallelic EXOSC5 variants and cardiac conduction abnormalities. These clinical and molecular data provide compelling evidence that cardiac conduction abnormalities and arrhythmias are part of the EXOSC5-related disease spectrum and argue for proactive screening due to potential risk of sudden cardiac death.

Phenotype of CM-AVM2 caused by variants in EPHB4: how much overlap with hereditary hemorrhagic telangiectasia (HHT)?

PURPOSE: EPHB4 variants were recently reported to cause capillary malformation-arteriovenous malformation 2 (CM-AVM2). CM-AVM2 mimics RASA1-related CM-AVM1 and hereditary hemorrhagic telangiectasia (HHT), as clinical features include capillary malformations (CMs), telangiectasia, and arteriovenous malformations (AVMs). Epistaxis, another clinical feature that overlaps with HHT, was reported in several cases. Based on the clinical overlap of CM-AVM2 and HHT, we hypothesized that patients considered clinically suspicious for HHT with no variant detected in an HHT gene (ENG, ACVRL1, or SMAD4) may have an EPHB4 variant. METHODS: Exome sequencing or a next-generation sequencing panel including EPHB4 was performed on individuals with previously negative molecular genetic testing for the HHT genes and/or RASA1. RESULTS: An EPHB4 variant was identified in ten unrelated cases. Seven cases had a pathogenic EPHB4 variant, including one with mosaicism. Three cases had an EPHB4 variant of uncertain significance. The majority had epistaxis (6/10 cases) and telangiectasia (8/10 cases), as well as CMs. Two of ten cases had a central nervous system AVM. CONCLUSIONS: Our results emphasize the importance of considering CM-AVM2 as part of the clinical differential for HHT and other vascular malformation syndromes. Yet, these cases highlight significant differences in the cutaneous presentations of CM-AVM2 versus HHT.

Genome sequencing reveals a deep intronic splicing ACVRL1 mutation hotspot in Hereditary Haemorrhagic Telangiectasia.

INTRODUCTION: Hereditary haemorrhagic telangiectasia (HHT) is a genetically heterogeneous disorder caused by mutations in the genes ENG, ACVRL1, and SMAD4. Yet the genetic cause remains unknown for some families even after exhaustive exome analysis. We hypothesised that non-coding regions of the known HHT genes may harbour variants that disrupt splicing in these cases. METHODS: DNA from 35 individuals with clinical findings of HHT and 2 healthy controls from 13 families underwent whole genome sequencing. Additionally, 87 unrelated cases suspected to have HHT were evaluated using a custom designed next-generation sequencing panel to capture the coding and non-coding regions of ENG, ACVRL1 and SMAD4. Individuals from both groups had tested negative previously for a mutation in the coding region of known HHT genes. Samples were sequenced on a HiSeq2500 instrument and data were analysed to identify novel and rare variants. RESULTS: Eight cases had a novel non-coding ACVRL1 variant that disrupted splicing. One family had an ACVRL1intron 9:chromosome 3 translocation, the first reported case of a translocation causing HHT. The other seven cases had a variant located within a ~300 bp CT-rich 'hotspot' region of ACVRL1intron 9 that disrupted splicing. CONCLUSIONS: Despite the difficulty of interpreting deep intronic variants, our study highlights the importance of non-coding regions in the disease mechanism of HHT, particularly the CT-rich hotspot region of ACVRL1intron 9. The addition of this region to HHT molecular diagnostic testing algorithms will improve clinical sensitivity.

Developmental genomics of limb malformations: Allelic series in association with gene dosage effects contribute to the clinical variability.

Genetic heterogeneity, reduced penetrance, and variable expressivity, the latter including asymmetric body axis plane presentations, have all been described in families with congenital limb malformations (CLMs). Interfamilial and intrafamilial heterogeneity highlight the complexity of the underlying genetic pathogenesis of these developmental anomalies. Family-based genomics by exome sequencing (ES) and rare variant analyses combined with whole-genome array-based comparative genomic hybridization were implemented to investigate 18 families with limb birth defects. Eleven of 18 (61%) families revealed explanatory variants, including 7 single-nucleotide variant alleles and 3 copy number variants (CNVs), at previously reported "disease trait associated loci": BHLHA9, GLI3, HOXD cluster, HOXD13, NPR2, and WNT10B. Breakpoint junction analyses for all three CNV alleles revealed mutational signatures consistent with microhomology-mediated break-induced replication, a mechanism facilitated by Alu/Alu-mediated rearrangement. Homozygous duplication of BHLHA9 was observed in one Turkish kindred and represents a novel contributory genetic mechanism to Gollop-Wolfgang Complex (MIM: 228250), where triplication of the locus has been reported in one family from Japan (i.e., 4n = 2n + 2n versus 4n = 3n + 1n allelic configurations). Genes acting on limb patterning are sensitive to a gene dosage effect and are often associated with an allelic series. We extend an allele-specific gene dosage model to potentially assist, in an adjuvant way, interpretations of interconnections among an allelic series, clinical severity, and reduced penetrance of the BHLHA9-related CLM spectrum.