Search

Nonnative invasive plant species are a major cause of ecosystem degradation and impairment of ecosystem service benefits in the United States. Forested riparian areas provide many ecosystem service benefits and are vital to maintaining water quality of streams and rivers. These systems are also vulnerable to natural disturbances and invasion by nonnative plants. We assessed whether planting native trees on disturbed riparian sites may increase biotic resistance to invasive plant establishment in central Vermont in the northeastern United States. The density (stems per square meter) of invasive stems was higher in non-planted sites (x̄ = 4.1 stems m−2) compared with planted sites (x̄ = 1.3 stems m−2). More than 90% of the invasive plants were Japanese knotweed [Fallopia japonica (Houtt.) Ronse Decr.; syn. Polygonum cuspidatum Siebold & Zucc.]. There were no significant differences in total stem density of native vegetation between planted and non-planted sites. Other measured response variables such as native tree regeneration, species diversity, soil properties, and soil function showed no significant differences or trends in the paired riparian study sites. The results of this case study indicate that tree planting in disturbed riparian forest areas may assist conservation efforts by minimizing the risk of invasive plant colonization.

Whole exome sequencing (WES) of matched tumor-normal pairs in rare tumors has the potential to identify genome-wide mutations and copy number alterations (CNAs). We evaluated 27 rare cancer patients with tumor-normal matching by WES and tumor-only next generation sequencing (NGS) as a comparator. Our goal was to: 1) identify known and novel variants and CNAs in rare cancers with comparison to common cancers, 2) examine differences between germline and somatic variants and how that functionally impacts rare tumors, 3) detect and characterize alleles in biologically relevant genes-pathways that may be of clinical importance but not represented in classical cancer genes. We identified 3343 germline single nucleotide variants (SNVs) and small indel variants&mdash, 1670 in oncogenes and 1673 in tumor suppressor genes&mdash, generating an average of 124 germline variants/case. The number of somatic SNVs and small indels detected in all cases was 523:306 in oncogenes and 217 in tumor suppressor genes. Of the germline variants, six were identified to be pathogenic or likely pathogenic. In the 27 analyzed rare cancer cases, CNAs are variable depending on tumor type, germline pathogenic variants are more common. Cell fate pathway mutations (e.g., Hippo, Notch, Wnt) dominate pathogenesis and double hit (mutation + CNV) represent ~18% cases.

Variants implicated in childhood epilepsy have been identified in all four voltage-gated sodium channels that initiate action potentials in the central nervous system. Previous research has focused on the functional effects of particular variants within the most studied of these channels (NaV1.1, NaV1.2 and NaV1.6); however, there have been few comparative studies across channels to infer the impact of mutations in patients with epilepsy. Here we compare patterns of variation in patient and public databases to test the hypothesis that regions of known functional significance within voltage-gated sodium (NaV) channels have an increased burden of deleterious variants. We assessed mutational burden in different regions of the Nav channels by (1) performing Fisher exact tests on odds ratios to infer excess variants in domains, segments, and loops of each channel in patient databases versus public "control" databases, and (2) comparing the cumulative distribution of variant sites along DNA sequences of each gene in patient and public databases (i.e., independent of protein structure). Patient variant density was concordant among channels in regions known to play a role in channel function, with statistically significant higher patient variant density in S4-S6 and DIII-DIV and an excess of public variants in SI-S3, DI-DII, DII-DIII. On the other hand, channel-specific patterns of patient burden were found in the NaV1.6 inactivation gate and NaV1.1 S5-S6 linkers, while NaV1.2 and NaV1.6 S4-S5 linkers and S5 segments shared patient variant patterns that contrasted with those in NaV1.1. These different patterns may reflect different roles played by the NaV1.6 inactivation gate in action potential propagation, and by NaV1.1 S5-S6 linkers in loss of function and haploinsufficiency. Interestingly, NaV1.2 and NaV1.6 both lack amino acid substitutions over significantly long stretches in both the patient and public databases suggesting that new mutations in these regions may cause embryonic lethality or a non-epileptic disease phenotype.

// Shariful Islam 1 , Andrew L. Paek 2 , Michael Hammer 3 , Savithri Rangarajan 4 , Robert Ruijtenbeek 4 , Laurence Cooke 5 , Eric Weterings 6 and Daruka Mahadevan 5 1 Cancer Biology GIDP, University of Arizona Cancer Center, Tucson, AZ, USA 2 Department of Molecular and Cellular Biology, University of Arizona, Tucson, AZ, USA 3 Division of Biotechnology, University of Arizona Cancer Center, Tucson, AZ, USA 4 PamGene International B.V., 's-Hertogenbosch, The Netherlands 5 Department of Medicine, University of Arizona Cancer Center, Tucson, AZ, USA 6 Department of Radiation Oncology, University of Arizona, Tucson, AZ, USA Correspondence to: Daruka Mahadevan, email: dmahadevan@uacc.arizona.edu Keywords: DLBCL; aurora kinase; aneuploidy-polyploidy; TPX2; RanGAP1 Received: August 22, 2018 Accepted: October 06, 2018 Published: November 13, 2018 ABSTRACT Double-hit (DH) or double-expresser (DE) lymphomas are high-grade diffuse large B-cell lymphomas (DLBCL) that are mostly incurable with standard chemo-immunotherapy due to treatment resistance. The generation of drug-induced aneuploid/polyploid (DIAP) cells is a common effect of anti-DLBCL therapies (e.g. vincristine, doxorubicin). DIAP cells are thought to be responsible for treatment resistance, as they are capable of re-entering the cell cycle during off-therapy periods. Previously we have shown that combination of alisertib plus ibrutinib plus rituximab can partially abrogate DIAP cells and induce cell death. Here, we provide evidence that DIAP cells can re-enter the cell cycle and escape cell death during anti-DLBCL treatment. We also discuss MYC/BCL2 mediated molecular mechanism that underlie treatment resistance. We isolated aneuploid/polyploid populations of DH/DE-DLBCL cells after treatment with the aurora kinase (AK) inhibitor alisertib. Time-lapse microscopy of single polyploid cells revealed that following drug removal, a subset of these DIAP cells divide and proliferate by reductive cell divisions, including multipolar mitosis, meiosis-like nuclear fission and budding. Genomic, proteomic, and kinomic profiling demonstrated that alisertib-induced aneuploid/polyploid cells up-regulate DNA damage, DNA replication and immune evasion pathways. In addition, we identified amplified receptor tyrosine kinase and T-cell receptor signaling, as well as MYC-mediated dysregulation of the spindle assembly checkpoints RanGAP1, TPX2 and KPNA2 . We infer that these factors contribute to treatment resistance of DIAP cells. These findings provide opportunities to develop novel DH/DE-DLBCL therapies, specifically targeting DIAP cells. Key Points MYC mediated upregulation of TPX2, KPNA2 and RanGAP1 dysregulate the spindle assembly checkpoint in drug-induced polyploid cells. Drug-induced polyploid cells re-enter the cell cycle via multipolar mitosis, fission or budding, a mechanism of disease relapse.

