Cell cycle networks link gene expression dysregulation, mutation, and brain maldevelopment in autistic toddlers

Embryonic origin of two ASD subtypes of social symptom severity: the larger the brain cortical organoid size, the more severe the social symptoms

BACKGROUND

Social affective and communication symptoms are central to autism spectrum disorder (ASD), yet their severity differs across toddlers: Some toddlers with ASD display improving abilities across early ages and develop good social and language skills, while others with "profound" autism have persistently low social, language and cognitive skills and require lifelong care. The biological origins of these opposite ASD social severity subtypes and developmental trajectories are not known.

METHODS

Because ASD involves early brain overgrowth and excess neurons, we measured size and growth in 4910 embryonic-stage brain cortical organoids (BCOs) from a total of 10 toddlers with ASD and 6 controls (averaging 196 individual BCOs measured/subject). In a 2021 batch, we measured BCOs from 10 ASD and 5 controls. In a 2022 batch, we  tested replicability of BCO size and growth effects by generating and measuring an independent batch of BCOs from 6 ASD and 4 control subjects. BCO size was analyzed within the context of our large, one-of-a-kind social symptom, social attention, social brain and social and language psychometric normative datasets ranging from N = 266 to N = 1902 toddlers. BCO growth rates were examined by measuring size changes between 1- and 2-months of organoid development. Neurogenesis markers at 2-months were examined at the cellular level. At the molecular level, we measured activity and expression of Ndel1; Ndel1 is a prime target for cell cycle-activated kinases; known to regulate cell cycle, proliferation, neurogenesis, and growth; and known to be involved in neuropsychiatric conditions.

RESULTS

At the BCO level, analyses showed BCO size was significantly enlarged by 39% and 41% in ASD in the 2021 and 2022 batches. The larger the embryonic BCO size, the more severe the ASD social symptoms. Correlations between BCO size and social symptoms were r = 0.719 in the 2021 batch and r = 0. 873 in the replication 2022 batch. ASD BCOs grew at an accelerated rate nearly 3 times faster than controls. At the cell level, the two largest ASD BCOs had accelerated neurogenesis. At the molecular level, Ndel1 activity was highly correlated with the growth rate and size of BCOs. Two BCO subtypes were found in ASD toddlers: Those in one subtype had very enlarged BCO size with accelerated rate of growth and neurogenesis; a profound autism clinical phenotype displaying severe social symptoms, reduced social attention, reduced cognitive, very low language and social IQ; and substantially altered growth in specific cortical social, language and sensory regions. Those in a second subtype had milder BCO enlargement and milder social, attention, cognitive, language and cortical differences.

LIMITATIONS

Larger samples of ASD toddler-derived BCO and clinical phenotypes may reveal additional ASD embryonic subtypes.

CONCLUSIONS

By embryogenesis, the biological bases of two subtypes of ASD social and brain development-profound autism and mild autism-are already present and measurable and involve dysregulated cell proliferation and accelerated neurogenesis and growth. The larger the embryonic BCO size in ASD, the more severe the toddler's social symptoms and the more reduced the social attention, language ability, and IQ, and the more atypical the growth of social and language brain regions.

A 3D approach to understanding heterogeneity in early developing autisms

Phenotypic heterogeneity in early language, intellectual, motor, and adaptive functioning (LIMA) features are amongst the most striking features that distinguish different types of autistic individuals. Yet the current diagnostic criteria uses a single label of autism and implicitly emphasizes what individuals have in common as core social-communicative and restricted repetitive behavior difficulties. Subtype labels based on the non-core LIMA features may help to more meaningfully distinguish types of autisms with differing developmental paths and differential underlying biology. Using relatively large (n=615) publicly available data from early developing (24-68 months) standardized clinical tests tapping LIMA features, we show that stability-based relative cluster validation analysis can identify two robust and replicable clusters in the autism population with high levels of generalization accuracy (98%). These clusters can be described as Type I versus Type II autisms differentiated by relatively high versus low scores on LIMA features. These two types of autisms are also distinguished by different developmental trajectories over the first decade of life. Finally, these two types of autisms reveal striking differences in functional and structural neuroimaging phenotypes and their relationships with gene expression. This work emphasizes the potential importance of stratifying autism by a Type I versus Type II distinction focused on LIMA features and which may be of high prognostic and biological significance.

BRCA1 and ELK-1 regulate neural progenitor cell fate in the optic tectum in response to visual experience in Xenopus laevis tadpoles

In developing Xenopus tadpoles, the optic tectum begins to receive patterned visual input while visuomotor circuits are still undergoing neurogenesis and circuit assembly. This visual input regulates neural progenitor cell fate decisions such that maintaining tadpoles in the dark increases proliferation, expanding the progenitor pool, while visual stimulation promotes neuronal differentiation. To identify regulators of activity-dependent neural progenitor cell fate, we profiled the transcriptomes of proliferating neural progenitor cells and newly differentiated neurons using RNA-Seq. We used advanced bioinformatic analysis of 1,130 differentially expressed transcripts to identify six differentially regulated transcriptional regulators, including Breast Cancer 1 (BRCA1) and the ETS-family transcription factor, ELK-1, which are predicted to regulate the majority of the other differentially expressed transcripts. BRCA1 is known for its role in cancers, but relatively little is known about its potential role in regulating neural progenitor cell fate. ELK-1 is a multifunctional transcription factor which regulates immediate early gene expression. We investigated the potential functions of BRCA1 and ELK-1 in activity-regulated neurogenesis in the tadpole visual system using in vivo time-lapse imaging to monitor the fate of GFP-expressing SOX2+ neural progenitor cells in the optic tectum. Our longitudinal in vivo imaging analysis showed that knockdown of either BRCA1 or ELK-1 altered the fates of neural progenitor cells and furthermore that the effects of visual experience on neurogenesis depend on BRCA1 and ELK-1 expression. These studies provide insight into the potential mechanisms by which neural activity affects neural progenitor cell fate.

Genomic insights and advanced machine learning: characterizing autism spectrum disorder biomarkers and genetic interactions

Autism Spectrum Disorder (ASD) is a complex neurodevelopmental condition characterized by altered brain connectivity and function. In this study, we employed advanced bioinformatics and explainable AI to analyze gene expression associated with ASD, using data from five GEO datasets. Among 351 neurotypical controls and 358 individuals with autism, we identified 3,339 Differentially Expressed Genes (DEGs) with an adjusted p-value (≤ 0.05). A subsequent meta-analysis pinpointed 342 DEGs (adjusted p-value ≤ 0.001), including 19 upregulated and 10 down-regulated genes across all datasets. Shared genes, pathogenic single nucleotide polymorphisms (SNPs), chromosomal positions, and their impact on biological pathways were examined. We identified potential biomarkers (HOXB3, NR2F2, MAPK8IP3, PIGT, SEMA4D, and SSH1) through text mining, meriting further investigation. Additionally, ‎we shed light on the roles of RPS4Y1 and KDM5D genes in neurogenesis and neurodevelopment. Our analysis detected 1,286 SNPs linked to ASD-related conditions, of which 14 high-risk SNPs were located on chromosomes 10 and X. We highlighted potential missense SNPs associated with FGFR inhibitors, suggesting that it may serve as a promising biomarker for responsiveness to targeted therapies. Our explainable AI model identified the MID2 gene as a potential ASD biomarker. This research unveils vital genes and potential biomarkers, providing a foundation for novel gene discovery in complex diseases.

The molecular genetic landscape of human brain size variation

