Genome-wide identification of tandem repeats associated with splicing variation across 49 tissues in humans.

Tandem repeats (TRs) are one of the largest sources of polymorphism, and their length is associated with gene regulation. Although previous studies reported several tandem repeats regulating gene splicing in cis (spl-TRs), no large-scale study has been conducted. In this study, we established a genome-wide catalog of 9537 spl-TRs with a total of 58,290 significant TR-splicing associations across 49 tissues (false discovery rate 5%) by using Genotype-Tissue expression (GTex) Project data. Regression models explaining splicing variation by using spl-TRs and other flanking variants suggest that at least some of the spl-TRs directly modulate splicing. In our catalog, two spl-TRs are known loci for repeat expansion diseases, spinocerebellar ataxia 6 (SCA6) and 12 (SCA12). Splicing alterations by these spl-TRs were compatible with those observed in SCA6 and SCA12. Thus, our comprehensive spl-TR catalog may help elucidate the pathomechanism of genetic diseases.

Complex chromosomal 6q rearrangements revealed by combined long-molecule genomics technologies.

Technologies for detecting structural variation (SV) have advanced with the advent of long-read sequencing, which enables the validation of SV at a nucleotide level. Optical genome mapping (OGM), a technology based on physical mapping, can also provide comprehensive SVs analysis. We applied long-read whole genome sequencing (LRWGS) to accurately reconstruct breakpoint (BP) segments in a patient with complex chromosome 6q rearrangements that remained elusive by conventional karyotyping. Although all BPs were precisely identified by LRWGS, there were two possible ways to construct the BP segments in terms of their orders and orientations. Thus, we also used OGM analysis. Notably, OGM recognized entire inversions exceeding 500 kb in size, which LRWGS could not characterize. Consequently, here we successfully unveil the full genomic structure of this complex chromosomal 6q rearrangement and cryptic SVs through combined long-molecule genomic analyses, showcasing how LRWGS and OGM can complement each other in SV analysis.

Detection of hidden intronic DDC variant in aromatic L-amino acid decarboxylase deficiency by adaptive sampling.

Aromatic l-amino acid decarboxylase (AADC) deficiency is an autosomal recessive neurotransmitter disorder caused by pathogenic DOPA decarboxylase (DDC) variants. We previously reported Japanese siblings with AADC deficiency, which was confirmed by the lack of enzyme activity; however, only a heterozygous missense variant was detected. We therefore performed targeted long-read sequencing by adaptive sampling to identify any missing variants. Haplotype phasing and variant calling identified a novel deep intronic variant (c.714+255 C > A), which was predicted to potentially activate the noncanonical splicing acceptor site. Minigene assay revealed that wild-type and c.714+255 C > A alleles had different impacts on splicing. Three transcripts, including the canonical transcript, were detected from the wild-type allele, but only the noncanonical cryptic exon was produced from the variant allele, indicating that c.714+255 C > A was pathogenic. Target long-read sequencing may be used to detect hidden pathogenic variants in unresolved autosomal recessive cases with only one disclosed hit variant.

Novel missense variants cause intermediate phenotypes in the phenotypic spectrum of SLC5A6-related disorders.

SLC5A6 encodes the sodium-dependent multivitamin transporter, a transmembrane protein that uptakes biotin, pantothenic acid, and lipoic acid. Biallelic SLC5A6 variants cause sodium-dependent multivitamin transporter deficiency (SMVTD) and childhood-onset biotin-responsive peripheral motor neuropathy (COMNB), which both respond well to replacement therapy with the above three nutrients. SMVTD usually presents with various symptoms in multiple organs, such as gastrointestinal hemorrhage, brain atrophy, and global developmental delay, at birth or in infancy. Without nutrient replacement therapy, SMVTD can be lethal in early childhood. COMNB is clinically milder and has a later onset than SMVTD, at approximately 10 years of age. COMNB symptoms are mostly limited to peripheral motor neuropathy. Here we report three patients from one Japanese family harboring novel compound heterozygous missense variants in SLC5A6, namely NM_021095.4:c.[221C>T];[642G>C] p.[(Ser74Phe)];[(Gln214His)]. Both variants were predicted to be deleterious through multiple lines of evidence, including amino acid conservation, in silico predictions of pathogenicity, and protein structure considerations. Drosophila analysis also showed c.221C>T to be pathogenic. All three patients had congenital brain cysts on neonatal cranial imaging, but no other morphological abnormalities. They also had a mild motor developmental delay that almost completely resolved despite no treatment. In terms of severity, their phenotypes were intermediate between SMVTD and COMNB. From these findings we propose a new SLC5A6-related disorder, spontaneously remitting developmental delay with brain cysts (SRDDBC) whose phenotypic severity is between that of SMVTD and COMNB. Further clinical and genetic evidence is needed to support our suggestion.

