A foundational atlas of autism protein interactions reveals molecular convergence

Translating high-confidence (hc) autism spectrum disorder (ASD) genes into viable treatment targets remains elusive. We constructed a foundational protein-protein interaction (PPI) network in HEK293T cells involving 100 hcASD risk genes, revealing over 1,800 PPIs (87% novel). Interactors, expressed in the human brain and enriched for ASD but not schizophrenia genetic risk, converged on protein complexes involved in neurogenesis, tubulin biology, transcriptional regulation, and chromatin modification. A PPI map of 54 patient-derived missense variants identified differential physical interactions, and we leveraged AlphaFold-Multimer predictions to prioritize direct PPIs and specific variants for interrogation in Xenopus tropicalis and human forebrain organoids. A mutation in the transcription factor FOXP1 led to reconfiguration of DNA binding sites and altered development of deep cortical layer neurons in forebrain organoids. This work offers new insights into molecular mechanisms underlying ASD and describes a powerful platform to develop and test therapeutic strategies for many genetically-defined conditions.

Autism genes converge on microtubule biology and RNA-binding proteins during excitatory neurogenesis

Recent studies have identified over one hundred high-confidence (hc) autism spectrum disorder (ASD) genes. Systems biological and functional analyses on smaller subsets of these genes have consistently implicated excitatory neurogenesis. However, the extent to which the broader set of hcASD genes are involved in this process has not been explored systematically nor have the biological pathways underlying this convergence been identified. Here, we leveraged CROP-Seq to repress 87 hcASD genes in a human in vitro model of cortical neurogenesis. We identified 17 hcASD genes whose repression significantly alters developmental trajectory and results in a common cellular state characterized by disruptions in proliferation, differentiation, cell cycle, microtubule biology, and RNA-binding proteins (RBPs). We also characterized over 3,000 differentially expressed genes, 286 of which had expression profiles correlated with changes in developmental trajectory. Overall, we uncovered transcriptional disruptions downstream of hcASD gene perturbations, correlated these disruptions with distinct differentiation phenotypes, and reinforced neurogenesis, microtubule biology, and RBPs as convergent points of disruption in ASD.

Rare X-linked variants carry predominantly male risk in autism, Tourette syndrome, and ADHD

Autism spectrum disorder (ASD), Tourette syndrome (TS), and attention-deficit/hyperactivity disorder (ADHD) display strong male sex bias, due to a combination of genetic and biological factors, as well as selective ascertainment. While the hemizygous nature of chromosome X (Chr X) in males has long been postulated as a key point of "male vulnerability", rare genetic variation on this chromosome has not been systematically characterized in large-scale whole exome sequencing studies of "idiopathic" ASD, TS, and ADHD. Here, we take advantage of informative recombinations in simplex ASD families to pinpoint risk-enriched regions on Chr X, within which rare maternally-inherited damaging variants carry substantial risk in males with ASD. We then apply a modified transmission disequilibrium test to 13,052 ASD probands and identify a novel high confidence ASD risk gene at exome-wide significance (MAGEC3). Finally, we observe that rare damaging variants within these risk regions carry similar effect sizes in males with TS or ADHD, further clarifying genetic mechanisms underlying male vulnerability in multiple neurodevelopmental disorders that can be exploited for systematic gene discovery.

Pleiotropy of autism-associated chromatin regulators

Gene ontology analyses of high-confidence autism spectrum disorder (ASD) risk genes highlight chromatin regulation and synaptic function as major contributors to pathobiology. Our recent functional work in vivo has additionally implicated tubulin biology and cellular proliferation. As many chromatin regulators, including the ASD risk genes ADNP and CHD3, are known to directly regulate both tubulins and histones, we studied the five chromatin regulators most strongly associated with ASD (ADNP, CHD8, CHD2, POGZ and KMT5B) specifically with respect to tubulin biology. We observe that all five localize to microtubules of the mitotic spindle in vitro in human cells and in vivo in Xenopus. Investigation of CHD2 provides evidence that mutations present in individuals with ASD cause a range of microtubule-related phenotypes, including disrupted localization of the protein at mitotic spindles, cell cycle stalling, DNA damage and cell death. Lastly, we observe that ASD genetic risk is significantly enriched among tubulin-associated proteins, suggesting broader relevance. Together, these results provide additional evidence that the role of tubulin biology and cellular proliferation in ASD warrants further investigation and highlight the pitfalls of relying solely on annotated gene functions in the search for pathological mechanisms.

Genome-wide Association Study Points to Novel Locus for Gilles de la Tourette Syndrome

BACKGROUND

Tourette syndrome (TS) is a childhood-onset neurodevelopmental disorder of complex genetic architecture and is characterized by multiple motor tics and at least one vocal tic persisting for more than 1 year.

METHODS

We performed a genome-wide meta-analysis integrating a novel TS cohort with previously published data, resulting in a sample size of 6133 individuals with TS and 13,565 ancestry-matched control participants.

RESULTS

We identified a genome-wide significant locus on chromosome 5q15. Integration of expression quantitative trait locus, Hi-C (high-throughput chromosome conformation capture), and genome-wide association study data implicated the NR2F1 gene and associated long noncoding RNAs within the 5q15 locus. Heritability partitioning identified statistically significant enrichment in brain tissue histone marks, while polygenic risk scoring of brain volume data identified statistically significant associations with right and left thalamus volumes and right putamen volume.

CONCLUSIONS

