African Pygmies practicing a mobile hunter-gatherer lifestyle are phenotypically and genetically diverged from other anatomically modern humans, and they likely experienced strong selective pressures due to their unique lifestyle in the Central African rainforest. To identify genomic targets of adaptation, we sequenced the genomes of four Biaka Pygmies from the Central African Republic and jointly analyzed these data with the genome sequences of three Baka Pygmies from Cameroon and nine Yoruba famers. To account for the complex demographic history of these populations that includes both isolation and gene flow, we fit models using the joint allele frequency spectrum and validated them using independent approaches. Our two best-fit models both suggest ancient divergence between the ancestors of the farmers and Pygmies, 90,000 or 150,000 yr ago. We also find that bidirectional asymmetric gene flow is statistically better supported than a single pulse of unidirectional gene flow from farmers to Pygmies, as previously suggested. We then applied complementary statistics to scan the genome for evidence of selective sweeps and polygenic selection. We found that conventional statistical outlier approaches were biased toward identifying candidates in regions of high mutation or low recombination rate. To avoid this bias, we assigned P-values for candidates using whole-genome simulations incorporating demography and variation in both recombination and mutation rates. We found that genes and gene sets involved in muscle development, bone synthesis, immunity, reproduction, cell signaling and development, and energy metabolism are likely to be targets of positive natural selection in Western African Pygmies or their recent ancestors.

Although population-level genomic sequence data have been gathered extensively for humans, similar data from our closest living relatives are just beginning to emerge. Examination of genomic variation within great apes offers many opportunities to increase our understanding of the forces that have differentially shaped the evolutionary history of hominid taxa. Here, we expand upon the work of the Great Ape Genome Project by analyzing medium to high coverage whole-genome sequences from 14 western lowland gorillas (Gorilla gorilla gorilla), 2 eastern lowland gorillas (G. beringei graueri), and a single Cross River individual (G. gorilla diehli). We infer that the ancestors of western and eastern lowland gorillas diverged from a common ancestor approximately 261 ka, and that the ancestors of the Cross River population diverged from the western lowland gorilla lineage approximately 68 ka. Using a diffusion approximation approach to model the genome-wide site frequency spectrum, we infer a history of western lowland gorillas that includes an ancestral population expansion of 1.4-fold around 970 ka and a recent 5.6-fold contraction in population size 23 ka. The latter may correspond to a major reduction in African equatorial forests around the Last Glacial Maximum. We also analyze patterns of variation among western lowland gorillas to identify several genomic regions with strong signatures of recent selective sweeps. We find that processes related to taste, pancreatic and saliva secretion, sodium ion transmembrane transport, and cardiac muscle function are overrepresented in genomic regions predicted to have experienced recent positive selection.

Gibbons are small arboreal apes that display an accelerated rate of evolutionary chromosomal rearrangement and occupy a key node in the primate phylogeny between Old World monkeys and great apes. Here we present the assembly and analysis of a northern white-cheeked gibbon (Nomascus leucogenys) genome. We describe the propensity for a gibbon-specific retrotransposon (LAVA) to insert into chromosome segregation genes and alter transcription by providing a premature termination site, suggesting a possible molecular mechanism for the genome plasticity of the gibbon lineage. We further show that the gibbon genera (Nomascus, Hylobates, Hoolock and Symphalangus) experienced a near-instantaneous radiation ∼5 million years ago, coincident with major geographical changes in southeast Asia that caused cycles of habitat compression and expansion. Finally, we identify signatures of positive selection in genes important for forelimb development (TBX5) and connective tissues (COL1A1) that may have been involved in the adaptation of gibbons to their arboreal habitat.

Individuals with severe, sporadic disorders of infantile onset represent an important class of disease for which discovery of the underlying genetic architecture is not amenable to traditional genetic analysis. Full-genome sequencing of affected individuals and their parents provides a powerful alternative strategy for gene discovery. We performed whole-genome sequencing (WGS) on a family quartet containing an affected proband and her unaffected parents and sibling. The 15-year-old female proband had a severe epileptic encephalopathy consisting of early-onset seizures, features of autism, intellectual disability, ataxia, and sudden unexplained death in epilepsy. We discovered a de novo heterozygous missense mutation (c.5302A>G [p.Asn1768Asp]) in the voltage-gated sodium-channel gene SCN8A in the proband. This mutation alters an evolutionarily conserved residue in Nav1.6, one of the most abundant sodium channels in the brain. Analysis of the biophysical properties of the mutant channel demonstrated a dramatic increase in persistent sodium current, incomplete channel inactivation, and a depolarizing shift in the voltage dependence of steady-state fast inactivation. Current-clamp analysis in hippocampal neurons transfected with p.Asn1768Asp channels revealed increased spontaneous firing, paroxysmal-depolarizing-shift-like complexes, and an increased firing frequency, consistent with a dominant gain-of-function phenotype in the heterozygous proband. This work identifies SCN8A as the fifth sodium-channel gene to be mutated in epilepsy and demonstrates the value of WGS for the identification of pathogenic mutations causing severe, sporadic neurological disorders.

