Topoisomerase 1 facilitates nucleosome reassembly at stress genes during recovery.

Chromatin remodeling is essential to allow full development of alternative gene expression programs in response to environmental changes. In fission yeast, oxidative stress triggers massive transcriptional changes including the activation of hundreds of genes, with the participation of histone modifying complexes and chromatin remodelers. DNA transcription is associated to alterations in DNA topology, and DNA topoisomerases facilitate elongation along gene bodies. Here, we test whether the DNA topoisomerase Top1 participates in the RNA polymerase II-dependent activation of the cellular response to oxidative stress. Cells lacking Top1 are resistant to H2O2 stress. The transcriptome of Δtop1 strain was not greatly affected in the absence of stress, but activation of the anti-stress gene expression program was more sustained than in wild-type cells. Top1 associated to stress open reading frames. While the nucleosomes of stress genes are partially and transiently evicted during stress, the chromatin configuration remains open for longer times in cells lacking Top1, facilitating RNA polymerase II progression. We propose that, by removing DNA tension arising from transcription, Top1 facilitates nucleosome reassembly and works in synergy with the chromatin remodeler Hrp1 as opposing forces to transcription and to Snf22 / Hrp3 opening remodelers.

Description of a novel splice site variant in UBA1 gene causing VEXAS syndrome.

OBJECTIVE: The vacuoles, E1 enzyme, X-linked, autoinflammatory, somatic (VEXAS) syndrome is a complex immune disorder consequence of somatic UBA1 variants. Most reported pathogenic UBA1 variants are missense or splice site mutations directly impairing the translational start site at p. Met41, with recent studies showing that these variants are frequent causes of recurrent inflammation in older individuals. Here we aimed to characterize a novel UBA1 variant found in two patients clinically presenting with VEXAS syndrome. METHODS: Patients' data were collected from direct assessments and from their medical charts. Genomics analyses were performed by both Sanger and amplicon-based deep sequencing, mRNA studies were performed by both cDNA subcloning and mRNA sequencing. RESULTS: We report a novel, somatic variant in a canonical splice site of the UBA1 gene (c.346-2A>G), which was identified in two unrelated adult male patients with late-onset, unexplained inflammatory manifestations including recurrent fever, Sweet syndrome-like neutrophilic dermatosis, and lung inflammation responsive only to glucocorticoids. RNA analysis from patients' samples demonstrated aberrant mRNA splicing leading to multiple in-frame transcripts, including a transcript retaining the full sequence of intron 4 and a different transcript with the deletion of the first 15 nucleotides of exon 5. CONCLUSION: Here we describe the abnormal UBA1 transcription as a consequence of the novel c.346-2A>G variant identified in two patients with clinical features compatible with VEXAS syndrome. Overall, these results further demonstrate the expanding spectrum of variants in UBA1 leading to pathology and support for a complete gene evaluation in those candidate patients for VEXAS syndrome.

Spanish cohort of VEXAS syndrome: clinical manifestations, outcome of treatments and novel evidences about UBA1 mosaicism.

BACKGROUND: The vacuoles, E1-enzyme, X linked, autoinflammatory and somatic (VEXAS) syndrome is an adult-onset autoinflammatory disease (AID) due to postzygotic UBA1 variants. OBJECTIVES: To investigate the presence of VEXAS syndrome among patients with adult-onset undiagnosed AID. Additional studies evaluated the mosaicism distribution and the circulating cytokines. METHODS: Gene analyses were performed by both Sanger and amplicon-based deep sequencing. Patients' data were collected from their medical charts. Cytokines were quantified by Luminex. RESULTS: Genetic analyses of enrolled patients (n=42) identified 30 patients carrying UBA1 pathogenic variants, with frequencies compatible for postzygotic variants. All patients were male individuals who presented with a late-onset disease (mean 67.5 years; median 67.0 years) characterised by cutaneous lesions (90%), fever (66.7%), pulmonary manifestations (66.7%) and arthritis (53.3%). Macrocytic anaemia and increased erythrocyte sedimentation rate and ferritin were the most relevant analytical abnormalities. Glucocorticoids ameliorated the inflammatory manifestations, but most patients became glucocorticoid-dependent. Positive responses were obtained when targeting the haematopoietic component of the disease with either decitabine or allogeneic haematopoietic stem cell transplantation. Additional analyses detected the UBA1 variants in both haematopoietic and non-haematopoietic tissues. Finally, analysis of circulating cytokines did not identify inflammatory mediators of the disease. CONCLUSION: Thirty patients with adult-onset AID were definitively diagnosed with VEXAS syndrome through genetic analyses. Despite minor interindividual differences, their main characteristics were in concordance with previous reports. We detected for the first time the UBA1 mosaicism in non-haematopoietic tissue, which questions the previous concept of myeloid-restricted mosaicism and may have conceptual consequences for the disease mechanisms.