Human brain size changes dynamically through early development, peaks in adolescence, and varies up to 2-fold among adults. However, the molecular genetic underpinnings of interindividual variation in brain size remain unknown. Here, we leveraged postmortem brain RNA sequencing and measurements of brain weight (BW) in 2,531 individuals across three independent datasets to identify 928 genome-wide significant associations with BW. Genes associated with higher or lower BW showed distinct neurodevelopmental trajectories and spatial patterns that mapped onto functional and cellular axes of brain organization. Expression of BW genes was predictive of interspecies differences in brain size, and bioinformatic annotation revealed enrichment for neurogenesis and cell-cell communication. Genome-wide, transcriptome-wide, and phenome-wide association analyses linked BW gene sets to neuroimaging measurements of brain size and brain-related clinical traits. Cumulatively, these results represent a major step toward delineating the molecular pathways underlying human brain size variation in health and disease.

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.

Artificial intelligence and bioinformatics analyze markers of children's transcriptional genome to predict autism spectrum disorder

Introduction

Autism spectrum disorder (ASD), characterized by difficulties in social interaction and communication as well as restricted interests and repetitive behaviors, is extremely challenging to diagnose in toddlers. Early diagnosis and intervention are crucial however.

Methods

In this study, we developed a machine learning classification model based on mRNA expression data from the peripheral blood of 128 toddlers with ASD and 126 controls. Differentially expressed genes (DEGs) between ASD and controls were identified.

Results

We identified genes such as UBE4B, SPATA2 and RBM3 as DEGs, mainly involved in immune-related pathways. 21 genes were screened as key biomarkers using LASSO regression, yielding an accuracy of 86%. A neural network model based on these 21 genes achieved an AUC of 0.88.

Discussion

Our findings suggest that the identified neurotransmitters and 21 immune-related biomarkers may facilitate the early diagnosis of ASD. The mRNA expression profile sheds light on the biological underpinnings of ASD in toddlers and potential biomarkers for early identification. Nevertheless, larger samples are needed to validate these biomarkers.

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

Background

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.

Objectives

It is necessary to integrate DNA and RNA testing together to identify true causal genes and propose useful biomarkers for ASD.

Methods

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 for genes identified in gene-based GWAS with a RNA-seq dataset (GSE30573: 3 cases and 3 controls) using the DESeq2 package.

Results

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; PINX1-DT, p=2.14×10-6). Among these 5 genes, gene SOX7 (p=0.00087), LOC101929229 (p=0.009), and KIZ-AS1 (p=0.059) were replicated in ASD 2017 data. KIZ (p=0.06) was close to the boundary of replication in ASD 2017 data. Genes SOX7 (p=0.0017, adjusted p=0.0085), LOC101929229 (also known as PINX1-DT, p=5.83×10-7, adjusted p=1.18×10-5), and KIZ (p=0.00099, adjusted p=0.0055) indicated significant expression differences between cases and controls in the RNA-seq data. 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.

Conclusion

Gene SOX7 in the transcription factor family could be associated with ASD. This finding may provide new diagnostic and therapeutic strategies for ASD.

Novel correlative analysis identifies multiple genomic variations impacting ASD with macrocephaly

Autism spectrum disorders (ASD) display both phenotypic and genetic heterogeneity, impeding the understanding of ASD and development of effective means of diagnosis and potential treatments. Genes affected by genomic variations for ASD converge in dozens of gene ontologies (GOs), but the relationship between the variations at the GO level have not been well elucidated. In the current study, multiple types of genomic variations were mapped to GOs and correlations among GOs were measured in ASD and control samples. Several ASD-unique GO correlations were found, suggesting the importance of co-occurrence of genomic variations in genes from different functional categories in ASD etiology. Combined with experimental data, several variations related to WNT signaling, neuron development, synapse morphology/function and organ morphogenesis were found to be important for ASD with macrocephaly, and novel co-occurrence patterns of them in ASD patients were found. Furthermore, we applied this gene ontology correlation analysis method to find genomic variations that contribute to ASD etiology in combination with changes in gene expression and transcription factor binding, providing novel insights into ASD with macrocephaly and a new methodology for the analysis of genomic variation.

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.

Top2a promotes the development of social behavior via PRC2 and H3K27me3

Little is understood about the embryonic development of sociality. We screened 1120 known drugs and found that embryonic inhibition of topoisomerase IIα (Top2a) resulted in lasting social deficits in zebrafish. In mice, prenatal Top2 inhibition caused defects in social interaction and communication, which are behaviors that relate to core symptoms of autism. Mutation of Top2a in zebrafish caused down-regulation of a set of genes highly enriched for genes associated with autism in humans. Both the Top2a-regulated and autism-associated gene sets have binding sites for polycomb repressive complex 2 (PRC2), a regulatory complex responsible for H3K27 trimethylation (H3K27me3). Moreover, both gene sets are highly enriched for H3K27me3. Inhibition of the PRC2 component Ezh2 rescued social deficits caused by Top2 inhibition. Therefore, Top2a is a key component of an evolutionarily conserved pathway that promotes the development of social behavior through PRC2 and H3K27me3.

Pre-treatment clinical and gene expression patterns predict developmental change in early intervention in autism

Early detection and intervention are believed to be key to facilitating better outcomes in children with autism, yet the impact of age at treatment start on the outcome is poorly understood. While clinical traits such as language ability have been shown to predict treatment outcome, whether or not and how information at the genomic level can predict treatment outcome is unknown. Leveraging a cohort of toddlers with autism who all received the same standardized intervention at a very young age and provided a blood sample, here we find that very early treatment engagement (i.e., <24 months) leads to greater gains while controlling for time in treatment. Pre-treatment clinical behavioral measures predict 21% of the variance in the rate of skill growth during early intervention. Pre-treatment blood leukocyte gene expression patterns also predict the rate of skill growth, accounting for 13% of the variance in treatment slopes. Results indicated that 295 genes can be prioritized as driving this effect. These treatment-relevant genes highly interact at the protein level, are enriched for differentially histone acetylated genes in autism postmortem cortical tissue, and are normatively highly expressed in a variety of subcortical and cortical areas important for social communication and language development. This work suggests that pre-treatment biological and clinical behavioral characteristics are important for predicting developmental change in the context of early intervention and that individualized pre-treatment biology related to histone acetylation may be key.

Postsynaptic autism spectrum disorder genes and synaptic dysfunction

This review provides an overview of the synaptic dysfunction of neuronal circuits and the ensuing behavioral alterations caused by mutations in autism spectrum disorder (ASD)-linked genes directly or indirectly affecting the postsynaptic neuronal compartment. There are plenty of ASD risk genes, that may be broadly grouped into those involved in gene expression regulation (epigenetic regulation and transcription) and genes regulating synaptic activity (neural communication and neurotransmission). Notably, the effects mediated by ASD-associated genes can vary extensively depending on the developmental time and/or subcellular site of expression. Therefore, in order to gain a better understanding of the mechanisms of disruptions in postsynaptic function, an effort to better model ASD in experimental animals is required to improve standardization and increase reproducibility within and among studies. Such an effort holds promise to provide deeper insight into the development of these disorders and to improve the translational value of preclinical studies.

Gene Regulation Analysis Reveals Perturbations of Autism Spectrum Disorder during Neural System Development

