Front Cover, Volume 40, Issue 8

RING1 missense variants reveal sensitivity of DNA damage repair to H2A monoubiquitination dosage during neurogenesis

Polycomb repressive complex 1 (PRC1) modifies chromatin through catalysis of histone H2A lysine 119 monoubiquitination (H2AK119ub1). RING1 and RNF2 interchangeably serve as the catalytic subunit within PRC1. Pathogenic missense variants in PRC1 core components reveal functions of these proteins that are obscured in knockout models. While Ring1a knockout models remain healthy, the microcephaly and neuropsychiatric phenotypes associated with a pathogenic RING1 missense variant implicate unappreciated functions. Using an in vitro model of neurodevelopment, we observe that RING1 contributes to the broad placement of H2AK119ub1, and that its targets overlap with those of RNF2. PRC1 complexes harboring hypomorphic RING1 bind target loci but do not catalyze H2AK119ub1, reducing H2AK119ub1 by preventing catalytically active complexes from accessing the locus. This results in delayed DNA damage repair and cell cycle progression in neural progenitor cells (NPCs). Conversely, reduced H2AK119ub1 due to hypomorphic RING1 does not generate differential expression that impacts NPC differentiation. In contrast, hypomorphic RNF2 generates a greater reduction in H2AK119ub1 that results in both delayed DNA repair and widespread transcriptional changes. These findings suggest that the DNA damage response is more sensitive to H2AK119ub1 dosage change than is regulation of gene expression.

Biallelic loss-of-function variations in BTD cause profound biotinidase deficiency in an Indian patient

BACKGROUND

Biotinidase deficiency (BD) is a rare, autosomal recessive metabolic disorder characterized by neurocutaneous symptoms. This study investigates a case of profound BD in an Indian infant and the underlying genetic basis.

METHODS

A 10-month-old male presenting with seizures, hypotonia, ataxia, visual impairments, and developmental delay underwent biochemical and genetic analysis. Biotinidase activity was measured using an ELISA kit. Sanger sequencing of the biotinidase (BTD) gene was performed to identify genetic variations. In silico analysis was employed to assess the potential impact of the identified variants.

RESULTS

The infant biotinidase activity was undetectable and its suggest profound biotinidase deficiency. Novel biallelic loss-of-function variations (c.903G > A and c.946 C > T) in the BTD gene were identified, leading to premature stop codons and truncated, non-functional protein fragments.

CONCLUSION

This case expands our knowledge of BD genetic diversity and underscores the critical role of early diagnosis and newborn screening programs in managing this treatable condition.

Exome sequencing reveals neurodevelopmental genes in simplex consanguineous Iranian families with syndromic autism

BACKGROUND AND OBJECTIVE

Autosomal recessive genetic disorders pose significant health challenges in regions where consanguineous marriages are prevalent. The utilization of exome sequencing as a frequently employed methodology has enabled a clear delineation of diagnostic efficacy and mode of inheritance within multiplex consanguineous families. However, these aspects remain less elucidated within simplex families.

METHODS

In this study involving 12 unrelated simplex Iranian families presenting syndromic autism, we conducted singleton exome sequencing. The identified genetic variants were validated using Sanger sequencing, and for the missense variants in FOXG1 and DMD, 3D protein structure modeling was carried out to substantiate their pathogenicity. To examine the expression patterns of the candidate genes in the fetal brain, adult brain, and muscle, RT-qPCR was employed.

RESULTS

In four families, we detected an autosomal dominant gene (FOXG1), an autosomal recessive gene (CHKB), and two X-linked autism genes (IQSEC2 and DMD), indicating diverse inheritance patterns. In the remaining eight families, we were unable to identify any disease-associated genes. As a result, our variant detection rate stood at 33.3% (4/12), surpassing rates reported in similar studies of smaller cohorts. Among the four newly identified coding variants, three are de novo (heterozygous variant p.Trp546Ter in IQSEC2, heterozygous variant p.Ala188Glu in FOXG1, and hemizygous variant p.Leu211Met in DMD), while the homozygous variant p.Glu128Ter in CHKB was inherited from both healthy heterozygous parents. 3D protein structure modeling was carried out for the missense variants in FOXG1 and DMD, which predicted steric hindrance and spatial inhibition, respectively, supporting the pathogenicity of these human mutants. Additionally, the nonsense variant in CHKB is anticipated to influence its dimerization - crucial for choline kinase function - and the nonsense variant in IQSEC2 is predicted to eliminate three functional domains. Consequently, these distinct variants found in four unrelated individuals with autism are likely indicative of loss-of-function mutations.

CONCLUSIONS

In our two syndromic autism families, we discovered variants in two muscular dystrophy genes, DMD and CHKB. Given that DMD and CHKB are recognized for their participation in the non-cognitive manifestations of muscular dystrophy, it indicates that some genes transcend the boundary of apparently unrelated clinical categories, thereby establishing a novel connection between ASD and muscular dystrophy. Our findings also shed light on the complex inheritance patterns observed in Iranian consanguineous simplex families and emphasize the connection between autism spectrum disorder and muscular dystrophy. This underscores a likely genetic convergence between neurodevelopmental and neuromuscular disorders.

Genetic constraint at single amino acid resolution in protein domains improves missense variant prioritisation and gene discovery

BACKGROUND

One of the major hurdles in clinical genetics is interpreting the clinical consequences associated with germline missense variants in humans. Recent significant advances have leveraged natural variation observed in large-scale human populations to uncover genes or genomic regions that show a depletion of natural variation, indicative of selection pressure. We refer to this as "genetic constraint". Although existing genetic constraint metrics have been demonstrated to be successful in prioritising genes or genomic regions associated with diseases, their spatial resolution is limited in distinguishing pathogenic variants from benign variants within genes.

METHODS

We aim to identify missense variants that are significantly depleted in the general human population. Given the size of currently available human populations with exome or genome sequencing data, it is not possible to directly detect depletion of individual missense variants, since the average expected number of observations of a variant at most positions is less than one. We instead focus on protein domains, grouping homologous variants with similar functional impacts to examine the depletion of natural variations within these comparable sets. To accomplish this, we develop the Homologous Missense Constraint (HMC) score. We utilise the Genome Aggregation Database (gnomAD) 125 K exome sequencing data and evaluate genetic constraint at quasi amino-acid resolution by combining signals across protein homologues.

RESULTS

We identify one million possible missense variants under strong negative selection within protein domains. Though our approach annotates only protein domains, it nonetheless allows us to assess 22% of the exome confidently. It precisely distinguishes pathogenic variants from benign variants for both early-onset and adult-onset disorders. It outperforms existing constraint metrics and pathogenicity meta-predictors in prioritising de novo mutations from probands with developmental disorders (DD). It is also methodologically independent of these, adding power to predict variant pathogenicity when used in combination. We demonstrate utility for gene discovery by identifying seven genes newly significantly associated with DD that could act through an altered-function mechanism.

CONCLUSIONS

Grouping variants of comparable functional impacts is effective in evaluating their genetic constraint. HMC is a novel and accurate predictor of missense consequence for improved variant interpretation.

Variants of NAV3, a neuronal morphogenesis protein, cause intellectual disability, developmental delay, and microcephaly

Microtubule associated proteins (MAPs) are widely expressed in the central nervous system, and have established roles in cell proliferation, myelination, neurite formation, axon specification, outgrowth, dendrite, and synapse formation. We report eleven individuals from seven families harboring predicted pathogenic biallelic, de novo, and heterozygous variants in the NAV3 gene, which encodes the microtubule positive tip protein neuron navigator 3 (NAV3). All affected individuals have intellectual disability (ID), microcephaly, skeletal deformities, ocular anomalies, and behavioral issues. In mouse brain, Nav3 is expressed throughout the nervous system, with more prominent signatures in postmitotic, excitatory, inhibiting, and sensory neurons. When overexpressed in HEK293T and COS7 cells, pathogenic variants impaired NAV3 ability to stabilize microtubules. Further, knocking-down nav3 in zebrafish led to severe morphological defects, microcephaly, impaired neuronal growth, and behavioral impairment, which were rescued with co-injection of WT NAV3 mRNA and not by transcripts encoding the pathogenic variants. Our findings establish the role of NAV3 in neurodevelopmental disorders, and reveal its involvement in neuronal morphogenesis, and neuromuscular responses.

