Multilocus pathogenic variants contribute to intrafamilial clinical heterogeneity: a retrospective study of sibling pairs with neurodevelopmental disorders

BACKGROUND

Multilocus pathogenic variants (MPVs) are genetic changes that affect multiple gene loci or regions of the genome, collectively leading to multiple molecular diagnoses. MPVs may also contribute to intrafamilial phenotypic variability between affected individuals within a nuclear family. In this study, we aim to gain further insights into the influence of MPVs on a disease manifestation in individual research subjects and explore the complexities of the human genome within a familial context.

METHODS

We conducted a systematic reanalysis of exome sequencing data and runs of homozygosity (ROH) regions of 47 sibling pairs previously diagnosed with various neurodevelopmental disorders (NDD).

RESULTS

We found siblings with MPVs driven by long ROH regions in 8.5% of families (4/47). The patients with MPVs exhibited significantly higher FROH values (p-value = 1.4e-2) and larger total ROH length (p-value = 1.8e-2). Long ROH regions mainly contribute to this pattern; the siblings with MPVs have a larger total size of long ROH regions than their siblings in all families (p-value = 6.9e-3). Whereas the short ROH regions in the siblings with MPVs are lower in total size compared to their sibling pairs with single locus pathogenic variants (p-value = 0.029), and there are no statistically significant differences in medium ROH regions between sibling pairs (p-value = 0.52).

CONCLUSION

This study sheds light on the significance of considering MPVs in families with affected sibling pairs and the role of ROH as an adjuvant tool in explaining clinical variability within families. Identifying individuals carrying MPVs may have implications for disease management, identification of possible disease risks to different family members, genetic counseling and exploring personalized treatment approaches.

Bi-allelic ACBD6 variants lead to a neurodevelopmental syndrome with progressive and complex movement disorders

The acyl-CoA-binding domain-containing protein 6 (ACBD6) is ubiquitously expressed, plays a role in the acylation of lipids and proteins and regulates the N-myristoylation of proteins via N-myristoyltransferase enzymes (NMTs). However, its precise function in cells is still unclear, as is the consequence of ACBD6 defects on human pathophysiology. Using exome sequencing and extensive international data sharing efforts, we identified 45 affected individuals from 28 unrelated families (consanguinity 93%) with bi-allelic pathogenic, predominantly loss-of-function (18/20) variants in ACBD6. We generated zebrafish and Xenopus tropicalis acbd6 knockouts by CRISPR/Cas9 and characterized the role of ACBD6 on protein N-myristoylation with myristic acid alkyne (YnMyr) chemical proteomics in the model organisms and human cells, with the latter also being subjected further to ACBD6 peroxisomal localization studies. The affected individuals (23 males and 22 females), aged 1-50 years, typically present with a complex and progressive disease involving moderate-to-severe global developmental delay/intellectual disability (100%) with significant expressive language impairment (98%), movement disorders (97%), facial dysmorphism (95%) and mild cerebellar ataxia (85%) associated with gait impairment (94%), limb spasticity/hypertonia (76%), oculomotor (71%) and behavioural abnormalities (65%), overweight (59%), microcephaly (39%) and epilepsy (33%). The most conspicuous and common movement disorder was dystonia (94%), frequently leading to early-onset progressive postural deformities (97%), limb dystonia (55%) and cervical dystonia (31%). A jerky tremor in the upper limbs (63%), a mild head tremor (59%), parkinsonism/hypokinesia developing with advancing age (32%) and simple motor and vocal tics were among other frequent movement disorders. Midline brain malformations including corpus callosum abnormalities (70%), hypoplasia/agenesis of the anterior commissure (66%), short midbrain and small inferior cerebellar vermis (38% each) as well as hypertrophy of the clava (24%) were common neuroimaging findings. Acbd6-deficient zebrafish and Xenopus models effectively recapitulated many clinical phenotypes reported in patients including movement disorders, progressive neuromotor impairment, seizures, microcephaly, craniofacial dysmorphism and midbrain defects accompanied by developmental delay with increased mortality over time. Unlike ACBD5, ACBD6 did not show a peroxisomal localization and ACBD6-deficiency was not associated with altered peroxisomal parameters in patient fibroblasts. Significant differences in YnMyr-labelling were observed for 68 co- and 18 post-translationally N-myristoylated proteins in patient-derived fibroblasts. N-myristoylation was similarly affected in acbd6-deficient zebrafish and X. tropicalis models, including Fus, Marcks and Chchd-related proteins implicated in neurological diseases. The present study provides evidence that bi-allelic pathogenic variants in ACBD6 lead to a distinct neurodevelopmental syndrome accompanied by complex and progressive cognitive and movement disorders.

Sensory experiences questionnaire unravels differences in sensory profiles between MECP2-related disorders

The methyl CpG-binding protein-2 (MECP2) gene is located on the Xq28 region. Loss of function mutations or increased copies of MECP2 result in Rett syndrome (RTT) and MECP2 duplication syndrome (MDS), respectively. Individuals with both disorders exhibit overlapping autism symptoms, yet few studies have dissected the differences between these gene dosage sensitive disorders. Further, research examining sensory processing patterns in persons with RTT and MDS is largely absent. Thus, the goal of this study was to analyze and compare sensory processing patterns in persons with RTT and MDS. Towards this goal, caregivers of 50 female individuals with RTT and 122 male individuals with MDS, between 1 and 46 years of age, completed a standardized measure of sensory processing, the Sensory Experiences Questionnaire. Patterns detected in both disorders were compared against each other and against normative values. We found sensory processing abnormalities for both hyper- and hypo-sensitivity in both groups. Interestingly, abnormalities in MDS were more pronounced compared with in RTT, particularly with items concerning hypersensitivity and sensory seeking, but not hyposensitivity. Individuals with MDS also exhibited greater sensory symptoms compared with RTT in the areas of tactile and vestibular sensory processing and for both social and nonsocial stimuli. This study provides a first description of sensory symptoms in individuals with RTT and individuals with MDS. Similar to other neurodevelopmental disorders, a variety of sensory processing abnormalities was found. These findings reveal a first insight into sensory processing abnormalities caused by a dosage sensitive gene and may ultimately help guide therapeutic approaches for these disorders.

A Comprehensive Examination of Clinical Characteristics and Determinants of Long-Term Outcomes in Pediatric Cerebral Sinus Venous Thrombosis

BACKGROUND

To analyze the clinical and neuroimaging features, risk factors, treatment choices, and long-term clinical outcomes in children with cerebral sinus venous thrombosis (CSVT).

METHODS

This is a retrospective cohort study of children diagnosed with CSVT between 2002 and 2018 at Texas Children's Hospital.

RESULTS

A total of 183 children (male: 62.3%) with CSVT were included. The average presenting age was 7.7 years (S.D.: 5.6). The mean follow-up duration was 33.7 months (S.D.: 38.6). The most common presenting clinical feature was headache (36.6%). Head and neck infections other than meningitis (36.6%) were the most common risk factors. Prevalent neurological examination findings included motor deficit (21.3%) and altered mental status (AMS, 20.2%). Neuroimaging features included hemorrhagic infarction (19.6%), ischemic infarction (8.2%), and intracranial hemorrhage without infarction (5.5%). The most common site of thrombosis was the superior sagittal sinus (37.2%), with 78.2% of patients demonstrating involvement of multiple sinuses. Treatment of choice was low-molecular-weight heparin in 69.4% of patients. Factors associated with worse clinical outcomes included head and neck infections, malignancy (other than hematologic), cardiac disease, and recent surgery; seizure and dehydration on initial presentation; motor abnormalities and AMS on initial examination; ischemic infarct only; and involvement of vein of Trolard on neuroimaging. Thrombus condition on repeat imaging, receiving any anticoagulant/antithrombotic treatment, treatment duration, or follow-up duration was not associated with severity of long-term outcome.

