Gene Variants Involved in Nonsense-Mediated mRNA Decay Suggest a Role in Autism Spectrum Disorder

SOX7: Autism Associated Gene Identified by Analysis of Multi-Omics Data

Genome-wide association studies and next generation sequencing data analyses based on DNA information have identified thousands of mutations associated with autism spectrum disorder (ASD). However, more than 99% of identified mutations are non-coding. Thus, it is unclear which of these mutations might be functional and thus potentially causal variants. Transcriptomic profiling using total RNA-sequencing has been one of the most utilized approaches to link protein levels to genetic information at the molecular level. The transcriptome captures molecular genomic complexity that the DNA sequence solely does not. Some mutations alter a gene's DNA sequence but do not necessarily change expression and/or protein function. To date, few common variants reliably associated with the diagnosis status of ASD despite consistently high estimates of heritability. In addition, reliable biomarkers used to diagnose ASD or molecular mechanisms to define the severity of ASD do not exist. Therefore, it is necessary to integrate DNA and RNA testing together to identify true causal genes and propose useful biomarkers for ASD. We performed gene-based association studies with adaptive test using genome-wide association studies (GWAS) summary statistics with two large GWAS datasets (ASD 2019 data: 18,382 ASD cases and 27,969 controls [discovery data]; ASD 2017 data: 6,197 ASD cases and 7,377 controls [replication data]) which were obtained from the Psychiatric Genomics Consortium (PGC). In addition, we investigated differential expression between ASD cases and controls for genes identified in gene-based GWAS with two RNA-seq datasets (GSE211154: 20 cases and 19 controls; GSE30573: 3 cases and 3 controls). We identified 5 genes significantly associated with ASD in ASD 2019 data (KIZ-AS1, p=8.67×10-10; KIZ, p=1.16×10-9; XRN2, p=7.73×10-9; SOX7, p=2.22×10-7; LOC101929229 also known as PINX1-DT, p=2.14×10-6). Among these 5 genes, gene SOX7 (p=0.00087) and LOC101929229 (p=0.009) were replicated in ASD 2017 data. KIZ-AS1 (p=0.059) and KIZ (p=0.06) were close to the boundary of replication in ASD 2017 data. Genes SOX7 (p=0.036 in all samples; p=0.044 in white samples) indicated significant expression differences between cases and controls in the GSE211154 RNA-seq data. Furthermore, gene SOX7 was upregulated in cases than in controls in the GSE30573 RNA-seq data (p=0.0017; Benjamini-Hochberg adjusted p=0.0085). SOX7 encodes a member of the SOX (SRY-related HMG-box) family of transcription factors pivotally contributing to determining of the cell fate and identity in many lineages. The encoded protein may act as a transcriptional regulator after forming a protein complex with other proteins leading to autism. Gene SOX7 in the transcription factor family could be associated with ASD. This finding may provide new diagnostic and therapeutic strategies for ASD.

Direct and indirect effects of spliceosome disruption compromise gene regulation by Nonsense-Mediated mRNA Decay

Pre-mRNA splicing, carried out in the nucleus by a large ribonucleoprotein machine known as the spliceosome, is functionally and physically coupled to the mRNA surveillance pathway in the cytoplasm called nonsense mediated mRNA decay (NMD). The NMD pathway monitors for premature translation termination signals, which can result from alternative splicing, by relying on the exon junction complex (EJC) deposited on exon-exon junctions by the spliceosome. Recently, multiple genetic screens in human cell lines have identified numerous spliceosome components as putative NMD factors. Using publicly available RNA-seq datasets from K562 and HepG2 cells depleted of 18 different spliceosome components, we find that natural NMD targeted mRNA isoforms are upregulated when members of the catalytic spliceosome are reduced. While some of this increase could be due to widespread pleiotropic effects of spliceosome dysfunction (e.g., reduced expression of NMD factors due to mis-splicing of their mRNAs), we identify that AQR, SF3B1, SF3B4 and CDC40 may have a more direct role in NMD. We also test the hypothesis that increased production of novel NMD substrates may overwhelm the pathway to find a direct correlation between the amount of novel NMD substrates detected and the degree of NMD inhibition observed. Finally, similar transcriptome alterations and NMD substrate upregulation are also observed in cells treated with spliceosome inhibitors and in cells derived from retinitis pigmentosa patients with mutations in PRPF8 and PRPF31. Overall, our results show that regardless of the cause, spliceosome disruption upregulates a broad set of NMD targets, which could contribute to cellular dysfunction in spliceosomopathies.