Characterization of Novel Pathogenic Variants Leading to Caspase-8 Cleavage-Resistant RIPK1-Induced Autoinflammatory Syndrome.

Pathogenic RIPK1 variants have been described as the cause of two different inborn errors of immunity. Biallelic loss-of-function variants cause the recessively inherited RIPK1 deficiency, while monoallelic variants impairing the caspase-8-mediated RIPK1 cleavage provoke a novel autoinflammatory disease (AID) called cleavage-resistant RIPK1-induced autoinflammatory (CRIA) syndrome. The aim of this study was to characterize the pathogenicity of two novel RIPK1 variants located at the cleavage site of caspase-8 detected in patients with dominantly-inherited, early-onset undefined AID. RIPK1 genotyping was performed by Sanger and next-generation sequencing. Clinical and analytical data were collected from medical charts, and in silico and in vitro assays were performed to evaluate the functional consequences. Genetic analyses identified two novel heterozygous RIPK1 variants at the caspase-8 cleavage site (p.Leu321Arg and p.Asp324Gly), which displayed a perfect intrafamilial phenotype-genotype segregation following a dominant inheritance pattern. Structural analyses suggested that these variants disrupt the normal RIPK1 structure, probably making it less accessible to and/or less cleavable by caspase-8. In vitro experiments confirmed that the p.Leu321Arg and p.Asp324Gly RIPK1 variants were resistant to caspase-8-mediated cleavage and induced a constitutive activation of necroptotic pathway in a similar manner that previously characterized RIPK1 variants causing CRIA syndrome. All these results strongly supported the pathogenicity of the two novel RIPK1 variants and the diagnosis of CRIA syndrome in all enrolled patients. Moreover, the evidences here collected expand the phenotypic and genetic diversity of this recently described AID, and provide interesting data about effectiveness of treatments that may benefit future patients.

Characterization of 58 STRs and 94 SNPs with the ForenSeq™ DNA signature prep kit in Mexican-Mestizos from the Monterrey city (Northeast, Mexico).

BACKGROUND: STR allele frequency databases from populations are necessary to take full advantage of the increased power of discrimination offered by massively parallel sequencing (MPS) platforms. MATERIAL AND METHODS: For this reason, we sequenced 58 STRs (aSTRs, X-STRs, and Y-STRs) and 94 identity informative SNPs (iiSNPs) on 105 Mestizo (admixed) individuals from Monterrey City (Northeast, Mexico), with the Primer Set-A of the ForenSeq™ DNA Signature Prep Kit. RESULTS: Most of the STR markers were in Hardy Weinberg equilibrium, with a few exceptions. We found 346 different length-based alleles for these 58 STRs; nevertheless, they became 528 alleles when the sequence was assessed. The combined power of discrimination from autosomal STRs (aSTRs) was -virtually- 100% in both length and sequence-based alleles, while the power of exclusion was 99.9999999976065 and 99.9999999999494%, respectively. Haplotypes based on X-STRs and Y-STRs showed 100% of discriminatory capacity. CONCLUSIONS: These results provide -for the first time- forensic genomic population data from Mexico necessary for interpretation in kinship and criminal analyses.

Severe Autoinflammatory Manifestations and Antibody Deficiency Due to Novel Hypermorphic PLCG2 Mutations.

