1
|
Villanea FA, Peede D, Kaufman EJ, Añorve-Garibay V, Witt KE, Villa-Islas V, Zeloni R, Marnetto D, Moorjani P, Jay F, Valdmanis PN, Ávila-Arcos MC, Huerta-Sánchez E. The MUC19 gene in Denisovans, Neanderthals, and Modern Humans: An Evolutionary History of Recurrent Introgression and Natural Selection. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.09.25.559202. [PMID: 37808839 PMCID: PMC10557577 DOI: 10.1101/2023.09.25.559202] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/10/2023]
Abstract
All humans carry a small fraction of archaic ancestry across the genome, the legacy of gene flow from Neanderthals, Denisovans, and other hominids into the ancestors of modern humans. While the effects of Neanderthal ancestry on human fitness and health have been explored more thoroughly, there are fewer examples of adaptive introgression of Denisovan variants. Here, we study the gene MUC19, for which some modern humans carry a Denisovan-like haplotype. MUC19 is a mucin, a glycoprotein that forms gels with various biological functions, from lubrication to immunity. We find the diagnostic variants for the Denisovan-like MUC19 haplotype at high frequencies in admixed Latin American individuals among global population, and at highest frequency in 23 ancient Indigenous American individuals, all predating population admixture with Europeans and Africans. We find that some Neanderthals--Vindija and Chagyrskaya--carry the Denisovan-like MUC19 haplotype, and that it was likely introgressed into human populations through Neanderthal introgression rather than Denisovan introgression. Finally, we find that the Denisovan-like MUC19 haplotype carries a higher copy number of a 30 base-pair variable number tandem repeat relative to the Human-like haplotype, and that copy numbers of this repeat are exceedingly high in American populations. Our results suggest that the Denisovan-like MUC19 haplotype served as the raw genetic material for positive selection as American populations adapted to novel environments during their movement from Beringia into North and then South America.
Collapse
Affiliation(s)
| | - David Peede
- Department of Ecology, Evolution, and Organismal Biology, Brown University
- Center for Computational Molecular Biology, Brown University
- Institute at Brown for Environment and Society, Brown University
| | - Eli J Kaufman
- Division of Medical Genetics, Department of Medicine, University of Washington School of Medicine
| | - Valeria Añorve-Garibay
- Center for Computational Molecular Biology, Brown University
- International Laboratory for Human Genome Research, Universidad Nacional Autónoma de México
| | - Kelsey E Witt
- Center for Human Genetics and Department of Genetics and Biochemistry, Clemson University
| | - Viridiana Villa-Islas
- International Laboratory for Human Genome Research, Universidad Nacional Autónoma de México
| | - Roberta Zeloni
- Department of Neurosciences "Rita Levi Montalcini", University of Turin
| | - Davide Marnetto
- Department of Neurosciences "Rita Levi Montalcini", University of Turin
| | - Priya Moorjani
- Department of Molecular and Cell Biology, University of California, Berkeley
- Center for Computational Biology, University of California, Berkeley
| | - Flora Jay
- Université Paris-Saclay, CNRS, INRIA, Laboratoire Interdisciplinaire des Sciences du Numérique, 91400, Orsay, France
| | - Paul N Valdmanis
- Division of Medical Genetics, Department of Medicine, University of Washington School of Medicine
| | - María C Ávila-Arcos
- International Laboratory for Human Genome Research, Universidad Nacional Autónoma de México
| | - Emilia Huerta-Sánchez
- Department of Ecology, Evolution, and Organismal Biology, Brown University
- Center for Computational Molecular Biology, Brown University
| |
Collapse
|
2
|
Thawornwattana Y, Huang J, Flouri T, Mallet J, Yang Z. Inferring the Direction of Introgression Using Genomic Sequence Data. Mol Biol Evol 2023; 40:msad178. [PMID: 37552932 PMCID: PMC10439365 DOI: 10.1093/molbev/msad178] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2023] [Revised: 08/01/2023] [Accepted: 08/02/2023] [Indexed: 08/10/2023] Open
Abstract
Genomic data are informative about the history of species divergence and interspecific gene flow, including the direction, timing, and strength of gene flow. However, gene flow in opposite directions generates similar patterns in multilocus sequence data, such as reduced sequence divergence between the hybridizing species. As a result, inference of the direction of gene flow is challenging. Here, we investigate the information about the direction of gene flow present in genomic sequence data using likelihood-based methods under the multispecies-coalescent-with-introgression model. We analyze the case of two species, and use simulation to examine cases with three or four species. We find that it is easier to infer gene flow from a small population to a large one than in the opposite direction, and easier to infer inflow (gene flow from outgroup species to an ingroup species) than outflow (gene flow from an ingroup species to an outgroup species). It is also easier to infer gene flow if there is a longer time of separate evolution between the initial divergence and subsequent introgression. When introgression is assumed to occur in the wrong direction, the time of introgression tends to be correctly estimated and the Bayesian test of gene flow is often significant, while estimates of introgression probability can be even greater than the true probability. We analyze genomic sequences from Heliconius butterflies to demonstrate that typical genomic datasets are informative about the direction of interspecific gene flow, as well as its timing and strength.
Collapse
Affiliation(s)
| | - Jun Huang
- School of Biomedical Engineering, Capital Medical University, Beijing 100069, P.R. China
| | - Tomáš Flouri
- Department of Genetics, Evolution and Environment, University College London, London WC1E 6BT, UK
| | - James Mallet
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA 02138, USA
- Department of Genetics, Evolution and Environment, University College London, London WC1E 6BT, UK
| | - Ziheng Yang
- Department of Genetics, Evolution and Environment, University College London, London WC1E 6BT, UK
| |
Collapse
|
3
|
Ji J, Jackson DJ, Leaché AD, Yang Z. Power of Bayesian and Heuristic Tests to Detect Cross-Species Introgression with Reference to Gene Flow in the Tamias quadrivittatus Group of North American Chipmunks. Syst Biol 2023; 72:446-465. [PMID: 36504374 PMCID: PMC10275556 DOI: 10.1093/sysbio/syac077] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2021] [Revised: 11/15/2022] [Accepted: 12/01/2022] [Indexed: 10/25/2023] Open
Abstract
In the past two decades, genomic data have been widely used to detect historical gene flow between species in a variety of plants and animals. The Tamias quadrivittatus group of North America chipmunks, which originated through a series of rapid speciation events, are known to undergo massive amounts of mitochondrial introgression. Yet in a recent analysis of targeted nuclear loci from the group, no evidence for cross-species introgression was detected, indicating widespread cytonuclear discordance. The study used the heuristic method HYDE to detect gene flow, which may suffer from low power. Here we use the Bayesian method implemented in the program BPP to re-analyze these data. We develop a Bayesian test of introgression, calculating the Bayes factor via the Savage-Dickey density ratio using the Markov chain Monte Carlo (MCMC) sample under the model of introgression. We take a stepwise approach to constructing an introgression model by adding introgression events onto a well-supported binary species tree. The analysis detected robust evidence for multiple ancient introgression events affecting the nuclear genome, with introgression probabilities reaching 63%. We estimate population parameters and highlight the fact that species divergence times may be seriously underestimated if ancient cross-species gene flow is ignored in the analysis. We examine the assumptions and performance of HYDE and demonstrate that it lacks power if gene flow occurs between sister lineages or if the mode of gene flow does not match the assumed hybrid-speciation model with symmetrical population sizes. Our analyses highlight the power of likelihood-based inference of cross-species gene flow using genomic sequence data. [Bayesian test; BPP; chipmunks; introgression; MSci; multispecies coalescent; Savage-Dickey density ratio.].
Collapse
Affiliation(s)
- Jiayi Ji
- Department of Genetics, Evolution and Environment, University College London, London WC1E 6BT, UK
| | - Donavan J Jackson
- Department of Biology and Burke Museum of Natural History and Culture, University of Washington, Box 351800, Seattle, WA 98195-1800, USA
| | - Adam D Leaché
- Department of Biology and Burke Museum of Natural History and Culture, University of Washington, Box 351800, Seattle, WA 98195-1800, USA
| | - Ziheng Yang
- Department of Genetics, Evolution and Environment, University College London, London WC1E 6BT, UK
| |
Collapse
|
4
|
Yan Z, Ogilvie HA, Nakhleh L. Comparing inference under the multispecies coalescent with and without recombination. Mol Phylogenet Evol 2023; 181:107724. [PMID: 36720421 DOI: 10.1016/j.ympev.2023.107724] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/28/2022] [Revised: 01/17/2023] [Accepted: 01/24/2023] [Indexed: 01/29/2023]
Abstract
Accurate inference of population parameters plays a pivotal role in unravelling evolutionary histories. While recombination has been universally accepted as a fundamental process in the evolution of sexually reproducing organisms, it remains challenging to model it exactly. Thus, existing coalescent-based approaches make different assumptions or approximations to facilitate phylogenetic inference, which can potentially bring about biases in estimates of evolutionary parameters when recombination is present. In this article, we evaluate the performance of population parameter estimation using three methods-StarBEAST2, SNAPP, and diCal2-that represent three different types of inference. We performed whole-genome simulations in which recombination rates, mutation rates, and levels of incomplete lineage sorting were varied. We show that StarBEAST2 using short or medium-sized loci is robust to realistic rates of recombination, which is in agreement with previous studies. SNAPP, as expected, is generally unaffected by recombination events. Most surprisingly, diCal2, a method that is designed to explicitly account for recombination, performs considerably worse than other methods under comparison.
Collapse
Affiliation(s)
- Zhi Yan
- Department of Computer Science, Rice University, 6100 Main Street, Houston 77005, TX, USA.
| | - Huw A Ogilvie
- Department of Computer Science, Rice University, 6100 Main Street, Houston 77005, TX, USA.
| | - Luay Nakhleh
- Department of Computer Science, Rice University, 6100 Main Street, Houston 77005, TX, USA.
| |
Collapse
|
5
|
Keeling BA, Quam R, Martínez I, Arsuaga JL, Maroto J. Reassessment of the human mandible from Banyoles (Girona, Spain). J Hum Evol 2023; 174:103291. [PMID: 36493597 DOI: 10.1016/j.jhevol.2022.103291] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/09/2021] [Revised: 10/18/2022] [Accepted: 10/21/2022] [Indexed: 12/12/2022]
Abstract
Since the discovery of a human mandible in 1887 near the present-day city of Banyoles, northeastern Spain, researchers have generally emphasized its archaic features, including the lack of chin structures, and suggested affinities with the Neandertals or European Middle Pleistocene (Chibanian) specimens. Uranium-series and electron spin resonance dating suggest the mandible dates to the Late Pleistocene (Tarantian), approximately ca. 45-66 ka. In this study, we reassessed the taxonomic affinities of the Banyoles mandible by comparing it to samples of Middle Pleistocene fossils from Africa and Europe, Neandertals, Early and Upper Paleolithic modern humans, and recent modern humans. We evaluated the frequencies and expressions of morphological features and performed a three-dimensional geometric morphometric analysis on a virtual reconstruction of Banyoles to capture overall mandibular shape. Our results revealed no derived Neandertal morphological features in Banyoles. While a principal component analysis based on Euclidean distances from the first two principal components clearly grouped Banyoles with both fossil and recent Homo sapiens individuals, an analysis of the Procrustes residuals demonstrated that Banyoles did not fit into any of the comparative groups. The lack of Neandertal features in Banyoles is surprising considering its Late Pleistocene age. A consideration of the Middle Pleistocene fossil record in Europe and southwest Asia suggests that Banyoles is unlikely to represent a late-surviving Middle Pleistocene population. The lack of chin structures also complicates an assignment to H. sapiens, although early fossil H. sapiens do show somewhat variable development of the chin structures. Thus, Banyoles represents a non-Neandertal Late Pleistocene European individual and highlights the continuing signal of diversity in the hominin fossil record. The present situation makes Banyoles a prime candidate for ancient DNA or proteomic analyses, which may shed additional light on its taxonomic affinities.
Collapse
Affiliation(s)
- Brian A Keeling
- Department of Anthropology, Binghamton University, SUNY, New York, USA.
| | - Rolf Quam
- Department of Anthropology, Binghamton University, SUNY, New York, USA; Centro UCM-ISCIII de Evolución y Comportamiento Humanos, Madrid, Spain; Division of Anthropology, American Museum of Natural History, New York, USA; Cátedra de Otoacústica Evolutiva y Paleoantropología (HM Hospitales-Universidad de Alcalá), Departamento de Ciencias de la Vida, Universidad de Alcalá, Madrid, Spain
| | - Ignacio Martínez
- Cátedra de Otoacústica Evolutiva y Paleoantropología (HM Hospitales-Universidad de Alcalá), Departamento de Ciencias de la Vida, Universidad de Alcalá, Madrid, Spain; Centro de Investigación Francisco Javier Muñiz, Universidad de Buenos Aires, Calle Paraguay 2155, Primer piso, Ciudad Autónoma de Buenos Aires, 1121, Argentina
| | - Juan Luis Arsuaga
- Centro UCM-ISCIII de Evolución y Comportamiento Humanos, Madrid, Spain; Departamento de Geodinámica, Estratigrafía y Paleontología, Facultad de Ciencias Geológicas, Universidad Complutense de Madrid, Madrid, Spain
| | - Julià Maroto
- Grup d'Arqueologia i Prehistòria, Universitat de Girona, pl. Ferrater Mora, 1, 17071 Girona, Spain
| |
Collapse
|
6
|
Yang Z, Flouri T. Estimation of Cross-Species Introgression Rates using Genomic Data Despite Model Unidentifiability. Mol Biol Evol 2022; 39:6568285. [PMID: 35417543 PMCID: PMC9087891 DOI: 10.1093/molbev/msac083] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/20/2022] Open
Abstract
Full-likelihood implementations of the multispecies coalescent with introgression (MSci) model treat genealogical fluctuations across the genome as a major source of information to infer the history of species divergence and gene flow using multilocus sequence data. However, MSci models are known to have unidentifiability issues, whereby different models or parameters make the same predictions about the data and cannot be distinguished by the data. Previous studies of unidentifiability have focused on heuristic methods based on gene trees and do not make an efficient use of the information in the data. Here we study the unidentifiability of MSci models under the full-likelihood methods. We characterize the unidentifiability of the bidirectional introgression (BDI) model, which assumes that gene flow occurs in both directions. We derive simple rules for arbitrary BDI models, which create unidentifiability of the label-switching type. In general, an MSci model with k BDI events has 2k unidentifiable modes or towers in the posterior, with each BDI event between sister species creating within-model parameter unidentifiability and each BDI event between nonsister species creating between-model unidentifiability. We develop novel algorithms for processing Markov chain Monte Carlo samples to remove label-switching problems and implement them in the bpp program. We analyze real and synthetic data to illustrate the utility of the BDI models and the new algorithms. We discuss the unidentifiability of heuristic methods and provide guidelines for the use of MSci models to infer gene flow using genomic data.
Collapse
Affiliation(s)
- Ziheng Yang
- Department of Genetics, Evolution and Environment, University College London, Gower Street, London WC1E6BT, UK
| | - Tomáš Flouri
- Department of Genetics, Evolution and Environment, University College London, Gower Street, London WC1E6BT, UK
| |
Collapse
|
7
|
Zhu T, Flouri T, Yang Z. A simulation study to examine the impact of recombination on phylogenomic inferences under the multispecies coalescent model. Mol Ecol 2022; 31:2814-2829. [PMID: 35313033 PMCID: PMC9321900 DOI: 10.1111/mec.16433] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2021] [Revised: 01/25/2022] [Accepted: 02/28/2022] [Indexed: 11/28/2022]
Affiliation(s)
- Tianqi Zhu
- Institute of Applied Mathematics Academy of Mathematics and Systems Science Chinese Academy of Sciences Beijing 100190 China
- Key Laboratory of Random Complex Structures and Data Science, Academy of Mathematics and Systems Science, Chinese Academy of Sciences Beijing 100190 China
| | - Tomáš Flouri
- Department of Genetics, Evolution and Environment University College London London WC1E 6BT UK
| | - Ziheng Yang
- Department of Genetics, Evolution and Environment University College London London WC1E 6BT UK
| |
Collapse
|
8
|
Tricou T, Tannier E, de Vienne DM. OUP accepted manuscript. Syst Biol 2022; 71:1147-1158. [PMID: 35169846 PMCID: PMC9366450 DOI: 10.1093/sysbio/syac011] [Citation(s) in RCA: 21] [Impact Index Per Article: 10.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/01/2021] [Revised: 02/01/2021] [Accepted: 02/08/2022] [Indexed: 11/29/2022] Open
Abstract
Most species are extinct, those that are not are often unknown. Sequenced and sampled species are often a minority of known ones. Past evolutionary events involving horizontal gene flow, such as horizontal gene transfer, hybridization, introgression, and admixture, are therefore likely to involve “ghosts,” that is extinct, unknown, or unsampled lineages. The existence of these ghost lineages is widely acknowledged, but their possible impact on the detection of gene flow and on the identification of the species involved is largely overlooked. It is generally considered as a possible source of error that, with reasonable approximation, can be ignored. We explore the possible influence of absent species on an evolutionary study by quantifying the effect of ghost lineages on introgression as detected by the popular D-statistic method. We show from simulated data that under certain frequently encountered conditions, the donors and recipients of horizontal gene flow can be wrongly identified if ghost lineages are not taken into account. In particular, having a distant outgroup, which is usually recommended, leads to an increase in the error probability and to false interpretations in most cases. We conclude that introgression from ghost lineages should be systematically considered as an alternative possible, even probable, scenario. [ABBA–BABA; D-statistic; gene flow; ghost lineage; introgression; simulation.]
Collapse
Affiliation(s)
- Théo Tricou
- Correspondence to be sent to: CNRS Université Claude Bernard Lyon 1, Laboratoire de Biométrie et Biologie Évolutive (LBBE), Bâtiment Mendel, 43 boulevard du 11 Novembre 1918, Villeurbanne, 69622 Cedex, France; E-mail:
| | - Eric Tannier
- Laboratoire de Biométrie et Biologie Évolutive UMR5558, Univ Lyon, Université Lyon 1, CNRS, F-69622 Villeurbanne, France
- Inria, Centre de Recherche de Lyon, F-69603 Villeurbanne, France
| | - Damien M de Vienne
- Laboratoire de Biométrie et Biologie Évolutive UMR5558, Univ Lyon, Université Lyon 1, CNRS, F-69622 Villeurbanne, France
| |
Collapse
|
9
|
Hibbins MS, Hahn MW. Phylogenomic approaches to detecting and characterizing introgression. Genetics 2021; 220:6425633. [PMID: 34788444 PMCID: PMC9208645 DOI: 10.1093/genetics/iyab173] [Citation(s) in RCA: 41] [Impact Index Per Article: 13.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/16/2021] [Accepted: 10/02/2021] [Indexed: 12/26/2022] Open
Abstract
Phylogenomics has revealed the remarkable frequency with which introgression occurs across the tree of life. These discoveries have been enabled by the rapid growth of methods designed to detect and characterize introgression from whole-genome sequencing data. A large class of phylogenomic methods makes use of data across species to infer and characterize introgression based on expectations from the multispecies coalescent. These methods range from simple tests, such as the D-statistic, to model-based approaches for inferring phylogenetic networks. Here, we provide a detailed overview of the various signals that different modes of introgression are expected leave in the genome, and how current methods are designed to detect them. We discuss the strengths and pitfalls of these approaches and identify areas for future development, highlighting the different signals of introgression, and the power of each method to detect them. We conclude with a discussion of current challenges in inferring introgression and how they could potentially be addressed.
Collapse
Affiliation(s)
- Mark S Hibbins
- Department of Biology, Indiana University, Bloomington, IN 47405, USA
| | - Matthew W Hahn
- Department of Biology, Indiana University, Bloomington, IN 47405, USA.,Department of Computer Science, Indiana University, Bloomington, IN 47405, USA
| |
Collapse
|
10
|
Gopalan S, Atkinson EG, Buck LT, Weaver TD, Henn BM. Inferring archaic introgression from hominin genetic data. Evol Anthropol 2021; 30:199-220. [PMID: 33951239 PMCID: PMC8360192 DOI: 10.1002/evan.21895] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2019] [Revised: 08/03/2020] [Accepted: 03/29/2021] [Indexed: 01/05/2023]
Abstract
Questions surrounding the timing, extent, and evolutionary consequences of archaic admixture into human populations have a long history in evolutionary anthropology. More recently, advances in human genetics, particularly in the field of ancient DNA, have shed new light on the question of whether or not Homo sapiens interbred with other hominin groups. By the late 1990s, published genetic work had largely concluded that archaic groups made no lasting genetic contribution to modern humans; less than a decade later, this conclusion was reversed following the successful DNA sequencing of an ancient Neanderthal. This reversal of consensus is noteworthy, but the reasoning behind it is not widely understood across all academic communities. There remains a communication gap between population geneticists and paleoanthropologists. In this review, we endeavor to bridge this gap by outlining how technological advancements, new statistical methods, and notable controversies ultimately led to the current consensus.
Collapse
Affiliation(s)
- Shyamalika Gopalan
- Department of Ecology and Evolution, Stony Brook University, Stony Brook, New York, USA.,Department of Evolutionary Anthropology, Duke University, Durham, North Carolina, USA
| | - Elizabeth G Atkinson
- Department of Ecology and Evolution, Stony Brook University, Stony Brook, New York, USA.,Analytic and Translational Genetics Unit, Massachusetts General Hospital and Stanley Center for Psychiatric Research, Broad Institute, Boston, Massachusetts, USA
| | - Laura T Buck
- Research Centre in Evolutionary Anthropology and Palaeoecology, Liverpool John Moores University, Liverpool, UK
| | - Timothy D Weaver
- Department of Anthropology, University of California, Davis, California, USA
| | - Brenna M Henn
- Department of Ecology and Evolution, Stony Brook University, Stony Brook, New York, USA.,Department of Anthropology, University of California, Davis, California, USA.,UC Davis Genome Center, University of California, Davis, California, USA
| |
Collapse
|
11
|
Inference of gene flow in the process of speciation: Efficient maximum-likelihood implementation of a generalised isolation-with-migration model. Theor Popul Biol 2021; 140:1-15. [PMID: 33736959 DOI: 10.1016/j.tpb.2021.03.001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2020] [Revised: 02/28/2021] [Accepted: 03/01/2021] [Indexed: 11/21/2022]
Abstract
The 'isolation with migration' (IM) model has been extensively used in the literature to detect gene flow during the process of speciation. In this model, an ancestral population split into two or more descendant populations which subsequently exchanged migrants at a constant rate until the present. Of course, the assumption of constant gene flow until the present is often over-simplistic in the context of speciation. In this paper, we consider a 'generalised IM' (GIM) model: a two-population IM model in which migration rates and population sizes are allowed to change at some point in the past. By developing a maximum-likelihood implementation of this model, we enable inference on both historical and contemporary rates of gene flow between two closely related populations or species. The GIM model encompasses both the standard two-population IM model and the 'isolation with initial migration' (IIM) model as special cases, as well as a model of secondary contact. We examine for simulated data how our method can be used, by means of likelihood ratio tests or AIC scores, to distinguish between the following scenarios of population divergence: (a) divergence in complete isolation; (b) divergence with a period of gene flow followed by isolation; (c) divergence with a period of isolation followed by secondary contact; (d) divergence with ongoing gene flow. Our method is based on the coalescent and is suitable for data sets consisting of the number of nucleotide differences between one pair of DNA sequences at each of a large number of independent loci. As our method relies on an explicit expression for the likelihood, it is computationally very fast.
Collapse
|
12
|
Forsythe ES, Nelson ADL, Beilstein MA. Biased Gene Retention in the Face of Introgression Obscures Species Relationships. Genome Biol Evol 2020; 12:1646-1663. [PMID: 33011798 PMCID: PMC7533067 DOI: 10.1093/gbe/evaa149] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 07/10/2020] [Indexed: 12/13/2022] Open
Abstract
Phylogenomic analyses are recovering previously hidden histories of hybridization, revealing the genomic consequences of these events on the architecture of extant genomes. We applied phylogenomic techniques and several complementary statistical tests to show that introgressive hybridization appears to have occurred between close relatives of Arabidopsis, resulting in cytonuclear discordance and impacting our understanding of species relationships in the group. The composition of introgressed and retained genes indicates that selection against incompatible cytonuclear and nuclear-nuclear interactions likely acted during introgression, whereas linkage also contributed to genome composition through the retention of ancient haplotype blocks. We also applied divergence-based tests to determine the species branching order and distinguish donor from recipient lineages. Surprisingly, these analyses suggest that cytonuclear discordance arose via extensive nuclear, rather than cytoplasmic, introgression. If true, this would mean that most of the nuclear genome was displaced during introgression whereas only a small proportion of native alleles were retained.
Collapse
|
13
|
Cooper L, Bunnefeld L, Hearn J, Cook JM, Lohse K, Stone GN. Low-coverage genomic data resolve the population divergence and gene flow history of an Australian rain forest fig wasp. Mol Ecol 2020; 29:3649-3666. [PMID: 32567765 DOI: 10.1111/mec.15523] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/21/2020] [Revised: 06/09/2020] [Accepted: 06/12/2020] [Indexed: 12/16/2022]
Abstract
Population divergence and gene flow are key processes in evolution and ecology. Model-based analysis of genome-wide data sets allows discrimination between alternative scenarios for these processes even in nonmodel taxa. We used two complementary approaches (one based on the blockwise site frequency spectrum [bSFS], the second on the pairwise sequentially Markovian coalescent [PSMC]) to infer the divergence history of a fig wasp, Pleistodontes nigriventris. Pleistodontes nigriventris and its fig tree mutualist Ficus watkinsiana are restricted to rain forest patches along the eastern coast of Australia and are separated into The Northern population is to the north of the Southern populations by two dry forest corridors (the Burdekin and St. Lawrence Gaps). We generated whole genome sequence data for two haploid males per population and used the bSFS approach to infer the timing of divergence between northern and southern populations of P. nigriventris, and to discriminate between alternative isolation with migration (IM) and instantaneous admixture (ADM) models of postdivergence gene flow. Pleistodontes nigriventris has low genetic diversity (π = 0.0008), to our knowledge one of the lowest estimates reported for a sexually reproducing arthropod. We find strongest support for an ADM model in which the two populations diverged ca. 196 kya in the late Pleistocene, with almost 25% of northern lineages introduced from the south during an admixture event ca. 57 kya. This divergence history is highly concordant with individual population demographies inferred from each pair of haploid males using PSMC. Our analysis illustrates the inferences possible with genome-level data for small population samples of tiny, nonmodel organisms and adds to a growing body of knowledge on the population structure of Australian rain forest taxa.
Collapse
Affiliation(s)
- Lisa Cooper
- Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, UK
| | - Lynsey Bunnefeld
- Biological and Environmental Sciences, University of Stirling, Stirling, UK
| | - Jack Hearn
- Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, UK.,Vector Biology Department, Liverpool School of Tropical Medicine, Liverpool, UK
| | - James M Cook
- Hawkesbury Institute for the Environment, Western Sydney University, Richmond, NSW, Australia
| | - Konrad Lohse
- Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, UK
| | - Graham N Stone
- Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, UK
| |
Collapse
|
14
|
Gokcumen O. Archaic hominin introgression into modern human genomes. AMERICAN JOURNAL OF PHYSICAL ANTHROPOLOGY 2019; 171 Suppl 70:60-73. [PMID: 31702050 DOI: 10.1002/ajpa.23951] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/02/2019] [Revised: 10/04/2019] [Accepted: 10/08/2019] [Indexed: 01/01/2023]
Abstract
Ancient genomes from multiple Neanderthal and the Denisovan individuals, along with DNA sequence data from diverse contemporary human populations strongly support the prevalence of gene flow among different hominins. Recent studies now provide evidence for multiple gene flow events that leave genetic signatures in extant and ancient human populations. These events include older gene flow from an unknown hominin in Africa predating out-of-Africa migrations, and in the last 50,000-100,000 years, multiple gene flow events from Neanderthals into ancestral Eurasian human populations, and at least three distinct introgression events from a lineage close to Denisovans into ancestors of extant Southeast Asian and Oceanic populations. Some of these introgression events may have happened as late as 20,000 years before present and reshaped the way in which we think about human evolution. In this review, I aim to answer anthropologically relevant questions with regard to recent research on ancient hominin introgression in the human lineage. How have genomic data from archaic hominins changed our view of human evolution? Is there any doubt about whether introgression from ancient hominins to the ancestors of present-day humans occurred? What is the current view of human evolutionary history from the genomics perspective? What is the impact of introgression on human phenotypes?
Collapse
Affiliation(s)
- Omer Gokcumen
- Department of Biological Sciences, North Campus, University at Buffalo, Buffalo, New York
| |
Collapse
|
15
|
Abstract
Abstract
Many methods exist for detecting introgression between nonsister species, but the most commonly used require either a single sequence from four or more taxa or multiple sequences from each of three taxa. Here, we present a test for introgression that uses only a single sequence from three taxa. This test, denoted D3, uses similar logic as the standard D-test for introgression, but by using pairwise distances instead of site patterns it is able to detect the same signal of introgression with fewer species. We use simulations to show that D3 has statistical power almost equal to D, demonstrating its use on a data set of wild bananas (Musa). The new test is easy to apply and easy to interpret, and should find wide use among currently available data sets.
Collapse
Affiliation(s)
- Matthew W Hahn
- Department of Biology, Indiana University, Bloomington, IN
- Department of Computer Science, Indiana University, Bloomington, IN
| | - Mark S Hibbins
- Department of Biology, Indiana University, Bloomington, IN
| |
Collapse
|
16
|
Shultz DR, Montrey M, Shultz TR. Comparing fitness and drift explanations of Neanderthal replacement. Proc Biol Sci 2019; 286:20190907. [PMID: 31185865 DOI: 10.1098/rspb.2019.0907] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022] Open
Abstract
There is a general consensus among archaeologists that replacement of Neanderthals by anatomically modern humans in Europe occurred around 40-35 ka. However, the causal mechanism for this replacement continues to be debated. Proposed models have featured either fitness advantages in favour of anatomically modern humans or invoked neutral drift under various preconditions. Searching for specific fitness advantages in the archaeological record has proven difficult, as these may be obscured, absent or subject to interpretation. To bridge this gap, we rigorously compare the system-level properties of fitness- and drift-based explanations of Neanderthal replacement. Our stochastic simulations and analytical predictions show that, although both fitness and drift can produce replacement, they present important differences in (i) required initial conditions, (ii) reliability, (iii) time to replacement, and (iv) path to replacement (population histories). These results present useful opportunities for comparison with archaeological and genetic data. We find greater agreement between the available empirical evidence and the system-level properties of replacement by differential fitness, rather than by neutral drift.
Collapse
Affiliation(s)
- Daniel R Shultz
- 1 Department of Anthropology, McGill University , Montreal, Quebec , Canada.,2 Department of History, McGill University , Montreal, Quebec , Canada
| | - Marcel Montrey
- 3 Department of Psychology, McGill University , Montreal, Quebec , Canada
| | - Thomas R Shultz
- 3 Department of Psychology, McGill University , Montreal, Quebec , Canada.,4 School of Computer Science, McGill University , Montreal, Quebec , Canada
| |
Collapse
|
17
|
The Timing and Direction of Introgression Under the Multispecies Network Coalescent. Genetics 2019; 211:1059-1073. [PMID: 30670542 DOI: 10.1534/genetics.118.301831] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2018] [Accepted: 01/21/2019] [Indexed: 12/26/2022] Open
Abstract
Introgression is a pervasive biological process, and many statistical methods have been developed to infer its presence from genomic data. However, many of the consequences and genomic signatures of introgression remain unexplored from a methodological standpoint. Here, we develop a model for the timing and direction of introgression based on the multispecies network coalescent, and from it suggest new approaches for testing introgression hypotheses. We suggest two new statistics, D 1 and D 2, which can be used in conjunction with other information to test hypotheses relating to the timing and direction of introgression, respectively. D 1 may find use in evaluating cases of homoploid hybrid speciation (HHS), while D 2 provides a four-taxon test for polarizing introgression. Although analytical expectations for our statistics require a number of assumptions to be met, we show how simulations can be used to test hypotheses about introgression when these assumptions are violated. We apply the D 1 statistic to genomic data from the wild yeast Saccharomyces paradoxus-a proposed example of HHS-demonstrating its use as a test of this model. These methods provide new and powerful ways to address questions relating to the timing and direction of introgression.
Collapse
|
18
|
Beeravolu CR, Hickerson MJ, Frantz LAF, Lohse K. ABLE: blockwise site frequency spectra for inferring complex population histories and recombination. Genome Biol 2018; 19:145. [PMID: 30253810 PMCID: PMC6156964 DOI: 10.1186/s13059-018-1517-y] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2018] [Accepted: 08/22/2018] [Indexed: 01/08/2023] Open
Abstract
We introduce ABLE (Approximate Blockwise Likelihood Estimation), a novel simulation-based composite likelihood method that uses the blockwise site frequency spectrum to jointly infer past demography and recombination. ABLE is explicitly designed for a wide variety of data from unphased diploid genomes to genome-wide multi-locus data (for example, RADSeq) and can also accommodate arbitrarily large samples. We use simulations to demonstrate the accuracy of this method to infer complex histories of divergence and gene flow and reanalyze whole genome data from two species of orangutan. ABLE is available for download at https://github.com/champost/ABLE.
Collapse
Affiliation(s)
- Champak R Beeravolu
- Biology Department, The City College of New York, New York, 10031, NY, USA. .,Department of Evolutionary Biology and Environmental Studies, University of Zurich, Zurich, 8057, Switzerland.
| | - Michael J Hickerson
- Biology Department, The City College of New York, New York, 10031, NY, USA.,The Graduate Center, The City University of New York, New York, 10016, NY, USA.,Division of Invertebrate Zoology, American Museum of Natural History, New York, 10024, NY, USA
| | - Laurent A F Frantz
- Paleogenomics and Bio-Archaeology Research Network, Research Laboratory for Archeology and History of Art, University of Oxford, Oxford, OX1 3QY, UK.,School of Biological and Chemical Sciences, Queen Mary University of London, London, E1 4NS, UK
| | - Konrad Lohse
- Institute of Evolutionary Biology, University of Edinburgh, King's Buildings, Edinburgh, EH9 3FL, UK
| |
Collapse
|
19
|
Barlow A, Cahill JA, Hartmann S, Theunert C, Xenikoudakis G, Fortes GG, Paijmans JLA, Rabeder G, Frischauf C, Grandal-d'Anglade A, García-Vázquez A, Murtskhvaladze M, Saarma U, Anijalg P, Skrbinšek T, Bertorelle G, Gasparian B, Bar-Oz G, Pinhasi R, Slatkin M, Dalén L, Shapiro B, Hofreiter M. Partial genomic survival of cave bears in living brown bears. Nat Ecol Evol 2018; 2:1563-1570. [PMID: 30150744 DOI: 10.1038/s41559-018-0654-8] [Citation(s) in RCA: 64] [Impact Index Per Article: 10.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2017] [Accepted: 07/27/2018] [Indexed: 02/06/2023]
Abstract
Although many large mammal species went extinct at the end of the Pleistocene epoch, their DNA may persist due to past episodes of interspecies admixture. However, direct empirical evidence of the persistence of ancient alleles remains scarce. Here, we present multifold coverage genomic data from four Late Pleistocene cave bears (Ursus spelaeus complex) and show that cave bears hybridized with brown bears (Ursus arctos) during the Pleistocene. We develop an approach to assess both the directionality and relative timing of gene flow. We find that segments of cave bear DNA still persist in the genomes of living brown bears, with cave bears contributing 0.9 to 2.4% of the genomes of all brown bears investigated. Our results show that even though extinction is typically considered as absolute, following admixture, fragments of the gene pool of extinct species can survive for tens of thousands of years in the genomes of extant recipient species.
Collapse
Affiliation(s)
- Axel Barlow
- Institute for Biochemistry and Biology, University of Potsdam, Potsdam, Germany.
| | - James A Cahill
- Department of Ecology and Evolutionary Biology, University of California, Santa Cruz, Santa Cruz, CA, USA
| | - Stefanie Hartmann
- Institute for Biochemistry and Biology, University of Potsdam, Potsdam, Germany
| | - Christoph Theunert
- Department of Integrative Biology, University of California, Berkeley, Berkeley, CA, USA.,Department of Evolutionary Genetics, Max Planck Institute for Evolutionary Anthropology, Leipzig, Germany
| | | | - Gloria G Fortes
- Institute for Biochemistry and Biology, University of Potsdam, Potsdam, Germany.,Department of Life Sciences and Biotechnology, University of Ferrara, Ferrara, Italy
| | | | - Gernot Rabeder
- Institute of Palaeontology, University of Vienna, Vienna, Austria
| | | | | | - Ana García-Vázquez
- Instituto Universitario de Xeoloxía, Universidade da Coruña, A Coruña, Spain
| | | | - Urmas Saarma
- Department of Zoology, Institute of Ecology and Earth Sciences, University of Tartu, Tartu, Estonia
| | - Peeter Anijalg
- Department of Zoology, Institute of Ecology and Earth Sciences, University of Tartu, Tartu, Estonia
| | - Tomaž Skrbinšek
- Department of Biology, Biotechnical Faculty, University of Ljubljana, Ljubljana, Slovenia
| | - Giorgio Bertorelle
- Department of Life Sciences and Biotechnology, University of Ferrara, Ferrara, Italy
| | - Boris Gasparian
- Institute of Archaeology and Ethnography, National Academy of Sciences of the Republic of Armenia, Yerevan, Armenia
| | - Guy Bar-Oz
- Zinman Institute of Archaeology, University of Haifa, Haifa, Israel
| | - Ron Pinhasi
- Earth Institute, University College Dublin, Dublin, Ireland.,Department of Evolutionary Anthropology, University of Vienna, Vienna, Austria
| | - Montgomery Slatkin
- Department of Integrative Biology, University of California, Berkeley, Berkeley, CA, USA
| | - Love Dalén
- Department of Bioinformatics and Genetics, Swedish Museum of Natural History, Stockholm, Sweden
| | - Beth Shapiro
- Department of Ecology and Evolutionary Biology, University of California, Santa Cruz, Santa Cruz, CA, USA
| | - Michael Hofreiter
- Institute for Biochemistry and Biology, University of Potsdam, Potsdam, Germany
| |
Collapse
|
20
|
Povysil G, Hochreiter S. IBD Sharing between Africans, Neandertals, and Denisovans. Genome Biol Evol 2018; 8:3406-3416. [PMID: 28158547 PMCID: PMC5381509 DOI: 10.1093/gbe/evw234] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 09/19/2016] [Indexed: 12/03/2022] Open
Abstract
Interbreeding between ancestors of humans and other hominins outside of Africa has been studied intensively, while their common history within Africa still lacks proper attention. However, shedding light on human evolution in this time period about which little is known, is essential for understanding subsequent events outside of Africa. We investigate the genetic relationships of humans, Neandertals, and Denisovans by identifying very short DNA segments in the 1000 Genomes Phase 3 data that these hominins share identical by descent (IBD). By focusing on low frequency and rare variants, we identify very short IBD segments with high confidence. These segments reveal events from a very distant past because shorter IBD segments are presumably older than longer ones. We extracted two types of very old IBD segments that are not only shared among humans, but also with Neandertals and/or Denisovans. The first type contains longer segments that are found primarily in Asians and Europeans where more segments are found in South Asians than in East Asians for both Neandertal and Denisovan. These longer segments indicate complex admixture events outside of Africa. The second type consists of shorter segments that are shared mainly by Africans and therefore may indicate events involving ancestors of humans and other ancient hominins within Africa. Our results from the autosomes are further supported by an analysis of chromosome X, on which segments that are shared by Africans and match the Neandertal and/or Denisovan genome were even more prominent. Our results indicate that interbreeding with other hominins was a common feature of human evolution starting already long before ancestors of modern humans left Africa.
Collapse
Affiliation(s)
- Gundula Povysil
- Institute of Bioinformatics, Johannes Kepler University Linz, Austria
| | - Sepp Hochreiter
- Institute of Bioinformatics, Johannes Kepler University Linz, Austria
| |
Collapse
|
21
|
Nadeau NJ, Kawakami T. Population Genomics of Speciation and Admixture. POPULATION GENOMICS 2018. [DOI: 10.1007/13836_2018_24] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
|
22
|
Abstract
Homo sapiens phylogeography begins with the species' origin nearly 200 kya in Africa. First signs of the species outside Africa (in Arabia) are from 125 kya. Earliest dates elsewhere are now 100 kya in China, 45 kya in Australia and southern Europe (maybe even 60 kya in Australia), 32 kya in northeast Siberia, and maybe 20 kya in the Americas. Humans reached arctic regions and oceanic islands last-arctic North America about 5 kya, mid- and eastern Pacific islands about 2-1 kya, and New Zealand about 700 y ago. Initial routes along coasts seem the most likely given abundant and easily harvested shellfish there as indicated by huge ancient oyster shell middens on all continents. Nevertheless, the effect of geographic barriers-mountains and oceans-is clear. The phylogeographic pattern of diasporas from several single origins-northeast Africa to Eurasia, southeast Eurasia to Australia, and northeast Siberia to the Americas-allows the equivalent of a repeat experiment on the relation between geography and phylogenetic and cultural diversity. On all continents, cultural diversity is high in productive low latitudes, presumably because such regions can support populations of sustainable size in a small area, therefore allowing a high density of cultures. Of course, other factors operate. South America has an unusually low density of cultures in its tropical latitudes. A likely factor is the phylogeographic movement of peoples from the Old World bringing novel and hence, lethal diseases to the New World, a foretaste, perhaps, of present day global transport of tropical diseases.
Collapse
|
23
|
Stone GN, White SC, Csóka G, Melika G, Mutun S, Pénzes Z, Sadeghi SE, Schönrogge K, Tavakoli M, Nicholls JA. Tournament ABC analysis of the western Palaearctic population history of an oak gall wasp,Synergus umbraculus. Mol Ecol 2017; 26:6685-6703. [DOI: 10.1111/mec.14372] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2017] [Revised: 09/16/2017] [Accepted: 09/18/2017] [Indexed: 01/17/2023]
Affiliation(s)
- Graham N. Stone
- Institute of Evolutionary Biology; University of Edinburgh; Edinburgh UK
| | - Sarah C. White
- Institute of Evolutionary Biology; University of Edinburgh; Edinburgh UK
| | - György Csóka
- National Agricultural Research and Innovation Centre; Forest Research Institute; Mátrafüred Hungary
| | - George Melika
- Plant Health and Molecular Biology Laboratory; Directorate of Plant Protection, Soil Conservation and Agri-environment; Budapest Hungary
| | - Serap Mutun
- Department of Biology; Faculty of Science and Arts; Abant İzzet Baysal University; Bolu Turkey
| | - Zsolt Pénzes
- Department of Ecology; Faculty of Science and Informatics; University of Szeged; Szeged Hungary
| | - S. Ebrahim Sadeghi
- Agricultural Research, Education and Extension Organization (AREEO); Research Institute of Forests and Rangelands of Iran; Tehran Iran
| | | | - Majid Tavakoli
- Lorestan Agriculture and Natural Resources Research Center; Khorramabad Lorestan Iran
| | - James A. Nicholls
- Institute of Evolutionary Biology; University of Edinburgh; Edinburgh UK
| |
Collapse
|
24
|
Harris SE, Xue AT, Alvarado-Serrano D, Boehm JT, Joseph T, Hickerson MJ, Munshi-South J. Urbanization shapes the demographic history of a native rodent (the white-footed mouse, Peromyscus leucopus) in New York City. Biol Lett 2017; 12:rsbl.2015.0983. [PMID: 27072402 DOI: 10.1098/rsbl.2015.0983] [Citation(s) in RCA: 19] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2015] [Accepted: 03/14/2016] [Indexed: 01/25/2023] Open
Abstract
How urbanization shapes population genomic diversity and evolution of urban wildlife is largely unexplored. We investigated the impact of urbanization on white-footed mice,Peromyscus leucopus,in the New York City (NYC) metropolitan area using coalescent-based simulations to infer demographic history from the site-frequency spectrum. We assigned individuals to evolutionary clusters and then inferred recent divergence times, population size changes and migration using genome-wide single nucleotide polymorphisms genotyped in 23 populations sampled along an urban-to-rural gradient. Both prehistoric climatic events and recent urbanization impacted these populations. Our modelling indicates that post-glacial sea-level rise led to isolation of mainland and Long Island populations. These models also indicate that several urban parks represent recently isolated P. leucopus populations, and the estimated divergence times for these populations are consistent with the history of urbanization in NYC.
Collapse
Affiliation(s)
- Stephen E Harris
- The Graduate Center, City University of New York (CUNY), New York, NY 10016, USA
| | - Alexander T Xue
- The Graduate Center, City University of New York (CUNY), New York, NY 10016, USA The City College of New York, City University of New York, New York, NY 10031, USA
| | | | - Joel T Boehm
- The City College of New York, City University of New York, New York, NY 10031, USA
| | - Tyler Joseph
- The City College of New York, City University of New York, New York, NY 10031, USA
| | - Michael J Hickerson
- The Graduate Center, City University of New York (CUNY), New York, NY 10016, USA The City College of New York, City University of New York, New York, NY 10031, USA
| | - Jason Munshi-South
- Louis Calder Center-Biological Field Station, Fordham University, Armonk, NY 10504, USA
| |
Collapse
|
25
|
Brother or Other: The Place of Neanderthals in Human Evolution. VERTEBRATE PALEOBIOLOGY AND PALEOANTHROPOLOGY 2017. [DOI: 10.1007/978-3-319-46646-0_19] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|
26
|
Distribution of coalescent histories under the coalescent model with gene flow. Mol Phylogenet Evol 2016; 105:177-192. [DOI: 10.1016/j.ympev.2016.08.024] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/13/2016] [Revised: 08/16/2016] [Accepted: 08/31/2016] [Indexed: 12/19/2022]
|
27
|
Caldararo N. Denisovans, Melanesians, Europeans, and Neandertals: The Confusion of DNA Assumptions and the Biological Species Concept. J Mol Evol 2016; 83:78-87. [DOI: 10.1007/s00239-016-9755-7] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/02/2016] [Accepted: 07/30/2016] [Indexed: 10/21/2022]
|
28
|
|
29
|
Morozova I, Flegontov P, Mikheyev AS, Bruskin S, Asgharian H, Ponomarenko P, Klyuchnikov V, ArunKumar G, Prokhortchouk E, Gankin Y, Rogaev E, Nikolsky Y, Baranova A, Elhaik E, Tatarinova TV. Toward high-resolution population genomics using archaeological samples. DNA Res 2016; 23:295-310. [PMID: 27436340 PMCID: PMC4991838 DOI: 10.1093/dnares/dsw029] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/22/2015] [Accepted: 05/22/2016] [Indexed: 12/30/2022] Open
Abstract
The term ‘ancient DNA’ (aDNA) is coming of age, with over 1,200 hits in the PubMed database, beginning in the early 1980s with the studies of ‘molecular paleontology’. Rooted in cloning and limited sequencing of DNA from ancient remains during the pre-PCR era, the field has made incredible progress since the introduction of PCR and next-generation sequencing. Over the last decade, aDNA analysis ushered in a new era in genomics and became the method of choice for reconstructing the history of organisms, their biogeography, and migration routes, with applications in evolutionary biology, population genetics, archaeogenetics, paleo-epidemiology, and many other areas. This change was brought by development of new strategies for coping with the challenges in studying aDNA due to damage and fragmentation, scarce samples, significant historical gaps, and limited applicability of population genetics methods. In this review, we describe the state-of-the-art achievements in aDNA studies, with particular focus on human evolution and demographic history. We present the current experimental and theoretical procedures for handling and analysing highly degraded aDNA. We also review the challenges in the rapidly growing field of ancient epigenomics. Advancement of aDNA tools and methods signifies a new era in population genetics and evolutionary medicine research.
Collapse
Affiliation(s)
- Irina Morozova
- Institute of Evolutionary Medicine, University of Zurich, Zurich, Switzerland
| | - Pavel Flegontov
- Department of Biology and Ecology, Faculty of Science, University of Ostrava, Ostrava, Czech Republic Bioinformatics Center, A.A. Kharkevich Institute for Information Transmission Problems, Russian Academy of Sciences, Moscow, Russian Federation
| | - Alexander S Mikheyev
- Ecology and Evolution Unit, Okinawa Institute of Science and Technology Graduate University, Okinawa, Japan
| | - Sergey Bruskin
- Vavilov Institute of General Genetics RAS, Moscow, Russia
| | - Hosseinali Asgharian
- Department of Computational and Molecular Biology, University of Southern California, Los Angeles, CA, USA
| | - Petr Ponomarenko
- Center for Personalized Medicine, Children's Hospital Los Angeles, Los Angeles, CA, USA Spatial Sciences Institute, University of Southern California, Los Angeles, CA, USA
| | | | | | - Egor Prokhortchouk
- Research Center of Biotechnology RAS, Moscow, Russia Department of Biology, Lomonosov Moscow State University, Russia
| | | | - Evgeny Rogaev
- Vavilov Institute of General Genetics RAS, Moscow, Russia University of Massachusetts Medical School, Worcester, MA, USA
| | - Yuri Nikolsky
- Vavilov Institute of General Genetics RAS, Moscow, Russia F1 Genomics, San Diego, CA, USA School of Systems Biology, George Mason University, VA, USA
| | - Ancha Baranova
- School of Systems Biology, George Mason University, VA, USA Research Centre for Medical Genetics, Moscow, Russia Atlas Biomed Group, Moscow, Russia
| | - Eran Elhaik
- Department of Animal & Plant Sciences, University of Sheffield, Sheffield, South Yorkshire, UK
| | - Tatiana V Tatarinova
- Bioinformatics Center, A.A. Kharkevich Institute for Information Transmission Problems, Russian Academy of Sciences, Moscow, Russian Federation Center for Personalized Medicine, Children's Hospital Los Angeles, Los Angeles, CA, USA Spatial Sciences Institute, University of Southern California, Los Angeles, CA, USA
| |
Collapse
|
30
|
Payseur BA, Rieseberg LH. A genomic perspective on hybridization and speciation. Mol Ecol 2016; 25:2337-60. [PMID: 26836441 PMCID: PMC4915564 DOI: 10.1111/mec.13557] [Citation(s) in RCA: 292] [Impact Index Per Article: 36.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/12/2015] [Revised: 01/18/2016] [Accepted: 01/25/2016] [Indexed: 12/13/2022]
Abstract
Hybridization among diverging lineages is common in nature. Genomic data provide a special opportunity to characterize the history of hybridization and the genetic basis of speciation. We review existing methods and empirical studies to identify recent advances in the genomics of hybridization, as well as issues that need to be addressed. Notable progress has been made in the development of methods for detecting hybridization and inferring individual ancestries. However, few approaches reconstruct the magnitude and timing of gene flow, estimate the fitness of hybrids or incorporate knowledge of recombination rate. Empirical studies indicate that the genomic consequences of hybridization are complex, including a highly heterogeneous landscape of differentiation. Inferred characteristics of hybridization differ substantially among species groups. Loci showing unusual patterns - which may contribute to reproductive barriers - are usually scattered throughout the genome, with potential enrichment in sex chromosomes and regions of reduced recombination. We caution against the growing trend of interpreting genomic variation in summary statistics across genomes as evidence of differential gene flow. We argue that converting genomic patterns into useful inferences about hybridization will ultimately require models and methods that directly incorporate key ingredients of speciation, including the dynamic nature of gene flow, selection acting in hybrid populations and recombination rate variation.
Collapse
Affiliation(s)
- Bret A. Payseur
- Laboratory of Genetics, University of Wisconsin-Madison, Madison, WI 53706, USA
| | - Loren H. Rieseberg
- Department of Botany, University of British Columbia, Vancouver, BC V6T 1Z4, Canada
- Department of Biology, Indiana University, Bloomington, IN 47405, USA
| |
Collapse
|
31
|
Wen D, Yu Y, Hahn MW, Nakhleh L. Reticulate evolutionary history and extensive introgression in mosquito species revealed by phylogenetic network analysis. Mol Ecol 2016; 25:2361-72. [PMID: 26808290 DOI: 10.1111/mec.13544] [Citation(s) in RCA: 64] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2015] [Revised: 12/15/2015] [Accepted: 01/06/2016] [Indexed: 12/27/2022]
Abstract
The role of hybridization and subsequent introgression has been demonstrated in an increasing number of species. Recently, Fontaine et al. (Science, 347, 2015, 1258524) conducted a phylogenomic analysis of six members of the Anopheles gambiae species complex. Their analysis revealed a reticulate evolutionary history and pointed to extensive introgression on all four autosomal arms. The study further highlighted the complex evolutionary signals that the co-occurrence of incomplete lineage sorting (ILS) and introgression can give rise to in phylogenomic analyses. While tree-based methodologies were used in the study, phylogenetic networks provide a more natural model to capture reticulate evolutionary histories. In this work, we reanalyse the Anopheles data using a recently devised framework that combines the multispecies coalescent with phylogenetic networks. This framework allows us to capture ILS and introgression simultaneously, and forms the basis for statistical methods for inferring reticulate evolutionary histories. The new analysis reveals a phylogenetic network with multiple hybridization events, some of which differ from those reported in the original study. To elucidate the extent and patterns of introgression across the genome, we devise a new method that quantifies the use of reticulation branches in the phylogenetic network by each genomic region. Applying the method to the mosquito data set reveals the evolutionary history of all the chromosomes. This study highlights the utility of 'network thinking' and the new insights it can uncover, in particular in phylogenomic analyses of large data sets with extensive gene tree incongruence.
Collapse
Affiliation(s)
- Dingqiao Wen
- Department of Computer Science, Rice University, Houston, TX, 77005, USA
| | - Yun Yu
- Department of Computer Science, Rice University, Houston, TX, 77005, USA
| | - Matthew W Hahn
- Department of Biology, Indiana University, Bloomington, IN, 47405, USA.,School of Informatics and Computing, Indiana University, Bloomington, IN, 47405, USA
| | - Luay Nakhleh
- Department of Computer Science, Rice University, Houston, TX, 77005, USA.,Department of BioSciences, Rice University, Houston, TX, 77005, USA
| |
Collapse
|
32
|
Zhang W, Dasmahapatra KK, Mallet J, Moreira GRP, Kronforst MR. Genome-wide introgression among distantly related Heliconius butterfly species. Genome Biol 2016; 17:25. [PMID: 26921238 PMCID: PMC4769579 DOI: 10.1186/s13059-016-0889-0] [Citation(s) in RCA: 83] [Impact Index Per Article: 10.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2015] [Accepted: 01/28/2016] [Indexed: 12/30/2022] Open
Abstract
Background Although hybridization is thought to be relatively rare in animals, the raw genetic material introduced via introgression may play an important role in fueling adaptation and adaptive radiation. The butterfly genus Heliconius is an excellent system to study hybridization and introgression but most studies have focused on closely related species such as H. cydno and H. melpomene. Here we characterize genome-wide patterns of introgression between H. besckei, the only species with a red and yellow banded ‘postman’ wing pattern in the tiger-striped silvaniform clade, and co-mimetic H. melpomene nanna. Results We find a pronounced signature of putative introgression from H. melpomene into H. besckei in the genomic region upstream of the gene optix, known to control red wing patterning, suggesting adaptive introgression of wing pattern mimicry between these two distantly related species. At least 39 additional genomic regions show signals of introgression as strong or stronger than this mimicry locus. Gene flow has been on-going, with evidence of gene exchange at multiple time points, and bidirectional, moving from the melpomene to the silvaniform clade and vice versa. The history of gene exchange has also been complex, with contributions from multiple silvaniform species in addition to H. besckei. We also detect a signature of ancient introgression of the entire Z chromosome between the silvaniform and melpomene/cydno clades. Conclusions Our study provides a genome-wide portrait of introgression between distantly related butterfly species. We further propose a comprehensive and efficient workflow for gene flow identification in genomic data sets. Electronic supplementary material The online version of this article (doi:10.1186/s13059-016-0889-0) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Wei Zhang
- Department of Ecology & Evolution, University of Chicago, Chicago, IL, 60637, USA.
| | | | - James Mallet
- Department of Organismic & Evolutionary Biology, Harvard University, Cambridge, MA, 02138, USA
| | - Gilson R P Moreira
- Departamento de Zoologia, Universidade Federal do Rio Grande do Sul, Porto Alegre, RS, 91501-970, Brazil
| | - Marcus R Kronforst
- Department of Ecology & Evolution, University of Chicago, Chicago, IL, 60637, USA.
| |
Collapse
|
33
|
Frantz L, Meijaard E, Gongora J, Haile J, Groenen MA, Larson G. The Evolution of Suidae. Annu Rev Anim Biosci 2016; 4:61-85. [DOI: 10.1146/annurev-animal-021815-111155] [Citation(s) in RCA: 63] [Impact Index Per Article: 7.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Affiliation(s)
- Laurent Frantz
- Palaeogenomics & Bio-Archaeology Research Network, Research Laboratory for Archaeology and History of Art, University of Oxford, Oxford OX1 3QY, United Kingdom;
| | - Erik Meijaard
- IUCN/SSC Wild Pig Specialist Group, Jakarta 15412, Indonesia
- School of Archaeology and Anthropology, The Australian National University, Canberra, ACT 0200, Australia
| | - Jaime Gongora
- Faculty of Veterinary Science, The University of Sydney, Sydney, New South Wales 2006, Australia
| | - James Haile
- Palaeogenomics & Bio-Archaeology Research Network, Research Laboratory for Archaeology and History of Art, University of Oxford, Oxford OX1 3QY, United Kingdom;
| | - Martien A.M. Groenen
- Animal Breeding and Genomics Centre, Wageningen University, 6708 PB Wageningen, The Netherlands
| | - Greger Larson
- Palaeogenomics & Bio-Archaeology Research Network, Research Laboratory for Archaeology and History of Art, University of Oxford, Oxford OX1 3QY, United Kingdom;
| |
Collapse
|
34
|
Efficient Strategies for Calculating Blockwise Likelihoods Under the Coalescent. Genetics 2015; 202:775-86. [PMID: 26715666 PMCID: PMC4788249 DOI: 10.1534/genetics.115.183814] [Citation(s) in RCA: 51] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/15/2015] [Accepted: 12/15/2015] [Indexed: 01/08/2023] Open
Abstract
The inference of demographic history from genome data is hindered by a lack of efficient computational approaches. In particular, it has proved difficult to exploit the information contained in the distribution of genealogies across the genome. We have previously shown that the generating function (GF) of genealogies can be used to analytically compute likelihoods of demographic models from configurations of mutations in short sequence blocks (Lohse et al. 2011). Although the GF has a simple, recursive form, the size of such likelihood calculations explodes quickly with the number of individuals and applications of this framework have so far been mainly limited to small samples (pairs and triplets) for which the GF can be written by hand. Here we investigate several strategies for exploiting the inherent symmetries of the coalescent. In particular, we show that the GF of genealogies can be decomposed into a set of equivalence classes that allows likelihood calculations from nontrivial samples. Using this strategy, we automated blockwise likelihood calculations for a general set of demographic scenarios in Mathematica. These histories may involve population size changes, continuous migration, discrete divergence, and admixture between multiple populations. To give a concrete example, we calculate the likelihood for a model of isolation with migration (IM), assuming two diploid samples without phase and outgroup information. We demonstrate the new inference scheme with an analysis of two individual butterfly genomes from the sister species Heliconius melpomene rosina and H. cydno.
Collapse
|
35
|
Xue AT, Hickerson MJ. The aggregate site frequency spectrum for comparative population genomic inference. Mol Ecol 2015; 24:6223-40. [PMID: 26769405 PMCID: PMC4717917 DOI: 10.1111/mec.13447] [Citation(s) in RCA: 46] [Impact Index Per Article: 5.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/13/2015] [Revised: 10/26/2015] [Accepted: 10/28/2015] [Indexed: 12/11/2022]
Abstract
Understanding how assemblages of species responded to past climate change is a central goal of comparative phylogeography and comparative population genomics, an endeavour that has increasing potential to integrate with community ecology. New sequencing technology now provides the potential to perform complex demographic inference at unprecedented resolution across assemblages of nonmodel species. To this end, we introduce the aggregate site frequency spectrum (aSFS), an expansion of the site frequency spectrum to use single nucleotide polymorphism (SNP) data sets collected from multiple, co-distributed species for assemblage-level demographic inference. We describe how the aSFS is constructed over an arbitrary number of independent population samples and then demonstrate how the aSFS can differentiate various multispecies demographic histories under a wide range of sampling configurations while allowing effective population sizes and expansion magnitudes to vary independently. We subsequently couple the aSFS with a hierarchical approximate Bayesian computation (hABC) framework to estimate degree of temporal synchronicity in expansion times across taxa, including an empirical demonstration with a data set consisting of five populations of the threespine stickleback (Gasterosteus aculeatus). Corroborating what is generally understood about the recent postglacial origins of these populations, the joint aSFS/hABC analysis strongly suggests that the stickleback data are most consistent with synchronous expansion after the Last Glacial Maximum (posterior probability = 0.99). The aSFS will have general application for multilevel statistical frameworks to test models involving assemblages and/or communities, and as large-scale SNP data from nonmodel species become routine, the aSFS expands the potential for powerful next-generation comparative population genomic inference.
Collapse
Affiliation(s)
- Alexander T. Xue
- Department of Biology: Subprogram in Ecology, Evolutionary Biology, and Behavior, City College and Graduate Center of City University of New York, 160 Convent Avenue, Marshak Science Building, Room 526, New York, NY 10031
| | - Michael J. Hickerson
- Department of Biology: Subprogram in Ecology, Evolutionary Biology, and Behavior, City College and Graduate Center of City University of New York, 160 Convent Avenue, Marshak Science Building, Room 526, New York, NY 10031
| |
Collapse
|
36
|
|
37
|
Demos TC, Kerbis Peterhans JC, Joseph TA, Robinson JD, Agwanda B, Hickerson MJ. Comparative Population Genomics of African Montane Forest Mammals Support Population Persistence across a Climatic Gradient and Quaternary Climatic Cycles. PLoS One 2015; 10:e0131800. [PMID: 26394036 PMCID: PMC4578943 DOI: 10.1371/journal.pone.0131800] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2014] [Accepted: 06/07/2015] [Indexed: 12/03/2022] Open
Abstract
The Eastern Afromontane biodiversity hotspot (EABH) has the highest concentration of biodiversity in tropical Africa, yet few studies have investigated recent historical diversification processes in EABH lineages. Herein, we analyze restriction-site associated DNA-sequences (RAD-Seq) to study recent historical processes in co-distributed mouse (Hylomyscus) and shrew (Sylvisorex) species complexes, with an aim to better determine how historical paleoenvironmental processes might have contributed to the EABH’s high diversity. We analyzed complete SNP matrices of > 50,000 RAD loci to delineate populations, reconstruct the history of isolation and admixture, and discover geographic patterns of genetic partitioning. These analyses demonstrate that persistently unsuitable habitat may have isolated multiple populations distributed across montane habitat islands in the Itombwe Massif and Albertine Rift to the west as well as Mt Elgon and Kenyan Highlands to the east. We detected low genetic diversity in Kenyan Highland populations of both genera, consistent with smaller historical population sizes in this region. We additionally tested predictions that Albertine Rift populations are older and more persistently isolated compared to the Kenyan Highlands. Phylogenetic analyses support greater historical isolation among Albertine Rift populations of both shrews and mice compared to the Kenyan Highlands and suggest that there are genetically isolated populations from both focal genera in the Itombwe Massif, Democratic Republic of Congo. The Albertine Rift ecoregion has the highest mammalian tropical forest species richness per unit area on earth. Our results clearly support accelerating efforts to conserve this diversity.
Collapse
Affiliation(s)
- Terrence C. Demos
- Biology Department, City College of New York, City University of New York, New York, New York, 10031, United States of America
- Biology Doctoral Program, The Graduate School and University Center, City University of New York, New York, New York, 10016, United States of America
- Integrative Research Center, Field Museum of Natural History, Chicago, Illinois, 60605, United States of America
- * E-mail:
| | - Julian C. Kerbis Peterhans
- Integrative Research Center, Field Museum of Natural History, Chicago, Illinois, 60605, United States of America
- College of Professional Studies, Roosevelt University, Chicago, Illinois, 60605, United States of America
| | - Tyler A. Joseph
- Biology Department, City College of New York, City University of New York, New York, New York, 10031, United States of America
| | - John D. Robinson
- Biology Department, City College of New York, City University of New York, New York, New York, 10031, United States of America
- South Carolina Department of Natural Resources, Charleston, South Carolina, 29412, United States of America
| | - Bernard Agwanda
- Mammalogy Section, National Museums of Kenya, Nairobi, 00100, Kenya
| | - Michael J. Hickerson
- Biology Department, City College of New York, City University of New York, New York, New York, 10031, United States of America
- Biology Doctoral Program, The Graduate School and University Center, City University of New York, New York, New York, 10016, United States of America
- Division of Invertebrate Zoology, American Museum of Natural History, New York, New York 10024, United States of America
| |
Collapse
|
38
|
Nicholls JA, Pennington RT, Koenen EJM, Hughes CE, Hearn J, Bunnefeld L, Dexter KG, Stone GN, Kidner CA. Using targeted enrichment of nuclear genes to increase phylogenetic resolution in the neotropical rain forest genus Inga (Leguminosae: Mimosoideae). FRONTIERS IN PLANT SCIENCE 2015; 6:710. [PMID: 26442024 PMCID: PMC4584976 DOI: 10.3389/fpls.2015.00710] [Citation(s) in RCA: 53] [Impact Index Per Article: 5.9] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/22/2015] [Accepted: 08/25/2015] [Indexed: 05/20/2023]
Abstract
Evolutionary radiations are prominent and pervasive across many plant lineages in diverse geographical and ecological settings; in neotropical rainforests there is growing evidence suggesting that a significant fraction of species richness is the result of recent radiations. Understanding the evolutionary trajectories and mechanisms underlying these radiations demands much greater phylogenetic resolution than is currently available for these groups. The neotropical tree genus Inga (Leguminosae) is a good example, with ~300 extant species and a crown age of 2-10 MY, yet over 6 kb of plastid and nuclear DNA sequence data gives only poor phylogenetic resolution among species. Here we explore the use of larger-scale nuclear gene data obtained though targeted enrichment to increase phylogenetic resolution within Inga. Transcriptome data from three Inga species were used to select 264 nuclear loci for targeted enrichment and sequencing. Following quality control to remove probable paralogs from these sequence data, the final dataset comprised 259,313 bases from 194 loci for 24 accessions representing 22 Inga species and an outgroup (Zygia). Bayesian phylogenies reconstructed using either all loci concatenated or a gene-tree/species-tree approach yielded highly resolved phylogenies. We used coalescent approaches to show that the same targeted enrichment data also have significant power to discriminate among alternative within-species population histories within the widespread species I. umbellifera. In either application, targeted enrichment simplifies the informatics challenge of identifying orthologous loci associated with de novo genome sequencing. We conclude that targeted enrichment provides the large volumes of phylogenetically-informative sequence data required to resolve relationships within recent plant species radiations, both at the species level and for within-species phylogeographic studies.
Collapse
Affiliation(s)
- James A. Nicholls
- Ashworth Labs, Institute of Evolutionary Biology, School of Biological Sciences, University of EdinburghEdinburgh, UK
- Royal Botanic Garden EdinburghEdinburgh, UK
| | | | - Erik J. M. Koenen
- Institute of Systematic Botany, University of ZurichZürich, Switzerland
| | - Colin E. Hughes
- Institute of Systematic Botany, University of ZurichZürich, Switzerland
| | - Jack Hearn
- Ashworth Labs, Institute of Evolutionary Biology, School of Biological Sciences, University of EdinburghEdinburgh, UK
| | - Lynsey Bunnefeld
- Ashworth Labs, Institute of Evolutionary Biology, School of Biological Sciences, University of EdinburghEdinburgh, UK
| | - Kyle G. Dexter
- School of Geosciences, University of EdinburghEdinburgh, UK
| | - Graham N. Stone
- Ashworth Labs, Institute of Evolutionary Biology, School of Biological Sciences, University of EdinburghEdinburgh, UK
| | - Catherine A. Kidner
- Royal Botanic Garden EdinburghEdinburgh, UK
- Institute of Molecular Plant Sciences, School of Biological Sciences, University of EdinburghEdinburgh, UK
| |
Collapse
|
39
|
Inferring Bottlenecks from Genome-Wide Samples of Short Sequence Blocks. Genetics 2015; 201:1157-69. [PMID: 26341659 PMCID: PMC4649642 DOI: 10.1534/genetics.115.179861] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2015] [Accepted: 09/01/2015] [Indexed: 01/02/2023] Open
Abstract
The advent of the genomic era has necessitated the development of methods capable of analyzing large volumes of genomic data efficiently. Being able to reliably identify bottlenecks—extreme population size changes of short duration—not only is interesting in the context of speciation and extinction but also matters (as a null model) when inferring selection. Bottlenecks can be detected in polymorphism data via their distorting effect on the shape of the underlying genealogy. Here, we use the generating function of genealogies to derive the probability of mutational configurations in short sequence blocks under a simple bottleneck model. Given a large number of nonrecombining blocks, we can compute maximum-likelihood estimates of the time and strength of the bottleneck. Our method relies on a simple summary of the joint distribution of polymorphic sites. We extend the site frequency spectrum by counting mutations in frequency classes in short sequence blocks. Using linkage information over short distances in this way gives greater power to detect bottlenecks than the site frequency spectrum and potentially opens up a wide range of demographic histories to blockwise inference. Finally, we apply our method to genomic data from a species of pig (Sus cebifrons) endemic to islands in the center and west of the Philippines to estimate whether a bottleneck occurred upon island colonization and compare our scheme to Li and Durbin’s pairwise sequentially Markovian coalescent (PSMC) both for the pig data and using simulations.
Collapse
|
40
|
Gash DM, Deane AS. Neuron-based heredity and human evolution. Front Neurosci 2015; 9:209. [PMID: 26136649 PMCID: PMC4469835 DOI: 10.3389/fnins.2015.00209] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2015] [Accepted: 05/26/2015] [Indexed: 11/22/2022] Open
Abstract
It is widely recognized that human evolution has been driven by two systems of heredity: one DNA-based and the other based on the transmission of behaviorally acquired information via nervous system functions. The genetic system is ancient, going back to the appearance of life on Earth. It is responsible for the evolutionary processes described by Darwin. By comparison, the nervous system is relatively newly minted and in its highest form, responsible for ideation and mind-to-mind transmission of information. Here the informational capabilities and functions of the two systems are compared. While employing quite different mechanisms for encoding, storing and transmission of information, both systems perform these generic hereditary functions. Three additional features of neuron-based heredity in humans are identified: the ability to transfer hereditary information to other members of their population, not just progeny; a selection process for the information being transferred; and a profoundly shorter time span for creation and dissemination of survival-enhancing information in a population. The mechanisms underlying neuron-based heredity involve hippocampal neurogenesis and memory and learning processes modifying and creating new neural assemblages changing brain structure and functions. A fundamental process in rewiring brain circuitry is through increased neural activity (use) strengthening and increasing the number of synaptic connections. Decreased activity in circuitry (disuse) leads to loss of synapses. Use and disuse modifying an organ to bring about new modes of living, habits and functions are processes in line with Neolamarckian concepts of evolution (Packard, 1901). Evidence is presented of bipartite evolutionary processes-Darwinian and Neolamarckian-driving human descent from a common ancestor shared with the great apes.
Collapse
Affiliation(s)
- Don M. Gash
- Department of Anatomy and Neurobiology, College of Medicine, University of KentuckyLexington, KY, USA
| | | |
Collapse
|
41
|
Racimo F, Sankararaman S, Nielsen R, Huerta-Sánchez E. Evidence for archaic adaptive introgression in humans. Nat Rev Genet 2015; 16:359-71. [PMID: 25963373 PMCID: PMC4478293 DOI: 10.1038/nrg3936] [Citation(s) in RCA: 314] [Impact Index Per Article: 34.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
As modern and ancient DNA sequence data from diverse human populations accumulate, evidence is increasing in support of the existence of beneficial variants acquired from archaic humans that may have accelerated adaptation and improved survival in new environments - a process known as adaptive introgression. Within the past few years, a series of studies have identified genomic regions that show strong evidence for archaic adaptive introgression. Here, we provide an overview of the statistical methods developed to identify archaic introgressed fragments in the genome sequences of modern humans and to determine whether positive selection has acted on these fragments. We review recently reported examples of adaptive introgression, grouped by selection pressure, and consider the level of supporting evidence for each. Finally, we discuss challenges and recommendations for inferring selection on introgressed regions.
Collapse
Affiliation(s)
- Fernando Racimo
- Department of Integrative Biology, UC Berkeley, Berkeley CA 97420
| | | | - Rasmus Nielsen
- Department of Integrative Biology, UC Berkeley, Berkeley CA 97420
- Department of Statistics, UC Berkeley, Berkeley CA 97420
| | | |
Collapse
|
42
|
Paudel Y, Madsen O, Megens HJ, Frantz LAF, Bosse M, Crooijmans RPMA, Groenen MAM. Copy number variation in the speciation of pigs: a possible prominent role for olfactory receptors. BMC Genomics 2015; 16:330. [PMID: 25896665 PMCID: PMC4413995 DOI: 10.1186/s12864-015-1449-9] [Citation(s) in RCA: 66] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/14/2014] [Accepted: 03/09/2015] [Indexed: 12/02/2022] Open
Abstract
Background Unraveling the genetic mechanisms associated with reduced gene flow between genetically differentiated populations is key to understand speciation. Different types of structural variations (SVs) have been found as a source of genetic diversity in a wide range of species. Previous studies provided detailed knowledge on the potential evolutionary role of SVs, especially copy number variations (CNVs), between well diverged species of e.g. primates. However, our understanding of their significance during ongoing speciation processes is limited due to the lack of CNV data from closely related species. The genus Sus (pig and its close relatives) which started to diverge ~4 Mya presents an excellent model for studying the role of CNVs during ongoing speciation. Results In this study, we identified 1408 CNV regions (CNVRs) across the genus Sus. These CNVRs encompass 624 genes and were found to evolve ~2.5 times faster than single nucleotide polymorphisms (SNPs). The majority of these copy number variable genes are olfactory receptors (ORs) known to play a prominent role in food foraging and mate recognition in Sus. Phylogenetic analyses, including novel Bayesian analysis, based on CNVRs that overlap ORs retain the well-accepted topology of the genus Sus whereas CNVRs overlapping genes other than ORs show evidence for random drift and/or admixture. Conclusion We hypothesize that inter-specific variation in copy number of ORs provided the means for rapid adaptation to different environments during the diversification of the genus Sus in the Pliocene. Furthermore, these regions might have acted as barriers preventing massive gene flow between these species during the multiple hybridization events that took place later in the Pleistocene suggesting a possible prominent role of ORs in the ongoing Sus speciation. Electronic supplementary material The online version of this article (doi:10.1186/s12864-015-1449-9) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Yogesh Paudel
- Animal Breeding and Genomics Centre, Wageningen University, 6700 AH, Wageningen, The Netherlands. .,Current address: Roche Pharma Research and Early Development, Roche Innovation Center Basel, F. Hoffmann-La Roche Ltd, 4070, Basel, Switzerland.
| | - Ole Madsen
- Animal Breeding and Genomics Centre, Wageningen University, 6700 AH, Wageningen, The Netherlands.
| | - Hendrik-Jan Megens
- Animal Breeding and Genomics Centre, Wageningen University, 6700 AH, Wageningen, The Netherlands.
| | - Laurent A F Frantz
- Animal Breeding and Genomics Centre, Wageningen University, 6700 AH, Wageningen, The Netherlands.
| | - Mirte Bosse
- Animal Breeding and Genomics Centre, Wageningen University, 6700 AH, Wageningen, The Netherlands.
| | - Richard P M A Crooijmans
- Animal Breeding and Genomics Centre, Wageningen University, 6700 AH, Wageningen, The Netherlands.
| | - Martien A M Groenen
- Animal Breeding and Genomics Centre, Wageningen University, 6700 AH, Wageningen, The Netherlands.
| |
Collapse
|
43
|
Pease JB, Hahn MW. Detection and Polarization of Introgression in a Five-Taxon Phylogeny. Syst Biol 2015; 64:651-62. [PMID: 25888025 DOI: 10.1093/sysbio/syv023] [Citation(s) in RCA: 151] [Impact Index Per Article: 16.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2015] [Accepted: 04/10/2015] [Indexed: 12/22/2022] Open
Abstract
When multiple speciation events occur rapidly in succession, discordant genealogies due to incomplete lineage sorting (ILS) can complicate the detection of introgression. A variety of methods, including the [Formula: see text]-statistic (a.k.a. the "ABBA-BABA test"), have been proposed to infer introgression in the presence of ILS for a four-taxon clade. However, no integrated method exists to detect introgression using allelic patterns for more complex phylogenies. Here we explore the issues associated with previous systems of applying [Formula: see text]-statistics to a larger tree topology, and propose new [Formula: see text] tests as an integrated framework to infer both the taxa involved in and the direction of introgression for a symmetric five-taxon phylogeny. Using theory and simulations, we show that the [Formula: see text] statistics correctly identify the introgression donor and recipient lineages, even at low rates of introgression. [Formula: see text] is also shown to have extremely low false-positive rates. The [Formula: see text] tests are computationally inexpensive to calculate and can easily be applied to phylogenomic data sets, both genome-wide and in windows of the genome. In addition, we explore both the principles and problems of introgression detection in even more complex phylogenies.
Collapse
Affiliation(s)
- James B Pease
- Department of Biology, Indiana University, 1001 E. Third St., Bloomington, IN 47405, USA and School of Informatics and Computing, Indiana University, Bloomington, IN 47405, USA
| | - Matthew W Hahn
- Department of Biology, Indiana University, 1001 E. Third St., Bloomington, IN 47405, USA and School of Informatics and Computing, Indiana University, Bloomington, IN 47405, USA Department of Biology, Indiana University, 1001 E. Third St., Bloomington, IN 47405, USA and School of Informatics and Computing, Indiana University, Bloomington, IN 47405, USA
| |
Collapse
|
44
|
Nicholls JA, Pennington RT, Koenen EJM, Hughes CE, Hearn J, Bunnefeld L, Dexter KG, Stone GN, Kidner CA. Using targeted enrichment of nuclear genes to increase phylogenetic resolution in the neotropical rain forest genus Inga (Leguminosae: Mimosoideae). FRONTIERS IN PLANT SCIENCE 2015. [PMID: 26442024 DOI: 10.5061/dryad.r9c12] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Subscribe] [Scholar Register] [Indexed: 05/02/2023]
Abstract
Evolutionary radiations are prominent and pervasive across many plant lineages in diverse geographical and ecological settings; in neotropical rainforests there is growing evidence suggesting that a significant fraction of species richness is the result of recent radiations. Understanding the evolutionary trajectories and mechanisms underlying these radiations demands much greater phylogenetic resolution than is currently available for these groups. The neotropical tree genus Inga (Leguminosae) is a good example, with ~300 extant species and a crown age of 2-10 MY, yet over 6 kb of plastid and nuclear DNA sequence data gives only poor phylogenetic resolution among species. Here we explore the use of larger-scale nuclear gene data obtained though targeted enrichment to increase phylogenetic resolution within Inga. Transcriptome data from three Inga species were used to select 264 nuclear loci for targeted enrichment and sequencing. Following quality control to remove probable paralogs from these sequence data, the final dataset comprised 259,313 bases from 194 loci for 24 accessions representing 22 Inga species and an outgroup (Zygia). Bayesian phylogenies reconstructed using either all loci concatenated or a gene-tree/species-tree approach yielded highly resolved phylogenies. We used coalescent approaches to show that the same targeted enrichment data also have significant power to discriminate among alternative within-species population histories within the widespread species I. umbellifera. In either application, targeted enrichment simplifies the informatics challenge of identifying orthologous loci associated with de novo genome sequencing. We conclude that targeted enrichment provides the large volumes of phylogenetically-informative sequence data required to resolve relationships within recent plant species radiations, both at the species level and for within-species phylogeographic studies.
Collapse
Affiliation(s)
- James A Nicholls
- Ashworth Labs, Institute of Evolutionary Biology, School of Biological Sciences, University of Edinburgh Edinburgh, UK ; Royal Botanic Garden Edinburgh Edinburgh, UK
| | | | - Erik J M Koenen
- Institute of Systematic Botany, University of Zurich Zürich, Switzerland
| | - Colin E Hughes
- Institute of Systematic Botany, University of Zurich Zürich, Switzerland
| | - Jack Hearn
- Ashworth Labs, Institute of Evolutionary Biology, School of Biological Sciences, University of Edinburgh Edinburgh, UK
| | - Lynsey Bunnefeld
- Ashworth Labs, Institute of Evolutionary Biology, School of Biological Sciences, University of Edinburgh Edinburgh, UK
| | - Kyle G Dexter
- School of Geosciences, University of Edinburgh Edinburgh, UK
| | - Graham N Stone
- Ashworth Labs, Institute of Evolutionary Biology, School of Biological Sciences, University of Edinburgh Edinburgh, UK
| | - Catherine A Kidner
- Royal Botanic Garden Edinburgh Edinburgh, UK ; Institute of Molecular Plant Sciences, School of Biological Sciences, University of Edinburgh Edinburgh, UK
| |
Collapse
|
45
|
Frantz LAF, Madsen O, Megens HJ, Groenen MAM, Lohse K. Testing models of speciation from genome sequences: divergence and asymmetric admixture in Island South-East Asian Sus species during the Plio-Pleistocene climatic fluctuations. Mol Ecol 2014; 23:5566-74. [PMID: 25294645 PMCID: PMC4245187 DOI: 10.1111/mec.12958] [Citation(s) in RCA: 25] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2014] [Revised: 10/01/2014] [Accepted: 10/03/2014] [Indexed: 12/19/2022]
Abstract
In many temperate regions, ice ages promoted range contractions into refugia resulting in divergence (and potentially speciation), while warmer periods led to range expansions and hybridization. However, the impact these climatic oscillations had in many parts of the tropics remains elusive. Here, we investigate this issue using genome sequences of three pig (Sus) species, two of which are found on islands of the Sunda-shelf shallow seas in Island South-East Asia (ISEA). A previous study revealed signatures of interspecific admixture between these Sus species (Genome biology, 14, 2013, R107). However, the timing, directionality and extent of this admixture remain unknown. Here, we use a likelihood-based model comparison to more finely resolve this admixture history and test whether it was mediated by humans or occurred naturally. Our analyses suggest that interspecific admixture between Sunda-shelf species was most likely asymmetric and occurred long before the arrival of humans in the region. More precisely, we show that these species diverged during the late Pliocene but around 23% of their genomes have been affected by admixture during the later Pleistocene climatic transition. In addition, we show that our method provides a significant improvement over D-statistics which are uninformative about the direction of admixture.
Collapse
Affiliation(s)
- Laurent A F Frantz
- Animal Breeding and Genomics Group, Wageningen University, De Elst 1, Wageningen, WD, 6708, The Netherlands
| | | | | | | | | |
Collapse
|
46
|
Robinson JD, Bunnefeld L, Hearn J, Stone GN, Hickerson MJ. ABC inference of multi-population divergence with admixture from unphased population genomic data. Mol Ecol 2014; 23:4458-71. [PMID: 25113024 PMCID: PMC4285295 DOI: 10.1111/mec.12881] [Citation(s) in RCA: 48] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/02/2013] [Revised: 08/04/2014] [Accepted: 08/06/2014] [Indexed: 01/13/2023]
Abstract
Rapidly developing sequencing technologies and declining costs have made it possible to collect genome-scale data from population-level samples in nonmodel systems. Inferential tools for historical demography given these data sets are, at present, underdeveloped. In particular, approximate Bayesian computation (ABC) has yet to be widely embraced by researchers generating these data. Here, we demonstrate the promise of ABC for analysis of the large data sets that are now attainable from nonmodel taxa through current genomic sequencing technologies. We develop and test an ABC framework for model selection and parameter estimation, given histories of three-population divergence with admixture. We then explore different sampling regimes to illustrate how sampling more loci, longer loci or more individuals affects the quality of model selection and parameter estimation in this ABC framework. Our results show that inferences improved substantially with increases in the number and/or length of sequenced loci, while less benefit was gained by sampling large numbers of individuals. Optimal sampling strategies given our inferential models included at least 2000 loci, each approximately 2 kb in length, sampled from five diploid individuals per population, although specific strategies are model and question dependent. We tested our ABC approach through simulation-based cross-validations and illustrate its application using previously analysed data from the oak gall wasp, Biorhiza pallida.
Collapse
Affiliation(s)
- John D Robinson
- Department of Biology, City College of New York, 160 Convent Ave., MR 526, New York, NY, 10031, USA
| | | | | | | | | |
Collapse
|
47
|
Johansson S. The thinking Neanderthals: What do we know about Neanderthal cognition? WILEY INTERDISCIPLINARY REVIEWS. COGNITIVE SCIENCE 2014; 5:613-620. [PMID: 26308868 DOI: 10.1002/wcs.1317] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/05/2014] [Revised: 07/28/2014] [Accepted: 08/04/2014] [Indexed: 01/12/2023]
Abstract
The study of Neanderthal cognition is difficult, because of the archaeological invisibility of cognition, and because of the methodological issues that arise both from that invisibility and from their being close to modern humans. Nevertheless, fair progress has been made in gathering relevant evidence. There is now good evidence that Neanderthals were cognitively sophisticated, displaying many of the cognitive traits that were traditionally regarded as proxies for modern human cognition, notably including language. It can neither be proven nor excluded that they were our cognitive equals, but they were close enough to us, biologically and cognitively, to interbreed successfully and leave a genetic legacy in our DNA. WIREs Cogn Sci 2014, 5:613-620. doi: 10.1002/wcs.1317 CONFLICT OF INTEREST: The author has declared no conflicts of interest for this article. For further resources related to this article, please visit the WIREs website.
Collapse
|
48
|
Higham T, Douka K, Wood R, Ramsey CB, Brock F, Basell L, Camps M, Arrizabalaga A, Baena J, Barroso-Ruíz C, Bergman C, Boitard C, Boscato P, Caparrós M, Conard NJ, Draily C, Froment A, Galván B, Gambassini P, Garcia-Moreno A, Grimaldi S, Haesaerts P, Holt B, Iriarte-Chiapusso MJ, Jelinek A, Jordá Pardo JF, Maíllo-Fernández JM, Marom A, Maroto J, Menéndez M, Metz L, Morin E, Moroni A, Negrino F, Panagopoulou E, Peresani M, Pirson S, de la Rasilla M, Riel-Salvatore J, Ronchitelli A, Santamaria D, Semal P, Slimak L, Soler J, Soler N, Villaluenga A, Pinhasi R, Jacobi R. The timing and spatiotemporal patterning of Neanderthal disappearance. Nature 2014; 512:306-9. [DOI: 10.1038/nature13621] [Citation(s) in RCA: 561] [Impact Index Per Article: 56.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2014] [Accepted: 06/27/2014] [Indexed: 11/09/2022]
|
49
|
Shriner D, Tekola-Ayele F, Adeyemo A, Rotimi CN. Genome-wide genotype and sequence-based reconstruction of the 140,000 year history of modern human ancestry. Sci Rep 2014; 4:6055. [PMID: 25116736 PMCID: PMC4131216 DOI: 10.1038/srep06055] [Citation(s) in RCA: 44] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/05/2014] [Accepted: 07/28/2014] [Indexed: 12/11/2022] Open
Abstract
We investigated ancestry of 3,528 modern humans from 163 samples. We identified 19 ancestral components, with 94.4% of individuals showing mixed ancestry. After using whole genome sequences to correct for ascertainment biases in genome-wide genotype data, we dated the oldest divergence event to 140,000 years ago. We detected an Out-of-Africa migration 100,000–87,000 years ago, leading to peoples of the Americas, east and north Asia, and Oceania, followed by another migration 61,000–44,000 years ago, leading to peoples of the Caucasus, Europe, the Middle East, and south Asia. We dated eight divergence events to 33,000–20,000 years ago, coincident with the Last Glacial Maximum. We refined understanding of the ancestry of several ethno-linguistic groups, including African Americans, Ethiopians, the Kalash, Latin Americans, Mozabites, Pygmies, and Uygurs, as well as the CEU sample. Ubiquity of mixed ancestry emphasizes the importance of accounting for ancestry in history, forensics, and health.
Collapse
Affiliation(s)
- Daniel Shriner
- Center for Research on Genomics and Global Health, National Human Genome Research Institute, Building 12A, Room 4047, 12 South Drive, Bethesda, Maryland 20892 USA
| | - Fasil Tekola-Ayele
- Center for Research on Genomics and Global Health, National Human Genome Research Institute, Building 12A, Room 4047, 12 South Drive, Bethesda, Maryland 20892 USA
| | - Adebowale Adeyemo
- Center for Research on Genomics and Global Health, National Human Genome Research Institute, Building 12A, Room 4047, 12 South Drive, Bethesda, Maryland 20892 USA
| | - Charles N Rotimi
- Center for Research on Genomics and Global Health, National Human Genome Research Institute, Building 12A, Room 4047, 12 South Drive, Bethesda, Maryland 20892 USA
| |
Collapse
|