Low-frequency and rare genetic variants associated with rheumatoid arthritis risk

Rheumatoid arthritis (RA) has an estimated heritability of nearly 50%, which is particularly high in seropositive RA. HLA alleles account for a large proportion of this heritability, in addition to many common single-nucleotide polymorphisms with smaller individual effects. Low-frequency and rare variants, such as those captured by next-generation sequencing, can also have a large role in heritability in some individuals. Rare variant discovery has informed the development of drugs such as inhibitors of PCSK9 and Janus kinases. Some 34 low-frequency and rare variants are currently associated with RA risk. One variant (19:10352442G>C in TYK2) was identified in five separate studies, and might therefore represent a promising therapeutic target. Following a set of best practices in future studies, including studying diverse populations, using large sample sizes, validating RA and serostatus, replicating findings, adjusting for other variants and performing functional assessment, could help to ensure the relevance of identified variants. Exciting opportunities are now on the horizon for genetics in RA, including larger datasets and consortia, whole-genome sequencing and direct applications of findings in the management, and especially treatment, of RA.

Genetic Analysis and Functional Assessment of a TGFBR2 Variant in Micrognathia and Cleft Palate

Cleft lip and cleft palate are among the most common congenital anomalies and are the result of incomplete fusion of embryonic craniofacial processes or palatal shelves, respectively. We know that genetics play a large role in these anomalies but the list of known causal genes is far from complete. As part of a larger sequencing effort of patients with micrognathia and cleft palate we identified a candidate variant in transforming growth factor beta receptor 2 ( TGFBR2 ) which is rare, changing a highly conserved amino acid, and predicted to be pathogenic by a number of metrics. The family history and population genetics would suggest this specific variant would be incompletely penetrant, but this gene has been convincingly implicated in craniofacial development. In order to test the hypothesis this might be a causal variant, we used genome editing to create the orthologous variant in a new mouse model. Surprisingly, Tgfbr2 V387M mice did not exhibit craniofacial anomalies or have reduced survival suggesting this is, in fact, not a causal variant for cleft palate/ micrognathia. The discrepancy between in silico predictions and mouse phenotypes highlights the complexity of translating human genetic findings to mouse models. We expect these findings will aid in interpretation of future variants seen in TGFBR2 from ongoing sequencing of patients with congenital craniofacial anomalies.

A unified analysis of evolutionary and population constraint in protein domains highlights structural features and pathogenic sites

Protein evolution is constrained by structure and function, creating patterns in residue conservation that are routinely exploited to predict structure and other features. Similar constraints should affect variation across individuals, but it is only with the growth of human population sequencing that this has been tested at scale. Now, human population constraint has established applications in pathogenicity prediction, but it has not yet been explored for structural inference. Here, we map 2.4 million population variants to 5885 protein families and quantify residue-level constraint with a new Missense Enrichment Score (MES). Analysis of 61,214 structures from the PDB spanning 3661 families shows that missense depleted sites are enriched in buried residues or those involved in small-molecule or protein binding. MES is complementary to evolutionary conservation and a combined analysis allows a new classification of residues according to a conservation plane. This approach finds functional residues that are evolutionarily diverse, which can be related to specificity, as well as family-wide conserved sites that are critical for folding or function. We also find a possible contrast between lethal and non-lethal pathogenic sites, and a surprising clinical variant hot spot at a subset of missense enriched positions.

Heterozygous loss-of-function variants in DOCK4 cause neurodevelopmental delay and microcephaly

Neurons form the basic anatomical and functional structure of the nervous system, and defects in neuronal differentiation or formation of neurites are associated with various psychiatric and neurodevelopmental disorders. Dynamic changes in the cytoskeleton are essential for this process, which is, inter alia, controlled by the dedicator of cytokinesis 4 (DOCK4) through the activation of RAC1. Here, we clinically describe 7 individuals (6 males and one female) with variants in DOCK4 and overlapping phenotype of mild to severe global developmental delay. Additional symptoms include coordination or gait abnormalities, microcephaly, nonspecific brain malformations, hypotonia and seizures. Four individuals carry missense variants (three of them detected de novo) and three individuals carry null variants (two of them maternally inherited). Molecular modeling of the heterozygous missense variants suggests that the majority of them affect the globular structure of DOCK4. In vitro functional expression studies in transfected Neuro-2A cells showed that all missense variants impaired neurite outgrowth. Furthermore, Dock4 knockout Neuro-2A cells also exhibited defects in promoting neurite outgrowth. Our results, including clinical, molecular and functional data, suggest that loss-of-function variants in DOCK4 probable cause a variable spectrum of a novel neurodevelopmental disorder with microcephaly.

Bi-allelic variants in CELSR3 are implicated in central nervous system and urinary tract anomalies

CELSR3 codes for a planar cell polarity protein. We describe twelve affected individuals from eleven independent families with bi-allelic variants in CELSR3. Affected individuals presented with an overlapping phenotypic spectrum comprising central nervous system (CNS) anomalies (7/12), combined CNS anomalies and congenital anomalies of the kidneys and urinary tract (CAKUT) (3/12) and CAKUT only (2/12). Computational simulation of the 3D protein structure suggests the position of the identified variants to be implicated in penetrance and phenotype expression. CELSR3 immunolocalization in human embryonic urinary tract and transient suppression and rescue experiments of Celsr3 in fluorescent zebrafish reporter lines further support an embryonic role of CELSR3 in CNS and urinary tract formation.

Next-generation sequencing and bioinformatics in rare movement disorders

The ability to sequence entire exomes and genomes has revolutionized molecular testing in rare movement disorders, and genomic sequencing is becoming an integral part of routine diagnostic workflows for these heterogeneous conditions. However, interpretation of the extensive genomic variant information that is being generated presents substantial challenges. In this Perspective, we outline multidimensional strategies for genetic diagnosis in patients with rare movement disorders. We examine bioinformatics tools and computational metrics that have been developed to facilitate accurate prioritization of disease-causing variants. Additionally, we highlight community-driven data-sharing and case-matchmaking platforms, which are designed to foster the discovery of new genotype-phenotype relationships. Finally, we consider how multiomic data integration might optimize diagnostic success by combining genomic, epigenetic, transcriptomic and/or proteomic profiling to enable a more holistic evaluation of variant effects. Together, the approaches that we discuss offer pathways to the improved understanding of the genetic basis of rare movement disorders.

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.

Variants in the WDR44 WD40-repeat domain cause a spectrum of ciliopathy by impairing ciliogenesis initiation

WDR44 prevents ciliogenesis initiation by regulating RAB11-dependent vesicle trafficking. Here, we describe male patients with missense and nonsense variants within the WD40 repeats (WDR) of WDR44, an X-linked gene product, who display ciliopathy-related developmental phenotypes that we can model in zebrafish. The patient phenotypic spectrum includes developmental delay/intellectual disability, hypotonia, distinct craniofacial features and variable presence of brain, renal, cardiac and musculoskeletal abnormalities. We demonstrate that WDR44 variants associated with more severe disease impair ciliogenesis initiation and ciliary signaling. Because WDR44 negatively regulates ciliogenesis, it was surprising that pathogenic missense variants showed reduced abundance, which we link to misfolding of WDR autonomous repeats and degradation by the proteasome. We discover that disease severity correlates with increased RAB11 binding, which we propose drives ciliogenesis initiation dysregulation. Finally, we discover interdomain interactions between the WDR and NH2-terminal region that contains the RAB11 binding domain (RBD) and show patient variants disrupt this association. This study provides new insights into WDR44 WDR structure and characterizes a new syndrome that could result from impaired ciliogenesis.

Novel inherited CDX2 variant segregating in a family with diverse congenital malformations of the genitourinary system