CONCLUSIONS

CSVT may lead to unfavorable long-term outcomes in a remarkable portion of pediatric patients. Thus, a high index of suspicion and early and appropriate management of pediatric CSVT is imperative.

Homozygous missense variants in YKT6 result in loss of function and are associated with developmental delay, with or without severe infantile liver disease and risk for hepatocellular carcinoma

PURPOSE

YKT6 plays important roles in multiple intracellular vesicle trafficking events but has not been associated with Mendelian diseases.

METHODS

We report three unrelated individuals with rare homozygous missense variants in YKT6 who exhibited neurological disease with or without a progressive infantile liver disease. We modeled the variants in Drosophila. We generated wild-type and variant genomic rescue constructs (GRs) of the fly ortholog dYkt6 and compared their ability in rescuing the loss-of-function phenotypes in mutant flies. We also generated a dYkt6KozakGAL4 allele to assess the expression pattern of dYkt6.

RESULTS

Two individuals are homozygous for YKT6 [NM_006555.3:c.554A>G p.(Tyr185Cys)] and exhibited normal prenatal course followed by failure to thrive, developmental delay and progressive liver disease. Haplotype analysis identified a shared homozygous region flanking the variant, suggesting a common ancestry. The third individual is homozygous for YKT6 [NM_006555.3:c.191A>G p.(Tyr64Cys)] and exhibited neurodevelopmental disorders and optic atrophy. Fly dYkt6 is essential and is expressed in the fat body (analogous to liver) and central nervous system. Wild-type GR can rescue the lethality and autophagic flux defects whereas the variants are less efficient in rescuing the phenotypes.

CONCLUSION

The YKT6 variants are partial loss-of-function alleles and the p.(Tyr185Cys) is more severe than p.(Tyr64Cys).

HMZDupFinder: a robust computational approach for detecting intragenic homozygous duplications from exome sequencing data

