Priya Moorjani

The complete genome sequence of a Neanderthal from the Altai Mountains

We present a high-quality genome sequence of a Neanderthal woman from Siberia. We show that her parents were related at the level of half-siblings and that mating among close relatives was common among her recent ancestors. We also sequenced the genome of a Neanderthal from the Caucasus to low coverage. An analysis of the relationships and population history of available archaic genomes and 25 present-day human genomes shows that several gene flow events occurred among Neanderthals, Denisovans and early modern humans, possibly including gene flow into Denisovans from an unknown archaic group. Thus, interbreeding, albeit of low magnitude, occurred among many hominin groups in the Late Pleistocene. In addition, the high-quality Neanderthal genome allows us to establish a definitive list of substitutions that became fixed in modern humans after their separation from the ancestors of Neanderthals and Denisovans.

Ancient Admixture in Human History

Population mixture is an important process in biology. We present a suite of methods for learning about population mixtures, implemented in a software package called ADMIXTOOLS, that support formal tests for whether mixture occurred and make it possible to infer proportions and dates of mixture. We also describe the development of a new single nucleotide polymorphism (SNP) array consisting of 629,433 sites with clearly documented ascertainment that was specifically designed for population genetic analyses and that we genotyped in 934 individuals from 53 diverse populations. To illustrate the methods, we give a number of examples that provide new insights about the history of human admixture. The most striking finding is a clear signal of admixture into northern Europe, with one ancestral population related to present-day Basques and Sardinians and the other related to present-day populations of northeast Asia and the Americas. This likely reflects a history of admixture between Neolithic migrants and the indigenous Mesolithic population of Europe, consistent with recent analyses of ancient bones from Sweden and the sequencing of the genome of the Tyrolean “Iceman.”

Complete Khoisan and Bantu genomes from southern Africa

The genetic structure of the indigenous hunter-gatherer peoples of southern Africa, the oldest known lineage of modern human, is important for understanding human diversity. Studies based on mitochondrial and small sets of nuclear markers have shown that these hunter-gatherers, known as Khoisan, San, or Bushmen, are genetically divergent from other humans. However, until now, fully sequenced human genomes have been limited to recently diverged populations. Here we present the complete genome sequences of an indigenous hunter-gatherer from the Kalahari Desert and a Bantu from southern Africa, as well as protein-coding regions from an additional three hunter-gatherers from disparate regions of the Kalahari. We characterize the extent of whole-genome and exome diversity among the five men, reporting 1.3 million novel DNA differences genome-wide, including 13,146 novel amino acid variants. In terms of nucleotide substitutions, the Bushmen seem to be, on average, more different from each other than, for example, a European and an Asian. Observed genomic differences between the hunter-gatherers and others may help to pinpoint genetic adaptations to an agricultural lifestyle. Adding the described variants to current databases will facilitate inclusion of southern Africans in medical research efforts, particularly when family and medical histories can be correlated with genome-wide data.

Inferring Admixture Histories of Human Populations Using Linkage Disequilibrium

Long-range migrations and the resulting admixtures between populations have been important forces shaping human genetic diversity. Most existing methods for detecting and reconstructing historical admixture events are based on allele frequency divergences or patterns of ancestry segments in chromosomes of admixed individuals. An emerging new approach harnesses the exponential decay of admixture-induced linkage disequilibrium (LD) as a function of genetic distance. Here, we comprehensively develop LD-based inference into a versatile tool for investigating admixture. We present a new weighted LD statistic that can be used to infer mixture proportions as well as dates with fewer constraints on reference populations than previous methods. We define an LD-based three-population test for admixture and identify scenarios in which it can detect admixture events that previous formal tests cannot. We further show that we can uncover phylogenetic relationships among populations by comparing weighted LD curves obtained using a suite of references. Finally, we describe several improvements to the computation and fitting of weighted LD curves that greatly increase the robustness and speed of the calculations. We implement all of these advances in a software package, ALDER, which we validate in simulations and apply to test for admixture among all populations from the Human Genome Diversity Project (HGDP), highlighting insights into the admixture history of Central African Pygmies, Sardinians, and Japanese.