Anorectal malformations (ARMs) constitute a group of congenital defects of the gastrointestinal and urogenital systems. They affect males and females, with an estimated worldwide prevalence of 1 in 5000 live births. These malformations are clinically heterogeneous and can be part of a syndromic presentation (syndromic ARM) or as a nonsyndromic entity (nonsyndromic ARM). Despite the well-recognized heritability of nonsyndromic ARM, the genetic etiology in most patients is unknown. In this study, we describe three siblings with diverse congenital anomalies of the genitourinary system, anemia, delayed milestones, and skeletal anomalies. Genome sequencing identified a novel, paternally inherited heterozygous Caudal type Homeobox 2 (CDX2) variant (c.722A > G (p.Glu241Gly)), that was present in all three affected siblings. The variant identified in this family is absent from population databases and predicted to be damaging by most in silico pathogenicity tools. So far, only two other reports implicate variants in CDX2 with ARMs. Remarkably, the individuals described in these studies had similar clinical phenotypes and genetic alterations in CDX2 CDX2 encodes a transcription factor and is considered the master regulator of gastrointestinal development. This variant maps to the homeobox domain of the encoded protein, which is critical for interaction with DNA targets. Our finding provides a potential molecular diagnosis for this family's condition and supports the role of CDX2 in anorectal anomalies. It also highlights the clinical heterogeneity and variable penetrance of ARM predisposition variants, another well-documented phenomenon. Finally, it underscores the diagnostic utility of genomic profiling of ARMs to identify the genetic etiology of these defects.

Skewed X-chromosome inactivation in unsolved neurodevelopmental disease cases can guide re-evaluation For X-linked genes

Despite major advances in genome technology and analysis, >50% of patients with a neurodevelopmental disorder (NDD) remain undiagnosed after extensive evaluation. A point in case is our clinically heterogeneous cohort of NDD patients that remained undiagnosed after FRAXA testing, chromosomal microarray analysis and trio exome sequencing (ES). In this study, we explored the frequency of non-random X chromosome inactivation (XCI) in the mothers of male patients and affected females, the rationale being that skewed XCI might be masking previously discarded genetic variants found on the X chromosome. A multiplex fluorescent PCR-based assay was used to analyse the pattern of XCI after digestion with HhaI methylation-sensitive restriction enzyme. In families with skewed XCI, we re-evaluated trio-based ES and identified pathogenic variants and a deletion on the X chromosome. Linkage analysis and RT-PCR were used to further study the inactive X chromosome allele, and Xdrop long-DNA technology was used to define chromosome deletion boundaries. We found skewed XCI (>90%) in 16/186 (8.6%) mothers of NDD males and in 12/90 (13.3%) NDD females, far beyond the expected rate of XCI in the normal population (3.6%, OR = 4.10; OR = 2.51). By re-analyzing ES and clinical data, we solved 7/28 cases (25%) with skewed XCI, identifying variants in KDM5C, PDZD4, PHF6, TAF1, OTUD5 and ZMYM3, and a deletion in ATRX. We conclude that XCI profiling is a simple assay that targets a subgroup of patients that can benefit from re-evaluation of X-linked variants, thus improving the diagnostic yield in NDD patients and identifying new X-linked disorders.

Genetic insights into PHARC syndrome: identification of a novel frameshift mutation in ABHD12

BACKGROUND

Mutations in ABHD12 (OMIM: 613,599) are associated with polyneuropathy, hearing loss, ataxia, retinitis pigmentosa, and cataract (PHARC) syndrome (OMIM: 612674), which is a rare autosomal recessive neurodegenerative disease. PHARC syndrome is easily misdiagnosed as other neurologic disorders, such as retinitis pigmentosa, Charcot-Marie-Tooth disease, and Refsum disease, due to phenotype variability and slow progression. This paper presents a novel mutation in ABHD12 in two affected siblings with PHARC syndrome phenotypes. In addition, we summarize genotype-phenotype information of the previously reported patients with ABHD12 mutation.

METHODS

Following a thorough medical evaluation, whole-exome sequencing was done on the proband to look for potential genetic causes. This was followed by confirmation of identified variant in the proband and segregation analysis in the family by Sanger sequencing. The variants were interpreted based on the American College of Medical Genetics and Genomics (ACMG) guidelines.

RESULTS

A novel pathogenic homozygous frameshift variant, NM_001042472.3:c.601dup, p.(Val201GlyfsTer4), was identified in exon 6 of ABHD12 (ACMG criteria: PVS1 and PM2, PM1, PM4, PP3, and PP4). Through Sanger sequencing, we showed that this variant is co-segregated with the disease in the family. Further medical evaluations confirmed the compatibility of the patients' phenotype with PHARC syndrome.

CONCLUSIONS

Our findings expand the spectrum of mutations in the ABHD12 and emphasize the significance of multidisciplinary diagnostic collaboration among clinicians and geneticists to solve the differential diagnosis of related disorders. Moreover, a summary based on mutations found so far in the ABHD12 gene did not suggest a clear genotype-phenotype correlation for PHARC syndrome.

ARF1-related disorder: phenotypic and molecular spectrum

PURPOSE

ARF1 was previously implicated in periventricular nodular heterotopia (PVNH) in only five individuals and systematic clinical characterisation was not available. The aim of this study is to provide a comprehensive description of the phenotypic and genotypic spectrum of ARF1-related neurodevelopmental disorder.

METHODS

We collected detailed phenotypes of an international cohort of individuals (n=17) with ARF1 variants assembled through the GeneMatcher platform. Missense variants were structurally modelled, and the impact of several were functionally validated.

RESULTS

De novo variants (10 missense, 1 frameshift, 1 splice altering resulting in 9 residues insertion) in ARF1 were identified among 17 unrelated individuals. Detailed phenotypes included intellectual disability (ID), microcephaly, seizures and PVNH. No specific facial characteristics were consistent across all cases, however microretrognathia was common. Various hearing and visual defects were recurrent, and interestingly, some inflammatory features were reported. MRI of the brain frequently showed abnormalities consistent with a neuronal migration disorder.

CONCLUSION

We confirm the role of ARF1 in an autosomal dominant syndrome with a phenotypic spectrum including severe ID, microcephaly, seizures and PVNH due to impaired neuronal migration.

Comparative yield of molecular diagnostic algorithms for autism spectrum disorder diagnosis in India: evidence supporting whole exome sequencing as first tier test

BACKGROUND

Autism spectrum disorder (ASD) affects 1 in 100 children globally with a rapidly increasing prevalence. To the best of our knowledge, no data exists on the genetic architecture of ASD in India. This study aimed to identify the genetic architecture of ASD in India and to assess the use of whole exome sequencing (WES) as a first-tier test instead of chromosomal microarray (CMA) for genetic diagnosis.

METHODS

Between 2020 and 2022, 101 patient-parent trios of Indian origin diagnosed with ASD according to the Diagnostic and Statistical Manual, 5th edition, were recruited. All probands underwent a sequential genetic testing pathway consisting of karyotyping, Fragile-X testing (in male probands only), CMA and WES. Candidate variant validation and parental segregation analysis was performed using orthogonal methods.

RESULTS

Of 101 trios, no probands were identified with a gross chromosomal anomaly or Fragile-X. Three (2.9%) and 30 (29.7%) trios received a confirmed genetic diagnosis from CMA and WES, respectively. Amongst diagnosis from WES, SNVs were detected in 27 cases (90%) and CNVs in 3 cases (10%), including the 3 CNVs detected from CMA. Segregation analysis showed 66.6% (n = 3 for CNVs and n = 17 for SNVs) and 16.6% (n = 5) of the cases had de novo and recessive variants respectively, which is in concordance with the distribution of variant types and mode of inheritance observed in ASD patients of non-Hispanic white/ European ethnicity. MECP2 gene was the most recurrently mutated gene (n = 6; 20%) in the present cohort. Majority of the affected genes identified in the study cohort are involved in synaptic formation, transcription and its regulation, ubiquitination and chromatin remodeling.

CONCLUSIONS

Our study suggests de novo variants as a major cause of ASD in the Indian population, with Rett syndrome as the most commonly detected disorder. Furthermore, we provide evidence of a significant difference in the diagnostic yield between CMA (3%) and WES (30%) which supports the implementation of WES as a first-tier test for genetic diagnosis of ASD in India.

Gain and loss of function variants in EZH1 disrupt neurogenesis and cause dominant and recessive neurodevelopmental disorders

