Towards a molecular characterization of autism spectrum disorders: an exome sequencing and systems approach

Towards understanding sex differences in autism spectrum disorders

Autism Spectrum Disorder (ASD) is a neurodevelopmental condition characterized by social deficits, repetitive behaviours and lack of empathy. Its significant genetic heritability and potential comorbidities often lead to diagnostic and therapeutic challenges. This review addresses the biological basis of ASD, focusing on the sex differences in gene expression and hormonal influences. ASD is more commonly diagnosed in males at a ratio of 4:1, indicating a potential oversight in female-specific ASD research and a risk of underdiagnosis in females. We consider how ASD manifests differently across sexes by exploring differential gene expression in female and male brains and consider how variations in steroid hormones influence ASD characteristics. Synaptic function, including excitation/inhibition ratio imbalance, is influenced by gene mutations and this is explored as a key factor in the cognitive and behavioural manifestations of ASD. We also discuss the role of micro RNAs (miRNAs) and highlight a novel mutation in miRNA-873, which affects a suite of key synaptic genes, neurexin, neuroligin, SHANK and post-synaptic density proteins, implicated in the pathology of ASD. Our review suggests that genetic predisposition, sex differences in brain gene expression, and hormonal factors significantly contribute to the presentation, identification and severity of ASD, necessitating sex-specific considerations in diagnosis and treatments. These findings advocate for personalized interventions to improve the outcomes for individuals with ASD.

Oligogenic inheritance in severe adult obesity

BACKGROUND/OBJECTIVE

The genetic architecture of extreme non-syndromic obesity in adults remains to be elucidated. A range of genes are known to cause monogenic obesity, but even when pathogenic mutations are present, there may be variable penetrance.

METHODS

Whole-exome sequencing (WES) was carried out on a 15-year-old male proband of Pakistani ancestry who had severe obesity. This was followed by family segregation analysis, using Sanger sequencing. We also undertook re-analysis of WES data from 91 unrelated adults with severe obesity (86% white European ancestry) from the Personalised Medicine for Morbid Obesity (PMMO) cohort, recruited from the UK National Health Service.

RESULTS

We identified an oligogenic mode of inheritance of obesity in the proband's family-this provided the impetus to reanalyze existing sequence data in a separate dataset. Analysis of PMMO participant data revealed two further patients who carried more than one rare, predicted-deleterious mutation in a known monogenic obesity gene. In all three cases, the genes involved had known autosomal dominant inheritance, with incomplete penetrance.

CONCLUSION

Oligogenic inheritance may explain some of the variable penetrance in Mendelian forms of obesity. We caution clinicians and researchers to avoid confining sequence analysis to individual genes and, in particular, not to stop looking when the first potentially-causative mutation is found.

Heterozygous deletion of the autism-associated gene CHD8 impairs synaptic function through widespread changes in gene expression and chromatin compaction

Whole-exome sequencing of autism spectrum disorder (ASD) probands and unaffected family members has identified many genes harboring de novo variants suspected to play a causal role in the disorder. Of these, chromodomain helicase DNA-binding protein 8 (CHD8) is the most recurrently mutated. Despite the prevalence of CHD8 mutations, we have little insight into how CHD8 loss affects genome organization or the functional consequences of these molecular alterations in neurons. Here, we engineered two isogenic human embryonic stem cell lines with CHD8 loss-of-function mutations and characterized differences in differentiated human cortical neurons. We identified hundreds of genes with altered expression, including many involved in neural development and excitatory synaptic transmission. Field recordings and single-cell electrophysiology revealed a 3-fold decrease in firing rates and synaptic activity in CHD8+/- neurons, as well as a similar firing-rate deficit in primary cortical neurons from Chd8+/- mice. These alterations in neuron and synapse function can be reversed by CHD8 overexpression. Moreover, CHD8+/- neurons displayed a large increase in open chromatin across the genome, where the greatest change in compaction was near autism susceptibility candidate 2 (AUTS2), which encodes a transcriptional regulator implicated in ASD. Genes with changes in chromatin accessibility and expression in CHD8+/- neurons have significant overlap with genes mutated in probands for ASD, intellectual disability, and schizophrenia but not with genes mutated in healthy controls or other disease cohorts. Overall, this study characterizes key molecular alterations in genome structure and expression in CHD8+/- neurons and links these changes to impaired neuronal and synaptic function.

A new neurodevelopmental disorder linked to heterozygous variants in UNC79

PURPOSE

The "NALCN channelosome" is an ion channel complex that consists of multiple proteins, including NALCN, UNC79, UNC80, and FAM155A. Only a small number of individuals with a neurodevelopmental syndrome have been reported with disease causing variants in NALCN and UNC80. However, no pathogenic UNC79 variants have been reported, and in vivo function of UNC79 in humans is largely unknown.

METHODS

We used international gene-matching efforts to identify patients harboring ultrarare heterozygous loss-of-function UNC79 variants and no other putative responsible genes. We used genetic manipulations in Drosophila and mice to test potential causal relationships between UNC79 variants and the pathology.

RESULTS

We found 6 unrelated and affected patients with UNC79 variants. Five patients presented with overlapping neurodevelopmental features, including mild to moderate intellectual disability and a mild developmental delay, whereas a single patient reportedly had normal cognitive and motor development but was diagnosed with epilepsy and autistic features. All displayed behavioral issues and 4 patients had epilepsy. Drosophila with UNC79 knocked down displayed induced seizure-like phenotype. Mice with a heterozygous loss-of-function variant have a developmental delay in body weight compared with wild type. In addition, they have impaired ability in learning and memory.

CONCLUSION

Our results demonstrate that heterozygous loss-of-function UNC79 variants are associated with neurologic pathologies.

Dysregulation of the Wnt/β-catenin signaling pathway via Rnf146 upregulation in a VPA-induced mouse model of autism spectrum disorder

Autism spectrum disorder (ASD) is a neurodevelopmental disorder associated with impaired social behavior and communication, repetitive behaviors, and restricted interests. In addition to genetic factors, environmental factors such as prenatal drug exposure contribute to the development of ASD. However, how those prenatal factors induce behavioral deficits in the adult stage is not clear. To elucidate ASD pathogenesis at the molecular level, we performed a high-resolution mass spectrometry-based quantitative proteomic analysis on the prefrontal cortex (PFC) of mice exposed to valproic acid (VPA) in utero, a widely used animal model of ASD. Differentially expressed proteins (DEPs) in VPA-exposed mice showed significant overlap with ASD risk genes, including differentially expressed genes from the postmortem cortex of ASD patients. Functional annotations of the DEPs revealed significant enrichment in the Wnt/β-catenin signaling pathway, which is dysregulated by the upregulation of Rnf146 in VPA-exposed mice. Consistently, overexpressing Rnf146 in the PFC impaired social behaviors and altered the Wnt signaling pathway in adult mice. Furthermore, Rnf146-overexpressing PFC neurons showed increased excitatory synaptic transmission, which may underlie impaired social behavior. These results demonstrate that Rnf146 is critical for social behavior and that dysregulation of Rnf146 underlies social deficits in VPA-exposed mice.

Mice Mutated in the Third Fibronectin Domain of L1 Show Enhanced Hippocampal Neuronal Cell Death, Astrogliosis and Alterations in Behavior

Adhesion molecules play major roles in cell proliferation, migration, survival, neurite outgrowth and synapse formation during nervous system development and in adulthood. The neural cell adhesion molecule L1 contributes to these functions during development and in synapse formation and synaptic plasticity after trauma in adulthood. Mutations of L1 in humans result in L1 syndrome, which is associated with mild-to-severe brain malformations and mental disabilities. Furthermore, mutations in the extracellular domain were shown to cause a severe phenotype more often than mutations in the intracellular domain. To explore the outcome of a mutation in the extracellular domain, we generated mice with disruption of the dibasic sequences RK and KR that localize to position 858RKHSKR863 in the third fibronectin type III domain of murine L1. These mice exhibit alterations in exploratory behavior and enhanced marble burying activity. Mutant mice display higher numbers of caspase 3-positive neurons, a reduced number of principle neurons in the hippocampus, and an enhanced number of glial cells. Experiments suggest that disruption of the dibasic sequence in L1 results in subtle impairments in brain structure and functions leading to obsessive-like behavior in males and reduced anxiety in females.

