De Novo and Inherited Loss-of-Function Variants in TLK2: Clinical and Genotype-Phenotype Evaluation of a Distinct Neurodevelopmental Disorder

De novo TLK1 and MDM1 mutations in a patient with a neurodevelopmental disorder and immunodeficiency

The Tousled-like kinases 1 and 2 (TLK1/TLK2) regulate DNA replication, repair and chromatin maintenance. TLK2 variants underlie the neurodevelopmental disorder (NDD) 'Intellectual Disability, Autosomal Dominant 57' (MRD57), characterized by intellectual disability and microcephaly. Several TLK1 variants have been reported in NDDs but their functional significance is unknown. A male patient presenting with ID, seizures, global developmental delay, hypothyroidism, and primary immunodeficiency was determined to have a heterozygous TLK1 variant (c.1435C>G, p.Q479E), as well as a mutation in MDM1 (c.1197dupT, p.K400∗). Cells expressing TLK1 p.Q479E exhibited reduced cytokine responses and elevated DNA damage, but not increased radiation sensitivity or DNA repair defects. The TLK1 p.Q479E variant impaired kinase activity but not proximal protein interactions. Our study provides the first functional characterization of NDD-associated TLK1 variants and suggests that, such as TLK2, TLK1 variants may impact development in multiple tissues and should be considered in the diagnosis of rare NDDs.

Discovery and optimization of narrow spectrum inhibitors of Tousled like kinase 2 (TLK2) using quantitative structure activity relationships

The oxindole scaffold has been the center of several kinase drug discovery programs, some of which have led to approved medicines. A series of two oxindole matched pairs from the literature were identified where TLK2 was potently inhibited as an off-target kinase. The oxindole has long been considered a promiscuous kinase inhibitor template, but across these four specific literature oxindoles TLK2 activity was consistent, while the kinome profile was radically different ranging from narrow to broad spectrum kinome coverage. We synthesized a large series of analogues, utilizing quantitative structure-activity relationship (QSAR) analysis, water mapping of the kinase ATP binding sites, kinome profiling, and small-molecule x-ray structural analysis to optimize TLK2 inhibition and kinome selectivity. This resulted in the identification of several narrow spectrum, sub-family selective, chemical tool compounds including 128 (UNC-CA2-103) that could enable elucidation of TLK2 biology.

Discovery and Optimization of Narrow Spectrum Inhibitors of Tousled Like Kinase 2 (TLK2) Using Quantitative Structure Activity Relationships

The oxindole scaffold has been the center of several kinase drug discovery programs, some of which have led to approved medicines. A series of two oxindole matched pairs from the literature were identified where TLK2 was a potent off-target kinase. The oxindole has long been considered a promiscuous inhibitor template, but across these 4 specific literature oxindoles TLK2 activity was consistent, while the kinome profile was radically different from narrow to broad spectrum coverage. We synthesized a large series of analogues and through quantitative structure-activity relationship (QSAR) analysis, water mapping of the kinase ATP binding sites, small-molecule x-ray structural analysis and kinome profiling, narrow spectrum, sub-family selective, chemical tool compounds were identified to enable elucidation of TLK2 biology.

Pathogenic variants in the paired-related homeobox 1 gene (PRRX1) cause craniosynostosis with incomplete penetrance

PURPOSE

Studies have previously implicated PRRX1 in craniofacial development, including demonstration of murine Prrx1 expression in the preosteogenic cells of the cranial sutures. We investigated the role of heterozygous missense and loss-of-function (LoF) variants in PRRX1 associated with craniosynostosis.

METHODS

Trio-based genome, exome, or targeted sequencing were used to screen PRRX1 in patients with craniosynostosis; immunofluorescence analyses were used to assess nuclear localization of wild-type and mutant proteins.

RESULTS

Genome sequencing identified 2 of 9 sporadically affected individuals with syndromic/multisuture craniosynostosis, who were heterozygous for rare/undescribed variants in PRRX1. Exome or targeted sequencing of PRRX1 revealed a further 9 of 1449 patients with craniosynostosis harboring deletions or rare heterozygous variants within the homeodomain. By collaboration, 7 additional individuals (4 families) were identified with putatively pathogenic PRRX1 variants. Immunofluorescence analyses showed that missense variants within the PRRX1 homeodomain cause abnormal nuclear localization. Of patients with variants considered likely pathogenic, bicoronal or other multisuture synostosis was present in 11 of 17 cases (65%). Pathogenic variants were inherited from unaffected relatives in many instances, yielding a 12.5% penetrance estimate for craniosynostosis.

CONCLUSION

This work supports a key role for PRRX1 in cranial suture development and shows that haploinsufficiency of PRRX1 is a relatively frequent cause of craniosynostosis.

Identification of a de novo mutation in TLK1 associated with a neurodevelopmental disorder and immunodeficiency

Background

