PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of plosgenPLoS GeneticsSubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)View this Article
 
PLoS Genet. 2017 February; 13(2): e1006618.
Published online 2017 February 10. doi:  10.1371/journal.pgen.1006618
PMCID: PMC5328401

Widespread signatures of positive selection in common risk alleles associated to autism spectrum disorder

Scott M. Williams, Editor

Abstract

The human brain is the outcome of innumerable evolutionary processes; the systems genetics of psychiatric disorders could bear their signatures. On this basis, we analyzed five psychiatric disorders, attention deficit hyperactivity disorder, autism spectrum disorder (ASD), bipolar disorder, major depressive disorder, and schizophrenia (SCZ), using GWAS summary statistics from the Psychiatric Genomics Consortium. Machine learning-derived scores were used to investigate two natural-selection scenarios: complete selection (loci where a selected allele reached fixation) and incomplete selection (loci where a selected allele has not yet reached fixation). ASD GWAS results positively correlated with incomplete-selection (p = 3.53*10−4). Variants with ASD GWAS p<0.1 were shown to have a 19%-increased probability to be in the top-5% for incomplete-selection score (OR = 1.19, 95%CI = 1.11–1.8, p = 9.56*10−7). Investigating the effect directions of minor alleles, we observed an enrichment for positive associations in SNPs with ASD GWAS p<0.1 and top-5% incomplete-selection score (permutation p<10−4). Considering the set of these ASD-positive-associated variants, we observed gene-expression enrichments for brain and pituitary tissues (p = 2.3*10−5 and p = 3*10−5, respectively) and 53 gene ontology (GO) enrichments, such as nervous system development (GO:0007399, p = 7.57*10−12), synapse organization (GO:0050808, p = 8.29*10−7), and axon guidance (GO:0007411, p = 1.81*10−7). Previous genetic studies demonstrated that ASD positively correlates with childhood intelligence, college completion, and years of schooling. Accordingly, we hypothesize that certain ASD risk alleles were under positive selection during human evolution due to their involvement in neurogenesis and cognitive ability.

Author summary

Predisposition to psychiatric disorders is due to the contribution of many genes involved in numerous molecular mechanisms. Since brain evolution has played a pivotal role in determining the success of the human species, the molecular pathways involved with the onset of mental illnesses are likely to be informative as we seek an understanding of the mechanisms involved in the evolution of human brain. Accordingly, we tested whether the genetics of psychiatric disorders is enriched for signatures of positive selection. We observed a strong finding related to the genetics of autism spectrum disorders (ASD): common risk alleles are enriched for genomic signatures of incomplete selection (loci where a selected allele has not yet reached fixation). The genes where these alleles map tend to be expressed in brain and pituitary tissues, to be involved in molecular mechanisms related to nervous system development, and surprisingly, to be associated with increased cognitive ability. Previous studies identified signatures of purifying selection in genes affected by ASD rare alleles. Accordingly, at least two different evolutionary mechanisms appear to be present in relation to ASD genetics: 1) rare disruptive alleles eliminated by purifying selection; 2) common alleles selected for their beneficial effects on cognitive skills. This scenario would explain ASD prevalence, which is higher than that expected for a trait under purifying selection, as the evolutionary cost of polygenic adaptation related to cognitive ability.

Introduction

The human brain is a uniquely complex organ; it is the outcome of numerous evolutionary processes that were necessary for the success of the human species [1]. The same mechanisms that contributed to the evolution of human brain are likely to be involved in the pathogenesis of mental illnesses [2, 3]. Numerous evolutionary hypotheses have been proposed to account for the observation that phenotypic traits with deleterious effects on fitness such as psychotic and mood disorders, have not been removed by natural selection [4]. Risk alleles with large effects on predisposition to mental illness should be under negative selection, at least insofar as they interfere with reproductive fitness. This does appear to be true for rare and de novo mutations strongly associated with psychiatric disorders [5], but most genetic risk for these disorders is attributable to a polygenic predisposition [6]. Specifically, heritability analyses have shown that major psychiatric diseases are highly polygenic: their genetic predisposition should be due the additive result of hundred to thousand variants with small effect size [79].