Cntnap2-dependent molecular networks in autism spectrum disorder revealed through an integrative multi-omics analysis

Autism spectrum disorder (ASD) is a major neurodevelopmental disorder in which patients present with core symptoms of social communication impairment, restricted interest, and repetitive behaviors. Although various studies have been performed to identify ASD-related mechanisms, ASD pathology is still poorly understood. CNTNAP2 genetic variants have been found that represent ASD genetic risk factors, and disruption of Cntnap2 expression has been associated with ASD phenotypes in mice. In this study, we performed an integrative multi-omics analysis by combining quantitative proteometabolomic data obtained with Cntnap2 knockout (KO) mice with multi-omics data obtained from ASD patients and forebrain organoids to elucidate Cntnap2-dependent molecular networks in ASD. To this end, a mass spectrometry-based proteometabolomic analysis of the medial prefrontal cortex in Cntnap2 KO mice led to the identification of Cntnap2-associated molecular features, and these features were assessed in combination with multi-omics data obtained on the prefrontal cortex in ASD patients to identify bona fide ASD cellular processes. Furthermore, a reanalysis of single-cell RNA sequencing data obtained from forebrain organoids derived from patients with CNTNAP2-associated ASD revealed that the aforementioned identified ASD processes were mainly linked to excitatory neurons. On the basis of these data, we constructed Cntnap2-associated ASD network models showing mitochondrial dysfunction, axonal impairment, and synaptic activity. Our results may shed light on the Cntnap2-dependent molecular networks in ASD.

A rigorous in silico genomic interrogation at 1p13.3 reveals 16 autosomal dominant candidate genes in syndromic neurodevelopmental disorders

Genome-wide chromosomal microarray is extensively used to detect copy number variations (CNVs), which can diagnose microdeletion and microduplication syndromes. These small unbalanced chromosomal structural rearrangements ranging from 1 kb to 10 Mb comprise up to 15% of human mutations leading to monogenic or contiguous genomic disorders. Albeit rare, CNVs at 1p13.3 cause a variety of neurodevelopmental disorders (NDDs) including development delay (DD), intellectual disability (ID), autism, epilepsy, and craniofacial anomalies (CFA). Most of the 1p13.3 CNV cases reported in the pre-microarray era encompassed a large number of genes and lacked the demarcating genomic coordinates, hampering the discovery of positional candidate genes within the boundaries. In this study, we present four subjects with 1p13.3 microdeletions displaying DD, ID, autism, epilepsy, and CFA. In silico comparative genomic mapping with three previously reported subjects with CNVs and 22 unreported DECIPHER CNV cases has resulted in the identification of four different sub-genomic loci harboring five positional candidate genes for DD, ID, and CFA at 1p13.3. Most of these genes have pathogenic variants reported, and their interacting genes are involved in NDDs. RT-qPCR in various human tissues revealed a high expression pattern in the brain and fetal brain, supporting their functional roles in NDDs. Interrogation of variant databases and interacting protein partners led to the identification of another set of 11 potential candidate genes, which might have been dysregulated by the position effect of these CNVs at 1p13.3. Our studies define 1p13.3 as a genomic region harboring 16 NDD candidate genes and underscore the critical roles of small CNVs in in silico comparative genomic mapping for disease gene discovery. Our candidate genes will help accelerate the isolation of pathogenic heterozygous variants from exome/genome sequencing (ES/GS) databases.

Mitochondrial developmental encephalopathy with bilateral optic neuropathy related to homozygous variants in IMMT gene

IMMT gene codes for mitofilin, a mitochondrial inner membrane protein that regulates the morphology of mitochondrial cristae. The phenotype associated with mutations in this gene has not been yet established, but functional studies carried out show that its loss causes a mitochondrial alteration, both in the morphology of the mitochondrial crests and in their function. We present two cousins from an extended highly consanguineous family with developmental encephalopathy, hypotonia, nystagmus due to optic neuropathy. The likely pathogenic homozygous c.895A>G (p.Lys299Glu) variant in the IMMT gene co-segregates with the disease and associates altered mitochondrial cristae observed by electron microscopy.

Developmental disabilities across the world: A scientometric review from 1936 to 2020

BACKGROUND

Developmental disabilities have been largely studied in the past years. Their etiological mechanisms have been underpinned to the interactions between genetic and environmental factors. These factors show variability across the world. Thus, it is important to understand where the set of knowledge obtained on developmental disabilities originates from and whether it is generalizable to low- and middle-income countries.

AIMS

This study aims to understand the origins of the available literature on developmental disabilities, keeping a focus on parenting, and identify the main trend of research.

METHODS AND PROCEDURE

A sample of 11,315 publications from 1936 to 2020 were collected from Scopus and a graphical country analysis was conducted. Furthermore, a qualitative approach enabled the clustering of references by keywords into four main areas: "Expression of the disorder", "Physiological Factors", "How it is studied" and "Environmental factors". For each area, a document co-citation analysis (DCA) on CiteSpace software was performed.

OUTCOMES AND RESULTS

Results highlight the leading role of North America in the study of developmental disabilities. Trends in the literature and the documents' scientific relevance are discussed in details.

CONCLUSIONS AND IMPLICATIONS

Results demand for investigation in different socio-economical settings to generalize our knowledge. What this paper adds? The current paper tries to provide insight into the origins of the literature on developmental disabilities with a focus on parenting, together with an analysis of the trends of research in the field. The paper consisted of a multi-disciplinary and multi-method review. In fact, the review tried to integrate the analysis of the relation between developmental disabilities with a closer look at the scientific contributions to the field across the world. Specifically, the paper integrates a total of 11,315 papers published on almost a century of research (from 1936 to 2020). An initial qualitative analysis on keywords was combined to a subsequent quantitative approach in order to maximize the comprehension of the impact of almost a century of scientific contributions. Specifically, documents were studied with temporal and structural metrics on a scientometric approach. This allowed the exploration of patterns within the literature available on Scopus in a quantitative way. This method not only assessed the importance of single documents within the network. As a matter of fact, the document co-citation analysis used on CiteSpace software provided insight into the relations existing between multiple documents in the field of research. As a result, the leading role of North America in the literature of developmental disabilities and parenting emerged. This was accompanied by the review of the main trends of research within the existing literature.

Extending the clinical phenotype of SPTAN1: From DEE5 to migraine, epilepsy, and subependymal heterotopias without intellectual disability

Mutations in SPTAN1 gene, encoding the nonerythrocyte αII-spectrin, are responsible for a severe developmental and epileptic encephalopathy (DEE5) and a wide spectrum of neurodevelopmental disorders, as epilepsy with or without intellectual disability (ID) or ID with cerebellar syndrome. A certain genotype-phenotype correlation has been proposed according to the type and location of the mutation. Herein, we report three novel cases with de novo SPTAN1 mutations, one of them associated to a mild phenotype not previously described. They range from (1) severe developmental encephalopathy with ataxia and a mild cerebellar atrophy, without epilepsy; (2) moderate intellectual disability, severe language delay, ataxia and tremor; (3) normal intelligence, chronic migraine, and generalized tonic-clonic seizures. Remarkably, all these patients showed brain MRI abnormalities, being of special interest the subependymal heterotopias detected in the latter patient. Thus we extend the SPTAN1-related phenotypic spectrum, both in its radiological and clinical involvement. Furthermore, after systematic analysis of all the patients so far reported, we noted an excess of male versus female patients (20:9, p = 0.04), more pronounced among the milder phenotypes. Consequently, some protection factor might be suspected among female carriers, which if confirmed should be considered when establishing the pathogenicity of milder genetic variants in this gene.

