Representation of Protein Dynamics Disentangled by Time-Structure-Based Prior.

Representation learning (RL) is a universal technique for deriving low-dimensional disentangled representations from high-dimensional observations, aiding in a multitude of downstream tasks. RL has been extensively applied to various data types, including images and natural language. Here, we analyze molecular dynamics (MD) simulation data of biomolecules in terms of RL. Currently, state-of-the-art RL techniques, mainly motivated by the variational principle, try to capture slow motions in the representation (latent) space. Here, we propose two methods based on an alternative perspective on the disentanglement in the latent space. By disentanglement, we here mean the separation of underlying factors in the simulation data, aiding in detecting physically important coordinates for conformational transitions. The proposed methods introduce a simple prior that imposes temporal constraints in the latent space, serving as a regularization term to facilitate the capture of disentangled representations of dynamics. Comparison with other methods via the analysis of MD simulation trajectories for alanine dipeptide and chignolin validates that the proposed methods construct Markov state models (MSMs) whose implied time scales are comparable to those of the state-of-the-art methods. Using a measure based on total variation, we quantitatively evaluated that the proposed methods successfully disentangle physically important coordinates, aiding the interpretation of folding/unfolding transitions of chignolin. Overall, our methods provide good representations of complex biomolecular dynamics for downstream tasks, allowing for better interpretations of the conformational transitions.

Effective mechanical potential of cell-cell interaction explains three-dimensional morphologies during early embryogenesis.

Mechanical forces are critical for the emergence of diverse three-dimensional morphologies of multicellular systems. However, it remains unclear what kind of mechanical parameters at cellular level substantially contribute to tissue morphologies. This is largely due to technical limitations of live measurements of cellular forces. Here we developed a framework for inferring and modeling mechanical forces of cell-cell interactions. First, by analogy to coarse-grained models in molecular and colloidal sciences, we approximated cells as particles, where mean forces (i.e. effective forces) of pairwise cell-cell interactions are considered. Then, the forces were statistically inferred by fitting the mathematical model to cell tracking data. This method was validated by using synthetic cell tracking data resembling various in vivo situations. Application of our method to the cells in the early embryos of mice and the nematode Caenorhabditis elegans revealed that cell-cell interaction forces can be written as a pairwise potential energy in a manner dependent on cell-cell distances. Importantly, the profiles of the pairwise potentials were quantitatively different among species and embryonic stages, and the quantitative differences correctly described the differences of their morphological features such as spherical vs. distorted cell aggregates, and tightly vs. non-tightly assembled aggregates. We conclude that the effective pairwise potential of cell-cell interactions is a live measurable parameter whose quantitative differences can be a parameter describing three-dimensional tissue morphologies.

Molecular diagnosis of 405 individuals with autism spectrum disorder.

Autism spectrum disorder (ASD) is caused by combined genetic and environmental factors. Genetic heritability in ASD is estimated as 60-90%, and genetic investigations have revealed many monogenic factors. We analyzed 405 patients with ASD using family-based exome sequencing to detect disease-causing single-nucleotide variants (SNVs), small insertions and deletions (indels), and copy number variations (CNVs) for molecular diagnoses. All candidate variants were validated by Sanger sequencing or quantitative polymerase chain reaction and were evaluated using the American College of Medical Genetics and Genomics/Association for Molecular Pathology guidelines for molecular diagnosis. We identified 55 disease-causing SNVs/indels in 53 affected individuals and 13 disease-causing CNVs in 13 affected individuals, achieving a molecular diagnosis in 66 of 405 affected individuals (16.3%). Among the 55 disease-causing SNVs/indels, 51 occurred de novo, 2 were compound heterozygous (in one patient), and 2 were X-linked hemizygous variants inherited from unaffected mothers. The molecular diagnosis rate in females was significantly higher than that in males. We analyzed affected sibling cases of 24 quads and 2 quintets, but only one pair of siblings shared an identical pathogenic variant. Notably, there was a higher molecular diagnostic rate in simplex cases than in multiplex families. Our simulation indicated that the diagnostic yield is increasing by 0.63% (range 0-2.5%) per year. Based on our simple simulation, diagnostic yield is improving over time. Thus, periodical reevaluation of ES data should be strongly encouraged in undiagnosed ASD patients.