Our work presents novel insights into the neurobiology of TS, thereby opening up new directions for future studies.

Genomics, convergent neuroscience and progress in understanding autism spectrum disorder

More than a hundred genes have been identified that, when disrupted, impart large risk for autism spectrum disorder (ASD). Current knowledge about the encoded proteins - although incomplete - points to a very wide range of developmentally dynamic and diverse biological processes. Moreover, the core symptoms of ASD involve distinctly human characteristics, presenting challenges to interpreting evolutionarily distant model systems. Indeed, despite a decade of striking progress in gene discovery, an actionable understanding of pathobiology remains elusive. Increasingly, convergent neuroscience approaches have been recognized as an important complement to traditional uses of genetics to illuminate the biology of human disorders. These methods seek to identify intersection among molecular-level, cellular-level and circuit-level functions across multiple risk genes and have highlighted developing excitatory neurons in the human mid-gestational prefrontal cortex as an important pathobiological nexus in ASD. In addition, neurogenesis, chromatin modification and synaptic function have emerged as key potential mediators of genetic vulnerability. The continued expansion of foundational 'omics' data sets, the application of higher-throughput model systems and incorporating developmental trajectories and sex differences into future analyses will refine and extend these results. Ultimately, a systems-level understanding of ASD genetic risk holds promise for clarifying pathobiology and advancing therapeutics.

Enhancing WNT Signaling Restores Cortical Neuronal Spine Maturation and Synaptogenesis in Tbr1 Mutants

Tbr1 is a high-confidence autism spectrum disorder (ASD) gene encoding a transcription factor with distinct pre- and postnatal functions. Postnatally, Tbr1 conditional knockout (CKO) mutants and constitutive heterozygotes have immature dendritic spines and reduced synaptic density. Tbr1 regulates expression of several genes that underlie synaptic defects, including a kinesin (Kif1a) and a WNT-signaling ligand (Wnt7b). Furthermore, Tbr1 mutant corticothalamic neurons have reduced thalamic axonal arborization. LiCl and a GSK3β inhibitor, two WNT-signaling agonists, robustly rescue the dendritic spines and the synaptic and axonal defects, suggesting that this could have relevance for therapeutic approaches in some forms of ASD.

Fazel Darbandi et al. demonstrate that TBR1 directly regulates transcriptional circuits in cortical layers 5 and 6, which promote dendritic spine and synaptic density. Enhancing WNT signaling rescues dendritic spine maturation and synaptogenesis defects in Tbr1 mutants. These results provide insights into mechanisms that underlie ASD pathophysiology.

Parallel in vivo analysis of large-effect autism genes implicates cortical neurogenesis and estrogen in risk and resilience

Gene Ontology analyses of autism spectrum disorders (ASD) risk genes have repeatedly highlighted synaptic function and transcriptional regulation as key points of convergence. However, these analyses rely on incomplete knowledge of gene function across brain development. Here we leverage Xenopus tropicalis to study in vivo ten genes with the strongest statistical evidence for association with ASD. All genes are expressed in developing telencephalon at time points mapping to human mid-prenatal development, and mutations lead to an increase in the ratio of neural progenitor cells to maturing neurons, supporting previous in silico systems biological findings implicating cortical neurons in ASD vulnerability, but expanding the range of convergent functions to include neurogenesis. Systematic chemical screening identifies that estrogen, via Sonic hedgehog signaling, rescues this convergent phenotype in Xenopus and human models of brain development, suggesting a resilience factor that may mitigate a range of ASD genetic risks.

Synaptic processes and immune-related pathways implicated in Tourette syndrome

Tourette syndrome (TS) is a neuropsychiatric disorder of complex genetic architecture involving multiple interacting genes. Here, we sought to elucidate the pathways that underlie the neurobiology of the disorder through genome-wide analysis. We analyzed genome-wide genotypic data of 3581 individuals with TS and 7682 ancestry-matched controls and investigated associations of TS with sets of genes that are expressed in particular cell types and operate in specific neuronal and glial functions. We employed a self-contained, set-based association method (SBA) as well as a competitive gene set method (MAGMA) using individual-level genotype data to perform a comprehensive investigation of the biological background of TS. Our SBA analysis identified three significant gene sets after Bonferroni correction, implicating ligand-gated ion channel signaling, lymphocytic, and cell adhesion and transsynaptic signaling processes. MAGMA analysis further supported the involvement of the cell adhesion and trans-synaptic signaling gene set. The lymphocytic gene set was driven by variants in FLT3, raising an intriguing hypothesis for the involvement of a neuroinflammatory element in TS pathogenesis. The indications of involvement of ligand-gated ion channel signaling reinforce the role of GABA in TS, while the association of cell adhesion and trans-synaptic signaling gene set provides additional support for the role of adhesion molecules in neuropsychiatric disorders. This study reinforces previous findings but also provides new insights into the neurobiology of TS.

The neurodevelopmental disorder risk gene DYRK1A is required for ciliogenesis and control of brain size in Xenopus embryos

