Loss-of-function of the Zinc Finger Homeobox 4 (ZFHX4) gene underlies a neurodevelopmental disorder.

8q21.11 microdeletions encompassing the gene encoding transcription factor ZFHX4, have previously been associated by us with a syndromic form of intellectual disability, hypotonia, decreased balance and hearing loss. Here, we report on 57 individuals, 52 probands and 5 affected family members, with protein truncating variants (n=36), (micro)deletions (n=20) or an inversion (n=1) affecting ZFHX4 with variable developmental delay and intellectual disability, distinctive facial characteristics, morphological abnormalities of the central nervous system, behavioral alterations, short stature, hypotonia, and occasionally cleft palate and anterior segment dysgenesis. The phenotypes associated with 8q21.11 microdeletions and ZFHX4 intragenic loss-of-function variants largely overlap, identifying ZFHX4 as the main driver for the microdeletion syndrome, although leukocyte-derived DNA shows a mild common methylation profile for (micro)deletions only. We identify ZFHX4 as a transcription factor that is increasingly expressed during human brain development and neuronal differentiation. Furthermore, ZFHX4 interacting factors identified via IP-MS in neural progenitor cells, suggest an important role for ZFHX4 in cellular and developmental pathways, especially during histone modifications, cytosolic transport and development. Additionally, using CUT&RUN, we observed that ZFHX4 binds with the promoter regions of genes with crucial roles in embryonic, neuron and axon development. Since loss-of-function variants in ZFHX4 are found with consistent dysmorphic facial features, we investigated whether the disruption of zfhx4 causes craniofacial abnormalities in zebrafish. First-generation (F0) zfhx4 crispant zebrafish, (mosaic) mutant for zfhx4 loss-of-function variants, have significantly shorter Meckel's cartilages and smaller ethmoid plates compared to control zebrafish. Furthermore, behavioral assays show a decreased movement frequency in the zfhx4 crispant zebrafish in comparison with control zebrafish larvae. Although further research is needed, our in vivo work suggests a role for zfhx4 in facial skeleton patterning, palatal development and behavior.

Comparative 3D genome analysis between neural retina and retinal pigment epithelium reveals differential cis-regulatory interactions at retinal disease loci.

BACKGROUND: Vision depends on the interplay between photoreceptor cells of the neural retina and the underlying retinal pigment epithelium (RPE). Most genes involved in inherited retinal diseases display specific spatiotemporal expression within these interconnected retinal components through the local recruitment of cis-regulatory elements (CREs) in 3D nuclear space. RESULTS: To understand the role of differential chromatin architecture in establishing tissue-specific expression at inherited retinal disease loci, we mapped genome-wide chromatin interactions using in situ Hi-C and H3K4me3 HiChIP on neural retina and RPE/choroid from human adult donor eyes. We observed chromatin looping between active promoters and 32,425 and 8060 candidate CREs in the neural retina and RPE/choroid, respectively. A comparative 3D genome analysis between these two retinal tissues revealed that 56% of 290 known inherited retinal disease genes were marked by differential chromatin interactions. One of these was ABCA4, which is implicated in the most common autosomal recessive inherited retinal disease. We zoomed in on retina- and RPE-specific cis-regulatory interactions at the ABCA4 locus using high-resolution UMI-4C. Integration with bulk and single-cell epigenomic datasets and in vivo enhancer assays in zebrafish revealed tissue-specific CREs interacting with ABCA4. CONCLUSIONS: Through comparative 3D genome mapping, based on genome-wide, promoter-centric, and locus-specific assays of human neural retina and RPE, we have shown that gene regulation at key inherited retinal disease loci is likely mediated by tissue-specific chromatin interactions. These findings do not only provide insight into tissue-specific regulatory landscapes at retinal disease loci, but also delineate the search space for non-coding genomic variation underlying unsolved inherited retinal diseases.

Haploinsufficiency of ZFHX3, encoding a key player in neuronal development, causes syndromic intellectual disability.

