Through a multi-center collaboration study, we here report six individuals from five unrelated families, with mutations in KAT6A/MOZ detected by whole-exome sequencing. All five different de novo heterozygous truncating mutations were located in the C-terminal transactivation domain of KAT6A: NM_001099412.1: c.3116_3117 delCT, p.(Ser1039∗); c.3830_3831insTT, p.(Arg1278Serfs∗17); c.3879 dupA, p.(Glu1294Argfs∗19); c.4108G>T p.(Glu1370∗) and c.4292 dupT, p.(Leu1431Phefs∗8). An additional subject with a 0.23 MB microdeletion including the entire KAT6A reading frame was identified with genome-wide array comparative genomic hybridization. Finally, by detailed clinical characterization we provide evidence that heterozygous mutations in KAT6A cause a distinct intellectual disability syndrome. The common phenotype includes hypotonia, intellectual disability, early feeding and oromotor difficulties, microcephaly and/or craniosynostosis, and cardiac defects in combination with subtle facial features such as bitemporal narrowing, broad nasal tip, thin upper lip, posteriorly rotated or low-set ears, and microretrognathia. The identification of human subjects complements previous work from mice and zebrafish where knockouts of Kat6a/kat6a lead to developmental defects.

We present a method that utilizes DNA methylation profiling for prediction of the cytogenetic subtypes of acute lymphoblastic leukemia (ALL) cells from pediatric ALL patients. The primary aim of our study was to improve risk stratification of ALL patients into treatment groups using DNA methylation as a complement to current diagnostic methods. A secondary aim was to gain insight into the functional role of DNA methylation in ALL.

Despite the high heritability of autism spectrum disorders (ASD), characterized by persistent deficits in social communication and interaction and restricted, repetitive patterns of behavior, interests or activities, a genetic diagnosis can be established in only a minority of patients. Known genetic causes include chromosomal aberrations, such as the duplication of the 15q11-13 region, and monogenic causes, as in Rett and fragile-X syndromes. The genetic heterogeneity within ASD is striking, with even the most frequent causes responsible for only 1% of cases at the most. Even with the recent developments in next-generation sequencing, for the large majority of cases no molecular diagnosis can be established. Here, we report ten patients with ASD and other shared clinical characteristics, including intellectual disability and facial dysmorphisms caused by a mutation in ADNP, a transcription factor involved in the SWI/SNF remodeling complex. We estimate this gene to be mutated in at least 0.17% of ASD cases, making it one of the most frequent ASD-associated genes known to date.

Mutations in more than a hundred genes have been reported to cause X-linked recessive intellectual disability (ID) mainly in males. In contrast, the number of identified X-linked genes in which de novo mutations specifically cause ID in females is limited. Here, we report 17 females with de novo loss-of-function mutations in USP9X, encoding a highly conserved deubiquitinating enzyme. The females in our study have a specific phenotype that includes ID/developmental delay (DD), characteristic facial features, short stature, and distinct congenital malformations comprising choanal atresia, anal abnormalities, post-axial polydactyly, heart defects, hypomastia, cleft palate/bifid uvula, progressive scoliosis, and structural brain abnormalities. Four females from our cohort were identified by targeted genetic testing because their phenotype was suggestive for USP9X mutations. In several females, pigment changes along Blaschko lines and body asymmetry were observed, which is probably related to differential (escape from) X-inactivation between tissues. Expression studies on both mRNA and protein level in affected-female-derived fibroblasts showed significant reduction of USP9X level, confirming the loss-of-function effect of the identified mutations. Given that some features of affected females are also reported in known ciliopathy syndromes, we examined the role of USP9X in the primary cilium and found that endogenous USP9X localizes along the length of the ciliary axoneme, indicating that its loss of function could indeed disrupt cilium-regulated processes. Absence of dysregulated ciliary parameters in affected female-derived fibroblasts, however, points toward spatiotemporal specificity of ciliary USP9X (dys-)function.