DYRK1A [dual specificity tyrosine-(Y)-phosphorylation-regulated kinase 1 A] is a high-confidence autism risk gene that encodes a conserved kinase. In addition to autism, individuals with putative loss-of-function variants in DYRK1A exhibit microcephaly, intellectual disability, developmental delay and/or congenital anomalies of the kidney and urinary tract. DYRK1A is also located within the critical region for Down syndrome; therefore, understanding the role of DYRK1A in brain development is crucial for understanding the pathobiology of multiple developmental disorders. To characterize the function of this gene, we used the diploid frog Xenopus tropicalis We discover that Dyrk1a is expressed in ciliated tissues, localizes to ciliary axonemes and basal bodies, and is required for ciliogenesis. We also demonstrate that Dyrk1a localizes to mitotic spindles and that its inhibition leads to decreased forebrain size, abnormal cell cycle progression and cell death during brain development. These findings provide hypotheses about potential mechanisms of pathobiology and underscore the utility of X. tropicalis as a model system for understanding neurodevelopmental disorders.

Human 3D cellular model of hypoxic brain injury of prematurity

Owing to recent medical and technological advances in neonatal care, infants born extremely premature have increased survival rates1,2. After birth, these infants are at high risk of hypoxic episodes because of lung immaturity, hypotension and lack of cerebral-flow regulation, and can develop a severe condition called encephalopathy of prematurity3. Over 80% of infants born before post-conception week 25 have moderate-to-severe long-term neurodevelopmental impairments4. The susceptible cell types in the cerebral cortex and the molecular mechanisms underlying associated gray-matter defects in premature infants remain unknown. Here we used human three-dimensional brain-region-specific organoids to study the effect of oxygen deprivation on corticogenesis. We identified specific defects in intermediate progenitors, a cortical cell type associated with the expansion of the human cerebral cortex, and showed that these are related to the unfolded protein response and changes. Moreover, we verified these findings in human primary cortical tissue and demonstrated that a small-molecule modulator of the unfolded protein response pathway can prevent the reduction in intermediate progenitors following hypoxia. We anticipate that this human cellular platform will be valuable for studying the environmental and genetic factors underlying injury in the developing human brain.

Genome-wide de novo risk score implicates promoter variation in autism spectrum disorder

Whole-genome sequencing (WGS) has facilitated the first genome-wide evaluations of the contribution of de novo noncoding mutations to complex disorders. Using WGS, we identified 255,106 de novo mutations among sample genomes from members of 1902 quartet families in which one child, but not a sibling or their parents, was affected by autism spectrum disorder (ASD). In contrast to coding mutations, no noncoding functional annotation category, analyzed in isolation, was significantly associated with ASD. Casting noncoding variation in the context of a de novo risk score across multiple annotation categories, however, did demonstrate association with mutations localized to promoter regions. We found that the strongest driver of this promoter signal emanates from evolutionarily conserved transcription factor binding sites distal to the transcription start site. These data suggest that de novo mutations in promoter regions, characterized by evolutionary and functional signatures, contribute to ASD.

Comparative Flavivirus-Host Protein Interaction Mapping Reveals Mechanisms of Dengue and Zika Virus Pathogenesis

Mosquito-borne flaviviruses, including dengue virus (DENV) and Zika virus (ZIKV), are a growing public health concern. Systems-level analysis of how flaviviruses hijack cellular processes through virus-host protein-protein interactions (PPIs) provides information about their replication and pathogenic mechanisms. We used affinity purification-mass spectrometry (AP-MS) to compare flavivirus-host interactions for two viruses (DENV and ZIKV) in two hosts (human and mosquito). Conserved virus-host PPIs revealed that the flavivirus NS5 protein suppresses interferon stimulated genes by inhibiting recruitment of the transcription complex PAF1C and that chemical modulation of SEC61 inhibits DENV and ZIKV replication in human and mosquito cells. Finally, we identified a ZIKV-specific interaction between NS4A and ANKLE2, a gene linked to hereditary microcephaly, and showed that ZIKV NS4A causes microcephaly in Drosophila in an ANKLE2-dependent manner. Thus, comparative flavivirus-host PPI mapping provides biological insights and, when coupled with in vivo models, can be used to unravel pathogenic mechanisms.

Neonatal Tbr1 Dosage Controls Cortical Layer 6 Connectivity

An understanding of how heterozygous loss-of-function mutations in autism spectrum disorder (ASD) risk genes, such as TBR1, contribute to ASD remains elusive. Conditional Tbr1 deletion during late mouse gestation in cortical layer 6 neurons (Tbr1^layer6 mutants) provides novel insights into its function, including dendritic patterning, synaptogenesis, and cell-intrinsic physiology. These phenotypes occur in heterozygotes, providing insights into mechanisms that may underlie ASD pathophysiology. Restoring expression of Wnt7b largely rescues the synaptic deficit in Tbr1^layer6 mutant neurons. Furthermore, Tbr1^layer6 heterozygotes have increased anxiety-like behavior, a phenotype seen ASD. Integrating TBR1 chromatin immunoprecipitation sequencing (ChIP-seq) and RNA sequencing (RNA-seq) data from layer 6 neurons and activity of TBR1-bound candidate enhancers provides evidence for how TBR1 regulates layer 6 properties. Moreover, several putative TBR1 targets are ASD risk genes, placing TBR1 in a central position both for ASD risk and for regulating transcriptional circuits that control multiple steps in layer 6 development essential for the assembly of neural circuits.

De Novo Sequence and Copy Number Variants Are Strongly Associated with Tourette Disorder and Implicate Cell Polarity in Pathogenesis

