Exome sequencing of 457 autism families recruited online provides evidence for autism risk genes

Characterizing Rare DNA Copy-Number Variants in Pediatric Obsessive-Compulsive Disorder

OBJECTIVE

Pediatric obsessive-compulsive disorder (OCD) is a common neuropsychiatric disorder for which genetic factors play an important role. Recent studies have demonstrated an enrichment of rare de novo DNA single nucleotide variants in OCD cases compared to controls, and larger studies have examined copy-number variants (CNVs) using microarray data. Our study examines rare de novo CNVs using whole-exome sequencing (WES) data to provide additional insight into genetic factors and biological processes underlying OCD.

METHOD

We detected CNVs using whole-exome DNA sequencing (WES) data from 183 OCD trio families (unaffected parents and children with OCD) and 771 control families to test the hypothesis that rare de novo CNVs are enriched in OCD cases compared to controls. Our primary analysis used the eXome-Hidden Markov Model (XHMM) to identify CNVs in silico. We performed burden analyses comparing individuals with OCD vs. controls and downstream biological systems analyses of CNVs in probands with OCD. We then used a second algorithm (GATK-gCNV) to confirm our primary analysis.

RESULTS

Our findings demonstrate a higher rate of rare de novo CNVs detected by WES in individuals with OCD (0.07 CNVs per proband) compared to controls (0.005) (corrected rate ratio = 11.7 95% CI, 3.6-50.0, p = 4.00x10-6). We confirmed this enrichment using GATK-gCNV. The majority of these rare de novo CNVs in OCD cases are predicted to be pathogenic or likely pathogenic, and an examination of genes disrupted by rare de novo CNVs in OCD cases finds enrichment of several gene-ontology sets.

CONCLUSION

This study shows for the first time an enrichment of rare de novo CNVs detected by WES in OCD, complementing previous larger CNV studies and providing additional insight into genetic factors underlying OCD risk.

Mosaic H3K9me3 at BREACHes predicts synaptic gene expression associated with fragile X syndrome cognitive severity

Diseases vary in clinical presentation across individuals despite the same molecular diagnosis. In fragile X syndrome (FXS), mutation-length expansion of a CGG short tandem repeat (STR) in FMR1 causes reduced gene expression and FMRP loss. Nevertheless, FMR1 and FMRP are limited predictors of adaptive functioning and cognition in FXS patients, suggesting that molecular correlates of clinical measures would add diagnostic value. We recently uncovered Megabase-scale domains of heterochromatin (BREACHes) in FXS patient-derived iPSCs. Here, we identify BREACHes in FXS brain tissue (N=4) and absent from sex/age-matched neurotypical controls (N=4). BREACHes span >250 genes and exhibit patient-specific H3K9me3 variation. Using N=4 FXS iPSC lines and N=7 single-cell isogenic FXS iPSC subclones, we observe a strong correlation between inter-sample H3K9me3 variation and heterogeneous BREACH gene repression. We demonstrate improved prediction of cognitive metrics in FXS patients with an additive model of blood FMRP and mRNA levels of H3K9me3-mosaic, but not H3K9me3-invariant, BREACH genes. Our results highlight the utility of H3K9me3 variation at BREACHes for identifying genes associated with FXS clinical metrics.

Uncovering convergence and divergence between autism and schizophrenia using genomic tools and patients' neurons

Autism spectrum disorders (ASDs) are highly heritable and result in abnormal repetitive behaviors and impairment in communication and cognitive skills. Previous studies have focused on the genetic correlation between ASDs and other neuropsychiatric disorders, but an in-depth understanding of the correlation to other disorders is required. We conducted an extensive meta-analysis of common variants identified in ASDs by genome-wide association studies (GWAS) and compared it to the consensus genes and single nucleotide polymorphisms (SNPs) of Schizophrenia (SCZ). We found approximately 75% of the GWAS genes that are associated with ASD are also associated with SCZ. We further investigated the cellular phenotypes of neurons derived from induced pluripotent stem cell (iPSC) models in ASD and SCZ. Our findings revealed that ASD and SCZ neurons initially follow divergent developmental trajectories compared to control neurons. However, despite these early diametrical differences, both ASD and SCZ neurons ultimately display similar deficits in synaptic activity as they mature. This significant genetic overlap between ASD and SCZ, coupled with the convergence towards similar synaptic deficits, highlights the intricate interplay of genetic and developmental factors in shaping the shared underlying mechanisms of these complex neurodevelopmental and neuropsychiatric disorders.

Review: Dopamine, Serotonin, and the Translational Neuroscience of Aggression in Autism Spectrum Disorder

Objective

Autism spectrum disorder (ASD) is a neurodevelopmental disorder with a 1% to 2% prevalence in children. In addition to social communication deficits and restricted or repetitive behavior, ASD is often characterized by a heightened propensity for aggression. In fact, aggressive behavior is the primary reason for hospitalization in children with ASD, and current treatment options, despite some efficacy, are often associated with prominent side effects. Despite such high clinical toll, the neurobiology of aggression in ASD remains poorly understood.

Method

The neural circuits linked to both ASD and aggression were reviewed, with the goal of identifying overlapping components to help guide future treatment development. In discussing the clinical phenotype of aggression in ASD, some of the triggers and risk factors were noted to differ from those that cause aggression in neurotypical children. Preclinical and clinical studies on the neurobiology of aggression and ASD were synthesized to combine evidence from genetics, neuroimaging, pharmacology, and circuit manipulations. Dopamine and serotonin, 2 neuromodulators that contribute to development and behavioral control, were specifically studied.

Results

The literature indicates that the intricate interplay of the dopamine and serotonin systems has a pivotal role in shaping behavior, including the expression of aggression.

Conclusion

Understanding the balance between dopamine as an accelerator and serotonin as a brake may provide insights into the mechanisms of aggression in children with ASD. Although much work remains to be done, new perspectives promise to bridge the gap between human and animal studies and pinpoint the neurobiology of aggression in ASD.

Diversity & Inclusion Statement

One or more of the authors of this paper self-identifies as a member of one or more historically underrepresented sexual and/or gender groups in science. We actively worked to promote sex and gender balance in our author group.

Genetic variants and phenotypic data curated for the CAGI6 intellectual disability panel challenge

Neurodevelopmental disorders (NDDs) are common conditions including clinically diverse and genetically heterogeneous diseases, such as intellectual disability, autism spectrum disorders, and epilepsy. The intricate genetic underpinnings of NDDs pose a formidable challenge, given their multifaceted genetic architecture and heterogeneous clinical presentations. This work delves into the intricate interplay between genetic variants and phenotypic manifestations in neurodevelopmental disorders, presenting a dataset curated for the Critical Assessment of Genome Interpretation (CAGI6) ID Panel Challenge. The CAGI6 competition serves as a platform for evaluating the efficacy of computational methods in predicting phenotypic outcomes from genetic data. In this study, a targeted gene panel sequencing has been used to investigate the genetic causes of NDDs in a cohort of 415 paediatric patients. We identified 60 pathogenic and 49 likely pathogenic variants in 102 individuals that accounted for 25% of NDD cases in the cohort. The most mutated genes were ANKRD11, MECP2, ARID1B, ASH1L, CHD8, KDM5C, MED12 and PTCHD1 The majority of pathogenic variants were de novo, with some inherited from mildly affected parents. Loss-of-function variants were the most common type of pathogenic variant. In silico analysis tools were used to assess the potential impact of variants on splicing and structural/functional effects of missense variants. The study highlights the challenges in variant interpretation especially in cases with atypical phenotypic manifestations. Overall, this study provides valuable insights into the genetic causes of NDDs and emphasises the importance of understanding the underlying genetic factors for accurate diagnosis, and intervention development in neurodevelopmental conditions.

Saliva Sample-Based Non-Invasive Carrier Screening for Spinal Muscular Atrophy, Hereditary Hearing Loss, and Thalassemia in 13,926 Women of Reproductive Age From South Zhejiang

BACKGROUND

Although spinal muscular atrophy (SMA), hereditary hearing loss (HL), and thalassemia are common monogenic genetic diseases, the carrier frequencies and variant spectrums of these diseases show regional differences, even within China. Their carrier frequencies and variant spectrums in Southern Zhejiang, China are unclear.

METHODS

Saliva was collected for carrier screening and amniotic fluid, villi, and peripheral blood were collected for prenatal diagnosis. Real-time quantitative polymerase chain reaction (PCR) and multiplex ligation-dependent probe amplification (MLPA) were used to detect the copy number of SMN1 exon 7. PCR coupled with flow-through hybridization, MLPA, and Sanger sequencing were used to detect common genes for HL and thalassemia.

RESULTS

Common variants were detected in 15.14% (2109/13926) of the 13,926 women of reproductive age from South Zhejiang who participated in this study. The carrier frequencies of SMA, HL, and thalassemia were 2.11% (294/13926), 4.87% (678/13926), and 8.82% (1228/13926), respectively. In total, 56.47% (1117/1978) of husbands were successfully recalled. The total number of at-risk couples was 111 (111/13926, 0.80%). Further, 47 families underwent prenatal diagnosis. A total of 13 (13/13926; 0.93‰) affected pregnancies were identified.