Oculocutaneous albinism (OCA) is a genetically heterogeneous disorder. There are four known types of OCA: OCA1-OCA4. The clinical manifestations of all types of OCA include skin and hair hypopigmentation and visual impairment. Although there are a few documented observations of high frequency of albinism among Native Americans, including the Hopi, Zuni, Kuna, Jemez, Laguna, San Juan, and Navajo, no causative molecular defect has been previously reported. In the present study, we show that albinism in one Native American population, the Navajo, is caused by a LINE-mediated 122.5-kilobase deletion of the P gene, thus demonstrating that albinism in this population is OCA2. This deletion appears to be Navajo specific, because this allele was not detected in 34 other individuals with albinism who listed other Native American origins, nor has it been reported in any other ethnic group. The molecular characterization of this deletion allele allowed us to design a three-primer polymerase chain reaction system to estimate the carrier frequency in the Navajo population by screening 134 unrelated normally pigmented Navajos. The carrier frequency was found to be approximately 4.5%. The estimated prevalence of OCA2 in Navajos is between approximately 1 per 1,500 and 1 per 2,000. We further estimate that this mutation originated 400-1,000 years ago from a single founder.

Signals of archaic admixture have been identified through comparisons of the draft Neanderthal and Denisova genomes with those of living humans. Studies of individual loci contributing to these genome-wide average signals are required for characterization of the introgression process and investigation of whether archaic variants conferred an adaptive advantage to the ancestors of contemporary human populations. However, no definitive case of adaptive introgression has yet been described. Here we provide a DNA sequence analysis of the innate immune gene STAT2 and show that a haplotype carried by many Eurasians (but not sub-Saharan Africans) has a sequence that closely matches that of the Neanderthal STAT2. This haplotype, referred to as N, was discovered through a resequencing survey of the entire coding region of STAT2 in a global sample of 90 individuals. Analyses of publicly available complete genome sequence data show that haplotype N shares a recent common ancestor with the Neanderthal sequence (~80 thousand years ago) and is found throughout Eurasia at an average frequency of ~5%. Interestingly, N is found in Melanesian populations at ~10-fold higher frequency (~54%) than in Eurasian populations. A neutrality test that controls for demography rejects the hypothesis that a variant of N rose to high frequency in Melanesia by genetic drift alone. Although we are not able to pinpoint the precise target of positive selection, we identify nonsynonymous mutations in ERBB3, ESYT1, and STAT2-all of which are part of the same 250 kb introgressive haplotype-as good candidates.

Among patients with intractable epilepsy, the most commonly performed surgical procedure is craniotomy for amygdalohippocampectomy (AH). Stereotactic laser amygdalohippocampotomy (SLAH) has also been recently employed as a minimally invasive treatment for intractable temporal lobe epilepsy (TLE). Among patients treated with AH and SLAH approximately 65% and 54% of patients become seizure-free, respectively. Therefore, selection criteria for surgical candidates with improved prognostic value for post-operative seizure-free outcome are greatly needed. In this study, we perform RNA sequencing (RNA-Seq) on whole blood leukocyte samples taken from 16 patients with intractable TLE prior to SLAH to test the hypothesis that pre-operative leukocyte RNA expression profiles are prognostic for post-operative seizure outcome. Multidimensional scaling analysis of the RNA expression data indicated separate clustering of patients with seizure free (SF) and non-seizure-free (NSF) outcomes. Differential expression (DE) analysis performed on SF versus NSF groups revealed 24 significantly differentially expressed genes (≥2.0-fold change, p-value < 0.05, FDR <0.05). Network and pathway analyses identified differential activation of pathways involved in lipid metabolism, morphology of oligodendrocytes, inflammatory response, and development of astrocytes. These results suggest that pre-operative leukocyte expression profiles have prognostic value for seizure outcome following SLAH.