Genetic architecture of autism spectrum disorder: Lessons from large-scale genomic studies

Autism spectrum disorder (ASD) is a neurodevelopmental disorder with a strong genetic component. Recently developed genomic technologies, including microarray and next-generation sequencing (NGS), have enabled researchers to genetic analyses aimed at identifying genetic variations associated with ASD and to elucidate the genetic architecture of the disorder. Large-scale microarray, exome sequencing analyses, and robust statistical methods have resulted in successful gene discovery and identification of high-confidence ASD genes from among de novo and inherited variants. Efforts have been made to understand the genetic architecture of ASD using whole-genome sequencing and genome-wide association studies aimed at identifying noncoding mutations and common variants associated with ASD. In addition, the development of systems biology approaches has resulted in the integration of genetic findings with functional genomic datasets, thereby providing a unique insight into the functional convergence of ASD risk genes and their neurobiology. In this review, we summarize the latest findings of ASD genetic studies involving large cohorts and discuss their implications in ASD neurobiology and in clinical practice.

Roles of HIF and 2-Oxoglutarate-Dependent Dioxygenases in Controlling Gene Expression in Hypoxia

Simple Summary

Hypoxia—reduction in oxygen availability—plays key roles in both physiological and pathological processes. Given the importance of oxygen for cell and organism viability, mechanisms to sense and respond to hypoxia are in place. A variety of enzymes utilise molecular oxygen, but of particular importance to oxygen sensing are the 2-oxoglutarate (2-OG) dependent dioxygenases (2-OGDs). Of these, Prolyl-hydroxylases have long been recognised to control the levels and function of Hypoxia Inducible Factor (HIF), a master transcriptional regulator in hypoxia, via their hydroxylase activity. However, recent studies are revealing that such dioxygenases are involved in almost all aspects of gene regulation, including chromatin organisation, transcription and translation.

Abstract

Hypoxia—reduction in oxygen availability—plays key roles in both physiological and pathological processes. Given the importance of oxygen for cell and organism viability, mechanisms to sense and respond to hypoxia are in place. A variety of enzymes utilise molecular oxygen, but of particular importance to oxygen sensing are the 2-oxoglutarate (2-OG) dependent dioxygenases (2-OGDs). Of these, Prolyl-hydroxylases have long been recognised to control the levels and function of Hypoxia Inducible Factor (HIF), a master transcriptional regulator in hypoxia, via their hydroxylase activity. However, recent studies are revealing that dioxygenases are involved in almost all aspects of gene regulation, including chromatin organisation, transcription and translation. We highlight the relevance of HIF and 2-OGDs in the control of gene expression in response to hypoxia and their relevance to human biology and health.

Dysfunction of NMDA receptors in neuronal models of an autism spectrum disorder patient with a DSCAM mutation and in Dscam-knockout mice

Heterogeneity in the etiopathology of autism spectrum disorders (ASD) limits the development of generic remedies, requires individualistic and patient-specific research. Recent progress in human-induced pluripotent stem cell (iPSC) technology provides a novel platform for modeling ASDs for studying complex neuronal phenotypes. In this study, we generated telencephalic induced neuronal (iN) cells from iPSCs derived from an ASD patient with a heterozygous point mutation in the DSCAM gene. The mRNA of DSCAM and the density of DSCAM in dendrites were significantly decreased in ASD compared to control iN cells. RNA sequencing analysis revealed that several synaptic function-related genes including NMDA receptor subunits were downregulated in ASD iN cells. Moreover, NMDA receptor (R)-mediated currents were significantly reduced in ASD compared to control iN cells. Normal NMDA-R-mediated current levels were rescued by expressing wild-type DSCAM in ASD iN cells, and reduced currents were observed by truncated DSCAM expression in control iN cells. shRNA-mediated DSCAM knockdown in control iN cells resulted in the downregulation of an NMDA-R subunit, which was rescued by the overexpression of shRNA-resistant DSCAM. Furthermore, DSCAM was co-localized with NMDA-R components in the dendritic spines of iN cells whereas their co-localizations were significantly reduced in ASD iN cells. Levels of phospho-ERK1/2 were significantly lower in ASD iN cells, suggesting a potential mechanism. A neural stem cell-specific Dscam heterozygous knockout mouse model, showing deficits in social interaction and social memory with reduced NMDA-R currents. These data suggest that DSCAM mutation causes pathological symptoms of ASD by dysregulating NMDA-R function.

Autism-associated miR-873 regulates ARID1B, SHANK3 and NRXN2 involved in neurodevelopment

Autism spectrum disorders (ASD) are highly heritable neurodevelopmental disorders with significant genetic heterogeneity. Noncoding microRNAs (miRNAs) are recognised as playing key roles in development of ASD albeit the function of these regulatory genes remains unclear. We previously conducted whole-exome sequencing of Australian families with ASD and identified four novel single nucleotide variations in mature miRNA sequences. A pull-down transcriptome analysis using transfected SH-SY5Y cells proposed a mechanistic model to examine changes in binding affinity associated with a unique mutation found in the conserved 'seed' region of miR-873-5p (rs777143952: T > A). Results suggested several ASD-risk genes were differentially targeted by wild-type and mutant miR-873 variants. In the current study, a dual-luciferase reporter assay confirmed miR-873 variants have a 20-30% inhibition/dysregulation effect on candidate autism risk genes ARID1B, SHANK3 and NRXN2 and also confirmed the affected expression with qPCR. In vitro mouse hippocampal neurons transfected with mutant miR-873 showed less morphological complexity and enhanced sodium currents and excitatory neurotransmission compared to cells transfected with wild-type miR-873. A second in vitro study showed CRISPR/Cas9 miR-873 disrupted SH-SY5Y neuroblastoma cells acquired a neuronal-like morphology and increased expression of ASD important genes ARID1B, SHANK3, ADNP2, ANK2 and CHD8. These results represent the first functional evidence that miR-873 regulates key neural genes involved in development and cell differentiation.

Identification of primary copy number variations reveal enrichment of Calcium, and MAPK pathways sensitizing secondary sites for autism

Background

Autism is a neurodevelopmental condition with genetic heterogeneity. It is characterized by difficulties in reciprocal social interactions with strong repetitive behaviors and stereotyped interests. Copy number variations (CNVs) are genomic structural variations altering the genomic structure either by duplication or deletion. De novo or inherited CNVs are found in 5–10% of autistic subjects with a size range of few kilobases to several megabases. CNVs predispose humans to various diseases by altering gene regulation, generation of chimeric genes, and disruption of the coding region or through position effect. Although, CNVs are not the initiating event in pathogenesis; additional preceding mutations might be essential for disease manifestation. The present study is aimed to identify the primary CNVs responsible for autism susceptibility in healthy cohorts to sensitize secondary-hits. In the current investigation, primary-hit autism gene CNVs are characterized in 1715 healthy cohorts of varying ethnicities across 12 populations using Affymetrix high-resolution array study. Thirty-eight individuals from twelve families residing in Karnataka, India, with the age group of 13–73 years are included for the comparative CNV analysis. The findings are validated against global 179 autism whole-exome sequence datasets derived from Simons Simplex Collection. These datasets are deposited at the Simons Foundation Autism Research Initiative (SFARI) database.

Results

The study revealed that 34.8% of the subjects carried 2% primary-hit CNV burden with 73 singleton-autism genes in different clusters. Of these, three conserved CNV breakpoints were identified with ARHGAP11B, DUSP22, and CHRNA7 as the target genes across 12 populations. Enrichment analysis of the population-specific autism genes revealed two signaling pathways—calcium and mitogen-activated protein kinases (MAPK) in the CNV identified regions. These impaired pathways affected the downstream cascades of neuronal function and physiology, leading to autism behavior. The pathway analysis of enriched genes unravelled complex protein interaction networks, which sensitized secondary sites for autism. Further, the identification of miRNA targets associated with autism gene CNVs added severity to the condition.

