Allelic contribution of Nrxn1α to autism-relevant behavioral phenotypes in mice.

Copy number variations (CNVs) in the Neurexin 1 (NRXN1) gene, which encodes a presynaptic protein involved in neurotransmitter release, are some of the most frequently observed single-gene variants associated with autism spectrum disorder (ASD). To address the functional contribution of NRXN1 CNVs to behavioral phenotypes relevant to ASD, we carried out systematic behavioral phenotyping of an allelic series of Nrxn1 mouse models: one carrying promoter and exon 1 deletion abolishing Nrxn1α transcription, one carrying exon 9 deletion disrupting Nrxn1α protein translation, and one carrying an intronic deletion with no observable effect on Nrxn1α expression. We found that homozygous loss of Nrxn1α resulted in enhanced aggression in males, reduced affiliative social behaviors in females, and significantly altered circadian activities in both sexes. Heterozygous or homozygous loss of Nrxn1α affected the preference for social novelty in male mice, and notably, enhanced repetitive motor skills and motor coordination in both sexes. In contrast, mice bearing an intronic deletion of Nrxn1 did not display alterations in any of the behaviors assessed. These findings demonstrate the importance of Nrxn1α gene dosage in regulating social, circadian, and motor functions, and the variables of sex and genomic positioning of CNVs in the expression of autism-related phenotypes. Importantly, mice with heterozygous loss of Nrxn1, as found in numerous autistic individuals, show an elevated propensity to manifest autism-related phenotypes, supporting the use of models with this genomic architecture to study ASD etiology and assess additional genetic variants associated with autism.

Genome-wide significant risk loci for mood disorders in the Old Order Amish founder population.

Genome-wide association studies (GWAS) of mood disorders in large case-control cohorts have identified numerous risk loci, yet pathophysiological mechanisms remain elusive, primarily due to the very small effects of common variants. We sought to discover risk variants with larger effects by conducting a genome-wide association study of mood disorders in a founder population, the Old Order Amish (OOA, n = 1,672). Our analysis revealed four genome-wide significant risk loci, all of which were associated with >2-fold relative risk. Quantitative behavioral and neurocognitive assessments (n = 314) revealed effects of risk variants on sub-clinical depressive symptoms and information processing speed. Network analysis suggested that OOA-specific risk loci harbor novel risk-associated genes that interact with known neuropsychiatry-associated genes via gene interaction networks. Annotation of the variants at these risk loci revealed population-enriched, non-synonymous variants in two genes encoding neurodevelopmental transcription factors, CUX1 and CNOT1. Our findings provide insight into the genetic architecture of mood disorders and a substrate for mechanistic and clinical studies.

Ancestral haplotype reconstruction in endogamous populations using identity-by-descent.

In this work we develop a novel algorithm for reconstructing the genomes of ancestral individuals, given genotype or sequence data from contemporary individuals and an extended pedigree of family relationships. A pedigree with complete genomes for every individual enables the study of allele frequency dynamics and haplotype diversity across generations, including deviations from neutrality such as transmission distortion. When studying heritable diseases, ancestral haplotypes can be used to augment genome-wide association studies and track disease inheritance patterns. The building blocks of our reconstruction algorithm are segments of Identity-By-Descent (IBD) shared between two or more genotyped individuals. The method alternates between identifying a source for each IBD segment and assembling IBD segments placed within each ancestral individual. Unlike previous approaches, our method is able to accommodate complex pedigree structures with hundreds of individuals genotyped at millions of SNPs. We apply our method to an Old Order Amish pedigree from Lancaster, Pennsylvania, whose founders came to North America from Europe during the early 18th century. The pedigree includes 1338 individuals from the past 12 generations, 394 with genotype data. The motivation for reconstruction is to understand the genetic basis of diseases segregating in the family through tracking haplotype transmission over time. Using our algorithm thread, we are able to reconstruct an average of 224 ancestral individuals per chromosome. For these ancestral individuals, on average we reconstruct 79% of their haplotypes. We also identify a region on chromosome 16 that is difficult to reconstruct-we find that this region harbors a short Amish-specific copy number variation and the gene HYDIN. thread was developed for endogamous populations, but can be applied to any extensive pedigree with the recent generations genotyped. We anticipate that this type of practical ancestral reconstruction will become more common and necessary to understand rare and complex heritable diseases in extended families.

Rediscovering the value of families for psychiatric genetics research.