Epilepsy is a common neurological disease; however, few if any of the currently marketed antiseizure medications prevent or cure epilepsy. Discovery of pathological processes in the early stages of epileptogenesis has been challenging given the common use of preclinical models that induce seizures in physiologically normal animals. Moreover, despite known sex dimorphism in neurological diseases, females are rarely included in preclinical epilepsy models.

Admixture between long-separated populations is a defining feature of the genomes of many species. The mosaic block structure of admixed genomes can provide information about past contact events, including the time and extent of admixture. Here, we describe an improved wavelet-based technique that better characterizes ancestry block structure from observed genomic patterns. principal components analysis is first applied to genomic data to identify the primary population structure, followed by wavelet decomposition to develop a new characterization of local ancestry information along the chromosomes. For testing purposes, this method is applied to human genome-wide genotype data from Indonesia, as well as virtual genetic data generated using genome-scale sequential coalescent simulations under a wide range of admixture scenarios. Time of admixture is inferred using an approximate Bayesian computation framework, providing robust estimates of both admixture times and their associated levels of uncertainty. Crucially, we demonstrate that this revised wavelet approach, which we have released as the R package adwave, provides improved statistical power over existing wavelet-based techniques and can be used to address a broad range of admixture questions.

Timor, an eastern Indonesian island linking mainland Asia with Australia and the Pacific world, had a complex history, including its role as a contact zone between two language families (Austronesian and Trans-New Guinean), as well as preserving elements of a rich Austronesian cultural heritage, such as matrilocal marriage practices. Using an array of biparental (autosomal and X-chromosome single-nucleotide polymorphisms) and uniparental markers (Y chromosome and mitochondrial DNA), we reconstruct a broad genetic profile of Timorese in the Belu regency of West Timor, including the traditional princedom of Wehali, focusing on the effects of cultural practices, such as language and social change, on patterns of genetic diversity. Sex-linked data highlight the different histories and social pressures experienced by women and men. Measures of diversity and population structure show that Timorese men had greater local mobility than women, as expected in matrilocal communities, where women remain in their natal village, whereas men move to the home village of their wife. Reaching further back in time, maternal loci (mitochondrial DNA and the X chromosome) are dominated by lineages with immigrant Asian origins, whereas paternal loci (Y chromosome) tend to exhibit lineages of the earliest settlers in the eastern Indonesian region. The dominance of Asian female lineages is especially apparent in the X chromosome compared with the autosomes, suggesting that women played a paramount role during and after the period of Asian immigration into Timor, perhaps driven by the matrilocal marriage practices of expanding Austronesian communities.

Recently, de novo SCN8A missense mutations have been identified as a rare dominant cause of epileptic encephalopathies (EIEE13). Functional studies on the first described case demonstrated gain-of-function effects of the mutation. We describe a novel de novo mutation of SCN8A in a patient with epileptic encephalopathy, and functional characterization of the mutant protein.

Single nucleotide polymorphisms (SNPs) in commercial arrays have often been discovered in a small number of samples from selected populations. This ascertainment skews patterns of nucleotide diversity and affects population genetic inferences. We propose a demographic inference pipeline that explicitly models the SNP discovery protocol in an Approximate Bayesian Computation (ABC) framework. We simulated genomic regions according to a demographic model incorporating parameters for the divergence of three well-characterized HapMap populations and recreated the SNP distribution of a commercial array by varying the number of haploid samples and the allele frequency cut-off in the given regions. We then calculated summary statistics obtained from both the ascertained and genomic data and inferred ascertainment and demographic parameters. We implemented our pipeline to study the admixture process that gave rise to the present-day Mexican population. Our estimate of the time of admixture is closer to the historical dates than those in previous works which did not consider ascertainment bias. Although the use of whole genome sequences for demographic inference is becoming the norm, there are still underrepresented areas of the world from where only SNP array data are available. Our inference framework is applicable to those cases and will help with the demographic inference.