The History of African Gene Flow into Southern Europeans, Levantines, and Jews

Previous genetic studies have suggested a history of sub-Saharan African gene flow into some West Eurasian populations after the initial dispersal out of Africa that occurred at least 45,000 years ago. However, there has been no accurate characterization of the proportion of mixture, or of its date. We analyze genome-wide polymorphism data from about 40 West Eurasian groups to show that almost all Southern Europeans have inherited 1%–3% African ancestry with an average mixture date of around 55 generations ago, consistent with North African gene flow at the end of the Roman Empire and subsequent Arab migrations. Levantine groups harbor 4%–15% African ancestry with an average mixture date of about 32 generations ago, consistent with close political, economic, and cultural links with Egypt in the late middle ages. We also detect 3%–5% sub-Saharan African ancestry in all eight of the diverse Jewish populations that we analyzed. For the Jewish admixture, we obtain an average estimated date of about 72 generations. This may reflect descent of these groups from a common ancestral population that already had some African ancestry prior to the Jewish Diasporas.

Genome-wide association study of Tourette's syndrome

Tourettes syndrome (TS) is a developmental disorder that has one of the highest familial recurrence rates among neuropsychiatric diseases with complex inheritance. However, the identification of definitive TS susceptibility genes remains elusive. Here, we report the first genome-wide association study (GWAS) of TS in 1285 cases and 4964 ancestry-matched controls of European ancestry, including two European-derived population isolates, Ashkenazi Jews from North America and Israel and French Canadians from Quebec, Canada. In a primary meta-analysis of GWAS data from these European ancestry samples, no markers achieved a genome-wide threshold of significance (P<5 × 10−8); the top signal was found in rs7868992 on chromosome 9q32 within COL27A1 (P=1.85 × 10−6). A secondary analysis including an additional 211 cases and 285 controls from two closely related Latin American population isolates from the Central Valley of Costa Rica and Antioquia, Colombia also identified rs7868992 as the top signal (P=3.6 × 10−7 for the combined sample of 1496 cases and 5249 controls following imputation with 1000 Genomes data). This study lays the groundwork for the eventual identification of common TS susceptibility variants in larger cohorts and helps to provide a more complete understanding of the full genetic architecture of this disorder.

Pharmacogenetic analysis of genes implicated in rodent models of antidepressant response: association of TREK1 and treatment resistance in the STAR(*)D study.

Recent rodent models of antidepressant response implicate a novel set of genes in mechanisms of antidepressant action. The authors examined variants in four such genes (KCNK2 (TREK1), SLC18A2 (VMAT2), S100A10, and HDAC5) for association with remission in a large effectiveness trial of antidepressant treatments. Subjects were drawn from the Sequenced Treatment Alternatives to Relieve Depression (STAR*D) study, a multicenter, prospective, effectiveness trial in major depressive disorder (MDD). Outpatients with nonpsychotic MDD were initially treated with citalopram for up to 14 weeks; those who did not remit with citalopram were sequentially randomized to a series of next-step treatments, each for up to 12 weeks. Single-nucleotide polymorphisms in four genes were examined for association with remission, defined as a clinician-rated Quick Inventory of Depressive Symptomatology (QIDS-C16) score ⩽5. Of 1554 participants for whom DNA was available, 565 (36%) reached remission with citalopram treatment. No association with any of the four genes was identified. However, among the 751 who entered next-step treatment, variants in KCNK2 were associated with treatment response (Bonferroni-corrected, gene-based empirical p<0.001). In follow-up analyses, KCNK2 was also associated with effects of similar magnitude for third-step treatment among those with unsatisfactory benefit to both citalopram and one next-step pharmacotherapy (n=225). These findings indicate that genetic variation in KCNK2 may identify individuals at risk for treatment resistance. More broadly, they indicate the utility of animal models in identifying genes for pharmacogenetic studies of antidepressant response.

