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

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

Reliability of human cortical organoid generation.

The differentiation of pluripotent stem cells in three-dimensional cultures can recapitulate key aspects of brain development, but protocols are prone to variable results. Here we differentiated multiple human pluripotent stem cell lines for over 100 d using our previously developed approach to generate brain-region-specific organoids called cortical spheroids and, using several assays, found that spheroid generation was highly reliable and consistent. We anticipate the use of this approach for large-scale differentiation experiments and disease modeling.

Transcriptomic networks implicate neuronal energetic abnormalities in three mouse models harboring autism and schizophrenia-associated mutations.

Genetic risk for psychiatric illness is complex, so identification of shared molecular pathways where distinct forms of genetic risk might coincide is of substantial interest. A growing body of genetic and genomic studies suggest that such shared molecular pathways exist across disorders with different clinical presentations, such as schizophrenia and autism spectrum disorder (ASD). But how this relates to specific genetic risk factors is unknown. Further, whether some of the molecular changes identified in brain relate to potentially confounding antemortem or postmortem factors are difficult to prove. We analyzed the transcriptome from the cortex and hippocampus of three mouse lines modeling human copy number variants (CNVs) associated with schizophrenia and ASD: Df(h15q13)/+, Df(h22q11)/+, and Df(h1q21)/+ which carry the 15q13.3 deletion, 22q11.2 deletion, and 1q21.1 deletion, respectively. Although we found very little overlap of differential expression at the level of individual genes, gene network analysis identified two cortical and two hippocampal modules of co-expressed genes that were dysregulated across all three mouse models. One cortical module was associated with neuronal energetics and firing rate, and overlapped with changes identified in postmortem human brain from SCZ and ASD patients. These data highlight aspects of convergent gene expression in mouse models harboring major risk alleles, and strengthen the connection between changes in neuronal energetics and neuropsychiatric disorders in humans.

sRNAnalyzer-a flexible and customizable small RNA sequencing data analysis pipeline.

Although many tools have been developed to analyze small RNA sequencing (sRNA-Seq) data, it remains challenging to accurately analyze the small RNA population, mainly due to multiple sequence ID assignment caused by short read length. Additional issues in small RNA analysis include low consistency of microRNA (miRNA) measurement results across different platforms, miRNA mapping associated with miRNA sequence variation (isomiR) and RNA editing, and the origin of those unmapped reads after screening against all endogenous reference sequence databases. To address these issues, we built a comprehensive and customizable sRNA-Seq data analysis pipeline-sRNAnalyzer, which enables: (i) comprehensive miRNA profiling strategies to better handle isomiRs and summarization based on each nucleotide position to detect potential SNPs in miRNAs, (ii) different sequence mapping result assignment approaches to simulate results from microarray/qRT-PCR platforms and a local probabilistic model to assign mapping results to the most-likely IDs, (iii) comprehensive ribosomal RNA filtering for accurate mapping of exogenous RNAs and summarization based on taxonomy annotation. We evaluated our pipeline on both artificial samples (including synthetic miRNA and Escherichia coli cultures) and biological samples (human tissue and plasma). sRNAnalyzer is implemented in Perl and available at: http://srnanalyzer.systemsbiology.net/.

Absolute binding free energies for the SAMPL6 cucurbit[8]uril host-guest challenge via the AMOEBA polarizable force field.

As part of the SAMPL6 host-guest blind challenge, the AMOEBA force field was applied to calculate the absolute binding free energy for a cucurbit[8]uril host complexed with 14 diverse guests, ranging from small, rigid structures to drug molecules. The AMOEBA results from the initial submission prompted an investigation into aspects of the methodology and parameterization employed. Lessons learned from the blind challenge include: a double annihilation scheme (electrostatics and van der Waals) is needed to obtain proper sampling of guest conformations, annihilation of key torsion parameters of the guest are recommended for flexible guests, and a more thorough analysis of torsion parameters is warranted. When put in to practice with the AMOEBA model, the lessons learned improved the MUE from 2.63 to 1.20 kcal/mol and the RMSE from 3.62 to 1.68 kcal/mol, respectively. Overall, the AMOEBA protocol for determining absolute binding free energies benefitted from participation in the SAMPL6 host-guest blind challenge and the results suggest the implementation of the methodology in future host-guest calculations.

