Non-coding de novo mutations in chromatin interactions are implicated in autism spectrum disorder

Mitochondrial DNA mosaicism in normal human somatic cells

Somatic cells accumulate genomic alterations with age; however, our understanding of mitochondrial DNA (mtDNA) mosaicism remains limited. Here we investigated the genomes of 2,096 clones derived from three cell types across 31 donors, identifying 6,451 mtDNA variants with heteroplasmy levels of ≳0.3%. While the majority of these variants were unique to individual clones, suggesting stochastic acquisition with age, 409 variants (6%) were shared across multiple embryonic lineages, indicating their origin from heteroplasmy in fertilized eggs. The mutational spectrum exhibited replication-strand bias, implicating mtDNA replication as a major mutational process. We evaluated the mtDNA mutation rate (5.0 × 10-8 per base pair) and a turnover frequency of 10-20 per year, which are fundamental components shaping the landscape of mtDNA mosaicism over a lifetime. The expansion of mtDNA-truncating mutations toward homoplasmy was substantially suppressed. Our findings provide comprehensive insights into the origins, dynamics and functional consequences of mtDNA mosaicism in human somatic cells.

Low-level brain somatic mutations in exonic regions are collectively implicated in autism with germline mutations in autism risk genes

Low-level somatic mutations in the human brain are implicated in various neurological disorders. The contribution of low-level brain somatic mutations to autism spectrum disorder (ASD), however, remains poorly understood. Here, we performed high-depth exome sequencing with an average read depth of 559.3x in 181 cortical, cerebellar, and peripheral tissue samples to identify brain somatic single nucleotide variants (SNVs) in 24 ASD subjects and 31 controls. We detected ~2.4 brain somatic SNVs per exome per single brain region, with a variant allele frequency (VAF) as low as 0.3%. The mutational profiles, including the number, signature, and type, were not significantly different between the ASD patients and controls. Intriguingly, when considering genes with low-level brain somatic SNVs and ASD risk genes with damaging germline SNVs together, the merged set of genes carrying either somatic or germline SNVs in ASD patients was significantly involved in ASD-associated pathophysiology, including dendrite spine morphogenesis (p = 0.025), mental retardation (p = 0.012), and intrauterine growth retardation (p = 0.012). Additionally, the merged gene set showed ASD-associated spatiotemporal expression in the early and mid-fetal cortex, striatum, and thalamus (all p < 0.05). Patients with damaging mutations in the merged gene set had a greater ASD risk than did controls (odds ratio = 3.92, p = 0.025, 95% confidence interval = 1.12-14.79). The findings of this study suggest that brain somatic SNVs and germline SNVs may collectively contribute to ASD-associated pathophysiology.

Noncoding de novo mutations in SCN2A are associated with autism spectrum disorders

Coding de novo mutations (DNMs) contribute to the risk for autism spectrum disorders (ASD), but the contribution of noncoding DNMs remains relatively unexplored. Here we use whole genome sequencing (WGS) data of 12,411 individuals (including 3,508 probands and 2,218 unaffected siblings) from 3,357 families collected in Simons Foundation Powering Autism Research for Knowledge (SPARK) to detect DNMs associated with ASD, while examining Simons Simplex Collection (SSC) with 6383 individuals from 2274 families to replicate the results. For coding DNMs, SCN2A reached exome-wide significance (p=2.06×10-11) in SPARK. The 618 known dominant ASD genes as a group are strongly enriched for coding DNMs in cases than sibling controls (fold change=1.51, p =1.13×10-5 for SPARK; fold change=1.86, p =2.06×10-9 for SSC). For noncoding DNMs, we used two methods to assess statistical significance: a point-based test that analyzes sites with a Combined Annotation Dependent Depletion (CADD) score ≥15, and a segment-based test that analyzes 1kb genomic segments with segment-specific background mutation rates (inferred from expected rare mutations in Gnocchi genome constraint scores). The point-based test identified SCN2A as marginally significant (p=6.12×10-4) in SPARK, yet segment-based test identified CSMD1, RBFOX1 and CHD13 as exome-wide significant. We did not identify significant enrichment of noncoding DNMs (in all 1kb segments or those with Gnocchi>4) in the 618 known ASD genes as a group in cases than sibling controls. When combining evidence from both coding and noncoding DNMs, we found that SCN2A with 11 coding and 5 noncoding DNMs exhibited the strongest significance (p=4.15×10-13). In summary, we identified both coding and noncoding DNMs in SCN2A associated with ASD, while nominating additional candidates for further examination in future studies.

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.

A foundational atlas of autism protein interactions reveals molecular convergence

Translating high-confidence (hc) autism spectrum disorder (ASD) genes into viable treatment targets remains elusive. We constructed a foundational protein-protein interaction (PPI) network in HEK293T cells involving 100 hcASD risk genes, revealing over 1,800 PPIs (87% novel). Interactors, expressed in the human brain and enriched for ASD but not schizophrenia genetic risk, converged on protein complexes involved in neurogenesis, tubulin biology, transcriptional regulation, and chromatin modification. A PPI map of 54 patient-derived missense variants identified differential physical interactions, and we leveraged AlphaFold-Multimer predictions to prioritize direct PPIs and specific variants for interrogation in Xenopus tropicalis and human forebrain organoids. A mutation in the transcription factor FOXP1 led to reconfiguration of DNA binding sites and altered development of deep cortical layer neurons in forebrain organoids. This work offers new insights into molecular mechanisms underlying ASD and describes a powerful platform to develop and test therapeutic strategies for many genetically-defined conditions.