Similarly, genome evolution also seems to operate mainly on gene networks rather than single genes [10]. Different authors have hypothesized that human adaptation in response to the selection of polygenic phenotypes may occur via subtle allele frequency shifts at many loci [1113]. Signatures of polygenic adaptation have been identified in the context of several phenotypic traits, including immune response [14, 15], anthropometric traits [1619], metabolic traits [15, 17], skin pigmentation [17], telomere length [20], bone mineral density [21], and dietary patterns [22]. Two previous studies investigated the genome-wide enrichment of evolutionary signatures among schizophrenia (SCZ) risk alleles using the results of the large genome-wide association study (GWAS) meta-analysis conducted by the Psychiatric Genomics Consortium (PGC) [23]. In 2015, Xu and colleagues reported that genes near human accelerated regions conserved in non-human primates (pHARs) are enriched for SCZ-associated loci, and they are particularly related to the GABA (gamma-Aminobutyric acid)-related co-expression module [24]. In 2016, Srinivasan and colleagues obtained consistent results using a different approach based on Neanderthal selective sweep (NSS) score, which is an indicator of positive selection in early humans based on the depletion of Neanderthal-derived alleles [25]. They observed loci associated with schizophrenia are more prevalent in regions with a low NSS score (i.e., this is evidence of positive selection) [25].

In the present study, we investigated polygenic adaptation signatures in the systems genetics of five psychiatric disorders: attention deficit hyperactivity disorder (ADHD), autism spectrum disorder (ASD), bipolar disorder (BP), major depressive disorder (MDD), and SCZ.

Results

Our investigation was conducted using summary statistics from GWAS of psychiatric disorders and testing for enrichment of positive selection signatures, using the same concept as is applied in high-resolution polygenic risk score analysis [26, 27]. The positive selection signatures were identified using the hierarchical boosting (HB) algorithm, which is a machine-learning classification framework that combines the functionality of several selection tests to uncover different genetic features that are expected under selective sweeps [28].