A critical review of the impact of candidate copy number variants on autism spectrum disorder

Autism spectrum disorder (ASD) is a complex neurodevelopmental disorder (NDD) influenced by genetic, epigenetic, and environmental factors. Recent advancements in genomic analysis have shed light on numerous genes associated with ASD, highlighting the significant role of both common and rare genetic mutations, as well as copy number variations (CNVs), single nucleotide polymorphisms (SNPs) and unique de novo variants. These genetic variations disrupt neurodevelopmental pathways, contributing to the disorder's complexity. Notably, CNVs are present in 10 %-20 % of individuals with autism, with 3 %-7 % detectable through cytogenetic methods. While the role of submicroscopic CNVs in ASD has been recently studied, their association with genomic loci and genes has not been thoroughly explored. In this review, we focus on 47 CNV regions linked to ASD, encompassing 1632 genes, including protein-coding genes and long non-coding RNAs (lncRNAs), of which 659 show significant brain expression. Using a list of ASD-associated genes from SFARI, we detect 17 regions harboring at least one known ASD-related protein-coding gene. Of the remaining 30 regions, we identify 24 regions containing at least one protein-coding gene with brain-enriched expression and a nervous system phenotype in mouse mutants, and one lncRNA with both brain-enriched expression and upregulation in iPSC to neuron differentiation. This review not only expands our understanding of the genetic diversity associated with ASD but also underscores the potential of lncRNAs in contributing to its etiology. Additionally, the discovered CNVs will be a valuable resource for future diagnostic, therapeutic, and research endeavors aimed at prioritizing genetic variations in ASD.

Spatial expression of the nonsense-mediated mRNA decay factors UPF3A and UPF3B among mouse tissues

无义介导的信使RNA（mRNA）降解途径（nonsense-mediated mRNA decay，简称为NMD)是真核生物细胞内一种重要的基因转录后表达调控机制，它积极参与一系列细胞生理和生化过程，控制细胞命运和生命体的组织稳态。NMD的缺陷会导致人类疾病，如神经发育障碍、肿瘤发生和自身免疫疾病等。UPF3 (Up-frameshift protein 3)是一个核心的NMD因子，它最早在酵母中被发现。UPF3A和UPF3B是UPF3在生物进化到脊椎动物阶段出现的两个旁系同源物，在NMD中具有激活或抑制的作用。以往研究发现，UPF3B蛋白几乎在所有哺乳动物器官中均有表达，而UPF3A蛋白在除睾丸外的大多数哺乳动物组织中难以被检测到。解释这一现象的假说为：在NMD途径中，UPF3B具有比UPF3A更高的竞争性结合UPF2的能力，UPF3B和UPF2的结合促使UPF3A成为游离状态，而游离的UPF3A蛋白不稳定且易被降解。此假说提示UPF3A和UPF3B在NMD中存在拮抗作用。在本研究中，我们重新定量评估了UPF3A和UPF3B在野生型成年雄性和雌性小鼠的9个主要组织和生殖器官中的mRNA和蛋白表达，结果证实UPF3A在雄性生殖细胞中表达量最高。令人惊讶的是，我们发现在包括大脑和胸腺在内的大多数组织中，UPF3A与UPF3B的蛋白水平相当，而在小鼠脾、肺组织中，UPF3A表达高于UPF3B。公共基因表达数据进一步支持了上述发现。因此，我们的研究表明了UPF3A是小鼠组织中普遍表达的NMD因子。同时，该研究结果推测：在生理条件下，UPF3A和UPF3B蛋白之间不存在竞争抑制，且UPF3A在多种哺乳动物组织的稳态中发挥重要作用。