We previously established the contribution of de novo damaging sequence variants to Tourette disorder (TD) through whole-exome sequencing of 511 trios. Here, we sequence an additional 291 TD trios and analyze the combined set of 802 trios. We observe an overrepresentation of de novo damaging variants in simplex, but not multiplex, families; we identify a high-confidence TD risk gene, CELSR3 (cadherin EGF LAG seven-pass G-type receptor 3); we find that the genes mutated in TD patients are enriched for those related to cell polarity, suggesting a common pathway underlying pathobiology; and we confirm a statistically significant excess of de novo copy number variants in TD. Finally, we identify significant overlap of de novo sequence variants between TD and obsessive-compulsive disorder and de novo copy number variants between TD and autism spectrum disorder, consistent with shared genetic risk.

The Psychiatric Cell Map Initiative: A Convergent Systems Biological Approach to Illuminating Key Molecular Pathways in Neuropsychiatric Disorders

Although gene discovery in neuropsychiatric disorders, including autism spectrum disorder, intellectual disability, epilepsy, schizophrenia, and Tourette disorder, has accelerated, resulting in a large number of molecular clues, it has proven difficult to generate specific hypotheses without the corresponding datasets at the protein complex and functional pathway level. Here, we describe one path forward-an initiative aimed at mapping the physical and genetic interaction networks of these conditions and then using these maps to connect the genomic data to neurobiology and, ultimately, the clinic. These efforts will include a team of geneticists, structural biologists, neurobiologists, systems biologists, and clinicians, leveraging a wide array of experimental approaches and creating a collaborative infrastructure necessary for long-term investigation. This initiative will ultimately intersect with parallel studies that focus on other diseases, as there is a significant overlap with genes implicated in cancer, infectious disease, and congenital heart defects.

Publisher Correction: Whole genome sequencing in psychiatric disorders: the WGSPD consortium

In the version of this article initially published, the consortium authorship and corresponding authors were not presented correctly. In the PDF and print versions, the Whole Genome Sequencing for Psychiatric Disorders (WGSPD) consortium was missing from the author list at the beginning of the paper, where it should have appeared as the seventh author; it was present in the author list at the end of the paper, but the footnote directing readers to the Supplementary Note for a list of members was missing. In the HTML version, the consortium was listed as the last author instead of as the seventh, and the line directing readers to the Supplementary Note for a list of members appeared at the end of the paper under Author Information but not in association with the consortium name itself. Also, this line stated that both member names and affiliations could be found in the Supplementary Note; in fact, only names are given. In all versions of the paper, the corresponding author symbols were attached to A. Jeremy Willsey, Steven E. Hyman, Anjene M. Addington and Thomas Lehner; they should have been attached, respectively, to Steven E. Hyman, Anjene M. Addington, Thomas Lehner and Nelson B. Freimer. As a result of this shift, the respective contact links in the HTML version did not lead to the indicated individuals. The errors have been corrected in the HTML and PDF versions of the article.

Whole genome sequencing in psychiatric disorders: the WGSPD consortium

As technology advances, whole genome sequencing (WGS) is likely to supersede other genotyping technologies. The rate of this change depends on its relative cost and utility. Variants identified uniquely through WGS may reveal novel biological pathways underlying complex disorders and provide high-resolution insight into when, where, and in which cell type these pathways are affected. Alternatively, cheaper and less computationally intensive approaches may yield equivalent insights. Understanding the role of rare variants in the noncoding gene-regulating genome, through pilot WGS projects, will be critical to determine which of these two extremes best represents reality. With large cohorts, well-defined risk loci, and a compelling need to understand the underlying biology, psychiatric disorders have a role to play in this preliminary WGS assessment. The WGSPD consortium will integrate data for 18,000 individuals with psychiatric disorders, beginning with autism spectrum disorder, schizophrenia, bipolar disorder, and major depressive disorder, along with over 150,000 controls.

Identification of Developmental and Behavioral Markers Associated With Genetic Abnormalities in Autism Spectrum Disorder

OBJECTIVE

Aside from features associated with risk of neurogenetic syndromes in general (e.g., cognitive impairment), limited progress has been made in identifying phenotype-genotype relationships in autism spectrum disorder (ASD). The objective of this study was to extend work in the Simons Simplex Collection by comparing the phenotypic profiles of ASD probands with or without identified de novo loss of function mutations or copy number variants in high-confidence ASD-associated genes or loci.

METHOD

Analyses preemptively accounted for documented differences in sex and IQ in affected individuals with de novo mutations by matching probands with and without these genetic events on sex, IQ, and age before comparing them on multiple behavioral domains.

RESULTS

Children with de novo mutations (N=112) had a greater likelihood of motor delay during early development (later age at walking), but they were less impaired on certain measures of ASD core symptoms (parent-rated social communication abnormalities and clinician-rated diagnostic certainty about ASD) in later childhood. These children also showed relative strengths in verbal and language abilities, including a smaller discrepancy between nonverbal and verbal IQ and a greater likelihood of having achieved fluent language (i.e., regular use of complex sentences).

CONCLUSIONS

Children with ASD with de novo mutations may exhibit a "muted" symptom profile with respect to social communication and language deficits relative to those with ASD with no identified genetic abnormalities. Such findings suggest that examining early milestone differences and standardized testing results may be helpful in etiologic efforts, and potentially in clinical differentiation of various subtypes of ASD, but only if developmental and demographic variables are properly accounted for first.

The Genetic Etiology of Tourette Syndrome: Large-Scale Collaborative Efforts on the Precipice of Discovery