Novel compound heterozygous ABCA2 variants cause IDPOGSA, a variable phenotypic syndrome with intellectual disability.

The gene for ATP binding cassette subfamily A member 2 (ABCA2) is located at chromosome 9q34.3. Biallelic ABCA2 variants lead to intellectual developmental disorder with poor growth and with or without seizures or ataxia (IDPOGSA). In this study, we identified novel compound heterozygous ABCA2 variants (NM_001606.5:c.[5300-17C>A];[6379C>T]) by whole exome sequencing in a 28-year-old Korean female patient with intellectual disability. These variants included intronic and nonsense variants of paternal and maternal origin, respectively, and are absent from gnomAD. SpliceAI predicted that the intron variant creates a cryptic acceptor site. Reverse transcription-PCR using RNA extracted from a lymphoblastoid cell line of the patient confirmed two aberrant transcripts. Her clinical features are compatible with those of IDPOGSA.

Complete nanopore repeat sequencing of SCA27B (GAA-FGF14 ataxia) in Japanese.

BACKGROUND: Although pure GAA expansion is considered pathogenic in SCA27B, non-GAA repeat motif is mostly mixed into longer repeat sequences. This study aimed to unravel the complete sequencing of FGF14 repeat expansion to elucidate its repeat motifs and pathogenicity. METHODS: We screened FGF14 repeat expansion in a Japanese cohort of 460 molecularly undiagnosed adult-onset cerebellar ataxia patients and 1022 controls, together with 92 non-Japanese controls, and performed nanopore sequencing of FGF14 repeat expansion. RESULTS: In the Japanese population, the GCA motif was predominantly observed as the non-GAA motif, whereas the GGA motif was frequently detected in non-Japanese controls. The 5'-common flanking variant was observed in all Japanese GAA repeat alleles within normal length, demonstrating its meiotic stability against repeat expansion. In both patients and controls, pure GAA repeat was up to 400 units in length, whereas non-pathogenic GAA-GCA repeat was larger, up to 900 units, but they evolved from different haplotypes, as rs534066520, located just upstream of the repeat sequence, completely discriminated them. Both (GAA)≥250 and (GAA)≥200 were enriched in patients, whereas (GAA-GCA)≥200 was similarly observed in patients and controls, suggesting the pathogenic threshold of (GAA)≥200 for cerebellar ataxia. We identified 14 patients with SCA27B (3.0%), but their single-nucleotide polymorphism genotype indicated different founder alleles between Japanese and Caucasians. The low prevalence of SCA27B in Japanese may be due to the lower allele frequency of (GAA)≥250 in the Japanese population than in Caucasians (0.15% vs 0.32%-1.26%). CONCLUSIONS: FGF14 repeat expansion has unique features of pathogenicity and allelic origin, as revealed by a single ethnic study.

De novo ARF3 variants cause neurodevelopmental disorder with brain abnormality.

An optimal Golgi transport system is important for mammalian cells. The adenosine diphosphate (ADP) ribosylation factors (ARF) are key proteins for regulating cargo sorting at the Golgi network. In this family, ARF3 mainly works at the trans-Golgi network (TGN), and no ARF3-related phenotypes have yet been described in humans. We here report the clinical and genetic evaluations of two unrelated children with de novo pathogenic variants in the ARF3 gene: c.200A > T (p.Asp67Val) and c.296G > T (p.Arg99Leu). Although the affected individuals presented commonly with developmental delay, epilepsy and brain abnormalities, there were differences in severity, clinical course and brain lesions. In vitro subcellular localization assays revealed that the p.Arg99Leu mutant localized to Golgi apparatus, similar to the wild-type, whereas the p.Asp67Val mutant tended to show a disperse cytosolic pattern together with abnormally dispersed Golgi localization, similar to that observed in a known dominant negative variant (p.Thr31Asn). Pull-down assays revealed that the p.Asp67Val had a loss-of-function effect and the p.Arg99Leu variant had increased binding of the adaptor protein, Golgi-localized, γ-adaptin ear-containing, ARF-binding protein 1 (GGA1), supporting the gain of function. Furthermore, in vivo studies revealed that p.Asp67Val transfection led to lethality in flies. In contrast, flies expressing p.Arg99Leu had abnormal rough eye, as observed in the gain-of-function variant p.Gln71Leu. These data indicate that two ARF3 variants, the possibly loss-of-function p.Asp67Val and the gain-of-function p.Arg99Leu, both impair the Golgi transport system. Therefore, it may not be unreasonable that they showed different clinical features like diffuse brain atrophy (p.Asp67Val) and cerebellar hypoplasia (p.Arg99Leu).