Autism spectrum disorder (ASD) is a neurodevelopmental disorder that impedes patients' cognition, social, speech and communication skills. ASD is highly heterogeneous with a variety of etiologies and clinical manifestations. The prevalence rate of ASD increased steadily in recent years. Presently, molecular mechanisms underlying ASD occurrence and development remain to be elucidated. Here, we integrated multi-layer genomics data to investigate the transcriptome and pathway dysregulations in ASD development. The RNA sequencing (RNA-seq) expression profiles of induced pluripotent stem cells (iPSCs), neural progenitor cells (NPCs) and neuron cells from ASD and normal samples were compared in our study. We found that substantially more genes were differentially expressed in the NPCs than the iPSCs. Consistently, gene set variation analysis revealed that the activity of the known ASD pathways in NPCs and neural cells were significantly different from the iPSCs, suggesting that ASD occurred at the early stage of neural system development. We further constructed comprehensive brain- and neural-specific regulatory networks by incorporating transcription factor (TF) and gene interactions with long 5 non-coding RNA(lncRNA) and protein interactions. We then overlaid the transcriptomes of different cell types on the regulatory networks to infer the regulatory cascades. The variations of the regulatory cascades between ASD and normal samples uncovered a set of novel disease-associated genes and gene interactions, particularly highlighting the functional roles of ELF3 and the interaction between STAT1 and lncRNA ELF3-AS 1 in the disease development. These new findings extend our understanding of ASD and offer putative new therapeutic targets for further studies.

Atypical genomic cortical patterning in autism with poor early language outcome

Cortical regionalization develops via genomic patterning along anterior-posterior (A-P) and dorsal-ventral (D-V) gradients. Here, we find that normative A-P and D-V genomic patterning of cortical surface area (SA) and thickness (CT), present in typically developing and autistic toddlers with good early language outcome, is absent in autistic toddlers with poor early language outcome. Autistic toddlers with poor early language outcome are instead specifically characterized by a secondary and independent genomic patterning effect on CT. Genes involved in these effects can be traced back to midgestational A-P and D-V gene expression gradients and different prenatal cell types (e.g., progenitor cells and excitatory neurons), are functionally important for vocal learning and human-specific evolution, and are prominent in prenatal coexpression networks enriched for high-penetrance autism risk genes. Autism with poor early language outcome may be explained by atypical genomic cortical patterning starting in prenatal development, which may detrimentally affect later regional functional specialization and circuit formation.

The Role of Hippo Signaling Pathway in the Development of the Nervous System

Hippo signaling pathway is a highly conserved and crucial signaling pathway that controls the size of tissues and organs by regulating the proliferation, differentiation, and apoptosis of cells. The nervous system is a complicated system that participates in information collection, integration, and procession. The balance of various aspects of the nervous system is vital for the normal regulation of physiological conditions of the body, like the population and distribution of nerve cells, nerve connections, and so on. Defects in these aspects may lead to cognitive, behavioral, and neurological dysfunction, resulting in various nervous system diseases. Recently, accumulating evidence proposes that Hippo pathway maintains numerous biological functions in the nervous system development, including modulating the proliferation and differentiation of nerve cells and promoting the development of synapse, corpus callosum, and cortex. In this review, we will summarize recent findings of Hippo pathway in the nervous system to improve our understanding on its function and to provide potential therapeutic strategies of nervous system diseases in the future.

Gene Set Enrichment Analysis of Interaction Networks Weighted by Node Centrality

Gene set enrichment analysis (GSEA) is a powerful tool to associate a disease phenotype to a group of genes/proteins. GSEA attributes a specific weight to each gene/protein in the input list that depends on a metric of choice, which is usually represented by quantitative expression data. However, expression data are not always available. Here, GSEA based on betweenness centrality of a protein–protein interaction (PPI) network is described and applied to two cases, where an expression metric is missing. First, personalized PPI networks were generated from genes displaying alterations (assessed by array comparative genomic hybridization and whole exome sequencing) in four probands bearing a 16p13.11 microdeletion in common and several other point variants. Patients showed disease phenotypes linked to neurodevelopment. All networks were assembled around a cluster of first interactors of altered genes with high betweenness centrality. All four clusters included genes known to be involved in neurodevelopmental disorders with different centrality. Moreover, the GSEA results pointed out to the evidence of “cell cycle” among enriched pathways. Second, a large interaction network obtained by merging proteomics studies on three neurodegenerative disorders was analyzed from the topological point of view. We observed that most central proteins are often linked to Parkinson’s disease. The selection of these proteins improved the specificity of GSEA, with “Metabolism of amino acids and derivatives” and “Cellular response to stress or external stimuli” as top-ranked enriched pathways. In conclusion, betweenness centrality revealed to be a suitable metric for GSEA. Thus, centrality-based GSEA represents an opportunity for precision medicine and network medicine.

Blood biomarker discovery for autism spectrum disorder: A proteomic analysis

Autism spectrum disorder (ASD) is a neurodevelopmental disorder characterized by deficits in social communication and social interaction and restricted, repetitive patterns of behavior, interests, or activities. Given the lack of specific pharmacological therapy for ASD and the clinical heterogeneity of the disorder, current biomarker research efforts are geared mainly toward identifying markers for determining ASD risk or for assisting with a diagnosis. A wide range of putative biological markers for ASD is currently being investigated. Proteomic analyses indicate that the levels of many proteins in plasma/serum are altered in ASD, suggesting that a panel of proteins may provide a blood biomarker for ASD. Serum samples from 76 boys with ASD and 78 typically developing (TD) boys, 18 months-8 years of age, were analyzed to identify possible early biological markers for ASD. Proteomic analysis of serum was performed using SomaLogic’s SOMAScan^TM assay 1.3K platform. A total of 1,125 proteins were analyzed. There were 86 downregulated proteins and 52 upregulated proteins in ASD (FDR < 0.05). Combining three different algorithms, we found a panel of 9 proteins that identified ASD with an area under the curve (AUC) = 0.8599±0.0640, with specificity and sensitivity of 0.8217±0.1178 and 0.835±0.1176, respectively. All 9 proteins were significantly different in ASD compared with TD boys, and were significantly correlated with ASD severity as measured by ADOS total scores. Using machine learning methods, a panel of serum proteins was identified that may be useful as a blood biomarker for ASD in boys. Further verification of the protein biomarker panel with independent test sets is warranted.

Coupling of autism genes to tissue-wide expression and dysfunction of synapse, calcium signalling and transcriptional regulation

Autism Spectrum Disorder (ASD) is a heterogeneous disorder that is often accompanied with many co-morbidities. Recent genetic studies have identified various pathways from hundreds of candidate risk genes with varying levels of association to ASD. However, it is unknown which pathways are specific to the core symptoms or which are shared by the co-morbidities. We hypothesised that critical ASD candidates should appear widely across different scoring systems, and that comorbidity pathways should be constituted by genes expressed in the relevant tissues. We analysed the Simons Foundation for Autism Research Initiative (SFARI) database and four independently published scoring systems and identified 292 overlapping genes. We examined their mRNA expression using the Genotype-Tissue Expression (GTEx) database and validated protein expression levels using the human protein atlas (HPA) dataset. This led to clustering of the overlapping ASD genes into 2 groups; one with 91 genes primarily expressed in the central nervous system (CNS geneset) and another with 201 genes expressed in both CNS and peripheral tissues (CNS+PT geneset). Bioinformatic analyses showed a high enrichment of CNS development and synaptic transmission in the CNS geneset, and an enrichment of synapse, chromatin remodelling, gene regulation and endocrine signalling in the CNS+PT geneset. Calcium signalling and the glutamatergic synapse were found to be highly interconnected among pathways in the combined geneset. Our analyses demonstrate that 2/3 of ASD genes are expressed beyond the brain, which may impact peripheral function and involve in ASD co-morbidities, and relevant pathways may be explored for the treatment of ASD co-morbidities.

Identification and Functional Analysis of Long Non-coding RNAs in Autism Spectrum Disorders