Gilles de la Tourette Syndrome (TS) is a childhood-onset neurodevelopmental disorder that is characterized by multiple motor and phonic tics. It has a complex etiology with multiple genes likely interacting with environmental factors to lead to the onset of symptoms. The genetic basis of the disorder remains elusive. However, multiple resources and large-scale projects are coming together, launching a new era in the field and bringing us on the verge of discovery. The large-scale efforts outlined in this report are complementary and represent a range of different approaches to the study of disorders with complex inheritance. The Tourette Syndrome Association International Consortium for Genetics (TSAICG) has focused on large families, parent-proband trios and cases for large case-control designs such as genomewide association studies (GWAS), copy number variation (CNV) scans, and exome/genome sequencing. TIC Genetics targets rare, large effect size mutations in simplex trios, and multigenerational families. The European Multicentre Tics in Children Study (EMTICS) seeks to elucidate gene-environment interactions including the involvement of infection and immune mechanisms in TS etiology. Finally, TS-EUROTRAIN, a Marie Curie Initial Training Network, aims to act as a platform to unify large-scale projects in the field and to educate the next generation of experts. Importantly, these complementary large-scale efforts are joining forces to uncover the full range of genetic variation and environmental risk factors for TS, holding great promise for identifying definitive TS susceptibility genes and shedding light into the complex pathophysiology of this disorder.

Localized JNK signaling regulates organ size during development

A fundamental question of biology is what determines organ size. Despite demonstrations that factors within organs determine their sizes, intrinsic size control mechanisms remain elusive. Here we show that Drosophila wing size is regulated by JNK signaling during development. JNK is active in a stripe along the center of developing wings, and modulating JNK signaling within this stripe changes organ size. This JNK stripe influences proliferation in a non-canonical, Jun-independent manner by inhibiting the Hippo pathway. Localized JNK activity is established by Hedgehog signaling, where Ci elevates dTRAF1 expression. As the dTRAF1 homolog, TRAF4, is amplified in numerous cancers, these findings provide a new mechanism for how the Hedgehog pathway could contribute to tumorigenesis, and, more importantly, provides a new strategy for cancer therapies. Finally, modulation of JNK signaling centers in developing antennae and legs changes their sizes, suggesting a more generalizable role for JNK signaling in developmental organ size control.

DOI:http://dx.doi.org/10.7554/eLife.11491.001

A key challenge in biology is to understand what determines size. As an animal grows, signals are produced that control the size of its organs. Many of the signaling pathways that regulate size during normal animal development also contribute to the formation of tumors. Therefore, it is important to find out exactly how the signaling molecules that regulate size are linked to those that regulate tumor growth.

A protein called JNK activates a signaling pathway that triggers tumor growth. JNK signaling also stimulates cells to multiply in tissues that need repair, but it is not known whether it also regulates the size of organs during animal development. Here, Willsey et al. investigate whether JNK is active in the developing wings of fruit flies, which are commonly used as models of animal development.

The experiments show that JNK is active in a stripe across the developing wing and is required for the wing to grow to its proper size. A master signal protein called Hedgehog is responsible for establishing this stripe of JNK activity. Unexpectedly, rather than acting through its usual signaling pathway, JNK activates the Hippo pathway in the wing to control organ size during development.

Willsey et al.’s findings highlight potential new targets for cancer therapies. A future challenge will be to find out whether small patches of JNK signaling are found in the developing organs of other animals, and whether they can help explain how size changes between species.

DOI:http://dx.doi.org/10.7554/eLife.11491.002

The female protective effect in autism spectrum disorder is not mediated by a single genetic locus

Background

A 4:1 male to female sex bias has consistently been observed in autism spectrum disorder (ASD). Epidemiological and genetic studies suggest a female protective effect (FPE) may account for part of this bias; however, the mechanism of such protection is unknown. Quantitative assessment of ASD symptoms using the Social Responsiveness Scale (SRS) shows a bimodal distribution unique to females in multiplex families. This leads to the hypothesis that a single, common genetic locus on chromosome X might mediate the FPE and produce the ASD sex bias. Such a locus would represent a major therapeutic target and is likely to have been missed by conventional genome-wide association study (GWAS) analysis.

Methods

To explore this possibility, we performed an association study in affected versus unaffected females, considering three tiers of single nucleotide polymorphisms (SNPs) as follows: 1) regions of chromosome X that escape X-inactivation, 2) all of chromosome X, and 3) genome-wide.

Results

No evidence of a SNP meeting the criteria for a single FPE locus was observed, despite the analysis being well powered to detect this effect.

Conclusions

The results do not support the hypothesis that the FPE is mediated by a single genetic locus; however, this does not exclude the possibility of multiple genetic loci playing a role in the FPE.

Electronic supplementary material

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

A genome-wide association study of autism using the Simons Simplex Collection: Does reducing phenotypic heterogeneity in autism increase genetic homogeneity?

BACKGROUND

Phenotypic heterogeneity in autism has long been conjectured to be a major hindrance to the discovery of genetic risk factors, leading to numerous attempts to stratify children based on phenotype to increase power of discovery studies. This approach, however, is based on the hypothesis that phenotypic heterogeneity closely maps to genetic variation, which has not been tested. Our study examines the impact of subphenotyping of a well-characterized autism spectrum disorder (ASD) sample on genetic homogeneity and the ability to discover common genetic variants conferring liability to ASD.

METHODS

Genome-wide genotypic data of 2576 families from the Simons Simplex Collection were analyzed in the overall sample and phenotypic subgroups defined on the basis of diagnosis, IQ, and symptom profiles. We conducted a family-based association study, as well as estimating heritability and evaluating allele scores for each phenotypic subgroup.