CONCLUSION

Our findings confirm that SMA, HL, and thalassemia are highly prevalent in Southern Zhejiang, with some regional specificity, as compared with recent large population-based surveys in China. Further, a rapid saliva sample-based non-invasive screening method was established, and its feasibility was demonstrated.

Sex-specific DNA methylation signatures of autism spectrum disorder in newborn blood

Background

Autism spectrum disorder (ASD) comprises a group of neurodevelopmental conditions currently diagnosed by behavioral assessment in childhood, although neuropathology begins during gestation. A poorly understood male bias for ASD diagnosis is thought to be due to both biological sex differences and cultural biases against female diagnosis of ASD. Identification of molecular biomarkers of ASD likelihood in newborns would provide more objective screening and early intervention. Epigenetic dysregulation has been reported in multiple tissues from newborns who are later diagnosed with ASD, but this is the first study to investigate sex-specific DNA methylation signatures for ASD in newborn blood, an accessible and widely banked tissue.

Methods

DNA methylation was assayed from ASD and typically developing (TD) newborn blood (discovery set n = 196, replication set n = 90) using whole genome bisulfite sequencing (WGBS). Sex-stratified differentially methylated regions (DMRs) were assessed for replication, comparisons by sex, overlaps with DMRs from other tissues, and enrichment for biological processes and SFARI ASD-risk genes.

Results

We found that newborn blood ASD DMRs from both sexes significantly replicated in an independent cohort and were enriched for hypomethylation in ASD compared to TD samples, as well as location in promoters, CpG islands and CpG shores. Comparing females and males, we found that most DMRs with sex differences amongst TD individuals were also found in ASD individuals, plus many additional DMRs with sex differences that were only found in those with ASD. Newborn blood DMRs from females were enriched for the X chromosome and both sexes showed significant overlap with DMRs from umbilical cord blood and placenta but not post-mortem cortex. DMRs from all tissues were enriched for neurodevelopmental processes (females) and SFARI ASD-risk genes (females and males).

Limitations

This study is primarily limited by sample sizes, particularly amongst females.

Conclusions

Overall, we found a highly replicated sex-specific DNA methylation signature of ASD in newborn blood that showed support for the female protective effect and convergence with epigenetic and genetic signatures of ASD in newborns. These results demonstrate the utility of newborn blood in ASD screening and emphasizes the importance of sex-stratification in future studies.

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.

Detailed phenotyping of Tbr1-2A-CreER knock-in mice demonstrates significant impacts on TBR1 protein levels and axon development

Cre recombinase knock-in mouse lines have served as invaluable genetic tools for understanding key developmental processes altered in autism. However, insertion of exogenous DNA into the genome can have unintended effects on local gene regulation or protein function that must be carefully considered. Here, we analyze a recently generated Tbr1-2A-CreER knock-in mouse line, where a 2A-CreER cassette was inserted in-frame before the stop codon of the transcription factor gene Tbr1. Heterozygous TBR1 mutations in humans and mice are known to cause autism or autism-like behavioral phenotypes accompanied by structural brain malformations, most frequently a reduction of the anterior commissure (AC). Thus, it is critical for modified versions of Tbr1 to exhibit true wild-type-like activity. We evaluated the Tbr1-2A-CreER allele for its potential impact on Tbr1 function and complementation to Tbr1 loss-of-function alleles. In mice with one copy of the Tbr1-2A-CreER allele, we identified reduction of TBR1 protein in early postnatal cortex along with thinning of the AC, suggesting hypersensitivity of this structure to TBR1 dosage. Comparing Tbr1-2A-CreER and Tbr1-null mice to Tbr1-null complementation crosses showed reductions of TBR1 dosage ranging from 20% to 100%. Using six combinatorial genotypes, we found that moderate to severe TBR1 reductions (≥44%) were associated with cortical layer 5 expansion, while only the complete absence of TBR1 was associated with reeler-like "inverted" cortical layering. In total, these results strongly support the conclusion that Tbr1-2A-CreER is a hypomorphic allele. We advise caution when interpreting experiments using this allele, considering the sensitivity of various corticogenic processes to TBR1 dosage and the association of heterozygous TBR1 mutations with complex neurodevelopmental disorders.

Large-scale metagenomic analysis of oral microbiomes reveals markers for autism spectrum disorders

The link between the oral microbiome and neurodevelopmental disorders remains a compelling hypothesis, still requiring confirmation in large-scale datasets. Leveraging over 7000 whole-genome sequenced salivary samples from 2025 US families with children diagnosed with autism spectrum disorders (ASD), our cross-sectional study shows that the oral microbiome composition can discriminate ASD subjects from neurotypical siblings (NTs, AUC = 0.66), with 108 differentiating species (q < 0.005). The relative abundance of these species is highly correlated with cognitive impairment as measured by Full-Scale Intelligence Quotient (IQ). ASD children with IQ < 70 also exhibit lower microbiome strain sharing with parents (p < 10-6) with respect to NTs. A two-pronged functional enrichment analysis suggests the contribution of enzymes from the serotonin, GABA, and dopamine degradation pathways to the distinct microbial community compositions observed between ASD and NT samples. Although measures of restrictive eating diet and proxies of oral hygiene show relatively minor effects on the microbiome composition, the observed associations with ASD and IQ may still represent unaccounted-for underlying differences in lifestyle among groups. While causal relationships could not be established, our study provides substantial support to the investigation of oral microbiome biomarkers in ASD.

A deep learning model for prediction of autism status using whole-exome sequencing data

Autism is a developmental disability. Research demonstrated that children with autism benefit from early diagnosis and early intervention. Genetic factors are considered major contributors to the development of autism. Machine learning (ML), including deep learning (DL), has been evaluated in phenotype prediction, but this method has been limited in its application to autism. We developed a DL model, the Separate Translated Autism Research Neural Network (STAR-NN) model to predict autism status. The model was trained and tested using whole exome sequencing data from 43,203 individuals (16,809 individuals with autism and 26,394 non-autistic controls). Polygenic scores from common variants and the aggregated count of rare variants on genes were used as input. In STAR-NN, protein truncating variants, possibly damaging missense variants and mild effect missense variants on the same gene were separated at the input level and merged to one gene node. In this way, rare variants with different level of pathogenic effects were treated separately. We further validated the performance of STAR-NN using an independent dataset, including 13,827 individuals with autism and 14,052 non-autistic controls. STAR-NN achieved a modest ROC-AUC of 0.7319 on the testing dataset and 0.7302 on the independent dataset. STAR-NN outperformed other traditional ML models. Gene Ontology analysis on the selected gene features showed an enrichment for potentially informative pathways including calcium ion transport.

Deciphering autism heterogeneity: a molecular stratification approach in four mouse models

Autism spectrum disorder (ASD) is a complex neurodevelopmental condition characterized by impairments in social interaction and communication, as well as restrained or stereotyped behaviors. The inherent heterogeneity within the autism spectrum poses challenges for developing effective pharmacological treatments targeting core features. Successful clinical trials require the identification of robust markers to enable patient stratification. In this study, we identified molecular markers within the oxytocin and immediate early gene families across five interconnected brain structures of the social circuit. We used wild-type and four heterogeneous mouse models, each exhibiting unique autism-like behaviors modeling the autism spectrum. While dysregulations in the oxytocin family were model-specific, immediate early genes displayed widespread alterations, reflecting global changes across the four models. Through integrative analysis, we identified Egr1, Foxp1, Homer1a, Oxt, and Oxtr as five robust and discriminant molecular markers that allowed the successful stratification of the four models. Importantly, our stratification demonstrated predictive values when challenged with a fifth mouse model or identifying subgroups of mice potentially responsive to oxytocin treatment. Beyond providing insights into oxytocin and immediate early gene mRNA dynamics, this proof-of-concept study represents a significant step toward the potential stratification of individuals with ASD. This work has implications for the success of clinical trials and the development of personalized medicine in autism.

Return of genetic research results in 21,532 individuals with autism

PURPOSE

The aim of this study is to identify likely pathogenic (LP) and pathogenic (P) genetic results for autism that can be returned to participants in SPARK (SPARKforAutism.org): a large recontactable cohort of people with autism in the United States. We also describe the process to return these clinically confirmed genetic findings.

METHODS

We present results from microarray genotyping and exome sequencing of 21,532 individuals with autism and 17,785 of their parents. We returned LP and P (American College of Medical Genetics criteria) copy-number variants, chromosomal aneuploidies, and variants in genes with strong evidence of association with autism and intellectual disability.

RESULTS

We identified 1903 returnable LP/P variants in 1861 individuals with autism (8.6%). 89.5% of these variants were not known to participants. The diagnostic genetic result was returned to 589 participants (53% of those contacted). Features associated with a higher probability of having a returnable result include cognitive and medically complex features, being female, being White (versus non-White) and being diagnosed more than 20 years ago. We also find results among autistics across the spectrum, as well as in transmitting parents with neuropsychiatric features but no autism diagnosis.

CONCLUSION