Genetic variants in chromatin regulators are frequently found in neurodevelopmental disorders, but their effect in disease etiology is rarely determined. Here, we uncover and functionally define pathogenic variants in the chromatin modifier EZH1 as the cause of dominant and recessive neurodevelopmental disorders in 19 individuals. EZH1 encodes one of the two alternative histone H3 lysine 27 methyltransferases of the PRC2 complex. Unlike the other PRC2 subunits, which are involved in cancers and developmental syndromes, the implication of EZH1 in human development and disease is largely unknown. Using cellular and biochemical studies, we demonstrate that recessive variants impair EZH1 expression causing loss of function effects, while dominant variants are missense mutations that affect evolutionarily conserved aminoacids, likely impacting EZH1 structure or function. Accordingly, we found increased methyltransferase activity leading to gain of function of two EZH1 missense variants. Furthermore, we show that EZH1 is necessary and sufficient for differentiation of neural progenitor cells in the developing chick embryo neural tube. Finally, using human pluripotent stem cell-derived neural cultures and forebrain organoids, we demonstrate that EZH1 variants perturb cortical neuron differentiation. Overall, our work reveals a critical role of EZH1 in neurogenesis regulation and provides molecular diagnosis for previously undefined neurodevelopmental disorders.

De novo MCM6 variants in neurodevelopmental disorders: a recognizable phenotype related to zinc binding residues

The minichromosome maintenance (MCM) complex acts as a DNA helicase during DNA replication, and thereby regulates cell cycle progression and proliferation. In addition, MCM-complex components localize to centrosomes and play an independent role in ciliogenesis. Pathogenic variants in genes coding for MCM components and other DNA replication factors have been linked to growth and developmental disorders as Meier-Gorlin syndrome and Seckel syndrome. Trio exome/genome sequencing identified the same de novo MCM6 missense variant p.(Cys158Tyr) in two unrelated individuals that presented with overlapping phenotypes consisting of intra-uterine growth retardation, short stature, congenital microcephaly, endocrine features, developmental delay and urogenital anomalies. The identified variant affects a zinc binding cysteine in the MCM6 zinc finger signature. This domain, and specifically cysteine residues, are essential for MCM-complex dimerization and the induction of helicase activity, suggesting a deleterious effect of this variant on DNA replication. Fibroblasts derived from the two affected individuals showed defects both in ciliogenesis and cell proliferation. We additionally traced three unrelated individuals with de novo MCM6 variants in the oligonucleotide binding (OB)-fold domain, presenting with variable (neuro)developmental features including autism spectrum disorder, developmental delay, and epilepsy. Taken together, our findings implicate de novo MCM6 variants in neurodevelopmental disorders. The clinical features and functional defects related to the zinc binding residue resemble those observed in syndromes related to other MCM components and DNA replication factors, while de novo OB-fold domain missense variants may be associated with more variable neurodevelopmental phenotypes. These data encourage consideration of MCM6 variants in the diagnostic arsenal of NDD.

Differential effects of Mendelian GDAP1 clinical variants on mitochondria-lysosome membrane contacts sites

GDAP1 pathogenic variants cause Charcot-Marie-Tooth (CMT) disease, the most common hereditary motor and sensory neuropathy. CMT-GDAP1 can be axonal or demyelinating, with autosomal dominant or recessive inheritance, leading to phenotypic heterogeneity. Recessive GDAP1 variants cause a severe phenotype, whereas dominant variants are associated with a milder disease course. GDAP1 is an outer mitochondrial membrane protein involved in mitochondrial membrane contact sites (MCSs) with the plasmatic membrane, the endoplasmic reticulum (ER), and lysosomes. In GDAP1-deficient models, the pathophysiology includes morphological defects in mitochondrial network and ER, impaired Ca2+ homeostasis, oxidative stress, and mitochondrial MCSs defects. Nevertheless, the underlying pathophysiology of dominant variants is less understood. Here, we study the effect upon mitochondria-lysosome MCSs of two GDAP1 clinical variants located in the α-loop interaction domain of the protein. p.Thr157Pro dominant variant causes the increase in these MCSs that correlates with a hyper-fissioned mitochondrial network. In contrast, p.Arg161His recessive variant, which is predicted to significantly change the contact surface of GDAP1, causes decreased contacts with more elongated mitochondria. Given that mitochondria-lysosome MCSs regulate Ca2+ transfer from the lysosome to mitochondria, our results support that GDAP1 clinical variants have different consequences for Ca2+ handling and that could be primary insults determining differences in severity between dominant and recessive forms of the disease.

Insights Into PROM1-Macular Disease Using Multimodal Imaging

Purpose

To describe the features of genetically confirmed PROM1-macular dystrophy in multimodal images.

Methods

Thirty-six (36) eyes of 18 patients (5-66 years; mean age, 42.4 years) were prospectively studied by clinical examination and multimodal imaging. Short-wavelength autofluorescence (SW-AF) and quantitative fundus autofluorescence (qAF) images were acquired with a scanning laser ophthalmoscope (HRA+OCT, Heidelberg Engineering) modified by insertion of an internal autofluorescent reference. Further clinical testing included near-infrared autofluorescence (NIR-AF; HRA2, Heidelberg Engineering) with semiquantitative analysis, spectral domain-optical coherence tomography (HRA+OCT) and full-field electroretinography. All patients were genetically confirmed by exome sequencing.

Results

All 18 patients presented with varying degrees of maculopathy. One family with individuals affected across two generations exhibited granular fleck-like deposits across the posterior pole. Areas of granular deposition in SW-AF and NIR-AF corresponded to intermittent loss of the ellipsoid zone, whereas discrete regions of hypoautofluorescence corresponded with a loss of outer retinal layers in spectral-domain optical coherence tomography scans. For 18 of the 20 eyes, qAF levels within the macula were within the 95% confidence intervals of healthy age-matched individuals; nor was the mean NIR-AF signal increased relative to healthy eyes.

Conclusions

Although PROM1-macular dystrophy (Stargardt disease 4) can exhibit phenotypic overlap with recessive Stargardt disease, significantly increased SW-AF levels were not detected. As such, elevated bisretinoid lipofuscin may not be a feature of the pathophysiology of PROM1 disease. The qAF approach could serve as a method of early differential diagnosis and may help to identify appropriate disease targets as therapeutics become available to treat inherited retinal disease.

Heterozygous and homozygous variants in STX1A cause a neurodevelopmental disorder with or without epilepsy

The neuronal SNARE complex drives synaptic vesicle exocytosis. Therefore, one of its core proteins syntaxin 1A (STX1A) has long been suspected to play a role in neurodevelopmental disorders. We assembled eight individuals harboring ultra rare variants in STX1A who present with a spectrum of intellectual disability, autism and epilepsy. Causative variants comprise a homozygous splice variant, three de novo missense variants and two inframe deletions of a single amino acid. We observed a phenotype mainly driven by epilepsy in the individuals with missense variants in contrast to intellectual disability and autistic behavior in individuals with single amino acid deletions and the splicing variant. In silico modeling of missense variants and single amino acid deletions show different impaired protein-protein interactions. We hypothesize the two phenotypic courses of affected individuals to be dependent on two different pathogenic mechanisms: (1) a weakened inhibitory STX1A-STXBP1 interaction due to missense variants results in an STX1A-related developmental epileptic encephalopathy and (2) a hampered SNARE complex formation due to inframe deletions causes an STX1A-related intellectual disability and autism phenotype. Our description of a STX1A-related neurodevelopmental disorder with or without epilepsy thus expands the group of rare diseases called SNAREopathies.

Mendelian inheritance revisited: dominance and recessiveness in medical genetics

Understanding the consequences of genotype for phenotype (which ranges from molecule-level effects to whole-organism traits) is at the core of genetic diagnostics in medicine. Many measures of the deleteriousness of individual alleles exist, but these have limitations for predicting the clinical consequences. Various mechanisms can protect the organism from the adverse effects of functional variants, especially when the variant is paired with a wild type allele. Understanding why some alleles are harmful in the heterozygous state - representing dominant inheritance - but others only with the biallelic presence of pathogenic variants - representing recessive inheritance - is particularly important when faced with the deluge of rare genetic alterations identified by high throughput DNA sequencing. Both awareness of the specific quantitative and/or qualitative effects of individual variants and the elucidation of allelic and non-allelic interactions are essential to optimize genetic diagnosis and counselling.

Case report: Novel compound heterozygosity for pathogenic variants in MED23 in a syndromic patient with postnatal microcephaly