Identification of pathogenic deep intronic variant and exonic LINE-1 insertion in a patient with Meckel syndrome.

Biallelic CC2D2A variants are associated with a wide range of neurodevelopmental disorders including Meckel syndrome. Here we report a Japanese girl with Meckel syndrome harboring a pathogenic deep intronic variant (NM_001378615.1:c.1149+3569A>G) and an exonic LINE-1 insertion, which was predicted to cause aberrant splicing by SpliceAI and was detected by TEMP2 program, respectively. RNA analysis using urine-derived cells (UDCs) showed retention of 149-bp intronic sequences, leading to frameshift. Immunoblotting showed marked reduction of CC2D2A protein in the patient. Our report demonstrated that utilization of transposon detection tool and functional analysis using UDCs will increase diagnostic yield of genome sequencing.

CDKL5 deficiency causes epileptic seizures independent of cellular mosaicism.

OBJECTIVE: In a study using a mouse model of CDKL5 deficiency disorder (CDD), seizures are specific to female mice heterozygous for Cdkl5 mutations and not observed in hemizygous knockout males or homozygous knockout females. The aim of this study was to examine whether the clinical phenotype of patients with CDD can be impacted by the type of genetic variant. METHODS: Eleven CDD patients (six females and five males) were included in this study. The molecular diagnosis of hemizygous male patients was performed using digital PCR and their clinical phenotypes were compared with those of patients with mosaic or heterozygous CDKL5 variants. The severity of clinical phenotypes was graded by using CDKL5 Developmental Score and the adapted version of the CDKL5 Clinical Severity Assessment. The effect of cellular mosaicism on the severity of CDD was studied by comparing the clinical characteristics and comorbidities between individuals with hemizygous and mosaic or heterozygous CDKL5 variants. RESULTS: One of the five male patients was mosaic for the CDKL5 variant. All patients developed seizures irrespective of their genetic status of the pathogenic variant. However, cellular mosaicism of CDKL5 deficiency was associated with lesser severity of other comorbidities such as feeding, respiratory, and visual functional impairments. SIGNIFICANCE: This study provided evidence that cellular mosaicism of CDKL5 deficiency was not necessarily required for developing epilepsy. CDD patients not only exhibited clinical features of epilepsy but also exhibited the developmental consequences arising directly from the effect of the CDKL5 pathogenic variant.

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.

Glycolytic oscillations in HeLa cervical cancer cell spheroids.

Previous studies have unravelled glycolytic oscillations in cancer cells, such as HeLa cervical and DU145 prostate cancer cells, using a monolayer culture system. Here, we demonstrate glycolytic oscillations in HeLa cervical cancer cell spheroids. Experiments revealed that a small number of HeLa cells in spheroids exhibited heterogeneous oscillations with a higher frequency than those in monolayers. Model analyses and our previous experiments indicated that the higher frequencies of oscillations in spheroids were mostly due to the increase in glycolytic enzyme activity in the cells, and to the decrease in glucose concentration by diffusional transport of glucose from the surface to inside the spheroids, as well as the increase in cell density through spheroid formation. These results and our previous studies imply that more malignant cancer cells tend to exhibit glycolytic oscillations with higher frequencies than less malignant cells. Adjacent cells in spheroids oscillated within a 10% difference in frequency, but did not synchronize with each other. This suggests that weak cell-to-cell interactions might exist among HeLa cells connected with cadherins in the spheroid microenvironment; however, the interactions were not strong enough to induce synchronization of glycolytic oscillations.

Missense and truncating variants in CHD5 in a dominant neurodevelopmental disorder with intellectual disability, behavioral disturbances, and epilepsy.