The Tousled-like kinases 1 and 2 (TLK1/TLK2) regulate DNA replication, repair and chromatin maintenance. TLK2 variants are associated with 'Intellectual Disability, Autosomal Dominant 57' (MRD57), a neurodevelopmental disorder (NDD) characterized by intellectual disability (ID), autism spectrum disorder (ASD) and microcephaly. Several TLK1 variants have been reported in NDDs but their functional significance is unknown.

Methods

A male patient presenting with ID, seizures, global developmental delay, hypothyroidism, and primary immunodeficiency was determined to have a novel, heterozygous variant in TLK1 (c.1435C>G, p.Q479E) by genome sequencing (GS). Single cell gel electrophoresis, western blot, flow cytometry and RNA-seq were performed in patient-derived lymphoblast cell lines. In silico, biochemical and proteomic analysis were used to determine the functional impact of the p.Q479E variant and previously reported NDD-associated TLK1 variant, p.M566T.

Results

Transcriptome sequencing in patient-derived cells confirmed expression of TLK1 transcripts carrying the p.Q479E variant and revealed alterations in genes involved in class switch recombination and cytokine signaling. Cells expressing the p.Q479E variant exhibited reduced cytokine responses and higher levels of spontaneous DNA damage but not increased sensitivity to radiation or DNA repair defects. The p.Q479E and p.M566T variants impaired kinase activity but did not strongly alter localization or proximal protein interactions.

Conclusion

Our study provides the first functional characterization of TLK1 variants associated with NDDs and suggests potential involvement in central nervous system and immune system development. Our results indicate that, like TLK2 variants, TLK1 variants may impact development in multiple tissues and should be considered in the diagnosis of rare NDDs.

A cryptic microdeletion del(12)(p11.21p11.23) within an unbalanced translocation t(7;12)(q21.13;q23.1) implicates new candidate loci for intellectual disability and Kallmann syndrome

In a patient diagnosed with both Kallmann syndrome (KS) and intellectual disability (ID), who carried an apparently balanced translocation t(7;12)(q22;q24)dn, array comparative genomic hybridization (aCGH) disclosed a cryptic heterozygous 4.7 Mb deletion del(12)(p11.21p11.23), unrelated to the translocation breakpoint. This novel discovery prompted us to consider the possibility that the combination of KS and neurological disorder in this patient could be attributed to gene(s) within this specific deletion at 12p11.21-12p11.23, rather than disrupted or dysregulated genes at the translocation breakpoints. To further support this hypothesis, we expanded our study by screening five candidate genes at both breakpoints of the chromosomal translocation in a cohort of 48 KS patients. However, no mutations were found, thus reinforcing our supposition. In order to delve deeper into the characterization of the 12p11.21-12p11.23 region, we enlisted six additional patients with small copy number variations (CNVs) and analyzed eight individuals carrying small CNVs in this region from the DECIPHER database. Our investigation utilized a combination of complementary approaches. Firstly, we conducted a comprehensive phenotypic-genotypic comparison of reported CNV cases. Additionally, we reviewed knockout animal models that exhibit phenotypic similarities to human conditions. Moreover, we analyzed reported variants in candidate genes and explored their association with corresponding phenotypes. Lastly, we examined the interacting genes associated with these phenotypes to gain further insights. As a result, we identified a dozen candidate genes: TSPAN11 as a potential KS candidate gene, TM7SF3, STK38L, ARNTL2, ERGIC2, TMTC1, DENND5B, and ETFBKMT as candidate genes for the neurodevelopmental disorder, and INTS13, REP15, PPFIBP1, and FAR2 as candidate genes for KS with ID. Notably, the high-level expression pattern of these genes in relevant human tissues further supported their candidacy. Based on our findings, we propose that dosage alterations of these candidate genes may contribute to sexual and/or cognitive impairments observed in patients with KS and/or ID. However, the confirmation of their causal roles necessitates further identification of point mutations in these candidate genes through next-generation sequencing.

The Genetic Landscape of Children Born Small for Gestational Age with Persistent Short Stature

INTRODUCTION

Among children born small for gestational age, 10-15% fail to catch up and remain short (SGA-SS). The underlying mechanisms are mostly unknown. We aimed to decipher genetic aetiologies of SGA-SS within a large single-centre cohort.

METHODS

Out of 820 patients treated with growth hormone (GH), 256 were classified as SGA-SS (birth length and/or birth weight <-2 SD for gestational age and life-minimum height <-2.5 SD). Those with the DNA triplet available (child and both parents) were included in the study (176/256). Targeted testing (karyotype/FISH/MLPA/specific Sanger sequencing) was performed if a specific genetic disorder was clinically suggestive. All remaining patients underwent MS-MLPA to identify Silver-Russell syndrome, and those with unknown genetic aetiology were subsequently examined using whole-exome sequencing or targeted panel of 398 growth-related genes. Genetic variants were classified using ACMG guidelines.

RESULTS