Homozygous duplications contribute to genetic disease by altering gene dosage or disrupting gene regulation and can be more deleterious to organismal biology than heterozygous duplications. Intragenic exonic duplications can result in loss-of-function (LoF) or gain-of-function (GoF) alleles that when homozygosed, i.e. brought to homozygous state at a locus by identity by descent or state, could potentially result in autosomal recessive (AR) rare disease traits. However, the detection and functional interpretation of homozygous duplications from exome sequencing data remains a challenge. We developed a framework algorithm, HMZDupFinder, that is designed to detect exonic homozygous duplications from exome sequencing (ES) data. The HMZDupFinder algorithm can efficiently process large datasets and accurately identifies small intragenic duplications, including those associated with rare disease traits. HMZDupFinder called 965 homozygous duplications with three or less exons from 8,707 ES with a recall rate of 70.9% and a precision of 16.1%. We experimentally confirmed 8/10 rare homozygous duplications. Pathogenicity assessment of these copy number variant alleles allowed clinical genomics contextualization for three homozygous duplications alleles, including two affecting known OMIM disease genes EDAR (MIM# 224900), TNNT1(MIM# 605355), and one variant in a novel candidate disease gene: PAAF1.

Biallelic variation in the choline and ethanolamine transporter FLVCR1 underlies a pleiotropic disease spectrum from adult neurodegeneration to severe developmental disorders

FLVCR1 encodes Feline leukemia virus subgroup C receptor 1 (FLVCR1), a solute carrier (SLC) transporter within the Major Facilitator Superfamily. FLVCR1 is a widely expressed transmembrane protein with plasma membrane and mitochondrial isoforms implicated in heme, choline, and ethanolamine transport. While Flvcr1 knockout mice die in utero with skeletal malformations and defective erythropoiesis reminiscent of Diamond-Blackfan anemia, rare biallelic pathogenic FLVCR1 variants are linked to childhood or adult-onset neurodegeneration of the retina, spinal cord, and peripheral nervous system. We ascertained from research and clinical exome sequencing 27 individuals from 20 unrelated families with biallelic ultra-rare missense and predicted loss-of-function (pLoF) FLVCR1 variant alleles. We characterize an expansive FLVCR1 phenotypic spectrum ranging from adult-onset retinitis pigmentosa to severe developmental disorders with microcephaly, reduced brain volume, epilepsy, spasticity, and premature death. The most severely affected individuals, including three individuals with homozygous pLoF variants, share traits with Flvcr1 knockout mice and Diamond-Blackfan anemia including macrocytic anemia and congenital skeletal malformations. Pathogenic FLVCR1 missense variants primarily lie within transmembrane domains and reduce choline and ethanolamine transport activity compared with wild-type FLVCR1 with minimal impact on FLVCR1 stability or subcellular localization. Several variants disrupt splicing in a mini-gene assay which may contribute to genotype-phenotype correlations. Taken together, these data support an allele-specific gene dosage model in which phenotypic severity reflects residual FLVCR1 activity. This study expands our understanding of Mendelian disorders of choline and ethanolamine transport and demonstrates the importance of choline and ethanolamine in neurodevelopment and neuronal homeostasis.

MECP2-related disorders while gene-based therapies are on the horizon

The emergence of new genetic tools has led to the discovery of the genetic bases of many intellectual and developmental disabilities. This creates exciting opportunities for research and treatment development, and a few genetic disorders (e.g., spinal muscular atrophy) have recently been treated with gene-based therapies. MECP2 is found on the X chromosome and regulates the transcription of thousands of genes. Loss of MECP2 gene product leads to Rett Syndrome, a disease found primarily in females, and is characterized by developmental regression, motor dysfunction, midline hand stereotypies, autonomic nervous system dysfunction, epilepsy, scoliosis, and autistic-like behavior. Duplication of MECP2 causes MECP2 Duplication Syndrome (MDS). MDS is found mostly in males and presents with developmental delay, hypotonia, autistic features, refractory epilepsy, and recurrent respiratory infections. While these two disorders share several characteristics, their differences (e.g., affected sex, age of onset, genotype/phenotype correlations) are important to distinguish in the light of gene-based therapy because they require opposite solutions. This review explores the clinical features of both disorders and highlights these important clinical differences.

Development and validation of parent-reported gastrointestinal health scale in MECP2 duplication syndrome

BACKGROUND/AIMS

We aimed to develop a validated patient-reported Gastrointestinal Health Scale (GHS) specific to MECP2 Duplication Syndrome (MDS) to be used in clinical trials.

METHODS

MDS parents completed a Gastrointestinal Health Questionnaire (GHQ) to investigate the most relevant and important items associated with gastrointestinal problems in MECP2-related disorders. Item reduction was executed according to EORTC guidelines. We performed reliability and validity studies for the finalized scale.

RESULTS

A total of 106 surveys were eligible for item reduction and validation processes. The initial 55 items were reduced to 38 items based on parent responses, expert opinion, and initial confirmatory factor analysis (CFA). The final MDS-specific GHS included 38 items and 7 factors that underwent further reliability and validity assessments. The power of the study was at least 0.982. The Cronbach's alphas of the instruments were General Health: 0.799, Eating-Chewing-Swallowing: 0.809, Reflux: 0.794, Motility: 0.762, Mood: 0.906, Medication: 0.595, Parenting: 0.942 and all items together: 0.928. The correlation coefficient between total and individual item scores ranged from 0.215 to 0.730. Because of the ordinal nature of the variables, the diagonal weighted least squares estimation (DWLS) method was used to execute the CFA and Structural Equation Modeling. The GHS had excellent model fit with the acceptable range of fit indices values.

CONCLUSIONS

We developed a parent-reported, reliable, and valid MDS-specific GHS. This scale can be utilized in clinical settings or as an outcome measure in translational and clinical research.

Generation of five induced pluripotent stem cell lines from patients with MECP2 Duplication Syndrome

MECP2 Duplication Syndrome (MDS) is a rare, severe neurodevelopmental disorder arising from duplications in the Xq28 region containing the MECP2 gene that predominantly affects males. We generated five human induced pluripotent stem cell (iPSC) lines from the fibroblasts of individuals carrying between 0.355 and 11.2 Mb size duplications in the chromosomal locus containing MECP2. All lines underwent extensive testing to confirm MECP2 duplication and iPSC-related features such as morphology, pluripotency markers, and trilineage differentiation potential. These lines are a valuable resource for molecular and functional studies of MDS as well as screening for a variety of therapeutic approaches.

RNA methyltransferase SPOUT1/CENP-32 links mitotic spindle organization with the neurodevelopmental disorder SpADMiSS

SPOUT1/CENP-32 encodes a putative SPOUT RNA methyltransferase previously identified as a mitotic chromosome associated protein. SPOUT1/CENP-32 depletion leads to centrosome detachment from the spindle poles and chromosome misalignment. Aided by gene matching platforms, we identified 24 individuals with neurodevelopmental delays from 18 families with bi-allelic variants in SPOUT1/CENP-32 detected by exome/genome sequencing. Zebrafish spout1/cenp-32 mutants showed reduction in larval head size with concomitant apoptosis likely associated with altered cell cycle progression. In vivo complementation assays in zebrafish indicated that SPOUT1/CENP-32 missense variants identified in humans are pathogenic. Crystal structure analysis of SPOUT1/CENP-32 revealed that most disease-associated missense variants mapped to the catalytic domain. Additionally, SPOUT1/CENP-32 recurrent missense variants had reduced methyltransferase activity in vitro and compromised centrosome tethering to the spindle poles in human cells. Thus, SPOUT1/CENP-32 pathogenic variants cause an autosomal recessive neurodevelopmental disorder: SpADMiSS ( SPOUT1 Associated Development delay Microcephaly Seizures Short stature) underpinned by mitotic spindle organization defects and consequent chromosome segregation errors.

Detection of mosaic and population-level structural variants with Sniffles2

Calling structural variations (SVs) is technically challenging, but using long reads remains the most accurate way to identify complex genomic alterations. Here we present Sniffles2, which improves over current methods by implementing a repeat aware clustering coupled with a fast consensus sequence and coverage-adaptive filtering. Sniffles2 is 11.8 times faster and 29% more accurate than state-of-the-art SV callers across different coverages (5-50×), sequencing technologies (ONT and HiFi) and SV types. Furthermore, Sniffles2 solves the problem of family-level to population-level SV calling to produce fully genotyped VCF files. Across 11 probands, we accurately identified causative SVs around MECP2, including highly complex alleles with three overlapping SVs. Sniffles2 also enables the detection of mosaic SVs in bulk long-read data. As a result, we identified multiple mosaic SVs in brain tissue from a patient with multiple system atrophy. The identified SV showed a remarkable diversity within the cingulate cortex, impacting both genes involved in neuron function and repetitive elements.

Biallelic missense variants in COG3 cause a congenital disorder of glycosylation with impairment of retrograde vesicular trafficking

Biallelic variants in genes for seven out of eight subunits of the conserved oligomeric Golgi complex (COG) are known to cause recessive congenital disorders of glycosylation (CDG) with variable clinical manifestations. COG3 encodes a constituent subunit of the COG complex that has not been associated with disease traits in humans. Herein, we report two COG3 homozygous missense variants in four individuals from two unrelated consanguineous families that co-segregated with COG3-CDG presentations. Clinical phenotypes of affected individuals include global developmental delay, severe intellectual disability, microcephaly, epilepsy, facial dysmorphism, and variable neurological findings. Biochemical analysis of serum transferrin from one family showed the loss of a single sialic acid. Western blotting on patient-derived fibroblasts revealed reduced COG3 and COG4. Further experiments showed delayed retrograde vesicular recycling in patient cells. This report adds to the knowledge of the COG-CDG network by providing collective evidence for a COG3-CDG rare disease trait and implicating a likely pathology of the disorder as the perturbation of Golgi trafficking.

Break-induced replication underlies formation of inverted triplications and generates unexpected diversity in haplotype structures

Background

The duplication-triplication/inverted-duplication (DUP-TRP/INV-DUP) structure is a type of complex genomic rearrangement (CGR) hypothesized to result from replicative repair of DNA due to replication fork collapse. It is often mediated by a pair of inverted low-copy repeats (LCR) followed by iterative template switches resulting in at least two breakpoint junctions in cis . Although it has been identified as an important mutation signature of pathogenicity for genomic disorders and cancer genomes, its architecture remains unresolved and is predicted to display at least four structural variation (SV) haplotypes.

Results

Here we studied the genomic architecture of DUP-TRP/INV-DUP by investigating the genomic DNA of 24 patients with neurodevelopmental disorders identified by array comparative genomic hybridization (aCGH) on whom we found evidence for the existence of 4 out of 4 predicted SV haplotypes. Using a combination of short-read genome sequencing (GS), long- read GS, optical genome mapping and StrandSeq the haplotype structure was resolved in 18 samples. This approach refined the point of template switching between inverted LCRs in 4 samples revealing a DNA segment of ∼2.2-5.5 kb of 100% nucleotide similarity. A prediction model was developed to infer the LCR used to mediate the non-allelic homology repair.

Conclusions

These data provide experimental evidence supporting the hypothesis that inverted LCRs act as a recombinant substrate in replication-based repair mechanisms. Such inverted repeats are particularly relevant for formation of copy-number associated inversions, including the DUP-TRP/INV-DUP structures. Moreover, this type of CGR can result in multiple conformers which contributes to generate diverse SV haplotypes in susceptible loci .

Novel LSS variants in alopecia and intellectual disability syndrome: New case report and clinical spectrum of LSS-related rare disease traits

Pathogenic biallelic variants in LSS are associated with three Mendelian rare disease traits including congenital cataract type 44, autosomal recessive hypotrichosis type 14, and alopecia-intellectual disability syndrome type 4 (APMR4). We performed trio research exome sequencing on a family with a four-year-old male with global developmental delay, epilepsy and striking alopecia, and identified novel compound heterozygous LSS splice site (c.14+2T>C) and missense (c.1357 G>A; p.V453L) variant alleles. Rare features associated with APMR4 such as cryptorchidism, micropenis, mild cortical brain atrophy and thin corpus callosum were detected. Previously unreported APMR4 findings including cerebellar involvement in the form of unsteady ataxic gait, small vermis with prominent folia, were noted. A review of all reported variants to date in 29 families with LSS-related phenotypes showed an emerging genotype-phenotype correlation. Our report potentially expands LSS-related phenotypic spectrum and highlights the importance of performing brain imaging in LSS-related conditions.

Monoallelic variation in DHX9, the gene encoding the DExH-box helicase DHX9, underlies neurodevelopment disorders and Charcot-Marie-Tooth disease

DExD/H-box RNA helicases (DDX/DHX) are encoded by a large paralogous gene family; in a subset of these human helicase genes, pathogenic variation causes neurodevelopmental disorder (NDD) traits and cancer. DHX9 encodes a BRCA1-interacting nuclear helicase regulating transcription, R-loops, and homologous recombination and exhibits the highest mutational constraint of all DDX/DHX paralogs but remains unassociated with disease traits in OMIM. Using exome sequencing and family-based rare-variant analyses, we identified 20 individuals with de novo, ultra-rare, heterozygous missense or loss-of-function (LoF) DHX9 variant alleles. Phenotypes ranged from NDDs to the distal symmetric polyneuropathy axonal Charcot-Marie-Tooth disease (CMT2). Quantitative Human Phenotype Ontology (HPO) analysis demonstrated genotype-phenotype correlations with LoF variants causing mild NDD phenotypes and nuclear localization signal (NLS) missense variants causing severe NDD. We investigated DHX9 variant-associated cellular phenotypes in human cell lines. Whereas wild-type DHX9 was restricted to the nucleus, NLS missense variants abnormally accumulated in the cytoplasm. Fibroblasts from an individual with an NLS variant also showed abnormal cytoplasmic DHX9 accumulation. CMT2-associated missense variants caused aberrant nucleolar DHX9 accumulation, a phenomenon previously associated with cellular stress. Two NDD-associated variants, p.Gly411Glu and p.Arg761Gln, altered DHX9 ATPase activity. The severe NDD-associated variant p.Arg141Gln did not affect DHX9 localization but instead increased R-loop levels and double-stranded DNA breaks. Dhx9-/- mice exhibited hypoactivity in novel environments, tremor, and sensorineural hearing loss. All together, these results establish DHX9 as a critical regulator of mammalian neurodevelopment and neuronal homeostasis.

The clinical and genetic spectrum of autosomal-recessive TOR1A-related disorders

In the field of rare diseases, progress in molecular diagnostics led to the recognition that variants linked to autosomal-dominant neurodegenerative diseases of later onset can, in the context of biallelic inheritance, cause devastating neurodevelopmental disorders and infantile or childhood-onset neurodegeneration. TOR1A-associated arthrogryposis multiplex congenita 5 (AMC5) is a rare neurodevelopmental disorder arising from biallelic variants in TOR1A, a gene that in the heterozygous state is associated with torsion dystonia-1 (DYT1 or DYT-TOR1A), an early-onset dystonia with reduced penetrance. While 15 individuals with AMC5-TOR1A have been reported (less than 10 in detail), a systematic investigation of the full disease-associated spectrum has not been conducted. Here, we assess the clinical, radiological and molecular characteristics of 57 individuals from 40 families with biallelic variants in TOR1A. Median age at last follow-up was 3 years (0-24 years). Most individuals presented with severe congenital flexion contractures (95%) and variable developmental delay (79%). Motor symptoms were reported in 79% and included lower limb spasticity and pyramidal signs, as well as gait disturbances. Facial dysmorphism was an integral part of the phenotype, with key features being a broad/full nasal tip, narrowing of the forehead and full cheeks. Analysis of disease-associated manifestations delineated a phenotypic spectrum ranging from normal cognition and mild gait disturbance to congenital arthrogryposis, global developmental delay, intellectual disability, absent speech and inability to walk. In a subset, the presentation was consistent with foetal akinesia deformation sequence with severe intrauterine abnormalities. Survival was 71%, with higher mortality in males. Death occurred at a median age of 1.2 months (1 week-9 years), due to respiratory failure, cardiac arrest or sepsis. Analysis of brain MRI studies identified non-specific neuroimaging features, including a hypoplastic corpus callosum (72%), foci of signal abnormality in the subcortical and periventricular white matter (55%), diffuse white matter volume loss (45%), mega cisterna magna (36%) and arachnoid cysts (27%). The molecular spectrum included 22 distinct variants, defining a mutational hotspot in the C-terminal domain of the Torsin-1A protein. Genotype-phenotype analysis revealed an association of missense variants in the 3-helix bundle domain to an attenuated phenotype, while missense variants near the Walker A/B motif as well as biallelic truncating variants were linked to early death. In summary, this systematic cross-sectional analysis of a large cohort of individuals with biallelic TOR1A variants across a wide age-range delineates the clinical and genetic spectrum of TOR1A-related autosomal-recessive disease and highlights potential predictors for disease severity and survival.

Exploring gastrointestinal health in MECP2 duplication syndrome

BACKGROUND

MECP2 duplication syndrome (MDS) is a rare neurogenetic syndrome caused by duplications of MECP2 at the Xq28 region. Although constipation and gastrointestinal reflux are reported in MDS, a comprehensive characterization of gastrointestinal health has not been fully explored.

METHODS

We conducted a parent survey to explore the characteristics of gastrointestinal health in individuals with MDS using a secure online registry and compared differences in gastrointestinal symptoms between individuals with MDS and those with Rett syndrome (RTT).

KEY RESULTS

One hundred six surveys were analyzed. Symptoms commonly associated with constipation occurred in 72% to 89% of MDS individuals. Eleven percent of MDS individuals underwent surgery for complications associated with constipation. We observed a bimodal distribution for gastroesophageal reflux disease (GERD) and gastrostomy feeding, with higher prevalence in 0-3 and >12-year-old MDS individuals. Constipation and GERD were significantly more common, and gas bloating was significantly less common in MDS than in RTT. Biliary tract disease requiring surgery was an unrecognized problem in 5% of MDS individuals. We determined that gastrointestinal problems in MDS individuals contribute to caretaker burden.

CONCLUSION AND INFERENCES

Our study is the first in-depth investigation that characterizes gastrointestinal health in MDS and enumerates differences in gastrointestinal symptoms between MDS and RTT. Strategies to reduce gastrointestinal symptoms will alleviate caregiver burden in MDS. Further studies are needed to examine the mechanisms that cause gastrointestinal problems in MDS.

Bi-allelic variants in HMGCR cause an autosomal-recessive progressive limb-girdle muscular dystrophy

Statins are a mainstay intervention for cardiovascular disease prevention, yet their use can cause rare severe myopathy. HMG-CoA reductase, an essential enzyme in the mevalonate pathway, is the target of statins. We identified nine individuals from five unrelated families with unexplained limb-girdle like muscular dystrophy and bi-allelic variants in HMGCR via clinical and research exome sequencing. The clinical features resembled other genetic causes of muscular dystrophy with incidental high CPK levels (>1,000 U/L), proximal muscle weakness, variable age of onset, and progression leading to impaired ambulation. Muscle biopsies in most affected individuals showed non-specific dystrophic changes with non-diagnostic immunohistochemistry. Molecular modeling analyses revealed variants to be destabilizing and affecting protein oligomerization. Protein activity studies using three variants (p.Asp623Asn, p.Tyr792Cys, and p.Arg443Gln) identified in affected individuals confirmed decreased enzymatic activity and reduced protein stability. In summary, we showed that individuals with bi-allelic amorphic (i.e., null and/or hypomorphic) variants in HMGCR display phenotypes that resemble non-genetic causes of myopathy involving this reductase. This study expands our knowledge regarding the mechanisms leading to muscular dystrophy through dysregulation of the mevalonate pathway, autoimmune myopathy, and statin-induced myopathy.

Bi-allelic variants in the ESAM tight-junction gene cause a neurodevelopmental disorder associated with fetal intracranial hemorrhage

The blood-brain barrier (BBB) is an essential gatekeeper for the central nervous system and incidence of neurodevelopmental disorders (NDDs) is higher in infants with a history of intracerebral hemorrhage (ICH). We discovered a rare disease trait in thirteen individuals, including four fetuses, from eight unrelated families associated with homozygous loss-of-function variant alleles of ESAM which encodes an endothelial cell adhesion molecule. The c.115del (p.Arg39Glyfs∗33) variant, identified in six individuals from four independent families of Southeastern Anatolia, severely impaired the in vitro tubulogenic process of endothelial colony-forming cells, recapitulating previous evidence in null mice, and caused lack of ESAM expression in the capillary endothelial cells of damaged brain. Affected individuals with bi-allelic ESAM variants showed profound global developmental delay/unspecified intellectual disability, epilepsy, absent or severely delayed speech, varying degrees of spasticity, ventriculomegaly, and ICH/cerebral calcifications, the latter being also observed in the fetuses. Phenotypic traits observed in individuals with bi-allelic ESAM variants overlap very closely with other known conditions characterized by endothelial dysfunction due to mutation of genes encoding tight junction molecules. Our findings emphasize the role of brain endothelial dysfunction in NDDs and contribute to the expansion of an emerging group of diseases that we propose to rename as "tightjunctionopathies."

Bi-allelic TTI1 variants cause an autosomal-recessive neurodevelopmental disorder with microcephaly

Telomere maintenance 2 (TELO2), Tel2 interacting protein 2 (TTI2), and Tel2 interacting protein 1 (TTI1) are the three components of the conserved Triple T (TTT) complex that modulates activity of phosphatidylinositol 3-kinase-related protein kinases (PIKKs), including mTOR, ATM, and ATR, by regulating the assembly of mTOR complex 1 (mTORC1). The TTT complex is essential for the expression, maturation, and stability of ATM and ATR in response to DNA damage. TELO2- and TTI2-related bi-allelic autosomal-recessive (AR) encephalopathies have been described in individuals with moderate to severe intellectual disability (ID), short stature, postnatal microcephaly, and a movement disorder (in the case of variants within TELO2). We present clinical, genomic, and functional data from 11 individuals in 9 unrelated families with bi-allelic variants in TTI1. All present with ID, and most with microcephaly, short stature, and a movement disorder. Functional studies performed in HEK293T cell lines and fibroblasts and lymphoblastoid cells derived from 4 unrelated individuals showed impairment of the TTT complex and of mTOR pathway activity which is improved by treatment with Rapamycin. Our data delineate a TTI1-related neurodevelopmental disorder and expand the group of disorders related to the TTT complex.

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 variants in HECT E3 paralogs, HECTD4 and UBE3C, encoding ubiquitin ligases cause neurodevelopmental disorders that overlap with Angelman syndrome

PURPOSE

Pathogenic variants in genes encoding ubiquitin E3 ligases are known to cause neurodevelopmental syndromes. Additional neurodevelopmental disorders associated with the other genes encoding E3 ligases are yet to be identified.

METHODS

Chromosomal analysis and exome sequencing were used to identify the genetic causes in 10 patients from 7 unrelated families with syndromic neurodevelopmental, seizure, and movement disorders and neurobehavioral phenotypes.

RESULTS

In total, 4 patients were found to have 3 different homozygous loss-of-function (LoF) variants, and 3 patients had 4 compound heterozygous missense variants in the candidate E3 ligase gene, HECTD4, that were rare, absent from controls as homozygous, and predicted to be deleterious in silico. In 3 patients from 2 families with Angelman-like syndrome, paralog-directed candidate gene approach detected 2 LoF variants in the other candidate E3 ligase gene, UBE3C, a paralog of the Angelman syndrome E3 ligase gene, UBE3A. The RNA studies in 4 patients with LoF variants in HECTD4 and UBE3C provided evidence for the LoF effect.

CONCLUSION

HECTD4 and UBE3C are novel biallelic rare disease genes, expand the association of the other HECT E3 ligase group with neurodevelopmental syndromes, and could explain some of the missing heritability in patients with a suggestive clinical diagnosis of Angelman syndrome.

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.

TCEAL1 loss-of-function results in an X-linked dominant neurodevelopmental syndrome and drives the neurological disease trait in Xq22.2 deletions

An Xq22.2 region upstream of PLP1 has been proposed to underly a neurological disease trait when deleted in 46,XX females. Deletion mapping revealed that heterozygous deletions encompassing the smallest region of overlap (SRO) spanning six Xq22.2 genes (BEX3, RAB40A, TCEAL4, TCEAL3, TCEAL1, and MORF4L2) associate with an early-onset neurological disease trait (EONDT) consisting of hypotonia, intellectual disability, neurobehavioral abnormalities, and dysmorphic facial features. None of the genes within the SRO have been associated with monogenic disease in OMIM. Through local and international collaborations facilitated by GeneMatcher and Matchmaker Exchange, we have identified and herein report seven de novo variants involving TCEAL1 in seven unrelated families: three hemizygous truncating alleles; one hemizygous missense allele; one heterozygous TCEAL1 full gene deletion; one heterozygous contiguous deletion of TCEAL1, TCEAL3, and TCEAL4; and one heterozygous frameshift variant allele. Variants were identified through exome or genome sequencing with trio analysis or through chromosomal microarray. Comparison with previously reported Xq22 deletions encompassing TCEAL1 identified a more-defined syndrome consisting of hypotonia, abnormal gait, developmental delay/intellectual disability especially affecting expressive language, autistic-like behavior, and mildly dysmorphic facial features. Additional features include strabismus, refractive errors, variable nystagmus, gastroesophageal reflux, constipation, dysmotility, recurrent infections, seizures, and structural brain anomalies. An additional maternally inherited hemizygous missense allele of uncertain significance was identified in a male with hypertonia and spasticity without syndromic features. These data provide evidence that TCEAL1 loss of function causes a neurological rare disease trait involving significant neurological impairment with features overlapping the EONDT phenotype in females with the Xq22 deletion.

Bi-allelic CAMSAP1 variants cause a clinically recognizable neuronal migration disorder

Non-centrosomal microtubules are essential cytoskeletal filaments that are important for neurite formation, axonal transport, and neuronal migration. They require stabilization by microtubule minus-end-targeting proteins including the CAMSAP family of molecules. Using exome sequencing on samples from five unrelated families, we show that bi-allelic CAMSAP1 loss-of-function variants cause a clinically recognizable, syndromic neuronal migration disorder. The cardinal clinical features of the syndrome include a characteristic craniofacial appearance, primary microcephaly, severe neurodevelopmental delay, cortical visual impairment, and seizures. The neuroradiological phenotype comprises a highly recognizable combination of classic lissencephaly with a posterior more severe than anterior gradient similar to PAFAH1B1(LIS1)-related lissencephaly and severe hypoplasia or absence of the corpus callosum; dysplasia of the basal ganglia, hippocampus, and midbrain; and cerebellar hypodysplasia, similar to the tubulinopathies, a group of monogenic tubulin-associated disorders of cortical dysgenesis. Neural cell rosette lineages derived from affected individuals displayed findings consistent with these phenotypes, including abnormal morphology, decreased cell proliferation, and neuronal differentiation. Camsap1-null mice displayed increased perinatal mortality, and RNAScope studies identified high expression levels in the brain throughout neurogenesis and in facial structures, consistent with the mouse and human neurodevelopmental and craniofacial phenotypes. Together our findings confirm a fundamental role of CAMSAP1 in neuronal migration and brain development and define bi-allelic variants as a cause of a clinically distinct neurodevelopmental disorder in humans and mice.

Loss of non-motor kinesin KIF26A causes congenital brain malformations via dysregulated neuronal migration and axonal growth as well as apoptosis

Kinesins are canonical molecular motors but can also function as modulators of intracellular signaling. KIF26A, an unconventional kinesin that lacks motor activity, inhibits growth-factor-receptor-bound protein 2 (GRB2)- and focal adhesion kinase (FAK)-dependent signal transduction, but its functions in the brain have not been characterized. We report a patient cohort with biallelic loss-of-function variants in KIF26A, exhibiting a spectrum of congenital brain malformations. In the developing brain, KIF26A is preferentially expressed during early- and mid-gestation in excitatory neurons. Combining mice and human iPSC-derived organoid models, we discovered that loss of KIF26A causes excitatory neuron-specific defects in radial migration, localization, dendritic and axonal growth, and apoptosis, offering a convincing explanation of the disease etiology in patients. Single-cell RNA sequencing in KIF26A knockout organoids revealed transcriptional changes in MAPK, MYC, and E2F pathways. Our findings illustrate the pathogenesis of KIF26A loss-of-function variants and identify the surprising versatility of this non-motor kinesin.

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.

Biallelic variants in ADAMTS15 cause a novel form of distal arthrogryposis

PURPOSE

We aimed to identify the underlying genetic cause for a novel form of distal arthrogryposis.

METHODS

Rare variant family-based genomics, exome sequencing, and disease-specific panel sequencing were used to detect ADAMTS15 variants in affected individuals. Adamts15 expression was analyzed at the single-cell level during murine embryogenesis. Expression patterns were characterized using in situ hybridization and RNAscope.

RESULTS

We identified homozygous rare variant alleles of ADAMTS15 in 5 affected individuals from 4 unrelated consanguineous families presenting with congenital flexion contractures of the interphalangeal joints and hypoplastic or absent palmar creases. Radiographic investigations showed physiological interphalangeal joint morphology. Additional features included knee, Achilles tendon, and toe contractures, spinal stiffness, scoliosis, and orthodontic abnormalities. Analysis of mouse whole-embryo single-cell sequencing data revealed a tightly regulated Adamts15 expression in the limb mesenchyme between embryonic stages E11.5 and E15.0. A perimuscular and peritendinous expression was evident in in situ hybridization in the developing mouse limb. In accordance, RNAscope analysis detected a significant coexpression with Osr1, but not with markers for skeletal muscle or joint formation.

CONCLUSION

In aggregate, our findings provide evidence that rare biallelic recessive trait variants in ADAMTS15 cause a novel autosomal recessive connective tissue disorder, resulting in a distal arthrogryposis syndrome.

A reverse genetics and genomics approach to gene paralog function and disease: Myokymia and the juxtaparanode

The leucine-rich glioma-inactivated (LGI) family consists of four highly conserved paralogous genes, LGI1-4, that are highly expressed in mammalian central and/or peripheral nervous systems. LGI1 antibodies are detected in subjects with autoimmune limbic encephalitis and peripheral nerve hyperexcitability syndromes (PNHSs) such as Isaacs and Morvan syndromes. Pathogenic variations of LGI1 and LGI4 are associated with neurological disorders as disease traits including familial temporal lobe epilepsy and neurogenic arthrogryposis multiplex congenita 1 with myelin defects, respectively. No human disease has been reported associated with either LGI2 or LGI3. We implemented exome sequencing and family-based genomics to identify individuals with deleterious variants in LGI3 and utilized GeneMatcher to connect practitioners and researchers worldwide to investigate the clinical and electrophysiological phenotype in affected subjects. We also generated Lgi3-null mice and performed peripheral nerve dissection and immunohistochemistry to examine the juxtaparanode LGI3 microarchitecture. As a result, we identified 16 individuals from eight unrelated families with loss-of-function (LoF) bi-allelic variants in LGI3. Deep phenotypic characterization showed LGI3 LoF causes a potentially clinically recognizable PNHS trait characterized by global developmental delay, intellectual disability, distal deformities with diminished reflexes, visible facial myokymia, and distinctive electromyographic features suggestive of motor nerve instability. Lgi3-null mice showed reduced and mis-localized Kv1 channel complexes in myelinated peripheral axons. Our data demonstrate bi-allelic LoF variants in LGI3 cause a clinically distinguishable disease trait of PNHS, most likely caused by disturbed Kv1 channel distribution in the absence of LGI3.

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.

Biallelic Variants in the Ectonucleotidase ENTPD1 Cause a Complex Neurodevelopmental Disorder with Intellectual Disability, Distinct White Matter Abnormalities, and Spastic Paraplegia

OBJECTIVE

Human genomics established that pathogenic variation in diverse genes can underlie a single disorder. For example, hereditary spastic paraplegia is associated with >80 genes, with frequently only few affected individuals described for each gene. Herein, we characterize a large cohort of individuals with biallelic variation in ENTPD1, a gene previously linked to spastic paraplegia 64 (Mendelian Inheritance in Man # 615683).

METHODS

Individuals with biallelic ENTPD1 variants were recruited worldwide. Deep phenotyping and molecular characterization were performed.

RESULTS

A total of 27 individuals from 17 unrelated families were studied; additional phenotypic information was collected from published cases. Twelve novel pathogenic ENTPD1 variants are described (NM 001776.6): c.398_399delinsAA; p.(Gly133Glu), c.540del; p.(Thr181Leufs*18), c.640del; p.(Gly216Glufs*75), c.185 T > G; p.(Leu62*), c.1531 T > C; p.(*511Glnext*100), c.967C > T; p.(Gln323*), c.414-2_414-1del, and c.146 A > G; p.(Tyr49Cys) including 4 recurrent variants c.1109 T > A; p.(Leu370*), c.574-6_574-3del, c.770_771del; p.(Gly257Glufs*18), and c.1041del; p.(Ile348Phefs*19). Shared disease traits include childhood onset, progressive spastic paraplegia, intellectual disability (ID), dysarthria, and white matter abnormalities. In vitro assays demonstrate that ENTPD1 expression and function are impaired and that c.574-6_574-3del causes exon skipping. Global metabolomics demonstrate ENTPD1 deficiency leads to impaired nucleotide, lipid, and energy metabolism.

INTERPRETATION

The ENTPD1 locus trait consists of childhood disease onset, ID, progressive spastic paraparesis, dysarthria, dysmorphisms, and white matter abnormalities, with some individuals showing neurocognitive regression. Investigation of an allelic series of ENTPD1 (1) expands previously described features of ENTPD1-related neurological disease, (2) highlights the importance of genotype-driven deep phenotyping, (3) documents the need for global collaborative efforts to characterize rare autosomal recessive disease traits, and (4) provides insights into disease trait neurobiology. ANN NEUROL 2022;92:304-321.

Assessing the Burden on Caregivers of MECP2 Duplication Syndrome

BACKGROUND

MECP2 duplication syndrome (MDS) is a rare neurogenetic disorder characterized by severe neurodevelopmental disorder, refractory epilepsy, recurrent infections, and functional gastrointestinal problems. Because of the significant clinical problems and lifelong disability of children with this disorder we hypothesized that the burden on parents/caregivers of these children is significant. However, there are no reports of the impact on caregivers of individuals with MDS.

METHODS

We developed and validated a burden scale to investigate the challenges of caregivers of children and adults with MDS and identified factors contributing to the burden on caregivers. We developed a Health Insurance Portability and Accountability Act-compliant patient registry for families with MDS and delivered caregiver burden survey through the registry.

RESULTS

Of 237 completed surveys, 101 were eligible for the study. We identified increased levels of self-perceived anxiety, depression, and emotional exhaustion in caregivers that correlated with higher burden scores. Epilepsy was the only clinical feature that caused a higher burden in caregivers of individuals with MDS. In addition, a higher burden was found in Hispanic caregivers. The duration of care negatively correlated with burden score.

CONCLUSIONS

This is the first study to investigate the burden on caregivers of individuals with MDS and identify several factors contributing to increased burden. Addressing these concerns has the potential to improve the health of individuals with MDS and contribute to the well-being of their caretakers.

Exploring the characteristics and most bothersome symptoms in MECP2 duplication syndrome to pave the path toward developing parent-oriented outcome measures

BACKGROUND

MECP2 Duplication Syndrome (MDS), resulting from the duplication of Xq28 region, including MECP2, is a rare disorder with a nascent understanding in clinical features and severity. Studies using antisense oligonucleotides revealed a broad phenotypic rescue in transgenic mice. With human clinical trials on the horizon, there is a need to develop clinical outcome measures for MDS.

METHODS

We surveyed caregivers of MDS individuals to explore the frequency and severity of MDS clinical features, and identify the most meaningful symptoms/domains that need to be included in the outcome measure scales.

RESULTS

A total of 101 responses were eligible for the survey. The top six most meaningful symptoms to caregivers in descending order included epilepsy, gross motor, fine motor, communication, infection, and constipation problems. Epilepsy was present in 58.4% of the subjects and 75% were drug-resistant, Furthermore, ~12% required intensive care unit (ICU) admission. Infections were present in 55% of the subjects, and one-fourth of them required ICU admission. Constipation was present in ~85% of the subjects and one-third required enemas/suppositories.

CONCLUSION

Our study is one of the largest cohorts conducted on MDS individuals characterizing the frequency and severity of MDS symptoms. Additionally, these study results will contribute to establishing a foundation to develop parent-reported outcomes in MDS.

Clinical Profile and Long-Term Outcomes in Pediatric Cavernous Sinus Thrombosis

BACKGROUND

Pediatric cavernous sinus thrombosis (CST) is a rare condition with limited data regarding its clinical characteristics and long-term outcomes. The objective of this study was to describe the clinical and radiologic features, diagnostic evaluation, management, and long-term prognosis and to identify clinical variables associated with long-term outcomes in pediatric CST.

METHODS

A retrospective chart review of patients younger than 18 years diagnosed with a CST between 2004 and 2018 at a single center was conducted.

RESULTS

We identified 16 (M:F = 10:6) children with CST with a mean age of 7.6 years (10 days to 15 years) and average follow-up duration of 29 months (3 weeks to 144 months). The most common symptom and examination finding at presentation was eyelid swelling (n = 8). Six patients had bilateral CST. The most common etiologies were sinusitis (n = 5) and orbital cellulitis (n = 5). Treatments included antibiotics (n = 14), anticoagulation (n = 11), and surgery (n = 5). Only one patient died due to intracranial complications. Twelve patients had a normal examination at follow-up. None of the clinical variables including age (P = 0.14), gender (P = 0.09), use of antibiotics (P = 1.00) or anticoagulation (P = 1.00), surgery (P = 0.28), parenchymal abnormalities (P = 0.30), additional cerebral venous thrombosis (P = 0.28), and early versus late commencing of anticoagulation (P = 1.00) were significant when comparing patients with full/partial resolution versus those with no resolution of thrombosis on follow-up neuroimaging.

CONCLUSIONS

Our study is one of the largest cohorts with the longest follow-up data for the pediatric CST. Most of our patients had favorable outcomes at follow-up. We found no statistical difference between clinical variables when comparing patients with full/partial resolution versus those with no resolution of thrombosis on follow-up neuroimaging.

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.

Biallelic variants in SLC38A3 encoding a glutamine transporter cause epileptic encephalopathy

The solute carrier (SLC) superfamily encompasses >400 transmembrane transporters involved in the exchange of amino acids, nutrients, ions, metals, neurotransmitters and metabolites across biological membranes. SLCs are highly expressed in the mammalian brain; defects in nearly 100 unique SLC-encoding genes (OMIM: https://www.omim.org) are associated with rare Mendelian disorders including developmental and epileptic encephalopathy and severe neurodevelopmental disorders. Exome sequencing and family-based rare variant analyses on a cohort with neurodevelopmental disorders identified two siblings with developmental and epileptic encephalopathy and a shared deleterious homozygous splicing variant in SLC38A3. The gene encodes SNAT3, a sodium-coupled neutral amino acid transporter and a principal transporter of the amino acids asparagine, histidine, and glutamine, the latter being the precursor for the neurotransmitters GABA and glutamate. Additional subjects with a similar developmental and epileptic encephalopathy phenotype and biallelic predicted-damaging SLC38A3 variants were ascertained through GeneMatcher and collaborations with research and clinical molecular diagnostic laboratories. Untargeted metabolomic analysis was performed to identify novel metabolic biomarkers. Ten individuals from seven unrelated families from six different countries with deleterious biallelic variants in SLC38A3 were identified. Global developmental delay, intellectual disability, hypotonia, and absent speech were common features while microcephaly, epilepsy, and visual impairment were present in the majority. Epilepsy was drug-resistant in half. Metabolomic analysis revealed perturbations of glutamate, histidine, and nitrogen metabolism in plasma, urine, and CSF of selected subjects, potentially representing biomarkers of disease. Our data support the contention that SLC38A3 is a novel disease gene for developmental and epileptic encephalopathy and illuminate the likely pathophysiology of the disease as perturbations in glutamine homeostasis.

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.

Quantitative dissection of multilocus pathogenic variation in an Egyptian infant with severe neurodevelopmental disorder resulting from multiple molecular diagnoses

Genomic sequencing and clinical genomics have demonstrated that substantial subsets of atypical and/or severe disease presentations result from multilocus pathogenic variation (MPV) causing blended phenotypes. In an infant with a severe neurodevelopmental disorder, four distinct molecular diagnoses were found by exome sequencing (ES). The blended phenotype that includes brain malformation, dysmorphism, and hypotonia was dissected using the Human Phenotype Ontology (HPO). ES revealed variants in CAPN3 (c.259C > G:p.L87V), MUSK (c.1781C > T:p.A594V), NAV2 (c.1996G > A:p.G666R), and ZC4H2 (c.595A > C:p.N199H). CAPN3, MUSK, and ZC4H2 are established disease genes linked to limb-girdle muscular dystrophy (OMIM# 253600), congenital myasthenia (OMIM# 616325), and Wieacker-Wolff syndrome (WWS; OMIM# 314580), respectively. NAV2 is a retinoic-acid responsive novel disease gene candidate with biological roles in neurite outgrowth and cerebellar dysgenesis in mouse models. Using semantic similarity, we show that no gene identified by ES individually explains the proband phenotype, but rather the totality of the clinically observed disease is explained by the combination of disease-contributing effects of the identified genes. These data reveal that multilocus pathogenic variation can result in a blended phenotype with each gene affecting a different part of the nervous system and nervous system-muscle connection. We provide evidence from this n = 1 study that in patients with MPV and complex blended phenotypes resulting from multiple molecular diagnoses, quantitative HPO analysis can allow for dissection of phenotypic contribution of both established disease genes and novel disease gene candidates not yet proven to cause human disease.

Molecular and clinical insights into complex genomic rearrangements related to MECP2 duplication syndrome

MECP2 duplication syndrome (MDS) is caused by copy number variation (CNV) spanning the MECP2 gene at Xq28 and is a major cause of intellectual disability (ID) in males. Herein, we describe two unrelated males harboring non-recurrent complex Xq28 rearrangements associated with MDS. Copy number gains were initially detected by quantitative real-time polymerase chain reaction and further delineated by high-resolution array comparative genomic hybridization, familial segregation, expression analysis and X-chromosome inactivation (XCI) evaluation in a carrier mother. SNVs within the rearrangements and/or fluorescent in situ hybridization (FISH) were used to assess the parental origin of the rearrangements. Patient 1 exhibited an intrachromosomal rearrangement, whose structure is consistent with a triplicated segment presumably embedded in an inverted orientation between two duplicated sequences (DUP-TRP/INV-DUP). The rearrangement was inherited from the carrier mother, who exhibits extreme XCI skewing and subtle psychiatric symptoms. Patient 2 presented a de novo (X;Y) unbalanced translocation resulting in duplication of Xq28 and deletion of Yp, originated in the paternal gametogenesis. Neurodevelopmental trajectory and non-neurological symptoms were consistent with previous reports, with the exception of cerebellar vermis hypoplasia in patient 2. Although both patients share the core MDS phenotype, patient 1 showed MECP2 transcript levels in blood similar to controls. Understanding the molecular mechanisms related to MDS is essential for designing targeted therapeutic strategies.

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.

Deep clinicopathological phenotyping identifies a previously unrecognized pathogenic EMD splice variant

Exome sequencing (ES) has revolutionized rare disease management, yet only ~25%-30% of patients receive a molecular diagnosis. A limiting factor is the quality of available phenotypic data. Here, we describe how deep clinicopathological phenotyping yielded a molecular diagnosis for a 19-year-old proband with muscular dystrophy and negative clinical ES. Deep phenotypic analysis identified two critical data points: (1) the absence of emerin protein in muscle biopsy and (2) clinical features consistent with Emery-Dreifuss muscular dystrophy. Sequencing data analysis uncovered an ultra-rare, intronic variant in EMD, the gene encoding emerin. The variant, NM_000117.3: c.188-6A > G, is predicted to impact splicing by in silico tools. This case thus illustrates how better integration of clinicopathologic data into ES analysis can enhance diagnostic yield with implications for clinical practice.

Biallelic and monoallelic variants in PLXNA1 are implicated in a novel neurodevelopmental disorder with variable cerebral and eye anomalies

PURPOSE

To investigate the effect of PLXNA1 variants on the phenotype of patients with autosomal dominant and recessive inheritance patterns and to functionally characterize the zebrafish homologs plxna1a and plxna1b during development.

METHODS

We assembled ten patients from seven families with biallelic or de novo PLXNA1 variants. We describe genotype-phenotype correlations, investigated the variants by structural modeling, and used Morpholino knockdown experiments in zebrafish to characterize the embryonic role of plxna1a and plxna1b.

RESULTS

Shared phenotypic features among patients include global developmental delay (9/10), brain anomalies (6/10), and eye anomalies (7/10). Notably, seizures were predominantly reported in patients with monoallelic variants. Structural modeling of missense variants in PLXNA1 suggests distortion in the native protein. Our zebrafish studies enforce an embryonic role of plxna1a and plxna1b in the development of the central nervous system and the eye.

CONCLUSION

We propose that different biallelic and monoallelic variants in PLXNA1 result in a novel neurodevelopmental syndrome mainly comprising developmental delay, brain, and eye anomalies. We hypothesize that biallelic variants in the extracellular Plexin-A1 domains lead to impaired dimerization or lack of receptor molecules, whereas monoallelic variants in the intracellular Plexin-A1 domains might impair downstream signaling through a dominant-negative effect.

Loss of C2orf69 defines a fatal autoinflammatory syndrome in humans and zebrafish that evokes a glycogen-storage-associated mitochondriopathy

Human C2orf69 is an evolutionarily conserved gene whose function is unknown. Here, we report eight unrelated families from which 20 children presented with a fatal syndrome consisting of severe autoinflammation and progredient leukoencephalopathy with recurrent seizures; 12 of these subjects, whose DNA was available, segregated homozygous loss-of-function C2orf69 variants. C2ORF69 bears homology to esterase enzymes, and orthologs can be found in most eukaryotic genomes, including that of unicellular phytoplankton. We found that endogenous C2ORF69 (1) is loosely bound to mitochondria, (2) affects mitochondrial membrane potential and oxidative respiration in cultured neurons, and (3) controls the levels of the glycogen branching enzyme 1 (GBE1) consistent with a glycogen-storage-associated mitochondriopathy. We show that CRISPR-Cas9-mediated inactivation of zebrafish C2orf69 results in lethality by 8 months of age due to spontaneous epileptic seizures, which is preceded by persistent brain inflammation. Collectively, our results delineate an autoinflammatory Mendelian disorder of C2orf69 deficiency that disrupts the development/homeostasis of the immune and central nervous systems.

Bi-allelic loss-of-function variants in BCAS3 cause a syndromic neurodevelopmental disorder

BCAS3 microtubule-associated cell migration factor (BCAS3) is a large, highly conserved cytoskeletal protein previously proposed to be critical in angiogenesis and implicated in human embryogenesis and tumorigenesis. Here, we established BCAS3 loss-of-function variants as causative for a neurodevelopmental disorder. We report 15 individuals from eight unrelated families with germline bi-allelic loss-of-function variants in BCAS3. All probands share a global developmental delay accompanied by pyramidal tract involvement, microcephaly, short stature, strabismus, dysmorphic facial features, and seizures. The human phenotype is less severe compared with the Bcas3 knockout mouse model and cannot be explained by angiogenic defects alone. Consistent with being loss-of-function alleles, we observed absence of BCAS3 in probands' primary fibroblasts. By comparing the transcriptomic and proteomic data based on probands' fibroblasts with those of the knockout mouse model, we identified similar dysregulated pathways resulting from over-representation analysis, while the dysregulation of some proposed key interactors could not be confirmed. Together with the results from a tissue-specific Drosophila loss-of-function model, we demonstrate a vital role for BCAS3 in neural tissue development.