1
|
Fraïsse C, Roux C, Gagnaire PA, Romiguier J, Faivre N, Welch JJ, Bierne N. The divergence history of European blue mussel species reconstructed from Approximate Bayesian Computation: the effects of sequencing techniques and sampling strategies. PeerJ 2018; 6:e5198. [PMID: 30083438 PMCID: PMC6071616 DOI: 10.7717/peerj.5198] [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/01/2018] [Accepted: 06/19/2018] [Indexed: 01/25/2023] Open
Abstract
Genome-scale diversity data are increasingly available in a variety of biological systems, and can be used to reconstruct the past evolutionary history of species divergence. However, extracting the full demographic information from these data is not trivial, and requires inferential methods that account for the diversity of coalescent histories throughout the genome. Here, we evaluate the potential and limitations of one such approach. We reexamine a well-known system of mussel sister species, using the joint site frequency spectrum (jSFS) of synonymous mutations computed either from exome capture or RNA-seq, in an Approximate Bayesian Computation (ABC) framework. We first assess the best sampling strategy (number of: individuals, loci, and bins in the jSFS), and show that model selection is robust to variation in the number of individuals and loci. In contrast, different binning choices when summarizing the jSFS, strongly affect the results: including classes of low and high frequency shared polymorphisms can more effectively reveal recent migration events. We then take advantage of the flexibility of ABC to compare more realistic models of speciation, including variation in migration rates through time (i.e., periodic connectivity) and across genes (i.e., genome-wide heterogeneity in migration rates). We show that these models were consistently selected as the most probable, suggesting that mussels have experienced a complex history of gene flow during divergence and that the species boundary is semi-permeable. Our work provides a comprehensive evaluation of ABC demographic inference in mussels based on the coding jSFS, and supplies guidelines for employing different sequencing techniques and sampling strategies. We emphasize, perhaps surprisingly, that inferences are less limited by the volume of data, than by the way in which they are analyzed.
Collapse
Affiliation(s)
- Christelle Fraïsse
- Institut des Sciences de l’Evolution UMR5554, University Montpellier, CNRS, IRD, EPHE, Montpellier, France
- Department of Genetics, University of Cambridge, Cambridge, UK
- Institute of Science and Technology Austria, Klosterneuburg, Austria
| | - Camille Roux
- Université de Lille, Unité Evo-Eco-Paléo (EEP), UMR 8198, Villeneuve d’Ascq, France
| | - Pierre-Alexandre Gagnaire
- Institut des Sciences de l’Evolution UMR5554, University Montpellier, CNRS, IRD, EPHE, Montpellier, France
| | - Jonathan Romiguier
- Institut des Sciences de l’Evolution UMR5554, University Montpellier, CNRS, IRD, EPHE, Montpellier, France
| | - Nicolas Faivre
- Institut des Sciences de l’Evolution UMR5554, University Montpellier, CNRS, IRD, EPHE, Montpellier, France
- Department of Genetics, University of Cambridge, Cambridge, UK
| | - John J. Welch
- Department of Genetics, University of Cambridge, Cambridge, UK
| | - Nicolas Bierne
- Institut des Sciences de l’Evolution UMR5554, University Montpellier, CNRS, IRD, EPHE, Montpellier, France
- Department of Genetics, University of Cambridge, Cambridge, UK
| |
Collapse
|
2
|
Fischer AHL, Mozzherin D, Eren AM, Lans KD, Wilson N, Cosentino C, Smith J. SeaBase: a multispecies transcriptomic resource and platform for gene network inference. Integr Comp Biol 2014; 54:250-63. [PMID: 24907201 DOI: 10.1093/icb/icu065] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Abstract
Marine and aquatic animals are extraordinarily useful as models for identifying mechanisms of development and evolution, regeneration, resistance to cancer, longevity and symbiosis, among many other areas of research. This is due to the great diversity of these organisms and their wide-ranging capabilities. Genomics tools are essential for taking advantage of these "free lessons" of nature. However, genomics and transcriptomics are challenging in emerging model systems. Here, we present SeaBase, a tool for helping to meet these needs. Specifically, SeaBase provides a platform for sharing and searching transcriptome data. More importantly, SeaBase will support a growing number of tools for inferring gene network mechanisms. The first dataset available on SeaBase is a developmental transcriptomic profile of the sea anemone Nematostella vectensis (Anthozoa, Cnidaria). Additional datasets are currently being prepared and we are aiming to expand SeaBase to include user-supplied data for any number of marine and aquatic organisms, thereby supporting many potentially new models for gene network studies. SeaBase can be accessed online at: http://seabase.core.cli.mbl.edu.
Collapse
Affiliation(s)
- Antje H L Fischer
- *Marine Biological Laboratory, Woods Hole, MA 02543, USA; Department of Molecular and Cellular Biology, Harvard University, Cambridge, MA 02138, USA; Systems & Control Engineering, University of Magna Graecia, 88100 Catanzaro, Italy*Marine Biological Laboratory, Woods Hole, MA 02543, USA; Department of Molecular and Cellular Biology, Harvard University, Cambridge, MA 02138, USA; Systems & Control Engineering, University of Magna Graecia, 88100 Catanzaro, Italy
| | - Dmitry Mozzherin
- *Marine Biological Laboratory, Woods Hole, MA 02543, USA; Department of Molecular and Cellular Biology, Harvard University, Cambridge, MA 02138, USA; Systems & Control Engineering, University of Magna Graecia, 88100 Catanzaro, Italy
| | - A Murat Eren
- *Marine Biological Laboratory, Woods Hole, MA 02543, USA; Department of Molecular and Cellular Biology, Harvard University, Cambridge, MA 02138, USA; Systems & Control Engineering, University of Magna Graecia, 88100 Catanzaro, Italy
| | - Kristen D Lans
- *Marine Biological Laboratory, Woods Hole, MA 02543, USA; Department of Molecular and Cellular Biology, Harvard University, Cambridge, MA 02138, USA; Systems & Control Engineering, University of Magna Graecia, 88100 Catanzaro, Italy
| | - Nathan Wilson
- *Marine Biological Laboratory, Woods Hole, MA 02543, USA; Department of Molecular and Cellular Biology, Harvard University, Cambridge, MA 02138, USA; Systems & Control Engineering, University of Magna Graecia, 88100 Catanzaro, Italy
| | - Carlo Cosentino
- *Marine Biological Laboratory, Woods Hole, MA 02543, USA; Department of Molecular and Cellular Biology, Harvard University, Cambridge, MA 02138, USA; Systems & Control Engineering, University of Magna Graecia, 88100 Catanzaro, Italy
| | - Joel Smith
- *Marine Biological Laboratory, Woods Hole, MA 02543, USA; Department of Molecular and Cellular Biology, Harvard University, Cambridge, MA 02138, USA; Systems & Control Engineering, University of Magna Graecia, 88100 Catanzaro, Italy
| |
Collapse
|
3
|
Abstract
Gamete-recognition proteins often evolve rapidly, but it is not known if their divergence occurs within species and corresponds with the evolution of reproductive isolation, or if divergence typically accumulates between already isolated lineages. We examined the evolution of a candidate gamete-recognition protein in several sympatric and allopatric populations of Mytilus blue mussels, species that hybridize in nature. Within a single species, Mytilus galloprovincialis, we found adaptive divergence of Lysin-M7, a sperm acrosomal protein that dissolves the egg vitelline envelope during fertilization. Mytilus galloprovincialis Lysin-M7 alleles group into two distinct clades (termed G and G(D)), and individual alleles in these clades are separated from each other by at least three and up to eleven amino-acid substitutions. Maximum-likelihood estimates of selective pressure (dN/dS =omega) implicate selection in the divergence between M. galloprovincialis Lysin-M7 clades, and within the G(D) clade. Exact tests of population differentiation indicate that the relative frequency of G and G(D) Lysin-M7 alleles differs significantly among M. galloprovincialis populations. Compared with allopatric Mediterranean samples, Lysin-M7 alleles in the G(D) clade are found at elevated frequency in samples from the East Atlantic and California, areas of secondary contact and hybridization between Mytilus species, and Australia, an area of unknown species composition. Adaptive divergence between the alleles most common in allopatry and those found at elevated frequency in samples from sympatry suggests that selection pressures acting in hybridizing populations, likely following Pleistocene secondary contact with M. edulis in the East Atlantic, drove the divergence of Lysin-M7 in M. galloprovincialis.
Collapse
Affiliation(s)
- Stevan A Springer
- Department of Biosciences, Behavioural Ecology Research Group, Simon Fraser University, Burnaby, British Columbia, V5A 1S6, Canada.
| | | |
Collapse
|
5
|
Gilg MR, Hilbish TJ. Patterns of larval dispersal and their effect on the maintenance of a blue mussel hybrid zone in southwestern England. Evolution 2003; 57:1061-77. [PMID: 12836823 DOI: 10.1111/j.0014-3820.2003.tb00316.x] [Citation(s) in RCA: 48] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
The blue mussels Mytilus edulis and M. galloprovincialis hybridize in southwestern England. Within this hybrid zone environmentally based directional selection favors individuals with alleles specific to M. galloprovincialis. What forces are countering this directional selection and allowing for the maintenance of a stable hybrid population are unknown. We used both the genetics of recently settled larvae and a fine-scale model of the physical oceanography of the region to determine the patterns of larval dispersal throughout the hybrid zone and the bordering parental populations. Evidence from both the model and the genetics suggests that the hybrid zone lies between two barriers to dispersal. Start Point separates the M. edulis population from the hybrid zone and allows minimal dispersal from the hybrid zone into the M. edulis population, but none in the other direction. Likewise, the M. galloprovincialis populations along the northern coast of Cornwall regularly receive immigrating larvae from the hybrid zone, but larvae from the M. galloprovincialis population do not enter the hybrid zone. However, larvae settling at hybrid zone sites have high frequencies of alleles specific to M. edulis, suggesting that reproductive barriers, selection in the larval stage, or gene flow from an undetermined source is effectively balancing the directional selection observed in the adults.
Collapse
Affiliation(s)
- Matthew R Gilg
- Department of Biological Sciences, University of South Carolina, Columbia, South Carolina 29208, USA.
| | | |
Collapse
|