The genetic aetiology was elucidated in 74/176 (42%) children. Of these, 12/74 (16%) had pathogenic or likely pathogenic (P/LP) gene variants affecting pituitary development (LHX4, OTX2, PROKR2, PTCH1, POU1F1), the GH-IGF-1 or IGF-2 axis (GHSR, IGFALS, IGF1R, STAT3, HMGA2), 2/74 (3%) the thyroid axis (TRHR, THRA), 17/74 (23%) the cartilaginous matrix (ACAN, various collagens, FLNB, MATN3), and 7/74 (9%) the paracrine chondrocyte regulation (FGFR3, FGFR2, NPR2). In 12/74 (16%), we revealed P/LP affecting fundamental intracellular/intranuclear processes (CDC42, KMT2D, LMNA, NSD1, PTPN11, SRCAP, SON, SOS1, SOX9, TLK2). SHOX deficiency was found in 7/74 (9%), Silver-Russell syndrome in 12/74 (16%) (11p15, UPD7), and miscellaneous chromosomal aberrations in 5/74 (7%) children.

CONCLUSIONS

The high diagnostic yield sheds a new light on the genetic landscape of SGA-SS, with a central role for the growth plate with substantial contributions from the GH-IGF-1 and thyroid axes and intracellular regulation and signalling.

A microdeletion del(12)(p11.21p11.23) with a cryptic unbalanced translocation t(7;12)(q21.13;q23.1) implicates new candidate loci for intellectual disability and Kallmann syndrome

In an apparently balanced translocation t(7;12)(q22;q24)dn exhibiting both Kallmann syndrome (KS) and intellectual disability (ID), we detected a cryptic heterozygous 4.7 Mb del(12)(p11.21p11.23) unrelated to the translocation breakpoint. This new finding raised the possibility that KS combined with neurological disorder in this patient could be caused by gene(s) within this deletion at 12p11.21-12p11.23 instead of disrupted or dysregulated genes at the genomic breakpoints. Screening of five candidate genes at both breakpoints in 48 KS patients we recruited found no mutation, corroborating our supposition. To substantiate this hypothesis further, we recruited six additional subjects with small CNVs and analyzed eight individuals carrying small CNVs in this region from DECIPHER to dissect 12p11.21-12p11.23. We used multiple complementary approaches including a phenotypic-genotypic comparison of reported cases, a review of knockout animal models recapitulating the human phenotypes, and analyses of reported variants in the interacting genes with corresponding phenotypes. The results identified one potential KS candidate gene ( TSPAN11 ), seven candidate genes for the neurodevelopmental disorder ( TM7SF3 , STK38L , ARNTL2 , ERGIC2 , TMTC1 , DENND5B , and ETFBKMT ), and four candidate genes for KS with ID ( INTS13 , REP15 , PPFIBP1 , and FAR2 ). The high-level expression pattern in the relevant human tissues further suggested the candidacy of these genes. We propose that the dosage alterations of the candidate genes may contribute to sexual and/or cognitive impairment in patients with KS and/or ID. Further identification of point mutations through next generation sequencing will be necessary to confirm their causal roles.

Comprehensive analysis of transcriptome characteristics and identification of TLK2 as a potential biomarker in dermatofibrosarcoma protuberans

Background: Dermatofibrosarcoma protuberans (DFSP) is a rare cutaneous sarcoma characterized by local invasion and recurrence. RNA sequencing (RNA-seq) allows the qualification of cellular RNA populations and provides information on the transcriptional state. However, few studies have comprehensively analyzed DFSP transcriptional data.

Methods: Fourteen DFSP samples with paired non-neoplastic soft tissue from Chinese patients undergoing Mohs micrographic surgery were used for RNA-seq analysis. Differential expression analysis and enrichment analysis for RNA-seq data were performed to identify fusion genes, biomarkers, and microenvironment characteristics of DFSP.

Results: This study systemically describes the transcriptomic characteristics of DFSP. First, we performed gene fusion analysis and identified a novel FBN1-CSAD fusion event in a DFSP patient with fibrosarcomatous transformation. Then, we identified TLK2 as a biomarker for DFSP based on functional enrichment analysis, and validated its accuracy for diagnosing DFSP by immunohistochemical staining and joint analysis with public data. Finally, microenvironment analysis described the infiltration characteristics of immune and stromal cells in DFSP.

Conclusion: This study demonstrates that RNA-seq can serve as a promising strategy for exploring molecular mechanisms in DFSP. Our results provide new insights into accurate diagnosis and therapeutic targets of DFSP.

A Diffusion Tensor Imaging Analysis of Frontal Lobe White Matter Microstructure in Trigonocephaly Patients

BACKGROUND

Children with trigonocephaly are at risk for neurodevelopmental disorders. The aim of this study is to investigate white matter properties of the frontal lobes in young, unoperated patients with metopic synostosis as compared to healthy controls using diffusion tension imaging (DTI).

METHODS

Preoperative DTI data sets of 46 patients with trigonocephaly with a median age of 0.49 (interquartile range: 0.38) years were compared with 21 controls with a median age of 1.44 (0.98) years. White matter metrics of the tracts in the frontal lobe were calculated using FMRIB Software Library (FSL). The mean value of tract-specific fractional anisotropy (FA) and mean diffusivity (MD) were estimated for each subject and compared to healthy controls. By linear regression, FA and MD values per tract were assessed by trigonocephaly, sex, and age.