Conclusion

These findings contribute to an atlas of primary-hit genes to detect autism susceptibility in healthy cohorts, indicating their impact on secondary sites for manifestation.

ORVAL: a novel platform for the prediction and exploration of disease-causing oligogenic variant combinations

A tremendous amount of DNA sequencing data is being produced around the world with the ambition to capture in more detail the mechanisms underlying human diseases. While numerous bioinformatics tools exist that allow the discovery of causal variants in Mendelian diseases, little to no support is provided to do the same for variant combinations, an essential task for the discovery of the causes of oligogenic diseases. ORVAL (the Oligogenic Resource for Variant AnaLysis), which is presented here, provides an answer to this problem by focusing on generating networks of candidate pathogenic variant combinations in gene pairs, as opposed to isolated variants in unique genes. This online platform integrates innovative machine learning methods for combinatorial variant pathogenicity prediction with visualization techniques, offering several interactive and exploratory tools, such as pathogenic gene and protein interaction networks, a ranking of pathogenic gene pairs, as well as visual mappings of the cellular location and pathway information. ORVAL is the first web-based exploration platform dedicated to identifying networks of candidate pathogenic variant combinations with the sole ambition to help in uncovering oligogenic causes for patients that cannot rely on the classical disease analysis tools. ORVAL is available at https://orval.ibsquare.be.

Spatio-Temporal Roles of ASD-Associated Variants in Human Brain Development

Transcriptional regulation of the genome arguably provides the basis for the anatomical elaboration and dynamic operation of the human brain. It logically follows that genetic variations affecting gene transcription contribute to mental health disorders, including autism spectrum disorder (ASD). A number of recent studies have shown the role of de novo variants (DNVs) in disrupting early neurodevelopment. However, there is limited knowledge concerning the role of inherited variants during the early brain development of ASD. In this study, we investigate the role of rare inherited variations in neurodevelopment. We conducted co-expression network analyses using an anatomically comprehensive atlas of the developing human brain and examined whether rare coding and regulatory variants, identified from our genetic screening of Australian families with ASD, work in different spatio-temporal functions.

A Caenorhabditis elegans Model for Integrating the Functions of Neuropsychiatric Risk Genes Identifies Components Required for Normal Dendritic Morphology

Analysis of patient-derived DNA samples has identified hundreds of variants that are likely involved in neuropsychiatric diseases such as autism spectrum disorder (ASD) and schizophrenia (SCZ). While these studies couple behavioral phenotypes to individual genotypes, the number and diversity of candidate genes implicated in these disorders highlights the fact that the mechanistic underpinnings of these disorders are largely unknown. Here, we describe a RNAi-based screening platform that uses C. elegans to screen candidate neuropsychiatric risk genes (NRGs) for roles in controlling dendritic arborization. To benchmark this approach, we queried published lists of NRGs whose variants in ASD and SCZ are predicted to result in complete or partial loss of gene function. We found that a significant fraction (>16%) of these candidate NRGs are essential for dendritic development. Furthermore, these gene sets are enriched for dendritic arbor phenotypes (>14 fold) when compared to control RNAi datasets of over 500 human orthologs. The diversity of PVD structural abnormalities observed in these assays suggests that the functions of diverse NRGs (encoding transcription factors, chromatin remodelers, molecular chaperones and cytoskeleton-related proteins) converge to regulate neuronal morphology and that individual NRGs may play distinct roles in dendritic branching. We also demonstrate that the experimental value of this platform by providing additional insights into the molecular frameworks of candidate NRGs. Specifically, we show that ANK2/UNC-44 function is directly integrated with known regulators of dendritic arborization and suggest that altering the dosage of ARID1B/LET-526 expression during development affects neuronal morphology without diminishing aspects of cell fate specification.

Breakpoint Mapping of Symptomatic Balanced Translocations Links the EPHA6, KLF13 and UBR3 Genes to Novel Disease Phenotype

De novo balanced chromosomal aberrations (BCAs), such as reciprocal translocations and inversions, are genomic aberrations that, in approximately 25% of cases, affect the human phenotype. Delineation of the exact structure of BCAs may provide a precise diagnosis and/or point to new disease loci. We report on six patients with de novo balanced chromosomal translocations (BCTs) and one patient with a de novo inversion, in whom we mapped breakpoints to a resolution of 1 bp, using shallow whole-genome mate pair sequencing. In all seven cases, a disruption of at least one gene was found. In two patients, the phenotypic impact of the disrupted genes is well known (NFIA, ATP7A). In five patients, the aberration damaged genes: PARD3, EPHA6, KLF13, STK24, UBR3, MLLT10 and TLE3, whose influence on the human phenotype is poorly understood. In particular, our results suggest novel candidate genes for retinal degeneration with anophthalmia (EPHA6), developmental delay with speech impairment (KLF13), and developmental delay with brain dysembryoplastic neuroepithelial tumor (UBR3). In conclusion, identification of the exact structure of symptomatic BCTs using next generation sequencing is a viable method for both diagnosis and finding novel disease candidate genes in humans.

Mutations of MACF1, Encoding Microtubule-Actin Crosslinking-Factor 1, Cause Spectraplakinopathy

As a member of spectraplakin family of cytoskeletal crosslinking proteins, microtubule-actin crosslinking factor 1 (MACF1) controls cytoskeleton network dynamics. Knockout of Macf1 in mice resulted in the developmental retardation and embryonic lethality. Spectraplakinopathy type I, a novel neuromuscular condition characterized by periodic hypotonia, lax muscles, joint contracture, and diminished motor skill, was reported to be associated with heterozygous genomic duplication involving the MACF1 loci, with incomplete penetrance and highly variable clinical presentation in a single pedigree. In this study, parental-derived compound heterozygous novel missense mutations of MACF1, c.1517C>T (p.Thr506Ile) and c.11654T>C (p.Ile3885Thr), were found to co-segregate with disease status in two affected brothers presenting with progressive spastic tetraplegia, dystonia, joint contracture, feeding difficulty and developmental delay. We speculated that MACF1 mutations cause spectraplakinopathy inherited in an autosomal recessive manner. Our clinical findings expanded the phenotype of this neuromuscular disorder and provided new insights into the function of MACF1.

An integrative analysis of non-coding regulatory DNA variations associated with autism spectrum disorder

A number of genetic studies have identified rare protein-coding DNA variations associated with autism spectrum disorder (ASD), a neurodevelopmental disorder with significant genetic etiology and heterogeneity. In contrast, the contributions of functional, regulatory genetic variations that occur in the extensive non-protein-coding regions of the genome remain poorly understood. Here we developed a genome-wide analysis to identify the rare single nucleotide variants (SNVs) that occur in non-coding regions and determined the regulatory function and evolutionary conservation of these variants. Using publicly available datasets and computational predictions, we identified SNVs within putative regulatory regions in promoters, transcription factor binding sites, and microRNA genes and their target sites. Overall, we found that the regulatory variants in ASD cases were enriched in ASD-risk genes and genes involved in fetal neurodevelopment. As with previously reported coding mutations, we found an enrichment of the regulatory variants associated with dysregulation of neurodevelopmental and synaptic signaling pathways. Among these were several rare inherited SNVs found in the mature sequence of microRNAs predicted to affect the regulation of ASD-risk genes. We show a paternally inherited miR-873-5p variant with altered binding affinity for several risk-genes including NRXN2 and CNTNAP2 putatively overlay maternally inherited loss-of-function coding variations in NRXN1 and CNTNAP2 to likely increase the genetic liability in an idiopathic ASD case. Our analysis pipeline provides a new resource for identifying loss-of-function regulatory DNA variations that may contribute to the genetic etiology of complex disorders.