Novel modelling was utilised in the present study to reveal significant relationships between the abundance of the Australian freshwater stream-specialist fish Galaxias olidus and metrics defining ...

Siberia is one of the coldest environments on Earth and has great seasonal temperature variation. Long-term settlement in northern Siberia undoubtedly required biological adaptation to severe cold stress, dramatic variation in photoperiod, and limited food resources. In addition, recent archeological studies show that humans first occupied Siberia at least 45,000 years ago; yet our understanding of the demographic history of modern indigenous Siberians remains incomplete. In this study, we use whole-exome sequencing data from the Nganasans and Yakuts to infer the evolutionary history of these two indigenous Siberian populations. Recognizing the complexity of the adaptive process, we designed a model-based test to systematically search for signatures of polygenic selection. Our approach accounts for stochasticity in the demographic process and the hitchhiking effect of classic selective sweeps, as well as potential biases resulting from recombination rate and mutation rate heterogeneity. Our demographic inference shows that the Nganasans and Yakuts diverged ∼12,000-13,000 years ago from East-Asian ancestors in a process involving continuous gene flow. Our polygenic selection scan identifies seven candidate gene sets with Siberian-specific signals. Three of these gene sets are related to diet, especially to fat metabolism, consistent with the hypothesis of adaptation to a fat-rich animal diet. Additional testing rejects the effect of hitchhiking and favors a model in which selection yields small allele frequency changes at multiple unlinked genes.

The Ashkenazim are a well-studied population with a dynamic and unique history. The spread of Jewish populations from the Middle East to locations in the Diaspora in the past 2000 years represents a microcosm of the evolutionary processes affecting modern humans as they dispersed out of Africa (e.g. founder effects, migrations, admixture and adaptation). Ashkenazim – Jews who trace their ancestry to founders who moved to the Rhine Valley in the fourth century – are no exception. Recent studies both support and refute long-standing hypotheses of Ashkenazi origins. Together, the results Strengthen the inferences that Ashkenazi Jews trace their ancestry to a genetically diverse population in the Middle East that underwent a series of founder events – reducing genetic diversity relative to the ancestral population – and that subsequently experienced significant gene flow during the period of Ashkenazi residency in Europe. Key Concepts Ashkenazi Jews are genetically more similar to each other than to individuals from other populations. Ashkenazi Jews are most closely related to other Jewish populations. The Ashkenazim likely arose from a genetically diverse population in the Middle East. Ashkenazi Jews appear equally close to both Middle Eastern and European populations. Ashkenazi Jews are not closely related to modern populations that best represent the Khazars. Ashkenazi Jews likely received 30–60% gene flow from Europeans. Ashkenazi Jews likely underwent multiple extreme events of population size reduction and expansion. About half of the genomic regions under putative selection in the Ashkenazim are not found to be under selection in Europeans. Most Ashkenazi disease loci are not under strong positive selection. Keywords: Ashkenazim; population genetics; genomics; demography; SNP array; admixture; founder effect; selection

Comparisons of whole-genome sequences from ancient and contemporary samples have pointed to several instances of archaic admixture through interbreeding between the ancestors of modern non-Africans and now extinct hominids such as Neanderthals and Denisovans. One implication of these findings is that some adaptive features in contemporary humans may have entered the population via gene flow with archaic forms in Eurasia. Within Africa, fossil evidence suggests that anatomically modern humans (AMH) and various archaic forms coexisted for much of the last 200,000 yr; however, the absence of ancient DNA in Africa has limited our ability to make a direct comparison between archaic and modern human genomes. Here, we use statistical inference based on high coverage whole-genome data (greater than 60×) from contemporary African Pygmy hunter-gatherers as an alternative means to study the evolutionary history of the genus Homo. Using whole-genome simulations that consider demographic histories that include both isolation and gene flow with neighboring farming populations, our inference method rejects the hypothesis that the ancestors of AMH were genetically isolated in Africa, thus providing model-based whole genome-level evidence of African archaic admixture. Our inferences also suggest a complex human evolutionary history in Africa, which involves at least a single admixture event from an unknown archaic population into the ancestors of AMH, likely within the last 30,000 yr.

OBJECTIVES We examined autosomal genome-wide SNPs and Y-chromosome data from 15 Siberian and 12 reference populations to study the affinities of Siberian populations, and to address hypotheses about the origin of the Samoyed peoples. METHODS Samples were genotyped for 567 096 autosomal SNPs and 147 Y-chromosome polymorphic sites. For several analyses, we used 281 093 SNPs from the intersection of our data with publicly available ancient Siberian samples. To examine genetic relatedness among populations, we applied PCA, FST , TreeMix, and ADMIXTURE analyses. To explore the potential effect of demography and evolutionary processes, the distribution of ROH and IBD sharing within population were studied. RESULTS Analyses of autosomal and Y-chromosome data reveal high differentiation of the Siberian groups. The Siberian populations have a large proportion of their genome in ROH and IBD segments. Several populations (ie, Nganasans, Evenks, Yukagirs, and Koryaks) do not appear to have experienced admixture with other Siberian populations (ie, producing only positive f3), while for the other tested populations the composition of mixing sources always included Nganasans or Evenks. The Nganasans from the Taymyr Peninsula demonstrate the greatest level of shared shorter ROH and IBD with nearly all other Siberian populations. CONCLUSIONS Autosomal SNP and Y-chromosome data demonstrate that Samoyedic populations differ significantly in their genetic composition. Genetic relationship is observed only between Forest and Tundra Nentsi. Selkups are affiliated with the Kets from the Yenisey River, while the Nganasans are separated from their linguistic neighbors, showing closer affinities with the Evenks and Yukagirs.

Evolutionary genetic studies have shown a positive correlation between levels of nucleotide diversity and either rates of recombination or genetic distance to genes. Both positive-directional and purifying selection have been offered as the source of these correlations via genetic hitchhiking and background selection, respectively. Phylogenetically conserved elements (CEs) are short (∼100 bp), widely distributed (comprising ∼5% of genome), sequences that are often found far from genes. While the function of many CEs is unknown, CEs also are associated with reduced diversity at linked sites. Using high coverage (>80×) whole genome data from two human populations, the Yoruba and the CEU, we perform fine scale evaluations of diversity, rates of recombination, and linkage to genes. We find that the local rate of recombination has a stronger effect on levels of diversity than linkage to genes, and that these effects of recombination persist even in regions far from genes. Our whole genome modeling demonstrates that, rather than recombination or GC-biased gene conversion, selection on sites within or linked to CEs better explains the observed genomic diversity patterns. A major implication is that very few sites in the human genome are predicted to be free of the effects of selection. These sites, which we refer to as the human "neutralome," comprise only 1.2% of the autosomes and 5.1% of the X chromosome. Demographic analysis of the neutralome reveals larger population sizes and lower rates of growth for ancestral human populations than inferred by previous analyses.

