PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of tpLink to Publisher's site
 
Transl Psychiatry. Sep 2013; 3(9): e301.
Published online Sep 3, 2013. doi:  10.1038/tp.2013.75
PMCID: PMC3784764
A rare WNT1 missense variant overrepresented in ASD leads to increased Wnt signal pathway activation
P-M Martin,1 X Yang,1 N Robin,2 E Lam,3 J S Rabinowitz,2 C A Erdman,1,3 J Quinn,1,3 L A Weiss,1,3 S P Hamilton,1,3 P-Y Kwok,3 R T Moon,2 and B N R Cheyette1*
1Department of Psychiatry, University of California, San Francisco, San Francisco, CA, USA
2Department of Pharmacology and Howard Hughes Medical Institute, University of Washington School of Medicine, University of Washington, Seattle, WA, USA
3Institute for Human Genetics, University of California San Francisco, San Francisco, CA, USA
*Department of Psychiatry, University of California, Rock Hall Room 284D, 1550 4th Street, San Francisco, CA 94158-2324, USA. E-mail: bc/at/ucsf.edu
Received April 25, 2013; Revised July 23, 2013; Accepted July 24, 2013.
Wnt signaling, which encompasses multiple biochemical pathways that regulate neural development downstream of extracellular Wnt glycoprotein ligands, has been suggested to contribute to major psychiatric disorders including autism spectrum disorders (ASD). We used next-generation sequencing and Sequenom genotyping technologies to resequence 10 Wnt signaling pathway genes in 198 ASD patients and 240 matched controls. Results for single-nucleotide polymorphisms (SNPs) of interest were confirmed in a second set of 91 ASD and 144 control samples. We found a significantly increased burden of extremely rare missense variants predicted to be deleterious by PolyPhen-2, distributed across seven genes in the ASD sample (3.5% in ASD vs 0.8% in controls; Fisher's exact test, odds ratio (OR)=4.37, P=0.04). We also found a missense variant in WNT1 (S88R) that was overrepresented in the ASD sample (8 A/T in 267 ASD (minor allele frequency (MAF)=1.69%) vs 1 A/T in 377 controls (MAF=0.13%), OR=13.0, Fisher's exact test, P=0.0048; OR=8.2 and P=0.053 after correction for population stratification). Functional analysis revealed that WNT1-S88R is more active than wild-type WNT1 in assays for the Wnt/β-catenin signaling pathway. Our findings of a higher burden in ASD of rare missense variants distributed across 7 of 10 Wnt signaling pathway genes tested, and of a functional variant at the WNT1 locus associated with ASD, support that dysfunction of this pathway contributes to ASD susceptibility. Given recent findings of common molecular mechanisms in ASD, schizophrenia and affective disorders, these loci merit scrutiny in other psychiatric conditions as well.
Keywords: association study, functional analysis, mutation burden, psychiatric disorders, rare mutations
Autism spectrum disorders (ASD) have a strong genetic etiology. The last decade of research in psychiatric genetics has led to groundbreaking discoveries, including the identification of ASD-associated sequence and structural variants in a variety of genes implicated in brain development—particularly in synapse formation, maintenance and function. However, the mutations identified so far explain only a fraction of ASD cases, underscoring the need to continue the search for missing genetic links in ASD.
Wnt intercellular signaling refers to a set of key developmental pathways that help to organize the nervous system. Among other roles, Wnt signaling influences subdivision of the anterior neural tube into major brain regions (such as forebrain, midbrain and hindbrain), neural precursor proliferation, neural cell fate and migration, axon guidance, dendrite development and synapse formation.1 Given this broad spectrum of roles, it is no surprise that dysregulation of Wnt signaling causes a similarly broad spectrum of deleterious effects on neural development, and largely on this basis it has long been considered a candidate pathway in psychiatric pathogenesis. Its candidacy in this regard is supported by evidence that Wnt signaling pathway components are pharmacological targets of lithium,2, 3 hallucinogenic drugs 4 and antipsychotic medications,5, 6 as well as by evidence derived from behavioral tests in animal models7, 8, 9, 10, 11 and from positive associations in scattered human genetic studies.12, 13, 14, 15, 16, 17
Here we used massively parallel sequencing technology to simultaneously search for rare variants in 10 Wnt signaling pathway genes in a sample of Caucasian patients with ASD and healthy controls. We selected these 10 genes out of a much larger set of pathway molecules based primarily on empirical support for their individual relevance to psychiatry in human genetic (DKK1, DKK4, NEUROG1, CTNNB1),18, 19, 20, 21 animal behavior (DVL1, DIXDC1),7, 10 biochemical (WNT1, NEUROD1)22, 23 or other types of studies (PRICKLE4, DACT2).24, 25 Compared with whole exome sequencing strategies that are now underway, this very targeted approach benefited from our focus on functionally linked loci for which we have substantial expertise, allowing us to immediately follow up the most promising variant in in vitro functional assays.
ASD samples (Caucasians)
Samples were lymphoblastoid cell lines from the Simons Simplex Collection. The diagnosis in these samples was based on the ADI-R (Autism Diagnostic Interview-Revised).26
Unscreened controls (Caucasians)
Healthy control subjects were recruited by the survey research company (Knowledge Networks, Menlo Park, CA, USA) from a nationally representative internet-based panel that was selected by random digit dialing.27 Subjects completed an online version of the CIDI-SF (Composite International Diagnostic Interview-Short Form) for lifetime history of anxiety, mood and substance use disorders. Subjects consented to anonymization and deposition of their clinical information and DNA in the National Institute of Mental Health (NIMH) repository for use in any medical research.
Pooling samples
Concentration of individual DNA samples was measured using Pico Green. Based on previous reports using an identical pooling strategy,28, 29 we expected sequencing error rates to be ~0.5%. To distinguish a single real mutant allele in 876 chromosomes (minor allele frequency (MAF)=0.11%) from random sequencing artifacts, we pooled samples into 11 sets of 38–40 individuals each. Therefore, the frequency of a single mutant allele in each pool was 1%, significantly above the detection threshold of 0.5%.
Target enrichment
Each pool of genomic DNA from 38–40 individuals was enriched by PCR amplification for 10 genes involved in Wnt signaling (WNT1, DVL1, DIXDC1, PRICKLE4, DACT2, NEUROG1, DKK1, DKK4, NEUROD1, and CTNNB1), plus 1 control gene, DACT1, during the validation step (Supplementary Table S1). For each gene, the entire locus was amplified in 1 or 2 ‘long-range' PCR reactions. The total amplification length for all 11 genes was 71.5 kb.
Sample preparation and high-throughput sequencing procedure
Each pool of target-enriched genomic DNA was used to build a prepped DNA library following Illumina sample preparation protocols. In brief, each pool was ligated with adapters containing a unique index allowing us to combine 11 pools in each sequencing lane of the HiSeq flow cell (HiSeq 2000, Illumina, San Diego, CA, USA). Finished flow cells carrying 438 individuals (1 lane containing 11 indexed pools of 38–40 individuals) were sequenced in the HiSeq 2000 for paired-end analysis.
Initial data processing (sequence assembly, SNP calling, quality assessment)
Bowtie30 was used for read alignment against the latest human genome reference. SAMtools,31 Picard (not published; http://picard.sourceforge.net/), BEDTools32 and the Genome Analysis Toolkit (GATK;33) were used for post-alignment processing. Multisample realignment around potential insertion/deletions (indels) and base quality score recalibration were performed before variant calling. For single-nucleotide polymorphism (SNP) calling, we used SyZyGy, an algorithm specifically designed to call bases in pooled samples.34
Analysis of conservation and prediction as deleterious
For each SNP, analysis by SeattleSeq (http://snp.gs.washington.edu/SeattleSeqAnnotation137/) provided information about base position conservation in the genome, whether it was intronic or exonic, and if exonic whether it was synonymous, missense or nonsense. Predicted deleteriousness for each SNP in exonic sequence was determined by PolyPhen-235 version v2.1.0r367 or more recent. If either version predicted that the SNP was deleterious, it was counted as deleterious.
Human cDNAs and mutagenesis
Human WNT1 (clone ID: 30915309) complementary DNA (cDNA) was purchased from Thermo Scientific (Rockford, IL, USA) and subcloned into a pcDNA 3.1(−) vector. The WNT1 S88R variant was generated using the QuikChange site-directed mutagenesis kit (Stratagene, La Jolla, CA, USA) and confirmed by Sanger sequencing. In subsequent functional assays, in order to avoid potential confounds arising from nonsystemic biases in measurements of plasmid yield or purity, each construct was prepared independently >2 times each by different experimenters in two different laboratories. Similar functional results were obtained from each independent preparation compared with wild type (WT).
Luciferase reporter assays
This was performed essentially as described previously.36 Briefly, both HEK293T and SH-SY5Y BAR-Luc cells were transfected with 1, 5, 20 and 50 ng pcDNA3.1(−) vector containing either WT WNT1 or WNT1-S88R. HEK293T cells were additionally transfected with 20 ng per well pBAR-Luc plasmid and 10 ng per well pRL-TK (Renilla luciferase control plasmid). Cells were cultured for 24 h, and then luciferase activity measured using the Dual-Luciferase Reporter Assay System (Promega, Madison, WI, USA) and an Envision multi-label plate reader (PerkinElmer, Waltham, MA, USA).
Quantitative reverse transcription-PCR
Total RNA from cultured HEK293T cells transfected with an empty pcDNA3.1(−) vector, or a pcDNA3.1(−) vector containing either WT WNT1 or WNT1-S88R, was isolated with the RNeasy Mini kit according to the manufacturer's instructions (Qiagen, Hilden, Germany). The cDNA was synthesized from 2 μg of total RNA with Superscript III Reverse Transcriptase (Invitrogen, Carlsbad, CA, USA) using random hexamer primers. Quantitative reverse transcription PCR was performed using the PerfeCTa SYBR Green FastMix, ROX (Quanta, Byfleet, UK) with gene-specific primer pairs (WNT1 forward (Fwd): 5′-GCAGCGACAACATTGACTTC-3′, WNT1 reverse (Rev): 5′-GTGGCACTTGCACTCCTG-3′ AXIN2 Fwd: 5′-TTATGCTTTGCACTACGTCCCTCCA-3′ AXIN2 Rev: 5′-CGCAACATGGTCAACCCTCAGAC-3′) on an ABI 7900HT real-time PCR system (ABI Advanced Technologies, Columbia, MD, USA). Data were analyzed using the comparative CT method (User Bulletin No. 2, PerkinElmer Life Sciences, Boston, MA, USA). To normalize for loaded cDNA, glyceraldehyde 3-phosphate dehydrogenase (GAPDH) was used as an endogenous control (GAPDH Fwd: 5′-CAAGGTCATCCATGACAACTTTG-3′, GAPDH Rev: 5′-GGCCATCCACAGTCTTCTGG-3′). Results shown reflect combined data from two independent experiments, each done in triplicate reactions.
Statistical analyses
Unless stated otherwise in Figure legends, the data are expressed as mean±s.e.m. Data were analyzed by Student's t-test, Fisher's exact test or Pearson's correlation test where appropriate, using GraphPad Prism software (GraphPad Software, La Jolla, CA, USA). PLINK software37 was used for Multidimensional Scaling and logistic regression. A P-value of <0.05 was considered to be significant.
Using a next-generation (massively parallel Illumina-based) strategy, we sequenced promoter, 3' untranslated, intronic and coding regions of 10 Wnt pathway genes (WNT1, DVL1, DIXDC1, PRICKLE4, DACT2, NEUROG1, DKK1, DKK4, NEUROD1 and CTNNB1) in genomic DNA obtained from a first set of 198 Caucasian ASD patients and 240 healthy Caucasian controls, arranged in pools of 38–40 samples each. The number of reads was homogeneously distributed over the different amplicons and pools (and therefore between ASD and controls), reads mapped 88% of the targeted regions and coverage averaged 865-fold/individual. To validate the next-generation sequencing strategy employed, we compared results obtained using this strategy to results obtained using Sanger sequencing of the exons in a control locus (DACT1) in a subset of samples (40 ASD). Every DACT1 variant (n=8), including several singleton variants (n=3), found by Sanger sequencing in this set was also detected by the next-generation sequencing strategy, and the allele frequencies of all SNPs were identical using each methodology (Pearson's correlation test; r=0.99; P<0.0001; Supplementary Table S1).
In the 10 genes under study, we identified 652 SNPs overall in this set of 438 combined ASD and control samples. We restricted most of our attention to nonsynonymous sequence variants found only in our own sample set (not referenced in dbSNP, 1000 Genomes, the Exome Variant Server (EVS) or the Exome Chip Design Consortium). To determine whether any of the putative SNPs found by this method were artifacts produced by the next-generation sequencing technology, we confirmed each of them by individual genotyping using either IPlex Sequenom or Sanger sequencing. Using this approach we confirmed 12 novel singleton missense variants distributed across 8 of the 10 loci included in our study (Supplementarty Table S2). Moreover, genotyping in a second set of 91 ASD and 144 control samples failed to identify any additional occurrences of each of these rare deleterious missense SNPs. Although each of these SNPs were found only once among the probands in our samples and had not been reported in any public database, genotyping of parents demonstrated that none of them had occurred de novo; all were transmitted from either the mother or father. Moreover, sequencing of family trios showed that each was also present in an unaffected sibling. This verifies that each of these variants is not an idiosyncratic artifact arising during lymphoblastoid cell line derivation, and suggests that they are extremely rare alleles (single-nucleotide variants) in the Caucasian population represented by our samples.
Although we found a greater number of extremely rare missense variants in our ASD sample compared with our control sample (4% in ASD vs 1.7% in controls; Figure 1), this difference in the ‘rare missense variant burden' did not reach statistical significance (Fisher's exact test, P=0.11). Nonetheless, when considering only those rare missense variants that were also predicted to be deleterious by PolyPhen-2, the difference in ASD versus control groups did become significant (3.5% in ASD vs 0.8% in controls; Fisher's exact test, OR=4.37, P=0.04; Figure 1). We also examined the distribution of less informative variants (that is, synonymous, untranslated regions and intronic variants) in each group, and found a very similar distribution for all of them (synonymous: 3.5% in ASD vs 2.9% in controls; untranslated regions: 11.1% in ASD vs 9.6% in controls; intronic: 24% in ASD vs 21% in controls; Figure 1). The similarity in the incidence of less informative variants argues against the possibility that the difference in the incidence of rare deleterious variants can be attributed to differences in recent ancestry between the ASD and control sample groups.
Figure 1
Figure 1
Burden of novel unique variants in the Wnt pathway in autism spectrum disorders (ASD). Percentage of individuals carrying a novel unique variant that was nonsynonymous and predicted to be deleterious by the PolyPhen-2 algorithm, compared with nonsynonymous, (more ...)
In addition to these novel SNPs, we found one previously referenced rare (<1%) missense SNP that was overrepresented in our ASD group versus our control group. This SNP (rs61758378) occurred in WNT1, the gene encoding the first-discovered mammalian Wnt ligand. To eliminate the possibility of a sequencing artifact we confirmed this result by genotyping all samples using IPlex Sequenom. We also genotyped an additional 91 ASD and 144 control samples to replicate and potentially extend these findings. In this way, we not only confirmed MAF results for rs61758378 obtained using the Illumina platform in the original sample set, we further found that the rs61758378 minor allele remained significantly overrepresented in the combined ASD sample compared with the combined control sample (8 A/T in 267 ASD (MAF=1.69%) vs 1 A/T in 377 controls (MAF=0.13%); Fisher's exact test, allelic P=0.0048; Table 1). Furthermore, when we compared the rs61758378 MAF in our combined ASD sample with the reported MAF in the EVS for a much larger group of individuals, the association also remained significant (8 A/T in 267 in our combined ASD sample (MAF=1.69%) vs 45 A/T in 4300 European/American individuals in the EVS (MAF=0.5%), Fisher's exact test, allelic P=0.011). Because different allele frequencies are reported in European populations (TSI, Tuscans in Italy, MAF=3.0% CEU, Northern Europeans in Utah, MAF=0.5%), we wanted to evaluate and address possible effects of subtle stratification in ancestry between cases and controls. We therefore performed Multidimensional Scaling using genome-wide SNP genotype data for a subset of the total samples used in this analysis (235 SSCs and 329 controls) using PLINK. Multidimensional Scaling using reference CEU and TSI data showed that the ratio of individuals clustering with CEU and TSI was not significantly different between ASD and controls. Nevertheless, we noted that one of the first four dimensions was significantly different between ASD and controls, indicating possible stratification. We sought to account for this using a logistic regression analysis including the first four dimensions as covariates. In this analysis, the odds ratio (OR) remained high (OR=8.17) and the P-value remained close to significance (P=0.053). We note that our analyses have low power for SNPs with such a low MAF and in a modest number of samples.
Table 1
Table 1
Association of WNT1 rs61758378 with ASD
The rs61758378 minor allele changes a conserved serine at position 88 into arginine. Although this change is not predicted to be deleterious by PolyPhen-2 (HumDiv score=0.174), serine 88 is nonetheless a very highly conserved residue across species (PhasCons score=0.998), and the change to arginine at this residue is predicted to be deleterious by Grantham score, which categorizes codon replacements into classes of increasing chemical dissimilarity (Grantham score=110, moderately radical change).38 Given the potential association between the WNT1 rs61758378 minor allele and ASD in our sample set, we therefore tested for functional effects of this variant on Wnt signaling pathway activation. To accomplish this, we engineered this variant into a human WNT1 cDNA (WNT1-S88R) using site-directed mutagenesis and then tested it alongside WT Wnt1 in several assays for Wnt signaling activity.
First, we performed a plasmid-based Wnt reporter assay in human embryonic kidney (HEK293T) tissue culture cells. This is a well-established method to quantify transcriptional activation downstream of the Wnt/β-catenin signaling pathway, based in this case on transient transfection of the pBAR-Luc reporter plasmid that expresses luciferase downstream of a β-catenin-responsive promoter.36 Cells were co-transfected with pBAR-Luc and a dose range (1–50 ng) of plasmids encoding WNT1-S88R or WT WNT1 to test for the ability of these proteins to activate Wnt/β-catenin pathway-dependent transcription. At all doses tested, WNT1-S88R showed increased activity in this assay compared with WT WNT1 (Figure 2a). Moreover, this difference was statistically significant at the 20 ng dose (WT WNT1, 1.00±0.09 vs WNT1-S88R, 1.37±0.087, t-test, P=0.02) and at the 50 ng dose (WT WNT1, 1.00±0.18 vs WNT1-S88R, 1.86±0.22, t-test, P=0.016; Figure 2a). Similar results were obtained when experiments were replicated in a different laboratory by a different investigator using independently prepared DNA constructs and reagents (pSuperTop; data not shown).
Figure 2
Figure 2
Functional analysis of S88R. (a) Plasmid-based Wnt/β-catenin pathway reporter assay in HEK293T cells. HEK293T cells were transfected with the pBAR-Luc reporter plasmid together with plasmids encoding either WNT1-S88R or wild-type (WT) WNT1 at (more ...)
The aforementioned assays rely on transient transfection of reporter plasmid, and on this basis they are conceptually more variable than assays based on read-out of an endogenous Wnt pathway target gene or of a single-copy Wnt pathway reporter integrated into the genome. Accordingly, we confirmed these results in the same (HEK293T) cell line using quantitative reverse transcription-PCR to directly measure transcriptional activation of the endogenous Wnt/β-catenin pathway target gene AXIN2.39 Mirroring results obtained with the Wnt reporter assay, compared with WT WNT1, WNT1-S88R showed significantly increased transcriptional activation of AXIN2 (WT Wnt1, 1.00±0.025 vs Wnt1-S88R, 2.16±0.43, t-test, P=0.02; Figure 2b). We confirmed these results via a third in vitro assay that relies on SH-SY5Y BAR-Luc cells. This is an independently derived immortalized cell line with a single copy of a Wnt/β-catenin pathway luciferase reporter stably integrated into its genome.36 Consistent with data obtained by the plasmid-based Wnt reporter assays and by quantitative reverse transcription-PCR in HEK293T cells, SH-SY5Y BAR-Luc cells showed significantly increased activation of their reporter after transfection of WNT1-S88R (WT Wnt1, 1.00±0.04 vs Wnt1-S88R, 1.28±0.03, t-test, P=0.0003; Figure 2c).
Taken together, our data from cell-based transcription target assays suggest that the rs61758378 minor allele associated with ASD in our sample (WNT1-S88R) encodes a modestly more active form of the WNT1 ligand compared with the most prevalent WNT1 allele found in Caucasians.
Our finding in ASD of a higher burden of rare missense variants distributed across 7 of 10 Wnt signaling pathway genes tested, and of an ASD-associated functional variant at the WNT1 locus, supports that dysfunction of this pathway contributes to ASD susceptibility.
In functional assays for the Wnt/β-catenin pathway, the biochemical cascade physiologically activated by the WNT1 gene product,40, 41 WNT1-S88R consistently showed increased activity. In animal models, hyperactivation of Wnt/β-catenin signaling can lead to dominant neurodevelopmental defects including in brain regionalization,42, 43 neural precursor proliferation, cell fate specification, neuron migration,44, 45, 46, 47 and in dendrite and synapse formation or function.48, 49 A mutation in a gene, such as WNT1, that modestly increases activation of this downstream signaling pathway would therefore be expected to result in neurodevelopmental and behavioral phenotypes. Moreover, as ASD and other major psychiatric disorders are generally not associated with gross neuroanatomical defects, changes in Wnt signaling activity contributing to susceptibility for these disorders should be modest. Anatomical changes resulting from a signaling variant such as WNT1-S88R might only be evident at a cellular level (for example, in neural subtype numbers and position, neural morphology, connectivity and so on) or at a subcellular level (for example, in dendritic spine morphology, synapse numbers, synapse strength, synapse plasticity and so on), yet could still be sufficient to cause behavioral symptoms. This of course conforms closely to the central features of major psychiatric disorders, which are distinguished from neurological disorders in part by the absence of focal or grossly evident structural brain defects. Moreover, it is not our intent to posit that this or the other rare coding variants we have identified here in WNT1 and in Wnt signaling pathway component genes are monogenic causes of psychiatric illness.50 Genetic susceptibility to these disorders is highly heterogeneous and thought to be frequently multigenic; that is, for most patients genetic loading is thought to be distributed across many loci that affect a common pathogenic pathway.51, 52 In this susceptibility model, an allelic variant at any one locus may only appreciably increase disease risk in the context of variants at other loci that affect the same pathway in the same individual. The rare nonsynonymous sequence variants found in both affected and unaffected individuals within the same family in this study may contribute to ASD in the context of other genetic and nongenetic factors that increase susceptibility.
Several lines of evidence support that Wnt ligands (of which there are 19 family members in humans) participate in ASD susceptibility. A recent whole exome sequencing study identified copy number variants affecting both WNT3 and WNT9B in an ASD patient,53 whereas another study showed that WNT3 expression is increased in the prefrontal cortex of young autistic patients compared with controls.54 Rare variants in WNT2 have been reported to segregate with ASD.16 Although other groups have not always replicated this linkage,55, 56 two more recent studies have found SNPs in WNT2 associated with either ASD or speech delay.15, 17 In a Valproate-exposure rat model of ASD, expression of both WNT1 and WNT2 was upregulated in the developing forebrain and this was associated with increased Wnt/β-catenin pathway activity in that brain region.57 Other recent sequencing studies have identified several ASD-susceptibility genes that encode functionally linked components of a protein network involved in the regulation of β-catenin-dependent transcription, as occurs in Wnt/β-catenin signaling.21, 58, 59 The findings of the present study add to this evidence by suggesting that minor alleles and extremely rare variants at the WNT1 locus and other Wnt pathway loci may contribute to risk for this psychiatric disorder, a possibility that merits further attention in ongoing human genetic studies focused on ASD. Given evidence that genetic etiologies are shared across several psychiatric disorders,60, 61, 62 as well as that Wnt signaling is involved more broadly in psychiatric pathogenic and therapeutic mechanisms,63 these loci should also be carefully scrutinized in ongoing human genetic studies of other major psychiatric disorders, such as schizophrenia and major affective disorders.
Acknowledgments
P-MM has received support from a Pilot Project award sponsored by Kristine Yaffe, the Roy and Marie Scola Endowed Chair in Psychiatry at UCSF; XYY was funded by NIH T32 MH089920; this research has also been supported by grants to BNRC and SPH together and to BNRC alone from the UCSF Research Allocation Program (RAP) including the Program for Breakthrough Biomedical Research (PBBR), by the Department of Psychiatry, and by the Center for Neurobiology and Psychiatry. LAW has support from the International Mental Health Research Organization. We acknowledge SFARI for ASD genomic DNA samples and for phenotypic data made available on SFARI Base. We also thank all the families who donated samples and the principal investigators involved in their collection (A Beaudet, R Bernier, J Constantino, E Cook, E Fombonne, D Geschwind, D Grice, A Klin, D Ledbetter, C Lord, C Martin, D Martin, R Maxim, J Miles, O Ousley, B Peterson, J Piggot, C Saulnier, M State, W Stone, J Sutcliffe, C Walsh and E Wijsman). Approved researchers can obtain the SSC population data set by applying at https://base.sfari.org. Collection of the control data set was supported by grants from the NIMH and the National Alliance for Research on Schizophrenia and Depression.
Notes
The authors declare no conflict of interest.
Footnotes
Supplementary Information accompanies the paper on the Translational Psychiatry website (http://www.nature.com/tp)
Supplementary Material
Supplementary Table S1
Supplementary Table S2
  • Mulligan KA, Cheyette BN. Wnt signaling in vertebrate neural development and function. J Neuroimmune Pharmacol. 2012;7:774–787. [PMC free article] [PubMed]
  • Wexler EM, Geschwind DH, Palmer TD. Lithium regulates adult hippocampal progenitor development through canonical Wnt pathway activation. Mol Psychiatry. 2008;13:285–292. [PubMed]
  • O'Brien WT, Klein PS. Validating GSK3 as an in vivo target of lithium action. Biochem Soc Trans. 2009;37:1133–1138. [PMC free article] [PubMed]
  • Svenningsson P, Tzavara ET, Carruthers R, Rachleff I, Wattler S, Nehls M, et al. Diverse psychotomimetics act through a common signaling pathway. Science. 2003;302:1412–1415. [PubMed]
  • Freyberg Z, Ferrando SJ, Javitch JA. Roles of the Akt/GSK-3 and Wnt signaling pathways in schizophrenia and antipsychotic drug action. Am J Psychiatry. 2010;167:388–396. [PMC free article] [PubMed]
  • Sutton LP, Honardoust D, Mouyal J, Rajakumar N, Rushlow W. Activation of the canonical Wnt pathway by the antipsychotics haloperidol and clozapine involves dishevelled-3. J Neurochem. 2007;102:153–169. [PubMed]
  • Long JM, LaPorte P, Paylor R, Wynshaw-Boris A. Expanded characterization of the social interaction abnormalities in mice lacking Dvl1. Genes Brain Behav. 2004;3:51–62. [PubMed]
  • Gould TD, O'Donnell KC, Picchini AM, Dow ER, Chen G, Manji HK. Generation and behavioral characterization of beta-catenin forebrain-specific conditional knock-out mice. Behav Brain Res. 2008;189:117–125. [PMC free article] [PubMed]
  • Mines MA, Yuskaitis CJ, King MK, Beurel E, Jope RS. GSK3 influences social preference and anxiety-related behaviors during social interaction in a mouse model of fragile X syndrome and autism. PLoS One. 2010;5:e9706. [PMC free article] [PubMed]
  • Kivimae S, Martin PM, Kapfhamer D, Ruan Y, Heberlein U, Rubenstein JL, et al. Abnormal behavior in mice mutant for the Disc1 binding partner, Dixdc1. Transl Psychiatry. 2011;1:e43. [PMC free article] [PubMed]
  • Sowers LP, Loo L, Wu Y, Campbell E, Ulrich JD, Wu S, et al. Disruption of the non-canonical Wnt gene PRICKLE2 leads to autism-like behaviors with evidence for hippocampal synaptic dysfunction Mol Psychiatryadvance online publication, 28 May 2013; doi: : 10.1038/mp.2013.71 (e-pub ahead of print). [PubMed]
  • Cui DH, Jiang KD, Jiang SD, Xu YF, Yao H. The tumor suppressor adenomatous polyposis coli gene is associated with susceptibility to schizophrenia. Mol Psychiatry. 2005;10:669–677. [PubMed]
  • Mao Y, Ge X, Frank CL, Madison JM, Koehler AN, Doud MK, et al. Disrupted in schizophrenia 1 regulates neuronal progenitor proliferation via modulation of GSK3beta/beta-catenin signaling. Cell. 2009;136:1017–1031. [PMC free article] [PubMed]
  • Katsu T, Ujike H, Nakano T, Tanaka Y, Nomura A, Nakata K, et al. The human frizzled-3 (FZD3) gene on chromosome 8p21, a receptor gene for Wnt ligands, is associated with the susceptibility to schizophrenia. Neurosci Lett. 2003;353:53–56. [PubMed]
  • Lin PI, Chien YL, Wu YY, Chen CH, Gau SS, Huang YS, et al. The WNT2 gene polymorphism associated with speech delay inherent to autism. Res Dev Disabil. Res Dev Disabil. 2012;33:1533–1540. [PubMed]
  • Wassink TH, Piven J, Vieland VJ, Huang J, Swiderski RE, Pietila J, et al. Evidence supporting WNT2 as an autism susceptibility gene. Am J Med Genet. 2001;105:406–413. [PubMed]
  • Marui T, Funatogawa I, Koishi S, Yamamoto K, Matsumoto H, Hashimoto O, et al. Association between autism and variants in the wingless-type MMTV integration site family member 2 ( WNT2) gene. Int J Neuropsychopharmacol. 2010;13:443–449. [PubMed]
  • Aleksic B, Kushima I, Ito Y, Nakamura Y, Ujike H, Suzuki M, et al. Genetic association study of KREMEN1 and DKK1 and schizophrenia in a Japanese population. Schizophr Res. 2010;118:113–117. [PubMed]
  • Fanous AH, Chen X, Wang X, Amdur RL, O'Neill FA, Walsh D, et al. Association between the 5q31.1 gene neurogenin1 and schizophrenia. Am J Med Genet B Neuropsychiatr Genet. 2007;144B:207–214. [PubMed]
  • Proitsi P, Li T, Hamilton G, Di Forti M, Collier D, Killick R, et al. Positional pathway screen of wnt signaling genes in schizophrenia: association with DKK4. Biol Psychiatry. 2008;63:13–16. [PubMed]
  • O'Roak BJ, Vives L, Girirajan S, Karakoc E, Krumm N, Coe BP, et al. Sporadic autism exomes reveal a highly interconnected protein network of de novo mutations. Nature. 2012;485:246–250. [PMC free article] [PubMed]
  • Miyaoka T, Seno H, Ishino H. Increased expression of Wnt-1 in schizophrenic brains. Schizophr Res. 1999;38:1–6. [PubMed]
  • Aston C, Jiang L, Sokolov BP. Microarray analysis of postmortem temporal cortex from patients with schizophrenia. J Neurosci Res. 2004;77:858–866. [PubMed]
  • Fisher DA, Kivimae S, Hoshino J, Suriben R, Martin PM, Baxter N, et al. Three Dact gene family members are expressed during embryonic development and in the adult brains of mice. Dev Dyn. 2006;235:2620–2630. [PubMed]
  • Veeman MT, Slusarski DC, Kaykas A, Louie SH, Moon RT. Zebrafish prickle, a modulator of noncanonical Wnt/Fz signaling, regulates gastrulation movements. Curr Biol. 2003;13:680–685. [PubMed]
  • Lord C, Rutter M, Le Couteur A. Autism Diagnostic Interview-Revised: a revised version of a diagnostic interview for caregivers of individuals with possible pervasive developmental disorders. J Autism Dev Disord. 1994;24:659–685. [PubMed]
  • Sanders AR, Levinson DF, Duan J, Dennis JM, Li R, Kendler KS, et al. The Internet-based MGS2 control sample: self report of mental illness. Am J Psychiatry. 2010;167:854–865. [PubMed]
  • Calvo SE, Tucker EJ, Compton AG, Kirby DM, Crawford G, Burtt NP, et al. High-throughput, pooled sequencing identifies mutations in NUBPL and FOXRED1 in human complex I deficiency. Nat Genet. 2010;42:851–858. [PMC free article] [PubMed]
  • Nejentsev S, Walker N, Riches D, Egholm M, Todd JA. Rare variants of IFIH1, a gene implicated in antiviral responses, protect against type 1 diabetes. Science. 2009;324:387–389. [PMC free article] [PubMed]
  • Langmead B, Trapnell C, Pop M, Salzberg SL. Ultrafast and memory-efficient alignment of short DNA sequences to the human genome. Genome Biol. 2009;10:R25. [PMC free article] [PubMed]
  • Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics. 2009;25:2078–2079. [PMC free article] [PubMed]
  • Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26:841–842. [PMC free article] [PubMed]
  • DePristo MA, Banks E, Poplin R, Garimella KV, Maguire JR, Hartl C, et al. A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat Genet. 2011;43:491–498. [PMC free article] [PubMed]
  • Rivas MA, Beaudoin M, Gardet A, Stevens C, Sharma Y, Zhang CK, et al. Deep resequencing of GWAS loci identifies independent rare variants associated with inflammatory bowel disease. Nat Genet. 2011;43:1066–1073. [PMC free article] [PubMed]
  • Adzhubei IA, Schmidt S, Peshkin L, Ramensky VE, Gerasimova A, Bork P, et al. A method and server for predicting damaging missense mutations. Nat Methods. 2010;7:248–249. [PMC free article] [PubMed]
  • Biechele TL, Moon RT. Assaying beta-catenin/TCF transcription with beta-catenin/TCF transcription-based reporter constructs. Methods Mol Biol. 2008;468:99–110. [PubMed]
  • Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MA, Bender D, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet. 2007;81:559–575. [PubMed]
  • Grantham R. Amino acid difference formula to help explain protein evolution. Science. 1974;185:862–864. [PubMed]
  • Jho EH, Zhang T, Domon C, Joo CK, Freund JN, Costantini F. Wnt/beta-catenin/Tcf signaling induces the transcription of Axin2, a negative regulator of the signaling pathway. Mol Cell Biol. 2002;22:1172–1183. [PMC free article] [PubMed]
  • Moon RT, Brown JD, Torres M. WNTs modulate cell fate and behavior during vertebrate development. Trends Genet. 1997;13:157–162. [PubMed]
  • Willert K, Nusse R. Beta-catenin: a key mediator of Wnt signaling. Curr Opin Genet Dev. 1998;8:95–102. [PubMed]
  • Fekany-Lee K, Gonzalez E, Miller-Bertoglio V, Solnica-Krezel L. The homeobox gene bozozok promotes anterior neuroectoderm formation in zebrafish through negative regulation of BMP2/4 and Wnt pathways. Development. 2000;127:2333–2345. [PubMed]
  • Takai A, Inomata H, Arakawa A, Yakura R, Matsuo-Takasaki M, Sasai Y. Anterior neural development requires Del1, a matrix-associated protein that attenuates canonical Wnt signaling via the Ror2 pathway. Development. 2010;137:3293–3302. [PubMed]
  • Tang K, Yang J, Gao X, Wang C, Liu L, Kitani H, et al. Wnt-1 promotes neuronal differentiation and inhibits gliogenesis in P19 cells. Biochem Biophys Res Commun. 2002;293:167–173. [PubMed]
  • Ligon KL, Echelard Y, Assimacopoulos S, Danielian PS, Kaing S, Grove EA, et al. Loss of Emx2 funtion leads to ectopic expression of Wnt1 in the developing telencephalon and cortical dysplasia. Development. 2003;130:2275–2287. [PubMed]
  • Bhat KM. Wingless activity in the precursor cells specifies neuronal migratory behavior in the Drosophila nerve cord. Dev Biol. 2007;311:613–622. [PMC free article] [PubMed]
  • Joksimovic M, Patel M, Taketo MM, Johnson R, Awatramani R. Ectopic Wnt/beta-catenin signaling induces neurogenesis in the spinal cord and hindbrain floor plate. PLoS One. 2012;7:e30266. [PMC free article] [PubMed]
  • Yu X, Malenka RC. Beta-catenin is critical for dendritic morphogenesis. Nat Neurosci. 2003;6:1169–1177. [PubMed]
  • Ataman B, Ashley J, Gorczyca M, Ramachandran P, Fouquet W, Sigrist SJ, et al. Rapid activity-dependent modifications in synaptic structure and function require bidirectional Wnt signaling. Neuron. 2008;57:705–718. [PMC free article] [PubMed]
  • De Ferrari GV, Moon RT. The ups and downs of Wnt signaling in prevalent neurological disorders. Oncogene. 2006;25:7545–7553. [PubMed]
  • Buxbaum JD, Daly MJ, Devlin B, Lehner T, Roeder K, State MW. The autism sequencing consortium: large-scale, high-throughput sequencing in autism spectrum disorders. Neuron. 2012;76:1052–1056. [PMC free article] [PubMed]
  • Devlin B, Scherer SW. Genetic architecture in autism spectrum disorder. Curr Opin Genet Dev. 2012;22:229–237. [PubMed]
  • Levy D, Ronemus M, Yamrom B, Lee YH, Leotta A, Kendall J, et al. Rare de novo and transmitted copy-number variation in autistic spectrum disorders. Neuron. 2011;70:886–897. [PubMed]
  • Chow ML, Pramparo T, Winn ME, Barnes CC, Li HR, Weiss L, et al. Age-dependent brain gene expression and copy number anomalies in autism suggest distinct pathological processes at young versus mature ages. PLoS Genet. 2012;8:e1002592. [PMC free article] [PubMed]
  • McCoy PA, Shao Y, Wolpert CM, Donnelly SL, Ashley-Koch A, Abel HL, et al. No association between the WNT2 gene and autistic disorder. Am J Med Genet. 2002;114:106–109. [PubMed]
  • Li J, Nguyen L, Gleason C, Lotspeich L, Spiker D, Risch N, et al. Lack of evidence for an association between WNT2 and RELN polymorphisms and autism. Am J Med Genet B Neuropsychiatr Genet. 2004;126B:51–57. [PubMed]
  • Wang Z, Xu L, Zhu X, Cui W, Sun Y, Nishijo H, et al. Demethylation of specific Wnt/beta-catenin pathway genes and its upregulation in rat brain induced by prenatal valproate exposure. Anat Rec (Hoboken) 2010;293:1947–1953. [PubMed]
  • O'Roak BJ, Vives L, Fu W, Egertson JD, Stanaway IB, Phelps IG, et al. Multiplex targeted sequencing identifies recurrently mutated genes in autism spectrum disorders. Science. 2012;338:1619–1622. [PMC free article] [PubMed]
  • Neale BM, Kou Y, Liu L, Ma'ayan A, Samocha KE, Sabo A, et al. Patterns and rates of exonic de novo mutations in autism spectrum disorders. Nature. 2012;485:242–245. [PMC free article] [PubMed]
  • State MW, Levitt P. The conundrums of understanding genetic risks for autism spectrum disorders. Nat Neurosci. 2011;14:1499–1506. [PMC free article] [PubMed]
  • Talkowski ME, Rosenfeld JA, Blumenthal I, Pillalamarri V, Chiang C, Heilbut A, et al. Sequencing chromosomal abnormalities reveals neurodevelopmental loci that confer risk across diagnostic boundaries. Cell. 2012;149:525–537. [PMC free article] [PubMed]
  • Smoller JW. Identification of risk loci with shared effects on five major psychiatric disorders: a genome-wide analysis. Lancet. 2013;381:1371–1379. [PMC free article] [PubMed]
  • Okerlund ND, Cheyette BN. Synaptic Wnt signaling-a contributor to major psychiatric disorders. J Neurodev Disord. 2011;3:162–174. [PMC free article] [PubMed]
Articles from Translational Psychiatry are provided here courtesy of
Nature Publishing Group