Initially, we verified whether the GWAS significance (-log10 p value) reported for the variants investigated correlate with the HB scores related to the corresponding genomic regions for incomplete and complete selection (loci where a selected allele has not yet reached fixation and loci where a selected allele has reached fixation, respectively) using a non-parametric test (Spearman's correlation). Significant positive correlations were observed between ASD GWAS results and HB scores for incomplete selection (p = 3.53*10−4): higher GWAS significance correlated with higher HB scores (S1 Table). We confirmed these results by conducting 10,000 random permutations of the GWAS results with respect to the corresponding HB scores and testing whether the observed correlations were significantly higher from the ones in the null distribution of the permuted datasets (ASD vs. incomplete selection–permutation p = 1*10−4, S1 Fig). Then, we investigated which GWAS significance threshold is more enriched for natural selection signatures. We considered the top-5% of HB scores as suggestive evidence of natural selection. Variants with ASD GWAS p < 0.1 have a 19%-increased probability to be in the top-5% of the HB scores for incomplete selection (Odds Ratio (OR) = 1.19, 95%CI = 1.11–1.8, p = 9.56*10−7; Fig 1; S2 Table). To evaluate the direction of these natural-selection enrichments, we tested whether the variants included within the significant thresholds (ASD GWAS p < 0.1) and with suggestive evidence of positive selection (top-5% HB scores for incomplete and complete selection, respectively) show an overrepresentation with respect to a specific effect direction (positive association [GWAS OR > 1] vs. negative association [GWAS OR < 1). Taking into account the effects of minor alleles, we compared the median OR of the variants identified with respect to the median value of ORs calculated on the basis of 10,000 permutations of the original ASD. The median value of OR for variants with ASD GWAS p < 0.1 and HB score for incomplete selection in the top-5% was 1.057; this resulted in higher than the median values of ORs from permuted datasets (p < 10−4; Fig 2). Considering the variants with ASD GWAS p < 0.1, top-5% incomplete-selection score, and ASD GWAS OR > 1, we conducted enrichment analysis for tissue-specific gene expression and for Gene Ontologies (GO) related to biological processes to gain insights regarding the molecular mechanisms involved. SNPs were assigned to the genes where they are located and/or to the nearest genes (±50KB). With a false discovery rate (FDR) < 5%, we observed significant enrichments for genes highly expressed in brain and pituitary tissues (p = 2.3*10−5 and p = 3*10−5, respectively; S3 Table). Applying a type I error rate at 5% after Bonferroni multiple testing correction, 53 GO enrichments (Table 1) were identified with the top GO result related to nervous system development (GO:0007399, p = 7.57*10−12). Gene sets related to the parent GO terms are reported in S4 Table. Considering the significant GO enrichments, we observed a large similarity-based network including several terms related to nervous system development (Fig 3).

Fig 1
Enrichments for scores related to incomplete selection in ASD GWAS considering different thresholds.
Fig 2
Null distribution of median OR (considering minor allele effect) generated from 10,000 permutations of ASD dataset.
Fig 3
Similarity network based on GO enrichment results.
Table 1
Gene Ontology (GO) enrichments for variants with ASD GWAS OR >1 and p < 0.1 and located in genomic regions with HB score for incomplete selection in top-5%.

To investigate our ASD results further, we considered ASD genetic correlation results available from LD Hub v1.3.1 [29] (available at http://ldsc.broadinstitute.org/). Applying a FDR at 5%, we identified eight significant correlations (S5 Table). Among these, we observed positive correlation with several advantageous traits: years of schooling (rg = 0.277, p = 2.9*10−13), college completion (rg = 0.339, p = 1*10−6), childhood intelligence (rg = 0.425, 5.74*10−5), openness to experience (rg = 0.421, p = 2*10−3).

We also observed nominally significant results in relation to SCZ: a positive correlation between SCZ GWAS results and HB scores for complete selection (p = 6.37*10−3, S1 Table; permutation p = 3*10−3, S2 Fig). Although this result did not survive after multiple testing correction (see Statistical Analysis), we nevertheless tested enrichment for GWAS significance thresholds because this nominal result replicated previous findings obtained by independent studies with different methods [24, 25]. We observed a 22%-increased probability to be in the top-5% of scores for incomplete selection for variants with SCZ GWAS p<0.001 (OR = 1.22, 95%CI = 1.05–1.41, p = 7.9*10−3; Fig 4). Because the result of this subsequent analysis also did not survive Bonferroni multiple testing correction, we did not investigate the SCZ results further. No significant result was observed for the correlation analysis conducted using ADHD, BD, and MDD GWAS summary statistics (S1 Table).

Fig 4
Genome-wide enrichments for scores related to incomplete selection considering different SCZ significance thresholds.

Discussion

To our knowledge, the current findings provide the first genome-wide evidence for the strong presence of natural-selection signatures in the systems genetics of ASD. Although many authors hypothesized a strong involvement of evolutionary mechanisms in ASD [4] and several studies investigated candidate genes [30]. no evidence of widespread positive selection on ASD-associated SNPs was previously reported. Using genome-wide data, we observed that common alleles associated with increased risk for ASD present a signature of positive selection in European populations. This strongly suggests that these variants have undergone positive selection during the course of human evolutionary history. Genetic correlation results support this hypothesis: ASD genetics (i.e., the set of risk variants that collectively–on a population level–influences ASD risk) strongly correlates with years of schooling, college completion, childhood intelligence, and openness to experience (S5 Table, data available at http://ldsc.broadinstitute.org/) [29]. Although these are robust genetic correlation results, further studies are needed to confirm the role of cognitive abilities in the evolutionary mechanisms involved in ASD genetics. However, multiple forms of evidence support that autism and high intelligence quotient share a diverse set of correlates, such as large brain size, fast brain growth, increased sensory and visual-spatial abilities, enhanced synaptic functions, increased attentional focus, high socioeconomic status, more deliberative decision-making, and high levels of positive assortative mating [31]. Accordingly, the genomic signatures observed in ASD risk alleles could be due to their positive associations with cognitive ability. In 2012, the Center for Disease Control and Prevention reported an ASD prevalence of 1.47% in the US population (1 in 68 children; data available at http://www.cdc.gov/ncbddd/autism/data.html). These epidemiological data can be interpreted in relation to the recent studies regarding the genetics of ASD where common genetic variation explains 49% of the ASD heritability and inherited rare and de novo mutations account for only 6% [9, 32]. Considering rare variation, the strongest contribution seems to be played by de novo single nucleotide variants (SNV) and copy number variants (CNV) with a reduced (but still significant) contribution from rare inherited variants [33]. In particular, rare variation associated with ASD is enrichened for disruptive truncating alleles [3336]. Accordingly, two different evolutionary mechanisms are likely to be present in relation to ASD genetics. Rare disruptive variants predisposing to ASD are under strong purifying selection (selective removal of deleterious alleles), as already observed in certain ASD genes [37, 38]. Conversely, genetic predisposition to ASD due to common variants is highly polygenic and, taken together on a population level, these alleles present beneficial effects with respect to cognitive ability. This positive selection for ASD risk alleles increased their occurrence in human populations, and this provides a possible explanation for the disease prevalence observed by epidemiological studies. A trait related to rare alleles under strong purifying selection should present much lower prevalence. Furthermore, we observed that common ASD risk alleles with evidence of positive selection are enriched for many biological processes related to developmental mechanisms and, in particular, to mechanisms related to nervous system development. This agrees with the strong evidence indicating that the processes related to human brain development are more responsible for distinctive human traits [39]. Accordingly, ASD risk alleles could positively affect these mechanisms, causing better cognitive ability in carriers as a consequence. However, an excessive burden of these risk variants is correlated with the onset of the developmental disorders included in the autism spectrum as the evolutionary cost.

Although they did not survive multiple testing correction, the results of our SCZ analysis are consistent with two previous studies that observed enrichments for selection signatures based on inter-species comparisons (pHARs: conserved regions in non-human primates vs. humans; NSS: selective sweep based Neanderthals vs. Humans comparisons) [24, 25]. Our results suggest that natural-selection processes related to SCZ genes were also present within human evolutionary history and not only in the divergence between Homo Sapiens and non-human species. However, these intra-species adaptation mechanisms were particularly strong because they caused fixation of the selected allele (complete selection). We can also speculate that different evolutionary mechanisms may have acted on the genetic architecture of these two traits. Indeed, it appears that SCZ risk alleles associated with positive selection can reach fixation (complete selection) more easily than ASD risk alleles (incomplete selection).

Although SCZ GWAS has so far been more powerful than ASD (over 100 genome-wide significant SCZ loci vs. 1 genome-wide significant ASD locus identified; this is largely due to the sample size difference; S3 Fig), the strongest evolutionary findings were observed in ASD analysis, suggesting that natural selection acted powerfully on the systems genetics of this disorder.

As mentioned above, signatures of polygenic adaptation have been identified in the context of several phenotypic traits. However, for some of them (e.g., immune response, skin pigmentation, and anthropometric traits) single-locus signals of selection are also present. Conversely, to our knowledge no single-locus evolutionary signature has been directly related to brain functions. This can support a speculative hypothesis about the genetic architecture of brain functions. Genetic mechanisms pivotal for neural activities are strongly conserved and alleles with large effect on gene functions have an extremely high probability of being deleterious (and of consequent elimination by purifying selection). Conversely, alleles with small effects could modify the brain systems more subtly and, in some cases, provide small beneficial effects. According to our interpretation of our data, such small-effect alleles were accumulated across the genome (polygenic adaptation) to the benefit of most but to the detriment of some.

In conclusion, the present study provides evidence regarding the role of human evolution in shaping the genetic architecture of psychiatry disorders, providing a hypothesis to explain the ASD prevalence as the evolutionary cost of the polygenic adaptation of the disease risk alleles.

Materials and methods

Ethics statement

Publicly available GWAS summary statistics from the PGC and computational methods were used and therefore no additional ethics approval was needed. The ethics approval of PGC studies can be found in the related articles [6, 23, 4042].

GWAS summary statistics

We used the summary statistics of the large GWAS meta-analyses conducted on these traits by the PGC (data available at https://www.med.unc.edu/pgc/results-and-downloads) [6, 23, 4042]. Details regarding the summary statistics used are reported in S6 Table. To exclude bias related to linkage disequilibrium (LD), we clumped the data considering the European reference panel from the 1,000 Genomes Project [43] and the following parameters: 500kb window, r2 < 0.25, imputation info score > 0.9, and minor allele frequency (MAF) > 0.02. For SCZ we included a single MHC single nucleotide polymorphism (SNP). To our knowledge, the PGC ASD GWAS was conducted on the largest cohort currently published and it includes samples from the Geschwind Autism Center of Excellence, the Autism Genome Project [44, 45], the Autism Genetic Resource Exchange[46, 47], the Montreal/Boston Collection [48], and the Simons Simplex Collection [49]. We were unable to identify any additional available ASD GWAS for replication.

Hierarchical boosting algorithm

We considered natural-selection scores derived from the HB algorithm, which is a machine-learning classification framework that combines the functionality of several selection tests to uncover different genetic features expected under selective sweeps [28].

Specifically, HB is based on a machine-learning algorithm called boosting (from the mboost R package), which is a supervised algorithm that estimates linear regressions of input variables (summary statistics of selection tests) to maximize the differences between two competing scenarios (complete vs. incomplete selective sweeps). The HB method sequentially applies different boosting functions to a hierarchical classification scheme to classify optimally genomic regions into different evolutive regimes. The HB method was successfully tested with respect to simulations in a period between 10 thousand years ago (Kya) and 45 Kya, which should include both ancient and recent selective sweeps [28]. In our study, we considered two different natural-selection scenarios: complete selection (loci where a selected allele reached fixation) and incomplete selection (loci where a selected allele has not yet reached fixation). The HB method assigns a score with respect a specific genomic region. Accordingly, we assigned the HB scores to the alleles investigated in GWAS based on their locations. Since the original PGC GWAS were conducted on European populations, we considered the HB scores calculated on a European (CEU) population from the 1,000 Genomes Project.

Statistical analysis

Statistical analysis was performed using the computing environment R (https://www.r-project.org/). For the initial correlation analysis (GWAS significance, -log10 p value, vs. HB score), we applied a Bonferroni correction accounting for the number of psychiatric disorders (ADHD, ASD, BP, MDD, SCZ) and the number of selection signals (incomplete-selection and complete-selection) to calculate the significance threshold adjusted for multiple testing correction (p = 5*10−3). For the subsequent testing of multiple GWAS significance thresholds, we applied a Bonferroni correction for the number of GWAS significance thresholds tested (ASD: 6, p = 8.3*10−3; SCZ: 9, p = 5.5*10−3). The difference in the number of GWAS significance thresholds tested is due to the difference of statistical power in the original GWAS (S3 Fig). The power analysis applied to the original ASD and SCZ GWAS was conducted using QUANTO software (http://biostats.usc.edu/Quanto.html).

To gain insights about the molecular mechanisms involved in the evolutionary pressures related to the genetics of psychiatry disorders, the results were further investigated to understand the direction of these evolutionary enrichments and to identify which tissues and biological processes are involved. SNPs were assigned to the genes where they are located and/or to the nearest genes (±50KB). DAVID v6.8 [50] (https://david.ncifcrf.gov/) was used to verify the enrichment for tissue-specific gene expression. DAVID v6.8. uses the following gene expression databases (DAVID Update Date: May 2016): Cancer Genome Anatomy Project (CGAP_EST_QUARTILE), Cancer Genome Anatomy Project Serial Analysis of Gene Expression (CGAP_SAGE_QUARTILE), U133A database of the Genomics Institute of the Novartis Research Foundation (GNF_U133A_QUARTILE), UniGene (UNIGENE_EST_QUARTILE), and Uniprot tissue (UP_TISSUE). The FDR was applied to correct the results for multiple testing [51] and q values < 0.05 were considered to be significant. Panther v11.0 [52] (http://www.pantherdb.org/) was used to conduct an overrepresentation test considering the biological processes included in the GO database for Homo Sapiens (released 2016-07-29). This statistical test calculates whether, considering a specific functional category, the genes included in the analysis are statistically different (overrepresented or underrepresented) from the chance expectation. Bonferroni correction for the number of tests conducted (i.e., GO for biological processes) was applied to adjust the results. To further validate our GO-enrichment results, we randomly permuted the original dataset 100 times and performed the GO-enrichment analysis on the permuted data. The number of GO enrichments observed was significantly higher than null distribution generated by the random permutations (ppermutation < 0.01; S4 Fig). The GO enrichment results were further investigated using REVIGO [53] (available at http://revigo.irb.hr/). Specifically, GO enrichments were used to make a graph-based visualization considering an allowed similarity of 0.7, UniProt as reference database, and Jian and Conrath method as the semantic similarity measure.

To further follow-up our findings, we used the information about genetic correlations provided by LD Hub v1.3.1 [29] (available at http://ldsc.broadinstitute.org/ldhub/). This web-tool provides information regarding genetic correlation results for 189 traits. The genetic correlations were calculated using the LD score method [54]. In our study, we applied the FDR [51] to correct the results for multiple testing and q values < 0.05 were considered to be significant.

Supporting information

S1 Table

Correlations (Spearman's rho, p value) between GWAS significance of psychiatric disorders and HB scores for incomplete and complete selection.

(DOCX)

S2 Table

Enrichment for incomplete-selection signals in ASD GWAS considering different thresholds.

(DOCX)

S3 Table

Details of the gene-expression enrichment results.

(DOCX)

S4 Table

Gene sets related to the parent GO terms.

(DOCX)

S5 Table

Genetic correlations between ASD and the phenotype traits available in the LD hub database v 1.3.1 (available at http://ldsc.broadinstitute.org/).

(DOCX)

S6 Table

Details of the GWAS summary statistics used in the present study.

(DOCX)

S1 Fig

Distribution of Spearman's rho generated from 10,000 permutations of ASD dataset.

Blue lines represent the Spearman's rho observed in the real data.

(DOCX)

S2 Fig

Distribution of Spearman's rho generated from 10,000 permutations of SCZ datasets.

Blue lines represent the Spearman's rho observed in the real data.

(DOCX)

S3 Fig

Statistical power of PGC ASD and SCZ GWAS calculated considering different thresholds, a minor allele frequency of 10%, and an effect size of OR = 1.1.

The sample size for ASD and SCZ cohorts are those reported in S6 Table.

(DOCX)

S4 Fig

GO-enrichment results in 100 random permutations.

We report the null distribution and the frequency the GOs observed. In the ASD data, we observed enrichment for 53 GOs (ppermutation<0.01).

(DOCX)

Acknowledgments

The authors are grateful to the research groups of the Psychiatric Genomics Consortium for their publicly available GWAS summary statistics.

Funding Statement

This study was supported by National Institutes of Health grants R01 DA12690, R01 AA017535, P50 AA012870, U01 MH109532, R21 AA024404, the VA Connecticut MIRECC, and a NARSAD Young Investigator Award (to RP) from the Brain & Behavior Research Foundation. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Data Availability

Data Availability

The summary statistics of the GWAS of psychiatric diseases are available at https://www.med.unc.edu/pgc/results-and-downloads.

References

1. Teffer K, Semendeferi K. Human prefrontal cortex: evolution, development, and pathology. Prog Brain Res. 2012;195:191–218. doi: 10.1016/B978-0-444-53860-4.00009-X [PubMed]
2. Agnati LF, et al. Possible genetic and epigenetic links between human inner speech, schizophrenia and altruism. Brain Res. 2012;1476:38–57. doi: 10.1016/j.brainres.2012.02.074 [PubMed]
3. Nettle D, Bateson M. The evolutionary origins of mood and its disorders. Curr Biol. 2012;22:R712–21. doi: 10.1016/j.cub.2012.06.020 [PubMed]
4. Daly M, et al. Natural Selection and Neuropsychiatric Disease: Theory, Observation, and Emerging Genetic Findings In: Lehner T, Miller BL, State MW, editors. Genomics, Circuits, and Pathways in Clinical Neuropsychiatry: Academic Press; 2016. p. 51–61.
5. Rees E, et al. De novo mutation in schizophrenia. Schizophr Bull. 2012;38:377–81. doi: 10.1093/schbul/sbs047 [PMC free article] [PubMed]
6. Robinson EB, et al. Genetic risk for autism spectrum disorders and neuropsychiatric variation in the general population. Nat Genet. 2016;48:552–5. doi: 10.1038/ng.3529 [PMC free article] [PubMed]
7. Loh PR, et al. Contrasting genetic architectures of schizophrenia and other complex diseases using fast variance-components analysis. Nat Genet. 2015;47:1385–92. doi: 10.1038/ng.3431 [PMC free article] [PubMed]
8. Kendler KS. A joint history of the nature of genetic variation and the nature of schizophrenia. Mol Psychiatry. 2015;20:77–83. doi: 10.1038/mp.2014.94 [PMC free article] [PubMed]
9. Gaugler T, et al. Most genetic risk for autism resides with common variation. Nat Genet. 2014;46:881–5. doi: 10.1038/ng.3039 [PMC free article] [PubMed]
10. Jain K, Stephan W. Response of Polygenic Traits Under Stabilizing Selection and Mutation When Loci Have Unequal Effects. G3 (Bethesda). 2015;5:1065–74. doi: 10.1534/g3.115.017970 . [PMC free article] [PubMed]
11. Pritchard JK, et al. The genetics of human adaptation: hard sweeps, soft sweeps, and polygenic adaptation. Curr Biol. 2010;20:R208–15. doi: 10.1016/j.cub.2009.11.055 [PMC free article] [PubMed]
12. Pritchard JK, Di Rienzo A. Adaptation—not by sweeps alone. Nat Rev Genet. 2010;11:665–7. doi: 10.1038/nrg2880 [PMC free article] [PubMed]
13. Stephan W. Signatures of positive selection: from selective sweeps at individual loci to subtle allele frequency changes in polygenic adaptation. Mol Ecol. 2016;25:79–88. doi: 10.1111/mec.13288 [PubMed]
14. Daub JT, et al. Evidence for polygenic adaptation to pathogens in the human genome. Mol Biol Evol. 2013;30:1544–58. doi: 10.1093/molbev/mst080 [PubMed]
15. Vatsiou AI, et al. Changes in selective pressures associated with human population expansion may explain metabolic and immune related pathways enriched for signatures of positive selection. BMC Genomics. 2016;17:504 doi: 10.1186/s12864-016-2783-2 [PMC free article] [PubMed]
16. Polimanti R, et al. Evidence of Polygenic Adaptation in the Systems Genetics of Anthropometric Traits. PLoS One. 2016;11:e0160654 doi: 10.1371/journal.pone.0160654 [PMC free article] [PubMed]
17. Berg JJ, Coop G. A population genetic signal of polygenic adaptation. PLoS Genet. 2014;10:e1004412 doi: 10.1371/journal.pgen.1004412 [PMC free article] [PubMed]
18. Turchin MC, et al. Evidence of widespread selection on standing variation in Europe at height-associated SNPs. Nat Genet. 2012;44:1015–9. doi: 10.1038/ng.2368 [PMC free article] [PubMed]
19. Robinson MR, et al. Population genetic differentiation of height and body mass index across Europe. Nat Genet. 2015;47:1357–62. doi: 10.1038/ng.3401 [PMC free article] [PubMed]
20. Hansen ME, et al. Shorter telomere length in Europeans than in Africans due to polygenetic adaptation. Hum Mol Genet. 2016. 25:2324–30. doi: 10.1093/hmg/ddw070 . [PMC free article] [PubMed]
21. Medina-Gomez C, et al. BMD Loci Contribute to Ethnic and Developmental Differences in Skeletal Fragility across Populations: Assessment of Evolutionary Selection Pressures. Mol Biol Evol. 2015;32:2961–72. doi: 10.1093/molbev/msv170 [PMC free article] [PubMed]
22. White L, et al. Genetic adaptation to levels of dietary selenium in recent human history. Mol Biol Evol. 2015;32:1507–18. doi: 10.1093/molbev/msv043 [PubMed]
23. Schizophrenia Working Group of the Psychiatric Genomics Consortium. Biological insights from 108 schizophrenia-associated genetic loci. Nature. 2014;511:421–7. doi: 10.1038/nature13595 [PMC free article] [PubMed]
24. Xu K, et al. Genomic and network patterns of schizophrenia genetic variation in human evolutionary accelerated regions. Mol Biol Evol. 2015;32:1148–60. doi: 10.1093/molbev/msv031 [PMC free article] [PubMed]
25. Srinivasan S, et al. Genetic Markers of Human Evolution Are Enriched in Schizophrenia. Biol Psychiatry. 2016;80:284–92. doi: 10.1016/j.biopsych.2015.10.009 [PubMed]
26. Euesden J, et al. PRSice: Polygenic Risk Score software. Bioinformatics. 2015;31:1466–8. doi: 10.1093/bioinformatics/btu848 [PMC free article] [PubMed]
27. Polimanti R, et al. Cross-Phenotype Polygenic Risk Score Analysis of Persistent Post-Concussive Symptoms in U.S. Army Soldiers with Deployment-Acquired Traumatic Brain Injury. J Neurotrauma. 2016. doi: 10.1089/neu.2016.4550 . [PubMed]
28. Pybus M, et al. Hierarchical boosting: a machine-learning framework to detect and classify hard selective sweeps in human populations. Bioinformatics. 2015;31:3946–52. doi: 10.1093/bioinformatics/btv493 [PubMed]
29. ZHENG J, et al. LD Hub: a centralized database and web interface to perform LD score regression that maximizes the potential of summary level GWAS data for SNP heritability and genetic correlation analysis. Bioinformatics. 2017; 33:272–9. doi: 10.1093/bioinformatics/btw613 . [PubMed]
30. Oksenberg N, et al. Function and regulation of AUTS2, a gene implicated in autism and human evolution. PLoS Genet. 2013;9:e1003221 doi: 10.1371/journal.pgen.1003221 [PMC free article] [PubMed]
31. Crespi BJ. Autism As a Disorder of High Intelligence. Front Neurosci. 2016;10:300 doi: 10.3389/fnins.2016.00300 [PMC free article] [PubMed]
32. de la Torre-Ubieta L, et al. Advancing the understanding of autism disease mechanisms through genetics. Nat Med. 2016;22:345–61. doi: 10.1038/nm.4071 [PMC free article] [PubMed]
33. Krumm N, et al. Excess of rare, inherited truncating mutations in autism. Nat Genet. 2015;47:582–8. doi: 10.1038/ng.3303 [PMC free article] [PubMed]
34. De Rubeis S, et al. Synaptic, transcriptional and chromatin genes disrupted in autism. Nature. 2014;515:209–15. doi: 10.1038/nature13772 [PMC free article] [PubMed]
35. Sanders SJ, et al. Insights into Autism Spectrum Disorder Genomic Architecture and Biology from 71 Risk Loci. Neuron. 2015;87:1215–33. doi: 10.1016/j.neuron.2015.09.016 [PMC free article] [PubMed]
36. Turner TN, et al. Genome Sequencing of Autism-Affected Families Reveals Disruption of Putative Noncoding Regulatory DNA. Am J Hum Genet. 2016;98:58–74. doi: 10.1016/j.ajhg.2015.11.023 [PubMed]
37. Iossifov I, et al. De novo gene disruptions in children on the autistic spectrum. Neuron. 2012;74:285–99. doi: 10.1016/j.neuron.2012.04.009 [PMC free article] [PubMed]
38. Tsur E, et al. The Unique Evolutionary Signature of Genes Associated with Autism Spectrum Disorder. Behav Genet. 2016. 46:754–62. doi: 10.1007/s10519-016-9804-4 . [PubMed]
39. Kuzawa CW, et al. Metabolic costs and evolutionary implications of human brain development. Proc Natl Acad Sci U S A. 2014;111:13010–5. doi: 10.1073/pnas.1323099111 [PubMed]
40. Psychiatric GWAS Consortium Bipolar Disorder Working Group. Large-scale genome-wide association analysis of bipolar disorder identifies a new susceptibility locus near ODZ4. Nat Genet. 2011;43:977–83. doi: 10.1038/ng.943 [PMC free article] [PubMed]
41. Neale BM, et al. Meta-analysis of genome-wide association studies of attention-deficit/hyperactivity disorder. J Am Acad Child Adolesc Psychiatry. 2010;49:884–97. doi: 10.1016/j.jaac.2010.06.008 [PMC free article] [PubMed]
42. Major Depressive Disorder Working Group of the Psychiatric GC, et al. A mega-analysis of genome-wide association studies for major depressive disorder. Mol Psychiatry. 2013;18:497–511. doi: 10.1038/mp.2012.21 [PMC free article] [PubMed]
43. 1000 Genomes Project Consortium, et al. A global reference for human genetic variation. Nature. 2015;526:68–74. doi: 10.1038/nature15393 [PMC free article] [PubMed]
44. Anney R, et al. A genome-wide scan for common alleles affecting risk for autism. Hum Mol Genet. 2010;19:4072–82. doi: 10.1093/hmg/ddq307 [PMC free article] [PubMed]
45. Anney R, et al. Individual common variants exert weak effects on the risk for autism spectrum disorders. Hum Mol Genet. 2012;21:4781–92. doi: 10.1093/hmg/dds301 [PMC free article] [PubMed]
46. Lajonchere CM, Consortium A. Changing the landscape of autism research: the autism genetic resource exchange. Neuron. 2010;68:187–91. doi: 10.1016/j.neuron.2010.10.009 [PMC free article] [PubMed]
47. Geschwind DH, et al. The autism genetic resource exchange: a resource for the study of autism and related neuropsychiatric conditions. Am J Hum Genet. 2001;69:463–6. doi: 10.1086/321292 [PubMed]
48. Weiss LA, et al. A genome-wide linkage and association scan reveals novel loci for autism. Nature. 2009;461:802–8. doi: 10.1038/nature08490 [PMC free article] [PubMed]
49. Chaste P, et al. A genome-wide association study of autism using the Simons Simplex Collection: Does reducing phenotypic heterogeneity in autism increase genetic homogeneity? Biol Psychiatry. 2015;77:775–84. doi: 10.1016/j.biopsych.2014.09.017 [PMC free article] [PubMed]
50. Huang da W, et al. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57. doi: 10.1038/nprot.2008.211 [PubMed]
51. Benjamini Y, Hochberg Y. Controlling the False Discovery Rate: A Practical and Powerful Approach to Multiple Testing. Journal of the Royal Statistical Society Series B (Methodological). 1995;57:289–300. doi: 10.2307/2346101
52. Mi H, et al. Large-scale gene function analysis with the PANTHER classification system. Nat Protoc. 2013;8:1551–66. doi: 10.1038/nprot.2013.092 [PubMed]
53. Supek F, et al. REVIGO summarizes and visualizes long lists of gene ontology terms. PLoS One. 2011;6:e21800 doi: 10.1371/journal.pone.0021800 [PMC free article] [PubMed]
54. Bulik-Sullivan B, et al. An atlas of genetic correlations across human diseases and traits. Nat Genet. 2015;47:1236–41. doi: 10.1038/ng.3406 [PMC free article] [PubMed]

Articles from PLoS Genetics are provided here courtesy of Public Library of Science