Biallelic loss-of-function variants in MED23 cause a recessive syndromic intellectual disability condition with or without epilepsy (MRT18). Due to the small number of reported individuals, the clinical phenotype of the disorder has not been fully delineated yet, and the spectrum and frequency of neurologic features have not been fully characterized. Here, we report a 5-year-old girl with compound heterozygous for two additional MED23 variants. Besides global developmental delay, axial hypotonia and peripheral increased muscular tone, absent speech, and generalized tonic seizures, which fit well MRT18, the occurrence of postnatal progressive microcephaly has been here documented. A retrospective assessment of the previously reported clinical data for these subjects confirms the occurrence of postnatal progressive microcephaly as a previously unappreciated feature of the phenotype of MED23-related disorder.

A novel case of two siblings harbouring homozygous variant in the NEUROG1 gene with autism as an additional phenotype: a case report

INTRODUCTION

NEUROG1 gene is yet to be associated with a set of human phenotypes in the OMIM database. Three cases have previously been diagnosed with cranial dysinnervation due to biallelic variants in the NEUROG1 gene. This is the fourth and a novel report of a sibling pair harboring a homozygous variant in the NEUROG1 gene with autism as an additional phenotype. A brief review of the literature in conjunction with a genotype-phenotype correlation has been described. A potential hypothesis for the presence of the autistic phenotype in the present case has also been elucidated.

CASE PRESENTATION

A female aged 6 years and 9 months born to endogamous and phenotypically healthy parents was diagnosed with global developmental delay, autism spectrum disorder, hearing loss, corneal opacity and no eye blinking. Her MRI of the brain revealed mild peritrigonal white matter hyperintensity, and MRI and CT scan of the temporal bones showed abnormal cranial nerves. The proband's younger sister, aged 4-years, was similarly affected. Whole exome sequencing was performed in the proband, which revealed a novel homozygous, likely pathogenic, truncating frameshift variant, c.228_231dup (p.Thr78ProfsTer122) in exon 1 of the NEUROG1 gene (ENST00000314744.4). Segregation analysis by Sanger sequencing showed the proband and her younger sister to be homozygotes and their parents to be heterozygous carriers.

CONCLUSION

This is the fourth report across the globe with a variant identified in the NEUROG1 gene to be associated with cranial dysinnervation phenotype. An additional phenotype of autism in two female siblings was a novel observation. We provide a hypothetical framework which could explain the pleiotropic effect of a dysfunctional NEUROG1 protein leading to autism and posit it as a candidate for diagnosis of autism spectrum disorder with congenital cranial dysinnervation disorder.

Mutation of the proline P81 into a serine modifies the tumour suppressor function of the von Hippel-Lindau gene in the ccRCC

BACKGROUND

The von Hippel-Lindau disease is an autosomal dominant syndrome associated with tumour formation in various tissues, such as retina, central nervous system, kidney, and adrenal glands. VHL gene deletion or mutations support the development of various cancers. Unclassified VHL variants also referred as "of unknown significance" result from gene mutations that have an unknown or unclear effect on protein functions. The P81S mutation has been linked to low penetrance Type 1 disease but its pathogenic function was not clearly determined.

METHODS

We established a stable cell line expressing the pVHL213 (c.241C>T, P81S) mutant. Using biochemical and physiological approaches, we herein analysed pVHL folding, stability and function in the context of this VHL single missense mutation.

RESULTS

The P81S mutation mostly affects the non-canonical function of the pVHL protein. The cells expressing the pVHL213P81S acquire invasive properties in relation with modified architecture network.

CONCLUSION

We demonstrated the pathogenic role of this mutation in tumour development in vhl patients and confirm a medical follow up of family carrying the c.241C>T, P81S.

Ribosomal protein S6 kinase beta-1 gene variants cause hypertrophic cardiomyopathy

BACKGROUND

Hypertrophic cardiomyopathy (HCM) is a genetic heart muscle disease with preserved or increased ejection fraction in the absence of secondary causes. Mutations in the sarcomeric protein-encoding genes predominantly cause HCM. However, relatively little is known about the genetic impact of signalling proteins on HCM.

METHODS AND RESULTS

Here, using exome and targeted sequencing methods, we analysed two independent cohorts comprising 401 Indian patients with HCM and 3521 Indian controls. We identified novel variants in ribosomal protein S6 kinase beta-1 (RPS6KB1 or S6K1) gene in two unrelated Indian families as a potential candidate gene for HCM. The two unrelated HCM families had the same heterozygous missense S6K1 variant (p.G47W). In a replication association study, we identified two S6K1 heterozygotes variants (p.Q49K and p.Y62H) in the UK Biobank cardiomyopathy cohort (n=190) compared with matched controls (n=16 479). These variants are neither detected in region-specific controls nor in the human population genome data. Additionally, we observed an S6K1 variant (p.P445S) in an Arab patient with HCM. Functional consequences were evaluated using representative S6K1 mutated proteins compared with wild type in cellular models. The mutated proteins activated the S6K1 and hyperphosphorylated the rpS6 and ERK1/2 signalling cascades, suggesting a gain-of-function effect.

CONCLUSIONS

Our study demonstrates for the first time that the variants in the S6K1 gene are associated with HCM, and early detection of the S6K1 variant carriers can help to identify family members at risk and subsequent preventive measures. Further screening in patients with HCM with different ethnic populations will establish the specificity and frequency of S6K1 gene variants.

GWAS meta-analysis of intrahepatic cholestasis of pregnancy implicates multiple hepatic genes and regulatory elements

Intrahepatic cholestasis of pregnancy (ICP) is a pregnancy-specific liver disorder affecting 0.5-2% of pregnancies. The majority of cases present in the third trimester with pruritus, elevated serum bile acids and abnormal serum liver tests. ICP is associated with an increased risk of adverse outcomes, including spontaneous preterm birth and stillbirth. Whilst rare mutations affecting hepatobiliary transporters contribute to the aetiology of ICP, the role of common genetic variation in ICP has not been systematically characterised to date. Here, we perform genome-wide association studies (GWAS) and meta-analyses for ICP across three studies including 1138 cases and 153,642 controls. Eleven loci achieve genome-wide significance and have been further investigated and fine-mapped using functional genomics approaches. Our results pinpoint common sequence variation in liver-enriched genes and liver-specific cis-regulatory elements as contributing mechanisms to ICP susceptibility.

Mitochondrial bioenergetic is impaired in Monocarboxylate transporter 1 deficiency: a new clinical case and review of the literature

Background

Monocarboxylate transporter 1 (MCT1) deficiency has recently been described as a rare cause of recurrent ketosis, the result of impaired ketone utilization in extrahepatic tissues. To date, only six patients with this condition have been identified, and clinical and biochemical details remain incomplete.

Results

The present work reports a patient suffering from severe, recurrent episodes of metabolic acidosis and psychomotor delay, showing a pathogenic loss-of-function variation c.747_750del in homozygosity in SLC16A1 (which codes for MCT1). Persistent ketotic and lactic acidosis was accompanied by an abnormal excretion of organic acids related to redox balance disturbances. Together with an altered bioenergetic profile detected in patient-derived fibroblasts, this suggests possible mitochondrial dysfunction. Brain MRI revealed extensive, diffuse bilateral, symmetric signal alterations for the subcortical white matter and basal ganglia, together with corpus callosum agenesia.

Conclusions

These findings suggest that the clinical spectrum of MCT1 deficiency not only involves recurrent atacks of ketoacidosis, but may also cause lactic acidosis and neuromotor delay with a distinctive neuroimaging pattern including agenesis of corpus callosum and other brain signal alterations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13023-022-02389-4.

Novel biallelic variants expand the SLC5A6-related phenotypic spectrum