Neurodevelopmental disorders (NDDs) result from impaired development and functioning of the brain. Here, we identify loss-of-function (LoF) variation in ZFHX3 as a cause for syndromic intellectual disability (ID). ZFHX3 is a zinc-finger homeodomain transcription factor involved in various biological processes, including cell differentiation and tumorigenesis. We describe 42 individuals with protein-truncating variants (PTVs) or (partial) deletions of ZFHX3, exhibiting variable intellectual disability and autism spectrum disorder, recurrent facial features, relative short stature, brachydactyly, and, rarely, cleft palate. ZFHX3 LoF associates with a specific methylation profile in whole blood extracted DNA. Nuclear abundance of ZFHX3 increases during human brain development and neuronal differentiation. ZFHX3 was found to interact with the chromatin remodeling BRG1/Brm-associated factor complex and the cleavage and polyadenylation complex, suggesting a function in chromatin remodeling and mRNA processing. Furthermore, ChIP-seq for ZFHX3 revealed that it predominantly binds promoters of genes involved in nervous system development. We conclude that loss-of-function variants in ZFHX3 are a cause of syndromic ID associating with a specific DNA methylation profile.

A distant global control region is essential for normal expression of anterior HOXA genes during mouse and human craniofacial development.

Craniofacial abnormalities account for approximately one third of birth defects. The regulatory programs that build the face require precisely controlled spatiotemporal gene expression, achieved through tissue-specific enhancers. Clusters of coactivated enhancers and their target genes, known as superenhancers, are important in determining cell identity but have been largely unexplored in development. In this study we identified superenhancer regions unique to human embryonic craniofacial tissue. To demonstrate the importance of such regions in craniofacial development and disease, we focused on an ~600 kb noncoding region located between NPVF and NFE2L3. We identified long range interactions with this region in both human and mouse embryonic craniofacial tissue with the anterior portion of the HOXA gene cluster. Mice lacking this superenhancer exhibit perinatal lethality, and present with highly penetrant skull defects and orofacial clefts phenocopying Hoxa2-/- mice. Moreover, we identified two cases of de novo copy number changes of the superenhancer in humans both with severe craniofacial abnormalities. This evidence suggests we have identified a critical noncoding locus control region that specifically regulates anterior HOXA genes and copy number changes are pathogenic in human patients.

A novel neurodevelopmental syndrome caused by loss-of-function of the Zinc Finger Homeobox 3 (ZFHX3) gene.

Neurodevelopmental disorders (NDDs) result from impaired development and functioning of the brain. Here, we identify loss-of-function variation in ZFHX3 as a novel cause for syndromic intellectual disability (ID). ZFHX3, previously known as ATBF1, is a zinc-finger homeodomain transcription factor involved in multiple biological processes including cell differentiation and tumorigenesis. Through international collaboration, we collected clinical and morphometric data (Face2Gene) of 41 individuals with protein truncating variants (PTVs) or (partial) deletions of ZFHX3 . We used data mining, RNA and protein analysis to identify the subcellular localization and spatiotemporal expression of ZFHX3 in multiple in vitro models. We identified the DNA targets of ZFHX3 using ChIP seq. Immunoprecipitation followed by mass spectrometry indicated potential binding partners of endogenous ZFHX3 in neural stem cells that were subsequently confirmed by reversed co-immunoprecipitation and western blot. We evaluated a DNA methylation profile associated with ZFHX3 haploinsufficiency using DNA methylation analysis on whole blood extracted DNA of six individuals with ZFHX3 PTVs and four with a (partial) deletion of ZFHX3 . A reversed genetic approach characterized the ZFHX3 orthologue in Drosophila melanogaster . Loss-of-function variation of ZFHX3 consistently associates with (mild) ID and/or behavioural problems, postnatal growth retardation, feeding difficulties, and recognizable facial characteristics, including the rare occurrence of cleft palate. Nuclear abundance of ZFHX3 increases during human brain development and neuronal differentiation in neural stem cells and SH-SY5Y cells, ZFHX3 interacts with the chromatin remodelling BRG1/Brm-associated factor complex and the cleavage and polyadenylation complex. In line with a role for chromatin remodelling, ZFHX3 haploinsufficiency associates with a specific DNA methylation profile in leukocyte-derived DNA. The target genes of ZFHX3 are implicated in neuron and axon development. In Drosophila melanogaster , z fh2, considered to be the ZFHX3 orthologue, is expressed in the third instar larval brain. Ubiquitous and neuron-specific knockdown of zfh2 results in adult lethality underscoring a key role for zfh2 in development and neurodevelopment. Interestingly, ectopic expression of zfh2 as well as ZFHX3 in the developing wing disc results in a thoracic cleft phenotype. Collectively, our data shows that loss-of-function variants in ZFHX3 are a cause of syndromic ID, that associates with a specific DNA methylation profile. Furthermore, we show that ZFHX3 participates in chromatin remodelling and mRNA processing.

