1
|
Nikolic N, Devloo-Delva F, Bailleul D, Noskova E, Rougeux C, Delord C, Borsa P, Liautard-Haag C, Hassan M, Marie AD, Feutry P, Grewe P, Davies C, Farley J, Fernando D, Biton-Porsmoguer S, Poisson F, Parker D, Leone A, Aulich J, Lansdell M, Marsac F, Arnaud-Haond S. Stepping up to genome scan allows stock differentiation in the worldwide distributed blue shark Prionace glauca. Mol Ecol 2023; 32:1000-1019. [PMID: 36511846 DOI: 10.1111/mec.16822] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2022] [Revised: 11/29/2022] [Accepted: 12/01/2022] [Indexed: 12/15/2022]
Abstract
The blue shark Prionace glauca is a top predator with one of the widest geographical distributions of any shark species. It is classified as Critically Endangered in the Mediterranean Sea, and Near Threatened globally. Previous genetic studies did not reject the null hypothesis of a single global population. The blue shark was proposed as a possible archetype of the "grey zone of population differentiation," coined to designate cases where population structure may be too recent or too faint to be detected using a limited set of markers. Here, blue shark samples collected throughout its global range were sequenced using a specific RAD method (DArTseq), which recovered 37,655 genome-wide single nucleotide polymorphisms (SNPs). Two main groups emerged, with Mediterranean Sea and northern Atlantic samples (Northern population) differentiated significantly from the Indo-west Pacific samples (Southern population). Significant pairwise FST values indicated further genetic differentiation within the Atlantic Ocean, and between the Atlantic Ocean and the Mediterranean Sea. Reconstruction of recent demographic history suggested divergence between Northern and Southern populations occurred about 500 generations ago and revealed a drastic reduction in effective population size from a large ancestral population. Our results illustrate the power of genome scans to detect population structure and reconstruct demographic history in highly migratory marine species. Given that the management plans of the blue shark (targeted or bycatch) fisheries currently assume panmictic regional stocks, we strongly recommend that the results presented here be considered in future stock assessments and conservation strategies.
Collapse
Affiliation(s)
- Natacha Nikolic
- UMR MARBEC, University of Montpellier, IRD, Ifremer, CNRS, Sète, France.,INRAE, Ecobiop, AQUA, Saint-Pée-sur-Nivelle, France.,ARBRE, Agence de Recherche pour la Biodiversité à la Réunion, Saint-Gilles, France
| | - Floriaan Devloo-Delva
- CSIRO Environment, Hobart, Tasmania, Australia.,School of Natural Sciences-Quantitative Marine Science, University of Tasmania, Hobart, Tasmania, Australia
| | - Diane Bailleul
- UMR MARBEC, University of Montpellier, IRD, Ifremer, CNRS, Sète, France
| | - Ekaterina Noskova
- Computer Technologies Laboratory, ITMO University, St Petersburg, Russia
| | | | - Chrystelle Delord
- UMR MARBEC, University of Montpellier, IRD, Ifremer, CNRS, Sète, France
| | - Philippe Borsa
- Institut de recherche pour le développement, UMR ENTROPIE, Montpellier, France
| | | | - Mohamad Hassan
- UMR MARBEC, University of Montpellier, IRD, Ifremer, CNRS, Sète, France.,Animal Production Department, Tishreen University, Latakia, Syria
| | - Amandine D Marie
- ARBRE, Agence de Recherche pour la Biodiversité à la Réunion, Saint-Gilles, France
| | | | - Peter Grewe
- CSIRO Environment, Hobart, Tasmania, Australia
| | | | | | | | | | - François Poisson
- UMR MARBEC, University of Montpellier, IRD, Ifremer, CNRS, Sète, France
| | - Denham Parker
- Department of Forestry, Fisheries and the Environment, (DFFE), Cape Town, South Africa.,Department of Biological Sciences, University of Cape Town, Rondebosch, South Africa
| | - Agostino Leone
- UMR MARBEC, University of Montpellier, IRD, Ifremer, CNRS, Sète, France
| | | | | | - Francis Marsac
- UMR MARBEC, University of Montpellier, IRD, Ifremer, CNRS, Sète, France
| | | |
Collapse
|
2
|
Thia JA, Korhonen PK, Young ND, Gasser RB, Umina PA, Yang Q, Edwards O, Walsh T, Hoffmann AA. The redlegged earth mite draft genome provides new insights into pesticide resistance evolution and demography in its invasive Australian range. J Evol Biol 2023; 36:381-398. [PMID: 36573922 PMCID: PMC10107102 DOI: 10.1111/jeb.14144] [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: 08/04/2022] [Revised: 10/13/2022] [Accepted: 11/03/2022] [Indexed: 12/28/2022]
Abstract
Genomic data provide valuable insights into pest management issues such as resistance evolution, historical patterns of pest invasions and ongoing population dynamics. We assembled the first reference genome for the redlegged earth mite, Halotydeus destructor (Tucker, 1925), to investigate adaptation to pesticide pressures and demography in its invasive Australian range using whole-genome pool-seq data from regionally distributed populations. Our reference genome comprises 132 autosomal contigs, with a total length of 48.90 Mb. We observed a large complex of ace genes, which has presumably evolved from a long history of organophosphate selection in H. destructor and may contribute towards organophosphate resistance through copy number variation, target-site mutations and structural variants. In the putative ancestral H. destructor ace gene, we identified three target-site mutations (G119S, A201S and F331Y) segregating in organophosphate-resistant populations. Additionally, we identified two new para sodium channel gene mutations (L925I and F1020Y) that may contribute to pyrethroid resistance. Regional structuring observed in population genomic analyses indicates that gene flow in H. destructor does not homogenize populations across large geographic distances. However, our demographic analyses were equivocal on the magnitude of gene flow; the short invasion history of H. destructor makes it difficult to distinguish scenarios of complete isolation vs. ongoing migration. Nonetheless, we identified clear signatures of reduced genetic diversity and smaller inferred effective population sizes in eastern vs. western populations, which is consistent with the stepping-stone invasion pathway of this pest in Australia. These new insights will inform development of diagnostic genetic markers of resistance, further investigation into the multifaceted organophosphate resistance mechanism and predictive modelling of resistance evolution and spread.
Collapse
Affiliation(s)
- Joshua A Thia
- Bio21 Institute, School of BioSciences, The University of Melbourne, Melbourne, Victoria, Australia
| | - Pasi K Korhonen
- Department of Veterinary Biosciences, Melbourne Veterinary School, The University of Melbourne, Melbourne, Victoria, Australia
| | - Neil D Young
- Department of Veterinary Biosciences, Melbourne Veterinary School, The University of Melbourne, Melbourne, Victoria, Australia
| | - Robin B Gasser
- Department of Veterinary Biosciences, Melbourne Veterinary School, The University of Melbourne, Melbourne, Victoria, Australia
| | | | - Qiong Yang
- Bio21 Institute, School of BioSciences, The University of Melbourne, Melbourne, Victoria, Australia
| | - Owain Edwards
- Land and Water, CSIRO, Floreat, Western Australia, Australia
| | - Tom Walsh
- CSIRO, Black Mountain Laboratories, Canberra, Australian Capital Territory, Australia.,Applied BioSciences, Macquarie University, Sydney, New South Wales, Australia
| | - Ary A Hoffmann
- Bio21 Institute, School of BioSciences, The University of Melbourne, Melbourne, Victoria, Australia
| |
Collapse
|
3
|
Martin Cerezo ML, Raval R, de Haro Reyes B, Kucka M, Chan FY, Bryk J. Identification and quantification of chimeric sequencing reads in a highly multiplexed RAD-seq protocol. Mol Ecol Resour 2022; 22:2860-2870. [PMID: 35668693 PMCID: PMC9796921 DOI: 10.1111/1755-0998.13661] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/15/2021] [Revised: 04/29/2022] [Accepted: 05/23/2022] [Indexed: 01/07/2023]
Abstract
Highly multiplexed approaches have become common in genomic studies. They have improved the cost-effectiveness of genotyping hundreds of individuals using combinatorially barcoded adapters. These strategies, however, can potentially misassigned reads to incorrect samples. Here, we used a modified quaddRAD protocol to analyse the occurrence of index hopping and PCR chimeras in a series of experiments with up to 100 multiplexed samples per sequencing lane (639 samples in total). We created two types of sequencing libraries: four libraries of type A, where PCRs were run on individual samples before multiplexing, and three libraries of type B, where PCRs were run on pooled samples. We used fixed pairs of inner barcodes to identify chimeric reads. Type B libraries show a higher percentage of misassigned reads (1.15%) than type A libraries (0.65%). We also quantify the commonly undetectable chimeric sequences that occur whenever multiplexed groups of samples with different outer barcodes are sequenced together on a single flow cell. Our results suggest that these types of chimeric sequences represent up to 1.56% and 1.29% of reads in type A and B libraries, respectively. We also show that increasing the number of mismatches allowed for barcode rescue to above 2 dramatically increases the number of recovered chimeric reads. We provide recommendations for developing highly multiplexed RAD-seq protocols and analysing the resulting data to minimize the generation of chimeric sequences, allowing their quantification and a finer control on the number of PCR cycles necessary to generate enough input DNA for library preparation.
Collapse
Affiliation(s)
- Maria Luisa Martin Cerezo
- Department of Biological and Geographical Sciences, School of Applied SciencesUniversity of HuddersfieldHuddersfieldUK,IFM BiologyLinköping UniversityLinköpingSweden
| | - Rohan Raval
- Department of Biological and Geographical Sciences, School of Applied SciencesUniversity of HuddersfieldHuddersfieldUK
| | - Bernardo de Haro Reyes
- Department of Biological and Geographical Sciences, School of Applied SciencesUniversity of HuddersfieldHuddersfieldUK,IFM BiologyLinköping UniversityLinköpingSweden
| | - Marek Kucka
- Friedrich Miescher Laboratory of the Max Planck SocietyTübingenGermany
| | | | - Jarosław Bryk
- Department of Biological and Geographical Sciences, School of Applied SciencesUniversity of HuddersfieldHuddersfieldUK
| |
Collapse
|
4
|
Jafari O, Zeinalabedini M, Robledo D, Fernandes JMO, Hedayati AA, Arefnezhad B. Genotyping-by-Sequencing Reveals the Impact of Restocking on Wild Common Carp Populations of the Southern Caspian Basin. Front Ecol Evol 2022. [DOI: 10.3389/fevo.2022.872176] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023] Open
Abstract
Understanding the population structure and level of genetic diversity of wild populations is fundamental for appropriate stock management and species conservation. The common carp (Cyprinus carpio) is one of the most important bony fish throughout the Southern coastline of the Caspian Sea, but captures of this species have seen a dramatic reduction during the last decade. As a consequence, a restocking program has been put in place to maintain C. carpio populations, but its impact is not clear. In the present study, the population structure and genetic diversity of C. carpio in the Southern Caspian basin was determined using 17,828 single-nucleotide polymorphism (SNP) markers. A total of 117 individuals collected from four different locations in the southern Caspian basin and a farm were genotyped by genotyping-by-sequencing. The overall Fst obtained was 0.04, indicating a low level of differentiation between populations, and most genetic diversity was attributed to within population variation (97%). The low Fst values suggest that frequent migration events between different locations occur, and three migration events were inferred in the present study. However, each population still showed a distinct genetic profile, which allowed distinguishing the origin of the fish. This indicates that the ongoing restocking program is maintaining the differences between populations to some extent. Nonetheless, high inbreeding and low heterozygosity were detected in all populations, suggesting that additional conservation efforts are required to protect C. carpio populations in the Southern coast of the Caspian Sea.
Collapse
|
5
|
Aguirre‐Sarabia I, Díaz‐Arce N, Pereda‐Agirre I, Mendibil I, Urtizberea A, Gerritsen HD, Burns F, Holmes I, Landa J, Coscia I, Quincoces I, Santurtún M, Zanzi A, Martinsohn JT, Rodríguez‐Ezpeleta N. Evidence of stock connectivity, hybridization, and misidentification in white anglerfish supports the need of a genetics-informed fisheries management framework. Evol Appl 2021; 14:2221-2230. [PMID: 34603494 PMCID: PMC8477593 DOI: 10.1111/eva.13278] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2021] [Revised: 06/29/2021] [Accepted: 07/02/2021] [Indexed: 11/27/2022] Open
Abstract
Understanding population connectivity within a species as well as potential interactions with its close relatives is crucial to define management units and to derive efficient management actions. However, although genetics can reveal mismatches between biological and management units and other relevant but hidden information such as species misidentification or hybridization, the uptake of genetic methods by the fisheries management process is far from having been consolidated. Here, we have assessed the power of genetics to better understand the population connectivity of white (Lophius piscatorius) and its interaction with its sister species, the black anglerfish (Lophius budegassa). Our analyses, based on thousands of genome-wide single nucleotide polymorphisms, show three findings that are crucial for white anglerfish management. We found (i) that white anglerfish is likely composed of a single panmictic population throughout the Northeast Atlantic, challenging the three-stock based management, (ii) that a fraction of specimens classified as white anglerfish using morphological characteristics are genetically identified as black anglerfish (L. budegassa), and iii) that the two Lophius species naturally hybridize leading to a population of hybrids of up to 20% in certain areas. Our results set the basics for a genetics-informed white anglerfish assessment framework that accounts for stock connectivity, revises and establishes new diagnostic characters for Lophius species identification, and evaluates the effect of hybrids in the current and future assessments of the white anglerfish. Furthermore, our study contributes to provide additional evidence of the potentially negative consequences of ignoring genetic data for assessing fisheries resources.
Collapse
Affiliation(s)
- Imanol Aguirre‐Sarabia
- Marine ResearchAZTI Basque Research and Technology Alliance (BRTA)SukarrietaBizkaiaSpain
| | - Natalia Díaz‐Arce
- Marine ResearchAZTI Basque Research and Technology Alliance (BRTA)SukarrietaBizkaiaSpain
| | - Iker Pereda‐Agirre
- Marine ResearchAZTI Basque Research and Technology Alliance (BRTA)SukarrietaBizkaiaSpain
| | - Iñaki Mendibil
- Marine ResearchAZTI Basque Research and Technology Alliance (BRTA)SukarrietaBizkaiaSpain
| | - Agurtzane Urtizberea
- Marine ResearchAZTI Basque Research and Technology Alliance (BRTA)SukarrietaBizkaiaSpain
| | | | - Finlay Burns
- Marine LaboratoryMarine Scotland ScienceAberdeenshireUK
| | - Ian Holmes
- Lowestoft LaboratoryCentre for Environment, Fisheries and Aquaculture ScienceLowestoftSuffolkUK
| | - Jorge Landa
- Centro Oceanográfico de SantanderInstituto Español de Oceanografía (IEO)SantanderSpain
| | - Ilaria Coscia
- School of Science, Engineering and EnvironmentUniversity of SalfordSalfordUK
| | - Iñaki Quincoces
- Marine ResearchAZTI Basque Research and Technology Alliance (BRTA)SukarrietaBizkaiaSpain
| | - Marina Santurtún
- Marine ResearchAZTI Basque Research and Technology Alliance (BRTA)SukarrietaBizkaiaSpain
| | | | | | | |
Collapse
|
6
|
Christiansen H, Heindler FM, Hellemans B, Jossart Q, Pasotti F, Robert H, Verheye M, Danis B, Kochzius M, Leliaert F, Moreau C, Patel T, Van de Putte AP, Vanreusel A, Volckaert FAM, Schön I. Facilitating population genomics of non-model organisms through optimized experimental design for reduced representation sequencing. BMC Genomics 2021; 22:625. [PMID: 34418978 PMCID: PMC8380342 DOI: 10.1186/s12864-021-07917-3] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/02/2021] [Accepted: 07/26/2021] [Indexed: 12/26/2022] Open
Abstract
BACKGROUND Genome-wide data are invaluable to characterize differentiation and adaptation of natural populations. Reduced representation sequencing (RRS) subsamples a genome repeatedly across many individuals. However, RRS requires careful optimization and fine-tuning to deliver high marker density while being cost-efficient. The number of genomic fragments created through restriction enzyme digestion and the sequencing library setup must match to achieve sufficient sequencing coverage per locus. Here, we present a workflow based on published information and computational and experimental procedures to investigate and streamline the applicability of RRS. RESULTS In an iterative process genome size estimates, restriction enzymes and size selection windows were tested and scaled in six classes of Antarctic animals (Ostracoda, Malacostraca, Bivalvia, Asteroidea, Actinopterygii, Aves). Achieving high marker density would be expensive in amphipods, the malacostracan target taxon, due to the large genome size. We propose alternative approaches such as mitogenome or target capture sequencing for this group. Pilot libraries were sequenced for all other target taxa. Ostracods, bivalves, sea stars, and fish showed overall good coverage and marker numbers for downstream population genomic analyses. In contrast, the bird test library produced low coverage and few polymorphic loci, likely due to degraded DNA. CONCLUSIONS Prior testing and optimization are important to identify which groups are amenable for RRS and where alternative methods may currently offer better cost-benefit ratios. The steps outlined here are easy to follow for other non-model taxa with little genomic resources, thus stimulating efficient resource use for the many pressing research questions in molecular ecology.
Collapse
Affiliation(s)
- Henrik Christiansen
- Laboratory of Biodiversity and Evolutionary Genomics, KU Leuven, Leuven, Belgium.
| | - Franz M Heindler
- Laboratory of Biodiversity and Evolutionary Genomics, KU Leuven, Leuven, Belgium
| | - Bart Hellemans
- Laboratory of Biodiversity and Evolutionary Genomics, KU Leuven, Leuven, Belgium
| | - Quentin Jossart
- Marine Biology Group, Vrije Universiteit Brussel (VUB), Brussels, Belgium
| | | | - Henri Robert
- OD Nature, Royal Belgian Institute of Natural Sciences, Brussels, Belgium
| | - Marie Verheye
- OD Nature, Royal Belgian Institute of Natural Sciences, Brussels, Belgium
| | - Bruno Danis
- Marine Biology Laboratory, Université Libre de Bruxelles (ULB), Brussels, Belgium
| | - Marc Kochzius
- Marine Biology Group, Vrije Universiteit Brussel (VUB), Brussels, Belgium
| | - Frederik Leliaert
- Marine Biology Research Group, Ghent University, Ghent, Belgium.,Meise Botanic Garden, Meise, Belgium
| | - Camille Moreau
- Marine Biology Laboratory, Université Libre de Bruxelles (ULB), Brussels, Belgium.,Université de Bourgogne Franche-Comté (UBFC) UMR CNRS 6282 Biogéosciences, Dijon, France
| | - Tasnim Patel
- OD Nature, Royal Belgian Institute of Natural Sciences, Brussels, Belgium
| | - Anton P Van de Putte
- Laboratory of Biodiversity and Evolutionary Genomics, KU Leuven, Leuven, Belgium.,OD Nature, Royal Belgian Institute of Natural Sciences, Brussels, Belgium.,Marine Biology Laboratory, Université Libre de Bruxelles (ULB), Brussels, Belgium
| | - Ann Vanreusel
- Marine Biology Research Group, Ghent University, Ghent, Belgium
| | - Filip A M Volckaert
- Laboratory of Biodiversity and Evolutionary Genomics, KU Leuven, Leuven, Belgium
| | - Isa Schön
- OD Nature, Royal Belgian Institute of Natural Sciences, Brussels, Belgium
| |
Collapse
|
7
|
Feron R, Pan Q, Wen M, Imarazene B, Jouanno E, Anderson J, Herpin A, Journot L, Parrinello H, Klopp C, Kottler VA, Roco AS, Du K, Kneitz S, Adolfi M, Wilson CA, McCluskey B, Amores A, Desvignes T, Goetz FW, Takanashi A, Kawaguchi M, Detrich HW, Oliveira MA, Nóbrega RH, Sakamoto T, Nakamoto M, Wargelius A, Karlsen Ø, Wang Z, Stöck M, Waterhouse RM, Braasch I, Postlethwait JH, Schartl M, Guiguen Y. RADSex: A computational workflow to study sex determination using restriction site-associated DNA sequencing data. Mol Ecol Resour 2021; 21:1715-1731. [PMID: 33590960 DOI: 10.1111/1755-0998.13360] [Citation(s) in RCA: 20] [Impact Index Per Article: 6.7] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2020] [Revised: 02/09/2021] [Accepted: 02/11/2021] [Indexed: 12/13/2022]
Abstract
The study of sex determination and sex chromosome organization in nonmodel species has long been technically challenging, but new sequencing methodologies now enable precise and high-throughput identification of sex-specific genomic sequences. In particular, restriction site-associated DNA sequencing (RAD-Seq) is being extensively applied to explore sex determination systems in many plant and animal species. However, software specifically designed to search for and visualize sex-biased markers using RAD-Seq data is lacking. Here, we present RADSex, a computational analysis workflow designed to study the genetic basis of sex determination using RAD-Seq data. RADSex is simple to use, requires few computational resources, makes no prior assumptions about the type of sex-determination system or structure of the sex locus, and offers convenient visualization through a dedicated R package. To demonstrate the functionality of RADSex, we re-analysed a published data set of Japanese medaka, Oryzias latipes, where we uncovered a previously unknown Y chromosome polymorphism. We then used RADSex to analyse new RAD-Seq data sets from 15 fish species spanning multiple taxonomic orders. We identified the sex determination system and sex-specific markers in six of these species, five of which had no known sex-markers prior to this study. We show that RADSex greatly facilitates the study of sex determination systems in nonmodel species thanks to its speed of analyses, low resource usage, ease of application and visualization options. Furthermore, our analysis of new data sets from 15 species provides new insights on sex determination in fish.
Collapse
Affiliation(s)
- Romain Feron
- INRAE, LPGP, Rennes, France.,Department of Ecology and Evolution, University of Lausanne, Lausanne, Switzerland.,Swiss Institute of Bioinformatics, Lausanne, Switzerland
| | - Qiaowei Pan
- INRAE, LPGP, Rennes, France.,Department of Ecology and Evolution, University of Lausanne, Lausanne, Switzerland
| | - Ming Wen
- INRAE, LPGP, Rennes, France.,State Key Laboratory of Developmental Biology of Freshwater Fish, College of Life Science, Hunan Normal University, Changsha, China
| | | | | | - Jennifer Anderson
- INRAE, LPGP, Rennes, France.,Department of Organismal Biology, Systematic Biology, Uppsala University, Uppsala, Sweden
| | | | - Laurent Journot
- Institut de Génomique Fonctionnelle, IGF, CNRS, INSERM, Univ. Montpellier, Montpellier, France
| | - Hugues Parrinello
- Institut de Génomique Fonctionnelle, IGF, CNRS, INSERM, Univ. Montpellier, Montpellier, France
| | - Christophe Klopp
- SIGENAE, Mathématiques et Informatique Appliquées de Toulouse, INRAE, Castanet Tolosan, France
| | - Verena A Kottler
- Physiological Chemistry, Biocenter, University of Wuerzburg, Wuerzburg, Germany
| | - Alvaro S Roco
- Physiological Chemistry, Biocenter, University of Wuerzburg, Wuerzburg, Germany
| | - Kang Du
- Department of Chemistry and Biochemistry, The Xiphophorus Genetic Stock Center, Texas State University, San Marcos, TX, USA.,Developmental Biochemistry, Biocenter, University of Wuerzburg, Wuerzburg, Germany
| | - Susanne Kneitz
- Physiological Chemistry, Biocenter, University of Wuerzburg, Wuerzburg, Germany
| | - Mateus Adolfi
- Developmental Biochemistry, Biocenter, University of Wuerzburg, Wuerzburg, Germany
| | | | | | - Angel Amores
- Institute of Neuroscience, University of Oregon, Eugene, OR, USA
| | - Thomas Desvignes
- Institute of Neuroscience, University of Oregon, Eugene, OR, USA
| | - Frederick W Goetz
- Environmental and Fisheries Sciences Division, Northwest Fisheries Science Center, National Marine Fisheries Service, National Oceanic and Atmospheric Administration, Seattle, WA, USA
| | - Ato Takanashi
- Department of Materials and Life Sciences, Faculty of Science and Technology, Sophia University, Tokyo, Japan
| | - Mari Kawaguchi
- Department of Materials and Life Sciences, Faculty of Science and Technology, Sophia University, Tokyo, Japan
| | - Harry William Detrich
- Department of Marine and Environmental Sciences, Marine Science Center, Northeastern University, Nahant, MA, USA
| | - Marcos A Oliveira
- Reproductive and Molecular Biology Group, Department of Structural and Functional Biology, Institute of Biosciences, São Paulo State University, Botucatu, Brazil
| | - Rafael H Nóbrega
- Reproductive and Molecular Biology Group, Department of Structural and Functional Biology, Institute of Biosciences, São Paulo State University, Botucatu, Brazil
| | - Takashi Sakamoto
- Department of Aquatic Marine Biosciences, Tokyo University of Marine Science and Technology, Tokyo, Japan
| | - Masatoshi Nakamoto
- Department of Aquatic Marine Biosciences, Tokyo University of Marine Science and Technology, Tokyo, Japan
| | | | | | - Zhongwei Wang
- Physiological Chemistry, Biocenter, University of Wuerzburg, Wuerzburg, Germany.,Institute of Hydrobiology, Chinese Academy of Sciences, Beijing, China
| | - Matthias Stöck
- Leibniz-Institute of Freshwater Ecology and Inland Fisheries, IGB, Berlin, Germany
| | - Robert M Waterhouse
- Department of Ecology and Evolution, University of Lausanne, Lausanne, Switzerland.,Swiss Institute of Bioinformatics, Lausanne, Switzerland
| | - Ingo Braasch
- Department of Integrative Biology, Ecology, Evolution, and Behavior Program, Michigan State University, East Lansing, MI, USA
| | | | - Manfred Schartl
- Department of Chemistry and Biochemistry, The Xiphophorus Genetic Stock Center, Texas State University, San Marcos, TX, USA.,Developmental Biochemistry, Biocenter, University of Wuerzburg, Wuerzburg, Germany
| | | |
Collapse
|
8
|
Benestan LM, Rougemont Q, Senay C, Normandeau E, Parent E, Rideout R, Bernatchez L, Lambert Y, Audet C, Parent GJ. Population genomics and history of speciation reveal fishery management gaps in two related redfish species ( Sebastes mentella and Sebastes fasciatus). Evol Appl 2021; 14:588-606. [PMID: 33664797 PMCID: PMC7896722 DOI: 10.1111/eva.13143] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2020] [Revised: 08/27/2020] [Accepted: 09/01/2020] [Indexed: 12/18/2022] Open
Abstract
Understanding the processes shaping population structure and reproductive isolation of marine organisms can improve their management and conservation. Using genomic markers combined with estimation of individual ancestries, assignment tests, spatial ecology, and demographic modeling, we (i) characterized the contemporary population structure, (ii) assessed the influence of space, fishing depth, and sampling years on contemporary distribution, and (iii) reconstructed the speciation history of two cryptic redfish species, Sebastes mentella and S. fasciatus. We genotyped 860 individuals in the Northwest Atlantic Ocean using 24,603 filtered single nucleotide polymorphisms (SNPs). Our results confirmed the clear genetic distinctiveness of the two species and identified three ecotypes within S. mentella and five populations in S. fasciatus. Multivariate analyses highlighted the influence of spatial distribution and depth on the overall genomic variation, while demographic modeling revealed that secondary contact models best explained inter- and intragenomic divergence. These species, ecotypes, and populations can be considered as a rare and wide continuum of genomic divergence in the marine environment. This acquired knowledge pertaining to the evolutionary processes driving population divergence and reproductive isolation will help optimizing the assessment of demographic units and possibly to refine fishery management units.
Collapse
Affiliation(s)
- Laura M. Benestan
- CEFEUniv Montpellier, CNRS, EPHE‐PSL UniversityIRD, Univ Paul Valéry Montpellier 3MontpellierFrance
| | - Quentin Rougemont
- Département de biologie, Institut de Biologie Intégrative et des Systèmes (IBIS)Université LavalQuébecQCCanada
| | - Caroline Senay
- Fisheries and Oceans CanadaMaurice‐Lamontagne InstituteMont‐JoliQCCanada
| | - Eric Normandeau
- Département de biologie, Institut de Biologie Intégrative et des Systèmes (IBIS)Université LavalQuébecQCCanada
| | - Eric Parent
- Fisheries and Oceans CanadaMaurice‐Lamontagne InstituteMont‐JoliQCCanada
| | - Rick Rideout
- Fisheries and Oceans CanadaNorthwest Atlantic Fisheries CentreN.L.St. John’sCanada
| | - Louis Bernatchez
- Département de biologie, Institut de Biologie Intégrative et des Systèmes (IBIS)Université LavalQuébecQCCanada
| | - Yvan Lambert
- Fisheries and Oceans CanadaMaurice‐Lamontagne InstituteMont‐JoliQCCanada
| | - Céline Audet
- Institut des sciences de la mer de RimouskiUniversité du Québec à RimouskiRimouskiQCCanada
| | | |
Collapse
|
9
|
Ewart KM, Johnson RN, Joseph L, Ogden R, Frankham GJ, Lo N. Phylogeography of the iconic Australian pink cockatoo, Lophochroa leadbeateri. Biol J Linn Soc Lond 2021. [DOI: 10.1093/biolinnean/blaa225] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/14/2022]
Abstract
Abstract
The pink cockatoo (Lophochroa leadbeateri; or Major Mitchell’s cockatoo) is one of Australia’s most iconic bird species. Two subspecies based on morphology are separated by a biogeographical divide, the Eyrean Barrier. Testing the genetic basis for this subspecies delineation, clarifying barriers to gene flow and identifying any cryptic genetic diversity will likely have important implications for conservation and management. Here, we used genome-wide single nucleotide polymorphisms (SNPs) and mitochondrial DNA data to conduct the first range-wide genetic assessment of the species. The aims were to investigate the phylogeography of the pink cockatoo, to characterize conservation units and to reassess subspecies boundaries. We found consistent but weak genetic structure between the two subspecies based on nuclear SNPs. However, phylogenetic analysis of nuclear SNPs and mitochondrial DNA sequence data did not recover reciprocally monophyletic groups, indicating incomplete evolutionary separation between the subspecies. Consequently, we have proposed that the two currently recognized subspecies be treated as separate management units rather than evolutionarily significant units. Given that poaching is suspected to be a threat to this species, we assessed the utility of our data for wildlife forensic applications. We demonstrated that a subspecies identification test could be designed using as few as 20 SNPs.
Collapse
Affiliation(s)
- Kyle M Ewart
- School of Life and Environmental Sciences, The University of Sydney, Sydney, NSW, Australia
- Australian Centre for Wildlife Genomics, Australian Museum Research Institute, Sydney, NSW, Australia
| | - Rebecca N Johnson
- School of Life and Environmental Sciences, The University of Sydney, Sydney, NSW, Australia
- Australian Centre for Wildlife Genomics, Australian Museum Research Institute, Sydney, NSW, Australia
| | - Leo Joseph
- Australian National Wildlife Collection, National Research Collections Australia, CSIRO, Canberra, ACT, Australia
| | - Rob Ogden
- Royal (Dick) School of Veterinary Studies and the Roslin Institute, University of Edinburgh, Edinburgh, UK
| | - Greta J Frankham
- Australian Centre for Wildlife Genomics, Australian Museum Research Institute, Sydney, NSW, Australia
- Centre for Forensic Science, University of Technology Sydney, Broadway, NSW, Australia
| | - Nathan Lo
- School of Life and Environmental Sciences, The University of Sydney, Sydney, NSW, Australia
| |
Collapse
|
10
|
Double-digest RAD-sequencing: do pre- and post-sequencing protocol parameters impact biological results? Mol Genet Genomics 2021; 296:457-471. [PMID: 33469716 DOI: 10.1007/s00438-020-01756-9] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2020] [Accepted: 12/14/2020] [Indexed: 02/06/2023]
Abstract
Next-generation sequencing technologies have opened a new era of research in population genetics. Following these new sequencing opportunities, the use of restriction enzyme-based genotyping techniques, such as restriction site-associated DNA sequencing (RAD-seq) or double-digest RAD-sequencing (ddRAD-seq), has dramatically increased in the last decade. From DNA sampling to SNP calling, the laboratory and bioinformatic parameters of enzyme-based techniques have been investigated in the literature. However, the impact of those parameters on downstream analyses and biological results remains less documented. In this study, we investigated the effects of sevral pre- and post-sequencing settings on ddRAD-seq results for two biological systems: a complex of butterfly species (Coenonympha sp.) and several populations of common beech (Fagus sylvatica). Our results suggest that pre-sequencing parameters (i.e., DNA quantity, number of PCR cycles during library preparation) have a significant impact on the number of recovered reads and SNPs, on the number of unique alleles and on individual heterozygosity. In the same way, we found that post-sequencing settings (i.e., clustering and minimum coverage thresholds) influenced loci reconstruction (e.g., number of loci, mean coverage) and SNP calling (e.g., number of SNPs; heterozygosity) but had only a marginal impact on downstream analyses (e.g., measure of genetic differentiation, estimation of individual admixture, and demographic inferences). In addition, replication analyses confirmed the reproducibility of the ddRAD-seq procedure. Overall, this study assesses the degree of sensitivity of ddRAD-seq data to pre- and post-sequencing protocols, and illustrates its robustness when studying population genetics.
Collapse
|
11
|
Nazareno AG, Knowles LL, Dick CW, Lohmann LG. By Animal, Water, or Wind: Can Dispersal Mode Predict Genetic Connectivity in Riverine Plant Species? FRONTIERS IN PLANT SCIENCE 2021; 12:626405. [PMID: 33643353 PMCID: PMC7907645 DOI: 10.3389/fpls.2021.626405] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/05/2020] [Accepted: 01/06/2021] [Indexed: 05/08/2023]
Abstract
Seed dispersal is crucial to gene flow among plant populations. Although the effects of geographic distance and barriers to gene flow are well studied in many systems, it is unclear how seed dispersal mediates gene flow in conjunction with interacting effects of geographic distance and barriers. To test whether distinct seed dispersal modes (i.e., hydrochory, anemochory, and zoochory) have a consistent effect on the level of genetic connectivity (i.e., gene flow) among populations of riverine plant species, we used unlinked single-nucleotide polymorphisms (SNPs) for eight co-distributed plant species sampled across the Rio Branco, a putative biogeographic barrier in the Amazon basin. We found that animal-dispersed plant species exhibited higher levels of genetic diversity and lack of inbreeding as a result of the stronger genetic connectivity than plant species whose seeds are dispersed by water or wind. Interestingly, our results also indicated that the Rio Branco facilitates gene dispersal for all plant species analyzed, irrespective of their mode of dispersal. Even at a small spatial scale, our findings suggest that ecology rather than geography play a key role in shaping the evolutionary history of plants in the Amazon basin. These results may help improve conservation and management policies in Amazonian riparian forests, where degradation and deforestation rates are high.
Collapse
Affiliation(s)
- Alison G. Nazareno
- Departamentos de Botânica, Universidade de São Paulo, São Paulo, Brazil
- Department of Genetics, Ecology and Evolution, Federal University of Minas Gerais, Belo Horizonte, Brazil
- Department of Ecology and Evolutionary Biology, University of Michigan, Ann Arbor, MI, United States
- *Correspondence: Alison G. Nazareno,
| | - L. Lacey Knowles
- Department of Ecology and Evolutionary Biology, University of Michigan, Ann Arbor, MI, United States
| | - Christopher W. Dick
- Department of Ecology and Evolutionary Biology, University of Michigan, Ann Arbor, MI, United States
- Smithsonian Tropical Research Institute, Panama City, Panama
| | - Lúcia G. Lohmann
- Departamentos de Botânica, Universidade de São Paulo, São Paulo, Brazil
- Lúcia G. Lohmann,
| |
Collapse
|
12
|
Vaux F, Gemmell MR, Hills SFK, Marshall BA, Beu AG, Crampton JS, Trewick SA, Morgan-Richards M. Lineage Identification Affects Estimates of Evolutionary Mode in Marine Snails. Syst Biol 2020; 69:1106-1121. [PMID: 32163159 DOI: 10.1093/sysbio/syaa018] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/01/2019] [Revised: 02/13/2020] [Accepted: 03/09/2020] [Indexed: 11/14/2022] Open
Abstract
In order to study evolutionary pattern and process, we need to be able to accurately identify species and the evolutionary lineages from which they are derived. Determining the concordance between genetic and morphological variation of living populations, and then directly comparing extant and fossil morphological data, provides a robust approach for improving our identification of lineages through time. We investigate genetic and shell morphological variation in extant species of Penion marine snails from New Zealand, and extend this analysis into deep time using fossils. We find that genetic and morphological variation identify similar patterns and support most currently recognized extant species. However, some taxonomic over-splitting is detected due to shell size being a poor trait for species delimitation, and we identify incorrect assignment of some fossil specimens. We infer that a single evolutionary lineage (Penion sulcatus) has existed for 22 myr, with most aspects of shell shape and shell size evolving under a random walk. However, by removing samples previously classified as the extinct species P. marwicki, we instead detect morphological stasis for one axis of shell shape variation. This result demonstrates how lineage identification can change our perception of evolutionary pattern and process. [Genotyping by sequencing; geometric morphometrics; morphological evolution; Neogastropoda; phenotype; speciation; stasis.].
Collapse
Affiliation(s)
- Felix Vaux
- Wildlife and Ecology Group, School of Agriculture and Environment, Massey University, Private Bag 11-222, Palmerston North 4410, New Zealand.,Department of Fisheries and Wildlife, Coastal Oregon Marine Experiment Station, Hatfield Marine Science Center, Oregon State University, 2030 SE Marine Science Dr, Newport, OR 97365, USA.,Department of Zoology, University of Otago, 340 Great King Street, Dunedin 9016, Otago, New Zealand
| | - Michael R Gemmell
- Wildlife and Ecology Group, School of Agriculture and Environment, Massey University, Private Bag 11-222, Palmerston North 4410, New Zealand
| | - Simon F K Hills
- Wildlife and Ecology Group, School of Agriculture and Environment, Massey University, Private Bag 11-222, Palmerston North 4410, New Zealand
| | - Bruce A Marshall
- Museum of New Zealand Te Papa Tongarewa, Wellington, New Zealand
| | - Alan G Beu
- GNS Science, PO Box 30-368, Lower Hutt 5011, New Zealand
| | - James S Crampton
- School of Geography, Environment & Earth Sciences, Victoria University of Wellington, PO Box 600, Wellington 6012, New Zealand
| | - Steven A Trewick
- Wildlife and Ecology Group, School of Agriculture and Environment, Massey University, Private Bag 11-222, Palmerston North 4410, New Zealand
| | - Mary Morgan-Richards
- Wildlife and Ecology Group, School of Agriculture and Environment, Massey University, Private Bag 11-222, Palmerston North 4410, New Zealand
| |
Collapse
|
13
|
Carrier E, Ferchaud AL, Normandeau E, Sirois P, Bernatchez L. Estimating the contribution of Greenland Halibut ( Reinhardtius hippoglossoides) stocks to nurseries by means of genotyping-by-sequencing: Sex and time matter. Evol Appl 2020; 13:2155-2167. [PMID: 33005216 PMCID: PMC7513701 DOI: 10.1111/eva.12979] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/03/2019] [Revised: 03/19/2020] [Accepted: 04/01/2020] [Indexed: 12/14/2022] Open
Abstract
Identification of stocks and quantification of their relative contribution to recruitment are major objectives toward improving the management and conservation of marine exploited species. Next-generation sequencing allows for thousands of genomic markers to be analyzed, which provides the resolution needed to address these questions in marine species with weakly differentiated populations. Greenland Halibut (Reinhardtius hippoglossoides) is one of the most important exploited demersal species throughout the North Atlantic, in particular in the Gulf of St. Lawrence, Canada. There, two nurseries are known, the St. Lawrence Estuary and the northern Anticosti Island, but their contribution to the renewal of stocks remains unknown. The goals of this study were (a) to document the genetic structure and (b) to estimate the contribution of the different identified breeding stocks to nurseries. We sampled 100 juveniles per nursery and 50 adults from seven sites ranging from Saguenay Fjord to offshore Newfoundland, with some sites sampled over two consecutive years in order to evaluate the temporal stability of the contribution. Our results show that after removing sex-linked markers, the Estuary/Gulf of St. Lawrence represents a single stock which is genetically distinct from the Atlantic around Newfoundland (F ST = 0.00146, p-value = .001). Population assignment showed that recruitment in both nurseries is predominantly associated with the St. Lawrence stock. However, we found that the relative contribution of both stocks to the nurseries is temporally variable with 1% contribution of the Newfoundland stock one year but up to 33% for the second year, which may be caused by year-to-year variation in larval transport into the Gulf of St. Lawrence. This study serves as a model for the identification of stocks for fisheries resources in a context where few barriers to dispersal occur, in addition to demonstrating the importance of considering sex-linked markers and temporal replicates in studies of population genomics.
Collapse
Affiliation(s)
- Emilie Carrier
- Institut de biologie intégrative et des systèmes (IBIS) Université Laval Québec City QC Canada
| | - Anne-Laure Ferchaud
- Institut de biologie intégrative et des systèmes (IBIS) Université Laval Québec City QC Canada
| | - Eric Normandeau
- Institut de biologie intégrative et des systèmes (IBIS) Université Laval Québec City QC Canada
| | - Pascal Sirois
- Département des sciences fondamentales Université du Québec à Chicoutimi Chicoutimi QC Canada
| | - Louis Bernatchez
- Institut de biologie intégrative et des systèmes (IBIS) Université Laval Québec City QC Canada
| |
Collapse
|
14
|
Ewart KM, Lo N, Ogden R, Joseph L, Ho SYW, Frankham GJ, Eldridge MDB, Schodde R, Johnson RN. Phylogeography of the iconic Australian red-tailed black-cockatoo (Calyptorhynchus banksii) and implications for its conservation. Heredity (Edinb) 2020; 125:85-100. [PMID: 32398870 PMCID: PMC7426920 DOI: 10.1038/s41437-020-0315-y] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/30/2019] [Revised: 04/08/2020] [Accepted: 04/08/2020] [Indexed: 01/31/2023] Open
Abstract
Advances in sequencing technologies have revolutionized wildlife conservation genetics. Analysis of genomic data sets can provide high-resolution estimates of genetic structure, genetic diversity, gene flow, and evolutionary history. These data can be used to characterize conservation units and to effectively manage the genetic health of species in a broad evolutionary context. Here we utilize thousands of genome-wide single-nucleotide polymorphisms (SNPs) and mitochondrial DNA to provide the first genetic assessment of the Australian red-tailed black-cockatoo (Calyptorhynchus banksii), a widespread bird species comprising populations of varying conservation concern. We identified five evolutionarily significant units, which are estimated to have diverged during the Pleistocene. These units are only partially congruent with the existing morphology-based subspecies taxonomy. Genetic clusters inferred from mitochondrial DNA differed from those based on SNPs and were less resolved. Our study has a range of conservation and taxonomic implications for this species. In particular, we provide advice on the potential genetic rescue of the Endangered and restricted-range subspecies C. b. graptogyne, and propose that the western C. b. samueli population is diagnosable as a separate subspecies. The results of our study highlight the utility of considering the phylogeographic relationships inferred from genome-wide SNPs when characterizing conservation units and management priorities, which is particularly relevant as genomic data sets become increasingly accessible.
Collapse
Affiliation(s)
- Kyle M Ewart
- School of Life and Environmental Sciences, University of Sydney, Sydney, NSW, Australia.
- Australian Centre for Wildlife Genomics, Australian Museum Research Institute, Sydney, NSW, Australia.
| | - Nathan Lo
- School of Life and Environmental Sciences, University of Sydney, Sydney, NSW, Australia
| | - Rob Ogden
- Royal (Dick) School of Veterinary Studies and the Roslin Institute, Easter Bush Campus, University of Edinburgh, Edinburgh, EH25 9RG, UK
| | - Leo Joseph
- Australian National Wildlife Collection, CSIRO, Canberra, ACT, Australia
| | - Simon Y W Ho
- School of Life and Environmental Sciences, University of Sydney, Sydney, NSW, Australia
| | - Greta J Frankham
- Australian Centre for Wildlife Genomics, Australian Museum Research Institute, Sydney, NSW, Australia
| | - Mark D B Eldridge
- Australian Centre for Wildlife Genomics, Australian Museum Research Institute, Sydney, NSW, Australia
| | - Richard Schodde
- Australian National Wildlife Collection, CSIRO, Canberra, ACT, Australia
| | - Rebecca N Johnson
- School of Life and Environmental Sciences, University of Sydney, Sydney, NSW, Australia
- Australian Centre for Wildlife Genomics, Australian Museum Research Institute, Sydney, NSW, Australia
| |
Collapse
|
15
|
Xiao JH, Ding X, Li L, Ma H, Ci XQ, van der Merwe M, Conran JG, Li J. Miocene diversification of a golden-thread nanmu tree species ( Phoebe zhennan, Lauraceae) around the Sichuan Basin shaped by the East Asian monsoon. Ecol Evol 2020; 10:10543-10557. [PMID: 33072279 PMCID: PMC7548194 DOI: 10.1002/ece3.6710] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2020] [Revised: 07/29/2020] [Accepted: 07/30/2020] [Indexed: 01/04/2023] Open
Abstract
Understanding the role of climate changes and geography as drivers of population divergence and speciation is a long‐standing goal of evolutionary biology and can inform conservation. In this study, we used restriction site‐associated DNA sequencing (RAD‐seq) to evaluate genetic diversity, population structure, and infer demographic history of the endangered tree, Phoebe zhennan which is distributed around the Sichuan Basin. Genomic patterns revealed two distinct clusters, each largely confined to the West and East. Despite sympatry of the two genomic clusters at some sites, individuals show little or no evidence of genomic introgression. Demographic modeling supported an initial divergence time between the West and East lineages at ~15.08 Ma with further diversification within the West lineage at ~7.12 Ma. These times largely coincide with the two independent intensifications of the East Asian monsoon that were initiated during the middle (Langhian) and late Miocene (Messinian), respectively. These results suggest that the Miocene intensification phases of the East Asian monsoon played a pivotal role in shaping the current landscape‐level patterns of genetic diversity within P. zhennan, as has been found for the interspecific divergence of other subtropical Chinese plants. Based on isolation‐by‐distance and species distribution modeling, we hypothesize that P. zhennan followed a ring diversification which was facilitated by the Sichuan Basin acting as barrier to gene flow. In situ and ex situ conservation management plans should consider the results obtained in this study to help secure the future of this beautiful and culturally significant endangered tree.
Collapse
Affiliation(s)
- Jian-Hua Xiao
- Plant Phylogenetics and Conservation Group Centre for Integrative Conservation Xishuangbanna Tropical Botanical Garden Chinese Academy of Sciences Kunming China.,University of Chinese Academy of Sciences Beijing China
| | - Xin Ding
- Department of Landscape Architecture Guangdong Eco-engineering Polytechnic Guangzhou China
| | - Lang Li
- Plant Phylogenetics and Conservation Group Centre for Integrative Conservation Xishuangbanna Tropical Botanical Garden Chinese Academy of Sciences Kunming China.,Center of Conservation Biology Core Botanical Gardens Chinese Academy of Sciences Mengla China
| | - Hui Ma
- Plant Phylogenetics and Conservation Group Centre for Integrative Conservation Xishuangbanna Tropical Botanical Garden Chinese Academy of Sciences Kunming China
| | - Xiu-Qin Ci
- Plant Phylogenetics and Conservation Group Centre for Integrative Conservation Xishuangbanna Tropical Botanical Garden Chinese Academy of Sciences Kunming China.,Center of Conservation Biology Core Botanical Gardens Chinese Academy of Sciences Mengla China
| | - Marlien van der Merwe
- Research Centre for Ecological Resilience Australian Institute of Botanical Science The Royal Botanic Garden Sydney Sydney NSW Australia
| | - John G Conran
- Australian Centre for Evolutionary Biology and Biodiversity (ACEBB) Sprigg Geobiology Centre (SGC) School of Biological Sciences The University of Adelaide Adelaide SA Australia
| | - Jie Li
- Plant Phylogenetics and Conservation Group Centre for Integrative Conservation Xishuangbanna Tropical Botanical Garden Chinese Academy of Sciences Kunming China.,Center of Conservation Biology Core Botanical Gardens Chinese Academy of Sciences Mengla China
| |
Collapse
|
16
|
Marandel F, Charrier G, Lamy J, Le Cam S, Lorance P, Trenkel VM. Estimating effective population size using RADseq: Effects of SNP selection and sample size. Ecol Evol 2020; 10:1929-1937. [PMID: 32128126 PMCID: PMC7042749 DOI: 10.1002/ece3.6016] [Citation(s) in RCA: 29] [Impact Index Per Article: 7.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2019] [Revised: 12/13/2019] [Accepted: 12/23/2019] [Indexed: 01/16/2023] Open
Abstract
Effective population size (Ne ) is a key parameter of population genetics. However, N e remains challenging to estimate for natural populations as several factors are likely to bias estimates. These factors include sampling design, sequencing method, and data filtering. One issue inherent to the restriction site-associated DNA sequencing (RADseq) protocol is missing data and SNP selection criteria (e.g., minimum minor allele frequency, number of SNPs). To evaluate the potential impact of SNP selection criteria on Ne estimates (Linkage Disequilibrium method) we used RADseq data for a nonmodel species, the thornback ray. In this data set, the inbreeding coefficient F IS was positively correlated with the amount of missing data, implying data were missing nonrandomly. The precision of Ne estimates decreased with the number of SNPs. Mean Ne estimates (averaged across 50 random data sets with2000 SNPs) ranged between 237 and 1784. Increasing the percentage of missing data from 25% to 50% increased Ne estimates between 82% and 120%, while increasing the minor allele frequency (MAF) threshold from 0.01 to 0.1 decreased estimates between 71% and 75%. Considering these effects is important when interpreting RADseq data-derived estimates of effective population size in empirical studies.
Collapse
Affiliation(s)
| | - Grégory Charrier
- Laboratoire des Sciences de l’Environnement Marin (LEMAR, UMR 6539 CNRS/IRD/UBO/Ifremer)Université de Bretagne OccidentaleInstitut Universitaire Européen de la MerPlouzanéFrance
| | - Jean‐Baptiste Lamy
- IfremerGénétique et Pathologie des Mollusques Marin (SG2M‐LGPMM)La TrembladeFrance
| | - Sabrina Le Cam
- Laboratoire des Sciences de l’Environnement Marin (LEMAR, UMR 6539 CNRS/IRD/UBO/Ifremer)Université de Bretagne OccidentaleInstitut Universitaire Européen de la MerPlouzanéFrance
- IfremerGénétique et Pathologie des Mollusques Marin (SG2M‐LGPMM)La TrembladeFrance
| | - Pascal Lorance
- IfremerEcologie et Modèles pour l’HalieutiqueNantesFrance
| | | |
Collapse
|
17
|
Graham CF, Boreham DR, Manzon RG, Stott W, Wilson JY, Somers CM. How "simple" methodological decisions affect interpretation of population structure based on reduced representation library DNA sequencing: A case study using the lake whitefish. PLoS One 2020; 15:e0226608. [PMID: 31978053 PMCID: PMC6980518 DOI: 10.1371/journal.pone.0226608] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/25/2019] [Accepted: 12/01/2019] [Indexed: 12/30/2022] Open
Abstract
Reduced representation (RRL) sequencing approaches (e.g., RADSeq, genotyping by sequencing) require decisions about how much to invest in genome coverage and sequencing depth, as well as choices of values for adjustable bioinformatics parameters. To empirically explore the importance of these “simple” methodological decisions, we generated two independent sequencing libraries for the same 142 individual lake whitefish (Coregonus clupeaformis) using a nextRAD RRL approach: (1) a larger number of loci at low sequencing depth based on a 9mer (library A); and (2) fewer loci at higher sequencing depth based on a 10mer (library B). The fish were selected from populations with different levels of expected genetic subdivision. Each library was analyzed using the STACKS pipeline followed by three types of population structure assessment (FST, DAPC and ADMIXTURE) with iterative increases in the stringency of sequencing depth and missing data requirements, as well as more specific a priori population maps. Library B was always able to resolve strong population differentiation in all three types of assessment regardless of the selected parameters, largely due to retention of more loci in analyses. In contrast, library A produced more variable results; increasing the minimum sequencing depth threshold (-m) resulted in a reduced number of retained loci, and therefore lost resolution at high -m values for FST and ADMIXTURE, but not DAPC. When detecting fine population differentiation, the population map influenced the number of loci and missing data, which generated artefacts in all downstream analyses tested. Similarly, when examining fine scale population subdivision, library B was robust to changing parameters but library A lost resolution depending on the parameter set. We used library B to examine actual subdivision in our study populations. All three types of analysis found complete subdivision among populations in Lake Huron, ON and Dore Lake, SK, Canada using 10,640 SNP loci. Weak population subdivision was detected in Lake Huron with fish from sites in the north-west, Search Bay, North Point and Hammond Bay, showing slight differentiation. Overall, we show that apparently simple decisions about library construction and bioinformatics parameters can have important impacts on the interpretation of population subdivision. Although potentially more costly on a per-locus basis, early investment in striking a balance between the number of loci and sequencing effort is well worth the reduced genomic coverage for population genetics studies. More conservative stringency settings on STACKS parameters lead to a final dataset that was more consistent and robust when examining both weak and strong population differentiation. Overall, we recommend that researchers approach “simple” methodological decisions with caution, especially when working on non-model species for the first time.
Collapse
Affiliation(s)
- Carly F. Graham
- Department of Biology, University of Regina, Regina, Saskatchewan, Canada
| | - Douglas R. Boreham
- Medical Sciences, Northern Ontario School of Medicine, Greater Sudbury, Ontario, Canada
| | - Richard G. Manzon
- Department of Biology, University of Regina, Regina, Saskatchewan, Canada
| | - Wendylee Stott
- Department of Fisheries and Wildlife, Michigan State University, East Lansing, Michigan, USA
| | - Joanna Y. Wilson
- Department of Biology, McMaster University, Hamilton, Ontario, Canada
| | | |
Collapse
|
18
|
Baltazar-Soares M, Blanchet S, Cote J, Tarkan AS, Záhorská E, Gozlan RE, Eizaguirre C. Genomic footprints of a biological invasion: Introduction from Asia and dispersal in Europe of the topmouth gudgeon (Pseudorasbora parva). Mol Ecol 2019; 29:71-85. [PMID: 31755610 PMCID: PMC7003831 DOI: 10.1111/mec.15313] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/01/2019] [Revised: 11/11/2019] [Accepted: 11/17/2019] [Indexed: 12/14/2022]
Abstract
Facilitated by the intensification of global trading, the introduction and dispersal of species to areas in which they are historically non-native is nowadays common. From an evolutionary standpoint, invasions are paradoxical: not only non-native environments could be different from native ones for which introduced individuals would be ill-adapted, but also small founding population size should be associated with reduced adaptive potential. As such, biological invasions are considered valuable real-time evolutionary experiments. Here, we investigated the population structure and adaptive potential of the highly invasive topmouth gudgeon (Pseudorasbora parva) across Europe and East Asia. We RAD-sequenced 301 specimens from sixteen populations and three distinct within-catchment invaded regions as well as two locations in the native range. With 13,785 single nucleotide polymorphisms, we provide conclusive evidence for a genome-wide signature of two distinct invasion events, in Slovakia and Turkey, each originating from a specific area in the native range. A third invaded area, in France, appears to be the result of dispersal within the invasive range. Few loci showed signs of selection, the vast majority of which being identified in the Slovakian region. Functional annotation suggests that faster early stage development, resistance to pollution and immunocompetence contribute to the invasion success of the local habitats. By showing that populations in the invasive range have different evolutionary histories, our study reinforces the idea that populations, rather than species, are the units to consider in invasion biology.
Collapse
Affiliation(s)
| | - Simon Blanchet
- CNRS, Station d'Ecologie Théorique et Expérimentale (SETE), Moulis, France
| | - Julien Cote
- UMR5174 (Laboratoire Evolution et Diversité Biologique), CNRS, University Toulouse III Paul Sabatier, Toulouse, France
| | - Ali S Tarkan
- Faculty of Fisheries, Muğla Sıtkı Koçman University, Kötekli, Muğla, Turkey.,Department of Ecology and Vertebrate Zoology, Faculty of Biology and Environmental Protection, University of Łódź, Łódź, Poland
| | - Eva Záhorská
- Faculty of Natural Sciences, Department of Ecology, Comenius University, Bratislava, Slovakia
| | - Rodolphe E Gozlan
- ISEM, Université de Montpellier, CNRS, IRD, EPHE, Montpellier, France
| | - Christophe Eizaguirre
- School of Chemical and Biological Sciences, Queen Mary University of London, London, UK
| |
Collapse
|
19
|
Mynhardt S, Bennett NC, Bloomer P. New insights from RADseq data on differentiation in the Hottentot golden mole species complex from South Africa. Mol Phylogenet Evol 2019; 143:106667. [PMID: 31676418 DOI: 10.1016/j.ympev.2019.106667] [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: 03/15/2019] [Revised: 10/28/2019] [Accepted: 10/28/2019] [Indexed: 10/25/2022]
Abstract
Golden moles (Family Chrysochloridae) are small subterranean mammals, endemic to sub-Saharan Africa, and many of the 21 species are listed as threatened on the IUCN Red List. Most species have highly restricted ranges; however two species, the Hottentot golden mole (Amblysomus hottentotus) and the Cape golden mole (Chrysochloris asiatica) have relatively wide ranges. We recently uncovered cryptic diversity within A. hottentotus, through a phylogeographic analysis of this taxon using two mitochondrial gene regions and a nuclear intron. To further investigate this cryptic diversity, we generated nuclear SNP data from across the genome of A. hottentotus, by means of double-digest restriction-site associated DNA sequencing (ddRADSeq), and mapped reads to the Cape golden mole genome. We conducted a phylogenetic analysis and investigated population differentiation. Our results support the distinctiveness of A. h. meesteri. Furthermore, we provide evidence from nuclear SNPs in support of our previous finding that Central coastal samples represent a unique cryptic lineage that is highly divergent from A. h. pondoliae farther south. Although mtDNA suggests that Umtata may represent a unique lineage sister to A. h. longiceps, mito-nuclear discordance from our RADseq data indicate that these samples may instead be closer to A. h. pondoliae, and therefore may not represent a distinct lineage. We stress the importance of recognizing that understudied populations, such as that of Umtata, may represent populations or ESUs under threat and in need of conservation attention. We present a high-quality filtered SNP dataset, comprising thousands of SNPs, which may serve as a useful resource for future golden mole studies. We have thus added to the growing body of research demonstrating the power and utility of RADseq to investigate population differentiation.
Collapse
Affiliation(s)
- Samantha Mynhardt
- Molecular Ecology and Evolution Programme, Department of Biochemistry, Genetics and Microbiology, University of Pretoria, Private Bag X20, Hatfield 0028, South Africa; Department of Zoology and Entomology, University of Pretoria, Private Bag X20, Hatfield 0028, South Africa.
| | - Nigel C Bennett
- Department of Zoology and Entomology, University of Pretoria, Private Bag X20, Hatfield 0028, South Africa; Mammal Research Institute, Department of Zoology and Entomology, University of Pretoria, Private Bag X20, Hatfield 0028, South Africa.
| | - Paulette Bloomer
- Molecular Ecology and Evolution Programme, Department of Biochemistry, Genetics and Microbiology, University of Pretoria, Private Bag X20, Hatfield 0028, South Africa; Mammal Research Institute, Department of Zoology and Entomology, University of Pretoria, Private Bag X20, Hatfield 0028, South Africa.
| |
Collapse
|
20
|
Ewart KM, Johnson RN, Ogden R, Joseph L, Frankham GJ, Lo N. Museum specimens provide reliable SNP data for population genomic analysis of a widely distributed but threatened cockatoo species. Mol Ecol Resour 2019; 19:1578-1592. [PMID: 31484222 DOI: 10.1111/1755-0998.13082] [Citation(s) in RCA: 28] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/27/2018] [Revised: 08/05/2019] [Accepted: 08/12/2019] [Indexed: 12/20/2022]
Abstract
Natural history museums harbour a plethora of biological specimens which are of potential use in population and conservation genetic studies. Although technical advancements in museum genomics have enabled genome-wide markers to be generated from aged museum specimens, the suitability of these data for robust biological inference is not well characterized. The aim of this study was to test the utility of museum specimens in population and conservation genomics by assessing the biological and technical validity of single nucleotide polymorphism (SNP) data derived from such samples. To achieve this, we generated thousands of SNPs from 47 red-tailed black cockatoo (Calyptorhychus banksii) traditional museum samples (i.e. samples that were not collected with the primary intent of DNA analysis) and 113 fresh tissue samples (cryopreserved liver/muscle) using a restriction site-associated DNA marker approach (DArTseq™ ). Thousands of SNPs were successfully generated from most of the traditional museum samples (with a mean age of 44 years, ranging from 5 to 123 years), although 38% did not provide useful data. These SNPs exhibited higher error rates and contained significantly more missing data compared with SNPs from fresh tissue samples, likely due to considerable DNA fragmentation. However, based on simulation results, the level of genotyping error had a negligible effect on inference of population structure in this species. We did identify a bias towards low diversity SNPs in older samples that appears to compromise temporal inferences of genetic diversity. This study demonstrates the utility of a RADseq-based method to produce reliable genome-wide SNP data from traditional museum specimens.
Collapse
Affiliation(s)
- Kyle M Ewart
- School of Life and Environmental Sciences, The University of Sydney, Sydney, NSW, Australia.,Australian Centre for Wildlife Genomics, Australian Museum Research Institute, Sydney, NSW, Australia
| | - Rebecca N Johnson
- School of Life and Environmental Sciences, The University of Sydney, Sydney, NSW, Australia.,Australian Centre for Wildlife Genomics, Australian Museum Research Institute, Sydney, NSW, Australia
| | - Rob Ogden
- Royal (Dick) School of Veterinary Studies and the Roslin Institute, University of Edinburgh, Edinburgh, UK
| | - Leo Joseph
- Australian National Wildlife Collection, CSIRO, Canberra, Australia
| | - Greta J Frankham
- Australian Centre for Wildlife Genomics, Australian Museum Research Institute, Sydney, NSW, Australia
| | - Nathan Lo
- School of Life and Environmental Sciences, The University of Sydney, Sydney, NSW, Australia
| |
Collapse
|
21
|
McCartney‐Melstad E, Gidiş M, Shaffer HB. An empirical pipeline for choosing the optimal clustering threshold in RADseq studies. Mol Ecol Resour 2019; 19:1195-1204. [DOI: 10.1111/1755-0998.13029] [Citation(s) in RCA: 29] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2017] [Revised: 04/18/2019] [Accepted: 04/22/2019] [Indexed: 11/27/2022]
Affiliation(s)
- Evan McCartney‐Melstad
- Department of Ecology and Evolutionary Biology, La Kretz Center for California Conservation Science, and Institute of the Environment and Sustainability University of California Los Angeles California
| | - Müge Gidiş
- Faculty of Arts and Science, Department of Biochemistry Kütahya Dumlupınar University Kutahya Turkey
| | - H. Bradley Shaffer
- Department of Ecology and Evolutionary Biology, La Kretz Center for California Conservation Science, and Institute of the Environment and Sustainability University of California Los Angeles California
| |
Collapse
|
22
|
Díaz-Arce N, Rodríguez-Ezpeleta N. Selecting RAD-Seq Data Analysis Parameters for Population Genetics: The More the Better? Front Genet 2019; 10:533. [PMID: 31191624 PMCID: PMC6549478 DOI: 10.3389/fgene.2019.00533] [Citation(s) in RCA: 29] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2019] [Accepted: 05/16/2019] [Indexed: 11/25/2022] Open
Abstract
Restriction site-associated DNA sequencing (RAD-seq) has become a powerful and widely used tool in molecular ecology studies as it allows to cost-effectively recover thousands of polymorphic sites across individuals of non-model organisms. However, its successful implementation in population genetics relies on correct data processing that would minimize potential loci-assembly biases and consequent genotyping error rates. RAD-seq data processing when no reference genome is available involves the assembly of hundreds of thousands high-throughput sequencing reads into orthologous loci, for which various key parameter values need to be selected by the researcher. Previous studies exploring the effect of these parameter values found or assumed that a larger number of recovered polymorphic loci is associated with a better assembly. Here, using three RAD-seq datasets from different species, we explore the effect of read filtering, loci assembly and polymorphic site selection on number of markers obtained and genetic differentiation inferred using the Stacks software. We find (i) that recovery of higher numbers of polymorphic loci is not necessarily associated with higher genetic differentiation, (ii) that the presence of PCR duplicates, selected loci assembly parameters and selected SNP filtering parameters affect the number of recovered polymorphic loci and degree of genetic differentiation, and (iii) that this effect is different in each dataset, meaning that defining a systematic universal protocol for RAD-seq data analysis may lead to missing relevant information about population differentiation.
Collapse
|
23
|
Nazareno AG, Dick CW, Lohmann LG. Tangled banks: A landscape genomic evaluation of Wallace's Riverine barrier hypothesis for three Amazon plant species. Mol Ecol 2019; 28:980-997. [PMID: 30450714 DOI: 10.1111/mec.14948] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/05/2017] [Revised: 08/18/2018] [Accepted: 10/02/2018] [Indexed: 12/20/2022]
Abstract
Wallace's Riverine Barrier hypothesis is one of the earliest biogeographic explanations for Amazon speciation, but it has rarely been tested in plants. In this study, we used three woody Amazonian plant species to evaluate Wallace's Hypothesis using tools of landscape genomics. We generated unlinked single-nucleotide polymorphism (SNP) data from the nuclear genomes of 234 individuals (78 for each plant species) across 13 sampling sites along the Rio Branco, Brazil, for Amphirrhox longifolia (8,075 SNPs), Psychotria lupulina (9,501 SNPs) and Passiflora spinosa (14,536 SNPs). Although significantly different migration rates were estimated between species, the population structure data do not support the hypothesis that the Rio Branco-an allopatric barrier for primates and birds-is a significant genetic barrier for Amphirrhox longifolia, Passiflora spinosa or Psychotria lupulina. Overall, we demonstrated that medium-sized rivers in the Amazon Basin, such as the Rio Branco, are permeable barriers to gene flow for animal-dispersed and animal-pollinated plant species.
Collapse
Affiliation(s)
- Alison G Nazareno
- Departamento de Botânica, Universidade de São Paulo, São Paulo, Brazil
| | - Christopher W Dick
- Department of Ecology and Evolutionary Biology, University of Michigan, Ann Arbor, Michigan.,Smithsonian Tropical Research Institute, Panama City, Panama
| | - Lúcia G Lohmann
- Departamento de Botânica, Universidade de São Paulo, São Paulo, Brazil
| |
Collapse
|
24
|
Population Genomics Applied to Fishery Management and Conservation. POPULATION GENOMICS 2019. [DOI: 10.1007/13836_2019_66] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/24/2022]
|
25
|
Nadukkalam Ravindran P, Bentzen P, Bradbury IR, Beiko RG. RADProc: A computationally efficient de novo locus assembler for population studies using RADseq data. Mol Ecol Resour 2018; 19:272-282. [PMID: 30312001 DOI: 10.1111/1755-0998.12954] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2018] [Revised: 09/09/2018] [Accepted: 09/26/2018] [Indexed: 11/30/2022]
Abstract
Restriction site-associated DNA sequencing (RADseq) is a powerful tool for genotyping of individuals, but the identification of loci and assignment of sequence reads is a crucial and often challenging step. The optimal parameter settings for a given de novo RADseq assembly vary between data sets and can be difficult and computationally expensive to determine. Here, we introduce RADProc, a software package that uses a graph data structure to represent all sequence reads and their similarity relationships. Storing sequence-comparison results in a graph eliminates unnecessary and redundant sequence similarity calculations. De novo locus formation for a given parameter set can be performed on the precomputed graph, making parameter sweeps far more efficient. RADProc also uses a clustering approach for faster nucleotide-distance calculation. The performance of RADProc compares favourably with that of the widely used Stacks software. The run-time comparisons between RADProc and Stacks for 32 different parameter settings using 20 green-crab (Carcinus maenas) samples showed that RADProc took as little as 2 hr 40 min compared to 78 hr by Stacks, while 16 brown trout (Salmo trutta L.) samples were processed by RADProc and Stacks in 23 and 263 hr, respectively. Comparisons of the de novo loci formed, and catalog built using both the methods demonstrate that the improvement in processing speeds achieved by RADProc does not affect much the actual loci formed and the results of downstream analyses based on those loci.
Collapse
Affiliation(s)
| | - Paul Bentzen
- Marine Gene Probe Laboratory, Department of Biology, Dalhousie University, Halifax, Nova Scotia, Canada
| | - Ian R Bradbury
- Salmonids Section, Science Branch, Department of Fisheries and Oceans Canada, St. John's, Newfoundland, Canada
| | - Robert G Beiko
- Faculty of Computer Science, Dalhousie University, Halifax, Nova Scotia, Canada
| |
Collapse
|
26
|
Catchen J. When structure leads to sex: Untangling signals in population genetic data sets. Mol Ecol 2018; 26:6763-6766. [PMID: 29282805 DOI: 10.1111/mec.14422] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/26/2017] [Revised: 11/20/2017] [Accepted: 11/04/2017] [Indexed: 11/28/2022]
Abstract
A robust signal of population structure often provides the first glimpse into the evolutionary history of a species and its populations. In this issue of Molecular Ecology, new work from Louis Bernatchez's group (Benestan et al., ) starts with an investigation of apparent structure in two marine species and concludes with an identification of sex-linked genes, and in the process provides a model for robust analysis. Structure is the genetic signal left by natural selection as well as by neutral processes like migration and gene flow. Neutral areas of the genome can reveal the geographical relationships and related gene flow between populations over time and space, while selection can resist the natural genomic turnover created by recombination and generate adaptive structure between populations that can be detected. However, artefacts in a data set can easily hide the true signal of structure; mutation, whether it is a true appearance of a recent, minor allele, or more commonly, an error in SNP calling or molecular library construction, can easily conceal patterns of population structure (e.g., geographical structure in mackerel, Rodriguez-Ezpeleta et al. ()). A demographic structure that results from the most "forceful" evolutionary processes can overwhelm another signal generated by other, unrelated phenotypes. For example, the structure among diverged freshwater and marine threespine stickleback populations results from such strong selection and linkage disequilibrium across the genome that it impairs the ability to disentangle the genetic basis of particular evolved morphological traits (e.g., opercle development, Alligood ()). Finally, there might be conflicting inferences for what underlies structure patterns. Structure may be created by differential patterns of meiotic recombination, and genetic maps are a reliable means for identifying genomic regions that resist recombination. But, without additional information (Anderson et al., ), it can be difficult to distinguish the recombination-suppressing effect of a segregating genomic inversion (Small et al., ) from that of sex-linked selection.
Collapse
Affiliation(s)
- Julian Catchen
- Department of Animal Biology, University of Illinois at Urbana-Champaign, Urbana, IL, USA
| |
Collapse
|
27
|
Nadukkalam Ravindran P, Bentzen P, Bradbury IR, Beiko RG. PMERGE: Computational filtering of paralogous sequences from RAD-seq data. Ecol Evol 2018; 8:7002-7013. [PMID: 30073062 PMCID: PMC6065343 DOI: 10.1002/ece3.4219] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/08/2018] [Revised: 04/23/2018] [Accepted: 04/24/2018] [Indexed: 01/01/2023] Open
Abstract
Restriction-site associated DNA sequencing (RAD-seq) can identify and score thousands of genetic markers from a group of samples for population-genetics studies. One challenge of de novo RAD-seq analysis is to distinguish paralogous sequence variants (PSVs) from true single-nucleotide polymorphisms (SNPs) associated with orthologous loci. In the absence of a reference genome, it is difficult to differentiate true SNPs from PSVs, and their impact on downstream analysis remains unclear. Here, we introduce a network-based approach, PMERGE that connects fragments based on their DNA sequence similarity to identify probable PSVs. Applying our method to de novo RAD-seq data from 150 Atlantic salmon (Salmo salar) samples collected from 15 locations across the Southern Newfoundland coast allowed the identification of 87% of total PSVs identified through alignment to the Atlantic salmon genome. Removal of these paralogs altered the inferred population structure, highlighting the potential impact of filtering in RAD-seq analysis. PMERGE is also applied to a green crab (Carcinus maenas) data set consisting of 242 samples from 11 different locations and was successfully able to identify and remove the majority of paralogous loci (62%). The PMERGE software can be run as part of the widely used Stacks analysis package.
Collapse
Affiliation(s)
| | - Paul Bentzen
- Marine Gene Probe LaboratoryDepartment of BiologyDalhousie UniversityHalifaxNSCanada
| | - Ian R. Bradbury
- Salmonids SectionScience BranchDepartment of Fisheries and Oceans CanadaSt. John'sNLCanada
| | - Robert G. Beiko
- Faculty of Computer ScienceDalhousie UniversityHalifaxNSCanada
| |
Collapse
|
28
|
Stobie CS, Oosthuizen CJ, Cunningham MJ, Bloomer P. Exploring the phylogeography of a hexaploid freshwater fish by RAD sequencing. Ecol Evol 2018; 8:2326-2342. [PMID: 29468047 PMCID: PMC5817159 DOI: 10.1002/ece3.3821] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/05/2017] [Accepted: 12/17/2017] [Indexed: 12/12/2022] Open
Abstract
The KwaZulu‐Natal yellowfish (Labeobarbus natalensis) is an abundant cyprinid, endemic to KwaZulu‐Natal Province, South Africa. In this study, we developed a single‐nucleotide polymorphism (SNP) dataset from double‐digest restriction site‐associated DNA (ddRAD) sequencing of samples across the distribution. We addressed several hidden challenges, primarily focusing on proper filtering of RAD data and selecting optimal parameters for data processing in polyploid lineages. We used the resulting high‐quality SNP dataset to investigate the population genetic structure of L. natalensis. A small number of mitochondrial markers present in these data had disproportionate influence on the recovered genetic structure. The presence of singleton SNPs also confounded genetic structure. We found a well‐supported division into northern and southern lineages, with further subdivision into five populations, one of which reflects north–south admixture. Approximate Bayesian Computation scenario testing supported a scenario where an ancestral population diverged into northern and southern lineages, which then diverged to yield the current five populations. All river systems showed similar levels of genetic diversity, which appears unrelated to drainage system size. Nucleotide diversity was highest in the smallest river system, the Mbokodweni, which, together with adjacent small coastal systems, should be considered as a key catchment for conservation.
Collapse
Affiliation(s)
- Cora Sabriel Stobie
- Molecular Ecology and Evolution Programme Department of Genetics University of Pretoria Pretoria South Africa
| | - Carel J Oosthuizen
- Molecular Ecology and Evolution Programme Department of Genetics University of Pretoria Pretoria South Africa
| | - Michael J Cunningham
- Molecular Ecology and Evolution Programme Department of Genetics University of Pretoria Pretoria South Africa
| | - Paulette Bloomer
- Molecular Ecology and Evolution Programme Department of Genetics University of Pretoria Pretoria South Africa
| |
Collapse
|
29
|
Luikart G, Kardos M, Hand BK, Rajora OP, Aitken SN, Hohenlohe PA. Population Genomics: Advancing Understanding of Nature. POPULATION GENOMICS 2018. [DOI: 10.1007/13836_2018_60] [Citation(s) in RCA: 30] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/07/2023]
|
30
|
Adding loci improves phylogeographic resolution in red mangroves despite increased missing data: comparing microsatellites and RAD-Seq and investigating loci filtering. Sci Rep 2017; 7:17598. [PMID: 29242627 PMCID: PMC5730610 DOI: 10.1038/s41598-017-16810-7] [Citation(s) in RCA: 70] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2017] [Accepted: 11/15/2017] [Indexed: 12/17/2022] Open
Abstract
The widespread adoption of RAD-Seq data in phylogeography means genealogical relationships previously evaluated using relatively few genetic markers can now be addressed with thousands of loci. One challenge, however, is that RAD-Seq generates complete genotypes for only a small subset of loci or individuals. Simulations indicate that loci with missing data can produce biased estimates of key population genetic parameters, although the influence of such biases in empirical studies is not well understood. Here we compare microsatellite data (8 loci) and RAD-Seq data (six datasets ranging from 239 to 25,198 loci) from red mangroves (Rhizophora mangle) in Florida to evaluate how different levels of data filtering influence phylogeographic inferences. For all datasets, we calculated population genetic statistics and evaluated population structure, and for RAD-Seq datasets, we additionally examined population structure using coalescence. We found higher FST using microsatellites, but that RAD-Seq-based estimates approached those based on microsatellites as more loci with more missing data were included. Analyses of RAD-Seq datasets resolved the classic Gulf-Atlantic coastal phylogeographic break, which was not significant in the microsatellite analyses. Applying multiple levels of filtering to RAD-Seq datasets can provide a more complete picture of potential biases in the data and elucidate subtle phylogeographic patterns.
Collapse
|
31
|
|
32
|
Rodriguez-Ezpeleta N, Álvarez P, Irigoien X. Genetic Diversity and Connectivity in Maurolicus muelleri in the Bay of Biscay Inferred from Thousands of SNP Markers. Front Genet 2017; 8:195. [PMID: 29234350 PMCID: PMC5712365 DOI: 10.3389/fgene.2017.00195] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2017] [Accepted: 11/16/2017] [Indexed: 12/30/2022] Open
Abstract
Mesopelagic fish are largely abundant poorly studied fish that are still intact, but which, due to their potentially great added value, will be imminently exploited by humans. Therefore, studies that provide information to anticipate the anthropogenic impact on this important resource are urgently needed. In particular, knowledge about their connectivity, potential adaptation and resilience are needed. This information can be obtained through the analysis of genome-wide markers which are now relatively easily and cost-efficiently discovered thanks to high-throughput sequencing technologies. Here, we have generated thousands of SNP markers in Maurolicus muelleri, based on the restriction-site associated DNA sequencing method, and preformed population connectivity and genetic diversity analyses in a subset of samples collected from the Bay of Biscay. Our study proves the method valid for obtaining genome-wide markers in this species and provides the first insights into the population genomics of M. muelleri. Importantly, the genomic resources developed here are made available for future studies and set the basics for additional endeavors on this issue.
Collapse
|
33
|
Benestan L, Moore JS, Sutherland BJG, Le Luyer J, Maaroufi H, Rougeux C, Normandeau E, Rycroft N, Atema J, Harris LN, Tallman RF, Greenwood SJ, Clark FK, Bernatchez L. Sex matters in massive parallel sequencing: Evidence for biases in genetic parameter estimation and investigation of sex determination systems. Mol Ecol 2017; 26:6767-6783. [DOI: 10.1111/mec.14217] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2016] [Revised: 03/23/2017] [Accepted: 03/29/2017] [Indexed: 12/26/2022]
Affiliation(s)
- Laura Benestan
- Département de Biologie; Université Laval; Québec QC Canada
- Institut de Biologie Intégrative et des Systèmes (IBIS); Université Laval; Québec QC Canada
| | - Jean-Sébastien Moore
- Département de Biologie; Université Laval; Québec QC Canada
- Institut de Biologie Intégrative et des Systèmes (IBIS); Université Laval; Québec QC Canada
| | - Ben J. G. Sutherland
- Département de Biologie; Université Laval; Québec QC Canada
- Institut de Biologie Intégrative et des Systèmes (IBIS); Université Laval; Québec QC Canada
| | - Jérémy Le Luyer
- Département de Biologie; Université Laval; Québec QC Canada
- Institut de Biologie Intégrative et des Systèmes (IBIS); Université Laval; Québec QC Canada
| | - Halim Maaroufi
- Institut de Biologie Intégrative et des Systèmes (IBIS); Université Laval; Québec QC Canada
| | - Clément Rougeux
- Département de Biologie; Université Laval; Québec QC Canada
- Institut de Biologie Intégrative et des Systèmes (IBIS); Université Laval; Québec QC Canada
| | - Eric Normandeau
- Département de Biologie; Université Laval; Québec QC Canada
- Institut de Biologie Intégrative et des Systèmes (IBIS); Université Laval; Québec QC Canada
| | | | - Jelle Atema
- Department of Biology; Boston University; Boston MA USA
| | - Les N. Harris
- Fisheries and Oceans Canada; Freshwater Institute; Winnipeg MB Canada
| | - Ross F. Tallman
- Fisheries and Oceans Canada; Freshwater Institute; Winnipeg MB Canada
| | - Spencer J. Greenwood
- Department of Biomedical Sciences & AVC Lobster Science Centre; Atlantic Veterinary College; University of Prince Edward Island; Charlottetown PE Canada
| | - Fraser K. Clark
- Department of Biomedical Sciences & AVC Lobster Science Centre; Atlantic Veterinary College; University of Prince Edward Island; Charlottetown PE Canada
| | - Louis Bernatchez
- Département de Biologie; Université Laval; Québec QC Canada
- Institut de Biologie Intégrative et des Systèmes (IBIS); Université Laval; Québec QC Canada
| |
Collapse
|
34
|
Xu S, Song N, Zhao L, Cai S, Han Z, Gao T. Genomic evidence for local adaptation in the ovoviviparous marine fish Sebastiscus marmoratus with a background of population homogeneity. Sci Rep 2017; 7:1562. [PMID: 28484228 PMCID: PMC5431535 DOI: 10.1038/s41598-017-01742-z] [Citation(s) in RCA: 27] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2017] [Accepted: 04/07/2017] [Indexed: 11/09/2022] Open
Abstract
Advances in next-generation sequencing techniques have allowed for the generation of genome-wide sequence data, to gain insight into the dynamics influencing genetic structure and the local adaptation of marine fish. Here, using genotyping-by-sequencing (GBS) technique, we identified 31,119 single nucleotide polymorphisms (SNPs) for Sebastiscus marmoratus in 59 individuals from three populations in Chinese coastal waters. Based on all SNPs, there was little evidence of genetic differentiation among populations. However, outlier tests revealed 329 SNPs putatively under divergent selection across populations. Structural and phylogenetic topology analyses based on the outliers showed clear genetic differentiation among populations. Gene Ontology (GO) annotation results revealed that most of these outliers are known or hypothesized to be involved in metabolic process. Together with previous work using mitochondrial cytochrome b sequences, the present results further suggest that the population structure is strongly influenced by locally adaptive pressure. Overall, adaptive evolution in a heterogeneous environment plays an important role in inducing genetic differentiation among local populations. This study increases understanding of the factors (including gene flow and local adaptation) promoting and constraining population genetic differentiation in marine organisms.
Collapse
Affiliation(s)
- Shengyong Xu
- Institute of Evolution & Marine Biodiversity, Ocean University of China, 5th Yushan Road, Qingdao, 266003, P.R. China
| | - Na Song
- Institute of Evolution & Marine Biodiversity, Ocean University of China, 5th Yushan Road, Qingdao, 266003, P.R. China
| | - Linlin Zhao
- The First Institute of Oceanography, State Oceanic Administration, 6th Xianxialing Road, Qingdao, 266061, P.R. China
| | - Shanshan Cai
- Fishery College, Zhejiang Ocean University, 1st Haidanan Road, Zhoushan, 316022, P.R. China
| | - Zhiqiang Han
- Fishery College, Zhejiang Ocean University, 1st Haidanan Road, Zhoushan, 316022, P.R. China.
| | - Tianxiang Gao
- Fishery College, Zhejiang Ocean University, 1st Haidanan Road, Zhoushan, 316022, P.R. China.
| |
Collapse
|
35
|
Affiliation(s)
- Josephine R. Paris
- Biosciences College of Life and Environmental Sciences University of Exeter Exeter UK
| | - Jamie R. Stevens
- Biosciences College of Life and Environmental Sciences University of Exeter Exeter UK
| | - Julian M. Catchen
- Department of Animal Biology University of Illinois at Urbana–Champaign Urbana IL 61801 USA
| |
Collapse
|
36
|
Van Wyngaarden M, Snelgrove PVR, DiBacco C, Hamilton LC, Rodríguez‐Ezpeleta N, Jeffery NW, Stanley RRE, Bradbury IR. Identifying patterns of dispersal, connectivity and selection in the sea scallop, Placopecten magellanicus, using RADseq-derived SNPs. Evol Appl 2017; 10:102-117. [PMID: 28035239 PMCID: PMC5192885 DOI: 10.1111/eva.12432] [Citation(s) in RCA: 48] [Impact Index Per Article: 6.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/10/2016] [Accepted: 08/23/2016] [Indexed: 12/29/2022] Open
Abstract
Understanding patterns of dispersal and connectivity among marine populations can directly inform fisheries conservation and management. Advances in high-throughput sequencing offer new opportunities for estimating marine connectivity. We used restriction-site-associated DNA sequencing to examine dispersal and realized connectivity in the sea scallop Placopecten magellanicus, an economically important marine bivalve. Based on 245 individuals sampled rangewide at 12 locations from Newfoundland to the Mid-Atlantic Bight, we identified and genotyped 7163 single nucleotide polymorphisms; 112 (1.6%) were identified as outliers potentially under directional selection. Bayesian clustering revealed a discontinuity between northern and southern samples, and latitudinal clines in allele frequencies were observed in 42.9% of the outlier loci and in 24.6% of neutral loci. Dispersal estimates derived using these clines and estimates of linkage disequilibrium imply limited dispersal; 373.1 ± 407.0 km (mean ± SD) for outlier loci and 641.0 ± 544.6 km (mean ± SD) for neutral loci. Our analysis suggests restricted dispersal compared to the species range (>2000 km) and that dispersal and effective connectivity differ. These observations support the hypothesis that limited effective dispersal structures scallop populations along eastern North America. These findings can help refine the appropriate scale of management and conservation in this commercially valuable species.
Collapse
Affiliation(s)
| | - Paul V. R. Snelgrove
- Department of BiologyMemorial University of NewfoundlandSt. John'sNLCanada
- Department of Ocean SciencesMemorial University of NewfoundlandSt. John'sNLCanada
| | | | | | | | - Nicholas W. Jeffery
- Fisheries and Oceans CanadaNorthwest Atlantic Fisheries CentreSt. John'sNLCanada
| | - Ryan R. E. Stanley
- Bedford Institute of OceanographyDartmouthNSCanada
- Faculty of Computer ScienceDalhousie UniversityHalifaxNSCanada
| | - Ian R. Bradbury
- Department of Ocean SciencesMemorial University of NewfoundlandSt. John'sNLCanada
- Fisheries and Oceans CanadaNorthwest Atlantic Fisheries CentreSt. John'sNLCanada
| |
Collapse
|
37
|
Tay YC, Chng MWP, Sew WWG, Rheindt FE, Tun KPP, Meier R. Beyond the Coral Triangle: high genetic diversity and near panmixia in Singapore's populations of the broadcast spawning sea star Protoreaster nodosus. ROYAL SOCIETY OPEN SCIENCE 2016; 3:160253. [PMID: 27853600 PMCID: PMC5108950 DOI: 10.1098/rsos.160253] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 04/12/2016] [Accepted: 07/19/2016] [Indexed: 06/06/2023]
Abstract
The Coral Triangle is widely considered the most important centre of marine biodiversity in Asia while areas on its periphery such as the South China Sea, have received much less interest. Here, we demonstrate that a small population of the knobbly sea star Protoreaster nodosus in Singapore has similarly high levels of genetic diversity as comparable Indonesian populations from the Coral Triangle. The high genetic diversity of this population is remarkable because it is maintained despite decades of continued anthropogenic disturbance. We postulate that it is probably due to broadcast spawning which is likely to maintain high levels of population connectivity. To test this, we analysed 6140 genome-wide single nucleotide polymorphism (SNP) loci for Singapore's populations and demonstrate a pattern of near panmixia. We here document a second case of high genetic diversity and low genetic structure for a broadcast spawner in Singapore, which suggests that such species have high resilience against anthropogenic disturbances. The study demonstrates the feasibility and power of using genome-wide SNPs for connectivity studies of marine invertebrates without a sequenced genome.
Collapse
Affiliation(s)
- Y. C. Tay
- Department of Biological Sciences, National University of Singapore, 14 Science Drive 4, Singapore 117543, Republic of Singapore
| | - M. W. P. Chng
- Centre for Liveable Cities, Ministry of National Development, The URA Centre, 45 Maxwell Road no. 07-01, Singapore 069118, Republic of Singapore
| | - W. W. G. Sew
- DHI Water and Environment (S) Pte Ltd, 1 Cleantech Loop, Cleantech One, no. 03-05, Singapore 637141, Republic of Singapore
| | - F. E. Rheindt
- Department of Biological Sciences, National University of Singapore, 14 Science Drive 4, Singapore 117543, Republic of Singapore
| | - K. P. P. Tun
- National Parks Board, 1 Cluny Road, Singapore 259569, Republic of Singapore
| | - R. Meier
- Department of Biological Sciences, National University of Singapore, 14 Science Drive 4, Singapore 117543, Republic of Singapore
- Lee Kong Chian Natural History Museum, Department of Biological Sciences, National University of Singapore, Singapore, Republic of Singapore
| |
Collapse
|
38
|
Berthouly-Salazar C, Mariac C, Couderc M, Pouzadoux J, Floc’h JB, Vigouroux Y. Genotyping-by-Sequencing SNP Identification for Crops without a Reference Genome: Using Transcriptome Based Mapping as an Alternative Strategy. FRONTIERS IN PLANT SCIENCE 2016; 7:777. [PMID: 27379109 PMCID: PMC4908121 DOI: 10.3389/fpls.2016.00777] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/29/2016] [Accepted: 05/19/2016] [Indexed: 05/26/2023]
Abstract
Next-generation sequencing opens the way for genomic studies of diversity even for non-model crops and animals. Genome reduction techniques are becoming progressively more popular as they allow a fraction of the genome to be sequenced for multiple individuals and/or populations. These techniques are an efficient way to explore genome diversity in non-model crops and animals for which no reference genome is available. Genome reduction techniques emerged with the development of specific pipelines such as UNEAK (Universal Network Enabled Analysis Kit) and Stacks. However, even for non-model crops and animals, transcriptomes are easier to obtain, thereby making it possible to directly map reads. We investigate the direct use of transcriptome as an alternative strategy. Our specific objective was to compare SNPs obtained from the UNEAK pipeline as well as SNPs obtained by directly mapping genotyping-by-sequencing reads on a transcriptome. We assessed the feasibility of both SNP datasets, UNEAK and transcriptome mapping, to investigate the diversity of 91 samples of wild pearl millet sampled across its distribution area. Both approaches produced several tens of thousands of single nucleotide variants, but differed in the way the variants were identified, leading to differences in the frequency spectrum associated with marked differences in the assessment of diversity. Difference in the frequency spectrum significantly biased a large set of diversity analyses as well as detection of selection approaches. However, whatever the approach, we found very similar inference of genetic structure, with three major genetic groups from West, Central, and East Africa. For non-model crops, using transcriptome data as a reference is thus a particularly promising way to obtain a more thorough analysis of datasets generated using genome reduction techniques.
Collapse
|