The sodium (Na+):multivitamin transporter (SMVT), encoded by SLC5A6, belongs to the sodium:solute symporter family and is required for the Na+-dependent uptake of biotin (vitamin B7), pantothenic acid (vitamin B5), the vitamin-like substance α-lipoic acid, and iodide. Compound heterozygous SLC5A6 variants have been reported in individuals with variable multisystemic disorder, including failure to thrive, developmental delay, seizures, cerebral palsy, brain atrophy, gastrointestinal problems, immunodeficiency, and/or osteopenia. We expand the phenotypic spectrum associated with biallelic SLC5A6 variants affecting function by reporting five individuals from three families with motor neuropathies. We identified the homozygous variant c.1285 A > G [p.(Ser429Gly)] in three affected siblings and a simplex patient and the maternally inherited c.280 C > T [p.(Arg94*)] variant and the paternally inherited c.485 A > G [p.(Tyr162Cys)] variant in the simplex patient of the third family. Both missense variants were predicted to affect function by in silico tools. 3D homology modeling of the human SMVT revealed 13 transmembrane helices (TMs) and Tyr162 and Ser429 to be located at the cytoplasmic facing region of TM4 and within TM11, respectively. The SLC5A6 missense variants p.(Tyr162Cys) and p.(Ser429Gly) did not affect plasma membrane localization of the ectopically expressed multivitamin transporter suggesting reduced but not abolished function, such as lower catalytic activity. Targeted therapeutic intervention yielded clinical improvement in four of the five patients. Early molecular diagnosis by exome sequencing is essential for timely replacement therapy in affected individuals.

ECNano: A cost-effective workflow for target enrichment sequencing and accurate variant calling on 4800 clinically significant genes using a single MinION flowcell

BACKGROUND

The application of long-read sequencing using the Oxford Nanopore Technologies (ONT) MinION sequencer is getting more diverse in the medical field. Having a high sequencing error of ONT and limited throughput from a single MinION flowcell, however, limits its applicability for accurate variant detection. Medical exome sequencing (MES) targets clinically significant exon regions, allowing rapid and comprehensive screening of pathogenic variants. By applying MES with MinION sequencing, the technology can achieve a more uniform capture of the target regions, shorter turnaround time, and lower sequencing cost per sample.

METHOD

We introduced a cost-effective optimized workflow, ECNano, comprising a wet-lab protocol and bioinformatics analysis, for accurate variant detection at 4800 clinically important genes and regions using a single MinION flowcell. The ECNano wet-lab protocol was optimized to perform long-read target enrichment and ONT library preparation to stably generate high-quality MES data with adequate coverage. The subsequent variant-calling workflow, Clair-ensemble, adopted a fast RNN-based variant caller, Clair, and was optimized for target enrichment data. To evaluate its performance and practicality, ECNano was tested on both reference DNA samples and patient samples.

RESULTS

ECNano achieved deep on-target depth of coverage (DoC) at average > 100× and > 98% uniformity using one MinION flowcell. For accurate ONT variant calling, the generated reads sufficiently covered 98.9% of pathogenic positions listed in ClinVar, with 98.96% having at least 30× DoC. ECNano obtained an average read length of 1000 bp. The long reads of ECNano also covered the adjacent splice sites well, with 98.5% of positions having ≥ 30× DoC. Clair-ensemble achieved > 99% recall and accuracy for SNV calling. The whole workflow from wet-lab protocol to variant detection was completed within three days.

CONCLUSION

We presented ECNano, an out-of-the-box workflow comprising (1) a wet-lab protocol for ONT target enrichment sequencing and (2) a downstream variant detection workflow, Clair-ensemble. The workflow is cost-effective, with a short turnaround time for high accuracy variant calling in 4800 clinically significant genes and regions using a single MinION flowcell. The long-read exon captured data has potential for further development, promoting the application of long-read sequencing in personalized disease treatment and risk prediction.

Consolidation of the clinical and genetic definition of a SOX4-related neurodevelopmental syndrome

BACKGROUND

A neurodevelopmental syndrome was recently reported in four patients with SOX4 heterozygous missense variants in the high-mobility-group (HMG) DNA-binding domain. The present study aimed to consolidate clinical and genetic knowledge of this syndrome.

METHODS

We newly identified 17 patients with SOX4 variants, predicted variant pathogenicity using in silico tests and in vitro functional assays and analysed the patients' phenotypes.

RESULTS

All variants were novel, distinct and heterozygous. Seven HMG-domain missense and five stop-gain variants were classified as pathogenic or likely pathogenic variant (L/PV) as they precluded SOX4 transcriptional activity in vitro. Five HMG-domain and non-HMG-domain missense variants were classified as of uncertain significance (VUS) due to negative results from functional tests. When known, inheritance was de novo or from a mosaic unaffected or non-mosaic affected parent for patients with L/PV, and from a non-mosaic asymptomatic or affected parent for patients with VUS. All patients had neurodevelopmental, neurological and dysmorphic features, and at least one cardiovascular, ophthalmological, musculoskeletal or other somatic anomaly. Patients with L/PV were overall more affected than patients with VUS. They resembled patients with other neurodevelopmental diseases, including the SOX11-related and Coffin-Siris (CSS) syndromes, but lacked the most specific features of CSS.

CONCLUSION

These findings consolidate evidence of a fairly non-specific neurodevelopmental syndrome due to SOX4 haploinsufficiency in neurogenesis and multiple other developmental processes.

Novel phenotype and genotype spectrum of WDR62 in two patients with associated primary autosomal recessive microcephaly

BACKGROUND

Microcephaly is a prominent feature of patients with primary autosomal recessive microcephaly 2 (MCPH2) caused by mutations in the WD Repeat Domain 62 (WDR62; OMIM: 613,583).

AIM

The study aimed to identify the underlying genetic factor(s) causing microcephaly in two patients in a consanguineous Iranian family.

METHODS

Two male patients (11 and 27 years old) were noticed due to microcephaly, neurodevelopmental delay, and occasional seizures. The younger patient (the proband) was subjected to paired-end whole-exome sequencing followed by Sanger sequencing to detect any underlying genetic factor.

RESULTS

Upon examination, both patients showed microcephaly as a prominent manifestation; they were under-weighted as well. The patients had a moderate gross motor impairment, severe cognitive disability and speech delay, increased deep tendon reflexes, flexible joint contractures, sensorineural hearing loss, and vertical nystagmus as a new ocular finding. The proband had more severe neurodevelopmental delay symptoms. The brain magnetic resonance imaging series revealed severe structural and cortical brain abnormalities in addition to hemiatrophy. Using Whole-exome Sequencing, a novel homozygous missense variant-NM_001083961.2; c.1598A > G: p.(His533Arg)-was identified in the WDR62. Subsequently, in silico analyses determined the possible impacts of the novel variant on the structure and function of WDR62 protein.

CONCLUSIONS

Herein, we identified a novel homozygous missense variant in the WDR62 in two patients with MCPH2. Vertical nystagmus and sensorineural hearing loss were detected as novel neurological findings. The present study expands the phenotype and genotype spectrum of MCPH2.

A novel deep intronic variant strongly associates with Alkaptonuria

Alkaptonuria is a rare autosomal recessive inherited disorder of tyrosine metabolism, which causes ochronosis, arthropathy, cardiac valvular calcification, and urolithiasis. The epidemiology of alkaptonuria in East Asia is not clear. In this study, patients diagnosed with alkaptonuria from January 2010 to June 2020 were reviewed. Their clinical and molecular features were further compared with those of patients from other countries. Three patients were found to have alkaptonuria. Mutation analyses of the homogentisate 1,2-dioxygenase gene (HGD) showed four novel variants c.16-2063 A > C, p.(Thr196Ile), p.(Gly344AspfsTer25), and p.(Gly362Arg) in six mutated alleles (83.3%). RNA sequencing revealed that c.16-2063 A > C activates a cryptic exon, causing protein truncation p.(Tyr5_Ile6insValTer17). A literature search identified another 6 patients with alkaptonuria in East Asia; including our cases, 13 of the 18 mutated alleles have not been reported elsewhere in the world. Alkaptonuria is rare in Taiwan and East Asia, with HGD variants being mostly novel and private.

Childhood-onset dystonia-causing KMT2B variants result in a distinctive genomic hypermethylation profile

Background

Dystonia is a clinically and genetically heterogeneous movement disorder characterized by sustained or intermittent muscle contractions causing abnormal, often repetitive, movements and/or postures. Heterozygous variants in lysine methyltransferase 2B (KMT2B), encoding a histone H3 methyltransferase, have been associated with a childhood-onset, progressive and complex form of dystonia (dystonia 28, DYT28). Since 2016, more than one hundred rare KMT2B variants have been reported, including frameshift, nonsense, splice site, missense and other in-frame changes, many having an uncertain clinical impact.