RESULTS

Association analyses revealed no genome-wide significant association signal. Subphenotyping did not increase power substantially. Moreover, allele scores built from the most associated single nucleotide polymorphisms, based on the odds ratio in the full sample, predicted case status in subsets of the sample equally well and heritability estimates were very similar for all subgroups.

CONCLUSIONS

In genome-wide association analysis of the Simons Simplex Collection sample, reducing phenotypic heterogeneity had at most a modest impact on genetic homogeneity. Our results are based on a relatively small sample, one with greater homogeneity than the entire population; if they apply more broadly, they imply that analysis of subphenotypes is not a productive path forward for discovering genetic risk variants in ASD.

The autism-associated chromatin modifier CHD8 regulates other autism risk genes during human neurodevelopment

Recent studies implicate chromatin modifiers in autism spectrum disorder (ASD) through the identification of recurrent de novo loss of function mutations in affected individuals. ASD risk genes are co-expressed in human midfetal cortex, suggesting that ASD risk genes converge in specific regulatory networks during neurodevelopment. To elucidate such networks, we identify genes targeted by CHD8, a chromodomain helicase strongly associated with ASD, in human midfetal brain, human neural stem cells (hNSCs) and embryonic mouse cortex. CHD8 targets are strongly enriched for other ASD risk genes in both human and mouse neurodevelopment, and converge in ASD-associated co-expression networks in human midfetal cortex. CHD8 knockdown in hNSCs results in dysregulation of ASD risk genes directly targeted by CHD8. Integration of CHD8-binding data into ASD risk models improves detection of risk genes. These results suggest loss of CHD8 contributes to ASD by perturbing an ancient gene regulatory network during human brain development.

Autism spectrum disorders: from genes to neurobiology

Advances in genome-wide technology, coupled with the availability of large cohorts, are finally yielding a steady stream of autism spectrum disorder (ASD) genes carrying mutations of large effect. These findings represent important molecular clues, but at the same time present notable challenges to traditional strategies for moving from genes to neurobiology. A remarkable degree of genetic heterogeneity, the biological pleiotropy of ASD genes, and the tremendous complexity of the human brain are prompting the development of new strategies for translating genetic discoveries into therapeutic targets. Recent developments in systems biology approaches that 'contextualize' these genetic findings along spatial, temporal, and cellular axes of human brain development are beginning to bridge the gap between high-throughput gene discovery and testable pathophysiological hypotheses.

The contribution of de novo coding mutations to autism spectrum disorder

Whole exome sequencing has proven to be a powerful tool for understanding the genetic architecture of human disease. Here we apply it to more than 2,500 simplex families, each having a child with an autistic spectrum disorder. By comparing affected to unaffected siblings, we show that 13% of de novo missense mutations and 43% of de novo likely gene-disrupting (LGD) mutations contribute to 12% and 9% of diagnoses, respectively. Including copy number variants, coding de novo mutations contribute to about 30% of all simplex and 45% of female diagnoses. Almost all LGD mutations occur opposite wild-type alleles. LGD targets in affected females significantly overlap the targets in males of lower intelligence quotient (IQ), but neither overlaps significantly with targets in males of higher IQ. We estimate that LGD mutation in about 400 genes can contribute to the joint class of affected females and males of lower IQ, with an overlapping and similar number of genes vulnerable to contributory missense mutation. LGD targets in the joint class overlap with published targets for intellectual disability and schizophrenia, and are enriched for chromatin modifiers, FMRP-associated genes and embryonically expressed genes. Most of the significance for the latter comes from affected females.

De novo insertions and deletions of predominantly paternal origin are associated with autism spectrum disorder

Whole-exome sequencing (WES) studies have demonstrated the contribution of de novo loss-of-function single-nucleotide variants (SNVs) to autism spectrum disorder (ASD). However, challenges in the reliable detection of de novo insertions and deletions (indels) have limited inclusion of these variants in prior analyses. By applying a robust indel detection method to WES data from 787 ASD families (2,963 individuals), we demonstrate that de novo frameshift indels contribute to ASD risk (OR = 1.6; 95% CI = 1.0-2.7; p = 0.03), are more common in female probands (p = 0.02), are enriched among genes encoding FMRP targets (p = 6 × 10(-9)), and arise predominantly on the paternal chromosome (p < 0.001). On the basis of mutation rates in probands versus unaffected siblings, we conclude that de novo frameshift indels contribute to risk in approximately 3% of individuals with ASD. Finally, by observing clustering of mutations in unrelated probands, we uncover two ASD-associated genes: KMT2E (MLL5), a chromatin regulator, and RIMS1, a regulator of synaptic vesicle release.

Increased frequency of de novo copy number variants in congenital heart disease by integrative analysis of single nucleotide polymorphism array and exome sequence data

RATIONALE

Congenital heart disease (CHD) is among the most common birth defects. Most cases are of unknown pathogenesis.

OBJECTIVE

To determine the contribution of de novo copy number variants (CNVs) in the pathogenesis of sporadic CHD.

METHODS AND RESULTS