As it is likely that both common and rare genetic variation are important for complex disease risk, studies that examine the full range of the allelic frequency distribution should be utilized to dissect the genetic influences on mental illness. The rate limiting factor for inferring an association between a variant and a phenotype is inevitably the total number of copies of the minor allele captured in the studied sample. For rare variation, with minor allele frequencies of 0.5% or less, very large samples of unrelated individuals are necessary to unambiguously associate a locus with an illness. Unfortunately, such large samples are often cost prohibitive. However, by using alternative analytic strategies and studying related individuals, particularly those from large multiplex families, it is possible to reduce the required sample size while maintaining statistical power. We contend that using whole genome sequence (WGS) in extended pedigrees provides a cost-effective strategy for psychiatric gene mapping that complements common variant approaches and WGS in unrelated individuals. This was our impetus for forming the "Pedigree-Based Whole Genome Sequencing of Affective and Psychotic Disorders" consortium. In this review, we provide a rationale for the use of WGS with pedigrees in modern psychiatric genetics research. We begin with a focused review of the current literature, followed by a short history of family-based research in psychiatry. Next, we describe several advantages of pedigrees for WGS research, including power estimates, methods for studying the environment, and endophenotypes. We conclude with a brief description of our consortium and its goals.

Functional significance of rare neuroligin 1 variants found in autism.

Genetic mutations contribute to the etiology of autism spectrum disorder (ASD), a common, heterogeneous neurodevelopmental disorder characterized by impairments in social interaction, communication, and repetitive and restricted patterns of behavior. Since neuroligin3 (NLGN3), a cell adhesion molecule at the neuronal synapse, was first identified as a risk gene for ASD, several additional variants in NLGN3 and NLGN4 were found in ASD patients. Moreover, synaptopathies are now known to cause several neuropsychiatric disorders including ASD. In humans, NLGNs consist of five family members, and neuroligin1 (NLGN1) is a major component forming a complex on excitatory glutamatergic synapses. However, the significance of NLGN1 in neuropsychiatric disorders remains unknown. Here, we systematically examine five missense variants of NLGN1 that were detected in ASD patients, and show molecular and cellular alterations caused by these variants. We show that a novel NLGN1 Pro89Leu (P89L) missense variant found in two ASD siblings leads to changes in cellular localization, protein degradation, and to the impairment of spine formation. Furthermore, we generated the knock-in P89L mice, and we show that the P89L heterozygote mice display abnormal social behavior, a core feature of ASD. These results, for the first time, implicate rare variants in NLGN1 as functionally significant and support that the NLGN synaptic pathway is of importance in the etiology of neuropsychiatric disorders.

Correction: Functional significance of rare neuroligin 1 variants found in autism.

[This corrects the article DOI: 10.1371/journal.pgen.1006940.].

Rare variants in the genetic background modulate cognitive and developmental phenotypes in individuals carrying disease-associated variants.

PURPOSE: To assess the contribution of rare variants in the genetic background toward variability of neurodevelopmental phenotypes in individuals with rare copy-number variants (CNVs) and gene-disruptive variants. METHODS: We analyzed quantitative clinical information, exome sequencing, and microarray data from 757 probands and 233 parents and siblings who carry disease-associated variants. RESULTS: The number of rare likely deleterious variants in functionally intolerant genes ("other hits") correlated with expression of neurodevelopmental phenotypes in probands with 16p12.1 deletion (n=23, p=0.004) and in autism probands carrying gene-disruptive variants (n=184, p=0.03) compared with their carrier family members. Probands with 16p12.1 deletion and a strong family history presented more severe clinical features (p=0.04) and higher burden of other hits compared with those with mild/no family history (p=0.001). The number of other hits also correlated with severity of cognitive impairment in probands carrying pathogenic CNVs (n=53) or de novo pathogenic variants in disease genes (n=290), and negatively correlated with head size among 80 probands with 16p11.2 deletion. These co-occurring hits involved known disease-associated genes such as SETD5, AUTS2, and NRXN1, and were enriched for cellular and developmental processes. CONCLUSION: Accurate genetic diagnosis of complex disorders will require complete evaluation of the genetic background even after a candidate disease-associated variant is identified.

Increased burden of deleterious variants in essential genes in autism spectrum disorder.

