1
|
Venu V, Harjunmaa E, Dreau A, Brady S, Absher D, Kingsley DM, Jones FC. Fine-scale contemporary recombination variation and its fitness consequences in adaptively diverging stickleback fish. Nat Ecol Evol 2024; 8:1337-1352. [PMID: 38839849 PMCID: PMC11239493 DOI: 10.1038/s41559-024-02434-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2023] [Accepted: 05/02/2024] [Indexed: 06/07/2024]
Abstract
Despite deep evolutionary conservation, recombination rates vary greatly across the genome and among individuals, sexes and populations. Yet the impact of this variation on adaptively diverging populations is not well understood. Here we characterized fine-scale recombination landscapes in an adaptively divergent pair of marine and freshwater populations of threespine stickleback from River Tyne, Scotland. Through whole-genome sequencing of large nuclear families, we identified the genomic locations of almost 50,000 crossovers and built recombination maps for marine, freshwater and hybrid individuals at a resolution of 3.8 kb. We used these maps to quantify the factors driving variation in recombination rates. We found strong heterochiasmy between sexes but also differences in recombination rates among ecotypes. Hybrids showed evidence of significant recombination suppression in overall map length and in individual loci. Recombination rates were lower not only within individual marine-freshwater-adaptive loci, but also between loci on the same chromosome, suggesting selection on linked gene 'cassettes'. Through temporal sampling along a natural hybrid zone, we found that recombinants showed traits associated with reduced fitness. Our results support predictions that divergence in cis-acting recombination modifiers, whose functions are disrupted in hybrids, may play an important role in maintaining differences among adaptively diverging populations.
Collapse
Affiliation(s)
- Vrinda Venu
- Friedrich Miescher Laboratory of the Max Planck Society, Tübingen, Germany.
- Los Alamos National Laboratory, New Mexico, NM, USA.
| | - Enni Harjunmaa
- Friedrich Miescher Laboratory of the Max Planck Society, Tübingen, Germany
- CeGAT GmbH, Tübingen, Germany
| | - Andreea Dreau
- Friedrich Miescher Laboratory of the Max Planck Society, Tübingen, Germany
- Evotec SE 'Campus Curie', Toulouse, France
| | - Shannon Brady
- Deptartment of Developmental Biology, Stanford University School of Medicine, Stanford, CA, USA
| | - Devin Absher
- HudsonAlpha Institute for Biotechnology, Huntsville, AL, USA
| | - David M Kingsley
- Deptartment of Developmental Biology, Stanford University School of Medicine, Stanford, CA, USA
- Howard Hughes Medical Institute, Chevy Chase, MD, USA
| | - Felicity C Jones
- Friedrich Miescher Laboratory of the Max Planck Society, Tübingen, Germany.
- Groningen Institute for Evolutionary Life Sciences (GELIFES), University of Groningen, Groningen, the Netherlands.
| |
Collapse
|
2
|
Chen Z, Zhou M, Sun Y, Tang X, Zhang Z, Huang L. Exploration of Genome-Wide Recombination Rate Variation Patterns at Different Scales in Pigs. Animals (Basel) 2024; 14:1345. [PMID: 38731349 PMCID: PMC11083071 DOI: 10.3390/ani14091345] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/05/2024] [Revised: 04/27/2024] [Accepted: 04/28/2024] [Indexed: 05/13/2024] Open
Abstract
Meiotic recombination is a prevalent process in eukaryotic sexual reproduction organisms that plays key roles in genetic diversity, breed selection, and species evolution. However, the recombination events differ across breeds and even within breeds. In this study, we initially computed large-scale population recombination rates of both sexes using approximately 52 K SNP genotypes in a total of 3279 pigs from four different Chinese and Western breeds. We then constructed a high-resolution historical recombination map using approximately 16 million SNPs from a sample of unrelated individuals. Comparative analysis of porcine recombination events from different breeds and at different resolutions revealed the following observations: Firstly, the 1Mb-scale pig recombination maps of the same sex are moderately conserved among different breeds, with the similarity of recombination events between Western pigs and Chinese indigenous pigs being lower than within their respective groups. Secondly, we identified 3861 recombination hotspots in the genome and observed medium- to high-level correlation between historical recombination rates (0.542~0.683) and estimates of meiotic recombination rates. Third, we observed that recombination hotspots are significantly far from the transcription start sites of pig genes, and the silico-predicted PRDM9 zinc finger domain DNA recognition motif is significantly enriched in the regions of recombination hotspots compared to recombination coldspots, highlighting the potential role of PRDM9 in regulating recombination hotspots in pigs. Our study analyzed the variation patterns of the pig recombination map at broad and fine scales, providing a valuable reference for genomic selection breeding and laying a crucial foundation for further understanding the molecular mechanisms of pig genome recombination.
Collapse
Affiliation(s)
| | | | | | | | - Zhiyan Zhang
- National Key Laboratory for Swine Genetic Improvement and Germplasm Innovation, Jiangxi Agricultural University, Nanchang 330045, China
| | | |
Collapse
|
3
|
Shaw DE, Naftaly AS, White MA. Positive Selection Drives cis-regulatory Evolution Across the Threespine Stickleback Y Chromosome. Mol Biol Evol 2024; 41:msae020. [PMID: 38306314 PMCID: PMC10899008 DOI: 10.1093/molbev/msae020] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2023] [Revised: 12/31/2023] [Accepted: 01/24/2024] [Indexed: 02/04/2024] Open
Abstract
Allele-specific gene expression evolves rapidly on heteromorphic sex chromosomes. Over time, the accumulation of mutations on the Y chromosome leads to widespread loss of gametolog expression, relative to the X chromosome. It remains unclear if expression evolution on degrading Y chromosomes is primarily driven by mutations that accumulate through processes of selective interference, or if positive selection can also favor the down-regulation of coding regions on the Y chromosome that contain deleterious mutations. Identifying the relative rates of cis-regulatory sequence evolution across Y chromosomes has been challenging due to the limited number of reference assemblies. The threespine stickleback (Gasterosteus aculeatus) Y chromosome is an excellent model to identify how regulatory mutations accumulate on Y chromosomes due to its intermediate state of divergence from the X chromosome. A large number of Y-linked gametologs still exist across 3 differently aged evolutionary strata to test these hypotheses. We found that putative enhancer regions on the Y chromosome exhibited elevated substitution rates and decreased polymorphism when compared to nonfunctional sites, like intergenic regions and synonymous sites. This suggests that many cis-regulatory regions are under positive selection on the Y chromosome. This divergence was correlated with X-biased gametolog expression, indicating the loss of expression from the Y chromosome may be favored by selection. Our findings provide evidence that Y-linked cis-regulatory regions exhibit signs of positive selection quickly after the suppression of recombination and allow comparisons with recent theoretical models that suggest the rapid divergence of regulatory regions may be favored to mask deleterious mutations on the Y chromosome.
Collapse
Affiliation(s)
- Daniel E Shaw
- Department of Genetics, University of Georgia, Athens, GA 30602, USA
| | | | - Michael A White
- Department of Genetics, University of Georgia, Athens, GA 30602, USA
| |
Collapse
|
4
|
Bascón-Cardozo K, Bours A, Manthey G, Durieux G, Dutheil JY, Pruisscher P, Odenthal-Hesse L, Liedvogel M. Fine-Scale Map Reveals Highly Variable Recombination Rates Associated with Genomic Features in the Eurasian Blackcap. Genome Biol Evol 2024; 16:evad233. [PMID: 38198800 PMCID: PMC10781513 DOI: 10.1093/gbe/evad233] [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] [Accepted: 12/12/2023] [Indexed: 01/12/2024] Open
Abstract
Recombination is responsible for breaking up haplotypes, influencing genetic variability, and the efficacy of selection. Bird genomes lack the protein PR domain-containing protein 9, a key determinant of recombination dynamics in most metazoans. Historical recombination maps in birds show an apparent stasis in positioning recombination events. This highly conserved recombination pattern over long timescales may constrain the evolution of recombination in birds. At the same time, extensive variation in recombination rate is observed across the genome and between different species of birds. Here, we characterize the fine-scale historical recombination map of an iconic migratory songbird, the Eurasian blackcap (Sylvia atricapilla), using a linkage disequilibrium-based approach that accounts for population demography. Our results reveal variable recombination rates among and within chromosomes, which associate positively with nucleotide diversity and GC content and negatively with chromosome size. Recombination rates increased significantly at regulatory regions but not necessarily at gene bodies. CpG islands are associated strongly with recombination rates, though their specific position and local DNA methylation patterns likely influence this relationship. The association with retrotransposons varied according to specific family and location. Our results also provide evidence of heterogeneous intrachromosomal conservation of recombination maps between the blackcap and its closest sister taxon, the garden warbler. These findings highlight the considerable variability of recombination rates at different scales and the role of specific genomic features in shaping this variation. This study opens the possibility of further investigating the impact of recombination on specific population-genomic features.
Collapse
Affiliation(s)
- Karen Bascón-Cardozo
- MPRG Behavioural Genomics, Max Planck Institute for Evolutionary Biology, Plön 24306, Germany
| | - Andrea Bours
- MPRG Behavioural Genomics, Max Planck Institute for Evolutionary Biology, Plön 24306, Germany
| | - Georg Manthey
- Institute of Avian Research “Vogelwarte Helgoland”, Wilhelmshaven 26386, Germany
| | - Gillian Durieux
- MPRG Behavioural Genomics, Max Planck Institute for Evolutionary Biology, Plön 24306, Germany
| | - Julien Y Dutheil
- Department for Theoretical Biology, Max Planck Institute for Evolutionary Biology, Plön 24306, Germany
| | - Peter Pruisscher
- MPRG Behavioural Genomics, Max Planck Institute for Evolutionary Biology, Plön 24306, Germany
- Department of Zoology, Stockholm University, Stockholm SE-106 91, Sweden
| | - Linda Odenthal-Hesse
- Department Evolutionary Genetics, Max Planck Institute for Evolutionary Biology, Plön 24306, Germany
| | - Miriam Liedvogel
- MPRG Behavioural Genomics, Max Planck Institute for Evolutionary Biology, Plön 24306, Germany
- Institute of Avian Research “Vogelwarte Helgoland”, Wilhelmshaven 26386, Germany
- Department of Biology and Environmental Sciences, Carl von Ossietzky University of Oldenburg, Oldenburg 26129, Germany
| |
Collapse
|
5
|
Dutreux F, Dutta A, Peltier E, Bibi-Triki S, Friedrich A, Llorente B, Schacherer J. Lessons from the meiotic recombination landscape of the ZMM deficient budding yeast Lachancea waltii. PLoS Genet 2023; 19:e1010592. [PMID: 36608114 PMCID: PMC9851511 DOI: 10.1371/journal.pgen.1010592] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2022] [Revised: 01/19/2023] [Accepted: 12/22/2022] [Indexed: 01/09/2023] Open
Abstract
Meiotic recombination is a driving force for genome evolution, deeply characterized in a few model species, notably in the budding yeast Saccharomyces cerevisiae. Interestingly, Zip2, Zip3, Zip4, Spo16, Msh4, and Msh5, members of the so-called ZMM pathway that implements the interfering meiotic crossover pathway in S. cerevisiae, have been lost in Lachancea yeast species after the divergence of Lachancea kluyveri from the rest of the clade. In this context, after investigating meiosis in L. kluyveri, we determined the meiotic recombination landscape of Lachancea waltii. Attempts to generate diploid strains with fully hybrid genomes invariably resulted in strains with frequent whole-chromosome aneuploidy and multiple extended regions of loss of heterozygosity (LOH), which mechanistic origin is so far unclear. Despite the lack of multiple ZMM pro-crossover factors in L. waltii, numbers of crossovers and noncrossovers per meiosis were higher than in L. kluyveri but lower than in S. cerevisiae, for comparable genome sizes. Similar to L. kluyveri but opposite to S. cerevisiae, L. waltii exhibits an elevated frequency of zero-crossover bivalents. Lengths of gene conversion tracts for both crossovers and non-crossovers in L. waltii were comparable to those observed in S. cerevisiae and shorter than in L. kluyveri despite the lack of Mlh2, a factor limiting conversion tract size in S. cerevisiae. L. waltii recombination hotspots were not shared with either S. cerevisiae or L. kluyveri, showing that meiotic recombination hotspots can evolve at a rather limited evolutionary scale within budding yeasts. Finally, L. waltii crossover interference was reduced relative to S. cerevisiae, with interference being detected only in the 25 kb distance range. Detection of positive inference only at short distance scales in the absence of multiple ZMM factors required for interference-sensitive crossovers in other systems likely reflects interference between early recombination precursors such as DSBs.
Collapse
Affiliation(s)
- Fabien Dutreux
- Université de Strasbourg, CNRS, GMGM UMR 7156, Strasbourg, France
| | - Abhishek Dutta
- Université de Strasbourg, CNRS, GMGM UMR 7156, Strasbourg, France
| | - Emilien Peltier
- Université de Strasbourg, CNRS, GMGM UMR 7156, Strasbourg, France
| | | | - Anne Friedrich
- Université de Strasbourg, CNRS, GMGM UMR 7156, Strasbourg, France
| | - Bertrand Llorente
- CNRS UMR7258, INSERM U1068, Aix Marseille Université UM105, Institut Paoli-Calmettes, CRCM, Marseille, France,* E-mail: (BL); (JS)
| | - Joseph Schacherer
- Université de Strasbourg, CNRS, GMGM UMR 7156, Strasbourg, France,Institut Universitaire de France (IUF), Paris, France,* E-mail: (BL); (JS)
| |
Collapse
|
6
|
Wooldridge LK, Dumont BL. Rapid Evolution of the Fine-scale Recombination Landscape in Wild House Mouse (Mus musculus) Populations. Mol Biol Evol 2022; 40:6889355. [PMID: 36508360 PMCID: PMC9825251 DOI: 10.1093/molbev/msac267] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2021] [Revised: 11/23/2022] [Accepted: 11/30/2022] [Indexed: 12/14/2022] Open
Abstract
Meiotic recombination is an important evolutionary force and an essential meiotic process. In many species, recombination events concentrate into hotspots defined by the site-specific binding of PRMD9. Rapid evolution of Prdm9's zinc finger DNA-binding array leads to remarkably abrupt shifts in the genomic distribution of hotspots between species, but the question of how Prdm9 allelic variation shapes the landscape of recombination between populations remains less well understood. Wild house mice (Mus musculus) harbor exceptional Prdm9 diversity, with >150 alleles identified to date, and pose a particularly powerful system for addressing this open question. We employed a coalescent-based approach to construct broad- and fine-scale sex-averaged recombination maps from contemporary patterns of linkage disequilibrium in nine geographically isolated wild house mouse populations, including multiple populations from each of three subspecies. Comparing maps between wild mouse populations and subspecies reveals several themes. First, we report weak fine- and broad-scale recombination map conservation across subspecies and populations, with genetic divergence offering no clear prediction for recombination map divergence. Second, most hotspots are unique to one population, an outcome consistent with minimal sharing of Prdm9 alleles between surveyed populations. Finally, by contrasting aggregate hotspot activity on the X versus autosomes, we uncover evidence for population-specific differences in the degree and direction of sex dimorphism for recombination. Overall, our findings illuminate the variability of both the broad- and fine-scale recombination landscape in M. musculus and underscore the functional impact of Prdm9 allelic variation in wild mouse populations.
Collapse
|
7
|
Winbush A, Singh ND. Variation in fine-scale recombination rate in temperature-evolved Drosophila melanogaster populations in response to selection. G3 GENES|GENOMES|GENETICS 2022; 12:6663992. [PMID: 35961026 PMCID: PMC9526048 DOI: 10.1093/g3journal/jkac208] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/09/2022] [Accepted: 08/05/2022] [Indexed: 11/16/2022]
Abstract
Meiotic recombination plays a critical evolutionary role in maintaining fitness in response to selective pressures due to changing environments. Variation in recombination rate has been observed amongst and between species and populations and within genomes across numerous taxa. Studies have demonstrated a link between changes in recombination rate and selection, but the extent to which fine-scale recombination rate varies between evolved populations during the evolutionary period in response to selection is under active research. Here, we utilize a set of 3 temperature-evolved Drosophila melanogaster populations that were shown to have diverged in several phenotypes, including recombination rate, based on the temperature regime in which they evolved. Using whole-genome sequencing data from these populations, we generated linkage disequilibrium-based fine-scale recombination maps for each population. With these maps, we compare recombination rates and patterns among the 3 populations and show that they have diverged at fine scales but are conserved at broader scales. We further demonstrate a correlation between recombination rates and genomic variation in the 3 populations. Lastly, we show variation in localized regions of enhanced recombination rates, termed warm spots, between the populations with these warm spots and associated genes overlapping areas previously shown to have diverged in the 3 populations due to selection. These data support the existence of recombination modifiers in these populations which are subject to selection during evolutionary change.
Collapse
Affiliation(s)
- Ari Winbush
- Department of Biology, Institute of Ecology and Evolution, University of Oregon , Eugene, OR 97403, USA
| | - Nadia D Singh
- Department of Biology, Institute of Ecology and Evolution, University of Oregon , Eugene, OR 97403, USA
| |
Collapse
|
8
|
Shi Y, Chen B, Kong S, Zeng Q, Li L, Liu B, Pu F, Xu P. Comparative genomics analysis and genome assembly integration with the recombination landscape contribute to Takifugu bimaculatus assembly refinement. Gene 2022; 849:146910. [PMID: 36167181 DOI: 10.1016/j.gene.2022.146910] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/03/2022] [Revised: 09/13/2022] [Accepted: 09/19/2022] [Indexed: 11/28/2022]
Abstract
Takifugu genus has been brought to the fore in scientific and practical research due to its compact genome, explosive speciation progress and economic value. Here we updated the chromosome-level genome of Takifugu bimaculatus by an ultra-high-density linkage map, a classic and accurate way of chromosome assembly. The map constituted a robust assembly frame, with 92.2% (372.77 Mb) of the draft genome cumulatively placed. With intraspecies and interspecies comparative genomic analysis, we developed a criterion to quantify the differences between assemblies and established a novel way to integrate information from multiple assemblies. The integrated assembly rectified potential mis-assemblies, greatly improving the genome contiguity and correctness. Our results rendered profound information on the genetic recombination of T. bimaculatus and provided new insights into effective genome assembly. The consolidated assembly will be a contributory tool of T. bimaculatus and broadly across the Takifugu by providing a convincing reference for genomic research.
Collapse
Affiliation(s)
- Yue Shi
- State Key Laboratory of Marine Environmental Science, College of Ocean and Earth Sciences, Xiamen University, Xiamen 361102, China; Fujian Key Laboratory of Genetics and Breeding of Marine Organisms, College of Ocean and Earth Sciences, Xiamen University, Xiamen 361102, China
| | - Baohua Chen
- State Key Laboratory of Marine Environmental Science, College of Ocean and Earth Sciences, Xiamen University, Xiamen 361102, China; Fujian Key Laboratory of Genetics and Breeding of Marine Organisms, College of Ocean and Earth Sciences, Xiamen University, Xiamen 361102, China
| | - Shengnan Kong
- State Key Laboratory of Marine Environmental Science, College of Ocean and Earth Sciences, Xiamen University, Xiamen 361102, China; Fujian Key Laboratory of Genetics and Breeding of Marine Organisms, College of Ocean and Earth Sciences, Xiamen University, Xiamen 361102, China
| | - Qingmin Zeng
- Fisheries Research Institute of Fujian, Xiamen 361000, China
| | - Leibin Li
- Fisheries Research Institute of Fujian, Xiamen 361000, China
| | - Bo Liu
- Fisheries Research Institute of Fujian, Xiamen 361000, China
| | - Fei Pu
- State Key Laboratory of Marine Environmental Science, College of Ocean and Earth Sciences, Xiamen University, Xiamen 361102, China; Fujian Key Laboratory of Genetics and Breeding of Marine Organisms, College of Ocean and Earth Sciences, Xiamen University, Xiamen 361102, China
| | - Peng Xu
- State Key Laboratory of Marine Environmental Science, College of Ocean and Earth Sciences, Xiamen University, Xiamen 361102, China; Fujian Key Laboratory of Genetics and Breeding of Marine Organisms, College of Ocean and Earth Sciences, Xiamen University, Xiamen 361102, China.
| |
Collapse
|
9
|
Protacio RU, Davidson MK, Wahls WP. Adaptive Control of the Meiotic Recombination Landscape by DNA Site-dependent Hotspots With Implications for Evolution. Front Genet 2022; 13:947572. [PMID: 35812747 PMCID: PMC9257126 DOI: 10.3389/fgene.2022.947572] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2022] [Accepted: 06/01/2022] [Indexed: 02/04/2023] Open
Abstract
Meiosis is an essential component of the sexual life cycle in eukaryotes. The independent assortment of chromosomes in meiosis increases genetic diversity at the level of whole chromosomes and meiotic recombination increases genetic diversity within chromosomes. The resulting variability fuels evolution. Interestingly, global mapping of recombination in diverse taxa revealed dramatic changes in its frequency distribution between closely related species, subspecies, and even isolated populations of the same species. New insight into mechanisms for these evolutionarily rapid changes has come from analyses of environmentally induced plasticity of recombination in fission yeast. Many different DNA sites, and where identified their binding/activator proteins, control the positioning of recombination at hotspots. Each different class of hotspots functions as an independently controlled rheostat that modulates rates of recombination over a broad dynamic range in response to changing conditions. Together, this independent modulation can rapidly and dramatically alter the global frequency distribution of recombination. This process likely contributes substantially to (i.e., can largely explain) evolutionarily rapid, Prdm9-independent changes in the recombination landscape. Moreover, the precise control mechanisms allow cells to dynamically favor or disfavor newly arising combinations of linked alleles in response to changing extracellular and intracellular conditions, which has striking implications for the impacts of meiotic recombination on evolution.
Collapse
|
10
|
Jenck CS, Lehto WR, Hunnicutt KE, Murphy SM, Quinn TW, Larson EL, Tinghitella RM. Genetic divergence among threespine stickleback that differ in nuptial coloration. J Evol Biol 2022; 35:934-947. [PMID: 35716364 DOI: 10.1111/jeb.14035] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2022] [Revised: 05/11/2022] [Accepted: 05/18/2022] [Indexed: 11/30/2022]
Abstract
Sexual signals are shaped by their intended and unintended receivers as well as the signalling environment. This interplay between sexual and natural selection can lead to divergence in signals in heterogeneous environments. Yet, the extent to which gene flow is restricted when signalling phenotypes vary across environments and over what spatial scales remains an outstanding question. In this study, we quantify gene flow between two colour morphs, red and black, of freshwater threespine stickleback fish (Gasterosteus aculeatus). We capitalize on the very recent divergence of signalling phenotypes in this system to characterize within-species and among-morph genetic variation and to test for levels of gene flow between colour morphs in Oregon and Washington. Despite limited evidence for assortative mating between allopatric red and black populations, we found that black populations are genetically distinct from nearby red populations and that the black morph appears to have evolved independently at least twice in Oregon and Washington. Surprisingly, we uncovered a group of stickleback in one small coastal stream, Connor Creek, which is genetically and morphologically distinct from the red and black colour morphs and from marine stickleback. Historically, both colour morphs have coexisted in this location and sometimes hybridized, raising new questions about the origins and history of these fish, which were first described as anadromous-black hybrids >50 years ago. Understanding how genetic variation is currently partitioned within and among populations and colour morphs in this system should prompt future studies to assess the relative roles of habitat, ecological and pre- and post-reproductive barriers in the genetic divergence and phenotypic patterns we observe in nature.
Collapse
Affiliation(s)
- Clara S Jenck
- Department of Biological Sciences, University of Denver, Denver, Colorado, USA
| | - Whitley R Lehto
- Department of Biological Sciences, University of Denver, Denver, Colorado, USA
| | - Kelsie E Hunnicutt
- Department of Biological Sciences, University of Denver, Denver, Colorado, USA
| | - Shannon M Murphy
- Department of Biological Sciences, University of Denver, Denver, Colorado, USA
| | - Thomas W Quinn
- Department of Biological Sciences, University of Denver, Denver, Colorado, USA
| | - Erica L Larson
- Department of Biological Sciences, University of Denver, Denver, Colorado, USA
| | - Robin M Tinghitella
- Department of Biological Sciences, University of Denver, Denver, Colorado, USA
| |
Collapse
|
11
|
Wall JD, Robinson JA, Cox LA. High-Resolution Estimates of Crossover and Noncrossover Recombination from a Captive Baboon Colony. Genome Biol Evol 2022; 14:evac040. [PMID: 35325119 PMCID: PMC9048888 DOI: 10.1093/gbe/evac040] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 03/02/2022] [Indexed: 11/17/2022] Open
Abstract
Homologous recombination has been extensively studied in humans and a handful of model organisms. Much less is known about recombination in other species, including nonhuman primates. Here, we present a study of crossovers (COs) and noncrossover (NCO) recombination in olive baboons (Papio anubis) from two pedigrees containing a total of 20 paternal and 17 maternal meioses, and compare these results to linkage disequilibrium (LD) based recombination estimates from 36 unrelated olive baboons. We demonstrate how COs, combined with LD-based recombination estimates, can be used to identify genome assembly errors. We also quantify sex-specific differences in recombination rates, including elevated male CO and reduced female CO rates near telomeres. Finally, we add to the increasing body of evidence suggesting that while most NCO recombination tracts in mammals are short (e.g., <500 bp), there is a non-negligible fraction of longer (e.g., >1 kb) NCO tracts. For NCO tracts shorter than 10 kb, we fit a mixture of two (truncated) geometric distributions model to the NCO tract length distribution and estimate that >99% of all NCO tracts are very short (mean 24 bp), but the remaining tracts can be quite long (mean 4.3 kb). A single geometric distribution model for NCO tract lengths is incompatible with the data, suggesting that LD-based methods for estimating NCO recombination rates that make this assumption may need to be modified.
Collapse
Affiliation(s)
- Jeffrey D. Wall
- Institute for Human Genetics, University of California San Francisco, USA
| | | | - Laura A. Cox
- Center for Precision Medicine, Department of Internal Medicine, Wake Forest School of Medicine, Winston-Salem, USA
| |
Collapse
|
12
|
Protacio RU, Mukiza TO, Davidson MK, Wahls WP. Molecular mechanisms for environmentally induced and evolutionarily rapid redistribution (plasticity) of meiotic recombination. Genetics 2022; 220:iyab212. [PMID: 34888655 PMCID: PMC9097252 DOI: 10.1093/genetics/iyab212] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2021] [Accepted: 11/15/2021] [Indexed: 02/04/2023] Open
Abstract
It has long been known (circa 1917) that environmental conditions, as well as speciation, can affect dramatically the frequency distribution of Spo11/Rec12-dependent meiotic recombination. Here, by analyzing DNA sequence-dependent meiotic recombination hotspots in the fission yeast Schizosaccharomyces pombe, we reveal a molecular basis for these phenomena. The impacts of changing environmental conditions (temperature, nutrients, and osmolarity) on local rates of recombination are mediated directly by DNA site-dependent hotspots (M26, CCAAT, and Oligo-C). This control is exerted through environmental condition-responsive signal transduction networks (involving Atf1, Pcr1, Php2, Php3, Php5, and Rst2). Strikingly, individual hotspots modulate rates of recombination over a very broad dynamic range in response to changing conditions. They can range from being quiescent to being highly proficient at promoting activity of the basal recombination machinery (Spo11/Rec12 complex). Moreover, each different class of hotspot functions as an independently controlled rheostat; a condition that increases the activity of one class can decrease the activity of another class. Together, the independent modulation of recombination rates by each different class of DNA site-dependent hotspots (of which there are many) provides a molecular mechanism for highly dynamic, large-scale changes in the global frequency distribution of meiotic recombination. Because hotspot-activating DNA sites discovered in fission yeast are conserved functionally in other species, this process can also explain the previously enigmatic, Prdm9-independent, evolutionarily rapid changes in hotspot usage between closely related species, subspecies, and isolated populations of the same species.
Collapse
Affiliation(s)
- Reine U Protacio
- Department of Biochemistry and Molecular Biology, University of Arkansas for Medical Sciences, Little Rock, AR 72205-7199, USA
| | - Tresor O Mukiza
- Department of Biochemistry and Molecular Biology, University of Arkansas for Medical Sciences, Little Rock, AR 72205-7199, USA
| | - Mari K Davidson
- Department of Biochemistry and Molecular Biology, University of Arkansas for Medical Sciences, Little Rock, AR 72205-7199, USA
| | - Wayne P Wahls
- Department of Biochemistry and Molecular Biology, University of Arkansas for Medical Sciences, Little Rock, AR 72205-7199, USA
| |
Collapse
|
13
|
Liu Z, Roesti M, Marques D, Hiltbrunner M, Saladin V, Peichel CL. Chromosomal fusions facilitate adaptation to divergent environments in threespine stickleback. Mol Biol Evol 2021; 39:6462204. [PMID: 34908155 PMCID: PMC8826639 DOI: 10.1093/molbev/msab358] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
Chromosomal fusions are hypothesized to facilitate adaptation to divergent environments, both by bringing together previously unlinked adaptive alleles and by creating regions of low recombination that facilitate the linkage of adaptive alleles; but, there is little empirical evidence to support this hypothesis. Here, we address this knowledge gap by studying threespine stickleback (Gasterosteus aculeatus), in which ancestral marine fish have repeatedly adapted to freshwater across the northern hemisphere. By comparing the threespine and ninespine stickleback (Pungitius pungitius) genomes to a de novo assembly of the fourspine stickleback (Apeltes quadracus) and an outgroup species, we find two chromosomal fusion events involving the same chromosomes have occurred independently in the threespine and ninespine stickleback lineages. On the fused chromosomes in threespine stickleback, we find an enrichment of quantitative trait loci underlying traits that contribute to marine versus freshwater adaptation. By comparing whole-genome sequences of freshwater and marine threespine stickleback populations, we also find an enrichment of regions under divergent selection on these two fused chromosomes. There is elevated genetic diversity within regions under selection in the freshwater population, consistent with a simulation study showing that gene flow can increase diversity in genomic regions associated with local adaptation and our demographic models showing gene flow between the marine and freshwater populations. Integrating our results with previous studies, we propose that these fusions created regions of low recombination that enabled the formation of adaptative clusters, thereby facilitating freshwater adaptation in the face of recurrent gene flow between marine and freshwater threespine sticklebacks.
Collapse
Affiliation(s)
- Zuyao Liu
- Division of Evolutionary Ecology, Institute of Ecology and Evolution, University of Bern, Bern, Switzerland
| | - Marius Roesti
- Division of Evolutionary Ecology, Institute of Ecology and Evolution, University of Bern, Bern, Switzerland
| | - David Marques
- Division of Aquatic Ecology and Evolution, Institute of Ecology and Evolution, University of Bern, Bern, Switzerland.,Department of Fish Ecology and Evolution, Centre for Ecology, Evolution, and Biogeochemistry, Swiss Federal Institute of Aquatic Science and Technology (EAWAG), Kastanienbaum, Switzerland.,Natural History Museum Basel, Basel, Switzerland
| | - Melanie Hiltbrunner
- Division of Evolutionary Ecology, Institute of Ecology and Evolution, University of Bern, Bern, Switzerland
| | - Verena Saladin
- Division of Evolutionary Ecology, Institute of Ecology and Evolution, University of Bern, Bern, Switzerland
| | - Catherine L Peichel
- Division of Evolutionary Ecology, Institute of Ecology and Evolution, University of Bern, Bern, Switzerland
| |
Collapse
|
14
|
Dahms C, Kemppainen P, Zanella LN, Zanella D, Carosi A, Merilä J, Momigliano P. Cast away in the Adriatic: Low degree of parallel genetic differentiation in three-spined sticklebacks. Mol Ecol 2021; 31:1234-1253. [PMID: 34843145 DOI: 10.1111/mec.16295] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2021] [Revised: 11/16/2021] [Accepted: 11/19/2021] [Indexed: 12/14/2022]
Abstract
The three-spined stickleback (Gasterosteus aculeatus) has repeatedly and independently adapted to freshwater habitats from standing genetic variation (SGV) following colonization from the sea. However, in the Mediterranean Sea G. aculeatus is believed to have gone extinct, and thus the spread of locally adapted alleles between different freshwater populations via the sea since then has been highly unlikely. This is expected to limit parallel evolution, that is the extent to which phylogenetically related alleles can be shared among independently colonized freshwater populations. Using whole genome and 2b-RAD sequencing data, we compared levels of genetic differentiation and genetic parallelism of 15 Adriatic stickleback populations to 19 Pacific, Atlantic and Caspian populations, where gene flow between freshwater populations across extant marine populations is still possible. Our findings support previous studies suggesting that Adriatic populations are highly differentiated (average FST ≈ 0.45), of low genetic diversity and connectivity, and likely to stem from multiple independent colonizations during the Pleistocene. Linkage disequilibrium network analyses in combination with linear mixed models nevertheless revealed several parallel marine-freshwater differentiated genomic regions, although still not to the extent observed elsewhere in the world. We hypothesize that current levels of genetic parallelism in the Adriatic lineages are a relic of freshwater adaptation from SGV prior to the extinction of marine sticklebacks in the Mediterranean that has persisted despite substantial genetic drift experienced by the Adriatic stickleback isolates.
Collapse
Affiliation(s)
- Carolin Dahms
- Ecological Genetics Research Unit, Organismal and Evolutionary Biology Research Programme, Faculty of Biological and Environmental Sciences, University of Helsinki, Helsinki, Finland
| | - Petri Kemppainen
- Ecological Genetics Research Unit, Organismal and Evolutionary Biology Research Programme, Faculty of Biological and Environmental Sciences, University of Helsinki, Helsinki, Finland
| | - Linda N Zanella
- Department of Zoology, Faculty of Science, University of Zagreb, Zagreb, Croatia
| | - Davor Zanella
- Department of Zoology, Faculty of Science, University of Zagreb, Zagreb, Croatia
| | - Antonella Carosi
- Department of Chemistry, Biology and Biotechnologies, University of Perugia, Perugia, Italy
| | - Juha Merilä
- Ecological Genetics Research Unit, Organismal and Evolutionary Biology Research Programme, Faculty of Biological and Environmental Sciences, University of Helsinki, Helsinki, Finland.,Division for Ecology and Biodiversity, School of Biological Sciences, Faculty of Science, The University of Hong Kong, Hong Kong SAR, Hong Kong
| | - Paolo Momigliano
- Ecological Genetics Research Unit, Organismal and Evolutionary Biology Research Programme, Faculty of Biological and Environmental Sciences, University of Helsinki, Helsinki, Finland
| |
Collapse
|
15
|
Turbek SP, Semenov GA, Enbody ED, Campagna L, Taylor SA. Variable Signatures of Selection Despite Conserved Recombination Landscapes Early in Speciation. J Hered 2021; 112:485-496. [PMID: 34499149 DOI: 10.1093/jhered/esab054] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2021] [Accepted: 09/08/2021] [Indexed: 11/14/2022] Open
Abstract
Recently diverged taxa often exhibit heterogeneous landscapes of genomic differentiation, characterized by regions of elevated differentiation on an otherwise homogeneous background. While divergence peaks are generally interpreted as regions responsible for reproductive isolation, they can also arise due to background selection, selective sweeps unrelated to speciation, and variation in recombination and mutation rates. To investigate the association between patterns of recombination and landscapes of genomic differentiation during the early stages of speciation, we generated fine-scale recombination maps for six southern capuchino seedeaters (Sporophila) and two subspecies of White Wagtail (Motacilla alba), two recent avian radiations in which divergent selection on pigmentation genes has likely generated peaks of differentiation. We compared these recombination maps to those of Collared (Ficedula albicollis) and Pied Flycatchers (Ficedula hypoleuca), non-sister taxa characterized by moderate genomic divergence and a heterogenous landscape of genomic differentiation shaped in part by background selection. Although recombination landscapes were conserved within all three systems, we documented a weaker negative correlation between recombination rate and genomic differentiation in the recent radiations. All divergence peaks between capuchinos, wagtails, and flycatchers were located in regions with lower-than-average recombination rates, and most divergence peaks in capuchinos and flycatchers fell in regions of exceptionally reduced recombination. Thus, co-adapted allelic combinations in these regions may have been protected early in divergence, facilitating rapid diversification. Despite largely conserved recombination landscapes, divergence peaks are specific to each focal comparison in capuchinos, suggesting that regions of elevated differentiation have not been generated by variation in recombination rate alone.
Collapse
Affiliation(s)
- Sheela P Turbek
- Department of Ecology and Evolutionary Biology, University of Colorado, Boulder, CO, USA
| | - Georgy A Semenov
- Department of Ecology and Evolutionary Biology, University of Colorado, Boulder, CO, USA
| | - Erik D Enbody
- Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden
| | - Leonardo Campagna
- Fuller Evolutionary Biology Program, Cornell Lab of Ornithology, Ithaca, NY, USA.,Department of Ecology and Evolutionary Biology, Cornell University, Ithaca, NY, USA
| | - Scott A Taylor
- Department of Ecology and Evolutionary Biology, University of Colorado, Boulder, CO, USA
| |
Collapse
|
16
|
Imai Y, Olaya I, Sakai N, Burgess SM. Meiotic Chromosome Dynamics in Zebrafish. Front Cell Dev Biol 2021; 9:757445. [PMID: 34692709 PMCID: PMC8531508 DOI: 10.3389/fcell.2021.757445] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/12/2021] [Accepted: 09/14/2021] [Indexed: 11/13/2022] Open
Abstract
Recent studies in zebrafish have revealed key features of meiotic chromosome dynamics, including clustering of telomeres in the bouquet configuration, biogenesis of chromosome axis structures, and the assembly and disassembly of the synaptonemal complex that aligns homologs end-to-end. The telomere bouquet stage is especially pronounced in zebrafish meiosis and sub-telomeric regions play key roles in mediating pairing and homologous recombination. In this review, we discuss the temporal progression of these events in meiosis prophase I and highlight the roles of proteins associated with meiotic chromosome architecture in homologous recombination. Finally, we discuss the interplay between meiotic mutants and gonadal sex differentiation and future research directions to study meiosis in living cells, including cell culture.
Collapse
Affiliation(s)
- Yukiko Imai
- Department of Gene Function and Phenomics, National Institute of Genetics, Mishima, Japan
| | - Ivan Olaya
- Department of Molecular and Cellular Biology, University of California, Davis, Davis, CA, United States.,Integrative Genetics and Genomics Graduate Group, University of California, Davis, Davis, CA, United States
| | - Noriyoshi Sakai
- Department of Gene Function and Phenomics, National Institute of Genetics, Mishima, Japan.,Department of Genetics, School of Life Sciences, SOKENDAI (The Graduate University for Advanced Studies), Mishima, Japan
| | - Sean M Burgess
- Department of Molecular and Cellular Biology, University of California, Davis, Davis, CA, United States
| |
Collapse
|
17
|
Hu B, Tian Y, Li Q, Liu S. Genomic signatures of artificial selection in the Pacific oyster,
Crassostrea gigas. Evol Appl 2021; 15:618-630. [PMID: 35505882 PMCID: PMC9046764 DOI: 10.1111/eva.13286] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/03/2021] [Revised: 07/06/2021] [Accepted: 07/30/2021] [Indexed: 01/01/2023] Open
Abstract
The Pacific oyster, Crassostrea gigas, is an important aquaculture shellfish around the world with great economic and ecological value. Selective breeding programs have been carried out globally to improve production and performance traits, while genomic signatures of artificial selection remain largely unexplored. In China, we performed selective breeding of C. gigas for over a decade, leading to production of several fast‐growing strains. In the present study, we conducted whole‐genome resequencing of 20 oysters from two fast‐growing strains that have been successively selected for 10 generations, and 20 oysters from the two corresponding wild populations. Sequencing depth of >10× was achieved for each sample, leading to identification of over 12.20 million SNPs. The population structures investigated with three independent methods (principal component analysis, phylogenetic tree, and structure) suggested distinct patterns among selected and wild oyster populations. Assessment of the linkage disequilibrium (LD) decay clearly indicated the changes in genetic diversity during selection. Fixation index (Fst) combined with cross‐population composite likelihood ratio (XP‐CLR) allowed for identification of 768 and 664 selective sweeps (encompassing 1042 and 872 genes) tightly linked to selection in the two fast‐growing strains. KEGG enrichment and functional analyses revealed that 33 genes are important for growth regulation, which act as key components of various signaling pathways with close connection and further take part in regulating the process of cell cycle. This work provides valuable information for the understanding of genomic signatures for long‐term selective breeding and will also be important for growth study and genome‐assisted breeding of the Pacific oyster in the future.
Collapse
Affiliation(s)
- Boyang Hu
- Key Laboratory of Mariculture (Ocean University of China) Ministry of Education, and College of Fisheries Ocean University of China Qingdao China
| | - Yuan Tian
- Key Laboratory of Mariculture (Ocean University of China) Ministry of Education, and College of Fisheries Ocean University of China Qingdao China
| | - Qi Li
- Key Laboratory of Mariculture (Ocean University of China) Ministry of Education, and College of Fisheries Ocean University of China Qingdao China
- Laboratory for Marine Fisheries Science and Food Production Processes Qingdao National Laboratory for Marine Science and Technology Qingdao China
| | - Shikai Liu
- Key Laboratory of Mariculture (Ocean University of China) Ministry of Education, and College of Fisheries Ocean University of China Qingdao China
- Laboratory for Marine Fisheries Science and Food Production Processes Qingdao National Laboratory for Marine Science and Technology Qingdao China
| |
Collapse
|
18
|
Abstract
The repeated adaptation of oceanic threespine sticklebacks to fresh water has made it a premier organism to study parallel evolution. These small fish have multiple distinct ecotypes that display a wide range of diverse phenotypic traits. Ecotypes are easily crossed in the laboratory, and families are large and develop quickly enough for quantitative trait locus analyses, positioning the threespine stickleback as a versatile model organism to address a wide range of biological questions. Extensive genomic resources, including linkage maps, a high-quality reference genome, and developmental genetics tools have led to insights into the genomic basis of adaptation and the identification of genomic changes controlling traits in vertebrates. Recently, threespine sticklebacks have been used as a model system to identify the genomic basis of highly complex traits, such as behavior and host-microbiome and host-parasite interactions. We review the latest findings and new avenues of research that have led the threespine stickleback to be considered a supermodel of evolutionary genomics.
Collapse
Affiliation(s)
- Kerry Reid
- Department of Ecology and Evolution, Stony Brook University, Stony Brook, New York 11794, USA;
| | - Michael A Bell
- University of California Museum of Paleontology, Berkeley, California 94720, USA
| | - Krishna R Veeramah
- Department of Ecology and Evolution, Stony Brook University, Stony Brook, New York 11794, USA;
| |
Collapse
|
19
|
Roberts Kingman GA, Vyas DN, Jones FC, Brady SD, Chen HI, Reid K, Milhaven M, Bertino TS, Aguirre WE, Heins DC, von Hippel FA, Park PJ, Kirch M, Absher DM, Myers RM, Di Palma F, Bell MA, Kingsley DM, Veeramah KR. Predicting future from past: The genomic basis of recurrent and rapid stickleback evolution. SCIENCE ADVANCES 2021; 7:7/25/eabg5285. [PMID: 34144992 PMCID: PMC8213234 DOI: 10.1126/sciadv.abg5285] [Citation(s) in RCA: 42] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/11/2021] [Accepted: 05/05/2021] [Indexed: 05/30/2023]
Abstract
Similar forms often evolve repeatedly in nature, raising long-standing questions about the underlying mechanisms. Here, we use repeated evolution in stickleback to identify a large set of genomic loci that change recurrently during colonization of freshwater habitats by marine fish. The same loci used repeatedly in extant populations also show rapid allele frequency changes when new freshwater populations are experimentally established from marine ancestors. Marked genotypic and phenotypic changes arise within 5 years, facilitated by standing genetic variation and linkage between adaptive regions. Both the speed and location of changes can be predicted using empirical observations of recurrence in natural populations or fundamental genomic features like allelic age, recombination rates, density of divergent loci, and overlap with mapped traits. A composite model trained on these stickleback features can also predict the location of key evolutionary loci in Darwin's finches, suggesting that similar features are important for evolution across diverse taxa.
Collapse
Affiliation(s)
- Garrett A Roberts Kingman
- Department of Developmental Biology, Stanford University School of Medicine, Stanford, CA 94305-5329, USA
| | - Deven N Vyas
- Department of Ecology and Evolution, Stony Brook University, Stony Brook, NY 11794-5245, USA
| | - Felicity C Jones
- Friedrich Miescher Laboratory of the Max Planck Society, Max-Planck-Ring, Tübingen, Germany
| | - Shannon D Brady
- Department of Developmental Biology, Stanford University School of Medicine, Stanford, CA 94305-5329, USA
| | - Heidi I Chen
- Department of Developmental Biology, Stanford University School of Medicine, Stanford, CA 94305-5329, USA
| | - Kerry Reid
- Department of Ecology and Evolution, Stony Brook University, Stony Brook, NY 11794-5245, USA
| | - Mark Milhaven
- Department of Ecology and Evolution, Stony Brook University, Stony Brook, NY 11794-5245, USA
- School of Life Sciences, Arizona State University, Tempe, AZ 85281, USA
| | - Thomas S Bertino
- Department of Ecology and Evolution, Stony Brook University, Stony Brook, NY 11794-5245, USA
| | - Windsor E Aguirre
- Department of Biological Sciences, DePaul University, Chicago, IL 60614-3207, USA
| | - David C Heins
- Department of Ecology and Evolutionary Biology, Tulane University, New Orleans, LA 70118, USA
| | - Frank A von Hippel
- Department of Community, Environment and Policy, Mel & Enid Zuckerman College of Public Health, University of Arizona, Tucson, AZ 85724, USA
| | - Peter J Park
- Department of Biology, Farmingdale State College, Farmingdale, NY 11735-1021, USA
| | - Melanie Kirch
- Friedrich Miescher Laboratory of the Max Planck Society, Max-Planck-Ring, Tübingen, Germany
| | - Devin M Absher
- HudsonAlpha Institute for Biotechnology, 601 Genome Way, Huntsville, AL 35806, USA
| | - Richard M Myers
- HudsonAlpha Institute for Biotechnology, 601 Genome Way, Huntsville, AL 35806, USA
| | - Federica Di Palma
- Broad Institute of MIT and Harvard, 7 Cambridge Center, Cambridge, MA 02142, USA
| | - Michael A Bell
- University of California Museum of Paleontology, University of California, Berkeley, Berkeley, CA 94720, USA.
| | - David M Kingsley
- Department of Developmental Biology, Stanford University School of Medicine, Stanford, CA 94305-5329, USA.
- Howard Hughes Medical Institute, Chevy Chase, MD 20815, USA
| | - Krishna R Veeramah
- Department of Ecology and Evolution, Stony Brook University, Stony Brook, NY 11794-5245, USA.
| |
Collapse
|
20
|
Schield DR, Pasquesi GIM, Perry BW, Adams RH, Nikolakis ZL, Westfall AK, Orton RW, Meik JM, Mackessy SP, Castoe TA. Snake Recombination Landscapes Are Concentrated in Functional Regions despite PRDM9. Mol Biol Evol 2021; 37:1272-1294. [PMID: 31926008 DOI: 10.1093/molbev/msaa003] [Citation(s) in RCA: 31] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
Meiotic recombination in vertebrates is concentrated in hotspots throughout the genome. The location and stability of hotspots have been linked to the presence or absence of PRDM9, leading to two primary models for hotspot evolution derived from mammals and birds. Species with PRDM9-directed recombination have rapid turnover of hotspots concentrated in intergenic regions (i.e., mammals), whereas hotspots in species lacking PRDM9 are concentrated in functional regions and have greater stability over time (i.e., birds). Snakes possess PRDM9, yet virtually nothing is known about snake recombination. Here, we examine the recombination landscape and test hypotheses about the roles of PRDM9 in rattlesnakes. We find substantial variation in recombination rate within and among snake chromosomes, and positive correlations between recombination rate and gene density, GC content, and genetic diversity. Like mammals, snakes appear to have a functional and active PRDM9, but rather than being directed away from genes, snake hotspots are concentrated in promoters and functional regions-a pattern previously associated only with species that lack a functional PRDM9. Snakes therefore provide a unique example of recombination landscapes in which PRDM9 is functional, yet recombination hotspots are associated with functional genic regions-a combination of features that defy existing paradigms for recombination landscapes in vertebrates. Our findings also provide evidence that high recombination rates are a shared feature of vertebrate microchromosomes. Our results challenge previous assumptions about the adaptive role of PRDM9 and highlight the diversity of recombination landscape features among vertebrate lineages.
Collapse
Affiliation(s)
- Drew R Schield
- Department of Biology, University of Texas at Arlington, Arlington, TX
| | | | - Blair W Perry
- Department of Biology, University of Texas at Arlington, Arlington, TX
| | - Richard H Adams
- Department of Biology, University of Texas at Arlington, Arlington, TX.,Department of Computer & Electrical Engineering and Computer Science, Florida Atlantic University, Boca Raton, FL
| | | | | | - Richard W Orton
- Department of Biology, University of Texas at Arlington, Arlington, TX
| | - Jesse M Meik
- Department of Biological Sciences, Tarleton State University, Stephenville, TX
| | - Stephen P Mackessy
- School of Biological Sciences, University of Northern Colorado, Greeley, CO
| | - Todd A Castoe
- Department of Biology, University of Texas at Arlington, Arlington, TX
| |
Collapse
|
21
|
Hu J, Wuitchik SJS, Barry TN, Jamniczky HA, Rogers SM, Barrett RDH. Heritability of DNA methylation in threespine stickleback (Gasterosteus aculeatus). Genetics 2021; 217:1-15. [PMID: 33683369 PMCID: PMC8045681 DOI: 10.1093/genetics/iyab001] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/06/2020] [Accepted: 12/30/2020] [Indexed: 12/13/2022] Open
Abstract
Epigenetic mechanisms underlying phenotypic change are hypothesized to contribute to population persistence and adaptation in the face of environmental change. To date, few studies have explored the heritability of intergenerationally stable methylation levels in natural populations, and little is known about the relative contribution of cis- and trans-regulatory changes to methylation variation. Here, we explore the heritability of DNA methylation, and conduct methylation quantitative trait loci (meQTLs) analysis to investigate the genetic architecture underlying methylation variation between marine and freshwater ecotypes of threespine stickleback (Gasterosteus aculeatus). We quantitatively measured genome-wide DNA methylation in fin tissue using reduced representation bisulfite sequencing of F1 and F2 crosses, and their marine and freshwater source populations. We identified cytosines (CpG sites) that exhibited stable methylation levels across generations. We found that additive genetic variance explained an average of 24-35% of the methylation variance, with a number of CpG sites possibly autonomous from genetic control. We also detected both cis- and trans-meQTLs, with only trans-meQTLs overlapping with previously identified genomic regions of high differentiation between marine and freshwater ecotypes. Finally, we identified the genetic architecture underlying two key CpG sites that were differentially methylated between ecotypes. These findings demonstrate a potential role for DNA methylation in facilitating adaptation to divergent environments and improve our understanding of the heritable basis of population epigenomic variation.
Collapse
Affiliation(s)
- Juntao Hu
- National Observation and Research Station for Yangtze Estuarine Wetland Ecosystems, and Ministry of Education Key Laboratory for Biodiversity Science and Ecological Engineering, School of Life Sciences, Institute of Biodiversity Science, Fudan University, Shanghai 200438, China
- Redpath Museum and Department of Biology, McGill University, Montreal, QC H3A 0C4, Canada
| | - Sara J S Wuitchik
- Informatics Group, Harvard University, Cambridge, MA 02138, USA
- Department of Biology, Boston University, Boston, MA 02215, USA
- Department of Biological Sciences, University of Calgary, Calgary, AB T2N 1N4, Canada
| | - Tegan N Barry
- Department of Biological Sciences, University of Calgary, Calgary, AB T2N 1N4, Canada
| | - Heather A Jamniczky
- Department of Cell Biology and Anatomy, McCaig Institute for Bone and Joint Health, University of Calgary, Calgary, AB T2N 1N4, Canada
| | - Sean M Rogers
- Department of Biological Sciences, University of Calgary, Calgary, AB T2N 1N4, Canada
| | - Rowan D H Barrett
- Redpath Museum and Department of Biology, McGill University, Montreal, QC H3A 0C4, Canada
| |
Collapse
|
22
|
Peichel CL, McCann SR, Ross JA, Naftaly AFS, Urton JR, Cech JN, Grimwood J, Schmutz J, Myers RM, Kingsley DM, White MA. Assembly of the threespine stickleback Y chromosome reveals convergent signatures of sex chromosome evolution. Genome Biol 2020; 21:177. [PMID: 32684159 PMCID: PMC7368989 DOI: 10.1186/s13059-020-02097-x] [Citation(s) in RCA: 63] [Impact Index Per Article: 15.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/20/2020] [Accepted: 07/08/2020] [Indexed: 01/15/2023] Open
Abstract
BACKGROUND Heteromorphic sex chromosomes have evolved repeatedly across diverse species. Suppression of recombination between X and Y chromosomes leads to degeneration of the Y chromosome. The progression of degeneration is not well understood, as complete sequence assemblies of heteromorphic Y chromosomes have only been generated across a handful of taxa with highly degenerate sex chromosomes. Here, we describe the assembly of the threespine stickleback (Gasterosteus aculeatus) Y chromosome, which is less than 26 million years old and at an intermediate stage of degeneration. Our previous work identified that the non-recombining region between the X and the Y spans approximately 17.5 Mb on the X chromosome. RESULTS We combine long-read sequencing with a Hi-C-based proximity guided assembly to generate a 15.87 Mb assembly of the Y chromosome. Our assembly is concordant with cytogenetic maps and Sanger sequences of over 90 Y chromosome BAC clones. We find three evolutionary strata on the Y chromosome, consistent with the three inversions identified by our previous cytogenetic analyses. The threespine stickleback Y shows convergence with more degenerate sex chromosomes in the retention of haploinsufficient genes and the accumulation of genes with testis-biased expression, many of which are recent duplicates. However, we find no evidence for large amplicons identified in other sex chromosome systems. We also report an excellent candidate for the master sex-determination gene: a translocated copy of Amh (Amhy). CONCLUSIONS Together, our work shows that the evolutionary forces shaping sex chromosomes can cause relatively rapid changes in the overall genetic architecture of Y chromosomes.
Collapse
Affiliation(s)
- Catherine L. Peichel
- Divisions of Human Biology and Basic Sciences, Fred Hutchinson Cancer Research Center, Seattle, WA 98109 USA
- Institute of Ecology and Evolution, University of Bern, Baltzerstrasse 6, 3012 Bern, Switzerland
| | - Shaugnessy R. McCann
- Divisions of Human Biology and Basic Sciences, Fred Hutchinson Cancer Research Center, Seattle, WA 98109 USA
| | - Joseph A. Ross
- Divisions of Human Biology and Basic Sciences, Fred Hutchinson Cancer Research Center, Seattle, WA 98109 USA
- Graduate Program in Molecular and Cellular Biology, University of Washington, Seattle, WA 98195 USA
| | | | - James R. Urton
- Divisions of Human Biology and Basic Sciences, Fred Hutchinson Cancer Research Center, Seattle, WA 98109 USA
- Graduate Program in Molecular and Cellular Biology, University of Washington, Seattle, WA 98195 USA
| | - Jennifer N. Cech
- Divisions of Human Biology and Basic Sciences, Fred Hutchinson Cancer Research Center, Seattle, WA 98109 USA
- Graduate Program in Molecular and Cellular Biology, University of Washington, Seattle, WA 98195 USA
| | - Jane Grimwood
- HudsonAlpha Institute for Biotechnology, Huntsville, AL 35806 USA
| | - Jeremy Schmutz
- HudsonAlpha Institute for Biotechnology, Huntsville, AL 35806 USA
| | - Richard M. Myers
- HudsonAlpha Institute for Biotechnology, Huntsville, AL 35806 USA
| | - David M. Kingsley
- Department of Developmental Biology and Howard Hughes Medical Institute, Stanford University School of Medicine, Stanford, CA 94305 USA
| | - Michael A. White
- Divisions of Human Biology and Basic Sciences, Fred Hutchinson Cancer Research Center, Seattle, WA 98109 USA
- Department of Genetics, University of Georgia, Athens, GA 30602 USA
| |
Collapse
|
23
|
Yamasaki YY, Kakioka R, Takahashi H, Toyoda A, Nagano AJ, Machida Y, Møller PR, Kitano J. Genome-wide patterns of divergence and introgression after secondary contact between Pungitius sticklebacks. Philos Trans R Soc Lond B Biol Sci 2020; 375:20190548. [PMID: 32654635 DOI: 10.1098/rstb.2019.0548] [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/11/2022] Open
Abstract
Speciation is a continuous process. Although it is known that differential adaptation can initiate divergence even in the face of gene flow, we know relatively little about the mechanisms driving complete reproductive isolation and the genomic patterns of divergence and introgression at the later stages of speciation. Sticklebacks contain many pairs of sympatric species differing in levels of reproductive isolation and divergence history. Nevertheless, most previous studies have focused on young species pairs. Here, we investigated two sympatric stickleback species, Pungitius pungitius and P. sinensis, whose habitats overlap in eastern Hokkaido; these species show hybrid male sterility, suggesting that they may be at a late stage of speciation. Our demographic analysis using whole-genome sequence data showed that these species split 1.73 Ma and came into secondary contact 37 200 years ago after a period of allopatry. This long period of allopatry might have promoted the evolution of intrinsic incompatibility. Although we detected on-going gene flow and signatures of introgression, overall genomic divergence was high, with considerable heterogeneity across the genome. The heterogeneity was significantly associated with variation in recombination rate. This sympatric pair provides new avenues to investigate the late stages of the stickleback speciation continuum. This article is part of the theme issue 'Towards the completion of speciation: the evolution of reproductive isolation beyond the first barriers'.
Collapse
Affiliation(s)
- Yo Y Yamasaki
- Ecological Genetics Laboratory, National Institute of Genetics, Yata 1111, Mishima, Shizuoka 411-8540, Japan
| | - Ryo Kakioka
- Ecological Genetics Laboratory, National Institute of Genetics, Yata 1111, Mishima, Shizuoka 411-8540, Japan
| | - Hiroshi Takahashi
- National Fisheries University, 2-7-1 Nagata-honmachi, Shimonoseki, Yamaguchi 759-6595, Japan
| | - Atsushi Toyoda
- Comparative Genomics Laboratory, National Institute of Genetics, Yata 1111, Mishima, Shizuoka 411-8540, Japan
| | - Atsushi J Nagano
- Faculty of Agriculture, Ryukoku University, Otsu, Shiga 520-2194, Japan
| | - Yoshiyasu Machida
- Bihoro Museum, Midori 253-4, Bihoro, Abashiri, Hokkaido 092-0002, Japan
| | - Peter R Møller
- Natural History Museum of Denmark, University of Copenhagen, Universitatetsparken 15, Copenhagen 2100, Denmark
| | - Jun Kitano
- Ecological Genetics Laboratory, National Institute of Genetics, Yata 1111, Mishima, Shizuoka 411-8540, Japan
| |
Collapse
|
24
|
Schwarzkopf EJ, Motamayor JC, Cornejo OE. Genetic differentiation and intrinsic genomic features explain variation in recombination hotspots among cocoa tree populations. BMC Genomics 2020; 21:332. [PMID: 32349675 PMCID: PMC7191684 DOI: 10.1186/s12864-020-6746-2] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/10/2020] [Accepted: 04/21/2020] [Indexed: 12/16/2022] Open
Abstract
Background Recombination plays an important evolutionary role by breaking up haplotypes and shuffling genetic variation. This process impacts the ability of selection to eliminate deleterious mutations or increase the frequency of beneficial mutations in a population. To understand the role of recombination generating and maintaining haplotypic variation in a population, we can construct fine-scale recombination maps. Such maps have been used to study a variety of model organisms and proven to be informative of how selection and demographics shape species-wide variation. Here we present a fine-scale recombination map for ten populations of Theobroma cacao – a non-model, long-lived, woody crop. We use this map to elucidate the dynamics of recombination rates in distinct populations of the same species, one of which is domesticated. Results Mean recombination rates in range between 2.5 and 8.6 cM/Mb for most populations of T. cacao with the exception of the domesticated Criollo (525 cM/Mb) and Guianna, a more recently established population (46.5 cM/Mb). We found little overlap in the location of hotspots of recombination across populations. We also found that hotspot regions contained fewer known retroelement sequences than expected and were overrepresented near transcription start and termination sites. We find mutations in FIGL-1, a protein shown to downregulate cross-over frequency in Arabidopsis, statistically associated to higher recombination rates in domesticated Criollo. Conclusions We generated fine-scale recombination maps for ten populations of Theobroma cacao and used them to understand what processes are associated with population-level variation in this species. Our results provide support to the hypothesis of increased recombination rates in domesticated plants (Criollo population). We propose a testable mechanistic hypothesis for the change in recombination rate in domesticated populations in the form of mutations to a previously identified recombination-suppressing protein. Finally, we establish a number of possible correlates of recombination hotspots that help explain general patterns of recombination in this species.
Collapse
Affiliation(s)
| | | | - Omar E Cornejo
- School of Biological Sciences, Washington State University, Pullman, WA, USA.
| |
Collapse
|
25
|
Yoshida K, Ravinet M, Makino T, Toyoda A, Kokita T, Mori S, Kitano J. Accumulation of Deleterious Mutations in Landlocked Threespine Stickleback Populations. Genome Biol Evol 2020; 12:479-492. [PMID: 32232440 PMCID: PMC7197494 DOI: 10.1093/gbe/evaa065] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 03/24/2020] [Indexed: 02/06/2023] Open
Abstract
Colonization of new habitats often reduces population sizes and may result in the accumulation of deleterious mutations by genetic drift. Compared with the genomic basis for adaptation to new environments, genome-wide analysis of deleterious mutations in isolated populations remains limited. In the present study, we investigated the accumulation of deleterious mutations in five endangered freshwater populations of threespine stickleback (Gasterosteus aculeatus) in the central part of the mainland of Japan. Using whole-genome resequencing data, we first conducted phylogenomic analysis and confirmed at least two independent freshwater colonization events in the central mainland from ancestral marine ecotypes. Next, analyses of single nucleotide polymorphisms showed a substantial reduction of heterozygosity in freshwater populations compared with marine populations. Reduction in heterozygosity was more apparent at the center of each chromosome than the peripheries and on X chromosomes compared with autosomes. Third, bioinformatic analysis of deleterious mutations showed increased accumulation of putatively deleterious mutations in the landlocked freshwater populations compared with marine populations. For the majority of populations examined, the frequencies of putatively deleterious mutations were higher on X chromosomes than on autosomes. The interpopulation comparison indicated that the majority of putatively deleterious mutations may have accumulated independently. Thus, whole-genome resequencing of endangered populations can help to estimate the accumulation of deleterious mutations and inform us of which populations are the most severely endangered. Furthermore, analysis of variation among chromosomes can give insights into whether any particular chromosomes are likely to accumulate deleterious mutations.
Collapse
Affiliation(s)
- Kohta Yoshida
- Ecological Genetics Laboratory, National Institute of Genetics, Mishima, Shizuoka, Japan.,Max Planck Institute for Developmental Biology, Tübingen, Germany
| | - Mark Ravinet
- Ecological Genetics Laboratory, National Institute of Genetics, Mishima, Shizuoka, Japan.,Department of Biosciences, Centre for Ecological and Evolutionary Synthesis, University of Oslo, Norway.,School of Life Sciences, University of Nottingham, University Park, Nottingham, United Kingdom
| | - Takashi Makino
- Graduate School of Life Sciences, Tohoku University, Sendai, Miyagi, Japan
| | - Atsushi Toyoda
- Comparative Genomics Laboratory, National Institute of Genetics, Mishima, Shizuoka, Japan
| | - Tomoyuki Kokita
- Department of Marine Bioscience, Fukui Prefectural University, Obama, Fukui, Japan
| | - Seiichi Mori
- Biological Laboratories, Gifu-kyoritsu University, Ogaki, Gifu, Japan
| | - Jun Kitano
- Ecological Genetics Laboratory, National Institute of Genetics, Mishima, Shizuoka, Japan
| |
Collapse
|
26
|
Inferring the Genomic Landscape of Recombination Rate Variation in European Aspen ( Populus tremula). G3-GENES GENOMES GENETICS 2020; 10:299-309. [PMID: 31744900 PMCID: PMC6945010 DOI: 10.1534/g3.119.400504] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
The rate of meiotic recombination is one of the central factors determining genome-wide levels of linkage disequilibrium which has important consequences for the efficiency of natural selection and for the dissection of quantitative traits. Here we present a new, high-resolution linkage map for Populus tremula that we use to anchor approximately two thirds of the P. tremula draft genome assembly on to the expected 19 chromosomes, providing us with the first chromosome-scale assembly for P. tremula (Table 2). We then use this resource to estimate variation in recombination rates across the P. tremula genome and compare these results to recombination rates based on linkage disequilibrium in a large number of unrelated individuals. We also assess how variation in recombination rates is associated with a number of genomic features, such as gene density, repeat density and methylation levels. We find that recombination rates obtained from the two methods largely agree, although the LD-based method identifies a number of genomic regions with very high recombination rates that the map-based method fails to detect. Linkage map and LD-based estimates of recombination rates are positively correlated and show similar correlations with other genomic features, showing that both methods can accurately infer recombination rate variation across the genome. Recombination rates are positively correlated with gene density and negatively correlated with repeat density and methylation levels, suggesting that recombination is largely directed toward gene regions in P. tremula.
Collapse
|
27
|
Spence JP, Song YS. Inference and analysis of population-specific fine-scale recombination maps across 26 diverse human populations. SCIENCE ADVANCES 2019; 5:eaaw9206. [PMID: 31681842 PMCID: PMC6810367 DOI: 10.1126/sciadv.aaw9206] [Citation(s) in RCA: 72] [Impact Index Per Article: 14.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/05/2019] [Accepted: 09/13/2019] [Indexed: 05/28/2023]
Abstract
Fine-scale rates of meiotic recombination vary by orders of magnitude across the genome and differ between species and even populations. Studying cross-population differences has been stymied by the confounding effects of demographic history. To address this problem, we developed a demography-aware method to infer fine-scale recombination rates and applied it to 26 diverse human populations, inferring population-specific recombination maps. These maps recapitulate many aspects of the history of these populations including signatures of the trans-Atlantic slave trade and the Iberian colonization of the Americas. We also investigated modulators of the local recombination rate, finding further evidence that Polycomb group proteins and the trimethylation of H3K27 elevate recombination rates. Further differences in the recombination landscape across the genome and between populations are driven by variation in the gene that encodes the DNA binding protein PRDM9, and we quantify the weak effect of meiotic drive acting to remove its binding sites.
Collapse
Affiliation(s)
- Jeffrey P. Spence
- Graduate Group in Computational Biology, University of California, Berkeley, Berkeley, CA, USA
| | - Yun S. Song
- Computer Science Division and Department of Statistics, University of California, Berkeley, Berkeley, CA, USA
- Chan Zuckerberg Biohub, San Francisco, CA, USA
| |
Collapse
|