RESULTS

The mean FA and MD values in the frontal lobe tracts of untreated trigonocephaly patients, younger than 3 years, were not significantly different in comparison to controls, where age showed to be a significant associated factor.

CONCLUSIONS

Microstructural parameters of white matter tracts of the frontal lobe of patients with trigonocephaly are comparable to those of controls aged 0-3 years.

Inherited variants in CHD3 show variable expressivity in Snijders Blok-Campeau syndrome

PURPOSE

Common diagnostic next-generation sequencing strategies are not optimized to identify inherited variants in genes associated with dominant neurodevelopmental disorders as causal when the transmitting parent is clinically unaffected, leaving a significant number of cases with neurodevelopmental disorders undiagnosed.

METHODS

We characterized 21 families with inherited heterozygous missense or protein-truncating variants in CHD3, a gene in which de novo variants cause Snijders Blok-Campeau syndrome.

RESULTS

Computational facial and Human Phenotype Ontology-based comparisons showed that the phenotype of probands with inherited CHD3 variants overlaps with the phenotype previously associated with de novo CHD3 variants, whereas heterozygote parents are mildly or not affected, suggesting variable expressivity. In addition, similarly reduced expression levels of CHD3 protein in cells of an affected proband and of healthy family members with a CHD3 protein-truncating variant suggested that compensation of expression from the wild-type allele is unlikely to be an underlying mechanism. Notably, most inherited CHD3 variants were maternally transmitted.

CONCLUSION

Our results point to a significant role of inherited variation in Snijders Blok-Campeau syndrome, a finding that is critical for correct variant interpretation and genetic counseling and warrants further investigation toward understanding the broader contributions of such variation to the landscape of human disease.

Report of two children with global developmental delay in association with de novo TLK2 variant and literature review

We describe clinical details, including novel findings, of two further children with the newly defined TLK2-related disorder. One patient was recruited to the Deciphering Developmental Delay (DDD) Study to identify underlying etiology of global developmental delay. The other was detected on whole-exome sequencing as part of second line investigations following normal microarray. Both patients were found to have de novo heterozygous pathogenic TLK2 variants. A novel c.6del p.(Glu3Lysfs*) loss-of-function frameshift variant was found in Patient 1. A c.1121+1G>A splice-donor variant was detected in Patient 2. TLK2-related neurodevelopmental disorder is a specific syndrome that has been recently described. Global developmental delay, behavioral problems, gastrointestinal disorders, and typical facial dysmorphism are common features. Neuropsychiatric disorders, ophthalmic, musculoskeletal and cranial abnormalities, as well as short stature, have also all been described. The novel findings we describe include sleep disturbance, nondifferentiation of lateral semi-circular canals (where asymmetric semi-circular canals were a feature in the previous cohort), vesico-ureteric reflux, and bilateral periauricular skin tags. Here, we report a novel TLK2 variant and previously undescribed features of TLK2-related disorder, to expand the clinical phenotype and provide further genotype-phenotype correlation.

The benefit of diagnostic whole genome sequencing in schizophrenia and other psychotic disorders

Schizophrenia has a multifactorial etiology, involving a polygenic architecture. The potential benefit of whole genome sequencing (WGS) in schizophrenia and other psychotic disorders is not well studied. We investigated the yield of clinical WGS analysis in 251 families with a proband diagnosed with schizophrenia (N = 190), schizoaffective disorder (N = 49), or other conditions involving psychosis (N = 48). Participants were recruited in Israel and USA, mainly of Jewish, Arab, and other European ancestries. Trio (parents and proband) WGS was performed for 228 families (90.8%); in the other families, WGS included parents and at least two affected siblings. In the secondary analyses, we evaluated the contribution of rare variant enrichment in particular gene sets, and calculated polygenic risk score (PRS) for schizophrenia. For the primary outcome, diagnostic rate was 6.4%; we found clinically significant, single nucleotide variants (SNVs) or small insertions or deletions (indels) in 14 probands (5.6%), and copy number variants (CNVs) in 2 (0.8%). Significant enrichment of rare loss-of-function variants was observed in a gene set of top schizophrenia candidate genes in affected individuals, compared with population controls (N = 6,840). The PRS for schizophrenia was significantly increased in the affected individuals group, compared to their unaffected relatives. Last, we were also able to provide pharmacogenomics information based on CYP2D6 genotype data for most participants, and determine their antipsychotic metabolizer status. In conclusion, our findings suggest that WGS may have a role in the setting of both research and genetic counseling for individuals with schizophrenia and other psychotic disorders and their families.

Cerebral Blood Flow of the Frontal Lobe in Untreated Children with Trigonocephaly versus Healthy Controls: An Arterial Spin Labeling Study

BACKGROUND