Expression of Cntnap2 (Caspr2) in multiple levels of sensory systems.

Genome-wide association studies and copy number variation analyses have linked contactin associated protein 2 (Caspr2, gene name Cntnap2) with autism spectrum disorder (ASD). In line with these findings, mice lacking Caspr2 (Cntnap2(-/-)) were shown to have core autism-like deficits including abnormal social behavior and communication, and behavior inflexibility. However the role of Caspr2 in ASD pathogenicity remains unclear. Here we have generated a new Caspr2:tau-LacZ knock-in reporter line (Cntnap2(tlacz/tlacz)), which enabled us to monitor the neuronal circuits in the brain expressing Caspr2. We show that Caspr2 is expressed in many brain regions and produced a comprehensive report of Caspr2 expression. Moreover, we found that Caspr2 marks all sensory modalities: it is expressed in distinct brain regions involved in different sensory processings and is present in all primary sensory organs. Olfaction-based behavioral tests revealed that mice lacking Caspr2 exhibit abnormal response to sensory stimuli and lack preference for novel odors. These results suggest that loss of Caspr2 throughout the sensory system may contribute to the sensory manifestations frequently observed in ASD.

Interaction proteomics of canonical Caspr2 (CNTNAP2) reveals the presence of two Caspr2 isoforms with overlapping interactomes.

Autism is a human developmental brain disorder characterized by impaired social interaction and communication. Contactin-associated protein-like 2 (Caspr2, CNTNAP2) is a known genetic risk factor of autism. However, how this protein might contribute to pathology is unclear. In this study, we demonstrate that Caspr2 is abundantly present in lipid raft and in the synaptic membrane but is highly depleted in the postsynaptic density. The Caspr2 protein level in hippocampus is present at a constant level during synapse formation and myelination from P0 to P84. Interaction proteomics revealed the interactors of Caspr2, including CNTN2, KCNAs, members of the ADAM family (ADAM22, ADAM23 and ADAM11), members of LGI family and MAGUKs (DLGs and MPPs). Interestingly, a short form of Caspr2 was detected, which lacks most of the extracellular domains, however, is still associated with ADAM22 and to a lesser extent LGI1 and Kv1 channels. The comprehensive Caspr2 interactome revealed here might aid in understanding the molecular mechanisms underlying autism. This article is part of a Special Issue titled Neuroproteomics: Applications in Neuroscience and Neurology.

Specific inhibition of secreted NRG1 types I-II by heparin enhances Schwann Cell myelination.

Primary cultures of mixed neuron and Schwann cells prepared from dorsal root ganglia (DRG) are extensively used as a model to study myelination. These dissociated DRG cultures have the particular advantage of bypassing the difficulty in purifying mouse Schwann cells, which is often required when using mutant mice. However, the drawback of this experimental system is that it yields low amounts of myelin. Here we report a simple and efficient method to enhance myelination in vitro. We show that the addition of heparin or low molecular weight heparin to mixed DRG cultures markedly increases Schwann cells myelination. The myelin promoting activity of heparin results from specific inhibition of the soluble immunoglobulin (Ig)-containing isoforms of neuregulin 1 (i.e., NRG1 types I and II) that negatively regulates myelination. Heparin supplement provides a robust and reproducible method to increase myelination in a simple and commonly used culture system. GLIA 2016;64:1227-1234.

Caspr and caspr2 are required for both radial and longitudinal organization of myelinated axons.

In myelinated peripheral axons, Kv1 potassium channels are clustered at the juxtaparanodal region and at an internodal line located along the mesaxon and below the Schmidt-Lanterman incisures. This polarized distribution is controlled by Schwann cells and requires specific cell adhesion molecules (CAMs). The accumulation of Kv1 channels at the juxtaparanodal region depends on the presence of Caspr2 at this site, as well as on the presence of Caspr at the adjacent paranodal junction. However, the localization of these channels along the mesaxonal internodal line still persists in the absence of each one of these CAMs. By generating mice lacking both Caspr and Caspr2 (caspr(-/-)/caspr2(-/-)), we now reveal compensatory functions of the two proteins in the organization of the axolemma. Although Kv1 channels are clustered along the inner mesaxon and in a circumferential ring below the incisures in the single mutants, in sciatic nerves of caspr(-/-)/caspr2(-/-) mice, these channels formed large aggregates that were dispersed along the axolemma, demonstrating that internodal localization of Kv1 channels requires either Caspr or Caspr2. Furthermore, deletion of both Caspr and Caspr2 also resulted in widening of the nodes of Ranvier, suggesting that Caspr2 (which is present at paranodes in the absence of Caspr) can partially compensate for the barrier function of Caspr at this site even without the formation of a distinct paranodal junction. Our results indicate that Caspr and Caspr2 are required for the organization of the axolemma both radially, manifested as the mesaxonal line, and longitudinally, demarcated by the nodal domains.

