1
|
Zheng HX, Yan S, Zhang M, Gu Z, Wang J, Jin L. Mitochondrial DNA Genomes Reveal Relaxed Purifying Selection During Human Population Expansion after the Last Glacial Maximum. Mol Biol Evol 2024; 41:msae175. [PMID: 39162340 PMCID: PMC11373649 DOI: 10.1093/molbev/msae175] [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/27/2023] [Revised: 07/31/2024] [Accepted: 08/06/2024] [Indexed: 08/21/2024] Open
Abstract
Modern humans have experienced explosive population growth in the past thousand years. We hypothesized that recent human populations have inhabited environments with relaxation of selective constraints, possibly due to the more abundant food supply after the Last Glacial Maximum. The ratio of nonsynonymous to synonymous mutations (N/S ratio) is a useful and common statistic for measuring selective constraints. In this study, we reconstructed a high-resolution phylogenetic tree using a total of 26,419 East Eurasian mitochondrial DNA genomes, which were further classified into expansion and nonexpansion groups on the basis of the frequencies of their founder lineages. We observed a much higher N/S ratio in the expansion group, especially for nonsynonymous mutations with moderately deleterious effects, indicating a weaker effect of purifying selection in the expanded clades. However, this observation on N/S ratio was unlikely in computer simulations where all individuals were under the same selective constraints. Thus, we argue that the expanded populations were subjected to weaker selective constraints than the nonexpanded populations were. The mildly deleterious mutations were retained during population expansion, which could have a profound impact on present-day disease patterns.
Collapse
Affiliation(s)
- Hong-Xiang Zheng
- State Key Laboratory of Genetic Engineering, School of Life Sciences, Human Phenome Institute and Center for Evolutionary Biology, Fudan University, Shanghai, China
- Collaborative Innovation Center for Genetics and Development, Fudan University, Shanghai, China
| | - Shi Yan
- Ministry of Education Key Laboratory of Contemporary Anthropology, Department of Anthropology and Human Genetics, School of Life Sciences, Fudan University, Shanghai, China
- School of Ethnology and Sociology, Minzu University of China, Beijing, China
| | - Menghan Zhang
- Institute of Modern Languages and Linguistics, Fudan University, Shanghai, China
- Research Institute of Intelligent Complex Systems, Fudan University, Shanghai, China
| | - Zhenglong Gu
- State Key Laboratory of Genetic Engineering, School of Life Sciences, Human Phenome Institute and Center for Evolutionary Biology, Fudan University, Shanghai, China
- Greater Bay Area Institute of Precision Medicine (Guangzhou), Fudan University, Guangzhou, China
| | - Jiucun Wang
- State Key Laboratory of Genetic Engineering, School of Life Sciences, Human Phenome Institute and Center for Evolutionary Biology, Fudan University, Shanghai, China
- Collaborative Innovation Center for Genetics and Development, Fudan University, Shanghai, China
- Research Unit of Dissecting Population Genetics and Developing New Technologies for Treatment and Prevention of Skin Phenotypes and Dermatological Diseases (2019RU058), Chinese Academy of Medical Sciences, Beijing, China
| | - Li Jin
- State Key Laboratory of Genetic Engineering, School of Life Sciences, Human Phenome Institute and Center for Evolutionary Biology, Fudan University, Shanghai, China
- Collaborative Innovation Center for Genetics and Development, Fudan University, Shanghai, China
- Research Unit of Dissecting Population Genetics and Developing New Technologies for Treatment and Prevention of Skin Phenotypes and Dermatological Diseases (2019RU058), Chinese Academy of Medical Sciences, Beijing, China
| |
Collapse
|
2
|
Kyriazis CC, Lohmueller KE. Constraining models of dominance for nonsynonymous mutations in the human genome. PLoS Genet 2024; 20:e1011198. [PMID: 39302992 PMCID: PMC11446423 DOI: 10.1371/journal.pgen.1011198] [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/25/2024] [Revised: 10/02/2024] [Accepted: 09/04/2024] [Indexed: 09/22/2024] Open
Abstract
Dominance is a fundamental parameter in genetics, determining the dynamics of natural selection on deleterious and beneficial mutations, the patterns of genetic variation in natural populations, and the severity of inbreeding depression in a population. Despite this importance, dominance parameters remain poorly known, particularly in humans or other non-model organisms. A key reason for this lack of information about dominance is that it is extremely challenging to disentangle the selection coefficient (s) of a mutation from its dominance coefficient (h). Here, we explore dominance and selection parameters in humans by fitting models to the site frequency spectrum (SFS) for nonsynonymous mutations. When assuming a single dominance coefficient for all nonsynonymous mutations, we find that numerous h values can fit the data, so long as h is greater than ~0.15. Moreover, we also observe that theoretically-predicted models with a negative relationship between h and s can also fit the data well, including models with h = 0.05 for strongly deleterious mutations. Finally, we use our estimated dominance and selection parameters to inform simulations revisiting the question of whether the out-of-Africa bottleneck has led to differences in genetic load between African and non-African human populations. These simulations suggest that the relative burden of genetic load in non-African populations depends on the dominance model assumed, with slight increases for more weakly recessive models and slight decreases shown for more strongly recessive models. Moreover, these results also demonstrate that models of partially recessive nonsynonymous mutations can explain the observed severity of inbreeding depression in humans, bridging the gap between molecular population genetics and direct measures of fitness in humans. Our work represents a comprehensive assessment of dominance and deleterious variation in humans, with implications for parameterizing models of deleterious variation in humans and other mammalian species.
Collapse
Affiliation(s)
- Christopher C. Kyriazis
- Department of Ecology and Evolutionary Biology, University of California, Los Angeles, California, United States of America
| | - Kirk E. Lohmueller
- Department of Ecology and Evolutionary Biology, University of California, Los Angeles, California, United States of America
- Interdepartmental Program in Bioinformatics, University of California, Los Angeles, California, United States of America
- Department of Human Genetics, David Geffen School of Medicine, University of California, Los Angeles, California, United States of America
| |
Collapse
|
3
|
Bieniek-Kobuszewska M, Panasiewicz G. Polymorphism Identification in the Coding Sequences (ORFs) of the Porcine Pregnancy-Associated Glycoprotein 2-like Gene Subfamily in Pigs. Genes (Basel) 2024; 15:1149. [PMID: 39336740 PMCID: PMC11431107 DOI: 10.3390/genes15091149] [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/19/2024] [Revised: 08/22/2024] [Accepted: 08/27/2024] [Indexed: 09/30/2024] Open
Abstract
Pregnancy-associated glycoproteins (PAGs) are a polygenic family with many scattered genes and pseudogenes resulting from the duplication or fusion of a pseudogene with expression beginning in the trophoblast during the peri-implantation period and continuing in the trophectoderm. In this study, single-nucleotide polymorphism (SNP) and insertion/deletion (InDels) in the open reading frame (nine exons) of crossbreed pigs are reported for the first time. Novel SNPs/InDels were researched using genomic DNA templates isolated from the leukocytes of crossbreed pigs (N = 25), which were amplified, gel-out-purified, and sequenced. Sixteen SNPs and one InDel (g.6961_6966 Ins TGCCAA) were identified in the crossbreed pigs. In silico analysis revealed that among 16 SNPs, only 10 SNPs cause amino acid (aa) substitutions, and InDel codes asparagine (N298) and alanine (A299). The results provide a novel broad-based database (main pattern) that will be critical for future research into the possible correlations between the SNP genotypes of the pPAG2-L subfamily in pigs of various breeds whose reproductive traits are known.
Collapse
Affiliation(s)
- Martyna Bieniek-Kobuszewska
- Voivodeship Sanitary-Epidemiological Station in Olsztyn, Laboratory of Epidemiological and Clinical Research, Department of Virology and Serology, Zolnierska Str. 16, 10-561 Olsztyn, Poland;
| | - Grzegorz Panasiewicz
- Department of Animal Anatomy and Physiology, Faculty of Biology and Biotechnology, University of Warmia and Mazury in Olsztyn, Oczapowskiego Str. 1A, 10-719 Olsztyn, Poland
| |
Collapse
|
4
|
Jiang J, Xu YC, Zhang ZQ, Chen JF, Niu XM, Hou XH, Li XT, Wang L, Zhang YE, Ge S, Guo YL. Forces driving transposable element load variation during Arabidopsis range expansion. THE PLANT CELL 2024; 36:840-862. [PMID: 38036296 PMCID: PMC10980350 DOI: 10.1093/plcell/koad296] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/17/2023] [Revised: 10/25/2023] [Accepted: 11/06/2023] [Indexed: 12/02/2023]
Abstract
Genetic load refers to the accumulated and potentially life-threatening deleterious mutations in populations. Understanding the mechanisms underlying genetic load variation of transposable element (TE) insertion, a major large-effect mutation, during range expansion is an intriguing question in biology. Here, we used 1,115 global natural accessions of Arabidopsis (Arabidopsis thaliana) to study the driving forces of TE load variation during its range expansion. TE load increased with range expansion, especially in the recently established Yangtze River basin population. Effective population size, which explains 62.0% of the variance in TE load, high transposition rate, and selective sweeps contributed to TE accumulation in the expanded populations. We genetically mapped and identified multiple candidate causal genes and TEs, and revealed the genetic architecture of TE load variation. Overall, this study reveals the variation in TE genetic load during Arabidopsis expansion and highlights the causes of TE load variation from the perspectives of both population genetics and quantitative genetics.
Collapse
Affiliation(s)
- Juan Jiang
- State Key Laboratory of Systematic and Evolutionary Botany, Institute of Botany, Chinese Academy of Sciences, Beijing 100093, China
- China National Botanical Garden, Beijing 100093, China
- College of Life Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
| | - Yong-Chao Xu
- State Key Laboratory of Systematic and Evolutionary Botany, Institute of Botany, Chinese Academy of Sciences, Beijing 100093, China
- China National Botanical Garden, Beijing 100093, China
| | - Zhi-Qin Zhang
- State Key Laboratory of Systematic and Evolutionary Botany, Institute of Botany, Chinese Academy of Sciences, Beijing 100093, China
- China National Botanical Garden, Beijing 100093, China
- College of Life Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
| | - Jia-Fu Chen
- State Key Laboratory of Systematic and Evolutionary Botany, Institute of Botany, Chinese Academy of Sciences, Beijing 100093, China
- China National Botanical Garden, Beijing 100093, China
- College of Life Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
| | - Xiao-Min Niu
- State Key Laboratory of Systematic and Evolutionary Botany, Institute of Botany, Chinese Academy of Sciences, Beijing 100093, China
- China National Botanical Garden, Beijing 100093, China
| | - Xing-Hui Hou
- State Key Laboratory of Systematic and Evolutionary Botany, Institute of Botany, Chinese Academy of Sciences, Beijing 100093, China
- China National Botanical Garden, Beijing 100093, China
| | - Xin-Tong Li
- State Key Laboratory of Systematic and Evolutionary Botany, Institute of Botany, Chinese Academy of Sciences, Beijing 100093, China
- China National Botanical Garden, Beijing 100093, China
- College of Life Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
| | - Li Wang
- Agricultural Synthetic Biology Center, Agricultural Genomics Institute at Shenzhen, Chinese Academy of Agricultural Sciences, Shenzhen 518000, China
| | - Yong E Zhang
- College of Life Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
- State Key Laboratory of Integrated Management of Pest Insects and Rodents & Key Laboratory of the Zoological Systematics and Evolution, Institute of Zoology, Chinese Academy of Sciences, Beijing 100101, China
| | - Song Ge
- State Key Laboratory of Systematic and Evolutionary Botany, Institute of Botany, Chinese Academy of Sciences, Beijing 100093, China
- China National Botanical Garden, Beijing 100093, China
- College of Life Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
| | - Ya-Long Guo
- State Key Laboratory of Systematic and Evolutionary Botany, Institute of Botany, Chinese Academy of Sciences, Beijing 100093, China
- China National Botanical Garden, Beijing 100093, China
- College of Life Sciences, University of Chinese Academy of Sciences, Beijing 100049, China
| |
Collapse
|
5
|
La Rocca LA, Frank J, Bentzen HB, Pantel JT, Gerischer K, Bovier A, Krawitz PM. Understanding recessive disease risk in multi-ethnic populations with different degrees of consanguinity. Am J Med Genet A 2024; 194:e63452. [PMID: 37921563 DOI: 10.1002/ajmg.a.63452] [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: 09/30/2023] [Accepted: 10/10/2023] [Indexed: 11/04/2023]
Abstract
Population medical genetics aims at translating clinically relevant findings from recent studies of large cohorts into healthcare for individuals. Genetic counseling concerning reproductive risks and options is still mainly based on family history, and consanguinity is viewed to increase the risk for recessive diseases regardless of the demographics. However, in an increasingly multi-ethnic society with diverse approaches to partner selection, healthcare professionals should also sharpen their intuition for the influence of different mating schemes in non-equilibrium dynamics. We, therefore, revisited the so-called out-of-Africa model and studied in forward simulations with discrete and not overlapping generations the effect of inbreeding on the average number of recessive lethals in the genome. We were able to reproduce in both frameworks the drop in the incidence of recessive disorders, which is a transient phenomenon during and after the growth phase of a population, and therefore showed their equivalence. With the simulation frameworks, we also provide the means to study and visualize the effect of different kin sizes and mating schemes on these parameters for educational purposes.
Collapse
Affiliation(s)
- Luis A La Rocca
- Institute for Applied Mathematics, University of Bonn, Bonn, Germany
| | - Julia Frank
- Institute for Genomic Statistics and Bioinformatics, University of Bonn, Bonn, Germany
| | - Heidi Beate Bentzen
- Centre for Medical Ethics, Faculty of Medicine, Univeristy of Oslo, Oslo, Norway
| | - Jean Tori Pantel
- Department of Digitalization and General Practice, University Hospital RWTH Aachen, Aachen, Germany
| | - Konrad Gerischer
- Institute for Genomic Statistics and Bioinformatics, University of Bonn, Bonn, Germany
| | - Anton Bovier
- Institute for Genomic Statistics and Bioinformatics, University of Bonn, Bonn, Germany
| | - Peter M Krawitz
- Institute for Applied Mathematics, University of Bonn, Bonn, Germany
| |
Collapse
|
6
|
Buffalo V, Kern AD. A quantitative genetic model of background selection in humans. PLoS Genet 2024; 20:e1011144. [PMID: 38507461 PMCID: PMC10984650 DOI: 10.1371/journal.pgen.1011144] [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: 09/17/2023] [Revised: 04/01/2024] [Accepted: 01/19/2024] [Indexed: 03/22/2024] Open
Abstract
Across the human genome, there are large-scale fluctuations in genetic diversity caused by the indirect effects of selection. This "linked selection signal" reflects the impact of selection according to the physical placement of functional regions and recombination rates along chromosomes. Previous work has shown that purifying selection acting against the steady influx of new deleterious mutations at functional portions of the genome shapes patterns of genomic variation. To date, statistical efforts to estimate purifying selection parameters from linked selection models have relied on classic Background Selection theory, which is only applicable when new mutations are so deleterious that they cannot fix in the population. Here, we develop a statistical method based on a quantitative genetics view of linked selection, that models how polygenic additive fitness variance distributed along the genome increases the rate of stochastic allele frequency change. By jointly predicting the equilibrium fitness variance and substitution rate due to both strong and weakly deleterious mutations, we estimate the distribution of fitness effects (DFE) and mutation rate across three geographically distinct human samples. While our model can accommodate weaker selection, we find evidence of strong selection operating similarly across all human samples. Although our quantitative genetic model of linked selection fits better than previous models, substitution rates of the most constrained sites disagree with observed divergence levels. We find that a model incorporating selective interference better predicts observed divergence in conserved regions, but overall our results suggest uncertainty remains about the processes generating fitness variation in humans.
Collapse
Affiliation(s)
- Vince Buffalo
- Department of Integrative Biology, University of California, Berkeley, Berkeley, California, United States of America
- Institute of Ecology and Evolution and Department of Biology, University of Oregon, Eugene, Oregon, United States of America
| | - Andrew D. Kern
- Institute of Ecology and Evolution and Department of Biology, University of Oregon, Eugene, Oregon, United States of America
| |
Collapse
|
7
|
Kyriazis CC, Lohmueller KE. Constraining models of dominance for nonsynonymous mutations in the human genome. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.02.25.582010. [PMID: 38463985 PMCID: PMC10925099 DOI: 10.1101/2024.02.25.582010] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 03/12/2024]
Abstract
Dominance is a fundamental parameter in genetics, determining the dynamics of natural selection on deleterious and beneficial mutations, the patterns of genetic variation in natural populations, and the severity of inbreeding depression in a population. Despite this importance, dominance parameters remain poorly known, particularly in humans or other non-model organisms. A key reason for this lack of information about dominance is that it is extremely challenging to disentangle the selection coefficient (s) of a mutation from its dominance coefficient (h). Here, we explore dominance and selection parameters in humans by fitting models to the site frequency spectrum (SFS) for nonsynonymous mutations. When assuming a single dominance coefficient for all nonsynonymous mutations, we find that numerous h values can fit the data, so long as h is greater than ~0.15. Moreover, we also observe that theoretically-predicted models with a negative relationship between h and s can also fit the data well, including models with h=0.05 for strongly deleterious mutations. Finally, we use our estimated dominance and selection parameters to inform simulations revisiting the question of whether the out-of-Africa bottleneck has led to differences in genetic load between African and non-African human populations. These simulations suggest that the relative burden of genetic load in non-African populations depends on the dominance model assumed, with slight increases for more weakly recessive models and slight decreases shown for more strongly recessive models. Moreover, these results also demonstrate that models of partially recessive nonsynonymous mutations can explain the observed severity of inbreeding depression in humans, bridging the gap between molecular population genetics and direct measures of fitness in humans. Our work represents a comprehensive assessment of dominance and deleterious variation in humans, with implications for parameterizing models of deleterious variation in humans and other mammalian species.
Collapse
Affiliation(s)
| | - Kirk E. Lohmueller
- Department of Ecology and Evolutionary Biology, University of California, Los Angeles, USA
- Interdepartmental Program in Bioinformatics, University of California, Los Angeles, USA
- Department of Human Genetics, David Geffen School of Medicine, Los Angeles, USA
| |
Collapse
|
8
|
Barili V, Ambrosini E, Bortesi B, Minari R, De Sensi E, Cannizzaro IR, Taiani A, Michiara M, Sikokis A, Boggiani D, Tommasi C, Serra O, Bonatti F, Adorni A, Luberto A, Caggiati P, Martorana D, Uliana V, Percesepe A, Musolino A, Pellegrino B. Genetic Basis of Breast and Ovarian Cancer: Approaches and Lessons Learnt from Three Decades of Inherited Predisposition Testing. Genes (Basel) 2024; 15:219. [PMID: 38397209 PMCID: PMC10888198 DOI: 10.3390/genes15020219] [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: 12/30/2023] [Revised: 02/02/2024] [Accepted: 02/05/2024] [Indexed: 02/25/2024] Open
Abstract
Germline variants occurring in BRCA1 and BRCA2 give rise to hereditary breast and ovarian cancer (HBOC) syndrome, predisposing to breast, ovarian, fallopian tube, and peritoneal cancers marked by elevated incidences of genomic aberrations that correspond to poor prognoses. These genes are in fact involved in genetic integrity, particularly in the process of homologous recombination (HR) DNA repair, a high-fidelity repair system for mending DNA double-strand breaks. In addition to its implication in HBOC pathogenesis, the impairment of HR has become a prime target for therapeutic intervention utilizing poly (ADP-ribose) polymerase (PARP) inhibitors. In the present review, we introduce the molecular roles of HR orchestrated by BRCA1 and BRCA2 within the framework of sensitivity to PARP inhibitors. We examine the genetic architecture underneath breast and ovarian cancer ranging from high- and mid- to low-penetrant predisposing genes and taking into account both germline and somatic variations. Finally, we consider higher levels of complexity of the genomic landscape such as polygenic risk scores and other approaches aiming to optimize therapeutic and preventive strategies for breast and ovarian cancer.
Collapse
Affiliation(s)
- Valeria Barili
- Department of Medicine and Surgery, University of Parma, 43126 Parma, Italy
| | - Enrico Ambrosini
- Medical Genetics, University Hospital of Parma, 43126 Parma, Italy
| | - Beatrice Bortesi
- Medical Oncology Unit, University Hospital of Parma, 43126 Parma, Italy
| | - Roberta Minari
- Medical Oncology Unit, University Hospital of Parma, 43126 Parma, Italy
| | - Erika De Sensi
- Department of Medicine and Surgery, University of Parma, 43126 Parma, Italy
| | | | - Antonietta Taiani
- Department of Medicine and Surgery, University of Parma, 43126 Parma, Italy
| | - Maria Michiara
- Medical Oncology Unit, University Hospital of Parma, 43126 Parma, Italy
- Breast Unit, University Hospital of Parma, 43126 Parma, Italy
| | - Angelica Sikokis
- Medical Oncology Unit, University Hospital of Parma, 43126 Parma, Italy
- Breast Unit, University Hospital of Parma, 43126 Parma, Italy
| | - Daniela Boggiani
- Medical Oncology Unit, University Hospital of Parma, 43126 Parma, Italy
- Breast Unit, University Hospital of Parma, 43126 Parma, Italy
| | - Chiara Tommasi
- Department of Medicine and Surgery, University of Parma, 43126 Parma, Italy
- Medical Oncology Unit, University Hospital of Parma, 43126 Parma, Italy
- Breast Unit, University Hospital of Parma, 43126 Parma, Italy
| | - Olga Serra
- Medical Oncology Unit, University Hospital of Parma, 43126 Parma, Italy
- Breast Unit, University Hospital of Parma, 43126 Parma, Italy
| | - Francesco Bonatti
- Medical Oncology Unit, University Hospital of Parma, 43126 Parma, Italy
| | - Alessia Adorni
- Medical Oncology Unit, University Hospital of Parma, 43126 Parma, Italy
| | - Anita Luberto
- Department of Medicine and Surgery, University of Parma, 43126 Parma, Italy
| | | | - Davide Martorana
- Medical Genetics, University Hospital of Parma, 43126 Parma, Italy
| | - Vera Uliana
- Medical Genetics, University Hospital of Parma, 43126 Parma, Italy
| | - Antonio Percesepe
- Department of Medicine and Surgery, University of Parma, 43126 Parma, Italy
- Medical Genetics, University Hospital of Parma, 43126 Parma, Italy
| | - Antonino Musolino
- Department of Medicine and Surgery, University of Parma, 43126 Parma, Italy
- Medical Oncology Unit, University Hospital of Parma, 43126 Parma, Italy
- Breast Unit, University Hospital of Parma, 43126 Parma, Italy
| | - Benedetta Pellegrino
- Medical Oncology Unit, University Hospital of Parma, 43126 Parma, Italy
- Breast Unit, University Hospital of Parma, 43126 Parma, Italy
| |
Collapse
|
9
|
Lucas-Sánchez M, Abdeli A, Bekada A, Calafell F, Benhassine T, Comas D. The Impact of Recent Demography on Functional Genetic Variation in North African Human Groups. Mol Biol Evol 2024; 41:msad283. [PMID: 38152862 PMCID: PMC10783648 DOI: 10.1093/molbev/msad283] [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: 05/22/2023] [Revised: 11/22/2023] [Accepted: 12/19/2023] [Indexed: 12/29/2023] Open
Abstract
The strategic location of North Africa has made the region the core of a wide range of human demographic events, including migrations, bottlenecks, and admixture processes. This has led to a complex and heterogeneous genetic and cultural landscape, which remains poorly studied compared to other world regions. Whole-exome sequencing is particularly relevant to determine the effects of these demographic events on current-day North Africans' genomes, since it allows to focus on those parts of the genome that are more likely to have direct biomedical consequences. Whole-exome sequencing can also be used to assess the effect of recent demography in functional genetic variation and the efficacy of natural selection, a long-lasting debate. In the present work, we use newly generated whole-exome sequencing and genome-wide array genotypes to investigate the effect of demography in functional variation in 7 North African populations, considering both cultural and demographic differences and with a special focus on Amazigh (plur. Imazighen) groups. We detect genetic differences among populations related to their degree of isolation and the presence of bottlenecks in their recent history. We find differences in the functional part of the genome that suggest a relaxation of purifying selection in the more isolated groups, allowing for an increase of putatively damaging variation. Our results also show a shift in mutational load coinciding with major demographic events in the region and reveal differences within and between cultural and geographic groups.
Collapse
Affiliation(s)
- Marcel Lucas-Sánchez
- Departament de Medicina i Ciències de la Vida, Institut de Biologia Evolutiva (CSIC-UPF), Universitat Pompeu Fabra, Barcelona, Spain
| | - Amine Abdeli
- Faculté des Sciences Biologiques, Laboratoire de Biologie Cellulaire et Moléculaire, Université des Sciences et de la Technologie Houari Boumediene, Alger, Algeria
| | - Asmahan Bekada
- Département de Biotechnologie, Faculté des Sciences de la Nature et de la Vie, Université Oran 1 (Ahmad Ben Bella), Oran, Algeria
| | - Francesc Calafell
- Departament de Medicina i Ciències de la Vida, Institut de Biologia Evolutiva (CSIC-UPF), Universitat Pompeu Fabra, Barcelona, Spain
| | - Traki Benhassine
- Faculté des Sciences Biologiques, Laboratoire de Biologie Cellulaire et Moléculaire, Université des Sciences et de la Technologie Houari Boumediene, Alger, Algeria
| | - David Comas
- Departament de Medicina i Ciències de la Vida, Institut de Biologia Evolutiva (CSIC-UPF), Universitat Pompeu Fabra, Barcelona, Spain
| |
Collapse
|
10
|
Kim AS, Kreiner JM, Hernández F, Bock DG, Hodgins KA, Rieseberg LH. Temporal collections to study invasion biology. Mol Ecol 2023; 32:6729-6742. [PMID: 37873879 DOI: 10.1111/mec.17176] [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: 08/22/2023] [Revised: 10/07/2023] [Accepted: 10/11/2023] [Indexed: 10/25/2023]
Abstract
Biological invasions represent an extraordinary opportunity to study evolution. This is because accidental or deliberate species introductions have taken place for centuries across large geographical scales, frequently prompting rapid evolutionary transitions in invasive populations. Until recently, however, the utility of invasions as evolutionary experiments has been hampered by limited information on the makeup of populations that were part of earlier invasion stages. Now, developments in ancient and historical DNA technologies, as well as the quickening pace of digitization for millions of specimens that are housed in herbaria and museums globally, promise to help overcome this obstacle. In this review, we first introduce the types of temporal data that can be used to study invasions, highlighting the timescale captured by each approach and their respective limitations. We then discuss how ancient and historical specimens as well as data available from prior invasion studies can be used to answer questions on mechanisms of (mal)adaptation, rates of evolution, or community-level changes during invasions. By bridging the gap between contemporary and historical invasive populations, temporal data can help us connect pattern to process in invasion science. These data will become increasingly important if invasions are to achieve their full potential as experiments of evolution in nature.
Collapse
Affiliation(s)
- Amy S Kim
- Department of Botany and Biodiversity Research Centre, University of British Columbia, Vancouver, British Columbia, Canada
| | - Julia M Kreiner
- Department of Botany and Biodiversity Research Centre, University of British Columbia, Vancouver, British Columbia, Canada
| | - Fernando Hernández
- Department of Botany and Biodiversity Research Centre, University of British Columbia, Vancouver, British Columbia, Canada
| | - Dan G Bock
- Department of Botany and Biodiversity Research Centre, University of British Columbia, Vancouver, British Columbia, Canada
| | - Kathryn A Hodgins
- School of Biological Sciences, Monash University, Melbourne, Victoria, Australia
| | - Loren H Rieseberg
- Department of Botany and Biodiversity Research Centre, University of British Columbia, Vancouver, British Columbia, Canada
| |
Collapse
|
11
|
Swinford NA, Prall SP, Gopalan S, Williams CM, Sheehama J, Scelza BA, Henn BM. Increased homozygosity due to endogamy results in fitness consequences in a human population. Proc Natl Acad Sci U S A 2023; 120:e2309552120. [PMID: 37847737 PMCID: PMC10614605 DOI: 10.1073/pnas.2309552120] [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: 06/08/2023] [Accepted: 09/18/2023] [Indexed: 10/19/2023] Open
Abstract
Recessive alleles have been shown to directly affect both human Mendelian disease phenotypes and complex traits. Pedigree studies also suggest that consanguinity results in increased childhood mortality and adverse health phenotypes, presumably through penetrance of recessive mutations. Here, we test whether the accumulation of homozygous, recessive alleles decreases reproductive success in a human population. We address this question among the Namibian Himba, an endogamous agro-pastoralist population, who until very recently practiced natural fertility. Using a sample of 681 individuals, we show that Himba exhibit elevated levels of "inbreeding," calculated as the fraction of the genome in runs of homozygosity (FROH). Many individuals contain multiple long segments of ROH in their genomes, indicating that their parents had high kinship coefficients. However, we do not find evidence that this is explained by first-cousin consanguinity, despite a reported social preference for cross-cousin marriages. Rather, we show that elevated haplotype sharing in the Himba is due to a bottleneck, likely in the past 60 generations. We test whether increased recessive mutation load results in observed fitness consequences by assessing the effect of FROH on completed fertility in a cohort of postreproductive women (n = 69). We find that higher FROH is significantly associated with lower fertility. Our data suggest a multilocus genetic effect on fitness driven by the expression of deleterious recessive alleles, especially those in long ROH. However, these effects are not the result of consanguinity but rather elevated background identity by descent.
Collapse
Affiliation(s)
- N. A. Swinford
- Department of Anthropology, University of California Davis, Davis, CA95616
| | - S. P. Prall
- Department of Anthropology, University of Missouri, Columbia, MO65211
| | - S. Gopalan
- Department of Evolutionary Anthropology, Duke University, Durham, NC27708
| | - C. M. Williams
- Center for Computational Molecular Biology, Brown University, Providence, RI02912
| | - J. Sheehama
- Department of Human, Biological and Translational Medical Sciences, School of Medicine University of Namibia, Oshakati10005, Namibia
| | - B. A. Scelza
- Department of Anthropology, University of California Los Angeles, Los Angeles, CA90095
| | - B. M. Henn
- Department of Anthropology, University of California Davis, Davis, CA95616
- Center for Population Biology, University of California Davis, Davis, CA95616
- Genome Center, University of California Davis, Davis, CA95616
| |
Collapse
|
12
|
Nigenda-Morales SF, Lin M, Nuñez-Valencia PG, Kyriazis CC, Beichman AC, Robinson JA, Ragsdale AP, Urbán R J, Archer FI, Viloria-Gómora L, Pérez-Álvarez MJ, Poulin E, Lohmueller KE, Moreno-Estrada A, Wayne RK. The genomic footprint of whaling and isolation in fin whale populations. Nat Commun 2023; 14:5465. [PMID: 37699896 PMCID: PMC10497599 DOI: 10.1038/s41467-023-40052-z] [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: 06/08/2022] [Accepted: 07/10/2023] [Indexed: 09/14/2023] Open
Abstract
Twentieth century industrial whaling pushed several species to the brink of extinction, with fin whales being the most impacted. However, a small, resident population in the Gulf of California was not targeted by whaling. Here, we analyzed 50 whole-genomes from the Eastern North Pacific (ENP) and Gulf of California (GOC) fin whale populations to investigate their demographic history and the genomic effects of natural and human-induced bottlenecks. We show that the two populations diverged ~16,000 years ago, after which the ENP population expanded and then suffered a 99% reduction in effective size during the whaling period. In contrast, the GOC population remained small and isolated, receiving less than one migrant per generation. However, this low level of migration has been crucial for maintaining its viability. Our study exposes the severity of whaling, emphasizes the importance of migration, and demonstrates the use of genome-based analyses and simulations to inform conservation strategies.
Collapse
Affiliation(s)
- Sergio F Nigenda-Morales
- Advanced Genomics Unit, National Laboratory of Genomics for Biodiversity (Langebio), Center for Research and Advanced Studies (Cinvestav), Irapuato, Guanajuato, 36824, Mexico.
- Department of Biological Sciences, California State University San Marcos, San Marcos, CA, 92096, USA.
| | - Meixi Lin
- Department of Ecology and Evolutionary Biology, University of California, Los Angeles, Los Angeles, CA, 90095, USA.
- Department of Plant Biology, Carnegie Institution for Science, Stanford, CA, 94305, USA.
| | - Paulina G Nuñez-Valencia
- Advanced Genomics Unit, National Laboratory of Genomics for Biodiversity (Langebio), Center for Research and Advanced Studies (Cinvestav), Irapuato, Guanajuato, 36824, Mexico
- Centro de Ciencias Genómicas, Universidad Nacional Autónoma de México (UNAM), Cuernavaca, Morelos, México
| | - Christopher C Kyriazis
- Department of Ecology and Evolutionary Biology, University of California, Los Angeles, Los Angeles, CA, 90095, USA
| | - Annabel C Beichman
- Department of Genome Sciences, University of Washington, Seattle, WA, 98195, USA
| | - Jacqueline A Robinson
- Institute for Human Genetics, University of California, San Francisco (UCSF), San Francisco, CA, 94143, USA
| | - Aaron P Ragsdale
- Advanced Genomics Unit, National Laboratory of Genomics for Biodiversity (Langebio), Center for Research and Advanced Studies (Cinvestav), Irapuato, Guanajuato, 36824, Mexico
- Department of Integrative Biology, University of Wisconsin, Madison, WI, 53706, USA
| | - Jorge Urbán R
- Departamento de Ciencias Marinas y Costeras, Universidad Autónoma de Baja California Sur (UABCS), La Paz, Baja California Sur, Mexico
| | - Frederick I Archer
- Marine Mammal and Turtle Division, Southwest Fisheries Science Center, La Jolla, CA, 92037, USA
| | - Lorena Viloria-Gómora
- Departamento de Ciencias Marinas y Costeras, Universidad Autónoma de Baja California Sur (UABCS), La Paz, Baja California Sur, Mexico
| | - María José Pérez-Álvarez
- Escuela de Medicina Veterinaria, Facultad de Medicina y Ciencias de la Salud, Universidad Mayor, Santiago, Chile
- Millennium Institute Biodiversity of Antarctic and Subantarctic Ecosystems (BASE), Universidad de Chile, Santiago, Chile
| | - Elie Poulin
- Millennium Institute Biodiversity of Antarctic and Subantarctic Ecosystems (BASE), Universidad de Chile, Santiago, Chile
| | - Kirk E Lohmueller
- Department of Ecology and Evolutionary Biology, University of California, Los Angeles, Los Angeles, CA, 90095, USA.
- Interdepartmental Program in Bioinformatics, University of California, Los Angeles, CA, 90095, USA.
- Department of Human Genetics, David Geffen School of Medicine, University of California, Los Angeles, CA, USA.
| | - Andrés Moreno-Estrada
- Advanced Genomics Unit, National Laboratory of Genomics for Biodiversity (Langebio), Center for Research and Advanced Studies (Cinvestav), Irapuato, Guanajuato, 36824, Mexico.
| | - Robert K Wayne
- Department of Ecology and Evolutionary Biology, University of California, Los Angeles, Los Angeles, CA, 90095, USA
| |
Collapse
|
13
|
Matheson J, Bertram J, Masel J. Human deleterious mutation rate implies high fitness variance, with declining mean fitness compensated by rarer beneficial mutations of larger effect. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.09.01.555871. [PMID: 37732183 PMCID: PMC10508744 DOI: 10.1101/2023.09.01.555871] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 09/22/2023]
Abstract
Each new human has an expected Ud = 2 - 10 new deleterious mutations. This deluge of deleterious mutations cannot all be purged, and therefore accumulate in a declining fitness ratchet. Using a novel simulation framework designed to efficiently handle genome-wide linkage disequilibria across many segregating sites, we find that rarer, beneficial mutations of larger effect are sufficient to compensate fitness declines due to the fixation of many slightly deleterious mutations. Drift barrier theory posits a similar asymmetric pattern of fixations to explain ratcheting genome size and complexity, but in our theory, the cause is Ud > 1 rather than small population size. In our simulations, Ud ~2 - 10 generates high within-population variance in relative fitness; two individuals will typically differ in fitness by 15-40%. Ud ~2 - 10 also slows net adaptation by ~13%-39%. Surprisingly, fixation rates are more sensitive to changes in the beneficial than the deleterious mutation rate, e.g. a 10% increase in overall mutation rate leads to faster adaptation; this puts to rest dysgenic fears about increasing mutation rates due to rising paternal age.
Collapse
Affiliation(s)
- Joseph Matheson
- Department of Ecology and Evolutionary Biology, University of Arizona, Tucson, AZ, 85721, USA
- Department of Ecology, Behavior, and Evolution, University of California San Diego, San Diego, CA, 92093, USA
| | - Jason Bertram
- Department of Ecology and Evolutionary Biology, University of Arizona, Tucson, AZ, 85721, USA
- Department of Mathematics, University of Western Ontario, London ON, Canada
| | - Joanna Masel
- Department of Ecology and Evolutionary Biology, University of Arizona, Tucson, AZ, 85721, USA
| |
Collapse
|
14
|
Antinucci M, Comas D, Calafell F. Population history modulates the fitness effects of Copy Number Variation in the Roma. Hum Genet 2023; 142:1327-1343. [PMID: 37311904 PMCID: PMC10449987 DOI: 10.1007/s00439-023-02579-5] [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: 04/17/2023] [Accepted: 06/02/2023] [Indexed: 06/15/2023]
Abstract
We provide the first whole genome Copy Number Variant (CNV) study addressing Roma, along with reference populations from South Asia, the Middle East and Europe. Using CNV calling software for short-read sequence data, we identified 3171 deletions and 489 duplications. Taking into account the known population history of the Roma, as inferred from whole genome nucleotide variation, we could discern how this history has shaped CNV variation. As expected, patterns of deletion variation, but not duplication, in the Roma followed those obtained from single nucleotide polymorphisms (SNPs). Reduced effective population size resulting in slightly relaxed natural selection may explain our observation of an increase in intronic (but not exonic) deletions within Loss of Function (LoF)-intolerant genes. Over-representation analysis for LoF-intolerant gene sets hosting intronic deletions highlights a substantial accumulation of shared biological processes in Roma, intriguingly related to signaling, nervous system and development features, which may be related to the known profile of private disease in the population. Finally, we show the link between deletions and known trait-related SNPs reported in the genome-wide association study (GWAS) catalog, which exhibited even frequency distributions among the studied populations. This suggests that, in general human populations, the strong association between deletions and SNPs associated to biomedical conditions and traits could be widespread across continental populations, reflecting a common background of potentially disease/trait-related CNVs.
Collapse
Affiliation(s)
- Marco Antinucci
- Institute of Evolutionary Biology (UPF-CSIC), Department of Medicine and Life Sciences, Universitat Pompeu Fabra, Barcelona, Spain
| | - David Comas
- Institute of Evolutionary Biology (UPF-CSIC), Department of Medicine and Life Sciences, Universitat Pompeu Fabra, Barcelona, Spain
| | - Francesc Calafell
- Institute of Evolutionary Biology (UPF-CSIC), Department of Medicine and Life Sciences, Universitat Pompeu Fabra, Barcelona, Spain.
| |
Collapse
|
15
|
Kumar M, Conroy G, Ogbourne S, Cairns K, Borburgh L, Subramanian S. Genomic signatures of bottleneck and founder effects in dingoes. Ecol Evol 2023; 13:e10525. [PMID: 37732287 PMCID: PMC10508967 DOI: 10.1002/ece3.10525] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/28/2023] [Revised: 07/29/2023] [Accepted: 08/30/2023] [Indexed: 09/22/2023] Open
Abstract
Dingoes arrived in Australia during the mid-Holocene and are the top-order terrestrial predator on the continent. Although dingoes subsequently spread across the continent, the initial founding population(s) could have been small. We investigated this hypothesis by sequencing the whole genomes of three dingoes and also obtaining the genome data from nine additional dingoes and 56 canines, including wolves, village dogs and breed dogs, and examined the signatures of bottlenecks and founder effects. We found that the nucleotide diversity of dingoes was low, 36% less than highly inbred breed dogs and 3.3 times lower than wolves. The number of runs of homozygosity (RoH) segments in dingoes was 1.6-4.7 times higher than in other canines. While examining deleterious mutational load, we observed that dingoes carried elevated ratios of nonsynonymous-to-synonymous diversities, significantly higher numbers of homozygous deleterious Single Nucleotide Variants (SNVs), and increased numbers of loss of function SNVs, compared to breed dogs, village dogs, and wolves. Our findings can be explained by bottlenecks and founder effects during the establishment of dingoes in mainland Australia. These findings highlight the need for conservation-based management of dingoes and the need for wildlife managers to be cognisant of these findings when considering the use of lethal control measures across the landscape.
Collapse
Affiliation(s)
- Manoharan Kumar
- School of Science, Technology, and EngineeringThe University of the Sunshine CoastMoreton BayQueenslandAustralia
| | - Gabriel Conroy
- Centre for BioinnovationThe University of the Sunshine CoastSippy DownsQueenslandAustralia
- School of Science, Technology, and EngineeringThe University of the Sunshine CoastSippy DownsQueenslandAustralia
| | - Steven Ogbourne
- Centre for BioinnovationThe University of the Sunshine CoastSippy DownsQueenslandAustralia
| | - Kylie Cairns
- Evolution and Ecology Research Centre, School of Biological, Earth and Environmental SciencesUNSW AustraliaSydneyNew South WalesAustralia
- Centre for Ecosystem Science, School of Biological, Earth and Environmental SciencesUNSW AustraliaSydneyNew South WalesAustralia
| | - Liesbeth Borburgh
- School of Science, Technology, and EngineeringThe University of the Sunshine CoastSippy DownsQueenslandAustralia
| | - Sankar Subramanian
- School of Science, Technology, and EngineeringThe University of the Sunshine CoastMoreton BayQueenslandAustralia
- Centre for BioinnovationThe University of the Sunshine CoastSippy DownsQueenslandAustralia
| |
Collapse
|
16
|
Ramasamy U, Elizur A, Subramanian S. Deleterious mutation load in the admixed mice population. Front Ecol Evol 2023. [DOI: 10.3389/fevo.2023.1084502] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/18/2023] Open
Abstract
Deleterious mutation loads are known to correlate negatively with effective population size (Ne). Due to this reason, previous studies observed a higher proportion of harmful mutations in small populations than that in large populations. However, the mutational load in an admixed population that derived from introgression between individuals from two populations with vastly different Ne is not known. We investigated this using the whole genome data from two subspecies of the mouse (Mus musculus castaneus and Mus musculus musculus) with significantly different Ne. We used the ratio of diversities at nonsynonymous and synonymous sites (dN/dS) to measure the harmful mutation load. Our results showed that this ratio observed for the admixed population was intermediate between those of the parental populations. The dN/dS ratio of the hybrid population was significantly higher than that of M. m. castaneus but lower than that of M. m. musculus. Our analysis revealed a significant positive correlation between the proportion of M. m. musculus ancestry in admixed individuals and their dN/dS ratio. This suggests that the admixed individuals with high proportions of M. m. musculus ancestry have large dN/dS ratios. We also used the proportion of deleterious nonsynonymous SNVs as a proxy for deleterious mutation load, which also produced similar results. The observed results were in concordance with those expected by theory. We also show a shift in the distribution of fitness effects of nonsynonymous SNVs in the admixed genomes compared to the parental populations. These findings suggest that the deleterious mutation load of the admixed population is determined by the proportion of the ancestries of the subspecies. Therefore, it is important to consider the status and the level of genetic admixture of the populations whilst estimating the mutation loads.
Collapse
|
17
|
Mooney JA, Marsden CD, Yohannes A, Wayne RK, Lohmueller KE. Long-term Small Population Size, Deleterious Variation, and Altitude Adaptation in the Ethiopian Wolf, a Severely Endangered Canid. Mol Biol Evol 2023; 40:msac277. [PMID: 36585842 PMCID: PMC9847632 DOI: 10.1093/molbev/msac277] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2022] [Revised: 11/07/2022] [Accepted: 12/22/2022] [Indexed: 01/01/2023] Open
Abstract
Ethiopian wolves, a canid species endemic to the Ethiopian Highlands, have been steadily declining in numbers for decades. Currently, out of 35 extant species, it is now one of the world's most endangered canids. Most conservation efforts have focused on preventing disease, monitoring movements and behavior, and assessing the geographic ranges of sub-populations. Here, we add an essential layer by determining the Ethiopian wolf's demographic and evolutionary history using high-coverage (∼40×) whole-genome sequencing from 10 Ethiopian wolves from the Bale Mountains. We observe exceptionally low diversity and enrichment of weakly deleterious variants in the Ethiopian wolves in comparison with two North American gray wolf populations and four dog breeds. These patterns are consequences of long-term small population size, rather than recent inbreeding. We infer the demographic history of the Ethiopian wolf and find it to be concordant with historic records and previous genetic analyses, suggesting Ethiopian wolves experienced a series of both ancient and recent bottlenecks, resulting in a census population size of fewer than 500 individuals and an estimated effective population size of approximately 100 individuals. Additionally, long-term small population size may have limited the accumulation of strongly deleterious recessive mutations. Finally, as the Ethiopian wolves have inhabited high-altitude areas for thousands of years, we searched for evidence of high-altitude adaptation, finding evidence of positive selection at a transcription factor in a hypoxia-response pathway [CREB-binding protein (CREBBP)]. Our findings are pertinent to continuing conservation efforts and understanding how demography influences the persistence of deleterious variation in small populations.
Collapse
Affiliation(s)
- Jazlyn A Mooney
- Department of Human Genetics, University of California Los Angeles, Los Angeles, CA 90095, USA
- Department of Biology, Stanford University, Stanford, CA, USA
- Department of Quantitative and Computational Biology, University of Southern California, Los Angeles, CA, USA
| | - Clare D Marsden
- Department of Ecology & Evolutionary Biology, University of California Los Angeles, Los Angeles, CA, USA
| | - Abigail Yohannes
- Department of Chemistry and Biochemistry, University of California Los Angeles, Los Angeles, CA, USA
| | - Robert K Wayne
- Department of Ecology & Evolutionary Biology, University of California Los Angeles, Los Angeles, CA, USA
| | - Kirk E Lohmueller
- Department of Human Genetics, University of California Los Angeles, Los Angeles, CA 90095, USA
- Department of Ecology & Evolutionary Biology, University of California Los Angeles, Los Angeles, CA, USA
| |
Collapse
|
18
|
Aguirre-Liguori JA, Morales-Cruz A, Gaut BS. Evaluating the persistence and utility of five wild Vitis species in the context of climate change. Mol Ecol 2022; 31:6457-6472. [PMID: 36197804 PMCID: PMC10092629 DOI: 10.1111/mec.16715] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2022] [Revised: 09/23/2022] [Accepted: 09/29/2022] [Indexed: 01/13/2023]
Abstract
Crop wild relatives (CWRs) have the capacity to contribute novel traits to agriculture. Given climate change, these contributions may be especially vital for the persistence of perennial crops, because perennials are often clonally propagated and consequently do not evolve rapidly. By studying the landscape genomics of samples from five Vitis CWRs (V. arizonica, V. mustangensis, V. riparia, V. berlandieri and V. girdiana) in the context of projected climate change, we addressed two goals. The first was to assess the relative potential of different CWR accessions to persist in the face of climate change. By integrating species distribution models with adaptive genetic variation, additional genetic features such as genomic load and a phenotype (resistance to Pierce's Disease), we predicted that accessions from one species (V. mustangensis) are particularly well-suited to persist in future climates. The second goal was to identify which CWR accessions may contribute to bioclimatic adaptation for grapevine (V. vinifera) cultivation. To do so, we evaluated whether CWR accessions have the allelic capacity to persist if moved to locations where grapevines are cultivated in the United States. We identified six candidates from V. mustangensis and hypothesized that they may prove useful for contributing alleles that can mitigate climate impacts on viticulture. By identifying candidate germplasm, this study takes a conceptual step toward assessing the genomic and bioclimatic characteristics of CWRs.
Collapse
Affiliation(s)
- Jonas A Aguirre-Liguori
- Department of Ecology and Evolutionary Biology, University of California, Irvine, California, USA
| | - Abraham Morales-Cruz
- Department of Ecology and Evolutionary Biology, University of California, Irvine, California, USA
| | - Brandon S Gaut
- Department of Ecology and Evolutionary Biology, University of California, Irvine, California, USA
| |
Collapse
|
19
|
Thorpe HA, Tourrette E, Yahara K, Vale FF, Liu S, Oleastro M, Alarcon T, Perets TT, Latifi-Navid S, Yamaoka Y, Martinez-Gonzalez B, Karayiannis I, Karamitros T, Sgouras DN, Elamin W, Pascoe B, Sheppard SK, Ronkainen J, Aro P, Engstrand L, Agreus L, Suerbaum S, Thorell K, Falush D. Repeated out-of-Africa expansions of Helicobacter pylori driven by replacement of deleterious mutations. Nat Commun 2022; 13:6842. [PMID: 36369175 PMCID: PMC9652371 DOI: 10.1038/s41467-022-34475-3] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/05/2022] [Accepted: 10/26/2022] [Indexed: 11/13/2022] Open
Abstract
Helicobacter pylori lives in the human stomach and has a population structure resembling that of its host. However, H. pylori from Europe and the Middle East trace substantially more ancestry from modern African populations than the humans that carry them. Here, we use a collection of Afro-Eurasian H. pylori genomes to show that this African ancestry is due to at least three distinct admixture events. H. pylori from East Asia, which have undergone little admixture, have accumulated many more non-synonymous mutations than African strains. European and Middle Eastern bacteria have elevated African ancestry at the sites of these mutations, implying selection to remove them during admixture. Simulations show that population fitness can be restored after bottlenecks by migration and subsequent admixture of small numbers of bacteria from non-bottlenecked populations. We conclude that recent spread of African DNA has been driven by deleterious mutations accumulated during the original out-of-Africa bottleneck.
Collapse
Affiliation(s)
- Harry A Thorpe
- Department of Biostatistics, University of Oslo, Oslo, Norway
| | - Elise Tourrette
- CAS Key Laboratory of Molecular Virology and Immunology, Institut Pasteur of Shanghai, Chinese Academy of Sciences, Shanghai, China
| | - Koji Yahara
- Antimicrobial Resistance Research Center, National Institute of Infectious Diseases, Tokyo, Japan
| | - Filipa F Vale
- Pathogen Genome Bioinformatics and Computational Biology, Research Institute for Medicines (iMed.ULisboa), Faculty of Pharmacy, Universidade de Lisboa, Lisbon, Portugal
| | - Siqi Liu
- CAS Key Laboratory of Molecular Virology and Immunology, Institut Pasteur of Shanghai, Chinese Academy of Sciences, Shanghai, China
- University of Chinese Academy of Sciences, Beijing, China
| | - Mónica Oleastro
- National Reference Laboratory for Gastrointestinal Infections, Department of Infectious Diseases, National Institute of Health Dr Ricardo Jorge, Lisbon, Portugal
| | - Teresa Alarcon
- Department of Microbiology, Hospital Universitario La Princesa, Instituto de Investigación Sanitaria Princesa, Madrid, Spain
| | - Tsachi-Tsadok Perets
- Gastroenterology Laboratory, Rabin Medical Center, Petah Tikva, Israel
- Department of Digital Medical Technologies, Holon Institute of Technology, Holon, Israel
| | - Saeid Latifi-Navid
- Department of Biology, Faculty of Sciences, University of Mohaghegh Ardabili, Ardabil, Iran
| | - Yoshio Yamaoka
- Department of Environmental and Preventive Medicine, Oita University Faculty of Medicine, Yufu, Oita, Japan
- Department of Medicine-Gastroenterology, Baylor College of Medicine, Houston, TX, USA
| | | | - Ioannis Karayiannis
- Laboratory of Medical Microbiology, Hellenic Pasteur Institute, Athens, Greece
| | | | | | - Wael Elamin
- G42 Healthcare, Abu Dhabi, UAE
- Elrazi University, Khartoum, Sudan
| | - Ben Pascoe
- Department of Biology, University of Oxford, Oxford, UK
| | - Samuel K Sheppard
- Ineos Oxford Institute, Department of Biology, University of Oxford, Oxford, UK
| | - Jukka Ronkainen
- Center for Life Course Health Research, University of Oulu, Oulu, Finland
- Primary Health Care Center, Tornio, Finland
| | | | - Lars Engstrand
- Center for Translational Microbiome Research, Department for Microbiology, Tumor, and Cell Biology, Karolinska Institutet, Stockholm, Sweden
| | - Lars Agreus
- Division of Family Medicine, Karolinska Institutet, Stockholm, Sweden
| | - Sebastian Suerbaum
- Department of Medical Microbiology and Hospital Epidemiology, Max von Pettenkofer Institute, Faculty of Medicine, LMU Munich, Munich, Germany
- Department of Medical Microbiology and Hospital Epidemiology, Hannover Medical School, Hanover, Germany
- DZIF German Center for Infection Research, Hannover-Braunschweig and Munich Partner Sites, Munich, Germany
| | - Kaisa Thorell
- Institute of Biomedicine, Department of Infectious Diseases, University of Gothenburg, Gothenburg, Sweden
- Department of Clinical Microbiology, Sahlgrenska University Hospital, Gothenburg, Sweden
| | - Daniel Falush
- CAS Key Laboratory of Molecular Virology and Immunology, Institut Pasteur of Shanghai, Chinese Academy of Sciences, Shanghai, China.
| |
Collapse
|
20
|
Halabian R, Makałowski W. A Map of 3' DNA Transduction Variants Mediated by Non-LTR Retroelements on 3202 Human Genomes. BIOLOGY 2022; 11:1032. [PMID: 36101413 PMCID: PMC9311842 DOI: 10.3390/biology11071032] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/07/2022] [Revised: 07/05/2022] [Accepted: 07/06/2022] [Indexed: 05/03/2023]
Abstract
As one of the major structural constituents, mobile elements comprise more than half of the human genome, among which Alu, L1, and SVA elements are still active and continue to generate new offspring. One of the major characteristics of L1 and SVA elements is their ability to co-mobilize adjacent downstream sequences to new loci in a process called 3' DNA transduction. Transductions influence the structure and content of the genome in different ways, such as increasing genome variation, exon shuffling, and gene duplication. Moreover, given their mutagenicity capability, 3' transductions are often involved in tumorigenesis or in the development of some diseases. In this study, we analyzed 3202 genomes sequenced at high coverage by the New York Genome Center to catalog and characterize putative 3' transduced segments mediated by L1s and SVAs. Here, we present a genome-wide map of inter/intrachromosomal 3' transduction variants, including their genomic and functional location, length, progenitor location, and allelic frequency across 26 populations. In total, we identified 7103 polymorphic L1s and 3040 polymorphic SVAs. Of these, 268 and 162 variants were annotated as high-confidence L1 and SVA 3' transductions, respectively, with lengths that ranged from 7 to 997 nucleotides. We found specific loci within chromosomes X, 6, 7, and 6_GL000253v2_alt as master L1s and SVAs that had yielded more transductions, among others. Together, our results demonstrate the dynamic nature of transduction events within the genome and among individuals and their contribution to the structural variations of the human genome.
Collapse
Affiliation(s)
| | - Wojciech Makałowski
- Institute of Bioinformatics, Faculty of Medicine, University of Münster, 48149 Münster, Germany;
| |
Collapse
|
21
|
Perrier A, Sánchez‐Castro D, Willi Y. Environment dependence of the expression of mutational load and species' range limits. J Evol Biol 2022; 35:731-741. [PMID: 35290676 PMCID: PMC9314787 DOI: 10.1111/jeb.13997] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2021] [Revised: 02/10/2022] [Accepted: 02/21/2022] [Indexed: 12/21/2022]
Abstract
Theoretical and empirical research on the causes of species' range limits suggest the contribution of several intrinsic and extrinsic factors, with potentially complex interactions among them. An intrinsic factor proposed by recent theory is mutational load increasing towards range edges because of genetic drift. Furthermore, environmental quality may decline towards range edges and enhance the expression of load. Here, we tested whether the expression of mutational load associated with range limits in the North American plant Arabidopsis lyrata was enhanced under stressful environmental conditions by comparing the performance of within- versus between-population crosses at common garden sites across the species' distribution and beyond. Heterosis, reflecting the expression of load, increased with heightened estimates of genomic load and with environmental stress caused by warming, but the interaction was not significant. We conclude that range-edge populations suffer from a twofold genetic Allee effect caused by increased mutational load and stress-dependent load linked to general heterozygote deficiency, but there is no synergistic effect between them.
Collapse
Affiliation(s)
- Antoine Perrier
- Department of BiologyUniversity of VirginiaCharlottesvilleVirginiaUSA
- Department of Environmental SciencesUniversity of BaselBaselSwitzerland
| | | | - Yvonne Willi
- Department of Environmental SciencesUniversity of BaselBaselSwitzerland
| |
Collapse
|
22
|
Niedbalski SD, Long JC. Novel alleles gained during the Beringian isolation period. Sci Rep 2022; 12:4289. [PMID: 35277570 PMCID: PMC8917172 DOI: 10.1038/s41598-022-08212-1] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/27/2021] [Accepted: 02/09/2022] [Indexed: 01/23/2023] Open
Abstract
During the Last Glacial Maximum, a small band of Siberians entered the Beringian corridor, where they persisted, isolated from gene flow, for several thousand years before expansion into the Americas. The ecological features of the Beringian environment, coupled with an extended period of isolation at small population size, would have provided evolutionary opportunity for novel genetic variation to arise as both rare standing variants and new mutations were driven to high frequency through both neutral and directed processes. Here we perform a full genome investigation of Native American populations in the Thousand Genomes Project Phase 3 to identify unique high frequency alleles that can be dated to an origin in Beringia. Our analyses demonstrate that descendant populations of Native Americans harbor 20,424 such variants, which is on a scale comparable only to Africa and the Out of Africa bottleneck. This is consistent with simulations of a serial founder effects model. Tests for selection reveal that some of these Beringian variants were likely driven to high frequency by adaptive processes, and bioinformatic analyses suggest possible phenotypic pathways that were under selection during the Beringian Isolation period. Specifically, pathways related to cardiac processes and melanocyte function appear to be enriched for selected Beringian variants.
Collapse
Affiliation(s)
- Sara D Niedbalski
- Human Evolutionary Genetics Unit, UMR 2000, CNRS, Institut Pasteur, Paris, France.,Department of Anthropology, University of New Mexico, Albuquerque, NM, 87131, USA
| | - Jeffrey C Long
- Department of Anthropology, University of New Mexico, Albuquerque, NM, 87131, USA.
| |
Collapse
|
23
|
Pediatric neurological cancer incidence and trends in the United States, 2000–2018. Cancer Causes Control 2022; 33:687-699. [DOI: 10.1007/s10552-021-01535-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2021] [Accepted: 12/01/2021] [Indexed: 10/19/2022]
|
24
|
The Myth of the Genetically Sick African. GENEALOGY 2022. [DOI: 10.3390/genealogy6010015] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/27/2023]
Abstract
Western medicine has an unfortunate history where it has been applied to address the health of African Americans. At its origins, it was aligned with the objectives of colonialism and chattel slavery. The degree to which medical “science” concerned itself with persons of African descent was to keep them alive for sale on the auction block, or to keep them healthy as they toiled to generate wealth for their European owners. Medicine in early America relied upon both dead and live African bodies to test its ideas to benefit Europeans. As medicine moved from quackery to a discipline based in science, its understanding of human biological variation was flawed. This was not a problem confined to medicine alone, but to the biological sciences in general. Biology had no solid theoretical basis until after 1859. As medicine further developed in the 20th century, it never doubted the difference between Europeans and Africans, and also asserted the innate inferiority of the latter. The genomic revolution in the latter 20th century produced tools that were deployed in a biomedical culture still mired in “racial” medicine. This lack of theoretical perspective still misdirects research associated with health disparity. In contrast to this is evolutionary medicine, which relies on a sound unification of evolutionary (ultimate) and physiological, cellular, and molecular (proximate) mechanisms. Utilizing the perspectives of evolutionary medicine is a prerequisite for an effective intervention in health disparity and finally dispelling the myth of the genetically sick African.
Collapse
|
25
|
Conover JL, Wendel JF. Deleterious Mutations Accumulate Faster in Allopolyploid than Diploid Cotton (Gossypium) and Unequally between Subgenomes. Mol Biol Evol 2022; 39:6517786. [PMID: 35099532 PMCID: PMC8841602 DOI: 10.1093/molbev/msac024] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022] Open
Abstract
Abstract
Whole genome duplication (polyploidization) is among the most dramatic mutational processes in nature, so understanding how natural selection differs in polyploids relative to diploids is an important goal. Population genetics theory predicts that recessive deleterious mutations accumulate faster in allopolyploids than diploids due to the masking effect of redundant gene copies, but this prediction is hitherto unconfirmed. Here, we use the cotton genus (Gossypium), which contains seven allopolyploids derived from a single polyploidization event 1-2 million years ago, to investigate deleterious mutation accumulation. We use two methods of identifying deleterious mutations at the nucleotide and amino acid level, along with whole-genome resequencing of 43 individuals spanning six allopolyploid species and their two diploid progenitors, to demonstrate that deleterious mutations accumulate faster in allopolyploids than in their diploid progenitors. We find that, unlike what would be expected under models of demographic changes alone, strongly deleterious mutations show the biggest difference between ploidy levels, and this effect diminishes for moderately and mildly deleterious mutations. We further show that the proportion of nonsynonymous mutations that are deleterious differs between the two co-resident subgenomes in the allopolyploids, suggesting that homoeologous masking acts unequally between subgenomes. Our results provide a genome-wide perspective on classic notions of the significance of gene duplication that likely are broadly applicable to allopolyploids, with implications for our understanding of the evolutionary fate of deleterious mutations. Finally, we note that some measures of selection (e.g. dN/dS, πN/πS) may be biased when species of different ploidy levels are compared.
Collapse
Affiliation(s)
- Justin L Conover
- Department of Ecology, Evolution, and Organismal Biology, Iowa State University, Ames, IA, 50011, USA
| | - Jonathan F Wendel
- Department of Ecology, Evolution, and Organismal Biology, Iowa State University, Ames, IA, 50011, USA
| |
Collapse
|
26
|
Stover DA, Housman G, Stone AC, Rosenberg MS, Verrelli BC. Evolutionary Genetic Signatures of Selection on Bone-Related Variation within Human and Chimpanzee Populations. Genes (Basel) 2022; 13:183. [PMID: 35205228 PMCID: PMC8871609 DOI: 10.3390/genes13020183] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2021] [Revised: 01/19/2022] [Accepted: 01/19/2022] [Indexed: 02/06/2023] Open
Abstract
Bone strength and the incidence and severity of skeletal disorders vary significantly among human populations, due in part to underlying genetic differentiation. While clinical models predict that this variation is largely deleterious, natural population variation unrelated to disease can go unnoticed, altering our perception of how natural selection has shaped bone morphologies over deep and recent time periods. Here, we conduct the first comparative population-based genetic analysis of the main bone structural protein gene, collagen type I α 1 (COL1A1), in clinical and 1000 Genomes Project datasets in humans, and in natural populations of chimpanzees. Contrary to predictions from clinical studies, we reveal abundant COL1A1 amino acid variation, predicted to have little association with disease in the natural population. We also find signatures of positive selection associated with intron haplotype structure, linkage disequilibrium, and population differentiation in regions of known gene expression regulation in humans and chimpanzees. These results recall how recent and deep evolutionary regimes can be linked, in that bone morphology differences that developed among vertebrates over 450 million years of evolution are the result of positive selection on subtle type I collagen functional variation segregating within populations over time.
Collapse
Affiliation(s)
- Daryn A. Stover
- School of Life Sciences, Arizona State University, Tempe, AZ 85287, USA;
- Arizona State University at Lake Havasu, Lake Havasu, AZ 86403, USA
| | - Genevieve Housman
- Section of Genetic Medicine, University of Chicago, Chicago, IL 60637, USA;
| | - Anne C. Stone
- School of Human Evolution and Social Change, Arizona State University, Tempe, AZ 85287, USA;
| | - Michael S. Rosenberg
- Center for Biological Data Science, Virginia Commonwealth University, Richmond, VA 23284, USA;
| | - Brian C. Verrelli
- School of Life Sciences, Arizona State University, Tempe, AZ 85287, USA;
- Center for Biological Data Science, Virginia Commonwealth University, Richmond, VA 23284, USA;
| |
Collapse
|
27
|
Soni V, Eyre-Walker A. OUP accepted manuscript. Genome Biol Evol 2022; 14:6528851. [PMID: 35166775 PMCID: PMC8882387 DOI: 10.1093/gbe/evac028] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 02/09/2022] [Indexed: 12/05/2022] Open
Abstract
The rate of amino acid substitution has been shown to be correlated to a number of factors including the rate of recombination, the age of the gene, the length of the protein, mean expression level, and gene function. However, the extent to which these correlations are due to adaptive and nonadaptive evolution has not been studied in detail, at least not in hominids. We find that the rate of adaptive evolution is significantly positively correlated to the rate of recombination, protein length and gene expression level, and negatively correlated to gene age. These correlations remain significant when each factor is controlled for in turn, except when controlling for expression in an analysis of protein length; and they also generally remain significant when biased gene conversion is taken into account. However, the positive correlations could be an artifact of population size contraction. We also find that the rate of nonadaptive evolution is negatively correlated to each factor, and all these correlations survive controlling for each other and biased gene conversion. Finally, we examine the effect of gene function on rates of adaptive and nonadaptive evolution; we confirm that virus-interacting proteins (VIPs) have higher rates of adaptive and lower rates of nonadaptive evolution, but we also demonstrate that there is significant variation in the rate of adaptive and nonadaptive evolution between GO categories when removing VIPs. We estimate that the VIP/non-VIP axis explains about 5–8 fold more of the variance in evolutionary rate than GO categories.
Collapse
Affiliation(s)
- Vivak Soni
- School of Life Sciences, University of Sussex, Brighton, United Kingdom
| | - Adam Eyre-Walker
- School of Life Sciences, University of Sussex, Brighton, United Kingdom
- Corresponding author: E-mail:
| |
Collapse
|
28
|
Ackiss AS, Magee MR, Sass GG, Turnquist K, McIntyre PB, Larson WA. Genomic and environmental influences on resilience in a cold-water fish near the edge of its range. Evol Appl 2021; 14:2794-2814. [PMID: 34950230 PMCID: PMC8674893 DOI: 10.1111/eva.13313] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2021] [Revised: 10/05/2021] [Accepted: 10/07/2021] [Indexed: 02/06/2023] Open
Abstract
Small, isolated populations present a challenge for conservation. The dueling effects of selection and drift in a limited pool of genetic diversity make the responses of small populations to environmental perturbations erratic and difficult to predict. This is particularly true at the edge of a species range, where populations often persist at the limits of their environmental tolerances. Populations of cisco, Coregonus artedi, in inland lakes have experienced numerous extirpations along the southern edge of their range in recent decades, which are thought to result from environmental degradation and loss of cold, well-oxygenated habitat as lakes warm. Yet, cisco extirpations do not show a clear latitudinal pattern, suggesting that local environmental factors and potentially local adaptation may influence resilience. Here, we used genomic tools to investigate the nature of this pattern of resilience. We used restriction site-associated DNA capture (Rapture) sequencing to survey genomic diversity and differentiation in southern inland lake cisco populations and compared the frequency of deleterious mutations that potentially influence fitness across lakes. We also examined haplotype diversity in a region of the major histocompatibility complex involved in stress and immune system response. We correlated these metrics to spatial and environmental factors including latitude, lake size, and measures of oxythermal habitat and found significant relationships between genetic metrics and broad and local factors. High levels of genetic differentiation among populations were punctuated by a phylogeographic break and residual patterns of isolation-by-distance. Although the prevalence of deleterious mutations and inbreeding coefficients was significantly correlated with latitude, neutral and non-neutral genetic diversity were most strongly correlated with lake surface area. Notably, differences among lakes in the availability of estimated oxythermal habitat left no clear population genomic signature. Our results shed light on the complex dynamics influencing these isolated populations and provide valuable information for their conservation.
Collapse
Affiliation(s)
- Amanda S. Ackiss
- Wisconsin Cooperative Fishery Research UnitCollege of Natural ResourcesUniversity of Wisconsin‐Stevens PointStevens PointWisconsinUSA
- U.S. Geological SurveyGreat Lakes Science CenterAnn ArborMichiganUSA
| | | | - Greg G. Sass
- Escanaba Lake Research StationWisconsin Department of Natural ResourcesBoulder JunctionWisconsinUSA
| | - Keith Turnquist
- Wisconsin Cooperative Fishery Research UnitCollege of Natural ResourcesUniversity of Wisconsin‐Stevens PointStevens PointWisconsinUSA
| | - Peter B. McIntyre
- Department of Natural Resources and the EnvironmentCornell UniversityIthacaNew YorkUSA
| | - Wesley A. Larson
- U.S. Geological SurveyWisconsin Cooperative Fishery Research UnitCollege of Natural ResourcesUniversity of Wisconsin‐Stevens PointStevens PointWisconsinUSA
- National Oceanographic and Atmospheric AdministrationNational Marine Fisheries ServiceAlaska Fisheries Science CenterAuke Bay LaboratoriesJuneauAlaskaUSA
| |
Collapse
|
29
|
Whole-exome analysis in Tunisian Imazighen and Arabs shows the impact of demography in functional variation. Sci Rep 2021; 11:21125. [PMID: 34702931 PMCID: PMC8548440 DOI: 10.1038/s41598-021-00576-0] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2021] [Accepted: 10/14/2021] [Indexed: 11/08/2022] Open
Abstract
Human populations are genetically affected by their demographic history, which shapes the distribution of their functional genomic variation. However, the genetic impact of recent demography is debated. This issue has been studied in different populations, but never in North Africans, despite their relevant cultural and demographic diversity. In this study we address the question by analyzing new whole-exome sequences from two culturally different Tunisian populations, an isolated Amazigh population and a close non-isolated Arab-speaking population, focusing on the distribution of functional variation. Both populations present clear differences in their variant frequency distribution, in general and for putatively damaging variation. This suggests a relevant effect in the Amazigh population of genetic isolation, drift, and inbreeding, pointing to relaxed purifying selection. We also discover the enrichment in Imazighen of variation associated to specific diseases or phenotypic traits, but the scarce genetic and biomedical data in the region limits further interpretation. Our results show the genomic impact of recent demography and reveal a clear genetic differentiation probably related to culture. These findings highlight the importance of considering cultural and demographic heterogeneity within North Africa when defining population groups, and the need for more data to improve knowledge on the region's health and disease landscape.
Collapse
|
30
|
Svensson E, Günther T, Hoischen A, Hervella M, Munters AR, Ioana M, Ridiche F, Edlund H, van Deuren RC, Soficaru A, de-la-Rua C, Netea MG, Jakobsson M. Genome of Peştera Muierii skull shows high diversity and low mutational load in pre-glacial Europe. Curr Biol 2021; 31:2973-2983.e9. [PMID: 34010592 DOI: 10.1016/j.cub.2021.04.045] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2020] [Revised: 04/09/2021] [Accepted: 04/19/2021] [Indexed: 11/24/2022]
Abstract
Few complete human genomes from the European Early Upper Palaeolithic (EUP) have been sequenced. Using novel sampling and DNA extraction approaches, we sequenced the genome of a woman from "Peştera Muierii," Romania who lived ∼34,000 years ago to 13.5× coverage. The genome shows similarities to modern-day Europeans, but she is not a direct ancestor. Although her cranium exhibits both modern human and Neanderthal features, the genome shows similar levels of Neanderthal admixture (∼3.1%) to most EUP humans but only half compared to the ∼40,000-year-old Peştera Oase 1. All EUP European hunter-gatherers display high genetic diversity, demonstrating that the severe loss of diversity occurred during and after the Last Glacial Maximum (LGM) rather than just during the out-of-Africa migration. The prevalence of genetic diseases is expected to increase with low diversity; however, pathogenic variant load was relatively constant from EUP to modern times, despite post-LGM hunter-gatherers having the lowest diversity ever observed among Europeans.
Collapse
Affiliation(s)
- Emma Svensson
- Human Evolution, Department of Organismal Biology, Uppsala University, 752 36 Uppsala, Sweden
| | - Torsten Günther
- Human Evolution, Department of Organismal Biology, Uppsala University, 752 36 Uppsala, Sweden.
| | - Alexander Hoischen
- Department of Human Genetics, Radboud University Medical Center, 6526 Nijmegen, the Netherlands; Department of Internal Medicine and Radboud Center for Infectious Diseases (RCI), Radboud University Medical Center, 6526 Nijmegen, the Netherlands
| | - Montserrat Hervella
- Department of Genetics, Physical Anthropology and Animal Physiology, Faculty of Science and Technology, University of the Basque Country (UPV/EHU), B° Sarriena s/n 48940 Leioa, Bizkaia, Spain
| | - Arielle R Munters
- Human Evolution, Department of Organismal Biology, Uppsala University, 752 36 Uppsala, Sweden
| | - Mihai Ioana
- Laboratory of Human Genetics, University of Medicine and Pharmacy, Craiova, Romania
| | | | - Hanna Edlund
- Human Evolution, Department of Organismal Biology, Uppsala University, 752 36 Uppsala, Sweden
| | - Rosanne C van Deuren
- Department of Internal Medicine and Radboud Center for Infectious Diseases (RCI), Radboud University Medical Center, 6526 Nijmegen, the Netherlands
| | - Andrei Soficaru
- "Francisc J. Rainer" Institute of Anthropology, Romanian Academy, 050474 Bucharest, Romania
| | - Concepción de-la-Rua
- Department of Genetics, Physical Anthropology and Animal Physiology, Faculty of Science and Technology, University of the Basque Country (UPV/EHU), B° Sarriena s/n 48940 Leioa, Bizkaia, Spain
| | - Mihai G Netea
- Department of Internal Medicine and Radboud Center for Infectious Diseases (RCI), Radboud University Medical Center, 6526 Nijmegen, the Netherlands; Laboratory of Human Genetics, University of Medicine and Pharmacy, Craiova, Romania
| | - Mattias Jakobsson
- Human Evolution, Department of Organismal Biology, Uppsala University, 752 36 Uppsala, Sweden.
| |
Collapse
|
31
|
Font-Porterias N, Caro-Consuegra R, Lucas-Sánchez M, Lopez M, Giménez A, Carballo-Mesa A, Bosch E, Calafell F, Quintana-Murci L, Comas D. The Counteracting Effects of Demography on Functional Genomic Variation: The Roma Paradigm. Mol Biol Evol 2021; 38:2804-2817. [PMID: 33713133 PMCID: PMC8233508 DOI: 10.1093/molbev/msab070] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Abstract
Demographic history plays a major role in shaping the distribution of genomic variation. Yet the interaction between different demographic forces and their effects in the genomes is not fully resolved in human populations. Here, we focus on the Roma population, the largest transnational ethnic minority in Europe. They have a South Asian origin and their demographic history is characterized by recent dispersals, multiple founder events, and extensive gene flow from non-Roma groups. Through the analyses of new high-coverage whole exome sequences and genome-wide array data for 89 Iberian Roma individuals together with forward simulations, we show that founder effects have reduced their genetic diversity and proportion of rare variants, gene flow has counteracted the increase in mutational load, runs of homozygosity show ancestry-specific patterns of accumulation of deleterious homozygotes, and selection signals primarily derive from preadmixture adaptation in the Roma population sources. The present study shows how two demographic forces, bottlenecks and admixture, act in opposite directions and have long-term balancing effects on the Roma genomes. Understanding how demography and gene flow shape the genome of an admixed population provides an opportunity to elucidate how genomic variation is modeled in human populations.
Collapse
Affiliation(s)
- Neus Font-Porterias
- Departament de Ciències Experimentals i de la Salut, Institut de Biologia Evolutiva (UPF-CSIC), Universitat Pompeu Fabra, Barcelona, Spain
| | - Rocio Caro-Consuegra
- Departament de Ciències Experimentals i de la Salut, Institut de Biologia Evolutiva (UPF-CSIC), Universitat Pompeu Fabra, Barcelona, Spain
| | - Marcel Lucas-Sánchez
- Departament de Ciències Experimentals i de la Salut, Institut de Biologia Evolutiva (UPF-CSIC), Universitat Pompeu Fabra, Barcelona, Spain
| | - Marie Lopez
- Human Evolutionary Genetics Unit, Institut Pasteur, UMR2000, CNRS, Paris, France
| | - Aaron Giménez
- Facultat de Sociologia, Universitat Autònoma de Barcelona, Barcelona, Spain
| | | | - Elena Bosch
- Departament de Ciències Experimentals i de la Salut, Institut de Biologia Evolutiva (UPF-CSIC), Universitat Pompeu Fabra, Barcelona, Spain.,Centro de Investigación Biomédica en Red de Salud Mental (CIBERSAM), Reus, Spain
| | - Francesc Calafell
- Departament de Ciències Experimentals i de la Salut, Institut de Biologia Evolutiva (UPF-CSIC), Universitat Pompeu Fabra, Barcelona, Spain
| | - Lluís Quintana-Murci
- Human Evolutionary Genetics Unit, Institut Pasteur, UMR2000, CNRS, Paris, France.,Human Genomics and Evolution, Collège de France, Paris, France
| | - David Comas
- Departament de Ciències Experimentals i de la Salut, Institut de Biologia Evolutiva (UPF-CSIC), Universitat Pompeu Fabra, Barcelona, Spain
| |
Collapse
|
32
|
Cuna A, Sampath V, Khashu M. Racial Disparities in Necrotizing Enterocolitis. Front Pediatr 2021; 9:633088. [PMID: 33681105 PMCID: PMC7930220 DOI: 10.3389/fped.2021.633088] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 11/24/2020] [Accepted: 01/28/2021] [Indexed: 12/20/2022] Open
Abstract
Necrotizing enterocolitis (NEC) is a serious disease of the intestinal tract affecting 5-10% of pre-term infants with up to 50% mortality in those that require surgery. There is wide variation in the rates and outcomes of NEC by race and ethnicity, and the reasons for this disparity are poorly understood. In this article, we review the epidemiology and discuss possible explanations for racial and ethnic differences in NEC. Most of the current evidence investigating the role of race in NEC comes from North America and suggests that Hispanic ethnicity and non-Hispanic Black race are associated with higher risk of NEC compared to non-Hispanic White populations. Differences in pre-term births, breastfeeding rates, and various sociodemographic factors does not fully account for the observed disparities in NEC incidence and outcomes. While genetic studies are beginning to identify candidate genes that may increase or decrease risk for NEC among racial populations, current data remain limited by small sample sizes and lack of validation. Complex interactions between social and biological determinants likely underly the differences in NEC outcomes among racial groups. Larger datasets with detailed social, phenotypic, and genotypic information, coupled with advanced bioinformatics techniques are needed to comprehensively understand racial disparities in NEC.
Collapse
Affiliation(s)
- Alain Cuna
- School of Medicine, University of Missouri Kansas City, Kansas City, MO, United States.,Division of Neonatology, Children's Mercy Kansas City, Kansas City, MO, United States
| | - Venkatesh Sampath
- School of Medicine, University of Missouri Kansas City, Kansas City, MO, United States.,Division of Neonatology, Children's Mercy Kansas City, Kansas City, MO, United States
| | - Minesh Khashu
- Neonatal Service, University Hospitals Dorset, Poole, United Kingdom.,Bournemouth University, Dorset, United Kingdom
| |
Collapse
|
33
|
Vy HMT, Jordan DM, Balick DJ, Do R. Probing the aggregated effects of purifying selection per individual on 1,380 medical phenotypes in the UK Biobank. PLoS Genet 2021; 17:e1009337. [PMID: 33493176 PMCID: PMC7861521 DOI: 10.1371/journal.pgen.1009337] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2020] [Revised: 02/04/2021] [Accepted: 01/04/2021] [Indexed: 11/26/2022] Open
Abstract
Understanding the relationship between natural selection and phenotypic variation has been a long-standing challenge in human population genetics. With the emergence of biobank-scale datasets, along with new statistical metrics to approximate strength of purifying selection at the variant level, it is now possible to correlate a proxy of individual relative fitness with a range of medical phenotypes. We calculated a per-individual deleterious load score by summing the total number of derived alleles per individual after incorporating a weight that approximates strength of purifying selection. We assessed four methods for the weight, including GERP, phyloP, CADD, and fitcons. By quantitatively tracking each of these scores with the site frequency spectrum, we identified phyloP as the most appropriate weight. The phyloP-weighted load score was then calculated across 15,129,142 variants in 335,161 individuals from the UK Biobank and tested for association on 1,380 medical phenotypes. After accounting for multiple test correction, we observed a strong association of the load score amongst coding sites only on 27 traits including body mass, adiposity and metabolic rate. We further observed that the association signals were driven by common variants (derived allele frequency > 5%) with high phyloP score (phyloP > 2). Finally, through permutation analyses, we showed that the load score amongst coding sites had an excess of nominally significant associations on many medical phenotypes. These results suggest a broad impact of deleterious load on medical phenotypes and highlight the deleterious load score as a tool to disentangle the complex relationship between natural selection and medical phenotypes.
Collapse
Affiliation(s)
- Ha My T. Vy
- The Charles Bronfman Institute for Personalized Medicine, Icahn School of Medicine at Mount Sinai, New York, New York, United States of America
- Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York, New York, United States of America
| | - Daniel M. Jordan
- The Charles Bronfman Institute for Personalized Medicine, Icahn School of Medicine at Mount Sinai, New York, New York, United States of America
- Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York, New York, United States of America
| | - Daniel J. Balick
- Division of Genetics, Brigham and Women’s Hospital, Harvard Medical School, Boston, Massachusetts, United States of America
- Department of Biomedical Informatics, Harvard Medical School, Boston, Massachusetts, United States of America
| | - Ron Do
- The Charles Bronfman Institute for Personalized Medicine, Icahn School of Medicine at Mount Sinai, New York, New York, United States of America
- Department of Genetics and Genomic Sciences, Icahn School of Medicine at Mount Sinai, New York, New York, United States of America
| |
Collapse
|
34
|
Pejaver V, Urresti J, Lugo-Martinez J, Pagel KA, Lin GN, Nam HJ, Mort M, Cooper DN, Sebat J, Iakoucheva LM, Mooney SD, Radivojac P. Inferring the molecular and phenotypic impact of amino acid variants with MutPred2. Nat Commun 2020; 11:5918. [PMID: 33219223 PMCID: PMC7680112 DOI: 10.1038/s41467-020-19669-x] [Citation(s) in RCA: 291] [Impact Index Per Article: 72.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2019] [Accepted: 10/23/2020] [Indexed: 01/02/2023] Open
Abstract
Identifying pathogenic variants and underlying functional alterations is challenging. To this end, we introduce MutPred2, a tool that improves the prioritization of pathogenic amino acid substitutions over existing methods, generates molecular mechanisms potentially causative of disease, and returns interpretable pathogenicity score distributions on individual genomes. Whilst its prioritization performance is state-of-the-art, a distinguishing feature of MutPred2 is the probabilistic modeling of variant impact on specific aspects of protein structure and function that can serve to guide experimental studies of phenotype-altering variants. We demonstrate the utility of MutPred2 in the identification of the structural and functional mutational signatures relevant to Mendelian disorders and the prioritization of de novo mutations associated with complex neurodevelopmental disorders. We then experimentally validate the functional impact of several variants identified in patients with such disorders. We argue that mechanism-driven studies of human inherited disease have the potential to significantly accelerate the discovery of clinically actionable variants.
Collapse
Affiliation(s)
- Vikas Pejaver
- Department of Computer Science, Indiana University, Bloomington, IN, USA
- Department of Biomedical Informatics and Medical Education, University of Washington, Seattle, WA, USA
| | - Jorge Urresti
- Department of Psychiatry, University of California San Diego, La Jolla, CA, USA
| | - Jose Lugo-Martinez
- Department of Computer Science, Indiana University, Bloomington, IN, USA
- Computational Biology Department, School of Computer Science, Carnegie Mellon University, 5000 Forbes Avenue, Pittsburgh, PA, 15213, USA
| | - Kymberleigh A Pagel
- Department of Computer Science, Indiana University, Bloomington, IN, USA
- Institute for Computational Medicine, Whiting School of Engineering, Johns Hopkins University, 220 Hackerman Hall, 3400 N Charles St, Baltimore, MD, 21218, USA
| | - Guan Ning Lin
- Department of Psychiatry, University of California San Diego, La Jolla, CA, USA
- School of Biomedical Engineering, Shanghai Jiao Tong University, Shanghai, 200030, People's Republic of China
| | - Hyun-Jun Nam
- Department of Psychiatry, University of California San Diego, La Jolla, CA, USA
| | - Matthew Mort
- Institute of Medical Genetics, School of Medicine, Cardiff University, Cardiff, UK
| | - David N Cooper
- Institute of Medical Genetics, School of Medicine, Cardiff University, Cardiff, UK
| | - Jonathan Sebat
- Department of Psychiatry, University of California San Diego, La Jolla, CA, USA
- Beyster Center for Genomics of Psychiatric Diseases, University of California San Diego, La Jolla, CA, USA
- Department of Cellular and Molecular Medicine, University of California San Diego, La Jolla, CA, USA
| | - Lilia M Iakoucheva
- Department of Psychiatry, University of California San Diego, La Jolla, CA, USA.
| | - Sean D Mooney
- Department of Biomedical Informatics and Medical Education, University of Washington, Seattle, WA, USA.
| | - Predrag Radivojac
- Department of Computer Science, Indiana University, Bloomington, IN, USA.
- Khoury College of Computer Sciences, Northeastern University, Boston, MA, USA.
| |
Collapse
|
35
|
Toffano AA, Chiarot G, Zamuner S, Marchi M, Salvi E, Waxman SG, Faber CG, Lauria G, Giacometti A, Simeoni M. Computational pipeline to probe NaV1.7 gain-of-function variants in neuropathic painful syndromes. Sci Rep 2020; 10:17930. [PMID: 33087732 PMCID: PMC7578092 DOI: 10.1038/s41598-020-74591-y] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2020] [Accepted: 09/30/2020] [Indexed: 01/09/2023] Open
Abstract
Applications of machine learning and graph theory techniques to neuroscience have witnessed an increased interest in the last decade due to the large data availability and unprecedented technology developments. Their employment to investigate the effect of mutational changes in genes encoding for proteins modulating the membrane of excitable cells, whose biological correlates are assessed at electrophysiological level, could provide useful predictive clues. We apply this concept to the analysis of variants in sodium channel NaV1.7 subunit found in patients with chronic painful syndromes, by the implementation of a dedicated computational pipeline empowering different and complementary techniques including homology modeling, network theory, and machine learning. By testing three templates of different origin and sequence identities, we provide an optimal condition for its use. Our findings reveal the usefulness of our computational pipeline in supporting the selection of candidates for cell electrophysiology assay and with potential clinical applications.
Collapse
Affiliation(s)
- Alberto A Toffano
- Dipartimento di Scienze Molecolari e Nanosistemi, Universitá Ca' Foscari Venezia, Venezia-Mestre, Italy
| | - Giacomo Chiarot
- Dipartimento di Scienze Ambientali, Informatica e Statistica, Universitá Ca' Foscari Venezia, Venezia-Mestre, Italy
| | - Stefano Zamuner
- Laboratory of Statistical Biophysics, Institute of Physics, School of Basic Sciences, Ècole Polytechnique Fèdèrale de Lausanne (EPFL), Lausanne, Switzerland
| | - Margherita Marchi
- Neuroalgology Unit, Fondazione IRCCS Istituto Neurologico "Carlo Besta", Milan, Italy
| | - Erika Salvi
- Neuroalgology Unit, Fondazione IRCCS Istituto Neurologico "Carlo Besta", Milan, Italy
| | - Stephen G Waxman
- Center for Neuroscience and Regeneration Research, VA Connecticut Healthcare System and Yale Medical School, West Haven, USA
| | - Catharina G Faber
- MHeNs school for Mental Health and Neuroscience, Maastricht University, Maastricht, The Netherlands.,Department of Neurology, Maastricht University Medical Center, Maastricht, The Netherlands
| | - Giuseppe Lauria
- Neuroalgology Unit, Fondazione IRCCS Istituto Neurologico "Carlo Besta", Milan, Italy.,Department of Biomedical and Clinical Sciences "Luigi Sacco", University of Milan, Milan, Italy
| | - Achille Giacometti
- Dipartimento di Scienze Molecolari e Nanosistemi, Universitá Ca' Foscari Venezia, Venezia-Mestre, Italy.,European Centre for Living Technology (ECLT), Venice, Italy
| | - Marta Simeoni
- Dipartimento di Scienze Ambientali, Informatica e Statistica, Universitá Ca' Foscari Venezia, Venezia-Mestre, Italy. .,European Centre for Living Technology (ECLT), Venice, Italy.
| |
Collapse
|
36
|
Rinker DC, Simonti CN, McArthur E, Shaw D, Hodges E, Capra JA. Neanderthal introgression reintroduced functional ancestral alleles lost in Eurasian populations. Nat Ecol Evol 2020; 4:1332-1341. [PMID: 32719451 PMCID: PMC7529911 DOI: 10.1038/s41559-020-1261-z] [Citation(s) in RCA: 21] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/26/2019] [Accepted: 06/29/2020] [Indexed: 12/31/2022]
Abstract
Neanderthal ancestry remains across modern Eurasian genomes and introgressed sequences influence diverse phenotypes. Here, we demonstrate that introgressed sequences reintroduced thousands of ancestral alleles that were lost in Eurasian populations before introgression. Our simulations and variant effect predictions argue that these reintroduced alleles (RAs) are more likely to be tolerated by modern humans than are introgressed Neanderthal-derived alleles (NDAs) due to their distinct evolutionary histories. Consistent with this, we show enrichment for RAs and depletion for NDAs on introgressed haplotypes with expression quantitative trait loci (eQTL) and phenotype associations. Analysis of available cross-population eQTLs and massively parallel reporter assay data show that RAs commonly influence gene expression independent of linked NDAs. We further validate these independent effects for one RA in vitro. Finally, we demonstrate that NDAs are depleted for regulatory activity compared to RAs, while RAs have activity levels similar to non-introgressed variants. In summary, our study reveals that Neanderthal introgression reintroduced thousands of lost ancestral variants with gene regulatory activity and that these RAs were more tolerated than NDAs. Thus, RAs and their distinct evolutionary histories must be considered when evaluating the effects of introgression.
Collapse
Affiliation(s)
- David C Rinker
- Department of Biological Sciences, Vanderbilt University, Nashville, TN, USA
| | - Corinne N Simonti
- Department of Biological Sciences, Georgia Institute of Technology, Atlanta, GA, USA
| | - Evonne McArthur
- Vanderbilt Genetics Institute, Vanderbilt University, Nashville, TN, USA
| | - Douglas Shaw
- Vanderbilt Genetics Institute, Vanderbilt University, Nashville, TN, USA
- Department of Biochemistry, Vanderbilt University, Nashville, TN, USA
| | - Emily Hodges
- Vanderbilt Genetics Institute, Vanderbilt University, Nashville, TN, USA
- Department of Biochemistry, Vanderbilt University, Nashville, TN, USA
| | - John A Capra
- Department of Biological Sciences, Vanderbilt University, Nashville, TN, USA.
- Vanderbilt Genetics Institute, Vanderbilt University, Nashville, TN, USA.
- Departments of Biomedical Informatics and Computer Science, Vanderbilt University, Nashville, TN, USA.
| |
Collapse
|
37
|
Frequency of genomic secondary findings among 21,915 eMERGE network participants. Genet Med 2020; 22:1470-1477. [PMID: 32546831 PMCID: PMC7713503 DOI: 10.1038/s41436-020-0810-9] [Citation(s) in RCA: 61] [Impact Index Per Article: 15.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/30/2019] [Revised: 02/24/2020] [Accepted: 04/06/2020] [Indexed: 12/31/2022] Open
Abstract
PURPOSE Discovering an incidental finding (IF) or secondary finding (SF) is a potential result of genomic testing, but few data exist describing types and frequencies of SFs likely to appear in broader clinical populations. METHODS The Electronic Medical Records and Genomics Network Phase III (eMERGE III) developed a CLIA-compliant sequencing panel of 109 genes and 1551 variants of clinical relevance or research interest and deployed this panel at ten clinical sites. We evaluated medically actionable SFs across 67 genes and 14 single-nucleotide variants (SNVs) in a diverse cohort of 21,915 participants drawn from a variety of settings (e.g., primary care, biobanks, specialty clinics). RESULTS Correcting for testing indication, we found a 3.02% overall frequency of SFs; 2.54% from 59 genes the American College of Medical Genetics and Genomics recommends for SF return, and 0.48% in other genes, primarily HFE and PALB2. SFs associated with cancer susceptibility were most frequent (1.38%), followed by cardiovascular diseases (0.87%), and lipid disorders (0.50%). After removing HFE, the frequency of SFs and proportion of pathogenic versus likely pathogenic SFs did not differ in those self-identifying as White versus others. CONCLUSION Here we present frequencies and types of medically actionable secondary findings to support informed decision making by patients, participants, and practitioners engaged in genomic medicine.
Collapse
|
38
|
Yue Z, Chu X, Xia J. PredCID: prediction of driver frameshift indels in human cancer. Brief Bioinform 2020; 22:5860690. [PMID: 32591774 DOI: 10.1093/bib/bbaa119] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2020] [Revised: 05/14/2020] [Accepted: 05/16/2020] [Indexed: 11/12/2022] Open
Abstract
The discrimination of driver from passenger mutations has been a hot topic in the field of cancer biology. Although recent advances have improved the identification of driver mutations in cancer genomic research, there is no computational method specific for the cancer frameshift indels (insertions or/and deletions) yet. In addition, existing pathogenic frameshift indel predictors may suffer from plenty of missing values because of different choices of transcripts during the variant annotation processes. In this study, we proposed a computational model, called PredCID (Predictor for Cancer driver frameshift InDels), for accurately predicting cancer driver frameshift indels. Gene, DNA, transcript and protein level features are combined together and selected for classification with eXtreme Gradient Boosting classifier. Benchmarking results on the cross-validation dataset and independent dataset showed that PredCID achieves better and robust performance compared with existing noncancer-specific methods in distinguishing cancer driver frameshift indels from passengers and is therefore a valuable method for deeper understanding of frameshift indels in human cancer. PredCID is freely available for academic research at http://bioinfo.ahu.edu.cn:8080/PredCID.
Collapse
Affiliation(s)
| | - Xinlu Chu
- Institutes of Physical Science and Information Technology, Anhui University
| | - Junfeng Xia
- Key Laboratory of Intelligent Computing and Signal Processing of Ministry of Education, Institutes of Physical Science and Information Technology, Anhui University
| |
Collapse
|
39
|
Perrier A, Sánchez-Castro D, Willi Y. Expressed mutational load increases toward the edge of a species' geographic range. Evolution 2020; 74:1711-1723. [PMID: 32538471 DOI: 10.1111/evo.14042] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2020] [Revised: 05/26/2020] [Accepted: 06/05/2020] [Indexed: 12/14/2022]
Abstract
There is no general explanation for why species have restricted geographic distributions. One hypothesis posits that range expansion or increasing scarcity of suitable habitat results in accumulation of mutational load due to enhanced genetic drift, which constrains population performance toward range limits and further expansion. We tested this hypothesis in the North American plant, Arabidopsis lyrata. We experimentally assessed mutational load by crossing plants of 20 populations from across the entire species range and by raising the offspring of within- and between-population crosses at five common garden sites within and beyond the range. Offspring performance was tracked over three growing seasons. The heterosis effect, depicting expressed mutational load, was increased in populations with heightened genomic estimates of load, longer expansion distance or long-term isolation, and a selfing mating system. The decline in performance of within-population crosses amounted to 80%. Mutation accumulation due to past range expansion and long-term isolation of populations in the area of range margins is therefore a strong determinant of population-mean performance, and the magnitude of effect may be sufficient to cause range limits.
Collapse
Affiliation(s)
- Antoine Perrier
- Department of Environmental Sciences, University of Basel, Basel, 4056, Switzerland
| | - Darío Sánchez-Castro
- Department of Environmental Sciences, University of Basel, Basel, 4056, Switzerland
| | - Yvonne Willi
- Department of Environmental Sciences, University of Basel, Basel, 4056, Switzerland
| |
Collapse
|
40
|
Gurdasani D, Carstensen T, Fatumo S, Chen G, Franklin CS, Prado-Martinez J, Bouman H, Abascal F, Haber M, Tachmazidou I, Mathieson I, Ekoru K, DeGorter MK, Nsubuga RN, Finan C, Wheeler E, Chen L, Cooper DN, Schiffels S, Chen Y, Ritchie GRS, Pollard MO, Fortune MD, Mentzer AJ, Garrison E, Bergström A, Hatzikotoulas K, Adeyemo A, Doumatey A, Elding H, Wain LV, Ehret G, Auer PL, Kooperberg CL, Reiner AP, Franceschini N, Maher D, Montgomery SB, Kadie C, Widmer C, Xue Y, Seeley J, Asiki G, Kamali A, Young EH, Pomilla C, Soranzo N, Zeggini E, Pirie F, Morris AP, Heckerman D, Tyler-Smith C, Motala AA, Rotimi C, Kaleebu P, Barroso I, Sandhu MS. Uganda Genome Resource Enables Insights into Population History and Genomic Discovery in Africa. Cell 2020; 179:984-1002.e36. [PMID: 31675503 DOI: 10.1016/j.cell.2019.10.004] [Citation(s) in RCA: 124] [Impact Index Per Article: 31.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2018] [Revised: 04/03/2019] [Accepted: 10/02/2019] [Indexed: 12/19/2022]
Abstract
Genomic studies in African populations provide unique opportunities to understand disease etiology, human diversity, and population history. In the largest study of its kind, comprising genome-wide data from 6,400 individuals and whole-genome sequences from 1,978 individuals from rural Uganda, we find evidence of geographically correlated fine-scale population substructure. Historically, the ancestry of modern Ugandans was best represented by a mixture of ancient East African pastoralists. We demonstrate the value of the largest sequence panel from Africa to date as an imputation resource. Examining 34 cardiometabolic traits, we show systematic differences in trait heritability between European and African populations, probably reflecting the differential impact of genes and environment. In a multi-trait pan-African GWAS of up to 14,126 individuals, we identify novel loci associated with anthropometric, hematological, lipid, and glycemic traits. We find that several functionally important signals are driven by Africa-specific variants, highlighting the value of studying diverse populations across the region.
Collapse
Affiliation(s)
- Deepti Gurdasani
- William Harvey Research Institute, Queen Mary's University of London, London, UK
| | | | - Segun Fatumo
- London School of Hygiene and Tropical Medicine, London, UK; Uganda Medical Informatics Centre (UMIC), MRC/UVRI and LSHTM (Uganda Research Unit), Entebbe, Uganda; H3Africa Bioinformatics Network (H3ABioNet) Node, Center for Genomics Research and Innovation (CGRI)/National Biotechnology Development Agency CGRI/NABDA, Abuja, Nigeria
| | - Guanjie Chen
- Center for Research on Genomics and Global Health, National Institute of Health, Bethesda, MD, USA
| | | | | | | | | | - Marc Haber
- Wellcome Sanger Institute, Hinxton, Cambridge, UK
| | - Ioanna Tachmazidou
- GSK Medicines Research Centre, Gunnels Wood Road, Stevenage Hertfordshire SG1 2NY, UK
| | - Iain Mathieson
- Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA
| | - Kenneth Ekoru
- Medical Research Council/Uganda Virus Research Institute (MRC/UVRI) and London School of Hygiene & Tropical Medicine Uganda Research Unit on AIDS, Entebbe, Uganda; Department of Medicine, University of Cambridge, Cambridge, UK
| | - Marianne K DeGorter
- Department of Pathology, Stanford University School of Medicine, Stanford, CA, USA
| | - Rebecca N Nsubuga
- Medical Research Council/Uganda Virus Research Institute (MRC/UVRI) and London School of Hygiene & Tropical Medicine Uganda Research Unit on AIDS, Entebbe, Uganda
| | - Chris Finan
- Wellcome Sanger Institute, Hinxton, Cambridge, UK
| | - Eleanor Wheeler
- Wellcome Sanger Institute, Hinxton, Cambridge, UK; MRC Epidemiology Unit, University of Cambridge, Cambridge, UK
| | - Li Chen
- Wellcome Sanger Institute, Hinxton, Cambridge, UK
| | - David N Cooper
- Institute of Medical Genetics, School of Medicine, Cardiff University, Cardiff, UK
| | - Stephan Schiffels
- Department of Archaeogenetics, Max Planck Institute for the Science of Human History, Jena, Germany
| | - Yuan Chen
- Wellcome Sanger Institute, Hinxton, Cambridge, UK
| | | | | | | | - Alex J Mentzer
- The Wellcome Centre for Human Genetics, University of Oxford, Oxford, UK
| | | | | | - Konstantinos Hatzikotoulas
- Wellcome Sanger Institute, Hinxton, Cambridge, UK; Institute of Translational Genomics, Helmholtz Zentrum München, German Research Center for Environmental Health, Neuherberg, Germany
| | - Adebowale Adeyemo
- Center for Research on Genomics and Global Health, National Institute of Health, Bethesda, MD, USA
| | - Ayo Doumatey
- Center for Research on Genomics and Global Health, National Institute of Health, Bethesda, MD, USA
| | | | - Louise V Wain
- Department of Health Sciences, University of Leicester, Leicester, UK; National Institute for Health Research, Leicester Biomedical Research Centre, Glenfield Hospital, Leicester, UK
| | - Georg Ehret
- McKusick-Nathans Institute of Genetic Medicine, Johns Hopkins University School of Medicine, Baltimore, MD, USA; Geneva University Hospitals, Rue Gabrielle-Perret-Gentil 4, 1211 Genève 14, Switzerland
| | - Paul L Auer
- Zilber School of Public Health, University of Wisconsin-Milwaukee, Milwaukee, WI, USA
| | - Charles L Kooperberg
- Division of Public Health Sciences, Fred Hutchinson Cancer Research Center, Seattle, WA, USA
| | - Alexander P Reiner
- Department of Epidemiology, University of Washington, Seattle, WA, USA; Division of Public Health Sciences, Fred Hutchinson Cancer Research Center, Seattle, WA, USA
| | - Nora Franceschini
- Department of Epidemiology, University of North Carolina, Chapel Hill, NC, USA
| | - Dermot Maher
- Medical Research Council/Uganda Virus Research Institute (MRC/UVRI) and London School of Hygiene & Tropical Medicine Uganda Research Unit on AIDS, Entebbe, Uganda
| | - Stephen B Montgomery
- Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, USA; Department of Genetics, Stanford University School of Medicine, Stanford, CA, USA
| | | | | | - Yali Xue
- Wellcome Sanger Institute, Hinxton, Cambridge, UK
| | - Janet Seeley
- London School of Hygiene and Tropical Medicine, London, UK; Medical Research Council/Uganda Virus Research Institute (MRC/UVRI) and London School of Hygiene & Tropical Medicine Uganda Research Unit on AIDS, Entebbe, Uganda
| | - Gershim Asiki
- Medical Research Council/Uganda Virus Research Institute (MRC/UVRI) and London School of Hygiene & Tropical Medicine Uganda Research Unit on AIDS, Entebbe, Uganda
| | - Anatoli Kamali
- Medical Research Council/Uganda Virus Research Institute (MRC/UVRI) and London School of Hygiene & Tropical Medicine Uganda Research Unit on AIDS, Entebbe, Uganda
| | - Elizabeth H Young
- Wellcome Sanger Institute, Hinxton, Cambridge, UK; Department of Medicine, University of Cambridge, Cambridge, UK
| | - Cristina Pomilla
- Wellcome Sanger Institute, Hinxton, Cambridge, UK; Department of Medicine, University of Cambridge, Cambridge, UK
| | - Nicole Soranzo
- Wellcome Sanger Institute, Hinxton, Cambridge, UK; Department of Haematology, University of Cambridge, Cambridge, UK; The National Institute for Health Research Blood and Transplant Unit (NIHR BTRU) in Donor Health and Genomics, University of Cambridge, Cambridge, UK
| | - Eleftheria Zeggini
- Institute of Translational Genomics, Helmholtz Zentrum München, German Research Center for Environmental Health, Neuherberg, Germany
| | - Fraser Pirie
- Department of Diabetes and Endocrinology, University of KwaZulu-Natal, Durban, South Africa
| | - Andrew P Morris
- The Wellcome Centre for Human Genetics, University of Oxford, Oxford, UK; Department of Biostatistics, University of Liverpool, Liverpool, UK
| | | | | | - Ayesha A Motala
- Department of Diabetes and Endocrinology, University of KwaZulu-Natal, Durban, South Africa.
| | - Charles Rotimi
- Center for Research on Genomics and Global Health, National Institute of Health, Bethesda, MD, USA.
| | - Pontiano Kaleebu
- London School of Hygiene and Tropical Medicine, London, UK; Uganda Medical Informatics Centre (UMIC), MRC/UVRI and LSHTM (Uganda Research Unit), Entebbe, Uganda; Medical Research Council/Uganda Virus Research Institute (MRC/UVRI) and London School of Hygiene & Tropical Medicine Uganda Research Unit on AIDS, Entebbe, Uganda.
| | - Inês Barroso
- Wellcome Sanger Institute, Hinxton, Cambridge, UK; MRC Epidemiology Unit, University of Cambridge, Cambridge, UK.
| | - Manj S Sandhu
- Department of Medicine, University of Cambridge, Cambridge, UK.
| |
Collapse
|
41
|
Ulaganathan VK. TraPS-VarI: Identifying genetic variants altering phosphotyrosine based signalling motifs. Sci Rep 2020; 10:8453. [PMID: 32439998 PMCID: PMC7242328 DOI: 10.1038/s41598-020-65146-2] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2019] [Accepted: 04/26/2020] [Indexed: 12/17/2022] Open
Abstract
Patient stratification and individualized therapeutic strategies rely on the established knowledge of genotype-specific molecular and cellular alterations of biological and therapeutic significance. Whilst almost all approved drugs have been developed based on the Reference Sequence protein database (RefSeq), the latest genome sequencing studies establish the substantial prevalence of non-synonymous genetic mutations in the general population, including stop-insertion and frame shift mutations within the coding regions of membrane proteins. While the availability of individual genotypes are becoming increasingly common, the biological and clinical interpretations of mutations among individual genomes is largely lagging behind. Lately, transmembrane proteins of haematopoietic (myeloid and lymphoid) derived immune cells have attracted much attention as important targets for cancer immunotherapies. As such, the signalling properties of haematological transmembrane receptors rely on the membrane-proximal phosphotyrosine based sequence motifs (TBSMs) such as ITAM (immunoreceptor tyrosine-based activation motif), ITIM (immunoreceptor tyrosine-based inhibition motif) and signal transducer and activator of transcription 3 (STAT3)-recruiting YxxQ motifs. However, mutations that alter the coding regions of transmembrane proteins, resulting in either insertion or deletion of crucial signal modulating TBSMs, remains unknown. To conveniently identify individual cell line-specific or patient-specific membrane protein altering mutations, we present the Transmembrane Protein Sequence Variant Identifier (TraPS-VarI). TraPS-VarI is an annotation tool for accurate mapping of the effect of an individual's mutation in the transmembrane protein sequence, and to identify the prevalence of TBSMs. TraPS-VarI is a biologist and clinician-friendly algorithm with a web interface and an associated database browser (https://www.traps-vari.org/).
Collapse
Affiliation(s)
- Vijay Kumar Ulaganathan
- Department of Molecular Biology, Max Planck Institute of Biochemistry, Am Klopferspitz 18, Martinsried, 82152, Germany.
- Department of Neuroimmunology, Universitätsmedizin Göttingen, Von-Siebold-Str. 3A, Göttingen, 37075, Germany.
| |
Collapse
|
42
|
de Manuel M, Barnett R, Sandoval-Velasco M, Yamaguchi N, Garrett Vieira F, Zepeda Mendoza ML, Liu S, Martin MD, Sinding MHS, Mak SST, Carøe C, Liu S, Guo C, Zheng J, Zazula G, Baryshnikov G, Eizirik E, Koepfli KP, Johnson WE, Antunes A, Sicheritz-Ponten T, Gopalakrishnan S, Larson G, Yang H, O'Brien SJ, Hansen AJ, Zhang G, Marques-Bonet T, Gilbert MTP. The evolutionary history of extinct and living lions. Proc Natl Acad Sci U S A 2020; 117:10927-10934. [PMID: 32366643 PMCID: PMC7245068 DOI: 10.1073/pnas.1919423117] [Citation(s) in RCA: 45] [Impact Index Per Article: 11.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/23/2022] Open
Abstract
Lions are one of the world's most iconic megafauna, yet little is known about their temporal and spatial demographic history and population differentiation. We analyzed a genomic dataset of 20 specimens: two ca. 30,000-y-old cave lions (Panthera leo spelaea), 12 historic lions (Panthera leo leo/Panthera leo melanochaita) that lived between the 15th and 20th centuries outside the current geographic distribution of lions, and 6 present-day lions from Africa and India. We found that cave and modern lions shared an ancestor ca. 500,000 y ago and that the 2 lineages likely did not hybridize following their divergence. Within modern lions, we found 2 main lineages that diverged ca. 70,000 y ago, with clear evidence of subsequent gene flow. Our data also reveal a nearly complete absence of genetic diversity within Indian lions, probably due to well-documented extremely low effective population sizes in the recent past. Our results contribute toward the understanding of the evolutionary history of lions and complement conservation efforts to protect the diversity of this vulnerable species.
Collapse
Affiliation(s)
- Marc de Manuel
- Institute of Evolutionary Biology (UPF-CSIC), PRBB, 08003 Barcelona, Spain
| | - Ross Barnett
- Section for Evolutionary Genomics, The GLOBE Institute, Faculty of Health and Medical Sciences, University of Copenhagen, 1353 Copenhagen, Denmark
| | - Marcela Sandoval-Velasco
- Section for Evolutionary Genomics, The GLOBE Institute, Faculty of Health and Medical Sciences, University of Copenhagen, 1353 Copenhagen, Denmark
| | - Nobuyuki Yamaguchi
- Institute of Tropical Biodiversity and Sustainable Development, University Malaysia Terengganu, 21030 Kuala Nerus, Terengganu, Malaysia;
| | - Filipe Garrett Vieira
- Section for Evolutionary Genomics, The GLOBE Institute, Faculty of Health and Medical Sciences, University of Copenhagen, 1353 Copenhagen, Denmark
| | - M Lisandra Zepeda Mendoza
- Section for Evolutionary Genomics, The GLOBE Institute, Faculty of Health and Medical Sciences, University of Copenhagen, 1353 Copenhagen, Denmark
- School of Medical and Dental Sciences, Institute of Microbiology and Infection, University of Birmingham, B15 2TT Edgbaston, Birmingham, United Kingdom
| | | | - Michael D Martin
- Norwegian University of Science and Technology (NTNU) University Museum, 7012 Trondheim, Norway
| | - Mikkel-Holger S Sinding
- Section for Evolutionary Genomics, The GLOBE Institute, Faculty of Health and Medical Sciences, University of Copenhagen, 1353 Copenhagen, Denmark
| | - Sarah S T Mak
- Section for Evolutionary Genomics, The GLOBE Institute, Faculty of Health and Medical Sciences, University of Copenhagen, 1353 Copenhagen, Denmark
| | - Christian Carøe
- Section for Evolutionary Genomics, The GLOBE Institute, Faculty of Health and Medical Sciences, University of Copenhagen, 1353 Copenhagen, Denmark
| | - Shanlin Liu
- Section for Evolutionary Genomics, The GLOBE Institute, Faculty of Health and Medical Sciences, University of Copenhagen, 1353 Copenhagen, Denmark
- BGI-Shenzhen, 518083 Shenzhen, China
| | | | - Jiao Zheng
- BGI-Shenzhen, 518083 Shenzhen, China
- BGI Education Center, University of Chinese Academy of Sciences, 518083 Shenzhen, China
| | - Grant Zazula
- Yukon Palaeontology Program, Department of Tourism and Culture, Government of Yukon, Y1A 2C6 Whitehorse, Yukon, Canada
| | - Gennady Baryshnikov
- Zoological Institute, Russian Academy of Sciences, 199034 St. Petersburg, Russia
| | - Eduardo Eizirik
- Laboratory of Genomics and Molecular Biology, Escola de Ciências da Saúde e da Vida, Pontifícia Universidade Católica do Rio Grande do Sul (PUCRS), Porto Alegre, RS 90619-900, Brazil
- Instituto Nacional de Ciência e Tecnologia - Ecologia Evolução e Conservação da Biodiversidade (INCT-EECBio), Goiânia, GO 74690-900, Brazil
- Instituto Pró-Carnívoros, Atibaia, SP 12945-010, Brazil
| | - Klaus-Peter Koepfli
- Center for Species Survival, Smithsonian Conservation Biology Institute, National Zoological Park, Front Royal, VA 22630
| | - Warren E Johnson
- Center for Species Survival, Smithsonian Conservation Biology Institute, National Zoological Park, Front Royal, VA 22630
- The Walter Reed Biosystematics Unit, Museum Support Center MRC-534, Smithsonian Institution, Suitland, MD 20746-2863
- Walter Reed Army Institute of Research, Silver Spring, MD 20910
| | - Agostinho Antunes
- Interdisciplinary Centre of Marine and Environmental Research (CIIMAR/CIMAR), University of Porto, 4450-208 Porto, Portugal
- Department of Biology, Faculty of Sciences, University of Porto, 4169-007 Porto, Portugal
| | - Thomas Sicheritz-Ponten
- Section for Evolutionary Genomics, The GLOBE Institute, Faculty of Health and Medical Sciences, University of Copenhagen, 1353 Copenhagen, Denmark
- Centre of Excellence for Omics-Driven Computational Biodiscovery (COMBio), Faculty of Applied Sciences, Asian Institute of Medicine, Science and Technology (AIMST), 08100 Bedong, Kedah, Malaysia
| | - Shyam Gopalakrishnan
- Section for Evolutionary Genomics, The GLOBE Institute, Faculty of Health and Medical Sciences, University of Copenhagen, 1353 Copenhagen, Denmark
| | - Greger Larson
- The Palaeogenomics and Bio-Archaeology Research Network, Research Laboratory for Archaeology and the History of Art, University of Oxford, OX1 3QY Oxford, United Kingdom
| | - Huanming Yang
- BGI-Shenzhen, 518083 Shenzhen, China
- James D. Watson Institute of Genome Science, 310008 Hangzhou, China
| | - Stephen J O'Brien
- Laboratory of Genomic Diversity, Center for Computer Technologies, ITMO (Information Technologies, Mechanics and Optics) University, 197101 St. Petersburg, Russia;
- Guy Harvey Oceanographic Center Halmos College of Natural Sciences and Oceanography, Nova Southeastern University, Ft. Lauderdale, FL 33004
| | - Anders J Hansen
- Section for GeoGenetics, The GLOBE Institute, Faculty of Health and Medical Sciences, University of Copenhagen, 1350 Copenhagen, Denmark
| | - Guojie Zhang
- BGI-Shenzhen, 518083 Shenzhen, China
- Section for Ecology and Evolution, Department of Biology, Faculty of Science, University of Copenhagen, 2100 Copenhagen, Denmark
- State Key Laboratory of Genetic Resources and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, 650223 Kunming, China
| | - Tomas Marques-Bonet
- Institute of Evolutionary Biology (UPF-CSIC), PRBB, 08003 Barcelona, Spain;
- Centre Nacional d'Anàlisi Genòmica, Centre for Genomic Regulation (CNAG-CRG), The Barcelona Institute of Science and Technology, 08028 Barcelona, Spain
- Institució Catalana de Recerca i Estudis Avançats (ICREA), 08003 Barcelona, Spain
- Institut Català de Paleontologia Miquel Crusafont, Universitat Autònoma de Barcelona, Edifici ICTA-ICP, 08193 Cerdanyola del Vallès, Barcelona, Spain
| | - M Thomas P Gilbert
- Section for Evolutionary Genomics, The GLOBE Institute, Faculty of Health and Medical Sciences, University of Copenhagen, 1353 Copenhagen, Denmark;
- Norwegian University of Science and Technology (NTNU) University Museum, 7012 Trondheim, Norway
| |
Collapse
|
43
|
Yoshida K, Ravinet M, Makino T, Toyoda A, Kokita T, Mori S, Kitano J. Accumulation of Deleterious Mutations in Landlocked Threespine Stickleback Populations. Genome Biol Evol 2020; 12:479-492. [PMID: 32232440 PMCID: PMC7197494 DOI: 10.1093/gbe/evaa065] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 03/24/2020] [Indexed: 02/06/2023] Open
Abstract
Colonization of new habitats often reduces population sizes and may result in the accumulation of deleterious mutations by genetic drift. Compared with the genomic basis for adaptation to new environments, genome-wide analysis of deleterious mutations in isolated populations remains limited. In the present study, we investigated the accumulation of deleterious mutations in five endangered freshwater populations of threespine stickleback (Gasterosteus aculeatus) in the central part of the mainland of Japan. Using whole-genome resequencing data, we first conducted phylogenomic analysis and confirmed at least two independent freshwater colonization events in the central mainland from ancestral marine ecotypes. Next, analyses of single nucleotide polymorphisms showed a substantial reduction of heterozygosity in freshwater populations compared with marine populations. Reduction in heterozygosity was more apparent at the center of each chromosome than the peripheries and on X chromosomes compared with autosomes. Third, bioinformatic analysis of deleterious mutations showed increased accumulation of putatively deleterious mutations in the landlocked freshwater populations compared with marine populations. For the majority of populations examined, the frequencies of putatively deleterious mutations were higher on X chromosomes than on autosomes. The interpopulation comparison indicated that the majority of putatively deleterious mutations may have accumulated independently. Thus, whole-genome resequencing of endangered populations can help to estimate the accumulation of deleterious mutations and inform us of which populations are the most severely endangered. Furthermore, analysis of variation among chromosomes can give insights into whether any particular chromosomes are likely to accumulate deleterious mutations.
Collapse
Affiliation(s)
- Kohta Yoshida
- Ecological Genetics Laboratory, National Institute of Genetics, Mishima, Shizuoka, Japan.,Max Planck Institute for Developmental Biology, Tübingen, Germany
| | - Mark Ravinet
- Ecological Genetics Laboratory, National Institute of Genetics, Mishima, Shizuoka, Japan.,Department of Biosciences, Centre for Ecological and Evolutionary Synthesis, University of Oslo, Norway.,School of Life Sciences, University of Nottingham, University Park, Nottingham, United Kingdom
| | - Takashi Makino
- Graduate School of Life Sciences, Tohoku University, Sendai, Miyagi, Japan
| | - Atsushi Toyoda
- Comparative Genomics Laboratory, National Institute of Genetics, Mishima, Shizuoka, Japan
| | - Tomoyuki Kokita
- Department of Marine Bioscience, Fukui Prefectural University, Obama, Fukui, Japan
| | - Seiichi Mori
- Biological Laboratories, Gifu-kyoritsu University, Ogaki, Gifu, Japan
| | - Jun Kitano
- Ecological Genetics Laboratory, National Institute of Genetics, Mishima, Shizuoka, Japan
| |
Collapse
|
44
|
Wang W, Chen C, Wang X, Zhang L, Shen D, Wang S, Gao B, Mao J, Song C. Development of Molecular Markers Based on the L1 Retrotransposon Insertion Polymorphisms in Pigs (Sus scrofa) and Their Association with Economic Traits. RUSS J GENET+ 2020. [DOI: 10.1134/s1022795420020131] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022]
|
45
|
Aris-Brosou S. Direct Evidence of an Increasing Mutational Load in Humans. Mol Biol Evol 2020; 36:2823-2829. [PMID: 31424543 DOI: 10.1093/molbev/msz192] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
The extent to which selection has shaped present-day human populations has attracted intense scrutiny, and examples of local adaptations abound. However, the evolutionary trajectory of alleles that, today, are deleterious has received much less attention. To address this question, the genomes of 2,062 individuals, including 1,179 ancient humans, were reanalyzed to assess how frequencies of risk alleles and their homozygosity changed through space and time in Europe over the past 45,000 years. Although the overall deleterious homozygosity has consistently decreased, risk alleles have steadily increased in frequency over that period of time. Those that increased most are associated with diseases such as asthma, Crohn disease, diabetes, and obesity, which are highly prevalent in present-day populations. These findings may not run against the existence of local adaptations but highlight the limitations imposed by drift and population dynamics on the strength of selection in purging deleterious mutations from human populations.
Collapse
Affiliation(s)
- Stéphane Aris-Brosou
- Department of Biology, University of Ottawa, Ottawa, ON, Canada.,Department of Mathematics and Statistics, University of Ottawa, Ottawa, ON, Canada
| |
Collapse
|
46
|
Beichman AC, Koepfli KP, Li G, Murphy W, Dobrynin P, Kliver S, Tinker MT, Murray MJ, Johnson J, Lindblad-Toh K, Karlsson EK, Lohmueller KE, Wayne RK. Aquatic Adaptation and Depleted Diversity: A Deep Dive into the Genomes of the Sea Otter and Giant Otter. Mol Biol Evol 2019; 36:2631-2655. [PMID: 31212313 PMCID: PMC7967881 DOI: 10.1093/molbev/msz101] [Citation(s) in RCA: 37] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Abstract
Despite its recent invasion into the marine realm, the sea otter (Enhydra lutris) has evolved a suite of adaptations for life in cold coastal waters, including limb modifications and dense insulating fur. This uniquely dense coat led to the near-extinction of sea otters during the 18th-20th century fur trade and an extreme population bottleneck. We used the de novo genome of the southern sea otter (E. l. nereis) to reconstruct its evolutionary history, identify genes influencing aquatic adaptation, and detect signals of population bottlenecks. We compared the genome of the southern sea otter with the tropical freshwater-living giant otter (Pteronura brasiliensis) to assess common and divergent genomic trends between otter species, and with the closely related northern sea otter (E. l. kenyoni) to uncover population-level trends. We found signals of positive selection in genes related to aquatic adaptations, particularly limb development and polygenic selection on genes related to hair follicle development. We found extensive pseudogenization of olfactory receptor genes in both the sea otter and giant otter lineages, consistent with patterns of sensory gene loss in other aquatic mammals. At the population level, the southern sea otter and the northern sea otter showed extremely low genomic diversity, signals of recent inbreeding, and demographic histories marked by population declines. These declines may predate the fur trade and appear to have resulted in an increase in putatively deleterious variants that could impact the future recovery of the sea otter.
Collapse
Affiliation(s)
- Annabel C Beichman
- Department of Ecology and Evolutionary Biology, University of California, Los Angeles, CA
| | - Klaus-Peter Koepfli
- Center for Species Survival, Smithsonian Conservation Biology Institute, National Zoological Park, Washington, DC
- Institute of Molecular and Cellular Biology, Siberian Branch of the Russian Academy of Sciences, Novosibirsk, Russian Federation
| | - Gang Li
- College of Life Science, Shaanxi Normal University, Xi’an, Shaanxi, China
| | - William Murphy
- Department of Veterinary Integrative Biosciences, Texas A&M University, College Station, TX
| | - Pasha Dobrynin
- Center for Species Survival, Smithsonian Conservation Biology Institute, National Zoological Park, Washington, DC
- Institute of Molecular and Cellular Biology, Siberian Branch of the Russian Academy of Sciences, Novosibirsk, Russian Federation
| | - Sergei Kliver
- Institute of Molecular and Cellular Biology, Siberian Branch of the Russian Academy of Sciences, Novosibirsk, Russian Federation
| | - Martin T Tinker
- Department of Ecology and Evolutionary Biology, University of California, Santa Cruz, CA
| | | | - Jeremy Johnson
- Vertebrate Genome Biology, Broad Institute of MIT and Harvard, Cambridge, MA
| | - Kerstin Lindblad-Toh
- Vertebrate Genome Biology, Broad Institute of MIT and Harvard, Cambridge, MA
- Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden
| | - Elinor K Karlsson
- Vertebrate Genome Biology, Broad Institute of MIT and Harvard, Cambridge, MA
- Bioinformatics and Integrative Biology, University of Massachusetts Medical School, Worcester, MA
| | - Kirk E Lohmueller
- Department of Ecology and Evolutionary Biology, University of California, Los Angeles, CA
- Interdepartmental Program in Bioinformatics, University of California, Los Angeles, CA
- Department of Human Genetics, David Geffen School of Medicine, University of California, Los Angeles, CA
| | - Robert K Wayne
- Department of Ecology and Evolutionary Biology, University of California, Los Angeles, CA
| |
Collapse
|
47
|
Mai TL, Chuang TJ. A-to-I RNA editing contributes to the persistence of predicted damaging mutations in populations. Genome Res 2019; 29:1766-1776. [PMID: 31515285 PMCID: PMC6836733 DOI: 10.1101/gr.246033.118] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/05/2018] [Accepted: 09/04/2019] [Indexed: 12/13/2022]
Abstract
Adenosine-to-inosine (A-to-I) RNA editing is a very common co-/posttranscriptional modification that can lead to A-to-G changes at the RNA level and compensate for G-to-A genomic changes to a certain extent. It has been shown that each healthy individual can carry dozens of missense variants predicted to be severely deleterious. Why strongly detrimental variants are preserved in a population and not eliminated by negative natural selection remains mostly unclear. Here, we ask if RNA editing correlates with the burden of deleterious A/G polymorphisms in a population. Integrating genome and transcriptome sequencing data from 447 human lymphoblastoid cell lines, we show that nonsynonymous editing activities (prevalence/level) are negatively correlated with the deleteriousness of A-to-G genomic changes and positively correlated with that of G-to-A genomic changes within the population. We find a significantly negative correlation between nonsynonymous editing activities and allele frequency of A within the population. This negative editing-allele frequency correlation is particularly strong when editing sites are located in highly important genes/loci. Examinations of deleterious missense variants from the 1000 Genomes Project further show a significantly higher proportion of rare missense mutations for G-to-A changes than for other types of changes. The proportion for G-to-A changes increases with increasing deleterious effects of the changes. Moreover, the deleteriousness of G-to-A changes is significantly positively correlated with the percentage of editing enzyme binding motifs at the variants. Overall, we show that nonsynonymous editing is associated with the increased burden of G-to-A missense mutations in healthy individuals, expanding RNA editing in pathogenomics studies.
Collapse
Affiliation(s)
- Te-Lun Mai
- Genomics Research Center, Academia Sinica, Taipei 11529, Taiwan
| | | |
Collapse
|
48
|
Fuselli S. Beyond drugs: the evolution of genes involved in human response to medications. Proc Biol Sci 2019; 286:20191716. [PMID: 31640517 DOI: 10.1098/rspb.2019.1716] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/29/2022] Open
Abstract
The genetic variation of our species reflects human demographic history and adaptation to diverse local environments. Part of this genetic variation affects individual responses to exogenous substances, such as food, pollutants and drugs, and plays an important role in drug efficacy and safety. This review provides a synthesis of the evolution of loci implicated in human pharmacological response and metabolism, interpreted within the theoretical framework of population genetics and molecular evolution. In particular, I review and discuss key evolutionary aspects of different pharmacogenes in humans and other species, such as the relationship between the type of substrates and rate of evolution; the selective pressure exerted by landscape variables or dietary habits; expected and observed patterns of rare genetic variation. Finally, I discuss how this knowledge can be translated directly or after the implementation of specific studies, into practical guidelines.
Collapse
Affiliation(s)
- Silvia Fuselli
- Department of Life Sciences and Biotechnology, University of Ferrara, Ferrara, Italy
| |
Collapse
|
49
|
Szpiech ZA, Mak ACY, White MJ, Hu D, Eng C, Burchard EG, Hernandez RD. Ancestry-Dependent Enrichment of Deleterious Homozygotes in Runs of Homozygosity. Am J Hum Genet 2019; 105:747-762. [PMID: 31543216 PMCID: PMC6817522 DOI: 10.1016/j.ajhg.2019.08.011] [Citation(s) in RCA: 23] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2019] [Accepted: 08/27/2019] [Indexed: 12/20/2022] Open
Abstract
Runs of homozygosity (ROH) are important genomic features that manifest when an individual inherits two haplotypes that are identical by descent. Their length distributions are informative about population history, and their genomic locations are useful for mapping recessive loci contributing to both Mendelian and complex disease risk. We have previously shown that ROH, and especially long ROH that are likely the result of recent parental relatedness, are enriched for homozygous deleterious coding variation in a worldwide sample of outbred individuals. However, the distribution of ROH in admixed populations and their relationship to deleterious homozygous genotypes is understudied. Here we analyze whole-genome sequencing data from 1,441 unrelated individuals from self-identified African American, Puerto Rican, and Mexican American populations. These populations are three-way admixed between European, African, and Native American ancestries and provide an opportunity to study the distribution of deleterious alleles partitioned by local ancestry and ROH. We re-capitulate previous findings that long ROH are enriched for deleterious variation genome-wide. We then partition by local ancestry and show that deleterious homozygotes arise at a higher rate when ROH overlap African ancestry segments than when they overlap European or Native American ancestry segments of the genome. These results suggest that, while ROH on any haplotype background are associated with an inflation of deleterious homozygous variation, African haplotype backgrounds may play a particularly important role in the genetic architecture of complex diseases for admixed individuals, highlighting the need for further study of these populations.
Collapse
Affiliation(s)
- Zachary A Szpiech
- Department of Bioengineering and Therapeutic Sciences, University of California San Francisco, San Francisco, CA 94158, USA; Department of Biological Sciences, Auburn University, Auburn, AL 36842, USA.
| | - Angel C Y Mak
- Department of Medicine, University of California San Francisco, San Francisco, CA 94158, USA
| | - Marquitta J White
- Department of Medicine, University of California San Francisco, San Francisco, CA 94158, USA
| | - Donglei Hu
- Department of Medicine, University of California San Francisco, San Francisco, CA 94158, USA
| | - Celeste Eng
- Department of Medicine, University of California San Francisco, San Francisco, CA 94158, USA
| | - Esteban G Burchard
- Department of Medicine, University of California San Francisco, San Francisco, CA 94158, USA
| | - Ryan D Hernandez
- Department of Bioengineering and Therapeutic Sciences, University of California San Francisco, San Francisco, CA 94158, USA; Institute for Human Genetics, University of California San Francisco, San Francisco, CA 94158, USA; Quantitative Biosciences Institute, University of California San Francisco, San Francisco, CA 94158, USA; Department of Human Genetics, McGill University, Montreal, QC H3A 0G1, Canada; Genome Quebec Innovation Center, McGill University, Montreal, QC H3A 0G1, Canada.
| |
Collapse
|
50
|
Han S, Andrés AM, Marques-Bonet T, Kuhlwilm M. Genetic Variation in Pan Species Is Shaped by Demographic History and Harbors Lineage-Specific Functions. Genome Biol Evol 2019; 11:1178-1191. [PMID: 30847478 PMCID: PMC6482415 DOI: 10.1093/gbe/evz047] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 03/05/2019] [Indexed: 01/08/2023] Open
Abstract
Chimpanzees (Pan troglodytes) and bonobos (Pan paniscus) are the closest living relatives of humans, but the two species show distinct behavioral and physiological differences, particularly regarding female reproduction. Despite their recent rapid decline, the demographic histories of the two species have been different during the past 1–2 Myr, likely having an impact on their genomic diversity. Here, we analyze the inferred functional consequences of genetic variation across 69 individuals, making use of the most complete data set of genomes in the Pan clade to date. We test to which extent the demographic history influences the efficacy of purifying selection in these species. We find that small historical effective population sizes (Ne) correlate not only with low levels of genetic diversity but also with a larger number of deleterious alleles in homozygosity and an increased proportion of deleterious changes at low frequencies. To investigate the putative genetic basis for phenotypic differences between chimpanzees and bonobos, we exploit the catalog of putatively deleterious protein-coding changes in each lineage. We show that bonobo-specific nonsynonymous changes are enriched in genes related to age at menarche in humans, suggesting that the prominent physiological differences in the female reproductive system between chimpanzees and bonobos might be explained, in part, by putatively adaptive changes on the bonobo lineage.
Collapse
Affiliation(s)
- Sojung Han
- Institut de Biologia Evolutiva, Consejo Superior de Investigaciones Científicas-Universitat Pompeu Fabra, Barcelona, Spain
| | - Aida M Andrés
- Department of Genetics, Max Planck Institute for Evolutionary Anthropology, Leipzig, Germany.,Department of Genetics, Evolution and Environment, UCL Genetics Institute, University College London, London, United Kingdom
| | - Tomas Marques-Bonet
- Institut de Biologia Evolutiva, Consejo Superior de Investigaciones Científicas-Universitat Pompeu Fabra, Barcelona, Spain.,Institucio Catalana de Recerca i Estudis Avançats (ICREA), 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, Edifici ICTA-ICP, Barcelona, Spain
| | - Martin Kuhlwilm
- Institut de Biologia Evolutiva, Consejo Superior de Investigaciones Científicas-Universitat Pompeu Fabra, Barcelona, Spain
| |
Collapse
|