SPARK offers an opportunity to assess returnable results among autistic people who have not been ascertained clinically. SPARK also provides practical experience returning genetic results for a behavioral condition at a large scale.

Barriers, motivators and strategies to increase participation in genetic research among Asian and Black families of autistic individuals

Genetic research can help advance our knowledge of autism and positively impact the progress of care for individuals with autism. Asian American and Pacific Islander (AAPI) and Black participants remain significantly underrepresented in genetic research in autism in the United States, including nationwide, multisite, genetic consortiums like Simons Foundation Powering Autism Research for Knowledge (SPARK). Few studies have explored the unique motivators and barriers that influence participation in genetics research across underrepresented groups with autism and strategies to increase participation. Therefore, the aim of this study was to understand the perspectives of AAPI and Black parents of individuals with autism about participating in genetic research, specifically motivators (e.g., desire to know more about the relationship between autism and genetics) and/or barriers (e.g., mistrust of research staff) that may impact their decision to participate in genetic research. Using a mixed-methods approach, we collected surveys (n = 134) across the United States and conducted three focus groups with parents of individuals with autism (n = 16) who identified as AAPI and Black from two large metropolitan cities. No significant differences were observed in the survey data but findings from the focus groups elucidate shared motivators for participation (e.g., to help advance the autism field for future generations) and nuanced differences in barriers that influence Black and AAPI parents' decision to participate (e.g., different beliefs about the source of autism). Practical suggestions to improve outreach and study engagement in genetic research in autism were identified and discussed.

Unveiling the role of IGF1R in autism spectrum disorder: a multi-omics approach to decipher common pathogenic mechanisms in the IGF signaling pathway

Autism spectrum disorder (ASD) is a complex neurodevelopmental condition marked by impairments in social interaction, communication, and repetitive behaviors. Emerging evidence suggests that the insulin-like growth factor (IGF) signaling pathway plays a critical role in ASD pathogenesis; however, the precise pathogenic mechanisms remain elusive. This study utilizes multi-omics approaches to investigate the pathogenic mechanisms of ASD susceptibility genes within the IGF pathway. Whole-exome sequencing (WES) revealed a significant enrichment of rare variants in key IGF signaling components, particularly the IGF receptor 1 (IGF1R), in a cohort of Chinese Han individuals diagnosed with ASD, as well as in ASD patients from the SFARI SPARK WES database. Subsequent single-cell RNA sequencing (scRNA-seq) of cortical tissues from children with ASD demonstrated elevated expression of IGF receptors in parvalbumin (PV) interneurons, suggesting a substantial impact on their development. Notably, IGF1R appears to mediate the effects of IGF2R on these neurons. Additionally, transcriptomic analysis of brain organoids derived from ASD patients indicated a significant association between IGF1R and ASD. Protein-protein interaction (PPI) and gene regulatory network (GRN) analyses further identified ASD susceptibility genes that interact with and regulate IGF1R expression. In conclusion, IGF1R emerges as a central node within the IGF signaling pathway, representing a potential common pathogenic mechanism and therapeutic target for ASD. These findings highlight the need for further investigation into the modulation of this pathway as a strategy for ASD intervention.

Beyond the Spectrum: Subtype-Specific Molecular Insights into Autism Spectrum Disorder Via Multimodal Data Integration

Some toddlers with autism spectrum disorder (ASD) have mild social symptoms and developmental improvement in skills, but for others, symptoms and abilities are moderately or even severely affected. Those with profound autism have the most severe social, language, and cognitive symptoms and are at the greatest risk of having a poor developmental outcome. The little that is known about the underlying biology of this important profound autism subtype, points clearly to embryonic dysregulation of proliferation, differentiation and neurogenesis. Because it is essential to gain foundational knowledge of the molecular biology associated with profound, moderate, and mild autism clinical subtypes, we used well-validated, data-driven patient subtyping methods to integrate clinical and molecular data at 1 to 3 years of age in a cohort of 363 ASD and controls representative of the general pediatric population in San Diego County. Clinical data were diagnostic, language, cognitive and adaptive ability scores. Molecular measures were 50 MSigDB Hallmark gene pathway activity scores derived from RNAseq gene expression. Subtyping identified four ASD, typical and mixed diagnostic clusters. 93% of subjects in one cluster were profound autism and 93% in a different cluster were control toddlers; a third cluster was 76% moderate ability ASD; and the last cluster was a mix of mild ASD and control toddlers. Among the four clusters, the profound autism subtype had the most severe social symptoms, language, cognitive, adaptive, social attention eye tracking, social fMRI activation, and age-related decline in abilities, while mild autism toddlers mixed within typical and delayed clusters had mild social symptoms, and neurotypical language, cognitive and adaptive scores that improved with age compared with profound and moderate autism toddlers in other clusters. In profound autism, 7 subtype-specific dysregulated gene pathways were found; they control embryonic proliferation, differentiation, neurogenesis, and DNA repair. To find subtype-common dysregulated pathways, we compared all ASD vs TD and found 17 ASD subtype-common dysregulated pathways. These common pathways showed a severity gradient with the greatest dysregulation in profound and least in mild. Collectively, results raise the new hypothesis that the continuum of ASD heterogeneity is moderated by subtype-common pathways and the distinctive nature of profound autism is driven by the differentially added profound subtype-specific embryonic pathways.

Loss-of-Function Variants in CUL3 Cause a Syndromic Neurodevelopmental Disorder

OBJECTIVE

De novo variants in cullin-3 ubiquitin ligase (CUL3) have been strongly associated with neurodevelopmental disorders (NDDs), but no large case series have been reported so far. Here, we aimed to collect sporadic cases carrying rare variants in CUL3, describe the genotype-phenotype correlation, and investigate the underlying pathogenic mechanism.

METHODS

Genetic data and detailed clinical records were collected via multicenter collaboration. Dysmorphic facial features were analyzed using GestaltMatcher. Variant effects on CUL3 protein stability were assessed using patient-derived T-cells.

RESULTS

We assembled a cohort of 37 individuals with heterozygous CUL3 variants presenting a syndromic NDD characterized by intellectual disability with or without autistic features. Of these, 35 have loss-of-function (LoF) and 2 have missense variants. CUL3 LoF variants in patients may affect protein stability leading to perturbations in protein homeostasis, as evidenced by decreased ubiquitin-protein conjugates in vitro. Notably, we show that 4E-BP1 (EIF4EBP1), a prominent substrate of CUL3, fails to be targeted for proteasomal degradation in patient-derived cells.

INTERPRETATION

Our study further refines the clinical and mutational spectrum of CUL3-associated NDDs, expands the spectrum of cullin RING E3 ligase-associated neuropsychiatric disorders, and suggests haploinsufficiency via LoF variants is the predominant pathogenic mechanism. ANN NEUROL 2024.

Impact of KDM6B mosaic brain knockout on synaptic function and behavior

Autism spectrum disorders (ASD) are complex neurodevelopmental conditions characterized by impairments in social communication, repetitive behaviors, and restricted interests. Epigenetic modifications serve as critical regulators of gene expression playing a crucial role in controlling brain function and behavior. Lysine (K)-specific demethylase 6B (KDM6B), a stress-inducible H3K27me3 demethylase, has emerged as one of the highest ASD risk genes, but the precise effects of KDM6B mutations on neuronal activity and behavioral function remain elusive. Here we show the impact of KDM6B mosaic brain knockout on the manifestation of different autistic-like phenotypes including repetitive behaviors, social interaction, and significant cognitive deficits. Moreover, KDM6B mosaic knockout display abnormalities in hippocampal excitatory synaptic transmission decreasing NMDA receptor mediated synaptic transmission and plasticity. Understanding the intricate interplay between epigenetic modifications and neuronal function may provide novel insights into the pathophysiology of ASD and potentially inform the development of targeted therapeutic interventions.

Association of genetic variants with autism spectrum disorder in Japanese children revealed by targeted sequencing

Introduction

Autism spectrum disorders (ASD) represent a heterogeneous group of neurodevelopmental disorders with strong genetic predispositions. Although an increasing number of genetic variants have been implicated in the pathogenesis of ASD, little is known about the relationship between ASD-associated genetic variants and individual ASD traits. Therefore, we aimed to investigate these relationships.

Methods

Here, we report a case-control association study of 32 Japanese children with ASD (mainly with high-functioning autism [HFA]) and 36 with typical development (TD). We explored previously established ASD-associated genes using a next-generation sequencing panel and determined the association between Social Responsiveness Scale (SRS) T-scores and intelligence quotient (IQ) scores.

Results

In the genotype-phenotype analyses, 40 variants of five genes (SCN1A, SHANK3, DYRK1A, CADPS, and SCN2A) were associated with ASD/TD phenotypes. In particular, 10 SCN1A variants passed permutation filtering (false discovery rate <0.05). In the quantitative association analyses, 49 variants of 12 genes (CHD8, SCN1A, SLC6A1, KMT5B, CNTNAP2, KCNQ3, SCN2A, ARID1B, SHANK3, DYRK1A, FOXP1, and GRIN2B) and 50 variants of 10 genes (DYRK1A, SCN2A, SLC6A1, ARID1B, CNTNAP2, SHANK3, FOXP1, PTEN, SCN1A, and CHD8) were associated with SRS T- and IQ-scores, respectively.