De Novo Truncating Variants in the Last Exon of SEMA6B Cause Progressive Myoclonic Epilepsy.

De novo variants (DNVs) cause many genetic diseases. When DNVs are examined in the whole coding regions of genes in next-generation sequencing analyses, pathogenic DNVs often cluster in a specific region. One such region is the last exon and the last 50 bp of the penultimate exon, where truncating DNVs cause escape from nonsense-mediated mRNA decay [NMD(-) region]. Such variants can have dominant-negative or gain-of-function effects. Here, we first developed a resource of rates of truncating DNVs in NMD(-) regions under the null model of DNVs. Utilizing this resource, we performed enrichment analysis of truncating DNVs in NMD(-) regions in 346 developmental and epileptic encephalopathy (DEE) trios. We observed statistically significant enrichment of truncating DNVs in semaphorin 6B (SEMA6B) (p value: 2.8 × 10-8; exome-wide threshold: 2.5 × 10-6). The initial analysis of the 346 individuals and additional screening of 1,406 and 4,293 independent individuals affected by DEE and developmental disorders collectively identified four truncating DNVs in the SEMA6B NMD(-) region in five individuals who came from unrelated families (p value: 1.9 × 10-13) and consistently showed progressive myoclonic epilepsy. RNA analysis of lymphoblastoid cells established from an affected individual showed that the mutant allele escaped NMD, indicating stable production of the truncated protein. Importantly, heterozygous truncating variants in the NMD(+) region of SEMA6B are observed in general populations, and SEMA6B is most likely loss-of-function tolerant. Zebrafish expressing truncating variants in the NMD(-) region of SEMA6B orthologs displayed defective development of brain neurons and enhanced pentylenetetrazole-induced seizure behavior. In summary, we show that truncating DNVs in the final exon of SEMA6B cause progressive myoclonic epilepsy.

Gain-of-Function MN1 Truncation Variants Cause a Recognizable Syndrome with Craniofacial and Brain Abnormalities.

MN1 was originally identified as a tumor-suppressor gene. Knockout mouse studies have suggested that Mn1 is associated with craniofacial development. However, no MN1-related phenotypes have been established in humans. Here, we report on three individuals who have de novo MN1 variants that lead to a protein lacking the carboxyl (C) terminus and who presented with severe developmental delay, craniofacial abnormalities with specific facial features, and structural abnormalities in the brain. An in vitro study revealed that the deletion of the C-terminal region led to increased protein stability, an inhibitory effect on cell proliferation, and enhanced MN1 aggregation in nuclei compared to what occurred in the wild type, suggesting that a gain-of-function mechanism is involved in this disease. Considering that C-terminal deletion increases the fraction of intrinsically disordered regions of MN1, it is possible that altered phase separation could be involved in the mechanism underlying the disease. Our data indicate that MN1 participates in transcriptional regulation of target genes through interaction with the transcription factors PBX1, PKNOX1, and ZBTB24 and that mutant MN1 impairs the binding with ZBTB24 and RING1, which is an E3 ubiquitin ligase. On the basis of our findings, we propose the model that C-terminal deletion interferes with MN1's interaction molecules related to the ubiquitin-mediated proteasome pathway, including RING1, and increases the amount of the mutant protein; this increase leads to the dysregulation of MN1 target genes by inhibiting rapid MN1 protein turnover.

Complete SAMD12 repeat expansion sequencing in a four-generation BAFME1 family with anticipation.