Genetic studies of intellectual disability and identification of monogenic causes of obesity in humans have made immense contribution toward the understanding of the brain and control of body mass. The leptin > melanocortin > SIM1 pathway is dysregulated in multiple monogenic human obesity syndromes but its downstream targets are still unknown. In ten individuals from six families, with overlapping 6q16.1 deletions, we describe a disorder of variable developmental delay, intellectual disability, and susceptibility to obesity and hyperphagia. The 6q16.1 deletions segregated with the phenotype in multiplex families and were shown to be de novo in four families, and there was dramatic phenotypic overlap among affected individuals who were independently ascertained without bias from clinical features. Analysis of the deletions revealed a ∼350 kb critical region on chromosome 6q16.1 that encompasses a gene for proneuronal transcription factor POU3F2, which is important for hypothalamic development and function. Using morpholino and mutant zebrafish models, we show that POU3F2 lies downstream of SIM1 and controls oxytocin expression in the hypothalamic neuroendocrine preoptic area. We show that this finding is consistent with the expression patterns of POU3F2 and related genes in the human brain. Our work helps to further delineate the neuro-endocrine control of energy balance/body mass and demonstrates that this molecular pathway is conserved across multiple species.

Clustered copy number variants (CNVs) as detected by chromosomal microarray analysis (CMA) are often reported as germline chromothripsis. However, such cases might need further investigations by massive parallel whole genome sequencing (WGS) in order to accurately define the underlying complex rearrangement, predict the occurrence mechanisms and identify additional complexities. Here, we utilized WGS to delineate the rearrangement structure of 21 clustered CNV carriers first investigated by CMA and identified a total of 83 breakpoint junctions (BPJs). The rearrangements were further sub-classified depending on the patterns observed: I) Cases with only deletions (n = 8) often had additional structural rearrangements, such as insertions and inversions typical to chromothripsis; II) cases with only duplications (n = 7) or III) combinations of deletions and duplications (n = 6) demonstrated mostly interspersed duplications and BPJs enriched with microhomology. In two cases the rearrangement mutational signatures indicated both a breakage-fusion-bridge cycle process and haltered formation of a ring chromosome. Finally, we observed two cases with Alu- and LINE-mediated rearrangements as well as two unrelated individuals with seemingly identical clustered CNVs on 2p25.3, possibly a rare European founder rearrangement. In conclusion, through detailed characterization of the derivative chromosomes we show that multiple mechanisms are likely involved in the formation of clustered CNVs and add further evidence for chromoanagenesis mechanisms in both "simple" and highly complex chromosomal rearrangements. Finally, WGS characterization adds positional information, important for a correct clinical interpretation and deciphering mechanisms involved in the formation of these rearrangements.

RNA binding proteins are key players in posttranscriptional regulation and have been implicated in neurodevelopmental and neuropsychiatric disorders. Here, we report a significant burden of heterozygous, likely gene-disrupting variants in CSDE1 (encoding a highly constrained RNA binding protein) among patients with autism and related neurodevelopmental disabilities. Analysis of 17 patients identifies common phenotypes including autism, intellectual disability, language and motor delay, seizures, macrocephaly, and variable ocular abnormalities. HITS-CLIP revealed that Csde1-binding targets are enriched in autism-associated gene sets, especially FMRP targets, and in neuronal development and synaptic plasticity-related pathways. Csde1 knockdown in primary mouse cortical neurons leads to an overgrowth of the neurites and abnormal dendritic spine morphology/synapse formation and impaired synaptic transmission, whereas mutant and knockdown experiments in Drosophila result in defects in synapse growth and synaptic transmission. Our study defines a new autism-related syndrome and highlights the functional role of CSDE1 in synapse development and synaptic transmission.

Most genes associated with neurodevelopmental disorders (NDDs) were identified with an excess of de novo mutations (DNMs) but the significance in case-control mutation burden analysis is unestablished. Here, we sequence 63 genes in 16,294 NDD cases and an additional 62 genes in 6,211 NDD cases. By combining these with published data, we assess a total of 125 genes in over 16,000 NDD cases and compare the mutation burden to nonpsychiatric controls from ExAC. We identify 48 genes (25 newly reported) showing significant burden of ultra-rare (MAF < 0.01%) gene-disruptive mutations (FDR 5%), six of which reach family-wise error rate (FWER) significance (p < 1.25E-06). Among these 125 targeted genes, we also reevaluate DNM excess in 17,426 NDD trios with 6,499 new autism trios. We identify 90 genes enriched for DNMs (FDR 5%; e.g., GABRG2 and UIMC1); of which, 61 reach FWER significance (p < 3.64E-07; e.g., CASZ1). In addition to doubling the number of patients for many NDD risk genes, we present phenotype-genotype correlations for seven risk genes (CTCF, HNRNPU, KCNQ3, ZBTB18, TCF12, SPEN, and LEO1) based on this large-scale targeted sequencing effort.