Transcriptional and functional consequences of alterations to MEF2C and its topological organization in neuronal models.

Point mutations and structural variants that directly disrupt the coding sequence of MEF2C have been associated with a spectrum of neurodevelopmental disorders (NDDs). However, the impact of MEF2C haploinsufficiency on neurodevelopmental pathways and synaptic processes is not well understood, nor are the complex mechanisms that govern its regulation. To explore the functional changes associated with structural variants that alter MEF2C expression and/or regulation, we generated an allelic series of 204 isogenic human induced pluripotent stem cell (hiPSC)-derived neural stem cells and glutamatergic induced neurons. These neuronal models harbored CRISPR-engineered mutations that involved direct deletion of MEF2C or deletion of the boundary points for topologically associating domains (TADs) and chromatin loops encompassing MEF2C. Systematic profiling of mutation-specific alterations, contrasted to unedited controls that were exposed to the same guide RNAs for each edit, revealed that deletion of MEF2C caused differential expression of genes associated with neurodevelopmental pathways and synaptic function. We also discovered significant reduction in synaptic activity measured by multielectrode arrays (MEAs) in neuronal cells. By contrast, we observed robust buffering against MEF2C regulatory disruption following deletion of a distal 5q14.3 TAD and loop boundary, whereas homozygous loss of a proximal loop boundary resulted in down-regulation of MEF2C expression and reduced electrophysiological activity on MEA that was comparable to direct gene disruption. Collectively, these studies highlight the considerable functional impact of MEF2C deletion in neuronal cells and systematically characterize the complex interactions that challenge a priori predictions of regulatory consequences from structural variants that disrupt three-dimensional genome organization.

Multi-omics approach dissects cis-regulatory mechanisms underlying North Carolina macular dystrophy, a retinal enhanceropathy.

North Carolina macular dystrophy (NCMD) is a rare autosomal-dominant disease affecting macular development. The disease is caused by non-coding single-nucleotide variants (SNVs) in two hotspot regions near PRDM13 and by duplications in two distinct chromosomal loci, overlapping DNase I hypersensitive sites near either PRDM13 or IRX1. To unravel the mechanisms by which these variants cause disease, we first established a genome-wide multi-omics retinal database, RegRet. Integration of UMI-4C profiles we generated on adult human retina then allowed fine-mapping of the interactions of the PRDM13 and IRX1 promoters and the identification of eighteen candidate cis-regulatory elements (cCREs), the activity of which was investigated by luciferase and Xenopus enhancer assays. Next, luciferase assays showed that the non-coding SNVs located in the two hotspot regions of PRDM13 affect cCRE activity, including two NCMD-associated non-coding SNVs that we identified herein. Interestingly, the cCRE containing one of these SNVs was shown to interact with the PRDM13 promoter, demonstrated in vivo activity in Xenopus, and is active at the developmental stage when progenitor cells of the central retina exit mitosis, suggesting that this region is a PRDM13 enhancer. Finally, mining of single-cell transcriptional data of embryonic and adult retina revealed the highest expression of PRDM13 and IRX1 when amacrine cells start to synapse with retinal ganglion cells, supporting the hypothesis that altered PRDM13 or IRX1 expression impairs interactions between these cells during retinogenesis. Overall, this study provides insight into the cis-regulatory mechanisms of NCMD and supports that this condition is a retinal enhanceropathy.