Craniofacial surgery is the standard treatment for children with moderate to severe trigonocephaly. The added value of surgery to release restriction of the frontal lobes is unproven, however. In this study, the authors aim to address the hypothesis that the frontal lobe perfusion is not restricted in trigonocephaly patients by investigating cerebral blood flow.

METHODS

Between 2018 and 2020, trigonocephaly patients for whom a surgical correction was considered underwent magnetic resonance imaging brain studies with arterial spin labeling to measure cerebral perfusion. The mean value of cerebral blood flow in the frontal lobe was calculated for each subject and compared to that of healthy controls.

RESULTS

Magnetic resonance imaging scans of 36 trigonocephaly patients (median age, 0.5 years; interquartile range, 0.3; 11 female patients) were included and compared to those of 16 controls (median age, 0.83 years; interquartile range, 0.56; 10 female patients). The mean cerebral blood flow values in the frontal lobe of the trigonocephaly patients (73.0 ml/100 g/min; SE, 2.97 ml/100 g/min) were not significantly different in comparison to control values (70.5 ml/100 g/min; SE, 4.45 ml/100 g/min; p = 0.65). The superior, middle, and inferior gyri of the frontal lobe showed no significant differences either.

CONCLUSIONS

The authors' findings suggest that the frontal lobes of trigonocephaly patients aged less than 18 months have a normal cerebral blood flow before surgery. In addition to the very low prevalence of papilledema or impaired skull growth previously reported, this finding further supports the authors' hypothesis that craniofacial surgery for trigonocephaly is rarely indicated for signs of raised intracranial pressure or restricted perfusion for patients younger than 18 months.

CLINICAL QUESTION/LEVEL OF EVIDENCE

Risk, II.

Tousled-like kinase 2 targets ASF1 histone chaperones through client mimicry

Tousled-like kinases (TLKs) are nuclear serine-threonine kinases essential for genome maintenance and proper cell division in animals and plants. A major function of TLKs is to phosphorylate the histone chaperone proteins ASF1a and ASF1b to facilitate DNA replication-coupled nucleosome assembly, but how TLKs selectively target these critical substrates is unknown. Here, we show that TLK2 selectivity towards ASF1 substrates is achieved in two ways. First, the TLK2 catalytic domain recognizes consensus phosphorylation site motifs in the ASF1 C-terminal tail. Second, a short sequence at the TLK2 N-terminus docks onto the ASF1a globular N-terminal domain in a manner that mimics its histone H3 client. Disrupting either catalytic or non-catalytic interactions through mutagenesis hampers ASF1 phosphorylation by TLK2 and cell growth. Our results suggest that the stringent selectivity of TLKs for ASF1 is enforced by an unusual interaction mode involving mutual recognition of a short sequence motifs by both kinase and substrate.

Tousled-like kinase 2 (TLK2) phosphorylates ASF1 histone chaperones to promote nucleosome assembly in S phase. Here, the authors show that TLK2 targets ASF1 by simulating its client protein histone H3, exploiting a primordial protein interaction surface for regulatory control.

Functional analysis of TLK2 variants and their proximal interactomes implicates impaired kinase activity and chromatin maintenance defects in their pathogenesis

INTRODUCTION