Individuals with Williams syndrome (WS) have an elevated risk for anxiety disorders throughout the life span, making it a research priority to identify the individual factors associated with anxiety. Most of the existing literature is based on questionnaire data and suggests that impaired executive functions (EF) increase the risk for anxiety in WS. The aim of this study was to use direct measures by trained clinicians to investigate the effects of general intelligence, inhibition, sustained attention, and working memory on anxiety in WS, to further elucidate potential underlying mechanisms.

Pathogenic variants in MYH3 cause distal arthrogryposis type 2A and type 2B3 as well as contractures, pterygia and spondylocarpotarsal fusion syndromes types 1A and 1B. These disorders are ultra-rare and their natural course and phenotypic variability are not well described. In this study, we summarize the clinical features and genetic findings of 17 patients from 10 unrelated families with vertebral malformations caused by dominant or recessive pathogenic variants in MYH3. Twelve novel pathogenic variants in MYH3 (NM_002470.4) were identified: three of them were de novo or inherited in autosomal dominant way and nine were inherited in autosomal recessive way. The patients had vertebral segmentation anomalies accompanied with variable joint contractures, short stature and dysmorphic facial features. There was a significant phenotypic overlap between dominant and recessive MYH3-associated conditions regarding the degree of short stature as well as the number of vertebral fusions. All monoallelic variants caused significantly decreased SMAD3 phosphorylation, which is consistent with the previously proposed pathogenic mechanism of impaired canonical TGF-β signaling. Most of the biallelic variants were predicted to be protein-truncating, while one missense variant c.4244T>G,p.(Leu1415Arg), which was inherited in an autosomal recessive way, was found to alter the phosphorylation level of p38, suggesting an inhibition of the non-canonical pathway of TGF-β signaling. In conclusion, the identification of 12 novel pathogenic variants and overlapping phenotypes in 17 affected individuals from 10 unrelated families expands the mutation and phenotype spectrum of MYH3-associated skeletal disorders. We show that disturbances of canonical or non-canonical TGF-β signaling pathways are involved in pathogenesis of MYH3-associated skeletal fusion (MASF) syndrome.

Spondyloepimetaphyseal dysplasia with severe short stature, RPL13-related (SEMD-RPL13), MIM#618728), is a rare autosomal dominant disorder characterized by short stature and skeletal changes such as mild spondylar and epimetaphyseal dysplasia affecting primarily the lower limbs. The genetic cause was first reported in 2019 by Le Caignec et al., and six disease-causing variants in the gene coding for a ribosomal protein, RPL13 (NM_000977.3) have been identified to date. This study presents clinical and radiographic data from 12 affected individuals aged 2-64 years from seven unrelated families, showing highly variable manifestations. The affected individuals showed a range from mild to severe short stature, retaining the same radiographic pattern of spondylar- and epi-metaphyseal dysplasia, but with varying severity of the hip and knee deformities. Two new missense variants, c.548 G>A, p.(Arg183His) and c.569 G>T, p.(Arg190Leu), and a previously known splice variant c.477+1G>A were identified, confirming mutational clustering in a highly specific RNA binding motif. Structural analysis and interpretation of the variants' impact on the protein suggests that disruption of extra-ribosomal functions of the protein through binding of mRNA may play a role in the skeletal phenotype of SEMD-RPL13. In addition, we present gonadal and somatic mosaicism for the condition.