Assessment of patient clinical descriptions and pathogenic variants from gene panel sequences in the CAGI-5 intellectual disability challenge

The Critical Assessment of Genome Interpretation-5 intellectual disability challenge asked to use computational methods to predict patient clinical phenotypes and the causal variant(s) based on an analysis of their gene panel sequence data. Sequence data for 74 genes associated with intellectual disability (ID) and/or autism spectrum disorders (ASD) from a cohort of 150 patients with a range of neurodevelopmental manifestations (i.e. ID, autism, epilepsy, microcephaly, macrocephaly, hypotonia, ataxia) have been made available for this challenge. For each patient, predictors had to report the causative variants and which of the seven phenotypes were present. Since neurodevelopmental disorders are characterized by strong comorbidity, tested individuals often present more than one pathological condition. Considering the overall clinical manifestation of each patient, the correct phenotype has been predicted by at least one group for 93 individuals (62%). ID and ASD were the best predicted among the seven phenotypic traits. Also, causative or potentially pathogenic variants were predicted correctly by at least one group. However, the prediction of the correct causative variant seems to be insufficient to predict the correct phenotype. In some cases, the correct prediction has been supported by rare or common variants in genes different from the causative one.

Family-based exome sequencing and case-control analysis implicate CEP41 as an ASD gene

Autism Spectrum Disorder (ASD) is a complex neurodevelopmental disorder with a strong genetic component. Although next-generation sequencing (NGS) technologies have been successfully applied to gene identification in de novo ASD, the genetic architecture of familial ASD remains largely unexplored. Our approach, which leverages the high specificity and sensitivity of NGS technology, has focused on rare variants in familial autism. We used NGS exome sequencing in 26 families with distantly related affected individuals to identify genes with private gene disrupting and missense variants of interest (VOI). We found that the genes carrying VOIs were enriched for biological processes related to cell projection organization and neuron development, which is consistent with the neurodevelopmental hypothesis of ASD. For a subset of genes carrying VOIs, we then used targeted NGS sequencing and gene-based variant burden case-control analysis to test for association with ASD. Missense variants in one gene, CEP41, associated significantly with ASD (p = 6.185^e-05). Homozygous gene-disrupting variants in CEP41 were initially found to be responsible for recessive Joubert syndrome. Using a zebrafish model, we evaluated the mechanism by which the CEP41 variants might contribute to ASD. We found that CEP41 missense variants affect development of the axonal tract, cranial neural crest migration and social behavior phenotype. Our work demonstrates the involvement of CEP41 heterozygous missense variants in ASD and that biological processes involved in cell projection organization and neuron development are enriched in ASD families we have studied.

MACF1 Controls Migration and Positioning of Cortical GABAergic Interneurons in Mice

GABAergic interneurons develop in the ganglionic eminence in the ventral telencephalon and tangentially migrate into the cortical plate during development. However, key molecules controlling interneuron migration remain poorly identified. Here, we show that microtubule-actin cross-linking factor 1 (MACF1) regulates GABAergic interneuron migration and positioning in the developing mouse brain. To investigate the role of MACF1 in developing interneurons, we conditionally deleted the MACF1 gene in mouse interneuron progenitors and their progeny using Dlx5/6-Cre-IRES-EGFP and Nkx2.1-Cre drivers. We found that MACF1 deletion results in a marked reduction and defective positioning of interneurons in the mouse cerebral cortex and hippocampus, suggesting abnormal interneuron migration. Indeed, the speed and mode of interneuron migration were abnormal in the MACF1-mutant brain, compared with controls. Additionally, MACF1-deleted interneurons showed a significant reduction in the length of their leading processes and dendrites in the mouse brain. Finally, loss of MACF1 decreased microtubule stability in cortical interneurons. Our findings suggest that MACF1 plays a critical role in cortical interneuron migration and positioning in the developing mouse brain.

Understanding Neurodevelopmental Disorders: The Promise of Regulatory Variation in the 3'UTRome

Neurodevelopmental disorders have a strong genetic component, but despite widespread efforts, the specific genetic factors underlying these disorders remain undefined for a large proportion of affected individuals. Given the accessibility of exome sequencing, this problem has thus far been addressed from a protein-centric standpoint; however, protein-coding regions only make up ∼1% to 2% of the human genome. With the advent of whole genome sequencing we are in the midst of a paradigm shift as it is now possible to interrogate the entire sequence of the human genome (coding and noncoding) to fill in the missing heritability of complex disorders. These new technologies bring new challenges, as the number of noncoding variants identified per individual can be overwhelming, making it prudent to focus on noncoding regions of known function, for which the effects of variation can be predicted and directly tested to assess pathogenicity. The 3'UTRome is a region of the noncoding genome that perfectly fulfills these criteria and is of high interest when searching for pathogenic variation related to complex neurodevelopmental disorders. Herein, we review the regulatory roles of the 3'UTRome as binding sites for microRNAs or RNA binding proteins, or during alternative polyadenylation. We detail existing evidence that these regions contribute to neurodevelopmental disorders and outline strategies for identification and validation of novel putatively pathogenic variation in these regions. This evidence suggests that studying the 3'UTRome will lead to the identification of new risk factors, new candidate disease genes, and a better understanding of the molecular mechanisms contributing to neurodevelopmental disorders.

Microarray analysis of gene expression in the cyclooxygenase knockout mice - a connection to autism spectrum disorder

The cellular and molecular events that take place during brain development play an important role in governing function of the mature brain. Lipid-signalling molecules such as prostaglandin E₂ (PGE₂ ) play an important role in healthy brain development. Abnormalities along the COX-PGE₂ signalling pathway due to genetic or environmental causes have been linked to autism spectrum disorder (ASD). This study aims to evaluate the effect of altered COX-PGE₂ signalling on development and function of the prenatal brain using male mice lacking cyclooxygenase-1 and cyclooxygenase-2 (COX-1^-/- and COX-2^-/- ) as potential model systems of ASD. Microarray analysis was used to determine global changes in gene expression during embryonic days 16 (E16) and 19 (E19). Gene Ontology: Biological Process (GO:BP) and Kyoto Encyclopedia of Genes and Genomes (KEGG) were implemented to identify affected developmental genes and cellular processes. We found that in both knockouts the brain at E16 had nearly twice as many differentially expressed genes, and affected biological pathways containing various ASD-associated genes important in neuronal function. Interestingly, using GeneMANIA and Cytoscape we also show that the ASD-risk genes identified in both COX-1^-/- and COX-2^-/- models belong to protein-interaction networks important for brain development despite of different cellular localization of these enzymes. Lastly, we identified eight genes that belong to the Wnt signalling pathways exclusively in the COX-2^-/- mice at E16. The level of PKA-phosphorylated β-catenin (S552), a major activator of the Wnt pathway, was increased in this model, suggesting crosstalk between the COX-2-PGE₂ and Wnt pathways during early brain development. Overall, these results provide further molecular insight into the contribution of the COX-PGE₂ pathways to ASD and demonstrate that COX-1^-/- and COX-2^-/- animals might be suitable new model systems for studying the disorders.

Family studies to find rare high risk variants in migraine

INTRODUCTION

Migraine has long been known as a common complex disease caused by genetic and environmental factors. The pathophysiology and the specific genetic susceptibility are poorly understood. Common variants only explain a small part of the heritability of migraine. It is thought that rare genetic variants with bigger effect size may be involved in the disease. Since migraine has a tendency to cluster in families, a family approach might be the way to find these variants. This is also indicated by identification of migraine-associated loci in classical linkage-analyses in migraine families. A single migraine study using a candidate-gene approach was performed in 2010 identifying a rare mutation in the TRESK potassium channel segregating in a large family with migraine with aura, but this finding has later become questioned. The technologies of next-generation sequencing (NGS) now provides an affordable tool to investigate the genetic variation in the entire exome or genome. The family-based study design using NGS is described in this paper. We also review family studies using NGS that have been successful in finding rare variants in other common complex diseases in order to argue the promising application of a family approach to migraine.