TrkB-dependent regulation of molecular signaling across septal cell types

The lateral septum (LS), a GABAergic structure located in the basal forebrain, is implicated in social behavior, learning, and memory. We previously demonstrated that expression of tropomyosin kinase receptor B (TrkB) in LS neurons is required for social novelty recognition. To better understand molecular mechanisms by which TrkB signaling controls behavior, we locally knocked down TrkB in LS and used bulk RNA-sequencing to identify changes in gene expression downstream of TrkB. TrkB knockdown induces upregulation of genes associated with inflammation and immune responses, and downregulation of genes associated with synaptic signaling and plasticity. Next, we generated one of the first atlases of molecular profiles for LS cell types using single nucleus RNA-sequencing (snRNA-seq). We identified markers for the septum broadly, and the LS specifically, as well as for all neuronal cell types. We then investigated whether the differentially expressed genes (DEGs) induced by TrkB knockdown map to specific LS cell types. Enrichment testing identified that downregulated DEGs are broadly expressed across neuronal clusters. Enrichment analyses of these DEGs demonstrated that downregulated genes are uniquely expressed in the LS, and associated with either synaptic plasticity or neurodevelopmental disorders. Upregulated genes are enriched in LS microglia, associated with immune response and inflammation, and linked to both neurodegenerative disease and neuropsychiatric disorders. In addition, many of these genes are implicated in regulating social behaviors. In summary, the findings implicate TrkB signaling in the LS as a critical regulator of gene networks associated with psychiatric disorders that display social deficits, including schizophrenia and autism, and with neurodegenerative diseases, including Alzheimer's.

Protocol for the development of joint attention-based subclassification of autism spectrum disorder and validation using multi-modal data

BACKGROUND

Heterogeneity in clinical manifestation and underlying neuro-biological mechanisms are major obstacles to providing personalized interventions for individuals with autism spectrum disorder (ASD). Despite various efforts to unify disparate data modalities and machine learning techniques for subclassification, replicable ASD clusters remain elusive. Our study aims to introduce a novel method, utilizing the objective behavioral biomarker of gaze patterns during joint attention, to subclassify ASD. We will assess whether behavior-based subgrouping yields clinically, genetically, and neurologically distinct ASD groups.

METHODS

We propose a study involving 60 individuals with ASD recruited from a specialized psychiatric clinic to perform joint attention tasks. Through the examination of gaze patterns in social contexts, we will conduct a semi-supervised clustering analysis, yielding two primary clusters: good gaze response group and poor gaze response group. Subsequent comparison will occur across these clusters, scrutinizing neuroanatomical structure and connectivity using structural as well as functional brain imaging studies, genetic predisposition through single nucleotide polymorphism data, and assorted socio-demographic and clinical information.

CONCLUSIONS

The aim of the study is to investigate the discriminative properties and the validity of the joint attention-based subclassification of ASD using multi-modality data.

TRIAL REGISTRATION

Clinical trial, KCT0008530, Registered 16 June 2023, https://cris.nih.go.kr/cris/index/index.do .

TrkB-dependent regulation of molecular signaling across septal cell types

The lateral septum (LS), a GABAergic structure located in the basal forebrain, is implicated in social behavior, learning and memory. We previously demonstrated that expression of tropomyosin kinase receptor B (TrkB) in LS neurons is required for social novelty recognition. To better understand molecular mechanisms by which TrkB signaling controls behavior, we locally knocked down TrkB in LS and used bulk RNA-sequencing to identify changes in gene expression downstream of TrkB. TrkB knockdown induces upregulation of genes associated with inflammation and immune responses, and downregulation of genes associated with synaptic signaling and plasticity. Next, we generated one of the first atlases of molecular profiles for LS cell types using single nucleus RNA-sequencing (snRNA-seq). We identified markers for the septum broadly, and the LS specifically, as well as for all neuronal cell types. We then investigated whether the differentially expressed genes (DEGs) induced by TrkB knockdown map to specific LS cell types. Enrichment testing identified that downregulated DEGs are broadly expressed across neuronal clusters. Enrichment analyses of these DEGs demonstrated that downregulated genes are uniquely expressed in the LS, and associated with either synaptic plasticity or neurodevelopmental disorders. Upregulated genes are enriched in LS microglia, associated with immune response and inflammation, and linked to both neurodegenerative disease and neuropsychiatric disorders. In addition, many of these genes are implicated in regulating social behaviors. In summary, the findings implicate TrkB signaling in the LS as a critical regulator of gene networks associated with psychiatric disorders that display social deficits, including schizophrenia and autism, and with neurodegenerative diseases, including Alzheimer's.

Estimation of intrafamilial DNA contamination in family trio genome sequencing using deviation from Mendelian inheritance

With the increasing number of sequencing projects involving families, quality control tools optimized for family genome sequencing are needed. However, accurately quantifying contamination in a DNA mixture is particularly difficult when genetically related family members are the sources. We developed TrioMix, a maximum likelihood estimation (MLE) framework based on Mendel's law of inheritance, to quantify DNA mixture between family members in genome sequencing data of parent-offspring trios. TrioMix can accurately deconvolute any intrafamilial DNA contamination, including parent-offspring, sibling-sibling, parent-parent, and even multiple familial sources. In addition, TrioMix can be applied to detect genomic abnormalities that deviate from Mendelian inheritance patterns, such as uniparental disomy (UPD) and chimerism. A genome-wide depth and variant allele frequency plot generated by TrioMix facilitates tracing the origin of Mendelian inheritance deviations. We showed that TrioMix could accurately deconvolute genomes in both simulated and real data sets.