HDAC9 structural variants disrupting TWIST1 transcriptional regulation lead to craniofacial and limb malformations.

Structural variants (SVs) can affect protein-coding sequences as well as gene regulatory elements. However, SVs disrupting protein-coding sequences that also function as cis-regulatory elements remain largely uncharacterized. Here, we show that craniosynostosis patients with SVs containing the histone deacetylase 9 (HDAC9) protein-coding sequence are associated with disruption of TWIST1 regulatory elements that reside within the HDAC9 sequence. Based on SVs within the HDAC9-TWIST1 locus, we defined the 3'-HDAC9 sequence as a critical TWIST1 regulatory region, encompassing craniofacial TWIST1 enhancers and CTCF sites. Deletions of either Twist1 enhancers (eTw5-7Δ/Δ) or CTCF site (CTCF-5Δ/Δ) within the Hdac9 protein-coding sequence led to decreased Twist1 expression and altered anterior/posterior limb expression patterns of SHH pathway genes. This decreased Twist1 expression results in a smaller sized and asymmetric skull and polydactyly that resembles Twist1+/- mouse phenotype. Chromatin conformation analysis revealed that the Twist1 promoter interacts with Hdac9 sequences that encompass Twist1 enhancers and a CTCF site, and that interactions depended on the presence of both regulatory regions. Finally, a large inversion of the entire Hdac9 sequence (Hdac9 INV/+) in mice that does not disrupt Hdac9 expression but repositions Twist1 regulatory elements showed decreased Twist1 expression and led to a craniosynostosis-like phenotype and polydactyly. Thus, our study elucidates essential components of TWIST1 transcriptional machinery that reside within the HDAC9 sequence. It suggests that SVs encompassing protein-coding sequences could lead to a phenotype that is not attributed to its protein function but rather to a disruption of the transcriptional regulation of a nearby gene.

Expanding the Phenotype of B3GALNT2-Related Disorders.

Dystroglycanopathies are a group of congenital muscular dystrophies (CMDs) that include a broad phenotypic spectrum ranging from late-onset limb-girdle muscular dystrophy to severe muscle-eye-brain disease, Walker-Warburg syndrome, and Fukuyama congenital muscular dystrophy. In addition to clinical heterogeneity, CMDs are characterized by genetic heterogeneity. To date, 18 genes have been associated with CMDs. One of them is B3GALNT2, which encodes the ß-1,3-N-acetylgalactosaminyltransferase 2 that glycosylates α-dystroglycan. In this study, using exome sequencing, we identify a homozygous frameshift variant in B3GALNT2 due to a mixed uniparental disomy of chromosome 1 in a 7-year-old girl with global developmental delay, severely delayed active language development, and autism spectrum disorder but without any symptoms of muscular dystrophy. In addition to this case, we also provide an overview of all previously reported cases, further expanding the phenotypic spectrum.

MYT1L-associated neurodevelopmental disorder: description of 40 new cases and literature review of clinical and molecular aspects.