Genetic and environmental factors, alone or in combination, contribute to the pathogenesis of autism spectrum disorder (ASD). Although many protein-coding genes have now been identified as disease risk genes for ASD, a detailed illustration of long non-coding RNAs (lncRNAs) associated with ASD remains elusive. In this study, we first identified ASD-related lncRNAs based on genomic variant data of individuals with ASD from a twin study. In total, 532 ASD-related lncRNAs were identified, and 86.7% of these ASD-related lncRNAs were further validated by an independent copy number variant (CNV) dataset. Then, the functions and associated biological pathways of ASD-related lncRNAs were explored by enrichment analysis of their three different types of functional neighbor genes (i.e., genomic neighbors, competing endogenous RNA (ceRNA) neighbors, and gene co-expression neighbors in the cortex). The results have shown that most of the functional neighbor genes of ASD-related lncRNAs were enriched in nervous system development, inflammatory responses, and transcriptional regulation. Moreover, we explored the differential functions of ASD-related lncRNAs in distinct brain regions by using gene co-expression network analysis based on tissue-specific gene expression profiles. As a set, ASD-related lncRNAs were mainly associated with nervous system development and dopaminergic synapse in the cortex, but associated with transcriptional regulation in the cerebellum. In addition, a functional network analysis was conducted for the highly reliable functional neighbor genes of ASD-related lncRNAs. We found that all the highly reliable functional neighbor genes were connected in a single functional network, which provided additional clues for the action mechanisms of ASD-related lncRNAs. Finally, we predicted several potential drugs based on the enrichment of drug-induced pathway sets in the ASD-altered biological pathway list. Among these drugs, several (e.g., amoxapine, piperine, and diflunisal) were partly supported by the previous reports. In conclusion, ASD-related lncRNAs participated in the pathogenesis of ASD through various known biological pathways, which may be differential in distinct brain regions. Detailed investigation into ASD-related lncRNAs can provide clues for developing potential ASD diagnosis biomarkers and therapy.

Cerebral organoids as tools to identify the developmental roots of autism

Some autism spectrum disorders (ASD) likely arise as a result of abnormalities during early embryonic development of the brain. Studying human embryonic brain development directly is challenging, mainly due to ethical and practical constraints. However, the recent development of cerebral organoids provides a powerful tool for studying both normal human embryonic brain development and, potentially, the origins of neurodevelopmental disorders including ASD. Substantial evidence now indicates that cerebral organoids can mimic normal embryonic brain development and neural cells found in organoids closely resemble their in vivo counterparts. However, with prolonged culture, significant differences begin to arise. We suggest that cerebral organoids, in their current form, are most suitable to model earlier neurodevelopmental events and processes such as neurogenesis and cortical lamination. Processes implicated in ASDs which occur at later stages of development, such as synaptogenesis and neural circuit formation, may also be modeled using organoids. The accuracy of such models will benefit from continuous improvements to protocols for organoid differentiation.

Integrated Systems Analysis Explores Dysfunctional Molecular Modules and Regulatory Factors in Children with Autism Spectrum Disorder

Autism spectrum disorder (ASD) is a genetic neurodevelopmental disorder involving multiple genes that occurs in early childhood, and a number of risk genes have been reported in previous studies. However, the molecular mechanism of the polygenic regulation leading to pathological changes in ASD remains unclear. First, we identified 8 dysregulated gene coexpression modules by analyzing blood transcriptome data from 96 children with ASD and 42 controls. These modules are rich in ASD risk genes and function related to metabolism, immunity, neurodevelopment, and signaling. The regulatory factors of each module including microRNA (miRNA) and transcription factors (TFs) were subsequently predicted based on transcriptional and posttranscriptional regulation. We identified a set of miRNAs that regulate metabolic and immune modules, as well as transcription factors that cause dysregulation of the modules, and we constructed a coregulatory network between the regulatory factors and modules. Our work reveals dysfunctional modules in children with ASD, elucidates the role of miRNA and transcription factor dysregulation in the pathophysiology of ASD, and helps us to further understand the underlying molecular mechanism of ASD.

Maternal Immune Activation Alters Fetal Brain Development and Enhances Proliferation of Neural Precursor Cells in Rats

Maternal immune activation (MIA) caused by exposure to pathogens or inflammation during critical periods of neurodevelopment is a major risk factor for behavioral deficits and psychiatric illness in offspring. A spectrum of behavioral abnormalities can be recapitulated in rodents by inducing MIA using the viral mimetic, PolyI:C. Many studies have focused on long-term changes in brain structure and behavioral outcomes in offspring following maternal PolyI:C exposure, but acute changes in prenatal development are not well-characterized. Using RNA-Sequencing, we profiled acute transcriptomic changes in rat conceptuses (decidua along with nascent embryo and placenta) after maternal PolyI:C exposure during early gestation, which enabled us to capture gene expression changes provoked by MIA inclusive to the embryonic milieu. We identified a robust increase in expression of genes related to antiviral inflammation following maternal PolyI:C exposure, and a corresponding decrease in transcripts associated with nervous system development. At mid-gestation, regions of the developing cortex were thicker in fetuses prenatally challenged with PolyI:C, with females displaying a thicker ventricular zone and males a thicker cortical mantle. Along these lines, neural precursor cells (NPCs) isolated from fetal brains prenatally challenged with PolyI:C exhibited a higher rate of self-renewal. Expression of Notch1 and the Notch ligand, delta-like ligand 1, which are both highly implicated in maintenance of NPCs and nervous system development, was increased following PolyI:C exposure. These results suggest that MIA elicits rapid gene expression changes within the conceptus, including repression of neurodevelopmental pathways, resulting in profound alterations in fetal brain development.

Prenatal Origins of ASD: The When, What, and How of ASD Development

Autism spectrum disorder (ASD) is a largely heritable, multistage prenatal disorder that impacts a child's ability to perceive and react to social information. Most ASD risk genes are expressed prenatally in many ASD-relevant brain regions and fall into two categories: broadly expressed regulatory genes that are expressed in the brain and other organs, and brain-specific genes. In trimesters one to three (Epoch-1), one set of broadly expressed (the majority) and brain-specific risk genes disrupts cell proliferation, neurogenesis, migration, and cell fate, while in trimester three and early postnatally (Epoch-2) another set (the majority being brain specific) disrupts neurite outgrowth, synaptogenesis, and the 'wiring' of the cortex. A proposed model is that upstream, highly interconnected regulatory ASD gene mutations disrupt transcriptional programs or signaling pathways resulting in dysregulation of downstream processes such as proliferation, neurogenesis, synaptogenesis, and neural activity. Dysregulation of signaling pathways is correlated with ASD social symptom severity. Since the majority of ASD risk genes are broadly expressed, many ASD individuals may benefit by being treated as having a broader medical disorder. An important future direction is the noninvasive study of ASD cell biology.

Human in vitro models for understanding mechanisms of autism spectrum disorder

Early brain development is a critical epoch for the development of autism spectrum disorder (ASD). In vivo animal models have, until recently, been the principal tool used to study early brain development and the changes occurring in neurodevelopmental disorders such as ASD. In vitro models of brain development represent a significant advance in the field. Here, we review the main methods available to study human brain development in vitro and the applications of these models for studying ASD and other psychiatric disorders. We discuss the main findings from stem cell models to date focusing on cell cycle and proliferation, cell death, cell differentiation and maturation, and neuronal signaling and synaptic stimuli. To be able to generalize the results from these studies, we propose a framework of experimental design and power considerations for using in vitro models to study ASD. These include both technical issues such as reproducibility and power analysis and conceptual issues such as the brain region and cell types being modeled.

