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.

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.

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.

Prostate cancer (PCa) is the most common male cancer and the second leading cause of cancer death in United States men. Controversy continues over the effectiveness of prostate-specific antigen (PSA) for distinguishing aggressive from indolent PCa. There is a critical need for more specific and sensitive biomarkers to detect and distinguish low- versus high-risk PCa cases. Discovery metabolomics were performed utilizing ultra-performance liquid chromatography-tandem mass spectrometry (UPLC-MS) on plasma samples from 159 men with treatment naïve prostate cancer participating in the North Carolina-Louisiana PCa Project to determine if there were metabolites associated with aggressive PCa. Thirty-five identifiable plasma small molecules were associated with PCa aggressiveness, 15 of which were sphingolipids; nine common molecules were present in both African-American and European-American men. The molecules most associated with PCa aggressiveness were glycosphingolipids; levels of trihexosylceramide and tetrahexosylceramide were most closely associated with high-aggressive PCa. The Cancer Genome Atlas was queried to determine gene alterations within glycosphingolipid metabolism that are associated with PCa and other cancers. Genes that encode enzymes associated with the metabolism of glycosphingolipids were altered in 12% of PCa and >30% of lung, uterine, and ovarian cancers. These data suggest that the identified plasma (glyco)sphingolipids should be further validated for their association with aggressive PCa, suggesting that specific sphingolipids may be included in a diagnostic signature for PCa.

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.

Widely considered probiotic organisms, Bifidobacteria are common inhabitants of the alimentary tract of animals including insects. Bifidobacteria identified from the honey bee are found in larval guts and throughout the alimentary tract, but attain their greatest abundance in the adult hind gut. To further understand the role of Bifidobacteria in honey bees, we sequenced two strains of Bifidobacterium cultured from different alimentary tract environments and life stages.

Microorganisms developing in the liner of the spent fuel pool (SFP) and the fuel transfer channel (FTC) of a Nuclear Power Plant (NPP) can form high radiation resistant biofilms and cause corrosion. Due to difficulties and limitations to obtain large samples from SFP and FTC, cotton swabs were used to collect the biofilm from the wall of these installations. Molecular characterization was performed using massively parallel sequencing to obtain a taxonomic and functional gene classification. Also, samples from the drainage system were evaluated because microorganisms may travel over the 12-meter column of the pool water of the Brazilian Nuclear Power Plant (Angra1), which has been functioning since 1985. Regardless of the treatment of the pool water, our data reveal the unexpected presence of Fungi (Basidiomycota and Ascomycota) as the main contaminators of the SFP and FTC. Ustilaginomycetes (Basidiomycota) was the major class contributor (70%) in the SFP and FTC reflecting the little diversity in these sites; nevertheless, Proteobacteria, Actinobacteria, Firmicutes (Bacilli) were present in small proportions. Mapping total reads against six fungal reference genomes indicate that there is, in fact, a high abundance of fungal sequences in samples collected from SFP and FTC. Analysis of the ribosomal internal transcribed spacer (ITS) 1 and 2 regions and the protein found in the mitochondria of eukaryotic cells, cytochrome b (cytb) grouped our sample fungi in the clade 7 as Ustilago and Pseudozyma. In contrast, in the drainage system, Alphaproteobacteria were present in high abundances (55%). The presence of Sphingopyxis, Mesorhizobium, Erythrobacter, Sphingomonas, Novosphingobium, Sphingobium, Chelativorans, Oceanicaulis, Acidovorax, and Cyanobacteria was observed. Based on genomic annotation data, the assessment of the biological function found a higher proportion of protein-coding sequences related to respiration and protein metabolism in SFP and FTC samples. The knowledge of this biological inventory present in the system may contribute to further studies of potential microorganisms that might be useful for bioremediation of nuclear waste.

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.

Despite the benefits of quantitative data generated by massively parallel sequencing, resolving mitotypes from mixtures occurring in certain ratios remains challenging. In this study, a bioinformatic mixture deconvolution method centered on population-based phasing was developed and validated. The method was first tested on 270 in silico two-person mixtures varying in mixture proportions. An assortment of external reference panels containing information on haplotypic variation (from similar and different haplogroups) was leveraged to assess the effect of panel composition on phasing accuracy. Building on these simulations, mitochondrial genomes from the Human Mitochondrial DataBase were sourced to populate the panels and key parameter values were identified by deconvolving an additional 7290 in silico two-person mixtures. Finally, employing an optimized reference panel and phasing parameters, the approach was validated with in vitro two-person mixtures with differing proportions. Deconvolution was most accurate when the haplotypes in the mixture were similar to haplotypes present in the reference panel and when the mixture ratios were neither highly imbalanced nor subequal (e.g., 4:1). Overall, errors in haplotype estimation were largely bounded by the accuracy of the mixture's genotype results. The proposed framework is the first available approach that automates the reconstruction of complete individual mitotypes from mixtures, even in ratios that have traditionally been considered problematic.

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.

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.

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.

Quantitative genomic and proteomic evaluation of human latent fingerprint depositions represents a challenge within the forensic field, due to the high variability in the amount of DNA and protein initially deposited. To better assess recovery techniques for touch depositions, we present a method to produce simple and customizable artificial fingerprints. These artificial fingerprint samples include the primary components of a typical latent fingerprint, specifically sebaceous fluid, eccrine perspiration, extracellular DNA, and proteinaceous epidermal skin material (i.e., shed skin cells). A commercially available emulsion of sebaceous and eccrine perspiration material provides a chemically-relevant suspension solution for fingerprint deposition, simplifying artificial fingerprint production. Extracted human genomic DNA is added to accurately mimic the extracellular DNA content of a typical latent print and comparable DNA yields are recovered from the artificial prints relative to human prints across surface types. Capitalizing on recent advancements in the use of protein sequence identification for human forensic analysis, these samples also contain a representative quantity of protein, originating from epidermal skin cells collected from the fingers and palms of volunteers. Proteomic sequencing by liquid chromatography-tandem mass spectrometry (LC-MS/MS) analysis indicates a high level of protein overlap between artificial and latent prints. Data are available via ProteomeXchange with identifier PXD015445. By including known quantities of DNA and protein into each artificial print, this method enables total DNA and protein recovery to be quantitatively assessed across different sample collection and extraction methods to better evaluate extraction efficiency. Collectively, these artificial fingerprint samples are simple to make, highly versatile and customizable, and accurately represent the biochemical composition and biological signatures of human fingerprints.

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.

Multi-locus genotype data are widely used in population genetics and disease studies. In evaluating the utility of multi-locus data, the independence of markers is commonly considered in many genomic assessments. Generally, pairwise non-random associations are tested by linkage disequilibrium; however, the dependence of one panel might be triplet, quartet, or other. Therefore, a compatible and user-friendly software is necessary for testing and assessing the global linkage disequilibrium among mixed genetic data.

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.