We investigated the genome of a 5-year-old male who presented with global developmental delay (motor, cognitive, and speech), hypotonia, possibly ataxia, and cerebellar hypoplasia of unknown origin. Whole genome sequencing (WGS) and mRNA sequencing (RNA-seq) were performed on a family having an affected proband, his unaffected parents, and maternal grandfather. To explore the molecular and functional consequences of the variant, we performed cell proliferation assays, quantitative real-time PCR (qRT-PCR) array, immunoblotting, calcium imaging, and neurite outgrowth experiments in SH-SY5Y neuroblastoma cells to compare the properties of the wild-type TATA-box-binding protein factor 1 (TAF1), deletion of TAF1, and TAF1 variant p.Ser1600Gly samples. The whole genome data identified several gene variants. However, the genome sequence data failed to implicate a candidate gene as many of the variants were of unknown significance. By combining genome sequence data with transcriptomic data, a probable candidate variant, p.Ser1600Gly, emerged in TAF1. Moreover, the RNA-seq data revealed a 90:10 extremely skewed X-chromosome inactivation (XCI) in the mother. Our results showed that neuronal ion channel genes were differentially expressed between TAF1 deletion and TAF1 variant p.Ser1600Gly cells, when compared with their respective controls, and that the TAF1 variant may impair neuronal differentiation and cell proliferation. Taken together, our data suggest that this novel variant in TAF1 plays a key role in the development of a recently described X-linked syndrome, TAF1 intellectual disability syndrome, and further extends our knowledge of a potential link between TAF1 deficiency and defects in neuronal cell function.

As a result of the combination of great linguistic and cultural diversity, the highland populations of Daghestan present an excellent opportunity to test the hypothesis of language-gene coevolution at a fine geographic scale. However, previous genetic studies generally have been restricted to uniparental markers and have not included many of the key populations of the region. To improve our understanding of the genetic structure of Daghestani populations and to investigate possible correlations between genetic and linguistic variation, we analyzed ~550,000 autosomal single nucleotide polymorphisms, phylogenetically informative Y chromosome markers and mtDNA haplotypes in 21 ethnic Daghestani groups. We found high levels of population structure in Daghestan consistent with the hypothesis of long-term isolation among populations of the highland Caucasus. Highland Daghestani populations exhibit extremely high levels of between-population diversity for all genetic systems tested, leading to some of the highest FST values observed for any region of the world. In addition, we find a significant positive correlation between gene and language diversity, suggesting that these two aspects of human diversity have coevolved as a result of historical patterns of social interaction among highland farmers at the community level. Finally, our data are consistent with the hypothesis that most Daghestanian-speaking groups descend from a common ancestral population (~6000-6500 years ago) that spread to the Caucasus by demic diffusion followed by population fragmentation and low levels of gene flow.

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.

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.

Summary Objective 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. Design Whole exome sequencing was used to discover the variant. We generated a mutant cDNA, transfected HEK293 cells, and performed Western blotting to assess protein stability. To study channel functional properties, patch-clamp experiments were carried out in transfected neuronal ND7/23 cells. Results The proband exhibited seizure onset at 6 months of age, diffuse brain atrophy, and more profound developmental impairment than the original case. The mutation p.Arg233Gly in the voltage sensing transmembrane segment D1S4 was present in the proband and absent in both parents. This mutation results in a temperature-sensitive reduction in protein expression as well as reduced sodium current amplitude and density and a relative increased response to a slow ramp stimulus, though this did not result in an absolute increased current at physiological temperatures. Conclusion The new de novo SCN8A mutation is clearly deleterious, resulting in an unstable protein with reduced channel activity. This differs from the gain-of-function attributes of the first SCN8A mutation in epileptic encephalopathy, pointing to heterogeneity of mechanisms. Since Nav1.6 is expressed in both excitatory and inhibitory neurons, a differential effect of a loss-of-function of Nav1.6 Arg223Gly on inhibitory interneurons may underlie the epilepsy phenotype in this patient.

Recently, Elhaik et al1 criticized several aspects of an analysis in which Mendez et al2 estimate the time to the most recent common ancestor (TMRCA) for the Y chromosome tree incorporating a newly identified basal branch called A00. Elhaik et al asserted that Mendez et al derived an inflated estimate of the TMRCA by applying incorrect assumptions and approximations, numerical miscalculations, and data manipulation. In particular they focused on (1) the method used to estimate the Y chromosome mutation rate, (2) the relative lengths of sequences that should be compared in order to estimate branch lengths in a tree, and (3) the implications of the Y chromosome TMRCA estimates reported by Mendez et al for human evolution. Here we show that these criticisms result from a misunderstanding of population genetic theory, as well as a misrepresentation of the methodology of Mendez et al. However, before addressing the various arguments put forward by Elhaik et al, we comment on conceptual and theoretical issues surrounding the significance of the Y chromosome TMRCA.

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.

Partially recessive variants under positive selection are expected to go to fixation more quickly on the X chromosome as a result of hemizygosity, an effect known as faster-X. Conversely, purifying selection is expected to reduce substitution rates more effectively on the X chromosome. Previous work in humans contrasted divergence on the autosomes and X chromosome, with results tending to support the faster-X effect. However, no study has yet incorporated both divergence and polymorphism to quantify the effects of both purifying and positive selection, which are opposing forces with respect to divergence. In this study, we develop a framework that integrates previously developed theory addressing differential rates of X and autosomal evolution with methods that jointly estimate the level of purifying and positive selection via modeling of the distribution of fitness effects (DFE). We then utilize this framework to estimate the proportion of nonsynonymous substitutions fixed by positive selection (α) using exome sequence data from a West African population. We find that varying the female to male breeding ratio (β) has minimal impact on the DFE for the X chromosome, especially when compared with the effect of varying the dominance coefficient of deleterious alleles (h). Estimates of α range from 46% to 51% and from 4% to 24% for the X chromosome and autosomes, respectively. While dependent on h, the magnitude of the difference between α values estimated for these two systems is highly statistically significant over a range of biologically realistic parameter values, suggesting faster-X has been operating in humans.