We report the discovery of an African American Y chromosome that carries the ancestral state of all SNPs that defined the basal portion of the Y chromosome phylogenetic tree. We sequenced ∼240 kb of this chromosome to identify private, derived mutations on this lineage, which we named A00. We then estimated the time to the most recent common ancestor (TMRCA) for the Y tree as 338 thousand years ago (kya) (95% confidence interval = 237-581 kya). Remarkably, this exceeds current estimates of the mtDNA TMRCA, as well as those of the age of the oldest anatomically modern human fossils. The extremely ancient age combined with the rarity of the A00 lineage, which we also find at very low frequency in central Africa, point to the importance of considering more complex models for the origin of Y chromosome diversity. These models include ancient population structure and the possibility of archaic introgression of Y chromosomes into anatomically modern humans. The A00 lineage was discovered in a large database of consumer samples of African Americans and has not been identified in traditional hunter-gatherer populations from sub-Saharan Africa. This underscores how the stochastic nature of the genealogical process can affect inference from a single locus and warrants caution during the interpretation of the geographic location of divergent branches of the Y chromosome phylogenetic tree for the elucidation of human origins.

Indonesia, an island nation as large as continental Europe, hosts a sizeable proportion of global human diversity, yet remains surprisingly undercharacterized genetically. Here, we substantially expand on existing studies by reporting genome-scale data for nearly 500 individuals from 25 populations in Island Southeast Asia, New Guinea, and Oceania, notably including previously unsampled islands across the Indonesian archipelago. We use high-resolution analyses of haplotype diversity to reveal fine detail of regional admixture patterns, with a particular focus on the Holocene. We find that recent population history within Indonesia is complex, and that populations from the Philippines made important genetic contributions in the early phases of the Austronesian expansion. Different, but interrelated processes, acted in the east and west. The Austronesian migration took several centuries to spread across the eastern part of the archipelago, where genetic admixture postdates the archeological signal. As with the Neolithic expansion further east in Oceania and in Europe, genetic mixing with local inhabitants in eastern Indonesia lagged behind the arrival of farming populations. In contrast, western Indonesia has a more complicated admixture history shaped by interactions with mainland Asian and Austronesian newcomers, which for some populations occurred more than once. Another layer of complexity in the west was introduced by genetic contact with South Asia and strong demographic events in isolated local groups.

Despite recent advances in our understanding of synaptic transmission associated with epileptogenesis, the molecular mechanisms that control seizure frequency in patients with temporal lobe epilepsy (TLE) remain obscure. RNA-Seq was performed on hippocampal tissue resected from 12 medically intractable TLE patients with pre-surgery seizure frequencies ranging from 0.33 to 120 seizures per month. Differential expression (DE) analysis of individuals with low (LSF, mean = 4 seizure/month) versus high (HSF, mean = 60 seizures/month) seizure frequency identified 979 genes with ≥2-fold change in transcript abundance (FDR-adjusted p-value ö0.05). Comparisons with post-mortem controls revealed a large number of downregulated genes in the HSF (1676) versus LSF (399) groups. More than 50 signaling pathways were inferred to be deactivated or activated, with Signal Transduction as the central hub in the pathway network. While neuroinflammation pathways were activated in both groups, key neuronal system pathways were systematically deactivated in the HSF group, including calcium, CREB and Opioid signaling. We also infer that enhanced expression of a signaling cascade promoting synaptic downscaling may have played a key role in maintaining a higher seizure threshold in the LSF cohort. These results suggest that therapeutic approaches targeting synaptic scaling pathways may aid in the treatment of seizures in TLE.

SCN8A encephalopathy is a severe, early-onset epilepsy disorder resulting from de novo gain-of-function mutations in the voltage-gated sodium channel Nav1.6. To identify the effects of this disorder on mRNA expression, RNA-seq was performed on brain tissue from a knock-in mouse expressing the patient mutation p.Asn1768Asp (N1768D). RNA was isolated from forebrain, cerebellum, and brainstem both before and after seizure onset, and from age-matched wildtype littermates. Altered transcript profiles were observed only in forebrain and only after seizures. The abundance of 50 transcripts increased more than 3-fold and 15 transcripts decreased more than 3-fold after seizures. The elevated transcripts included two anti-convulsant neuropeptides and more than a dozen genes involved in reactive astrocytosis and response to neuronal damage. There was no change in the level of transcripts encoding other voltage-gated sodium, potassium or calcium channels. Reactive astrocytosis was observed in the hippocampus of mutant mice after seizures. There is considerable overlap between the genes affected in this genetic model of epilepsy and those altered by chemically induced seizures, traumatic brain injury, ischemia, and inflammation. The data support the view that gain-of-function mutations of SCN8A lead to pathogenic alterations in brain function contributing to encephalopathy.