Autoinflammatory diseases (AIDs) were first described as clinical disorders characterized by recurrent episodes of seemingly unprovoked sterile inflammation. In the past few years, the identification of novel AIDs expanded their phenotypes toward more complex clinical pictures associating vasculopathy, autoimmunity, or immunodeficiency. Herein, we describe two unrelated patients suffering since the neonatal period from a complex disease mainly characterized by severe sterile inflammation, recurrent bacterial infections, and marked humoral immunodeficiency. Whole-exome sequencing detected a novel, de novo heterozygous PLCG2 variant in each patient (p.Ala708Pro and p.Leu845_Leu848del). A clear enhanced PLCγ2 activity for both variants was demonstrated by both ex vivo calcium responses of the patient's B cells to IgM stimulation and in vitro assessment of PLC activity. These data supported the autoinflammation and PLCγ2-associated antibody deficiency and immune dysregulation (APLAID) diagnosis in both patients. Immunological evaluation revealed a severe decrease of immunoglobulins and B cells, especially class-switched memory B cells, with normal T and NK cell counts. Analysis of bone marrow of one patient revealed a reduced immature B cell fraction compared with controls. Additional investigations showed that both PLCG2 variants activate the NLRP3-inflammasome through the alternative pathway instead of the canonical pathway. Collectively, the evidences here shown expand APLAID diversity toward more severe phenotypes than previously reported including dominantly inherited agammaglobulinemia, add novel data about its genetic basis, and implicate the alternative NLRP3-inflammasome activation pathway in the basis of sterile inflammation.

Historical human remains identification through maternal and paternal genetic signatures in a founder population with extensive genealogical record.

OBJECTIVES: We describe a method to identify human remains excavated from unmarked graves in historical Québec cemeteries by combining parental-lineage genetic markers with the whole-population genealogy of Québec contained in the BALSAC database. MATERIALS AND METHODS: The remains of six men were exhumed from four historical cemeteries in the province of Québec, Canada. DNA was extracted from the remains and genotyped to reveal their mitochondrial and Y-chromosome haplotypes, which were compared to a collection of haplotypes of genealogically-anchored modern volunteers. Maternal and paternal genealogies were searched in the BALSAC genealogical record for parental couples matching the mitochondrial and the Y-chromosome haplotypic signatures, to identify candidate sons from whom the remains could have originated. RESULTS: Analysis of the matching genealogies identified the parents of one man inhumed in the cemetery of the investigated parish during its operating time. The candidate individual died in 1833 at the age of 58, a plausible age at death in light of osteological analysis of the remains. DISCUSSION: This study demonstrates the promising potential of coupling genetic information from living individuals to genealogical data in BALSAC to identify historical human remains. If genetic coverage is increased, the genealogical information in BALSAC could enable the identification of 87% of the men (n = 178,435) married in Québec before 1850, with high discriminatory power in most cases since >75% of the parental couples have unique biparental signatures in most regions. Genotyping and identifying Québec's historical human remains are a key to reconstructing the genomes of the founders of Québec and reinhuming archeological remains with a marked grave.

Unexpected relevant role of gene mosaicism in patients with primary immunodeficiency diseases.

BACKGROUND: Postzygotic de novo mutations lead to the phenomenon of gene mosaicism. The 3 main types are called somatic, gonadal, and gonosomal mosaicism, which differ in terms of the body distribution of postzygotic mutations. Mosaicism has been reported occasionally in patients with primary immunodeficiency diseases (PIDs) since the early 1990s, but its real involvement has not been systematically addressed. OBJECTIVE: We sought to investigate the incidence of gene mosaicism in patients with PIDs. METHODS: The amplicon-based deep sequencing method was used in the 3 parts of the study that establish (1) the allele frequency of germline variants (n = 100), (2) the incidence of parental gonosomal mosaicism in families with PIDs with de novo mutations (n = 92), and (3) the incidence of mosaicism in families with PIDs with moderate-to-high suspicion of gene mosaicism (n = 36). Additional investigations evaluated body distribution of postzygotic mutations, their stability over time, and their characteristics. RESULTS: The range of allele frequency (44.1% to 55.6%) was established for germline variants. Those with minor allele frequencies of less than 44.1% were assumed to be postzygotic. Mosaicism was detected in 30 (23.4%) of 128 families with PIDs, with a variable minor allele frequency (0.8% to 40.5%). Parental gonosomal mosaicism was detected in 6 (6.5%) of 92 families with de novo mutations, and a high incidence of mosaicism (63.9%) was detected among families with moderate-to-high suspicion of gene mosaicism. In most analyzed cases mosaicism was found to be both uniformly distributed and stable over time. CONCLUSION: This study represents the largest performed to date to investigate mosaicism in patients with PIDs, revealing that it affects approximately 25% of enrolled families. Our results might have serious consequences regarding treatment and genetic counseling and reinforce the use of next-generation sequencing-based methods in the routine analyses of PIDs.

Great ape genetic diversity and population history.