Benign adult familial myoclonic epilepsy type 1 (BAFME1) is an autosomal dominant, adult-onset neurological disease caused by SAMD12 repeat expansion. In BAFME1, anticipation, such as the earlier onset of tremor and/or seizures in the next generation, was reported. This could be explained by intergenerational repeat instability, leading to larger expansions in successive generations. We report a four-generation BAFME1-affected family with anticipation. Using Nanopore long-read sequencing, detailed information regarding the sizes, configurations, and compositions of the expanded SAMD12 repeats across generations was obtained. Unexpectedly, a grandmother-mother-daughter triad showed similar repeat structures but with slight repeat expansions, despite quite variable age of onset of seizures (range: 52-14 years old), implying a complex relationship between the SAMD12 repeat expansion sequence and anticipation. This study suggests that different factor(s) from repeat expansion could modify the anticipation in BAFME1.

Distal arthrogryposis in a girl arising from a novel TNNI2 variant inherited from paternal somatic mosaicism.

TNNI2 at 11p15.5 encodes troponin I2, fast skeletal type, which is a member of the troponin I gene family and a component of the troponin complex. Distal arthrogryposis (DA) is characterized by congenital limb contractures without primary neurological or muscular effects. DA is inherited in an autosomal dominant fashion and is clinically and genetically heterogeneous. Exome sequencing identified a causative variant in TNNI2 [NM_003282.4:c.532T>C p.(Phe178Leu)] in a Japanese girl with typical DA2b. Interestingly, the familial study using Sanger sequencing suggested a mosaic variant in her healthy father. Subsequent targeted amplicon-based deep sequencing detected the TNNI2 variant with variant allele frequencies of 9.4-17.7% in genomic DNA derived from peripheral blood leukocytes, saliva, hair, and nails in the father. We confirmed a disease-causing variant in TNNI2 in the proband inherited from her asymptomatic father with its somatic variant. Our case demonstrates that careful clinical and genetic evaluation is required in DA.

A novel homozygous CHMP1A variant arising from segmental uniparental disomy causes pontocerebellar hypoplasia type 8.

Pontocerebellar hypoplasia (PCH) is currently classified into 16 subgroups. Using mostly next-generation sequencing, pathogenic variants have been identified in as many as 24 PCH-associated genes. PCH type 8 (PCH8) is a rare heterogeneous disorder. Its clinical presentation includes severe development delay, increased muscle tone, microcephaly, and magnetic resonance imaging (MRI) abnormalities such as reduced cerebral white matter, a thin corpus callosum, and brainstem and cerebellar hypoplasia. To date, only two variants in the CHMP1A gene (MIM: 164010), NM_002768.5: c.88 C > T (p.Glu30*) and c.28-13 G > A, have been identified homozygously in seven patients with PCH8 from four families (MIM: 614961). CHMP1A is a subunit of the endosomal sorting complex required for transport III (ESCRT-III), which regulates the formation and release of extracellular vesicles. Biallelic CHMP1A loss of function impairs the ESCRT-III-mediated release of extracellular vesicles, which causes impaired progenitor proliferation in the developing brain. Herein, we report a patient with PCH8 who had a homozygous CHMP1A variant, c.122delA (p.Asn41Metfs*2), which arose from segmental uniparental disomy. Although our patient had similar MRI findings to those of previously reported patients, with no progression, we report some novel neurological and developmental findings that expand our knowledge of the clinical consequences associated with CHMP1A variants.

Long-read sequencing revealing intragenic deletions in exome-negative spastic paraplegias.

Hereditary spastic paraplegias (HSPs) are a heterogeneous group of neurodegenerative disorders characterized by progressive spasticity and weakness in the lower extremities. To date, a total of 88 types of SPG are known. To diagnose HSP, multiple technologies, including microarray, direct sequencing, multiplex ligation-dependent probe amplification, and short-read next-generation sequencing, are often chosen based on the frequency of HSP subtypes. Exome sequencing (ES) is commonly used. We used ES to analyze ten cases of HSP from eight families. We identified pathogenic variants in three cases (from three different families); however, we were unable to determine the cause of the other seven cases using ES. We therefore applied long-read sequencing to the seven undetermined HSP cases (from five families). We detected intragenic deletions within the SPAST gene in four families, and a deletion within PSEN1 in the remaining family. The size of the deletion ranged from 4.7 to 12.5 kb and involved 1-7 exons. All deletions were entirely included in one long read. We retrospectively performed an ES-based copy number variation analysis focusing on pathogenic deletions, but were not able to accurately detect these deletions. This study demonstrated the efficiency of long-read sequencing in detecting intragenic pathogenic deletions in ES-negative HSP patients.

