1
|
Kamitaki N, Hujoel MLA, Mukamel RE, Gebara E, McCarroll SA, Loh PR. A sequence of SVA retrotransposon insertions in ASIP shaped human pigmentation. Nat Genet 2024; 56:1583-1591. [PMID: 39048794 PMCID: PMC11319198 DOI: 10.1038/s41588-024-01841-4] [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: 08/08/2023] [Accepted: 06/21/2024] [Indexed: 07/27/2024]
Abstract
Retrotransposons comprise about 45% of the human genome1, but their contributions to human trait variation and evolution are only beginning to be explored2,3. Here, we find that a sequence of SVA retrotransposon insertions in an early intron of the ASIP (agouti signaling protein) gene has probably shaped human pigmentation several times. In the UK Biobank (n = 169,641), a recent 3.3-kb SVA insertion polymorphism associated strongly with lighter skin pigmentation (0.22 [0.21-0.23] s.d.; P = 2.8 × 10-351) and increased skin cancer risk (odds ratio = 1.23 [1.18-1.27]; P = 1.3 × 10-28), appearing to underlie one of the strongest common genetic influences on these phenotypes within European populations4-6. ASIP expression in skin displayed the same association pattern, with the SVA insertion allele exhibiting 2.2-fold (1.9-2.6) increased expression. This effect had an unusual apparent mechanism: an earlier, nonpolymorphic, human-specific SVA retrotransposon 3.9 kb upstream appeared to have caused ASIP hypofunction by nonproductive splicing, which the new (polymorphic) SVA insertion largely eliminated. Extended haplotype homozygosity indicated that the insertion allele has risen to allele frequencies up to 11% in European populations over the past several thousand years. These results indicate that a sequence of retrotransposon insertions contributed to a species-wide increase, then a local decrease, of human pigmentation.
Collapse
Affiliation(s)
- Nolan Kamitaki
- Division of Genetics, Department of Medicine, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, USA.
- Center for Data Sciences, Brigham and Women's Hospital, Boston, MA, USA.
- Program in Medical and Population Genetics, Broad Institute of MIT and Harvard, Cambridge, MA, USA.
- Stanley Center for Psychiatric Research, Broad Institute of MIT and Harvard, Cambridge, MA, USA.
- Department of Genetics, Harvard Medical School, Boston, MA, USA.
- Department of Biomedical Informatics, Harvard Medical School, Boston, MA, USA.
| | - Margaux L A Hujoel
- Division of Genetics, Department of Medicine, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, USA
- Center for Data Sciences, Brigham and Women's Hospital, Boston, MA, USA
- Program in Medical and Population Genetics, Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - Ronen E Mukamel
- Division of Genetics, Department of Medicine, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, USA
- Center for Data Sciences, Brigham and Women's Hospital, Boston, MA, USA
- Program in Medical and Population Genetics, Broad Institute of MIT and Harvard, Cambridge, MA, USA
| | - Edward Gebara
- Program in Medical and Population Genetics, Broad Institute of MIT and Harvard, Cambridge, MA, USA
- Stanley Center for Psychiatric Research, Broad Institute of MIT and Harvard, Cambridge, MA, USA
- Department of Genetics, Harvard Medical School, Boston, MA, USA
| | - Steven A McCarroll
- Program in Medical and Population Genetics, Broad Institute of MIT and Harvard, Cambridge, MA, USA.
- Stanley Center for Psychiatric Research, Broad Institute of MIT and Harvard, Cambridge, MA, USA.
- Department of Genetics, Harvard Medical School, Boston, MA, USA.
| | - Po-Ru Loh
- Division of Genetics, Department of Medicine, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, USA.
- Center for Data Sciences, Brigham and Women's Hospital, Boston, MA, USA.
- Program in Medical and Population Genetics, Broad Institute of MIT and Harvard, Cambridge, MA, USA.
| |
Collapse
|
2
|
Pfennig A, Lachance J. The evolutionary fate of Neanderthal DNA in 30,780 admixed genomes with recent African-like ancestry. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.07.25.605203. [PMID: 39091830 PMCID: PMC11291122 DOI: 10.1101/2024.07.25.605203] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Subscribe] [Scholar Register] [Indexed: 08/04/2024]
Abstract
Following introgression, Neanderthal DNA was initially purged from non-African genomes, but the evolutionary fate of remaining introgressed DNA has not been explored yet. To fill this gap, we analyzed 30,780 admixed genomes with African-like ancestry from the All of Us research program, in which Neanderthal alleles encountered novel genetic backgrounds during the last 15 generations. Observed amounts of Neanderthal DNA approximately match expectations based on ancestry proportions, suggesting neutral evolution. Nevertheless, we identified genomic regions that have significantly less or more Neanderthal ancestry than expected and are associated with spermatogenesis, innate immunity, and other biological processes. We also identified three novel introgression desert-like regions in recently admixed genomes, whose genetic features are compatible with hybrid incompatibilities and intrinsic negative selection. Overall, we find that much of the remaining Neanderthal DNA in human genomes is not under strong selection, and complex evolutionary dynamics have shaped introgression landscapes in our species.
Collapse
Affiliation(s)
- Aaron Pfennig
- School of Biological Sciences, Georgia Institute of Technology, 950 Atlantic Dr, Atlanta, 30332, GA, USA
| | - Joseph Lachance
- School of Biological Sciences, Georgia Institute of Technology, 950 Atlantic Dr, Atlanta, 30332, GA, USA
| |
Collapse
|
3
|
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
|
4
|
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
|
5
|
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
|
6
|
Aneli S, Ceccatelli Berti C, Gilea AI, Birolo G, Mutti G, Pavesi A, Baruffini E, Goffrini P, Capelli C. Functional characterization of archaic-specific variants in mitonuclear genes: insights from comparative analysis in S. cerevisiae. Hum Mol Genet 2024; 33:1152-1163. [PMID: 38558123 DOI: 10.1093/hmg/ddae057] [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: 01/17/2024] [Revised: 02/29/2024] [Accepted: 03/14/2024] [Indexed: 04/04/2024] Open
Abstract
Neanderthal and Denisovan hybridisation with modern humans has generated a non-random genomic distribution of introgressed regions, the result of drift and selection dynamics. Cross-species genomic incompatibility and more efficient removal of slightly deleterious archaic variants have been proposed as selection-based processes involved in the post-hybridisation purge of archaic introgressed regions. Both scenarios require the presence of functionally different alleles across Homo species onto which selection operated differently according to which populations hosted them, but only a few of these variants have been pinpointed so far. In order to identify functionally divergent archaic variants removed in humans, we focused on mitonuclear genes, which are underrepresented in the genomic landscape of archaic humans. We searched for non-synonymous, fixed, archaic-derived variants present in mitonuclear genes, rare or absent in human populations. We then compared the functional impact of archaic and human variants in the model organism Saccharomyces cerevisiae. Notably, a variant within the mitochondrial tyrosyl-tRNA synthetase 2 (YARS2) gene exhibited a significant decrease in respiratory activity and a substantial reduction of Cox2 levels, a proxy for mitochondrial protein biosynthesis, coupled with the accumulation of the YARS2 protein precursor and a lower amount of mature enzyme. Our work suggests that this variant is associated with mitochondrial functionality impairment, thus contributing to the purging of archaic introgression in YARS2. While different molecular mechanisms may have impacted other mitonuclear genes, our approach can be extended to the functional screening of mitonuclear genetic variants present across species and populations.
Collapse
Affiliation(s)
- Serena Aneli
- Department of Public Health Sciences and Pediatrics, University of Turin, C.so Galileo Galilei 22, Turin 10126, Italy
| | - Camilla Ceccatelli Berti
- Department of Chemistry, Life Sciences and Environmental Sustainability, University of Parma, Parco Area delle Scienze 11/a, Parma 43124, Italy
| | - Alexandru Ionut Gilea
- Department of Chemistry, Life Sciences and Environmental Sustainability, University of Parma, Parco Area delle Scienze 11/a, Parma 43124, Italy
| | - Giovanni Birolo
- Department of Medical Sciences, University of Turin, Via Santena 5, Turin 10126, Italy
| | - Giacomo Mutti
- Barcelona Supercomputing Centre (BSC-CNS), Department of Life Sciences, Plaça Eusebi Güell, 1-3, Barcelona 08034, Spain
- Institute for Research in Biomedicine (IRB Barcelona), Department of Mechanisms of Disease, The Barcelona Institute of Science and Technology, Baldiri Reixac, 10, Barcelona 08028, Spain
| | - Angelo Pavesi
- Department of Chemistry, Life Sciences and Environmental Sustainability, University of Parma, Parco Area delle Scienze 11/a, Parma 43124, Italy
| | - Enrico Baruffini
- Department of Chemistry, Life Sciences and Environmental Sustainability, University of Parma, Parco Area delle Scienze 11/a, Parma 43124, Italy
| | - Paola Goffrini
- Department of Chemistry, Life Sciences and Environmental Sustainability, University of Parma, Parco Area delle Scienze 11/a, Parma 43124, Italy
| | - Cristian Capelli
- Department of Chemistry, Life Sciences and Environmental Sustainability, University of Parma, Parco Area delle Scienze 11/a, Parma 43124, Italy
- Department of Biology, University of Oxford, 11a Mansfield Rd, Oxford OX1 3SZ, United Kingdom
| |
Collapse
|
7
|
Bolognini D, Halgren A, Lou RN, Raveane A, Rocha JL, Guarracino A, Soranzo N, Chin J, Garrison E, Sudmant PH. Global diversity, recurrent evolution, and recent selection on amylase structural haplotypes in humans. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.02.07.579378. [PMID: 38370750 PMCID: PMC10871346 DOI: 10.1101/2024.02.07.579378] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/20/2024]
Abstract
The adoption of agriculture, first documented ~12,000 years ago in the Fertile Crescent, triggered a rapid shift toward starch-rich diets in human populations. Amylase genes facilitate starch digestion and increased salivary amylase copy number has been observed in some modern human populations with high starch intake, though evidence of recent selection is lacking. Here, using 52 long-read diploid assemblies and short read data from ~5,600 contemporary and ancient humans, we resolve the diversity, evolutionary history, and selective impact of structural variation at the amylase locus. We find that amylase genes have higher copy numbers in populations with agricultural subsistence compared to fishing, hunting, and pastoral groups. We identify 28 distinct amylase structural architectures and demonstrate that nearly identical structures have arisen recurrently on different haplotype backgrounds throughout recent human history. AMY1 and AMY2A genes each exhibit multiple duplications/deletions with mutation rates >10,000-fold the SNP mutation rate, whereas AMY2B gene duplications share a single origin. Using a pangenome graph-based approach to infer structural haplotypes across thousands of humans, we identify extensively duplicated haplotypes present at higher frequencies in modern day populations with traditionally agricultural diets. Leveraging 533 ancient human genomes we find that duplication-containing haplotypes (i.e. haplotypes with more amylase gene copies than the ancestral haplotype) have increased in frequency more than seven-fold over the last 12,000 years providing evidence for recent selection in West Eurasians. Together, our study highlights the potential impacts of the agricultural revolution on human genomes and the importance of long-read sequencing in identifying signatures of selection at structurally complex loci.
Collapse
Affiliation(s)
| | - Alma Halgren
- Department of Integrative Biology, University of California Berkeley, Berkeley, USA
| | - Runyang Nicolas Lou
- Department of Integrative Biology, University of California Berkeley, Berkeley, USA
| | | | - Joana L Rocha
- Department of Integrative Biology, University of California Berkeley, Berkeley, USA
| | - Andrea Guarracino
- Department of Genetics, Genomics, and Informatics, University of Tennessee Health Science Center, Memphis, USA
| | | | - Jason Chin
- Foundation for Biological Data Science, Belmont, USA
| | - Erik Garrison
- Department of Genetics, Genomics, and Informatics, University of Tennessee Health Science Center, Memphis, USA
| | - Peter H Sudmant
- Department of Integrative Biology, University of California Berkeley, Berkeley, USA
- Center for Computational Biology, University of California Berkeley, Berkeley, USA
| |
Collapse
|
8
|
Brázda V, Šislerová L, Cucchiarini A, Mergny JL. G-quadruplex propensity in H. neanderthalensis, H. sapiens and Denisovans mitochondrial genomes. NAR Genom Bioinform 2024; 6:lqae060. [PMID: 38817800 PMCID: PMC11137754 DOI: 10.1093/nargab/lqae060] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2023] [Revised: 04/18/2024] [Accepted: 05/17/2024] [Indexed: 06/01/2024] Open
Abstract
Current methods of processing archaeological samples combined with advances in sequencing methods lead to disclosure of a large part of H. neanderthalensis and Denisovans genetic information. It is hardly surprising that the genome variability between modern humans, Denisovans and H. neanderthalensis is relatively limited. Genomic studies may provide insight on the metabolism of extinct human species or lineages. Detailed analysis of G-quadruplex sequences in H. neanderthalensis and Denisovans mitochondrial DNA showed us interesting features. Relatively similar patterns in mitochondrial DNA are found compared to modern humans, with one notable exception for H. neanderthalensis. An interesting difference between H. neanderthalensis and H. sapiens corresponds to a motif found in the D-loop region of mtDNA, which is responsible for mitochondrial DNA replication. This area is directly responsible for the number of mitochondria and consequently for the efficient energy metabolism of cell. H. neanderthalensis harbor a long uninterrupted run of guanines in this region, which may cause problems for replication, in contrast with H. sapiens, for which this run is generally shorter and interrupted. One may propose that the predominant H. sapiens motif provided a selective advantage for modern humans regarding mtDNA replication and function.
Collapse
Affiliation(s)
- Václav Brázda
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 00 Brno, Czech Republic
- Brno University of Technology, Faculty of Chemistry, Purkyňova 118, 612 00 Brno, Czech Republic
| | - Lucie Šislerová
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 00 Brno, Czech Republic
- Brno University of Technology, Faculty of Chemistry, Purkyňova 118, 612 00 Brno, Czech Republic
| | - Anne Cucchiarini
- Laboratoire d’Optique et Biosciences (LOB), Ecole Polytechnique, CNRS, INSERM, Institut Polytechnique de Paris, 91120 Palaiseau, France
| | - Jean-Louis Mergny
- Institute of Biophysics of the Czech Academy of Sciences, Královopolská 135, 612 00 Brno, Czech Republic
- Laboratoire d’Optique et Biosciences (LOB), Ecole Polytechnique, CNRS, INSERM, Institut Polytechnique de Paris, 91120 Palaiseau, France
| |
Collapse
|
9
|
Stankey CT, Bourges C, Haag LM, Turner-Stokes T, Piedade AP, Palmer-Jones C, Papa I, Silva Dos Santos M, Zhang Q, Cameron AJ, Legrini A, Zhang T, Wood CS, New FN, Randzavola LO, Speidel L, Brown AC, Hall A, Saffioti F, Parkes EC, Edwards W, Direskeneli H, Grayson PC, Jiang L, Merkel PA, Saruhan-Direskeneli G, Sawalha AH, Tombetti E, Quaglia A, Thorburn D, Knight JC, Rochford AP, Murray CD, Divakar P, Green M, Nye E, MacRae JI, Jamieson NB, Skoglund P, Cader MZ, Wallace C, Thomas DC, Lee JC. A disease-associated gene desert directs macrophage inflammation through ETS2. Nature 2024; 630:447-456. [PMID: 38839969 PMCID: PMC11168933 DOI: 10.1038/s41586-024-07501-1] [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/17/2023] [Accepted: 05/01/2024] [Indexed: 06/07/2024]
Abstract
Increasing rates of autoimmune and inflammatory disease present a burgeoning threat to human health1. This is compounded by the limited efficacy of available treatments1 and high failure rates during drug development2, highlighting an urgent need to better understand disease mechanisms. Here we show how functional genomics could address this challenge. By investigating an intergenic haplotype on chr21q22-which has been independently linked to inflammatory bowel disease, ankylosing spondylitis, primary sclerosing cholangitis and Takayasu's arteritis3-6-we identify that the causal gene, ETS2, is a central regulator of human inflammatory macrophages and delineate the shared disease mechanism that amplifies ETS2 expression. Genes regulated by ETS2 were prominently expressed in diseased tissues and more enriched for inflammatory bowel disease GWAS hits than most previously described pathways. Overexpressing ETS2 in resting macrophages reproduced the inflammatory state observed in chr21q22-associated diseases, with upregulation of multiple drug targets, including TNF and IL-23. Using a database of cellular signatures7, we identified drugs that might modulate this pathway and validated the potent anti-inflammatory activity of one class of small molecules in vitro and ex vivo. Together, this illustrates the power of functional genomics, applied directly in primary human cells, to identify immune-mediated disease mechanisms and potential therapeutic opportunities.
Collapse
Affiliation(s)
- C T Stankey
- Genetic Mechanisms of Disease Laboratory, The Francis Crick Institute, London, UK
- Department of Immunology and Inflammation, Imperial College London, London, UK
- Washington University School of Medicine, St Louis, MO, USA
| | - C Bourges
- Genetic Mechanisms of Disease Laboratory, The Francis Crick Institute, London, UK
| | - L M Haag
- Division of Gastroenterology, Infectious Diseases and Rheumatology, Charité-Universitätsmedizin Berlin, Berlin, Germany
| | - T Turner-Stokes
- Genetic Mechanisms of Disease Laboratory, The Francis Crick Institute, London, UK
- Department of Immunology and Inflammation, Imperial College London, London, UK
| | - A P Piedade
- Genetic Mechanisms of Disease Laboratory, The Francis Crick Institute, London, UK
| | - C Palmer-Jones
- Department of Gastroenterology, Royal Free Hospital, London, UK
- Institute for Liver and Digestive Health, Division of Medicine, University College London, London, UK
| | - I Papa
- Genetic Mechanisms of Disease Laboratory, The Francis Crick Institute, London, UK
| | | | - Q Zhang
- Genomics of Inflammation and Immunity Group, Human Genetics Programme, Wellcome Sanger Institute, Hinxton, UK
| | - A J Cameron
- Wolfson Wohl Cancer Centre, School of Cancer Sciences, University of Glasgow, Glasgow, UK
| | - A Legrini
- Wolfson Wohl Cancer Centre, School of Cancer Sciences, University of Glasgow, Glasgow, UK
| | - T Zhang
- Wolfson Wohl Cancer Centre, School of Cancer Sciences, University of Glasgow, Glasgow, UK
| | - C S Wood
- Wolfson Wohl Cancer Centre, School of Cancer Sciences, University of Glasgow, Glasgow, UK
| | - F N New
- NanoString Technologies, Seattle, WA, USA
| | - L O Randzavola
- Department of Immunology and Inflammation, Imperial College London, London, UK
| | - L Speidel
- Ancient Genomics Laboratory, The Francis Crick Institute, London, UK
- Genetics Institute, University College London, London, UK
| | - A C Brown
- Wellcome Centre for Human Genetics, University of Oxford, Oxford, UK
| | - A Hall
- The Sheila Sherlock Liver Centre, Royal Free Hospital, London, UK
- Department of Cellular Pathology, Royal Free Hospital, London, UK
| | - F Saffioti
- Institute for Liver and Digestive Health, Division of Medicine, University College London, London, UK
- The Sheila Sherlock Liver Centre, Royal Free Hospital, London, UK
| | - E C Parkes
- Genetic Mechanisms of Disease Laboratory, The Francis Crick Institute, London, UK
| | - W Edwards
- Cambridge Institute of Therapeutic Immunology and Infectious Disease, University of Cambridge, Cambridge, UK
| | - H Direskeneli
- Department of Internal Medicine, Division of Rheumatology, Marmara University, Istanbul, Turkey
| | - P C Grayson
- Systemic Autoimmunity Branch, NIAMS, National Institutes of Health, Bethesda, MD, USA
| | - L Jiang
- Department of Rheumatology, Zhongshan Hospital, Fudan University, Shanghai, China
| | - P A Merkel
- Division of Rheumatology, Department of Medicine, University of Pennsylvania, Philadelphia, PA, USA
- Division of Epidemiology, Department of Biostatistics, Epidemiology and Informatics, University of Pennsylvania, Philadelphia, PA, USA
| | - G Saruhan-Direskeneli
- Department of Physiology, Istanbul University, Istanbul Faculty of Medicine, Istanbul, Turkey
| | - A H Sawalha
- Division of Rheumatology, Department of Pediatrics, University of Pittsburgh, Pittsburgh, PA, USA
- Division of Rheumatology and Clinical Immunology, Department of Medicine, University of Pittsburgh, Pittsburgh, PA, USA
- Lupus Center of Excellence, University of Pittsburgh, Pittsburgh, PA, USA
- Department of Immunology, University of Pittsburgh, Pittsburgh, PA, USA
| | - E Tombetti
- Department of Biomedical and Clinical Sciences, Milan University, Milan, Italy
- Internal Medicine and Rheumatology, ASST FBF-Sacco, Milan, Italy
| | - A Quaglia
- Department of Cellular Pathology, Royal Free Hospital, London, UK
- UCL Cancer Institute, London, UK
| | - D Thorburn
- Institute for Liver and Digestive Health, Division of Medicine, University College London, London, UK
- The Sheila Sherlock Liver Centre, Royal Free Hospital, London, UK
| | - J C Knight
- Wellcome Centre for Human Genetics, University of Oxford, Oxford, UK
- Chinese Academy of Medical Sciences Institute, Nuffield Department of Medicine, University of Oxford, Oxford, UK
- NIHR Comprehensive Biomedical Research Centre, Oxford, UK
| | - A P Rochford
- Department of Gastroenterology, Royal Free Hospital, London, UK
- Institute for Liver and Digestive Health, Division of Medicine, University College London, London, UK
| | - C D Murray
- Department of Gastroenterology, Royal Free Hospital, London, UK
- Institute for Liver and Digestive Health, Division of Medicine, University College London, London, UK
| | - P Divakar
- NanoString Technologies, Seattle, WA, USA
| | - M Green
- Experimental Histopathology STP, The Francis Crick Institute, London, UK
| | - E Nye
- Experimental Histopathology STP, The Francis Crick Institute, London, UK
| | - J I MacRae
- Metabolomics STP, The Francis Crick Institute, London, UK
| | - N B Jamieson
- Wolfson Wohl Cancer Centre, School of Cancer Sciences, University of Glasgow, Glasgow, UK
| | - P Skoglund
- Ancient Genomics Laboratory, The Francis Crick Institute, London, UK
| | - M Z Cader
- Cambridge Institute of Therapeutic Immunology and Infectious Disease, University of Cambridge, Cambridge, UK
- Department of Medicine, University of Cambridge, Cambridge, UK
| | - C Wallace
- Cambridge Institute of Therapeutic Immunology and Infectious Disease, University of Cambridge, Cambridge, UK
- MRC Biostatistics Unit, Cambridge Institute of Public Health, Cambridge, UK
| | - D C Thomas
- Cambridge Institute of Therapeutic Immunology and Infectious Disease, University of Cambridge, Cambridge, UK
- Department of Medicine, University of Cambridge, Cambridge, UK
| | - J C Lee
- Genetic Mechanisms of Disease Laboratory, The Francis Crick Institute, London, UK.
- Department of Gastroenterology, Royal Free Hospital, London, UK.
- Institute for Liver and Digestive Health, Division of Medicine, University College London, London, UK.
| |
Collapse
|
10
|
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
|
11
|
Yee SW, Ferrández-Peral L, Alentorn-Moron P, Fontsere C, Ceylan M, Koleske ML, Handin N, Artegoitia VM, Lara G, Chien HC, Zhou X, Dainat J, Zalevsky A, Sali A, Brand CM, Wolfreys FD, Yang J, Gestwicki JE, Capra JA, Artursson P, Newman JW, Marquès-Bonet T, Giacomini KM. Illuminating the function of the orphan transporter, SLC22A10, in humans and other primates. Nat Commun 2024; 15:4380. [PMID: 38782905 PMCID: PMC11116522 DOI: 10.1038/s41467-024-48569-7] [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: 08/14/2023] [Accepted: 05/06/2024] [Indexed: 05/25/2024] Open
Abstract
SLC22A10 is an orphan transporter with unknown substrates and function. The goal of this study is to elucidate its substrate specificity and functional characteristics. In contrast to orthologs from great apes, human SLC22A10, tagged with green fluorescent protein, is not expressed on the plasma membrane. Cells expressing great ape SLC22A10 orthologs exhibit significant accumulation of estradiol-17β-glucuronide, unlike those expressing human SLC22A10. Sequence alignments reveal a proline at position 220 in humans, which is a leucine in great apes. Replacing proline with leucine in SLC22A10-P220L restores plasma membrane localization and uptake function. Neanderthal and Denisovan genomes show proline at position 220, akin to modern humans, indicating functional loss during hominin evolution. Human SLC22A10 is a unitary pseudogene due to a fixed missense mutation, P220, while in great apes, its orthologs transport sex steroid conjugates. Characterizing SLC22A10 across species sheds light on its biological role, influencing organism development and steroid homeostasis.
Collapse
Affiliation(s)
- Sook Wah Yee
- Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, CA, USA
| | - Luis Ferrández-Peral
- Institute of Evolutionary Biology (UPF-CSIC), PRBB, Dr. Aiguader 88, 08003, Barcelona, Spain
| | - Pol Alentorn-Moron
- Institute of Evolutionary Biology (UPF-CSIC), PRBB, Dr. Aiguader 88, 08003, Barcelona, Spain
| | - Claudia Fontsere
- Institute of Evolutionary Biology (UPF-CSIC), PRBB, Dr. Aiguader 88, 08003, Barcelona, Spain
- Center for Evolutionary Hologenomics, The Globe Institute, University of Copenhagen, Øster Farimagsgade 5A, 1352, Copenhagen, Denmark
| | - Merve Ceylan
- Department of Pharmacy, Uppsala University, Uppsala, Sweden
| | - Megan L Koleske
- Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, CA, USA
| | - Niklas Handin
- Department of Pharmacy, Uppsala University, Uppsala, Sweden
| | - Virginia M Artegoitia
- United States Department of Agriculture, Agricultural Research Service, Western Human Nutrition Research Center, Davis, CA, 95616, USA
| | - Giovanni Lara
- Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, CA, USA
| | - Huan-Chieh Chien
- Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, CA, USA
| | - Xujia Zhou
- Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, CA, USA
| | - Jacques Dainat
- Joint Research Unit for Infectious Diseases and Vectors Ecology Genetics Evolution and Control (MIVEGEC), University of Montpellier, French National Center for Scientific Research (CNRS 5290), French National Research Institute for Sustainable Development (IRD 224), 911 Avenue Agropolis, BP 64501, 34394, Montpellier Cedex 5, France
| | - Arthur Zalevsky
- Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, CA, USA
| | - Andrej Sali
- Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, CA, USA
- Department of Pharmaceutical Chemistry, University of California, San Francisco, CA, USA
- Quantitative Biosciences Institute (QBI), University of California, San Francisco, San Francisco, CA, US
| | - 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
| | - Finn D Wolfreys
- Department of Ophthalmology, University of California, San Francisco, CA, USA
| | - Jia Yang
- Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, CA, USA
| | - Jason E Gestwicki
- Department of Pharmaceutical Chemistry, University of California, San Francisco, CA, USA
- Institute for Neurodegenerative Diseases, University of California, San Francisco, CA, USA
| | - John A Capra
- Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, CA, USA
- Bakar Computational Health Sciences Institute, University of California, San Francisco, CA, USA
- Department of Epidemiology and Biostatistics, University of California, San Francisco, CA, USA
- Institute for Human Genetics, University of California San Francisco, San Francisco, CA, USA
| | - Per Artursson
- Department of Pharmacy, Uppsala University, Uppsala, Sweden
- Science for Life Laboratories, Uppsala University, Uppsala, Sweden
| | - John W Newman
- United States Department of Agriculture, Agricultural Research Service, Western Human Nutrition Research Center, Davis, CA, 95616, USA
- Department of Nutrition, University of California, Davis, Davis, CA, 95616, USA
| | - Tomàs Marquès-Bonet
- Institute of Evolutionary Biology (UPF-CSIC), PRBB, Dr. Aiguader 88, 08003, Barcelona, Spain
- Catalan Institution of Research and Advanced Studies (ICREA), Passeig de Lluís Companys, 23, 08010, Barcelona, Spain
- CNAG, Centro Nacional de Analisis Genomico, Barcelona Institute of Science and Technology (BIST), Baldiri i Reixac 4, 08028, Barcelona, Spain
- Institut Català de Paleontologia Miquel Crusafont, Universitat Autònoma de Barcelona, Edifici ICTA-ICP, c/ Columnes s/n, 08193, Cerdanyola del Vallès, Barcelona, Spain
| | - Kathleen M Giacomini
- Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, CA, USA.
- Institute for Human Genetics, University of California San Francisco, San Francisco, CA, USA.
| |
Collapse
|
12
|
Moya R, Wang X, Tsien RW, Maurano MT. Structural characterization of a polymorphic repeat at the CACNA1C schizophrenia locus. MEDRXIV : THE PREPRINT SERVER FOR HEALTH SCIENCES 2024:2024.03.05.24303780. [PMID: 38798557 PMCID: PMC11118589 DOI: 10.1101/2024.03.05.24303780] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/29/2024]
Abstract
Genetic variation within intron 3 of the CACNA1C calcium channel gene is associated with schizophrenia and bipolar disorder, but analysis of the causal variants and their effect is complicated by a nearby variable-number tandem repeat (VNTR). Here, we used 155 long-read genome assemblies from 78 diverse individuals to delineate the structure and population variability of the CACNA1C intron 3 VNTR. We categorized VNTR sequences into 7 Types of structural alleles using sequence differences among repeat units. Only 12 repeat units at the 5' end of the VNTR were shared across most Types, but several Types were related through a series of large and small duplications. The most diverged Types were rare and present only in individuals with African ancestry, but the multiallelic structural polymorphism Variable Region 2 was present across populations at different frequencies, consistent with expansion of the VNTR preceding the emergence of early hominins. VR2 was in complete linkage disequilibrium with fine-mapped schizophrenia variants (SNPs) from genome-wide association studies (GWAS). This risk haplotype was associated with decreased CACNA1C gene expression in brain tissues profiled by the GTEx project. Our work suggests that sequence variation within a human-specific VNTR affects gene expression, and provides a detailed characterization of new alleles at a flagship neuropsychiatric locus.
Collapse
Affiliation(s)
- Raquel Moya
- Institute for Systems Genetics, NYU School of Medicine, New York, NY 10016, USA
| | - Xiaohan Wang
- Neuroscience Institute, NYU School of Medicine, New York, NY 10016, USA
- Department of Neuroscience and Physiology, New York University, New York, NY 10016
| | - Richard W. Tsien
- Neuroscience Institute, NYU School of Medicine, New York, NY 10016, USA
- Department of Neuroscience and Physiology, New York University, New York, NY 10016
| | - Matthew T. Maurano
- Institute for Systems Genetics, NYU School of Medicine, New York, NY 10016, USA
- Department of Pathology, NYU School of Medicine, New York, NY 10016, USA
| |
Collapse
|
13
|
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
|
14
|
André M, Brucato N, Hudjasov G, Pankratov V, Yermakovich D, Montinaro F, Kreevan R, Kariwiga J, Muke J, Boland A, Deleuze JF, Meyer V, Evans N, Cox MP, Leavesley M, Dannemann M, Org T, Metspalu M, Mondal M, Ricaut FX. Positive selection in the genomes of two Papua New Guinean populations at distinct altitude levels. Nat Commun 2024; 15:3352. [PMID: 38688933 PMCID: PMC11061283 DOI: 10.1038/s41467-024-47735-1] [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: 01/10/2023] [Accepted: 04/08/2024] [Indexed: 05/02/2024] Open
Abstract
Highlanders and lowlanders of Papua New Guinea have faced distinct environmental stress, such as hypoxia and environment-specific pathogen exposure, respectively. In this study, we explored the top genomics regions and the candidate driver SNPs for selection in these two populations using newly sequenced whole-genomes of 54 highlanders and 74 lowlanders. We identified two candidate SNPs under selection - one in highlanders, associated with red blood cell traits and another in lowlanders, which is associated with white blood cell count - both potentially influencing the heart rate of Papua New Guineans in opposite directions. We also observed four candidate driver SNPs that exhibit linkage disequilibrium with an introgressed haplotype, highlighting the need to explore the possibility of adaptive introgression within these populations. This study reveals that the signatures of positive selection in highlanders and lowlanders of Papua New Guinea align closely with the challenges they face, which are specific to their environments.
Collapse
Affiliation(s)
- Mathilde André
- Estonian Biocentre, Institute of Genomics, University of Tartu, Riia 23b, 51010, Tartu, Tartumaa, Estonia
- Centre for Genomics, Evolution & Medicine, Institute of Genomics, University of Tartu, Riia 23b, 51010, Tartu, Tartumaa, Estonia
| | - Nicolas Brucato
- Centre de Recherche sur la Biodiversité et l'Environnement (CRBE), Université de Toulouse, CNRS, IRD, Toulouse INP, Université Toulouse 3 - Paul Sabatier (UT3), Toulouse, France
| | - Georgi Hudjasov
- Centre for Genomics, Evolution & Medicine, Institute of Genomics, University of Tartu, Riia 23b, 51010, Tartu, Tartumaa, Estonia
| | - Vasili Pankratov
- Centre for Genomics, Evolution & Medicine, Institute of Genomics, University of Tartu, Riia 23b, 51010, Tartu, Tartumaa, Estonia
| | - Danat Yermakovich
- Centre for Genomics, Evolution & Medicine, Institute of Genomics, University of Tartu, Riia 23b, 51010, Tartu, Tartumaa, Estonia
| | - Francesco Montinaro
- Estonian Biocentre, Institute of Genomics, University of Tartu, Riia 23b, 51010, Tartu, Tartumaa, Estonia
- Department of Biosciences, Biotechnology and the Environment, University of Bari, Bari, Italy
| | - Rita Kreevan
- Centre for Genomics, Evolution & Medicine, Institute of Genomics, University of Tartu, Riia 23b, 51010, Tartu, Tartumaa, Estonia
| | - Jason Kariwiga
- Strand of Anthropology, Sociology and Archaeology, School of Humanities and Social Sciences, University of Papua New Guinea, University 134, PO Box 320, National Capital District, Papua New Guinea
- School of Social Science, University of Queensland, St Lucia, QLD, Australia
| | - John Muke
- Social Research Institute Ltd, Port Moresby, Papua New Guinea
| | - Anne Boland
- Université Paris-Saclay, CEA, Centre National de Recherche en Génomique Humaine (CNRGH), 91057, Evry, France
| | - Jean-François Deleuze
- Université Paris-Saclay, CEA, Centre National de Recherche en Génomique Humaine (CNRGH), 91057, Evry, France
| | - Vincent Meyer
- Université Paris-Saclay, CEA, Centre National de Recherche en Génomique Humaine (CNRGH), 91057, Evry, France
| | - Nicholas Evans
- ARC Centre of Excellence for the Dynamics of Language, Coombs Building, Fellows Road, CHL, CAP, Australian National University, Canberra, ACT, Australia
| | - Murray P Cox
- School of Natural Sciences, Massey University, Palmerston North, New Zealand
- Department of Statistics, University of Auckland, Auckland, New Zealand
| | - Matthew Leavesley
- Strand of Anthropology, Sociology and Archaeology, School of Humanities and Social Sciences, University of Papua New Guinea, University 134, PO Box 320, National Capital District, Papua New Guinea
- College of Arts, Society and Education, James Cook University, P.O. Box 6811, Cairns, QLD, 4870, Australia
- ARC Centre of Excellence for Australian Biodiversity and Heritage, University of Wollongong, Wollongong, NSW, 2522, Australia
| | - Michael Dannemann
- Centre for Genomics, Evolution & Medicine, Institute of Genomics, University of Tartu, Riia 23b, 51010, Tartu, Tartumaa, Estonia
| | - Tõnis Org
- Centre for Genomics, Evolution & Medicine, Institute of Genomics, University of Tartu, Riia 23b, 51010, Tartu, Tartumaa, Estonia
| | - Mait Metspalu
- Estonian Biocentre, Institute of Genomics, University of Tartu, Riia 23b, 51010, Tartu, Tartumaa, Estonia
| | - Mayukh Mondal
- Centre for Genomics, Evolution & Medicine, Institute of Genomics, University of Tartu, Riia 23b, 51010, Tartu, Tartumaa, Estonia.
- Institute of Clinical Molecular Biology, Christian-Albrechts-Universität zu Kiel, 24118, Kiel, Germany.
| | - François-Xavier Ricaut
- Centre de Recherche sur la Biodiversité et l'Environnement (CRBE), Université de Toulouse, CNRS, IRD, Toulouse INP, Université Toulouse 3 - Paul Sabatier (UT3), Toulouse, France.
| |
Collapse
|
15
|
Voinescu CD, Mozere M, Genovese G, Downie ML, Gupta S, Gale DP, Bockenhauer D, Kleta R, Arcos-Burgos M, Stanescu HC. A Neanderthal haplotype introgressed into the human genome confers protection against membranous nephropathy. Kidney Int 2024; 105:791-798. [PMID: 38367960 DOI: 10.1016/j.kint.2024.01.017] [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: 01/19/2023] [Revised: 12/27/2023] [Accepted: 01/08/2024] [Indexed: 02/19/2024]
Abstract
Class 2 HLA and PLA2R1 alleles are exceptionally strong genetic risk factors for membranous nephropathy (MN), leading, through an unknown mechanism, to a targeted autoimmune response. Introgressed archaic haplotypes (introduced from an archaic human genome into the modern human genome) might influence phenotypes through gene dysregulation. Here, we investigated the genomic region surrounding the PLA2R1 gene. We reconstructed the phylogeny of Neanderthal and modern haplotypes in this region and calculated the probability of the observed clustering being the result of introgression or common descent. We imputed variants for the participants in our previous genome-wide association study and we compared the distribution of Neanderthal variants between MN cases and controls. The region associated with the lead MN risk locus in the PLA2R1 gene was confirmed and showed that, within a 507 kb region enriched in introgressed sequence, a stringently defined 105 kb haplotype, intersecting the coding regions for PLA2R1 and ITGB6, is inherited from Neanderthals. Thus, introgressed Neanderthal haplotypes overlapping PLA2R1 are differentially represented in MN cases and controls, with enrichment In controls suggesting a protective effect.
Collapse
Affiliation(s)
- Cătălin D Voinescu
- Centre for Genetics and Genomics, Department of Renal Medicine, UCL Division of Medicine, University College London, London, UK
| | - Monika Mozere
- Department of Human and Medical Genetics, Faculty of Medicine, Vilnius University, Vilnius, Lithuania
| | - Giulio Genovese
- Program in Medical and Population Genetics, Broad Institute of MIT and Harvard, Cambridge, Massachusetts, USA; Stanley Center for Psychiatric Research, Broad Institute of MIT and Harvard, Cambridge, Massachusetts, USA; Department of Genetics, Harvard Medical School, Boston, Massachusetts, USA
| | - Mallory L Downie
- Centre for Genetics and Genomics, Department of Renal Medicine, UCL Division of Medicine, University College London, London, UK
| | - Sanjana Gupta
- Centre for Genetics and Genomics, Department of Renal Medicine, UCL Division of Medicine, University College London, London, UK
| | - Daniel P Gale
- Centre for Genetics and Genomics, Department of Renal Medicine, UCL Division of Medicine, University College London, London, UK
| | - Detlef Bockenhauer
- Centre for Genetics and Genomics, Department of Renal Medicine, UCL Division of Medicine, University College London, London, UK
| | - Robert Kleta
- Centre for Genetics and Genomics, Department of Renal Medicine, UCL Division of Medicine, University College London, London, UK
| | - Mauricio Arcos-Burgos
- Grupo de Investigación en Psiquiatría, Departamento de Psiquiatría, Instituto de Investigaciones Médicas, Facultad de Medicina, Universidad de Antioquia, Medellín, Colombia
| | - Horia C Stanescu
- Centre for Genetics and Genomics, Department of Renal Medicine, UCL Division of Medicine, University College London, London, UK.
| |
Collapse
|
16
|
Vallini L, Zampieri C, Shoaee MJ, Bortolini E, Marciani G, Aneli S, Pievani T, Benazzi S, Barausse A, Mezzavilla M, Petraglia MD, Pagani L. The Persian plateau served as hub for Homo sapiens after the main out of Africa dispersal. Nat Commun 2024; 15:1882. [PMID: 38528002 DOI: 10.1038/s41467-024-46161-7] [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/05/2023] [Accepted: 02/16/2024] [Indexed: 03/27/2024] Open
Abstract
A combination of evidence, based on genetic, fossil and archaeological findings, indicates that Homo sapiens spread out of Africa between ~70-60 thousand years ago (kya). However, it appears that once outside of Africa, human populations did not expand across all of Eurasia until ~45 kya. The geographic whereabouts of these early settlers in the timeframe between ~70-60 to 45 kya has been difficult to reconcile. Here we combine genetic evidence and palaeoecological models to infer the geographic location that acted as the Hub for our species during the early phases of colonisation of Eurasia. Leveraging on available genomic evidence we show that populations from the Persian Plateau carry an ancestry component that closely matches the population that settled the Hub outside Africa. With the paleoclimatic data available to date, we built ecological models showing that the Persian Plateau was suitable for human occupation and that it could sustain a larger population compared to other West Asian regions, strengthening this claim.
Collapse
Affiliation(s)
| | - Carlo Zampieri
- Department of Biology, University of Padova, Padova, Italy
| | - Mohamed Javad Shoaee
- Department of Archaeology, Max Planck Institute for Geoanthropology, Jena, Germany
| | - Eugenio Bortolini
- Department of Cultural Heritage, University of Bologna, Bologna, Italy
| | - Giulia Marciani
- Department of Cultural Heritage, University of Bologna, Bologna, Italy
- Research Unit Prehistory and Anthropology, Department of Physical Sciences, Earth and Environment, University of Siena, Siena, Italy
| | - Serena Aneli
- Department of Public Health Sciences and Pediatrics, University of Turin, Turin, Italy
| | - Telmo Pievani
- Department of Biology, University of Padova, Padova, Italy
| | - Stefano Benazzi
- Department of Cultural Heritage, University of Bologna, Bologna, Italy
| | - Alberto Barausse
- Department of Biology, University of Padova, Padova, Italy
- Department of Industrial Engineering, University of Padova, Padova, Italy
| | | | - Michael D Petraglia
- Human Origins Program, Smithsonian Institution, Washington, DC, 20560, USA
- School of Social Science, The University of Queensland, Brisbane, QLD, Australia
- Australian Research Centre for Human Evolution, Griffith University, Brisbane, QLD, Australia
| | - Luca Pagani
- Department of Biology, University of Padova, Padova, Italy.
- Institute of Genomics, University of Tartu, Tartu, Estonia.
| |
Collapse
|
17
|
Chen Y, Yu XY, Xu SJ, Shi XQ, Zhang XX, Sun C. An indel introduced by Neanderthal introgression, rs3835124:ATTTATT > ATT, might contribute to prostate cancer risk by regulating PDK1 expression. Ann Hum Genet 2024; 88:126-137. [PMID: 37846608 DOI: 10.1111/ahg.12533] [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/28/2023] [Revised: 07/25/2023] [Accepted: 08/04/2023] [Indexed: 10/18/2023]
Abstract
INTRODUCTION Prostate cancer is one of the most common cancer types in males and rs12621278:A > G has been suggested to be associated with this disease by previous genome-wide association studies. One thousand genomes project data analysis indicated that rs12621278:A > G is within two long-core haplotypes. However, the origin, causal variant(s), and molecular function of these haplotypes were remaining unclear. MATERIALS AND METHODS Population genetics analysis and functional genomics work was performed for this locus. RESULTS Phylogeny analysis verified that the rare haplotype is derived from Neanderthal introgression. Genome annotation suggested that three genetic variants in the core haplotypes, rs116108611:G > A, rs139972066:AAAAAAAA > AAAAAAAAA, and rs3835124:ATTTATT > ATT, are located in functional regions. Luciferase assay indicated that rs139972066:AAAAAAAA > AAAAAAAAA and rs116108611:G > A are not able to alter ITGA6 (integrin alpha 6) and ITGA6 antisense RNA 1 expression, respectively. In contrast, rs3835124:ATTTATT > ATT can significantly influence PDK1 (pyruvate dehydrogenase kinase 1) expression, which was verified by expression quantitative trait locus analysis. This genetic variant can alter transcription factor cut like homeobox 1 interaction efficiency. The introgressed haplotype was observed to be subject to positive selection in East Asian populations. The molecular function of the haplotype suggested that Neanderthal should be with lower PDK1 expression and further different energy homeostasis from modern human. CONCLUSION This study provided new insight into the contribution of Neanderthal introgression to human phenotypes.
Collapse
Affiliation(s)
- Ying Chen
- College of Life Sciences, Shaanxi Normal University, Xi'an, Shaanxi, P. R. China
| | - Xin-Yi Yu
- College of Life Sciences, Shaanxi Normal University, Xi'an, Shaanxi, P. R. China
| | - Shuang-Jia Xu
- College of Life Sciences, Shaanxi Normal University, Xi'an, Shaanxi, P. R. China
| | - Xiao-Qian Shi
- College of Life Sciences, Shaanxi Normal University, Xi'an, Shaanxi, P. R. China
| | - Xin-Xin Zhang
- College of Life Sciences, Shaanxi Normal University, Xi'an, Shaanxi, P. R. China
| | - Chang Sun
- College of Life Sciences, Shaanxi Normal University, Xi'an, Shaanxi, P. R. China
| |
Collapse
|
18
|
Goltermann O, Alagöz G, Molz B, Fisher SE. Neuroimaging genomics as a window into the evolution of human sulcal organization. Cereb Cortex 2024; 34:bhae078. [PMID: 38466113 PMCID: PMC10926775 DOI: 10.1093/cercor/bhae078] [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: 10/23/2023] [Revised: 02/03/2024] [Accepted: 02/04/2024] [Indexed: 03/12/2024] Open
Abstract
Primate brain evolution has involved prominent expansions of the cerebral cortex, with largest effects observed in the human lineage. Such expansions were accompanied by fine-grained anatomical alterations, including increased cortical folding. However, the molecular bases of evolutionary alterations in human sulcal organization are not yet well understood. Here, we integrated data from recently completed large-scale neuroimaging genetic analyses with annotations of the human genome relevant to various periods and events in our evolutionary history. These analyses identified single-nucleotide polymorphism (SNP) heritability enrichments in fetal brain human-gained enhancer (HGE) elements for a number of sulcal structures, including the central sulcus, which is implicated in human hand dexterity. We zeroed in on a genomic region that harbors DNA variants associated with left central sulcus shape, an HGE element, and genetic loci involved in neurogenesis including ZIC4, to illustrate the value of this approach for probing the complex factors contributing to human sulcal evolution.
Collapse
Affiliation(s)
- Ole Goltermann
- Max Planck School of Cognition, Stephanstrasse 1a, 04103 Leipzig, Germany
- Language and Genetics Department, Max Planck Institute for Psycholinguistics, Wundtlaan 1, 6525 XD Nijmegen, The Netherlands
- Institute of Systems Neuroscience, University Medical Center Hamburg-Eppendorf, Martinistrasse 52, 20246 Hamburg, Germany
| | - Gökberk Alagöz
- Language and Genetics Department, Max Planck Institute for Psycholinguistics, Wundtlaan 1, 6525 XD Nijmegen, The Netherlands
| | - Barbara Molz
- Language and Genetics Department, Max Planck Institute for Psycholinguistics, Wundtlaan 1, 6525 XD Nijmegen, The Netherlands
| | - Simon E Fisher
- Language and Genetics Department, Max Planck Institute for Psycholinguistics, Wundtlaan 1, 6525 XD Nijmegen, The Netherlands
- Donders Institute for Brain, Cognition & Behaviour, Radboud University, Thomas van Aquinostraat 4, 6525 GD Nijmegen, The Netherlands
| |
Collapse
|
19
|
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
|
20
|
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
|
21
|
Mallick S, Micco A, Mah M, Ringbauer H, Lazaridis I, Olalde I, Patterson N, Reich D. The Allen Ancient DNA Resource (AADR) a curated compendium of ancient human genomes. Sci Data 2024; 11:182. [PMID: 38341426 PMCID: PMC10858950 DOI: 10.1038/s41597-024-03031-7] [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: 08/10/2023] [Accepted: 01/31/2024] [Indexed: 02/12/2024] Open
Abstract
More than two hundred papers have reported genome-wide data from ancient humans. While the raw data for the vast majority are fully publicly available testifying to the commitment of the paleogenomics community to open data, formats for both raw data and meta-data differ. There is thus a need for uniform curation and a centralized, version-controlled compendium that researchers can download, analyze, and reference. Since 2019, we have been maintaining the Allen Ancient DNA Resource (AADR), which aims to provide an up-to-date, curated version of the world's published ancient human DNA data, represented at more than a million single nucleotide polymorphisms (SNPs) at which almost all ancient individuals have been assayed. The AADR has gone through six public releases at the time of writing and review of this manuscript, and crossed the threshold of >10,000 individuals with published genome-wide ancient DNA data at the end of 2022. This note is intended as a citable descriptor of the AADR.
Collapse
Affiliation(s)
- Swapan Mallick
- Department of Genetics, Harvard Medical School, Boston, MA, 02115, USA.
- Broad Institute of MIT and Harvard, Cambridge, MA, 02142, USA.
- Howard Hughes Medical Institute, Boston, MA, 02115, USA.
| | - Adam Micco
- Department of Genetics, Harvard Medical School, Boston, MA, 02115, USA
- Howard Hughes Medical Institute, Boston, MA, 02115, USA
| | - Matthew Mah
- Department of Genetics, Harvard Medical School, Boston, MA, 02115, USA
- Howard Hughes Medical Institute, Boston, MA, 02115, USA
| | - Harald Ringbauer
- Department of Human Evolutionary Biology, Harvard University, Cambridge, MA, 02138, USA
- Max Planck Institute for Evolutionary Anthropology, Leipzig, 04103, Germany
| | - Iosif Lazaridis
- Department of Genetics, Harvard Medical School, Boston, MA, 02115, USA
- Department of Human Evolutionary Biology, Harvard University, Cambridge, MA, 02138, USA
| | - Iñigo Olalde
- Department of Genetics, Harvard Medical School, Boston, MA, 02115, USA
- BIOMICs Research Group, University of the Basque Country, 01006, Vitoria-Gasteiz, Spain
| | - Nick Patterson
- Broad Institute of MIT and Harvard, Cambridge, MA, 02142, USA
- Department of Human Evolutionary Biology, Harvard University, Cambridge, MA, 02138, USA
| | - David Reich
- Department of Genetics, Harvard Medical School, Boston, MA, 02115, USA.
- Broad Institute of MIT and Harvard, Cambridge, MA, 02142, USA.
- Howard Hughes Medical Institute, Boston, MA, 02115, USA.
- Department of Human Evolutionary Biology, Harvard University, Cambridge, MA, 02138, USA.
| |
Collapse
|
22
|
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
|
23
|
Kubat J, Paterson R, Patramanis I, Barker G, Demeter F, Filoux A, Kullmer O, Mackie M, Marques-Bonet T, Huong NTM, Tuan NA, Pheng S, Rippengal J, Schrenk F, Souksavatdy V, Tshen LT, Wattanapituksakul A, Wang W, Zanolli C, Cappellini E, Bacon AM. Geometric morphometrics and paleoproteomics enlighten the paleodiversity of Pongo. PLoS One 2023; 18:e0291308. [PMID: 38100471 PMCID: PMC10723683 DOI: 10.1371/journal.pone.0291308] [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/16/2023] [Accepted: 08/27/2023] [Indexed: 12/17/2023] Open
Abstract
Pleistocene Pongo teeth show substantial variation in size and morphology, fueling taxonomic debates about the paleodiversity of the genus. We investigated prominent features of the enamel-dentine-junction junction (EDJ)-phylogenetically informative internal structures-of 71 fossil Pongo lower molars from various sites by applying geometric morphometrics and conducted paleoproteomic analyses from enamel proteins to attempt to identify extinct orangutan species. Forty-three orangutan lower molars representing Pongo pygmaeus and Pongo abelii were included for comparison. The shape of the EDJ was analyzed by placing five landmarks on the tip of the main dentine horns, and 142 semilandmarks along the marginal ridges connecting the dentine horns. Paleoproteomic analyses were conducted on 15 teeth of Late Pleistocene Pongo using high-resolution tandem mass spectrometry. The geometric morphometric results show variations in EDJ shape regarding aspects of the height and position of the dentine horns and connecting ridges. Despite the issue of molar position and sample size, modern molars are distinguished from fossil counterparts by their elongated tooth outline and narrowly positioned dentine horns. Proteomic results show that neither a distinction of P. pygmaeus and P. abelii, nor a consistent allocation of fossil specimens to extant species is feasible. Based on the EDJ shape, the (late) Middle to Late Pleistocene Pongo samples from Vietnam share the same morphospace, supporting the previous allocation to P. devosi, although substantial overlap with Chinese fossils could also indicate close affinities with P. weidenreichi. The hypothesis that both species represent one chronospecies cannot be ruled out. Two fossil specimens, one from Tam Hay Marklot (Laos, Late Pleistocene), and another from Sangiran (Java, Early to Middle Pleistocene), along with some specimens within the Punung sample (Java), exhibit affinities with Pongo abelii. The Punung fossils might represent a mix of early Late Pleistocene and later specimens (terminal Pleistocene to Holocene) related to modern Pongo. The taxonomy and phylogeny of the complete Punung sample needs to be further investigated.
Collapse
Affiliation(s)
- Jülide Kubat
- CNRS, BABEL, Université Paris Cité, Paris, France
- Department of Palaeoanthropology, Senckenberg Research Institute and Natural History Museum Frankfurt, Frankfurt am Main, Germany
| | - Ryan Paterson
- Globe Institute, University of Copenhagen, Copenhagen, Denmark
| | | | - Graeme Barker
- McDonald Institute for Archaeological Research, University of Cambridge, Cambridge, United Kingdom
| | - Fabrice Demeter
- Lundbeck Foundation GeoGenetics Centre, Globe Institute, University of Copenhagen, Copenhagen, Denmark
- UMR 7206 Eco Anthropologie, Muséum National d’Histoire Naturelle, CNRS, Paris, France
| | - Arnaud Filoux
- Palaeontological Research and Education Centre, Mahasarakham University, Mahasarakham, Thailand
| | - Ottmar Kullmer
- Department of Palaeoanthropology, Senckenberg Research Institute and Natural History Museum Frankfurt, Frankfurt am Main, Germany
| | - Meaghan Mackie
- Globe Institute, University of Copenhagen, Copenhagen, Denmark
- The Novo Nordisk Foundation Center for Protein Research, University of Copenhagen, Copenhagen, Denmark
| | - Tomas Marques-Bonet
- Institute of Evolutionary Biology (UPF-CSIC), PRBB, Barcelona, Spain
- Catalan Institution of Research and Advanced Studies (ICREA), Passeig de Lluís Companys, Barcelona, Spain
- CNAG-CRG, Centre for Genomic Regulation (CRG), Barcelona Institute of Science and Technology (BIST), Barcelona, Spain
- Institut Català de Paleontologia Miquel Crusafont, Universitat Autònoma de Barcelona, Barcelona, Spain
| | - Nguyen Thi Mai Huong
- Anthropological and Palaeoenvironmental Department, Institute of Archaeology, Ha Noi, Vietnam
| | - Nguyen Anh Tuan
- Anthropological and Palaeoenvironmental Department, Institute of Archaeology, Ha Noi, Vietnam
| | - Sytha Pheng
- Royal University of Fine Arts, Phnom Penh, Cambodia
| | - Jessica Rippengal
- McDonald Institute for Archaeological Research, University of Cambridge, Cambridge, United Kingdom
| | - Friedemann Schrenk
- Department of Palaeoanthropology, Senckenberg Research Institute and Natural History Museum Frankfurt, Frankfurt am Main, Germany
| | - Viengkeo Souksavatdy
- Department of Heritage, Ministry of Information, Culture and Tourism, Vientiane, Lao People’s Democratic Republic
| | - Lim Tze Tshen
- Department of Geology, Faculty of Science, Universiti Malaya, Kuala Lumpur, Malaysia
| | | | - Wei Wang
- Institute of Cultural Heritage, Shandong University, Qingdao, China
| | | | | | | |
Collapse
|
24
|
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
|
25
|
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
|
26
|
Houldcroft CJ, Underdown S. Infectious disease in the Pleistocene: Old friends or old foes? AMERICAN JOURNAL OF BIOLOGICAL ANTHROPOLOGY 2023; 182:513-531. [PMID: 38006200 DOI: 10.1002/ajpa.24737] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/12/2022] [Revised: 03/01/2023] [Accepted: 03/14/2023] [Indexed: 11/26/2023]
Abstract
The impact of endemic and epidemic disease on humans has traditionally been seen as a comparatively recent historical phenomenon associated with the Neolithisation of human groups, an increase in population size led by sedentarism, and increasing contact with domesticated animals as well as species occupying opportunistic symbiotic and ectosymbiotic relationships with humans. The orthodox approach is that Neolithisation created the conditions for increasing population size able to support a reservoir of infectious disease sufficient to act as selective pressure. This orthodoxy is the result of an overly simplistic reliance on skeletal data assuming that no skeletal lesions equated to a healthy individual, underpinned by the assumption that hunter-gatherer groups were inherently healthy while agricultural groups acted as infectious disease reservoirs. The work of van Blerkom, Am. J. Phys. Anthropol., vol. suppl 37 (2003), Wolfe et al., Nature, vol. 447 (2007) and Houldcroft and Underdown, Am. J. Phys. Anthropol., vol. 160, (2016) has changed this landscape by arguing that humans and pathogens have long been fellow travelers. The package of infectious diseases experienced by our ancient ancestors may not be as dissimilar to modern infectious diseases as was once believed. The importance of DNA, from ancient and modern sources, to the study of the antiquity of infectious disease, and its role as a selective pressure cannot be overstated. Here we consider evidence of ancient epidemic and endemic infectious diseases with inferences from modern and ancient human and hominin DNA, and from circulating and extinct pathogen genomes. We argue that the pandemics of the past are a vital tool to unlock the weapons needed to fight pandemics of the future.
Collapse
Affiliation(s)
| | - Simon Underdown
- Human Origins and Palaeoenvironmental Research Group, School of Social Sciences, Oxford Brookes University, Oxford, UK
- Center for Microbial Ecology and Genomics, Department of Biochemistry, Genetics and Microbiology, University of Pretoria, Pretoria, South Africa
| |
Collapse
|
27
|
Harris DN, Platt A, Hansen MEB, Fan S, McQuillan MA, Nyambo T, Mpoloka SW, Mokone GG, Belay G, Fokunang C, Njamnshi AK, Tishkoff SA. Diverse African genomes reveal selection on ancient modern human introgressions in Neanderthals. Curr Biol 2023; 33:4905-4916.e5. [PMID: 37837965 PMCID: PMC10841429 DOI: 10.1016/j.cub.2023.09.066] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2023] [Revised: 07/18/2023] [Accepted: 09/26/2023] [Indexed: 10/16/2023]
Abstract
Comparisons of Neanderthal genomes to anatomically modern human (AMH) genomes show a history of Neanderthal-to-AMH introgression stemming from interbreeding after the migration of AMHs from Africa to Eurasia. All non-sub-Saharan African AMHs have genomic regions genetically similar to Neanderthals that descend from this introgression. Regions of the genome with Neanderthal similarities have also been identified in sub-Saharan African populations, but their origins have been unclear. To better understand how these regions are distributed across sub-Saharan Africa, the source of their origin, and what their distribution within the genome tells us about early AMH and Neanderthal evolution, we analyzed a dataset of high-coverage, whole-genome sequences from 180 individuals from 12 diverse sub-Saharan African populations. In sub-Saharan African populations with non-sub-Saharan African ancestry, as much as 1% of their genomes can be attributed to Neanderthal sequence introduced by recent migration, and subsequent admixture, of AMH populations originating from the Levant and North Africa. However, most Neanderthal homologous regions in sub-Saharan African populations originate from migration of AMH populations from Africa to Eurasia ∼250 kya, and subsequent admixture with Neanderthals, resulting in ∼6% AMH ancestry in Neanderthals. These results indicate that there have been multiple migration events of AMHs out of Africa and that Neanderthal and AMH gene flow has been bi-directional. Observing that genomic regions where AMHs show a depletion of Neanderthal introgression are also regions where Neanderthal genomes show a depletion of AMH introgression points to deleterious interactions between introgressed variants and background genomes in both groups-a hallmark of incipient speciation.
Collapse
Affiliation(s)
- Daniel N Harris
- Department of Genetics, University of Pennsylvania, Philadelphia, PA 19104, USA
| | - Alexander Platt
- Department of Genetics, University of Pennsylvania, Philadelphia, PA 19104, USA
| | - Matthew E B Hansen
- Department of Genetics, University of Pennsylvania, Philadelphia, PA 19104, USA
| | - Shaohua Fan
- State Key Laboratory of Genetic Engineering, Human Phenome Institute, Zhangjiang Fudan International Innovation Center, School of Life Science, Fudan University, Shanghai 200438, China
| | - Michael A McQuillan
- Department of Genetics, University of Pennsylvania, Philadelphia, PA 19104, USA
| | - Thomas Nyambo
- Department of Biochemistry and Molecular Biology, Hubert Kairuki Memorial University, Dar es Salaam, Tanzania
| | - Sununguko Wata Mpoloka
- Department of Biological Sciences, Faculty of Science, University of Botswana, Private Bag UB 0022, Gaborone, Botswana
| | - Gaonyadiwe George Mokone
- Department of Biomedical Sciences, Faculty of Medicine, University of Botswana, Private Bag UB 0022, Gaborone, Botswana
| | - Gurja Belay
- Department of Microbial Cellular and Molecular Biology, Addis Ababa University, P.O. Box 1176, Addis Ababa, Ethiopia
| | - Charles Fokunang
- Department of Pharmacotoxicology and Pharmacokinetics, Faculty of Medicine and Biomedical Sciences, The University of Yaoundé I, P.O. Box 337, Yaoundé, Cameroon
| | - Alfred K Njamnshi
- Brain Research Africa Initiative (BRAIN), P.O. Box 25625, Yaoundé, Cameroon; Neuroscience Lab, Faculty of Medicine and Biomedical Sciences, The University of Yaoundé I, Yaoundé, Cameroon
| | - Sarah A Tishkoff
- Department of Genetics, University of Pennsylvania, Philadelphia, PA 19104, USA; Department of Biology, University of Pennsylvania, Philadelphia, PA 19104, USA.
| |
Collapse
|
28
|
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
|
29
|
Di Pietro L, Boroumand M, Lattanzi W, Manconi B, Salvati M, Cabras T, Olianas A, Flore L, Serrao S, Calò CM, Francalacci P, Parolini O, Castagnola M. A Catalog of Coding Sequence Variations in Salivary Proteins' Genes Occurring during Recent Human Evolution. Int J Mol Sci 2023; 24:15010. [PMID: 37834461 PMCID: PMC10573131 DOI: 10.3390/ijms241915010] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/17/2023] [Revised: 10/04/2023] [Accepted: 10/06/2023] [Indexed: 10/15/2023] Open
Abstract
Saliva houses over 2000 proteins and peptides with poorly clarified functions, including proline-rich proteins, statherin, P-B peptides, histatins, cystatins, and amylases. Their genes are poorly conserved across related species, reflecting an evolutionary adaptation. We searched the nucleotide substitutions fixed in these salivary proteins' gene loci in modern humans compared with ancient hominins. We mapped 3472 sequence variants/nucleotide substitutions in coding, noncoding, and 5'-3' untranslated regions. Despite most of the detected variations being within noncoding regions, the frequency of coding variations was far higher than the general rate found throughout the genome. Among the various missense substitutions, specific substitutions detected in PRB1 and PRB2 genes were responsible for the introduction/abrogation of consensus sequences recognized by convertase enzymes that cleave the protein precursors. Overall, these changes that occurred during the recent human evolution might have generated novel functional features and/or different expression ratios among the various components of the salivary proteome. This may have influenced the homeostasis of the oral cavity environment, possibly conditioning the eating habits of modern humans. However, fixed nucleotide changes in modern humans represented only 7.3% of all the substitutions reported in this study, and no signs of evolutionary pressure or adaptative introgression from archaic hominins were found on the tested genes.
Collapse
Affiliation(s)
- Lorena Di Pietro
- Dipartimento Scienze della Vita e Sanità Pubblica, Università Cattolica del Sacro Cuore, 00168 Rome, Italy; (L.D.P.)
- Fondazione Policlinico Universitario A. Gemelli IRCCS, 00168 Rome, Italy
| | - Mozhgan Boroumand
- Laboratorio di Proteomica, Centro Europeo di Ricerca sul Cervello, IRCCS Fondazione Santa Lucia, 00179 Rome, Italy
| | - Wanda Lattanzi
- Dipartimento Scienze della Vita e Sanità Pubblica, Università Cattolica del Sacro Cuore, 00168 Rome, Italy; (L.D.P.)
- Fondazione Policlinico Universitario A. Gemelli IRCCS, 00168 Rome, Italy
| | - Barbara Manconi
- Dipartimento di Scienze della Vita e Dell’ambiente, Università di Cagliari, 09042 Monserrato, Italy
| | - Martina Salvati
- Dipartimento Scienze della Vita e Sanità Pubblica, Università Cattolica del Sacro Cuore, 00168 Rome, Italy; (L.D.P.)
| | - Tiziana Cabras
- Dipartimento di Scienze della Vita e Dell’ambiente, Università di Cagliari, 09042 Monserrato, Italy
| | - Alessandra Olianas
- Dipartimento di Scienze della Vita e Dell’ambiente, Università di Cagliari, 09042 Monserrato, Italy
| | - Laura Flore
- Dipartimento di Scienze della Vita e Dell’ambiente, Università di Cagliari, 09042 Monserrato, Italy
| | - Simone Serrao
- Department of Medicine and Surgery, Proteomics and Metabolomics Unit, University of Milano-Bicocca, 20854 Vedano al Lambro, Italy
| | - Carla M. Calò
- Dipartimento di Scienze della Vita e Dell’ambiente, Università di Cagliari, 09042 Monserrato, Italy
| | - Paolo Francalacci
- Dipartimento di Scienze della Vita e Dell’ambiente, Università di Cagliari, 09042 Monserrato, Italy
| | - Ornella Parolini
- Dipartimento Scienze della Vita e Sanità Pubblica, Università Cattolica del Sacro Cuore, 00168 Rome, Italy; (L.D.P.)
- Fondazione Policlinico Universitario A. Gemelli IRCCS, 00168 Rome, Italy
| | - Massimo Castagnola
- Laboratorio di Proteomica, Centro Europeo di Ricerca sul Cervello, IRCCS Fondazione Santa Lucia, 00179 Rome, Italy
| |
Collapse
|
30
|
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
|
31
|
Yee SW, Ferrández-Peral L, Alentorn P, Fontsere C, Ceylan M, Koleske ML, Handin N, Artegoitia VM, Lara G, Chien HC, Zhou X, Dainat J, Zalevsky A, Sali A, Brand CM, Capra JA, Artursson P, Newman JW, Marques-Bonet T, Giacomini KM. Illuminating the Function of the Orphan Transporter, SLC22A10 in Humans and Other Primates. RESEARCH SQUARE 2023:rs.3.rs-3263845. [PMID: 37790518 PMCID: PMC10543398 DOI: 10.21203/rs.3.rs-3263845/v1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/05/2023]
Abstract
SLC22A10 is classified as an orphan transporter with unknown substrates and function. Here we describe the discovery of the substrate specificity and functional characteristics of SLC22A10. The human SLC22A10 tagged with green fluorescent protein was found to be absent from the plasma membrane, in contrast to the SLC22A10 orthologs found in great apes. Estradiol-17β-glucuronide accumulated in cells expressing great ape SLC22A10 orthologs (over 4-fold, p<0.001). In contrast, human SLC22A10 displayed no uptake function. Sequence alignments revealed two amino acid differences including a proline at position 220 of the human SLC22A10 and a leucine at the same position of great ape orthologs. Site-directed mutagenesis yielding the human SLC22A10-P220L produced a protein with excellent plasma membrane localization and associated uptake function. Neanderthal and Denisovan genomes show human-like sequences at proline 220 position, corroborating that SLC22A10 were rendered nonfunctional during hominin evolution after the divergence from the pan lineage (chimpanzees and bonobos). These findings demonstrate that human SLC22A10 is a unitary pseudogene and was inactivated by a missense mutation that is fixed in humans, whereas orthologs in great apes transport sex steroid conjugates.
Collapse
Affiliation(s)
- Sook Wah Yee
- Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, California, USA
| | | | - Pol Alentorn
- Institute of Evolutionary Biology (UPF-CSIC), PRBB, 08003 Barcelona, Spain
| | - Claudia Fontsere
- Institute of Evolutionary Biology (UPF-CSIC), PRBB, 08003 Barcelona, Spain; Center for Evolutionary Hologenomics, The Globe Institute, University of Copenhagen, Øster Farimagsgade 5A, 1352 Copenhagen, Denmark
| | - Merve Ceylan
- Department of Pharmacy and Science for Life Laboratory, Uppsala University, P.O. Box 580, 75123, Uppsala, Sweden
| | - Megan L. Koleske
- Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, California, USA
| | - Niklas Handin
- Department of Pharmacy and Science for Life Laboratory, Uppsala University, P.O. Box 580, 75123, Uppsala, Sweden
| | - Virginia M. Artegoitia
- United States Department of Agriculture, Agricultural Research Service, Western Human Nutrition Research Center, Davis, CA 95616, USA
| | - Giovanni Lara
- Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, California, USA
| | - Huan-Chieh Chien
- Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, California, USA
| | - Xujia Zhou
- Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, California, USA
| | - Jacques Dainat
- Joint Research Unit for Infectious Diseases and Vectors Ecology Genetics Evolution and Control (MIVEGEC), University of Montpellier, French National Center for Scientific Research (CNRS 5290), French National Research Institute for Sustainable Development (IRD 224), 911 Avenue Agropolis, BP 64501, 34394 Montpellier Cedex 5, France
| | - Arthur Zalevsky
- Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, California, USA
| | - Andrej Sali
- Department of Bioengineering and Therapeutic Sciences, UCSF Box 0775 1700 4th St, University of California, San Francisco, San Francisco, CA 94158, United States; Department of Pharmaceutical Chemistry, University of California, San Francisco, UCSF Box 2880 600 16th St, San Francisco, CA 94143, United States; Quantitative Biosciences Institute (QBI), University of California, San Francisco, 1700 4th St, San Francisco, CA, United States
| | - 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
| | - 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
| | - Per Artursson
- Department of Pharmacy and Science for Life Laboratory, Uppsala University, P.O. Box 580, 75123, Uppsala, Sweden
| | - John W. Newman
- United States Department of Agriculture, Agricultural Research Service, Western Human Nutrition Research Center, Davis, CA 95616, USA; Department of Nutrition, University of California, Davis, Davis, CA 95616, USA; UC Davis West Coast Metabolomics Center, Davis, CA 95616, USA
| | - Tomas Marques-Bonet
- Institute of Evolutionary Biology (UPF-CSIC), PRBB, 08003 Barcelona, Spain; Institute of Evolutionary Biology (UPF-CSIC), PRBB, Dr. Aiguader 88, 08003 Barcelona, Spain; Catalan Institution of Research and Advanced Studies (ICREA), Passeig de Lluís Companys, 23, 08010, Barcelona, Spain; CNAG, Centro Nacional de Analisis Genomico, Barcelona Institute of Science and Technology (BIST), Baldiri i Reixac 4, 08028 Barcelona, Spain; Institut Català de Paleontologia Miquel Crusafont, Universitat Autònoma de Barcelona, Edifici ICTA-ICP, c/ Columnes s/n, 08193 Cerdanyola del Vallès, Barcelona, Spain
| | - Kathleen M. Giacomini
- Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, California, USA
- Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, California, USA; Institute for Human Genetics, University of California San Francisco, San Francisco, CA, USA
| |
Collapse
|
32
|
Glunčić M, Vlahović I, Rosandić M, Paar V. Tandem NBPF 3mer HORs (Olduvai triplets) in Neanderthal and two novel HOR tandem arrays in human chromosome 1 T2T-CHM13 assembly. Sci Rep 2023; 13:14420. [PMID: 37660151 PMCID: PMC10475015 DOI: 10.1038/s41598-023-41517-3] [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: 05/24/2023] [Accepted: 08/28/2023] [Indexed: 09/04/2023] Open
Abstract
It is known that the ~ 1.6 kb Neuroblastoma BreakPoint Family (NBPF) repeats are human specific and contributing to cognitive capabilities, with increasing frequency in higher order repeat 3mer HORs (Olduvai triplets). From chimpanzee to modern human there is a discontinuous jump from 0 to ~ 50 tandemly organized 3mer HORs. Here we investigate the structure of NBPF 3mer HORs in the Neanderthal genome assembly of Pääbo et al., comparing it to the results obtained for human hg38.p14 chromosome 1. Our findings reveal corresponding NBPF 3mer HOR arrays in Neanderthals with slightly different monomer structures and numbers of HOR copies compared to humans. Additionally, we compute the NBPF 3mer HOR pattern for the complete telomere-to-telomere human genome assembly (T2T-CHM13) by Miga et al., identifying two novel tandem arrays of NBPF 3mer HOR repeats with 5 and 9 NBPF 3mer HOR copies. We hypothesize that these arrays correspond to novel NBPF genes (here referred to as NBPFA1 and NBPFA2). Further improving the quality of the Neanderthal genome using T2T-CHM13 as a reference would be of great interest in determining the presence of such distant novel NBPF genes in the Neanderthal genome and enhancing our understanding of human evolution.
Collapse
Affiliation(s)
- Matko Glunčić
- Faculty of Science, University of Zagreb, 10000, Zagreb, Croatia.
| | | | - Marija Rosandić
- University Hospital Centre Zagreb (Ret.), 10000, Zagreb, Croatia
- Croatian Academy of Sciences and Arts, 10000, Zagreb, Croatia
| | - Vladimir Paar
- Faculty of Science, University of Zagreb, 10000, Zagreb, Croatia
- Croatian Academy of Sciences and Arts, 10000, Zagreb, Croatia
| |
Collapse
|
33
|
Flegontov P, Işıldak U, Maier R, Yüncü E, Changmai P, Reich D. Modeling of African population history using f-statistics is biased when applying all previously proposed SNP ascertainment schemes. PLoS Genet 2023; 19:e1010931. [PMID: 37676865 PMCID: PMC10508636 DOI: 10.1371/journal.pgen.1010931] [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: 01/22/2023] [Revised: 09/19/2023] [Accepted: 08/21/2023] [Indexed: 09/09/2023] Open
Abstract
f-statistics have emerged as a first line of analysis for making inferences about demographic history from genome-wide data. Not only are they guaranteed to allow robust tests of the fits of proposed models of population history to data when analyzing full genome sequencing data-that is, all single nucleotide polymorphisms (SNPs) in the individuals being analyzed-but they are also guaranteed to allow robust tests of models for SNPs ascertained as polymorphic in a population that is an outgroup in a phylogenetic sense to all groups being analyzed. True "outgroup ascertainment" is in practice impossible in humans because our species has arisen from a substructured ancestral population that does not descend from a homogeneous ancestral population going back many hundreds of thousands of years into the past. However, initial studies suggested that non-outgroup-ascertainment schemes might produce robust enough results using f-statistics, and that motivated widespread fitting of models to data using non-outgroup-ascertained SNP panels such as the "Affymetrix Human Origins array" which has been genotyped on thousands of modern individuals from hundreds of populations, or the "1240k" in-solution enrichment reagent which has been the source of about 70% of published genome-wide data for ancient humans. In this study, we show that while analyses of population history using such panels work well for studies of relationships among non-African populations and one African outgroup, when co-modeling more than one sub-Saharan African and/or archaic human groups (Neanderthals and Denisovans), fitting of f-statistics to such SNP sets is expected to frequently lead to false rejection of true demographic histories, and failure to reject incorrect models. Analyzing panels of SNPs polymorphic in archaic humans, which has been suggested as a solution for the ascertainment problem, has limited statistical power and retains important biases. However, by carrying out simulations of diverse demographic histories, we show that bias in inferences based on f-statistics can be minimized by ascertaining on variants common in a union of diverse African groups; such ascertainment retains high statistical power while allowing co-analysis of archaic and modern groups.
Collapse
Affiliation(s)
- Pavel Flegontov
- Department of Human Evolutionary Biology, Harvard University, Cambridge, Massachusetts, United States of America
- Department of Biology and Ecology, Faculty of Science, University of Ostrava, Ostrava, Czechia
- Kalmyk Research Center of the Russian Academy of Sciences, Elista, Russia
| | - Ulaş Işıldak
- Department of Biology and Ecology, Faculty of Science, University of Ostrava, Ostrava, Czechia
| | - Robert Maier
- Department of Human Evolutionary Biology, Harvard University, Cambridge, Massachusetts, United States of America
| | - Eren Yüncü
- Department of Biology and Ecology, Faculty of Science, University of Ostrava, Ostrava, Czechia
| | - Piya Changmai
- Department of Biology and Ecology, Faculty of Science, University of Ostrava, Ostrava, Czechia
| | - David Reich
- Department of Human Evolutionary Biology, Harvard University, Cambridge, Massachusetts, United States of America
- Department of Genetics, Harvard Medical School, Boston, Massachusetts, United States of America
- Howard Hughes Medical Institute, Harvard Medical School, Boston, Massachusetts, United States of America
- Broad Institute of Harvard and MIT, Cambridge, Massachusetts, United States of America
| |
Collapse
|
34
|
Yee SW, Ferrández-Peral L, Alentorn P, Fontsere C, Ceylan M, Koleske ML, Handin N, Artegoitia VM, Lara G, Chien HC, Zhou X, Dainat J, Zalevsky A, Sali A, Brand CM, Capra JA, Artursson P, Newman JW, Marques-Bonet T, Giacomini KM. Illuminating the Function of the Orphan Transporter, SLC22A10 in Humans and Other Primates. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.08.08.552553. [PMID: 37609337 PMCID: PMC10441401 DOI: 10.1101/2023.08.08.552553] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 08/24/2023]
Abstract
SLC22A10 is classified as an orphan transporter with unknown substrates and function. Here we describe the discovery of the substrate specificity and functional characteristics of SLC22A10. The human SLC22A10 tagged with green fluorescent protein was found to be absent from the plasma membrane, in contrast to the SLC22A10 orthologs found in great apes. Estradiol-17β-glucuronide accumulated in cells expressing great ape SLC22A10 orthologs (over 4-fold, p<0.001). In contrast, human SLC22A10 displayed no uptake function. Sequence alignments revealed two amino acid differences including a proline at position 220 of the human SLC22A10 and a leucine at the same position of great ape orthologs. Site-directed mutagenesis yielding the human SLC22A10-P220L produced a protein with excellent plasma membrane localization and associated uptake function. Neanderthal and Denisovan genomes show human-like sequences at proline 220 position, corroborating that SLC22A10 were rendered nonfunctional during hominin evolution after the divergence from the pan lineage (chimpanzees and bonobos). These findings demonstrate that human SLC22A10 is a unitary pseudogene and was inactivated by a missense mutation that is fixed in humans, whereas orthologs in great apes transport sex steroid conjugates.
Collapse
Affiliation(s)
- Sook Wah Yee
- Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, California, USA
| | | | - Pol Alentorn
- Institute of Evolutionary Biology (UPF-CSIC), PRBB, 08003 Barcelona, Spain
| | - Claudia Fontsere
- Institute of Evolutionary Biology (UPF-CSIC), PRBB, 08003 Barcelona, Spain; Center for Evolutionary Hologenomics, The Globe Institute, University of Copenhagen, Øster Farimagsgade 5A, 1352 Copenhagen, Denmark
| | - Merve Ceylan
- Department of Pharmacy and Science for Life Laboratory, Uppsala University, P.O. Box 580, 75123, Uppsala, Sweden
| | - Megan L. Koleske
- Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, California, USA
| | - Niklas Handin
- Department of Pharmacy and Science for Life Laboratory, Uppsala University, P.O. Box 580, 75123, Uppsala, Sweden
| | - Virginia M. Artegoitia
- United States Department of Agriculture, Agricultural Research Service, Western Human Nutrition Research Center, Davis, CA 95616, USA
| | - Giovanni Lara
- Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, California, USA
| | - Huan-Chieh Chien
- Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, California, USA
| | - Xujia Zhou
- Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, California, USA
| | - Jacques Dainat
- Joint Research Unit for Infectious Diseases and Vectors Ecology Genetics Evolution and Control (MIVEGEC), University of Montpellier, French National Center for Scientific Research (CNRS 5290), French National Research Institute for Sustainable Development (IRD 224), 911 Avenue Agropolis, BP 64501, 34394 Montpellier Cedex 5, France
| | - Arthur Zalevsky
- Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, California, USA
| | - Andrej Sali
- Department of Bioengineering and Therapeutic Sciences, UCSF Box 0775 1700 4th St, University of California, San Francisco, San Francisco, CA 94158, United States; Department of Pharmaceutical Chemistry, University of California, San Francisco, UCSF Box 2880 600 16th St, San Francisco, CA 94143, United States; Quantitative Biosciences Institute (QBI), University of California, San Francisco, 1700 4th St, San Francisco, CA, United States
| | - 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
| | - 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
| | - Per Artursson
- Department of Pharmacy and Science for Life Laboratory, Uppsala University, P.O. Box 580, 75123, Uppsala, Sweden
| | - John W. Newman
- United States Department of Agriculture, Agricultural Research Service, Western Human Nutrition Research Center, Davis, CA 95616, USA; Department of Nutrition, University of California, Davis, Davis, CA 95616, USA; UC Davis West Coast Metabolomics Center, Davis, CA 95616, USA
| | - Tomas Marques-Bonet
- Institute of Evolutionary Biology (UPF-CSIC), PRBB, 08003 Barcelona, Spain; Institute of Evolutionary Biology (UPF-CSIC), PRBB, Dr. Aiguader 88, 08003 Barcelona, Spain; Catalan Institution of Research and Advanced Studies (ICREA), Passeig de Lluís Companys, 23, 08010, Barcelona, Spain; CNAG, Centro Nacional de Analisis Genomico, Barcelona Institute of Science and Technology (BIST), Baldiri i Reixac 4, 08028 Barcelona, Spain; Institut Català de Paleontologia Miquel Crusafont, Universitat Autònoma de Barcelona, Edifici ICTA-ICP, c/ Columnes s/n, 08193 Cerdanyola del Vallès, Barcelona, Spain
| | - Kathleen M. Giacomini
- Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, California, USA
- Department of Bioengineering and Therapeutic Sciences, University of California, San Francisco, California, USA; Institute for Human Genetics, University of California San Francisco, San Francisco, CA, USA
| |
Collapse
|
35
|
Caglayan E, Ayhan F, Liu Y, Vollmer RM, Oh E, Sherwood CC, Preuss TM, Yi SV, Konopka G. Molecular features driving cellular complexity of human brain evolution. Nature 2023; 620:145-153. [PMID: 37468639 PMCID: PMC11161302 DOI: 10.1038/s41586-023-06338-4] [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/24/2022] [Accepted: 06/16/2023] [Indexed: 07/21/2023]
Abstract
Human-specific genomic changes contribute to the unique functionalities of the human brain1-5. The cellular heterogeneity of the human brain6,7 and the complex regulation of gene expression highlight the need to characterize human-specific molecular features at cellular resolution. Here we analysed single-nucleus RNA-sequencing and single-nucleus assay for transposase-accessible chromatin with sequencing datasets for human, chimpanzee and rhesus macaque brain tissue from posterior cingulate cortex. We show a human-specific increase of oligodendrocyte progenitor cells and a decrease of mature oligodendrocytes across cortical tissues. Human-specific regulatory changes were accelerated in oligodendrocyte progenitor cells, and we highlight key biological pathways that may be associated with the proportional changes. We also identify human-specific regulatory changes in neuronal subtypes, which reveal human-specific upregulation of FOXP2 in only two of the neuronal subtypes. We additionally identify hundreds of new human accelerated genomic regions associated with human-specific chromatin accessibility changes. Our data also reveal that FOS::JUN and FOX motifs are enriched in the human-specifically accessible chromatin regions of excitatory neuronal subtypes. Together, our results reveal several new mechanisms underlying the evolutionary innovation of human brain at cell-type resolution.
Collapse
Affiliation(s)
- Emre Caglayan
- Department of Neuroscience, UT Southwestern Medical Center, Dallas, TX, USA
- Peter O'Donnell Jr. Brain Institute, UT Southwestern Medical Center, Dallas, TX, USA
| | - Fatma Ayhan
- Department of Neuroscience, UT Southwestern Medical Center, Dallas, TX, USA
- Peter O'Donnell Jr. Brain Institute, UT Southwestern Medical Center, Dallas, TX, USA
| | - Yuxiang Liu
- Department of Neuroscience, UT Southwestern Medical Center, Dallas, TX, USA
- Peter O'Donnell Jr. Brain Institute, UT Southwestern Medical Center, Dallas, TX, USA
| | - Rachael M Vollmer
- Department of Neuroscience, UT Southwestern Medical Center, Dallas, TX, USA
- Peter O'Donnell Jr. Brain Institute, UT Southwestern Medical Center, Dallas, TX, USA
| | - Emily Oh
- Department of Neuroscience, UT Southwestern Medical Center, Dallas, TX, USA
- Peter O'Donnell Jr. Brain Institute, UT Southwestern Medical Center, Dallas, TX, USA
| | - Chet C Sherwood
- Center for the Advanced Study of Human Paleobiology, Department of Anthropology, The George Washington University, Washington, DC, USA
| | - Todd M Preuss
- Division of Neuropharmacology and Neurologic Diseases, Yerkes National Primate Research Center, Emory University, Atlanta, GA, USA
- Department of Pathology, Emory University School of Medicine, Atlanta, GA, USA
| | - Soojin V Yi
- Department of Ecology, Evolution, and Marine Biology, University of California, Santa Barbara, Santa Barbara, CA, USA.
- Department of Molecular, Cellular and Developmental Biology, University of California, Santa Barbara, Santa Barbara, CA, USA.
- Neuroscience Research Institute, University of California, Santa Barbara, Santa Barbara, CA, USA.
| | - Genevieve Konopka
- Department of Neuroscience, UT Southwestern Medical Center, Dallas, TX, USA.
- Peter O'Donnell Jr. Brain Institute, UT Southwestern Medical Center, Dallas, TX, USA.
| |
Collapse
|
36
|
Baumann M, Plisson H, Maury S, Renou S, Coqueugniot H, Vanderesse N, Kolobova K, Shnaider S, Rots V, Guérin G, Rendu W. On the Quina side: A Neanderthal bone industry at Chez-Pinaud site, France. PLoS One 2023; 18:e0284081. [PMID: 37315040 DOI: 10.1371/journal.pone.0284081] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/05/2022] [Accepted: 03/23/2023] [Indexed: 06/16/2023] Open
Abstract
Did Neanderthal produce a bone industry? The recent discovery of a large bone tool assemblage at the Neanderthal site of Chagyrskaya (Altai, Siberia, Russia) and the increasing discoveries of isolated finds of bone tools in various Mousterian sites across Eurasia stimulate the debate. Assuming that the isolate finds may be the tip of the iceberg and that the Siberian occurrence did not result from a local adaptation of easternmost Neanderthals, we looked for evidence of a similar industry in the Western side of their spread area. We assessed the bone tool potential of the Quina bone-bed level currently under excavation at chez Pinaud site (Jonzac, Charente-Maritime, France) and found as many bone tools as flint ones: not only the well-known retouchers but also beveled tools, retouched artifacts and a smooth-ended rib. Their diversity opens a window on a range of activities not expected in a butchering site and not documented by the flint tools, all involved in the carcass processing. The re-use of 20% of the bone blanks, which are mainly from large ungulates among faunal remains largely dominated by reindeer, raises the question of blank procurement and management. From the Altai to the Atlantic shore, through a multitude of sites where only a few objects have been reported so far, evidence of a Neanderthal bone industry is emerging which provides new insights on Middle Paleolithic subsistence strategies.
Collapse
Affiliation(s)
| | - Hugues Plisson
- PACEA UMR 5199, CNRS, University of Bordeaux, Pessac, France
| | | | | | - Hélène Coqueugniot
- PACEA UMR 5199, CNRS, University of Bordeaux, Pessac, France
- Ecole Pratique des Hautes Etudes-PSL University, Paris, France
| | | | - Ksenyia Kolobova
- Paleolithic Department, Institute of archaeology and Ethnography, Siberian Branch of the Russian Academy of Sciences, Novosibirsk, Russian Federation
| | - Svetlana Shnaider
- ZooSCAN, International Research Laboratory 2013, CNRS, Novosibirsk, Russian Federation
| | | | - Guillaume Guérin
- Géosciences Rennes UMR 6118, CNRS, University of Rennes, Rennes, France
| | - William Rendu
- ZooSCAN, International Research Laboratory 2013, CNRS, Novosibirsk, Russian Federation
| |
Collapse
|
37
|
Gisladottir RS, Helgason A, Halldorsson BV, Helgason H, Borsky M, Chien YR, Gudnason J, Gudjonsson SA, Moisik S, Dediu D, Thorleifsson G, Tragante V, Bustamante M, Jonsdottir GA, Stefansdottir L, Rutsdottir G, Magnusson SH, Hardarson M, Ferkingstad E, Halldorsson GH, Rognvaldsson S, Skuladottir A, Ivarsdottir EV, Norddahl G, Thorgeirsson G, Jonsdottir I, Ulfarsson MO, Holm H, Stefansson H, Thorsteinsdottir U, Gudbjartsson DF, Sulem P, Stefansson K. Sequence variants affecting voice pitch in humans. SCIENCE ADVANCES 2023; 9:eabq2969. [PMID: 37294764 PMCID: PMC10256171 DOI: 10.1126/sciadv.abq2969] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/30/2022] [Accepted: 05/04/2023] [Indexed: 06/11/2023]
Abstract
The genetic basis of the human vocal system is largely unknown, as are the sequence variants that give rise to individual differences in voice and speech. Here, we couple data on diversity in the sequence of the genome with voice and vowel acoustics in speech recordings from 12,901 Icelanders. We show how voice pitch and vowel acoustics vary across the life span and correlate with anthropometric, physiological, and cognitive traits. We found that voice pitch and vowel acoustics have a heritable component and discovered correlated common variants in ABCC9 that associate with voice pitch. The ABCC9 variants also associate with adrenal gene expression and cardiovascular traits. By showing that voice and vowel acoustics are influenced by genetics, we have taken important steps toward understanding the genetics and evolution of the human vocal system.
Collapse
Affiliation(s)
- Rosa S. Gisladottir
- deCODE Genetics/Amgen Inc., Sturlugata 8, 101 Reykjavik, Iceland
- Department of Icelandic and Comparative Cultural Studies, University of Iceland, Saemundargata 2, 102 Reykjavik, Iceland
| | - Agnar Helgason
- deCODE Genetics/Amgen Inc., Sturlugata 8, 101 Reykjavik, Iceland
- Department of Anthropology, University of Iceland, Saemundargata 10, 102 Reykjavik, Iceland
| | - Bjarni V. Halldorsson
- deCODE Genetics/Amgen Inc., Sturlugata 8, 101 Reykjavik, Iceland
- Department of Engineering, Reykjavik University, Menntavegur 1, 101 Reykjavik, Iceland
| | - Hannes Helgason
- deCODE Genetics/Amgen Inc., Sturlugata 8, 101 Reykjavik, Iceland
| | - Michal Borsky
- Department of Engineering, Reykjavik University, Menntavegur 1, 101 Reykjavik, Iceland
| | - Yu-Ren Chien
- Department of Engineering, Reykjavik University, Menntavegur 1, 101 Reykjavik, Iceland
| | - Jon Gudnason
- Department of Engineering, Reykjavik University, Menntavegur 1, 101 Reykjavik, Iceland
| | | | - Scott Moisik
- Division of Linguistics and Multilingual Studies, Nanyang Technological University, 50 Nanyang Avenue, Singapore 639798, Singapore
| | - Dan Dediu
- Department of Catalan Philology and General Linguistics, University of Barcelona, Gran Via 585, Barcelona 08007, Spain
- University of Barcelona Institute for Complex Systems (UBICS), Martí Franquès 1, Barcelona 08028, Spain
- Catalan Institute for Research and Advanced Studies (ICREA), Passeig Lluís Companys 23, Barcelona 08010, Spain
| | | | | | | | | | | | | | | | | | - Egil Ferkingstad
- deCODE Genetics/Amgen Inc., Sturlugata 8, 101 Reykjavik, Iceland
| | - Gisli H. Halldorsson
- deCODE Genetics/Amgen Inc., Sturlugata 8, 101 Reykjavik, Iceland
- School of Engineering and Natural Sciences, University of Iceland, Dunhagi 5, 107 Reykjavik, Iceland
| | | | | | | | | | - Gudmundur Thorgeirsson
- deCODE Genetics/Amgen Inc., Sturlugata 8, 101 Reykjavik, Iceland
- Faculty of Medicine, University of Iceland, Vatnsmyrarvegur 16, 101 Reykjavik, Iceland
| | - Ingileif Jonsdottir
- deCODE Genetics/Amgen Inc., Sturlugata 8, 101 Reykjavik, Iceland
- Faculty of Medicine, University of Iceland, Vatnsmyrarvegur 16, 101 Reykjavik, Iceland
| | - Magnus O. Ulfarsson
- deCODE Genetics/Amgen Inc., Sturlugata 8, 101 Reykjavik, Iceland
- School of Engineering and Natural Sciences, University of Iceland, Dunhagi 5, 107 Reykjavik, Iceland
| | - Hilma Holm
- deCODE Genetics/Amgen Inc., Sturlugata 8, 101 Reykjavik, Iceland
| | | | - Unnur Thorsteinsdottir
- deCODE Genetics/Amgen Inc., Sturlugata 8, 101 Reykjavik, Iceland
- Faculty of Medicine, University of Iceland, Vatnsmyrarvegur 16, 101 Reykjavik, Iceland
| | - Daniel F. Gudbjartsson
- deCODE Genetics/Amgen Inc., Sturlugata 8, 101 Reykjavik, Iceland
- School of Engineering and Natural Sciences, University of Iceland, Dunhagi 5, 107 Reykjavik, Iceland
| | - Patrick Sulem
- deCODE Genetics/Amgen Inc., Sturlugata 8, 101 Reykjavik, Iceland
| | - Kari Stefansson
- deCODE Genetics/Amgen Inc., Sturlugata 8, 101 Reykjavik, Iceland
- Faculty of Medicine, University of Iceland, Vatnsmyrarvegur 16, 101 Reykjavik, Iceland
| |
Collapse
|
38
|
Kuderna LFK, Gao H, Janiak MC, Kuhlwilm M, Orkin JD, Bataillon T, Manu S, Valenzuela A, Bergman J, Rousselle M, Silva FE, Agueda L, Blanc J, Gut M, de Vries D, Goodhead I, Harris RA, Raveendran M, Jensen A, Chuma IS, Horvath JE, Hvilsom C, Juan D, Frandsen P, Schraiber JG, de Melo FR, Bertuol F, Byrne H, Sampaio I, Farias I, Valsecchi J, Messias M, da Silva MNF, Trivedi M, Rossi R, Hrbek T, Andriaholinirina N, Rabarivola CJ, Zaramody A, Jolly CJ, Phillips-Conroy J, Wilkerson G, Abee C, Simmons JH, Fernandez-Duque E, Kanthaswamy S, Shiferaw F, Wu D, Zhou L, Shao Y, Zhang G, Keyyu JD, Knauf S, Le MD, Lizano E, Merker S, Navarro A, Nadler T, Khor CC, Lee J, Tan P, Lim WK, Kitchener AC, Zinner D, Gut I, Melin AD, Guschanski K, Schierup MH, Beck RMD, Umapathy G, Roos C, Boubli JP, Rogers J, Farh KKH, Marques Bonet T. A global catalog of whole-genome diversity from 233 primate species. Science 2023; 380:906-913. [PMID: 37262161 DOI: 10.1126/science.abn7829] [Citation(s) in RCA: 32] [Impact Index Per Article: 32.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/19/2021] [Accepted: 02/06/2023] [Indexed: 06/03/2023]
Abstract
The rich diversity of morphology and behavior displayed across primate species provides an informative context in which to study the impact of genomic diversity on fundamental biological processes. Analysis of that diversity provides insight into long-standing questions in evolutionary and conservation biology and is urgent given severe threats these species are facing. Here, we present high-coverage whole-genome data from 233 primate species representing 86% of genera and all 16 families. This dataset was used, together with fossil calibration, to create a nuclear DNA phylogeny and to reassess evolutionary divergence times among primate clades. We found within-species genetic diversity across families and geographic regions to be associated with climate and sociality, but not with extinction risk. Furthermore, mutation rates differ across species, potentially influenced by effective population sizes. Lastly, we identified extensive recurrence of missense mutations previously thought to be human specific. This study will open a wide range of research avenues for future primate genomic research.
Collapse
Affiliation(s)
- Lukas F K Kuderna
- IBE, Institute of Evolutionary Biology (UPF-CSIC), Department of Medicine and Life Sciences, Universitat Pompeu Fabra. PRBB, C. Doctor Aiguader N88, 08003 Barcelona, Spain
- Illumina Artificial Intelligence Laboratory, Illumina Inc., Foster City, CA 94404, USA
| | - Hong Gao
- Illumina Artificial Intelligence Laboratory, Illumina Inc., Foster City, CA 94404, USA
| | - Mareike C Janiak
- School of Science, Engineering & Environment, University of Salford, Salford M5 4WT, UK
| | - Martin Kuhlwilm
- IBE, Institute of Evolutionary Biology (UPF-CSIC), Department of Medicine and Life Sciences, Universitat Pompeu Fabra. PRBB, C. Doctor Aiguader N88, 08003 Barcelona, Spain
- Department of Evolutionary Anthropology, University of Vienna, Djerassiplatz 1, 1030 Vienna, Austria
- Human Evolution and Archaeological Sciences (HEAS), University of Vienna, Austria
| | - Joseph D Orkin
- IBE, Institute of Evolutionary Biology (UPF-CSIC), Department of Medicine and Life Sciences, Universitat Pompeu Fabra. PRBB, C. Doctor Aiguader N88, 08003 Barcelona, Spain
- Département d'anthropologie, Université de Montréal, 3150 Jean-Brillant, Montréal, QC H3T 1N8, Canada
| | - Thomas Bataillon
- Bioinformatics Research Centre, Aarhus University, Aarhus, Denmark
| | - Shivakumara Manu
- Academy of Scientific and Innovative Research (AcSIR), Ghaziabad 201002, India
- Laboratory for the Conservation of Endangered Species, CSIR-Centre for Cellular and Molecular Biology, Hyderabad 500007, India
| | - Alejandro Valenzuela
- IBE, Institute of Evolutionary Biology (UPF-CSIC), Department of Medicine and Life Sciences, Universitat Pompeu Fabra. PRBB, C. Doctor Aiguader N88, 08003 Barcelona, Spain
| | - Juraj Bergman
- Bioinformatics Research Centre, Aarhus University, Aarhus, Denmark
- Section for Ecoinformatics and Biodiversity, Department of Biology, Aarhus University, Aarhus, Denmark
| | | | - Felipe Ennes Silva
- Research Group on Primate Biology and Conservation, Mamirauá Institute for Sustainable Development, Estrada da Bexiga 2584, CEP 69553-225, Tefé, Amazonas, Brazil
- Evolutionary Biology and Ecology (EBE), Département de Biologie des Organismes, Université libre de Bruxelles (ULB), Av. Franklin D. Roosevelt 50, CP 160/12, B-1050 Brussels Belgium
| | - Lidia Agueda
- CNAG-CRG, Centre for Genomic Regulation (CRG), Barcelona Institute of Science and Technology (BIST), Baldiri I Reixac 4, 08028 Barcelona, Spain
| | - Julie Blanc
- CNAG-CRG, Centre for Genomic Regulation (CRG), Barcelona Institute of Science and Technology (BIST), Baldiri I Reixac 4, 08028 Barcelona, Spain
| | - Marta Gut
- CNAG-CRG, Centre for Genomic Regulation (CRG), Barcelona Institute of Science and Technology (BIST), Baldiri I Reixac 4, 08028 Barcelona, Spain
| | - Dorien de Vries
- School of Science, Engineering & Environment, University of Salford, Salford M5 4WT, UK
| | - Ian Goodhead
- School of Science, Engineering & Environment, University of Salford, Salford M5 4WT, UK
| | - R Alan Harris
- Human Genome Sequencing Center and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX 77030, USA
| | - Muthuswamy Raveendran
- Human Genome Sequencing Center and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX 77030, USA
| | - Axel Jensen
- Department of Ecology and Genetics, Animal Ecology, Uppsala University, SE-75236 Uppsala, Sweden
| | | | - Julie E Horvath
- North Carolina Museum of Natural Sciences, Raleigh, NC 27601, USA
- Department of Biological and Biomedical Sciences, North Carolina Central University, Durham, NC 27707, USA
- Department of Biological Sciences, North Carolina State University, Raleigh, NC 27695, USA
- Department of Evolutionary Anthropology, Duke University, Durham, NC 27708, USA
- Renaissance Computing Institute, University of North Carolina at Chapel Hill, Chapel Hill, NC 27599, USA
| | | | - David Juan
- IBE, Institute of Evolutionary Biology (UPF-CSIC), Department of Medicine and Life Sciences, Universitat Pompeu Fabra. PRBB, C. Doctor Aiguader N88, 08003 Barcelona, Spain
| | | | - Joshua G Schraiber
- Illumina Artificial Intelligence Laboratory, Illumina Inc., Foster City, CA 94404, USA
| | | | - Fabrício Bertuol
- Universidade Federal do Amazonas, Departamento de Genética, Laboratório de Evolução e Genética Animal (LEGAL), Manaus, Amazonas 69080-900, Brazil
| | - Hazel Byrne
- Department of Anthropology, University of Utah, Salt Lake City. UT 84102, USA
| | | | - Izeni Farias
- Universidade Federal do Amazonas, Departamento de Genética, Laboratório de Evolução e Genética Animal (LEGAL), Manaus, Amazonas 69080-900, Brazil
| | - João Valsecchi
- Research Group on Terrestrial Vertebrate Ecology, Mamirauá Institute for Sustainable Development, Tefé, Amazonas, Brazil
- Rede de Pesquisa para Estudos sobre Diversidade, Conservação e Uso da Fauna na Amazônia - RedeFauna, Manaus, Amazonas, Brazil
- Comunidad de Manejo de Fauna Silvestre en la Amazonía y en Latinoamérica - ComFauna, Iquitos, Loreto, Peru
| | - Malu Messias
- Universidade Federal de Rondônia, Porto Velho, Rondônia, Brazil
| | | | - Mihir Trivedi
- Laboratory for the Conservation of Endangered Species, CSIR-Centre for Cellular and Molecular Biology, Hyderabad 500007, India
| | - Rogerio Rossi
- Instituto de Biociências, Universidade Federal do Mato Grosso, Cuiabá, MT, Brazil
| | - Tomas Hrbek
- Universidade Federal do Amazonas, Departamento de Genética, Laboratório de Evolução e Genética Animal (LEGAL), Manaus, Amazonas 69080-900, Brazil
- Department of Biology, Trinity University, San Antonio, TX 78212, USA
| | - Nicole Andriaholinirina
- Life Sciences and Environment, Technology and Environment of Mahajanga, University of Mahajanga, Mahajanga, Madagascar
| | - Clément J Rabarivola
- Life Sciences and Environment, Technology and Environment of Mahajanga, University of Mahajanga, Mahajanga, Madagascar
| | - Alphonse Zaramody
- Life Sciences and Environment, Technology and Environment of Mahajanga, University of Mahajanga, Mahajanga, Madagascar
| | - Clifford J Jolly
- Department of Anthropology, New York University, New York, NY 10003, USA
| | - Jane Phillips-Conroy
- Department of Neuroscience, Washington University School of Medicine in St. Louis, St. Louis, MO 63110, USA
| | - Gregory Wilkerson
- Keeling Center for Comparative Medicine and Research, MD Anderson Cancer Center, Bastrop TX 78602, USA
| | - Christian Abee
- Keeling Center for Comparative Medicine and Research, MD Anderson Cancer Center, Bastrop TX 78602, USA
| | - Joe H Simmons
- Keeling Center for Comparative Medicine and Research, MD Anderson Cancer Center, Bastrop TX 78602, USA
| | | | - Sree Kanthaswamy
- School of Mathematical and Natural Sciences, Arizona State University, Phoenix, AZ 85004, USA
| | - Fekadu Shiferaw
- Guinea Worm Eradication Program, The Carter Center Ethiopia, Addis Ababa, Ethiopia
| | - Dongdong Wu
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, Yunnan 650223, China
| | - Long Zhou
- Center for Evolutionary and Organismal Biology, Zhejiang University School of Medicine, Hangzhou 310058, China
| | - Yong Shao
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, Yunnan 650223, China
| | - Guojie Zhang
- Center for Evolutionary and Organismal Biology, Zhejiang University School of Medicine, Hangzhou 310058, China
- Villum Centre for Biodiversity Genomics, Section for Ecology and Evolution, Department of Biology, University of Copenhagen, DK-2100 Copenhagen, Denmark
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, Yunnan 650223, China
- Liangzhu Laboratory, Zhejiang University Medical Center, 1369 West Wenyi Road, Hangzhou 311121, China
- Women's Hospital, School of Medicine, Zhejiang University, 1 Xueshi Road, Shangcheng District, Hangzhou 310006, China
| | - Julius D Keyyu
- Tanzania Wildlife Research Institute (TAWIRI), Head Office, P.O. Box 661, Arusha, Tanzania
| | - Sascha Knauf
- Institute of International Animal Health/One Health, Friedrich-Loeffler-Institut, Federal Research Institute for Animal Health, 17493 Greifswald-Insel Riems, Germany
| | - Minh D Le
- Department of Environmental Ecology, Faculty of Environmental Sciences, University of Science and Central Institute for Natural Resources and Environmental Studies, Vietnam National University, Hanoi, Vietnam
| | - Esther Lizano
- IBE, Institute of Evolutionary Biology (UPF-CSIC), Department of Medicine and Life Sciences, Universitat Pompeu Fabra. PRBB, C. Doctor Aiguader N88, 08003 Barcelona, Spain
- Institut Català de Paleontologia Miquel Crusafont, Universitat Autònoma de Barcelona, Barcelona, Spain
| | - Stefan Merker
- Department of Zoology, State Museum of Natural History Stuttgart, Stuttgart, Germany
| | - Arcadi Navarro
- IBE, Institute of Evolutionary Biology (UPF-CSIC), Department of Medicine and Life Sciences, Universitat Pompeu Fabra. PRBB, C. Doctor Aiguader N88, 08003 Barcelona, Spain
- Institució Catalana de Recerca i Estudis Avançats (ICREA) and Universitat Pompeu Fabra. Pg. Luís Companys 23, 08010 Barcelona, Spain
- Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology, Av. Doctor Aiguader, N88, 08003 Barcelona, Spain
- BarcelonaBeta Brain Research Center, Pasqual Maragall Foundation, C. Wellington 30, 08005 Barcelona, Spain
| | - Tilo Nadler
- Cuc Phuong Commune, Nho Quan District, Ninh Binh Province, Vietnam
| | - Chiea Chuen Khor
- Genome Institute of Singapore, Agency for Science, Technology and Research, Singapore
| | - Jessica Lee
- Mandai Nature, 80 Mandai Lake Road, Singapore
| | - Patrick Tan
- Genome Institute of Singapore, Agency for Science, Technology and Research, Singapore
- SingHealth Duke-NUS Institute of Precision Medicine (PRISM), Singapore
- Cancer and Stem Cell Biology Program, Duke-NUS Medical School, Singapore
| | - Weng Khong Lim
- SingHealth Duke-NUS Institute of Precision Medicine (PRISM), Singapore
- Cancer and Stem Cell Biology Program, Duke-NUS Medical School, Singapore
- SingHealth Duke-NUS Genomic Medicine Centre, Singapore
| | - Andrew C Kitchener
- Department of Natural Sciences, National Museums Scotland, Chambers Street, Edinburgh EH1 1JF, UK, and School of Geosciences, Drummond Street, Edinburgh EH8 9XP, UK
| | - Dietmar Zinner
- Cognitive Ethology Laboratory, Germany Primate Center, Leibniz Institute for Primate Research, 37077 Göttingen, Germany
- Department of Primate Cognition, Georg-August-Universität Göttingen, 37077 Göttingen, Germany
- Leibniz ScienceCampus Primate Cognition, 37077 Göttingen, Germany
| | - Ivo Gut
- CNAG-CRG, Centre for Genomic Regulation (CRG), Barcelona Institute of Science and Technology (BIST), Baldiri I Reixac 4, 08028 Barcelona, Spain
| | - Amanda D Melin
- Department of Anthropology and Archaeology, University of Calgary, 2500 University Dr NW, Calgary, AB T2N 1N4, Canada
- Department of Medical Genetics, University of Calgary, 3330 Hospital Drive NW, HMRB 202, Calgary, AB T2N 4N1, Canada
- Alberta Children's Hospital Research Institute, University of Calgary, 3330 Hospital Drive NW, HMRB 202, Calgary, AB T2N 4N1, Canada
| | - Katerina Guschanski
- Department of Ecology and Genetics, Animal Ecology, Uppsala University, SE-75236 Uppsala, Sweden
- Institute of Ecology and Evolution, School of Biological Sciences, University of Edinburgh, Edinburgh, UK
| | | | - Robin M D Beck
- School of Science, Engineering & Environment, University of Salford, Salford M5 4WT, UK
| | - Govindhaswamy Umapathy
- Academy of Scientific and Innovative Research (AcSIR), Ghaziabad 201002, India
- Laboratory for the Conservation of Endangered Species, CSIR-Centre for Cellular and Molecular Biology, Hyderabad 500007, India
| | - Christian Roos
- Gene Bank of Primates and Primate Genetics Laboratory, German Primate Center, Leibniz Institute for Primate Research, Kellnerweg 4, 37077 Göttingen, Germany
| | - Jean P Boubli
- School of Science, Engineering & Environment, University of Salford, Salford M5 4WT, UK
| | - Jeffrey Rogers
- Human Genome Sequencing Center and Department of Molecular and Human Genetics, Baylor College of Medicine, Houston, TX 77030, USA
| | - Kyle Kai-How Farh
- Illumina Artificial Intelligence Laboratory, Illumina Inc., Foster City, CA 94404, USA
| | - Tomas Marques Bonet
- IBE, Institute of Evolutionary Biology (UPF-CSIC), Department of Medicine and Life Sciences, Universitat Pompeu Fabra. PRBB, C. Doctor Aiguader N88, 08003 Barcelona, Spain
- CNAG-CRG, Centre for Genomic Regulation (CRG), Barcelona Institute of Science and Technology (BIST), Baldiri I Reixac 4, 08028 Barcelona, Spain
- Institut Català de Paleontologia Miquel Crusafont, Universitat Autònoma de Barcelona, Barcelona, Spain
- Institució Catalana de Recerca i Estudis Avançats (ICREA) and Universitat Pompeu Fabra. Pg. Luís Companys 23, 08010 Barcelona, Spain
| |
Collapse
|
39
|
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
|
40
|
Ågren R, Patil S, Zhou X, Sahlholm K, Pääbo S, Zeberg H. Major Genetic Risk Factors for Dupuytren's Disease Are Inherited From Neandertals. Mol Biol Evol 2023; 40:msad130. [PMID: 37315093 DOI: 10.1093/molbev/msad130] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/16/2023] Open
Abstract
Dupuytren's disease is characterized by fingers becoming permanently bent in a flexed position. Whereas people of African ancestry are rarely afflicted by Dupuytren's disease, up to ∼30% of men over 60 years suffer from this condition in northern Europe. Here, we meta-analyze 3 biobanks comprising 7,871 cases and 645,880 controls and find 61 genome-wide significant variants associated with Dupuytren's disease. We show that 3 of the 61 loci harbor alleles of Neandertal origin, including the second and third most strongly associated ones (P = 6.4 × 10-132 and P = 9.2 × 10-69, respectively). For the most strongly associated Neandertal variant, we identify EPDR1 as the causal gene. Dupuytren's disease is an example of how admixture with Neandertals has shaped regional differences in disease prevalence.
Collapse
Affiliation(s)
- Richard Ågren
- Department of Physiology and Pharmacology, Karolinska Institutet, Stockholm, Sweden
| | - Snehal Patil
- Department of Biostatistics, University of Michigan, Ann Arbor, MI, USA
| | - Xiang Zhou
- Department of Biostatistics, University of Michigan, Ann Arbor, MI, USA
| | - Kristoffer Sahlholm
- Department of Physiology and Pharmacology, Karolinska Institutet, Stockholm, Sweden
- Department of Integrative Medical Biology, Wallenberg Centre for Molecular Medicine, Umeå University, Umeå, Sweden
| | - Svante Pääbo
- Department of Evolutionary Genetics, Max Planck Institute for Evolutionary Anthropology, Leipzig, Germany
- Human Evolutionary Genomics Unit, Okinawa Institute of Science and Technology, Okinawa, Japan
| | - Hugo Zeberg
- Department of Physiology and Pharmacology, Karolinska Institutet, Stockholm, Sweden
- Department of Evolutionary Genetics, Max Planck Institute for Evolutionary Anthropology, Leipzig, Germany
- Human Evolutionary Genomics Unit, Okinawa Institute of Science and Technology, Okinawa, Japan
| |
Collapse
|
41
|
Petersen J, Englmaier L, Artemov AV, Poverennaya I, Mahmoud R, Bouderlique T, Tesarova M, Deviatiiarov R, Szilvásy-Szabó A, Akkuratov EE, Pajuelo Reguera D, Zeberg H, Kaucka M, Kastriti ME, Krivanek J, Radaszkiewicz T, Gömöryová K, Knauth S, Potesil D, Zdrahal Z, Ganji RS, Grabowski A, Buhl ME, Zikmund T, Kavkova M, Axelson H, Lindgren D, Kramann R, Kuppe C, Erdélyi F, Máté Z, Szabó G, Koehne T, Harkany T, Fried K, Kaiser J, Boor P, Fekete C, Rozman J, Kasparek P, Prochazka J, Sedlacek R, Bryja V, Gusev O, Adameyko I. A previously uncharacterized Factor Associated with Metabolism and Energy (FAME/C14orf105/CCDC198/1700011H14Rik) is related to evolutionary adaptation, energy balance, and kidney physiology. Nat Commun 2023; 14:3092. [PMID: 37248239 DOI: 10.1038/s41467-023-38663-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/24/2021] [Accepted: 05/11/2023] [Indexed: 05/31/2023] Open
Abstract
In this study we use comparative genomics to uncover a gene with uncharacterized function (1700011H14Rik/C14orf105/CCDC198), which we hereby name FAME (Factor Associated with Metabolism and Energy). We observe that FAME shows an unusually high evolutionary divergence in birds and mammals. Through the comparison of single nucleotide polymorphisms, we identify gene flow of FAME from Neandertals into modern humans. We conduct knockout experiments on animals and observe altered body weight and decreased energy expenditure in Fame knockout animals, corresponding to genome-wide association studies linking FAME with higher body mass index in humans. Gene expression and subcellular localization analyses reveal that FAME is a membrane-bound protein enriched in the kidneys. Although the gene knockout results in structurally normal kidneys, we detect higher albumin in urine and lowered ferritin in the blood. Through experimental validation, we confirm interactions between FAME and ferritin and show co-localization in vesicular and plasma membranes.
Collapse
Affiliation(s)
- Julian Petersen
- Department of Orthodontics, University Leipzig Medical Center, Leipzig, Germany.
| | - Lukas Englmaier
- CeMM Research Center for Molecular Medicine of the Austrian Academy of Sciences, 1090, Vienna, Austria
- Ludwig Boltzmann Institute for Rare and Undiagnosed Diseases, 1090, Vienna, Austria
| | - Artem V Artemov
- Department of Neuroimmunology, Center for Brain Research, Medical University Vienna, Vienna, Austria
| | - Irina Poverennaya
- Department of Neuroimmunology, Center for Brain Research, Medical University Vienna, Vienna, Austria
| | - Ruba Mahmoud
- Department of Orthodontics, University Leipzig Medical Center, Leipzig, Germany
| | - Thibault Bouderlique
- Department of Neuroimmunology, Center for Brain Research, Medical University Vienna, Vienna, Austria
| | - Marketa Tesarova
- Central European Institute of Technology, Brno University of Technology, Brno, Czech Republic
| | - Ruslan Deviatiiarov
- Regulatory Genomics Research Center, Institute of Fundamental Medicine and Biology, Kazan Federal University, Kazan, Russia
- Endocrinology Research Center, Moscow, Russia
| | - Anett Szilvásy-Szabó
- Laboratory of Integrative Neuroendocrinology, Institute of Experimental Medicine, 1083, Budapest, Hungary
| | - Evgeny E Akkuratov
- Department of Applied Physics, Royal Institute of Technology, Science for Life Laboratory, 171 65, Stockholm, Sweden
- University of Oxford, MRC Weatherall Institute of Molecular Medicine, Radcliffe Department of Medicine, Oxford, OX3 9DS, UK
| | - David Pajuelo Reguera
- Institute of Molecular Genetics of the Czech Academy of Science, Czech Centre for Phenogenomics, Vestec, Czech Republic
| | - Hugo Zeberg
- Department of Neuroscience, Karolinska Institutet, Stockholm, Sweden
- Department of Physiology and Pharmacology, Karolinska Institutet, Stockholm, Sweden
| | - Marketa Kaucka
- Max Planck Institute for Evolutionary Biology, Plön, 24306, Germany
| | - Maria Eleni Kastriti
- Department of Neuroimmunology, Center for Brain Research, Medical University Vienna, Vienna, Austria
- Department of Physiology and Pharmacology, Karolinska Institutet, Stockholm, Sweden
| | - Jan Krivanek
- Department of Histology and Embryology, Faculty of Medicine, Masaryk University, Brno, Czech Republic
| | - Tomasz Radaszkiewicz
- Institute of Experimental Biology, Faculty of Science, Masaryk University, Brno, Czech Republic
| | - Kristína Gömöryová
- Institute of Experimental Biology, Faculty of Science, Masaryk University, Brno, Czech Republic
| | - Sarah Knauth
- Department of Orthodontics, University Leipzig Medical Center, Leipzig, Germany
| | - David Potesil
- Central European Institute of Technology, Masaryk University, Brno, Czech Republic
| | - Zbynek Zdrahal
- Central European Institute of Technology, Masaryk University, Brno, Czech Republic
| | - Ranjani Sri Ganji
- Central European Institute of Technology, Masaryk University, Brno, Czech Republic
| | - Anna Grabowski
- Department of Neuroimmunology, Center for Brain Research, Medical University Vienna, Vienna, Austria
| | - Miriam E Buhl
- Institute of Pathology & Electron Microscopy Facility, RWTH Aachen University Hospital, Aachen, Germany
| | - Tomas Zikmund
- Central European Institute of Technology, Brno University of Technology, Brno, Czech Republic
| | - Michaela Kavkova
- Central European Institute of Technology, Brno University of Technology, Brno, Czech Republic
- Department of Histology and Embryology, Faculty of Medicine, Masaryk University, Brno, Czech Republic
| | - Håkan Axelson
- Translational Cancer Research, Department of Laboratory Medicine, Lund University, Medicon Village, Scheelevägen 2, Lund, Sweden
| | - David Lindgren
- Translational Cancer Research, Department of Laboratory Medicine, Lund University, Medicon Village, Scheelevägen 2, Lund, Sweden
| | - Rafael Kramann
- Institute of Experimental Medicine and Systems Biology, RWTH Aachen University, Aachen, Germany
| | - Christoph Kuppe
- Institute of Experimental Medicine and Systems Biology, RWTH Aachen University, Aachen, Germany
| | - Ferenc Erdélyi
- Medical Gene Technology Unit, Institute of Experimental Medicine, Budapest, Hungary
| | - Zoltán Máté
- Medical Gene Technology Unit, Institute of Experimental Medicine, Budapest, Hungary
| | - Gábor Szabó
- Medical Gene Technology Unit, Institute of Experimental Medicine, Budapest, Hungary
| | - Till Koehne
- Department of Orthodontics, University Leipzig Medical Center, Leipzig, Germany
| | - Tibor Harkany
- Department of Molecular Neurosciences, Center for Brain Research, Medical University Vienna, Vienna, Austria
| | - Kaj Fried
- Department of Physiology and Pharmacology, Karolinska Institutet, Stockholm, Sweden
| | - Jozef Kaiser
- Central European Institute of Technology, Brno University of Technology, Brno, Czech Republic
| | - Peter Boor
- Institute of Pathology & Electron Microscopy Facility, RWTH Aachen University Hospital, Aachen, Germany
| | - Csaba Fekete
- Laboratory of Integrative Neuroendocrinology, Institute of Experimental Medicine, 1083, Budapest, Hungary
| | - Jan Rozman
- Institute of Molecular Genetics of the Czech Academy of Science, Czech Centre for Phenogenomics, Vestec, Czech Republic
- Luxembourg Centre for Systems Biomedicine, University of Luxembourg, 6, avenue du Swing, 4367, Belvaux, Luxembourg
| | - Petr Kasparek
- Institute of Molecular Genetics of the Czech Academy of Science, Czech Centre for Phenogenomics, Vestec, Czech Republic
| | - Jan Prochazka
- Institute of Molecular Genetics of the Czech Academy of Science, Czech Centre for Phenogenomics, Vestec, Czech Republic
| | - Radislav Sedlacek
- Institute of Molecular Genetics of the Czech Academy of Science, Czech Centre for Phenogenomics, Vestec, Czech Republic
| | - Vitezslav Bryja
- Institute of Experimental Biology, Faculty of Science, Masaryk University, Brno, Czech Republic
| | - Oleg Gusev
- Endocrinology Research Center, Moscow, Russia
- Intractable Disease Research Center, Graduate School of Medicine, Juntendo University, Tokyo, Japan
- Center for Integrative Medical Sciences, RIKEN, Yokohama City, Kanagawa, Japan
| | - Igor Adameyko
- Department of Neuroimmunology, Center for Brain Research, Medical University Vienna, Vienna, Austria.
- Department of Physiology and Pharmacology, Karolinska Institutet, Stockholm, Sweden.
| |
Collapse
|
42
|
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
|
43
|
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
|
44
|
Jackson A, Lin SJ, Jones EA, Chandler KE, Orr D, Moss C, Haider Z, Ryan G, Holden S, Harrison M, Burrows N, Jones WD, Loveless M, Petree C, Stewart H, Low K, Donnelly D, Lovell S, Drosou K, Varshney GK, Banka S. Clinical, genetic, epidemiologic, evolutionary, and functional delineation of TSPEAR-related autosomal recessive ectodermal dysplasia 14. HGG ADVANCES 2023; 4:100186. [PMID: 37009414 PMCID: PMC10064225 DOI: 10.1016/j.xhgg.2023.100186] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/04/2022] [Accepted: 02/27/2023] [Indexed: 06/11/2023] Open
Abstract
TSPEAR variants cause autosomal recessive ectodermal dysplasia (ARED) 14. The function of TSPEAR is unknown. The clinical features, the mutation spectrum, and the underlying mechanisms of ARED14 are poorly understood. Combining data from new and previously published individuals established that ARED14 is primarily characterized by dental anomalies such as conical tooth cusps and hypodontia, like those seen in individuals with WNT10A-related odontoonychodermal dysplasia. AlphaFold-predicted structure-based analysis showed that most of the pathogenic TSPEAR missense variants likely destabilize the β-propeller of the protein. Analysis of 100000 Genomes Project (100KGP) data revealed multiple founder TSPEAR variants across different populations. Mutational and recombination clock analyses demonstrated that non-Finnish European founder variants likely originated around the end of the last ice age, a period of major climatic transition. Analysis of gnomAD data showed that the non-Finnish European population TSPEAR gene-carrier rate is ∼1/140, making it one of the commonest AREDs. Phylogenetic and AlphaFold structural analyses showed that TSPEAR is an ortholog of drosophila Closca, an extracellular matrix-dependent signaling regulator. We, therefore, hypothesized that TSPEAR could have a role in enamel knot, a structure that coordinates patterning of developing tooth cusps. Analysis of mouse single-cell RNA sequencing (scRNA-seq) data revealed highly restricted expression of Tspear in clusters representing enamel knots. A tspeara -/-;tspearb -/- double-knockout zebrafish model recapitulated the clinical features of ARED14 and fin regeneration abnormalities of wnt10a knockout fish, thus suggesting interaction between tspear and wnt10a. In summary, we provide insights into the role of TSPEAR in ectodermal development and the evolutionary history, epidemiology, mechanisms, and consequences of its loss of function variants.
Collapse
Affiliation(s)
- Adam Jackson
- Division of Evolution, Infection and Genomics, School of Biological Sciences, Faculty of Biology, Medicine and Health, University of Manchester, Manchester, UK
- Manchester Centre for Genomic Medicine, St Mary’s Hospital, Manchester University NHS Foundation Trust, Health Innovation Manchester, Manchester, UK
| | - Sheng-Jia Lin
- Genes and Human Disease Research Program, Oklahoma Medical Research Foundation, Oklahoma City, OK, USA
| | - Elizabeth A. Jones
- Division of Evolution, Infection and Genomics, School of Biological Sciences, Faculty of Biology, Medicine and Health, University of Manchester, Manchester, UK
- Manchester Centre for Genomic Medicine, St Mary’s Hospital, Manchester University NHS Foundation Trust, Health Innovation Manchester, Manchester, UK
| | - Kate E. Chandler
- Manchester Centre for Genomic Medicine, St Mary’s Hospital, Manchester University NHS Foundation Trust, Health Innovation Manchester, Manchester, UK
| | - David Orr
- Manchester Centre for Genomic Medicine, St Mary’s Hospital, Manchester University NHS Foundation Trust, Health Innovation Manchester, Manchester, UK
| | - Celia Moss
- Department of Dermatology, Birmingham Children’s Hospital, Birmingham Women’s and Children’s NHS Foundation Trust, Birmingham, UK
| | - Zahra Haider
- Department of Dermatology, Birmingham Children’s Hospital, Birmingham Women’s and Children’s NHS Foundation Trust, Birmingham, UK
| | - Gavin Ryan
- West Midlands Regional Genetics Laboratory, Birmingham Women’s and Children’s NHS Foundation Trust, Birmingham, UK
| | - Simon Holden
- Clinical Genetics, Addenbrooke’s Hospital, Cambridge, UK
| | - Mike Harrison
- Department of Pediatric Dentistry, Guy’s and St Thomas' Dental Institute, London, UK
| | - Nigel Burrows
- Department of Dermatology, Cambridge University Hospitals NHS Foundation Trust, Cambridge, UK
| | - Wendy D. Jones
- North East Thames Regional Genetics Service, Great Ormond Street Hospital for Children, Great Ormond Street NHS Foundation Trust, London, UK
| | - Mary Loveless
- Genes and Human Disease Research Program, Oklahoma Medical Research Foundation, Oklahoma City, OK, USA
| | - Cassidy Petree
- Genes and Human Disease Research Program, Oklahoma Medical Research Foundation, Oklahoma City, OK, USA
| | - Helen Stewart
- Oxford Centre for Genomic Medicine, Oxford University Hospitals NHS Foundation Trust, Oxford, UK
| | - Karen Low
- Department of Clinical Genetics, St Michael’s Hospital, Bristol, UK
| | - Deirdre Donnelly
- Department of Genetic Medicine, Belfast HSC Trust, Lisburn Road, Belfast, UK
| | - Simon Lovell
- Division of Evolution, Infection and Genomics, School of Biological Sciences, Faculty of Biology, Medicine and Health, University of Manchester, Manchester, UK
| | - Konstantina Drosou
- Department of Earth and Environmental Sciences, Manchester Institute of Biotechnology, University of Manchester, Manchester, UK
- Division of Cell Matrix Biology and Regenerative Medicine, Faculty of Biology, Medicine and Health, University of Manchester, 99 Oxford Road, Manchester, UK
| | - Gaurav K. Varshney
- Genes and Human Disease Research Program, Oklahoma Medical Research Foundation, Oklahoma City, OK, USA
| | - Siddharth Banka
- Division of Evolution, Infection and Genomics, School of Biological Sciences, Faculty of Biology, Medicine and Health, University of Manchester, Manchester, UK
- Manchester Centre for Genomic Medicine, St Mary’s Hospital, Manchester University NHS Foundation Trust, Health Innovation Manchester, Manchester, UK
| |
Collapse
|
45
|
Pathogenic Variants Associated with Rare Monogenic Diseases Established in Ancient Neanderthal and Denisovan Genome-Wide Data. Genes (Basel) 2023; 14:genes14030727. [PMID: 36980999 PMCID: PMC10048696 DOI: 10.3390/genes14030727] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2023] [Revised: 03/08/2023] [Accepted: 03/13/2023] [Indexed: 03/18/2023] Open
Abstract
Ancient anatomically modern humans (AMHs) encountered other archaic human species, most notably Neanderthals and Denisovans, when they left Africa and spread across Europe and Asia ~60,000 years ago. They interbred with them, and modern human genomes retain DNA inherited from these interbreeding events. High quality (high coverage) ancient human genomes have recently been sequenced allowing for a direct estimation of individual heterozygosity, which has shown that genetic diversity in these archaic human groups was very low, indicating low population sizes. In this study, we analyze ten ancient human genome-wide data, including four sequenced with high-coverage. We screened these ancient genome-wide data for pathogenic mutations associated with monogenic diseases, and established unusual aggregation of pathogenic mutations in individual subjects, including quadruple homozygous cases of pathogenic variants in the PAH gene associated with the condition phenylketonuria in a ~120,000 years old Neanderthal. Such aggregation of pathogenic mutations is extremely rare in contemporary populations, and their existence in ancient humans could be explained by less significant clinical manifestations coupled with small community sizes, leading to higher inbreeding levels. Our results suggest that pathogenic variants associated with rare diseases might be the result of introgression from other archaic human species, and archaic admixture thus could have influenced disease risk in modern humans.
Collapse
|
46
|
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
|
47
|
Wang PY, Yang Y, Shi XQ, Chen Y, Liu SD, Wang HY, Peng T, Shi Q, Zhang W, Sun C. Distilling functional variations for human UGT2B4 upstream region based on selection signals and implications for phenotypes of Neanderthal and Denisovan. Sci Rep 2023; 13:3134. [PMID: 36823244 PMCID: PMC9950360 DOI: 10.1038/s41598-023-29682-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2022] [Accepted: 02/08/2023] [Indexed: 02/25/2023] Open
Abstract
Our previous work identified one region upstream human UGT2B4 (UDP glucuronosyltransferase family 2 member B4) which is associated with breast cancer and under balancing selection. However, the distribution, functional variation and molecular mechanism underlying breast cancer and balancing selection remain unclear. In current study, the two haplotypes with deep divergence are described by analyzing 1000 genomes project data and observed to be with high frequencies in all human populations. Through population genetics analysis and genome annotation, the potential functional region is identified and verified by reporter gene assay. Further mutagenesis indicates that the functional mutations are rs66862535 and rs68096061. Both SNPs can alter the interaction efficiency of transcription factor POU2F1 (POU class 2 homeobox 1). Through chromosome conformation capture, it is identified that the enhancer containing these two SNPs can interact with UGT2B4 promoter. Expression quantitative trait loci analysis indicates that UGT2B4 expression is dependent on the genotype of this locus. The common haplotype in human is lost in four genomes of archaic hominins, which suggests that Neanderthal and Denisovan should present relatively lower UGT2B4 expression and further higher steroid hormone level. This study provides new insight into the contribution of ancient population structure to human phenotypes.
Collapse
Affiliation(s)
- Pin-Yi Wang
- grid.412498.20000 0004 1759 8395College of Life Sciences, Shaanxi Normal University, Xi’an, 710119 Shaanxi People’s Republic of China ,grid.440773.30000 0000 9342 2456State Key Laboratory for Conservation and Utilization of Bio-Resources in Yunnan, Yunnan University, Kunming, 650091 Yunnan People’s Republic of China
| | - Yuan Yang
- grid.412498.20000 0004 1759 8395College of Life Sciences, Shaanxi Normal University, Xi’an, 710119 Shaanxi People’s Republic of China
| | - Xiao-Qian Shi
- grid.412498.20000 0004 1759 8395College of Life Sciences, Shaanxi Normal University, Xi’an, 710119 Shaanxi People’s Republic of China
| | - Ying Chen
- grid.412498.20000 0004 1759 8395College of Life Sciences, Shaanxi Normal University, Xi’an, 710119 Shaanxi People’s Republic of China
| | - Shao-Dong Liu
- grid.412498.20000 0004 1759 8395College of Life Sciences, Shaanxi Normal University, Xi’an, 710119 Shaanxi People’s Republic of China
| | - Hong-Yan Wang
- grid.412498.20000 0004 1759 8395College of Life Sciences, Shaanxi Normal University, Xi’an, 710119 Shaanxi People’s Republic of China
| | - Tao Peng
- grid.440773.30000 0000 9342 2456State Key Laboratory for Conservation and Utilization of Bio-Resources in Yunnan, Yunnan University, Kunming, 650091 Yunnan People’s Republic of China
| | - Qiang Shi
- grid.412498.20000 0004 1759 8395College of Life Sciences, Shaanxi Normal University, Xi’an, 710119 Shaanxi People’s Republic of China
| | - Wei Zhang
- grid.16753.360000 0001 2299 3507Department of Preventive Medicine, Northwestern University Feinberg School of Medicine, Chicago, IL 60611 USA ,grid.449428.70000 0004 1797 7280Institute of Precision Medicine, Jining Medical University, Jining, 272067 Shandong People’s Republic of China
| | - Chang Sun
- College of Life Sciences, Shaanxi Normal University, Xi'an, 710119, Shaanxi, People's Republic of China.
| |
Collapse
|
48
|
Boeckx C. What made us "hunter-gatherers of words". Front Neurosci 2023; 17:1080861. [PMID: 36845441 PMCID: PMC9947416 DOI: 10.3389/fnins.2023.1080861] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/26/2022] [Accepted: 01/19/2023] [Indexed: 02/11/2023] Open
Abstract
This paper makes three interconnected claims: (i) the "human condition" cannot be captured by evolutionary narratives that reduce it to a recent 'cognitive modernity', nor by narratives that eliminates all cognitive differences between us and out closest extinct relatives, (ii) signals from paleogenomics, especially coming from deserts of introgression but also from signatures of positive selection, point to the importance of mutations that impact neurodevelopment, plausibly leading to temperamental differences, which may impact cultural evolutionary trajectories in specific ways, and (iii) these trajectories are expected to affect the language phenotypes, modifying what is being learned and how it is put to use. In particular, I hypothesize that these different trajectories influence the development of symbolic systems, the flexible ways in which symbols combine, and the size and configurations of the communities in which these systems are put to use.
Collapse
Affiliation(s)
- Cedric Boeckx
- Section of General Linguistics, Universitat de Barcelona, Barcelona, Spain
- Institute of Complex Systems, Universitat de Barcelona, Barcelona, Spain
- Catalan Institute for Research and Advanced Studies (ICREA), Barcelona, Spain
| |
Collapse
|
49
|
Flegontov P, Işıldak U, Maier R, Yüncü E, Changmai P, Reich D. Modeling of African population history using f -statistics can be highly biased and is not addressed by previously suggested SNP ascertainment schemes. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.01.22.525077. [PMID: 36711923 PMCID: PMC9882349 DOI: 10.1101/2023.01.22.525077] [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/24/2023]
Abstract
f -statistics have emerged as a first line of analysis for making inferences about demographic history from genome-wide data. These statistics can provide strong evidence for either admixture or cladality, which can be robust to substantial rates of errors or missing data. f -statistics are guaranteed to be unbiased under "SNP ascertainment" (analyzing non-randomly chosen subsets of single nucleotide polymorphisms) only if it relies on a population that is an outgroup for all groups analyzed. However, ascertainment on a true outgroup that is not co-analyzed with other populations is often impractical and uncommon in the literature. In this study focused on practical rather than theoretical aspects of SNP ascertainment, we show that many non-outgroup ascertainment schemes lead to false rejection of true demographic histories, as well as to failure to reject incorrect models. But the bias introduced by common ascertainments such as the 1240K panel is mostly limited to situations when more than one sub-Saharan African and/or archaic human groups (Neanderthals and Denisovans) or non-human outgroups are co-modelled, for example, f 4 -statistics involving one non-African group, two African groups, and one archaic group. Analyzing panels of SNPs polymorphic in archaic humans, which has been suggested as a solution for the ascertainment problem, cannot fix all these problems since for some classes of f -statistics it is not a clean outgroup ascertainment, and in other cases it demonstrates relatively low power to reject incorrect demographic models since it provides a relatively small number of variants common in anatomically modern humans. And due to the paucity of high-coverage archaic genomes, archaic individuals used for ascertainment often act as sole representatives of the respective groups in an analysis, and we show that this approach is highly problematic. By carrying out large numbers of simulations of diverse demographic histories, we find that bias in inferences based on f -statistics introduced by non-outgroup ascertainment can be minimized if the derived allele frequency spectrum in the population used for ascertainment approaches the spectrum that existed at the root of all groups being co-analyzed. Ascertaining on sites with variants common in a diverse group of African individuals provides a good approximation to such a set of SNPs, addressing the great majority of biases and also retaining high statistical power for studying population history. Such a "pan-African" ascertainment, although not completely problem-free, allows unbiased exploration of demographic models for the widest set of archaic and modern human populations, as compared to the other ascertainment schemes we explored.
Collapse
Affiliation(s)
- Pavel Flegontov
- Department of Human Evolutionary Biology, Harvard University, Cambridge, MA, USA
- Department of Biology and Ecology, Faculty of Science, University of Ostrava, Ostrava, Czechia
- Kalmyk Research Center of the Russian Academy of Sciences, Elista, Russia
| | - Ulaş Işıldak
- Department of Biology and Ecology, Faculty of Science, University of Ostrava, Ostrava, Czechia
| | - Robert Maier
- Department of Human Evolutionary Biology, Harvard University, Cambridge, MA, USA
| | - Eren Yüncü
- Department of Biology and Ecology, Faculty of Science, University of Ostrava, Ostrava, Czechia
| | - Piya Changmai
- Department of Biology and Ecology, Faculty of Science, University of Ostrava, Ostrava, Czechia
| | - David Reich
- Department of Human Evolutionary Biology, Harvard University, Cambridge, MA, USA
- Department of Genetics, Harvard Medical School, Boston, MA 02115, USA
- Howard Hughes Medical Institute, Harvard Medical School, Boston, MA, USA
- Broad Institute of Harvard and MIT, Cambridge, MA, USA
| |
Collapse
|
50
|
Popli D, Peyrégne S, Peter BM. KIN: a method to infer relatedness from low-coverage ancient DNA. Genome Biol 2023; 24:10. [PMID: 36650598 PMCID: PMC9843908 DOI: 10.1186/s13059-023-02847-7] [Citation(s) in RCA: 15] [Impact Index Per Article: 15.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2022] [Accepted: 01/04/2023] [Indexed: 01/19/2023] Open
Abstract
Genetic kinship of ancient individuals can provide insights into their culture and social hierarchy, and is relevant for downstream genetic analyses. However, estimating relatedness from ancient DNA is difficult due to low-coverage, ascertainment bias, or contamination from various sources. Here, we present KIN, a method to estimate the relatedness of a pair of individuals from the identical-by-descent segments they share. KIN accurately classifies up to 3rd-degree relatives using at least 0.05x sequence coverage and differentiates siblings from parent-child pairs. It incorporates additional models to adjust for contamination and detect inbreeding, which improves classification accuracy.
Collapse
Affiliation(s)
- Divyaratan Popli
- grid.419518.00000 0001 2159 1813 Department of Evolutionary Genetics, Max Planck Institute for Evolutionary Anthropology, Leipzig, Germany
| | - Stéphane Peyrégne
- grid.419518.00000 0001 2159 1813 Department of Evolutionary Genetics, Max Planck Institute for Evolutionary Anthropology, Leipzig, Germany
| | - Benjamin M. Peter
- grid.419518.00000 0001 2159 1813 Department of Evolutionary Genetics, Max Planck Institute for Evolutionary Anthropology, Leipzig, Germany
| |
Collapse
|