METHOD

PubMed was searched to find studies that looked for rare genetic variants in common complex diseases through a family-based design using NGS, excluding studies looking for de-novo mutations, or using a candidate-gene approach and studies on cancer. All issues from Nature Genetics and PLOS genetics 2014, 2015 and 2016 (UTAI June) were screened for relevant papers. Reference lists from included and other relevant papers were also searched. For the description of the family-based study design using NGS an in-house protocol was used.

RESULTS

Thirty-two successful studies, which covered 16 different common complex diseases, were included in this paper. We also found a single migraine study. Twenty-three studies found one or a few family specific variants (less than five), while other studies found several possible variants. Not all of them were genome wide significant. Four studies performed follow-up analyses in unrelated cases and controls and calculated odds ratios that supported an association between detected variants and risk of disease. Studies of 11 diseases identified rare variants that segregated fully or to a large degree with the disease in the pedigrees.

CONCLUSION

It is possible to find rare high risk variants for common complex diseases through a family-based approach. One study using a family approach and NGS to find rare variants in migraine has already been published but with strong limitations. More studies are under way.

Delineating SPTAN1 associated phenotypes: from isolated epilepsy to encephalopathy with progressive brain atrophy

De novo in-frame deletions and duplications in the SPTAN1 gene, encoding the non-erythrocyte αII spectrin, have been associated with severe West syndrome with hypomyelination and pontocerebellar atrophy. We aimed at comprehensively delineating the phenotypic spectrum associated with SPTAN1 mutations. Using different molecular genetic techniques, we identified 20 patients with a pathogenic or likely pathogenic SPTAN1 variant and reviewed their clinical, genetic and imaging data. SPTAN1 de novo alterations included seven unique missense variants and nine in-frame deletions/duplications of which 12 were novel. The recurrent three-amino acid duplication p.(Asp2303_Leu2305dup) occurred in five patients. Our patient cohort exhibited a broad spectrum of neurodevelopmental phenotypes, comprising six patients with mild to moderate intellectual disability, with or without epilepsy and behavioural disorders, and 14 patients with infantile epileptic encephalopathy, of which 13 had severe neurodevelopmental impairment and four died in early childhood. Imaging studies suggested that the severity of neurological impairment and epilepsy correlates with that of structural abnormalities as well as the mutation type and location. Out of seven patients harbouring mutations outside the α/β spectrin heterodimerization domain, four had normal brain imaging and three exhibited moderately progressive brain and/or cerebellar atrophy. Twelve of 13 patients with mutations located within the spectrin heterodimer contact site exhibited severe and progressive brain, brainstem and cerebellar atrophy, with hypomyelination in most. We used fibroblasts from five patients to study spectrin aggregate formation by Triton-X extraction and immunocytochemistry followed by fluorescence microscopy. αII/βII aggregates and αII spectrin in the insoluble protein fraction were observed in fibroblasts derived from patients with the mutations p.(Glu2207del), p.(Asp2303_Leu2305dup) and p.(Arg2308_Met2309dup), all falling in the nucleation site of the α/β spectrin heterodimer region. Molecular modelling of the seven SPTAN1 amino acid changes provided preliminary evidence for structural alterations of the A-, B- and/or C-helices within each of the mutated spectrin repeats. We conclude that SPTAN1-related disorders comprise a wide spectrum of neurodevelopmental phenotypes ranging from mild to severe and progressive. Spectrin aggregate formation in fibroblasts with mutations in the α/β heterodimerization domain seems to be associated with a severe neurodegenerative course and suggests that the amino acid stretch from Asp2303 to Met2309 in the α20 repeat is important for α/β spectrin heterodimer formation and/or αII spectrin function.

Hypermasculinised facial morphology in boys and girls with Autism Spectrum Disorder and its association with symptomatology

Elevated prenatal testosterone exposure has been associated with Autism Spectrum Disorder (ASD) and facial masculinity. By employing three-dimensional (3D) photogrammetry, the current study investigated whether prepubescent boys and girls with ASD present increased facial masculinity compared to typically-developing controls. There were two phases to this research. 3D facial images were obtained from a normative sample of 48 boys and 53 girls (3.01-12.44 years old) to determine typical facial masculinity/femininity. The sexually dimorphic features were used to create a continuous 'gender score', indexing degree of facial masculinity. Gender scores based on 3D facial images were then compared for 54 autistic and 54 control boys (3.01-12.52 years old), and also for 20 autistic and 60 control girls (4.24-11.78 years). For each sex, increased facial masculinity was observed in the ASD group relative to control group. Further analyses revealed that increased facial masculinity in the ASD group correlated with more social-communication difficulties based on the Social Affect score derived from the Autism Diagnostic Observation Scale-Generic (ADOS-G). There was no association between facial masculinity and the derived Restricted and Repetitive Behaviours score. This is the first study demonstrating facial hypermasculinisation in ASD and its relationship to social-communication difficulties in prepubescent children.

Whole genome sequencing resource identifies 18 new candidate genes for autism spectrum disorder

We are performing whole-genome sequencing of families with autism spectrum disorder (ASD) to build a resource (MSSNG) for subcategorizing the phenotypes and underlying genetic factors involved. Here we report sequencing of 5,205 samples from families with ASD, accompanied by clinical information, creating a database accessible on a cloud platform and through a controlled-access internet portal. We found an average of 73.8 de novo single nucleotide variants and 12.6 de novo insertions and deletions or copy number variations per ASD subject. We identified 18 new candidate ASD-risk genes and found that participants bearing mutations in susceptibility genes had significantly lower adaptive ability (P = 6 × 10^-4). In 294 of 2,620 (11.2%) of ASD cases, a molecular basis could be determined and 7.2% of these carried copy number variations and/or chromosomal abnormalities, emphasizing the importance of detecting all forms of genetic variation as diagnostic and therapeutic targets in ASD.

Anatomy and Cell Biology of Autism Spectrum Disorder: Lessons from Human Genetics

Until recently autism spectrum disorder (ASD) was regarded as a neurodevelopmental condition with unknown causes and pathogenesis. In the footsteps of the revolution of genome technologies and genetics, and with its high degree of heritability, ASD became the first neuropsychiatric disorder for which clues towards molecular and cellular pathogenesis were uncovered by genetic identification of susceptibility genes. Currently several hundreds of risk genes have been assigned, with a recurrence below 1% in the ASD population. The multitude and diversity of known ASD genes has extended the clinical notion that ASD comprises very heterogeneous conditions ranging from severe intellectual disabilities to mild high-functioning forms. The results of genetics have allowed to pinpoint a limited number of cellular and molecular processes likely involved in ASD including protein synthesis, signal transduction, transcription/chromatin remodelling and synaptic function all playing an essential role in the regulation of synaptic homeostasis during brain development. In this context, we highlight the role of protein synthesis as a key process in ASD pathogenesis as it might be central in synaptic deregulation and a potential target for intervention. These current insights should lead to a rational design of interventions in molecular and cellular pathways of ASD pathogenesis that may be applied to affected individuals in the future.

A Subset of Autism-Associated Genes Regulate the Structural Stability of Neurons