Conclusion

Our data suggest that these identified variants are essential for the genetic architecture of HFA.

hnRNP R promotes O-GlcNAcylation of eIF4G and facilitates axonal protein synthesis

Motoneurons critically depend on precise spatial and temporal control of translation for axon growth and the establishment and maintenance of neuromuscular connections. While defects in local translation have been implicated in the pathogenesis of motoneuron disorders, little is known about the mechanisms regulating axonal protein synthesis. Here, we report that motoneurons derived from Hnrnpr knockout mice show reduced axon growth accompanied by lowered synthesis of cytoskeletal and synaptic components in axons. Mutant mice display denervated neuromuscular junctions and impaired motor behavior. In axons, hnRNP R is a component of translation initiation complexes and, through interaction with O-linked β-N-acetylglucosamine (O-GlcNAc) transferase (Ogt), modulates O-GlcNAcylation of eIF4G. Restoring axonal O-GlcNAc levels rescued local protein synthesis and axon growth defects of hnRNP R knockout motoneurons. Together, these findings demonstrate a function of hnRNP R in controlling the local production of key factors required for axon growth and formation of neuromuscular innervations.

Rare de novo damaging DNA variants are enriched in attention-deficit/hyperactivity disorder and implicate risk genes

Research demonstrates the important role of genetic factors in attention-deficit/hyperactivity disorder (ADHD). DNA sequencing of families provides a powerful approach for identifying de novo (spontaneous) variants, leading to the discovery of hundreds of clinically informative risk genes for other childhood neurodevelopmental disorders. This approach has yet to be extensively leveraged in ADHD. We conduct whole-exome DNA sequencing in 152 families, comprising a child with ADHD and both biological parents, and demonstrate a significant enrichment of rare and ultra-rare de novo gene-damaging mutations in ADHD cases compared to unaffected controls. Combining these results with a large independent case-control DNA sequencing cohort (3206 ADHD cases and 5002 controls), we identify lysine demethylase 5B (KDM5B) as a high-confidence risk gene for ADHD and estimate that 1057 genes contribute to ADHD risk. Using our list of genes harboring ultra-rare de novo damaging variants, we show that these genes overlap with previously reported risk genes for other neuropsychiatric conditions and are enriched in several canonical biological pathways, suggesting early neurodevelopmental underpinnings of ADHD. This work provides insight into the biology of ADHD and demonstrates the discovery potential of DNA sequencing in larger parent-child trio cohorts.

Single nucleotide variations encoding missense mutations in G protein-coupled receptors may contribute to autism

Autism is a neurodevelopmental condition with a range of symptoms that vary in intensity and severity from person to person. Genetic sequencing has identified thousands of genes containing mutations in autistic individuals, which may contribute to the development of autistic symptoms. Several of these genes encode G protein-coupled receptors (GPCRs), which are cell surface expressed proteins that transduce extracellular messages to the intracellular space. Mutations in GPCRs can impact their function, resulting in aberrant signalling within cells and across neurotransmitter systems in the brain. This review summarises the current knowledge on autism-associated single nucleotide variations encoding missense mutations in GPCRs and the impact of these genetic mutations on GPCR function. For some autism-associated mutations, changes in GPCR expression levels, ligand affinity, potency and efficacy have been observed. However, for many the functional consequences remain unknown. Thus, further work to characterise the functional impacts of the genetically identified mutations is required. LINKED ARTICLES: This article is part of a themed issue Therapeutic Targeting of G Protein-Coupled Receptors: hot topics from the Australasian Society of Clinical and Experimental Pharmacologists and Toxicologists 2021 Virtual Annual Scientific Meeting. To view the other articles in this section visit http://onlinelibrary.wiley.com/doi/10.1111/bph.v181.14/issuetoc.

Rare predicted deleterious FEZF2 variants are associated with a neurodevelopmental phenotype

FEZF2 encodes a transcription factor critical to neurodevelopment that regulates other neurodevelopment genes. Rare variants in FEZF2 have previously been suggested to play a role in autism, and cases of 3p14 microdeletions that include FEZF2 share a neurodevelopmental phenotype including mild dysmorphic features and intellectual disability. We identified seven heterozygous predicted deleterious variants in FEZF2 (three frameshifts, one recurrent missense in two independent cases, one nonsense, and one complete gene deletion) in unrelated individuals with neurodevelopmental disorders including developmental delay/intellectual disability, autism, and/or attention-deficit/hyperactivity. Variants were confirmed to be de novo in five of seven cases and paternally inherited from an affected father in one. Predicted deleterious variants in FEZF2 may affect the expression of genes that are involved in fate choice pathways in developing neurons, and thus contribute to the neurodevelopmental phenotype. Future studies are needed to clarify the mechanism by which FEZF2 leads to this neurodevelopmental disorder.

The Importance of Large-Scale Genomic Studies to Unravel Genetic Risk Factors for Autism

Autism spectrum disorder (ASD) is a common and highly heritable neurodevelopmental disorder. During the last 15 years, advances in genomic technologies and the availability of increasingly large patient cohorts have greatly expanded our knowledge of the genetic architecture of ASD and its neurobiological mechanisms. Over two hundred risk regions and genes carrying rare de novo and transmitted high-impact variants have been identified. Additionally, common variants with small individual effect size are also important, and a number of loci are now being uncovered. At the same time, these new insights have highlighted ongoing challenges. In this perspective article, we summarize developments in ASD genetic research and address the enormous impact of large-scale genomic initiatives on ASD gene discovery.

NR4A2 as a Novel Target Gene for Developmental and Epileptic Encephalopathy: A Systematic Review of Related Disorders and Therapeutic Strategies

The NR4A2 gene encodes an orphan transcription factor of the steroid-thyroid hormone-retinoid receptor superfamily. This review focuses on the clinical findings associated with the pathogenic variants so far reported, including three unreported cases. Also, its role in neurodegenerative diseases, such as Parkinson's or Alzheimer's disease, is examined, as well as a brief exploration on recent proposals to develop novel therapies for these neurological diseases based on small molecules that could modulate NR4A2 transcriptional activity. The main characteristic shared by all patients is mild to severe developmental delay/intellectual disability. Moderate to severe disorder of the expressive and receptive language is present in at least 42%, while neuro-psychiatric issues were reported in 53% of patients. Movement disorders, including dystonia, chorea or ataxia, are described in 37% patients, although probably underestimated because of its frequent onset in late adolescence-young adulthood. Finally, epilepsy was surprisingly present in 42% of patients, being drug-resistant in three of them. The age at onset varied widely, from five months to twenty-six years, as did the classification of epilepsy, which ranged from focal epilepsy to infantile spasms or Lennox-Gastaut syndrome. Accordingly, we propose that NR4A2 should be considered as a first-tier target gene for the genetic diagnosis of developmental and epileptic encephalopathy.

Detailed phenotyping of Tbr1-2A-CreER knock-in mice demonstrates significant impacts on TBR1 protein levels and axon development

Spatiotemporal control of Cre-mediated recombination has been an invaluable tool for understanding key developmental processes. For example, knock-in of Cre into cell type marker gene loci drives Cre expression under endogenous promoter and enhancer sequences, greatly facilitating the study of diverse neuronal subtypes in the cerebral cortex. However, insertion of exogenous DNA into the genome can have unintended effects on local gene regulation or protein function that must be carefully considered. Here, we analyze a recently generated Tbr1-2A-CreER knock-in mouse line, where a 2A-CreER cassette was inserted in-frame just before the stop codon of the transcription factor gene Tbr1 . Heterozygous TBR1 mutations in humans and mice are known to cause autism or autism-like behavioral phenotypes accompanied by structural brain malformations, most frequently a reduction of the anterior commissure. Thus, it is critical for modified versions of Tbr1 to exhibit true wild-type-like activity. We evaluated the Tbr1-2A-CreER allele for its potential impact on Tbr1 function and complementation to Tbr1 loss-of-function alleles. In mice with one copy of the Tbr1-2A-CreER allele, we identified reduction of TBR1 protein in early postnatal cortex along with thinning of the anterior commissure, suggesting hypersensitivity of this structure to TBR1 dosage. Comparing Tbr1-2A-CreER and Tbr1 -null heterozygous and homozygous mice to Tbr1 -null complementation crosses showed reductions of TBR1 dosage ranging from 28.4% to 95.9%. Using these combinatorial genotypes, we found that low levels of TBR1 protein (∼16%) are sufficient to establish cortical layer positioning, while greater levels (>50%) are required for normal suppression of layer 5 identity. In total, these results strongly support the conclusion that Tbr1-2A-CreER is a hypomorphic allele. We advise caution when interpreting experiments using this allele, such as transcriptomic studies, considering the sensitivity of various corticogenic processes to TBR1 dosage and the association of heterozygous TBR1 mutations with complex neurodevelopmental disorders.

Detection of autism spectrum disorder-related pathogenic trio variants by a novel structure-based approach

BACKGROUND