Examining patterns of molecular genetic variation in both modern-day and ancient humans has proved to be a powerful approach to learn about our origins. Rapid advances in DNA sequencing technology have allowed us to characterize increasing amounts of genomic information. Although this clearly provides unprecedented power for inference, it also introduces more complexity into the way we use and interpret such data. Here, we review ongoing debates that have been influenced by improvements in our ability to sequence DNA and discuss some of the analytical challenges that need to be overcome in order to fully exploit the rich historical information that is contained in the entirety of the human genome.

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.

Analyses of ancient DNA from extinct humans reveal signals of at least two independent hybridization events in the history of non-African populations. To date, there are very few examples of specific genetic variants that have been rigorously identified as introgressive. Here, we survey DNA sequence variation in the OAS gene cluster on chromosome 12 and provide strong evidence that a haplotype extending for ~185 kb introgressed from Neandertals. This haplotype is nearly restricted to Eurasians and is estimated to have diverged from the Neandertal sequence ~125 kya. Despite the potential for novel functional variation, the observed frequency of this haplotype is consistent with neutral introgression. This is the second locus in the human genome, after STAT2, carrying distinct haplotypes that appear to have introgressed separately from both Neandertals and Denisova.

SummaryObjective Two major classes of SCN1A variants are associated with Dravet syndrome (DS): those that result in haploinsufficiency (truncating) and those that result in an amino acid substitution (missense). The aim of this retrospective study was to describe the first large cohort of Japanese patients with SCN1A mutation–positive DS (n = 285), and investigate the relationship between variant (type and position) and clinical expression and response to treatment. Methods We sequenced all exons and intron–exon boundaries of SCN1A in our cohort, investigated differences in the distribution of truncating and missense variants, tested for associations between variant type and phenotype, and compared these patterns with those of cohorts with milder epilepsy and healthy individuals. Results Unlike truncation variants, missense variants are found at higher density in the S4 voltage sensor and pore loops and at lower density in the domain I–II and II–III linkers and the first three segments of domain II. Relative to healthy individuals, there is an increased frequency of truncating (but not missense) variants in the noncoding C-terminus. The rate of cognitive decline is more rapid for patients with truncation variants regardless of age at seizure onset, whereas age at onset is a predictor of the rate of cognitive decline for patients with missense variants. Significance We found significant differences in the distribution of truncating and missense variants across the SCN1A sequence among healthy individuals, patients with DS, and those with milder forms of SCN1A-variant positive epilepsy. Testing for associations with phenotype revealed that variant type can be predictive of rate of cognitive decline. Analysis of descriptive medication data suggests that in addition to conventional drug therapy in DS, bromide, clonazepam and topiramate may reduce seizure frequency.

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.

Here we report the Simons Genome Diversity Project data set: high quality genomes from 300 individuals from 142 diverse populations. These genomes include at least 5.8 million base pairs that are not present in the human reference genome. Our analysis reveals key features of the landscape of human genome variation, including that the rate of accumulation of mutations has accelerated by about 5% in non-Africans compared to Africans since divergence. We show that the ancestors of some pairs of present-day human populations were substantially separated by 100,000 years ago, well before the archaeologically attested onset of behavioural modernity. We also demonstrate that indigenous Australians, New Guineans and Andamanese do not derive substantial ancestry from an early dispersal of modern humans; instead, their modern human ancestry is consistent with coming from the same source as that of other non-Africans., Facultad de Ciencias Naturales y Museo

Although some mitochondrial, X chromosome, and autosomal sequence diversity data are available for our closest relatives,Pan troglodytesandPan paniscus, data from the nonrecombining portion of the Y chromosome (NRY) are more limited. We examined ≈3 kb of NRY DNA from 101 chimpanzees, seven bonobos, and 42 humans to investigate: (i) relative levels of intraspecific diversity; (ii) the degree of paternal lineage sorting among species and subspecies of the genusPan; and (iii) the date of the chimpanzee/bonobo divergence. We identified 10 informative sequence-tagged sites associated with 23 polymorphisms on the NRY from the genusPan. Nucleotide diversity was significantly higher on the NRY of chimpanzees and bonobos than on the human NRY. Similar to mtDNA, but unlike X-linked and autosomal loci, lineages defined by mutations on the NRY were not shared among subspecies ofP. troglodytes. Comparisons with mtDNAND2sequences from some of the same individuals revealed a larger female versus male effective population size for chimpanzees. The NRY-based divergence time between chimpanzees and bonobos was estimated at ≈1.8 million years ago. In contrast to human populations who appear to have had a low effective size and a recent origin with subsequent population growth, some taxa within the genusPanmay be characterized by large populations of relatively constant size, more ancient origins, and high levels of subdivision.

Prado-Martinez, Javier et al., Most great ape genetic variation remains uncharacterized1,2; however, its study is critical for understanding population history3–6, recombination7, selection8 and susceptibility to disease9,10.Herewe 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.Wefind extensive inbreeding in almost all wild populations, with eastern gorillas being the most extreme. Inferred effective population sizes have varied radically over timein 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 genomediversity provides a framework for understanding evolution and a resource for more effective management of wild and captive great ape populations., We thank the following funding agencies: ERC Starting Grant (260372) to T.M.-B.; NIH grants HG002385 to E.E.E., R01_HG005226 to K.R.V., A.E.W., M.F.H., L.S. and J.D.W., GM100233 and NSF HOMINID grant 1032255 to D.R. and He.Li.;MICINN(Spain)BFU2011-28549to T.M.-B.,BFU2010-19443to Ja.Be.,Spanish Government and FEDER for grants BFU2009-13409-C02-02 and BFU2012-38236 to A.N. and J.P.-M., Direccio´ General de Recerca, Generalitat de Catalunya (Grup de Recerca Consolidat 2009 SGR 1101) to Ja.Be., D.C., A.N. and T.M.-B.; ERC Advanced Grant (233297) and Max Planck Society to S. Paabo; Danish Council for Independent Research Natural Sciences to H.S.; Spanish Grant (CGL-2010-20170) and Zoo de Barcelona (Beca PRIC) to A.R.-H.; EUPRIM-Net to BPRC; DP1ES022577-04 NIH grant to S.A.T.; NSF Grant 0755823 to M.K.G.; P.G. is supported by the G. Harold and Leila Y. Mathers Foundation. A.N. and T.M.-B. are ICREA Research Investigators (Institut Catala d’Estudis i RecercaAvancats delaGeneralitat deCatalunya). J.P.-M. is supported by the Zoo de Barcelona and l’Ajuntament de Barcelona.