SOX7: Novel Autistic Gene Identified by Analysis of Multi-Omics Data

Background

Despite thousands of variants identified by genome-wide association studies (GWAS) to be associated with autism spectrum disorder (ASD), it is unclear which mutations are causal because most are noncoding. Consequently, reliable diagnostic biomarkers are lacking. RNA-seq analysis captures biomolecular complexity that GWAS cannot by considering transcriptomic patterns. Therefore, integrating DNA and RNA testing may reveal causal genes and useful biomarkers for ASD.

Methods

We performed gene-based association studies using an adaptive test method with GWAS summary statistics from two large Psychiatric Genomics Consortium (PGC) datasets (ASD2019: 18,382 cases and 27,969 controls; ASD2017: 6,197 cases and 7,377 controls). We also investigated differential expression for genes identified with the adaptive test using an RNA-seq dataset (GSE30573: 3 cases and 3 controls) and DESeq2.

Results

We identified 5 genes significantly associated with ASD in ASD2019 (KIZ-AS1, p = 8.67×10- 10; KIZ, p = 1.16×10- 9; XRN2, p = 7.73×10- 9; SOX7, p = 2.22×10- 7; LOC101929229 (also known as PINX1-DT), p = 2.14×10- 6). Two of the five genes were replicated in ASD2017: SOX7 (p = 0.00087) and LOC101929229 (p = 0.009), and KIZ was close to the replication boundary of replication (p = 0.06). We identified significant expression differences for SOX7 (p = 0.0017, adjusted p = 0.0085), LOC101929229 (p = 5.83×10- 7, adjusted p = 1.18×10- 5), and KIZ (p = 0.00099, adjusted p = 0.0055). SOX7 encodes a transcription factor that regulates developmental pathways, alterations in which may contribute to ASD.

Limitations

The limitation of the gene-based analysis is the reliance on a reference population for estimating linkage disequilibrium between variants. The similarity of this reference population to the population of study is crucial to the accuracy of many gene-based analyses, including those performed in this study. As a result, the extent of our findings is limited to European populations, as this was our reference of choice. Future work includes a tighter integration of DNA and RNA information as well as extensions to non-European populations that have been under-researched.

Conclusions

These findings suggest that SOX7 and its related SOX family genes encode transcription factors that are critical to the downregulation of the canonical Wnt/β-catenin signaling pathway, an important developmental signaling pathway, providing credence to the biologic plausibility of the association between gene SOX7 and autism spectrum disorder.

Genetic ablation of metabotropic glutamate receptor 5 in rats results in an autism-like behavioral phenotype

Autism spectrum disorder (ASD) is a neurodevelopmental disorder characterized by deficits in communication, and social skills, as well as repetitive and/or restrictive interests and behaviors. The severity of ASD varies from mild to severe, drastically interfering with the quality of life of affected individuals. The current occurrence of ASD in the United States is about 1 in 44 children. The precise pathophysiology of ASD is still unknown, but it is believed that ASD is heterogeneous and can arise due to genetic etiology. Although various genes have been implicated in predisposition to ASD, metabotropic glutamate receptor 5 (mGluR5) is one of the most common downstream targets, which may be involved in autism. mGluR5 signaling has been shown to play a crucial role in neurodevelopment and neural transmission making it a very attractive target for understanding the pathogenesis of ASD. In the present study, we determined the effect of genetic ablation of mGluR5 (Grm5) on an ASD-like phenotype using a rat model to better understand the role of mGluR5 signaling in behavior patterns and clinical manifestations of ASD. We observed that mGluR5 Ko rats exhibited exaggerated self-grooming and increased marble burying, as well as deficits in social novelty. Our results suggest that mGluR5 Ko rats demonstrate an ASD-like phenotype, specifically impaired social interaction as well as repetitive and anxiety-like behavior, which are correlates of behavior symptoms observed in individuals with ASD. The mGluR5 Ko rat model characterized in this study may be explored to understand the molecular mechanisms underlying ASD and for developing effective therapeutic modalities.