Glutamatergic synapse dysfunction is believed to underlie the development of Autism Spectrum Disorder (ASD) and Intellectual Disability (ID) in many individuals. However, identification of genetic markers that contribute to synaptic dysfunction in these individuals is notoriously difficult. Based on genomic analysis, structural modeling, and functional data, we recently established the involvement of the TRIO-RAC1 pathway in ASD and ID. Furthermore, we identified a pathological de novo missense mutation hotspot in TRIO's GEF1 domain. ASD/ID-related missense mutations within this domain compromise glutamatergic synapse function and likely contribute to the development of ASD/ID. The number of ASD/ID cases with mutations identified within TRIO's GEF1 domain is increasing. However, tools for accurately predicting whether such mutations are detrimental to protein function are lacking.

METHODS

Here we deployed advanced protein structural modeling techniques to predict potential de novo pathogenic and benign mutations within TRIO's GEF1 domain. Mutant TRIO-9 constructs were generated and expressed in CA1 pyramidal neurons of organotypic cultured hippocampal slices. AMPA receptor-mediated postsynaptic currents were examined in these neurons using dual whole-cell patch clamp electrophysiology. We also validated these findings using orthogonal co-immunoprecipitation and fluorescence lifetime imaging (FLIM-FRET) experiments to assay TRIO mutant overexpression effects on TRIO-RAC1 binding and on RAC1 activity in HEK293/T cells.

RESULTS

Missense mutations in TRIO's GEF1 domain that were predicted to disrupt TRIO-RAC1 binding or stability were tested experimentally and found to greatly impair TRIO-9's influence on glutamatergic synapse function. In contrast, missense mutations in TRIO's GEF1 domain that were predicted to have minimal effect on TRIO-RAC1 binding or stability did not impair TRIO-9's influence on glutamatergic synapse function in our experimental assays. In orthogonal assays, we find most of the mutations predicted to disrupt binding display loss of function but mutants predicted to disrupt stability do not reflect our results from neuronal electrophysiological data.

LIMITATIONS

We present a method to predict missense mutations in TRIO's GEF1 domain that may compromise TRIO function and test for effects in a limited number of assays. Possible limitations arising from the model systems employed here can be addressed in future studies. Our method does not provide evidence for whether these mutations confer ASD/ID risk or the likelihood that such mutations will result in the development of ASD/ID.

CONCLUSIONS

Here we show that a combination of structure-based computational predictions and experimental validation can be employed to reliably predict whether missense mutations in the human TRIO gene impede TRIO protein function and compromise TRIO's role in glutamatergic synapse regulation. With the growing accessibility of genome sequencing, the use of such tools in the accurate identification of pathological mutations will be instrumental in diagnostics of ASD/ID.

Genomic analysis of 116 autism families strengthens known risk genes and highlights promising candidates

Autism spectrum disorder (ASD) is a complex neurodevelopmental condition with a strong genetic component in which rare variants contribute significantly to risk. We performed whole genome and/or exome sequencing (WGS and WES) and SNP-array analysis to identify both rare sequence and copy number variants (SNVs and CNVs) in 435 individuals from 116 ASD families. We identified 37 rare potentially damaging de novo SNVs (pdSNVs) in the cases (n = 144). Interestingly, two of them (one stop-gain and one missense variant) occurred in the same gene, BRSK2. Moreover, the identification of 8 severe de novo pdSNVs in genes not previously implicated in ASD (AGPAT3, IRX5, MGAT5B, RAB8B, RAP1A, RASAL2, SLC9A1, YME1L1) highlighted promising candidates. Potentially damaging CNVs (pdCNVs) provided support to the involvement of inherited variants in PHF3, NEGR1, TIAM1 and HOMER1 in neurodevelopmental disorders (NDD), although mostly acting as susceptibility factors with incomplete penetrance. Interpretation of identified pdSNVs/pdCNVs according to the ACMG guidelines led to a molecular diagnosis in 19/144 cases, although this figure represents a lower limit and is expected to increase thanks to further clarification of the role of likely pathogenic variants in ASD/NDD candidate genes not yet established. In conclusion, our study highlights promising ASD candidate genes and contributes to characterize the allelic diversity, mode of inheritance and phenotypic impact of de novo and inherited risk variants in ASD/NDD genes.

Identifying associations of de novo noncoding variants with autism through integration of gene expression, sequence and sex information

Whole-genome sequencing (WGS) data is facilitating genome-wide identification of rare noncoding variants, while elucidating their roles in disease remains challenging. Towards this end, we first revisit a reported significant brain-related association signal of autism spectrum disorder (ASD) detected from de novo noncoding variants attributed to deep-learning and show that local GC content can capture similar association signals. We further show that the association signal appears driven by variants from male proband-female sibling pairs that are upstream of assigned genes. We then develop Expression Neighborhood Sequence Association Study (ENSAS), which utilizes gene expression correlations and sequence information, to more systematically identify phenotype-associated variant sets. Applying ENSAS to the same set of de novo variants, we identify gene expression-based neighborhoods showing significant ASD association signal, enriched for synapse-related gene ontology terms. For these top neighborhoods, we also identify chromatin states annotations of variants that are predictive of the proband-sibling local GC content differences. Our work provides new insights into associations of non-coding de novo mutations in ASD and presents an analytical framework applicable to other phenotypes.

Case report: A familial B-acute lymphoblastic leukemia associated with a new germline pathogenic variant in PAX5. The first report in Mexico

B-cell acute lymphoblastic leukemia (B-ALL) is one of the most common childhood cancers worldwide. Although most cases are sporadic, some familial forms, inherited as autosomal dominant traits with incomplete penetrance, have been described over the last few years. Germline pathogenic variants in transcription factors such as PAX5, IKZF1, and ETV6 have been identified as causal in familial forms. The proband was a 7-year-old Mexican girl diagnosed with high-risk B-ALL at five years and 11 months of age. Family history showed that the proband's mother had high-risk B-ALL at 16 months of age. She received chemotherapy and was discharged at nine years of age without any evidence of recurrence of leukemia. The proband's father was outside the family nucleus, but no history of leukemia or cancer was present up to the last contact with the mother. We performed exome sequencing on the proband and the proband's mother and identified the PAX5 variant NM_016734.3:c.963del: p.(Ala322LeufsTer11), located in the transactivation domain of the PAX5 protein. The variant was classified as probably pathogenic according to the ACMG criteria. To the best of our knowledge, this is the first Mexican family with an inherited increased risk of childhood B-ALL caused by a novel germline pathogenic variant of PAX5. Identifying individuals with a hereditary predisposition to cancer is essential for modern oncological practice. Individuals at high risk of leukemia would benefit from hematopoietic stem cell transplantation, but family members carrying the pathogenic variant should be excluded as hematopoietic stem cell donors.

Topologically associating domains define the impact of de novo promoter variants on autism spectrum disorder risk

Whole-genome sequencing (WGS) studies of autism spectrum disorder (ASD) have demonstrated the roles of rare promoter de novo variants (DNVs). However, most promoter DNVs in ASD are not located immediately upstream of known ASD genes. In this study analyzing WGS data of 5,044 ASD probands, 4,095 unaffected siblings, and their parents, we show that promoter DNVs within topologically associating domains (TADs) containing ASD genes are significantly and specifically associated with ASD. An analysis considering TADs as functional units identified specific TADs enriched for promoter DNVs in ASD and indicated that common variants in these regions also confer ASD heritability. Experimental validation using human induced pluripotent stem cells (iPSCs) showed that likely deleterious promoter DNVs in ASD can influence multiple genes within the same TAD, resulting in overall dysregulation of ASD-associated genes. These results highlight the importance of TADs and gene-regulatory mechanisms in better understanding the genetic architecture of ASD.

The Genotype and Phenotypes in Families (GPF) platform manages the large and complex data at SFARI