Comparisons of whole-genome sequences from ancient and contemporary samples have pointed to several instances of archaic admixture through interbreeding between the ancestors of modern non-Africans and now extinct hominids such as Neanderthals and Denisovans. One implication of these findings is that some adaptive features in contemporary humans may have entered the population via gene flow with archaic forms in Eurasia. Within Africa, fossil evidence suggests that anatomically modern humans (AMH) and various archaic forms coexisted for much of the last 200,000 yr; however, the absence of ancient DNA in Africa has limited our ability to make a direct comparison between archaic and modern human genomes. Here, we use statistical inference based on high coverage whole-genome data (greater than 60×) from contemporary African Pygmy hunter-gatherers as an alternative means to study the evolutionary history of the genus Homo. Using whole-genome simulations that consider demographic histories that include both isolation and gene flow with neighboring farming populations, our inference method rejects the hypothesis that the ancestors of AMH were genetically isolated in Africa, thus providing the first whole genome-level evidence of African archaic admixture. Our inferences also suggest a complex human evolutionary history in Africa, which involves at least a single admixture event from an unknown archaic population into the ancestors of AMH, likely within the last 30,000 yr.

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.

Recent analysis of DNA extracted from two Eurasian forms of archaic human shows that more genetic variants are shared with humans currently living in Eurasia than with anatomically modern humans in sub-Saharan Africa. Although these genome-wide average measures of genetic similarity are consistent with the hypothesis of archaic admixture in Eurasia, analyses of individual loci exhibiting the signal of archaic introgression are needed to test alternative hypotheses and investigate the admixture process. Here, we provide a detailed sequence analysis of the innate immune gene OAS1, a locus with a divergent Melanesian haplotype that is very similar to the Denisova sequence from the Altai region of Siberia. We resequenced a 7-kb region encompassing the OAS1 gene in 88 individuals from six Old World populations (San, Biaka, Mandenka, French Basque, Han Chinese, and Papua New Guineans) and discovered previously unknown and ancient genetic variation. The 5' region of this gene has unusual patterns of diversity, including 1) higher levels of nucleotide diversity in Papuans than in sub-Saharan Africans, 2) very deep ancestry with an estimated time to the most recent common ancestor of >3 myr, and 3) a basal branching pattern with Papuan individuals on either side of the rooted network. A global geographic survey of >1,500 individuals showed that the divergent Papuan haplotype is nearly restricted to populations from eastern Indonesia and Melanesia. Polymorphic sites within this haplotype are shared with the draft Denisova genome over a span of ∼90 kb and are associated with an extended block of linkage disequilibrium, supporting the hypothesis that this haplotype introgressed from an archaic source that likely lived in Eurasia.

A long-debated question concerns the fate of archaic forms of the genus Homo : did they go extinct without interbreeding with anatomically modern humans, or are their genes present in contemporary populations? This question is typically focused on the genetic contribution of archaic forms outside of Africa. Here we use DNA sequence data gathered from 61 noncoding autosomal regions in a sample of three sub-Saharan African populations (Mandenka, Biaka, and San) to test models of African archaic admixture. We use two complementary approximate-likelihood approaches and a model of human evolution that involves recent population structure, with and without gene flow from an archaic population. Extensive simulation results reject the null model of no admixture and allow us to infer that contemporary African populations contain a small proportion of genetic material (≈2%) that introgressed ≈35 kya from an archaic population that split from the ancestors of anatomically modern humans ≈700 kya. Three candidate regions showing deep haplotype divergence, unusual patterns of linkage disequilibrium, and small basal clade size are identified and the distributions of introgressive haplotypes surveyed in a sample of populations from across sub-Saharan Africa. One candidate locus with an unusual segment of DNA that extends for >31 kb on chromosome 4 seems to have introgressed into modern Africans from a now-extinct taxon that may have lived in central Africa. Taken together our results suggest that polymorphisms present in extant populations introgressed via relatively recent interbreeding with hominin forms that diverged from the ancestors of modern humans in the Lower-Middle Pleistocene.

Modern day Latin America resulted from the encounter of Europeans with the indigenous peoples of the Americas in 1492, followed by waves of migration from Europe and Africa. As a result, the genomic structure of present day Latin Americans was determined both by the genetic structure of the founding populations and the numbers of migrants from these different populations. Here, we analyzed DNA collected from two well-established communities in Colorado (33 unrelated individuals) and Ecuador (20 unrelated individuals) with a measurable prevalence of the BRCA1 c.185delAG and the GHR c.E180 mutations, respectively, using Affymetrix Genome-wide Human SNP 6.0 arrays to identify their ancestry. These mutations are thought to have been brought to these communities by Sephardic Jewish progenitors. Principal component analysis and clustering methods were employed to determine the genome-wide patterns of continental ancestry within both populations using single nucleotide polymorphisms, complemented by determination of Y-chromosomal and mitochondrial DNA haplotypes. When examining the presumed European component of these two communities, we demonstrate enrichment for Sephardic Jewish ancestry not only for these mutations, but also for other segments as well. Although comparison of both groups to a reference Hispanic/Latino population of Mexicans demonstrated proximity and similarity to other modern day communities derived from a European and Native American two-way admixture, identity-by-descent and Y-chromosome mapping demonstrated signatures of Sephardim in both communities. These findings are consistent with historical accounts of Jewish migration from the realms that comprise modern Spain and Portugal during the Age of Discovery. More importantly, they provide a rationale for the occurrence of mutations typically associated with the Jewish Diaspora in Latin American communities.

Bacteria play an important role in the pathogenesis of inflammatory bowel disease (IBD), its complications and its symptoms. Antibiotics can decrease tissue invasion and eliminate aggressive bacterial species. They are used in IBD to treat infective complications and for altering bacterial flora, which may result in specific anti-inflammatory effects. In addition, suppression of bacterial metabolic activities or direct effects of antibiotics on intestinal structures and functions may result in symptoms which cannot be differentiated from symptoms caused by inflammation. Although current clinical trials do not fulfill criteria of evidence-based treatment, a few placebo- or standard treatment-controlled studies suggest that metronidazole and ciprofloxacin are effective in Crohn’s colitis and ileocolitis, perianal fistulae and pouchitis. Administration of probiotics, prebiotics and synbiotics can restore a predominance of beneficial species. However, beneficial effects of probiotics in IBD are modest, strain-specific and limited to certain manifestations of disease and duration of use of the probiotic. For probiotics there is reasonable evidence of efficacy in relapse prevention in chronic pouchitis and ulcerative colitis, and suggestive evidence for postoperative prevention in pouchitis. Therapeutic manipulation of the intestinal flora offers considerable promise for treating IBD, but must be supported by large controlled therapeutic trials before widespread clinical acceptance. These agents may become a component of treating IBD in combination with traditional anti-inflammatory and immunosuppressive agents. Probiotic strategies, based on metagenomic or metabonomic analyses, and new classes of probiotics might play an important role in the future management of IBD.