Located in the critical 1p36 microdeletion region, the chromodomain helicase DNA-binding protein 5 (CHD5) gene encodes a subunit of the nucleosome remodeling and deacetylation (NuRD) complex required for neuronal development. Pathogenic variants in six of nine chromodomain (CHD) genes cause autosomal dominant neurodevelopmental disorders, while CHD5-related disorders are still unknown. Thanks to GeneMatcher and international collaborations, we assembled a cohort of 16 unrelated individuals harboring heterozygous CHD5 variants, all identified by exome sequencing. Twelve patients had de novo CHD5 variants, including ten missense and two splice site variants. Three familial cases had nonsense or missense variants segregating with speech delay, learning disabilities, and/or craniosynostosis. One patient carried a frameshift variant of unknown inheritance due to unavailability of the father. The most common clinical features included language deficits (81%), behavioral symptoms (69%), intellectual disability (64%), epilepsy (62%), and motor delay (56%). Epilepsy types were variable, with West syndrome observed in three patients, generalized tonic-clonic seizures in two, and other subtypes observed in one individual each. Our findings suggest that, in line with other CHD-related disorders, heterozygous CHD5 variants are associated with a variable neurodevelopmental syndrome that includes intellectual disability with speech delay, epilepsy, and behavioral problems as main features.

Proteomic Analysis of Hepatocellular Carcinoma Tissues With Encapsulation Shows Up-regulation of Leucine Aminopeptidase 3 and Phosphoenolpyruvate Carboxykinase 2.

BACKGROUND/AIM: Cancer is the most fatal disease worldwide whose most lethal characteristics are invasion and metastasis. Hepatocellular carcinoma (HCC) is one of the most fatal cancers worldwide. HCC often shows encapsulation, which is related to better prognosis. In this study, proteomic analysis of HCC tissues with and without encapsulation was performed, in order to elucidate the factors which play important roles in encapsulation. MATERIALS AND METHODS: Five HCC tissues surrounded by a capsule and five HCC tissues which broke the capsule were obtained from patients diagnosed with HCC who underwent surgical liver resection. Protein samples from these tissues were separated by two-dimensional gel electrophoresis (2-DE), and the protein spots whose expression was different between encapsulated and non-encapsulated HCC tissues were identified through gel imaging analysis software. The selected protein spots were analyzed and identified by liquid chromatography-tandem mass spectrometry (LC-MS/MS). RESULTS: Two-DE analysis showed 14 spots whose expression was different between encapsulated and non-encapsulated HCC tissues. Of these, 9 were up-regulated and 5 were down-regulated in HCC tissues without encapsulation. The validation by Western blot confirmed that leucine aminopeptidase 3 (LAP3) and phosphoenolpyruvate carboxykinase mitochondrial (PCK2) were up-regulated significantly in HCC tissues with a capsule, compared to HCC tissues that broke the capsule. CONCLUSION: These findings suggest that LAP3 and PCK2 could be factors responsible for the maintenance of encapsulation in HCC tissues.

De novo ATP1A3 variants cause polymicrogyria.

Polymicrogyria is a common malformation of cortical development whose etiology remains elusive. We conducted whole-exome sequencing for 124 patients with polymicrogyria and identified de novo ATP1A3 variants in eight patients. Mutated ATP1A3 causes functional brain diseases, including alternating hemiplegia of childhood (AHC), rapid-onset dystonia parkinsonism (RDP), and cerebellar ataxia, areflexia, pes cavus, optic nerve atrophy, and sensorineural deafness (CAPOS). However, our patients showed no clinical features of AHC, RDP, or CAPOS and had a completely different phenotype: a severe form of polymicrogyria with epilepsy and developmental delay. Detected variants had different locations in ATP1A3 and different functional properties compared with AHC-, RDP-, or CAPOS-associated variants. In the developing cerebral cortex of mice, radial neuronal migration was impaired in neurons overexpressing the ATP1A3 variant of the most severe patients, suggesting that this variant is involved in cortical malformation pathogenesis. We propose a previously unidentified category of polymicrogyria associated with ATP1A3 abnormalities.