Autism spectrum disorder (ASD) is a heterogeneous, highly heritable neurodevelopmental syndrome characterized by impaired social interaction, communication, and repetitive behavior. It is estimated that hundreds of genes contribute to ASD. We asked if genes with a strong effect on survival and fitness contribute to ASD risk. Human orthologs of genes with an essential role in pre- and postnatal development in the mouse [essential genes (EGs)] are enriched for disease genes and under strong purifying selection relative to human orthologs of mouse genes with a known nonlethal phenotype [nonessential genes (NEGs)]. This intolerance to deleterious mutations, commonly observed haploinsufficiency, and the importance of EGs in development suggest a possible cumulative effect of deleterious variants in EGs on complex neurodevelopmental disorders. With a comprehensive catalog of 3,915 mammalian EGs, we provide compelling evidence for a stronger contribution of EGs to ASD risk compared with NEGs. By examining the exonic de novo and inherited variants from 1,781 ASD quartet families, we show a significantly higher burden of damaging mutations in EGs in ASD probands compared with their non-ASD siblings. The analysis of EGs in the developing brain identified clusters of coexpressed EGs implicated in ASD. Finally, we suggest a high-priority list of 29 EGs with potential ASD risk as targets for future functional and behavioral studies. Overall, we show that large-scale studies of gene function in model organisms provide a powerful approach for prioritization of genes and pathogenic variants identified by sequencing studies of human disease.

Genomic view of bipolar disorder revealed by whole genome sequencing in a genetic isolate.

Bipolar disorder is a common, heritable mental illness characterized by recurrent episodes of mania and depression. Despite considerable effort to elucidate the genetic underpinnings of bipolar disorder, causative genetic risk factors remain elusive. We conducted a comprehensive genomic analysis of bipolar disorder in a large Old Order Amish pedigree. Microsatellite genotypes and high-density SNP-array genotypes of 388 family members were combined with whole genome sequence data for 50 of these subjects, comprising 18 parent-child trios. This study design permitted evaluation of candidate variants within the context of haplotype structure by resolving the phase in sequenced parent-child trios and by imputation of variants into multiple unsequenced siblings. Non-parametric and parametric linkage analysis of the entire pedigree as well as on smaller clusters of families identified several nominally significant linkage peaks, each of which included dozens of predicted deleterious variants. Close inspection of exonic and regulatory variants in genes under the linkage peaks using family-based association tests revealed additional credible candidate genes for functional studies and further replication in population-based cohorts. However, despite the in-depth genomic characterization of this unique, large and multigenerational pedigree from a genetic isolate, there was no convergence of evidence implicating a particular set of risk loci or common pathways. The striking haplotype and locus heterogeneity we observed has profound implications for the design of studies of bipolar and other related disorders.

A population-based study of KCNH7 p.Arg394His and bipolar spectrum disorder.

We conducted blinded psychiatric assessments of 26 Amish subjects (52 ± 11 years) from four families with prevalent bipolar spectrum disorder, identified 10 potentially pathogenic alleles by exome sequencing, tested association of these alleles with clinical diagnoses in the larger Amish Study of Major Affective Disorder (ASMAD) cohort, and studied mutant potassium channels in neurons. Fourteen of 26 Amish had bipolar spectrum disorder. The only candidate allele shared among them was rs78247304, a non-synonymous variant of KCNH7 (c.1181G>A, p.Arg394His). KCNH7 c.1181G>A and nine other potentially pathogenic variants were subsequently tested within the ASMAD cohort, which consisted of 340 subjects grouped into controls subjects and affected subjects from overlapping clinical categories (bipolar 1 disorder, bipolar spectrum disorder and any major affective disorder). KCNH7 c.1181G>A had the highest enrichment among individuals with bipolar spectrum disorder (χ(2) = 7.3) and the strongest family-based association with bipolar 1 (P = 0.021), bipolar spectrum (P = 0.031) and any major affective disorder (P = 0.016). In vitro, the p.Arg394His substitution allowed normal expression, trafficking, assembly and localization of HERG3/Kv11.3 channels, but altered the steady-state voltage dependence and kinetics of activation in neuronal cells. Although our genome-wide statistical results do not alone prove association, cumulative evidence from multiple independent sources (parallel genome-wide study cohorts, pharmacological studies of HERG-type potassium channels, electrophysiological data) implicates neuronal HERG3/Kv11.3 potassium channels in the pathophysiology of bipolar spectrum disorder. Such a finding, if corroborated by future studies, has implications for mental health services among the Amish, as well as development of drugs that specifically target HERG3/Kv11.3.

From mouse to human: evolutionary genomics analysis of human orthologs of essential genes.