Autism spectrum disorder (ASD) comprises a range of neurological conditions that affect individuals’ ability to communicate and interact with others. People with ASD often exhibit marked qualitative difficulties in social interaction, communication, and behavior. Alterations in neurite arborization and dendritic spine morphology, including size, shape, and number, are hallmarks of almost all neurological conditions, including ASD. As experimental evidence emerges in recent years, it becomes clear that although there is broad heterogeneity of identified autism risk genes, many of them converge into similar cellular pathways, including those regulating neurite outgrowth, synapse formation and spine stability, and synaptic plasticity. These mechanisms together regulate the structural stability of neurons and are vulnerable targets in ASD. In this review, we discuss the current understanding of those autism risk genes that affect the structural connectivity of neurons. We sub-categorize them into (1) cytoskeletal regulators, e.g., motors and small RhoGTPase regulators; (2) adhesion molecules, e.g., cadherins, NCAM, and neurexin superfamily; (3) cell surface receptors, e.g., glutamatergic receptors and receptor tyrosine kinases; (4) signaling molecules, e.g., protein kinases and phosphatases; and (5) synaptic proteins, e.g., vesicle and scaffolding proteins. Although the roles of some of these genes in maintaining neuronal structural stability are well studied, how mutations contribute to the autism phenotype is still largely unknown. Investigating whether and how the neuronal structure and function are affected when these genes are mutated will provide insights toward developing effective interventions aimed at improving the lives of people with autism and their families.

Recent Advances in Autism Spectrum Disorders: Applications of Whole Exome Sequencing Technology

Autism spectrum disorders (ASD) is characterized by three core symptoms with impaired reciprocal social interaction and communication, a pattern of repetitive behavior and/or restricted interests in early childhood. The prevalence is higher in male children than in female children. As a complex neurodevelopmental disorder, the phenotype and severity of autism are extremely heterogeneous with differences from one patient to another. Genetics has a key role in the etiology of autism. Environmental factors are also interacting with the genetic profile and cause abnormal changes in neuronal development, brain growth, and functional connectivity. The term of exome represents less than 1% of the human genome, but contains 85% of known disease-causing variants. Whole-exome sequencing (WES) is an application of the next generation sequencing technology to determine the variations of all coding regions, or exons of known genes. For this reason, WES has been extensively used for clinical studies in the recent years. WES has achieved great success in the past years for identifying Mendelian disease genes. This review evaluates the potential of current findings in ASD for application in next generation sequencing technology, particularly WES. WES and whole-genome sequencing (WGS) approaches may lead to the discovery of underlying genetic factors for ASD and may thereby identify novel therapeutic targets for this disorder.

Variation in the Profile of Anxiety Disorders in Boys with an ASD According to Method and Source of Assessment

To determine any variation that might occur due to the type of assessment and source used to assess them, the prevalence of 7 anxiety disorders were investigated in a sample of 140 boys with an Autism spectrum disorder (ASD) and 50 non-ASD (NASD) boys via the Child and Adolescent Symptom Inventory and the KIDSCID Clinical Interview. Boys with an ASD were significantly more anxious than their NASD peers. Data collected from the boys with an ASD themselves showed differences in the severity and diagnostic criterion of anxiety disorders to data collected from the boys' parents. There were age-related variations to the pattern of anxiety disorder differences across reports from the boys with an ASD and reports from their parents.

Replicate exome-sequencing in a multiple-generation family: improved interpretation of next-generation sequencing data

Background

Whole-exome sequencing (WES) is rapidly evolving into a tool of choice for rapid, and inexpensive identification of molecular genetic lesions within targeted regions of the human genome. While biases in WES coverage of nucleotides in targeted regions are recognized, it is not well understood how repetition of WES improves the interpretation of sequencing results in a clinical diagnostic setting.

Method

To address this, we compared independently generated exome-capture of six individuals from three-generations sequenced in triplicate. This generated between 48x-86x mean target depth of high-quality mapped bases (>Q20) for each technical replicate library. Cumulatively, we achieved 179 - 208x average target coverage for each individual in the pedigree. Using this experimental design, we evaluated stochastics in WES interpretation, genotyping sensitivity, and accuracy to detect de novo variants.

Results

In this study, we show that repetition of WES improved the interpretation of the capture target regions after aggregating the data (93.5 - 93.9 %). Compared to 81.2 - 89.6 % (50.2-55.4 Mb of 61.7 M) coverage of targeted bases at ≥20x in the individual technical replicates, the aggregated data covered 93.5 - 93.9 % of targeted bases (57.7 – 58.0 of 61.7 M) at ≥20x threshold, suggesting a 4.3 – 12.7 % improvement in coverage. Each individual’s aggregate dataset recovered 3.4 – 6.4 million bases within variable targeted regions. We uncovered technical variability (2-5 %) inherent to WES technique. We also show improved interpretation in assessing clinically important regions that lack interpretation under current conditions, affecting 12–16 of the 56 genes recommended for secondary analysis by American College of Medical Genetics (ACMG). We demonstrate that comparing technical replicate WES datasets and their derived aggregate data can effectively address overall WES genotyping discrepancies.

Conclusion

We describe a method to evaluate the reproducibility and stochastics in exome library preparation, and delineate the advantages of aggregating the data derived from technical replicates. The implications of this study are directly applicable to improved experimental design and provide an opportunity to rapidly, efficiently, and accurately arrive at reliable candidate nucleotide variants.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2107-y) contains supplementary material, which is available to authorized users.

A neuroligin-3 mutation implicated in autism causes abnormal aggression and increases repetitive behavior in mice

BACKGROUND

Aggression is common in patients with autism spectrum disorders (ASD) along with the core symptoms of impairments in social communication and repetitive behavior. Risperidone, an atypical antipsychotic, is widely used to treat aggression in ASD. In order to understand the neurobiological underpinnings of these challenging behaviors, a thorough characterisation of behavioral endophenotypes in animal models is required.

METHODS

We investigated aggression in mice containing the ASD-associated R451C (arginine to cysteine residue 451 substitution) mutation in neuroligin-3 (NL3). Furthermore, we sought to verify social interaction impairments and assess olfaction, anxiety, and repetitive and restrictive behavior in NL3(R451C) mutant mice.

RESULTS

We show a pronounced elevation in aggressive behavior in NL3(R451C) mutant mice. Treatment with risperidone reduced this aggression to wild-type (WT) levels. Juvenile and adult social interactions were also investigated, and subtle differences in initiation of interaction were seen in juvenile NL3(R451C) mice. No genotype differences in olfactory discrimination or anxiety were observed indicating that aggression was not dependent on altered olfaction, stress response, or social preference. We also describe repetitive behavior in NL3(R451C) mice as assessed by a clinically relevant object exploration task.

CONCLUSIONS

The presence of aberrant aggression and other behavioral phenotypes in NL3(R451C) mice consistent with clinical traits strengthen face validity of this model of ASD. Furthermore, we demonstrate predictive validity in this model through the reversal of the aggressive phenotype with risperidone. This is the first demonstration that risperidone can ameliorate aggression in an animal model of ASD and will inform mechanistic and therapeutic research into the neurobiology underlying abnormal behaviors in ASD.

Understanding the molecular basis of autism in a dish using hiPSCs-derived neurons from ASD patients

Autism spectrum disorder (ASD) is a complex neurodevelopmental disorder characterized by deficits in social cognition, language development, and repetitive/restricted behaviors. Due to the complexity and heterogeneity of ASD and lack of a proper human cellular model system, the pathophysiological mechanism of ASD during the developmental process is largely unknown. However, recent progress in induced pluripotent stem cell (iPSC) technology as well as in vitro neural differentiation techniques have allowed us to functionally characterize neurons and analyze cortical development during neural differentiation. These technical advances will increase our understanding of the pathogenic mechanisms of heterogeneous ASD and help identify molecular biomarkers for patient stratification as well as personalized medicine. In this review, we summarize our current knowledge of iPSC generation, differentiation of specific neuronal subtypes from iPSCs, and phenotypic characterizations of human ASD patient-derived iPSC models. Finally, we discuss the current limitations of iPSC technology and future directions of ASD pathophysiology studies using iPSCs.

Neurexin-1 regulates sleep and synaptic plasticity in Drosophila melanogaster