Transcriptome signatures from discordant sibling pairs reveal changes in peripheral blood immune cell composition in Autism Spectrum Disorder

Notwithstanding several research efforts in the past years, robust and replicable molecular signatures for autism spectrum disorders from peripheral blood remain elusive. The available literature on blood transcriptome in ASD suggests that through accurate experimental design it is possible to extract important information on the disease pathophysiology at the peripheral level. Here we exploit the availability of a resource for molecular biomarkers in ASD, the Italian Autism Network (ITAN) collection, for the investigation of transcriptomic signatures in ASD based on a discordant sibling pair design. Whole blood samples from 75 discordant sibling pairs selected from the ITAN network where submitted to RNASeq analysis and data analyzed by complementary approaches. Overall, differences in gene expression between affected and unaffected siblings were small. In order to assess the contribution of differences in the relative proportion of blood cells between discordant siblings, we have applied two different cell deconvolution algorithms, showing that the observed molecular signatures mainly reflect changes in peripheral blood immune cell composition, in particular NK cells. The results obtained by the cell deconvolution approach are supported by the analysis performed by WGCNA. Our report describes the largest differential gene expression profiling in peripheral blood of ASD subjects and controls conducted by RNASeq. The observed signatures are consistent with the hypothesis of immune alterations in autism and an increased risk of developing autism in subjects exposed to prenatal infections or stress. Our study also points to a potential role of NMUR1, HMGB3, and PTPRN2 in ASD.

Transcriptomic Analysis Reveals Abnormal Expression of Prion Disease Gene Pathway in Brains from Patients with Autism Spectrum Disorders

The role of infections in the pathogenesis of autism spectrum disorder (ASD) is still controversial. In this study, we aimed to evaluate markers of infections and immune activation in ASD by performing a meta-analysis of publicly available whole-genome transcriptomic datasets of brain samples from autistic patients and otherwise normal people. Among the differentially expressed genes, no significant enrichment was observed for infectious diseases previously associated with ASD, including herpes simplex virus-1 (HSV-1), cytomegalovirus and Epstein–Barr virus in brain samples, nor was it found in peripheral blood from ASD patients. Interestingly, a significant number of genes belonging to the “prion diseases” pathway were found to be modulated in our ASD brain meta-analysis. Overall, our data do not support an association between infection and ASD. However, the data do provide support for the involvement of pathways related to other neurodegenerative diseases and give input to uncover novel pathogenetic mechanisms underlying ASD.

New Horizons for Molecular Genetics Diagnostic and Research in Autism Spectrum Disorder

Autism spectrum disorder (ASD) is a highly heritable, heterogeneous, and complex pervasive neurodevelopmental disorder (PND) characterized by distinctive abnormalities of human cognitive functions, social interaction, and speech development.Nowadays, several genetic changes including chromosome abnormalities, genetic variations, transcriptional epigenetics, and noncoding RNA have been identified in ASD. However, the association between these genetic modifications and ASDs has not been confirmed yet.The aim of this review is to summarize the key findings in ASD from genetic viewpoint that have been identified from the last few decades of genetic and molecular research.

Default mode-visual network hypoconnectivity in an autism subtype with pronounced social visual engagement difficulties

Social visual engagement difficulties are hallmark early signs of autism (ASD) and are easily quantified using eye tracking methods. However, it is unclear how these difficulties are linked to atypical early functional brain organization in ASD. With resting state fMRI data in a large sample of ASD toddlers and other non-ASD comparison groups, we find ASD-related functional hypoconnnectivity between ‘social brain’ circuitry such as the default mode network (DMN) and visual and attention networks. An eye tracking-identified ASD subtype with pronounced early social visual engagement difficulties (GeoPref ASD) is characterized by marked DMN-occipito-temporal cortex (OTC) hypoconnectivity. Increased DMN-OTC hypoconnectivity is also related to increased severity of social-communication difficulties, but only in GeoPref ASD. Early and pronounced social-visual circuit hypoconnectivity is a key underlying neurobiological feature describing GeoPref ASD and may be critical for future social-communicative development and represent new treatment targets for early intervention in these individuals.

Many parents of children with autism spectrum disorder (ASD) spot the first signs when their child is still a toddler, by noticing that their child is less interested than other toddlers in people and in social play. These early differences in behavior can have long-term implications for brain development. The brains of toddlers with little interest in social stimuli will receive less social input than those of other toddlers. This will make it even harder for the brain to develop the circuits required to support social skills.

But even among children with ASD, there are large differences in children's interest in the social world. One way of measuring these differences is to track eye movements. Lombardo et al. presented toddlers with and without ASD with images of moving colorful geometric shapes next to videos of dancing children. The majority of toddlers, including most of those with ASD, spent more time looking at the children than the shapes. But about 20% of the toddlers with ASD spent most of their time looking at the shapes. These toddlers also had the most severe social symptoms.

To find out why, Lombardo et al. measured the toddlers' brain activity while they slept. During sleep, or when at rest, the brain shows stereotyped patterns of activity. Groups of brain regions that work together – such as those involved in vision – fire in synchrony. Lombardo et al. found that toddlers who preferred looking at shapes over people showed different patterns of brain activity while asleep compared to other children. In the toddlers who preferred shapes, brain networks involved in social skills were less likely to coordinate their activity with networks that support vision and attention.

These findings suggest there may be multiple subtypes of ASD, with different symptoms resulting from different patterns of brain activity. At present, all children who receive a diagnosis of ASD receive much the same behavioral therapy. But in the future, studies of brain networks could allow children to receive more specific diagnoses. This could in turn lead to more effective and personalized treatments.

Multiscale Neuroscience of Psychiatric Disorders

The human brain comprises a multiscale network with multiple levels of organization. Neurons with dendritic and axonal connections form the microscale fabric of brain circuitry, and macroscale brain regions and white matter connections form the infrastructure for system-level brain communication and information integration. In this review, we discuss the emerging trend of multiscale neuroscience, the multidisciplinary field that brings together data from these different levels of nervous system organization to form a better understanding of between-scale relationships of brain structure, function, and behavior in health and disease. We provide a broad overview of this developing field and discuss recent findings of exemplary multiscale neuroscience studies that illustrate the importance of studying cross-scale interactions among the genetic, molecular, cellular, and macroscale levels of brain circuitry and connectivity and behavior. We particularly consider a central, overarching goal of these multiscale neuroscience studies of human brain connectivity: to obtain insight into how disease-related alterations at one level of organization may underlie alterations observed at other scales of brain network organization in mental disorders. We conclude by discussing the current limitations, challenges, and future directions of the field.

A perturbed gene network containing PI3K-AKT, RAS-ERK and WNT-β-catenin pathways in leukocytes is linked to ASD genetics and symptom severity

Hundreds of genes are implicated in autism spectrum disorder (ASD) but the mechanisms through which they contribute to ASD pathophysiology remain elusive. Here, we analyzed leukocyte transcriptomics from 1–4 year-old male toddlers with ASD or typical development from the general population. We discovered a perturbed gene network that includes genes that are highly expressed during fetal brain development and which is dysregulated in hiPSC-derived neuron models of ASD. High-confidence ASD risk genes emerge as upstream regulators of the network, and many risk genes may impact the network by modulating RAS/ERK, PI3K/AKT, and WNT/β-catenin signaling pathways. We found that the degree of dysregulation in this network correlated with the severity of ASD symptoms in the toddlers. These results demonstrate how the heterogeneous genetics of ASD may dysregulate a core network to influence brain development at prenatal and very early postnatal ages and, thereby, the severity of later ASD symptoms.