Three KINSSHIP syndrome patients with mosaic and germline AFF3 variants.

AFF3 at 2q11.2 encodes the nuclear transcriptional activator AF4/FMR2 Family Member 3. AFF3 constitutes super elongation complex like 3, which plays a role in promoting the expression of genes involved in neurogenesis and development. The degron motif in AFF3 with nine highly conserved amino acids is recognized by E3 ubiquitin ligase to induce protein degradation. Recently, AFF3 missense variants in this region and variants featuring deletion including this region were identified and shown to cause KINSSHIP syndrome. In this study, we identified two novel and one previously reported missense variants in the degron of AFF3 in three unrelated Japanese patients. Notably, two of these three variants exhibited mosaicism in the examined tissues. This study suggests that mosaic variants also cause KINSSHIP syndrome, showing various phenotypes.

Repeat conformation heterogeneity in cerebellar ataxia, neuropathy, vestibular areflexia syndrome.

Cerebellar ataxia, neuropathy, vestibular areflexia syndrome (CANVAS) is a late-onset, slow-progressing multisystem neurodegenerative disorder. Biallelic AAGGG repeat expansion in RFC1 has been identified as causative of this disease, and repeat conformation heterogeneity (ACAGG repeat) was also recently implied. To molecularly characterize this disease in Japanese patients with adult-onset ataxia, we accumulated and screened 212 candidate families by an integrated approach consisting of flanking PCR, repeat-primed PCR, Southern blotting and long-read sequencing using Sequel II, GridION or PromethION. We identified 16 patients from 11 families, of whom seven had ACAGG expansions [(ACAGG)exp/(ACAGG)exp] (ACAGG homozygotes), two had ACAGG and AAGGG expansions [(ACAGG)exp/(AAGGG)exp] (ACAGG/AAGGG compound heterozygotes) and seven had AAGGG expansions [(AAGGG)exp/(AAGGG)exp] (AAGGG homozygotes). The overall detection rate was 5.2% (11/212 families including one family having two expansion genotypes). Long-read sequencers revealed the entire sequence of both AAGGG and ACAGG repeat expansions at the nucleotide level of resolution. Clinical assessment and neuropathology results suggested that patients with ACAGG expansions have similar clinical features to previously reported patients with homozygous AAGGG expansions, although motor neuron involvement was more notable in patients with ACAGG expansions (even if one allele was involved). Furthermore, a later age of onset and slower clinical progression were implied in patients with ACAGG/AAGGG compound heterozygous expansions compared with either ACAGG or AAGGG homozygotes in our very limited cohort. Our study clearly shows the occurrence of repeat conformation heterogeneity, with possible different impacts on the affected nervous systems. The difference in disease onset and progression between compound heterozygotes and homozygotes might also be suspected but with very limited certainty due to the small sample number of cases in our study. Studies of additional patients are needed to confirm this.

Patients with biallelic GGC repeat expansions in NOTCH2NLC exhibiting a typical neuronal intranuclear inclusion disease phenotype.

We report two patients with autosomal dominant neuronal intranuclear inclusion disease (NIID) harboring the biallelic GGC repeat expansion in NOTCH2NLC to uncover the impact of repeat expansion zygosity on the clinical phenotype. The zygosity of the entire NOTCH2NLC GGC repeat expansion and DNA methylation were comprehensively evaluated using fluorescent amplicon length PCR (AL-PCR), Southern blotting and targeted long-read sequencing, and detailed genetic/epigenetic and clinical features were described. In AL-PCR, we could not recognize the wild-type allele in both patients. Targeted long-read sequencing revealed that one patient harbored a homozygous repeat expansion. The other patient harbored compound heterozygous repeat expansions. The GGC repeats and the nearest CpG island were hypomethylated in all expanded alleles in both patients. Both patients harboring the biallelic GGC repeat expansion showed a typical dementia-dominant NIID phenotype. In conclusion, the biallelic GGC repeat expansion in two typical NIID patients indicated that NOTCH2NLC-related diseases could be completely dominant.

Pathogenic variants detected by RNA sequencing in Cornelia de Lange syndrome.