Results

We characterize the genome-wide peripheral blood DNA methylation profiles of a cohort of 18 patients with pathogenic and unclassified KMT2B variants. We resolve the “episignature” associated with KMT2B haploinsufficiency, proving that this approach is robust in diagnosing clinically unsolved cases, properly classifying them with respect to other partially overlapping dystonic phenotypes, other rare neurodevelopmental disorders and healthy controls. Notably, defective KMT2B function in DYT28 causes a non-random DNA hypermethylation across the genome, selectively involving promoters and other regulatory regions positively controlling gene expression.

Conclusions

We demonstrate a distinctive DNA hypermethylation pattern associated with DYT28, provide an epigenetic signature for this disorder enabling accurate diagnosis and reclassification of ambiguous genetic findings and suggest potential therapeutic approaches.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13148-021-01145-y.

Broadening INPP5E phenotypic spectrum: detection of rare variants in syndromic and non-syndromic IRD

Pathogenic variants in INPP5E cause Joubert syndrome (JBTS), a ciliopathy with retinal involvement. However, despite sporadic cases in large cohort sequencing studies, a clear association with non-syndromic inherited retinal degenerations (IRDs) has not been made. We validate this association by reporting 16 non-syndromic IRD patients from ten families with bi-allelic mutations in INPP5E. Additional two patients showed early onset IRD with limited JBTS features. Detailed phenotypic description for all probands is presented. We report 14 rare INPP5E variants, 12 of which have not been reported in previous studies. We present tertiary protein modeling and analyze all INPP5E variants for deleteriousness and phenotypic correlation. We observe that the combined impact of INPP5E variants in JBTS and non-syndromic IRD patients does not reveal a clear genotype-phenotype correlation, suggesting the involvement of genetic modifiers. Our study cements the wide phenotypic spectrum of INPP5E disease, adding proof that sequence defects in this gene can lead to early-onset non-syndromic IRD.

Variants in GCNA, X-linked germ-cell genome integrity gene, identified in men with primary spermatogenic failure

Male infertility impacts millions of couples yet, the etiology of primary infertility remains largely unknown. A critical element of successful spermatogenesis is maintenance of genome integrity. Here, we present a genomic study of spermatogenic failure (SPGF). Our initial analysis (n = 176) did not reveal known gene-candidates but identified a potentially significant single-nucleotide variant (SNV) in X-linked germ-cell nuclear antigen (GCNA). Together with a larger follow-up study (n = 2049), 7 likely clinically relevant GCNA variants were identified. GCNA is critical for genome integrity in male meiosis and knockout models exhibit impaired spermatogenesis and infertility. Single-cell RNA-seq and immunohistochemistry confirm human GCNA expression from spermatogonia to elongated spermatids. Five identified SNVs were located in key functional regions, including N-terminal SUMO-interacting motif and C-terminal Spartan-like protease domain. Notably, variant p.Ala115ProfsTer7 results in an early frameshift, while Spartan-like domain missense variants p.Ser659Trp and p.Arg664Cys change conserved residues, likely affecting 3D structure. For variants within GCNA's intrinsically disordered region, we performed computational modeling for consensus motifs. Two SNVs were predicted to impact the structure of these consensus motifs. All identified variants have an extremely low minor allele frequency in the general population and 6 of 7 were not detected in > 5000 biological fathers. Considering evidence from animal models, germ-cell-specific expression, 3D modeling, and computational predictions for SNVs, we propose that identified GCNA variants disrupt structure and function of the respective protein domains, ultimately arresting germ-cell division. To our knowledge, this is the first study implicating GCNA, a key genome integrity factor, in human male infertility.

Rare deleterious mutations of HNRNP genes result in shared neurodevelopmental disorders

BACKGROUND

With the increasing number of genomic sequencing studies, hundreds of genes have been implicated in neurodevelopmental disorders (NDDs). The rate of gene discovery far outpaces our understanding of genotype-phenotype correlations, with clinical characterization remaining a bottleneck for understanding NDDs. Most disease-associated Mendelian genes are members of gene families, and we hypothesize that those with related molecular function share clinical presentations.

METHODS

We tested our hypothesis by considering gene families that have multiple members with an enrichment of de novo variants among NDDs, as determined by previous meta-analyses. One of these gene families is the heterogeneous nuclear ribonucleoproteins (hnRNPs), which has 33 members, five of which have been recently identified as NDD genes (HNRNPK, HNRNPU, HNRNPH1, HNRNPH2, and HNRNPR) and two of which have significant enrichment in our previous meta-analysis of probands with NDDs (HNRNPU and SYNCRIP). Utilizing protein homology, mutation analyses, gene expression analyses, and phenotypic characterization, we provide evidence for variation in 12 HNRNP genes as candidates for NDDs. Seven are potentially novel while the remaining genes in the family likely do not significantly contribute to NDD risk.

RESULTS

We report 119 new NDD cases (64 de novo variants) through sequencing and international collaborations and combined with published clinical case reports. We consider 235 cases with gene-disruptive single-nucleotide variants or indels and 15 cases with small copy number variants. Three hnRNP-encoding genes reach nominal or exome-wide significance for de novo variant enrichment, while nine are candidates for pathogenic mutations. Comparison of HNRNP gene expression shows a pattern consistent with a role in cerebral cortical development with enriched expression among radial glial progenitors. Clinical assessment of probands (n = 188-221) expands the phenotypes associated with HNRNP rare variants, and phenotypes associated with variation in the HNRNP genes distinguishes them as a subgroup of NDDs.

CONCLUSIONS

Overall, our novel approach of exploiting gene families in NDDs identifies new HNRNP-related disorders, expands the phenotypes of known HNRNP-related disorders, strongly implicates disruption of the hnRNPs as a whole in NDDs, and supports that NDD subtypes likely have shared molecular pathogenesis. To date, this is the first study to identify novel genetic disorders based on the presence of disorders in related genes. We also perform the first phenotypic analyses focusing on related genes. Finally, we show that radial glial expression of these genes is likely critical during neurodevelopment. This is important for diagnostics, as well as developing strategies to best study these genes for the development of therapeutics.

Pathogenic variants in CDH11 impair cell adhesion and cause Teebi hypertelorism syndrome

Teebi hypertelorism syndrome (THS; OMIM 145420) is a rare craniofacial disorder characterized by hypertelorism, prominent forehead, short nose with broad or depressed nasal root. Some cases of THS have been attributed to SPECC1L variants. Homozygous variants in CDH11 truncating the transmembrane and intracellular domains have been implicated in Elsahy-Waters syndrome (EWS; OMIM 211380) with hypertelorism. We report THS due to CDH11 heterozygous missense variants on 19 subjects from 9 families. All affected residues in the extracellular region of Cadherin-11 (CHD11) are highly conserved across vertebrate species and classical cadherins. Six of the variants that cluster around the EC2-EC3 and EC3-EC4 linker regions are predicted to affect Ca²⁺ binding that is required for cadherin stability. Two of the additional variants [c.164G > C, p.(Trp55Ser) and c.418G > A, p.(Glu140Lys)] are also notable as they are predicted to directly affect trans-homodimer formation. Immunohistochemical study demonstrates that CDH11 is strongly expressed in human facial mesenchyme. Using multiple functional assays, we show that five variants from the EC1, EC2-EC3 linker, and EC3 regions significantly reduced the cell-substrate trans adhesion activity and one variant from EC3-EC4 linker results in changes in cell morphology, focal adhesion, and migration, suggesting dominant negative effect. Characteristic features in this cohort included depressed nasal root, cardiac and umbilical defects. These features distinguished this phenotype from that seen in SPECC1L-related hypertelorism syndrome and CDH11-related EWS. Our results demonstrate heterozygous variants in CDH11, which decrease cell-cell adhesion and increase cell migratory behavior, cause a form of THS, as termed CDH11-related THS.

Association of SLC32A1 Missense Variants With Genetic Epilepsy With Febrile Seizures Plus

OBJECTIVE

To identify the causative gene in a large unsolved family with genetic epilepsy with febrile seizures plus (GEFS+), we sequenced the genomes of family members, and then determined the contribution of the identified gene to the pathogenicity of epilepsies by examining sequencing data from 2,772 additional patients.

