1
|
Wang X, Heckel G. Genome-wide relaxation of selection and the evolution of the island syndrome in Orkney voles. Genome Res 2024; 34:851-862. [PMID: 38955466 PMCID: PMC11293545 DOI: 10.1101/gr.278487.123] [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/07/2023] [Accepted: 05/14/2024] [Indexed: 07/04/2024]
Abstract
Island populations often experience different ecological and demographic conditions than their counterparts on the continent, resulting in divergent evolutionary forces affecting their genomes. Random genetic drift and selection both may leave their imprints on island populations, although the relative impact depends strongly on the specific conditions. Here we address their contributions to the island syndrome in a rodent with an unusually clear history of isolation. Common voles (Microtus arvalis) were introduced by humans on the Orkney archipelago north of Scotland >5000 years ago and rapidly evolved to exceptionally large size. Our analyses show that the genomes of Orkney voles were dominated by genetic drift, with extremely low diversity, variable Tajima's D, and very high divergence from continental conspecifics. Increased d N/d S ratios over a wide range of genes in Orkney voles indicated genome-wide relaxation of purifying selection. We found evidence of hard sweeps on key genes of the lipid metabolism pathway only in continental voles. The marked increase of body size in Orkney-a typical phenomenon of the island syndrome-may thus be associated to the relaxation of positive selection on genes related to this pathway. On the other hand, a hard sweep on immune genes of Orkney voles likely reflects the divergent ecological conditions and possibly the history of human introduction. The long-term isolated Orkney voles show that adaptive changes may still impact the evolutionary trajectories of such populations despite the pervasive consequences of genetic drift at the genome level.
Collapse
Affiliation(s)
- Xuejing Wang
- Institute of Ecology and Evolution, University of Bern, 3012 Bern, Switzerland
| | - Gerald Heckel
- Institute of Ecology and Evolution, University of Bern, 3012 Bern, Switzerland;
- Swiss Institute of Bioinformatics, 1015 Lausanne, Switzerland
| |
Collapse
|
2
|
Hemstrom W, Grummer JA, Luikart G, Christie MR. Next-generation data filtering in the genomics era. Nat Rev Genet 2024:10.1038/s41576-024-00738-6. [PMID: 38877133 DOI: 10.1038/s41576-024-00738-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 04/25/2024] [Indexed: 06/16/2024]
Abstract
Genomic data are ubiquitous across disciplines, from agriculture to biodiversity, ecology, evolution and human health. However, these datasets often contain noise or errors and are missing information that can affect the accuracy and reliability of subsequent computational analyses and conclusions. A key step in genomic data analysis is filtering - removing sequencing bases, reads, genetic variants and/or individuals from a dataset - to improve data quality for downstream analyses. Researchers are confronted with a multitude of choices when filtering genomic data; they must choose which filters to apply and select appropriate thresholds. To help usher in the next generation of genomic data filtering, we review and suggest best practices to improve the implementation, reproducibility and reporting standards for filter types and thresholds commonly applied to genomic datasets. We focus mainly on filters for minor allele frequency, missing data per individual or per locus, linkage disequilibrium and Hardy-Weinberg deviations. Using simulated and empirical datasets, we illustrate the large effects of different filtering thresholds on common population genetics statistics, such as Tajima's D value, population differentiation (FST), nucleotide diversity (π) and effective population size (Ne).
Collapse
Affiliation(s)
- William Hemstrom
- Department of Biological Sciences, Purdue University, West Lafayette, IN, USA.
| | - Jared A Grummer
- Flathead Lake Biological Station, Wildlife Biology Program and Division of Biological Sciences, University of Montana, Missoula, MT, USA
| | - Gordon Luikart
- Flathead Lake Biological Station, Wildlife Biology Program and Division of Biological Sciences, University of Montana, Missoula, MT, USA
| | - Mark R Christie
- Department of Biological Sciences, Purdue University, West Lafayette, IN, USA.
- Department of Forestry and Natural Resources, Purdue University, West Lafayette, IN, USA.
| |
Collapse
|
3
|
Imamoto M, Nakamura H, Aibara M, Hatashima R, Kimirei IA, Kashindye BB, Itoh T, Nikaido M. Severe Bottleneck Impacted the Genomic Structure of Egg-Eating Cichlids in Lake Victoria. Mol Biol Evol 2024; 41:msae093. [PMID: 38782570 PMCID: PMC11166178 DOI: 10.1093/molbev/msae093] [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/04/2023] [Revised: 05/03/2024] [Accepted: 05/06/2024] [Indexed: 05/25/2024] Open
Abstract
Within 15,000 years, the explosive adaptive radiation of haplochromine cichlids in Lake Victoria, East Africa, generated 500 endemic species. In the 1980s, the upsurge of Nile perch, a carnivorous fish artificially introduced to the lake, drove the extinction of more than 200 endemic cichlids. The Nile perch predation particularly harmed piscivorous cichlids, including paedophages, cichlids eat eggs and fries, which is an example of the unique trophic adaptation seen in African cichlids. Here, aiming to investigate past demographic events possibly triggered by the invasion of Nile perch and the subsequent impacts on the genetic structure of cichlids, we conducted large-scale comparative genomics. We discovered evidence of recent bottleneck events in 4 species, including 2 paedophages, which began during the 1970s to 1980s, and population size rebounded during the 1990s to 2000s. The timing of the bottleneck corresponded to the historical records of endemic haplochromines" disappearance and later resurgence, which is likely associated with the introduction of Nile perch by commercial demand to Lake Victoria in the 1950s. Interestingly, among the 4 species that likely experienced bottleneck, Haplochromis sp. "matumbi hunter," a paedophagous cichlid, showed the most severe bottleneck signatures. The components of shared ancestry inferred by ADMIXTURE suggested a high genetic differentiation between matumbi hunter and other species. In contrast, our phylogenetic analyses highly supported the monophyly of the 5 paedophages, consistent with the results of previous studies. We conclude that high genetic differentiation of matumbi hunter occurred due to the loss of shared genetic components among haplochromines in Lake Victoria caused by the recent severe bottleneck.
Collapse
Affiliation(s)
- Minami Imamoto
- Department of Life Science and Technology, School of Life Science and Technology, Tokyo Institute of Technology, Tokyo, Japan
| | - Haruna Nakamura
- Department of Life Science and Technology, School of Life Science and Technology, Tokyo Institute of Technology, Tokyo, Japan
- Research Center for Integrative Evolutionary Science, The Graduate University for Advanced Studies, SOKENDAI, Kanagawa, Japan
| | - Mitsuto Aibara
- Department of Life Science and Technology, School of Life Science and Technology, Tokyo Institute of Technology, Tokyo, Japan
| | - Ryo Hatashima
- Department of Life Science and Technology, School of Life Science and Technology, Tokyo Institute of Technology, Tokyo, Japan
| | - Ismael A Kimirei
- Tanzania Fisheries Research Institute (TAFIRI), Dar es Salaam, Tanzania
| | | | - Takehiko Itoh
- Department of Life Science and Technology, School of Life Science and Technology, Tokyo Institute of Technology, Tokyo, Japan
| | - Masato Nikaido
- Department of Life Science and Technology, School of Life Science and Technology, Tokyo Institute of Technology, Tokyo, Japan
| |
Collapse
|
4
|
Lau SCY, Wilson NG, Golledge NR, Naish TR, Watts PC, Silva CNS, Cooke IR, Allcock AL, Mark FC, Linse K, Strugnell JM. Genomic evidence for West Antarctic Ice Sheet collapse during the Last Interglacial. Science 2023; 382:1384-1389. [PMID: 38127761 DOI: 10.1126/science.ade0664] [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: 07/22/2022] [Accepted: 10/30/2023] [Indexed: 12/23/2023]
Abstract
The marine-based West Antarctic Ice Sheet (WAIS) is considered vulnerable to irreversible collapse under future climate trajectories, and its tipping point may lie within the mitigated warming scenarios of 1.5° to 2°C of the United Nations Paris Agreement. Knowledge of ice loss during similarly warm past climates could resolve this uncertainty, including the Last Interglacial when global sea levels were 5 to 10 meters higher than today and global average temperatures were 0.5° to 1.5°C warmer than preindustrial levels. Using a panel of genome-wide, single-nucleotide polymorphisms of a circum-Antarctic octopus, we show persistent, historic signals of gene flow only possible with complete WAIS collapse. Our results provide the first empirical evidence that the tipping point of WAIS loss could be reached even under stringent climate mitigation scenarios.
Collapse
Affiliation(s)
- Sally C Y Lau
- Centre for Sustainable Tropical Fisheries and Aquaculture and College of Science and Engineering, James Cook University, Townsville, Qld, Australia
- Securing Antarctica's Environmental Future, James Cook University, Townsville, Qld, Australia
| | - Nerida G Wilson
- Collections & Research, Western Australian Museum, Welshpool, WA, Australia
- School of Biological Sciences, University of Western Australia, Perth, WA, Australia
- Securing Antarctica's Environmental Future, Western Australian Museum, Welshpool, WA, Australia
| | - Nicholas R Golledge
- Antarctic Research Centre, Victoria University of Wellington, Wellington, New Zealand
| | - Tim R Naish
- Antarctic Research Centre, Victoria University of Wellington, Wellington, New Zealand
| | - Phillip C Watts
- Department of Biological and Environmental Science, University of Jyväskylä, Jyväskylä, Finland
| | - Catarina N S Silva
- Centre for Sustainable Tropical Fisheries and Aquaculture and College of Science and Engineering, James Cook University, Townsville, Qld, Australia
- Centre for Functional Ecology - Science for People & the Planet (CFE), Associate Laboratory TERRA, Department of Life Sciences, University of Coimbra, Portugal
| | - Ira R Cooke
- Centre for Tropical Bioinformatics and Molecular Biology, James Cook University, Townsville, Qld, Australia
| | - A Louise Allcock
- School of Natural Sciences and Ryan Institute, University of Galway, Galway, Ireland
| | - Felix C Mark
- Alfred Wegener Institute, Helmholtz Centre for Polar and Marine Research, Bremerhaven, Germany
| | | | - Jan M Strugnell
- Centre for Sustainable Tropical Fisheries and Aquaculture and College of Science and Engineering, James Cook University, Townsville, Qld, Australia
- Securing Antarctica's Environmental Future, James Cook University, Townsville, Qld, Australia
| |
Collapse
|
5
|
Qin HT, Mӧller M, Milne R, Luo YH, Zhu GF, Li DZ, Liu J, Gao LM. Multiple paternally inherited chloroplast capture events associated with Taxus speciation in the Hengduan Mountains. Mol Phylogenet Evol 2023; 189:107915. [PMID: 37666379 DOI: 10.1016/j.ympev.2023.107915] [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: 10/27/2022] [Revised: 06/16/2023] [Accepted: 09/01/2023] [Indexed: 09/06/2023]
Abstract
Mountainous regions provide a multitude of habitats and opportunities for complex speciation scenarios. Hybridization leading to chloroplast capture, which can be revealed by incongruent phylogenetic trees, is one possible outcome. Four allopatric Taxus lineages (three species and an undescribed lineage) from the Hengduan Mountains, southwest China, exhibit conflicting phylogenetic relationships between nuclear and chloroplast phylogenies. Here, we use multi-omic data at the population level to investigate their historical speciation processes. Population genomic analysis based on ddRAD-seq data revealed limited contemporary inter-specific gene flow involving only populations located close to another species. In a historical context, chloroplast and nuclear data (transcriptome) consistently showed conflicting phylogenetic relationships for T. florinii and the Emei type lineage. ILS and chloroplast recombination were excluded as possible causes, and transcriptome and ddRAD-seq data revealed an absence of the mosaic nuclear genomes that characterize hybrid origin scenarios. Therefore, T. florinii appears to have originated when a lineage of T. florinii captured the T. chinensis plastid type, whereas plastid introgression in the opposite direction generated the Emei Type. All four species have distinct ecological niche based on community investigations and ecological niche analyses. We propose that the origins of both species represent very rare examples of chloroplast capture events despite the paternal cpDNA inheritance of gymnosperms. Specifically, allopatrically and/or ecologically diverged parental species experienced a rare secondary contact, subsequent hybridization and reciprocal chloroplast capture, generating two new lineages, each of which acquired a unique ecological niche. These events might have been triggered by orogenic activities of the Hengduan Mountains and an intensification of the Asian monsoon in the late Miocene, and may represent a scenario more common in these mountains than presently known.
Collapse
Affiliation(s)
- Han-Tao Qin
- CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming 650201, Yunnan, China; Germplasm Bank of Wild Species, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming 650201, Yunnan, China; University of Chinese Academy of Sciences, Beijing 100049, China
| | - Michael Mӧller
- Royal Botanic Garden Edinburgh, Edinburgh EH3 5LR, United Kingdom
| | - Richard Milne
- Institute of Molecular Plant Sciences, School of Biological Sciences, University of Edinburgh, Edinburgh EH9 3JH, United Kingdom
| | - Ya-Huang Luo
- CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming 650201, Yunnan, China; Lijiang Forest Biodiversity National Observation and Research Station, Kunming Institute of Botany, Chinese Academy of Sciences, Lijiang 674100, Yunnan, China
| | - Guang-Fu Zhu
- CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming 650201, Yunnan, China; Germplasm Bank of Wild Species, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming 650201, Yunnan, China; University of Chinese Academy of Sciences, Beijing 100049, China
| | - De-Zhu Li
- CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming 650201, Yunnan, China; Germplasm Bank of Wild Species, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming 650201, Yunnan, China; University of Chinese Academy of Sciences, Beijing 100049, China; Lijiang Forest Biodiversity National Observation and Research Station, Kunming Institute of Botany, Chinese Academy of Sciences, Lijiang 674100, Yunnan, China.
| | - Jie Liu
- CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming 650201, Yunnan, China; Germplasm Bank of Wild Species, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming 650201, Yunnan, China.
| | - Lian-Ming Gao
- CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming 650201, Yunnan, China; Lijiang Forest Biodiversity National Observation and Research Station, Kunming Institute of Botany, Chinese Academy of Sciences, Lijiang 674100, Yunnan, China.
| |
Collapse
|
6
|
Arana A, Esteves J, Ramírez R, Galetti PM, Pérez Z J, Ramirez JL. Population genomics reveals how 5 ka of human occupancy led the Lima leaf-toed gecko (Phyllodactylus sentosus) to the brink of extinction. Sci Rep 2023; 13:18465. [PMID: 37891335 PMCID: PMC10611785 DOI: 10.1038/s41598-023-45715-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2023] [Accepted: 10/23/2023] [Indexed: 10/29/2023] Open
Abstract
Small species with high home fidelity, high ecological specialization or low vagility are particularly prone to suffer from habitat modification and fragmentation. The Lima leaf-toed gecko (Phyllodactylus sentosus) is a critically endangered Peruvian species that shelters mostly in pre-Incan archeological areas called huacas, where the original environmental conditions are maintained. We used genotyping by sequencing to understand the population genomic history of P. sentosus. We found low genetic diversity (He 0.0406-0.134 and nucleotide diversity 0.0812-0.145) and deviations of the observed heterozygosity relative to the expected heterozygosity in some populations (Fis - 0.0202 to 0.0187). In all analyses, a clear population structuring was observed that cannot be explained by isolation by distance alone. Also, low levels of historical gene flow were observed between most populations, which decreased as shown in contemporary migration rate analysis. Demographic inference suggests these populations experienced bottleneck events during the last 5 ka. These results indicate that habitat modification since pre-Incan civilizations severely affected these populations, which currently face even more drastic urbanization threats. Finally, our predictions show that this species could become extinct in a decade without further intervention, which calls for urgent conservation actions being undertaken.
Collapse
Affiliation(s)
- Alejandra Arana
- Facultad de Ciencias Biológicas, Universidad Nacional Mayor de San Marcos, Lima, Peru
| | - Juan Esteves
- Facultad de Ciencias Biológicas, Universidad Nacional Mayor de San Marcos, Lima, Peru
| | - Rina Ramírez
- Facultad de Ciencias Biológicas, Universidad Nacional Mayor de San Marcos, Lima, Peru
| | - Pedro M Galetti
- Departamento de Genética e Evolução, Universidade Federal de São Carlos, São Carlos, SP, 13565-905, Brazil
| | - José Pérez Z
- Facultad de Ciencias Biológicas, Universidad Nacional Mayor de San Marcos, Lima, Peru
| | - Jorge L Ramirez
- Facultad de Ciencias Biológicas, Universidad Nacional Mayor de San Marcos, Lima, Peru.
| |
Collapse
|
7
|
Genomic and ecological evidence shed light on the recent demographic history of two related invasive insects. Sci Rep 2022; 12:19629. [PMID: 36385480 PMCID: PMC9669014 DOI: 10.1038/s41598-022-21548-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/05/2022] [Accepted: 09/28/2022] [Indexed: 11/17/2022] Open
Abstract
Hypogeococcus pungens is a species complex native to southern South America that is composed of at least five putative species, each one specialized in the use of different host plants. Two of these undescribed species were registered as invasive in Central and North America: Hyp-C is a cactophagous mealybug that became an important pest that threatens endemic cactus species in Puerto Rico, and Hyp-AP feeds on Amaranthaceae and Portulacaceae hosts, but does not produce severe damage to the host plants. We quantified genomic variation and investigated the demographic history of both invasive species by means of coalescent-based simulations using high throughput sequencing data. We also evaluated the incidence of host plant infestation produced by both species and used an ecological niche modeling approach to assess potential distribution under current and future climatic scenarios. Our genetic survey evinced the footprints of strong effective population size reduction and signals of genetic differentiation among populations within each species. Incidence of plant attacks varied between species and among populations within species, with some host plant species preferred over others. Ecological niche modeling suggested that under future climatic scenarios both species would expand their distribution ranges in Puerto Rico. These results provide valuable information for the design of efficient management and control strategies of the Puerto Rican cactus pest and shed light on the evolutionary pathways of biological invasions.
Collapse
|
8
|
Hoey JA, Able KW, Pinsky ML. Genetic decline and recovery of a demographically rebuilt fishery species. Mol Ecol 2022; 31:5684-5698. [PMID: 36114805 PMCID: PMC9828022 DOI: 10.1111/mec.16697] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/12/2020] [Revised: 09/09/2022] [Accepted: 09/15/2022] [Indexed: 01/13/2023]
Abstract
The demographic history of a population is important for conservation and evolution, but this history is unknown for many populations. Methods that use genomic data have been developed to infer demography, but they can be challenging to implement and interpret, particularly for large populations. Thus, understanding if and when genetic estimates of demography correspond to true population history is important for assessing the performance of these genetic methods. Here, we used double-digest restriction-site associated DNA (ddRAD) sequencing data from archived collections of larval summer flounder (Paralichthys dentatus, n = 279) from three cohorts (1994-1995, 1997-1998 and 2008-2009) along the U.S. East coast to examine how contemporary effective population size and genetic diversity responded to changes in abundance in a natural population. Despite little to no detectable change in genetic diversity, coalescent-based demographic modelling from site frequency spectra revealed that summer flounder effective population size declined dramatically in the early 1980s. The timing and direction of change corresponded well with the observed decline in spawning stock census abundance in the late 1980s from independent fish surveys. Census abundance subsequently recovered and achieved the prebottleneck size. Effective population size also grew following the bottleneck. Our results for summer flounder demonstrate that genetic sampling and site frequency spectra can be useful for detecting population dynamics, even in species with large effective sizes.
Collapse
Affiliation(s)
- Jennifer A. Hoey
- Ecology, Evolution, & Natural ResourcesRutgers UniversityNew BrunswickNew JerseyUSA,Institute for Biodiversity Science and SustainabilityCalifornia Academy of SciencesSan FranciscoCaliforniaUSA
| | - Kenneth W. Able
- Marine Field Station, Department of Marine and Coastal Sciences, Rutgers UniversityTuckertonNew JerseyUSA
| | - Malin L. Pinsky
- Ecology, Evolution, & Natural ResourcesRutgers UniversityNew BrunswickNew JerseyUSA
| |
Collapse
|
9
|
Zhang J, Richards ZT, Adam AAS, Chan CX, Shinzato C, Gilmour J, Thomas L, Strugnell JM, Miller DJ, Cooke I. Evolutionary responses of a reef-building coral to climate change at the end of the last glacial maximum. Mol Biol Evol 2022; 39:msac201. [PMID: 36219871 PMCID: PMC9578555 DOI: 10.1093/molbev/msac201] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2022] [Revised: 09/04/2022] [Accepted: 09/13/2022] [Indexed: 11/13/2022] Open
Abstract
Climate change threatens the survival of coral reefs on a global scale, primarily through mass bleaching and mortality as a result of marine heatwaves. While these short-term effects are clear, predicting the fate of coral reefs over the coming century is a major challenge. One way to understand the longer-term effects of rapid climate change is to examine the response of coral populations to past climate shifts. Coastal and shallow-water marine ecosystems such as coral reefs have been reshaped many times by sea-level changes during the Pleistocene, yet, few studies have directly linked this with its consequences on population demographics, dispersal, and adaptation. Here we use powerful analytical techniques, afforded by haplotype phased whole-genomes, to establish such links for the reef-building coral, Acropora digitifera. We show that three genetically distinct populations are present in northwestern Australia, and that their rapid divergence since the last glacial maximum (LGM) can be explained by a combination of founder-effects and restricted gene flow. Signatures of selective sweeps, too strong to be explained by demographic history, are present in all three populations and overlap with genes that show different patterns of functional enrichment between inshore and offshore habitats. In contrast to rapid divergence in the host, we find that photosymbiont communities are largely undifferentiated between corals from all three locations, spanning almost 1000 km, indicating that selection on host genes and not acquisition of novel symbionts, has been the primary driver of adaptation for this species in northwestern Australia.
Collapse
Affiliation(s)
- Jia Zhang
- Department of Molecular and Cell Biology, James Cook University, Townsville, QLD, 4811, Australia
- Centre for Tropical Bioinformatics and Molecular Biology, James Cook University, Townsville, QLD, 4811, Australia
- ARC Centre of Excellence for Coral Reef Studies, James Cook University, Townsville, QLD, 4811, Australia
| | - Zoe T Richards
- Coral Conservation and Research Group, Trace and Environmental DNA Laboratory, School of Molecular and Life Sciences, Curtin University, Bentley, WA 6102, Australia
- Collections and Research, Western Australian Museum, 49 Kew Street Welshpool, WA 6106, Australia
| | - Arne A S Adam
- Coral Conservation and Research Group, Trace and Environmental DNA Laboratory, School of Molecular and Life Sciences, Curtin University, Bentley, WA 6102, Australia
| | - Cheong Xin Chan
- The University of Queensland, School of Chemistry and Molecular Biosciences, Australian Centre for Ecogenomics, Brisbane, QLD 4072, Australia
| | - Chuya Shinzato
- Atmosphere and Ocean Research Institute, The University of Tokyo277-8564, Chiba, Japan
| | - James Gilmour
- Australia Institute of Marine Science, Indian Oceans Marine Research Centre, Crawley, WA, 6009, Australia
| | - Luke Thomas
- Australia Institute of Marine Science, Indian Oceans Marine Research Centre, Crawley, WA, 6009, Australia
- Oceans Graduate School, The UWA Oceans Institute, The University of Western Australia, Perth, WA, 6009, Australia
| | - Jan M Strugnell
- Department of Marine Biology and Aquaculture, James Cook University, Townsville, QLD, 4811, Australia
- Centre for Sustainable Fisheries and Aquaculture, James Cook University, Townsville, QLD, 4811, Australia
| | - David J Miller
- Department of Molecular and Cell Biology, James Cook University, Townsville, QLD, 4811, Australia
- Centre for Tropical Bioinformatics and Molecular Biology, James Cook University, Townsville, QLD, 4811, Australia
- ARC Centre of Excellence for Coral Reef Studies, James Cook University, Townsville, QLD, 4811, Australia
- Marine Climate Change Unit, Okinawa Institute of Science and Technology, Onna-son, Okinawa, Japan 904-0495
| | - Ira Cooke
- Department of Molecular and Cell Biology, James Cook University, Townsville, QLD, 4811, Australia
- Centre for Tropical Bioinformatics and Molecular Biology, James Cook University, Townsville, QLD, 4811, Australia
| |
Collapse
|
10
|
Attwood SW, Hill SC, Aanensen DM, Connor TR, Pybus OG. Phylogenetic and phylodynamic approaches to understanding and combating the early SARS-CoV-2 pandemic. Nat Rev Genet 2022; 23:547-562. [PMID: 35459859 PMCID: PMC9028907 DOI: 10.1038/s41576-022-00483-8] [Citation(s) in RCA: 45] [Impact Index Per Article: 22.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 03/23/2022] [Indexed: 01/05/2023]
Abstract
Determining the transmissibility, prevalence and patterns of movement of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infections is central to our understanding of the impact of the pandemic and to the design of effective control strategies. Phylogenies (evolutionary trees) have provided key insights into the international spread of SARS-CoV-2 and enabled investigation of individual outbreaks and transmission chains in specific settings. Phylodynamic approaches combine evolutionary, demographic and epidemiological concepts and have helped track virus genetic changes, identify emerging variants and inform public health strategy. Here, we review and synthesize studies that illustrate how phylogenetic and phylodynamic techniques were applied during the first year of the pandemic, and summarize their contributions to our understanding of SARS-CoV-2 transmission and control.
Collapse
Affiliation(s)
- Stephen W Attwood
- Department of Zoology, University of Oxford, Oxford, UK.
- Pathogen Genomics Unit, Public Health Wales NHS Trust, Cardiff, UK.
| | - Sarah C Hill
- Department of Pathobiology and Population Sciences, Royal Veterinary College, University of London, London, UK
| | - David M Aanensen
- Centre for Genomic Pathogen Surveillance, Wellcome Genome Campus, Hinxton, UK
- Big Data Institute, Li Ka Shing Centre for Health Information and Discovery, Nuffield Department of Medicine, University of Oxford, Oxford, UK
| | - Thomas R Connor
- Pathogen Genomics Unit, Public Health Wales NHS Trust, Cardiff, UK
- School of Biosciences, Cardiff University, Cardiff, UK
| | - Oliver G Pybus
- Department of Zoology, University of Oxford, Oxford, UK.
- Department of Pathobiology and Population Sciences, Royal Veterinary College, University of London, London, UK.
| |
Collapse
|
11
|
Chak STC, Harris SE, Hultgren KM, Duffy JE, Rubenstein DR. Demographic inference provides insights into the extirpation and ecological dominance of eusocial snapping shrimps. J Hered 2022; 113:552-562. [PMID: 35921239 DOI: 10.1093/jhered/esac035] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/08/2021] [Accepted: 07/27/2022] [Indexed: 11/14/2022] Open
Abstract
Although eusocial animals often achieve ecological dominance in the ecosystems where they occur, many populations are unstable, resulting in local extinction. Both patterns may be linked to the characteristic demography of eusocial species-high reproductive skew and reproductive division of labor support stable effective population sizes that make eusocial groups more competitive in some species, but also lower effective population sizes that increase susceptibility to population collapse in others. Here, we examine the relationship between demography and social organization in Synalpheus snapping shrimps, a group in which eusociality has evolved recently and repeatedly. We show using coalescent demographic modelling that eusocial species have had lower but more stable effective population sizes across 100,000 generations. Our results are consistent with the idea that stable population sizes may enable competitive dominance in eusocial shrimps, but they also suggest that recent population declines are likely caused by eusocial shrimps' heightened sensitivity to environmental changes, perhaps as a result of their low effective population sizes and localized dispersal. Thus, although the unique life histories and demography of eusocial shrimps have likely contributed to their persistence and ecological dominance over evolutionary timescales, these social traits may also make them vulnerable to contemporary environmental change.
Collapse
Affiliation(s)
- Solomon T C Chak
- Department of Ecology, Evolution and Environmental Biology, Columbia University, New York, NY, USA.,Department of Biological Sciences, New Jersey Institute of Technology, Newark, NJ, USA.,Department of Biological Sciences, SUNY College at Old Westbury, Old Westbury, NY, USA
| | - Stephen E Harris
- Department of Ecology, Evolution and Environmental Biology, Columbia University, New York, NY, USA.,Biology Department, SUNY Purchase College, Purchase, NY, USA
| | | | - J Emmett Duffy
- Tennenbaum Marine Observatories Network, Smithsonian Institution, Edgewater, MD, USA
| | - Dustin R Rubenstein
- Department of Ecology, Evolution and Environmental Biology, Columbia University, New York, NY, USA
| |
Collapse
|
12
|
Salloum PM, Santure AW, Lavery SD, de Villemereuil P. Finding the adaptive needles in a population-structured haystack: a case study in a New Zealand mollusc. J Anim Ecol 2022; 91:1209-1221. [PMID: 35318661 PMCID: PMC9311215 DOI: 10.1111/1365-2656.13692] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/04/2021] [Accepted: 03/09/2022] [Indexed: 11/30/2022]
Abstract
Genetic adaptation to future environmental conditions is crucial to help species persist as the climate changes. Genome scans are powerful tools to understand adaptive landscapes, enabling us to correlate genetic diversity with environmental gradients while disentangling neutral from adaptive variation. However, low gene flow can lead to both local adaptation and highly structured populations, and is a major confounding factor for genome scans, resulting in an inflated number of candidate loci. Here, we compared candidate locus detection in a marine mollusc (Onithochiton neglectus), taking advantage of a natural geographical contrast in the levels of genetic structure between its populations. O. neglectus is endemic to New Zealand and distributed throughout an environmental gradient from the subtropical north to the subantarctic south. Due to a brooding developmental mode, populations tend to be locally isolated. However, adult hitchhiking on rafting kelp increases connectivity among southern populations. We applied two genome scans for outliers (Bayescan and PCAdapt) and two genotype–environment association (GEA) tests (BayeScEnv and RDA). To limit issues with false positives, we combined results using the geometric mean of q‐values and performed association tests with random environmental variables. This novel approach is a compromise between stringent and relaxed approaches widely used before, and allowed us to classify candidate loci as low confidence or high confidence. Genome scans for outliers detected a large number of significant outliers in strong and moderately structured populations. No high‐confidence GEA loci were detected in the context of strong population structure. However, 86 high‐confidence loci were associated predominantly with latitudinally varying abiotic factors in the less structured southern populations. This suggests that the degree of connectivity driven by kelp rafting over the southern scale may be insufficient to counteract local adaptation in this species. Our study supports the expectation that genome scans may be prone to errors in highly structured populations. Nonetheless, it also empirically demonstrates that careful statistical controls enable the identification of candidate loci that invite more detailed investigations. Ultimately, genome scans are valuable tools to help guide further research aiming to determine the potential of non‐model species to adapt to future environments.
Collapse
Affiliation(s)
- P M Salloum
- School of Biological Sciences, University of Auckland, Auckland, New Zealand
| | - A W Santure
- School of Biological Sciences, University of Auckland, Auckland, New Zealand
| | - S D Lavery
- School of Biological Sciences, University of Auckland, Auckland, New Zealand.,Institute of Marine Science, Leigh Marine Laboratory, University of Auckland, Warkworth, New Zealand
| | - P de Villemereuil
- Institut de Systématique, Évolution, Biodiversité (ISYEB), École Pratique des Hautes Études
- PSL, MNHN, CNRS, SU, UA, Paris, France
| |
Collapse
|
13
|
Petretto E, Dettori ML, Pazzola M, Manca F, Amills M, Vacca GM. Mitochondrial DNA diversity of the Sardinian local cattle stock. Sci Rep 2022; 12:2486. [PMID: 35169207 PMCID: PMC8847569 DOI: 10.1038/s41598-022-06420-3] [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: 10/01/2021] [Accepted: 01/25/2022] [Indexed: 12/14/2022] Open
Abstract
The aim of this research was to characterize the genetic diversity of the Sarda (Sa, n = 131), Sardo Bruna (SB, n = 44) and Sardo Modicana (SM, n = 26) cattle breeds, reared in the island of Sardinia (Italy). A portion of the mitochondrial DNA hypervariable region was sequenced, in order to identify a potential signature of African introgression. The FST coefficients among populations ranged between 0.056 for Sa vs SB and 0.167 for SB vs SM. AMOVA analysis indicated there was a significant differentiation of the three breeds, although most of diversity was gathered at the within-breed level. The Median Joining Network of the Sardinian sequences showed a potential founder effect signature. A MJ network including Sardinian cattle plus African, Italian, Iberian and Asian sequences, revealed the presence of haplogroup T3, already detected in Sa cattle, and the presence of Hg T1 and Hg T1′2′3, in Sa and SB. The presence of a private haplotype belonging to haplogroup T1, which is characteristic of African taurine breeds, may be due to the introgression of Sardinian breeds with African cattle, either directly (most probable source: North African cattle) or indirectly (through a Mediterranean intermediary already introgressed with African blood).
Collapse
Affiliation(s)
- Elena Petretto
- Department of Veterinary Medicine, University of Sassari, Via Vienna 2, 07100, Sassari, Italy.,Departament de Ciència Animal i dels Aliments, Universitat Autònoma de Barcelona, 08193, Bellaterra, Spain
| | - Maria Luisa Dettori
- Department of Veterinary Medicine, University of Sassari, Via Vienna 2, 07100, Sassari, Italy.
| | - Michele Pazzola
- Department of Veterinary Medicine, University of Sassari, Via Vienna 2, 07100, Sassari, Italy
| | - Fabio Manca
- Department of Veterinary Medicine, University of Sassari, Via Vienna 2, 07100, Sassari, Italy
| | - Marcel Amills
- Departament de Ciència Animal i dels Aliments, Universitat Autònoma de Barcelona, 08193, Bellaterra, Spain.,Centre for Research in Agricultural Genomics (CRAG), CSIC-IRTA-UAB-UB, Campus Universitat Autònoma de Barcelona, 08193, Bellaterra, Spain
| | - Giuseppe Massimo Vacca
- Department of Veterinary Medicine, University of Sassari, Via Vienna 2, 07100, Sassari, Italy
| |
Collapse
|
14
|
Zee MJ, Whiting JR, Paris JR, Bassar RD, Travis J, Weigel D, Reznick DN, Fraser BA. Rapid genomic convergent evolution in experimental populations of Trinidadian guppies (
Poecilia reticulata
). Evol Lett 2022; 6:149-161. [PMID: 35386829 PMCID: PMC8966473 DOI: 10.1002/evl3.272] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/10/2021] [Revised: 12/08/2021] [Accepted: 12/14/2021] [Indexed: 01/14/2023] Open
Affiliation(s)
- Mijke J. Zee
- Biosciences University of Exeter Exeter EX4 4QD United Kingdom
| | | | | | - Ron D. Bassar
- Department of Biology Williams College Williamstown Massachusetts 01267
| | - Joseph Travis
- Department of Biological Science Florida State University Tallahassee Florida 32306
| | - Detlef Weigel
- Department of Molecular Biology Max Planck Institute for Developmental Biology Tübingen 72076 Germany
| | - David N. Reznick
- Department of Biology University of California, Riverside Riverside California 92521
| | | |
Collapse
|
15
|
von Thaden A, Cocchiararo B, Mueller SA, Reiners TE, Reinert K, Tuchscherer I, Janke A, Nowak C. Informing conservation strategies with museum genomics: Long-term effects of past anthropogenic persecution on the elusive European wildcat. Ecol Evol 2021; 11:17932-17951. [PMID: 35003648 PMCID: PMC8717334 DOI: 10.1002/ece3.8385] [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: 06/16/2021] [Revised: 11/04/2021] [Accepted: 11/09/2021] [Indexed: 12/13/2022] Open
Abstract
Like many carnivore species, European wildcats (Felis silvestris) have suffered severe anthropogenic population declines in the past, resulting in a strong population bottleneck at the beginning of the 20th century. In Germany, the species has managed to survive its near extinction in small isolated areas and is currently recolonizing former habitats owing to legal protection and concerted conservation efforts. Here, we SNP-genotyped and mtDNA-sequenced 56 historical and 650 contemporary samples to assess the impact of massive persecution on genetic diversity, population structure, and hybridization dynamics of wildcats. Spatiotemporal analyses suggest that the presumed postglacial differentiation between two genetically distinct metapopulations in Germany is in fact the result of the anthropogenic bottleneck followed by re-expansion from few secluded refugia. We found that, despite the bottleneck, populations experienced no severe genetic erosion, nor suffered from elevated inbreeding or showed signs of increased hybridization with domestic cats. Our findings have significant implications for current wildcat conservation strategies, as the data analyses show that the two presently recognized wildcat population clusters should be treated as a single conservation unit. Although current populations appear under no imminent threat from genetic factors, fostering connectivity through the implementation of forest corridors will facilitate the preservation of genetic diversity and promote long-term viability. The present study documents how museum collections can be used as essential resource for assessing long-term anthropogenic effects on natural populations, for example, regarding population structure and the delineation of appropriate conservation units, potentially informing todays' species conservation.
Collapse
Affiliation(s)
- Alina von Thaden
- Conservation Genetics GroupSenckenberg Research Institute and Natural History Museum FrankfurtGelnhausenGermany
- Institute of Ecology, Evolution & DiversityJohann Wolfgang Goethe‐University, BiologicumFrankfurt am MainGermany
| | - Berardino Cocchiararo
- Conservation Genetics GroupSenckenberg Research Institute and Natural History Museum FrankfurtGelnhausenGermany
- LOEWE Centre for Translational Biodiversity Genomics (LOEWE‐TBG)Frankfurt am MainGermany
| | - Sarah Ashley Mueller
- Conservation Genetics GroupSenckenberg Research Institute and Natural History Museum FrankfurtGelnhausenGermany
- Institute of Ecology, Evolution & DiversityJohann Wolfgang Goethe‐University, BiologicumFrankfurt am MainGermany
| | - Tobias Erik Reiners
- Conservation Genetics GroupSenckenberg Research Institute and Natural History Museum FrankfurtGelnhausenGermany
| | - Katharina Reinert
- Conservation Genetics GroupSenckenberg Research Institute and Natural History Museum FrankfurtGelnhausenGermany
- Department of Physical GeographyJohann Wolfgang Goethe‐UniversityFrankfurt am MainGermany
| | - Iris Tuchscherer
- Conservation Genetics GroupSenckenberg Research Institute and Natural History Museum FrankfurtGelnhausenGermany
- Institute of Ecology, Evolution & DiversityJohann Wolfgang Goethe‐University, BiologicumFrankfurt am MainGermany
| | - Axel Janke
- Institute of Ecology, Evolution & DiversityJohann Wolfgang Goethe‐University, BiologicumFrankfurt am MainGermany
- LOEWE Centre for Translational Biodiversity Genomics (LOEWE‐TBG)Frankfurt am MainGermany
- Senckenberg Biodiversity and Climate Research CentreSenckenberg Gesellschaft für NaturforschungFrankfurt am MainGermany
| | - Carsten Nowak
- Conservation Genetics GroupSenckenberg Research Institute and Natural History Museum FrankfurtGelnhausenGermany
- LOEWE Centre for Translational Biodiversity Genomics (LOEWE‐TBG)Frankfurt am MainGermany
| |
Collapse
|
16
|
Li SH, Liu Y, Yeh CF, Fu Y, Yeung CKL, Lee CC, Chiu CC, Kuo TH, Chan FT, Chen YC, Ko WY, Yao CT. Not out of the woods yet: Signatures of the prolonged negative genetic consequences of a population bottleneck in a rapidly re-expanding wader, the black-faced spoonbill Platalea minor. Mol Ecol 2021; 31:529-545. [PMID: 34726290 DOI: 10.1111/mec.16260] [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: 06/30/2021] [Revised: 09/27/2021] [Accepted: 10/20/2021] [Indexed: 11/30/2022]
Abstract
The long-term persistence of a population which has suffered a bottleneck partly depends on how historical demographic dynamics impacted its genetic diversity and the accumulation of deleterious mutations. Here we provide genomic evidence for the genetic effect of a recent population bottleneck in the endangered black-faced spoonbill (Platalea minor) after its rapid population recovery. Our data suggest that the bird's effective population size, Ne , had been relatively stable (7500-9000) since 22,000 years ago; however, a recent brief yet severe bottleneck (Ne = 20) which we here estimated to occur around the 1940s wiped out >99% of its historical Ne in roughly three generations. Despite a >15-fold population recovery since 1988, we found that black-faced spoonbill population has higher levels of inbreeding (7.4 times more runs of homozygosity) than its sister species, the royal spoonbill (P. regia), which is not thought to have undergone a marked population contraction. Although the two spoonbills have similar levels of genome-wide genetic diversity, our results suggest that selection on more genes was relaxed in the black-faced spoonbill; moreover individual black-faced spoonbills carry more putatively deleterious mutations (Grantham's score > 50), and may therefore express more deleterious phenotypic effects than royal spoonbills. Here we demonstrate the value of using genomic indices to monitor levels of genetic erosion, inbreeding and mutation load in species with conservation concerns. To mitigate the prolonged negative genetic effect of a population bottleneck, we recommend that all possible measures should be employed to maintain population growth of a threatened species.
Collapse
Affiliation(s)
- Shou-Hsien Li
- School of Life Science, National Taiwan Normal University, Taipei, Taiwan
| | - Yang Liu
- State Key Laboratory of Biocontrol, School of Ecology, Sun Yat-Sen University, Guangzhou, China
| | - Chia-Fen Yeh
- School of Life Science, National Taiwan Normal University, Taipei, Taiwan
| | - Yuchen Fu
- School of Life Science, National Taiwan Normal University, Taipei, Taiwan
| | | | - Chun-Cheng Lee
- School of Life Science, National Taiwan Normal University, Taipei, Taiwan
| | - Chi-Cheng Chiu
- School of Life Science, National Taiwan Normal University, Taipei, Taiwan
| | | | - Fang-Tse Chan
- Division of Zoology, Taiwan Endemic Species Research Institute, Nantou, Taiwan
| | - Yu-Chia Chen
- Department of Life Sciences, National Yanming Medical University, Taipei, Taiwan
| | - Wen-Ya Ko
- Department of Life Sciences, National Yanming Medical University, Taipei, Taiwan
| | - Cheng-Te Yao
- High Altitude Research Station, Taiwan Endemic Species Research Institute, Nantou, Taiwan
| |
Collapse
|
17
|
Trejo-Salazar RE, Castellanos-Morales G, Hernández-Rosales D, Gámez N, Gasca-Pineda J, Morales Garza MR, Medellin R, Eguiarte LE. Discordance in maternal and paternal genetic markers in lesser long-nosed bat Leptonycteris yerbabuenae, a migratory bat: recent expansion to the North and male phylopatry. PeerJ 2021; 9:e12168. [PMID: 34703665 PMCID: PMC8487242 DOI: 10.7717/peerj.12168] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2021] [Accepted: 08/26/2021] [Indexed: 11/20/2022] Open
Abstract
Leptonycteris yerbabuenae, the lesser long-nosed bat is an abundant migratory nectar-feeding bat found in most of Mexico, and in some areas of northern Central America and small sections of southwestern USA. We analyzed the distribution of the maternal and paternal lineages of this species with phylogeographic methods based on two mitochondrial markers, Cyt-b and D-loop, and a marker located in the Y chromosome, DBY. We obtained tissue samples from 220 individuals from 23 localities. Levels of genetic diversity (haplotype diversity, Hd ) were high (Cyt-b = 0.757; D-loop = 0.8082; DBY = 0.9137). No clear patterns of population genetic structure were found for mitochondrial markers, while male genetic differentiation suggested the presence of two lineages: one from Mexican Pacific coast states and another from central-southern Mexico; in accordance to strong male philopatry and higher female migration. We used genealogical reconstructions based on Bayesian tools to calculate divergence times, and to test coalescent models to explain changes in L. yerbabuenae historical demography. Our results show that recent demographic changes were consistent with global climatic changes (∼130,000 kyr ago for Cyt-b and ∼160,000 kyr for D-loop) and divergence times dated from molecular genealogies exhibited older divergence times, Cyt-b (4.03 mya), D-loop (10.26 mya) and DBY (12.23 mya). Accordingly, the female lineage underwent demographic expansion associated to Pleistocene climate change, whereas the male lineage remained constant.
Collapse
Affiliation(s)
- Roberto-Emiliano Trejo-Salazar
- Pograma de Doctorado en Ciencias Biomédicas, Instituto de Ecología, Universidad Nacional Autónoma de México, Ciudad Universitaria, Ciudad de Mexico, México
- Ecología Evolutiva, Instituto de Ecología, Universidad Nacional Autónoma de México, Ciudad Universitaria, Ciudad de México, México
| | | | - DulceCarolina Hernández-Rosales
- Ecología Evolutiva, Instituto de Ecología, Universidad Nacional Autónoma de México, Ciudad Universitaria, Ciudad de México, México
| | - Niza Gámez
- Facultad de Estudios Superiores Zaragoza, Universidad Nacional Autónoma de México, Ciudad de Mexico, Ciudad de Mexico, Mexico
| | - Jaime Gasca-Pineda
- Ecología Evolutiva, Instituto de Ecología, Universidad Nacional Autónoma de México, Ciudad Universitaria, Ciudad de México, México
| | - Miguel Rene Morales Garza
- Facultad de Ciencia y Tecnología, Universidad Simón Bolívar, Ciudad de Mexico, Ciudad de Mexico, Mexico
| | - Rodrigo Medellin
- Ecología de la Biodiversidad, Instituto de Ecología, Universidad Nacional Autónoma de México, Ciudad de Mexico, Ciudad de Mexico, Mexico
| | - Luis E. Eguiarte
- Ecología Evolutiva, Instituto de Ecología, Universidad Nacional Autónoma de México, Ciudad Universitaria, Ciudad de México, México
| |
Collapse
|
18
|
Berger DJ, Crellen T, Lamberton PHL, Allan F, Tracey A, Noonan JD, Kabatereine NB, Tukahebwa EM, Adriko M, Holroyd N, Webster JP, Berriman M, Cotton JA. Whole-genome sequencing of Schistosoma mansoni reveals extensive diversity with limited selection despite mass drug administration. Nat Commun 2021; 12:4776. [PMID: 34362894 PMCID: PMC8346512 DOI: 10.1038/s41467-021-24958-0] [Citation(s) in RCA: 23] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/06/2020] [Accepted: 07/06/2021] [Indexed: 02/07/2023] Open
Abstract
Control and elimination of the parasitic disease schistosomiasis relies on mass administration of praziquantel. Whilst these programmes reduce infection prevalence and intensity, their impact on parasite transmission and evolution is poorly understood. Here we examine the genomic impact of repeated mass drug administration on Schistosoma mansoni populations with documented reduced praziquantel efficacy. We sequenced whole-genomes of 198 S. mansoni larvae from 34 Ugandan children from regions with contrasting praziquantel exposure. Parasites infecting children from Lake Victoria, a transmission hotspot, form a diverse panmictic population. A single round of treatment did not reduce this diversity with no apparent population contraction caused by long-term praziquantel use. We find evidence of positive selection acting on members of gene families previously implicated in praziquantel action, but detect no high frequency functionally impactful variants. As efforts to eliminate schistosomiasis intensify, our study provides a foundation for genomic surveillance of this major human parasite.
Collapse
Affiliation(s)
- Duncan J Berger
- Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton, Cambridgeshire, UK.
- Department of Pathology and Pathogen Biology, Centre for Emerging, Endemic and Exotic Diseases, Royal Veterinary College, University of London, Herts, UK.
| | - Thomas Crellen
- Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton, Cambridgeshire, UK
- Imperial College London, Department of Infectious Disease Epidemiology, London, UK
- Big Data Institute, Nuffield Department of Medicine, University of Oxford, Oxford, UK
| | - Poppy H L Lamberton
- Imperial College London, Department of Infectious Disease Epidemiology, London, UK
- Institute for Biodiversity, Animal Health, and Comparative Medicine, and Wellcome Centre for Integrative Parasitology, University of Glasgow, Glasgow, UK
| | - Fiona Allan
- The Natural History Museum, Department of Life Sciences, London, UK
| | - Alan Tracey
- Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton, Cambridgeshire, UK
| | - Jennifer D Noonan
- Institute of Parasitology, Faculty of Agricultural and Environmental Sciences, McGill University, Montreal, Quebec, Canada
| | - Narcis B Kabatereine
- Vector Borne & Neglected Tropical Disease Control Division, Ministry of Health, Kampala, Uganda
| | - Edridah M Tukahebwa
- Vector Borne & Neglected Tropical Disease Control Division, Ministry of Health, Kampala, Uganda
| | - Moses Adriko
- Vector Borne & Neglected Tropical Disease Control Division, Ministry of Health, Kampala, Uganda
| | - Nancy Holroyd
- Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton, Cambridgeshire, UK
| | - Joanne P Webster
- Department of Pathology and Pathogen Biology, Centre for Emerging, Endemic and Exotic Diseases, Royal Veterinary College, University of London, Herts, UK.
- Imperial College London, Department of Infectious Disease Epidemiology, London, UK.
| | - Matthew Berriman
- Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton, Cambridgeshire, UK.
| | - James A Cotton
- Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton, Cambridgeshire, UK.
| |
Collapse
|
19
|
Walton W, Stone GN, Lohse K. Discordant Pleistocene population size histories in a guild of hymenopteran parasitoids. Mol Ecol 2021; 30:4538-4550. [PMID: 34252238 DOI: 10.1111/mec.16074] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/23/2019] [Revised: 06/23/2021] [Accepted: 07/06/2021] [Indexed: 01/03/2023]
Abstract
Signatures of past changes in population size have been detected in genome-wide variation in many species. However, the causes of such demographic changes and the extent to which they are shared across co-distributed species remain poorly understood. During Pleistocene glacial maxima, many temperate European species were confined to southern refugia. While vicariance and range expansion processes associated with glacial cycles have been widely documented, it is unclear whether refugial populations of co-distributed species have experienced shared histories of population size change. We analyse whole-genome sequence data to reconstruct and compare demographic histories during the Quaternary for Iberian refuge populations in a single ecological guild (seven species of chalcid parasitoid wasps associated with oak cynipid galls). For four of these species, we find support for large changes in effective population size (Ne ) through the Pleistocene that coincide with major climate events. However, there is little evidence that the timing, direction and magnitude of demographic change are shared across species, suggesting that demographic histories in this guild are largely idiosyncratic, even at the scale of a single glacial refugium.
Collapse
Affiliation(s)
- William Walton
- Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, UK
| | - Graham N Stone
- Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, UK
| | - Konrad Lohse
- Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, UK
| |
Collapse
|
20
|
Genetic structure and population history in two critically endangered Kaua‘i honeycreepers. CONSERV GENET 2021. [DOI: 10.1007/s10592-021-01382-x] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
|
21
|
North HL, McGaughran A, Jiggins CD. Insights into invasive species from whole-genome resequencing. Mol Ecol 2021; 30:6289-6308. [PMID: 34041794 DOI: 10.1111/mec.15999] [Citation(s) in RCA: 31] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2020] [Revised: 03/12/2021] [Accepted: 04/30/2021] [Indexed: 12/12/2022]
Abstract
Studies of invasive species can simultaneously inform management strategies and quantify rapid evolution in the wild. The role of genomics in invasion science is increasingly recognised, and the growing availability of reference genomes for invasive species is paving the way for whole-genome resequencing studies in a wide range of systems. Here, we survey the literature to assess the application of whole-genome resequencing data in invasion biology. For some applications, such as the reconstruction of invasion routes in time and space, sequencing the whole genome of many individuals can increase the accuracy of existing methods. In other cases, population genomic approaches such as haplotype analysis can permit entirely new questions to be addressed and new technologies applied. To date whole-genome resequencing has only been used in a handful of invasive systems, but these studies have confirmed the importance of processes such as balancing selection and hybridization in allowing invasive species to reuse existing adaptations and rapidly overcome the challenges of a foreign ecosystem. The use of genomic data does not constitute a paradigm shift per se, but by leveraging new theory, tools, and technologies, population genomics can provide unprecedented insight into basic and applied aspects of invasion science.
Collapse
Affiliation(s)
- Henry L North
- Department of Zoology, University of Cambridge, Cambridge, UK
| | - Angela McGaughran
- Te Aka Mātuatua/School of Science, University of Waikato, Hamilton, New Zealand
| | - Chris D Jiggins
- Department of Zoology, University of Cambridge, Cambridge, UK
| |
Collapse
|
22
|
Sanderson BJ, Feng G, Hu N, Carlson CH, Smart LB, Keefover-Ring K, Yin T, Ma T, Liu J, DiFazio SP, Olson MS. Sex determination through X-Y heterogamety in Salix nigra. Heredity (Edinb) 2021; 126:630-639. [PMID: 33510464 PMCID: PMC8115673 DOI: 10.1038/s41437-020-00397-3] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/21/2020] [Revised: 12/01/2020] [Accepted: 12/14/2020] [Indexed: 01/30/2023] Open
Abstract
The development of non-recombining sex chromosomes has radical effects on the evolution of discrete sexes and sexual dimorphism. Although dioecy is rare in plants, sex chromosomes have evolved repeatedly throughout the diversification of angiosperms, and many of these sex chromosomes are relatively young compared to those found in vertebrates. In this study, we designed and used a sequence capture array to identify a novel sex-linked region (SLR) in Salix nigra, a basal species in the willow clade, and demonstrated that this species has XY heterogamety. We did not detect any genetic overlap with the previously characterized ZW SLRs in willows, which map to a different chromosome. The S. nigra SLR is characterized by strong recombination suppression across a 2 MB region and an excess of low-frequency alleles, resulting in a low Tajima's D compared to the remainder of the genome. We speculate that either a recent bottleneck in population size or factors related to positive or background selection generated this differential pattern of Tajima's D on the X and autosomes. This discovery provides insights into factors that may influence the evolution of sex chromosomes in plants and contributes to a large number of recent observations that underscore their dynamic nature.
Collapse
Affiliation(s)
- Brian J. Sanderson
- grid.264784.b0000 0001 2186 7496Department of Biological Sciences, Texas Tech University, Lubbock, TX 79409-3131 USA ,grid.268154.c0000 0001 2156 6140Present Address: Department of Biology, West Virginia University, Morgantown, WV 26506-6057 USA
| | - Guanqiao Feng
- grid.264784.b0000 0001 2186 7496Department of Biological Sciences, Texas Tech University, Lubbock, TX 79409-3131 USA
| | - Nan Hu
- grid.264784.b0000 0001 2186 7496Department of Biological Sciences, Texas Tech University, Lubbock, TX 79409-3131 USA
| | - Craig H. Carlson
- grid.5386.8000000041936877XHorticulture Section, School of Integrative Plant Science, Cornell University, Cornell AgriTech, Geneva, NY 14456 USA
| | - Lawrence B. Smart
- grid.5386.8000000041936877XHorticulture Section, School of Integrative Plant Science, Cornell University, Cornell AgriTech, Geneva, NY 14456 USA
| | - Ken Keefover-Ring
- grid.14003.360000 0001 2167 3675Departments of Botany and Geography, University of Wisconsin-Madison, Madison, WI 53706 USA
| | - Tongming Yin
- grid.410625.40000 0001 2293 4910Key Laboratory of Tree Genetics and Biotechnology of Jiangsu Province and Education Department of China, Nanjing Forestry University, Nanjing, China
| | - Tao Ma
- grid.13291.380000 0001 0807 1581Key Laboratory of Bio-Resource and Eco-Environment of Ministry of Education & College of Life Sciences, Sichuan University, Chengdu, 610065 China
| | - Jianquan Liu
- grid.13291.380000 0001 0807 1581Key Laboratory of Bio-Resource and Eco-Environment of Ministry of Education & College of Life Sciences, Sichuan University, Chengdu, 610065 China ,grid.32566.340000 0000 8571 0482State Key Laboratory of Grassland Agro-Ecosystem, Institute of Innovation Ecology & College of Life Sciences, Lanzhou University, Lanzhou, 730000 China
| | - Stephen P. DiFazio
- grid.268154.c0000 0001 2156 6140Department of Biology, West Virginia University, Morgantown, WV 26506-6057 USA
| | - Matthew S. Olson
- grid.264784.b0000 0001 2186 7496Department of Biological Sciences, Texas Tech University, Lubbock, TX 79409-3131 USA
| |
Collapse
|
23
|
Schmidt TL, Swan T, Chung J, Karl S, Demok S, Yang Q, Field MA, Muzari MO, Ehlers G, Brugh M, Bellwood R, Horne P, Burkot TR, Ritchie S, Hoffmann AA. Spatial population genomics of a recent mosquito invasion. Mol Ecol 2021; 30:1174-1189. [PMID: 33421231 DOI: 10.1111/mec.15792] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2020] [Revised: 12/20/2020] [Accepted: 01/04/2021] [Indexed: 02/06/2023]
Abstract
Population genomic approaches can characterize dispersal across a single generation through to many generations in the past, bridging the gap between individual movement and intergenerational gene flow. These approaches are particularly useful when investigating dispersal in recently altered systems, where they provide a way of inferring long-distance dispersal between newly established populations and their interactions with existing populations. Human-mediated biological invasions represent such altered systems which can be investigated with appropriate study designs and analyses. Here we apply temporally restricted sampling and a range of population genomic approaches to investigate dispersal in a 2004 invasion of Aedes albopictus (the Asian tiger mosquito) in the Torres Strait Islands (TSI) of Australia. We sampled mosquitoes from 13 TSI villages simultaneously and genotyped 373 mosquitoes at genome-wide single nucleotide polymorphisms (SNPs): 331 from the TSI, 36 from Papua New Guinea (PNG) and four incursive mosquitoes detected in uninvaded regions. Within villages, spatial genetic structure varied substantially but overall displayed isolation by distance and a neighbourhood size of 232-577. Close kin dyads revealed recent movement between islands 31-203 km apart, and deep learning inferences showed incursive Ae. albopictus had travelled to uninvaded regions from both adjacent and nonadjacent islands. Private alleles and a co-ancestry matrix indicated direct gene flow from PNG into nearby islands. Outlier analyses also detected four linked alleles introgressed from PNG, with the alleles surrounding 12 resistance-associated cytochrome P450 genes. By treating dispersal as both an intergenerational process and a set of discrete events, we describe a highly interconnected invasive system.
Collapse
Affiliation(s)
- Thomas L Schmidt
- Pest and Environmental Adaptation Research Group, School of BioSciences, Bio21 Institute, University of Melbourne, Parkville, VIC, Australia
| | - Tom Swan
- Australian Institute of Tropical Health and Medicine, James Cook University, Cairns, QLD, Australia.,College of Public Health, Medical and Veterinary Sciences, James Cook University, Cairns, QLD, Australia
| | - Jessica Chung
- Pest and Environmental Adaptation Research Group, School of BioSciences, Bio21 Institute, University of Melbourne, Parkville, VIC, Australia.,Melbourne Bioinformatics, University of Melbourne, Parkville, VIC, Australia
| | - Stephan Karl
- Australian Institute of Tropical Health and Medicine, James Cook University, Cairns, QLD, Australia.,Vector Borne Diseases Unit, Papua New Guinea Institute of Medical Research, Madang, Papua New Guinea
| | - Samuel Demok
- Vector Borne Diseases Unit, Papua New Guinea Institute of Medical Research, Madang, Papua New Guinea
| | - Qiong Yang
- Pest and Environmental Adaptation Research Group, School of BioSciences, Bio21 Institute, University of Melbourne, Parkville, VIC, Australia
| | - Matt A Field
- Australian Institute of Tropical Health and Medicine, James Cook University, Cairns, QLD, Australia.,John Curtin School of Medical Research, Australian National University, Canberra, ACT, Australia
| | - Mutizwa Odwell Muzari
- Medical Entomology, Tropical Public Health Services Cairns, Cairns and Hinterland Hospital & Health Services, Cairns, QLD, Australia
| | - Gerhard Ehlers
- Medical Entomology, Tropical Public Health Services Cairns, Cairns and Hinterland Hospital & Health Services, Cairns, QLD, Australia
| | - Mathew Brugh
- Medical Entomology, Tropical Public Health Services Cairns, Cairns and Hinterland Hospital & Health Services, Cairns, QLD, Australia
| | - Rodney Bellwood
- Medical Entomology, Tropical Public Health Services Cairns, Cairns and Hinterland Hospital & Health Services, Cairns, QLD, Australia
| | - Peter Horne
- Medical Entomology, Tropical Public Health Services Cairns, Cairns and Hinterland Hospital & Health Services, Cairns, QLD, Australia
| | - Thomas R Burkot
- Australian Institute of Tropical Health and Medicine, James Cook University, Cairns, QLD, Australia
| | - Scott Ritchie
- College of Public Health, Medical and Veterinary Sciences, James Cook University, Cairns, QLD, Australia.,Institute of Vector-Borne Disease, Monash University, Clayton, VIC, Australia
| | - Ary A Hoffmann
- Pest and Environmental Adaptation Research Group, School of BioSciences, Bio21 Institute, University of Melbourne, Parkville, VIC, Australia
| |
Collapse
|
24
|
Klingler KB, Jahner JP, Parchman TL, Ray C, Peacock MM. Genomic variation in the American pika: signatures of geographic isolation and implications for conservation. BMC Ecol Evol 2021; 21:2. [PMID: 33514306 PMCID: PMC7853312 DOI: 10.1186/s12862-020-01739-9] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2020] [Accepted: 12/17/2020] [Indexed: 01/12/2023] Open
Abstract
Background Distributional responses by alpine taxa to repeated, glacial-interglacial cycles throughout the last two million years have significantly influenced the spatial genetic structure of populations. These effects have been exacerbated for the American pika (Ochotona princeps), a small alpine lagomorph constrained by thermal sensitivity and a limited dispersal capacity. As a species of conservation concern, long-term lack of gene flow has important consequences for landscape genetic structure and levels of diversity within populations. Here, we use reduced representation sequencing (ddRADseq) to provide a genome-wide perspective on patterns of genetic variation across pika populations representing distinct subspecies. To investigate how landscape and environmental features shape genetic variation, we collected genetic samples from distinct geographic regions as well as across finer spatial scales in two geographically proximate mountain ranges of eastern Nevada. Results Our genome-wide analyses corroborate range-wide, mitochondrial subspecific designations and reveal pronounced fine-scale population structure between the Ruby Mountains and East Humboldt Range of eastern Nevada. Populations in Nevada were characterized by low genetic diversity (π = 0.0006–0.0009; θW = 0.0005–0.0007) relative to populations in California (π = 0.0014–0.0019; θW = 0.0011–0.0017) and the Rocky Mountains (π = 0.0025–0.0027; θW = 0.0021–0.0024), indicating substantial genetic drift in these isolated populations. Tajima’s D was positive for all sites (D = 0.240–0.811), consistent with recent contraction in population sizes range-wide. Conclusions Substantial influences of geography, elevation and climate variables on genetic differentiation were also detected and may interact with the regional effects of anthropogenic climate change to force the loss of unique genetic lineages through continued population extirpations in the Great Basin and Sierra Nevada.
Collapse
Affiliation(s)
| | - Joshua P Jahner
- Department of Biology, University of Nevada, Reno, 89557, USA
| | - Thomas L Parchman
- Department of Biology, University of Nevada, Reno, 89557, USA.,Program in Ecology, Evolution, and Conservation Biology, University of Nevada, Reno, NV, 89557, USA
| | - Chris Ray
- Department of Ecology and Evolutionary Biology, University of Colorado, Boulder, CO, 80309-0334, USA
| | - Mary M Peacock
- Department of Biology, University of Nevada, Reno, 89557, USA. .,Program in Ecology, Evolution, and Conservation Biology, University of Nevada, Reno, NV, 89557, USA.
| |
Collapse
|
25
|
Jay F, Boitard S, Austerlitz F. An ABC Method for Whole-Genome Sequence Data: Inferring Paleolithic and Neolithic Human Expansions. Mol Biol Evol 2020; 36:1565-1579. [PMID: 30785202 DOI: 10.1093/molbev/msz038] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/24/2022] Open
Abstract
Species generally undergo a complex demographic history consisting, in particular, of multiple changes in population size. Genome-wide sequencing data are potentially highly informative for reconstructing this demographic history. A crucial point is to extract the relevant information from these very large data sets. Here, we design an approach for inferring past demographic events from a moderate number of fully sequenced genomes. Our new approach uses Approximate Bayesian Computation, a simulation-based statistical framework that allows 1) identifying the best demographic scenario among several competing scenarios and 2) estimating the best-fitting parameters under the chosen scenario. Approximate Bayesian Computation relies on the computation of summary statistics. Using a cross-validation approach, we show that statistics such as the lengths of haplotypes shared between individuals, or the decay of linkage disequilibrium with distance, can be combined with classical statistics (e.g., heterozygosity and Tajima's D) to accurately infer complex demographic scenarios including bottlenecks and expansion periods. We also demonstrate the importance of simultaneously estimating the genotyping error rate. Applying our method on genome-wide human-sequence databases, we finally show that a model consisting in a bottleneck followed by a Paleolithic and a Neolithic expansion is the most relevant for Eurasian populations.
Collapse
Affiliation(s)
- Flora Jay
- Laboratoire EcoAnthropologie et Ethnobiologie, CNRS/MNHN/Université Paris Diderot, Paris, France.,Laboratoire de Recherche en Informatique, CNRS/Université Paris-Sud/Université Paris-Saclay, Orsay, France
| | - Simon Boitard
- GenPhySE, Université de Toulouse, INRA, INPT, INP-ENVT, Castanet Tolosan, France
| | - Frédéric Austerlitz
- Laboratoire EcoAnthropologie et Ethnobiologie, CNRS/MNHN/Université Paris Diderot, Paris, France
| |
Collapse
|
26
|
Chattopadhyay B, Garg KM, Ray R, Rheindt FE. Fluctuating fortunes: genomes and habitat reconstructions reveal global climate-mediated changes in bats' genetic diversity. Proc Biol Sci 2019; 286:20190304. [PMID: 31530139 PMCID: PMC6784725 DOI: 10.1098/rspb.2019.0304] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2019] [Accepted: 08/23/2019] [Indexed: 12/21/2022] Open
Abstract
Over the last approximately 2.6 Myr, Earth's climate has been dominated by cyclical ice ages that have profoundly affected species' population sizes, but the impact of impending anthropogenic climate change on species' extinction potential remains a worrying problem. We investigated 11 bat species from different taxonomic, ecological and geographical backgrounds using combined information from palaeoclimatic habitat reconstructions and genomes to analyse biotic impacts of historic climate change. We discover tightly correlated fluctuations between species' historic distribution and effective population size, identify frugivores as particularly susceptible to global warming, pinpoint large insectivores as having overall low effective population size and flag the onset of the Holocene (approx. 10-12 000 years ago) as the period with the generally lowest effective population sizes across the last approximately 1 Myr. Our study shows that combining genomic and palaeoclimatological approaches reveals effects of climatic shifts on genetic diversity and may help predict impacts of future climate change.
Collapse
Affiliation(s)
| | - Kritika M. Garg
- Department of Biological Sciences, National University of Singapore, Singapore
| | - Rajasri Ray
- Center for Ecological Sciences, Indian Institute of Science, Bangalore, 560012 Karnataka, India
- Centre for Studies in Ethnobiology, Biodiversity and Sustainability (CEiBa), BG Road, Mokdumpur, Malda-732103 West Bengal, India
| | - Frank E. Rheindt
- Department of Biological Sciences, National University of Singapore, Singapore
| |
Collapse
|
27
|
Linck E, Battey CJ. Minor allele frequency thresholds strongly affect population structure inference with genomic data sets. Mol Ecol Resour 2019; 19:639-647. [DOI: 10.1111/1755-0998.12995] [Citation(s) in RCA: 184] [Impact Index Per Article: 36.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/23/2018] [Revised: 12/14/2018] [Accepted: 01/04/2019] [Indexed: 01/25/2023]
Affiliation(s)
- Ethan Linck
- Department of Biology and Burke Museum of Natural History and Culture, University of Washington Seattle Washington
| | - C. J. Battey
- Department of Biology and Institute of Ecology and Evolution University of Oregon Eugene Oregon
| |
Collapse
|
28
|
Population Genomics and Structure of the Critically Endangered Mariana Crow ( Corvus kubaryi). Genes (Basel) 2019; 10:genes10030187. [PMID: 30832245 PMCID: PMC6471520 DOI: 10.3390/genes10030187] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/07/2018] [Revised: 02/16/2019] [Accepted: 02/19/2019] [Indexed: 11/29/2022] Open
Abstract
The Mariana Crow, or Åga (Corvus kubaryi), is a critically endangered species (IUCN -International Union for Conservation of Nature), endemic to the islands of Guam and Rota in the Mariana Archipelago. It is locally extinct on Guam, and numbers have declined dramatically on Rota to a historical low of less than 55 breeding pairs throughout the island in 2013. Because of its extirpation on Guam and population decline on Rota, it is of critical importance to assess the genetic variation among individuals to assist ongoing recovery efforts. We conducted a population genomics analysis comparing the Guam and Rota populations and studied the genetic structure of the Rota population. We used blood samples from five birds from Guam and 78 birds from Rota. We identified 145,552 candidate single nucleotide variants (SNVs) from a genome sequence of an individual from Rota and selected a subset of these to develop an oligonucleotide in-solution capture assay. The Guam and Rota populations were genetically differentiated from each other. Crow populations sampled broadly across their range on Rota showed significant genetic structuring – a surprising result given the small size of this island and the good flight capabilities of the species. Knowledge of its genetic structure will help improve management strategies to help with its recovery.
Collapse
|
29
|
Sainudiin R, Véber A. Full likelihood inference from the site frequency spectrum based on the optimal tree resolution. Theor Popul Biol 2018; 124:1-15. [DOI: 10.1016/j.tpb.2018.07.002] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2017] [Revised: 06/11/2018] [Accepted: 07/09/2018] [Indexed: 10/28/2022]
|
30
|
Xue AT, Ruggiero RP, Hickerson MJ, Boissinot S. Differential Effect of Selection against LINE Retrotransposons among Vertebrates Inferred from Whole-Genome Data and Demographic Modeling. Genome Biol Evol 2018; 10:1265-1281. [PMID: 29688421 PMCID: PMC5963298 DOI: 10.1093/gbe/evy083] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 04/20/2018] [Indexed: 12/30/2022] Open
Abstract
Variation in LINE composition is one of the major determinants for the substantial size and structural differences among vertebrate genomes. In particular, the larger genomes of mammals are characterized by hundreds of thousands of copies from a single LINE clade, L1, whereas nonmammalian vertebrates possess a much greater diversity of LINEs, yet with orders of magnitude less in copy number. It has been proposed that such variation in copy number among vertebrates is due to differential effect of LINE insertions on host fitness. To investigate LINE selection, we deployed a framework of demographic modeling, coalescent simulations, and probabilistic inference against population-level whole-genome data sets for four model species: one population each of threespine stickleback, green anole, and house mouse, as well as three human populations. Specifically, we inferred a null demographic background utilizing SNP data, which was then exploited to simulate a putative null distribution of summary statistics that was compared with LINE data. Subsequently, we applied the inferred null demographic model with an additional exponential size change parameter, coupled with model selection, to test for neutrality as well as estimate the strength of either negative or positive selection. We found a robust signal for purifying selection in anole and mouse, but a lack of clear evidence for selection in stickleback and human. Overall, we demonstrated LINE insertion dynamics that are not in accordance to a mammalian versus nonmammalian dichotomy, and instead the degree of existing LINE activity together with host-specific demographic history may be the main determinants of LINE abundance.
Collapse
Affiliation(s)
- Alexander T Xue
- Department of Biology: Subprogram in Ecology, Evolutionary Biology, and Behavior, City College and Graduate Center of City University of New York
- Human Genetics Institute of New Jersey and Department of Genetics, Rutgers University, Piscataway
| | - Robert P Ruggiero
- New York University Abu Dhabi, Saadiyat Island Campus, United Arab Emirates
| | - Michael J Hickerson
- Department of Biology: Subprogram in Ecology, Evolutionary Biology, and Behavior, City College and Graduate Center of City University of New York
- Division of Invertebrate Zoology, American Museum of Natural History, New York, New York
| | - Stéphane Boissinot
- New York University Abu Dhabi, Saadiyat Island Campus, United Arab Emirates
| |
Collapse
|
31
|
Bennett KL, Kaddumukasa M, Shija F, Djouaka R, Misinzo G, Lutwama J, Linton YM, Walton C. Comparative phylogeography of Aedes mosquitoes and the role of past climatic change for evolution within Africa. Ecol Evol 2018; 8:3019-3036. [PMID: 29531714 PMCID: PMC5838080 DOI: 10.1002/ece3.3668] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/19/2017] [Revised: 10/26/2017] [Accepted: 10/27/2017] [Indexed: 01/01/2023] Open
Abstract
The study of demographic processes involved in species diversification and evolution ultimately provides explanations for the complex distribution of biodiversity on earth, indicates regions important for the maintenance and generation of biodiversity, and identifies biological units important for conservation or medical consequence. African and forest biota have both received relatively little attention with regard to understanding their diversification, although one possible mechanism is that this has been driven by historical climate change. To investigate this, we implemented a standard population genetics approach along with Approximate Bayesian Computation, using sequence data from two exon-primed intron-crossing (EPIC) nuclear loci and mitochondrial cytochrome oxidase subunit I, to investigate the evolutionary history of five medically important and inherently forest dependent mosquito species of the genus Aedes. By testing different demographic hypotheses, we show that Aedes bromeliae and Aedes lilii fit the same model of lineage diversification, admixture, expansion, and recent population structure previously inferred for Aedes aegypti. In addition, analyses of population structure show that Aedes africanus has undergone lineage diversification and expansion while Aedes hansfordi has been impacted by population expansion within Uganda. This congruence in evolutionary history is likely to relate to historical climate-driven habitat change within Africa during the late Pleistocene and Holocene epoch. We find differences in the population structure of mosquitoes from Tanzania and Uganda compared to Benin and Uganda which could relate to differences in the historical connectivity of forests across the continent. Our findings emphasize the importance of recent climate change in the evolution of African forest biota.
Collapse
Affiliation(s)
- Kelly Louise Bennett
- Faculty of Life SciencesComputational Evolutionary Biology GroupUniversity of ManchesterManchesterUK
| | - Martha Kaddumukasa
- Department of Arbovirology, Emerging and Re‐emerging InfectionsUganda Virus Research InstituteEntebbeUganda
- WITS Institute for Malaria ResearchSchool of Pathology Faculty of Health SciencesUniversity of WitwatersrandParktownJohannesburg
| | - Fortunate Shija
- Faculty of Life SciencesComputational Evolutionary Biology GroupUniversity of ManchesterManchesterUK
- Department of Veterinary Microbiology and ParasitologySokoine University of AgricultureMorogoroTanzania
| | - Rousseau Djouaka
- Agro‐Eco‐Health Platform for West and Central AfricaInternational Institute for Tropical AgricultureCotonouRepublic of Benin
| | - Gerald Misinzo
- Agro‐Eco‐Health Platform for West and Central AfricaInternational Institute for Tropical AgricultureCotonouRepublic of Benin
| | - Julius Lutwama
- Department of Arbovirology, Emerging and Re‐emerging InfectionsUganda Virus Research InstituteEntebbeUganda
| | - Yvonne Marie Linton
- Department of EntomologyNational Museum of Natural HistorySmithsonian InstitutionWashingtonDCUSA
- Walter Reed Biosystematics UnitSmithsonian Institution Museum Support CenterSuitlandMDUSA
- Walter Reed Army Institute of ResearchSilver SpringMDUSA
- Uniformed Services University of Health SciencesBethesdaMDUSA
| | - Catherine Walton
- Faculty of Life SciencesComputational Evolutionary Biology GroupUniversity of ManchesterManchesterUK
| |
Collapse
|
32
|
Ahn J, Conkright B, Boca SM, Madhavan S. POPSTR: Inference of Admixed Population Structure Based on Single-Nucleotide Polymorphisms and Copy Number Variations. J Comput Biol 2018; 25:417-429. [PMID: 29293371 DOI: 10.1089/cmb.2017.0127] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/12/2022] Open
Abstract
Statistical approaches for population structure estimation have been predominantly driven by a particular data type, single-nucleotide polymorphisms (SNPs). However, in the presence of weak identifiability in SNPs, population structure estimation can suffer from undesirable accuracy loss. Copy number variations (CNVs) are genomic structural variants with loci that are commonly shared within a specific population and thus provide valuable information for estimation of the ancestry of sampled populations. We develop a Bayesian joint modeling framework of SNPs and CNVs, called POPSTR, to better understand population structure than approaches that use SNPs solely. To deal with the increased data volume, we use the Metropolis Adjusted Langevin algorithm (MALA) that guides the target distribution in a computationally efficient way. We illustrate applications of our approach using the HapMap 2005 project data. We carry out simulation studies and show that the performance of our approach is comparable or better than that of popular benchmarks, STRUCTURE and ADMIXTURE. We also observe that using only CNVs can be remarkably efficient if SNP data are not available.
Collapse
Affiliation(s)
- Jaeil Ahn
- 1 Department of Biostatistics, Bioinformatics, and Biomathematics, Georgetown University , Washington, District of Columbia
| | - Brian Conkright
- 2 Innovation Center for Biomedical Informatics, Georgetown University , Washington, District of Columbia
| | - Simina M Boca
- 1 Department of Biostatistics, Bioinformatics, and Biomathematics, Georgetown University , Washington, District of Columbia.,2 Innovation Center for Biomedical Informatics, Georgetown University , Washington, District of Columbia.,3 Department of Oncology, Georgetown University , Washington, District of Columbia
| | - Subha Madhavan
- 2 Innovation Center for Biomedical Informatics, Georgetown University , Washington, District of Columbia.,3 Department of Oncology, Georgetown University , Washington, District of Columbia
| |
Collapse
|
33
|
Meyer D, C Aguiar VR, Bitarello BD, C Brandt DY, Nunes K. A genomic perspective on HLA evolution. Immunogenetics 2018; 70:5-27. [PMID: 28687858 PMCID: PMC5748415 DOI: 10.1007/s00251-017-1017-3] [Citation(s) in RCA: 114] [Impact Index Per Article: 19.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2016] [Accepted: 06/16/2017] [Indexed: 12/20/2022]
Abstract
Several decades of research have convincingly shown that classical human leukocyte antigen (HLA) loci bear signatures of natural selection. Despite this conclusion, many questions remain regarding the type of selective regime acting on these loci, the time frame at which selection acts, and the functional connections between genetic variability and natural selection. In this review, we argue that genomic datasets, in particular those generated by next-generation sequencing (NGS) at the population scale, are transforming our understanding of HLA evolution. We show that genomewide data can be used to perform robust and powerful tests for selection, capable of identifying both positive and balancing selection at HLA genes. Importantly, these tests have shown that natural selection can be identified at both recent and ancient timescales. We discuss how findings from genomewide association studies impact the evolutionary study of HLA genes, and how genomic data can be used to survey adaptive change involving interaction at multiple loci. We discuss the methodological developments which are necessary to correctly interpret genomic analyses involving the HLA region. These developments include adapting the NGS analysis framework so as to deal with the highly polymorphic HLA data, as well as developing tools and theory to search for signatures of selection, quantify differentiation, and measure admixture within the HLA region. Finally, we show that high throughput analysis of molecular phenotypes for HLA genes-namely transcription levels-is now a feasible approach and can add another dimension to the study of genetic variation.
Collapse
Affiliation(s)
- Diogo Meyer
- Department of Genetics and Evolutionary Biology, University of São Paulo, 05508-090, São Paulo, SP, Brazil.
| | - Vitor R C Aguiar
- Department of Genetics and Evolutionary Biology, University of São Paulo, 05508-090, São Paulo, SP, Brazil
| | - Bárbara D Bitarello
- Department of Genetics and Evolutionary Biology, University of São Paulo, 05508-090, São Paulo, SP, Brazil
- Department of Evolutionary Genetics, Max Planck Institute for Evolutionary Anthropology, Leipzig, Germany
| | - Débora Y C Brandt
- Department of Genetics and Evolutionary Biology, University of São Paulo, 05508-090, São Paulo, SP, Brazil
- Department of Integrative Biology, University of California, Berkeley, CA, USA
| | - Kelly Nunes
- Department of Genetics and Evolutionary Biology, University of São Paulo, 05508-090, São Paulo, SP, Brazil
| |
Collapse
|
34
|
Nunziata SO, Weisrock DW. Estimation of contemporary effective population size and population declines using RAD sequence data. Heredity (Edinb) 2017; 120:196-207. [PMID: 29269932 DOI: 10.1038/s41437-017-0037-y] [Citation(s) in RCA: 44] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2017] [Revised: 10/12/2017] [Accepted: 11/09/2017] [Indexed: 12/12/2022] Open
Abstract
Large genomic data sets generated with restriction site-associated DNA sequencing (RADseq), in combination with demographic inference methods, are improving our ability to gain insights into the population history of species. We used a simulation approach to examine the potential for RADseq data sets to accurately estimate effective population size (N e) over the course of stable and declining population trends, and we compare the ability of two methods of analysis to accurately distinguish stable from steadily declining populations over a contemporary time scale (20 generations). Using a linkage disequilibrium-based analysis, individual sampling (i.e., n ≥ 30) had the greatest effect on N e estimation and the detection of population size declines, with declines reliably detected across scenarios ~10 generations after they began. Coalescent-based inference required fewer sampled individuals (i.e., n = 15), and instead was most influenced by the size of the SNP data set, with 25,000-50,000 SNPs required for accurate detection of population trends and at least 20 generations after decline began. The number of samples available and targeted number of RADseq loci are important criteria when choosing between these methods. Neither method suffered any apparent bias due to the effects of allele dropout typical of RAD data. With an understanding of the limitations and biases of these approaches, researchers can make more informed decisions when designing their sampling and analyses. Overall, our results reveal that demographic inference using RADseq data can be successfully applied to infer recent population size change and may be an important tool for population monitoring and conservation biology.
Collapse
Affiliation(s)
- Schyler O Nunziata
- Department of Biology, University of Kentucky, Lexington, KY, 40506, USA.
| | - David W Weisrock
- Department of Biology, University of Kentucky, Lexington, KY, 40506, USA
| |
Collapse
|
35
|
Comparison of Single Genome and Allele Frequency Data Reveals Discordant Demographic Histories. G3-GENES GENOMES GENETICS 2017; 7:3605-3620. [PMID: 28893846 PMCID: PMC5677151 DOI: 10.1534/g3.117.300259] [Citation(s) in RCA: 40] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 02/03/2023]
Abstract
Inference of demographic history from genetic data is a primary goal of population genetics of model and nonmodel organisms. Whole genome-based approaches such as the pairwise/multiple sequentially Markovian coalescent methods use genomic data from one to four individuals to infer the demographic history of an entire population, while site frequency spectrum (SFS)-based methods use the distribution of allele frequencies in a sample to reconstruct the same historical events. Although both methods are extensively used in empirical studies and perform well on data simulated under simple models, there have been only limited comparisons of them in more complex and realistic settings. Here we use published demographic models based on data from three human populations (Yoruba, descendants of northwest-Europeans, and Han Chinese) as an empirical test case to study the behavior of both inference procedures. We find that several of the demographic histories inferred by the whole genome-based methods do not predict the genome-wide distribution of heterozygosity, nor do they predict the empirical SFS. However, using simulated data, we also find that the whole genome methods can reconstruct the complex demographic models inferred by SFS-based methods, suggesting that the discordant patterns of genetic variation are not attributable to a lack of statistical power, but may reflect unmodeled complexities in the underlying demography. More generally, our findings indicate that demographic inference from a small number of genomes, routine in genomic studies of nonmodel organisms, should be interpreted cautiously, as these models cannot recapitulate other summaries of the data.
Collapse
|
36
|
Kamdem C, Fouet C, White BJ. Chromosome arm-specific patterns of polymorphism associated with chromosomal inversions in the major African malaria vector, Anopheles funestus. Mol Ecol 2017; 26:5552-5566. [PMID: 28833796 PMCID: PMC5927613 DOI: 10.1111/mec.14335] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/20/2016] [Revised: 08/08/2017] [Accepted: 08/14/2017] [Indexed: 02/02/2023]
Abstract
Chromosomal inversions facilitate local adaptation of beneficial mutations and modulate genetic polymorphism, but the extent of their effects within the genome is still insufficiently understood. The genome of Anopheles funestus, a malaria mosquito endemic to sub-Saharan Africa, contains an impressive number of paracentric polymorphic inversions, which are unevenly distributed among chromosomes and provide an excellent framework for investigating the genomic impacts of chromosomal rearrangements. Here, we present results of a fine-scale analysis of genetic variation within the genome of two weakly differentiated populations of Anopheles funestus inhabiting contrasting moisture conditions in Cameroon. Using population genomic analyses, we found that genetic divergence between the two populations is centred on regions of the genome corresponding to three inversions, which are characterized by high values of FST , absolute sequence divergence and fixed differences. Importantly, in contrast to the 2L chromosome arm, which is collinear, nucleotide diversity is significantly reduced along the entire length of three autosome arms bearing multiple overlapping chromosomal rearrangements. These findings support the idea that interactions between reduced recombination and natural selection within inversions contribute to sculpt nucleotide polymorphism across chromosomes in An. funestus.
Collapse
Affiliation(s)
- Colince Kamdem
- Department of Entomology, University of California, Riverside, CA 92521
| | - Caroline Fouet
- Department of Entomology, University of California, Riverside, CA 92521
| | - Bradley J. White
- Department of Entomology, University of California, Riverside, CA 92521
| |
Collapse
|
37
|
Simons ND, Eick GN, Ruiz-Lopez MJ, Omeja PA, Chapman CA, Goldberg TL, Ting N, Sterner KN. Cis-regulatory evolution in a wild primate: Infection-associated genetic variation drives differential expression of MHC-DQA1 in vitro. Mol Ecol 2017; 26:4523-4535. [PMID: 28665019 PMCID: PMC5570663 DOI: 10.1111/mec.14221] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/27/2016] [Revised: 06/10/2017] [Accepted: 06/13/2017] [Indexed: 12/14/2022]
Abstract
Few studies have combined genetic association analyses with functional characterization of infection-associated SNPs in natural populations of nonhuman primates. Here, we investigate the relationship between host genetic variation, parasitism and natural selection in a population of red colobus (Procolobus rufomitratus tephrosceles) in Kibale National Park, Uganda. We collected parasitological, cellular and genomic data to test the following hypotheses: (i) MHC-DQA1 regulatory genetic variation is associated with control of whipworm (Trichuris) infection in a natural population of red colobus; (ii) infection-associated SNPs are functional in driving differential gene expression in vitro; and (iii) balancing selection has shaped patterns of variation in the MHC-DQA1 promoter. We identified two SNPs in the MHC-DQA1 promoter, both in transcription factor binding sites, and both of which are associated with decreased control of Trichuris infection. We characterized the function of both SNPs by testing differences in gene expression between the two alleles of each SNP in two mammalian cell lines. Alleles of one of the SNPs drove differential gene expression in both cell lines, while the other SNP drove differences in expression in one of the cell lines. Additionally, we found evidence of balancing selection acting on the MHC-DQA1 promoter, including extensive trans-species polymorphisms between red colobus and other primates, and an excess of intermediate-frequency alleles relative to genome-wide, coding and noncoding RADseq data. Our data suggest that balancing selection provides adaptive regulatory flexibility that outweighs the consequences of increased parasite infection intensity in heterozygotes.
Collapse
Affiliation(s)
- Noah D. Simons
- Department of Anthropology, University of Oregon, Eugene, OR 97403, USA
| | - Geeta N. Eick
- Department of Anthropology, University of Oregon, Eugene, OR 97403, USA
| | | | - Patrick A. Omeja
- Makerere University Biological Field Station, P.O Box 967, Fort Portal, Uganda
| | - Colin A. Chapman
- Makerere University Biological Field Station, P.O Box 967, Fort Portal, Uganda
- Department of Anthropology and McGill School of Environment, Montreal, Quebec, H3A 2T7, Canada
| | - Tony L. Goldberg
- Department of Pathobiological Sciences, University of Wisconsin-Madison, Madison, WI 53725, USA
- Global Health Institute, University of Wisconsin-Madison, Madison, WI 53715, USA
| | - Nelson Ting
- Department of Anthropology, University of Oregon, Eugene, OR 97403, USA
- Institute of Ecology and Evolution, University of Oregon, Eugene, OR 97403, USA
| | | |
Collapse
|
38
|
Affiliation(s)
- Sean M. Rovito
- Unidad de Genómica Avanzada (Langebio), Centro de Investigación y de Estudios Avanzados del Instituto Politécnico Nacional, km 9.6 Libramiento Norte Carretera Irapuato-León, Irapuato, Gto CP36821, México
| |
Collapse
|
39
|
Afonso Silva AC, Bragg JG, Potter S, Fernandes C, Coelho MM, Moritz C. Tropical specialist vs. climate generalist: Diversification and demographic history of sister species of
Carlia
skinks from northwestern Australia. Mol Ecol 2017; 26:4045-4058. [DOI: 10.1111/mec.14185] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2016] [Revised: 04/29/2016] [Accepted: 05/02/2017] [Indexed: 01/05/2023]
Affiliation(s)
- Ana C. Afonso Silva
- Research School of Biology and Centre for Biodiversity Analysis Australian National University Acton ACT Australia
- cE3c ‐ Centre for Ecology, Evolution and Environmental Changes Departamento de Biologia Animal Faculdade de Ciências Universidade de Lisboa Lisboa Portugal
| | - Jason G. Bragg
- Research School of Biology and Centre for Biodiversity Analysis Australian National University Acton ACT Australia
- Royal Botanic Garden Sydney NSW Australia
| | - Sally Potter
- Research School of Biology and Centre for Biodiversity Analysis Australian National University Acton ACT Australia
| | - Carlos Fernandes
- cE3c ‐ Centre for Ecology, Evolution and Environmental Changes Departamento de Biologia Animal Faculdade de Ciências Universidade de Lisboa Lisboa Portugal
| | - Maria Manuela Coelho
- cE3c ‐ Centre for Ecology, Evolution and Environmental Changes Departamento de Biologia Animal Faculdade de Ciências Universidade de Lisboa Lisboa Portugal
| | - Craig Moritz
- Research School of Biology and Centre for Biodiversity Analysis Australian National University Acton ACT Australia
| |
Collapse
|
40
|
|
41
|
Genome-culture coevolution promotes rapid divergence of killer whale ecotypes. Nat Commun 2016; 7:11693. [PMID: 27243207 PMCID: PMC4895049 DOI: 10.1038/ncomms11693] [Citation(s) in RCA: 135] [Impact Index Per Article: 16.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2015] [Accepted: 04/18/2016] [Indexed: 12/22/2022] Open
Abstract
Analysing population genomic data from killer whale ecotypes, which we estimate have globally radiated within less than 250,000 years, we show that genetic structuring including the segregation of potentially functional alleles is associated with socially inherited ecological niche. Reconstruction of ancestral demographic history revealed bottlenecks during founder events, likely promoting ecological divergence and genetic drift resulting in a wide range of genome-wide differentiation between pairs of allopatric and sympatric ecotypes. Functional enrichment analyses provided evidence for regional genomic divergence associated with habitat, dietary preferences and post-zygotic reproductive isolation. Our findings are consistent with expansion of small founder groups into novel niches by an initial plastic behavioural response, perpetuated by social learning imposing an altered natural selection regime. The study constitutes an important step towards an understanding of the complex interaction between demographic history, culture, ecological adaptation and evolution at the genomic level. Killer whales have evolved into specialized ecotypes based on hunting strategies and ecological niches. Here, Andrew Foote and colleagues sequenced the whole genome of individual killer whales representing 5 different ecotypes from North Pacific and Antarctic, and show expansion of small founder groups to adapt to specific ecological niches.
Collapse
|
42
|
Boitard S, Rodríguez W, Jay F, Mona S, Austerlitz F. Inferring Population Size History from Large Samples of Genome-Wide Molecular Data - An Approximate Bayesian Computation Approach. PLoS Genet 2016; 12:e1005877. [PMID: 26943927 PMCID: PMC4778914 DOI: 10.1371/journal.pgen.1005877] [Citation(s) in RCA: 99] [Impact Index Per Article: 12.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/31/2015] [Accepted: 01/27/2016] [Indexed: 12/02/2022] Open
Abstract
Inferring the ancestral dynamics of effective population size is a long-standing question in population genetics, which can now be tackled much more accurately thanks to the massive genomic data available in many species. Several promising methods that take advantage of whole-genome sequences have been recently developed in this context. However, they can only be applied to rather small samples, which limits their ability to estimate recent population size history. Besides, they can be very sensitive to sequencing or phasing errors. Here we introduce a new approximate Bayesian computation approach named PopSizeABC that allows estimating the evolution of the effective population size through time, using a large sample of complete genomes. This sample is summarized using the folded allele frequency spectrum and the average zygotic linkage disequilibrium at different bins of physical distance, two classes of statistics that are widely used in population genetics and can be easily computed from unphased and unpolarized SNP data. Our approach provides accurate estimations of past population sizes, from the very first generations before present back to the expected time to the most recent common ancestor of the sample, as shown by simulations under a wide range of demographic scenarios. When applied to samples of 15 or 25 complete genomes in four cattle breeds (Angus, Fleckvieh, Holstein and Jersey), PopSizeABC revealed a series of population declines, related to historical events such as domestication or modern breed creation. We further highlight that our approach is robust to sequencing errors, provided summary statistics are computed from SNPs with common alleles. Molecular data sampled from extant individuals contains considerable information about their demographic history. In particular, one classical question in population genetics is to reconstruct past population size changes from such data. Relating these changes to various climatic, geological or anthropogenic events allows characterizing the main factors driving genetic diversity and can have major outcomes for conservation. Until recently, mostly very simple histories, including one or two population size changes, could be estimated from genetic data. This has changed with the sequencing of entire genomes in many species, and several methods allow now inferring complex histories consisting of several tens of population size changes. However, analyzing entire genomes, while accounting for recombination, remains a statistical and numerical challenge. These methods, therefore, can only be applied to small samples with a few diploid genomes. We overcome this limitation by using an approximate estimation approach, where observed genomes are summarized using a small number of statistics related to allele frequencies and linkage disequilibrium. In contrast to previous approaches, we show that our method allows us to reconstruct also the most recent part (the last 100 generations) of the population size history. As an illustration, we apply it to large samples of whole-genome sequences in four cattle breeds.
Collapse
Affiliation(s)
- Simon Boitard
- Institut de Systématique, Évolution, Biodiversité ISYEB - UMR 7205 - CNRS & MNHN & UPMC & EPHE, Ecole Pratique des Hautes Etudes, Sorbonne Universités, Paris, France
- GABI, INRA, AgroParisTech, Université Paris-Saclay, Jouy-en-Josas, France
- * E-mail:
| | - Willy Rodríguez
- UMR CNRS 5219, Institut de Mathématiques de Toulouse, Université de Toulouse, Toulouse, France
| | - Flora Jay
- UMR 7206 Eco-anthropologie et Ethnobiologie, Muséum National d’Histoire Naturelle, CNRS, Université Paris Diderot, Paris, France
- LRI, Paris-Sud University, CNRS UMR 8623, Orsay, France
| | - Stefano Mona
- Institut de Systématique, Évolution, Biodiversité ISYEB - UMR 7205 - CNRS & MNHN & UPMC & EPHE, Ecole Pratique des Hautes Etudes, Sorbonne Universités, Paris, France
| | - Frédéric Austerlitz
- UMR 7206 Eco-anthropologie et Ethnobiologie, Muséum National d’Histoire Naturelle, CNRS, Université Paris Diderot, Paris, France
| |
Collapse
|
43
|
Browning SR, Browning BL. Accurate Non-parametric Estimation of Recent Effective Population Size from Segments of Identity by Descent. Am J Hum Genet 2015; 97:404-18. [PMID: 26299365 PMCID: PMC4564943 DOI: 10.1016/j.ajhg.2015.07.012] [Citation(s) in RCA: 182] [Impact Index Per Article: 20.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/13/2015] [Accepted: 07/28/2015] [Indexed: 10/23/2022] Open
Abstract
Existing methods for estimating historical effective population size from genetic data have been unable to accurately estimate effective population size during the most recent past. We present a non-parametric method for accurately estimating recent effective population size by using inferred long segments of identity by descent (IBD). We found that inferred segments of IBD contain information about effective population size from around 4 generations to around 50 generations ago for SNP array data and to over 200 generations ago for sequence data. In human populations that we examined, the estimates of effective size were approximately one-third of the census size. We estimate the effective population size of European-ancestry individuals in the UK four generations ago to be eight million and the effective population size of Finland four generations ago to be 0.7 million. Our method is implemented in the open-source IBDNe software package.
Collapse
Affiliation(s)
- Sharon R Browning
- Department of Biostatistics, University of Washington, Seattle, WA 98195, USA.
| | - Brian L Browning
- Division of Medical Genetics, Department of Medicine, University of Washington, Seattle, WA 98195, USA
| |
Collapse
|
44
|
Inferring Bottlenecks from Genome-Wide Samples of Short Sequence Blocks. Genetics 2015; 201:1157-69. [PMID: 26341659 PMCID: PMC4649642 DOI: 10.1534/genetics.115.179861] [Citation(s) in RCA: 34] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/24/2015] [Accepted: 09/01/2015] [Indexed: 01/02/2023] Open
Abstract
The advent of the genomic era has necessitated the development of methods capable of analyzing large volumes of genomic data efficiently. Being able to reliably identify bottlenecks—extreme population size changes of short duration—not only is interesting in the context of speciation and extinction but also matters (as a null model) when inferring selection. Bottlenecks can be detected in polymorphism data via their distorting effect on the shape of the underlying genealogy. Here, we use the generating function of genealogies to derive the probability of mutational configurations in short sequence blocks under a simple bottleneck model. Given a large number of nonrecombining blocks, we can compute maximum-likelihood estimates of the time and strength of the bottleneck. Our method relies on a simple summary of the joint distribution of polymorphic sites. We extend the site frequency spectrum by counting mutations in frequency classes in short sequence blocks. Using linkage information over short distances in this way gives greater power to detect bottlenecks than the site frequency spectrum and potentially opens up a wide range of demographic histories to blockwise inference. Finally, we apply our method to genomic data from a species of pig (Sus cebifrons) endemic to islands in the center and west of the Philippines to estimate whether a bottleneck occurred upon island colonization and compare our scheme to Li and Durbin’s pairwise sequentially Markovian coalescent (PSMC) both for the pig data and using simulations.
Collapse
|
45
|
Stewart JE, Brooks K, Brannen PM, Cline WO, Brewer MT. Elevated Genetic Diversity in the Emerging Blueberry Pathogen Exobasidium maculosum. PLoS One 2015; 10:e0132545. [PMID: 26207812 PMCID: PMC4514876 DOI: 10.1371/journal.pone.0132545] [Citation(s) in RCA: 7] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2015] [Accepted: 06/17/2015] [Indexed: 01/26/2023] Open
Abstract
Emerging diseases caused by fungi are increasing at an alarming rate. Exobasidium leaf and fruit spot of blueberry, caused by the fungus Exobasidium maculosum, is an emerging disease that has rapidly increased in prevalence throughout the southeastern USA, severely reducing fruit quality in some plantings. The objectives of this study were to determine the genetic diversity of E. maculosum in the southeastern USA to elucidate the basis of disease emergence and to investigate if populations of E. maculosum are structured by geography, host species, or tissue type. We sequenced three conserved loci from 82 isolates collected from leaves and fruit of rabbiteye blueberry (Vaccinium virgatum), highbush blueberry (V. corymbosum), and southern highbush blueberry (V. corymbosum hybrids) from commercial fields in Georgia and North Carolina, USA, and 6 isolates from lowbush blueberry (V. angustifolium) from Maine, USA, and Nova Scotia, Canada. Populations of E. maculosum from the southeastern USA and from lowbush blueberry in Maine and Nova Scotia are distinct, but do not represent unique species. No difference in genetic structure was detected between different host tissues or among different host species within the southeastern USA; however, differentiation was detected between populations in Georgia and North Carolina. Overall, E. maculosum showed extreme genetic diversity within the conserved loci with 286 segregating sites among the 1,775 sequenced nucleotides and each isolate representing a unique multilocus haplotype. However, 94% of the nucleotide substitutions were silent, so despite the high number of mutations, selective constraints have limited changes to the amino acid sequences of the housekeeping genes. Overall, these results suggest that the emergence of Exobasidium leaf and fruit spot is not due to a recent introduction or host shift, or the recent evolution of aggressive genotypes of E. maculosum, but more likely as a result of an increasing host population or an environmental change.
Collapse
Affiliation(s)
- Jane E. Stewart
- Department of Plant Pathology, University of Georgia, Athens, Georgia, United States of America
| | - Kyle Brooks
- Department of Plant Pathology, University of Georgia, Athens, Georgia, United States of America
| | - Phillip M. Brannen
- Department of Plant Pathology, University of Georgia, Athens, Georgia, United States of America
| | - William O. Cline
- Department of Plant Pathology, North Carolina State University, Horticultural Crops Research Station, Castle Hayne, North Carolina, United States of America
| | - Marin T. Brewer
- Department of Plant Pathology, University of Georgia, Athens, Georgia, United States of America
| |
Collapse
|
46
|
Morin PA, Parsons KM, Archer FI, Ávila-Arcos MC, Barrett-Lennard LG, Dalla Rosa L, Duchêne S, Durban JW, Ellis GM, Ferguson SH, Ford JK, Ford MJ, Garilao C, Gilbert MTP, Kaschner K, Matkin CO, Petersen SD, Robertson KM, Visser IN, Wade PR, Ho SYW, Foote AD. Geographic and temporal dynamics of a global radiation and diversification in the killer whale. Mol Ecol 2015; 24:3964-79. [PMID: 26087773 DOI: 10.1111/mec.13284] [Citation(s) in RCA: 60] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2014] [Revised: 06/09/2015] [Accepted: 06/17/2015] [Indexed: 02/05/2023]
Abstract
Global climate change during the Late Pleistocene periodically encroached and then released habitat during the glacial cycles, causing range expansions and contractions in some species. These dynamics have played a major role in geographic radiations, diversification and speciation. We investigate these dynamics in the most widely distributed of marine mammals, the killer whale (Orcinus orca), using a global data set of over 450 samples. This marine top predator inhabits coastal and pelagic ecosystems ranging from the ice edge to the tropics, often exhibiting ecological, behavioural and morphological variation suggestive of local adaptation accompanied by reproductive isolation. Results suggest a rapid global radiation occurred over the last 350 000 years. Based on habitat models, we estimated there was only a 15% global contraction of core suitable habitat during the last glacial maximum, and the resources appeared to sustain a constant global effective female population size throughout the Late Pleistocene. Reconstruction of the ancestral phylogeography highlighted the high mobility of this species, identifying 22 strongly supported long-range dispersal events including interoceanic and interhemispheric movement. Despite this propensity for geographic dispersal, the increased sampling of this study uncovered very few potential examples of ancestral dispersal among ecotypes. Concordance of nuclear and mitochondrial data further confirms genetic cohesiveness, with little or no current gene flow among sympatric ecotypes. Taken as a whole, our data suggest that the glacial cycles influenced local populations in different ways, with no clear global pattern, but with secondary contact among lineages following long-range dispersal as a potential mechanism driving ecological diversification.
Collapse
Affiliation(s)
- Phillip A Morin
- Southwest Fisheries Science Center, National Marine Fisheries Service, NOAA, 8901 La Jolla Shores Dr., La Jolla, CA, 92037, USA
| | - Kim M Parsons
- Alaska Fisheries Science Center, National Marine Fisheries Service, NOAA, 7600 Sand Point Way NE, Seattle, WA, 98115, USA
| | - Frederick I Archer
- Southwest Fisheries Science Center, National Marine Fisheries Service, NOAA, 8901 La Jolla Shores Dr., La Jolla, CA, 92037, USA
| | - María C Ávila-Arcos
- Centre for GeoGenetics, Natural History Museum of Denmark, University of Copenhagen, Øster Voldgade 5-7, 1350, Copenhagen, Denmark
| | - Lance G Barrett-Lennard
- Vancouver Aquarium Marine Science Centre, 845 Avison Way, Vancouver, British Columbia, V6G 3E2, Canada
| | - Luciano Dalla Rosa
- Laboratório de Ecologia e Conservação da Megafauna Marinha, Instituto de Oceanografia, Universidade Federal do Rio Grande, Av. Itália km. 8 s/n, Campus Carreiros, Rio Grande, RS, 96201-900, Brazil
| | - Sebastián Duchêne
- School of Biological Sciences, University of Sydney, Sydney, NSW, 2006, Australia
| | - John W Durban
- Southwest Fisheries Science Center, National Marine Fisheries Service, NOAA, 8901 La Jolla Shores Dr., La Jolla, CA, 92037, USA.,Alaska Fisheries Science Center, National Marine Fisheries Service, NOAA, 7600 Sand Point Way NE, Seattle, WA, 98115, USA
| | - Graeme M Ellis
- Fisheries and Oceans Canada, Pacific Biological Station, 3190 Hammond Bay Rd, Nanaimo, British Columbia, Canada
| | - Steven H Ferguson
- Fisheries & Oceans Canada, 501 University Crescent, Winnipeg, Manitoba, R3T 2N6, Canada
| | - John K Ford
- Fisheries and Oceans Canada, Pacific Biological Station, 3190 Hammond Bay Rd, Nanaimo, British Columbia, Canada
| | - Michael J Ford
- Northwest Fisheries Science Center, National Marine Fisheries Service, NOAA 2725 Montlake Blvd E, Seattle, WA, USA
| | - Cristina Garilao
- GEOMAR Helmholtz-Zentrum für Ozeanforschung Kiel Düsternbrooker Weg 2, 24105, Kiel, Germany
| | - M Thomas P Gilbert
- Centre for GeoGenetics, Natural History Museum of Denmark, University of Copenhagen, Øster Voldgade 5-7, 1350, Copenhagen, Denmark.,Trace and Environmental DNA laboratory, Department of Environment and Agriculture, Curtin University, Perth, Western Australia, 6845, Australia
| | - Kristin Kaschner
- Department of Biometry and Environmental System Analysis, Albert-Ludwigs-University of Freiburg, Tennenbacher Strasse 4, 79106, Freiburg, Germany
| | - Craig O Matkin
- North Gulf Oceanic Society, 3430 Main St. Ste. B1, Homer, AK, 99603, USA
| | - Stephen D Petersen
- Assiniboine Park Zoo, 2595 Roblin Blvd, Winnipeg, Manitoba, R3P 2N7, Canada
| | - Kelly M Robertson
- Southwest Fisheries Science Center, National Marine Fisheries Service, NOAA, 8901 La Jolla Shores Dr., La Jolla, CA, 92037, USA
| | - Ingrid N Visser
- Orca Research Trust, P.O. Box 402043, Tutukaka, Northland, 0153, New Zealand
| | - Paul R Wade
- Alaska Fisheries Science Center, National Marine Fisheries Service, NOAA, 7600 Sand Point Way NE, Seattle, WA, 98115, USA
| | - Simon Y W Ho
- School of Biological Sciences, University of Sydney, Sydney, NSW, 2006, Australia
| | - Andrew D Foote
- Centre for GeoGenetics, Natural History Museum of Denmark, University of Copenhagen, Øster Voldgade 5-7, 1350, Copenhagen, Denmark.,Department of Evolutionary Biology, Evolutionary Biology Centre, Uppsala University, Norbyvägen 18D, SE-752 36, Uppsala, Sweden
| |
Collapse
|
47
|
Ferrando A, Manunza A, Jordana J, Capote J, Pons A, Pais J, Delgado T, Atoche P, Cabrera B, Martínez A, Landi V, Delgado JV, Argüello A, Vidal O, Lalueza-Fox C, Ramírez O, Amills M. A mitochondrial analysis reveals distinct founder effect signatures in Canarian and Balearic goats. Anim Genet 2015; 46:452-6. [DOI: 10.1111/age.12302] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 03/25/2015] [Indexed: 11/28/2022]
Affiliation(s)
- A. Ferrando
- Departament de Ciència Animal i dels Aliments; Universitat Autònoma de Barcelona; Bellaterra 08193 Spain
| | - A. Manunza
- Department of Animal Genetics; Center for Research in Agricultural Genomics (CSIC-IRTA-UAB-UB); Campus Universitat Autònoma de Barcelona; Bellaterra 08193 Spain
| | - J. Jordana
- Departament de Ciència Animal i dels Aliments; Universitat Autònoma de Barcelona; Bellaterra 08193 Spain
| | - J. Capote
- Instituto Canario de Investigaciones Agrarias; La Laguna 38108 Tenerife Spain
| | - A. Pons
- Àrea Tècnica Agrària; Servei de Millora Agrària i Pesquera (SEMILLA); Son Ferriol 07198 Spain
| | - J. Pais
- Museo Arqueológico Benahorita; Camino de Las Adelfas; 38760 Los Llanos de Aridane La Palma Spain
| | - T. Delgado
- Museo Canario; Doctor Verneau; 2, 35001 Las Palmas de Gran Canaria Spain
| | - P. Atoche
- Departamento de Ciencias Históricas; Universidad de Las Palmas de Gran Canaria; 35001 Las Palmas de Gran Canaria
| | - B. Cabrera
- Department of Animal Genetics; Center for Research in Agricultural Genomics (CSIC-IRTA-UAB-UB); Campus Universitat Autònoma de Barcelona; Bellaterra 08193 Spain
| | - A. Martínez
- Departamento de Genética; Universidad de Córdoba; Córdoba 14071 Spain
| | - V. Landi
- Departamento de Genética; Universidad de Córdoba; Córdoba 14071 Spain
| | - J. V. Delgado
- Departamento de Genética; Universidad de Córdoba; Córdoba 14071 Spain
| | - A. Argüello
- Department of Animal Science; Universidad de Las Palmas de Gran Canaria; Arucas 35413 Spain
| | - O. Vidal
- Departament de Biologia; Universitat de Girona; Girona 17071 Spain
| | - C. Lalueza-Fox
- Institut de Biologia Evolutiva (CSIC - Universitat Pompeu Fabra); Barcelona 08003 Spain
| | - O. Ramírez
- Institut de Biologia Evolutiva (CSIC - Universitat Pompeu Fabra); Barcelona 08003 Spain
| | - M. Amills
- Departament de Ciència Animal i dels Aliments; Universitat Autònoma de Barcelona; Bellaterra 08193 Spain
- Department of Animal Genetics; Center for Research in Agricultural Genomics (CSIC-IRTA-UAB-UB); Campus Universitat Autònoma de Barcelona; Bellaterra 08193 Spain
| |
Collapse
|
48
|
Hoffmann A, Griffin P, Dillon S, Catullo R, Rane R, Byrne M, Jordan R, Oakeshott J, Weeks A, Joseph L, Lockhart P, Borevitz J, Sgrò C. A framework for incorporating evolutionary genomics into biodiversity conservation and management. ACTA ACUST UNITED AC 2015. [DOI: 10.1186/s40665-014-0009-x] [Citation(s) in RCA: 126] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
|
49
|
Ding X, Li M, Gu H, Peng X, Zhang Z, Wu F. Detecting SNP combinations discriminating human populations from HapMap data. IEEE Trans Nanobioscience 2015; 14:220-8. [PMID: 25594978 DOI: 10.1109/tnb.2015.2391134] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/06/2022]
Abstract
The genomes of different human beings are similar. There are only a relatively small number of genetic differences between people. The genetic differences between people are very worthy of study. Researchers have proposed the fixation index FST measurement to find the single nucleotide polymorphisms (SNPs) which can reflect human population differences. However, most SNPs have interactions and they work together, which leads to the differences among human populations. The number of all possible m-locus combinations chosen from n SNPs grows exponentially. Most methods concern on 2-locus interactions. In this paper, we propose a novel method to find a new coordinate system under which the energy distributions of different populations are quite different. We select out candidate SNPs from n SNPs by using the information of the axes in the coordinate system. The number of candidate SNPs is small, thus SNP-SNP interactions can be searched efficiently. The method can also find interactions of more than two loci. These interactions should be able to reflect the evolution of human populations from another way. The numbers of SNP-SNP interactions are regarded as the differences between pairwise populations and a hierarchical clustering algorithm is used to construct the evolutionary tree. In the experiments, we apply the method to SNP data of four chromosomes separately and the trees constructed on these four chromosomes are highly consistent. Furthermore, the trees are also consistent with previous studies, which indicates that evolutionary information is well mined. The method provides a new insight to analyze the human population differences.
Collapse
|
50
|
Amount of information needed for model choice in Approximate Bayesian Computation. PLoS One 2014; 9:e99581. [PMID: 24959900 PMCID: PMC4069000 DOI: 10.1371/journal.pone.0099581] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2014] [Accepted: 05/16/2014] [Indexed: 11/19/2022] Open
Abstract
Approximate Bayesian Computation (ABC) has become a popular technique in evolutionary genetics for elucidating population structure and history due to its flexibility. The statistical inference framework has benefited from significant progress in recent years. In population genetics, however, its outcome depends heavily on the amount of information in the dataset, whether that be the level of genetic variation or the number of samples and loci. Here we look at the power to reject a simple constant population size coalescent model in favor of a bottleneck model in datasets of varying quality. Not only is this power dependent on the number of samples and loci, but it also depends strongly on the level of nucleotide diversity in the observed dataset. Whilst overall model choice in an ABC setting is fairly powerful and quite conservative with regard to false positives, detecting weaker bottlenecks is problematic in smaller or less genetically diverse datasets and limits the inferences possible in non-model organism where the amount of information regarding the two models is often limited. Our results show it is important to consider these limitations when performing an ABC analysis and that studies should perform simulations based on the size and nature of the dataset in order to fully assess the power of the study.
Collapse
|