It is commonly thought that human genetic diversity in non-African populations was shaped primarily by an out-of-Africa dispersal 50-100 thousand yr ago (kya). Here, we present a study of 456 geographically diverse high-coverage Y chromosome sequences, including 299 newly reported samples. Applying ancient DNA calibration, we date the Y-chromosomal most recent common ancestor (MRCA) in Africa at 254 (95% CI 192-307) kya and detect a cluster of major non-African founder haplogroups in a narrow time interval at 47-52 kya, consistent with a rapid initial colonization model of Eurasia and Oceania after the out-of-Africa bottleneck. In contrast to demographic reconstructions based on mtDNA, we infer a second strong bottleneck in Y-chromosome lineages dating to the last 10 ky. We hypothesize that this bottleneck is caused by cultural changes affecting variance of reproductive success among males.

A major unanswered question in the evolution of Homo sapiens is when anatomically modern human populations began to expand: was demographic growth associated with the invention of particular technologies or behavioral innovations by hunter-gatherers in the Late Pleistocene, or with the acquisition of farming in the Neolithic?

Dravet syndrome (DS) is a rare, devastating form of childhood epilepsy that is often associated with mutations in the voltage-gated sodium channel gene, SCN1A. There is considerable variability in expressivity within families, as well as among individuals carrying the same primary mutation, suggesting that clinical outcome is modulated by variants at other genes. To identify modifier gene variants that contribute to clinical outcome, we sequenced the exomes of 22 individuals at both ends of a phenotype distribution (i.e., mild and severe cognitive condition). We controlled for variation associated with different mutation types by limiting inclusion to individuals with a de novo truncation mutation resulting in SCN1A haploinsufficiency. We performed tests aimed at identifying 1) single common variants that are enriched in either phenotypic group, 2) sets of common or rare variants aggregated in and around genes associated with clinical outcome, and 3) rare variants in 237 candidate genes associated with neuronal excitability. While our power to identify enrichment of a common variant in either phenotypic group is limited as a result of the rarity of mild phenotypes in individuals with SCN1A truncation variants, our top candidates did not map to functional regions of genes, or in genes that are known to be associated with neurological pathways. In contrast, we found a statistically-significant excess of rare variants predicted to be damaging and of small effect size in genes associated with neuronal excitability in severely affected individuals. A KCNQ2 variant previously associated with benign neonatal seizures is present in 3 of 12 individuals in the severe category. To compare our results with the healthy population, we performed a similar analysis on whole exome sequencing data from 70 Japanese individuals in the 1000 genomes project. Interestingly, the frequency of rare damaging variants in the same set of neuronal excitability genes in healthy individuals is nearly as high as in severely affected individuals. Rather than a single common gene/variant modifying clinical outcome in SCN1A-related epilepsies, our results point to the cumulative effect of rare variants with little to no measurable phenotypic effect (i.e., typical genetic background) unless present in combination with a disease-causing truncation mutation in SCN1A.

De novo mutations of the voltage-gated sodium channel gene SCN8A have recently been recognized as a cause of epileptic encephalopathy, which is characterized by refractory seizures with developmental delay and cognitive disability. We previously described the heterozygous SCN8A missense mutation p.Asn1768Asp in a child with epileptic encephalopathy that included seizures, ataxia, and sudden unexpected death in epilepsy (SUDEP). The mutation results in increased persistent sodium current and hyperactivity of transfected neurons. We have characterized a knock-in mouse model expressing this dominant gain-of-function mutation to investigate the pathology of the altered channel in vivo. The mutant channel protein is stable in vivo. Heterozygous Scn8a(N1768D/+) mice exhibit seizures and SUDEP, confirming the causality of the de novo mutation in the proband. Using video/EEG analysis, we detect ictal discharges that coincide with convulsive seizures and myoclonic jerks. Prior to seizure onset, heterozygous mutants are not defective in motor learning or fear conditioning, but do exhibit mild impairment of motor coordination and social discrimination. Homozygous mutant mice exhibit earlier seizure onset than heterozygotes and more rapid progression to death. Analysis of the intermediate phenotype of functionally hemizygous Scn8a(N1768D/-) mice indicates that severity is increased by a double dose of mutant protein and reduced by the presence of wild-type protein. Scn8a(N1768D) mutant mice provide a model of epileptic encephalopathy that will be valuable for studying the in vivo effects of hyperactive Nav1.6 and the response to therapeutic interventions.