The exploration of genotypic variants impacting phenotypes is a cornerstone in genetics research. The emergence of vast collections containing deeply genotyped and phenotyped families has made it possible to pursue the search for variants associated with complex diseases. However, managing these large-scale datasets requires specialized computational tools tailored to organize and analyze the extensive data. GPF (Genotypes and Phenotypes in Families) is an open-source platform ( https://github.com/iossifovlab/gpf ) that manages genotypes and phenotypes derived from collections of families. The GPF interface allows interactive exploration of genetic variants, enrichment analysis for de novo mutations, and phenotype/genotype association tools. In addition, GPF allows researchers to share their data securely with the broader scientific community. GPF is used to disseminate two large-scale family collection datasets (SSC, SPARK) for the study of autism funded by the SFARI foundation. However, GPF is versatile and can manage genotypic data from other small or large family collections. Our GPF-SFARI GPF instance ( https://gpf.sfari.org/ ) provides protected access to comprehensive genotypic and phenotypic data for the SSC and SPARK. In addition, GPF-SFARI provides public access to an extensive collection of de novo mutations identified in individuals with autism and related disorders and to gene-level statistics of the protected datasets characterizing the genes' roles in autism. Here, we highlight the primary features of GPF within the context of GPF-SFARI.

Factors Associated With Underutilization of Genetic Testing in Autism Spectrum Disorders

BACKGROUND

We sought to identify patient and provider factors associated with low completion of genetic testing, specifically chromosomal microarray (CMA), for autism spectrum disorder (ASD).

METHODS

Medical record review was conducted of children newly diagnosed with ASD without prior genetic testing at a single academic medical center from February 2015 through January 2016.

RESULTS

Only 41.9% of individuals with ASD completed CMA testing over at least 18 months from diagnosis (n = 140 of 334). Time to CMA completion varied, with a median of 86.5 days (interquartile range 2 to 214.5 days). Provider recommendation of genetic testing at the diagnostic visit and greater number of follow-up visits were associated with CMA completion. On multivariate regression, CMA completion was inversely associated with age (odds ratio [OR] = 0.8 for each year older, 95% confidence interval [CI] 0.7, 0.9; P = 0.001) and directly associated with intellectual disability or global developmental delay (OR = 2.2, 95% CI 1.3, 3.8; P = 0.004), first-degree relative with ASD (OR = 2.5, 95% CI 1.0, 6.0; P = 0.044), and public insurance (OR = 1.7, 95% CI 1.0, 2.9; P = 0.037). Parental concern and cost/insurance coverage were the most frequently documented barriers.

CONCLUSIONS

Workflows to support early genetic testing recommendation and ordering soon after diagnosis may increase utilization, incorporating both family and provider perspectives. Genetic counseling highlighting the utility of genetic testing across the life span, phenotypic variability of genetic disorders, and possibility of de novo variants in ASD may also improve utilization.

Autism Spectrum Disorder Detection by Hybrid Convolutional Recurrent Neural Networks from Structural and Resting State Functional MRI Images

This study aims to increase the accuracy of autism spectrum disorder (ASD) diagnosis based on cognitive and behavioral phenotypes through multiple neuroimaging modalities. We apply machine learning (ML) algorithms to classify ASD patients and healthy control (HC) participants using structural magnetic resonance imaging (s-MRI) together with resting state functional MRI (rs-f-MRI and f-MRI) data from the large multisite data repository ABIDE (autism brain imaging data exchange) and identify important brain connectivity features. The 2D f-MRI images were converted into 3D s-MRI images, and datasets were preprocessed using the Montreal Neurological Institute (MNI) atlas. The data were then denoised to remove any confounding factors. We show, by using three fusion strategies such as early fusion, late fusion, and cross fusion, that, in this implementation, hybrid convolutional recurrent neural networks achieve better performance in comparison to either convolutional neural networks (CNNs) or recurrent neural networks (RNNs). The proposed model classifies subjects as autistic or not according to how functional and anatomical connectivity metrics provide an overall diagnosis based on the autism diagnostic observation schedule (ADOS) standard. Our hybrid network achieved an accuracy of 96% by fusing s-MRI and f-MRI together, which outperforms the methods used in previous studies.

Development and Psychometric Evaluation of the Open-Source Challenging Behavior Scale (OS-CBS)

At present, there are no brief, freely-available, informant-report measures that evaluate key challenging behaviors relevant to youth with autism spectrum disorder (ASD) or other developmental disabilities (DD). This paper describes the development, refinement, and initial psychometric evaluation of a new 18-item measure, the Open-Source Challenging Behavior Scale (OS-CBS). In a large sample (n = 2004, 169 with ASD, ages 2-17), results of psychometric analyses indicated a clear factor structure (property destruction, aggression, elopement, conduct problems, and self-injury and a general factor with high loadings from all items) based on exploratory structural equation modeling, good scale reliability (α = .66-.83 for subscales, α = .91 total scale), measurement invariance across demographics, and good construct validity. The OS-CBS is a psychometrically-sound instrument for screening and monitoring intervention progress.

CellGO: a novel deep learning-based framework and webserver for cell-type-specific gene function interpretation

Interpreting the function of genes and gene sets identified from omics experiments remains a challenge, as current pathway analysis tools often fail to consider the critical biological context, such as tissue or cell-type specificity. To address this limitation, we introduced CellGO. CellGO tackles this challenge by leveraging the visible neural network (VNN) and single-cell gene expressions to mimic cell-type-specific signaling propagation along the Gene Ontology tree within a cell. This design enables a novel scoring system to calculate the cell-type-specific gene-pathway paired active scores, based on which, CellGO is able to identify cell-type-specific active pathways associated with single genes. In addition, by aggregating the activities of single genes, CellGO extends its capability to identify cell-type-specific active pathways for a given gene set. To enhance biological interpretation, CellGO offers additional features, including the identification of significantly active cell types and driver genes and community analysis of pathways. To validate its performance, CellGO was assessed using a gene set comprising mixed cell-type markers, confirming its ability to discern active pathways across distinct cell types. Subsequent benchmarking analyses demonstrated CellGO's superiority in effectively identifying cell types and their corresponding cell-type-specific pathways affected by gene knockouts, using either single genes or sets of genes differentially expressed between knockout and control samples. Moreover, CellGO demonstrated its ability to infer cell-type-specific pathogenesis for disease risk genes. Accessible as a Python package, CellGO also provides a user-friendly web interface, making it a versatile and accessible tool for researchers in the field.

Dynamic Regulation of brsk2 in the Social and Motor Development of Zebrafish: A Developmental Behavior Analysis

Both social and motor development play an essential role in an individual's physical, psychological, and social well-being. It is essential to conduct a dynamic analysis at multiple time points during the developmental process as it helps us better understand and evaluate the trajectory and changes in individual development. Recently, some studies found that mutations in the BRSK2 gene may contribute to motor impairments, delays in achieving motor milestones, and deficits in social behavior and communication skills in patients. However, little is known about the dynamic analysis of social and motor development at multiple time points during the development of the brsk2 gene. We generated a novel brsk2-deficient (brsk2ab-/-) zebrafish model through CRISPR/Cas9 editing and conducted comprehensive morphological and neurobehavioral evaluations, including that of locomotor behaviors, social behaviors, and anxiety behaviors from the larval to adult stages of development. Compared to wild-type zebrafish, brsk2ab-/- zebrafish exhibited a catch-up growth pattern of body length and gradually improved locomotor activities during the developmental process. In contrast, multimodal behavior tests showed that the brsk2ab-/- zebrafish displayed escalating social deficiency and anxiety-like behaviors over time. We reported for the first time that the brsk2 gene had dynamic regulatory effects on motor and social development. It helps us understand developmental trends, capture changes, facilitate early interventions, and provide personalized support and development opportunities for individuals.

NGS Custom Panel Implementation in Patients with Non-Syndromic Autism Spectrum Disorders in the Clinical Routine of a Tertiary Hospital

Autism spectrum disorder (ASD) is a set of neurodevelopmental disorders characterized by deficiencies in communication, social interaction, and repetitive and restrictive behaviors. The discovery of genetic involvement in the etiology of ASD has made this condition a strong candidate for genome-based diagnostic tests. Next-generation sequencing (NGS) is useful for the detection of variants in the sequence of different genes in ASD patients. Herein, we present the implementation of a personalized NGS panel for autism (AutismSeq) for patients with essential ASD over a prospective period of four years in the clinical routine of a tertiary hospital. The cohort is composed of 48 individuals, older than 3 years, who met the DSM-5 (The Diagnostic and Statistical Manual of Mental Disorders) diagnostic criteria for ASD. The NGS customized panel (AutismSeq) turned out to be a tool with good diagnostic efficacy in routine clinical care, where we detected 12 "pathogenic" (including pathogenic, likely pathogenic, and VUS (variant of uncertain significance) possibly pathogenic variations) in 11 individuals, and 11 VUS in 10 individuals, which had previously been negative for chromosomal microarray analysis and other previous genetic studies, such as karyotype, fragile-X, or MLPA/FISH (Multiplex Ligation dependent Probe Amplification/Fluorescence in situ hybridization) analysis. Our results demonstrate the high genetic and clinical heterogeneity of individuals with ASD and the current difficulty of molecular diagnosis. Our study also shows that an NGS-customized panel might be useful for diagnosing patients with essential/primary autism and that it is cost-effective for most genetic laboratories.

Genetic Ablation of GIGYF1, Associated With Autism, Causes Behavioral and Neurodevelopmental Defects in Zebrafish and Mice

BACKGROUND

Autism spectrum disorder is characterized by deficits in social communication and restricted or repetitive behaviors. Due to the extremely high genetic and phenotypic heterogeneity, it is critical to pinpoint the genetic factors for understanding the pathology of these disorders.

METHODS

We analyzed the exomes generated by the SPARK (Simons Powering Autism Research) project and performed a meta-analysis with previous data. We then generated 1 zebrafish knockout model and 3 mouse knockout models to examine the function of GIGYF1 in neurodevelopment and behavior. Finally, we performed whole tissue and single-nucleus transcriptome analysis to explore the molecular and cellular function of GIGYF1.

RESULTS

GIGYF1 variants are significantly associated with various neurodevelopmental disorder phenotypes, including autism, global developmental delay, intellectual disability, and sleep disturbance. Loss of GIGYF1 causes similar behavioral effects in zebrafish and mice, including elevated levels of anxiety and reduced social engagement, which is reminiscent of the behavioral deficits in human patients carrying GIGYF1 variants. Moreover, excitatory neuron-specific Gigyf1 knockout mice recapitulate the increased repetitive behaviors and impaired social memory, suggesting a crucial role of Gigyf1 in excitatory neurons, which correlates with the observations in single-nucleus RNA sequencing. We also identified a series of downstream target genes of GIGYF1 that affect many aspects of the nervous system, especially synaptic transmission.

CONCLUSIONS

De novo variants of GIGYF1 are associated with neurodevelopmental disorders, including autism spectrum disorder. GIGYF1 is involved in neurodevelopment and animal behavior, potentially through regulating hippocampal CA2 neuronal numbers and disturbing synaptic transmission.

Whole genome analysis of rare deleterious variants adds further evidence to BRSK2 and other risk genes in Autism Spectrum Disorder

Autism spectrum disorder (ASD) is a complex neurodevelopmental condition with a strong genetic component in which rare variants contribute significantly to risk. We have performed whole genome and/or exome sequencing (WGS and WES) and SNP-array analysis to identify both rare sequence and copy number variants (SNVs and CNVs) in 435 individuals from 116 ASD families. We identified 37 rare potentially damaging de novo SNVs (pdSNVs) in cases (n = 144). Interestingly, two of them (one stop-gain and one missense variant) occurred in the same gene, BRSK2. Moreover, the identification of 9 severe de novo pdSNVs in genes not previously implicated in ASD (RASAL2, RAP1A, IRX5, SLC9A1, AGPAT3, MGAT3, RAB8B, MGAT5B, YME1L1), highlighted novel candidates. Potentially damaging CNVs (pdCNVs) provided support to the involvement of inherited variants in PHF3, NEGR1, TIAM1 and HOMER1 in neurodevelopmental disorders (NDD), although mostly acting as susceptibility factors with incomplete penetrance. Interpretation of identified pdSNVs/pdCNVs according to the ACMG guidelines led to a molecular diagnosis in 19/144 cases, but this figure represents a lower limit and is expected to increase thanks to further clarification of the role of likely pathogenic variants in new ASD/NDD candidates. In conclusion, our study strengthens the role of BRSK2 and other neurodevelopmental genes in ASD risk, highlights novel candidates and contributes to characterize the allelic diversity, mode of inheritance and phenotypic impact of de novo and inherited risk variants in ASD/NDD genes.

Identification of de novo Mutations in the Chinese Autism Spectrum Disorder Cohort via Whole-Exome Sequencing Unveils Brain Regions Implicated in Autism

Autism spectrum disorder (ASD) is a highly heritable neurodevelopmental disorder characterized by deficits in social interactions and repetitive behaviors. Although hundreds of ASD risk genes, implicated in synaptic formation and transcriptional regulation, have been identified through human genetic studies, the East Asian ASD cohorts are still under-represented in genome-wide genetic studies. Here, we applied whole-exome sequencing to 369 ASD trios including probands and unaffected parents of Chinese origin. Using a joint-calling analytical pipeline based on GATK toolkits, we identified numerous de novo mutations including 55 high-impact variants and 165 moderate-impact variants, as well as de novo copy number variations containing known ASD-related genes. Importantly, combined with single-cell sequencing data from the developing human brain, we found that the expression of genes with de novo mutations was specifically enriched in the pre-, post-central gyrus (PRC, PC) and banks of the superior temporal (BST) regions in the human brain. By further analyzing the brain imaging data with ASD and healthy controls, we found that the gray volume of the right BST in ASD patients was significantly decreased compared to healthy controls, suggesting the potential structural deficits associated with ASD. Finally, we found a decrease in the seed-based functional connectivity between BST/PC/PRC and sensory areas, the insula, as well as the frontal lobes in ASD patients. This work indicated that combinatorial analysis with genome-wide screening, single-cell sequencing, and brain imaging data reveal the brain regions contributing to the etiology of ASD.

Case report: A novel frameshift mutation in BRSK2 causes autism in a 16-year old Chinese boy

Serine/threonine protein kinases are involved in axon formation and neuronal polarization and have recently been implicated in autism spectrum disorder (ASD) and neurodevelopmental disorders (NDD). Here, we focus on BRSK2, which encodes brain-specific serine/threonine protein kinase 2. Although previous studies have reported 19 unrelated patients with BRSK2 pathogenic variation, only 15 of 19 patients have detailed clinical data. Therefore, more case reports are needed to enrich the phenotype associated with BRSK2 mutations. In this study, we report a novel de novo frameshift variant (c.442del, p.L148Cfs*39) identified by exome sequencing in a 16 year-old Chinese boy with ASD. The proband presented with attention-deficit, auditory hallucinations, limb tremor, and abnormal brain electrical activity mapping. This study expands the phenotypic spectrum of BRSK2-related cases and reveals the highly variable severity of disorders associated with BRSK2.

The contributions of rare inherited and polygenic risk to ASD in multiplex families

Autism spectrum disorder (ASD) has a complex genetic architecture involving contributions from both de novo and inherited variation. Few studies have been designed to address the role of rare inherited variation or its interaction with common polygenic risk in ASD. Here, we performed whole-genome sequencing of the largest cohort of multiplex families to date, consisting of 4,551 individuals in 1,004 families having two or more autistic children. Using this study design, we identify seven previously unrecognized ASD risk genes supported by a majority of rare inherited variants, finding support for a total of 74 genes in our cohort and a total of 152 genes after combined analysis with other studies. Autistic children from multiplex families demonstrate an increased burden of rare inherited protein-truncating variants in known ASD risk genes. We also find that ASD polygenic score (PGS) is overtransmitted from nonautistic parents to autistic children who also harbor rare inherited variants, consistent with combinatorial effects in the offspring, which may explain the reduced penetrance of these rare variants in parents. We also observe that in addition to social dysfunction, language delay is associated with ASD PGS overtransmission. These results are consistent with an additive complex genetic risk architecture of ASD involving rare and common variation and further suggest that language delay is a core biological feature of ASD.

Phenotypic effects of genetic variants associated with autism

While over 100 genes have been associated with autism, little is known about the prevalence of variants affecting them in individuals without a diagnosis of autism. Nor do we fully appreciate the phenotypic diversity beyond the formal autism diagnosis. Based on data from more than 13,000 individuals with autism and 210,000 undiagnosed individuals, we estimated the odds ratios for autism associated to rare loss-of-function (LoF) variants in 185 genes associated with autism, alongside 2,492 genes displaying intolerance to LoF variants. In contrast to autism-centric approaches, we investigated the correlates of these variants in individuals without a diagnosis of autism. We show that these variants are associated with a small but significant decrease in fluid intelligence, qualification level and income and an increase in metrics related to material deprivation. These effects were larger for autism-associated genes than in other LoF-intolerant genes. Using brain imaging data from 21,040 individuals from the UK Biobank, we could not detect significant differences in the overall brain anatomy between LoF carriers and non-carriers. Our results highlight the importance of studying the effect of the genetic variants beyond categorical diagnosis and the need for more research to understand the association between these variants and sociodemographic factors, to best support individuals carrying these variants.

Exploratory analysis of L1 retrotransposons expression in autism

BACKGROUND

Autism spectrum disorder (ASD) is a set of highly heterogeneous neurodevelopmental diseases whose genetic etiology is not completely understood. Several investigations have relied on transcriptome analysis from peripheral tissues to dissect ASD into homogenous molecular phenotypes. Recently, analysis of changes in gene expression from postmortem brain tissues has identified sets of genes that are involved in pathways previously associated with ASD etiology. In addition to protein-coding transcripts, the human transcriptome is composed by a large set of non-coding RNAs and transposable elements (TEs). Advancements in sequencing technologies have proven that TEs can be transcribed in a regulated fashion, and their dysregulation might have a role in brain diseases.

METHODS

We exploited published datasets comprising RNA-seq data from (1) postmortem brain of ASD subjects, (2) in vitro cell cultures where ten different ASD-relevant genes were knocked out and (3) blood of discordant siblings. We measured the expression levels of evolutionarily young full-length transposable L1 elements and characterized the genomic location of deregulated L1s assessing their potential impact on the transcription of ASD-relevant genes. We analyzed every sample independently, avoiding to pool together the disease subjects to unmask the heterogeneity of the molecular phenotypes.

RESULTS

We detected a strong upregulation of intronic full-length L1s in a subset of postmortem brain samples and in in vitro differentiated neurons from iPSC knocked out for ATRX. L1 upregulation correlated with an high number of deregulated genes and retained introns. In the anterior cingulate cortex of one subject, a small number of significantly upregulated L1s overlapped with ASD-relevant genes that were significantly downregulated, suggesting the possible existence of a negative effect of L1 transcription on host transcripts.

LIMITATIONS

Our analyses must be considered exploratory and will need to be validated in bigger cohorts. The main limitation is given by the small sample size and by the lack of replicates for postmortem brain samples. Measuring the transcription of locus-specific TEs is complicated by the repetitive nature of their sequence, which reduces the accuracy in mapping sequencing reads to the correct genomic locus.

CONCLUSIONS

L1 upregulation in ASD appears to be limited to a subset of subjects that are also characterized by a general deregulation of the expression of canonical genes and an increase in intron retention. In some samples from the anterior cingulate cortex, L1s upregulation seems to directly impair the expression of some ASD-relevant genes by a still unknown mechanism. L1s upregulation may therefore identify a group of ASD subjects with common molecular features and helps stratifying individuals for novel strategies of therapeutic intervention.

Loss-of-function variants in CUL3 cause a syndromic neurodevelopmental disorder

Purpose

De novo variants in CUL3 (Cullin-3 ubiquitin ligase) have been strongly associated with neurodevelopmental disorders (NDDs), but no large case series have been reported so far. Here we aimed to collect sporadic cases carrying rare variants in CUL3, describe the genotype-phenotype correlation, and investigate the underlying pathogenic mechanism.

Methods

Genetic data and detailed clinical records were collected via multi-center collaboration. Dysmorphic facial features were analyzed using GestaltMatcher. Variant effects on CUL3 protein stability were assessed using patient-derived T-cells.

Results

We assembled a cohort of 35 individuals with heterozygous CUL3 variants presenting a syndromic NDD characterized by intellectual disability with or without autistic features. Of these, 33 have loss-of-function (LoF) and two have missense variants. CUL3 LoF variants in patients may affect protein stability leading to perturbations in protein homeostasis, as evidenced by decreased ubiquitin-protein conjugates in vitro . Specifically, we show that cyclin E1 (CCNE1) and 4E-BP1 (EIF4EBP1), two prominent substrates of CUL3, fail to be targeted for proteasomal degradation in patient-derived cells.

Conclusion

Our study further refines the clinical and mutational spectrum of CUL3 -associated NDDs, expands the spectrum of cullin RING E3 ligase-associated neuropsychiatric disorders, and suggests haploinsufficiency via LoF variants is the predominant pathogenic mechanism.

Examining indicators of psychosocial risk and resilience in parents of autistic children

Background

Parents of autistic children experience increased levels of caregiver strain and adverse mental health outcomes, even in comparison to parents of children with other neurodevelopmental disabilities. Previous studies have largely attributed these increased levels of mental health concerns to their child behavioral concerns and autism symptomatology, but less attention has been given to other potential child factors, such as child adaptive functioning. Additionally, little is known about potential protective factors, such as parents' emotion regulation (ER) abilities, that may ameliorate the experience of caregiver strain, anxiety, and depression.

Objective

The current study examined the impact of child characteristics (restricted and repetitive behaviors, adaptive functioning and behavioral concerns) on parent mental health outcomes (caregiver strain, anxiety, depression and wellbeing). Additionally, we explore parents' ER abilities as a moderator of the impact of child characteristic on parents' mental health outcomes.

Results

Results of linear mixed effect models indicated a significant relationship between parents' ER abilities and all four parent outcomes. Additionally, children's adaptive functioning abilities and repetitive behaviors (RRBs) were significant predictors of caregiving strain. Parents' ER abilities were a significant moderator of the effect of children's repetitive behaviors and adaptive functioning challenges on caregiver strain, such that better ER abilities mitigated the impact of child clinical factors on caregiver strain. Finally, a significant difference was detected for mothers' and fathers' mental health, with mothers reporting higher caregiver strain, and more symptoms of anxiety and depression than did fathers.

Conclusion

This study leveraged a large sample of autistic children and their biological parents to examine the relationship between children's clinical characteristics and parents' psychological wellbeing. Results indicate that, although parents of autistic children do experience high rates of internalizing mental health concerns that relate to child adaptive functioning and RRBs, parent ER abilities act as a protective factor against parents' adverse mental health outcomes. Further, mothers in our sample reported significantly higher rates of depression, anxiety, and caregiver strain, as compared with fathers.

Targeted sequencing and clinical strategies in children with autism spectrum disorder: A cohort study

Objectives: Autism spectrum disorder (ASD) is a neurodevelopmental disorder with genetic and clinical heterogeneity. Owing to the advancement of sequencing technologies, an increasing number of ASD-related genes have been reported. We designed a targeted sequencing panel (TSP) for ASD based on next-generation sequencing (NGS) to provide clinical strategies for genetic testing of ASD and its subgroups.

Methods: TSP comprised 568 ASD-related genes and analyzed both single nucleotide variations (SNVs) and copy number variations (CNVs). The Autism Diagnostic Observation Schedule (ADOS) and the Griffiths Mental Development Scales (GMDS) were performed with the consent of ASD parents. Additional medical information of the selected cases was recorded.

Results: A total of 160 ASD children were enrolled in the cohort (male to female ratio 3.6:1). The total detection yield was 51.3% for TSP (82/160), among which SNVs and CNVs accounted for 45.6% (73/160) and 8.1% (13/160), respectively, with 4 children having both SNVs and CNV variants (2.5%). The detection rate of disease-associated variants in females (71.4%) was significantly higher than that in males (45.6%, p = 0.007). Pathogenic and likely pathogenic variants were detected in 16.9% (27/160) of the cases. SHANK3, KMT2A, and DLGAP2 were the most frequent variants among these patients. Eleven children had de novo SNVs, 2 of whom had de novo ASXL3 variants with mild global developmental delay (DD) and minor dysmorphic facial features besides autistic symptoms. Seventy-one children completed both ADOS and GMDS, of whom 51 had DD/intellectual disability (ID). In this subgroup of ASD children with DD/ID, we found that children with genetic abnormalities had lower language competence than those without positive genetic findings (p = 0.028). There was no correlation between the severity of ASD and positive genetic findings.

Conclusion: Our study revealed the potential of TSP, with lower cost and more efficient genetic diagnosis. We recommended that ASD children with DD or ID, especially those with lower language competence, undergo genetic testing. More precise clinical phenotypes may help in the decision-making of patients with genetic testing.

A predictive ensemble classifier for the gene expression diagnosis of ASD at ages 1 to 4 years

Autism Spectrum Disorder (ASD) diagnosis remains behavior-based and the median age of diagnosis is ~52 months, nearly 5 years after its first-trimester origin. Accurate and clinically-translatable early-age diagnostics do not exist due to ASD genetic and clinical heterogeneity. Here we collected clinical, diagnostic, and leukocyte RNA data from 240 ASD and typically developing (TD) toddlers (175 toddlers for training and 65 for test). To identify gene expression ASD diagnostic classifiers, we developed 42,840 models composed of 3570 gene expression feature selection sets and 12 classification methods. We found that 742 models had AUC-ROC ≥ 0.8 on both Training and Test sets. Weighted Bayesian model averaging of these 742 models yielded an ensemble classifier model with accurate performance in Training and Test gene expression datasets with ASD diagnostic classification AUC-ROC scores of 85-89% and AUC-PR scores of 84-92%. ASD toddlers with ensemble scores above and below the overall ASD ensemble mean of 0.723 (on a scale of 0 to 1) had similar diagnostic and psychometric scores, but those below this ASD ensemble mean had more prenatal risk events than TD toddlers. Ensemble model feature genes were involved in cell cycle, inflammation/immune response, transcriptional gene regulation, cytokine response, and PI3K-AKT, RAS and Wnt signaling pathways. We additionally collected targeted DNA sequencing smMIPs data on a subset of ASD risk genes from 217 of the 240 ASD and TD toddlers. This DNA sequencing found about the same percentage of SFARI Level 1 and 2 ASD risk gene mutations in TD (12 of 105) as in ASD (13 of 112) toddlers, and classification based only on the presence of mutation in these risk genes performed at a chance level of 49%. By contrast, the leukocyte ensemble gene expression classifier correctly diagnostically classified 88% of TD and ASD toddlers with ASD risk gene mutations. Our ensemble ASD gene expression classifier is diagnostically predictive and replicable across different toddler ages, races, and ethnicities; out-performs a risk gene mutation classifier; and has potential for clinical translation.

SHINE: protein language model-based pathogenicity prediction for short inframe insertion and deletion variants

Accurate variant pathogenicity predictions are important in genetic studies of human diseases. Inframe insertion and deletion variants (indels) alter protein sequence and length, but not as deleterious as frameshift indels. Inframe indel Interpretation is challenging due to limitations in the available number of known pathogenic variants for training. Existing prediction methods largely use manually encoded features including conservation, protein structure and function, and allele frequency to infer variant pathogenicity. Recent advances in deep learning modeling of protein sequences and structures provide an opportunity to improve the representation of salient features based on large numbers of protein sequences. We developed a new pathogenicity predictor for SHort Inframe iNsertion and dEletion (SHINE). SHINE uses pretrained protein language models to construct a latent representation of an indel and its protein context from protein sequences and multiple protein sequence alignments, and feeds the latent representation into supervised machine learning models for pathogenicity prediction. We curated training data from ClinVar and gnomAD, and created two test datasets from different sources. SHINE achieved better prediction performance than existing methods for both deletion and insertion variants in these two test datasets. Our work suggests that unsupervised protein language models can provide valuable information about proteins, and new methods based on these models can improve variant interpretation in genetic analyses.