Whole genome sequencing of 45 Japanese patients with intellectual disability.

Intellectual disability (ID) is characterized by significant limitations in both intellectual functioning and adaptive behaviors, originating before the age of 18 years. However, the genetic etiologies of ID are still incompletely elucidated due to the wide range of clinical and genetic heterogeneity. Whole genome sequencing (WGS) has been applied as a single-step clinical diagnostic tool for ID because it detects genetic variations with a wide range of resolution from single nucleotide variants (SNVs) to structural variants (SVs). To explore the causative genes for ID, we employed WGS in 45 patients from 44 unrelated Japanese families and performed a stepwise screening approach focusing on the coding variants in the genes. Here, we report 12 pathogenic and likely pathogenic variants: seven heterozygous variants of ADNP, SATB2, ANKRD11, PTEN, TCF4, SPAST, and KCNA2, three hemizygous variants of SMS, SLC6A8, and IQSEC2, and one homozygous variant in AGTPBP1. Of these, four were considered novel. Furthermore, a novel 76 kb deletion containing exons 1 and 2 in DYRK1A was identified. We confirmed the clinical and genetic heterogeneity and high frequency of de novo causative variants (8/12, 66.7%). This is the first report of WGS analysis in Japanese patients with ID. Our results would provide insight into the correlation between novel variants and expanded phenotypes of the disease.

Clinical manifestations and epilepsy treatment in Japanese patients with pathogenic CDKL5 variants.

OBJECTIVE: Patients with pathogenic cyclin-dependent kinase-like-5 gene (CDKL5) variants are designated CDKL5 deficiency disorder (CDD). This study aimed to delineate the clinical characteristics of Japanese patients with CDD and elucidate possible appropriate treatments. METHODS: We recruited patients with pathogenic or likely pathogenic CDKL5 variants from a cohort of approximately 1,100 Japanese patients with developmental and epileptic encephalopathies, who underwent genetic analysis. We retrospectively reviewed clinical, electroencephalogram, neuroimaging, and genetic information. RESULTS: We identified 29 patients (21 females, eight males). All patients showed severe developmental delay, especially in males. Involuntary movements were observed in 15 patients. No antiepileptic drugs (AEDs) achieved seizure freedom by monotherapy. AEDs achieving ≥ 50% reduction in seizure frequency were sodium valproate in two patients, vigabatrin in one, and lamotrigine in one. Seizure aggravation was observed during the use of lamotrigine, potassium bromide, and levetiracetam. Adrenocorticotrophic hormone (ACTH) was the most effective treatment. The ketogenic diet (KD), corpus callosotomy and vagus nerve stimulation did not improve seizure frequency in most patients, but KD was remarkably effective in one. The degree of brain atrophy on magnetic resonance imaging (MRI) reflected disease severity. Compared with females, males had lower levels of attained motor development and more severe cerebral atrophy on MRI. CONCLUSION: Our patients showed more severe global developmental delay than those in previous studies and had intractable epilepsy, likely because previous studies had lower numbers of males. Further studies are needed to investigate appropriate therapy for CDD, such as AED polytherapy or combination treatment involving ACTH, KD, and AEDs.

Resolving the data asynchronicity in high-speed atomic force microscopy measurement via the Kalman Smoother.

