1
|
Frei D, De-Kayne R, Selz OM, Seehausen O, Feulner PGD. Genomic variation from an extinct species is retained in the extant radiation following speciation reversal. Nat Ecol Evol 2022; 6:461-468. [PMID: 35210577 DOI: 10.1038/s41559-022-01665-7] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2021] [Accepted: 01/10/2022] [Indexed: 11/09/2022]
Abstract
Ecosystem degradation and biodiversity loss are major global challenges. When reproductive isolation between species is contingent on the interaction of intrinsic lineage traits with features of the environment, environmental change can weaken reproductive isolation and result in extinction through hybridization. By this process called speciation reversal, extinct species can leave traces in genomes of extant species through introgressive hybridization. Using historical and contemporary samples, we sequenced all four species of an Alpine whitefish radiation before and after anthropogenic lake eutrophication and the associated loss of one species through speciation reversal. Despite the extinction of this taxon, substantial fractions of its genome, including regions shaped by positive selection before eutrophication, persist within surviving species as a consequence of introgressive hybridization during eutrophication. Given the prevalence of environmental change, studying speciation reversal and its genomic consequences provides fundamental insights into evolutionary processes and informs biodiversity conservation.
Collapse
Affiliation(s)
- David Frei
- Department of Fish Ecology and Evolution, Center for Ecology, Evolution and Biogeochemistry, Eawag-Swiss Federal Institute of Aquatic Science and Technology, Kastanienbaum, Switzerland.,Division of Aquatic Ecology and Evolution, Institute of Ecology and Evolution, University of Bern, Bern, Switzerland
| | - Rishi De-Kayne
- Department of Fish Ecology and Evolution, Center for Ecology, Evolution and Biogeochemistry, Eawag-Swiss Federal Institute of Aquatic Science and Technology, Kastanienbaum, Switzerland.,Division of Aquatic Ecology and Evolution, Institute of Ecology and Evolution, University of Bern, Bern, Switzerland
| | - Oliver M Selz
- Department of Fish Ecology and Evolution, Center for Ecology, Evolution and Biogeochemistry, Eawag-Swiss Federal Institute of Aquatic Science and Technology, Kastanienbaum, Switzerland.,Division of Aquatic Ecology and Evolution, Institute of Ecology and Evolution, University of Bern, Bern, Switzerland
| | - Ole Seehausen
- Department of Fish Ecology and Evolution, Center for Ecology, Evolution and Biogeochemistry, Eawag-Swiss Federal Institute of Aquatic Science and Technology, Kastanienbaum, Switzerland.,Division of Aquatic Ecology and Evolution, Institute of Ecology and Evolution, University of Bern, Bern, Switzerland
| | - Philine G D Feulner
- Department of Fish Ecology and Evolution, Center for Ecology, Evolution and Biogeochemistry, Eawag-Swiss Federal Institute of Aquatic Science and Technology, Kastanienbaum, Switzerland. .,Division of Aquatic Ecology and Evolution, Institute of Ecology and Evolution, University of Bern, Bern, Switzerland.
| |
Collapse
|
2
|
Otte KA, Schlötterer C. Detecting selected haplotype blocks in evolve and resequence experiments. Mol Ecol Resour 2020; 21:93-109. [PMID: 32810339 PMCID: PMC7754423 DOI: 10.1111/1755-0998.13244] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/24/2020] [Revised: 07/30/2020] [Accepted: 08/04/2020] [Indexed: 12/15/2022]
Abstract
Shifting from the analysis of single nucleotide polymorphisms to the reconstruction of selected haplotypes greatly facilitates the interpretation of evolve and resequence (E&R) experiments. Merging highly correlated hitchhiker SNPs into haplotype blocks reduces thousands of candidates to few selected regions. Current methods of haplotype reconstruction from Pool‐seq data need a variety of data‐specific parameters that are typically defined ad hoc and require haplotype sequences for validation. Here, we introduce haplovalidate, a tool which detects selected haplotypes in Pool‐seq time series data without the need for sequenced haplotypes. Haplovalidate makes data‐driven choices of two key parameters for the clustering procedure, the minimum correlation between SNPs constituting a cluster and the window size. Applying haplovalidate to simulated E&R data reliably detects selected haplotype blocks with low false discovery rates. Importantly, our analyses identified a restriction of the haplotype block‐based approach to describe the genomic architecture of adaptation. We detected a substantial fraction of haplotypes containing multiple selection targets. These blocks were considered as one region of selection and therefore led to underestimation of the number of selection targets. We demonstrate that the separate analysis of earlier time points can significantly increase the separation of selection targets into individual haplotype blocks. We conclude that the analysis of selected haplotype blocks has great potential for the characterization of the adaptive architecture with E&R experiments.
Collapse
Affiliation(s)
- Kathrin A Otte
- Institut für Populationsgenetik, Vetmeduni Vienna, Vienna, Austria
| | | |
Collapse
|
3
|
Parejo M, Wragg D, Henriques D, Charrière JD, Estonba A. Digging into the Genomic Past of Swiss Honey Bees by Whole-Genome Sequencing Museum Specimens. Genome Biol Evol 2020; 12:2535-2551. [PMID: 32877519 PMCID: PMC7720081 DOI: 10.1093/gbe/evaa188] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 08/24/2020] [Indexed: 11/20/2022] Open
Abstract
Historical specimens in museum collections provide opportunities to gain insights into the genomic past. For the Western honey bee, Apis mellifera L., this is particularly important because its populations are currently under threat worldwide and have experienced many changes in management and environment over the last century. Using Swiss Apis mellifera mellifera as a case study, our research provides important insights into the genetic diversity of native honey bees prior to the industrial-scale introductions and trade of non-native stocks during the 20th century—the onset of intensive commercial breeding and the decline of wild honey bees following the arrival of Varroa destructor. We sequenced whole-genomes of 22 honey bees from the Natural History Museum in Bern collected in Switzerland, including the oldest A. mellifera sample ever sequenced. We identify both, a historic and a recent migrant, natural or human-mediated, which corroborates with the population history of honey bees in Switzerland. Contrary to what we expected, we find no evidence for a significant genetic bottleneck in Swiss honey bees, and find that genetic diversity is not only maintained, but even slightly increased, most probably due to modern apicultural practices. Finally, we identify signals of selection between historic and modern honey bee populations associated with genes enriched in functions linked to xenobiotics, suggesting a possible selective pressure from the increasing use and diversity of chemicals used in agriculture and apiculture over the last century.
Collapse
Affiliation(s)
- Melanie Parejo
- Agroscope, Swiss Bee Research Center, Bern, Switzerland.,Lab. Genetics, Department of Genetics, Faculty of Science and Technology, University of the Basque Country (UPV/EHU), Leioa, Spain
| | - David Wragg
- The Roslin Institute, University of Edinburgh, Edinburgh, United Kingdom
| | - Dora Henriques
- Instituto Politécnico de Bragança, Centro de Investigação de Montanha (CIMO), Bragança, Portugal
| | | | - Andone Estonba
- Lab. Genetics, Department of Genetics, Faculty of Science and Technology, University of the Basque Country (UPV/EHU), Leioa, Spain
| |
Collapse
|
4
|
Barghi N, Schlötterer C. Shifting the paradigm in Evolve and Resequence studies: From analysis of single nucleotide polymorphisms to selected haplotype blocks. Mol Ecol 2020; 28:521-524. [PMID: 30793868 PMCID: PMC6850332 DOI: 10.1111/mec.14992] [Citation(s) in RCA: 16] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/19/2018] [Revised: 12/16/2018] [Accepted: 12/18/2018] [Indexed: 12/18/2022]
Abstract
For almost a decade the combination of whole genome sequencing with experimental evolution (Evolve and Resequence, E&R; Turner, Stewart, Fields, Rice, & Tarone, 2011) has been used to study adaptation in outcrossing organisms. However, complications caused by inversions and hitchhiking variants have prevented this powerful approach from living up to its potential. In this issue of Molecular Ecology, Michalak, Kang, Schou, Garner, and Loeschke (2018), provide an important step ahead by using a population of Drosophila melanogaster devoid of segregating inversions to identify the genetic basis of resistance to five environmental stressors. They further address the challenge of hitchhiking variants by reconstructing selected haplotype blocks. While it is apparent that the haplotype block reconstruction needs further refinements, their work underpins the potential of E&R studies in Drosophila to address fundamental questions in evolutionary biology.
Collapse
Affiliation(s)
- Neda Barghi
- Institut für Populationsgenetik, Vetmeduni Vienna, Vienna, Austria
| | | |
Collapse
|
5
|
Bijma P, Wientjes YCJ, Calus MPL. Breeding Top Genotypes and Accelerating Response to Recurrent Selection by Selecting Parents with Greater Gametic Variance. Genetics 2020; 214:91-107. [PMID: 31772074 PMCID: PMC6944402 DOI: 10.1534/genetics.119.302643] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2019] [Accepted: 11/21/2019] [Indexed: 11/18/2022] Open
Abstract
Because of variation in linkage phase and heterozygosity among individuals, some individuals produce genetically more variable gametes than others. With the availability of genomic EBVs (GEBVs) or estimates of SNP-effects together with phased genotypes, differences in gametic variability can be quantified by simulating a set of virtual gametes of each selection candidate. Previous results in dairy cattle show that gametic variance can be large. Here, we show that breeders can increase the probability of breeding a top-ranking genotype and response to recurrent selection by selecting parents that produce more variable gametes, using the index [Formula: see text] where [Formula: see text] is the standardized normal truncation point belonging to selected proportion p, and SDgGEBV is the SD of the GEBV of an individual's gametes. Benefits of the index were considerably larger in an ongoing selection program with equilibrium genetic parameters than in an initially unselected population. Superiority of the index over selection on GEBV increased strongly with the magnitude of the [Formula: see text] indicating that benefits of the index may vary considerably among populations. Compared to selection on ordinary GEBV, the probability of breeding a top-ranking individual can be increased by ∼36%, and response to selection by ∼3.6% when selection is strong (P = 0.001) based on values for the Holstein-Friesian dairy cattle population. Two-stage selection, with a preselection on GEBV and a final selection on the index, considerably reduced computational requirements with little loss of benefits. Response to multiple generations of selection and inheritance of the SDgEBV require further study.
Collapse
Affiliation(s)
- Piter Bijma
- Wageningen University and Research, Animal Breeding and Genomics, 6708PB Wageningen, The Netherlands
| | - Yvonne C J Wientjes
- Wageningen University and Research, Animal Breeding and Genomics, 6708PB Wageningen, The Netherlands
| | - Mario P L Calus
- Wageningen University and Research, Animal Breeding and Genomics, 6708PB Wageningen, The Netherlands
| |
Collapse
|
6
|
A Multireference-Based Whole Genome Assembly for the Obligate Ant-Following Antbird, Rhegmatorhina melanosticta (Thamnophilidae). DIVERSITY-BASEL 2019. [DOI: 10.3390/d11090144] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/16/2022]
Abstract
Current generation high-throughput sequencing technology has facilitated the generation of more genomic-scale data than ever before, thus greatly improving our understanding of avian biology across a range of disciplines. Recent developments in linked-read sequencing (Chromium 10×) and reference-based whole-genome assembly offer an exciting prospect of more accessible chromosome-level genome sequencing in the near future. We sequenced and assembled a genome of the Hairy-crested Antbird (Rhegmatorhina melanosticta), which represents the first publicly available genome for any antbird (Thamnophilidae). Our objectives were to (1) assemble scaffolds to chromosome level based on multiple reference genomes, and report on differences relative to other genomes, (2) assess genome completeness and compare content to other related genomes, and (3) assess the suitability of linked-read sequencing technology for future studies in comparative phylogenomics and population genomics studies. Our R. melanosticta assembly was both highly contiguous (de novo scaffold N50 = 3.3 Mb, reference based N50 = 53.3 Mb) and relatively complete (contained close to 90% of evolutionarily conserved single-copy avian genes and known tetrapod ultraconserved elements). The high contiguity and completeness of this assembly enabled the genome to be successfully mapped to the chromosome level, which uncovered a consistent structural difference between R. melanosticta and other avian genomes. Our results are consistent with the observation that avian genomes are structurally conserved. Additionally, our results demonstrate the utility of linked-read sequencing for non-model genomics. Finally, we demonstrate the value of our R. melanosticta genome for future researchers by mapping reduced representation sequencing data, and by accurately reconstructing the phylogenetic relationships among a sample of thamnophilid species.
Collapse
|
7
|
Medina P, Thornlow B, Nielsen R, Corbett-Detig R. Estimating the Timing of Multiple Admixture Pulses During Local Ancestry Inference. Genetics 2018; 210:1089-1107. [PMID: 30206187 PMCID: PMC6218234 DOI: 10.1534/genetics.118.301411] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2018] [Accepted: 08/31/2018] [Indexed: 11/18/2022] Open
Abstract
Admixture, the mixing of genetically distinct populations, is increasingly recognized as a fundamental biological process. One major goal of admixture analyses is to estimate the timing of admixture events. Whereas most methods today can only detect the most recent admixture event, here, we present coalescent theory and associated software that can be used to estimate the timing of multiple admixture events in an admixed population. We extensively validate this approach and evaluate the conditions under which it can successfully distinguish one- from two-pulse admixture models. We apply our approach to real and simulated data of Drosophila melanogaster We find evidence of a single very recent pulse of cosmopolitan ancestry contributing to African populations, as well as evidence for more ancient admixture among genetically differentiated populations in sub-Saharan Africa. These results suggest our method can quantify complex admixture histories involving genetic material introduced by multiple discrete admixture pulses. The new method facilitates the exploration of admixture and its contribution to adaptation, ecological divergence, and speciation.
Collapse
Affiliation(s)
- Paloma Medina
- Department of Biomolecular Engineering, Genomics Institute, UC Santa Cruz, California 95064
| | - Bryan Thornlow
- Department of Biomolecular Engineering, Genomics Institute, UC Santa Cruz, California 95064
| | - Rasmus Nielsen
- Departments of Integrative Biology and Statistics and Museum of Natural History, University of Copenhagen, 2100 Denmark
| | - Russell Corbett-Detig
- Department of Biomolecular Engineering, Genomics Institute, UC Santa Cruz, California 95064
| |
Collapse
|
8
|
Linkage Disequilibrium Estimation in Low Coverage High-Throughput Sequencing Data. Genetics 2018; 209:389-400. [PMID: 29588288 PMCID: PMC5972415 DOI: 10.1534/genetics.118.300831] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/16/2018] [Accepted: 03/22/2018] [Indexed: 12/31/2022] Open
Abstract
High-throughput sequencing methods provide a cost-effective approach for genotyping and are commonly used in population genetics studies. A drawback of these methods, however, is that sequencing and genotyping errors can arise... High-throughput sequencing methods that multiplex a large number of individuals have provided a cost-effective approach for discovering genome-wide genetic variation in large populations. These sequencing methods are increasingly being utilized in population genetic studies across a diverse range of species. Two side-effects of these methods, however, are (1) sequencing errors and (2) heterozygous genotypes called as homozygous due to only one allele at a particular locus being sequenced, which occurs when the sequencing depth is insufficient. Both of these errors have a profound effect on the estimation of linkage disequilibrium (LD) and, if not taken into account, lead to inaccurate estimates. We developed a new likelihood method, GUS-LD, to estimate pairwise linkage disequilibrium using low coverage sequencing data that accounts for undercalled heterozygous genotypes and sequencing errors. Our findings show that accurate estimates were obtained using GUS-LD, whereas underestimation of LD results if no adjustment is made for the errors.
Collapse
|
9
|
Wit J, Gilleard JS. Resequencing Helminth Genomes for Population and Genetic Studies. Trends Parasitol 2017; 33:388-399. [DOI: 10.1016/j.pt.2017.01.009] [Citation(s) in RCA: 17] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2016] [Revised: 01/11/2017] [Accepted: 01/12/2017] [Indexed: 10/20/2022]
|
10
|
Exploring Evolutionary Relationships Across the Genome Using Topology Weighting. Genetics 2017; 206:429-438. [PMID: 28341652 PMCID: PMC5419486 DOI: 10.1534/genetics.116.194720] [Citation(s) in RCA: 126] [Impact Index Per Article: 18.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2016] [Accepted: 03/14/2017] [Indexed: 12/30/2022] Open
Abstract
We introduce the concept of topology weighting, a method for quantifying relationships between taxa that are not necessarily monophyletic, and visualizing how these relationships change across the genome. A given set of taxa can be related in a limited number of ways, but if each taxon is represented by multiple sequences, the number of possible topologies becomes very large. Topology weighting reduces this complexity by quantifying the contribution of each taxon topology to the full tree. We describe our method for topology weighting by iterative sampling of subtrees (Twisst), and test it on both simulated and real genomic data. Overall, we show that this is an informative and versatile approach, suitable for exploring relationships in almost any genomic dataset. Scripts to implement the method described are available at http://github.com/simonhmartin/twisst.
Collapse
|