The detection of recurrent somatic chromosomal rearrangements is standard of care for most leukemia types. Even though karyotype analysis-a low-resolution genome-wide chromosome analysis-is still the gold standard, it often needs to be complemented with other methods to increase resolution. To evaluate the feasibility and applicability of mate pair whole genome sequencing (MP-WGS) to detect structural chromosomal rearrangements in the diagnostic setting, we sequenced ten bone marrow samples from leukemia patients with recurrent rearrangements. Samples were selected based on cytogenetic and FISH results at leukemia diagnosis to include common rearrangements of prognostic relevance. Using MP-WGS and in-house bioinformatic analysis all sought rearrangements were successfully detected. In addition, unexpected complexity or additional, previously undetected rearrangements was unraveled in three samples. Finally, the MP-WGS analysis pinpointed the location of chromosome junctions at high resolution and we were able to identify the exact exons involved in the resulting fusion genes in all samples and the specific junction at the nucleotide level in half of the samples. The results show that our approach combines the screening character from karyotype analysis with the specificity and resolution of cytogenetic and molecular methods. As a result of the straightforward analysis and high-resolution detection of clinically relevant rearrangements, we conclude that MP-WGS is a feasible method for routine leukemia diagnostics of structural chromosomal rearrangements.

N-alpha-acetylation is a common co-translational protein modification that is essential for normal cell function in humans. We previously identified the genetic basis of an X-linked infantile lethal Mendelian disorder involving a c.109T>C (p.Ser37Pro) missense variant in NAA10, which encodes the catalytic subunit of the N-terminal acetyltransferase A (NatA) complex. The auxiliary subunit of the NatA complex, NAA15, is the dimeric binding partner for NAA10. Through a genotype-first approach with whole-exome or genome sequencing (WES/WGS) and targeted sequencing analysis, we identified and phenotypically characterized 38 individuals from 33 unrelated families with 25 different de novo or inherited, dominantly acting likely gene disrupting (LGD) variants in NAA15. Clinical features of affected individuals with LGD variants in NAA15 include variable levels of intellectual disability, delayed speech and motor milestones, and autism spectrum disorder. Additionally, mild craniofacial dysmorphology, congenital cardiac anomalies, and seizures are present in some subjects. RNA analysis in cell lines from two individuals showed degradation of the transcripts with LGD variants, probably as a result of nonsense-mediated decay. Functional assays in yeast confirmed a deleterious effect for two of the LGD variants in NAA15. Further supporting a mechanism of haploinsufficiency, individuals with copy-number variant (CNV) deletions involving NAA15 and surrounding genes can present with mild intellectual disability, mild dysmorphic features, motor delays, and decreased growth. We propose that defects in NatA-mediated N-terminal acetylation (NTA) lead to variable levels of neurodevelopmental disorders in humans, supporting the importance of the NatA complex in normal human development.

In this study, we take advantage of human induced pluripotent stem (iPS) cell-derived neural stem cells and brain organoids to study the role of p53 during human brain development. We knocked down (KD) p53 in human neuroepithelial stem (NES) cells derived from iPS cells. Upon p53KD, NES cells rapidly show centrosome amplification and genomic instability. Furthermore, a reduced proliferation rate, downregulation of genes involved in oxidative phosphorylation (OXPHOS), and an upregulation of glycolytic capacity was apparent upon loss of p53. In addition, p53KD neural stem cells display an increased pace of differentiating into neurons and exhibit a phenotype corresponding to more mature neurons compared to control neurons. Using brain organoids, we modeled more specifically cortical neurogenesis. Here we found that p53 loss resulted in brain organoids with disorganized stem cell layer and reduced cortical progenitor cells and neurons. Similar to NES cells, neural progenitors isolated from brain organoids also show a downregulation in several OXPHOS genes. Taken together, this demonstrates an important role for p53 in controlling genomic stability of neural stem cells and regulation of neuronal differentiation, as well as maintaining structural organization and proper metabolic gene profile of neural progenitors in human brain organoids.

Gastrointestinal atresias and urological defects are main causes of pediatric surgery in infants. As copy number variants (CNVs) have been shown to be involved in the development of congenital malformations, the aim of our study was to investigate the presence of CNVs in patients with gastrointestinal and urological malformations as well as the possibility of tissue-specific mosaicism for CNVs in the cohort.

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.