Twins discordant for myositis and systemic lupus erythematosus show markedly enriched autoantibodies in the affected twin supporting environmental influences in pathogenesis.

BACKGROUND: Studies of twin pairs discordant for autoimmune conditions provide a unique opportunity to explore contributing factors triggered by complex gene-environment interactions. METHODS: In this cross-sectional study, thirty-one monozygotic or dizygotic twin pairs discordant for myositis or systemic lupus erythematosus (SLE), along with matched healthy controls were evaluated for antibodies against a panel of 21 autoantigens. RESULTS: Autoantibody profiling revealed that 42% of the affected twins showed significant seropositivity against autoantigens in the panel. In many of these affected twins, but none of healthy controls, there were high levels of autoantibodies detected against two or more autoantigens commonly seen in systemic autoimmune diseases including Ro52, Ro60, RNP-70 K and/or RNP-A. In contrast, only 10% (3/31) of the unaffected twins showed seropositivity and these immunoreactivities were against single autoantigens not seen in systemic autoimmune diseases. While no significant differences in autoantibodies were detected between the affected or unaffected twins against thyroid peroxidase, transglutaminase and several cytokines, 23% of the affected twins with myositis showed autoantibodies against the gastric ATPase. Analysis of the monozygotic twins separately also revealed a higher frequencies of autoantibodies in the affected twins compared to the unaffected twins (P = 0.046). Lastly, clinical analysis of both the affected monozygotic and dizygotic twins revealed that the autoantibody seropositive affected twins had a greater global disease activity score compared to seronegative affected twins (P = 0.019). CONCLUSION: The findings of significantly more autoantibodies in the affected twins with myositis and SLE compared to the unaffected twins are consistent with potential non-genetic factors playing a role in autoantibody production and pathogenesis of these autoimmune disorders.

Dissecting the molecular basis of human interneuron migration in forebrain assembloids from Timothy syndrome.

Defects in interneuron migration can disrupt the assembly of cortical circuits and lead to neuropsychiatric disease. Using forebrain assembloids derived by integration of cortical and ventral forebrain organoids, we have previously discovered a cortical interneuron migration defect in Timothy syndrome (TS), a severe neurodevelopmental disease caused by a mutation in the L-type calcium channel (LTCC) Cav1.2. Here, we find that acute pharmacological modulation of Cav1.2 can regulate the saltation length, but not the frequency, of interneuron migration in TS. Interestingly, the defect in saltation length is related to aberrant actomyosin and myosin light chain (MLC) phosphorylation, while the defect in saltation frequency is driven by enhanced γ-aminobutyric acid (GABA) sensitivity and can be restored by GABA-A receptor antagonism. Finally, we describe hypersynchronous hCS network activity in TS that is exacerbated by interneuron migration. Taken together, these studies reveal a complex role of LTCC function in human cortical interneuron migration and strategies to restore deficits in the context of disease.

Long-term maturation of human cortical organoids matches key early postnatal transitions.

Human stem-cell-derived models provide the promise of accelerating our understanding of brain disorders, but not knowing whether they possess the ability to mature beyond mid- to late-fetal stages potentially limits their utility. We leveraged a directed differentiation protocol to comprehensively assess maturation in vitro. Based on genome-wide analysis of the epigenetic clock and transcriptomics, as well as RNA editing, we observe that three-dimensional human cortical organoids reach postnatal stages between 250 and 300 days, a timeline paralleling in vivo development. We demonstrate the presence of several known developmental milestones, including switches in the histone deacetylase complex and NMDA receptor subunits, which we confirm at the protein and physiological levels. These results suggest that important components of an intrinsic in vivo developmental program persist in vitro. We further map neurodevelopmental and neurodegenerative disease risk genes onto in vitro gene expression trajectories to provide a resource and webtool (Gene Expression in Cortical Organoids, GECO) to guide disease modeling.

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.