Most great ape genetic variation remains uncharacterized; however, its study is critical for understanding population history, recombination, selection and susceptibility to disease. Here we sequence to high coverage a total of 79 wild- and captive-born individuals representing all six great ape species and seven subspecies and report 88.8 million single nucleotide polymorphisms. Our analysis provides support for genetically distinct populations within each species, signals of gene flow, and the split of common chimpanzees into two distinct groups: Nigeria-Cameroon/western and central/eastern populations. We find extensive inbreeding in almost all wild populations, with eastern gorillas being the most extreme. Inferred effective population sizes have varied radically over time in different lineages and this appears to have a profound effect on the genetic diversity at, or close to, genes in almost all species. We discover and assign 1,982 loss-of-function variants throughout the human and great ape lineages, determining that the rate of gene loss has not been different in the human branch compared to other internal branches in the great ape phylogeny. This comprehensive catalogue of great ape genome diversity provides a framework for understanding evolution and a resource for more effective management of wild and captive great ape populations.

Natural Selection in the Great Apes.

Natural selection is crucial for the adaptation of populations to their environments. Here, we present the first global study of natural selection in the Hominidae (humans and great apes) based on genome-wide information from population samples representing all extant species (including most subspecies). Combining several neutrality tests we create a multi-species map of signatures of natural selection covering all major types of natural selection. We find that the estimated efficiency of both purifying and positive selection varies between species and is significantly correlated with their long-term effective population size. Thus, even the modest differences in population size among the closely related Hominidae lineages have resulted in differences in their ability to remove deleterious alleles and to adapt to changing environments. Most signatures of balancing and positive selection are species-specific, with signatures of balancing selection more often being shared among species. We also identify loci with evidence of positive selection across several lineages. Notably, we detect signatures of positive selection in several genes related to brain function, anatomy, diet and immune processes. Our results contribute to a better understanding of human evolution by putting the evidence of natural selection in humans within its larger evolutionary context. The global map of natural selection in our closest living relatives is available as an interactive browser at http://tinyurl.com/nf8qmzh.

Mendelian genes for Parkinson's disease contribute to the sporadic forms of the disease.

Parkinson's disease (PD) can be divided into familial (Mendelian) and sporadic forms. A number of causal genes have been discovered for the Mendelian form, which constitutes 10-20% of the total cases. Genome-wide association studies have successfully uncovered a number of susceptibility loci for sporadic cases but those only explain a small fraction (6-7%) of PD heritability. It has been observed that some genes that confer susceptibility to PD through common risk variants also contain rare causing mutations for the Mendelian forms of the disease. These results suggest a possible functional link between Mendelian and sporadic PD and led us to investigate the role that rare and low-frequency variants could have on the sporadic form. Through a targeting approach, we have resequenced at 49× coverage the exons and regulatory regions of 38 genes (including Mendelian and susceptibility PD genes) in 249 sporadic PD patients and 145 unrelated controls of European origin. Unlike susceptibility genes, Mendelian genes show a clear general enrichment of rare functional variants in PD cases, observed directly as well as with Tajima's D statistic and several collapsing methods. Our findings suggest that rare variation on PD Mendelian genes may have a role in the sporadic forms of the disease.

Y-chromosomal sequences of diverse Indian populations and the ancestry of the Andamanese.

We present 42 new Y-chromosomal sequences from diverse Indian tribal and non-tribal populations, including the Jarawa and Onge from the Andaman Islands, which are analysed within a calibrated Y-chromosomal phylogeny incorporating South Asian (in total 305 individuals) and worldwide (in total 1286 individuals) data from the 1000 Genomes Project. In contrast to the more ancient ancestry in the South than in the North that has been claimed, we detected very similar coalescence times within Northern and Southern non-tribal Indian populations. A closest neighbour analysis in the phylogeny showed that Indian populations have an affinity towards Southern European populations and that the time of divergence from these populations substantially predated the Indo-European migration into India, probably reflecting ancient shared ancestry rather than the Indo-European migration, which had little effect on Indian male lineages. Among the tribal populations, the Birhor (Austro-Asiatic-speaking) and Irula (Dravidian-speaking) are the nearest neighbours of South Asian non-tribal populations, with a common origin in the last few millennia. In contrast, the Riang (Tibeto-Burman-speaking) and Andamanese have their nearest neighbour lineages in East Asia. The Jarawa and Onge shared haplogroup D lineages with each other within the last ~7000 years, but had diverged from Japanese haplogroup D Y-chromosomes ~53000 years ago, most likely by a split from a shared ancestral population. This analysis suggests that Indian populations have complex ancestry which cannot be explained by a single expansion model.