METHODS

We performed whole genome sequencing of 3 members of a GEFS+ family. Subsequently, whole exome sequencing data from 1,165 patients with epilepsy from the Epi4K dataset and 1,329 Australian patients with epilepsy from the Epi25 dataset were interrogated. Targeted resequencing was performed on 278 patients with febrile seizures or GEFS+ phenotypes. Variants were validated and familial segregation examined by Sanger sequencing.

RESULTS

Eight previously unreported missense variants were identified in SLC32A1, coding for the vesicular inhibitory amino acid cotransporter VGAT. Two variants cosegregated with the phenotype in 2 large GEFS+ families containing 8 and 10 affected individuals, respectively. Six further variants were identified in smaller families with GEFS+ or idiopathic generalized epilepsy (IGE).

CONCLUSION

Missense variants in SLC32A1 cause GEFS+ and IGE. These variants are predicted to alter γ-aminobutyric acid (GABA) transport into synaptic vesicles, leading to altered neuronal inhibition. Examination of further epilepsy cohorts will determine the full genotype-phenotype spectrum associated with SLC32A1 variants.

Identification of novel, clonally stable, somatic mutations targeting transcription factors PAX5 and NKX2-3, the epigenetic regulator LRIF1, and BRAF in a case of atypical B-cell chronic lymphocytic leukemia harboring a t(14;18)(q32;q21)

Diagnosis of B-cell chronic lymphocytic leukemia (B-CLL) is usually straightforward, involving clinical, immunophenotypic (Matutes score), and (immuno)genetic analyses (to refine patient prognosis for treatment). CLL cases with atypical presentation (e.g., Matutes ≤ 3) are also encountered, and for these diseases, biology and prognostic impact are less clear. Here we report the genomic characterization of a case of atypical B-CLL in a 70-yr-old male patient; B-CLL cells showed a Matutes score of 3, chromosomal translocation t(14;18)(q32;q21) (BCL2/IGH), mutated IGHV, deletion 17p, and mutations in BCL2, NOTCH1 (subclonal), and TP53 (subclonal). Quite strikingly, a novel PAX5 mutation that was predicted to be loss of function was also seen. Exome sequencing identified, in addition, a potentially actionable BRAF mutation, together with novel somatic mutations affecting the homeobox transcription factor NKX2-3, known to control B-lymphocyte development and homing, and the epigenetic regulator LRIF1, which is implicated in chromatin compaction and gene silencing. Neither NKX2-3 nor LRIF1 mutations, predicted to be loss of function, have previously been reported in B-CLL. Sequencing confirmed the presence of these mutations together with BCL2, NOTCH1, and BRAF mutations, with the t(14;18)(q32;q21) translocation, in the initial diagnostic sample obtained 12 yr prior. This is suggestive of a role for these novel mutations in B-CLL initiation and stable clonal evolution, including upon treatment withdrawal. This case extends the spectrum of atypical B-CLL with t(14;18)(q32;q21) and highlights the value of more global precision genomics for patient follow-up and treatment in these patients.

Normal activated partial thromboplastin time in Chinese patients with mild hemophilia B

OBJECTIVES

Hemophilia B (HB, OMIM: 300746) is one of the most common bleeding disorders with an X-linked recessive inheritance pattern, caused by the deficiency of coagulation factor IX (FIX). FIX is encoded by the F9 gene located on Xq27.1. Diagnosis of HB is primarily suspected by prolonged activated partial thromboplastin time (APTT), decreased FIX activity (FIX:C) or genetic test of the F9 gene. We herein described a Chinese family with patients of mild HB.

METHODS

Sanger sequencing of the F9 gene was applied to identify mutation. Coagulation tests were performed.

RESULTS

The proband was a 5-year-old boy. He suffered prolonged bleeding after tonsillectomy recently and circumcision last year as well. His grandfather experienced prolonged bleeding after gastric surgery. Both patients showed normal APTT, though they had significantly decreased FIX:C. Sanger sequencing of the F9 gene revealed a novel hemizygous F9 c.639C > A (p.Asn213Lys) missense mutation in both patients. The proband's mother carried heterozygous mutation. This mutation was located in the activation peptide domain of FIX.

CONCLUSION

In conclusion, we confirmed that APTT could be normal in mild HB patients. Highly sensitive APTT for mild HB and molecular genetic test could confirm the diagnosis of mild HB.

A novel missense variant in the EML1 gene associated with bilateral ribbon-like subcortical heterotopia leads to ciliary defects

Heterotopia is a brain malformation caused by a failed migration of cortical neurons during development. Clinical symptoms of heterotopia vary in severity of intellectual disability and may be associated with epileptic disorders. Abnormal neuronal migration is known to be associated with mutations in the doublecortin gene (DCX), the platelet-activating factor acetylhydrolase gene (PAFAH1B1), or tubulin alpha-1A gene (TUBA1A). Recently, a new gene encoding echinoderm microtubule-associated protein-like 1 (EML1) was reported to cause a particular form of subcortical heterotopia, the ribbon-like subcortical heterotopia (RSH). EML1 mutations are inherited in an autosomal recessive manner. Only six unrelated EML1-associated heterotopia-affected families were reported so far. The EML1 protein is a member of the microtubule-associated proteins family, playing an important role in microtubule assembly and stabilization as well as in mitotic spindle formation in interphase. Herein, we present a novel homozygous missense variant in EML1 (NM_004434.2: c.692G>A, NP_004425.2: p.Gly231Asp) identified in a male RSH-affected patient. Our clinical and molecular findings confirm the genotype-phenotype associations of EML1 mutations and RSH. Analyses of patient-derived fibroblasts showed the significantly reduced length of primary cilia. In addition, our results presented, that the mutated EML1 protein did not change binding capacities with tubulin. The data described herein will expand the mutation spectrum of the EML1 gene and provide further insight into molecular and cellular bases of the pathogenic mechanisms underlying RSH.

MobiDetails: online DNA variants interpretation

MobiDetails is an expert tool, online application which gathers useful data for the interpretation of DNA variants in the context of molecular diagnosis. It brings together in a single tool many sources of data, such as population genetics, various kinds of predictors, Human Genome Variation Society (HGVS) nomenclatures, curated databases, and access to various annotations. Accurate interpretation of DNA variants is crucial and can impact the patient care or have familial outcomes (prenatal diagnosis). Its importance will increase in the coming years with the expansion of the personalized medicine. MobiDetails is specifically designed to help with this task. Exonic or intronic substitutions and small insertions/deletions related to more than 18,000 human genes are easily submitted and annotated in real-time. It is a responsive website that can be accessed using mobiles or tablets during medical staff meetings. MobiDetails is based on publicly available resources, does not include any specific data on patients or phenotypes, and is freely available for academic use at https://mobidetails.iurc.montp.inserm.fr/MD/ .

De novo CLTC variants are associated with a variable phenotype from mild to severe intellectual disability, microcephaly, hypoplasia of the corpus callosum, and epilepsy

PURPOSE

To delineate the genotype-phenotype correlation in individuals with likely pathogenic variants in the CLTC gene.

METHODS

We describe 13 individuals with de novo CLTC variants. Causality of variants was determined by using the tolerance landscape of CLTC and computer-assisted molecular modeling where applicable. Phenotypic abnormalities observed in the individuals identified with missense and in-frame variants were compared with those with nonsense or frameshift variants in CLTC.

RESULTS

All de novo variants were judged to be causal. Combining our data with that of 14 previously reported affected individuals (n = 27), all had intellectual disability (ID), ranging from mild to moderate/severe, with or without additional neurologic, behavioral, craniofacial, ophthalmologic, and gastrointestinal features. Microcephaly, hypoplasia of the corpus callosum, and epilepsy were more frequently observed in individuals with missense and in-frame variants than in those with nonsense and frameshift variants. However, this difference was not significant.

CONCLUSIONS

The wide phenotypic variability associated with likely pathogenic CLTC variants seems to be associated with allelic heterogeneity. The detailed clinical characterization of a larger cohort of individuals with pathogenic CLTC variants is warranted to support the hypothesis that missense and in-frame variants exert a dominant-negative effect, whereas the nonsense and frameshift variants would result in haploinsufficiency.