Pathogenic variants of the myelin transcription factor-1 like (MYT1L) gene include heterozygous missense, truncating variants and 2p25.3 microdeletions and cause a syndromic neurodevelopmental disorder (OMIM#616,521). Despite enrichment in de novo mutations in several developmental disorders and autism studies, the data on clinical characteristics and genotype-phenotype correlations are scarce, with only 22 patients with single nucleotide pathogenic variants reported. We aimed to further characterize this disorder at both the clinical and molecular levels by gathering a large series of patients with MYT1L-associated neurodevelopmental disorder. We collected genetic information on 40 unreported patients with likely pathogenic/pathogenic MYT1L variants and performed a comprehensive review of published data (total = 62 patients). We confirm that the main phenotypic features of the MYT1L-related disorder are developmental delay with language delay (95%), intellectual disability (ID, 70%), overweight or obesity (58%), behavioral disorders (98%) and epilepsy (23%). We highlight novel clinical characteristics, such as learning disabilities without ID (30%) and feeding difficulties during infancy (18%). We further describe the varied dysmorphic features (67%) and present the changes in weight over time of 27 patients. We show that patients harboring highly clustered missense variants in the 2-3-ZNF domains are not clinically distinguishable from patients with truncating variants. We provide an updated overview of clinical and genetic data of the MYT1L-associated neurodevelopmental disorder, hence improving diagnosis and clinical management of these patients.

A Reassessment of Copy Number Variations in Congenital Heart Defects: Picturing the Whole Genome.

Copy number variations (CNVs) can modulate phenotypes by affecting protein-coding sequences directly or through interference of gene expression. Recent studies in cancer and limb defects pinpointed the relevance of non-coding gene regulatory elements such as long non-coding RNAs (lncRNAs) and topologically associated domain (TAD)-related gene-enhancer interactions. The contribution of such non-coding elements is largely unexplored in congenital heart defects (CHD). We performed a retrospective analysis of CNVs reported in a cohort of 270 CHD patients. We reviewed the diagnostic yield of pathogenic CNVs, and performed a comprehensive reassessment of 138 CNVs of unknown significance (CNV-US), evaluating protein-coding genes, lncRNA genes, and potential interferences with TAD-related gene-enhancer interactions. Fifty-two of the 138 CNV-US may relate to CHD, revealing three candidate CHD regions, 19 candidate CHD genes, 80 lncRNA genes of interest, and six potentially CHD-related TAD interferences. Our study thus indicates a potential relevance of non-coding gene regulatory elements in CNV-related CHD pathogenesis. Shortcomings in our current knowledge on genomic variation call for continuous reporting of CNV-US in international databases, careful patient counseling, and additional functional studies to confirm these preliminary findings.

Expanding the molecular spectrum and the neurological phenotype related to CAMTA1 variants.

The CAMTA1-associated phenotype was initially defined in patients with intragenic deletions and duplications who showed nonprogressive congenital ataxia, with or without intellectual disability. Here, we describe 10 individuals with CAMTA1 variants: nine previously unreported (likely) pathogenic variants comprising one missense, four frameshift and four nonsense variants, and one missense variant of unknown significance. Six patients were diagnosed following whole exome sequencing and four individuals with exome-based targeted panel analysis. Most of them present with developmental delay, manifesting in speech and motor delay. Other frequent findings are hypotonia, cognitive impairment, cerebellar dysfunction, oculomotor abnormalities, and behavioral problems. Feeding problems occur more frequently than previously observed. In addition, we present a systematic review of 19 previously published individuals with causal variants, including copy number, truncating, and missense variants. We note a tendency of more severe cognitive impairment and recurrent dysmorphic features in individuals with a copy number variant. Pathogenic variants are predominantly observed in and near the N- and C- terminal functional domains. Clinical heterogeneity is observed, but 3'-terminal variants seem to associate with less pronounced cerebellar dysfunction.

Interpreting the impact of noncoding structural variation in neurodevelopmental disorders.

The emergence of novel sequencing technologies has greatly improved the identification of structural variation, revealing that a human genome harbors tens of thousands of structural variants (SVs). Since these SVs primarily impact noncoding DNA sequences, the next challenge is one of interpretation, not least to improve our understanding of human disease etiology. However, this task is severely complicated by the intricacy of the gene regulatory landscapes embedded within these noncoding regions, their incomplete annotation, as well as their dependence on the three-dimensional (3D) conformation of the genome. Also in the context of neurodevelopmental disorders (NDDs), reports of putatively causal, noncoding SVs are accumulating and understanding their impact on transcriptional regulation is presenting itself as the next step toward improved genetic diagnosis.

Phenotypic spectrum of the RBM10-mediated intellectual disability and congenital malformation syndrome beyond classic TARP syndrome features.

Pathogenic variants in the RBM10 gene cause a rare X-linked disorder described as TARP (Talipes equinovarus, Atrial septal defect, Robin sequence, and Persistent left vena cava superior) syndrome. We report two novel patients with truncating RBM10 variants in view of the literature, presenting a total of 26 patients from 15 unrelated families. Our results illustrate the highly pleiotropic nature of RBM10 pathogenic variants, beyond the classic TARP syndrome features. Major clinical characteristics include severe developmental delay, failure to thrive, brain malformations, neurological symptoms, respiratory issues, and facial dysmorphism. Minor features are growth retardation, cardiac, gastrointestinal, limb, and skeletal abnormalities. Additional recurrent features include genital and renal abnormalities as well as hearing and visual impairment. Thus, RBM10 loss of function variants typically cause an intellectual disability and congenital malformation syndrome that requires assessment of multiple organ systems at diagnosis and for which provided clinical features might simplify diagnostic assessment. Furthermore, evidence for an RBM10-related genotype-phenotype correlation is emerging, which can be important for prognosis.

A neuronal enhancer network upstream of MEF2C is compromised in patients with Rett-like characteristics.

Mutations in myocyte enhancer factor 2C (MEF2C), an important transcription factor in neurodevelopment, are associated with a Rett-like syndrome. Structural variants (SVs) upstream of MEF2C, which do not disrupt the gene itself, have also been found in patients with a similar phenotype, suggesting that disruption of MEF2C regulatory elements can also cause a Rett-like phenotype. To characterize those elements that regulate MEF2C during neural development and that are affected by these SVs, we used genomic tools coupled with both in vitro and in vivo functional assays. Through circularized chromosome conformation capture sequencing (4C-seq) and the assay for transposase-accessible chromatin using sequencing (ATAC-seq), we revealed a complex interaction network in which the MEF2C promoter physically contacts several distal enhancers that are deleted or translocated by disease-associated SVs. A total of 16 selected candidate regulatory sequences were tested for enhancer activity in vitro, with 14 found to be functional enhancers. Further analyses of their in vivo activity in zebrafish showed that each of these enhancers has a distinct activity pattern during development, with eight enhancers displaying neuronal activity. In summary, our results disentangle a complex regulatory network governing neuronal MEF2C expression that involves multiple distal enhancers. In addition, the characterized neuronal enhancers pose as novel candidates to screen for mutations in neurodevelopmental disorders, such as Rett-like syndrome.

De novo missense variants in RAC3 cause a novel neurodevelopmental syndrome.

PURPOSE: RAC3 is an underexamined member of the Rho GTPase gene family that is expressed in the developing brain and linked to key cellular functions. De novo missense variants in the homolog RAC1 were recently associated with developmental disorders. In the RAC subfamily, transforming missense changes at certain shared residues have been observed in human cancers and previously characterized in experimental studies. The purpose of this study was to determine whether constitutional dysregulation of RAC3 is associated with human disease. METHODS: We discovered a RAC3 variant in the index case using genome sequencing, and searched for additional variants using international data-sharing initiatives. Functional effects of the variants were assessed using a multifaceted approach generalizable to most clinical laboratory settings. RESULTS: We rapidly identified five individuals with de novo monoallelic missense variants in RAC3, including one recurrent change. Every participant had severe intellectual disability and brain malformations. In silico protein modeling, and prior in vivo and in situ experiments, supported a transforming effect for each of the three different RAC3 variants. All variants were observed in databases of somatic variation in cancer. CONCLUSIONS: Missense variants in RAC3 cause a novel brain disorder, likely through a mechanism of constitutive protein activation.

De novo mutations in MSL3 cause an X-linked syndrome marked by impaired histone H4 lysine 16 acetylation.

The etiological spectrum of ultra-rare developmental disorders remains to be fully defined. Chromatin regulatory mechanisms maintain cellular identity and function, where misregulation may lead to developmental defects. Here, we report pathogenic variations in MSL3, which encodes a member of the chromatin-associated male-specific lethal (MSL) complex responsible for bulk histone H4 lysine 16 acetylation (H4K16ac) in flies and mammals. These variants cause an X-linked syndrome affecting both sexes. Clinical features of the syndrome include global developmental delay, progressive gait disturbance, and recognizable facial dysmorphism. MSL3 mutations affect MSL complex assembly and activity, accompanied by a pronounced loss of H4K16ac levels in vivo. Patient-derived cells display global transcriptome alterations of pathways involved in morphogenesis and cell migration. Finally, we use histone deacetylase inhibitors to rebalance acetylation levels, alleviating some of the molecular and cellular phenotypes of patient cells. Taken together, we characterize a syndrome that allowed us to decipher the developmental importance of MSL3 in humans.

Severe hepatopathy and neurological deterioration after start of valproate treatment in a 6-year-old child with mitochondrial tryptophanyl-tRNA synthetase deficiency.

BACKGROUND: The first subjects with deficiency of mitochondrial tryptophanyl-tRNA synthetase (WARS2) were reported in 2017. Their clinical characteristics can be subdivided into three phenotypes (neonatal phenotype, severe infantile onset phenotype, Parkinson-like phenotype). RESULTS: Here, we report on a subject who presented with early developmental delay, motor weakness and intellectual disability and who was considered during several years as having a non-progressive encephalopathy. At the age of six years, she had an epileptic seizure which was treated with sodium valproate. In the months after treatment was started, she developed acute liver failure and severe progressive encephalopathy. Although valproate was discontinued, she died six months later. Spectrophotometric analysis of the oxidative phosphorylation complexes in liver revealed a deficient activity of complex III and low normal activities of the complexes I and IV. Activity staining in the BN-PAGE gel confirmed the low activities of complex I, III and IV and, in addition, showed the presence of a subcomplex of complex V. Histochemically, a mosaic pattern was seen in hepatocytes after cytochrome c oxidase staining. Using Whole Exome Sequencing two known pathogenic variants were detected in WARS2 (c.797delC, p.Pro266ArgfsTer10/ c.938 A > T, p.Lys313Met). CONCLUSION: This is the first report of severe hepatopathy in a subject with WARS2 deficiency. The hepatopathy occurred soon after start of sodium valproate treatment. In the literature, valproate-induced hepatotoxicity was reported in the subjects with pathogenic mutations in POLG and TWNK. This case report illustrates that the course of the disease in the subjects with a mitochondrial defect can be non-progressive during several years. The subject reported here was first diagnosed as having cerebral palsy. Only after a mitochondriotoxic medication was started, the disease became progressive, and the diagnosis of a mitochondrial defect was made.

CRISPR/Cas9-mediated genome editing in naïve human embryonic stem cells.

The combination of genome-edited human embryonic stem cells (hESCs) and subsequent neural differentiation is a powerful tool to study neurodevelopmental disorders. Since the naïve state of pluripotency has favourable characteristics for efficient genome-editing, we optimized a workflow for the CRISPR/Cas9 system in these naïve stem cells. Editing efficiencies of respectively 1.3-8.4% and 3.8-19% were generated with the Cas9 nuclease and the D10A Cas9 nickase mutant. Next to this, wildtype and genome-edited naïve hESCs were successfully differentiated to neural progenitor cells. As a proof-of-principle of our workflow, two monoclonal genome-edited naïve hESCs colonies were obtained for TUNA, a long non-coding RNA involved in pluripotency and neural differentiation. In these genome-edited hESCs, an effect was seen on expression of TUNA, although not on neural differentiation potential. In conclusion, we optimized a genome-editing workflow in naïve hESCs that can be used to study candidate genes involved in neural differentiation and/or functioning.