Recent studies suggest that transcript isoforms significantly overlap (approximately 60%) between brain tissue and Epstein-Barr virus-transformed lymphoblastoid cell lines (LCLs). Interestingly, 14 cohesion-related genes with variants that cause Cornelia de Lange Syndrome (CdLS) are highly expressed in the brain and LCLs. In this context, we first performed RNA sequencing of LCLs from 22 solved (with pathogenic variants) and 19 unsolved (with no confirmed variants) CdLS cases. Next, an RNA sequencing pipeline was developed using solved cases with two different methods: short variant analysis (for single-nucleotide and indel variants) and aberrant splicing detection analysis. Then, 19 unsolved cases were subsequently applied to our pipeline, and four pathogenic variants in NIPBL (one inframe deletion and three intronic variants) were newly identified. Two of three intronic variants were located at Alu elements in deep-intronic regions, creating cryptic exons. RNA sequencing with LCLs was useful for identifying hidden variants in exome-negative cases.

Amelioration of a neurodevelopmental disorder by carbamazepine in a case having a gain-of-function GRIA3 variant.

GRIA3 at Xq25 encodes glutamate ionotropic receptor AMPA type 3 (GluA3), a subunit of postsynaptic glutamate-gated ion channels mediating neurotransmission. Hemizygous loss-of-function (LOF) variants in GRIA3 cause a neurodevelopmental disorder (NDD) in male individuals. Here, we report a gain-of-function (GOF) variant at GRIA3 in a male patient. We identified a hemizygous de novo missense variant in GRIA3 in a boy with an NDD: c.1844C > T (p.Ala615Val) using whole-exome sequencing. His neurological signs, such as hypertonia and hyperreflexia, were opposite to those in previous cases having LOF GRIA3 variants. His seizures and hypertonia were ameliorated by carbamazepine, inhibiting glutamate release from presynapses. Patch-clamp recordings showed that the human GluA3 mutant (p.Ala615Val) had slower desensitization and deactivation kinetics. A fly line expressing a human GluA3 mutant possessing our variant and the Lurcher variant, which makes ion channels leaky, showed developmental defects, while one expressing a mutant possessing either of them did not. Collectively, these results suggest that p.Ala615Val has GOF effects. GRIA3 GOF variants may cause an NDD phenotype distinctive from that of LOF variants, and drugs suppressing glutamatergic neurotransmission may ameliorate this phenotype. This study should help in refining the clinical management of GRIA3-related NDDs.

Monogenic causes of pigmentary mosaicism.

Pigmentary mosaicism of the Ito type, also known as hypomelanosis of Ito, is a neurocutaneous syndrome considered to be predominantly caused by somatic chromosomal mosaicism. However, a few monogenic causes of pigmentary mosaicism have been recently reported. Eleven unrelated individuals with pigmentary mosaicism (mostly hypopigmented skin) were recruited for this study. Skin punch biopsies of the probands and trio-based blood samples (from probands and both biological parents) were collected, and genomic DNA was extracted and analyzed by exome sequencing. In all patients, plausible monogenic causes were detected with somatic and germline variants identified in five and six patients, respectively. Among the somatic variants, four patients had MTOR variant (36%) and another had an RHOA variant. De novo germline variants in USP9X, TFE3, and KCNQ5 were detected in two, one, and one patients, respectively. A maternally inherited PHF6 variant was detected in one patient with hyperpigmented skin. Compound heterozygous GTF3C5 variants were highlighted as strong candidates in the remaining patient. Exome sequencing, using patients' blood and skin samples is highly recommended as the first choice for detecting causative genetic variants of pigmentary mosaicism.

Biallelic null variants in ZNF142 cause global developmental delay with familial epilepsy and dysmorphic features.

Biallelic variants in ZNF142 at 2q35, which encodes zinc-finger protein 142, cause neurodevelopmental disorder with seizures or dystonia. We identified compound heterozygous null variants in ZNF142, NM_001105537.4:c.[1252C>T];[1274-2A>G],p.[Arg418*];[Glu426*], in Malaysian siblings suffering from global developmental delay with epilepsy and dysmorphism. cDNA analysis showed the marked reduction of ZNF142 transcript level through nonsense-mediated mRNA decay by these novel biallelic variants. The affected siblings present with global developmental delay and epilepsy in common, which were previously described, as well as dysmorphism, which was not recognized. It is important to collect patients with ZNF142 abnormality to define its phenotypic spectrum.