Understanding the core set of genes that are necessary for basic developmental functions is one of the central goals in biology. Studies in model organisms identified a significant fraction of essential genes through the analysis of null-mutations that lead to lethality. Recent large-scale next-generation sequencing efforts have provided unprecedented data on genetic variation in human. However, evolutionary and genomic characteristics of human essential genes have never been directly studied on a genome-wide scale. Here we use detailed phenotypic resources available for the mouse and deep genomics sequencing data from human populations to characterize patterns of genetic variation and mutational burden in a set of 2,472 human orthologs of known essential genes in the mouse. Consistent with the action of strong, purifying selection, these genes exhibit comparatively reduced levels of sequence variation, skew in allele frequency towards more rare, and exhibit increased conservation across the primate and rodent lineages relative to the remainder of genes in the genome. In individual genomes we observed ~12 rare mutations within essential genes predicted to be damaging. Consistent with the hypothesis that mutations in essential genes are risk factors for neurodevelopmental disease, we show that de novo variants in patients with Autism Spectrum Disorder are more likely to occur in this collection of genes. While incomplete, our set of human orthologs shows characteristics fully consistent with essential function in human and thus provides a resource to inform and facilitate interpretation of sequence data in studies of human disease.

Mapping autism risk loci using genetic linkage and chromosomal rearrangements.

Autism spectrum disorders (ASDs) are common, heritable neurodevelopmental conditions. The genetic architecture of ASDs is complex, requiring large samples to overcome heterogeneity. Here we broaden coverage and sample size relative to other studies of ASDs by using Affymetrix 10K SNP arrays and 1,181 [corrected] families with at least two affected individuals, performing the largest linkage scan to date while also analyzing copy number variation in these families. Linkage and copy number variation analyses implicate chromosome 11p12-p13 and neurexins, respectively, among other candidate loci. Neurexins team with previously implicated neuroligins for glutamatergic synaptogenesis, highlighting glutamate-related genes as promising candidates for contributing to ASDs.

Copy number variants encompassing Mendelian disease genes in a large multigenerational family segregating bipolar disorder.

BACKGROUND: Bipolar affective disorder (BP) is a common, highly heritable psychiatric disorder characterized by periods of depression and mania. Using dense SNP genotype data, we characterized CNVs in 388 members of an Old Order Amish Pedigree with bipolar disorder. We identified CNV regions arising from common ancestral mutations by utilizing the pedigree information. By combining this analysis with whole genome sequence data in the same individuals, we also explored the role of compound heterozygosity. RESULTS: Here we describe 541 inherited CNV regions, of which 268 are rare in a control population of European origin but present in a large number of Amish individuals. In addition, we highlight a set of CNVs found at higher frequencies in BP individuals, and within genes known to play a role in human development and disease. As in prior reports, we find no evidence for an increased burden of CNVs in BP individuals, but we report a trend towards a higher burden of CNVs in known Mendelian disease loci in bipolar individuals (BPI and BPII, p = 0.06). CONCLUSIONS: We conclude that CNVs may be contributing factors in the phenotypic presentation of mood disorders and co-morbid medical conditions in this family. These results reinforce the hypothesis of a complex genetic architecture underlying BP disorder, and suggest that the role of CNVs should continue to be investigated in BP data sets.

Common genetic variants on 5p14.1 associate with autism spectrum disorders.

Autism spectrum disorders (ASDs) represent a group of childhood neurodevelopmental and neuropsychiatric disorders characterized by deficits in verbal communication, impairment of social interaction, and restricted and repetitive patterns of interests and behaviour. To identify common genetic risk factors underlying ASDs, here we present the results of genome-wide association studies on a cohort of 780 families (3,101 subjects) with affected children, and a second cohort of 1,204 affected subjects and 6,491 control subjects, all of whom were of European ancestry. Six single nucleotide polymorphisms between cadherin 10 (CDH10) and cadherin 9 (CDH9)-two genes encoding neuronal cell-adhesion molecules-revealed strong association signals, with the most significant SNP being rs4307059 (P = 3.4 x 10(-8), odds ratio = 1.19). These signals were replicated in two independent cohorts, with combined P values ranging from 7.4 x 10(-8) to 2.1 x 10(-10). Our results implicate neuronal cell-adhesion molecules in the pathogenesis of ASDs, and represent, to our knowledge, the first demonstration of genome-wide significant association of common variants with susceptibility to ASDs.