We studied 538 CHD trios using genome-wide dense single nucleotide polymorphism arrays and whole exome sequencing. Results were experimentally validated using digital droplet polymerase chain reaction. We compared validated CNVs in CHD cases with CNVs in 1301 healthy control trios. The 2 complementary high-resolution technologies identified 63 validated de novo CNVs in 51 CHD cases. A significant increase in CNV burden was observed when comparing CHD trios with healthy trios, using either single nucleotide polymorphism array (P=7×10(-5); odds ratio, 4.6) or whole exome sequencing data (P=6×10(-4); odds ratio, 3.5) and remained after removing 16% of de novo CNV loci previously reported as pathogenic (P=0.02; odds ratio, 2.7). We observed recurrent de novo CNVs on 15q11.2 encompassing CYFIP1, NIPA1, and NIPA2 and single de novo CNVs encompassing DUSP1, JUN, JUP, MED15, MED9, PTPRE SREBF1, TOP2A, and ZEB2, genes that interact with established CHD proteins NKX2-5 and GATA4. Integrating de novo variants in whole exome sequencing and CNV data suggests that ETS1 is the pathogenic gene altered by 11q24.2-q25 deletions in Jacobsen syndrome and that CTBP2 is the pathogenic gene in 10q subtelomeric deletions.

CONCLUSIONS

We demonstrate a significantly increased frequency of rare de novo CNVs in CHD patients compared with healthy controls and suggest several novel genetic loci for CHD.

Modest impact on risk for autism spectrum disorder of rare copy number variants at 15q11.2, specifically breakpoints 1 to 2

The proximal region of chromosome 15 is one of the genomic hotspots for copy number variants (CNVs). Among the rearrangements observed in this region, CNVs from the interval between the common breakpoints 1 and 2 (BP1 and BP2) have been reported cosegregating with autism spectrum disorder (ASD). Although evidence supporting an association between BP1-BP2 CNVs and autism accumulates, the magnitude of the effect of BP1-BP2 CNVs remains elusive, posing a great challenge to recurrence-risk counseling. To gain further insight into their pathogenicity for ASD, we estimated the penetrance of the BP1-BP2 CNVs for ASD as well as their effects on ASD-related phenotypes in a well-characterized ASD sample (n = 2525 families). Transmission disequilibrium test revealed significant preferential transmission only for the duplicated chromosome in probands (20T:9NT). The penetrance of the BP1-BP2 CNVs for ASD was low, conferring additional risks of 0.3% (deletion) and 0.8% (duplication). Stepwise regression analyses suggest a greater effect of the CNVs on ASD-related phenotype in males and when maternally inherited. Taken together, the results are consistent with BP1-BP2 CNVs as risk factors for autism. However, their effect is modest, more akin to that seen for common variants. To be consistent with the current American College of Medical Genetics guidelines for interpretation of postnatal CNV, the BP1-BP2 deletion and duplication CNVs would probably best be classified as variants of uncertain significance (VOUS): they appear to have an impact on risk, but one so modest that these CNVs do not merit pathogenic status.

Coexpression networks implicate human midfetal deep cortical projection neurons in the pathogenesis of autism

Autism spectrum disorder (ASD) is a complex developmental syndrome of unknown etiology. Recent studies employing exome- and genome-wide sequencing have identified nine high-confidence ASD (hcASD) genes. Working from the hypothesis that ASD-associated mutations in these biologically pleiotropic genes will disrupt intersecting developmental processes to contribute to a common phenotype, we have attempted to identify time periods, brain regions, and cell types in which these genes converge. We have constructed coexpression networks based on the hcASD "seed" genes, leveraging a rich expression data set encompassing multiple human brain regions across human development and into adulthood. By assessing enrichment of an independent set of probable ASD (pASD) genes, derived from the same sequencing studies, we demonstrate a key point of convergence in midfetal layer 5/6 cortical projection neurons. This approach informs when, where, and in what cell types mutations in these specific genes may be productively studied to clarify ASD pathophysiology.

Adjusting head circumference for covariates in autism: clinical correlates of a highly heritable continuous trait

BACKGROUND

Brain development follows a different trajectory in children with autism spectrum disorders (ASD) than in typically developing children. A proxy for neurodevelopment could be head circumference (HC), but studies assessing HC and its clinical correlates in ASD have been inconsistent. This study investigates HC and clinical correlates in the Simons Simplex Collection cohort.

METHODS

We used a mixed linear model to estimate effects of covariates and the deviation from the expected HC given parental HC (genetic deviation). After excluding individuals with incomplete data, 7225 individuals in 1891 families remained for analysis. We examined the relationship between HC/genetic deviation of HC and clinical parameters.

RESULTS

Gender, age, height, weight, genetic ancestry, and ASD status were significant predictors of HC (estimate of the ASD effect = .2 cm). HC was approximately normally distributed in probands and unaffected relatives, with only a few outliers. Genetic deviation of HC was also normally distributed, consistent with a random sampling of parental genes. Whereas larger HC than expected was associated with ASD symptom severity and regression, IQ decreased with the absolute value of the genetic deviation of HC.

CONCLUSIONS

Measured against expected values derived from covariates of ASD subjects, statistical outliers for HC were uncommon. HC is a strongly heritable trait, and population norms for HC would be far more accurate if covariates including genetic ancestry, height, and age were taken into account. The association of diminishing IQ with absolute deviation from predicted HC values suggests HC could reflect subtle underlying brain development and warrants further investigation.

Using large clinical data sets to infer pathogenicity for rare copy number variants in autism cohorts

