1
|
Li L, Comi TJ, Bierman RF, Akey JM. Recurrent gene flow between Neanderthals and modern humans over the past 200,000 years. Science 2024; 385:eadi1768. [PMID: 38991054 DOI: 10.1126/science.adi1768] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/05/2023] [Accepted: 05/14/2024] [Indexed: 07/13/2024]
Abstract
Although it is well known that the ancestors of modern humans and Neanderthals admixed, the effects of gene flow on the Neanderthal genome are not well understood. We develop methods to estimate the amount of human-introgressed sequences in Neanderthals and apply it to whole-genome sequence data from 2000 modern humans and three Neanderthals. We estimate that Neanderthals have 2.5 to 3.7% human ancestry, and we leverage human-introgressed sequences in Neanderthals to revise estimates of Neanderthal ancestry in modern humans, show that Neanderthal population sizes were significantly smaller than previously estimated, and identify two distinct waves of modern human gene flow into Neanderthals. Our data provide insights into the genetic legacy of recurrent gene flow between modern humans and Neanderthals.
Collapse
Affiliation(s)
- Liming Li
- Department of Medical Genetics and Developmental Biology, School of Medicine, The Key Laboratory of Developmental Genes and Human Diseases, Ministry of Education, Southeast University, Nanjing 210009, China
- The Lewis-Sigler Institute for Integrative Genomics, Princeton University, Princeton, NJ 08540, USA
| | - Troy J Comi
- The Lewis-Sigler Institute for Integrative Genomics, Princeton University, Princeton, NJ 08540, USA
| | - Rob F Bierman
- The Lewis-Sigler Institute for Integrative Genomics, Princeton University, Princeton, NJ 08540, USA
| | - Joshua M Akey
- The Lewis-Sigler Institute for Integrative Genomics, Princeton University, Princeton, NJ 08540, USA
| |
Collapse
|
2
|
Ma X, Lu Y, Xu S. Adaptive Evolution of Two Distinct Adaptive Haplotypes of Neanderthal Origin at the Immunoglobulin Heavy-chain Locus in East Asian and European Populations. Mol Biol Evol 2024; 41:msae147. [PMID: 39011558 DOI: 10.1093/molbev/msae147] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2023] [Revised: 07/05/2024] [Accepted: 07/11/2024] [Indexed: 07/17/2024] Open
Abstract
Immunoglobulins (Igs) have a crucial role in humoral immunity. Two recent studies have reported a high-frequency Neanderthal-introgressed haplotype throughout Eurasia and a high-frequency Neanderthal-introgressed haplotype specific to southern East Asia at the immunoglobulin heavy-chain (IGH) gene locus on chromosome 14q32.33. Surprisingly, we found the previously reported high-frequency Neanderthal-introgressed haplotype does not exist throughout Eurasia. Instead, our study identified two distinct high-frequency haplotypes of putative Neanderthal origin in East Asia and Europe, although they shared introgressed alleles. Notably, the alleles of putative Neanderthal origin reduced the expression of IGHG1 and increased the expression of IGHG2 and IGHG3 in various tissues. These putatively introgressed alleles also affected the production of IgG1 upon antigen stimulation and increased the risk of systemic lupus erythematosus. Additionally, the greatest genetic differentiation across the whole genome between southern and northern East Asians was observed for the East Asian haplotype of putative Neanderthal origin. The frequency decreased from southern to northern East Asia and correlated positively with the genome-wide proportion of southern East Asian ancestry, indicating that this putative positive selection likely occurred in the common ancestor of southern East Asian populations before the admixture with northern East Asian populations.
Collapse
Affiliation(s)
- Xixian Ma
- Key Laboratory of Computational Biology, Shanghai Institute of Nutrition and Health, University of Chinese Academy of Sciences, Chinese Academy of Sciences, Shanghai 200031, China
- Institute of Systems and Physical Biology, Shenzhen Bay Laboratory, Shenzhen 518055, China
| | - Yan Lu
- State Key Laboratory of Genetic Engineering, Human Phenome Institute, Zhangjiang Fudan International Innovation Center, Center for Evolutionary Biology, School of Life Sciences, Department of Liver Surgery and Transplantation Liver Cancer Institute, Zhongshan Hospital, Fudan University, Shanghai, China
- Ministry of Education Key Laboratory of Contemporary Anthropology, Fudan University, Shanghai China
| | - Shuhua Xu
- State Key Laboratory of Genetic Engineering, Human Phenome Institute, Zhangjiang Fudan International Innovation Center, Center for Evolutionary Biology, School of Life Sciences, Department of Liver Surgery and Transplantation Liver Cancer Institute, Zhongshan Hospital, Fudan University, Shanghai, China
| |
Collapse
|
3
|
Reyna-Blanco CS, Caduff M, Galimberti M, Leuenberger C, Wegmann D. Inference of Locus-Specific Population Mixtures from Linked Genome-Wide Allele Frequencies. Mol Biol Evol 2024; 41:msae137. [PMID: 38958167 PMCID: PMC11255385 DOI: 10.1093/molbev/msae137] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2023] [Revised: 06/26/2024] [Accepted: 06/27/2024] [Indexed: 07/04/2024] Open
Abstract
Admixture between populations and species is common in nature. Since the influx of new genetic material might be either facilitated or hindered by selection, variation in mixture proportions along the genome is expected in organisms undergoing recombination. Various graph-based models have been developed to better understand these evolutionary dynamics of population splits and mixtures. However, current models assume a single mixture rate for the entire genome and do not explicitly account for linkage. Here, we introduce TreeSwirl, a novel method for inferring branch lengths and locus-specific mixture proportions by using genome-wide allele frequency data, assuming that the admixture graph is known or has been inferred. TreeSwirl builds upon TreeMix that uses Gaussian processes to estimate the presence of gene flow between diverged populations. However, in contrast to TreeMix, our model infers locus-specific mixture proportions employing a hidden Markov model that accounts for linkage. Through simulated data, we demonstrate that TreeSwirl can accurately estimate locus-specific mixture proportions and handle complex demographic scenarios. It also outperforms related D- and f-statistics in terms of accuracy and sensitivity to detect introgressed loci.
Collapse
Affiliation(s)
- Carlos S Reyna-Blanco
- Department of Biology, University of Fribourg, Fribourg 1700, Switzerland
- Swiss Institute of Bioinformatics, Fribourg 1700, Switzerland
| | - Madleina Caduff
- Department of Biology, University of Fribourg, Fribourg 1700, Switzerland
- Swiss Institute of Bioinformatics, Fribourg 1700, Switzerland
| | - Marco Galimberti
- Department of Biology, University of Fribourg, Fribourg 1700, Switzerland
- Swiss Institute of Bioinformatics, Fribourg 1700, Switzerland
- Department of Psychiatry, Yale University School of Medicine, New Haven, CT, USA
- Veterans Affairs Connecticut Healthcare System, West Haven, CT, USA
| | | | - Daniel Wegmann
- Department of Biology, University of Fribourg, Fribourg 1700, Switzerland
- Swiss Institute of Bioinformatics, Fribourg 1700, Switzerland
| |
Collapse
|
4
|
Yermakovich D, André M, Brucato N, Kariwiga J, Leavesley M, Pankratov V, Mondal M, Ricaut FX, Dannemann M. Denisovan admixture facilitated environmental adaptation in Papua New Guinean populations. Proc Natl Acad Sci U S A 2024; 121:e2405889121. [PMID: 38889149 PMCID: PMC11214076 DOI: 10.1073/pnas.2405889121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2024] [Accepted: 05/16/2024] [Indexed: 06/20/2024] Open
Abstract
Neandertals and Denisovans, having inhabited distinct regions in Eurasia and possibly Oceania for over 200,000 y, experienced ample time to adapt to diverse environmental challenges these regions presented. Among present-day human populations, Papua New Guineans (PNG) stand out as one of the few carrying substantial amounts of both Neandertal and Denisovan DNA, a result of past admixture events with these archaic human groups. This study investigates the distribution of introgressed Denisovan and Neandertal DNA within two distinct PNG populations, residing in the highlands of Mt Wilhelm and the lowlands of Daru Island. These locations exhibit unique environmental features, some of which may parallel the challenges that archaic humans once confronted and adapted to. Our results show that PNG highlanders carry higher levels of Denisovan DNA compared to PNG lowlanders. Among the Denisovan-like haplotypes with higher frequencies in highlander populations, those exhibiting the greatest frequency difference compared to lowlander populations also demonstrate more pronounced differences in population frequencies than frequency-matched nonarchaic variants. Two of the five most highly differentiated of those haplotypes reside in genomic areas linked to brain development genes. Conversely, Denisovan-like haplotypes more frequent in lowlanders overlap with genes associated with immune response processes. Our findings suggest that Denisovan DNA has provided genetic variation associated with brain biology and immune response to PNG genomes, some of which might have facilitated adaptive processes to environmental challenges.
Collapse
Affiliation(s)
- Danat Yermakovich
- Center of Genomics, Evolution and Medicine, Institute of Genomics, University of Tartu, Tartu51010, Estonia
| | - Mathilde André
- Center of Genomics, Evolution and Medicine, Institute of Genomics, University of Tartu, Tartu51010, Estonia
| | - Nicolas Brucato
- Centre de Recherche sur la Biodiversité et l'Environnement, Université de Toulouse, Centre National de la Recherche Scientifique, Institut de Recherche pour le Développement, Toulouse Institut National Polytechnique, Université Toulouse 3–Paul Sabatier, cedex 9, Toulouse31062, France
| | - Jason Kariwiga
- Strand of Anthropology, Sociology and Archaeology, School of Humanities and Social Sciences, University of Papua New Guinea, PO Box 320, University 134, National Capital District, Papua New Guinea
- School of Social Science, University of Queensland, St. Lucia, QLD4072, Australia
| | - Matthew Leavesley
- Strand of Anthropology, Sociology and Archaeology, School of Humanities and Social Sciences, University of Papua New Guinea, PO Box 320, University 134, National Capital District, Papua New Guinea
- The Australian Research Council Centre of Excellence for Australian Biodiversity and Heritage & College of Arts, Society and Education, James Cook University, Cairns, QLD4870, Australia
| | - Vasili Pankratov
- Center of Genomics, Evolution and Medicine, Institute of Genomics, University of Tartu, Tartu51010, Estonia
| | - Mayukh Mondal
- Center of Genomics, Evolution and Medicine, Institute of Genomics, University of Tartu, Tartu51010, Estonia
- Institute of Clinical Molecular Biology, Christian-Albrechts-Universität zu Kiel, Kiel24118, Germany
| | - François-Xavier Ricaut
- Centre de Recherche sur la Biodiversité et l'Environnement, Université de Toulouse, Centre National de la Recherche Scientifique, Institut de Recherche pour le Développement, Toulouse Institut National Polytechnique, Université Toulouse 3–Paul Sabatier, cedex 9, Toulouse31062, France
| | - Michael Dannemann
- Center of Genomics, Evolution and Medicine, Institute of Genomics, University of Tartu, Tartu51010, Estonia
| |
Collapse
|
5
|
Roca-Ayats N, Maceda I, Bruque CD, Martínez-Gil N, Garcia-Giralt N, Cozar M, Mellibovsky L, Van Hul W, Lao O, Grinberg D, Balcells S. Evolutionary and functional analyses of LRP5 in archaic and extant modern humans. Hum Genomics 2024; 18:53. [PMID: 38802968 PMCID: PMC11131306 DOI: 10.1186/s40246-024-00616-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2024] [Accepted: 05/07/2024] [Indexed: 05/29/2024] Open
Abstract
BACKGROUND The human lineage has undergone a postcranial skeleton gracilization (i.e. lower bone mass and strength relative to body size) compared to other primates and archaic populations such as the Neanderthals. This gracilization has been traditionally explained by differences in the mechanical load that our ancestors exercised. However, there is growing evidence that gracilization could also be genetically influenced. RESULTS We have analyzed the LRP5 gene, which is known to be associated with high bone mineral density conditions, from an evolutionary and functional point of view. Taking advantage of the published genomes of archaic Homo populations, our results suggest that this gene has a complex evolutionary history both between archaic and living humans and within living human populations. In particular, we identified the presence of different selective pressures in archaics and extant modern humans, as well as evidence of positive selection in the African and South East Asian populations from the 1000 Genomes Project. Furthermore, we observed a very limited evidence of archaic introgression in this gene (only at three haplotypes of East Asian ancestry out of the 1000 Genomes), compatible with a general erasing of the fingerprint of archaic introgression due to functional differences in archaics compared to extant modern humans. In agreement with this hypothesis, we observed private mutations in the archaic genomes that we experimentally validated as putatively increasing bone mineral density. In particular, four of five archaic missense mutations affecting the first β-propeller of LRP5 displayed enhanced Wnt pathway activation, of which two also displayed reduced negative regulation. CONCLUSIONS In summary, these data suggest a genetic component contributing to the understanding of skeletal differences between extant modern humans and archaic Homo populations.
Collapse
Affiliation(s)
- Neus Roca-Ayats
- Departament de Genètica, Microbiologia i Estadística and IBUB, Universitat de Barcelona, Barcelona, Spain
- Centro de Investigación Biomédica en Red de Enfermedades Raras (CIBERER) ISCIII, Barcelona, Spain
- Institut de Recerca Sant Joan de Déu (IRSJD), Barcelona, Spain
| | - Iago Maceda
- CNAG, Centre Nacional d'Analisi Genòmic, C/ Baldiri I Reixach 4, 08028, Barcelona, Spain
- Barcelona Institute of Science and Technology (BIST), Barcelona, Spain
- Universitat Pompeu Fabra (UPF), Barcelona, Spain
| | - Carlos David Bruque
- Unidad de Conocimiento Traslacional Hospitalaria Patagónica, Hospital de Alta Complejidad El Calafate - S.A.M.I.C., Santa Cruz, Argentina
| | - Núria Martínez-Gil
- Departament de Genètica, Microbiologia i Estadística and IBUB, Universitat de Barcelona, Barcelona, Spain
- Centro de Investigación Biomédica en Red de Enfermedades Raras (CIBERER) ISCIII, Barcelona, Spain
- Institut de Recerca Sant Joan de Déu (IRSJD), Barcelona, Spain
| | - Natàlia Garcia-Giralt
- Musculoskeletal Research Group, IMIM (Hospital del Mar Medical Research Institute), Centro de Investigación Biomédica en Red en Fragilidad y Envejecimiento Saludable (CIBERFES), ISCIII, Departament de Genètica, Microbiologia i Estadística, UB, Barcelona, Spain
| | - Mónica Cozar
- Departament de Genètica, Microbiologia i Estadística and IBUB, Universitat de Barcelona, Barcelona, Spain
- Centro de Investigación Biomédica en Red de Enfermedades Raras (CIBERER) ISCIII, Barcelona, Spain
- Institut de Recerca Sant Joan de Déu (IRSJD), Barcelona, Spain
| | - Leonardo Mellibovsky
- Musculoskeletal Research Group, IMIM (Hospital del Mar Medical Research Institute), Centro de Investigación Biomédica en Red en Fragilidad y Envejecimiento Saludable (CIBERFES), ISCIII, Barcelona, Spain
| | - Wim Van Hul
- Center of Medical Genetics, University of Antwerp, 2650, Antwerp, Belgium
| | - Oscar Lao
- Institute of Evolutionary Biology, CSIC-Universitat Pompeu Fabra, 08003, Barcelona, Spain.
| | - Daniel Grinberg
- Departament de Genètica, Microbiologia i Estadística and IBUB, Universitat de Barcelona, Barcelona, Spain
- Centro de Investigación Biomédica en Red de Enfermedades Raras (CIBERER) ISCIII, Barcelona, Spain
- Institut de Recerca Sant Joan de Déu (IRSJD), Barcelona, Spain
| | - Susanna Balcells
- Departament de Genètica, Microbiologia i Estadística and IBUB, Universitat de Barcelona, Barcelona, Spain
- Centro de Investigación Biomédica en Red de Enfermedades Raras (CIBERER) ISCIII, Barcelona, Spain
- Institut de Recerca Sant Joan de Déu (IRSJD), Barcelona, Spain
| |
Collapse
|
6
|
Iasi LNM, Chintalapati M, Skov L, Mesa AB, Hajdinjak M, Peter BM, Moorjani P. Neandertal ancestry through time: Insights from genomes of ancient and present-day humans. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.05.13.593955. [PMID: 38798350 PMCID: PMC11118355 DOI: 10.1101/2024.05.13.593955] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/29/2024]
Abstract
Gene flow from Neandertals has shaped the landscape of genetic and phenotypic variation in modern humans. We identify the location and size of introgressed Neandertal ancestry segments in more than 300 genomes spanning the last 50,000 years. We study how Neandertal ancestry is shared among individuals to infer the time and duration of the Neandertal gene flow. We find the correlation of Neandertal segment locations across individuals and their divergence to sequenced Neandertals, both support a model of single major Neandertal gene flow. Our catalog of introgressed segments through time confirms that most natural selection-positive and negative-on Neandertal ancestry variants occurred immediately after the gene flow, and provides new insights into how the contact with Neandertals shaped human origins and adaptation.
Collapse
Affiliation(s)
- Leonardo N. M. Iasi
- Department for Evolutionary Genetics, Max-Planck-Institute for Evolutionary Anthropology; Leipzig, 04301, Germany
| | - Manjusha Chintalapati
- Department of Molecular and Cell Biology, University of California Berkeley; Berkeley, CA 94720, USA
| | - Laurits Skov
- Department of Molecular and Cell Biology, University of California Berkeley; Berkeley, CA 94720, USA
| | - Alba Bossoms Mesa
- Department for Evolutionary Genetics, Max-Planck-Institute for Evolutionary Anthropology; Leipzig, 04301, Germany
| | - Mateja Hajdinjak
- Department for Evolutionary Genetics, Max-Planck-Institute for Evolutionary Anthropology; Leipzig, 04301, Germany
- The Francis Crick Institute; London, NW1 1AT, UK
| | - Benjamin M. Peter
- Department for Evolutionary Genetics, Max-Planck-Institute for Evolutionary Anthropology; Leipzig, 04301, Germany
- Department of Biology, University of Rochester; Rochester NY, 14620,USA
| | - Priya Moorjani
- Department of Molecular and Cell Biology, University of California Berkeley; Berkeley, CA 94720, USA
- Center for Computational Biology, University of California Berkeley; Berkeley, CA 94720, USA
| |
Collapse
|
7
|
Carrión-Castillo A, Boeckx C. Insights into the genetic architecture of cerebellar lobules derived from the UK Biobank. Sci Rep 2024; 14:9488. [PMID: 38664414 DOI: 10.1038/s41598-024-59699-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2023] [Accepted: 04/15/2024] [Indexed: 06/19/2024] Open
Abstract
In this work we endeavor to further understand the genetic architecture of the cerebellum by examining the genetic underpinnings of the different cerebellar lob(ul)es, identifying their genetic relation to cortical and subcortical regions, as well as to psychiatric disorders, as well as traces of their evolutionary trajectories. We confirm the moderate heritability of cerebellar volumes, and reveal genetic clustering and variability across their different substructures, which warranted a detailed analysis using this higher structural resolution. We replicated known genetic correlations with several subcortical volumes, and report new cortico-cerebellar genetic correlations, including negative genetic correlations between anterior cerebellar lobules and cingulate, and positive ones between lateral Crus I and lobule VI with cortical measures in the fusiform region. Heritability partitioning for evolutionary annotations highlighted that the vermis of Crus II has depleted heritability in genomic regions of "archaic introgression deserts", but no enrichment/depletion of heritability in any other cerebellar regions. Taken together, these findings reveal novel insights into the genetic underpinnings of the different cerebellar lobules.
Collapse
Affiliation(s)
- Amaia Carrión-Castillo
- Basque Center on Cognition, Brain and Language (BCBL), Donostia-San Sebastián, Spain.
- Ikerbasque, Basque Foundation for Science, Bilbao, Spain.
| | - Cedric Boeckx
- Universitat de Barcelona, Barcelona, Spain.
- Universitat de Barcelona Institute of Complex Systems, Barcelona, Spain.
- Universitat de Barcelona Institute of Neurosciences, Barcelona, Spain.
- Catalan Institute for Research and Advanced Studies (ICREA), Barcelona, Spain.
| |
Collapse
|
8
|
Liu X, Koyama S, Tomizuka K, Takata S, Ishikawa Y, Ito S, Kosugi S, Suzuki K, Hikino K, Koido M, Koike Y, Horikoshi M, Gakuhari T, Ikegawa S, Matsuda K, Momozawa Y, Ito K, Kamatani Y, Terao C. Decoding triancestral origins, archaic introgression, and natural selection in the Japanese population by whole-genome sequencing. SCIENCE ADVANCES 2024; 10:eadi8419. [PMID: 38630824 PMCID: PMC11023554 DOI: 10.1126/sciadv.adi8419] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/21/2023] [Accepted: 03/07/2024] [Indexed: 04/19/2024]
Abstract
We generated Japanese Encyclopedia of Whole-Genome/Exome Sequencing Library (JEWEL), a high-depth whole-genome sequencing dataset comprising 3256 individuals from across Japan. Analysis of JEWEL revealed genetic characteristics of the Japanese population that were not discernible using microarray data. First, rare variant-based analysis revealed an unprecedented fine-scale genetic structure. Together with population genetics analysis, the present-day Japanese can be decomposed into three ancestral components. Second, we identified unreported loss-of-function (LoF) variants and observed that for specific genes, LoF variants appeared to be restricted to a more limited set of transcripts than would be expected by chance, with PTPRD as a notable example. Third, we identified 44 archaic segments linked to complex traits, including a Denisovan-derived segment at NKX6-1 associated with type 2 diabetes. Most of these segments are specific to East Asians. Fourth, we identified candidate genetic loci under recent natural selection. Overall, our work provided insights into genetic characteristics of the Japanese population.
Collapse
Affiliation(s)
- Xiaoxi Liu
- Laboratory for Statistical and Translational Genetics, RIKEN Center for Integrative Medical Sciences, Yokohama, Japan
- Clinical Research Center, Shizuoka General Hospital, Shizuoka, Japan
| | - Satoshi Koyama
- Laboratory for Cardiovascular Genomics and Informatics, RIKEN Center for Integrative Medical Sciences, Yokohama, Japan
- Medical and Population Genetics and Cardiovascular Disease Initiative, Broad Institute of Harvard and MIT, Boston, MA, USA
- Cardiovascular Research Center, Massachusetts General Hospital, Boston, MA, USA
| | - Kohei Tomizuka
- Laboratory for Statistical and Translational Genetics, RIKEN Center for Integrative Medical Sciences, Yokohama, Japan
| | - Sadaaki Takata
- Laboratory for Genotyping Development, RIKEN Center for Integrative Medical Sciences, Yokohama, Japan
| | - Yuki Ishikawa
- Laboratory for Statistical and Translational Genetics, RIKEN Center for Integrative Medical Sciences, Yokohama, Japan
| | - Shuji Ito
- Laboratory for Statistical and Translational Genetics, RIKEN Center for Integrative Medical Sciences, Yokohama, Japan
- Laboratory for Bone and Joint Diseases, RIKEN Center for Medical Sciences, Tokyo, Japan
- Department of Orthopedic Surgery, Faculty of Medicine, Shimane University, Izumo, Japan
| | - Shunichi Kosugi
- Laboratory for Statistical and Translational Genetics, RIKEN Center for Integrative Medical Sciences, Yokohama, Japan
| | - Kunihiko Suzuki
- Laboratory for Genotyping Development, RIKEN Center for Integrative Medical Sciences, Yokohama, Japan
| | - Keiko Hikino
- Laboratory for Pharmacogenomics, RIKEN Center for Integrative Medical Sciences, Yokohama, Japan
| | - Masaru Koido
- Laboratory for Statistical and Translational Genetics, RIKEN Center for Integrative Medical Sciences, Yokohama, Japan
- Laboratory of Complex Trait Genomics, Department of Computational Biology and Medical Sciences, Graduate School of Frontier Sciences, The University of Tokyo, Tokyo, Japan
| | - Yoshinao Koike
- Laboratory for Statistical and Translational Genetics, RIKEN Center for Integrative Medical Sciences, Yokohama, Japan
- Laboratory for Bone and Joint Diseases, RIKEN Center for Medical Sciences, Tokyo, Japan
- Department of Orthopedic Surgery, Hokkaido University Graduate School of Medicine, Sapporo, Japan
| | - Momoko Horikoshi
- Laboratory for Genomics of Diabetes and Metabolism, RIKEN Center for Integrative Medical Sciences, Yokohama, Japan
| | - Takashi Gakuhari
- Institute for the Study of Ancient Civilizations and Cultural Resources, College of Human and Social Sciences, Kanazawa University, Kanazawa, Japan
| | - Shiro Ikegawa
- Laboratory for Bone and Joint Diseases, RIKEN Center for Medical Sciences, Tokyo, Japan
| | - Kochi Matsuda
- Laboratory of Genome Technology, Human Genome Center, Institute of Medical Science, The University of Tokyo, Tokyo, Japan
- Laboratory of Clinical Genome Sequencing, Department of Computational Biology and Medical Sciences, Graduate School of Frontier Sciences, The University of Tokyo, Tokyo, Japan
| | - Yukihide Momozawa
- Laboratory for Genotyping Development, RIKEN Center for Integrative Medical Sciences, Yokohama, Japan
| | - Kaoru Ito
- Laboratory for Cardiovascular Genomics and Informatics, RIKEN Center for Integrative Medical Sciences, Yokohama, Japan
| | - Yoichiro Kamatani
- Laboratory for Statistical and Translational Genetics, RIKEN Center for Integrative Medical Sciences, Yokohama, Japan
- Laboratory of Complex Trait Genomics, Department of Computational Biology and Medical Sciences, Graduate School of Frontier Sciences, The University of Tokyo, Tokyo, Japan
| | - Chikashi Terao
- Laboratory for Statistical and Translational Genetics, RIKEN Center for Integrative Medical Sciences, Yokohama, Japan
- Clinical Research Center, Shizuoka General Hospital, Shizuoka, Japan
- The Department of Applied Genetics, The School of Pharmaceutical Sciences, University of Shizuoka, Shizuoka, Japan
| |
Collapse
|
9
|
Zeberg H, Jakobsson M, Pääbo S. The genetic changes that shaped Neandertals, Denisovans, and modern humans. Cell 2024; 187:1047-1058. [PMID: 38367615 DOI: 10.1016/j.cell.2023.12.029] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2023] [Revised: 11/20/2023] [Accepted: 12/20/2023] [Indexed: 02/19/2024]
Abstract
Modern human ancestors diverged from the ancestors of Neandertals and Denisovans about 600,000 years ago. Until about 40,000 years ago, these three groups existed in parallel, occasionally met, and exchanged genes. A critical question is why modern humans, and not the other two groups, survived, became numerous, and developed complex cultures. Here, we discuss genetic differences among the groups and some of their functional consequences. As more present-day genome sequences become available from diverse groups, we predict that very few, if any, differences will distinguish all modern humans from all Neandertals and Denisovans. We propose that the genetic basis of what constitutes a modern human is best thought of as a combination of genetic features, where perhaps none of them is present in each and every present-day individual.
Collapse
Affiliation(s)
- Hugo Zeberg
- Max Planck Institute for Evolutionary Anthropology, Deutscher Platz 6, 04103 Leipzig, Germany; Department of Physiology and Pharmacology, Karolinska Institutet, 17165 Stockholm, Sweden.
| | - Mattias Jakobsson
- Department of Organismal Biology, Uppsala University, 75236 Uppsala, Sweden
| | - Svante Pääbo
- Max Planck Institute for Evolutionary Anthropology, Deutscher Platz 6, 04103 Leipzig, Germany; Okinawa Institute of Science and Technology, Onnason 904-0495, Okinawa, Japan.
| |
Collapse
|
10
|
Kerdoncuff E, Skov L, Patterson N, Zhao W, Lueng YY, Schellenberg GD, Smith JA, Dey S, Ganna A, Dey AB, Kardia SL, Lee J, Moorjani P. 50,000 years of Evolutionary History of India: Insights from ~2,700 Whole Genome Sequences. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.02.15.580575. [PMID: 38405782 PMCID: PMC10888882 DOI: 10.1101/2024.02.15.580575] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/27/2024]
Abstract
India has been underrepresented in whole genome sequencing studies. We generated 2,762 high coverage genomes from India-including individuals from most geographic regions, speakers of all major languages, and tribal and caste groups-providing a comprehensive survey of genetic variation in India. With these data, we reconstruct the evolutionary history of India through space and time at fine scales. We show that most Indians derive ancestry from three ancestral groups related to ancient Iranian farmers, Eurasian Steppe pastoralists and South Asian hunter-gatherers. We uncover a common source of Iranian-related ancestry from early Neolithic cultures of Central Asia into the ancestors of Ancestral South Indians (ASI), Ancestral North Indians (ANI), Austro-asiatic-related and East Asian-related groups in India. Following these admixtures, India experienced a major demographic shift towards endogamy, resulting in extensive homozygosity and identity-by-descent sharing among individuals. At deep time scales, Indians derive around 1-2% of their ancestry from gene flow from archaic hominins, Neanderthals and Denisovans. By assembling the surviving fragments of archaic ancestry in modern Indians, we recover ~1.5 Gb (or 50%) of the introgressing Neanderthal and ~0.6 Gb (or 20%) of the introgressing Denisovan genomes, more than any other previous archaic ancestry study. Moreover, Indians have the largest variation in Neanderthal ancestry, as well as the highest amount of population-specific Neanderthal segments among worldwide groups. Finally, we demonstrate that most of the genetic variation in Indians stems from a single major migration out of Africa that occurred around 50,000 years ago, with minimal contribution from earlier migration waves. Together, these analyses provide a detailed view of the population history of India and underscore the value of expanding genomic surveys to diverse groups outside Europe.
Collapse
Affiliation(s)
- Elise Kerdoncuff
- Department of Molecular and Cell Biology, University of California, Berkeley, United States of America
| | - Laurits Skov
- Department of Molecular and Cell Biology, University of California, Berkeley, United States of America
| | - Nick Patterson
- Department of Human Evolutionary Biology, Harvard University, Cambridge, Massachusetts, United States of America
| | - Wei Zhao
- Survey Research Center, Institute for Social Research, University of Michigan, Ann Arbor, Michigan, United States of America
- Department of Epidemiology, School of Public Health, University of Michigan, Ann Arbor, Michigan, United States of America
| | - Yuk Yee Lueng
- Neurodegeneration Genomics Center, Department of Pathology and Laboratory Medicine, Perelman School of Medicine, University of Pennsylvania, United States of America
| | - Gerard D. Schellenberg
- Neurodegeneration Genomics Center, Department of Pathology and Laboratory Medicine, Perelman School of Medicine, University of Pennsylvania, United States of America
| | - Jennifer A. Smith
- Survey Research Center, Institute for Social Research, University of Michigan, Ann Arbor, Michigan, United States of America
- Department of Epidemiology, School of Public Health, University of Michigan, Ann Arbor, Michigan, United States of America
| | - Sharmistha Dey
- Department of Biophysics, All India Institute of Medical Sciences, New Delhi, India
| | - Andrea Ganna
- Institute for Molecular Medicine Finland, Helsinki, Finland
| | - AB Dey
- Department of Geriatric Medicine, All India Institute of Medical Sciences, New Delhi, India
| | - Sharon L.R. Kardia
- Department of Epidemiology, School of Public Health, University of Michigan, Ann Arbor, Michigan, United States of America
| | - Jinkook Lee
- Department of Economics, and Center for Economic & Social Research, University of Southern California, Los Angeles, California, United States of America
| | - Priya Moorjani
- Department of Molecular and Cell Biology, University of California, Berkeley, United States of America
- Center for Computational Biology, University of California, Berkeley, United States of America
| |
Collapse
|
11
|
Peyrégne S, Slon V, Kelso J. More than a decade of genetic research on the Denisovans. Nat Rev Genet 2024; 25:83-103. [PMID: 37723347 DOI: 10.1038/s41576-023-00643-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 07/19/2023] [Indexed: 09/20/2023]
Abstract
Denisovans, a group of now extinct humans who lived in Eastern Eurasia in the Middle and Late Pleistocene, were first identified from DNA sequences just over a decade ago. Only ten fragmentary remains from two sites have been attributed to Denisovans based entirely on molecular information. Nevertheless, there has been great interest in using genetic data to understand Denisovans and their place in human history. From the reconstruction of a single high-quality genome, it has been possible to infer their population history, including events of admixture with other human groups. Additionally, the identification of Denisovan DNA in the genomes of present-day individuals has provided insights into the timing and routes of dispersal of ancient modern humans into Asia and Oceania, as well as the contributions of archaic DNA to the physiology of present-day people. In this Review, we synthesize more than a decade of research on Denisovans, reconcile controversies and summarize insights into their population history and phenotype. We also highlight how our growing knowledge about Denisovans has provided insights into our own evolutionary history.
Collapse
Affiliation(s)
- Stéphane Peyrégne
- Department of Evolutionary Genetics, Max-Planck-Institute for Evolutionary Anthropology, Leipzig, Germany.
| | - Viviane Slon
- Department of Evolutionary Genetics, Max-Planck-Institute for Evolutionary Anthropology, Leipzig, Germany
- Department of Anatomy and Anthropology, Faculty of Medicine, Tel Aviv University, Tel Aviv, Israel
- Department of Human Molecular Genetics and Biochemistry, Faculty of Medicine, Tel Aviv University, Tel Aviv, Israel
- The Dan David Center for Human Evolution and Biohistory Research, Tel Aviv University, Tel Aviv, Israel
| | - Janet Kelso
- Department of Evolutionary Genetics, Max-Planck-Institute for Evolutionary Anthropology, Leipzig, Germany.
| |
Collapse
|
12
|
Lopez Fang L, Peede D, Ortega-Del Vecchyo D, McTavish EJ, Huerta-Sánchez E. Leveraging shared ancestral variation to detect local introgression. PLoS Genet 2024; 20:e1010155. [PMID: 38190420 PMCID: PMC10798638 DOI: 10.1371/journal.pgen.1010155] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2022] [Revised: 01/19/2024] [Accepted: 12/04/2023] [Indexed: 01/10/2024] Open
Abstract
Introgression is a common evolutionary phenomenon that results in shared genetic material across non-sister taxa. Existing statistical methods such as Patterson's D statistic can detect introgression by measuring an excess of shared derived alleles between populations. The D statistic is effective to detect genome-wide patterns of introgression but can give spurious inferences of introgression when applied to local regions. We propose a new statistic, D+, that leverages both shared ancestral and derived alleles to infer local introgressed regions. Incorporating both shared derived and ancestral alleles increases the number of informative sites per region, improving our ability to identify local introgression. We use a coalescent framework to derive the expected value of this statistic as a function of different demographic parameters under an instantaneous admixture model and use coalescent simulations to compute the power and precision of D+. While the power of D and D+ is comparable, D+ has better precision than D. We apply D+ to empirical data from the 1000 Genome Project and Heliconius butterflies to infer local targets of introgression in humans and in butterflies.
Collapse
Affiliation(s)
- Lesly Lopez Fang
- Department of Life & Environmental Sciences, University of California, Merced, Merced, California, United States of America
- Quantitative & Systems Biology Graduate Group, University of California, Merced, Merced, California, United States of America
| | - David Peede
- Department of Ecology, Evolution and Organismal Biology, Brown University, Providence, Rhode Island, United States of America
- Center for Computational Biology, Brown University, Providence, Rhode Island, United States of America
- Institute at Brown for Environment and Society, Brown University, Providence, Rhode Island, United States of America
| | - Diego Ortega-Del Vecchyo
- Laboratorio Internacional de Investigación sobre el Genoma Humano, Universidad Nacional Autónoma de México, Santiago de Querétaro, Querétaro, México
| | - Emily Jane McTavish
- Department of Life & Environmental Sciences, University of California, Merced, Merced, California, United States of America
- Quantitative & Systems Biology Graduate Group, University of California, Merced, Merced, California, United States of America
| | - Emilia Huerta-Sánchez
- Department of Ecology, Evolution and Organismal Biology, Brown University, Providence, Rhode Island, United States of America
- Center for Computational Biology, Brown University, Providence, Rhode Island, United States of America
| |
Collapse
|
13
|
Huang X, Rymbekova A, Dolgova O, Lao O, Kuhlwilm M. Harnessing deep learning for population genetic inference. Nat Rev Genet 2024; 25:61-78. [PMID: 37666948 DOI: 10.1038/s41576-023-00636-3] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 07/11/2023] [Indexed: 09/06/2023]
Abstract
In population genetics, the emergence of large-scale genomic data for various species and populations has provided new opportunities to understand the evolutionary forces that drive genetic diversity using statistical inference. However, the era of population genomics presents new challenges in analysing the massive amounts of genomes and variants. Deep learning has demonstrated state-of-the-art performance for numerous applications involving large-scale data. Recently, deep learning approaches have gained popularity in population genetics; facilitated by the advent of massive genomic data sets, powerful computational hardware and complex deep learning architectures, they have been used to identify population structure, infer demographic history and investigate natural selection. Here, we introduce common deep learning architectures and provide comprehensive guidelines for implementing deep learning models for population genetic inference. We also discuss current challenges and future directions for applying deep learning in population genetics, focusing on efficiency, robustness and interpretability.
Collapse
Affiliation(s)
- Xin Huang
- Department of Evolutionary Anthropology, University of Vienna, Vienna, Austria.
- Human Evolution and Archaeological Sciences (HEAS), University of Vienna, Vienna, Austria.
| | - Aigerim Rymbekova
- Department of Evolutionary Anthropology, University of Vienna, Vienna, Austria
- Human Evolution and Archaeological Sciences (HEAS), University of Vienna, Vienna, Austria
| | - Olga Dolgova
- Integrative Genomics Laboratory, CIC bioGUNE - Centro de Investigación Cooperativa en Biociencias, Derio, Biscaya, Spain
| | - Oscar Lao
- Institute of Evolutionary Biology, CSIC-Universitat Pompeu Fabra, Barcelona, Spain.
| | - Martin Kuhlwilm
- Department of Evolutionary Anthropology, University of Vienna, Vienna, Austria.
- Human Evolution and Archaeological Sciences (HEAS), University of Vienna, Vienna, Austria.
| |
Collapse
|
14
|
Di Santo LN, Quilodrán CS, Currat M. Temporal Variation in Introgressed Segments' Length Statistics Computed from a Limited Number of Ancient Genomes Sheds Light on Past Admixture Pulses. Mol Biol Evol 2023; 40:msad252. [PMID: 37992125 PMCID: PMC10715198 DOI: 10.1093/molbev/msad252] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2023] [Revised: 10/16/2023] [Accepted: 11/09/2023] [Indexed: 11/24/2023] Open
Abstract
Hybridization is recognized as an important evolutionary force, but identifying and timing admixture events between divergent lineages remain a major aim of evolutionary biology. While this has traditionally been done using inferential tools on contemporary genomes, the latest advances in paleogenomics have provided a growing wealth of temporally distributed genomic data. Here, we used individual-based simulations to generate chromosome-level genomic data for a 2-population system and described temporal neutral introgression patterns under a single- and 2-pulse admixture model. We computed 6 summary statistics aiming to inform the timing and number of admixture pulses between interbreeding entities: lengths of introgressed sequences and their variance within genomes, as well as genome-wide introgression proportions and related measures. The first 2 statistics could confidently be used to infer interlineage hybridization history, peaking at the beginning and shortly after an admixture pulse. Temporal variation in introgression proportions and related statistics provided more limited insights, particularly when considering their application to ancient genomes still scant in number. Lastly, we computed these statistics on Homo sapiens paleogenomes and successfully inferred the hybridization pulse from Neanderthal that occurred approximately 40 to 60 kya. The scarce number of genomes dating from this period prevented more precise inferences, but the accumulation of paleogenomic data opens promising perspectives as our approach only requires a limited number of ancient genomes.
Collapse
Affiliation(s)
- Lionel N Di Santo
- Department of Genetics and Evolution, University of Geneva, Geneva CH-1205
| | | | - Mathias Currat
- Department of Genetics and Evolution, University of Geneva, Geneva CH-1205
- Institute of Genetics and Genomics in Geneva (IGE3), University of Geneva, Geneva CH-1205
| |
Collapse
|
15
|
Wroblewski TH, Witt KE, Lee SB, Malhi RS, Peede D, Huerta-Sánchez E, Villanea FA, Claw KG. Pharmacogenetic Variation in Neanderthals and Denisovans and Implications for Human Health and Response to Medications. Genome Biol Evol 2023; 15:evad222. [PMID: 38051947 PMCID: PMC10727477 DOI: 10.1093/gbe/evad222] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2023] [Revised: 11/08/2023] [Accepted: 11/16/2023] [Indexed: 12/07/2023] Open
Abstract
Modern humans carry both Neanderthal and Denisovan (archaic) genome elements that are part of the human gene pool and affect the life and health of living individuals. The impact of archaic DNA may be particularly evident in pharmacogenes-genes responsible for the processing of exogenous substances such as food, pollutants, and medications-as these can relate to changing environmental effects, and beneficial variants may have been retained as modern humans encountered new environments. However, the health implications and contribution of archaic ancestry in pharmacogenes of modern humans remain understudied. Here, we explore 11 key cytochrome P450 genes (CYP450) involved in 75% of all drug metabolizing reactions in three Neanderthal and one Denisovan individuals and examine archaic introgression in modern human populations. We infer the metabolizing efficiency of these 11 CYP450 genes in archaic individuals and find important predicted phenotypic differences relative to modern human variants. We identify several single nucleotide variants shared between archaic and modern humans in each gene, including some potentially function-altering mutations in archaic CYP450 genes, which may result in altered metabolism in living people carrying these variants. We also identified several variants in the archaic CYP450 genes that are novel and unique to archaic humans as well as one gene, CYP2B6, that shows evidence for a gene duplication found only in Neanderthals and modern Africans. Finally, we highlight CYP2A6, CYP2C9, and CYP2J2, genes which show evidence for archaic introgression into modern humans and posit evolutionary hypotheses that explain their allele frequencies in modern populations.
Collapse
Affiliation(s)
- Tadeusz H Wroblewski
- Department of Biomedical Informatics, Colorado Center for Personalized Medicine, University of Colorado Anschutz Medical Campus, Aurora, Colorado, USA
| | - Kelsey E Witt
- Center for Human Genetics and Department of Genetics and Biochemistry, Clemson University, South Carolina, USA
| | - Seung-been Lee
- Precision Medicine Institute, Macrogen Inc., Seoul, Republic of Korea
| | - Ripan S Malhi
- Department of Anthropology and Carl R. Woese Institute for Genomic Biology, University of Illinois Urbana-Champaign, Illinois, USA
| | - David Peede
- Department of Ecology, Evolution, and Organismal Biology and Center for Computational and Molecular Biology, Brown University, Providence, Rhode Island, USA
- Institute at Brown for Environment and Society, Brown University, Providence, Rhode Island, USA
| | - Emilia Huerta-Sánchez
- Department of Ecology, Evolution, and Organismal Biology and Center for Computational and Molecular Biology, Brown University, Providence, Rhode Island, USA
| | | | - Katrina G Claw
- Department of Biomedical Informatics, Colorado Center for Personalized Medicine, University of Colorado Anschutz Medical Campus, Aurora, Colorado, USA
| |
Collapse
|
16
|
Velazquez-Arcelay K, Colbran LL, McArthur E, Brand CM, Rinker DC, Siemann JK, McMahon DG, Capra JA. Archaic Introgression Shaped Human Circadian Traits. Genome Biol Evol 2023; 15:evad203. [PMID: 38095367 PMCID: PMC10719892 DOI: 10.1093/gbe/evad203] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 11/07/2023] [Indexed: 12/17/2023] Open
Abstract
When the ancestors of modern Eurasians migrated out of Africa and interbred with Eurasian archaic hominins, namely, Neanderthals and Denisovans, DNA of archaic ancestry integrated into the genomes of anatomically modern humans. This process potentially accelerated adaptation to Eurasian environmental factors, including reduced ultraviolet radiation and increased variation in seasonal dynamics. However, whether these groups differed substantially in circadian biology and whether archaic introgression adaptively contributed to human chronotypes remain unknown. Here, we traced the evolution of chronotype based on genomes from archaic hominins and present-day humans. First, we inferred differences in circadian gene sequences, splicing, and regulation between archaic hominins and modern humans. We identified 28 circadian genes containing variants with potential to alter splicing in archaics (e.g., CLOCK, PER2, RORB, and RORC) and 16 circadian genes likely divergently regulated between present-day humans and archaic hominins, including RORA. These differences suggest the potential for introgression to modify circadian gene expression. Testing this hypothesis, we found that introgressed variants are enriched among expression quantitative trait loci for circadian genes. Supporting the functional relevance of these regulatory effects, we found that many introgressed alleles have associations with chronotype. Strikingly, the strongest introgressed effects on chronotype increase morningness, consistent with adaptations to high latitude in other species. Finally, we identified several circadian loci with evidence of adaptive introgression or latitudinal clines in allele frequency. These findings identify differences in circadian gene regulation between modern humans and archaic hominins and support the contribution of introgression via coordinated effects on variation in human chronotype.
Collapse
Affiliation(s)
| | - Laura L Colbran
- Department of Genetics, Perelman School of Medicine, University of Pennsylvania, Philadelphia, Pennsylvania, USA
| | - Evonne McArthur
- Vanderbilt Genetics Institute, Vanderbilt University Medical Center, Nashville, Tennessee, USA
| | - Colin M Brand
- Department of Epidemiology and Biostatistics, University of California, SanFrancisco, California, USA
- Bakar Computational Health Sciences Institute, University of California, SanFrancisco, California, USA
| | - David C Rinker
- Department of Biological Sciences, Vanderbilt University, Nashville, Tennessee, USA
| | - Justin K Siemann
- Department of Biological Sciences, Vanderbilt University, Nashville, Tennessee, USA
| | - Douglas G McMahon
- Department of Biological Sciences, Vanderbilt University, Nashville, Tennessee, USA
| | - John A Capra
- Department of Epidemiology and Biostatistics, University of California, SanFrancisco, California, USA
- Bakar Computational Health Sciences Institute, University of California, SanFrancisco, California, USA
| |
Collapse
|
17
|
Zhang Y, Zhu Q, Shao Y, Jiang Y, Ouyang Y, Zhang L, Zhang W. Inferring Historical Introgression with Deep Learning. Syst Biol 2023; 72:1013-1038. [PMID: 37257491 DOI: 10.1093/sysbio/syad033] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2022] [Revised: 05/28/2023] [Accepted: 05/30/2023] [Indexed: 06/02/2023] Open
Abstract
Resolving phylogenetic relationships among taxa remains a challenge in the era of big data due to the presence of genetic admixture in a wide range of organisms. Rapidly developing sequencing technologies and statistical tests enable evolutionary relationships to be disentangled at a genome-wide level, yet many of these tests are computationally intensive and rely on phased genotypes, large sample sizes, restricted phylogenetic topologies, or hypothesis testing. To overcome these difficulties, we developed a deep learning-based approach, named ERICA, for inferring genome-wide evolutionary relationships and local introgressed regions from sequence data. ERICA accepts sequence alignments of both population genomic data and multiple genome assemblies, and efficiently identifies discordant genealogy patterns and exchanged regions across genomes when compared with other methods. We further tested ERICA using real population genomic data from Heliconius butterflies that have undergone adaptive radiation and frequent hybridization. Finally, we applied ERICA to characterize hybridization and introgression in wild and cultivated rice, revealing the important role of introgression in rice domestication and adaptation. Taken together, our findings demonstrate that ERICA provides an effective method for teasing apart evolutionary relationships using whole genome data, which can ultimately facilitate evolutionary studies on hybridization and introgression.
Collapse
Affiliation(s)
- Yubo Zhang
- State Key Laboratory of Protein and Plant Gene Research, Peking-Tsinghua Center for Life Sciences, Academy for Advanced Interdisciplinary Studies, Peking University, Beijing 100871, China
| | - Qingjie Zhu
- Chinese Institute for Brain Research, Beijing 102206, China
| | - Yi Shao
- Chinese Institute for Brain Research, Beijing 102206, China
| | - Yanchen Jiang
- State Key Laboratory of Protein and Plant Gene Research, Peking-Tsinghua Center for Life Sciences, Academy for Advanced Interdisciplinary Studies, Peking University, Beijing 100871, China
- State Key Laboratory of Protein and Plant Gene Research, School of Life Sciences, Peking University, Beijing 100871, China
| | - Yidan Ouyang
- National Key Laboratory of Crop Genetic Improvement and National Centre of Plant Gene Research (Wuhan), Hubei Hongshan Laboratory, Huazhong Agricultural University, Wuhan 430070, China
| | - Li Zhang
- Chinese Institute for Brain Research, Beijing 102206, China
| | - Wei Zhang
- State Key Laboratory of Protein and Plant Gene Research, Peking-Tsinghua Center for Life Sciences, Academy for Advanced Interdisciplinary Studies, Peking University, Beijing 100871, China
- State Key Laboratory of Protein and Plant Gene Research, School of Life Sciences, Peking University, Beijing 100871, China
| |
Collapse
|
18
|
Pivirotto AM, Platt A, Patel R, Kumar S, Hey J. Analyses of allele age and fitness impact reveal human beneficial alleles to be older than neutral controls. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.10.09.561569. [PMID: 37873438 PMCID: PMC10592680 DOI: 10.1101/2023.10.09.561569] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/25/2023]
Abstract
A classic population genetic prediction is that alleles experiencing directional selection should swiftly traverse allele frequency space, leaving detectable reductions in genetic variation in linked regions. However, despite this expectation, identifying clear footprints of beneficial allele passage has proven to be surprisingly challenging. We addressed the basic premise underlying this expectation by estimating the ages of large numbers of beneficial and deleterious alleles in a human population genomic data set. Deleterious alleles were found to be young, on average, given their allele frequency. However, beneficial alleles were older on average than non-coding, non-regulatory alleles of the same frequency. This finding is not consistent with directional selection and instead indicates some type of balancing selection. Among derived beneficial alleles, those fixed in the population show higher local recombination rates than those still segregating, consistent with a model in which new beneficial alleles experience an initial period of balancing selection due to linkage disequilibrium with deleterious recessive alleles. Alleles that ultimately fix following a period of balancing selection will leave a modest 'soft' sweep impact on the local variation, consistent with the overall paucity of species-wide 'hard' sweeps in human genomes.
Collapse
Affiliation(s)
| | - Alexander Platt
- Temple University, Department of Biology, Philadelphia PA 19122, USA
- University of Pennsylvania, Department of Genetics, Philadelphia PA 19104, USA
| | - Ravi Patel
- Temple University, Department of Biology, Philadelphia PA 19122, USA
- Institute for Genomics and Evolutionary Medicine, Temple University, PA 19122, USA
| | - Sudhir Kumar
- Temple University, Department of Biology, Philadelphia PA 19122, USA
- Institute for Genomics and Evolutionary Medicine, Temple University, PA 19122, USA
| | - Jody Hey
- Temple University, Department of Biology, Philadelphia PA 19122, USA
| |
Collapse
|
19
|
Teng SN, Svenning JC, Xu C. Large mammals and trees in eastern monsoonal China: anthropogenic losses since the Late Pleistocene and restoration prospects in the Anthropocene. Biol Rev Camb Philos Soc 2023; 98:1607-1632. [PMID: 37102332 DOI: 10.1111/brv.12968] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2022] [Revised: 04/16/2023] [Accepted: 04/18/2023] [Indexed: 04/28/2023]
Abstract
Massive human-induced declines of large-sized animals and trees (megabiota) from the Late Pleistocene to the Anthropocene have resulted in downsized ecosystems across the globe, in which components and functions have been greatly simplified. In response, active restoration projects of extant large-sized species or functional substitutes are needed at large scales to promote ecological processes that are important for ecosystem self-regulation and biodiversity maintenance. Despite the desired global scope of such projects, they have received little attention in East Asia. Here, we synthesise the biogeographical and ecological knowledge of megabiota in ancient and modern China, with relevant data mostly located in eastern monsoonal China (EMC), aiming to assess its potential for restoring functionally intact ecosystems modulated by megabiota. We found that during the Late Pleistocene, 12 mammalian megafaunal (carnivores ≥15 kg and herbivores ≥500 kg) species disappeared from EMC: one carnivore Crocuta ultima (East Asian spotted hyena) and 11 herbivores including six megaherbivores (≥1000 kg). The relative importance of climate change and humans in driving these losses remains debated, despite accumulating evidence in favour of the latter. Later massive depletion of megafauna and large-sized (45-500 kg) herbivores has been closely associated with agricultural expansion and societal development, especially during the late Holocene. While forests rich in large timber trees (33 taxa in written records) were common in the region 2000-3000 years ago, millennial-long logging has resulted in considerable range contractions and at least 39 threatened species. The wide distribution of C. ultima, which likely favoured open or semi-open habitats (like extant spotted hyenas), suggests the existence of mosaic open and closed vegetation in the Late Pleistocene across EMC, in line with a few pollen-based vegetation reconstructions and potentially, or at least partially, reflecting herbivory by herbivorous megafauna. The widespread loss of megaherbivores may have strongly compromised seed dispersal for both megafruit (fleshy fruits with widths ≥40 mm) and non-megafruit plant species in EMC, especially in terms of extra-long-distance (>10 km) dispersal, which is critical for plant species that rely on effective biotic agents to track rapid climate change. The former occurrence of large mammals and trees have translated into rich material and non-material heritages passed down across generations. Several reintroduction projects have been implemented or are under consideration, with the case of Elaphurus davidianus a notable success in recovering wild populations in the middle reaches of the Yangtze River, although trophic interactions with native carnivorous megafauna have not yet been restored. Lessons of dealing with human-wildlife conflicts are key to public support for maintaining landscapes shared with megafauna and large herbivores in the human-dominated Anthropocene. Meanwhile, potential human-wildlife conflicts, e.g. public health risks, need to be scientifically informed and effectively reduced. The Chinese government's strong commitment to improved policies of ecological protection and restoration (e.g. ecological redlines and national parks) provides a solid foundation for a scaling-up contribution to the global scope needed for solving the crisis of biotic downsizing and ecosystem degradation.
Collapse
Affiliation(s)
- Shuqing N Teng
- School of Life Sciences, Nanjing University, Nanjing, 210023, China
| | - Jens-Christian Svenning
- Center for Ecological Dynamics in a Novel Biosphere (ECONOVO) & Center for Biodiversity Dynamics in a Changing World (BIOCHANGE), Department of Biology, Aarhus University, Aarhus, 8000, Denmark
| | - Chi Xu
- School of Life Sciences, Nanjing University, Nanjing, 210023, China
- Key Laboratory of Restoration and Reconstruction of Degraded Ecosystems in northwestern China of Ministry of Education, Ningxia University, Yinchuan, 750021, China
| |
Collapse
|
20
|
Yang C, Zhou Y, Song Y, Wu D, Zeng Y, Nie L, Liu P, Zhang S, Chen G, Xu J, Zhou H, Zhou L, Qian X, Liu C, Tan S, Zhou C, Dai W, Xu M, Qi Y, Wang X, Guo L, Fan G, Wang A, Deng Y, Zhang Y, Jin J, He Y, Guo C, Guo G, Zhou Q, Xu X, Yang H, Wang J, Xu S, Mao Y, Jin X, Ruan J, Zhang G. The complete and fully-phased diploid genome of a male Han Chinese. Cell Res 2023; 33:745-761. [PMID: 37452091 PMCID: PMC10542383 DOI: 10.1038/s41422-023-00849-5] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2023] [Accepted: 06/29/2023] [Indexed: 07/18/2023] Open
Abstract
Since the release of the complete human genome, the priority of human genomic study has now been shifting towards closing gaps in ethnic diversity. Here, we present a fully phased and well-annotated diploid human genome from a Han Chinese male individual (CN1), in which the assemblies of both haploids achieve the telomere-to-telomere (T2T) level. Comparison of this diploid genome with the CHM13 haploid T2T genome revealed significant variations in the centromere. Outside the centromere, we discovered 11,413 structural variations, including numerous novel ones. We also detected thousands of CN1 alleles that have accumulated high substitution rates and a few that have been under positive selection in the East Asian population. Further, we found that CN1 outperforms CHM13 as a reference genome in mapping and variant calling for the East Asian population owing to the distinct structural variants of the two references. Comparison of SNP calling for a large cohort of 8869 Chinese genomes using CN1 and CHM13 as reference respectively showed that the reference bias profoundly impacts rare SNP calling, with nearly 2 million rare SNPs miss-called with different reference genomes. Finally, applying the CN1 as a reference, we discovered 5.80 Mb and 4.21 Mb putative introgression sequences from Neanderthal and Denisovan, respectively, including many East Asian specific ones undetected using CHM13 as the reference. Our analyses reveal the advances of using CN1 as a reference for population genomic studies and paleo-genomic studies. This complete genome will serve as an alternative reference for future genomic studies on the East Asian population.
Collapse
Affiliation(s)
- Chentao Yang
- Center for Genomic Research, International Institutes of Medicine, The Fourth Affiliated Hospital, Zhejiang University School of Medicine, Yiwu, Zhejiang, China
- Center for Evolutionary & Organismal Biology, & Women's Hospital, Zhejiang University School of Medicine, Hangzhou, Zhejiang, China
- BGI-Shenzhen, Shenzhen, Guangdong, China
| | - Yang Zhou
- BGI-Shenzhen, Shenzhen, Guangdong, China
- BGI Research-Wuhan, BGI, Wuhan, Hubei, China
| | - Yanni Song
- Shenzhen Branch, Guangdong Laboratory for Lingnan Modern Agriculture, Genome Analysis Laboratory of the Ministry of Agriculture and Rural Affairs, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, Guangdong, China
| | - Dongya Wu
- Center for Genomic Research, International Institutes of Medicine, The Fourth Affiliated Hospital, Zhejiang University School of Medicine, Yiwu, Zhejiang, China
- Center for Evolutionary & Organismal Biology, & Women's Hospital, Zhejiang University School of Medicine, Hangzhou, Zhejiang, China
- Liangzhu Laboratory, Zhejiang University Medical Center, Hangzhou, Zhejiang, China
- Institute of Crop Science & Institute of Bioinformatics, Zhejiang University, Hangzhou, Zhejiang, China
| | - Yan Zeng
- BGI-Shenzhen, Shenzhen, Guangdong, China
| | - Lei Nie
- BGI-Shenzhen, Shenzhen, Guangdong, China
| | | | - Shilong Zhang
- Bio-X Institutes, Key Laboratory for the Genetics of Developmental and Neuropsychiatric Disorders, Ministry of Education, Shanghai Jiao Tong University, Shanghai, China
| | - Guangji Chen
- BGI-Shenzhen, Shenzhen, Guangdong, China
- College of Life Sciences, University of Chinese Academy of Sciences, Beijing, China
| | - Jinjin Xu
- BGI-Shenzhen, Shenzhen, Guangdong, China
| | - Hongling Zhou
- Shenzhen Branch, Guangdong Laboratory for Lingnan Modern Agriculture, Genome Analysis Laboratory of the Ministry of Agriculture and Rural Affairs, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, Guangdong, China
| | - Long Zhou
- Center for Evolutionary & Organismal Biology, & Women's Hospital, Zhejiang University School of Medicine, Hangzhou, Zhejiang, China
- Liangzhu Laboratory, Zhejiang University Medical Center, Hangzhou, Zhejiang, China
- Innovation Center of Yangtze River Delta, Zhejiang University, Hangzhou, Zhejiang, China
| | - Xiaobo Qian
- BGI-Shenzhen, Shenzhen, Guangdong, China
- College of Life Sciences, University of Chinese Academy of Sciences, Beijing, China
| | - Chenlu Liu
- Life Sciences Institute, Zhejiang University, Hangzhou, Zhejiang, China
| | | | | | - Wei Dai
- BGI-Shenzhen, Shenzhen, Guangdong, China
| | - Mengyang Xu
- BGI-Shenzhen, Shenzhen, Guangdong, China
- BGI-Qingdao, BGI-Shenzhen, Qingdao, Shandong, China
| | - Yanwei Qi
- BGI-Qingdao, BGI-Shenzhen, Qingdao, Shandong, China
| | - Xiaobo Wang
- Shenzhen Branch, Guangdong Laboratory for Lingnan Modern Agriculture, Genome Analysis Laboratory of the Ministry of Agriculture and Rural Affairs, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, Guangdong, China
| | - Lidong Guo
- College of Life Sciences, University of Chinese Academy of Sciences, Beijing, China
- BGI-Qingdao, BGI-Shenzhen, Qingdao, Shandong, China
| | - Guangyi Fan
- BGI-Qingdao, BGI-Shenzhen, Qingdao, Shandong, China
| | - Aijun Wang
- BGI-Qingdao, BGI-Shenzhen, Qingdao, Shandong, China
| | - Yuan Deng
- BGI-Shenzhen, Shenzhen, Guangdong, China
| | - Yong Zhang
- BGI-Shenzhen, Shenzhen, Guangdong, China
| | | | - Yunqiu He
- Center for Genomic Research, International Institutes of Medicine, The Fourth Affiliated Hospital, Zhejiang University School of Medicine, Yiwu, Zhejiang, China
- Center for Evolutionary & Organismal Biology, & Women's Hospital, Zhejiang University School of Medicine, Hangzhou, Zhejiang, China
| | - Chunxue Guo
- BGI-Shenzhen, Shenzhen, Guangdong, China
- BGI-Hangzhou, Hangzhou, Zhejiang, China
| | - Guoji Guo
- School of Medicine, Zhejiang University, Hangzhou, Zhejiang, China
| | - Qing Zhou
- Liangzhu Laboratory, Zhejiang University Medical Center, Hangzhou, Zhejiang, China
- Life Sciences Institute, Zhejiang University, Hangzhou, Zhejiang, China
| | - Xun Xu
- BGI-Shenzhen, Shenzhen, Guangdong, China
| | | | - Jian Wang
- BGI-Shenzhen, Shenzhen, Guangdong, China
| | - Shuhua Xu
- State Key Laboratory of Genetic Engineering, Center for Evolutionary Biology, Collaborative Innovation Center for Genetics and Development, School of Life Sciences, Fudan University, Shanghai, China
- Human Phenome Institute, Zhangjiang Fudan International Innovation Center, and Ministry of Education Key Laboratory of Contemporary Anthropology, Fudan University, Shanghai, China
- Jiangsu Key Laboratory of Phylogenomics & Comparative Genomics, International Joint Center of Genomics of Jiangsu Province School of Life Sciences, Jiangsu Normal University, Xuzhou, Jiangsu, China
- Department of Liver Surgery and Transplantation Liver Cancer Institute, Zhongshan Hospital, Fudan University, Shanghai, China
- Center for Excellence in Animal Evolution and Genetics, Chinese Academy of Sciences, Kunming, Yunnan, China
| | - Yafei Mao
- Bio-X Institutes, Key Laboratory for the Genetics of Developmental and Neuropsychiatric Disorders, Ministry of Education, Shanghai Jiao Tong University, Shanghai, China
| | - Xin Jin
- BGI-Shenzhen, Shenzhen, Guangdong, China
| | - Jue Ruan
- Shenzhen Branch, Guangdong Laboratory for Lingnan Modern Agriculture, Genome Analysis Laboratory of the Ministry of Agriculture and Rural Affairs, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen, Guangdong, China.
| | - Guojie Zhang
- Center for Genomic Research, International Institutes of Medicine, The Fourth Affiliated Hospital, Zhejiang University School of Medicine, Yiwu, Zhejiang, China.
- Center for Evolutionary & Organismal Biology, & Women's Hospital, Zhejiang University School of Medicine, Hangzhou, Zhejiang, China.
- Liangzhu Laboratory, Zhejiang University Medical Center, Hangzhou, Zhejiang, China.
- Innovation Center of Yangtze River Delta, Zhejiang University, Hangzhou, Zhejiang, China.
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, Yunnan, China.
| |
Collapse
|
21
|
Pollen AA, Kilik U, Lowe CB, Camp JG. Human-specific genetics: new tools to explore the molecular and cellular basis of human evolution. Nat Rev Genet 2023; 24:687-711. [PMID: 36737647 PMCID: PMC9897628 DOI: 10.1038/s41576-022-00568-4] [Citation(s) in RCA: 18] [Impact Index Per Article: 18.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 12/08/2022] [Indexed: 02/05/2023]
Abstract
Our ancestors acquired morphological, cognitive and metabolic modifications that enabled humans to colonize diverse habitats, develop extraordinary technologies and reshape the biosphere. Understanding the genetic, developmental and molecular bases for these changes will provide insights into how we became human. Connecting human-specific genetic changes to species differences has been challenging owing to an abundance of low-effect size genetic changes, limited descriptions of phenotypic differences across development at the level of cell types and lack of experimental models. Emerging approaches for single-cell sequencing, genetic manipulation and stem cell culture now support descriptive and functional studies in defined cell types with a human or ape genetic background. In this Review, we describe how the sequencing of genomes from modern and archaic hominins, great apes and other primates is revealing human-specific genetic changes and how new molecular and cellular approaches - including cell atlases and organoids - are enabling exploration of the candidate causal factors that underlie human-specific traits.
Collapse
Affiliation(s)
- Alex A Pollen
- Eli and Edythe Broad Center of Regeneration Medicine and Stem Cell Research, University of California, San Francisco, San Francisco, CA, USA.
- Department of Neurology, University of California, San Francisco, San Francisco, CA, USA.
| | - Umut Kilik
- Institute of Human Biology (IHB), Roche Pharma Research and Early Development, Roche Innovation Center Basel, Basel, Switzerland
- University of Basel, Basel, Switzerland
| | - Craig B Lowe
- Department of Molecular Genetics and Microbiology, Duke University School of Medicine, Durham, NC, USA.
| | - J Gray Camp
- Institute of Human Biology (IHB), Roche Pharma Research and Early Development, Roche Innovation Center Basel, Basel, Switzerland.
- University of Basel, Basel, Switzerland.
| |
Collapse
|
22
|
Bacon AM, Bourgon N, Dufour E, Demeter F, Zanolli C, Westaway KE, Joannes-Boyau R, Duringer P, Ponche JL, Morley MW, Suzzoni E, Frangeul S, Boesch Q, Antoine PO, Boualaphane S, Sichanthongtip P, Sihanam D, Huong NTM, Tuan NA, Fiorillo D, Tombret O, Patole-Edoumba E, Zachwieja A, Luangkhoth T, Souksavatdy V, Dunn TE, Shackelford L, Hublin JJ. Palaeoenvironments and hominin evolutionary dynamics in southeast Asia. Sci Rep 2023; 13:16165. [PMID: 37758744 PMCID: PMC10533506 DOI: 10.1038/s41598-023-43011-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2023] [Accepted: 09/18/2023] [Indexed: 09/29/2023] Open
Abstract
Secure environmental contexts are crucial for hominin interpretation and comparison. The discovery of a Denisovan individual and associated fauna at Tam Ngu Hao 2 (Cobra) Cave, Laos, dating back to 164-131 ka, allows for environmental comparisons between this (sub)tropical site and the Palearctic Denisovan sites of Denisova Cave (Russia) and Baishiya Karst Cave (China). Denisovans from northern latitudes foraged in a mix of forested and open landscapes, including tundra and steppe. Using stable isotope values from the Cobra Cave assemblage, we demonstrate that, despite the presence of nearby canopy forests, the Denisovan individual from Cobra Cave primarily consumed plants and/or animals from open forests and savannah. Using faunal evidence and proxy indicators of climates, results herein highlight a local expansion of rainforest at ~ 130 ka, raising questions about how Denisovans responded to this local climate change. Comparing the diet and habitat of the archaic hominin from Cobra Cave with those of early Homo sapiens from Tam Pà Ling Cave (46-43 ka), Laos, it appears that only our species was able to exploit rainforest resources.
Collapse
Affiliation(s)
- Anne-Marie Bacon
- Université Paris Cité, CNRS, BABEL UMR 8045, 75012, Paris, France.
| | - Nicolas Bourgon
- IsoTROPIC Research Group, Max Planck Institute for Geoanthropology, 07745, Jena, Germany
- Max Planck Institute for Evolutionary Anthropology, Department of Human Evolution, 04103, Leipzig, Germany
| | - Elise Dufour
- UMR 7209 Archéozoologie, Archéobotanique, Sociétés, Pratiques, Environnements, MNHN, CNRS, Paris, France
| | - Fabrice Demeter
- Lundbeck Foundation GeoGenetics Centre, Globe Institute, University of Copenhagen, Copenhagen, Denmark
- Eco-anthropologie (EA), MNHN, CNRS, Université Paris Cité, Musée de l'Homme, 75016, Paris, France
| | - Clément Zanolli
- Univ. Bordeaux, CNRS, MCC, PACEA, UMR 5199, 33600, Pessac, France
| | - Kira E Westaway
- 'Traps' Luminescence Dating Facility, School of Natural Sciences, Macquarie University, Sydney, Australia
| | - Renaud Joannes-Boyau
- Geoarchaeology and Archaeometry Research Group (GARG), Southern Cross University, Lismore, NSW, Australia
| | - Philippe Duringer
- Ecole et Observatoire des Sciences de la Terre, Institut de Physique du Globe de Strasbourg, UMR 7516 CNRS, Université de Strasbourg, Strasbourg, France
| | - Jean-Luc Ponche
- Laboratoire Image, Ville Environnement, UMR 7362 UdS CNRS, Université de Strasbourg, Strasbourg, France
| | - Mike W Morley
- Flinders Microarchaeology Laboratory, Archaeology, College of Humanities and Social Sciences, Flinders University, Sturt Road, Bedford Park, Adelaide, SA, 5042, Australia
| | - Eric Suzzoni
- Spitteurs Pan, Technical Cave Supervision and Exploration, La Chapelle en Vercors, France
| | - Sébastien Frangeul
- Spitteurs Pan, Technical Cave Supervision and Exploration, La Chapelle en Vercors, France
| | - Quentin Boesch
- Ecole et Observatoire des Sciences de la Terre, Institut de Physique du Globe de Strasbourg, UMR 7516 CNRS, Université de Strasbourg, Strasbourg, France
| | - Pierre-Olivier Antoine
- Institut des Sciences de l'Évolution de Montpellier, Univ Montpellier, CNRS, IRD, Montpellier, France
| | | | | | - Daovee Sihanam
- Ministry of Information, Culture and Tourism, Vientiane, Lao PDR
| | | | | | - Denis Fiorillo
- UMR 7209 Archéozoologie, Archéobotanique, Sociétés, Pratiques, Environnements, MNHN, CNRS, Paris, France
| | - Olivier Tombret
- UMR 7209 Archéozoologie, Archéobotanique, Sociétés, Pratiques, Environnements, MNHN, CNRS, Paris, France
| | - Elise Patole-Edoumba
- Muséum d'histoire naturelle de La Rochelle, UMRU 24140 Dynamiques, interactions, interculturalité asiatiques (UBM, LRUniv), La Rochelle, France
| | - Alexandra Zachwieja
- Department of Biomedical Sciences, University of Minnesota Medical School Duluth, Duluth, USA
| | | | | | - Tyler E Dunn
- Anatomical Sciences Education Center, Oregon Health & Sciences University, Portland, OR, USA
| | - Laura Shackelford
- Department of Anthropology, University of Illinois at Urbana-Champaign, Urbana, IL, USA
- Carle Illinois College of Medicine, University of Illinois at Urbana-Champaign, Urbana, IL, USA
| | - Jean-Jacques Hublin
- Max Planck Institute for Evolutionary Anthropology, Department of Human Evolution, 04103, Leipzig, Germany
- Chaire de Paléoanthropologie, CIRB (UMR 7241-U1050), Collège de France, Paris, France
| |
Collapse
|
23
|
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
|
24
|
Velazquez-Arcelay K, Colbran LL, McArthur E, Brand C, Rinker D, Siemann J, McMahon D, Capra JA. Archaic Introgression Shaped Human Circadian Traits. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.02.03.527061. [PMID: 36778254 PMCID: PMC9915721 DOI: 10.1101/2023.02.03.527061] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/09/2023]
Abstract
Introduction When the ancestors of modern Eurasians migrated out of Africa and interbred with Eurasian archaic hominins, namely Neanderthals and Denisovans, DNA of archaic ancestry integrated into the genomes of anatomically modern humans. This process potentially accelerated adaptation to Eurasian environmental factors, including reduced ultra-violet radiation and increased variation in seasonal dynamics. However, whether these groups differed substantially in circadian biology, and whether archaic introgression adaptively contributed to human chronotypes remains unknown. Results Here we traced the evolution of chronotype based on genomes from archaic hominins and present-day humans. First, we inferred differences in circadian gene sequences, splicing, and regulation between archaic hominins and modern humans. We identified 28 circadian genes containing variants with potential to alter splicing in archaics (e.g., CLOCK, PER2, RORB, RORC), and 16 circadian genes likely divergently regulated between present-day humans and archaic hominins, including RORA. These differences suggest the potential for introgression to modify circadian gene expression. Testing this hypothesis, we found that introgressed variants are enriched among eQTLs for circadian genes. Supporting the functional relevance of these regulatory effects, we found that many introgressed alleles have associations with chronotype. Strikingly, the strongest introgressed effects on chronotype increase morningness, consistent with adaptations to high latitude in other species. Finally, we identified several circadian loci with evidence of adaptive introgression or latitudinal clines in allele frequency. Conclusions These findings identify differences in circadian gene regulation between modern humans and archaic hominins and support the contribution of introgression via coordinated effects on variation in human chronotype.
Collapse
Affiliation(s)
| | - Laura L. Colbran
- Department of Genetics, Perelman School of Medicine, University of Pennsylvania
| | | | - Colin Brand
- Department of Epidemiology and Biostatistics, University of California, San Francisco
- Bakar Computational Health Sciences Institute, University of California, San Francisco
| | - David Rinker
- Department of Biological Sciences, Vanderbilt University
| | - Justin Siemann
- Department of Biological Sciences, Vanderbilt University
| | | | - John A. Capra
- Department of Epidemiology and Biostatistics, University of California, San Francisco
- Bakar Computational Health Sciences Institute, University of California, San Francisco
| |
Collapse
|
25
|
Roca-Umbert A, Garcia-Calleja J, Vogel-González M, Fierro-Villegas A, Ill-Raga G, Herrera-Fernández V, Bosnjak A, Muntané G, Gutiérrez E, Campelo F, Vicente R, Bosch E. Human genetic adaptation related to cellular zinc homeostasis. PLoS Genet 2023; 19:e1010950. [PMID: 37747921 PMCID: PMC10553801 DOI: 10.1371/journal.pgen.1010950] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2023] [Revised: 10/05/2023] [Accepted: 08/31/2023] [Indexed: 09/27/2023] Open
Abstract
SLC30A9 encodes a ubiquitously zinc transporter (ZnT9) and has been consistently suggested as a candidate for positive selection in humans. However, no direct adaptive molecular phenotype has been demonstrated. Our results provide evidence for directional selection operating in two major complementary haplotypes in Africa and East Asia. These haplotypes are associated with differential gene expression but also differ in the Met50Val substitution (rs1047626) in ZnT9, which we show is found in homozygosis in the Denisovan genome and displays accompanying signatures suggestive of archaic introgression. Although we found no significant differences in systemic zinc content between individuals with different rs1047626 genotypes, we demonstrate that the expression of the derived isoform (ZnT9 50Val) in HEK293 cells shows a gain of function when compared with the ancestral (ZnT9 50Met) variant. Notably, the ZnT9 50Val variant was found associated with differences in zinc handling by the mitochondria and endoplasmic reticulum, with an impact on mitochondrial metabolism. Given the essential role of the mitochondria in skeletal muscle and since the derived allele at rs1047626 is known to be associated with greater susceptibility to several neuropsychiatric traits, we propose that adaptation to cold may have driven this selection event, while also impacting predisposition to neuropsychiatric disorders in modern humans.
Collapse
Affiliation(s)
- Ana Roca-Umbert
- Institut de Biologia Evolutiva (UPF-CSIC), Departament de Medicina i Ciències de la Vida, Universitat Pompeu Fabra, Parc de Recerca Biomèdica de Barcelona, Barcelona, Spain
| | - Jorge Garcia-Calleja
- Institut de Biologia Evolutiva (UPF-CSIC), Departament de Medicina i Ciències de la Vida, Universitat Pompeu Fabra, Parc de Recerca Biomèdica de Barcelona, Barcelona, Spain
| | - Marina Vogel-González
- Laboratory of Molecular Physiology, Department of Medicine and Life Sciences (MELIS), Universitat Pompeu Fabra, Barcelona, Spain
| | - Alejandro Fierro-Villegas
- Laboratory of Molecular Physiology, Department of Medicine and Life Sciences (MELIS), Universitat Pompeu Fabra, Barcelona, Spain
| | - Gerard Ill-Raga
- Laboratory of Molecular Physiology, Department of Medicine and Life Sciences (MELIS), Universitat Pompeu Fabra, Barcelona, Spain
| | - Víctor Herrera-Fernández
- Laboratory of Molecular Physiology, Department of Medicine and Life Sciences (MELIS), Universitat Pompeu Fabra, Barcelona, Spain
| | - Anja Bosnjak
- Laboratory of Molecular Physiology, Department of Medicine and Life Sciences (MELIS), Universitat Pompeu Fabra, Barcelona, Spain
| | - Gerard Muntané
- Institut de Biologia Evolutiva (UPF-CSIC), Departament de Medicina i Ciències de la Vida, Universitat Pompeu Fabra, Parc de Recerca Biomèdica de Barcelona, Barcelona, Spain
- Hospital Universitari Institut Pere Mata, IISPV, Universitat Rovira i Virgili, Reus, Spain
- Centro de Investigación Biomédica en Red de Salud Mental, Instituto de Salud Carlos III, Madrid, Spain
| | - Esteban Gutiérrez
- Laboratory of Molecular Physiology, Department of Medicine and Life Sciences (MELIS), Universitat Pompeu Fabra, Barcelona, Spain
| | - Felix Campelo
- ICFO-Institut de Ciencies Fotoniques, The Barcelona Institute of Science and Technology, Barcelona, Spain
| | - Rubén Vicente
- Laboratory of Molecular Physiology, Department of Medicine and Life Sciences (MELIS), Universitat Pompeu Fabra, Barcelona, Spain
| | - Elena Bosch
- Institut de Biologia Evolutiva (UPF-CSIC), Departament de Medicina i Ciències de la Vida, Universitat Pompeu Fabra, Parc de Recerca Biomèdica de Barcelona, Barcelona, Spain
- Centro de Investigación Biomédica en Red de Salud Mental, Instituto de Salud Carlos III, Madrid, Spain
| |
Collapse
|
26
|
Gautney JR. A new approach to exploratory data analysis in hominin phylogenetic reconstruction. J Hum Evol 2023; 182:103412. [PMID: 37499423 DOI: 10.1016/j.jhevol.2023.103412] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/23/2022] [Revised: 06/30/2023] [Accepted: 07/01/2023] [Indexed: 07/29/2023]
Abstract
The phylogenetic relationships between fossil hominin taxa have been a contentious topic for decades. Recent discoveries of new taxa, rather than resolving the issue, have only further confused it. Compounding this problem are the limitations of some of the tools frequently used by paleoanthropologists to analyze these relationships. Most commonly, phylogenetic questions are investigated using analytical methods such as maximum parsimony and Bayesian analysis. While these are useful analytical tools, these tree-building methods can have limitations when investigating taxa that may have complex evolutionary histories. Exploratory data analysis can provide information about patterns in a dataset that are obscured by tree-based methods. These patterns include phylogenetic signal conflict, which is not depicted in tree-based methods. Signal conflict can have a number of sources, including methodological issues with character choice, taxonomic issues, homoplasy, and gene flow between taxa. In this study, an exploratory data analysis of fossil hominin morphological data is conducted using the tree-based analytical method neighbor-joining and the network-based analytical method neighbor-net with the goal of visualizing phylogenetic signal conflict within a hominin morphological data set. The data set is divided into cranial regions, and each cranial region is analyzed individually to investigate which regions of the skull contain the highest levels of signal conflict. Results of this analysis show that conflicting phylogenetic signals are present in the hominin fossil record during the relatively speciose period between 3 and 1 Ma, and they also indicate that levels of signal conflict vary by cranial region. Possible sources of these conflicting signals are then explored. Exploratory data analyses such as this can be a useful tool in generating phylogenetic hypotheses and in refining character choice. This study also highlights the value network-based approaches can bring to the hominin phylogenetic analysis toolkit.
Collapse
Affiliation(s)
- Joanna R Gautney
- Department of Sociology and Anthropology, Weber State University, 1299 Edvalson St., Ogden, UT, USA.
| |
Collapse
|
27
|
Aquino Y, Bisiaux A, Li Z, O'Neill M, Mendoza-Revilla J, Merkling SH, Kerner G, Hasan M, Libri V, Bondet V, Smith N, de Cevins C, Ménager M, Luca F, Pique-Regi R, Barba-Spaeth G, Pietropaoli S, Schwartz O, Leroux-Roels G, Lee CK, Leung K, Wu JT, Peiris M, Bruzzone R, Abel L, Casanova JL, Valkenburg SA, Duffy D, Patin E, Rotival M, Quintana-Murci L. Dissecting human population variation in single-cell responses to SARS-CoV-2. Nature 2023; 621:120-128. [PMID: 37558883 PMCID: PMC10482701 DOI: 10.1038/s41586-023-06422-9] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2022] [Accepted: 07/11/2023] [Indexed: 08/11/2023]
Abstract
Humans display substantial interindividual clinical variability after SARS-CoV-2 infection1-3, the genetic and immunological basis of which has begun to be deciphered4. However, the extent and drivers of population differences in immune responses to SARS-CoV-2 remain unclear. Here we report single-cell RNA-sequencing data for peripheral blood mononuclear cells-from 222 healthy donors of diverse ancestries-that were stimulated with SARS-CoV-2 or influenza A virus. We show that SARS-CoV-2 induces weaker, but more heterogeneous, interferon-stimulated gene activity compared with influenza A virus, and a unique pro-inflammatory signature in myeloid cells. Transcriptional responses to viruses display marked population differences, primarily driven by changes in cell abundance including increased lymphoid differentiation associated with latent cytomegalovirus infection. Expression quantitative trait loci and mediation analyses reveal a broad effect of cell composition on population disparities in immune responses, with genetic variants exerting a strong effect on specific loci. Furthermore, we show that natural selection has increased population differences in immune responses, particularly for variants associated with SARS-CoV-2 response in East Asians, and document the cellular and molecular mechanisms by which Neanderthal introgression has altered immune functions, such as the response of myeloid cells to viruses. Finally, colocalization and transcriptome-wide association analyses reveal an overlap between the genetic basis of immune responses to SARS-CoV-2 and COVID-19 severity, providing insights into the factors contributing to current disparities in COVID-19 risk.
Collapse
Affiliation(s)
- Yann Aquino
- Human Evolutionary Genetics Unit, Institut Pasteur, Université Paris Cité, CNRS UMR2000, Paris, France
- Collège Doctoral, Sorbonne Université, Paris, France
| | - Aurélie Bisiaux
- Human Evolutionary Genetics Unit, Institut Pasteur, Université Paris Cité, CNRS UMR2000, Paris, France
| | - Zhi Li
- Human Evolutionary Genetics Unit, Institut Pasteur, Université Paris Cité, CNRS UMR2000, Paris, France
| | - Mary O'Neill
- Human Evolutionary Genetics Unit, Institut Pasteur, Université Paris Cité, CNRS UMR2000, Paris, France
| | - Javier Mendoza-Revilla
- Human Evolutionary Genetics Unit, Institut Pasteur, Université Paris Cité, CNRS UMR2000, Paris, France
| | - Sarah Hélène Merkling
- Insect-Virus Interactions Unit, Institut Pasteur, Université Paris Cité, CNRS UMR2000, Paris, France
| | - Gaspard Kerner
- Human Evolutionary Genetics Unit, Institut Pasteur, Université Paris Cité, CNRS UMR2000, Paris, France
| | - Milena Hasan
- Cytometry and Biomarkers UTechS, Institut Pasteur, Université Paris Cité, Paris, France
| | - Valentina Libri
- Cytometry and Biomarkers UTechS, Institut Pasteur, Université Paris Cité, Paris, France
| | - Vincent Bondet
- Translational Immunology Unit, Institut Pasteur, Université Paris Cité, Paris, France
| | - Nikaïa Smith
- Translational Immunology Unit, Institut Pasteur, Université Paris Cité, Paris, France
| | - Camille de Cevins
- Université Paris Cité, Imagine Institute, Laboratory of Inflammatory Responses and Transcriptomic Networks in Diseases, Atip-Avenir Team, INSERM UMR1163, Paris, France
| | - Mickaël Ménager
- Université Paris Cité, Imagine Institute, Laboratory of Inflammatory Responses and Transcriptomic Networks in Diseases, Atip-Avenir Team, INSERM UMR1163, Paris, France
- Labtech Single-Cell@Imagine, Imagine Institute, INSERM UMR1163, Paris, France
| | - Francesca Luca
- Center for Molecular Medicine and Genetics, Wayne State University, Detroit, MI, USA
- Department of Obstetrics and Gynecology, Wayne State University, Detroit, MI, USA
- Department of Biology, University of Rome Tor Vergata, Rome, Italy
| | - Roger Pique-Regi
- Center for Molecular Medicine and Genetics, Wayne State University, Detroit, MI, USA
- Department of Obstetrics and Gynecology, Wayne State University, Detroit, MI, USA
| | - Giovanna Barba-Spaeth
- Structural Virology Unit, Institut Pasteur, Université Paris Cité, CNRS UMR3569, Paris, France
| | - Stefano Pietropaoli
- Structural Virology Unit, Institut Pasteur, Université Paris Cité, CNRS UMR3569, Paris, France
| | - Olivier Schwartz
- Virus and Immunity Unit, Institut Pasteur, Université Paris Cité, CNRS UMR3569, Paris, France
| | | | - Cheuk-Kwong Lee
- Hong Kong Red Cross Blood Transfusion Service, Hospital Authority, Hong Kong SAR, China
| | - Kathy Leung
- WHO Collaborating Centre for Infectious Disease Epidemiology and Control, School of Public Health, Li Ka Shing Faculty of Medicine, The University of Hong Kong, Hong Kong SAR, China
- Laboratory of Data Discovery for Health (D24H), Hong Kong Science Park, Hong Kong SAR, China
| | - Joseph T Wu
- WHO Collaborating Centre for Infectious Disease Epidemiology and Control, School of Public Health, Li Ka Shing Faculty of Medicine, The University of Hong Kong, Hong Kong SAR, China
- Laboratory of Data Discovery for Health (D24H), Hong Kong Science Park, Hong Kong SAR, China
| | - Malik Peiris
- Division of Public Health Laboratory Sciences, School of Public Health, Li Ka Shing Faculty of Medicine, The University of Hong Kong, Hong Kong SAR, China
- HKU-Pasteur Research Pole, School of Public Health, The University of Hong Kong, Hong Kong SAR, China
- Centre for Immunology and Infection, Hong Kong Science Park, Hong Kong SAR, China
| | - Roberto Bruzzone
- HKU-Pasteur Research Pole, School of Public Health, The University of Hong Kong, Hong Kong SAR, China
- Centre for Immunology and Infection, Hong Kong Science Park, Hong Kong SAR, China
| | - Laurent Abel
- St Giles Laboratory of Human Genetics of Infectious Diseases, The Rockefeller University, New York, NY, USA
- Laboratory of Human Genetics of Infectious Diseases, INSERM UMR1163, Necker Hospital for Sick Children, Paris, France
- Université Paris Cité, Imagine Institute, Paris, France
| | - Jean-Laurent Casanova
- St Giles Laboratory of Human Genetics of Infectious Diseases, The Rockefeller University, New York, NY, USA
- Laboratory of Human Genetics of Infectious Diseases, INSERM UMR1163, Necker Hospital for Sick Children, Paris, France
- Université Paris Cité, Imagine Institute, Paris, France
- Department of Pediatrics, Necker Hospital for Sick Children, Paris, France
- Howard Hughes Medical Institute, New York, NY, USA
| | - Sophie A Valkenburg
- HKU-Pasteur Research Pole, School of Public Health, The University of Hong Kong, Hong Kong SAR, China
- Department of Microbiology and Immunology, Peter Doherty Institute for Infection and Immunity, University of Melbourne, Melbourne, Victoria, Australia
| | - Darragh Duffy
- Translational Immunology Unit, Institut Pasteur, Université Paris Cité, Paris, France
- Centre for Immunology and Infection, Hong Kong Science Park, Hong Kong SAR, China
| | - Etienne Patin
- Human Evolutionary Genetics Unit, Institut Pasteur, Université Paris Cité, CNRS UMR2000, Paris, France
| | - Maxime Rotival
- Human Evolutionary Genetics Unit, Institut Pasteur, Université Paris Cité, CNRS UMR2000, Paris, France.
| | - Lluis Quintana-Murci
- Human Evolutionary Genetics Unit, Institut Pasteur, Université Paris Cité, CNRS UMR2000, Paris, France.
- Chair Human Genomics and Evolution, Collège de France, Paris, France.
| |
Collapse
|
28
|
Ruan J, Timmermann A, Raia P, Yun KS, Zeller E, Mondanaro A, Di Febbraro M, Lemmon D, Castiglione S, Melchionna M. Climate shifts orchestrated hominin interbreeding events across Eurasia. Science 2023; 381:699-704. [PMID: 37561879 DOI: 10.1126/science.add4459] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/29/2022] [Accepted: 04/19/2023] [Indexed: 08/12/2023]
Abstract
When, where, and how often hominin interbreeding happened is largely unknown. We study the potential for Neanderthal-Denisovan admixture using species distribution models that integrate extensive fossil, archaeological, and genetic data with transient coupled general circulation model simulations of global climate and biomes. Our Pleistocene hindcast of past hominins' habitat suitability reveals pronounced climate-driven zonal shifts in the main overlap region of Denisovans and Neanderthals in central Eurasia. These shifts, which influenced the timing and intensity of potential interbreeding events, can be attributed to the response of climate and vegetation to past variations in atmospheric carbon dioxide and Northern Hemisphere ice-sheet volume. Therefore, glacial-interglacial climate swings likely played an important role in favoring gene flow between archaic humans.
Collapse
Affiliation(s)
- Jiaoyang Ruan
- Center for Climate Physics, Institute for Basic Science, Busan, South Korea
- Center for Climate Physics, Pusan National University, Busan, South Korea
| | - Axel Timmermann
- Center for Climate Physics, Institute for Basic Science, Busan, South Korea
- Center for Climate Physics, Pusan National University, Busan, South Korea
| | - Pasquale Raia
- DiSTAR, Monte Sant'Angelo, Napoli Università di Napoli Federico II, Naples, Italy
| | - Kyung-Sook Yun
- Center for Climate Physics, Institute for Basic Science, Busan, South Korea
- Center for Climate Physics, Pusan National University, Busan, South Korea
| | - Elke Zeller
- Center for Climate Physics, Institute for Basic Science, Busan, South Korea
- Department of Climate System, Pusan National University, Busan, South Korea
| | | | - Mirko Di Febbraro
- Department of Biosciences and Territory, University of Molise, C. da Fonte Lappone, Pesche, Italy
| | - Danielle Lemmon
- Center for Climate Physics, Institute for Basic Science, Busan, South Korea
- Center for Climate Physics, Pusan National University, Busan, South Korea
| | - Silvia Castiglione
- DiSTAR, Monte Sant'Angelo, Napoli Università di Napoli Federico II, Naples, Italy
| | - Marina Melchionna
- DiSTAR, Monte Sant'Angelo, Napoli Università di Napoli Federico II, Naples, Italy
| |
Collapse
|
29
|
Maasch JRMA, Torres MDT, Melo MCR, de la Fuente-Nunez C. Molecular de-extinction of ancient antimicrobial peptides enabled by machine learning. Cell Host Microbe 2023; 31:1260-1274.e6. [PMID: 37516110 DOI: 10.1016/j.chom.2023.07.001] [Citation(s) in RCA: 17] [Impact Index Per Article: 17.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2022] [Revised: 05/12/2023] [Accepted: 07/06/2023] [Indexed: 07/31/2023]
Abstract
Molecular de-extinction could offer avenues for drug discovery by reintroducing bioactive molecules that are no longer encoded by extant organisms. To prospect for antimicrobial peptides encrypted within extinct and extant human proteins, we introduce the panCleave random forest model for proteome-wide cleavage site prediction. Our model outperformed multiple protease-specific cleavage site classifiers for three modern human caspases, despite its pan-protease design. Antimicrobial activity was observed in vitro for modern and archaic protein fragments identified with panCleave. Lead peptides showed resistance to proteolysis and exhibited variable membrane permeabilization. Additionally, representative modern and archaic protein fragments showed anti-infective efficacy against A. baumannii in both a skin abscess infection model and a preclinical murine thigh infection model. These results suggest that machine-learning-based encrypted peptide prospection can identify stable, nontoxic peptide antibiotics. Moreover, we establish molecular de-extinction through paleoproteome mining as a framework for antibacterial drug discovery.
Collapse
Affiliation(s)
- Jacqueline R M A Maasch
- Department of Computer and Information Science, School of Engineering and Applied Science, University of Pennsylvania, Philadelphia, PA 19104, USA; Machine Biology Group, Departments of Psychiatry and Microbiology, Institute for Biomedical Informatics, Institute for Translational Medicine and Therapeutics, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA 19104, USA; Department of Bioengineering, Department of Chemical and Biomolecular Engineering, School of Engineering and Applied Science, University of Pennsylvania, Philadelphia, PA 19104, USA; Penn Institute for Computational Science, University of Pennsylvania, Philadelphia, PA 19104, USA
| | - Marcelo D T Torres
- Machine Biology Group, Departments of Psychiatry and Microbiology, Institute for Biomedical Informatics, Institute for Translational Medicine and Therapeutics, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA 19104, USA; Department of Bioengineering, Department of Chemical and Biomolecular Engineering, School of Engineering and Applied Science, University of Pennsylvania, Philadelphia, PA 19104, USA; Penn Institute for Computational Science, University of Pennsylvania, Philadelphia, PA 19104, USA
| | - Marcelo C R Melo
- Machine Biology Group, Departments of Psychiatry and Microbiology, Institute for Biomedical Informatics, Institute for Translational Medicine and Therapeutics, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA 19104, USA; Department of Bioengineering, Department of Chemical and Biomolecular Engineering, School of Engineering and Applied Science, University of Pennsylvania, Philadelphia, PA 19104, USA; Penn Institute for Computational Science, University of Pennsylvania, Philadelphia, PA 19104, USA
| | - Cesar de la Fuente-Nunez
- Machine Biology Group, Departments of Psychiatry and Microbiology, Institute for Biomedical Informatics, Institute for Translational Medicine and Therapeutics, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA 19104, USA; Department of Bioengineering, Department of Chemical and Biomolecular Engineering, School of Engineering and Applied Science, University of Pennsylvania, Philadelphia, PA 19104, USA; Penn Institute for Computational Science, University of Pennsylvania, Philadelphia, PA 19104, USA.
| |
Collapse
|
30
|
Graves JL. Favored Races in the Struggle for Life: Racism and the Speciation Concept. Cold Spring Harb Perspect Biol 2023; 15:a041454. [PMID: 37463717 PMCID: PMC10411861 DOI: 10.1101/cshperspect.a041454] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/20/2023]
Abstract
Evolutionary speciation, whether it be cladistic or phyletic, has always been associated with race concepts. Biological races are conceived as definable stages of divergence from a common ancestor. However, the species concept in Western science began within a special creationist framework. The sixteenth century European voyages of discovery resulted in special creationist schemes explaining the origin of the new peoples encountered. These were designed to provide the moral justification for their colonization and enslavement. By the seventeenth century, European naturalists were beginning to seriously question the meaning of the variation within the animals and plants they observed within the context of God's role in creation. By the middle of the nineteenth century, "the species question" was the most important intellectual enterprise within biology. Here I discuss how notions of speciation influenced and were influenced by conceptions of race within Homo sapiens.
Collapse
Affiliation(s)
- Joseph L Graves
- Department of Biological Sciences, North Carolina A&T State University, Greensboro, North Carolina 27410, USA
| |
Collapse
|
31
|
Kun E, Javan EM, Smith O, Gulamali F, de la Fuente J, Flynn BI, Vajrala K, Trutner Z, Jayakumar P, Tucker-Drob EM, Sohail M, Singh T, Narasimhan VM. The genetic architecture and evolution of the human skeletal form. Science 2023; 381:eadf8009. [PMID: 37471560 PMCID: PMC11075689 DOI: 10.1126/science.adf8009] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2022] [Accepted: 06/20/2023] [Indexed: 07/22/2023]
Abstract
The human skeletal form underlies bipedalism, but the genetic basis of skeletal proportions (SPs) is not well characterized. We applied deep-learning models to 31,221 x-rays from the UK Biobank to extract a comprehensive set of SPs, which were associated with 145 independent loci genome-wide. Structural equation modeling suggested that limb proportions exhibited strong genetic sharing but were independent of width and torso proportions. Polygenic score analysis identified specific associations between osteoarthritis and hip and knee SPs. In contrast to other traits, SP loci were enriched in human accelerated regions and in regulatory elements of genes that are differentially expressed between humans and great apes. Combined, our work identifies specific genetic variants that affect the skeletal form and ties a major evolutionary facet of human anatomical change to pathogenesis.
Collapse
Affiliation(s)
- Eucharist Kun
- Department of Integrative Biology, The University of Texas at Austin, Austin, TX, USA
| | - Emily M. Javan
- Department of Integrative Biology, The University of Texas at Austin, Austin, TX, USA
| | - Olivia Smith
- Department of Integrative Biology, The University of Texas at Austin, Austin, TX, USA
| | - Faris Gulamali
- Icahn School of Medicine at Mount Sinai, New York, NY, USA
| | - Javier de la Fuente
- Department of Psychology, The University of Texas at Austin, Austin, TX, USA
| | - Brianna I. Flynn
- Department of Integrative Biology, The University of Texas at Austin, Austin, TX, USA
| | - Kushal Vajrala
- Department of Integrative Biology, The University of Texas at Austin, Austin, TX, USA
| | - Zoe Trutner
- Department of Surgery and Perioperative Care, The University of Texas at Austin, Austin, TX, USA
| | - Prakash Jayakumar
- Department of Surgery and Perioperative Care, The University of Texas at Austin, Austin, TX, USA
| | | | - Mashaal Sohail
- Centro de Ciencias Genómicas (CCG), Universidad Nacional Autónoma de México (UNAM), 62209 Cuernavaca, Mexico
| | - Tarjinder Singh
- Department of Psychiatry, Columbia University Irving Medical Center, New York, NY, USA
- The New York Genome Center, New York, NY, USA
- Mortimer B. Zuckerman Mind Brain Behavior Institute at Columbia University, New York, NY, USA
| | - Vagheesh M. Narasimhan
- Department of Integrative Biology, The University of Texas at Austin, Austin, TX, USA
- Department of Statistics and Data Science, The University of Texas at Austin, Austin, TX, USA
| |
Collapse
|
32
|
Purcell RA, Theisen RM, Arnold KB, Chung AW, Selva KJ. Polyfunctional antibodies: a path towards precision vaccines for vulnerable populations. Front Immunol 2023; 14:1183727. [PMID: 37600816 PMCID: PMC10433199 DOI: 10.3389/fimmu.2023.1183727] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2023] [Accepted: 05/30/2023] [Indexed: 08/22/2023] Open
Abstract
Vaccine efficacy determined within the controlled environment of a clinical trial is usually substantially greater than real-world vaccine effectiveness. Typically, this results from reduced protection of immunologically vulnerable populations, such as children, elderly individuals and people with chronic comorbidities. Consequently, these high-risk groups are frequently recommended tailored immunisation schedules to boost responses. In addition, diverse groups of healthy adults may also be variably protected by the same vaccine regimen. Current population-based vaccination strategies that consider basic clinical parameters offer a glimpse into what may be achievable if more nuanced aspects of the immune response are considered in vaccine design. To date, vaccine development has been largely empirical. However, next-generation approaches require more rational strategies. We foresee a generation of precision vaccines that consider the mechanistic basis of vaccine response variations associated with both immunogenetic and baseline health differences. Recent efforts have highlighted the importance of balanced and diverse extra-neutralising antibody functions for vaccine-induced protection. However, in immunologically vulnerable populations, significant modulation of polyfunctional antibody responses that mediate both neutralisation and effector functions has been observed. Here, we review the current understanding of key genetic and inflammatory modulators of antibody polyfunctionality that affect vaccination outcomes and consider how this knowledge may be harnessed to tailor vaccine design for improved public health.
Collapse
Affiliation(s)
- Ruth A. Purcell
- Department of Microbiology and Immunology, Peter Doherty Institute for Infection and Immunity, University of Melbourne, Melbourne, VIC, Australia
| | - Robert M. Theisen
- Department of Biomedical Engineering, University of Michigan, Ann Arbor, MI, United States
| | - Kelly B. Arnold
- Department of Biomedical Engineering, University of Michigan, Ann Arbor, MI, United States
| | - Amy W. Chung
- Department of Microbiology and Immunology, Peter Doherty Institute for Infection and Immunity, University of Melbourne, Melbourne, VIC, Australia
| | - Kevin J. Selva
- Department of Microbiology and Immunology, Peter Doherty Institute for Infection and Immunity, University of Melbourne, Melbourne, VIC, Australia
| |
Collapse
|
33
|
Wang J, Wu J, Sun Q, Wu Q, Li Y, Duan S, Yang L, Wu W, Wang Z, Liu Y, Tang R, Yang J, Wang C, Liu C, Xu J, Wang M, He G. Extensive genetic admixture between Tai-Kadai-speaking people and their neighbours in the northeastern region of the Yungui Plateau inferred from genome-wide variations. BMC Genomics 2023; 24:317. [PMID: 37308851 DOI: 10.1186/s12864-023-09412-3] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2022] [Accepted: 05/27/2023] [Indexed: 06/14/2023] Open
Abstract
BACKGROUND Yungui Plateau in Southwest China is characterized by multi-language and multi-ethnic communities and is one of the regions with the wealthiest ethnolinguistic, cultural and genetic diversity in East Asia. There are numerous Tai-Kadai (TK)-speaking populations, but their detailed evolutionary history and biological adaptations are still unclear. RESULTS Here, we genotyped genome-wide SNP data of 77 unrelated TK-speaking Zhuang and Dong individuals from the Yungui Plateau and explored their detailed admixture history and adaptive features using clustering patterns, allele frequency differentiation and sharing haplotype patterns. TK-speaking Zhuang and Dong people in Guizhou are closely related to geographically close TK and Hmong-Mien (HM)-speaking populations. Besides, we identified that Guizhou TK-speaking people have a close genetic relationship with Austronesian (AN)-speaking Atayal and Paiwan people, which is supported by the common origin of the ancient Baiyue tribe. We additionally found subtle genetic differences among the newly studied TK people and previously reported Dais via the fine-scale genetic substructure analysis based on the shared haplotype chunks. Finally, we identified specific selection candidate signatures associated with several essential human immune systems and neurological disorders, which could provide evolutionary evidence for the allele frequency distribution pattern of genetic risk loci. CONCLUSIONS Our comprehensive genetic characterization of TK people suggested the strong genetic affinity within TK groups and extensive gene flow with geographically close HM and Han people. We also provided genetic evidence that supported the common origin hypothesis of TK and AN people. The best-fitted admixture models further suggested that ancestral sources from northern millet farmers and southern inland and coastal people contributed to the formation of the gene pool of the Zhuang and Dong people.
Collapse
Affiliation(s)
- Jiawen Wang
- School of Forensic Medicine, Guizhou Medical University, Guiyang, 550004, China.
| | - Jun Wu
- School of Forensic Medicine, Guizhou Medical University, Guiyang, 550004, China
| | - Qiuxia Sun
- Institute of Rare Diseases, West China Hospital of Sichuan University, Sichuan University, Chengdu, 610000, China
- Department of Forensic Medicine, College of Basic Medicine, Chongqing Medical University, Chongqing, 400331, China
| | - Qian Wu
- Qiannan Prefecture People's Hospital, Buyi and Miao Autonomous Prefecture of QianNan, Buyi and Miao Autonomous Prefecture of QianNan, 558000, China
| | - Youjing Li
- Congjiang People's Hospital, Congjiang, 557499, China
| | - Shuhan Duan
- Institute of Rare Diseases, West China Hospital of Sichuan University, Sichuan University, Chengdu, 610000, China
- School of Basic Medical Sciences, North Sichuan Medical College, Nanchong, 637000, China
| | - Lin Yang
- School of Forensic Medicine, Guizhou Medical University, Guiyang, 550004, China
| | - Wenxin Wu
- School of Forensic Medicine, Guizhou Medical University, Guiyang, 550004, China
| | - Zheng Wang
- Institute of Forensic Medicine, West China School of Basic Medical Sciences & Forensic Medicine, Sichuan University, Chengdu, 610041, China
| | - Yan Liu
- Institute of Rare Diseases, West China Hospital of Sichuan University, Sichuan University, Chengdu, 610000, China
- School of Basic Medical Sciences, North Sichuan Medical College, Nanchong, 637000, China
| | - Renkuan Tang
- Department of Forensic Medicine, College of Basic Medicine, Chongqing Medical University, Chongqing, 400331, China
| | - Junbao Yang
- School of Basic Medical Sciences, North Sichuan Medical College, Nanchong, 637000, China
| | - Chuanchao Wang
- Department of Anthropology and Ethnology, School of Sociology and Anthropology, Institute of Anthropology, Xiamen University, Xiamen, 361000, China
| | - Chao Liu
- Faculty of Forensic Medicine, Zhongshan School of Medicine, Sun Yat-Sen University, Guangzhou, 510275, China
| | - Jianwei Xu
- Department of Pharmacology, School of Basic Medicine, Guizhou Medical University, Guiyang, 550004, China.
| | - Mengge Wang
- Faculty of Forensic Medicine, Zhongshan School of Medicine, Sun Yat-Sen University, Guangzhou, 510275, China.
| | - Guanglin He
- Institute of Rare Diseases, West China Hospital of Sichuan University, Sichuan University, Chengdu, 610000, China.
- Center for Archaeological Science, Sichuan University, Chengdu, 610000, China.
| |
Collapse
|
34
|
Brand CM, Colbran LL, Capra JA. Resurrecting the alternative splicing landscape of archaic hominins using machine learning. Nat Ecol Evol 2023; 7:939-953. [PMID: 37142741 DOI: 10.1038/s41559-023-02053-5] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2022] [Accepted: 03/29/2023] [Indexed: 05/06/2023]
Abstract
Alternative splicing contributes to adaptation and divergence in many species. However, it has not been possible to directly compare splicing between modern and archaic hominins. Here, we unmask the recent evolution of this previously unobservable regulatory mechanism by applying SpliceAI, a machine-learning algorithm that identifies splice-altering variants (SAVs), to high-coverage genomes from three Neanderthals and a Denisovan. We discover 5,950 putative archaic SAVs, of which 2,186 are archaic-specific and 3,607 also occur in modern humans via introgression (244) or shared ancestry (3,520). Archaic-specific SAVs are enriched in genes that contribute to traits potentially relevant to hominin phenotypic divergence, such as the epidermis, respiration and spinal rigidity. Compared to shared SAVs, archaic-specific SAVs occur in sites under weaker selection and are more common in genes with tissue-specific expression. Further underscoring the importance of negative selection on SAVs, Neanderthal lineages with low effective population sizes are enriched for SAVs compared to Denisovan and shared SAVs. Finally, we find that nearly all introgressed SAVs in humans were shared across the three Neanderthals, suggesting that older SAVs were more tolerated in human genomes. Our results reveal the splicing landscape of archaic hominins and identify potential contributions of splicing to phenotypic differences among hominins.
Collapse
Affiliation(s)
- Colin M Brand
- Bakar Computational Health Sciences Institute, University of California, San Francisco, CA, USA
- Department of Epidemiology and Biostatistics, University of California, San Francisco, CA, USA
| | - Laura L Colbran
- Department of Genetics, Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| | - John A Capra
- Bakar Computational Health Sciences Institute, University of California, San Francisco, CA, USA.
- Department of Epidemiology and Biostatistics, University of California, San Francisco, CA, USA.
| |
Collapse
|
35
|
Rong S, Neil CR, Welch A, Duan C, Maguire S, Meremikwu IC, Meyerson M, Evans BJ, Fairbrother WG. Large-scale functional screen identifies genetic variants with splicing effects in modern and archaic humans. Proc Natl Acad Sci U S A 2023; 120:e2218308120. [PMID: 37192163 PMCID: PMC10214146 DOI: 10.1073/pnas.2218308120] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2022] [Accepted: 04/12/2023] [Indexed: 05/18/2023] Open
Abstract
Humans coexisted and interbred with other hominins which later became extinct. These archaic hominins are known to us only through fossil records and for two cases, genome sequences. Here, we engineer Neanderthal and Denisovan sequences into thousands of artificial genes to reconstruct the pre-mRNA processing patterns of these extinct populations. Of the 5,169 alleles tested in this massively parallel splicing reporter assay (MaPSy), we report 962 exonic splicing mutations that correspond to differences in exon recognition between extant and extinct hominins. Using MaPSy splicing variants, predicted splicing variants, and splicing quantitative trait loci, we show that splice-disrupting variants experienced greater purifying selection in anatomically modern humans than that in Neanderthals. Adaptively introgressed variants were enriched for moderate-effect splicing variants, consistent with positive selection for alternative spliced alleles following introgression. As particularly compelling examples, we characterized a unique tissue-specific alternative splicing variant at the adaptively introgressed innate immunity gene TLR1, as well as a unique Neanderthal introgressed alternative splicing variant in the gene HSPG2 that encodes perlecan. We further identified potentially pathogenic splicing variants found only in Neanderthals and Denisovans in genes related to sperm maturation and immunity. Finally, we found splicing variants that may contribute to variation among modern humans in total bilirubin, balding, hemoglobin levels, and lung capacity. Our findings provide unique insights into natural selection acting on splicing in human evolution and demonstrate how functional assays can be used to identify candidate causal variants underlying differences in gene regulation and phenotype.
Collapse
Affiliation(s)
- Stephen Rong
- Center for Computational Molecular Biology, Brown University, Providence, RI02912
- Department of Molecular Biology, Cell Biology, and Biochemistry, Brown University, Providence, RI02912
| | - Christopher R. Neil
- Department of Molecular Biology, Cell Biology, and Biochemistry, Brown University, Providence, RI02912
| | - Anastasia Welch
- Department of Molecular Biology, Cell Biology, and Biochemistry, Brown University, Providence, RI02912
| | - Chaorui Duan
- Department of Molecular Biology, Cell Biology, and Biochemistry, Brown University, Providence, RI02912
| | - Samantha Maguire
- Department of Molecular Biology, Cell Biology, and Biochemistry, Brown University, Providence, RI02912
| | - Ijeoma C. Meremikwu
- Department of Molecular Biology, Cell Biology, and Biochemistry, Brown University, Providence, RI02912
| | - Malcolm Meyerson
- Department of Molecular Biology, Cell Biology, and Biochemistry, Brown University, Providence, RI02912
| | - Ben J. Evans
- Department of Biology, McMaster University, Hamilton, ONL8S 4K1, Canada
| | - William G. Fairbrother
- Center for Computational Molecular Biology, Brown University, Providence, RI02912
- Department of Molecular Biology, Cell Biology, and Biochemistry, Brown University, Providence, RI02912
- Hassenfeld Child Health Innovation Institute of Brown University, Providence, RI02912
| |
Collapse
|
36
|
Witt KE, Funk A, Añorve-Garibay V, Fang LL, Huerta-Sánchez E. The Impact of Modern Admixture on Archaic Human Ancestry in Human Populations. Genome Biol Evol 2023; 15:evad066. [PMID: 37103242 PMCID: PMC10194819 DOI: 10.1093/gbe/evad066] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2022] [Revised: 03/07/2023] [Accepted: 04/17/2023] [Indexed: 04/28/2023] Open
Abstract
Admixture, the genetic merging of parental populations resulting in mixed ancestry, has occurred frequently throughout the course of human history. Numerous admixture events have occurred between human populations across the world, which have shaped genetic ancestry in modern humans. For example, populations in the Americas are often mosaics of different ancestries due to recent admixture events as part of European colonization. Admixed individuals also often have introgressed DNA from Neanderthals and Denisovans that may have come from multiple ancestral populations, which may affect how archaic ancestry is distributed across an admixed genome. In this study, we analyzed admixed populations from the Americas to assess whether the proportion and location of admixed segments due to recent admixture impact an individual's archaic ancestry. We identified a positive correlation between non-African ancestry and archaic alleles, as well as a slight increase of Denisovan alleles in Indigenous American segments relative to European segments in admixed genomes. We also identify several genes as candidates for adaptive introgression, based on archaic alleles present at high frequency in admixed American populations but low frequency in East Asian populations. These results provide insights into how recent admixture events between modern humans redistributed archaic ancestry in admixed genomes.
Collapse
Affiliation(s)
- Kelsey E Witt
- Ecology, Evolution, and Organismal Biology, Brown University, Providence, Rhode Island
- Center for Computational Molecular Biology, Brown University, Providence, Rhode Island
| | - Alyssa Funk
- Center for Computational Molecular Biology, Brown University, Providence, Rhode Island
- Molecular Biology, Cell Biology, & Biochemistry, Brown University, Providence, Rhode Island
| | - Valeria Añorve-Garibay
- Center for Computational Molecular Biology, Brown University, Providence, Rhode Island
- Licenciatura en Ciencias Genómicas, Escuela Nacional de Estudios Superiores Unidad Juriquilla, Universidad Nacional Autónoma de México, Querétaro, Mexico
- Laboratorio Internacional de Investigación sobre el Genoma Humano, Universidad Nacional Autónoma de México, Querétaro, Mexico
| | - Lesly Lopez Fang
- Department of Life & Environmental Sciences, University of California, Merced, California, United States of America
| | - Emilia Huerta-Sánchez
- Ecology, Evolution, and Organismal Biology, Brown University, Providence, Rhode Island
- Center for Computational Molecular Biology, Brown University, Providence, Rhode Island
| |
Collapse
|
37
|
Zheng W, He Y, Guo Y, Yue T, Zhang H, Li J, Zhou B, Zeng X, Li L, Wang B, Cao J, Chen L, Li C, Li H, Cui C, Bai C, Qi X, Su B. Large-scale genome sequencing redefines the genetic footprints of high-altitude adaptation in Tibetans. Genome Biol 2023; 24:73. [PMID: 37055782 PMCID: PMC10099689 DOI: 10.1186/s13059-023-02912-1] [Citation(s) in RCA: 12] [Impact Index Per Article: 12.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2022] [Accepted: 03/29/2023] [Indexed: 04/15/2023] Open
Abstract
BACKGROUND Tibetans are genetically adapted to high-altitude environments. Though many studies have been conducted, the genetic basis of the adaptation remains elusive due to the poor reproducibility for detecting selective signatures in the Tibetan genomes. RESULTS Here, we present whole-genome sequencing (WGS) data of 1001 indigenous Tibetans, covering the major populated areas of the Qinghai-Tibetan Plateau in China. We identify 35 million variants, and more than one-third of them are novel variants. Utilizing the large-scale WGS data, we construct a comprehensive map of allele frequency and linkage disequilibrium and provide a population-specific genome reference panel, referred to as 1KTGP. Moreover, with the use of a combined approach, we redefine the signatures of Darwinian-positive selection in the Tibetan genomes, and we characterize a high-confidence list of 4320 variants and 192 genes that have undergone selection in Tibetans. In particular, we discover four new genes, TMEM132C, ATP13A3, SANBR, and KHDRBS2, with strong signals of selection, and they may account for the adaptation of cardio-pulmonary functions in Tibetans. Functional annotation and enrichment analysis indicate that the 192 genes with selective signatures are likely involved in multiple organs and physiological systems, suggesting polygenic and pleiotropic effects. CONCLUSIONS Overall, the large-scale Tibetan WGS data and the identified adaptive variants/genes can serve as a valuable resource for future genetic and medical studies of high-altitude populations.
Collapse
Affiliation(s)
- Wangshan Zheng
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, 650223, China
- Kunming College of Life Science, University of Chinese Academy of Sciences, Beijing, 100101, China
| | - Yaoxi He
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, 650223, China
| | - Yongbo Guo
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, 650223, China
- Kunming College of Life Science, University of Chinese Academy of Sciences, Beijing, 100101, China
| | - Tian Yue
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, 650223, China
- Kunming College of Life Science, University of Chinese Academy of Sciences, Beijing, 100101, China
| | - Hui Zhang
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, 650223, China
| | - Jun Li
- Fukang Obstetrics, Gynecology and Children Branch Hospital, Tibetan Fukang Hospital, Lhasa, 850000, China
| | - Bin Zhou
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, 650223, China
- Kunming College of Life Science, University of Chinese Academy of Sciences, Beijing, 100101, China
| | - Xuerui Zeng
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, 650223, China
- Kunming College of Life Science, University of Chinese Academy of Sciences, Beijing, 100101, China
| | - Liya Li
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, 650223, China
| | - Bin Wang
- Fukang Obstetrics, Gynecology and Children Branch Hospital, Tibetan Fukang Hospital, Lhasa, 850000, China
| | - Jingxin Cao
- Fukang Obstetrics, Gynecology and Children Branch Hospital, Tibetan Fukang Hospital, Lhasa, 850000, China
| | - Li Chen
- Fukang Obstetrics, Gynecology and Children Branch Hospital, Tibetan Fukang Hospital, Lhasa, 850000, China
| | - Chunxia Li
- Fukang Obstetrics, Gynecology and Children Branch Hospital, Tibetan Fukang Hospital, Lhasa, 850000, China
| | - Hongyan Li
- Fukang Obstetrics, Gynecology and Children Branch Hospital, Tibetan Fukang Hospital, Lhasa, 850000, China
| | - Chaoying Cui
- High Altitude Medical Research Center, School of Medicine, Tibetan University, Lhasa, 850000, China
| | - Caijuan Bai
- High Altitude Medical Research Center, School of Medicine, Tibetan University, Lhasa, 850000, China
| | - Xuebin Qi
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, 650223, China.
- Fukang Obstetrics, Gynecology and Children Branch Hospital, Tibetan Fukang Hospital, Lhasa, 850000, China.
| | - Bing Su
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, 650223, China.
- Center for Excellence in Animal Evolution and Genetics, Chinese Academy of Sciences, Kunming, 650223, China.
| |
Collapse
|
38
|
Wei X, Robles CR, Pazokitoroudi A, Ganna A, Gusev A, Durvasula A, Gazal S, Loh PR, Reich D, Sankararaman S. The lingering effects of Neanderthal introgression on human complex traits. eLife 2023; 12:e80757. [PMID: 36939312 PMCID: PMC10076017 DOI: 10.7554/elife.80757] [Citation(s) in RCA: 7] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/03/2022] [Accepted: 03/17/2023] [Indexed: 03/21/2023] Open
Abstract
The genetic variants introduced into the ancestors of modern humans from interbreeding with Neanderthals have been suggested to contribute an unexpected extent to complex human traits. However, testing this hypothesis has been challenging due to the idiosyncratic population genetic properties of introgressed variants. We developed rigorous methods to assess the contribution of introgressed Neanderthal variants to heritable trait variation and applied these methods to analyze 235,592 introgressed Neanderthal variants and 96 distinct phenotypes measured in about 300,000 unrelated white British individuals in the UK Biobank. Introgressed Neanderthal variants make a significant contribution to trait variation (explaining 0.12% of trait variation on average). However, the contribution of introgressed variants tends to be significantly depleted relative to modern human variants matched for allele frequency and linkage disequilibrium (about 59% depletion on average), consistent with purifying selection on introgressed variants. Different from previous studies (McArthur et al., 2021), we find no evidence for elevated heritability across the phenotypes examined. We identified 348 independent significant associations of introgressed Neanderthal variants with 64 phenotypes. Previous work (Skov et al., 2020) has suggested that a majority of such associations are likely driven by statistical association with nearby modern human variants that are the true causal variants. Applying a customized fine-mapping led us to identify 112 regions across 47 phenotypes containing 4303 unique genetic variants where introgressed variants are highly likely to have a phenotypic effect. Examination of these variants reveals their substantial impact on genes that are important for the immune system, development, and metabolism.
Collapse
Affiliation(s)
- Xinzhu Wei
- Department of Computational Biology, Cornell UniversityNew YorkUnited States
| | - Christopher R Robles
- Department of Human Genetics, University of California, Los AngelesLos AngelesUnited States
| | - Ali Pazokitoroudi
- Department of Computer Science, University of California, Los AngelesLos AngelesUnited States
| | - Andrea Ganna
- Analytical and Translational Genetics Unit, Center for Genomic Medicine, Massachusetts General HospitalBostonUnited States
- Program in Medical and Population Genetics, Broad Institute of MIT and HarvardCambridgeUnited States
- Stanley Center for Psychiatric Research, Broad Institute of MIT and HarvardCambridgeUnited States
| | - Alexander Gusev
- Dana-Farber Cancer Institute, Harvard Medical SchoolBostonUnited States
| | - Arun Durvasula
- Department of Genetics, Harvard Medical SchoolBostonUnited States
- Department of Human Evolutionary Biology, Harvard UniversityCambridgeUnited States
| | - Steven Gazal
- Center for Genetic Epidemiology, Department of Public and Population Health Sciences, University of Southern CaliforniaLos AngelesUnited States
- Division of Genetics,Department of Medicine, Brigham and Women’s Hospital, Harvard Medical SchoolBostonUnited States
| | - Po-Ru Loh
- Program in Medical and Population Genetics, Broad Institute of MIT and HarvardCambridgeUnited States
| | - David Reich
- Program in Medical and Population Genetics, Broad Institute of MIT and HarvardCambridgeUnited States
- Department of Genetics, Harvard Medical SchoolBostonUnited States
- Department of Human Evolutionary Biology, Harvard UniversityCambridgeUnited States
- Howard Hughes Medical Institute, Harvard Medical SchoolBostonUnited States
| | - Sriram Sankararaman
- Department of Human Genetics, University of California, Los AngelesLos AngelesUnited States
- Department of Computer Science, University of California, Los AngelesLos AngelesUnited States
- Department of Computational Medicine, University of California, Los AngelesLos AngelesUnited States
| |
Collapse
|
39
|
Watanabe Y, Ohashi J. Modern Japanese ancestry-derived variants reveal the formation process of the current Japanese regional gradations. iScience 2023; 26:106130. [PMID: 36879818 PMCID: PMC9984562 DOI: 10.1016/j.isci.2023.106130] [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: 09/26/2022] [Revised: 12/02/2022] [Accepted: 01/31/2023] [Indexed: 02/05/2023] Open
Abstract
Modern Japanese people have two major ancestral populations: indigenous Jomon hunter-gatherers and continental East Asian farmers. To determine the formation process of the current Japanese population, we developed a detection method for variants derived from ancestral populations using a summary statistic, the ancestry marker index (AMI). We applied AMI to modern Japanese population samples and identified 208,648 single nucleotide polymorphisms (SNPs) that were likely derived from the Jomon people (Jomon-derived variants). Analysis of Jomon-derived variants in 10,842 modern Japanese individuals recruited from all over Japan revealed that the admixture proportions of the Jomon people varied between prefectures, probably owing to the prehistoric population size difference. The estimated allele frequencies of genome-wide SNPs in the ancestral populations of the modern Japanese suggested their adaptive phenotypic characteristics to their respective livelihoods. Based on our findings, we propose a formation model for the genotypic and phenotypic gradations of the current Japanese archipelago populations.
Collapse
Affiliation(s)
- Yusuke Watanabe
- Department of Biological Sciences, Graduate School of Science, The University of Tokyo, Tokyo 113-0033, Japan.,Genome Medical Science Project Toyama Project, National Center for Global Health and Medicine, Tokyo 162-8655, Japan
| | - Jun Ohashi
- Department of Biological Sciences, Graduate School of Science, The University of Tokyo, Tokyo 113-0033, Japan
| |
Collapse
|
40
|
Corcoran M, Chernyshev M, Mandolesi M, Narang S, Kaduk M, Ye K, Sundling C, Färnert A, Kreslavsky T, Bernhardsson C, Larena M, Jakobsson M, Karlsson Hedestam GB. Archaic humans have contributed to large-scale variation in modern human T cell receptor genes. Immunity 2023; 56:635-652.e6. [PMID: 36796364 DOI: 10.1016/j.immuni.2023.01.026] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/12/2022] [Revised: 11/21/2022] [Accepted: 01/23/2023] [Indexed: 02/18/2023]
Abstract
Human T cell receptors (TCRs) are critical for mediating immune responses to pathogens and tumors and regulating self-antigen recognition. Yet, variations in the genes encoding TCRs remain insufficiently defined. Detailed analysis of expressed TCR alpha, beta, gamma, and delta genes in 45 donors from four human populations-African, East Asian, South Asian, and European-revealed 175 additional TCR variable and junctional alleles. Most of these contained coding changes and were present at widely differing frequencies in the populations, a finding confirmed using DNA samples from the 1000 Genomes Project. Importantly, we identified three Neanderthal-derived, introgressed TCR regions including a highly divergent TRGV4 variant, which mediated altered butyrophilin-like molecule 3 (BTNL3) ligand reactivity and was frequent in all modern Eurasian population groups. Our results demonstrate remarkable variation in TCR genes in both individuals and populations, providing a strong incentive for including allelic variation in studies of TCR function in human biology.
Collapse
Affiliation(s)
- Martin Corcoran
- Department of Microbiology, Tumor and Cell Biology, Karolinska Institutet, 171 77 Stockholm, Sweden.
| | - Mark Chernyshev
- Department of Microbiology, Tumor and Cell Biology, Karolinska Institutet, 171 77 Stockholm, Sweden
| | - Marco Mandolesi
- Department of Microbiology, Tumor and Cell Biology, Karolinska Institutet, 171 77 Stockholm, Sweden
| | - Sanjana Narang
- Department of Microbiology, Tumor and Cell Biology, Karolinska Institutet, 171 77 Stockholm, Sweden
| | - Mateusz Kaduk
- Department of Microbiology, Tumor and Cell Biology, Karolinska Institutet, 171 77 Stockholm, Sweden
| | - Kewei Ye
- Department of Medicine, Solna, Karolinska Institutet and Karolinska University Hospital, Stockholm, Sweden; Center for Molecular Medicine, Karolinska Institutet, Stockholm, Sweden
| | - Christopher Sundling
- Department of Medicine, Solna, Karolinska Institutet and Karolinska University Hospital, Stockholm, Sweden; Center for Molecular Medicine, Karolinska Institutet, Stockholm, Sweden; Department of Infectious Diseases, Karolinska University Hospital, 171 76 Stockholm, Sweden
| | - Anna Färnert
- Department of Medicine, Solna, Karolinska Institutet and Karolinska University Hospital, Stockholm, Sweden; Center for Molecular Medicine, Karolinska Institutet, Stockholm, Sweden; Department of Infectious Diseases, Karolinska University Hospital, 171 76 Stockholm, Sweden
| | - Taras Kreslavsky
- Department of Medicine, Solna, Karolinska Institutet and Karolinska University Hospital, Stockholm, Sweden; Center for Molecular Medicine, Karolinska Institutet, Stockholm, Sweden
| | - Carolina Bernhardsson
- Department of Organismal Biology, Human Evolution, Norbyvägen 18C, 752 63 Uppsala, Sweden
| | - Maximilian Larena
- Department of Organismal Biology, Human Evolution, Norbyvägen 18C, 752 63 Uppsala, Sweden
| | - Mattias Jakobsson
- Department of Organismal Biology, Human Evolution, Norbyvägen 18C, 752 63 Uppsala, Sweden
| | | |
Collapse
|
41
|
Jagoda E, Marnetto D, Senevirathne G, Gonzalez V, Baid K, Montinaro F, Richard D, Falzarano D, LeBlanc EV, Colpitts CC, Banerjee A, Pagani L, Capellini TD. Regulatory dissection of the severe COVID-19 risk locus introgressed by Neanderthals. eLife 2023; 12:e71235. [PMID: 36763080 PMCID: PMC9917435 DOI: 10.7554/elife.71235] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/12/2021] [Accepted: 01/26/2023] [Indexed: 02/11/2023] Open
Abstract
Individuals infected with the SARS-CoV-2 virus present with a wide variety of symptoms ranging from asymptomatic to severe and even lethal outcomes. Past research has revealed a genetic haplotype on chromosome 3 that entered the human population via introgression from Neanderthals as the strongest genetic risk factor for the severe response to COVID-19. However, the specific variants along this introgressed haplotype that contribute to this risk and the biological mechanisms that are involved remain unclear. Here, we assess the variants present on the risk haplotype for their likelihood of driving the genetic predisposition to severe COVID-19 outcomes. We do this by first exploring their impact on the regulation of genes involved in COVID-19 infection using a variety of population genetics and functional genomics tools. We then perform a locus-specific massively parallel reporter assay to individually assess the regulatory potential of each allele on the haplotype in a multipotent immune-related cell line. We ultimately reduce the set of over 600 linked genetic variants to identify four introgressed alleles that are strong functional candidates for driving the association between this locus and severe COVID-19. Using reporter assays in the presence/absence of SARS-CoV-2, we find evidence that these variants respond to viral infection. These variants likely drive the locus' impact on severity by modulating the regulation of two critical chemokine receptor genes: CCR1 and CCR5. These alleles are ideal targets for future functional investigations into the interaction between host genomics and COVID-19 outcomes.
Collapse
Affiliation(s)
- Evelyn Jagoda
- Department of Human Evolutionary Biology, Harvard UniversityCambridgeUnited States
| | - Davide Marnetto
- Estonian Biocentre, Institute of Genomics, University of TartuTartuEstonia
| | - Gayani Senevirathne
- Department of Human Evolutionary Biology, Harvard UniversityCambridgeUnited States
| | - Victoria Gonzalez
- Department of Veterinary Microbiology, University of SaskatchewanSaskatoonCanada
- Vaccine and Infectious Disease Organization, University of SaskatchewanSaskatoonCanada
| | - Kaushal Baid
- Vaccine and Infectious Disease Organization, University of SaskatchewanSaskatoonCanada
| | - Francesco Montinaro
- Estonian Biocentre, Institute of Genomics, University of TartuTartuEstonia
- Department of Biology, University of BariBariItaly
| | - Daniel Richard
- Department of Human Evolutionary Biology, Harvard UniversityCambridgeUnited States
| | - Darryl Falzarano
- Department of Veterinary Microbiology, University of SaskatchewanSaskatoonCanada
- Vaccine and Infectious Disease Organization, University of SaskatchewanSaskatoonCanada
| | - Emmanuelle V LeBlanc
- Department of Biomedical and Molecular Sciences, Queen’s UniversityKingstonCanada
| | - Che C Colpitts
- Department of Biomedical and Molecular Sciences, Queen’s UniversityKingstonCanada
| | - Arinjay Banerjee
- Department of Veterinary Microbiology, University of SaskatchewanSaskatoonCanada
- Vaccine and Infectious Disease Organization, University of SaskatchewanSaskatoonCanada
- Department of Biology, University of WaterlooWaterlooCanada
- Department of Laboratory Medicine and Pathobiology, University of TorontoTorontoCanada
| | - Luca Pagani
- Estonian Biocentre, Institute of Genomics, University of TartuTartuEstonia
- Department of Biology, University of PadovaPadovaItaly
| | - Terence D Capellini
- Department of Human Evolutionary Biology, Harvard UniversityCambridgeUnited States
- Broad Institute of MIT and HarvardCambridgeUnited States
| |
Collapse
|
42
|
Zhang X, Kim B, Singh A, Sankararaman S, Durvasula A, Lohmueller KE. MaLAdapt Reveals Novel Targets of Adaptive Introgression From Neanderthals and Denisovans in Worldwide Human Populations. Mol Biol Evol 2023; 40:msad001. [PMID: 36617238 PMCID: PMC9887621 DOI: 10.1093/molbev/msad001] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/20/2022] [Revised: 12/25/2022] [Accepted: 12/28/2022] [Indexed: 01/09/2023] Open
Abstract
Adaptive introgression (AI) facilitates local adaptation in a wide range of species. Many state-of-the-art methods detect AI with ad-hoc approaches that identify summary statistic outliers or intersect scans for positive selection with scans for introgressed genomic regions. Although widely used, approaches intersecting outliers are vulnerable to a high false-negative rate as the power of different methods varies, especially for complex introgression events. Moreover, population genetic processes unrelated to AI, such as background selection or heterosis, may create similar genomic signals to AI, compromising the reliability of methods that rely on neutral null distributions. In recent years, machine learning (ML) methods have been increasingly applied to population genetic questions. Here, we present a ML-based method called MaLAdapt for identifying AI loci from genome-wide sequencing data. Using an Extra-Trees Classifier algorithm, our method combines information from a large number of biologically meaningful summary statistics to capture a powerful composite signature of AI across the genome. In contrast to existing methods, MaLAdapt is especially well-powered to detect AI with mild beneficial effects, including selection on standing archaic variation, and is robust to non-AI selective sweeps, heterosis from deleterious mutations, and demographic misspecification. Furthermore, MaLAdapt outperforms existing methods for detecting AI based on the analysis of simulated data and the validation of empirical signals through visual inspection of haplotype patterns. We apply MaLAdapt to the 1000 Genomes Project human genomic data and discover novel AI candidate regions in non-African populations, including genes that are enriched in functionally important biological pathways regulating metabolism and immune responses.
Collapse
Affiliation(s)
- Xinjun Zhang
- Department of Ecology and Evolutionary Biology, UCLA, Los Angeles, CA
| | - Bernard Kim
- Department of Biology, Stanford University, Palo Alto, CA
| | - Armaan Singh
- Department of Computer Science, UCLA, Los Angeles, CA
| | - Sriram Sankararaman
- Department of Computer Science, UCLA, Los Angeles, CA
- Department of Computational Medicine, UCLA, Los Angeles, CA
- Department of Human Genetics, David Geffen School of Medicine, UCLA, Los Angeles, CA
| | - Arun Durvasula
- Department of Genetics, Harvard Medical School, Boston, MA
| | - Kirk E Lohmueller
- Department of Ecology and Evolutionary Biology, UCLA, Los Angeles, CA
- Department of Human Genetics, David Geffen School of Medicine, UCLA, Los Angeles, CA
| |
Collapse
|
43
|
Kun E, Javan EM, Smith O, Gulamali F, de la Fuente J, Flynn BI, Vajrala K, Trutner Z, Jayakumar P, Tucker-Drob EM, Sohail M, Singh T, Narasimhan VM. The genetic architecture of the human skeletal form. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.01.03.521284. [PMID: 36712136 PMCID: PMC9881884 DOI: 10.1101/2023.01.03.521284] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
The human skeletal form underlies our ability to walk on two legs, but unlike standing height, the genetic basis of limb lengths and skeletal proportions is less well understood. Here we applied a deep learning model to 31,221 whole body dual-energy X-ray absorptiometry (DXA) images from the UK Biobank (UKB) to extract 23 different image-derived phenotypes (IDPs) that include all long bone lengths as well as hip and shoulder width, which we analyzed while controlling for height. All skeletal proportions are highly heritable (∼40-50%), and genome-wide association studies (GWAS) of these traits identified 179 independent loci, of which 102 loci were not associated with height. These loci are enriched in genes regulating skeletal development as well as associated with rare human skeletal diseases and abnormal mouse skeletal phenotypes. Genetic correlation and genomic structural equation modeling indicated that limb proportions exhibited strong genetic sharing but were genetically independent of width and torso proportions. Phenotypic and polygenic risk score analyses identified specific associations between osteoarthritis (OA) of the hip and knee, the leading causes of adult disability in the United States, and skeletal proportions of the corresponding regions. We also found genomic evidence of evolutionary change in arm-to-leg and hip-width proportions in humans consistent with striking anatomical changes in these skeletal proportions in the hominin fossil record. In contrast to cardiovascular, auto-immune, metabolic, and other categories of traits, loci associated with these skeletal proportions are significantly enriched in human accelerated regions (HARs), and regulatory elements of genes differentially expressed through development between humans and the great apes. Taken together, our work validates the use of deep learning models on DXA images to identify novel and specific genetic variants affecting the human skeletal form and ties a major evolutionary facet of human anatomical change to pathogenesis.
Collapse
Affiliation(s)
- Eucharist Kun
- Department of Integrative Biology, The University of Texas at Austin
| | - Emily M Javan
- Department of Integrative Biology, The University of Texas at Austin
| | - Olivia Smith
- Department of Integrative Biology, The University of Texas at Austin
| | | | | | - Brianna I Flynn
- Department of Integrative Biology, The University of Texas at Austin
| | - Kushal Vajrala
- Department of Integrative Biology, The University of Texas at Austin
| | - Zoe Trutner
- Department of Surgery and Perioperative Care, The University of Texas at Austin
| | - Prakash Jayakumar
- Department of Surgery and Perioperative Care, The University of Texas at Austin
| | | | - Mashaal Sohail
- Centro de Ciencias Genómicas (CCG), Universidad Nacional Autónoma de México (UNAM)
| | - Tarjinder Singh
- The Department of Psychiatry at Columbia University Irving Medical Center
- The New York Genome Center
- Mortimer B. Zuckerman Mind Brain Behavior Institute at Columbia University
| | - Vagheesh M Narasimhan
- Department of Integrative Biology, The University of Texas at Austin
- Department of Statistics and Data Science, The University of Texas at Austin
| |
Collapse
|
44
|
de March CA, Matsunami H, Abe M, Cobb M, Hoover KC. Genetic and functional odorant receptor variation in the Homo lineage. iScience 2022; 26:105908. [PMID: 36691623 PMCID: PMC9860384 DOI: 10.1016/j.isci.2022.105908] [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] [Received: 05/20/2022] [Revised: 10/07/2022] [Accepted: 12/26/2022] [Indexed: 12/30/2022] Open
Abstract
Humans, Neanderthals, and Denisovans independently adapted to a wide range of geographic environments and their associated food odors. Using ancient DNA sequences, we explored the in vitro function of thirty odorant receptor genes in the genus Homo. Our extinct relatives had highly conserved olfactory receptor sequence, but humans did not. Variations in odorant receptor protein sequence and structure may have produced variation in odor detection and perception. Variants led to minimal changes in specificity but had more influence on functional sensitivity. The few Neanderthal variants disturbed function, whereas Denisovan variants increased sensitivity to sweet and sulfur odors. Geographic adaptations may have produced greater functional variation in our lineage, increasing our olfactory repertoire and expanding our adaptive capacity. Our survey of olfactory genes and odorant receptors suggests that our genus has a shared repertoire with possible local ecological adaptations.
Collapse
Affiliation(s)
- Claire A. de March
- Institut de Chimie des Substances Naturelles, UPR2301 CNRS, Université Paris-Saclay, Gif-sur-Yvette 91190, France,Department of Molecular Genetics and Microbiology, Duke University Medical Center, Durham, NC 27710, USA,Department of Neurobiology, Duke Institute for Brain Sciences, Duke University, Durham, NC 27710, USA,Corresponding author
| | - Hiroaki Matsunami
- Department of Molecular Genetics and Microbiology, Duke University Medical Center, Durham, NC 27710, USA,Department of Neurobiology, Duke Institute for Brain Sciences, Duke University, Durham, NC 27710, USA
| | - Masashi Abe
- Department of Molecular Genetics and Microbiology, Duke University Medical Center, Durham, NC 27710, USA,Department of Biotechnology and Life Science, Tokyo University of Agriculture and Technology, Tokyo 184-8588, Japan
| | - Matthew Cobb
- Faculty of Life Sciences, The University of Manchester, Oxford Road, Manchester M13 9PL, UK
| | - Kara C. Hoover
- Department of Anthropology, University of Alaska Fairbanks, Fairbanks, AK 99775, USA,Corresponding author
| |
Collapse
|
45
|
Taufik L, Teixeira JC, Llamas B, Sudoyo H, Tobler R, Purnomo GA. Human Genetic Research in Wallacea and Sahul: Recent Findings and Future Prospects. Genes (Basel) 2022; 13:genes13122373. [PMID: 36553640 PMCID: PMC9778601 DOI: 10.3390/genes13122373] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/11/2022] [Revised: 12/08/2022] [Accepted: 12/13/2022] [Indexed: 12/23/2022] Open
Abstract
Genomic sequence data from worldwide human populations have provided a range of novel insights into our shared ancestry and the historical migrations that have shaped our global genetic diversity. However, a comprehensive understanding of these fundamental questions has been impeded by the lack of inclusion of many Indigenous populations in genomic surveys, including those from the Wallacean archipelago (which comprises islands of present-day Indonesia located east and west of Wallace's and Lydekker's Lines, respectively) and the former continent of Sahul (which once combined New Guinea and Australia during lower sea levels in the Pleistocene). Notably, these regions have been important areas of human evolution throughout the Late Pleistocene, as documented by diverse fossil and archaeological records which attest to the regional presence of multiple hominin species prior to the arrival of anatomically modern human (AMH) migrants. In this review, we collate and discuss key findings from the past decade of population genetic and phylogeographic literature focussed on the hominin history in Wallacea and Sahul. Specifically, we examine the evidence for the timing and direction of the ancient AMH migratory movements and subsequent hominin mixing events, emphasising several novel but consistent results that have important implications for addressing these questions. Finally, we suggest potentially lucrative directions for future genetic research in this key region of human evolution.
Collapse
Affiliation(s)
- Leonard Taufik
- Australian Centre for Ancient DNA, School of Biological Sciences, University of Adelaide, Adelaide, SA 5005, Australia
- Centre of Excellence for Australian Biodiversity and Heritage, University of Adelaide, Adelaide, SA 5005, Australia
- Mochtar Riady Institute for Nanotechnology, Tangerang 15810, Indonesia
- Correspondence: (L.T.); (G.A.P.)
| | - João C. Teixeira
- Australian Centre for Ancient DNA, School of Biological Sciences, University of Adelaide, Adelaide, SA 5005, Australia
- Centre of Excellence for Australian Biodiversity and Heritage, University of Adelaide, Adelaide, SA 5005, Australia
- Evolution of Cultural Diversity Initiative, Australian National University, Canberra, ACT 2601, Australia
- Centre for Interdisciplinary Studies, University of Coimbra, 3004-531 Coimbra, Portugal
| | - Bastien Llamas
- Australian Centre for Ancient DNA, School of Biological Sciences, University of Adelaide, Adelaide, SA 5005, Australia
- Centre of Excellence for Australian Biodiversity and Heritage, University of Adelaide, Adelaide, SA 5005, Australia
- Environment Institute, University of Adelaide, Adelaide, SA 5005, Australia
- National Centre for Indigenous Genomics, Australian National University, Canberra, ACT 2601, Australia
- Indigenous Genomics Research Group, Telethon Kids Institute, Adelaide, SA 5001, Australia
| | - Herawati Sudoyo
- Mochtar Riady Institute for Nanotechnology, Tangerang 15810, Indonesia
| | - Raymond Tobler
- Australian Centre for Ancient DNA, School of Biological Sciences, University of Adelaide, Adelaide, SA 5005, Australia
- Centre of Excellence for Australian Biodiversity and Heritage, University of Adelaide, Adelaide, SA 5005, Australia
- Evolution of Cultural Diversity Initiative, Australian National University, Canberra, ACT 2601, Australia
| | - Gludhug A. Purnomo
- Australian Centre for Ancient DNA, School of Biological Sciences, University of Adelaide, Adelaide, SA 5005, Australia
- Centre of Excellence for Australian Biodiversity and Heritage, University of Adelaide, Adelaide, SA 5005, Australia
- Correspondence: (L.T.); (G.A.P.)
| |
Collapse
|
46
|
Taravella Oill AM, Buetow KH, Wilson MA. The role of Neanderthal introgression in liver cancer. BMC Med Genomics 2022; 15:255. [PMID: 36503519 PMCID: PMC9743633 DOI: 10.1186/s12920-022-01405-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2022] [Accepted: 11/25/2022] [Indexed: 12/14/2022] Open
Abstract
BACKGROUND Neanderthal introgressed DNA has been linked to different normal and disease traits including immunity and metabolism-two important functions that are altered in liver cancer. However, there is limited understanding of the relationship between Neanderthal introgression and liver cancer risk. The aim of this study was to investigate the relationship between Neanderthal introgression and liver cancer risk. METHODS Using germline and somatic DNA and tumor RNA from liver cancer patients from The Cancer Genome Atlas, along with ancestry-match germline DNA from unaffected individuals from the 1000 Genomes Resource, and allele specific expression data from normal liver tissue from The Genotype-Tissue Expression project we investigated whether Neanderthal introgression impacts cancer etiology. Using a previously generated set of Neanderthal alleles, we identified Neanderthal introgressed haplotypes. We then tested whether somatic mutations are enriched or depleted on Neanderthal introgressed haplotypes compared to modern haplotypes. We also computationally assessed whether somatic mutations have a functional effect or show evidence of regulating expression of Neanderthal haplotypes. Finally, we compared patterns of Neanderthal introgression in liver cancer patients and the general population. RESULTS We find Neanderthal introgressed haplotypes exhibit an excess of somatic mutations compared to modern haplotypes. Variant Effect Predictor analysis revealed that most of the somatic mutations on these Neanderthal introgressed haplotypes are not functional. We did observe expression differences of Neanderthal alleles between tumor and normal for four genes that also showed a pattern of enrichment of somatic mutations on Neanderthal haplotypes. However, gene expression was similar between liver cancer patients with modern ancestry and liver cancer patients with Neanderthal ancestry at these genes. Provocatively, when analyzing all genes, we find evidence of Neanderthal introgression regulating expression in tumor from liver cancer patients in two genes, ARK1C4 and OAS1. Finally, we find that most genes do not show a difference in the proportion of Neanderthal introgression between liver cancer patients and the general population. CONCLUSION Our results suggest that Neanderthal introgression provides opportunity for somatic mutations to accumulate, and that some Neanderthal introgression may impact liver cancer risk.
Collapse
Affiliation(s)
- Angela M Taravella Oill
- Center for Evolution and Medicine, Arizona State University, Tempe, AZ, USA.
- School of Life Sciences, Arizona State University, Tempe, AZ, USA.
| | - Kenneth H Buetow
- Center for Evolution and Medicine, Arizona State University, Tempe, AZ, USA
- School of Life Sciences, Arizona State University, Tempe, AZ, USA
| | - Melissa A Wilson
- Center for Evolution and Medicine, Arizona State University, Tempe, AZ, USA
- School of Life Sciences, Arizona State University, Tempe, AZ, USA
| |
Collapse
|
47
|
Vespasiani DM, Jacobs GS, Cook LE, Brucato N, Leavesley M, Kinipi C, Ricaut FX, Cox MP, Gallego Romero I. Denisovan introgression has shaped the immune system of present-day Papuans. PLoS Genet 2022; 18:e1010470. [PMID: 36480515 PMCID: PMC9731433 DOI: 10.1371/journal.pgen.1010470] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2022] [Accepted: 10/10/2022] [Indexed: 12/13/2022] Open
Abstract
Modern humans have admixed with multiple archaic hominins. Papuans, in particular, owe up to 5% of their genome to Denisovans, a sister group to Neanderthals whose remains have only been identified in Siberia and Tibet. Unfortunately, the biological and evolutionary significance of these introgression events remain poorly understood. Here we investigate the function of both Denisovan and Neanderthal alleles characterised within a set of 56 genomes from Papuan individuals. By comparing the distribution of archaic and non-archaic variants we assess the consequences of archaic admixture across a multitude of different cell types and functional elements. We observe an enrichment of archaic alleles within cis-regulatory elements and transcribed regions of the genome, with Denisovan variants strongly affecting elements active within immune-related cells. We identify 16,048 and 10,032 high-confidence Denisovan and Neanderthal variants that fall within annotated cis-regulatory elements and with the potential to alter the affinity of multiple transcription factors to their cognate DNA motifs, highlighting a likely mechanism by which introgressed DNA can impact phenotypes. Lastly, we experimentally validate these predictions by testing the regulatory potential of five Denisovan variants segregating within Papuan individuals, and find that two are associated with a significant reduction of transcriptional activity in plasmid reporter assays. Together, these data provide support for a widespread contribution of archaic DNA in shaping the present levels of modern human genetic diversity, with different archaic ancestries potentially affecting multiple phenotypic traits within non-Africans.
Collapse
Affiliation(s)
- Davide M. Vespasiani
- Melbourne Integrative Genomics, University of Melbourne, Parkville, Australia
- School of Biosciences, University of Melbourne, Parkville, Australia
| | - Guy S. Jacobs
- Department of Archaeology, University of Cambridge, Cambridge, Uniteed Kingdom
| | - Laura E. Cook
- Melbourne Integrative Genomics, University of Melbourne, Parkville, Australia
- School of Biosciences, University of Melbourne, Parkville, Australia
| | - Nicolas Brucato
- Laboratoire de Evolution et Diversite Biologique, Université de Toulouse Midi-Pyrénées, Toulouse, France
| | - Matthew Leavesley
- School of Humanities and Social Sciences, University of Papua New Guinea, Port Moresby, Papua New Guinea
- College of Arts, Society and Education, James Cook University, Cairns, Australia
- ARC Centre of Excellence for Australian Biodiversity and Heritage, University of Wollongong, Wollongong, Australia
| | - Christopher Kinipi
- School of Humanities and Social Sciences, University of Papua New Guinea, Port Moresby, Papua New Guinea
| | - François-Xavier Ricaut
- Laboratoire de Evolution et Diversite Biologique, Université de Toulouse Midi-Pyrénées, Toulouse, France
| | - Murray P. Cox
- School of Natural Sciences, Massey University, Palmerston North, New Zealand
| | - Irene Gallego Romero
- Melbourne Integrative Genomics, University of Melbourne, Parkville, Australia
- School of Biosciences, University of Melbourne, Parkville, Australia
- Center for Stem Cell Systems, University of Melbourne, Parkville, Australia
- Center for Genomics, Evolution and Medicine, University of Tartu, Tartu, Estonia
- * E-mail:
| |
Collapse
|
48
|
Askapuli A, Vilar M, Garcia-Ortiz H, Zhabagin M, Sabitov Z, Akilzhanova A, Ramanculov E, Schamiloglu U, Martinez-Hernandez A, Contreras-Cubas C, Barajas-Olmos F, Schurr TG, Zhumadilov Z, Flores-Huacuja M, Orozco L, Hawks J, Saitou N. Kazak mitochondrial genomes provide insights into the human population history of Central Eurasia. PLoS One 2022; 17:e0277771. [PMID: 36445929 PMCID: PMC9707748 DOI: 10.1371/journal.pone.0277771] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/27/2022] [Accepted: 11/03/2022] [Indexed: 11/30/2022] Open
Abstract
As a historical nomadic group in Central Asia, Kazaks have mainly inhabited the steppe zone from the Altay Mountains in the East to the Caspian Sea in the West. Fine scale characterization of the genetic profile and population structure of Kazaks would be invaluable for understanding their population history and modeling prehistoric human expansions across the Eurasian steppes. With this mind, we characterized the maternal lineages of 200 Kazaks from Jetisuu at mitochondrial genome level. Our results reveal that Jetisuu Kazaks have unique mtDNA haplotypes including those belonging to the basal branches of both West Eurasian (R0, H, HV) and East Eurasian (A, B, C, D) lineages. The great diversity observed in their maternal lineages may reflect pivotal geographic location of Kazaks in Eurasia and implies a complex history for this population. Comparative analyses of mitochondrial genomes of human populations in Central Eurasia reveal a common maternal genetic ancestry for Turko-Mongolian speakers and their expansion being responsible for the presence of East Eurasian maternal lineages in Central Eurasia. Our analyses further indicate maternal genetic affinity between the Sherpas from the Tibetan Plateau with the Turko-Mongolian speakers.
Collapse
Affiliation(s)
- Ayken Askapuli
- School of Sciences and Humanities, Nazarbayev University, Astana, Kazakhstan
- National Center for Biotechnology, Astana, Kazakhstan
- National Laboratory Astana, Nazarbayev University, Astana, Kazakhstan
- Department of Integrative Biology, University of Wisconsin-Madison, Madison, Wisconsin, United States of America
| | - Miguel Vilar
- The Genographic Project, National Geographic Society, Washington, DC, United States of America
- Department of Anthropology, University of Maryland, College Park, Maryland, United States of America
| | - Humberto Garcia-Ortiz
- Immunogenomics and Metabolic Diseases Laboratory, National Institute of Genomic Medicine, Mexico City, Mexico
| | - Maxat Zhabagin
- School of Sciences and Humanities, Nazarbayev University, Astana, Kazakhstan
- National Center for Biotechnology, Astana, Kazakhstan
- National Laboratory Astana, Nazarbayev University, Astana, Kazakhstan
| | | | - Ainur Akilzhanova
- National Laboratory Astana, Nazarbayev University, Astana, Kazakhstan
| | - Erlan Ramanculov
- School of Sciences and Humanities, Nazarbayev University, Astana, Kazakhstan
- National Center for Biotechnology, Astana, Kazakhstan
| | - Uli Schamiloglu
- School of Sciences and Humanities, Nazarbayev University, Astana, Kazakhstan
| | - Angelica Martinez-Hernandez
- Immunogenomics and Metabolic Diseases Laboratory, National Institute of Genomic Medicine, Mexico City, Mexico
| | - Cecilia Contreras-Cubas
- Immunogenomics and Metabolic Diseases Laboratory, National Institute of Genomic Medicine, Mexico City, Mexico
| | - Francisco Barajas-Olmos
- Immunogenomics and Metabolic Diseases Laboratory, National Institute of Genomic Medicine, Mexico City, Mexico
| | - Theodore G. Schurr
- Department of Anthropology, University of Pennsylvania, Philadelphia, Pennsylvania, United States of America
| | - Zhaxybay Zhumadilov
- National Laboratory Astana, Nazarbayev University, Astana, Kazakhstan
- School of Medicine, Nazarbayev University, Astana, Kazakhstan
| | - Marlen Flores-Huacuja
- Immunogenomics and Metabolic Diseases Laboratory, National Institute of Genomic Medicine, Mexico City, Mexico
| | - Lorena Orozco
- Immunogenomics and Metabolic Diseases Laboratory, National Institute of Genomic Medicine, Mexico City, Mexico
| | - John Hawks
- Department of Integrative Biology, University of Wisconsin-Madison, Madison, Wisconsin, United States of America
- Department of Anthropology, University of Wisconsin-Madison, Madison, Wisconsin, United States of America
| | - Naruya Saitou
- Population Genetics Laboratory, National Institute of Genetics, Mishima, Shizuoka, Japan
- Department of Biological Sciences, Graduate School of Science, University of Tokyo, Tokyo, Japan
- Advanced Medical Research Center, Faculty of Medicine, University of the Ryukyus, Okinawa Ken, Japan
| |
Collapse
|
49
|
Zhu K, Du P, Li J, Zhang J, Hu X, Meng H, Chen L, Zhou B, Yang X, Xiong J, Allen E, Ren X, Ding Y, Xu Y, Chang X, Yu Y, Han S, Dong G, Wang CC, Wen S. Cultural and demic co-diffusion of Tubo Empire on Tibetan Plateau. iScience 2022; 25:105636. [PMID: 36582485 PMCID: PMC9792914 DOI: 10.1016/j.isci.2022.105636] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/14/2022] [Revised: 10/15/2022] [Accepted: 11/17/2022] [Indexed: 11/23/2022] Open
Abstract
A high point of Tibetan Plateau (TP) civilization, the expansive Tubo Empire (618-842 AD) wielded great influence across ancient western China. However, whether the Tubo expansion was cultural or demic remains unclear due to sparse ancient DNA sampling. Here, we reported ten ancient genomes at 0.017- to 0.867-fold coverages from the Dulan site with typical Tubo archaeological culture dating to 1308-1130 BP. Nine individuals from three different grave types have close relationship with previously reported ancient highlanders from the southwestern Himalayas and modern core-Tibetan populations. A Dulan-related Tubo ancestry contributed overwhelmingly (95%-100%) to the formation of modern Tibetans. A genetic outlier with dominant Eurasian steppe-related ancestry suggesting a potential population movement into the Tubo-controlled regions from Central Asia. Together with archeological evidence from burial styles and customs, our study suggested the impact of the Tubo empire on the northeast edge of the TP involved both cultural and demic diffusion.
Collapse
Affiliation(s)
- Kongyang Zhu
- Institute of Archaeological Science, Fudan University, Shanghai 200433, China,State Key Laboratory of Cellular Stress Biology, School of Life Sciences, Xiamen University, Xiamen 361102, China
| | - Panxin Du
- Ministry of Education Key Laboratory of Contemporary Anthropology, Department of Anthropology and Human Genetics, School of Life Sciences, Fudan University, Shanghai 200433, China,State Key Laboratory of Genetic Engineering, Collaborative Innovation Center for Genetics and Development, School of Life Sciences, and Human Phenome Institute, Fudan University, Shanghai 200433, China
| | - Jiyuan Li
- Qinghai Provincial Cultural Relics and Archaeology Institute, Xining 810007, China
| | | | - Xiaojun Hu
- Qinghai Provincial Cultural Relics and Archaeology Institute, Xining 810007, China
| | - Hailiang Meng
- Ministry of Education Key Laboratory of Contemporary Anthropology, Department of Anthropology and Human Genetics, School of Life Sciences, Fudan University, Shanghai 200433, China
| | - Liang Chen
- China-Central Asia Human Environmental “the Belt and Road” Joint Laboratory, Northwest University, Xi’an 710127, China
| | - Boyan Zhou
- Division of Biostatistics, Department of Population Health, School of Medicine, New York University, New York, NY 10016, USA
| | - Xiaomin Yang
- Department of Anthropology and Ethnology, Institute of Anthropology, School of Sociology and Anthropology, Xiamen University, Xiamen 361005, China
| | - Jianxue Xiong
- Institute of Archaeological Science, Fudan University, Shanghai 200433, China
| | - Edward Allen
- Institute of Archaeological Science, Fudan University, Shanghai 200433, China
| | - Xiaoying Ren
- Institute of Archaeological Science, Fudan University, Shanghai 200433, China
| | - Yi Ding
- Institute of Archaeological Science, Fudan University, Shanghai 200433, China
| | - Yiran Xu
- Institute of Archaeological Science, Fudan University, Shanghai 200433, China
| | - Xin Chang
- Institute of Archaeological Science, Fudan University, Shanghai 200433, China
| | - Yao Yu
- Institute of Archaeological Science, Fudan University, Shanghai 200433, China
| | - Sheng Han
- Ministry of Education Key Laboratory of Contemporary Anthropology, Department of Anthropology and Human Genetics, School of Life Sciences, Fudan University, Shanghai 200433, China
| | - Guanghui Dong
- Key Laboratory of Western China’s Environmental Systems (Ministry of Education), College of Earth and Environmental Sciences, Lanzhou University, Lanzhou 730000, China
| | - Chuan-Chao Wang
- State Key Laboratory of Cellular Stress Biology, School of Life Sciences, Xiamen University, Xiamen 361102, China,Ministry of Education Key Laboratory of Contemporary Anthropology, Department of Anthropology and Human Genetics, School of Life Sciences, Fudan University, Shanghai 200433, China,Department of Anthropology and Ethnology, Institute of Anthropology, School of Sociology and Anthropology, Xiamen University, Xiamen 361005, China,State Key Laboratory of Marine Environmental Science, Xiamen University, Xiamen 361102, China,Institute of Artificial Intelligence, Xiamen University, Xiamen 361005, China,Corresponding author
| | - Shaoqing Wen
- Institute of Archaeological Science, Fudan University, Shanghai 200433, China,Ministry of Education Key Laboratory of Contemporary Anthropology, Department of Anthropology and Human Genetics, School of Life Sciences, Fudan University, Shanghai 200433, China,MOE Laboratory for National Development and Intelligent Governance, Fudan University, Shanghai 200433, China,Center for the Belt and Road Archaeology and Ancient Civilizations, Shanghai 200433, China,Corresponding author
| |
Collapse
|
50
|
Koller D, Wendt FR, Pathak GA, De Lillo A, De Angelis F, Cabrera-Mendoza B, Tucci S, Polimanti R. Denisovan and Neanderthal archaic introgression differentially impacted the genetics of complex traits in modern populations. BMC Biol 2022; 20:249. [PMID: 36344982 PMCID: PMC9641937 DOI: 10.1186/s12915-022-01449-2] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2022] [Accepted: 10/24/2022] [Indexed: 11/09/2022] Open
Abstract
BACKGROUND Introgression from extinct Neanderthal and Denisovan human species has been shown to contribute to the genetic pool of modern human populations and their phenotypic spectrum. Evidence of how Neanderthal introgression shaped the genetics of human traits and diseases has been extensively studied in populations of European descent, with signatures of admixture reported for instance in genes associated with pigmentation, immunity, and metabolic traits. However, limited information is currently available about the impact of archaic introgression on other ancestry groups. Additionally, to date, no study has been conducted with respect to the impact of Denisovan introgression on the health and disease of modern populations. Here, we compare the way evolutionary pressures shaped the genetics of complex traits in East Asian and European populations, and provide evidence of the impact of Denisovan introgression on the health of East Asian and Central/South Asian populations. RESULTS Leveraging genome-wide association statistics from the Biobank Japan and UK Biobank, we assessed whether Denisovan and Neanderthal introgression together with other evolutionary genomic signatures were enriched for the heritability of physiological and pathological conditions in populations of East Asian and European descent. In EAS, Denisovan-introgressed loci were enriched for coronary artery disease heritability (1.69-fold enrichment, p=0.003). No enrichment for archaic introgression was observed in EUR. We also performed a phenome-wide association study of Denisovan and Neanderthal alleles in six ancestry groups available in the UK Biobank. In EAS, the Denisovan-introgressed SNP rs62391664 in the major histocompatibility complex region was associated with albumin/globulin ratio (beta=-0.17, p=3.57×10-7). Neanderthal-introgressed alleles were associated with psychiatric and cognitive traits in EAS (e.g., "No Bipolar or Depression"-rs79043717 beta=-1.5, p=1.1×10-7), and with blood biomarkers (e.g., alkaline phosphatase-rs11244089 beta=0.1, p=3.69×10-116) and red hair color (rs60733936 beta=-0.86, p=4.49×10-165) in EUR. In the other ancestry groups, Neanderthal alleles were associated with several traits, also including the use of certain medications (e.g., Central/South East Asia: indapamide - rs732632 beta=-2.38, p=5.22×10-7). CONCLUSIONS Our study provides novel evidence regarding the impact of archaic introgression on the genetics of complex traits in worldwide populations, highlighting the specific contribution of Denisovan introgression in EAS populations.
Collapse
Affiliation(s)
- Dora Koller
- Department of Psychiatry, Yale University School of Medicine, West Haven, CT, 06516, USA
- VA CT Healthcare Center, West Haven, CT, 06516, USA
- Department of Genetics, Microbiology and Statistics, Faculty of Biology, University of Barcelona, Barcelona, Catalonia, 08028, Spain
| | - Frank R Wendt
- Department of Psychiatry, Yale University School of Medicine, West Haven, CT, 06516, USA
- VA CT Healthcare Center, West Haven, CT, 06516, USA
| | - Gita A Pathak
- Department of Psychiatry, Yale University School of Medicine, West Haven, CT, 06516, USA
- VA CT Healthcare Center, West Haven, CT, 06516, USA
| | - Antonella De Lillo
- Department of Psychiatry, Yale University School of Medicine, West Haven, CT, 06516, USA
- Department of Biology, University of Rome "Tor Vergata", Rome, 00133, Italy
| | - Flavio De Angelis
- Department of Psychiatry, Yale University School of Medicine, West Haven, CT, 06516, USA
- VA CT Healthcare Center, West Haven, CT, 06516, USA
- Department of Biology, University of Rome "Tor Vergata", Rome, 00133, Italy
| | - Brenda Cabrera-Mendoza
- Department of Psychiatry, Yale University School of Medicine, West Haven, CT, 06516, USA
- VA CT Healthcare Center, West Haven, CT, 06516, USA
| | - Serena Tucci
- Department of Anthropology, Yale University, New Haven, CT, 06511, USA
| | - Renato Polimanti
- Department of Psychiatry, Yale University School of Medicine, West Haven, CT, 06516, USA.
- VA CT Healthcare Center, West Haven, CT, 06516, USA.
| |
Collapse
|