Key Points MYC mediated upregulation of TPX2, KPNA2 and RanGAP1 dysregulate the spindle assembly checkpoint in drug-induced polyploid cells. Drug-induced polyploid cells re-enter cell cycle via multipolar mitosis, fission or budding, a mechanism of disease relapse. Abstract Double-hit (DH) or double-expresser (DE) diffuse large B-cell lymphomas (DLBCL) are high-grade B-cell lymphomas that are mostly incurable with standard chemo-immunotherapy due to treatment resistance. The generation of drug-induced aneuploid/polyploid (DIAP) cells is a common effect of anti-DLBCL therapies (e.g. vincristine, doxorubicin). DIAP cells are thought to be responsible for treatment resistance, as they are capable of re-entering the cell cycle during off-therapy periods. Previously (DOI: 10.1158/1535-7163.MCT-17-0089) we demonstrated that combination of alisertib plus ibrutinib plus rituximab can partially abrogate DIAP cells and induce cell death. Here, we provide evidence that DIAP cells can re-enter the cell cycle and escape apoptosis during anti-DLBCL treatment. We also discuss MYC/BCL2 mediated molecular mechanisms that underlie treatment resistance. We isolated aneuploid/polyploid populations (2n-similar, 4n and 8n) of DH/DE-DLBCL cells after treatment with the aurora kinase (AK) inhibitor alisertib. Time-lapse microscopy of single polyploid cells (8n) revealed that following drug removal, these cells divide and proliferate by reductive cell divisions, including multipolar mitosis, meiosis-like nuclear fission and budding. Genomic, proteomic, and kinomic profiling of DIAP cells demonstrated up-regulation DNA damage, DNA replication and immune evasion pathways. In addition, the immune-kinome, AKT/MTorc1/2, and ERK/MAPK pathway are over-represented in 8n cells. Further, MYC-mediates dysregulation of the spindle assembly checkpoints by over-expression of RanGAP1, TPX2 and KPNA2. Multiple mechanisms contribute to DIAP which are amenable to novel therapeutic intervention(s) in MYC/BCL2-addicted high-grade B-cell lymphomas. Disclosures No relevant conflicts of interest to declare.

Hispanic/Latino populations possess a complex genetic structure that reflects recent admixture among and potentially ancient substructure within Native American, European, and West African source populations. Here, we quantify genome-wide patterns of SNP and haplotype variation among 100 individuals with ancestry from Ecuador, Colombia, Puerto Rico, and the Dominican Republic genotyped on the Illumina 610-Quad arrays and 112 Mexicans genotyped on Affymetrix 500K platform. Intersecting these data with previously collected high-density SNP data from 4,305 individuals, we use principal component analysis and clustering methods FRAPPE and STRUCTURE to investigate genome-wide patterns of African, European, and Native American population structure within and among Hispanic/Latino populations. Comparing autosomal, X and Y chromosome, and mtDNA variation, we find evidence of a significant sex bias in admixture proportions consistent with disproportionate contribution of European male and Native American female ancestry to present-day populations. We also find that patterns of linkage-disequilibria in admixed Hispanic/Latino populations are largely affected by the admixture dynamics of the populations, with faster decay of LD in populations of higher African ancestry. Finally, using the locus-specific ancestry inference method LAMP , we reconstruct fine-scale chromosomal patterns of admixture. We document moderate power to differentiate among potential subcontinental source populations within the Native American, European, and African segments of the admixed Hispanic/Latino genomes. Our results suggest future genome-wide association scans in Hispanic/Latino populations may require correction for local genomic ancestry at a subcontinental scale when associating differences in the genome with disease risk, progression, and drug efficacy, as well as for admixture mapping.

The early history of island Southeast Asia is often characterized as the story of two major population dispersals: the initial Paleolithic colonization of Sahul approximately 45 ka ago and the much later Neolithic expansion of Austronesian-speaking farmers approximately 4 ka ago. Here, in the largest survey of Indonesian Y chromosomes to date, we present evidence for multiple genetic strata that likely arose through a series of distinct migratory processes. We genotype an extensive battery of Y chromosome markers, including 85 single-nucleotide polymorphisms/indels and 12 short tandem repeats, in a sample of 1,917 men from 32 communities located across Indonesia. We find that the paternal gene pool is sharply subdivided between western and eastern locations, with a boundary running between the islands of Bali and Flores. Analysis of molecular variance reveals one of the highest levels of between-group variance yet reported for human Y chromosome data (e.g., Phi(ST) = 0.47). Eastern Y chromosome haplogroups are closely related to Melanesian lineages (i.e., within the C, M, and S subclades) and likely reflect the initial wave of colonization of the region, whereas the majority of western Y chromosomes (i.e., O-M119*, O-P203, and O-M95*) are related to haplogroups that may have entered Indonesia during the Paleolithic from mainland Asia. In addition, two novel markers (P201 and P203) provide significantly enhanced phylogenetic resolution of two key haplogroups (O-M122 and O-M119) that are often associated with the Austronesian expansion. This more refined picture leads us to put forward a four-phase colonization model in which Paleolithic migrations of hunter-gatherers shape the primary structure of current Indonesian Y chromosome diversity, and Neolithic incursions make only a minor impact on the paternal gene pool, despite the large cultural impact of the Austronesian expansion.

The geographical region between mainland Asia and New Guinea is characterized by numerous small islands with isolated human populations. Phenotypically, groups in the west are similar to their neighbours in mainland Southeast Asia, eastern groups near New Guinea are similar to Melanesians, and intervening populations are intermediate in appearance. A long-standing question is whether this pattern primarily reflects mixing between groups with distinct origins or whether natural selection has shaped this range of variation by acting differentially on populations across the region. To address this question, we genotyped a set of 37 single nucleotide polymorphisms that are evolutionarily independent, putatively neutral and highly informative for Asian–Melanesian ancestry in 1430 individuals from 60 populations spanning mainland Asia to Melanesia. Admixture analysis reveals a sharp transition from Asian to Melanesian genetic variants over a narrow geographical region in eastern Indonesia. Interestingly, this admixture cline roughly corresponds to the human phenotypic boundary noted by Alfred Russell Wallace in 1869. We conclude that this phenotypic gradient probably reflects mixing of two long-separated ancestral source populations—one descended from the initial Melanesian-like inhabitants of the region, and the other related to Asian groups that immigrated during the Paleolithic and/or with the spread of agriculture. A higher frequency of Asian X-linked markers relative to autosomal markers throughout the transition zone suggests that the admixture process was sex-biased, either favouring a westward expansion of patrilocal Melanesian groups or an eastward expansion of matrilocal Asian immigrants. The matrilocal marriage practices that dominated early Austronesian societies may be one factor contributing to this observed sex bias in admixture rates.