An integrated transcriptomic analysis of autism spectrum disorder

Autism spectrum disorder (ASD) is not a single disease but a set of disorders. To find clues of ASD pathogenesis in transcriptomic data, we performed an integrated transcriptomic analysis of ASD. After screening based on several standards in Gene Expression Omnibus (GEO) database, we obtained 11 series of transcriptomic data of different human tissues of ASD patients and healthy controls. Multidimensional scaling analysis revealed that datasets from the same tissue had bigger similarity than from different tissues. Functional enrichment analysis demonstrated that differential expressed genes were significantly enriched in inflammation/immune response, mitochondrion-related function and oxidative phosphorylation. Interestingly, genes enriched in inflammation/immune response were up-regulated in the brain tissues and down-regulated in the blood. In addition, drug prediction provided several compounds which might reverse gene expression profiles of ASD patients. And we also replicated the methods and criteria of transcriptomic analysis with datasets of ASD animal models and healthy controls, the results from animal models consolidated the results of transcriptomic analysis of ASD human tissues. In general, the results of our study may provide researchers a new sight of understanding the etiology of ASD and clinicians the possibilities of developing medical therapies.

Network-Based Integrative Analysis of Genomics, Epigenomics and Transcriptomics in Autism Spectrum Disorders

Current studies suggest that autism spectrum disorders (ASDs) may be caused by many genetic factors. In fact, collectively considering multiple studies aimed at characterizing the basic pathophysiology of ASDs, a large number of genes has been proposed. Addressing the problem of molecular data interpretation using gene networks helps to explain genetic heterogeneity in terms of shared pathways. Besides, the integrative analysis of multiple omics has emerged as an approach to provide a more comprehensive view of a disease. In this work, we carry out a network-based meta-analysis of the genes reported as associated with ASDs by studies that involved genomics, epigenomics, and transcriptomics. Collectively, our analysis provides a prioritization of the large number of genes proposed to be associated with ASDs, based on genes' relevance within the intracellular circuits, the strength of the supporting evidence of association with ASDs, and the number of different molecular alterations affecting genes. We discuss the presence of the prioritized genes in the SFARI (Simons Foundation Autism Research Initiative) database and in gene networks associated with ASDs by other investigations. Lastly, we provide the full results of our analyses to encourage further studies on common targets amenable to therapy.

Downregulation of an Evolutionary Young miR-1290 in an iPSC-Derived Neural Stem Cell Model of Autism Spectrum Disorder

The identification of several evolutionary young miRNAs, which arose in primates, raised several possibilities for the role of such miRNAs in human-specific disease processes. We previously have identified an evolutionary young miRNA, miR-1290, to be essential in neural stem cell proliferation and neuronal differentiation. Here, we show that miR-1290 is significantly downregulated during neuronal differentiation in reprogrammed induced pluripotent stem cell- (iPSC-) derived neurons obtained from idiopathic autism spectrum disorder (ASD) patients. Further, we identified that miR-1290 is actively released into extracellular vesicles. Supplementing ASD patient-derived neural stem cells (NSCs) with conditioned media from differentiated control-NSCs spiked with "artificial EVs" containing synthetic miR-1290 oligonucleotides significantly rescued differentiation deficits in ASD cell lines. Based on our earlier published study and the observations from the data presented here, we conclude that miR-1290 regulation could play a critical role during neuronal differentiation in early brain development.

Species-specific maturation profiles of human, chimpanzee and bonobo neural cells

Comparative analyses of neuronal phenotypes in closely related species can shed light on neuronal changes occurring during evolution. The study of post-mortem brains of nonhuman primates (NHPs) has been limited and often does not recapitulate important species-specific developmental hallmarks. We utilize induced pluripotent stem cell (iPSC) technology to investigate the development of cortical pyramidal neurons following migration and maturation of cells grafted in the developing mouse cortex. Our results show differential migration patterns in human neural progenitor cells compared to those of chimpanzees and bonobos both in vitro and in vivo, suggesting heterochronic changes in human neurons. The strategy proposed here lays the groundwork for further comparative analyses between humans and NHPs and opens new avenues for understanding the differences in the neural underpinnings of cognition and neurological disease susceptibility between species.

The ASD Living Biology: from cell proliferation to clinical phenotype

Autism spectrum disorder (ASD) has captured the attention of scientists, clinicians and the lay public because of its uncertain origins and striking and unexplained clinical heterogeneity. Here we review genetic, genomic, cellular, postmortem, animal model, and cell model evidence that shows ASD begins in the womb. This evidence leads to a new theory that ASD is a multistage, progressive disorder of brain development, spanning nearly all of prenatal life. ASD can begin as early as the 1st and 2nd trimester with disruption of cell proliferation and differentiation. It continues with disruption of neural migration, laminar disorganization, altered neuron maturation and neurite outgrowth, disruption of synaptogenesis and reduced neural network functioning. Among the most commonly reported high-confidence ASD (hcASD) genes, 94% express during prenatal life and affect these fetal processes in neocortex, amygdala, hippocampus, striatum and cerebellum. A majority of hcASD genes are pleiotropic, and affect proliferation/differentiation and/or synapse development. Proliferation and subsequent fetal stages can also be disrupted by maternal immune activation in the 1st trimester. Commonly implicated pathways, PI3K/AKT and RAS/ERK, are also pleiotropic and affect multiple fetal processes from proliferation through synapse and neural functional development. In different ASD individuals, variation in how and when these pleiotropic pathways are dysregulated, will lead to different, even opposing effects, producing prenatal as well as later neural and clinical heterogeneity. Thus, the pathogenesis of ASD is not set at one point in time and does not reside in one process, but rather is a cascade of prenatal pathogenic processes in the vast majority of ASD toddlers. Despite this new knowledge and theory that ASD biology begins in the womb, current research methods have not provided individualized information: What are the fetal processes and early-age molecular and cellular differences that underlie ASD in each individual child? Without such individualized knowledge, rapid advances in biological-based diagnostic, prognostic, and precision medicine treatments cannot occur. Missing, therefore, is what we call ASD Living Biology. This is a conceptual and paradigm shift towards a focus on the abnormal prenatal processes underlying ASD within each living individual. The concept emphasizes the specific need for foundational knowledge of a living child's development from abnormal prenatal beginnings to early clinical stages. The ASD Living Biology paradigm seeks this knowledge by linking genetic and in vitro prenatal molecular, cellular and neural measurements with in vivo post-natal molecular, neural and clinical presentation and progression in each ASD child. We review the first such study, which confirms the multistage fetal nature of ASD and provides the first in vitro fetal-stage explanation for in vivo early brain overgrowth. Within-child ASD Living Biology is a novel research concept we coin here that advocates the integration of in vitro prenatal and in vivo early post-natal information to generate individualized and group-level explanations, clinically useful prognoses, and precision medicine approaches that are truly beneficial for the individual infant and toddler with ASD.

Large-scale associations between the leukocyte transcriptome and BOLD responses to speech differ in autism early language outcome subtypes

Heterogeneity in early language development in autism spectrum disorders (ASD) is clinically important and may reflect neurobiologically distinct subtypes. Here we identify a large-scale association between multiple coordinated blood leukocyte gene co-expression modules and multivariate functional neuroimaging (fMRI) response to speech. Gene co-expression modules associated with multivariate fMRI response to speech are different for all pairwise comparisons between typically developing toddlers and toddlers with ASD and either poor versus good early language outcome. Associated co-expression modules are enriched in genes that are broadly expressed in the brain and many other tissues. These co-expression modules are also enriched for ASD, prenatal, human-specific and language-relevant genes. This work highlights distinctive neurobiology in ASD subtypes with different early language outcomes that is present well before such outcomes are known. Associations between neuroimaging measures and gene expression levels in blood leukocytes may offer a unique in-vivo window into identifying brain-relevant molecular mechanisms in ASD.