Lack of association between SLITRK1var321 and Tourette syndrome in a large family-based sample.

Tourette syndrome (TS) has a significant genetic component, yet no TS susceptibility genes have been identified definitively. Several studies have determined that first-degree relatives of patients with TS have at least a 5- to 15-fold increased risk of developing the disorder compared with the general population, an increase that represents one of the highest familial recurrence risks among neuropsychiatric diseases that are inherited in a non-Mendelian fashion.1 Recently, Slit- and Trk-like 1 (SLITRK1) was proposed as a candidate TS susceptibility gene, and a noncoding polymorphism in the 3′ untranslated region of this gene (var321) was reported to be associated with TS in a case–control sample.2 Additional studies in small samples or population isolates have failed to replicate this association.3,4 As part of a 20-year collaborative effort, the Tourette Syndrome Association International Consortium for Genetics (TSAICG) has systematically collected a clinical sample of over 1,000 patients with TS and their family members.5 We chose to screen these individuals for SLITRK1 var321 to determine a more accurate estimate of the prevalence of this variant in the white TS clinic population and to test for any association between var321 and TS. ### Methods. A total of 2,300 individuals from 646 independently ascertained nuclear families were recruited from tic disorder specialty clinics from the United States, Canada, Great Britain, and The Netherlands. A total of 1,048 individuals (172 parents and 876 offspring) were diagnosed with either TS (989 subjects) or chronic tics (CT) (59 subjects) (e-Methods on the Neurology ® Web site at www.neurology.org). A total of 440 …

Reconstructing Austronesian population history in Island Southeast Asia

Austronesian languages are spread across half the globe, from Easter Island to Madagascar. Evidence from linguistics and archaeology indicates that the ‘Austronesian expansion,’ which began 4,000–5,000 years ago, likely had roots in Taiwan, but the ancestry of present-day Austronesian-speaking populations remains controversial. Here, we analyse genome-wide data from 56 populations using new methods for tracing ancestral gene flow, focusing primarily on Island Southeast Asia. We show that all sampled Austronesian groups harbour ancestry that is more closely related to aboriginal Taiwanese than to any present-day mainland population. Surprisingly, western Island Southeast Asian populations have also inherited ancestry from a source nested within the variation of present-day populations speaking Austro-Asiatic languages, which have historically been nearly exclusive to the mainland. Thus, either there was once a substantial Austro-Asiatic presence in Island Southeast Asia, or Austronesian speakers migrated to and through the mainland, admixing there before continuing to western Indonesia.

Corticotrophin-releasing hormone type 1 receptor gene (CRHR1) variants predict posttraumatic stress disorder onset and course in pediatric injury patients

Posttraumatic stress disorder (PTSD) is a common and disabling anxiety disorder that may occur in the aftermath of exposure to potentially traumatic life events. PTSD is moderately heritable, but few specific molecular variants accounting for this heritability have been identified. Genes regulating the hypothalamic-pituitary-adrenal (HPA) axis, such as corticotrophin-releasing hormone type 1 receptor gene (CRHR1), have been implicated in traumatic-stress related phenotypes but have yet to be studied in relation to PTSD. The present study sought to examine the relation between 9 single nucleotide polymorphisms (SNPs) in the CRHR1 gene and posttraumatic stress symptoms in a prospective study of pediatric injury patients (n = 103) who were first assessed in the acute aftermath of their injury at the hospital. Results indicated that multiple SNPs were associated with acute symptoms at a univariate level, and after correction for multiple testing, rs12944712 was significantly related to acute PTSD symptoms. Longitudinal latent growth curve analyses suggest that rs12944712 is also related to both acute symptom level and trajectory of symptoms over time. The present study adds support for the role of CRHR1 in the stress response following potentially traumatic event exposure in youth. It should be noted that the sample size in this study was small, and therefore statistical power was low; following, results from this study should be considered preliminary. Although results are not definitive, the findings from this study warrant future replication studies on how variation in this gene relates to response to traumatic event exposure in youth.