It has been known for over a decade that a majority of men who self report as members of the Jewish priesthood (Cohanim) carry a characteristic Y chromosome haplotype termed the Cohen Modal Haplotype (CMH). The CMH has since been used to trace putative Jewish ancestral origins of various populations. However, the limited number of binary and STR Y chromosome markers used previously did not provide the phylogenetic resolution needed to infer the number of independent paternal lineages that are encompassed within the Cohanim or their coalescence times. Accordingly, we have genotyped 75 binary markers and 12 Y-STRs in a sample of 215 Cohanim from diverse Jewish communities, 1,575 Jewish men from across the range of the Jewish Diaspora, and 2,099 non-Jewish men from the Near East, Europe, Central Asia, and India. While Cohanim from diverse backgrounds carry a total of 21 Y chromosome haplogroups, 5 haplogroups account for 79.5% of Cohanim Y chromosomes. The most frequent Cohanim lineage (46.1%) is marked by the recently reported P58 T->C mutation, which is prevalent in the Near East. Based on genotypes at 12 Y-STRs, we identify an extended CMH on the J-P58* background that predominates in both Ashkenazi and non-Ashkenazi Cohanim and is remarkably absent in non-Jews. The estimated divergence time of this lineage based on 17 STRs is 3,190 ± 1,090 years. Notably, the second most frequent Cohanim lineage (J-M410*, 14.4%) contains an extended modal haplotype that is also limited to Ashkenazi and non-Ashkenazi Cohanim and is estimated to be 4.2 ± 1.3 ky old. These results support the hypothesis of a common origin of the CMH in the Near East well before the dispersion of the Jewish people into separate communities, and indicate that the majority of contemporary Jewish priests descend from a limited number of paternal lineages. Electronic supplementary material The online version of this article (doi:10.1007/s00439-009-0727-5) contains supplementary material, which is available to authorized users.

A central tenet of evolutionary social science holds that behaviors, such as those associated with social dominance, produce fitness effects that are subject to cultural selection. However, evidence for such selection is inconclusive because it is based on short-term statistical associations between behavior and fertility. Here, we show that the evolutionary effects of dominance at the population level can be detected using noncoding regions of DNA. Highly variable polymorphisms on the nonrecombining portion of the Y chromosome can be used to trace lines of descent from a common male ancestor. Thus, it is possible to test for the persistence of differential fertility among patrilines. We examine haplotype distributions defined by 12 short tandem repeats in a sample of 1269 men from 41 Indonesian communities and test for departures from neutral mutation-drift equilibrium based on the Ewens sampling formula. Our tests reject the neutral model in only 5 communities. Analysis and simulations show that we have sufficient power to detect such departures under varying demographic conditions, including founder effects, bottlenecks, and migration, and at varying levels of social dominance. We conclude that patrilines seldom are dominant for more than a few generations, and thus traits or behaviors that are strictly paternally inherited are unlikely to be under strong cultural selection.

Markers on the non-recombining portion of the human Y chromosome continue to have applications in many fields including evolutionary biology, forensics, medical genetics, and genealogical reconstruction. In 2002, the Y Chromosome Consortium published a single parsimony tree showing the relationships among 153 haplogroups based on 243 binary markers and devised a standardized nomenclature system to name lineages nested within this tree. Here we present an extensively revised Y chromosome tree containing 311 distinct haplogroups, including two new major haplogroups (S and T), and incorporating approximately 600 binary markers. We describe major changes in the topology of the parsimony tree and provide names for new and rearranged lineages within the tree following the rules presented by the Y Chromosome Consortium in 2002. Several changes in the tree topology have important implications for studies of human ancestry. We also present demography-independent age estimates for 11 of the major clades in the new Y chromosome tree.

Sudmant, Peter H. et al., In order to explore the diversity and selective signatures of duplication and deletion human copy-number variants (CNVs), we sequenced 236 individuals from 125 distinct human populations. We observed that duplications exhibit fundamentally different population genetic and selective signatures than deletions and are more likely to be stratified between human populations. Through reconstruction of the ancestral human genome, we identify megabases of DNA lost in different human lineages and pinpoint large duplications that introgressed from the extinct Denisova lineage now found at high frequency exclusively in Oceanic populations. We find that the proportion of CNV base pairs to single-nucleotide–variant base pairs is greater among non-Africans than it is among African populations, but we conclude that this difference is likely due to unique aspects of non-African population history as opposed to differences in CNV load., This project has been funded in part with federal funds from the National Cancer Institute, NIH, under contract HHSN26120080001E. The content of this publication does not necessarily reflect the views or policies of the Department of Health and Human Services, nor does the mention of trade names, commercial products, or organizations imply endorsement by the U.S. government. This research was supported in part by the Intramural Research Program of the NIH, National Cancer Institute, Center for Cancer Research. This work was also partly supported by NIH grant 2R01HG002385 and a grant (11631) from the Paul G. Allen Family Foundation to E.E.E. The sequencing for this study was supported by a grant from the Simons Foundation to D.R. (SFARI 280376) and by a HOMINID grant from the NSF to D.R. (BCS-1032255). T.K. is supported by a European Research Council Starting Investigator grant (FP7 - 26213). R.S. and S.D. received support from the Ministry of Education and Science, Russian Federation (14.Z50.31.0010). H.S., E.M., R.V., and M.M. are supported by Institutional Research Funding from the Estonian Research Council IUT24-1 and by the European Regional Development Fund (European Union) through the Centre of Excellence in Genomics to Estonian Biocentre and University of Tartu. S.A.T. is supported by NIH grants 5DP1ES022577 05, 1R01DK104339-01, and 1R01GM113657-01. C.T.-S. is supported by Wellcome Trust grant 098051. C.M.B. is supported by the NSF (award numbers 0924726 and 1153911). E.E.E. and D.R. are investigators of the Howard Hughes Medical Institute. Data are deposited into ENA (PRJEB9586 or ERP010710), and variant calls are deposited in dbVar (PRJNA285786). E.E.E. is on the scientific advisory board of DNAnexus, Incorporated, and is a consultant for Kunming University of Science and Technology (KUST) as part of the 1000 China Talent Program.

It is commonly thought that human genetic diversity in non-African populations was primarily shaped by a recent out-of-Africa dispersal. Here we present a large geographical Y chromosome study using 459 high coverage sequences, including 302 newly reported here. We date the Y chromosomal Most Recent Common Ancestor (MRCA) in Africa at 254 (95% CI 192- 307) kya and differentiation of African and non-African lineages 52-121 kya. The age estimates for major non-African founder haplogroups cluster closely 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. We find that the highest basal Y chromosome diversity outside Africa has been preserved in South and Southeast Asians while extant Y chromosome diversity in Europe and the Near East stems from a small number of mid-Holocene founders. In contrast to demographic reconstructions based on mtDNA, we infer a second strong bottleneck, followed by a fast recovery in Old World populations dating to the last 10 ky. We hypothesize that this bottleneck is caused by cultural changes affecting variance of reproductive success among males.