High-speed atomic force microscopy (HS-AFM) is a scanning probe microscopy that can capture structural dynamics of biomolecules in real time at single molecule level near physiological condition. Albeit much improvement, while scanning one frame of HS-AFM movies, biomolecules often change their conformations largely. Thus, the obtained frame images can be hampered by the time-difference, the asynchronicity, in the data acquisition. Here, to resolve this data asynchronicity in the HS-AFM movie, we developed Kalman filter and smoother methods, some of the sequential Bayesian filtering approaches. The Kalman filter/smoother methods use alternative steps of a short time-propagation by a linear dynamical system and a correction by the likelihood of AFM data acquired pixel by pixel. We first tested the method using a toy model of a diffusing cone, showing that the Kalman smoother method outperforms to reproduce the ground-truth movie. We then applied the Kalman smoother to a synthetic movie for conformational change dynamics of a motor protein, i.e., dynein, confirming the superiority of the Kalman smoother. Finally, we applied the Kalman smoother to two real HS-AFM movies, FlhAC and centralspindlin, reducing distortion and noise in the AFM movies. The method is general and can be applied to any HS-AFM movies.

Clinical and genetic characteristics of patients with Doose syndrome.

OBJECTIVE: To elucidate the genetic background and genotype-phenotype correlations for epilepsy with myoclonic-atonic seizures, also known as myoclonic-astatic epilepsy (MAE) or Doose syndrome. METHODS: We collected clinical information and blood samples from 29 patients with MAE. We performed whole-exome sequencing for all except one MAE case in whom custom capture sequencing identified a variant. RESULTS: We newly identified four variants: SLC6A1 and HNRNPU missense variants and microdeletions at 2q24.2 involving SCN1A and Xp22.31 involving STS. Febrile seizures preceded epileptic or afebrile seizures in four patients, of which two patients had gene variants. Myoclonic-atonic seizures occurred at onset in four patients, of which two had variants, and during the course of disease in three patients. Variants were more commonly identified in patients with a developmental delay or intellectual disability (DD/ID), but genetic status was not associated with the severity of DD/ID. Attention-deficit/hyperactivity disorder and autistic spectrum disorder were less frequently observed in patients with variants than in those with unknown etiology. SIGNIFICANCE: MAE patients had genetic heterogeneity, and HNRNPU and STS emerged as possible candidate causative genes. Febrile seizures prior to epileptic seizures and myoclonic-atonic seizure at onset indicate a genetic predisposition to MAE. Comorbid conditions were not related to genetic predisposition to MAE.

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.

Comprehensive analysis of coding variants highlights genetic complexity in developmental and epileptic encephalopathy.

Although there are many known Mendelian genes linked to epileptic or developmental and epileptic encephalopathy (EE/DEE), its genetic architecture is not fully explained. Here, we address this incompleteness by analyzing exomes of 743 EE/DEE cases and 2366 controls. We observe that damaging ultra-rare variants (dURVs) unique to an individual are significantly overrepresented in EE/DEE, both in known EE/DEE genes and the other non-EE/DEE genes. Importantly, enrichment of dURVs in non-EE/DEE genes is significant, even in the subset of cases with diagnostic dURVs (P = 0.000215), suggesting oligogenic contribution of non-EE/DEE gene dURVs. Gene-based analysis identifies exome-wide significant (P = 2.04 × 10-6) enrichment of damaging de novo mutations in NF1, a gene primarily linked to neurofibromatosis, in infantile spasm. Together with accumulating evidence for roles of oligogenic or modifier variants in severe neurodevelopmental disorders, our results highlight genetic complexity in EE/DEE, and indicate that EE/DEE is not an aggregate of simple Mendelian disorders.

Gene therapy improves motor and mental function of aromatic l-amino acid decarboxylase deficiency.