Parallel Development of Chromatin Patterns, Neuron Morphology, and Connections: Potential for Disruption in Autism

The phenotype of neurons and their connections depend on complex genetic and epigenetic processes that regulate the expression of genes in the nucleus during development and throughout life. Here we examined the distribution of nuclear chromatin patters in relation to the epigenetic landscape, phenotype and connections of neurons with a focus on the primate cerebral cortex. We show that nuclear patterns of chromatin in cortical neurons are related to neuron size and cortical connections. Moreover, we point to evidence that reveals an orderly sequence of events during development, linking chromatin and gene expression patterns, neuron morphology, function, and connections across cortical areas and layers. Based on this synthesis, we posit that systematic studies of changes in chromatin patterns and epigenetic marks across cortical areas will provide novel insights on the development and evolution of cortical networks, and their disruption in connectivity disorders of developmental origin, like autism. Achieving this requires embedding and interpreting genetic, transcriptional, and epigenetic studies within a framework that takes into consideration distinct types of neurons, local circuit interactions, and interareal pathways. These features vary systematically across cortical areas in parallel with laminar structure and are differentially affected in disorders. Finally, based on evidence that autism-associated genetic polymorphisms are especially prominent in excitatory neurons and connectivity disruption affects mostly limbic cortices, we employ this systematic approach to propose novel, targeted studies of projection neurons in limbic areas to elucidate the emergence and time-course of developmental disruptions in autism.

Investigation of epigenetic regulatory networks associated with autism spectrum disorder (ASD) by integrated global LINE-1 methylation and gene expression profiling analyses

BACKGROUND

The exact cause and mechanisms underlying the pathobiology of autism spectrum disorder (ASD) remain unclear. Dysregulation of long interspersed element-1 (LINE-1) has been reported in the brains of ASD-like mutant mice and ASD brain tissues. However, the role and methylation of LINE-1 in individuals with ASD remain unclear. In this study, we aimed to investigate whether LINE-1 insertion is associated with differentially expressed genes (DEGs) and to assess LINE-1 methylation in ASD.

METHODS

To identify DEGs associated with LINE-1 in ASD, we reanalyzed previously published transcriptome profiles and overlapped them with the list of LINE-1-containing genes from the TranspoGene database. An Ingenuity Pathway Analysis (IPA) of DEGs associated with LINE-1 insertion was conducted. DNA methylation of LINE-1 was assessed via combined bisulfite restriction analysis (COBRA) of lymphoblastoid cell lines from ASD individuals and unaffected individuals, and the methylation levels were correlated with the expression levels of LINE-1 and two LINE-1-inserted DEGs, C1orf27 and ARMC8.

RESULTS

We found that LINE-1 insertion was significantly associated with DEGs in ASD. The IPA showed that LINE-1-inserted DEGs were associated with ASD-related mechanisms, including sex hormone receptor signaling and axon guidance signaling. Moreover, we observed that the LINE-1 methylation level was significantly reduced in lymphoblastoid cell lines from ASD individuals with severe language impairment and was inversely correlated with the transcript level. The methylation level of LINE-1 was also correlated with the expression of the LINE-1-inserted DEG C1orf27 but not ARMC8.

CONCLUSIONS

In ASD individuals with severe language impairment, LINE-1 methylation was reduced and correlated with the expression levels of LINE-1 and the LINE-1-inserted DEG C1orf27. Our findings highlight the association of LINE-1 with DEGs in ASD blood samples and warrant further investigation. The molecular mechanisms of LINE-1 and the effects of its methylation in ASD pathobiology deserve further study.

Integrated genome-wide Alu methylation and transcriptome profiling analyses reveal novel epigenetic regulatory networks associated with autism spectrum disorder

Background

Alu elements are a group of repetitive elements that can influence gene expression through CpG residues and transcription factor binding. Altered gene expression and methylation profiles have been reported in various tissues and cell lines from individuals with autism spectrum disorder (ASD). However, the role of Alu elements in ASD remains unclear. We thus investigated whether Alu elements are associated with altered gene expression profiles in ASD.

Methods

We obtained five blood-based gene expression profiles from the Gene Expression Omnibus database and human Alu-inserted gene lists from the TranspoGene database. Differentially expressed genes (DEGs) in ASD were identified from each study and overlapped with the human Alu-inserted genes. The biological functions and networks of Alu-inserted DEGs were then predicted by Ingenuity Pathway Analysis (IPA). A combined bisulfite restriction analysis of lymphoblastoid cell lines (LCLs) derived from 36 ASD and 20 sex- and age-matched unaffected individuals was performed to assess the global DNA methylation levels within Alu elements, and the Alu expression levels were determined by quantitative RT-PCR.

Results

In ASD blood or blood-derived cells, 320 Alu-inserted genes were reproducibly differentially expressed. Biological function and pathway analysis showed that these genes were significantly associated with neurodevelopmental disorders and neurological functions involved in ASD etiology. Interestingly, estrogen receptor and androgen signaling pathways implicated in the sex bias of ASD, as well as IL-6 signaling and neuroinflammation signaling pathways, were also highlighted. Alu methylation was not significantly different between the ASD and sex- and age-matched control groups. However, significantly altered Alu methylation patterns were observed in ASD cases sub-grouped based on Autism Diagnostic Interview-Revised scores compared with matched controls. Quantitative RT-PCR analysis of Alu expression also showed significant differences between ASD subgroups. Interestingly, Alu expression was correlated with methylation status in one phenotypic ASD subgroup.

Conclusion

Alu methylation and expression were altered in LCLs from ASD subgroups. Our findings highlight the association of Alu elements with gene dysregulation in ASD blood samples and warrant further investigation. Moreover, the classification of ASD individuals into subgroups based on phenotypes may be beneficial and could provide insights into the still unknown etiology and the underlying mechanisms of ASD.

Electronic supplementary material

The online version of this article (10.1186/s13229-018-0213-9) contains supplementary material, which is available to authorized users.

Maternal immune activation dysregulation of the fetal brain transcriptome and relevance to the pathophysiology of autism spectrum disorder

Maternal immune activation (MIA) via infection during pregnancy is known to increase risk for autism spectrum disorder (ASD). However, it is unclear how MIA disrupts fetal brain gene expression in ways that may explain this increased risk. Here we examine how MIA dysregulates rat fetal brain gene expression (at a time point analogous to the end of the first trimester of human gestation) in ways relevant to ASD-associated pathophysiology. MIA downregulates expression of ASD-associated genes, with the largest enrichments in genes known to harbor rare highly penetrant mutations. MIA also downregulates expression of many genes also known to be persistently downregulated in the ASD cortex later in life and which are canonically known for roles in affecting prenatally late developmental processes at the synapse. Transcriptional and translational programs that are downstream targets of highly ASD-penetrant FMR1 and CHD8 genes are also heavily affected by MIA. MIA strongly upregulates expression of a large number of genes involved in translation initiation, cell cycle, DNA damage and proteolysis processes that affect multiple key neural developmental functions. Upregulation of translation initiation is common to and preserved in gene network structure with the ASD cortical transcriptome throughout life and has downstream impact on cell cycle processes. The cap-dependent translation initiation gene, EIF4E, is one of the most MIA-dysregulated of all ASD-associated genes and targeted network analyses demonstrate prominent MIA-induced transcriptional dysregulation of mTOR and EIF4E-dependent signaling. This dysregulation of translation initiation via alteration of the Tsc2-mTor-Eif4e axis was further validated across MIA rodent models. MIA may confer increased risk for ASD by dysregulating key aspects of fetal brain gene expression that are highly relevant to pathophysiology affecting ASD.