Detection of genomic rearrangements from targeted resequencing data in Parkinson's disease patients.

BACKGROUND: The analysis of coverage depth in next-generation sequencing data allows the detection of gene dose alterations. We explore the frequency of such structural events in a Spanish cohort of sporadic PD cases. METHODS: Gene dose alterations were detected with the eXome-Hidden Markov Model (XHMM) software from depth of coverage in resequencing data available for 38 Mendelian and other risk PD loci in 394 individuals (249 cases and 145 controls) and subsequently validated by quantitative PCR. RESULTS: We identified 10 PD patients with exon dosage alterations in PARK2, GBA-GBAP1, and DJ1. Additional functional variants, including 2 novel nonsense mutations (p.Arg1552Ter in LRRK2 and p.Trp90Ter in PINK1), were confirmed by Sanger sequencing. This combined approach disclosed the genetic cause of 12 PD cases. CONCLUSIONS: Gene dose alterations related to PD can be correctly identified from targeting resequencing data. This approach substantially improves the detection rate of cases with causal genetic alterations. © 2016 The Authors. Movement Disorders published by Wiley Periodicals, Inc. on behalf of International Parkinson and Movement Disorder Society.

Rare allelic forms of PRDM9 associated with childhood leukemogenesis.

One of the most rapidly evolving genes in humans, PRDM9, is a key determinant of the distribution of meiotic recombination events. Mutations in this meiotic-specific gene have previously been associated with male infertility in humans and recent studies suggest that PRDM9 may be involved in pathological genomic rearrangements. In studying genomes from families with children affected by B-cell precursor acute lymphoblastic leukemia (B-ALL), we characterized meiotic recombination patterns within a family with two siblings having hyperdiploid childhood B-ALL and observed unusual localization of maternal recombination events. The mother of the family carries a rare PRDM9 allele, potentially explaining the unusual patterns found. From exomes sequenced in 44 additional parents of children affected with B-ALL, we discovered a substantial and significant excess of rare allelic forms of PRDM9. The rare PRDM9 alleles are transmitted to the affected children in half the cases; nonetheless there remains a significant excess of rare alleles among patients relative to controls. We successfully replicated this latter observation in an independent cohort of 50 children with B-ALL, where we found an excess of rare PRDM9 alleles in aneuploid and infant B-ALL patients. PRDM9 variability in humans is thought to influence genomic instability, and these data support a potential role for PRDM9 variation in risk of acquiring aneuploidies or genomic rearrangements associated with childhood leukemogenesis.

An assessment of a massively parallel sequencing approach for the identification of individuals from mass graves of the Spanish Civil War (1936-1939).

Next-generation sequencing technologies have opened new opportunities in forensic genetics. Here, we assess the applicability and performance of the MiSeq FGx™ & ForenSeq™ DNA Signature Prep Kit (Illumina) for the identification of individuals from the mass graves of the Spanish Civil War (1936-1939). The main limitations for individual identification are the low number of possible first-degree living relatives and the high levels of DNA degradation reported in previous studies. Massively parallel sequencing technologies enabling the analysis of hundreds of regions and prioritizing short length amplicons constitute a promising tool for this kind of approaches. In this study, we first explore the power of this new technology to detect first- and second-degree kinship given different scenarios of DNA degradation. Second, we specifically assess its performance in a set of low DNA input samples previously analyzed with CE technologies. We conclude that this methodology will allow identification of up to second-degree relatives, even in situations with low sequencing performance and important levels of allele drop-out; it is thus a technology that resolves previous drawbacks and that will allow a successful approximation to the identification of remains.

Whole-exome sequencing reveals a rapid change in the frequency of rare functional variants in a founding population of humans.

Whole-exome or gene targeted resequencing in hundreds to thousands of individuals has shown that the majority of genetic variants are at low frequency in human populations. Rare variants are enriched for functional mutations and are expected to explain an important fraction of the genetic etiology of human disease, therefore having a potential medical interest. In this work, we analyze the whole-exome sequences of French-Canadian individuals, a founder population with a unique demographic history that includes an original population bottleneck less than 20 generations ago, followed by a demographic explosion, and the whole exomes of French individuals sampled from France. We show that in less than 20 generations of genetic isolation from the French population, the genetic pool of French-Canadians shows reduced levels of diversity, higher homozygosity, and an excess of rare variants with low variant sharing with Europeans. Furthermore, the French-Canadian population contains a larger proportion of putatively damaging functional variants, which could partially explain the increased incidence of genetic disease in the province. Our results highlight the impact of population demography on genetic fitness and the contribution of rare variants to the human genetic variation landscape, emphasizing the need for deep cataloguing of genetic variants by resequencing worldwide human populations in order to truly assess disease risk.