In patients with aromatic l-amino acid decarboxylase (AADC) deficiency, a decrease in catecholamines and serotonin levels in the brain leads to developmental delay and movement disorders. The beneficial effects of gene therapy in patients from 1 to 8 years of age with homogeneous severity of disease have been reported from Taiwan. We conducted an open-label phase 1/2 study of population including adolescent patients with different degrees of severity. Six patients were enrolled: four males (ages 4, 10, 15 and 19 years) and one female (age 12 years) with a severe phenotype who were not capable of voluntary movement or speech, and one female (age 5 years) with a moderate phenotype who could walk with support. The patients received a total of 2 × 1011 vector genomes of adeno-associated virus vector harbouring DDC via bilateral intraputaminal infusions. At up to 2 years after gene therapy, the motor function was remarkably improved in all patients. Three patients with the severe phenotype were able to stand with support, and one patient could walk with a walker, while the patient with the moderate phenotype could run and ride a bicycle. This moderate-phenotype patient also showed improvement in her mental function, being able to converse fluently and perform simple arithmetic. Dystonia disappeared and oculogyric crisis was markedly decreased in all patients. The patients exhibited transient choreic dyskinesia for a couple of months, but no adverse events caused by vector were observed. PET with 6-[18F]fluoro-l-m-tyrosine, a specific tracer for AADC, showed a persistently increased uptake in the broad areas of the putamen. In our study, older patients (>8 years of age) also showed improvement, although treatment was more effective in younger patients. The genetic background of our patients was heterogeneous, and some patients suspected of having remnant enzyme activity showed better improvement than the Taiwanese patients. In addition to the alleviation of motor symptoms, the cognitive and verbal functions were improved in a patient with the moderate phenotype. The restoration of dopamine synthesis in the putamen via gene transfer provides transformative medical benefit across all patient ages, genotypes, and disease severities included in this study, with the most pronounced improvements noted in moderate patients.10.1093/brain/awy331_video1awy331media15991361892001.

Optimization of fixative solution for retinal morphology: a comparison with Davidson's fixative and other fixation solutions.

PURPOSE: Numerous fixative solutions are available but many are not amenable to the histomorphological preservation of retinae. The investigators specifically focused on retinal histological studies, which rather than 4% formaldehyde (FA), often use Davidson's fixative. However the latter has its limitations. The purpose of this study was to produce a new fixative which maintains retinae closer to the in vivo conditions. STUDY DESIGN: Experimental design. METHODS: Four fixative formulations (4% paraformaldehyde, Davidson's fixative, modified Davidson's fixative and an in-house fixative - TB-Fix) were tested on retinae and the outcomes on histomorphology and immunohistochemical staining for selected antigenic markers was compared. RESULTS: TB-Fix markedly improved morphological detail following hematoxylin and eosin staining, most importantly eliminating the spongiform appearance in the plexiform layer and the swelling of somata (including Müller cells), when compared to FA, Davidson's fixative and its modified version. Retinal samples fixed with TB-Fix or FA showed comparable results in immunohistological staining for neurons and glia in the retina. Importantly, while the whole eye fixed with FA collapsed in shape and induced artificial retinal detachment, the eye fixed with TB-Fix avoided deformation and detachment. Furthermore, we found that TB-Fix also prevented detachment from the culture plate when used to fix HEK293 cells, which are known to detach from the plate easily. CONCLUSION: It was demonstrated that TB-Fix provides an overall improvement in the preservation of retinal morphology and chemical composition.

Integrative Analyses of De Novo Mutations Provide Deeper Biological Insights into Autism Spectrum Disorder.

Recent studies have established important roles of de novo mutations (DNMs) in autism spectrum disorders (ASDs). Here, we analyze DNMs in 262 ASD probands of Japanese origin and confirm the "de novo paradigm" of ASDs across ethnicities. Based on this consistency, we combine the lists of damaging DNMs in our and published ASD cohorts (total number of trios, 4,244) and perform integrative bioinformatics analyses. Besides replicating the findings of previous studies, our analyses highlight ATP-binding genes and fetal cerebellar/striatal circuits. Analysis of individual genes identified 61 genes enriched for damaging DNMs, including ten genes for which our dataset now contributes to statistical significance. Screening of compounds altering the expression of genes hit by damaging DNMs reveals a global downregulating effect of valproic acid, a known risk factor for ASDs, whereas cardiac glycosides upregulate these genes. Collectively, our integrative approach provides deeper biological and potential medical insights into ASDs.