SummaryWe present the largest exome sequencing study of autism spectrum disorder (ASD) to date (n=35,584 total samples, 11,986 with ASD). Using an enhanced Bayesian framework to integrate de novo and case-control rare variation, we identify 102 risk genes at a false discovery rate ≤ 0.1. Of these genes, 49 show higher frequencies of disruptive de novo variants in individuals ascertained for severe neurodevelopmental delay, while 53 show higher frequencies in individuals ascertained for ASD; comparing ASD cases with mutations in these groups reveals phenotypic differences. Expressed early in brain development, most of the risk genes have roles in regulation of gene expression or neuronal communication (i.e., mutations effect neurodevelopmental and neurophysiological changes), and 13 fall within loci recurrently hit by copy number variants. In human cortex single-cell gene expression data, expression of risk genes is enriched in both excitatory and inhibitory neuronal lineages, consistent with multiple paths to an excitatory/inhibitory imbalance underlying ASD.

No abstract

Recent research on disparate psychiatric disorders has implicated rare variants in genes involved in global gene regulation and chromatin modification, as well as many common variants located primarily in regulatory regions of the genome. Understanding precisely how these variants contribute to disease will require a deeper appreciation for the mechanisms of gene regulation in the developing and adult human brain. The PsychENCODE project aims to produce a public resource of multidimensional genomic data using tissue- and cell type–specific samples from approximately 1,000 phenotypically well-characterized, high-quality healthy and disease-affected human post-mortem brains, as well as functionally characterize disease-associated regulatory elements and variants in model systems. We are beginning with a focus on autism spectrum disorder, bipolar disorder and schizophrenia, and expect that this knowledge will apply to a wide variety of psychiatric disorders. This paper outlines the motivation and design of PsychENCODE.

Neuropsychiatric disorders have a complex genetic architecture. Human genetic population-based studies have identified numerous heritable sequence and structural genomic variants associated with susceptibility to neuropsychiatric disease. However, these germline variants do not fully account for disease risk. During brain development, progenitor cells undergo billions of cell divisions to generate the ~80 billion neurons in the brain. The failure to accurately repair DNA damage arising during replication, transcription, and cellular metabolism amid this dramatic cellular expansion can lead to somatic mutations. Somatic mutations that alter subsets of neuronal transcriptomes and proteomes can, in turn, affect cell proliferation and survival and lead to neurodevelopmental disorders. The long life span of individual neurons and the direct relationship between neural circuits and behavior suggest that somatic mutations in small populations of neurons can significantly affect individual neurodevelopment. The Brain Somatic Mosaicism Network has been founded to study somatic mosaicism both in neurotypical human brains and in the context of complex neuropsychiatric disorders.

A GAG deletion in the DYT1 gene is a major cause of early-onset dystonia, but clinical disease expression occurs in only 30% of mutation carriers. To gain insight into genetic factors that may influence penetrance, we evaluated three DYT1 single-nucleotide polymorphisms, including D216H, a coding-sequence variation that moderates the effects of the DYT1 GAG deletion in cellular models. We tested DYT1 GAG-deletion carriers with (n=119) and without (n=113) clinical signs of dystonia and control individuals (n=197) and found the frequency of the 216H allele to be increased in GAG-deletion carriers without dystonia and to be decreased in carriers with dystonia, compared with the control individuals. Analysis of haplotypes demonstrated a highly protective effect of the H allele in trans with the GAG deletion; there was also suggestive evidence that the D216 allele in cis is required for the disease to be penetrant. Our findings establish, for the first time, a clinically relevant gene modifier of DYT1.

The DYT6 gene for primary torsion dystonia (PTD) was mapped to chromosome 8p21-q22 in two Amish-Mennonite families who shared a haplotype of marker alleles across a 40 cM linked region. The objective of this study was to narrow the DYT6 region, clinically characterize DYT6 dystonia in a larger cohort, and to determine whether DYT6 is associated with dystonia in newly ascertained multiplex families. We systematically examined familial Amish-Mennonite dystonia cases, identifying five additional members from the original families, as well as three other multiplex Amish-Mennonite families, and evaluated the known DYT6 haplotype and recombination events. One of the three new families carried the shared haplotype, whereas the region was excluded in the two other families, suggesting genetic heterogeneity for PTD in the Amish-Mennonites. Clinical features in the five newly identified DYT6 carriers were similar to those initially described. In contrast, affected individuals from the excluded families had a later age of onset (46.9 years vs. 16.1 years in the DYT6), and the dystonia was both more likely to be of focal distribution and begin in the cervical muscles. Typing of additional markers in the DYT6-linked families revealed recombinations that now place the gene in a 23 cM region surrounding the centromere. In summary, the DYT6 gene is in a 23 cM region on chromosome 8q21-22 and does not account for all familial PTD in Amish-Mennonites.

The availability of high-quality RNA-sequencing and genotyping data of post-mortem brain collections from consortia such as CommonMind Consortium (CMC) and the Accelerating Medicines Partnership for Alzheimer’s Disease (AMP-AD) Consortium enable the generation of a large-scale brain cis-eQTL meta-analysis. Here we generate cerebral cortical eQTL from 1433 samples available from four cohorts (identifying >4.1 million significant eQTL for >18,000 genes), as well as cerebellar eQTL from 261 samples (identifying 874,836 significant eQTL for >10,000 genes). We find substantially improved power in the meta-analysis over individual cohort analyses, particularly in comparison to the Genotype-Tissue Expression (GTEx) Project eQTL. Additionally, we observed differences in eQTL patterns between cerebral and cerebellar brain regions. We provide these brain eQTL as a resource for use by the research community. As a proof of principle for their utility, we apply a colocalization analysis to identify genes underlying the GWAS association peaks for schizophrenia and identify a potentially novel gene colocalization with lncRNA RP11-677M14.2 (posterior probability of colocalization 0.975).

Microarrays derived from Solanum tuberosum expressed sequence tags were used to test the hypothesis that genetically distinct enveloped viruses elicit unique changes in Nicotiana benthamiana gene expression. The results of our study, which included Sonchus yellow net virus (SYNV), a plant rhabdovirus that replicates in the nucleus of infected cells, and Impatiens necrotic spot virus (INSV), a plant bunyavirus that replicates in the cytoplasm, were consistent with this hypothesis. Statistically significant changes (P¡0?01) in the expression of 275, 2646 and 4165 genes were detected in response to INSV at 2, 4 and 5 days post-inoculation (d.p.i.), respectively. In contrast, 35, 665 and 1458 genes were expressed differentially in response to SYNV at 5, 11 and 14 d.p.i., respectively. The microarray results were verified by Northern hybridization using a subset of these genes as probes. Notably, INSV, but not SYNV, induced expression of small heat-shock protein genes to high levels. In contrast to SYNV, infection by INSV resulted in downregulation of all histone genes, of which the downregulation of histone 2b expression to very low levels was confirmed by Northern hybridization. The expression of a putative WRKY transcription factor at 11 d.p.i., but not at 5 or 14 d.p.i., in SYNV-infected tissue suggested that the temporal response to virus infection was identified readily using our experimental design. Overall, infection by INSV resulted in larger fold changes in host gene expression relative to infection by SYNV. Taken together, the present data demonstrate differential responses of a common host to two genetically distinct viruses.