Neurexins are cell adhesion molecules that are important for synaptic plasticity and homeostasis, although links to sleep have not yet been investigated. We examined the effects of neurexin-1 perturbation on sleep in Drosophila, showing that neurexin-1 nulls displayed fragmented sleep and altered circadian rhythm. Conversely, the over-expression of neurexin-1 could increase and consolidate night-time sleep. This was not solely due to developmental effects as it could be induced acutely in adulthood, and was coupled with evidence of synaptic growth. The timing of over-expression could differentially impact sleep patterns, with specific night-time effects. These results show that neurexin-1 was dynamically involved in synaptic plasticity and sleep in Drosophila. Neurexin-1 and a number of its binding partners have been repeatedly associated with mental health disorders, including autism spectrum disorders, schizophrenia and Tourette syndrome, all of which are also linked to altered sleep patterns. How and when plasticity-related proteins such as neurexin-1 function during sleep can provide vital information on the interaction between synaptic homeostasis and sleep, paving the way for more informed treatments of human disorders.

Neuroligin 1 modulates striatal glutamatergic neurotransmission in a pathway and NMDAR subunit-specific manner

Together with its presynaptic partner Neurexin 1 (Nxn1), Neuroligin 1 (NL1) participates in synapse specification and synapse maintenance. We and others have shown that NL1 can also modulate glutamatergic synaptic function in the central nervous system of rodent models. These molecular/cellular changes can translate into altered animal behaviors that are thought to be analogous to symptomatology of neuropsychiatric disorders. For example, in dorsal striatum of NL1 deletion mice, we previously reported that the ratio N-methyl-D-aspartate receptor (NMDAR) mediated synaptic currents to α-amino-3-hydroxyl-5-methyl-4-isoxazole-propionate receptor (AMPAR) mediated synaptic currents (NMDA/AMPA) is reduced in medium spiny neuron (MSNs). Importantly, this reduction in NMDA/AMPA ratio correlated with increased repetitive grooming. The striatum is the input nucleus of the basal ganglia (BG). Classical models of this circuitry imply that there are two principal pathways that render distinct and somewhat opposite striatal outputs critical to the function of these nuclei in modulating motor behavior. Thus, we set out to better characterize the effects of NL1 deletion on direct and indirect pathways of the dorsal striatum by genetically labeling MSNs participating in the direct and indirect pathways. We demonstrate that a decrease in NMDAR-mediated currents is limited to MSNs of the direct pathway. Furthermore, the decrease in NMDAR-mediated currents is largely due to a reduction in function of NMDARs containing the GluN2A subunit. In contrast, indirect pathway MSNs in NL1 knockout (KO) mice showed a reduction in the frequency of miniature excitatory neurotransmission not observed in the direct pathway. Thus, NL1 deletion differentially affects direct and indirect pathway MSNs in dorsal striatum. These findings have potential implications for striatal function in NL1 KO mice.

The Role of Epigenetic Change in Autism Spectrum Disorders

Autism spectrum disorders (ASD) are a heterogeneous group of neurodevelopmental disorders characterized by problems with social communication, social interaction, and repetitive or restricted behaviors. ASD are comorbid with other disorders including attention deficit hyperactivity disorder, epilepsy, Rett syndrome, and Fragile X syndrome. Neither the genetic nor the environmental components have been characterized well enough to aid diagnosis or treatment of non-syndromic ASD. However, genome-wide association studies have amassed evidence suggesting involvement of hundreds of genes and a variety of associated genetic pathways. Recently, investigators have turned to epigenetics, a prime mediator of environmental effects on genomes and phenotype, to characterize changes in ASD that constitute a molecular level on top of DNA sequence. Though in their infancy, such studies have the potential to increase our understanding of the etiology of ASD and may assist in the development of biomarkers for its prediction, diagnosis, prognosis, and eventually in its prevention and intervention. This review focuses on the first few epigenome-wide association studies of ASD and discusses future directions.

CACNA1D de novo mutations in autism spectrum disorders activate Cav1.3 L-type calcium channels

BACKGROUND

Cav1.3 voltage-gated L-type calcium channels (LTCCs) are part of postsynaptic neuronal signaling networks. They play a key role in brain function, including fear memory and emotional and drug-taking behaviors. A whole-exome sequencing study identified a de novo mutation, p.A749G, in Cav1.3 α1-subunits (CACNA1D), the second main LTCC in the brain, as 1 of 62 high risk-conferring mutations in a cohort of patients with autism and intellectual disability. We screened all published genetic information available from whole-exome sequencing studies and identified a second de novo CACNA1D mutation, p.G407R. Both mutations are present only in the probands and not in their unaffected parents or siblings.

METHODS

We functionally expressed both mutations in tsA-201 cells to study their functional consequences using whole-cell patch-clamp.

RESULTS

The mutations p.A749G and p.G407R caused dramatic changes in channel gating by shifting (~15 mV) the voltage dependence for steady-state activation and inactivation to more negative voltages (p.A749G) or by pronounced slowing of current inactivation during depolarizing stimuli (p.G407R). In both cases, these changes are compatible with a gain-of-function phenotype.

CONCLUSIONS

Our data, together with the discovery that Cav1.3 gain-of-function causes primary aldosteronism with seizures, neurologic abnormalities, and intellectual disability, suggest that Cav1.3 gain-of-function mutations confer a major part of the risk for autism in the two probands and may even cause the disease. Our findings have immediate clinical relevance because blockers of LTCCs are available for therapeutic attempts in affected individuals. Patients should also be explored for other symptoms likely resulting from Cav1.3 hyperactivity, in particular, primary aldosteronism.

Contribution of common and rare variants of the PTCHD1 gene to autism spectrum disorders and intellectual disability

Recent findings revealed rare copy number variants and missense changes in the X-linked gene PTCHD1 in autism spectrum disorder (ASD) and intellectual disability (ID). Here, we aim to explore the contribution of common PTCHD1 variants in ASD and gain additional evidence for the role of rare variants of this gene in ASD and ID. A two-stage case-control association study investigated 28 tag single nucleotide polymorphisms (SNPs) in 994 ASD cases and 1035 controls from four European populations. Mutation screening was performed in 673 individuals who included 240 ASD cases, 183 ID patients and 250 controls. The case-control association study showed a significant association with rs7052177 (P=6.13E-4) in the ASD discovery sample that was replicated in an independent sample (P=0.03). A Mantel-Haenszel meta-analysis for rs7052177T considering the four European populations showed an odds ratio of 0.58 (P=7E-05). This SNP is predicted to be located in a transcription factor binding site. No rare missense PTCHD1 variants were found in our ASD cohort and only one was identified in the ID sample. A duplication (27 bp) in the promoter region, absent from 590 controls, was found in three ASD patients (Fisher exact test, P=0.024). A gene reporter assay showed a significant decrease in the transcriptional activity (26%) driven by this variant. Moreover, we found that the longest allele of a trinucleotide repeat located upstream from PTCHD1 was associated with ASD (P=0.003, permP=0.0186). Our results further support the involvement of PTCHD1 in ASD, suggesting that both common and rare variants contribute to the disorder.

Integrated systems analysis reveals a molecular network underlying autism spectrum disorders

Autism is a complex disease whose etiology remains elusive. We integrated previously and newly generated data and developed a systems framework involving the interactome, gene expression and genome sequencing to identify a protein interaction module with members strongly enriched for autism candidate genes. Sequencing of 25 patients confirmed the involvement of this module in autism, which was subsequently validated using an independent cohort of over 500 patients. Expression of this module was dichotomized with a ubiquitously expressed subcomponent and another subcomponent preferentially expressed in the corpus callosum, which was significantly affected by our identified mutations in the network center. RNA-sequencing of the corpus callosum from patients with autism exhibited extensive gene mis-expression in this module, and our immunochemical analysis showed that the human corpus callosum is predominantly populated by oligodendrocyte cells. Analysis of functional genomic data further revealed a significant involvement of this module in the development of oligodendrocyte cells in mouse brain. Our analysis delineates a natural network involved in autism, helps uncover novel candidate genes for this disease and improves our understanding of its molecular pathology.