Population and genomic lessons from genetic analysis of two Indian populations.

Indian demographic history includes special features such as founder effects, interpopulation segregation, complex social structure with a caste system and elevated frequency of consanguineous marriages. It also presents a higher frequency for some rare mendelian disorders and in the last two decades increased prevalence of some complex disorders. Despite the fact that India represents about one-sixth of the human population, deep genetic studies from this terrain have been scarce. In this study, we analyzed high-density genotyping and whole-exome sequencing data of a North and a South Indian population. Indian populations show higher differentiation levels than those reported between populations of other continents. In this work, we have analyzed its consequences, by specifically assessing the transferability of genetic markers from or to Indian populations. We show that there is limited genetic marker portability from available genetic resources such as HapMap or the 1,000 Genomes Project to Indian populations, which also present an excess of private rare variants. Conversely, tagSNPs show a high level of portability between the two Indian populations, in contrast to the common belief that North and South Indian populations are genetically very different. By estimating kinship from mates and consanguinity in our data from trios, we also describe different patterns of assortative mating and inbreeding in the two populations, in agreement with distinct mating preferences and social structures. In addition, this analysis has allowed us to describe genomic regions under recent adaptive selection, indicating differential adaptive histories for North and South Indian populations. Our findings highlight the importance of considering demography for design and analysis of genetic studies, as well as the need for extending human genetic variation catalogs to new populations and particularly to those with particular demographic histories.

Exome sequencing identifies mutations in the gene TTC7A in French-Canadian cases with hereditary multiple intestinal atresia.

BACKGROUND: Congenital multiple intestinal atresia (MIA) is a severe, fatal neonatal disorder, involving the occurrence of obstructions in the small and large intestines ultimately leading to organ failure. Surgical interventions are palliative but do not provide long-term survival. Severe immunodeficiency may be associated with the phenotype. A genetic basis for MIA is likely. We had previously ascertained a cohort of patients of French-Canadian origin, most of whom were deceased as infants or in utero. The goal of the study was to identify the molecular basis for the disease in the patients of this cohort. METHODS: We performed whole exome sequencing on samples from five patients of four families. Validation of mutations and familial segregation was performed using standard Sanger sequencing in these and three additional families with deceased cases. Exon skipping was assessed by reverse transcription-PCR and Sanger sequencing. RESULTS: Five patients from four different families were each homozygous for a four base intronic deletion in the gene TTC7A, immediately adjacent to a consensus GT splice donor site. The deletion was demonstrated to have deleterious effects on splicing causing the skipping of the attendant upstream coding exon, thereby leading to a predicted severe protein truncation. Parents were heterozygous carriers of the deletion in these families and in two additional families segregating affected cases. In a seventh family, an affected case was compound heterozygous for the same 4bp deletion and a second missense mutation p.L823P, also predicted as pathogenic. No other sequenced genes possessed deleterious variants explanatory for all patients in the cohort. Neither mutation was seen in a large set of control chromosomes. CONCLUSIONS: Based on our genetic results, TTC7A is the likely causal gene for MIA.

A population genetic approach to mapping neurological disorder genes using deep resequencing.

Deep resequencing of functional regions in human genomes is key to identifying potentially causal rare variants for complex disorders. Here, we present the results from a large-sample resequencing (n  =  285 patients) study of candidate genes coupled with population genetics and statistical methods to identify rare variants associated with Autism Spectrum Disorder and Schizophrenia. Three genes, MAP1A, GRIN2B, and CACNA1F, were consistently identified by different methods as having significant excess of rare missense mutations in either one or both disease cohorts. In a broader context, we also found that the overall site frequency spectrum of variation in these cases is best explained by population models of both selection and complex demography rather than neutral models or models accounting for complex demography alone. Mutations in the three disease-associated genes explained much of the difference in the overall site frequency spectrum among the cases versus controls. This study demonstrates that genes associated with complex disorders can be mapped using resequencing and analytical methods with sample sizes far smaller than those required by genome-wide association studies. Additionally, our findings support the hypothesis that rare mutations account for a proportion of the phenotypic variance of these complex disorders.