The Tousled-like kinases 1 and 2 (TLK1 and TLK2) are involved in many fundamental processes, including DNA replication, cell cycle checkpoint recovery and chromatin remodelling. Mutations in TLK2 were recently associated with 'Mental Retardation Autosomal Dominant 57' (MRD57, MIM# 618050), a neurodevelopmental disorder characterised by a highly variable phenotype, including mild-to-moderate intellectual disability, behavioural abnormalities, facial dysmorphisms, microcephaly, epilepsy and skeletal anomalies.

METHODS

We re-evaluate whole exome sequencing and array-CGH data from a large cohort of patients affected by neurodevelopmental disorders. Using spatial proteomics (BioID) and single-cell gel electrophoresis, we investigated the proximity interaction landscape of TLK2 and analysed the effects of p.(Asp551Gly) and a previously reported missense variant (c.1850C>T; p.(Ser617Leu)) on TLK2 interactions, localisation and activity.

RESULTS

We identified three new unrelated MRD57 families. Two were sporadic and caused by a missense change (c.1652A>G; p.(Asp551Gly)) or a 39 kb deletion encompassing TLK2, and one was familial with three affected siblings who inherited a nonsense change from an affected mother (c.1423G>T; p.(Glu475Ter)). The clinical phenotypes were consistent with those of previously reported cases. The tested mutations strongly impaired TLK2 kinase activity. Proximal interactions between TLK2 and other factors implicated in neurological disorders, including CHD7, CHD8, BRD4 and NACC1, were identified. Finally, we demonstrated a more relaxed chromatin state in lymphoblastoid cells harbouring the p.(Asp551Gly) variant compared with control cells, conferring susceptibility to DNA damage.

CONCLUSION

Our study identified novel TLK2 pathogenic variants, confirming and further expanding the MRD57-related phenotype. The molecular characterisation of missense variants increases our knowledge about TLK2 function and provides new insights into its role in neurodevelopmental disorders.

Family-Based Genome-Wide Association Study of Autism Spectrum Disorder in Middle Eastern Families

Autism spectrum disorder (ASD) is a neurodevelopmental disease characterized by abnormalities in language and social communication with substantial clinical heterogeneity. Genetic factors play an important role in ASD with heritability estimated between 70% to 80%. Genome-wide association studies (GWAS) have identified multiple loci associated with ASD. However, most studies were performed on European populations and little is known about the genetic architecture of ASD in Middle Eastern populations. Here, we report the first GWAS of ASD in the Middle eastern population of Qatar. We analyzed 171 families with ASD, using linear mixed models adjusting for relatedness and other confounders. Results showed that common single nucleotide polymorphisms (SNP) in seven loci are associated with ASD (p < 1 × 10-5). Although the identified loci did not reach genome-wide significance, many of the top associated SNPs are located within or near genes that have been implicated in ASD or related neurodevelopmental disorders. These include GORASP2, GABBR2, ANKS6, THSD4, ERCC6L, ARHGEF6, and HDAC8. Additionally, three of the top associated SNPs were significantly associated with gene expression. We also found evidence of association signals in two previously reported ASD-susceptibility loci (rs10099100 and rs4299400). Our results warrant further functional studies and replication to provide further insights into the genetic architecture of ASD.

Tousled-Like Kinases Suppress Innate Immune Signaling Triggered by Alternative Lengthening of Telomeres

The Tousled-like kinases 1 and 2 (TLK1/2) control histone deposition through the ASF1 histone chaperone and influence cell cycle progression and genome maintenance, yet the mechanisms underlying TLK-mediated genome stability remain uncertain. Here, we show that TLK loss results in severe chromatin decompaction and altered genome accessibility, particularly affecting heterochromatic regions. Failure to maintain heterochromatin increases spurious transcription of repetitive elements and induces features of alternative lengthening of telomeres (ALT). TLK depletion culminates in a cGAS-STING-TBK1-mediated innate immune response that is independent of replication-stress signaling and attenuated by the depletion of factors required to produce extra-telomeric DNA. Analysis of human cancers reveals that chromosomal instability correlates with high TLK2 and low STING levels in many cohorts. Based on these findings, we propose that high TLK levels contribute to immune evasion in chromosomally unstable and ALT+ cancers.

•

TLK-deficient cells have increased accessibility at heterochromatin regions

•

TLK1/2 suppress spurious transcription and telomere hyper-recombination

•

Extra-telomeric DNA generated upon TLK loss promotes innate immune signaling

•

cGAS-STING-TBK1 signaling in TLK-deficient cells is independent of replication stress

A DNA repair disorder caused by de novo monoallelic DDB1 variants is associated with a neurodevelopmental syndrome

The DNA damage-binding protein 1 (DDB1) is part of the CUL4-DDB1 ubiquitin E3 ligase complex (CRL4), which is essential for DNA repair, chromatin remodeling, DNA replication, and signal transduction. Loss-of-function variants in genes encoding the complex components CUL4 and PHIP have been reported to cause syndromic intellectual disability with hypotonia and obesity, but no phenotype has been reported in association with DDB1 variants. Here, we report eight unrelated individuals, identified through Matchmaker Exchange, with de novo monoallelic variants in DDB1, including one recurrent variant in four individuals. The affected individuals have a consistent phenotype of hypotonia, mild to moderate intellectual disability, and similar facies, including horizontal or slightly bowed eyebrows, deep-set eyes, full cheeks, a short nose, and large, fleshy and forward-facing earlobes, demonstrated in the composite face generated from the cohort. Digital anomalies, including brachydactyly and syndactyly, were common. Three older individuals have obesity. We show that cells derived from affected individuals have altered DDB1 function resulting in abnormal DNA damage signatures and histone methylation following UV-induced DNA damage. Overall, our study adds to the growing family of neurodevelopmental phenotypes mediated by disruption of the CRL4 ubiquitin ligase pathway and begins to delineate the phenotypic and molecular effects of DDB1 misregulation.

Mutation-specific pathophysiological mechanisms define different neurodevelopmental disorders associated with SATB1 dysfunction

Whereas large-scale statistical analyses can robustly identify disease-gene relationships, they do not accurately capture genotype-phenotype correlations or disease mechanisms. We use multiple lines of independent evidence to show that different variant types in a single gene, SATB1, cause clinically overlapping but distinct neurodevelopmental disorders. Clinical evaluation of 42 individuals carrying SATB1 variants identified overt genotype-phenotype relationships, associated with different pathophysiological mechanisms, established by functional assays. Missense variants in the CUT1 and CUT2 DNA-binding domains result in stronger chromatin binding, increased transcriptional repression, and a severe phenotype. In contrast, variants predicted to result in haploinsufficiency are associated with a milder clinical presentation. A similarly mild phenotype is observed for individuals with premature protein truncating variants that escape nonsense-mediated decay, which are transcriptionally active but mislocalized in the cell. Our results suggest that in-depth mutation-specific genotype-phenotype studies are essential to capture full disease complexity and to explain phenotypic variability.

Dosage-sensitive functions in embryonic development drove the survival of genes on sex-specific chromosomes in snakes, birds, and mammals

Different ancestral autosomes independently evolved into sex chromosomes in snakes, birds, and mammals. In snakes and birds, females are ZW and males are ZZ; in mammals, females are XX and males are XY. Although X and Z Chromosomes retain nearly all ancestral genes, sex-specific W and Y Chromosomes suffered extensive genetic decay. In both birds and mammals, the genes that survived on sex-specific chromosomes are enriched for broadly expressed, dosage-sensitive regulators of gene expression, subject to strong purifying selection. To gain deeper insight into the processes that govern survival on sex-specific chromosomes, we carried out a meta-analysis of survival across 41 species-three snakes, 24 birds, and 14 mammals-doubling the number of ancestral genes under investigation and increasing our power to detect enrichments among survivors relative to nonsurvivors. Of 2564 ancestral genes, representing an eighth of the ancestral amniote genome, only 324 survive on present-day sex-specific chromosomes. Survivors are enriched for dosage-sensitive developmental processes, particularly development of neural crest-derived structures, such as the face. However, there was no enrichment for expression in sex-specific tissues, involvement in sex determination or gonadogenesis pathways, or conserved sex-biased expression. Broad expression and dosage sensitivity contributed independently to gene survival, suggesting that pleiotropy imposes additional constraints on the evolution of dosage compensation. We propose that maintaining the viability of the heterogametic sex drove gene survival on amniote sex-specific chromosomes, and that subtle modulation of the expression of survivor genes and their autosomal orthologs has disproportionately large effects on development and disease.

amplimap: a versatile tool to process and analyze targeted NGS data

SUMMARY

amplimap is a command-line tool to automate the processing and analysis of data from targeted next-generation sequencing experiments with PCR-based amplicons or capture-based enrichment systems. From raw sequencing reads, amplimap generates output such as read alignments, annotated variant calls, target coverage statistics and variant allele counts and frequencies for each target base pair. In addition to its focus on user-friendliness and reproducibility, amplimap supports advanced features such as consensus base calling for read families based on unique molecular identifiers and filtering false positive variant calls caused by amplification of off-target loci.

AVAILABILITY AND IMPLEMENTATION

amplimap is available as a free Python package under the open-source Apache 2.0 License. Documentation, source code and installation instructions are available at https://github.com/koelling/amplimap.

Molecular causes of primary microcephaly and related diseases: a report from the UNIA Workshop

The International University of Andalucía (UNIA) Current Trends in Biomedicine Workshop on Molecular Causes of Primary Microcephaly and Related Diseases took place in Baeza, Spain, November 18-20, 2019. This meeting brought together scientists from Europe, the USA and China to discuss recent advances in our molecular and genetic understanding of a group of rare neurodevelopmental diseases characterised by primary microcephaly, a condition in which head circumference is smaller than normal at birth. Microcephaly can be caused by inherited mutations that affect key cellular processes, or environmental exposure to radiation or other toxins. It can also result from viral infection, as exemplified by the recent Zika virus outbreak in South America. Here we summarise a number of the scientific advances presented and topics discussed at the meeting.

Severe neurodevelopmental disease caused by a homozygous TLK2 variant

A distinct neurodevelopmental phenotype characterised mainly by mild motor and language delay and facial dysmorphism, caused by heterozygous de novo or dominant variants in the TLK2 gene has recently been described. All cases reported carried either truncating variants located throughout the gene, or missense changes principally located at the C-terminal end of the protein mostly resulting in haploinsufficiency of TLK2. Through whole exome sequencing, we identified a homozygous missense variant in TLK2 in a patient showing more severe symptoms than those previously described, including cerebellar vermis hypoplasia and West syndrome. Both parents are heterozygous for the variant and clinically unaffected highlighting that recessive variants in TLK2 can also be disease causing and may act through a different pathomechanism.

The Tousled-like kinases regulate genome and epigenome stability: implications in development and disease

The Tousled-like kinases (TLKs) are an evolutionarily conserved family of serine-threonine kinases that have been implicated in DNA replication, DNA repair, transcription, chromatin structure, viral latency, cell cycle checkpoint control and chromosomal stability in various organisms. The functions of the TLKs appear to depend largely on their ability to regulate the H3/H4 histone chaperone ASF1, although numerous TLK substrates have been proposed. Over the last few years, a clearer picture of TLK function has emerged through the identification of new partners, the definition of specific roles in development and the elucidation of their structural and biochemical properties. In addition, the TLKs have been clearly linked to human disease; both TLK1 and TLK2 are frequently amplified in human cancers and TLK2 mutations have been identified in patients with neurodevelopmental disorders characterized by intellectual disability (ID), autism spectrum disorder (ASD) and microcephaly. A better understanding of the substrates, regulation and diverse roles of the TLKs is needed to understand their functions in neurodevelopment and determine if they are viable targets for cancer therapy. In this review, we will summarize current knowledge of TLK biology and its potential implications in development and disease.

Enabling Global Clinical Collaborations on Identifiable Patient Data: The Minerva Initiative

The clinical utility of computational phenotyping for both genetic and rare diseases is increasingly appreciated; however, its true potential is yet to be fully realized. Alongside the growing clinical and research availability of sequencing technologies, precise deep and scalable phenotyping is required to serve unmet need in genetic and rare diseases. To improve the lives of individuals affected with rare diseases through deep phenotyping, global big data interrogation is necessary to aid our understanding of disease biology, assist diagnosis, and develop targeted treatment strategies. This includes the application of cutting-edge machine learning methods to image data. As with most digital tools employed in health care, there are ethical and data governance challenges associated with using identifiable personal image data. There are also risks with failing to deliver on the patient benefits of these new technologies, the biggest of which is posed by data siloing. The Minerva Initiative has been designed to enable the public good of deep phenotyping while mitigating these ethical risks. Its open structure, enabling collaboration and data sharing between individuals, clinicians, researchers and private enterprise, is key for delivering precision public health.

Mutations in the Neuronal Vesicular SNARE VAMP2 Affect Synaptic Membrane Fusion and Impair Human Neurodevelopment

VAMP2 encodes the vesicular SNARE protein VAMP2 (also called synaptobrevin-2). Together with its partners syntaxin-1A and synaptosomal-associated protein 25 (SNAP25), VAMP2 mediates fusion of synaptic vesicles to release neurotransmitters. VAMP2 is essential for vesicular exocytosis and activity-dependent neurotransmitter release. Here, we report five heterozygous de novo mutations in VAMP2 in unrelated individuals presenting with a neurodevelopmental disorder characterized by axial hypotonia (which had been present since birth), intellectual disability, and autistic features. In total, we identified two single-amino-acid deletions and three non-synonymous variants affecting conserved residues within the C terminus of the VAMP2 SNARE motif. Affected individuals carrying de novo non-synonymous variants involving the C-terminal region presented a more severe phenotype with additional neurological features, including central visual impairment, hyperkinetic movement disorder, and epilepsy or electroencephalography abnormalities. Reconstituted fusion involving a lipid-mixing assay indicated impairment in vesicle fusion as one of the possible associated disease mechanisms. The genetic synaptopathy caused by VAMP2 de novo mutations highlights the key roles of this gene in human brain development and function.

Whole genome paired-end sequencing elucidates functional and phenotypic consequences of balanced chromosomal rearrangement in patients with developmental disorders

BACKGROUND

Balanced chromosomal rearrangements associated with abnormal phenotype are rare events, but may be challenging for genetic counselling, since molecular characterisation of breakpoints is not performed routinely. We used next-generation sequencing to characterise breakpoints of balanced chromosomal rearrangements at the molecular level in patients with intellectual disability and/or congenital anomalies.

METHODS

Breakpoints were characterised by a paired-end low depth whole genome sequencing (WGS) strategy and validated by Sanger sequencing. Expression study of disrupted and neighbouring genes was performed by RT-qPCR from blood or lymphoblastoid cell line RNA.

RESULTS

Among the 55 patients included (41 reciprocal translocations, 4 inversions, 2 insertions and 8 complex chromosomal rearrangements), we were able to detect 89% of chromosomal rearrangements (49/55). Molecular signatures at the breakpoints suggested that DNA breaks arose randomly and that there was no major influence of repeated elements. Non-homologous end-joining appeared as the main mechanism of repair (55% of rearrangements). A diagnosis could be established in 22/49 patients (44.8%), 15 by gene disruption (KANSL1, FOXP1, SPRED1, TLK2, MBD5, DMD, AUTS2, MEIS2, MEF2C, NRXN1, NFIX, SYNGAP1, GHR, ZMIZ1) and 7 by position effect (DLX5, MEF2C, BCL11B, SATB2, ZMIZ1). In addition, 16 new candidate genes were identified. Systematic gene expression studies further supported these results. We also showed the contribution of topologically associated domain maps to WGS data interpretation.

CONCLUSION

Paired-end WGS is a valid strategy and may be used for structural variation characterisation in a clinical setting.

Host and Microbiome Genome-Wide Association Studies: Current State and Challenges

The involvement of the microbiome in health and disease is well established. Microbiome genome-wide association studies (mGWAS) are used to elucidate the interaction of host genetic variation with the microbiome. The emergence of this relatively new field has been facilitated by the advent of next generation sequencing technologies that enable the investigation of the complex interaction between host genetics and microbial communities. In this paper, we review recent studies investigating host-microbiome interactions using mGWAS. Additionally, we highlight the marked disparity in the sampling population of mGWAS carried out to date and draw attention to the critical need for inclusion of diverse populations.