Network‐based method for detecting dysregulated pathways in glioblastoma cancer

The knowledge on the biological molecular mechanisms underlying cancer is important for the precise diagnosis and treatment of cancer patients. Detecting dysregulated pathways in cancer can provide insights into the mechanism of cancer and help to detect novel drug targets. Based on the wide existing mutual exclusivity among mutated genes and the interrelationship between gene mutations and expression changes, this study presents a network‐based method to detect the dysregulated pathways from gene mutations and expression data of the glioblastoma cancer. First, the authors construct a gene network based on mutual exclusivity between each pair of genes and the interaction between gene mutations and expression changes. Then they detect all complete subgraphs using CFinder clustering algorithm in the constructed gene network. Next, the two gene sets whose overlapping scores are above a specific threshold are merged. Finally, they obtain two dysregulated pathways in which there are glioblastoma‐related multiple genes which are closely related to the two subtypes of glioblastoma. The results show that one dysregulated pathway revolving around epidermal growth factor receptor is likely to be associated with the primary subtype of glioblastoma, and the other dysregulated pathway revolving around TP53 is likely to be associated with the secondary subtype of glioblastoma.

Enriched expression of genes associated with autism spectrum disorders in human inhibitory neurons

Autism spectrum disorder (ASD) is highly heritable but genetically heterogeneous. The affected neural circuits and cell types remain unclear and may vary at different developmental stages. By analyzing multiple sets of human single cell transcriptome profiles, we found that ASD candidates showed relatively enriched gene expression in neurons, especially in inhibitory neurons. ASD candidates were also more likely to be the hubs of the co-expression gene module that is highly expressed in inhibitory neurons, a feature not detected for excitatory neurons. In addition, we found that upregulated genes in multiple ASD cortex samples were enriched with genes highly expressed in inhibitory neurons, suggesting a potential increase of inhibitory neurons and an imbalance in the ratio between excitatory and inhibitory neurons in ASD brains. Furthermore, the downstream targets of several ASD candidates, such as CHD8, EHMT1 and SATB2, also displayed enriched expression in inhibitory neurons. Taken together, our analyses of single cell transcriptomic data suggest that inhibitory neurons may be a major neuron subtype affected by the disruption of ASD gene networks, providing single cell functional evidence to support the excitatory/inhibitory (E/I) imbalance hypothesis.

Toddlers later diagnosed with autism exhibit multiple structural abnormalities in temporal corpus callosum fibers

Interhemispheric functional connectivity abnormalities are often reported in autism and it is thus not surprising that structural defects of the corpus callosum (CC) are consistently found using both traditional MRI and DTI techniques. Past DTI studies however, have subdivided the CC into 2 or 3 segments without regard for where fibers may project to within the cortex, thus placing limitations on our ability to understand the nature, timing and neurobehavioral impact of early CC abnormalities in autism. Leveraging a unique cohort of 97 toddlers (68 autism; 29 typical) we utilized a novel technique that identified seven CC tracts according to their cortical projections. Results revealed that younger (<2.5 years old), but not older toddlers with autism exhibited abnormally low mean, radial, and axial diffusivity values in the CC tracts connecting the occipital lobes and the temporal lobes. Fractional anisotropy and the cross sectional area of the temporal CC tract were significantly larger in young toddlers with autism. These findings indicate that water diffusion is more restricted and unidirectional in the temporal CC tract of young toddlers who develop autism. Such results may be explained by a potential overabundance of small caliber axons generated by excessive prenatal neural proliferation as proposed by previous genetic, animal model, and postmortem studies of autism. Furthermore, early diffusion measures in the temporal CC tract of the young toddlers were correlated with outcome measures of autism severity at later ages. These findings regarding the potential nature, timing, and location of early CC abnormalities in autism add to accumulating evidence, which suggests that altered inter-hemispheric connectivity, particularly across the temporal lobes, is a hallmark of the disorder.

Hierarchical cortical transcriptome disorganization in autism

Background

Autism spectrum disorders (ASD) are etiologically heterogeneous and complex. Functional genomics work has begun to identify a diverse array of dysregulated transcriptomic programs (e.g., synaptic, immune, cell cycle, DNA damage, WNT signaling, cortical patterning and differentiation) potentially involved in ASD brain abnormalities during childhood and adulthood. However, it remains unclear whether such diverse dysregulated pathways are independent of each other or instead reflect coordinated hierarchical systems-level pathology.

Methods

Two ASD cortical transcriptome datasets were re-analyzed using consensus weighted gene co-expression network analysis (WGCNA) to identify common co-expression modules across datasets. Linear mixed-effect models and Bayesian replication statistics were used to identify replicable differentially expressed modules. Eigengene network analysis was then utilized to identify between-group differences in how co-expression modules interact and cluster into hierarchical meta-modular organization. Protein-protein interaction analyses were also used to determine whether dysregulated co-expression modules show enhanced interactions.

Results

We find replicable evidence for 10 gene co-expression modules that are differentially expressed in ASD cortex. Rather than being independent non-interacting sources of pathology, these dysregulated co-expression modules work in synergy and physically interact at the protein level. These systems-level transcriptional signals are characterized by downregulation of synaptic processes coordinated with upregulation of immune/inflammation, response to other organism, catabolism, viral processes, translation, protein targeting and localization, cell proliferation, and vasculature development. Hierarchical organization of meta-modules (clusters of highly correlated modules) is also highly affected in ASD.

Conclusions

These findings highlight that dysregulation of the ASD cortical transcriptome is characterized by the dysregulation of multiple coordinated transcriptional programs producing synergistic systems-level effects that cannot be fully appreciated by studying the individual component biological processes in isolation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13229-017-0147-7) contains supplementary material, which is available to authorized users.

Cross-disorder comparative analysis of comorbid conditions reveals novel autism candidate genes

BACKGROUND

Numerous studies have highlighted the elevated degree of comorbidity associated with autism spectrum disorder (ASD). These comorbid conditions may add further impairments to individuals with autism and are substantially more prevalent compared to neurotypical populations. These high rates of comorbidity are not surprising taking into account the overlap of symptoms that ASD shares with other pathologies. From a research perspective, this suggests common molecular mechanisms involved in these conditions. Therefore, identifying crucial genes in the overlap between ASD and these comorbid disorders may help unravel the common biological processes involved and, ultimately, shed some light in the understanding of autism etiology.

RESULTS

In this work, we used a two-fold systems biology approach specially focused on biological processes and gene networks to conduct a comparative analysis of autism with 31 frequently comorbid disorders in order to define a multi-disorder subcomponent of ASD and predict new genes of potential relevance to ASD etiology. We validated our predictions by determining the significance of our candidate genes in high throughput transcriptome expression profiling studies. Using prior knowledge of disease-related biological processes and the interaction networks of the disorders related to autism, we identified a set of 19 genes not previously linked to ASD that were significantly differentially regulated in individuals with autism. In addition, these genes were of potential etiologic relevance to autism, given their enriched roles in neurological processes crucial for optimal brain development and function, learning and memory, cognition and social behavior.

CONCLUSIONS

Taken together, our approach represents a novel perspective of autism from the point of view of related comorbid disorders and proposes a model by which prior knowledge of interaction networks may enlighten and focus the genome-wide search for autism candidate genes to better define the genetic heterogeneity of ASD.