Nuclear receptor corepressor and histone deacetylase 3 govern circadian metabolic physiology.

Rhythmic changes in histone acetylation at circadian clock genes suggest that temporal modulation of gene expression is regulated by chromatin modifications. Furthermore, recent studies demonstrate a critical relationship between circadian and metabolic physiology. The nuclear receptor corepressor 1 (Ncor1) functions as an activating subunit for the chromatin modifying enzyme histone deacetylase 3 (Hdac3). Lack of Ncor1 is incompatible with life, and hence it is unknown whether Ncor1, and particularly its regulation of Hdac3, is critical for adult mammalian physiology. Here we show that specific, genetic disruption of the Ncor1-Hdac3 interaction in mice causes aberrant regulation of clock genes and results in abnormal circadian behaviour. These mice are also leaner and more insulin-sensitive owing to increased energy expenditure. Unexpectedly, loss of a functional Ncor1-Hdac3 complex in vivo does not lead to sustained increases in known catabolic genes, but instead significantly alters the oscillatory patterns of several metabolic genes, demonstrating that circadian regulation of metabolism is critical for normal energy balance. These findings indicate that activation of Hdac3 by Ncor1 is a nodal point in the epigenetic regulation of circadian and metabolic physiology.

Genotype, haplotype and copy-number variation in worldwide human populations.

Genome-wide patterns of variation across individuals provide a powerful source of data for uncovering the history of migration, range expansion, and adaptation of the human species. However, high-resolution surveys of variation in genotype, haplotype and copy number have generally focused on a small number of population groups. Here we report the analysis of high-quality genotypes at 525,910 single-nucleotide polymorphisms (SNPs) and 396 copy-number-variable loci in a worldwide sample of 29 populations. Analysis of SNP genotypes yields strongly supported fine-scale inferences about population structure. Increasing linkage disequilibrium is observed with increasing geographic distance from Africa, as expected under a serial founder effect for the out-of-Africa spread of human populations. New approaches for haplotype analysis produce inferences about population structure that complement results based on unphased SNPs. Despite a difference from SNPs in the frequency spectrum of the copy-number variants (CNVs) detected--including a comparatively large number of CNVs in previously unexamined populations from Oceania and the Americas--the global distribution of CNVs largely accords with population structure analyses for SNP data sets of similar size. Our results produce new inferences about inter-population variation, support the utility of CNVs in human population-genetic research, and serve as a genomic resource for human-genetic studies in diverse worldwide populations.

Phenotypic and ancestry-related assortative mating in autism.

BACKGROUND: Positive assortative mating (AM) in several neuropsychiatric traits, including autism, has been noted. However, it is unknown whether the pattern of AM is different in phenotypically defined autism subgroups [e.g., autism with and without intellectually disability (ID)]. It is also unclear what proportion of the phenotypic AM can be explained by the genetic similarity between parents of children with an autism diagnosis, and the consequences of AM on the genetic structure of the population. METHODS: To address these questions, we analyzed two family-based autism collections: the Simons Foundation Powering Autism Research for Knowledge (SPARK) (1575 families) and the Simons Simplex Collection (SSC) (2283 families). RESULTS: We found a similar degree of phenotypic and ancestry-related AM in parents of children with an autism diagnosis regardless of the presence of ID. We did not find evidence of AM for autism based on autism polygenic scores (PGS) (at a threshold of |r|> 0.1). The adjustment of ancestry-related AM or autism PGS accounted for only 0.3-4% of the fractional change in the estimate of the phenotypic AM. The ancestry-related AM introduced higher long-range linkage disequilibrium (LD) between single nucleotide polymorphisms (SNPs) on different chromosomes that are highly ancestry-informative compared to SNPs that are less ancestry-informative (D2 on the order of 1 × 10-5). LIMITATIONS: We only analyzed participants of European ancestry, limiting the generalizability of our results to individuals of non-European ancestry. SPARK and SSC were both multicenter studies. Therefore, there could be ancestry-related AM in SPARK and SSC due to geographic stratification. The study participants from each site were unknown, so we were unable to evaluate for geographic stratification. CONCLUSIONS: This study showed similar patterns of AM in autism with and without ID, and demonstrated that the common genetic influences of autism are likely relevant to both autism groups. The adjustment of ancestry-related AM and autism PGS accounted for < 5% of the fractional change in the estimate of the phenotypic AM. Future studies are needed to evaluate if the small increase of long-range LD induced by ancestry-related AM has impact on the downstream analysis.