Copy number variants (CNVs) have a major role in the etiology of autism spectrum disorders (ASD), and several of these have reached statistical significance in case-control analyses. Nevertheless, current ASD cohorts are not large enough to detect very rare CNVs that may be causative or contributory (that is, risk alleles). Here, we use a tiered approach, in which clinically significant CNVs are first identified in large clinical cohorts of neurodevelopmental disorders (including but not specific to ASD), after which these CNVs are then systematically identified within well-characterized ASD cohorts. We focused our initial analysis on 48 recurrent CNVs (segmental duplication-mediated 'hotspots') from 24 loci in 31 516 published clinical cases with neurodevelopmental disorders and 13 696 published controls, which yielded a total of 19 deletion CNVs and 11 duplication CNVs that reached statistical significance. We then investigated the overlap of these 30 CNVs in a combined sample of 3955 well-characterized ASD cases from three published studies. We identified 73 deleterious recurrent CNVs, including 36 deletions from 11 loci and 37 duplications from seven loci, for a frequency of 1 in 54; had we considered the ASD cohorts alone, only 58 CNVs from eight loci (24 deletions from three loci and 34 duplications from five loci) would have reached statistical significance. In conclusion, until there are sufficiently large ASD research cohorts with enough power to detect very rare causative or contributory CNVs, data from larger clinical cohorts can be used to infer the likely clinical significance of CNVs in ASD.

Complete haplotype sequence of the human immunoglobulin heavy-chain variable, diversity, and joining genes and characterization of allelic and copy-number variation

The immunoglobulin heavy-chain locus (IGH) encodes variable (IGHV), diversity (IGHD), joining (IGHJ), and constant (IGHC) genes and is responsible for antibody heavy-chain biosynthesis, which is vital to the adaptive immune response. Programmed V-(D)-J somatic rearrangement and the complex duplicated nature of the locus have impeded attempts to reconcile its genomic organization based on traditional B-lymphocyte derived genetic material. As a result, sequence descriptions of germline variation within IGHV are lacking, haplotype inference using traditional linkage disequilibrium methods has been difficult, and the human genome reference assembly is missing several expressed IGHV genes. By using a hydatidiform mole BAC clone resource, we present the most complete haplotype of IGHV, IGHD, and IGHJ gene regions derived from a single chromosome, representing an alternate assembly of ∼1 Mbp of high-quality finished sequence. From this we add 101 kbp of previously uncharacterized sequence, including functional IGHV genes, and characterize four large germline copy-number variants (CNVs). In addition to this germline reference, we identify and characterize eight CNV-containing haplotypes from a panel of nine diploid genomes of diverse ethnic origin, discovering previously unmapped IGHV genes and an additional 121 kbp of insertion sequence. We genotype four of these CNVs by using PCR in 425 individuals from nine human populations. We find that all four are highly polymorphic and show considerable evidence of stratification (Fst = 0.3-0.5), with the greatest differences observed between African and Asian populations. These CNVs exhibit weak linkage disequilibrium with SNPs from two commercial arrays in most of the populations tested.

Common genetic variants, acting additively, are a major source of risk for autism

Background

Autism spectrum disorders (ASD) are early onset neurodevelopmental syndromes typified by impairments in reciprocal social interaction and communication, accompanied by restricted and repetitive behaviors. While rare and especially de novo genetic variation are known to affect liability, whether common genetic polymorphism plays a substantial role is an open question and the relative contribution of genes and environment is contentious. It is probable that the relative contributions of rare and common variation, as well as environment, differs between ASD families having only a single affected individual (simplex) versus multiplex families who have two or more affected individuals.

Methods

By using quantitative genetics techniques and the contrast of ASD subjects to controls, we estimate what portion of liability can be explained by additive genetic effects, known as narrow-sense heritability. We evaluate relatives of ASD subjects using the same methods to evaluate the assumptions of the additive model and partition families by simplex/multiplex status to determine how heritability changes with status.

Results

By analyzing common variation throughout the genome, we show that common genetic polymorphism exerts substantial additive genetic effects on ASD liability and that simplex/multiplex family status has an impact on the identified composition of that risk. As a fraction of the total variation in liability, the estimated narrow-sense heritability exceeds 60% for ASD individuals from multiplex families and is approximately 40% for simplex families. By analyzing parents, unaffected siblings and alleles not transmitted from parents to their affected children, we conclude that the data for simplex ASD families follow the expectation for additive models closely. The data from multiplex families deviate somewhat from an additive model, possibly due to parental assortative mating.

Conclusions

Our results, when viewed in the context of results from genome-wide association studies, demonstrate that a myriad of common variants of very small effect impacts ASD liability.

De novo mutations revealed by whole-exome sequencing are strongly associated with autism

Multiple studies have confirmed the contribution of rare de novo copy number variations to the risk for autism spectrum disorders. But whereas de novo single nucleotide variants have been identified in affected individuals, their contribution to risk has yet to be clarified. Specifically, the frequency and distribution of these mutations have not been well characterized in matched unaffected controls, and such data are vital to the interpretation of de novo coding mutations observed in probands. Here we show, using whole-exome sequencing of 928 individuals, including 200 phenotypically discordant sibling pairs, that highly disruptive (nonsense and splice-site) de novo mutations in brain-expressed genes are associated with autism spectrum disorders and carry large effects. On the basis of mutation rates in unaffected individuals, we demonstrate that multiple independent de novo single nucleotide variants in the same gene among unrelated probands reliably identifies risk alleles, providing a clear path forward for gene discovery. Among a total of 279 identified de novo coding mutations, there is a single instance in probands, and none in siblings, in which two independent nonsense variants disrupt the same gene, SCN2A (sodium channel, voltage-gated, type II, α subunit), a result that is highly unlikely by chance.