Loss-of-function variants in ANKRD11 were identified as the cause of KBG syndrome, an autosomal dominant syndrome with specific dental, neurobehavioural, craniofacial and skeletal anomalies. We present the largest cohort of KBG syndrome cases confirmed by ANKRD11 variants reported so far, consisting of 20 patients from 13 families. Sixteen patients were molecularly diagnosed by Sanger sequencing of ANKRD11, one familial case and three sporadic patients were diagnosed through whole-exome sequencing and one patient was identified through genomewide array analysis. All patients were evaluated by a clinical geneticist. Detailed orofacial phenotyping, including orthodontic evaluation, intra-oral photographs and orthopantomograms, was performed in 10 patients and revealed besides the hallmark feature of macrodontia of central upper incisors, several additional dental anomalies as oligodontia, talon cusps and macrodontia of other teeth. Three-dimensional (3D) stereophotogrammetry was performed in 14 patients and 3D analysis of patients compared with controls showed consistent facial dysmorphisms comprising a bulbous nasal tip, upturned nose with a broad base and a round or triangular face. Many patients exhibited neurobehavioural problems, such as autism spectrum disorder or hyperactivity. One-third of patients presented with (conductive) hearing loss. Congenital heart defects, velopharyngeal insufficiency and hip anomalies were less frequent. On the basis of our observations, we recommend cardiac assessment in children and regular hearing tests in all individuals with a molecular diagnosis of KBG syndrome. As ANKRD11 is a relatively common gene in which sequence variants have been identified in individuals with neurodevelopmental disorders, it seems an important contributor to the aetiology of both sporadic and familial cases.

In order to identify genetic causes of VACTERL association (V vertebral defects, A anorectal malformations, C cardiac defects, T tracheoesofageal fistula, E esophageal atresia, R renal anomalies, L limb deformities), we have collected DNA samples from 20 patients diagnosed with VACTERL or with a VACTERL-like phenotype as well as samples from 19 aborted fetal cases with VACTERL. To investigate the importance of gene dose alterations in the genetic etiology of VACTERL association we have performed a systematic analysis of this cohort using a 180K array comparative genomic hybridization (array-CGH) platform. In addition, to further clarify the significance of PCSK5, HOXD13 and CHD7 genes in the VACTERL phenotype, mutation screening has been performed. We identified pathogenic gene dose imbalances in two fetal cases; a hemizygous deletion of the FANCB gene and a (9;18)(p24;q12) unbalanced translocation. In addition, one pathogenic mutation in CHD7 was detected, while no apparent disease-causing mutations were found in HOXD13 or PCSK5. Our study shows that although large gene dose alterations do not seem to be a common cause in VACTERL association, array-CGH is still important in clinical diagnostics to identify disease cause in individual cases.

Neu-Laxova syndrome (NLS) is a rare autosomal-recessive disorder characterized by a recognizable pattern of severe malformations leading to prenatal or early postnatal lethality. Homozygous mutations in PHGDH, a gene involved in the first and limiting step in L-serine biosynthesis, were recently identified as the cause of the disease in three families. By studying a cohort of 12 unrelated families affected by NLS, we provide evidence that NLS is genetically heterogeneous and can be caused by mutations in all three genes encoding enzymes of the L-serine biosynthesis pathway. Consistent with recently reported findings, we could identify PHGDH missense mutations in three unrelated families of our cohort. Furthermore, we mapped an overlapping homozygous chromosome 9 region containing PSAT1 in four consanguineous families. This gene encodes phosphoserine aminotransferase, the enzyme for the second step in L-serine biosynthesis. We identified six families with three different missense and frameshift PSAT1 mutations fully segregating with the disease. In another family, we discovered a homozygous frameshift mutation in PSPH, the gene encoding phosphoserine phosphatase, which catalyzes the last step of L-serine biosynthesis. Interestingly, all three identified genes have been previously implicated in serine-deficiency disorders, characterized by variable neurological manifestations. Our findings expand our understanding of NLS as a disorder of the L-serine biosynthesis pathway and suggest that NLS represents the severe end of serine-deficiency disorders, demonstrating that certain complex syndromes characterized by early lethality could indeed be the extreme end of the phenotypic spectrum of already known disorders.

We report the findings from 4437 individuals (3219 patients and 1218 relatives) who have been analyzed by whole genome sequencing (WGS) at the Genomic Medicine Center Karolinska-Rare Diseases (GMCK-RD) since mid-2015. GMCK-RD represents a long-term collaborative initiative between Karolinska University Hospital and Science for Life Laboratory to establish advanced, genomics-based diagnostics in the Stockholm healthcare setting.