1
|
Grupstra CGB, Gómez-Corrales M, Fifer JE, Aichelman HE, Meyer-Kaiser KS, Prada C, Davies SW. Integrating cryptic diversity into coral evolution, symbiosis and conservation. Nat Ecol Evol 2024; 8:622-636. [PMID: 38351091 DOI: 10.1038/s41559-023-02319-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2023] [Accepted: 12/12/2023] [Indexed: 04/13/2024]
Abstract
Understanding how diversity evolves and is maintained is critical to predicting the future trajectories of ecosystems under climate change; however, our understanding of these processes is limited in marine systems. Corals, which engineer reef ecosystems, are critically threatened by climate change, and global efforts are underway to conserve and restore populations as attempts to mitigate ocean warming continue. Recently, sequencing efforts have uncovered widespread undescribed coral diversity, including 'cryptic lineages'-genetically distinct but morphologically similar coral taxa. Such cryptic lineages have been identified in at least 24 coral genera spanning the anthozoan phylogeny and across ocean basins. These cryptic lineages co-occur in many reef systems, but their distributions often differ among habitats. Research suggests that cryptic lineages are ecologically specialized and several examples demonstrate differences in thermal tolerance, highlighting the critical implications of this diversity for predicting coral responses to future warming. Here, we draw attention to recent discoveries, discuss how cryptic diversity affects the study of coral adaptation and acclimation to future environments, explore how it shapes symbiotic partnerships, and highlight challenges and opportunities for conservation and restoration efforts.
Collapse
Affiliation(s)
| | | | - James E Fifer
- Department of Biology, Boston University, Boston, MA, USA
| | | | | | - Carlos Prada
- Department of Biological Sciences, University of Rhode Island, Kingston, RI, USA
| | - Sarah W Davies
- Department of Biology, Boston University, Boston, MA, USA.
| |
Collapse
|
2
|
Schmidt TL, Thia JA, Hoffmann AA. How Can Genomics Help or Hinder Wildlife Conservation? Annu Rev Anim Biosci 2024; 12:45-68. [PMID: 37788416 DOI: 10.1146/annurev-animal-021022-051810] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/05/2023]
Abstract
Genomic data are becoming increasingly affordable and easy to collect, and new tools for their analysis are appearing rapidly. Conservation biologists are interested in using this information to assist in management and planning but are typically limited financially and by the lack of genomic resources available for non-model taxa. It is therefore important to be aware of the pitfalls as well as the benefits of applying genomic approaches. Here, we highlight recent methods aimed at standardizing population assessments of genetic variation, inbreeding, and forms of genetic load and methods that help identify past and ongoing patterns of genetic interchange between populations, including those subjected to recent disturbance. We emphasize challenges in applying some of these methods and the need for adequate bioinformatic support. We also consider the promises and challenges of applying genomic approaches to understand adaptive changes in natural populations to predict their future adaptive capacity.
Collapse
Affiliation(s)
- Thomas L Schmidt
- School of BioSciences, Bio21 Institute, University of Melbourne, Parkville, Victoria, Australia;
| | - Joshua A Thia
- School of BioSciences, Bio21 Institute, University of Melbourne, Parkville, Victoria, Australia;
| | - Ary A Hoffmann
- School of BioSciences, Bio21 Institute, University of Melbourne, Parkville, Victoria, Australia;
| |
Collapse
|
3
|
Puritz JB, Guo X, Hare M, He Y, Hillier LW, Jin S, Liu M, Lotterhos KE, Minx P, Modak T, Proestou D, Rice ES, Tomlinson C, Warren WC, Witkop E, Zhao H, Gomez-Chiarri M. A second unveiling: Haplotig masking of the eastern oyster genome improves population-level inference. Mol Ecol Resour 2024; 24:e13801. [PMID: 37186213 DOI: 10.1111/1755-0998.13801] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2022] [Revised: 12/16/2022] [Accepted: 03/20/2023] [Indexed: 05/17/2023]
Abstract
Genome assembly can be challenging for species that are characterized by high amounts of polymorphism, heterozygosity, and large effective population sizes. High levels of heterozygosity can result in genome mis-assemblies and a larger than expected genome size due to the haplotig versions of a single locus being assembled as separate loci. Here, we describe the first chromosome-level genome for the eastern oyster, Crassostrea virginica. Publicly released and annotated in 2017, the assembly has a scaffold N50 of 54 mb and is over 97.3% complete based on BUSCO analysis. The genome assembly for the eastern oyster is a critical resource for foundational research into molluscan adaptation to a changing environment and for selective breeding for the aquaculture industry. Subsequent resequencing data suggested the presence of haplotigs in the original assembly, and we developed a post hoc method to break up chimeric contigs and mask haplotigs in published heterozygous genomes and evaluated improvements to the accuracy of downstream analysis. Masking haplotigs had a large impact on SNP discovery and estimates of nucleotide diversity and had more subtle and nuanced effects on estimates of heterozygosity, population structure analysis, and outlier detection. We show that haplotig masking can be a powerful tool for improving genomic inference, and we present an open, reproducible resource for the masking of haplotigs in any published genome.
Collapse
Affiliation(s)
- Jonathan B Puritz
- Department of Biological Sciences, University of Rhode Island, Kingston, Rhode Island, USA
| | - Ximing Guo
- Haskin Shellfish Research Laboratory, Department of Marine and Coastal Sciences, Rutgers University, Port Norris, New Jersey, USA
| | - Matthew Hare
- Department of Natural Resources and the Environment, Cornell University, Ithaca, New York, USA
| | - Yan He
- Haskin Shellfish Research Laboratory, Department of Marine and Coastal Sciences, Rutgers University, Port Norris, New Jersey, USA
| | - LaDeana W Hillier
- Department of Genome Sciences, University of Washington, Seattle, Washington, USA
| | - Shubo Jin
- Haskin Shellfish Research Laboratory, Department of Marine and Coastal Sciences, Rutgers University, Port Norris, New Jersey, USA
| | - Ming Liu
- Haskin Shellfish Research Laboratory, Department of Marine and Coastal Sciences, Rutgers University, Port Norris, New Jersey, USA
| | - Katie E Lotterhos
- Department of Marine and Environmental Sciences, Northeastern University Marine Science Center, Nahant, Massachusetts, USA
| | - Pat Minx
- Donald Danforth Plant Science Center, Olivette, Missouri, USA
| | - Tejashree Modak
- Department of Cell and Molecular Biology, University of Rhode Island, Kingston, Rhode Island, USA
| | - Dina Proestou
- USDA Agricultural Research Service, National Cold Water Marine Aquaculture Center, Kingston, Rhode Island, USA
| | - Edward S Rice
- Bond Life Sciences Center, University of Missouri, Columbia, Missouri, USA
| | - Chad Tomlinson
- McDonnell Genome Institute, Washington University School of Medicine, St Louis, Missouri, USA
| | - Wesley C Warren
- Departments of Animal Sciences and Surgery, Institute of Informatics and Data Sciences, Bond Life Sciences Center, University of Missouri, Columbia, Missouri, USA
| | - Erin Witkop
- Department of Fisheries, Animal and Veterinary Sciences, University of Rhode Island, Kingston, Rhode Island, USA
| | - Honggang Zhao
- Department of Natural Resources and the Environment, Cornell University, Ithaca, New York, USA
| | - Marta Gomez-Chiarri
- Department of Fisheries, Animal and Veterinary Sciences, University of Rhode Island, Kingston, Rhode Island, USA
| |
Collapse
|
4
|
Coimbra RTF, Winter S, Muneza A, Fennessy S, Otiende M, Mijele D, Masiaine S, Stacy-Dawes J, Fennessy J, Janke A. Genomic analysis reveals limited hybridization among three giraffe species in Kenya. BMC Biol 2023; 21:215. [PMID: 37833744 PMCID: PMC10576358 DOI: 10.1186/s12915-023-01722-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/27/2023] [Accepted: 10/02/2023] [Indexed: 10/15/2023] Open
Abstract
BACKGROUND In the speciation continuum, the strength of reproductive isolation varies, and species boundaries are blurred by gene flow. Interbreeding among giraffe (Giraffa spp.) in captivity is known, and anecdotal reports of natural hybrids exist. In Kenya, Nubian (G. camelopardalis camelopardalis), reticulated (G. reticulata), and Masai giraffe sensu stricto (G. tippelskirchi tippelskirchi) are parapatric, and thus, the country might be a melting pot for these taxa. We analyzed 128 genomes of wild giraffe, 113 newly sequenced, representing these three taxa. RESULTS We found varying levels of Nubian ancestry in 13 reticulated giraffe sampled across the Laikipia Plateau most likely reflecting historical gene flow between these two lineages. Although comparatively weaker signs of ancestral gene flow and potential mitochondrial introgression from reticulated into Masai giraffe were also detected, estimated admixture levels between these two lineages are minimal. Importantly, contemporary gene flow between East African giraffe lineages was not statistically significant. Effective population sizes have declined since the Late Pleistocene, more severely for Nubian and reticulated giraffe. CONCLUSIONS Despite historically hybridizing, these three giraffe lineages have maintained their overall genomic integrity suggesting effective reproductive isolation, consistent with the previous classification of giraffe into four species.
Collapse
Affiliation(s)
- Raphael T F Coimbra
- Senckenberg Biodiversity and Climate Research Centre, Frankfurt, Germany.
- Institute for Ecology, Evolution and Diversity, Goethe University, Frankfurt, Germany.
| | - Sven Winter
- Senckenberg Biodiversity and Climate Research Centre, Frankfurt, Germany
- Research Institute of Wildlife Ecology, University of Veterinary Medicine Vienna, Vienna, Austria
| | | | | | | | | | | | | | - Julian Fennessy
- Giraffe Conservation Foundation, Windhoek, Namibia
- School of Biology and Environmental Science, University College Dublin, Dublin, Ireland
| | - Axel Janke
- Senckenberg Biodiversity and Climate Research Centre, Frankfurt, Germany.
- Institute for Ecology, Evolution and Diversity, Goethe University, Frankfurt, Germany.
- LOEWE Centre for Translational Biodiversity Genomics, Frankfurt, Germany.
| |
Collapse
|
5
|
Kersten O, Star B, Krabberød AK, Atmore LM, Tørresen OK, Anker-Nilssen T, Descamps S, Strøm H, Johansson US, Sweet PR, Jakobsen KS, Boessenkool S. Hybridization of Atlantic puffins in the Arctic coincides with 20th-century climate change. SCIENCE ADVANCES 2023; 9:eadh1407. [PMID: 37801495 PMCID: PMC10558128 DOI: 10.1126/sciadv.adh1407] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/17/2023] [Accepted: 09/06/2023] [Indexed: 10/08/2023]
Abstract
The Arctic is experiencing the fastest rates of global warming, leading to shifts in the distribution of its biota and increasing the potential for hybridization. However, genomic evidence of recent hybridization events in the Arctic remains unexpectedly rare. Here, we use whole-genome sequencing of contemporary and 122-year-old historical specimens to investigate the origin of an Arctic hybrid population of Atlantic puffins (Fratercula arctica) on Bjørnøya, Norway. We show that the hybridization between the High Arctic, large-bodied subspecies F. a. naumanni and the temperate, smaller-sized subspecies F. a. arctica began as recently as six generations ago due to an unexpected southward range expansion of F. a. naumanni. Moreover, we find a significant temporal loss of genetic diversity across Arctic and temperate puffin populations. Our observations provide compelling genomic evidence of the impacts of recent distributional shifts and loss of diversity in Arctic communities during the 20th century.
Collapse
Affiliation(s)
- Oliver Kersten
- Centre for Ecological and Evolutionary Synthesis (CEES), Department of Biosciences, University of Oslo, Oslo, Norway
| | - Bastiaan Star
- Centre for Ecological and Evolutionary Synthesis (CEES), Department of Biosciences, University of Oslo, Oslo, Norway
| | - Anders K. Krabberød
- Section for Genetics and Evolutionary Biology (Evogene), Department of Biosciences, University of Oslo, Oslo, Norway
| | - Lane M. Atmore
- Centre for Ecological and Evolutionary Synthesis (CEES), Department of Biosciences, University of Oslo, Oslo, Norway
| | - Ole K. Tørresen
- Centre for Ecological and Evolutionary Synthesis (CEES), Department of Biosciences, University of Oslo, Oslo, Norway
| | | | | | - Hallvard Strøm
- Norwegian Polar Institute, Fram Centre, Langnes, Tromsø, Norway
| | | | - Paul R. Sweet
- American Museum of Natural History, New York, NY, USA
| | - Kjetill S. Jakobsen
- Centre for Ecological and Evolutionary Synthesis (CEES), Department of Biosciences, University of Oslo, Oslo, Norway
| | - Sanne Boessenkool
- Centre for Ecological and Evolutionary Synthesis (CEES), Department of Biosciences, University of Oslo, Oslo, Norway
| |
Collapse
|
6
|
Smith SK, Frazel PW, Khodadadi-Jamayran A, Zappile P, Marier C, Okhovat M, Brown S, Long MA, Heguy A, Phelps SM. De novo assembly and annotation of the singing mouse genome. BMC Genomics 2023; 24:569. [PMID: 37749493 PMCID: PMC10521431 DOI: 10.1186/s12864-023-09678-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/24/2023] [Accepted: 09/14/2023] [Indexed: 09/27/2023] Open
Abstract
BACKGROUND Developing genomic resources for a diverse range of species is an important step towards understanding the mechanisms underlying complex traits. Specifically, organisms that exhibit unique and accessible phenotypes-of-interest allow researchers to address questions that may be ill-suited to traditional model organisms. We sequenced the genome and transcriptome of Alston's singing mouse (Scotinomys teguina), an emerging model for social cognition and vocal communication. In addition to producing advertisement songs used for mate attraction and male-male competition, these rodents are diurnal, live at high-altitudes, and are obligate insectivores, providing opportunities to explore diverse physiological, ecological, and evolutionary questions. RESULTS Using PromethION, Illumina, and PacBio sequencing, we produced an annotated genome and transcriptome, which were validated using gene expression and functional enrichment analyses. To assess the usefulness of our assemblies, we performed single nuclei sequencing on cells of the orofacial motor cortex, a brain region implicated in song coordination, identifying 12 cell types. CONCLUSIONS These resources will provide the opportunity to identify the molecular basis of complex traits in singing mice as well as to contribute data that can be used for large-scale comparative analyses.
Collapse
Affiliation(s)
- Samantha K Smith
- Department of Integrative Biology, University of Texas at Austin, Austin, TX, 78712, USA.
| | - Paul W Frazel
- Department of Neuroscience and Physiology, New York University Grossman School of Medicine, New York, NY, 10016, USA
| | - Alireza Khodadadi-Jamayran
- Applied Bioinformatics Laboratory, New York University Grossman School of Medicine, New York, NY, 10016, USA
| | - Paul Zappile
- Genome Technology Center, New York University Grossman School of Medicine, New York, NY, 10016, USA
| | - Christian Marier
- Genome Technology Center, New York University Grossman School of Medicine, New York, NY, 10016, USA
| | - Mariam Okhovat
- Department of Integrative Biology, University of Texas at Austin, Austin, TX, 78712, USA
- Present Address: Oregon Health & Science University, Portland, OR, USA
| | - Stuart Brown
- NYU Center for Health Informatics and Bioinformatics, New York University Grossman School of Medicine, New York, NY, 10016, USA
- Present Address: Exxon Mobil Corporate, Houston, TX, USA
| | - Michael A Long
- Department of Neuroscience and Physiology, New York University Grossman School of Medicine, New York, NY, 10016, USA
| | - Adriana Heguy
- Genome Technology Center, New York University Grossman School of Medicine, New York, NY, 10016, USA
| | - Steven M Phelps
- Department of Integrative Biology, University of Texas at Austin, Austin, TX, 78712, USA
| |
Collapse
|
7
|
Smith SK, Aglyamova G, Oberg EW, Fuiman LA, Matz MV. Gene expression underlying variation in the survival skills of red drum larvae (Sciaenops ocellatus). JOURNAL OF FISH BIOLOGY 2023; 103:704-714. [PMID: 37300518 DOI: 10.1111/jfb.15480] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/22/2022] [Accepted: 06/09/2023] [Indexed: 06/12/2023]
Abstract
Mortality rates of marine fish larvae are incredibly high and can determine year-class strength. The major causes of larval mortality are predation and starvation, and the performance of larvae in survival skills that can mitigate this mortality (predator evasion, foraging) varies among individuals and cohorts, but the causes of the variation are not known. Transcriptomics can link gene expression variation to phenotypic variation at the whole-system level to investigate the molecular basis of behavioural variation. We used tag-based RNA-sequencing to examine the molecular basis of variation in predator evasion and routine swimming (trait related to foraging efficiency) in the larval red drum, Sciaenops ocellatus. We looked for functional gene networks in which interindividual variation would explain variation in larval behavioural performance. We identified co-expressed gene groups ("modules") associated with predator evasion traits and found enrichment of motor, neural and energy metabolism pathways. These functional associations and pattern of correlations between modules and traits suggest that energy availability and allocation were responsible for the magnitude of startle responses, while differential neural and motor activation were associated with differences in response latency.
Collapse
Affiliation(s)
- Samantha K Smith
- Department of Integrative Biology, University of Texas at Austin, Austin, Texas, USA
| | - Galina Aglyamova
- Department of Integrative Biology, University of Texas at Austin, Austin, Texas, USA
| | - Erik W Oberg
- Marine Science Institute, University of Texas at Austin, Port Aransas, Texas, USA
| | - Lee A Fuiman
- Marine Science Institute, University of Texas at Austin, Port Aransas, Texas, USA
| | - Mikhail V Matz
- Department of Integrative Biology, University of Texas at Austin, Austin, Texas, USA
| |
Collapse
|
8
|
Leiva C, Pérez-Portela R, Lemer S. Genomic signatures suggesting adaptation to ocean acidification in a coral holobiont from volcanic CO 2 seeps. Commun Biol 2023; 6:769. [PMID: 37481685 PMCID: PMC10363134 DOI: 10.1038/s42003-023-05103-7] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2023] [Accepted: 07/06/2023] [Indexed: 07/24/2023] Open
Abstract
Ocean acidification, caused by anthropogenic CO2 emissions, is predicted to have major consequences for reef-building corals, jeopardizing the scaffolding of the most biodiverse marine habitats. However, whether corals can adapt to ocean acidification and how remains unclear. We addressed these questions by re-examining transcriptome and genome data of Acropora millepora coral holobionts from volcanic CO2 seeps with end-of-century pH levels. We show that adaptation to ocean acidification is a wholistic process involving the three main compartments of the coral holobiont. We identified 441 coral host candidate adaptive genes involved in calcification, response to acidification, and symbiosis; population genetic differentiation in dinoflagellate photosymbionts; and consistent transcriptional microbiome activity despite microbial community shifts. Coral holobionts from natural analogues to future ocean conditions harbor beneficial genetic variants with far-reaching rapid adaptation potential. In the face of climate change, these populations require immediate conservation strategies as they could become key to coral reef survival.
Collapse
Affiliation(s)
- Carlos Leiva
- University of Guam Marine Laboratory, 303 University Drive, 96923, Mangilao, Guam, USA.
| | - Rocío Pérez-Portela
- Departament de Biologia Evolutiva, Ecologia i Ciències Ambientals, Facultat de Biologia, Universitat de Barcelona, Av. Diagonal 643, 08028, Barcelona, Spain
- Institut de Recerca de la Biodiversitat (IRBio), Universitat de Barcelona, Barcelona, Spain
| | - Sarah Lemer
- University of Guam Marine Laboratory, 303 University Drive, 96923, Mangilao, Guam, USA
| |
Collapse
|
9
|
Chambers EA, Tarvin RD, Santos JC, Ron SR, Betancourth-Cundar M, Hillis DM, Matz MV, Cannatella DC. 2b or not 2b? 2bRAD is an effective alternative to ddRAD for phylogenomics. Ecol Evol 2023; 13:e9842. [PMID: 36911313 PMCID: PMC9994478 DOI: 10.1002/ece3.9842] [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: 01/24/2023] [Revised: 02/02/2023] [Accepted: 02/03/2023] [Indexed: 03/10/2023] Open
Abstract
Restriction-site-associated DNA sequencing (RADseq) has become an accessible way to obtain genome-wide data in the form of single-nucleotide polymorphisms (SNPs) for phylogenetic inference. Nonetheless, how differences in RADseq methods influence phylogenetic estimation is poorly understood because most comparisons have largely relied on conceptual predictions rather than empirical tests. We examine how differences in ddRAD and 2bRAD data influence phylogenetic estimation in two non-model frog groups. We compare the impact of method choice on phylogenetic information, missing data, and allelic dropout, considering different sequencing depths. Given that researchers must balance input (funding, time) with output (amount and quality of data), we also provide comparisons of laboratory effort, computational time, monetary costs, and the repeatability of library preparation and sequencing. Both 2bRAD and ddRAD methods estimated well-supported trees, even at low sequencing depths, and had comparable amounts of missing data, patterns of allelic dropout, and phylogenetic signal. Compared to ddRAD, 2bRAD produced more repeatable datasets, had simpler laboratory protocols, and had an overall faster bioinformatics assembly. However, many fewer parsimony-informative sites per SNP were obtained from 2bRAD data when using native pipelines, highlighting a need for further investigation into the effects of each pipeline on resulting datasets. Our study underscores the importance of comparing RADseq methods, such as expected results and theoretical performance using empirical datasets, before undertaking costly experiments.
Collapse
Affiliation(s)
- E Anne Chambers
- Department of Integrative Biology and Biodiversity Center University of Texas at Austin Austin Texas USA.,Department of Environmental Science, Policy, and Management and Museum of Vertebrate Zoology University of California Berkeley Berkeley California USA
| | - Rebecca D Tarvin
- Department of Integrative Biology and Biodiversity Center University of Texas at Austin Austin Texas USA.,Department of Integrative Biology and Museum of Vertebrate Zoology University of California Berkeley Berkeley California USA
| | - Juan C Santos
- Department of Biological Sciences St John's University New York New York USA
| | - Santiago R Ron
- Museo de Zoología, Escuela de Ciencias Biológicas Pontificia Universidad Católica del Ecuador Quito Ecuador
| | | | - David M Hillis
- Department of Integrative Biology and Biodiversity Center University of Texas at Austin Austin Texas USA
| | - Mikhail V Matz
- Department of Integrative Biology and Biodiversity Center University of Texas at Austin Austin Texas USA
| | - David C Cannatella
- Department of Integrative Biology and Biodiversity Center University of Texas at Austin Austin Texas USA
| |
Collapse
|
10
|
Matias AMA, Popovic I, Thia JA, Cooke IR, Torda G, Lukoschek V, Bay LK, Kim SW, Riginos C. Cryptic diversity and spatial genetic variation in the coral Acropora tenuis and its endosymbionts across the Great Barrier Reef. Evol Appl 2023; 16:293-310. [PMID: 36793689 PMCID: PMC9923489 DOI: 10.1111/eva.13435] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/06/2022] [Revised: 05/20/2022] [Accepted: 05/29/2022] [Indexed: 11/26/2022] Open
Abstract
Genomic studies are uncovering extensive cryptic diversity within reef-building corals, suggesting that evolutionarily and ecologically relevant diversity is highly underestimated in the very organisms that structure coral reefs. Furthermore, endosymbiotic algae within coral host species can confer adaptive responses to environmental stress and may represent additional axes of coral genetic variation that are not constrained by taxonomic divergence of the cnidarian host. Here, we examine genetic variation in a common and widespread, reef-building coral, Acropora tenuis, and its associated endosymbiotic algae along the entire expanse of the Great Barrier Reef (GBR). We use SNPs derived from genome-wide sequencing to characterize the cnidarian coral host and organelles from zooxanthellate endosymbionts (genus Cladocopium). We discover three distinct and sympatric genetic clusters of coral hosts, whose distributions appear associated with latitude and inshore-offshore reef position. Demographic modelling suggests that the divergence history of the three distinct host taxa ranges from 0.5 to 1.5 million years ago, preceding the GBR's formation, and has been characterized by low-to-moderate ongoing inter-taxon gene flow, consistent with occasional hybridization and introgression typifying coral evolution. Despite this differentiation in the cnidarian host, A. tenuis taxa share a common symbiont pool, dominated by the genus Cladocopium (Clade C). Cladocopium plastid diversity is not strongly associated with host identity but varies with reef location relative to shore: inshore colonies contain lower symbiont diversity on average but have greater differences between colonies as compared with symbiont communities from offshore colonies. Spatial genetic patterns of symbiont communities could reflect local selective pressures maintaining coral holobiont differentiation across an inshore-offshore environmental gradient. The strong influence of environment (but not host identity) on symbiont community composition supports the notion that symbiont community composition responds to habitat and may assist in the adaptation of corals to future environmental change.
Collapse
Affiliation(s)
- Ambrocio Melvin A Matias
- Institute of Biology University of the Philippines Diliman Quezon City Philippines.,School of Biological Sciences The University of Queensland St. Lucia Queensland Australia
| | - Iva Popovic
- School of Biological Sciences The University of Queensland St. Lucia Queensland Australia
| | - Joshua A Thia
- Bio21 Institute, School of BioSciences The University of Melbourne Parkeville Victoria Australia
| | - Ira R Cooke
- College of Public Health, Medical and Veterinary Sciences James Cook University Townsville Queensland Australia
| | - Gergely Torda
- ARC Centre of Excellence for Coral Reef Studies James Cook University Townsville Queensland Australia
| | - Vimoksalehi Lukoschek
- ARC Centre of Excellence for Coral Reef Studies James Cook University Townsville Queensland Australia.,Gold Coast University Hospital QLD Health Southport Queensland Australia
| | - Line K Bay
- Australian Institute of Marine Science Townsville Queensland Australia
| | - Sun W Kim
- School of Biological Sciences The University of Queensland St. Lucia Queensland Australia
| | - Cynthia Riginos
- School of Biological Sciences The University of Queensland St. Lucia Queensland Australia
| |
Collapse
|
11
|
Martín-Blázquez R, Bakkali M, Ruiz-Estévez M, Garrido-Ramos MA. Comparison between the Gametophyte and the Sporophyte Transcriptomes of the Endangered Fern Vandenboschia speciosa. Genes (Basel) 2023; 14:genes14010166. [PMID: 36672907 PMCID: PMC9859580 DOI: 10.3390/genes14010166] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2022] [Revised: 01/04/2023] [Accepted: 01/05/2023] [Indexed: 01/10/2023] Open
Abstract
Genomic resources are essential to understanding the evolution and functional biology of organisms. Nevertheless, generating genomic resources from endangered species may be challenging due to the scarcity of available specimens and sampling difficulties. In this study, we compare the transcriptomes of the sporophyte and the gametophyte of the endangered fern Vandenboschia speciosa. After Illumina sequencing and de novo transcriptome assembly of the gametophyte, annotation proved the existence of cross-species contamination in the gametophyte sample. Thus, we developed an in silico decontamination step for the gametophyte sequences. Once the quality check of the decontaminated reads passed, we produced a de novo assembly with the decontaminated gametophyte reads (with 43,139 contigs) and another combining the sporophyte and in silico decontaminated gametophyte reads (with 42,918 contigs). A comparison of the enriched GO terms from the top 1000 most expressed transcripts from both tissues showed that the gametophyte GO term set was enriched in sequences involved in development, response to stress, and plastid organization, while the sporophyte GO term set had a larger representation of more general metabolic functions. This study complements the available genomic resources on the life cycle of the endangered fern Vandenboschia speciosa.
Collapse
Affiliation(s)
- Rubén Martín-Blázquez
- Department of Evolutionary Ecology, Estación Biológica de Doñana, Consejo Superior de Investigaciones Científicas (CSIC), Isla de la Cartuja, 41092 Sevilla, Spain
| | - Mohammed Bakkali
- Departamento de Genética, Universidad de Granada, 18071 Granada, Spain
| | | | | |
Collapse
|
12
|
Alvarado-Cerón V, Muñiz-Castillo AI, León-Pech MG, Prada C, Arias-González JE. A decade of population genetics studies of scleractinian corals: A systematic review. MARINE ENVIRONMENTAL RESEARCH 2023; 183:105781. [PMID: 36371949 DOI: 10.1016/j.marenvres.2022.105781] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/09/2022] [Revised: 10/11/2022] [Accepted: 10/12/2022] [Indexed: 06/16/2023]
Abstract
Coral reefs are the most diverse marine ecosystems. However, coral cover has decreased worldwide due to natural disturbances, climate change, and local anthropogenic drivers. In recent decades, various genetic methods and molecular markers have been developed to assess genetic diversity, structure, and connectivity in different coral species to determine the vulnerability of their populations. This review aims to identify population genetic studies of scleractinian corals in the last decade (2010-2020), and the techniques and molecular markers used. Bibliometric analysis was conducted to identify journals and authors working in this field. We then calculated the number of genetic studies by species and ecoregion based on data obtained from 178 studies found in Scopus and Web of Science. Coral Reefs and Molecular Ecology were the main journals published population genetics studies, and microsatellites are the most widely used molecular markers. The Caribbean, Australian Barrier Reef, and South Kuroshio in Japan are among the ecoregions with the most population genetics data. In contrast, we found limited information about the Coral Triangle, a region with the highest biodiversity and key to coral reef conservation. Notably, only 117 (out of 1500 described) scleractinian coral species have genetic studies. This review emphasizes which coral species have been studied and highlights remaining gaps and locations where such data is critical for coral conservation.
Collapse
Affiliation(s)
- Viridiana Alvarado-Cerón
- Departamento de Recursos del Mar, Centro de Investigación y de Estudios Avanzados del I.P.N., Unidad Mérida. Km. 6 Antigua carretera a Progreso, Cordemex, 97310, Mérida, Yucatán, Mexico.
| | - Aarón Israel Muñiz-Castillo
- Departamento de Recursos del Mar, Centro de Investigación y de Estudios Avanzados del I.P.N., Unidad Mérida. Km. 6 Antigua carretera a Progreso, Cordemex, 97310, Mérida, Yucatán, Mexico.
| | - María Geovana León-Pech
- Department of Biological Science, University of Rhode Island, 120 Flag Road, Kingston, RI, 02881, USA.
| | - Carlos Prada
- Department of Biological Science, University of Rhode Island, 120 Flag Road, Kingston, RI, 02881, USA.
| | - Jesús Ernesto Arias-González
- Departamento de Recursos del Mar, Centro de Investigación y de Estudios Avanzados del I.P.N., Unidad Mérida. Km. 6 Antigua carretera a Progreso, Cordemex, 97310, Mérida, Yucatán, Mexico.
| |
Collapse
|
13
|
Hoey JA, Able KW, Pinsky ML. Genetic decline and recovery of a demographically rebuilt fishery species. Mol Ecol 2022; 31:5684-5698. [PMID: 36114805 PMCID: PMC9828022 DOI: 10.1111/mec.16697] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/12/2020] [Revised: 09/09/2022] [Accepted: 09/15/2022] [Indexed: 01/13/2023]
Abstract
The demographic history of a population is important for conservation and evolution, but this history is unknown for many populations. Methods that use genomic data have been developed to infer demography, but they can be challenging to implement and interpret, particularly for large populations. Thus, understanding if and when genetic estimates of demography correspond to true population history is important for assessing the performance of these genetic methods. Here, we used double-digest restriction-site associated DNA (ddRAD) sequencing data from archived collections of larval summer flounder (Paralichthys dentatus, n = 279) from three cohorts (1994-1995, 1997-1998 and 2008-2009) along the U.S. East coast to examine how contemporary effective population size and genetic diversity responded to changes in abundance in a natural population. Despite little to no detectable change in genetic diversity, coalescent-based demographic modelling from site frequency spectra revealed that summer flounder effective population size declined dramatically in the early 1980s. The timing and direction of change corresponded well with the observed decline in spawning stock census abundance in the late 1980s from independent fish surveys. Census abundance subsequently recovered and achieved the prebottleneck size. Effective population size also grew following the bottleneck. Our results for summer flounder demonstrate that genetic sampling and site frequency spectra can be useful for detecting population dynamics, even in species with large effective sizes.
Collapse
Affiliation(s)
- Jennifer A. Hoey
- Ecology, Evolution, & Natural ResourcesRutgers UniversityNew BrunswickNew JerseyUSA,Institute for Biodiversity Science and SustainabilityCalifornia Academy of SciencesSan FranciscoCaliforniaUSA
| | - Kenneth W. Able
- Marine Field Station, Department of Marine and Coastal Sciences, Rutgers UniversityTuckertonNew JerseyUSA
| | - Malin L. Pinsky
- Ecology, Evolution, & Natural ResourcesRutgers UniversityNew BrunswickNew JerseyUSA
| |
Collapse
|
14
|
Kenkel CD, Smith J, Hubbard KA, Chadwick C, Lorenzen N, Tatters AO, Caron DA. Reduced representation sequencing accurately quantifies relative abundance and reveals population-level variation in Pseudo-nitzschia spp. HARMFUL ALGAE 2022; 118:102314. [PMID: 36195429 PMCID: PMC9869635 DOI: 10.1016/j.hal.2022.102314] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/09/2022] [Revised: 08/16/2022] [Accepted: 08/22/2022] [Indexed: 06/16/2023]
Abstract
Certain species within the genus Pseudo-nitzschia are able to produce the neurotoxin domoic acid (DA), which can cause illness in humans, mass-mortality of marine animals, and closure of commercial and recreational shellfisheries during toxic events. Understanding and forecasting blooms of these harmful species is a primary management goal. However, accurately predicting the onset and severity of bloom events remains difficult, in part because the underlying drivers of bloom formation have not been fully resolved. Furthermore, Pseudo-nitzschia species often co-occur, and recent work suggests that the genetic composition of a Pseudo-nitzschia bloom may be a better predictor of toxicity than prevailing environmental conditions. We developed a novel next-generation sequencing assay using restriction site-associated DNA (2b-RAD) genotyping and applied it to mock Pseudo-nitzschia communities generated by mixing cultures of different species in known abundances. On average, 94% of the variance in observed species abundance was explained by the expected abundance. In addition, the false positive rate was low (0.45% on average) and unrelated to read depth, and false negatives were never observed. Application of this method to environmental DNA samples collected during natural Pseudo-nitzschia spp. bloom events in Southern California revealed that increases in DA were associated with increases in the relative abundance of P. australis. Although the absolute correlation across time-points was weak, an independent species fingerprinting assay (Automated Ribosomal Intergenic Spacer Analysis) supported this and identified other potentially toxic species. Finally, we assessed population-level genomic variation by mining SNPs from the environmental 2bRAD dataset. Consistent shifts in allele frequencies in P. pungens and P. subpacifica were detected between high and low DA years, suggesting that different intraspecific variants may be associated with prevailing environmental conditions or the presence of DA. Taken together, this method presents a potentially cost-effective and high-throughput approach for studies aiming to evaluate both population and species dynamics in mixed samples.
Collapse
Affiliation(s)
- Carly D Kenkel
- Department of Biological Sciences, University of Southern California, 3616 Trousdale Parkway, Los Angeles, CA 90089, USA.
| | - Jayme Smith
- Southern California Coastal Water Research Project, 3535 Harbor Boulevard, Suite 110, Costa Mesa, CA, 92626, USA
| | - Katherine A Hubbard
- Florida Fish and Wildlife Conservation Commission-Fish and Wildlife Research Institute (FWC-FWRI), 100 8th Ave. SE, St. Petersburg, FL 33701, USA; Biology Department, Woods Hole Oceanographic Institution, 266 Woods Hole Road, Woods Hole, MA 02543, USA
| | - Christina Chadwick
- Florida Fish and Wildlife Conservation Commission-Fish and Wildlife Research Institute (FWC-FWRI), 100 8th Ave. SE, St. Petersburg, FL 33701, USA
| | - Nico Lorenzen
- Department of Biological Sciences, University of Southern California, 3616 Trousdale Parkway, Los Angeles, CA 90089, USA
| | - Avery O Tatters
- U.S. Environmental Protection Agency, Gulf Ecosystem Measurement and Modeling Division, 1 Sabine Island Drive, Gulf Breeze, FL, 32561, USA
| | - David A Caron
- Department of Biological Sciences, University of Southern California, 3616 Trousdale Parkway, Los Angeles, CA 90089, USA
| |
Collapse
|
15
|
Lutier M, Di Poi C, Gazeau F, Appolis A, Le Luyer J, Pernet F. Revisiting tolerance to ocean acidification: Insights from a new framework combining physiological and molecular tipping points of Pacific oyster. GLOBAL CHANGE BIOLOGY 2022; 28:3333-3348. [PMID: 35092108 DOI: 10.1111/gcb.16101] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/14/2021] [Revised: 12/02/2021] [Accepted: 01/06/2022] [Indexed: 06/14/2023]
Abstract
Studies on the impact of ocean acidification on marine organisms involve exposing organisms to future acidification scenarios, which has limited relevance for coastal calcifiers living in a mosaic of habitats. Identification of tipping points beyond which detrimental effects are observed is a widely generalizable proxy of acidification susceptibility at the population level. This approach is limited to a handful of studies that focus on only a few macro-physiological traits, thus overlooking the whole organism response. Here we develop a framework to analyze the broad macro-physiological and molecular responses over a wide pH range in juvenile oyster. We identify low tipping points for physiological traits at pH 7.3-6.9 that coincide with a major reshuffling in membrane lipids and transcriptome. In contrast, a drop in pH affects shell parameters above tipping points, likely impacting animal fitness. These findings were made possible by the development of an innovative methodology to synthesize and identify the main patterns of variations in large -omic data sets, fitting them to pH and identifying molecular tipping points. We propose the broad application of our framework to the assessment of effects of global change on other organisms.
Collapse
Affiliation(s)
| | - Carole Di Poi
- Univ Brest, Ifremer, CNRS, IRD, LEMAR, Plouzané, France
| | - Frédéric Gazeau
- Laboratoire d'Océanographie de Villefranche, LOV Sorbonne Université, CNRS, Villefranche-sur-Mer, France
| | | | - Jérémy Le Luyer
- EIO UPF/IRD/ILM/Ifremer, Labex CORAIL, Unité RMPF, Centre Océanologique du Pacifique, Vairao, French Polynesia
| | | |
Collapse
|
16
|
Conservation Genomics of Two Threatened Subspecies of Northern Giraffe: The West African and the Kordofan Giraffe. Genes (Basel) 2022; 13:genes13020221. [PMID: 35205265 PMCID: PMC8872558 DOI: 10.3390/genes13020221] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/30/2021] [Revised: 01/22/2022] [Accepted: 01/23/2022] [Indexed: 11/17/2022] Open
Abstract
Three of the four species of giraffe are threatened, particularly the northern giraffe (Giraffa camelopardalis), which collectively have the smallest known wild population estimates. Among the three subspecies of the northern giraffe, the West African giraffe (Giraffa camelopardalis peralta) had declined to 49 individuals by 1996 and only recovered due to conservation efforts undertaken in the past 25 years, while the Kordofan giraffe (Giraffa camelopardalis antiquorum) remains at <2300 individuals distributed in small, isolated populations over a large geographical range in Central Africa. These combined factors could lead to genetically depauperated populations. We analyzed 119 mitochondrial sequences and 26 whole genomes of northern giraffe individuals to investigate their population structure and assess the recent demographic history and current genomic diversity of West African and Kordofan giraffe. Phylogenetic and population structure analyses separate the three subspecies of northern giraffe and suggest genetic differentiation between populations from eastern and western areas of the Kordofan giraffe’s range. Both West African and Kordofan giraffe show a gradual decline in effective population size over the last 10 ka and have moderate genome-wide heterozygosity compared to other giraffe species. Recent inbreeding levels are higher in the West African giraffe and in Kordofan giraffe from Garamba National Park, Democratic Republic of Congo. Although numbers for both West African and some populations of Kordofan giraffe have increased in recent years, the threat of habitat loss, climate change impacts, and illegal hunting persists. Thus, future conservation actions should consider close genetic monitoring of populations to detect and, where practical, counteract negative trends that might develop.
Collapse
|
17
|
McCulloch GA, Guhlin J, Dutoit L, Harrop TWR, Dearden PK, Waters JM. Genomic signatures of parallel alpine adaptation in recently evolved flightless insects. Mol Ecol 2021; 30:6677-6686. [PMID: 34592029 DOI: 10.1111/mec.16204] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/22/2021] [Revised: 09/16/2021] [Accepted: 09/27/2021] [Indexed: 12/01/2022]
Abstract
Natural selection along elevational gradients has potential to drive predictable adaptations across distinct lineages, but the extent of such repeated evolution remains poorly studied for many widespread alpine taxa. We present parallel genomic analyses of two recently evolved flightless alpine insect lineages to test for molecular signatures of repeated alpine adaptation. Specifically, we compare low-elevation vs. alpine stonefly ecotypes from parallel stream populations in which flightless upland ecotypes have been independently derived. We map 67,922 polymorphic genetic markers, generated across 176 Zelandoperla fenestrata specimens from two independent alpine stream populations in New Zealand's Rock and Pillar Range, to a newly developed plecopteran reference genome. Genome-wide scans revealed 31 regions with outlier single nucleotide polymorphisms (SNPs) differentiating lowland vs. alpine ecotypes in Lug Creek, and 37 regions with outliers differentiating ecotypes in Six Mile Creek. Of these regions, 13% (8/60) yielded outlier SNPs across both within-stream ecotype comparisons, implying comparable genomic shifts contribute to this repeated alpine adaptation. Candidate genes closely linked to repeated outlier regions include several with documented roles in insect wing-development (e.g., dishevelled), suggesting that they may contribute to repeated alpine wing reduction. Additional candidate genes have been shown to influence insect fecundity (e.g., ovo) and lifespan (e.g., Mrp4), implying that they might contribute to life history differentiation between upland and lowland ecotypes. Additional outlier genes have potential roles in the evolution of reproductive isolation among ecotypes (hedgehog and Desaturase 1). These results demonstrate how replicated outlier tests across independent lineages can potentially contribute to the discovery of genes underpinning repeated adaptation.
Collapse
Affiliation(s)
| | - Joseph Guhlin
- Genomics Aotearoa and Department of Biochemistry, University of Otago, Dunedin, New Zealand
| | - Ludovic Dutoit
- Department of Zoology, University of Otago, Dunedin, New Zealand
| | - Thomas W R Harrop
- Genomics Aotearoa and Department of Biochemistry, University of Otago, Dunedin, New Zealand
| | - Peter K Dearden
- Genomics Aotearoa and Department of Biochemistry, University of Otago, Dunedin, New Zealand
| | | |
Collapse
|
18
|
Kostanjšek R, Diderichsen B, Recknagel H, Gunde-Cimerman N, Gostinčar C, Fan G, Kordiš D, Trontelj P, Jiang H, Bolund L, Luo Y. Toward the massive genome of Proteus anguinus-illuminating longevity, regeneration, convergent evolution, and metabolic disorders. Ann N Y Acad Sci 2021; 1507:5-11. [PMID: 34480358 DOI: 10.1111/nyas.14686] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2021] [Revised: 08/13/2021] [Accepted: 08/17/2021] [Indexed: 12/27/2022]
Abstract
Deciphering the genetic code of organisms with unusual phenotypes can help answer fundamental biological questions and provide insight into mechanisms relevant to human biomedical research. The cave salamander Proteus anguinus (Urodela: Proteidae), also known as the olm, is an example of a species with unique morphological and physiological adaptations to its subterranean environment, including regenerative abilities, resistance to prolonged starvation, and a life span of more than 100 years. However, the structure and sequence of the olm genome is still largely unknown owing to its enormous size, estimated at nearly 50 gigabases. An international Proteus Genome Research Consortium has been formed to decipher the olm genome. This perspective provides the scientific and biomedical rationale for exploring the olm genome and outlines potential outcomes, challenges, and methodological approaches required to analyze and annotate the genome of this unique amphibian.
Collapse
Affiliation(s)
- Rok Kostanjšek
- Department of Biology, Biotechnical Faculty, University of Ljubljana, Ljubljana, Slovenia
| | - Børge Diderichsen
- Department of Molecular Biology and Genetics, Aarhus University, Aarhus, Denmark
| | - Hans Recknagel
- Department of Biology, Biotechnical Faculty, University of Ljubljana, Ljubljana, Slovenia
| | - Nina Gunde-Cimerman
- Department of Biology, Biotechnical Faculty, University of Ljubljana, Ljubljana, Slovenia
| | - Cene Gostinčar
- Department of Biology, Biotechnical Faculty, University of Ljubljana, Ljubljana, Slovenia.,Lars Bolund Institute of Regenerative Medicine, Qingdao-Europe Advanced Institute for Life Sciences, BGI-Qingdao, BGI-Shenzhen, Qingdao, China
| | - Guangyi Fan
- Lars Bolund Institute of Regenerative Medicine, Qingdao-Europe Advanced Institute for Life Sciences, BGI-Qingdao, BGI-Shenzhen, Qingdao, China
| | - Dušan Kordiš
- Department of Molecular and Biomedical Sciences, Jožef Stefan Institute, Ljubljana, Slovenia
| | - Peter Trontelj
- Department of Biology, Biotechnical Faculty, University of Ljubljana, Ljubljana, Slovenia
| | | | - Lars Bolund
- Lars Bolund Institute of Regenerative Medicine, Qingdao-Europe Advanced Institute for Life Sciences, BGI-Qingdao, BGI-Shenzhen, Qingdao, China.,Department of Biomedicine, Aarhus University, Aarhus, Denmark
| | - Yonglun Luo
- Lars Bolund Institute of Regenerative Medicine, Qingdao-Europe Advanced Institute for Life Sciences, BGI-Qingdao, BGI-Shenzhen, Qingdao, China.,Department of Biomedicine, Aarhus University, Aarhus, Denmark.,Steno Diabetes Center Aarhus, Aarhus University Hospital, Aarhus, Denmark
| |
Collapse
|
19
|
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
|
20
|
Rippe JP, Dixon G, Fuller ZL, Liao Y, Matz M. Environmental specialization and cryptic genetic divergence in two massive coral species from the Florida Keys Reef Tract. Mol Ecol 2021; 30:3468-3484. [PMID: 33894013 DOI: 10.1111/mec.15931] [Citation(s) in RCA: 16] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/18/2020] [Revised: 03/22/2021] [Accepted: 04/14/2021] [Indexed: 01/02/2023]
Abstract
Broadcast-spawning coral species have wide geographical ranges spanning strong environmental gradients, but it is unclear how much spatially varying selection these gradients actually impose. Strong divergent selection might present a considerable barrier for demographic exchange between disparate reef habitats. We investigated whether the cross-shelf gradient is associated with spatially varying selection in two common coral species, Montastraea cavernosa and Siderastrea siderea, in the Florida Keys. To this end, we generated a de novo genome assembly for M. cavernosa and used 2bRAD to genotype 20 juveniles and 20 adults of both species from each of the three reef zones to identify signatures of selection occurring within a single generation. Unexpectedly, each species was found to be composed of four genetically distinct lineages, with gene flow between them still ongoing but highly reduced in 13.0%-54.7% of the genome. Each species includes two sympatric lineages that are only found in the deep (20 m) habitat, while the other lineages are found almost exclusively on the shallower reefs (3-10 m). The two "shallow" lineages of M. cavernosa are also specialized for either nearshore or offshore: comparison between adult and juvenile cohorts indicates that cross-shelf migrants are more than twice as likely to die before reaching adulthood than local recruits. S. siderea and M. cavernosa are among the most ecologically successful species on the Florida Keys Reef Tract, and this work offers important insight into the genomic background of divergent selection and environmental specialization that may in part explain their resilience and broad environmental range.
Collapse
Affiliation(s)
- John P Rippe
- Department of Integrative Biology, University of Texas at Austin, Austin, TX, USA
| | - Groves Dixon
- Department of Integrative Biology, University of Texas at Austin, Austin, TX, USA
| | - Zachary L Fuller
- Department of Biological Sciences, Columbia University, New York, NY, USA
| | - Yi Liao
- Department of Integrative Biology, University of Texas at Austin, Austin, TX, USA.,Department of Ecology and Evolutionary Biology, University of California Irvine, Irvine, CA, USA
| | - Mikhail Matz
- Department of Integrative Biology, University of Texas at Austin, Austin, TX, USA
| |
Collapse
|
21
|
Sun J, Li R, Chen C, Sigwart JD, Kocot KM. Benchmarking Oxford Nanopore read assemblers for high-quality molluscan genomes. Philos Trans R Soc Lond B Biol Sci 2021; 376:20200160. [PMID: 33813888 PMCID: PMC8059532 DOI: 10.1098/rstb.2020.0160] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 12/31/2020] [Indexed: 12/14/2022] Open
Abstract
Choosing the optimum assembly approach is essential to achieving a high-quality genome assembly suitable for comparative and evolutionary genomic investigations. Significant recent progress in long-read sequencing technologies such as PacBio and Oxford Nanopore Technologies (ONT) has also brought about a large variety of assemblers. Although these have been extensively tested on model species such as Homo sapiens and Drosophila melanogaster, such benchmarking has not been done in Mollusca, which lacks widely adopted model species. Molluscan genomes are notoriously rich in repeats and are often highly heterozygous, making their assembly challenging. Here, we benchmarked 10 assemblers based on ONT raw reads from two published molluscan genomes of differing properties, the gastropod Chrysomallon squamiferum (356.6 Mb, 1.59% heterozygosity) and the bivalve Mytilus coruscus (1593 Mb, 1.94% heterozygosity). By optimizing the assembly pipeline, we greatly improved both genomes from previously published versions. Our results suggested that 40-50X of ONT reads are sufficient for high-quality genomes, with Flye being the recommended assembler for compact and less heterozygous genomes exemplified by C. squamiferum, while NextDenovo excelled for more repetitive and heterozygous molluscan genomes exemplified by M. coruscus. A phylogenomic analysis using the two updated genomes with 32 other published high-quality lophotrochozoan genomes resulted in maximum support across all nodes, and we show that improved genome quality also leads to more complete matrices for phylogenomic inferences. Our benchmarking will ensure efficiency in future assemblies for molluscs and perhaps also for other marine phyla with few genomes available. This article is part of the Theo Murphy meeting issue 'Molluscan genomics: broad insights and future directions for a neglected phylum'.
Collapse
Affiliation(s)
- Jin Sun
- Institute of Evolution and Marine Biodiversity, Key Laboratory of Mariculture (Ministry of Education), Ocean University of China, Qingdao 266003, People's Republic of China
| | - Runsheng Li
- Department of Infectious Diseases and Public Health, Jockey Club College of Veterinary Medicine and Life Sciences, City University of Hong Kong, Kowloon, Hong Kong, People's Republic of China
| | - Chong Chen
- X-STAR, Japan Agency for Marine-Earth Science and Technology (JAMSTEC), 2–15 Natsushima-cho, Yokosuka, Kanagawa Prefecture 237-0061, Japan
| | - Julia D. Sigwart
- Senckenberg Museum, 60325 Frankfurt, Germany
- Marine Laboratory Queen's University Belfast, Portaferry, BT22 1PF, Northern Ireland
| | - Kevin M. Kocot
- Department of Biological Sciences and Alabama Museum of Natural History, University of Alabama, Tuscaloosa, AL 35487, USA
| |
Collapse
|
22
|
RNA sequencing describes both population structure and plasticity-selection dynamics in a non-model fish. BMC Genomics 2021; 22:273. [PMID: 33858341 PMCID: PMC8048188 DOI: 10.1186/s12864-021-07592-4] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2021] [Accepted: 04/05/2021] [Indexed: 01/03/2023] Open
Abstract
Background Messenger RNA sequencing is becoming more common in studies of non-model species and is most often used for gene expression-based investigations. However, the method holds potential for numerous other applications as well—including analyses of alternative splicing, population structure, and signatures of selection. To maximize the utility of mRNA data sets, distinct analyses may be combined such as by exploring dynamics between gene expression with signatures of selection in the context of population structure. Here, we compare two published data sets describing two populations of a minnow species endemic to the San Francisco Estuary (Sacramento splittail, Pogonichthys macrolepidotus): a microsatellite data set showing population structure, and an mRNA whole transcriptome data set obtained after the two populations were exposed to a salinity challenge. We compared measures of population structure and genetic variation using single nucleotide polymorphisms (SNPs) called from mRNA from the whole transcriptome sequencing study with those patterns determined from microsatellites. For investigating plasticity and evolution, intra- and inter-population transcriptome plasticity was investigated with differential gene expression, differential exon usage, and gene expression variation. Outlier SNP analysis was also performed on the mRNA data set and signatures of selection and phenotypic plasticity were investigated on an individual-gene basis. Results We found that mRNA sequencing revealed patterns of population structure consistent with those found with microsatellites, but with lower magnitudes of genetic variation and population differentiation consistent with widespread purifying selection expected when using mRNA. In addition, within individual genes, phenotypic plasticity or signatures of selection were found in almost mutual exclusion (except heatr6, nfu1, slc22a6, sya, and mmp13). Conclusions These results show that an mRNA sequencing data set may have multiple uses, including describing population structure and for investigating the mechanistic interplay of evolution and plasticity in adaptation. MRNA sequencing thus complements traditional sequencing methods used for population genetics, in addition to its utility for describing phenotypic plasticity. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-021-07592-4.
Collapse
|
23
|
Liu P, Ewald J, Galvez JH, Head J, Crump D, Bourque G, Basu N, Xia J. Ultrafast functional profiling of RNA-seq data for nonmodel organisms. Genome Res 2021; 31:713-720. [PMID: 33731361 PMCID: PMC8015844 DOI: 10.1101/gr.269894.120] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2020] [Accepted: 02/18/2021] [Indexed: 12/02/2022]
Abstract
Computational time and cost remain a major bottleneck for RNA-seq data analysis of nonmodel organisms without reference genomes. To address this challenge, we have developed Seq2Fun, a novel, all-in-one, ultrafast tool to directly perform functional quantification of RNA-seq reads without transcriptome de novo assembly. The pipeline starts with raw read quality control: sequencing error correction, removing poly(A) tails, and joining overlapped paired-end reads. It then conducts a DNA-to-protein search by translating each read into all possible amino acid fragments and subsequently identifies possible homologous sequences in a well-curated protein database. Finally, the pipeline generates several informative outputs including gene abundance tables, pathway and species hit tables, an HTML report to visualize the results, and an output of clean reads annotated with mapped genes ready for downstream analysis. Seq2Fun does not have any intermediate steps of file writing and loading, making I/O very efficient. Seq2Fun is written in C++ and can run on a personal computer with a limited number of CPUs and memory. It can process >2,000,000 reads/min and is >120 times faster than conventional workflows based on de novo assembly, while maintaining high accuracy in our various test data sets.
Collapse
Affiliation(s)
- Peng Liu
- Faculty of Agricultural and Environmental Sciences, McGill University, Montreal, Quebec H9X 3V9, Canada
| | - Jessica Ewald
- Faculty of Agricultural and Environmental Sciences, McGill University, Montreal, Quebec H9X 3V9, Canada
| | - Jose Hector Galvez
- Department of Human Genetics, McGill University, Montreal, Quebec H3A 0C7, Canada.,Canadian Center for Computational Genomics, McGill University, Montreal, Quebec H3A 0G1, Canada
| | - Jessica Head
- Faculty of Agricultural and Environmental Sciences, McGill University, Montreal, Quebec H9X 3V9, Canada
| | - Doug Crump
- Environment and Climate Change Canada, National Wildlife Research Centre, Ottawa, Ontario K1A 0H3, Canada
| | - Guillaume Bourque
- Department of Human Genetics, McGill University, Montreal, Quebec H3A 0C7, Canada.,Canadian Center for Computational Genomics, McGill University, Montreal, Quebec H3A 0G1, Canada
| | - Niladri Basu
- Faculty of Agricultural and Environmental Sciences, McGill University, Montreal, Quebec H9X 3V9, Canada
| | - Jianguo Xia
- Faculty of Agricultural and Environmental Sciences, McGill University, Montreal, Quebec H9X 3V9, Canada.,Department of Human Genetics, McGill University, Montreal, Quebec H3A 0C7, Canada
| |
Collapse
|
24
|
Marx HE, Scheidt S, Barker MS, Dlugosch KM. TagSeq for gene expression in non-model plants: A pilot study at the Santa Rita Experimental Range NEON core site. APPLICATIONS IN PLANT SCIENCES 2020; 8:e11398. [PMID: 33304661 PMCID: PMC7705334 DOI: 10.1002/aps3.11398] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 04/11/2020] [Accepted: 08/20/2020] [Indexed: 05/12/2023]
Abstract
PREMISE TagSeq is a cost-effective approach for gene expression studies requiring a large number of samples. To date, TagSeq studies in plants have been limited to those with a high-quality reference genome. We tested the suitability of reference transcriptomes for TagSeq in non-model plants, as part of a study of natural gene expression variation at the Santa Rita Experimental Range National Ecological Observatory Network (NEON) core site. METHODS Tissue for TagSeq was sampled from multiple individuals of four species (Bouteloua aristidoides and Eragrostis lehmanniana [Poaceae], Tidestromia lanuginosa [Amaranthaceae], and Parkinsonia florida [Fabaceae]) at two locations on three dates (56 samples total). One sample per species was used to create a reference transcriptome via standard RNA-seq. TagSeq performance was assessed by recovery of reference loci, specificity of tag alignments, and variation among samples. RESULTS A high fraction of tags aligned to each reference and mapped uniquely. Expression patterns were quantifiable for tens of thousands of loci, which revealed consistent spatial differentiation in expression for all species. DISCUSSION TagSeq using de novo reference transcriptomes was an effective approach to quantifying gene expression in this study. Tags were highly locus specific and generated biologically informative profiles for four non-model plant species.
Collapse
Affiliation(s)
- Hannah E. Marx
- Department of Ecology and Evolutionary BiologyUniversity of ArizonaTucsonArizona85721USA
- Department of Ecology and Evolutionary BiologyUniversity of MichiganAnn ArborMichigan48109‐1048USA
| | - Stephen Scheidt
- Howard University2400 6th Street NWWashingtonD.C.20059USA
- Solar System Exploration DivisionNASA Goddard Space Flight CenterGreenbeltMaryland20771USA
- Center for Research and Exploration in Space Science and TechnologyNASA Goddard Space Flight CenterGreenbeltMaryland20771USA
| | - Michael S. Barker
- Department of Ecology and Evolutionary BiologyUniversity of ArizonaTucsonArizona85721USA
| | - Katrina M. Dlugosch
- Department of Ecology and Evolutionary BiologyUniversity of ArizonaTucsonArizona85721USA
| |
Collapse
|
25
|
Blommaert J. Genome size evolution: towards new model systems for old questions. Proc Biol Sci 2020; 287:20201441. [PMID: 32842932 PMCID: PMC7482279 DOI: 10.1098/rspb.2020.1441] [Citation(s) in RCA: 39] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/18/2020] [Accepted: 07/29/2020] [Indexed: 12/20/2022] Open
Abstract
Genome size (GS) variation is a fundamental biological characteristic; however, its evolutionary causes and consequences are the topic of ongoing debate. Whether GS is a neutral trait or one subject to selective pressures, and how strong these selective pressures are, may remain open questions. Fundamentally, the genomic sequences responsible for this variation directly impact the potential evolutionary outcomes and, equally, are the targets of different evolutionary pressures. For example, duplications and deletions of genic regions (large or small) can have immediate and drastic phenotypic effects, while an expansion or contraction of non-coding DNA is less likely to cause catastrophic phenotypic effects. However, in the long term, the accumulation or deletion of ncDNA is likely to have larger effects. Modern sequencing technologies are allowing for the dissection of these proximate causes, but a combination of these new technologies with more traditional evolutionary experiments and approaches could revolutionize this debate and potentially resolve many of these arguments. Here, I discuss an ambitious way forward for GS research, putting it in context of historical debates, theories and sometimes contradictory evidence, and highlighting the promise of combining new sequencing technologies and analytical developments with more traditional experimental evolution approaches.
Collapse
Affiliation(s)
- Julie Blommaert
- Department of Organismal Biology, Uppsala University, Uppsala, Sweden
| |
Collapse
|
26
|
Muir AP, Dubois SF, Ross RE, Firth LB, Knights AM, Lima FP, Seabra R, Corre E, Le Corguillé G, Nunes FLD. Seascape genomics reveals population isolation in the reef-building honeycomb worm, Sabellaria alveolata (L.). BMC Evol Biol 2020; 20:100. [PMID: 32778052 PMCID: PMC7418442 DOI: 10.1186/s12862-020-01658-9] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/05/2019] [Accepted: 07/17/2020] [Indexed: 12/11/2022] Open
Abstract
BACKGROUND Under the threat of climate change populations can disperse, acclimatise or evolve in order to avoid fitness loss. In light of this, it is important to understand neutral gene flow patterns as a measure of dispersal potential, but also adaptive genetic variation as a measure of evolutionary potential. In order to assess genetic variation and how this relates to environment in the honeycomb worm (Sabellaria alveolata (L.)), a reef-building polychaete that supports high biodiversity, we carried out RAD sequencing using individuals from along its complete latitudinal range. Patterns of neutral population genetic structure were compared to larval dispersal as predicted by ocean circulation modelling, and outlier analyses and genotype-environment association tests were used to attempt to identify loci under selection in relation to local temperature data. RESULTS We genotyped 482 filtered SNPs, from 68 individuals across nine sites, 27 of which were identified as outliers using BAYESCAN and ARLEQUIN. All outlier loci were potentially under balancing selection, despite previous evidence of local adaptation in the system. Limited gene flow was observed among reef-sites (FST = 0.28 ± 0.10), in line with the low dispersal potential identified by the larval dispersal models. The North Atlantic reef emerged as a distinct population and this was linked to high local larval retention and the effect of the North Atlantic Current on dispersal. CONCLUSIONS As an isolated population, with limited potential for natural genetic or demographic augmentation from other reefs, the North Atlantic site warrants conservation attention in order to preserve not only this species, but above all the crucial functional ecological roles that are associated with their bioconstructions. Our study highlights the utility of using seascape genomics to identify populations of conservation concern.
Collapse
Affiliation(s)
- Anna P Muir
- Conservation Biology Research Group, Department of Biological Sciences, University of Chester, Parkgate Road, Chester, CH1 4BJ, UK.
- Laboratoire des Sciences de l'Environnement Marin, LEMAR UMR 6539 CNRS/UBO/IRD/Ifremer, Université de Brest (UBO), Université Européenne de Bretagne (UEB), Institut Universitaire Européen de la Mer (IUEM), 29280, Plouzané, France.
| | - Stanislas F Dubois
- Ifremer, DYNECO, Laboratory of Coastal Benthic Ecology, F-29280, Plouzané, France
| | - Rebecca E Ross
- Marine Biology and Ecology Research Centre, School of Biological and Marine Sciences, University of Plymouth, Plymouth, PL4 8AA, UK
- Institute of Marine Research, 1870 Nordnes, 5817, Bergen, Norway
| | - Louise B Firth
- School of Biological and Marine Sciences, University of Plymouth, Plymouth, PL4 8AA, UK
| | - Antony M Knights
- Marine Biology and Ecology Research Centre, School of Biological and Marine Sciences, University of Plymouth, Plymouth, PL4 8AA, UK
| | - Fernando P Lima
- CIBIO/InBIO, Centro de Investigação em Biodiversidade e Recursos Genéticos, Universidade do Porto, Vairão, Portugal
| | - Rui Seabra
- CIBIO/InBIO, Centro de Investigação em Biodiversidade e Recursos Genéticos, Universidade do Porto, Vairão, Portugal
| | - Erwan Corre
- CNRS, Sorbonne Université, FR2424, ABiMS, Station Biologique de Roscoff, 29680, Roscoff, France
| | - Gildas Le Corguillé
- CNRS, Sorbonne Université, FR2424, ABiMS, Station Biologique de Roscoff, 29680, Roscoff, France
| | - Flavia L D Nunes
- Laboratoire des Sciences de l'Environnement Marin, LEMAR UMR 6539 CNRS/UBO/IRD/Ifremer, Université de Brest (UBO), Université Européenne de Bretagne (UEB), Institut Universitaire Européen de la Mer (IUEM), 29280, Plouzané, France
- Ifremer, DYNECO, Laboratory of Coastal Benthic Ecology, F-29280, Plouzané, France
| |
Collapse
|
27
|
Studivan MS, Voss JD. Transcriptomic plasticity of mesophotic corals among natural populations and transplants of
Montastraea cavernosa
in the Gulf of Mexico and Belize. Mol Ecol 2020; 29:2399-2415. [DOI: 10.1111/mec.15495] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2019] [Revised: 05/14/2020] [Accepted: 05/28/2020] [Indexed: 02/06/2023]
Affiliation(s)
- Michael S. Studivan
- Harbor Branch Oceanographic Institute Florida Atlantic University Fort Pierce FL USA
- Cooperative Institute for Marine and Atmospheric Studies University of Miami Rosenstiel School of Marine and Atmospheric Sciences Miami FL USA
| | - Joshua D. Voss
- Harbor Branch Oceanographic Institute Florida Atlantic University Fort Pierce FL USA
| |
Collapse
|
28
|
Garrison NL, Brewer MS, Bond JE. Shifting evolutionary sands: transcriptome characterization of the Aptostichus atomarius species complex. BMC Evol Biol 2020; 20:68. [PMID: 32539685 PMCID: PMC7294663 DOI: 10.1186/s12862-020-01606-7] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2019] [Accepted: 03/20/2020] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND Mygalomorph spiders represent a diverse, yet understudied lineage for which genomic level data has only recently become accessible through high-throughput genomic and transcriptomic sequencing methods. The Aptostichus atomarius species complex (family Euctenizidae) includes two coastal dune endemic members, each with inland sister species - affording exploration of dune adaptation associated patterns at the transcriptomic level. We apply an RNAseq approach to examine gene family conservation across the species complex and test for patterns of positive selection along branches leading to dune endemic species. RESULTS An average of ~ 44,000 contigs were assembled for eight spiders representing dune (n = 2), inland (n = 4), and atomarius species complex outgroup taxa (n = 2). Transcriptomes were estimated to be 64% complete on average with 77 spider reference orthologs missing from all taxa. Over 18,000 orthologous gene clusters were identified within the atomarius complex members, > 5000 were detected in all species, and ~ 4700 were shared between species complex members and outgroup Aptostichus species. Gene family analysis with the FUSTr pipeline identified 47 gene families appearing to be under selection in the atomarius ingroup; four of the five top clusters include sequences strongly resembling other arthropod venom peptides. The COATS pipeline identified six gene clusters under positive selection on branches leading to dune species, three of which reflected the preferred species tree. Genes under selection were identified as Cytochrome P450 2c15 (also recovered in the FUSTr analysis), Niemann 2 Pick C1-like, and Kainate 2 isoform X1. CONCLUSIONS We have generated eight draft transcriptomes for a closely related and ecologically diverse group of trapdoor spiders, identifying venom gene families potentially under selection across the Aptostichus atomarius complex and chemosensory-associated gene families under selection in dune endemic lineages.
Collapse
Affiliation(s)
- Nicole L. Garrison
- School of Fisheries, Aquaculture, and Aquatic Sciences, Auburn University, 203 Swingle Hall, Auburn, AL 36849 USA
| | - Michael S. Brewer
- Department of Biology, East Carolina University, Howell Science Complex N407, 1000 E 5th St, Greenville, NC 27858 USA
| | - Jason E. Bond
- Department of Entomology and Nematology, University of California Davis, Academic Surge Building 1282, Davis, CA 95616-5270 USA
| |
Collapse
|
29
|
Scalzitti N, Jeannin-Girardon A, Collet P, Poch O, Thompson JD. A benchmark study of ab initio gene prediction methods in diverse eukaryotic organisms. BMC Genomics 2020; 21:293. [PMID: 32272892 PMCID: PMC7147072 DOI: 10.1186/s12864-020-6707-9] [Citation(s) in RCA: 36] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2019] [Accepted: 03/30/2020] [Indexed: 02/02/2023] Open
Abstract
Background The draft genome assemblies produced by new sequencing technologies present important challenges for automatic gene prediction pipelines, leading to less accurate gene models. New benchmark methods are needed to evaluate the accuracy of gene prediction methods in the face of incomplete genome assemblies, low genome coverage and quality, complex gene structures, or a lack of suitable sequences for evidence-based annotations. Results We describe the construction of a new benchmark, called G3PO (benchmark for Gene and Protein Prediction PrOgrams), designed to represent many of the typical challenges faced by current genome annotation projects. The benchmark is based on a carefully validated and curated set of real eukaryotic genes from 147 phylogenetically disperse organisms, and a number of test sets are defined to evaluate the effects of different features, including genome sequence quality, gene structure complexity, protein length, etc. We used the benchmark to perform an independent comparative analysis of the most widely used ab initio gene prediction programs and identified the main strengths and weaknesses of the programs. More importantly, we highlight a number of features that could be exploited in order to improve the accuracy of current prediction tools. Conclusions The experiments showed that ab initio gene structure prediction is a very challenging task, which should be further investigated. We believe that the baseline results associated with the complex gene test sets in G3PO provide useful guidelines for future studies.
Collapse
Affiliation(s)
- Nicolas Scalzitti
- Department of Computer Science, ICube, CNRS, University of Strasbourg, Strasbourg, France
| | - Anne Jeannin-Girardon
- Department of Computer Science, ICube, CNRS, University of Strasbourg, Strasbourg, France
| | - Pierre Collet
- Department of Computer Science, ICube, CNRS, University of Strasbourg, Strasbourg, France
| | - Olivier Poch
- Department of Computer Science, ICube, CNRS, University of Strasbourg, Strasbourg, France
| | - Julie D Thompson
- Department of Computer Science, ICube, CNRS, University of Strasbourg, Strasbourg, France.
| |
Collapse
|
30
|
Jackson JM, Pimsler ML, Oyen KJ, Strange JP, Dillon ME, Lozier JD. Local adaptation across a complex bioclimatic landscape in two montane bumble bee species. Mol Ecol 2020; 29:920-939. [PMID: 32031739 DOI: 10.1111/mec.15376] [Citation(s) in RCA: 26] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2019] [Revised: 01/21/2020] [Accepted: 02/02/2020] [Indexed: 12/21/2022]
Abstract
Understanding evolutionary responses to variation in temperature and precipitation across species ranges is of fundamental interest given ongoing climate change. The importance of temperature and precipitation for multiple aspects of bumble bee (Bombus) biology, combined with large geographic ranges that expose populations to diverse environmental pressures, make these insects well-suited for studying local adaptation. Here, we analyzed genome-wide sequence data from two widespread bumble bees, Bombus vosnesenskii and Bombus vancouverensis, using multiple environmental association analysis methods to investigate climate adaptation across latitude and altitude. The strongest signatures of selection were observed in B. vancouverensis, but despite unique responses between species for most loci, we detected several shared responses. Genes relating to neural and neuromuscular function and ion transport were especially evident with respect to temperature variables, while genes relating to cuticle formation, tracheal and respiratory system development, and homeostasis were associated with precipitation variables. Our data thus suggest that adaptive responses for tolerating abiotic variation are likely to be complex, but that several parallels among species can emerge even for these complex traits and landscapes. Results provide the framework for future work into mechanisms of thermal and desiccation tolerance in bumble bees and a set of genomic targets that might be monitored for future conservation efforts.
Collapse
Affiliation(s)
- Jason M Jackson
- Department of Biological Sciences, The University of Alabama, Tuscaloosa, AL, USA
| | - Meaghan L Pimsler
- Department of Biological Sciences, The University of Alabama, Tuscaloosa, AL, USA
| | - Kennan J Oyen
- Department of Zoology & Physiology and Program in Ecology, University of Wyoming, Laramie, WY, USA.,Department of Biological Sciences, University of Cincinnati, Cincinnati, OH, USA
| | - James P Strange
- Department of Entomology, The Ohio State University, Columbus, OH, USA
| | - Michael E Dillon
- Department of Zoology & Physiology and Program in Ecology, University of Wyoming, Laramie, WY, USA
| | - Jeffrey D Lozier
- Department of Biological Sciences, The University of Alabama, Tuscaloosa, AL, USA
| |
Collapse
|
31
|
Strader ME, Wong JM, Hofmann GE. Ocean acidification promotes broad transcriptomic responses in marine metazoans: a literature survey. Front Zool 2020; 17:7. [PMID: 32095155 PMCID: PMC7027112 DOI: 10.1186/s12983-020-0350-9] [Citation(s) in RCA: 41] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2018] [Accepted: 01/06/2020] [Indexed: 01/16/2023] Open
Abstract
For nearly a decade, the metazoan-focused research community has explored the impacts of ocean acidification (OA) on marine animals, noting that changes in ocean chemistry can impact calcification, metabolism, acid-base regulation, stress response and behavior in organisms that hold high ecological and economic value. Because OA interacts with several key physiological processes in marine organisms, transcriptomics has become a widely-used method to characterize whole organism responses on a molecular level as well as inform mechanisms that explain changes in phenotypes observed in response to OA. In the past decade, there has been a notable rise in studies that examine transcriptomic responses to OA in marine metazoans, and here we attempt to summarize key findings across these studies. We find that organisms vary dramatically in their transcriptomic responses to pH although common patterns are often observed, including shifts in acid-base ion regulation, metabolic processes, calcification and stress response mechanisms. We also see a rise in transcriptomic studies examining organismal response to OA in a multi-stressor context, often reporting synergistic effects of OA and temperature. In addition, there is an increase in studies that use transcriptomics to examine the evolutionary potential of organisms to adapt to OA conditions in the future through population and transgenerational experiments. Overall, the literature reveals complex organismal responses to OA, in which some organisms will face more dramatic consequences than others. This will have wide-reaching impacts on ocean communities and ecosystems as a whole.
Collapse
Affiliation(s)
- Marie E Strader
- 1Department of Ecology, Evolution and Marine Biology, University of California Santa Barbara, Santa Barbara, CA 93106 USA.,2Department of Biological Sciences, Auburn University, Auburn, AL 36849 USA
| | - Juliet M Wong
- 1Department of Ecology, Evolution and Marine Biology, University of California Santa Barbara, Santa Barbara, CA 93106 USA.,3Present address: Department of Biological Sciences, Florida International University, North Miami, FL 33181 USA
| | - Gretchen E Hofmann
- 1Department of Ecology, Evolution and Marine Biology, University of California Santa Barbara, Santa Barbara, CA 93106 USA
| |
Collapse
|
32
|
Patton AH, Margres MJ, Stahlke AR, Hendricks S, Lewallen K, Hamede RK, Ruiz-Aravena M, Ryder O, McCallum HI, Jones ME, Hohenlohe PA, Storfer A. Contemporary Demographic Reconstruction Methods Are Robust to Genome Assembly Quality: A Case Study in Tasmanian Devils. Mol Biol Evol 2020; 36:2906-2921. [PMID: 31424552 PMCID: PMC6878949 DOI: 10.1093/molbev/msz191] [Citation(s) in RCA: 55] [Impact Index Per Article: 13.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
Reconstructing species’ demographic histories is a central focus of molecular ecology and evolution. Recently, an expanding suite of methods leveraging either the sequentially Markovian coalescent (SMC) or the site-frequency spectrum has been developed to reconstruct population size histories from genomic sequence data. However, few studies have investigated the robustness of these methods to genome assemblies of varying quality. In this study, we first present an improved genome assembly for the Tasmanian devil using the Chicago library method. Compared with the original reference genome, our new assembly reduces the number of scaffolds (from 35,975 to 10,010) and increases the scaffold N90 (from 0.101 to 2.164 Mb). Second, we assess the performance of four contemporary genomic methods for inferring population size history (PSMC, MSMC, SMC++, Stairway Plot), using the two devil genome assemblies as well as simulated, artificially fragmented genomes that approximate the hypothesized demographic history of Tasmanian devils. We demonstrate that each method is robust to assembly quality, producing similar estimates of Ne when simulated genomes were fragmented into up to 5,000 scaffolds. Overall, methods reliant on the SMC are most reliable between ∼300 generations before present (gbp) and 100 kgbp, whereas methods exclusively reliant on the site-frequency spectrum are most reliable between the present and 30 gbp. Our results suggest that when used in concert, genomic methods for reconstructing species’ effective population size histories 1) can be applied to nonmodel organisms without highly contiguous reference genomes, and 2) are capable of detecting independently documented effects of historical geological events.
Collapse
Affiliation(s)
- Austin H Patton
- School of Biological Sciences, Washington State University, Pullman, WA
| | - Mark J Margres
- School of Biological Sciences, Washington State University, Pullman, WA.,Department of Organismic and Evolutionary Biology, Harvard University, MA
| | - Amanda R Stahlke
- Institute for Bioinformatics and Evolutionary Studies, University of Idaho, Moscow, ID
| | - Sarah Hendricks
- Institute for Bioinformatics and Evolutionary Studies, University of Idaho, Moscow, ID
| | - Kevin Lewallen
- Institute for Bioinformatics and Evolutionary Studies, University of Idaho, Moscow, ID
| | - Rodrigo K Hamede
- School of Natural Sciences, University of Tasmania, Hobart, Australia
| | | | - Oliver Ryder
- Institute for Conservation Research, San Diego, CA
| | | | - Menna E Jones
- School of Natural Sciences, University of Tasmania, Hobart, Australia
| | - Paul A Hohenlohe
- Institute for Bioinformatics and Evolutionary Studies, University of Idaho, Moscow, ID
| | - Andrew Storfer
- School of Biological Sciences, Washington State University, Pullman, WA
| |
Collapse
|
33
|
Dillon ME, Lozier JD. Adaptation to the abiotic environment in insects: the influence of variability on ecophysiology and evolutionary genomics. CURRENT OPINION IN INSECT SCIENCE 2019; 36:131-139. [PMID: 31698151 DOI: 10.1016/j.cois.2019.09.003] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/09/2019] [Revised: 09/10/2019] [Accepted: 09/10/2019] [Indexed: 06/10/2023]
Abstract
Advances in tools to gather environmental, phenotypic, and molecular data have accelerated our ability to detect abiotic drivers of variation across the genome-to-phenome spectrum in model and non-model insects. However, differences in the spatial and temporal resolution of these data sets may create gaps in our understanding of linkages between environment, genotype, and phenotype that yield missed or misleading results about adaptive variation. In this review we highlight sources of variability that might impact studies of phenotypic and 'omic environmental adaptation, challenges to collecting data at relevant scales, and possible solutions that link intensive fine-scale reductionist studies of mechanisms to large-scale biogeographic patterns.
Collapse
Affiliation(s)
- Michael E Dillon
- Department of Zoology & Physiology and Program in Ecology, The University of Wyoming, Laramie, Wyoming 82071, USA.
| | - Jeffrey D Lozier
- Department of Biological Sciences, The University of Alabama, Box 870344, Tuscaloosa, Alabama 35487, USA
| |
Collapse
|
34
|
Parkinson JE, Baker AC, Baums IB, Davies SW, Grottoli AG, Kitchen SA, Matz MV, Miller MW, Shantz AA, Kenkel CD. Molecular tools for coral reef restoration: Beyond biomarker discovery. Conserv Lett 2019. [DOI: 10.1111/conl.12687] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/01/2023] Open
Affiliation(s)
- John Everett Parkinson
- SECORE International Miami Florida
- Department of Integrative BiologyUniversity of South Florida Tampa Florida
| | - Andrew C. Baker
- Department of Marine Biology and Ecology, Rosenstiel School of Marine and Atmospheric ScienceUniversity of Miami Miami Florida
| | - Iliana B. Baums
- Department of BiologyPennsylvania State University University Park Pennsylvania
| | | | | | - Sheila A. Kitchen
- Department of BiologyPennsylvania State University University Park Pennsylvania
| | - Mikhail V. Matz
- Department of Integrative BiologyUniversity of Texas at Austin Austin Texas
| | | | - Andrew A. Shantz
- Department of Marine Biology and Ecology, Rosenstiel School of Marine and Atmospheric ScienceUniversity of Miami Miami Florida
| | - Carly D. Kenkel
- Department of Biological SciencesUniversity of Southern California Los Angeles California
| |
Collapse
|
35
|
Conservation Genomics in a Changing Arctic. Trends Ecol Evol 2019; 35:149-162. [PMID: 31699414 DOI: 10.1016/j.tree.2019.09.008] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2019] [Revised: 09/13/2019] [Accepted: 09/17/2019] [Indexed: 12/25/2022]
Abstract
Although logistically challenging to study, the Arctic is a bellwether for global change and is becoming a model for questions pertinent to the persistence of biodiversity. Disruption of Arctic ecosystems is accelerating, with impacts ranging from mixing of biotic communities to individual behavioral responses. Understanding these changes is crucial for conservation and sustainable economic development. Genomic approaches are providing transformative insights into biotic responses to environmental change, but have seen limited application in the Arctic due to a series of limitations. To meet the promise of genome analyses, we urge rigorous development of biorepositories from high latitudes to provide essential libraries to improve the conservation, monitoring, and management of Arctic ecosystems through genomic approaches.
Collapse
|
36
|
MARTIN CHRISTOPHERH, RICHARDS EMILIEJ. The paradox behind the pattern of rapid adaptive radiation: how can the speciation process sustain itself through an early burst? ANNUAL REVIEW OF ECOLOGY, EVOLUTION, AND SYSTEMATICS 2019; 50:569-593. [PMID: 36237480 PMCID: PMC9555815 DOI: 10.1146/annurev-ecolsys-110617-062443] [Citation(s) in RCA: 41] [Impact Index Per Article: 8.2] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 08/27/2023]
Abstract
Rapid adaptive radiation poses a distinct question apart from speciation and adaptation: what happens after one speciation event? That is, how are some lineages able to continue speciating through a rapid burst? This question connects global macroevolutionary patterns to microevolutionary processes. Here we review major features of rapid radiations in nature and their mismatch with theoretical models and what is currently known about speciation mechanisms. Rapid radiations occur on three major diversification axes - species richness, phenotypic disparity, and ecological diversity - with exceptional outliers on each axis. The paradox is that the hallmark early stage of adaptive radiation, a rapid burst of speciation and niche diversification, is contradicted by most existing speciation models which instead predict continuously decelerating speciation rates and niche subdivision through time. Furthermore, while speciation mechanisms such as magic traits, phenotype matching, and physical linkage of co-adapted alleles promote speciation, it is often not discussed how these mechanisms could promote multiple speciation events in rapid succession. Additional mechanisms beyond ecological opportunity are needed to understand how rapid radiations occur. We review the evidence for five emerging theories: 1) the 'transporter' hypothesis: introgression and the ancient origins of adaptive alleles, 2) the 'signal complexity' hypothesis: the dimensionality of sexual traits, 3) the connectivity of fitness landscapes, 4) 'diversity begets diversity', and 5) flexible stem/'plasticity first'. We propose new questions and predictions to guide future work on the mechanisms underlying the rare origins of rapid radiation.
Collapse
Affiliation(s)
- CHRISTOPHER H. MARTIN
- Department of Biology, University of North Carolina at Chapel Hill, NC, USA
- Integrative Biology and Museum of Vertebrate Zoology, University of California, Berkeley, CA, USA
| | - EMILIE J. RICHARDS
- Department of Biology, University of North Carolina at Chapel Hill, NC, USA
- Integrative Biology and Museum of Vertebrate Zoology, University of California, Berkeley, CA, USA
| |
Collapse
|
37
|
Laskowski KL, Doran C, Bierbach D, Krause J, Wolf M. Naturally clonal vertebrates are an untapped resource in ecology and evolution research. Nat Ecol Evol 2019; 3:161-169. [PMID: 30692622 DOI: 10.1038/s41559-018-0775-0] [Citation(s) in RCA: 16] [Impact Index Per Article: 3.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2018] [Accepted: 11/29/2018] [Indexed: 11/09/2022]
Abstract
Science requires replication. The development of many cloned or isogenic model organisms is a testament to this. But researchers are reluctant to use these traditional animal model systems for certain questions in evolution or ecology research, because of concerns over relevance or inbreeding. It has largely been overlooked that there are a substantial number of vertebrate species that reproduce clonally in nature. Here we highlight how use of these naturally evolved, phenotypically complex animals can push the boundaries of traditional experimental design and contribute to answering fundamental questions in the fields of ecology and evolution.
Collapse
Affiliation(s)
- Kate L Laskowski
- Department of Biology & Ecology of Fishes, Leibniz-Institute of Freshwater Ecology & Inland Fisheries, Berlin, Germany.
| | - Carolina Doran
- Department of Biology & Ecology of Fishes, Leibniz-Institute of Freshwater Ecology & Inland Fisheries, Berlin, Germany
| | - David Bierbach
- Department of Biology & Ecology of Fishes, Leibniz-Institute of Freshwater Ecology & Inland Fisheries, Berlin, Germany
| | - Jens Krause
- Department of Biology & Ecology of Fishes, Leibniz-Institute of Freshwater Ecology & Inland Fisheries, Berlin, Germany.,Faculty of Life Sciences, Humboldt-Universität zu Berlin, Berlin, Germany
| | - Max Wolf
- Department of Biology & Ecology of Fishes, Leibniz-Institute of Freshwater Ecology & Inland Fisheries, Berlin, Germany
| |
Collapse
|
38
|
De-Kayne R, Feulner PGD. A European Whitefish Linkage Map and Its Implications for Understanding Genome-Wide Synteny Between Salmonids Following Whole Genome Duplication. G3 (BETHESDA, MD.) 2018; 8:3745-3755. [PMID: 30297382 PMCID: PMC6288842 DOI: 10.1534/g3.118.200552] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/29/2018] [Accepted: 10/03/2018] [Indexed: 12/11/2022]
Abstract
Genomic datasets continue to increase in number due to the ease of production for a wider selection of species including non-model organisms. For many of these species, especially those with large or polyploid genomes, highly contiguous and well-annotated genomes are still rare due to the complexity and cost involved in their assembly. As a result, a common starting point for genomic work in non-model species is the production of a linkage map. Dense linkage maps facilitate the analysis of genomic data in a variety of ways, from broad scale observations regarding genome structure e.g., chromosome number and type or sex-related structural differences, to fine scale patterns e.g., recombination rate variation and co-localization of differentiated regions. Here we present both sex-averaged and sex-specific linkage maps for Coregonus sp. "Albock", a member of the European whitefish lineage (C. lavaretus spp. complex), containing 5395 single nucleotide polymorphism (SNP) loci across 40 linkage groups to facilitate future investigation into the genomic basis of whitefish adaptation and speciation. The map was produced using restriction-site associated digestion (RAD) sequencing data from two wild-caught parents and 156 F1 offspring. We discuss the differences between our sex-averaged and sex-specific maps and identify genome-wide synteny between C. sp. "Albock" and Atlantic Salmon (Salmo salar), which have diverged following the salmonid-specific whole genome duplication. Our analysis confirms that many patterns of synteny observed between Atlantic Salmon and Oncorhynchus and Salvelinus species are also shared by members of the Coregoninae subfamily. We also show that regions known for their species-specific rediploidization history can pose challenges for synteny identification since these regions have diverged independently in each salmonid species following the salmonid-specific whole genome duplication. The European whitefish map provided here will enable future studies to understand the distribution of loci of interest, e.g., FST outliers, along the whitefish genome as well as assisting with the de novo assembly of a whitefish reference genome.
Collapse
Affiliation(s)
- Rishi De-Kayne
- Department of Fish Ecology and Evolution, Centre of Ecology, Evolution and Biogeochemistry, EAWAG Swiss Federal Institute of Aquatic Science and Technology, Switzerland
- Division of Aquatic Ecology and Evolution, Institute of Ecology and Evolution, University of Bern, Switzerland
| | - Philine G D Feulner
- Department of Fish Ecology and Evolution, Centre of Ecology, Evolution and Biogeochemistry, EAWAG Swiss Federal Institute of Aquatic Science and Technology, Switzerland
- Division of Aquatic Ecology and Evolution, Institute of Ecology and Evolution, University of Bern, Switzerland
| |
Collapse
|