To investigate the paternal population history of populations in Northern Island Melanesia, 685 paternally unrelated males from 36 populations in this region and New Guinea were analyzed at 14 regionally informative binary markers and 7 short tandem repeat (STR) loci from the nonrecombining portion of the Y chromosome. Three newly defined binary markers (K6-P79, K7-P117, and M2-P87) aided in identifying considerable heterozygosity that would have otherwise gone undetected. Judging from their geographic distributions and network analyses of their associated STR profiles, 4 lineages appear to have developed in this region and to be of considerable age: K6-P79, K7-P117, M2-P87, and M2a-P22. The origins of K5-M230 and M-M4 are also confirmed as being located further west, probably in New Guinea. In the 25 adequately sampled populations, the number of different haplogroups ranged from 2 in the single most isolated group (the Aita of Bougainville), to 9, and measures of molecular diversity were generally not particularly low. The resulting pattern contradicts earlier findings that suggested far lower male-mediated diversity and gene exchange rates in the region. However, these earlier studies had not included the newly defined haplogroups. We could only identify a very weak signal of recent male Southeast Asian genetic influence (< 10%), which was almost entirely restricted to Austronesian (Oceanic)-speaking groups. This contradicts earlier assumptions on the ancestral composition of these groups and requires a revision of hypotheses concerning the settlement of the islands of the central Pacific, which commenced from this region.

To investigate associations between genetic, linguistic, and geographic variation in Africa, we type 50 Y chromosome SNPs in 1122 individuals from 40 populations representing African geographic and linguistic diversity. We compare these patterns of variation with those that emerge from a similar analysis of published mtDNA HVS1 sequences from 1918 individuals from 39 African populations. For the Y chromosome, Mantel tests reveal a strong partial correlation between genetic and linguistic distances (r=0.33, P=0.001) and no correlation between genetic and geographic distances (r=-0.08, P>0.10). In contrast, mtDNA variation is weakly correlated with both language (r=0.16, P=0.046) and geography (r=0.17, P=0.035). AMOVA indicates that the amount of paternal among-group variation is much higher when populations are grouped by linguistics (Phi(CT)=0.21) than by geography (Phi(CT)=0.06). Levels of maternal genetic among-group variation are low for both linguistics and geography (Phi(CT)=0.03 and 0.04, respectively). When Bantu speakers are removed from these analyses, the correlation with linguistic variation disappears for the Y chromosome and strengthens for mtDNA. These data suggest that patterns of differentiation and gene flow in Africa have differed for men and women in the recent evolutionary past. We infer that sex-biased rates of admixture and/or language borrowing between expanding Bantu farmers and local hunter-gatherers played an important role in influencing patterns of genetic variation during the spread of African agriculture in the last 4000 years.

Global-scale patterns of human population structure may be influenced by the rate of migration among populations that is nearly eight times higher for females than for males. This difference is attributed mainly to the widespread practice of patrilocality, in which women move into their mates' residences after marriage. Here we directly test this hypothesis by comparing global patterns of DNA sequence variation on the Y chromosome and mitochondrial DNA (mtDNA) in the same panel of 389 individuals from ten populations (four from Africa and two each from Europe, Asia and Oceania). We introduce a new strategy to assay Y-chromosome variation that identifies a high density of single-nucleotide polymorphisms, allows complete sequencing of all individuals rather than relying on predetermined markers and provides direct sequence comparisons with mtDNA. We found the overall proportion of between-group variation (Phi(ST)) to be 0.334 for the Y chromosome and 0.382 for mtDNA. Genetic differentiation between populations was similar for the Y chromosome and mtDNA at all geographic scales that we tested. Although patrilocality may be important at the local scale, patterns of genetic structure on the continental and global scales are not shaped by the higher rate of migration among females than among males.

Studies of human DNA sequence polymorphism reveal a range of diversity patterns throughout the genome. This variation among loci may be due to natural selection, demographic influences, and/or different sampling strategies. Here we build on a continuing study of noncoding regions on the X chromosome in a panel of 41 globally sampled humans representing African and non-African populations by examining patterns of DNA sequence variation at four loci (APXL, AMELX, TNFSF5, and RRM2P4) and comparing these patterns with those previously reported at six loci in the same panel of 41 individuals. We also include comparisons with patterns of noncoding variation seen at five additional X-linked loci that were sequenced in similar global panels. We find that, while almost all loci show a reduction in non-African diversity, the magnitude of the reduction varies substantially across loci. The large observed variance in non-African levels of diversity results in the rejection of a neutral model of molecular evolution with a multi-locus HKA test under both a constant size and a bottleneck model. In non-Africans, some loci harbor an excess of rare mutations over neutral equilibrium predictions, while other loci show no such deviation in the distribution of mutation frequencies. We also observe a positive relationship between recombination rate and frequency spectra in our non-African, but not in our African, sample. These results indicate that a simple out-of-Africa bottleneck model is not sufficient to explain the observed patterns of sequence variation and that diversity-reducing selection acting at a subset of loci and/or a more complex neutral model must be invoked.

The time to the most recent common ancestor (TMRCA) of the human mitochondria (mtDNA) is estimated to be older than that of the nonrecombining portion of the Y chromosome (NRY). Surveys of variation in globally distributed humans typically result in mtDNA TMRCA values just under 200 thousand years ago (kya), whereas those for the NRY range between 46 and 110 kya. A favored hypothesis for this finding is that natural selection has acted on the NRY, leading to a recent selective sweep. An alternate hypothesis is that sex-biased demographic processes are responsible. Here, we re-examine the disparity between NRY and mtDNA TMRCAs using data collected from individual human populations--a sampling strategy that minimizes the confounding influence of population subdivision in global data sets. We survey variation at 782 bp of the mitochondrial cytochrome c oxidase subunit 3 gene as well as at 26.5 kb of noncoding DNA from the NRY in a sample of 25 Khoisan, 24 Mongolians, and 24 Papua New Guineans. Data from both loci in all populations are best described by a model of constant population size, with the exception of Mongolian mtDNA, which appears to be experiencing rapid population growth. Taking these demographic models into account, we estimate the TMRCAs for each locus in each population. A pattern that is remarkably consistent across all three populations is an approximately twofold deeper coalescence for mtDNA than for the NRY. The oldest TMRCAs are observed for the Khoisan (73.6 kya for the NRY and 176.5 kya for mtDNA), whereas those in the non-African populations are consistently lower (averaging 47.7 kya for the NRY and 92.8 kya for mtDNA). Our data do not suggest that differential natural selection is the cause of this difference in TMRCAs. Rather, these results are most consistent with a higher female effective population size.