Neuronal defects in a human cellular model of 22q11.2 deletion syndrome.

22q11.2 deletion syndrome (22q11DS) is a highly penetrant and common genetic cause of neuropsychiatric disease. Here we generated induced pluripotent stem cells from 15 individuals with 22q11DS and 15 control individuals and differentiated them into three-dimensional (3D) cerebral cortical organoids. Transcriptional profiling across 100 days showed high reliability of differentiation and revealed changes in neuronal excitability-related genes. Using electrophysiology and live imaging, we identified defects in spontaneous neuronal activity and calcium signaling in both organoid- and 2D-derived cortical neurons. The calcium deficit was related to resting membrane potential changes that led to abnormal inactivation of voltage-gated calcium channels. Heterozygous loss of DGCR8 recapitulated the excitability and calcium phenotypes and its overexpression rescued these defects. Moreover, the 22q11DS calcium abnormality could also be restored by application of antipsychotics. Taken together, our study illustrates how stem cell derived models can be used to uncover and rescue cellular phenotypes associated with genetic forms of neuropsychiatric disease.

CYFIP1 overexpression increases fear response in mice but does not affect social or repetitive behavioral phenotypes.

Background: CYFIP1, a protein that interacts with FMRP and regulates protein synthesis and actin dynamics, is overexpressed in Dup15q syndrome as well as autism spectrum disorder (ASD). While CYFIP1 heterozygosity has been rigorously studied due to its loss in 15q11.2 deletion, Prader-Willi and Angelman syndrome, the effects of CYFIP1 overexpression, as is observed in patients with CYFIP1 duplication, are less well understood. Methods: We developed and validated a mouse model of human CYFIP1 overexpression (CYFIP1 OE) using qPCR and western blot analysis. We performed a large battery of behavior testing on these mice, including ultrasonic vocalizations, three-chamber social assay, home-cage behavior, Y-maze, elevated plus maze, open field test, Morris water maze, fear conditioning, prepulse inhibition, and the hot plate assay. We also performed RNA sequencing and analysis on the basolateral amygdala. Results: Extensive behavioral testing in CYFIP1 OE mice reveals no changes in the core behaviors related to ASD: social interactions and repetitive behaviors. However, we did observe mild learning deficits and an exaggerated fear response. Using RNA sequencing of the basolateral amygdala, a region associated with fear response, we observed changes in pathways related to cytoskeletal regulation, oligodendrocytes, and myelination. We also identified GABA-A subunit composition changes in basolateral amygdala neurons, which are essential components of the neural fear conditioning circuit. Conclusion: Overall, this research identifies the behavioral and molecular consequences of CYFIP1 overexpression and how they contribute to the variable phenotype seen in Dup15q syndrome and in ASD patients with excess CYFIP1.

Exogenous and evoked oxytocin restores social behavior in the Cntnap2 mouse model of autism.

Mouse models of neuropsychiatric diseases provide a platform for mechanistic understanding and development of new therapies. We previously demonstrated that knockout of the mouse homolog of CNTNAP2 (contactin-associated protein-like 2), in which mutations cause cortical dysplasia and focal epilepsy (CDFE) syndrome, displays many features that parallel those of the human disorder. Because CDFE has high penetrance for autism spectrum disorder (ASD), we performed an in vivo screen for drugs that ameliorate abnormal social behavior in Cntnap2 mutant mice and found that acute administration of the neuropeptide oxytocin improved social deficits. We found a decrease in the number of oxytocin immunoreactive neurons in the paraventricular nucleus (PVN) of the hypothalamus in mutant mice and an overall decrease in brain oxytocin levels. Administration of a selective melanocortin receptor 4 agonist, which causes endogenous oxytocin release, also acutely rescued the social deficits, an effect blocked by an oxytocin antagonist. We confirmed that oxytocin neurons mediated the behavioral improvement by activating endogenous oxytocin neurons in the paraventricular hypothalamus with Designer Receptors Exclusively Activated by Designer Drugs (DREADD). Last, we showed that chronic early postnatal treatment with oxytocin led to more lasting behavioral recovery and restored oxytocin immunoreactivity in the PVN. These data demonstrate dysregulation of the oxytocin system in Cntnap2 knockout mice and suggest that there may be critical developmental windows for optimal treatment to rectify this deficit.