1
|
Schwaner C, Farhat S, Boutet I, Tanguy A, Barbosa M, Grouzdev D, Pales Espinosa E, Allam B. Combination of RNAseq and RADseq to Identify Physiological and Adaptive Responses to Acidification in the Eastern Oyster (Crassostrea virginica). MARINE BIOTECHNOLOGY (NEW YORK, N.Y.) 2023; 25:997-1019. [PMID: 37864760 DOI: 10.1007/s10126-023-10255-y] [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: 02/27/2023] [Accepted: 09/29/2023] [Indexed: 10/23/2023]
Abstract
Ocean acidification (OA) is a major stressor threatening marine calcifiers, including the eastern oyster (Crassostrea virginica). In this paper, we provide insight into the molecular mechanisms associated with resilience to OA, with the dual intentions of probing both acclimation and adaptation potential in this species. C. virginica were spawned, and larvae were reared in control or acidified conditions immediately after fertilization. RNA samples were collected from larvae and juveniles, and DNA samples were collected from juveniles after undergoing OA-induced mortality and used to contrast gene expression (RNAseq) and SNP (ddRADseq) profiles from animals reared under both conditions. Results showed convergence of evidence from both approaches, particularly in genes involved in biomineralization that displayed significant changes in variant frequencies and gene expression levels among juveniles that survived acidification as compared to controls. Downregulated genes were related to immune processes, supporting previous studies demonstrating a reduction in immunity from exposure to OA. Acclimation to OA via regulation of gene expression might confer short-term resilience to immediate threats; however, the costs may not be sustainable, underscoring the importance of selection of resilient genotypes. Here, we identified SNPs associated with survival under OA conditions, suggesting that this commercially and ecologically important species might have the genetic variation needed for adaptation to future acidification. The identification of genetic features associated with OA resilience is a highly-needed step for the development of marker-assisted selection of oyster stocks for aquaculture and restoration activities.
Collapse
Affiliation(s)
- Caroline Schwaner
- School of Marine and Atmospheric Sciences, Stony Brook University, Stony Brook, NY, 11790, USA
| | - Sarah Farhat
- School of Marine and Atmospheric Sciences, Stony Brook University, Stony Brook, NY, 11790, USA
- Institut Systématique Evolution Biodiversité (ISYEB), Muséum national d'Histoire naturelle, CNRS, Sorbonne Université, EPHE, Université des Antilles, Paris, France
| | - Isabelle Boutet
- Station Biologique de Roscoff, CNRS/Sorbonne Université, Place Georges Teissier 29680, Roscoff, France
| | - Arnaud Tanguy
- Station Biologique de Roscoff, CNRS/Sorbonne Université, Place Georges Teissier 29680, Roscoff, France
| | - Michelle Barbosa
- School of Marine and Atmospheric Sciences, Stony Brook University, Stony Brook, NY, 11790, USA
| | - Denis Grouzdev
- School of Marine and Atmospheric Sciences, Stony Brook University, Stony Brook, NY, 11790, USA
| | | | - Bassem Allam
- School of Marine and Atmospheric Sciences, Stony Brook University, Stony Brook, NY, 11790, USA.
| |
Collapse
|
2
|
Maduna SN, Jónsdóttir ÓDB, Imsland AKD, Gíslason D, Reynolds P, Kapari L, Hangstad TA, Meier K, Hagen SB. Genomic Signatures of Local Adaptation under High Gene Flow in Lumpfish-Implications for Broodstock Provenance Sourcing and Larval Production. Genes (Basel) 2023; 14:1870. [PMID: 37895225 PMCID: PMC10606024 DOI: 10.3390/genes14101870] [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] [Received: 08/21/2023] [Revised: 09/20/2023] [Accepted: 09/23/2023] [Indexed: 10/29/2023] Open
Abstract
Aquaculture of the lumpfish (Cyclopterus lumpus L.) has become a large, lucrative industry owing to the escalating demand for "cleaner fish" to minimise sea lice infestations in Atlantic salmon mariculture farms. We used over 10K genome-wide single nucleotide polymorphisms (SNPs) to investigate the spatial patterns of genomic variation in the lumpfish along the coast of Norway and across the North Atlantic. Moreover, we applied three genome scans for outliers and two genotype-environment association tests to assess the signatures and patterns of local adaptation under extensive gene flow. With our 'global' sampling regime, we found two major genetic groups of lumpfish, i.e., the western and eastern Atlantic. Regionally in Norway, we found marginal evidence of population structure, where the population genomic analysis revealed a small portion of individuals with a different genetic ancestry. Nevertheless, we found strong support for local adaption under high gene flow in the Norwegian lumpfish and identified over 380 high-confidence environment-associated loci linked to gene sets with a key role in biological processes associated with environmental pressures and embryonic development. Our results bridge population genetic/genomics studies with seascape genomics studies and will facilitate genome-enabled monitoring of the genetic impacts of escapees and allow for genetic-informed broodstock selection and management in Norway.
Collapse
Affiliation(s)
- Simo Njabulo Maduna
- Department of Ecosystems in the Barents Region, Svanhovd Research Station, Norwegian Institute of Bioeconomy Research, 9925 Svanvik, Norway;
| | | | - Albert Kjartan Dagbjartarson Imsland
- Akvaplan-Niva Iceland Office, Akralind 6, 201 Kópavogur, Iceland; (Ó.D.B.J.); (A.K.D.I.)
- Department of Biological Sciences, High Technology Centre, University of Bergen, 5020 Bergen, Norway
| | | | | | - Lauri Kapari
- Akvaplan-Niva, Framsenteret, 9296 Tromsø, Norway;
| | | | | | - Snorre B. Hagen
- Department of Ecosystems in the Barents Region, Svanhovd Research Station, Norwegian Institute of Bioeconomy Research, 9925 Svanvik, Norway;
| |
Collapse
|
3
|
Orland C, Escalona M, Sahasrabudhe R, Marimuthu MPA, Nguyen O, Beraut E, Marshman B, Moore J, Raimondi P, Shapiro B. A Draft Reference Genome Assembly of the Critically Endangered Black Abalone, Haliotis cracherodii. J Hered 2022; 113:665-672. [PMID: 35567593 PMCID: PMC9709981 DOI: 10.1093/jhered/esac024] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/18/2022] [Accepted: 05/11/2022] [Indexed: 12/16/2022] Open
Abstract
The once abundant black abalone, Haliotis cracherodii, is a large, long-lived grazing marine mollusk that inhabits the rocky intertidal along the coast of California. The species has experienced dramatic declines since the mid-1980s largely due to the fatal bacterial disease called withering syndrome, leading to the collapse of an economically important fishery and to its inclusion into the IUCN listing as a critically endangered species. In some places impacted by the disease, populations of black abalone have declined by more than 90%, prompting population crashes associated with very little recruitment of new individuals and changes to intertidal communities. Habitats that were dominated by crustose coralline algae and bare rock have become dominated instead by fleshy algae and sessile invertebrates. Here, we present the first high-quality black abalone reference genome, assembled with PacBio HiFi long-reads and assembled with Dovetail Omni-C data to generate a scaffold-level assembly. The black abalone reference genome will be an essential resource in understanding the evolutionary history of this species as well as for exploring its current levels of genetic diversity and establishing future management and restoration plans.
Collapse
Affiliation(s)
- Chloé Orland
- Address correspondence to C. Orland at the address above, or e-mail:
| | | | - Ruta Sahasrabudhe
- UC Davis Genome Center, DNA Technologies and Expression Analysis Cores, University of California, Davis, Davis, CA, USA
| | - Mohan P A Marimuthu
- UC Davis Genome Center, DNA Technologies and Expression Analysis Cores, University of California, Davis, Davis, CA, USA
| | - Oanh Nguyen
- UC Davis Genome Center, DNA Technologies and Expression Analysis Cores, University of California, Davis, Davis, CA, USA
| | - Eric Beraut
- Ecology and Evolutionary Biology Department, University of California Santa Cruz, Santa Cruz, CA, USA
| | - Blythe Marshman
- Marine Laboratory, University of California Davis, Davis, CA, USA
| | - James Moore
- Marine Laboratory, University of California Davis, Davis, CA, USA
| | - Peter Raimondi
- Ecology and Evolutionary Biology Department, University of California Santa Cruz, Santa Cruz, CA, USA
| | - Beth Shapiro
- Ecology and Evolutionary Biology Department, University of California Santa Cruz, Santa Cruz, CA, USA
- Howard Hughes Medical Institute, University of California Santa Cruz, Santa Cruz, CA, USA
| |
Collapse
|
4
|
Vargas-Peralta CE, Araneda C, Galindo-Sánchez CE, Larraín MA, Del Río-Portilla MA, Lafarga-De la Cruz F. Species identification in Haliotis genus from the northeastern Pacific Ocean using genome-wide RAD-SNPs. Food Control 2022. [DOI: 10.1016/j.foodcont.2022.108979] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
|
5
|
Baratti M, Pinosio S, Gori M, Biricolti S, Chini G, Fratini S, Cannicci S, Caliani I, Oliva M, De Marchi L, Pretti C. Differential gene expression and chemical patterns of an intertidal crab inhabiting a polluted port and an adjacent marine protected area. THE SCIENCE OF THE TOTAL ENVIRONMENT 2022; 822:153463. [PMID: 35101492 DOI: 10.1016/j.scitotenv.2022.153463] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/15/2021] [Revised: 12/14/2021] [Accepted: 01/23/2022] [Indexed: 06/14/2023]
Abstract
The acquisition of data to safeguard marine protected areas located close to ports is important in order to develop plans that allow effective protection from pollution as well as sustainable development of the port. The area Secche della Meloria is a Marine Protected Area (MPA-MEL) three miles from Livorno Harbour (LH), which is characterized by a long history of pollution. Here we studied the bioaccumulation and transcriptomic patterns of the marbled crab, Pachygrapsus marmoratus (Fabricius, 1787) (Crustacea; Brachyura, Grapsidae), inhabiting the two selected sites. Results showed that the two crab populations are significantly different in their chemical composition of trace elements and Polyciclic Aromatic Hydrocarbons (PAHs), and gene expression patterns (1280 DEGs). Enrichment analysis indicated that crabs at LH had the highest stress response genes, and they were associated with higher levels of bioaccumulation detected in body tissues. We are confident that the significant differential gene expression profiles observed between crabs, characterized by significant chemical differences, is associated with responses to contaminant exposure.
Collapse
Affiliation(s)
- Mariella Baratti
- Institute of Biosciences and Bioresources, IBBR-CNR, Sesto Fiorentino I-50019, Italy.
| | - Sara Pinosio
- Institute of Biosciences and Bioresources, IBBR-CNR, Sesto Fiorentino I-50019, Italy; Institute of Applied Genomics IGA-CNR, Udine I-33100, Italy
| | - Massimo Gori
- Department of Agriculture, Food, Environment and Forestry (DAGRI), University of Florence, Firenze I-50144, Italy
| | - Stefano Biricolti
- Department of Agriculture, Food, Environment and Forestry (DAGRI), University of Florence, Firenze I-50144, Italy
| | - Giacomo Chini
- Department of Biology, University of Florence, Sesto Fiorentino I-5001, Italy
| | - Sara Fratini
- Department of Biology, University of Florence, Sesto Fiorentino I-5001, Italy
| | - Stefano Cannicci
- Department of Biology, University of Florence, Sesto Fiorentino I-5001, Italy; The Swire Institute of Marine Science and Area of Ecology and Biodiversity, School of Biological Sciences, The University of Hong Kong, Hong Kong SAR, PR China
| | - Ilaria Caliani
- Department of Physical, Earth and Environmental Sciences, University of Siena, Siena I-53100, Italy
| | - Matteo Oliva
- Interuniversity Consortium of Marine Biology "G. Bacci, CIBM ", Livorno I-57128, Italy
| | - Lucia De Marchi
- Interuniversity Consortium of Marine Biology "G. Bacci, CIBM ", Livorno I-57128, Italy; Departamento de Biologia e CESAM, Universidade de Aveiro, Aveiro 3810-193, Portugal
| | - Carlo Pretti
- Interuniversity Consortium of Marine Biology "G. Bacci, CIBM ", Livorno I-57128, Italy; Department of Veterinary Sciences, University of Pisa, Pisa I-56124, Italy
| |
Collapse
|
6
|
Sun Y, Zhang X, Wang Y, Zhang Z. Long-read RNA sequencing of Pacific abalone Haliotis discus hannai reveals innate immune system responses to environmental stress. FISH & SHELLFISH IMMUNOLOGY 2022; 122:131-145. [PMID: 35122948 DOI: 10.1016/j.fsi.2022.01.042] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/20/2021] [Revised: 01/28/2022] [Accepted: 01/31/2022] [Indexed: 06/14/2023]
Abstract
Haliotis discus hannai is a commercially important mollusk species, and the abalone aquaculture sector has been jeopardized by deteriorating environmental circumstances such as bacterial infection and thermal stress during the hot summers. However, due to a paucity of genetic information, such as transcriptome resources, our understanding of their stress adaptation is restricted. In this research, using single-molecule long-read (SMRT) sequencing technology, a library composed of ten tissues (i.e., haemocytes, gills, muscle, hepatopancreas, digestive tract, mantle, mucous gland, ovary, testis and head) was constructed and sequenced. In all, 41,855 high-quality unique transcripts, among which 24,778 were successfully annotated. Additionally, 13,463 SSRs, 1,169 transcription factors, and 18,124 lncRNAs were identified in H. discus hannai transcriptome. Furthermore, multiple immune-related transcripts were identified according to KEGG annotation, and a portion of these transcripts were mapped into several classical immune-related pathways, including the PI3K-AKT signaling pathway and Toll-like receptor signaling pathway. Additionally, 24 typical sequences related to the immunity pathway were detected by RT-PCR; the results showed that most of the immune-related genes showed significantly high expression at 72 h after bacterial challenges and thermal stress, especially the expression level of genes in gills was significantly higher than that in haemocytes under V. parahaemolyticus stress at 24 h. At the same time. The analysis of alternative splicing identified several innate immunity-related functions genes, including CD109 and caspase 2. These results suggest that the complex immune system, particularly the powerful innate immunity system, was crucial for H. discus hannai response to numerous environmental challenges.
Collapse
Affiliation(s)
- Yulong Sun
- College of Marine Science, Fujian Agriculture and Forestry University, Fuzhou, 350002, China; Fisheries College, Jimei University, Xiamen, 361021, China
| | - Xin Zhang
- College of Marine Science, Fujian Agriculture and Forestry University, Fuzhou, 350002, China; Fisheries College, Jimei University, Xiamen, 361021, China
| | - Yilei Wang
- Fisheries College, Jimei University, Xiamen, 361021, China.
| | - Ziping Zhang
- College of Marine Science, Fujian Agriculture and Forestry University, Fuzhou, 350002, China; Key Laboratory of Marine Biotechnology of Fujian Province, Institute of Oceanology, Fujian Agriculture and Forestry University, Fuzhou, 350002, China.
| |
Collapse
|
7
|
Tamagawa K, Yoshida K, Ohrui S, Takahashi Y. Population transcriptomics reveals the effect of gene flow on the evolution of range limits. Sci Rep 2022; 12:1318. [PMID: 35079049 PMCID: PMC8789792 DOI: 10.1038/s41598-022-05248-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2021] [Accepted: 01/10/2022] [Indexed: 11/17/2022] Open
Abstract
One of the most important questions in evolutionary biology is how the spatial distribution of species is limited. Asymmetric gene flow from core populations is suggested to increase the number of poorly adapted immigrants in the populations at the range edge. Genetic load due to migration, i.e., migration load, should prevent adaptation to the local habitat, leading to decreases in distribution range via local extinction or the limiting range expansion. However, few experimental studies have examined the effects of immigration on fitness and natural selection within recipient populations. To investigate the influence of migration load on the evolution of distribution range, we performed field and laboratory observations as well as population transcriptomics for the common river snail, Semisulcospira reiniana. This species meets the conditions that migration from source populations can prevent local adaptation in a sink population because they inhabit the broader range of environments, including middle/upper reaches of a river and estuaries within a single river and they may be more vulnerable to being swept away by water currents due to lowered spontaneous (upward) locomotion activity. We found that river steepness was related to the lower distribution limit of S. reiniana, with a narrower distribution range in the steeper river. Population transcriptomic analysis showed that gene flow was heavily asymmetric from the upstream populations to downstream ones in the steep river, suggesting a greater migration load in the steep river. The number of genes putatively involved in adaptation to the local habitat was lower in the steep river than in the gentle river. Gene expression profiles suggested that individuals achieve better local adaptation in the gentle river. Laboratory experiments suggested that evolutionary differences in salinity tolerance among local populations were only found in the gentle river. Our results consistent with the hypothesis that migration load owing to asymmetric gene flow disturbs local adaptation and restricts the distribution range of river snails.
Collapse
Affiliation(s)
- Katsunori Tamagawa
- Graduate School of Science, Chiba University, 1-33, Yayoi, Inage, Chiba, 263-8522, Japan.
- Graduate School of Life Sciences, Tohoku University, 6-3, Aoba, Aramaki, Sendai, Miyagi, 980-8578, Japan.
| | - Kotone Yoshida
- Graduate School of Science and Engineering, Chiba University, 1-33, Yayoi, Inage, Chiba, 263-8522, Japan
| | - Shiori Ohrui
- Graduate School of Science and Engineering, Chiba University, 1-33, Yayoi, Inage, Chiba, 263-8522, Japan
| | - Yuma Takahashi
- Graduate School of Science, Chiba University, 1-33, Yayoi, Inage, Chiba, 263-8522, Japan.
| |
Collapse
|
8
|
Hirase S, Yamasaki YY, Sekino M, Nishisako M, Ikeda M, Hara M, Merilä J, Kikuchi K. Genomic Evidence for Speciation with Gene Flow in Broadcast Spawning Marine Invertebrates. Mol Biol Evol 2021; 38:4683-4699. [PMID: 34311468 PMCID: PMC8557453 DOI: 10.1093/molbev/msab194] [Citation(s) in RCA: 15] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
How early stages of speciation in free-spawning marine invertebrates proceed is poorly understood. The Western Pacific abalones, Haliotis discus, H. madaka, and H. gigantea, occur in sympatry with shared breeding season and are capable of producing viable F1 hybrids in spite of being ecologically differentiated. Population genomic analyses revealed that although the three species are genetically distinct, there is evidence for historical and ongoing gene flow among these species. Evidence from demographic modeling suggests that reproductive isolation among the three species started to build in allopatry and has proceeded with gene flow, possibly driven by ecological selection. We identified 27 differentiation islands between the closely related H. discus and H. madaka characterized by high FST and dA, but not high dXY values, as well as high genetic diversity in one H. madaka population. These genomic signatures suggest differentiation driven by recent ecological divergent selection in presence of gene flow outside of the genomic islands of differentiation. The differentiation islands showed low polymorphism in H. gigantea, and both high FST, dXY, and dA values between H. discus and H. gigantea, as well as between H. madaka and H. gigantea. Collectively, the Western Pacific abalones appear to occupy the early stages speciation continuum, and the differentiation islands associated with ecological divergence among the abalones do not appear to have acted as barrier loci to gene flow in the younger divergences but appear to do so in older divergences.
Collapse
Affiliation(s)
- Shotaro Hirase
- Fisheries Laboratory, Graduate School of Agricultural and Life Sciences, The University of Tokyo, Maisaka, Hamamatsu, Shizuoka, Japan
| | - Yo Y Yamasaki
- Ecological Genetics Laboratory, Department of Genomics and Evolutionary Biology, National Institute of Genetics, Mishima, Shizuoka, Japan
| | - Masashi Sekino
- Bioinformatics and Biosciences Division, Fisheries Resources Institute, Japan Fisheries Research and Education Agency, Yokohama, Kanagawa, Japan
| | - Masato Nishisako
- Laboratory of Integrative Aquatic Biology, Graduate School of Agricultural Sciences, Tohoku University, Onagawa, Miyagi, Japan
| | - Minoru Ikeda
- Laboratory of Integrative Aquatic Biology, Graduate School of Agricultural Sciences, Tohoku University, Onagawa, Miyagi, Japan
| | - Motoyuki Hara
- Tohoku Ecosystem-Associated Marine Sciences, Tohoku University, Sendai, Miyagi, Japan
| | - Juha Merilä
- Ecological Genetics Research Unit, Organismal and Evolutionary Biology Research Programme, Faculty of Biological and Environmental Sciences, University of Helsinki, Helsinki, Finland
- Research Division of Ecology and Biodiversity, Faculty of Sciences, The University of Hong Kong, Hong Kong SAR, China
| | - Kiyoshi Kikuchi
- Fisheries Laboratory, Graduate School of Agricultural and Life Sciences, The University of Tokyo, Maisaka, Hamamatsu, Shizuoka, Japan
| |
Collapse
|
9
|
Clark RD, Aardema ML, Andolfatto P, Barber PH, Hattori A, Hoey JA, Montes HR, Pinsky ML. Genomic signatures of spatially divergent selection at clownfish range margins. Proc Biol Sci 2021; 288:20210407. [PMID: 34102891 PMCID: PMC8187997 DOI: 10.1098/rspb.2021.0407] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2021] [Accepted: 05/11/2021] [Indexed: 01/25/2023] Open
Abstract
Understanding how evolutionary forces interact to drive patterns of selection and distribute genetic variation across a species' range is of great interest in ecology and evolution, especially in an era of global change. While theory predicts how and when populations at range margins are likely to undergo local adaptation, empirical evidence testing these models remains sparse. Here, we address this knowledge gap by investigating the relationship between selection, gene flow and genetic drift in the yellowtail clownfish, Amphiprion clarkii, from the core to the northern periphery of the species range. Analyses reveal low genetic diversity at the range edge, gene flow from the core to the edge and genomic signatures of local adaptation at 56 single nucleotide polymorphisms in 25 candidate genes, most of which are significantly correlated with minimum annual sea surface temperature. Several of these candidate genes play a role in functions that are upregulated during cold stress, including protein turnover, metabolism and translation. Our results illustrate how spatially divergent selection spanning the range core to the periphery can occur despite the potential for strong genetic drift at the range edge and moderate gene flow from the core populations.
Collapse
Affiliation(s)
- René D. Clark
- Department of Ecology, Evolution and Natural Resources, Rutgers University, 14 College Farm Road, New Brunswick, NJ 08901, USA
| | - Matthew L. Aardema
- Department of Biology, Montclair State University, 1 Normal Avenue, Montclair, NJ 07043, USA
- Sackler Institute for Comparative Genomics, American Museum of Natural History, 200 Central Park West, New York, NY 10024-5102, USA
| | - Peter Andolfatto
- Department of Biological Sciences, Columbia University, New York, NY 10026, USA
| | - Paul H. Barber
- Department of Ecology and Evolutionary Biology, University of California-Los Angeles, Los Angeles, CA 90095, USA
| | - Akihisa Hattori
- Faculty of Liberal Arts and Education, Shiga University, 2-5-1 Hiratsu, Otsu, Shiga 520-0862, Japan
| | - Jennifer A. Hoey
- Department of Ecology, Evolution and Natural Resources, Rutgers University, 14 College Farm Road, New Brunswick, NJ 08901, USA
- Department of Ecology and Evolutionary Biology, University of California-Santa Cruz, 130 McAllister Way, Santa Cruz, CA 95060, USA
| | | | - Malin L. Pinsky
- Department of Ecology, Evolution and Natural Resources, Rutgers University, 14 College Farm Road, New Brunswick, NJ 08901, USA
| |
Collapse
|
10
|
Santos CA, Sonoda GG, Cortez T, Coutinho LL, Andrade SCS. Transcriptome Expression of Biomineralization Genes in Littoraria flava Gastropod in Brazilian Rocky Shore Reveals Evidence of Local Adaptation. Genome Biol Evol 2021; 13:6171147. [PMID: 33720344 PMCID: PMC8070887 DOI: 10.1093/gbe/evab050] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/04/2020] [Revised: 02/09/2021] [Accepted: 03/11/2021] [Indexed: 12/11/2022] Open
Abstract
Understanding how selection shapes population differentiation and local adaptation in marine species remains one of the greatest challenges in the field of evolutionary biology. The selection of genes in response to environment-specific factors and microenvironmental variation often results in chaotic genetic patchiness, which is commonly observed in rocky shore organisms. To identify these genes, the expression profile of the marine gastropod Littoraria flava collected from four Southeast Brazilian locations in ten rocky shore sites was analyzed. In this first L. flava transcriptome, 250,641 unigenes were generated, and 24% returned hits after functional annotation. Independent paired comparisons between 1) transects, 2) sites within transects, and 3) sites from different transects were performed for differential expression, detecting 8,622 unique differentially expressed genes. Araçá (AR) and São João (SJ) transect comparisons showed the most divergent gene products. For local adaptation, fitness-related differentially expressed genes were chosen for selection tests. Nine and 24 genes under adaptative and purifying selection, respectively, were most related to biomineralization in AR and chaperones in SJ. The biomineralization-genes perlucin and gigasin-6 were positively selected exclusively in the site toward the open ocean in AR, with sequence variants leading to pronounced protein structure changes. Despite an intense gene flow among L. flava populations due to its planktonic larva, gene expression patterns within transects may be the result of selective pressures. Our findings represent the first step in understanding how microenvironmental genetic variation is maintained in rocky shore populations and the mechanisms underlying local adaptation in marine species.
Collapse
Affiliation(s)
- Camilla A Santos
- Departamento de Genética e Biologia Evolutiva, Instituto de Biociências, Universidade de São Paulo, São Paulo, SP, Brazil
| | - Gabriel G Sonoda
- Departamento de Genética e Biologia Evolutiva, Instituto de Biociências, Universidade de São Paulo, São Paulo, SP, Brazil
| | - Thainá Cortez
- Departamento de Genética e Biologia Evolutiva, Instituto de Biociências, Universidade de São Paulo, São Paulo, SP, Brazil
| | - Luiz L Coutinho
- Departamento de Ciência Animal, Escola Superior de Agricultura Luiz de Queiroz (ESALQ), Universidade de São Paulo, Piracicaba, São Paulo, SP, Brazil
| | - Sónia C S Andrade
- Departamento de Genética e Biologia Evolutiva, Instituto de Biociências, Universidade de São Paulo, São Paulo, SP, Brazil
| |
Collapse
|
11
|
Tempestini A, Pinchuk AI, Dufresne F. Spatial genetic structure in Themisto libellula (Amphipoda: Hyperiidae) from the coastal Gulf of Alaska, Bering and Chukchi seas. Polar Biol 2020. [DOI: 10.1007/s00300-020-02745-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
|
12
|
Thorstensen MJ, Jeffrey JD, Treberg JR, Watkinson DA, Enders EC, Jeffries KM. Genomic signals found using RNA sequencing show signatures of selection and subtle population differentiation in walleye ( Sander vitreus) in a large freshwater ecosystem. Ecol Evol 2020; 10:7173-7188. [PMID: 32760520 PMCID: PMC7391302 DOI: 10.1002/ece3.6418] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2020] [Revised: 05/04/2020] [Accepted: 05/07/2020] [Indexed: 12/29/2022] Open
Abstract
RNA sequencing is an effective approach for studying aquatic species yielding both physiological and genomic data. However, its population genetic applications are not well-characterized. We investigate this possible role for RNA sequencing for population genomics in Lake Winnipeg, Manitoba, Canada, walleye (Sander vitreus). Lake Winnipeg walleye represent the largest component of the second-largest freshwater fishery in Canada. In the present study, large female walleye were sampled via nonlethal gill biopsy over two years at three spawning sites representing a latitudinal gradient in the lake. Genetic variation from sequenced mRNA was analyzed for neutral and adaptive markers to investigate population structure and possible adaptive variation. We find low population divergence (F ST = 0.0095), possible northward gene flow, and outlier loci that vary latitudinally in transcripts associated with cell membrane proteins and cytoskeletal function. These results indicate that Lake Winnipeg walleye may be effectively managed as a single demographically connected metapopulation with contributing subpopulations and suggest genomic differences possibly underlying observed phenotypic differences. Despite its high cost relative to other genotyping methods, RNA sequencing data can yield physiological in addition to genetic information discussed here. We therefore argue that it is useful for addressing diverse molecular questions in the conservation of freshwater species.
Collapse
Affiliation(s)
| | | | - Jason R. Treberg
- Department of Biological SciencesUniversity of ManitobaWinnipegMBCanada
| | | | - Eva C. Enders
- Freshwater Institute, Fisheries and Oceans CanadaWinnipegMBCanada
| | - Ken M. Jeffries
- Department of Biological SciencesUniversity of ManitobaWinnipegMBCanada
| |
Collapse
|
13
|
Zimmerman SJ, Aldridge CL, Oyler-McCance SJ. An empirical comparison of population genetic analyses using microsatellite and SNP data for a species of conservation concern. BMC Genomics 2020; 21:382. [PMID: 32487020 PMCID: PMC7268520 DOI: 10.1186/s12864-020-06783-9] [Citation(s) in RCA: 56] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/27/2019] [Accepted: 05/14/2020] [Indexed: 12/11/2022] Open
Abstract
BACKGROUND Use of genomic tools to characterize wildlife populations has increased in recent years. In the past, genetic characterization has been accomplished with more traditional genetic tools (e.g., microsatellites). The explosion of genomic methods and the subsequent creation of large SNP datasets has led to the promise of increased precision in population genetic parameter estimates and identification of demographically and evolutionarily independent groups, as well as questions about the future usefulness of the more traditional genetic tools. At present, few empirical comparisons of population genetic parameters and clustering analyses performed with microsatellites and SNPs have been conducted. RESULTS Here we used microsatellite and SNP data generated from Gunnison sage-grouse (Centrocercus minimus) samples to evaluate concordance of the results obtained from each dataset for common metrics of genetic diversity (HO, HE, FIS, AR) and differentiation (FST, GST, DJost). Additionally, we evaluated clustering of individuals using putatively neutral (SNPs and microsatellites), putatively adaptive, and a combined dataset of putatively neutral and adaptive loci. We took particular interest in the conservation implications of any differences. Generally, we found high concordance between microsatellites and SNPs for HE, FIS, AR, and all differentiation estimates. Although there was strong correlation between metrics from SNPs and microsatellites, the magnitude of the diversity and differentiation metrics were quite different in some cases. Clustering analyses also showed similar patterns, though SNP data was able to cluster individuals into more distinct groups. Importantly, clustering analyses with SNP data suggest strong demographic independence among the six distinct populations of Gunnison sage-grouse with some indication of evolutionary independence in two or three populations; a finding that was not revealed by microsatellite data. CONCLUSION We demonstrate that SNPs have three main advantages over microsatellites: more precise estimates of population-level diversity, higher power to identify groups in clustering methods, and the ability to consider local adaptation. This study adds to a growing body of work comparing the use of SNPs and microsatellites to evaluate genetic diversity and differentiation for a species of conservation concern with relatively high population structure and using the most common method of obtaining SNP genotypes for non-model organisms.
Collapse
Affiliation(s)
- Shawna J Zimmerman
- U.S. Geological Survey, Fort Collins Science Center, 2150 Centre Avenue, Bldg. C, Fort Collins, CO, 80526, USA.
- Department of Ecosystem Science and Sustainability and Natural Resource Ecology Laboratory, Colorado State University, Fort Collins, CO, 80526, USA.
| | - Cameron L Aldridge
- U.S. Geological Survey, Fort Collins Science Center, 2150 Centre Avenue, Bldg. C, Fort Collins, CO, 80526, USA
- Department of Ecosystem Science and Sustainability and Natural Resource Ecology Laboratory, Colorado State University, Fort Collins, CO, 80526, USA
| | - Sara J Oyler-McCance
- U.S. Geological Survey, Fort Collins Science Center, 2150 Centre Avenue, Bldg. C, Fort Collins, CO, 80526, USA
| |
Collapse
|
14
|
Klein AH, Ballard KR, Storey KB, Motti CA, Zhao M, Cummins SF. Multi-omics investigations within the Phylum Mollusca, Class Gastropoda: from ecological application to breakthrough phylogenomic studies. Brief Funct Genomics 2020; 18:377-394. [PMID: 31609407 DOI: 10.1093/bfgp/elz017] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2019] [Revised: 07/06/2019] [Accepted: 07/15/2019] [Indexed: 12/22/2022] Open
Abstract
Gastropods are the largest and most diverse class of mollusc and include species that are well studied within the areas of taxonomy, aquaculture, biomineralization, ecology, microbiome and health. Gastropod research has been expanding since the mid-2000s, largely due to large-scale data integration from next-generation sequencing and mass spectrometry in which transcripts, proteins and metabolites can be readily explored systematically. Correspondingly, the huge data added a great deal of complexity for data organization, visualization and interpretation. Here, we reviewed the recent advances involving gastropod omics ('gastropodomics') research from hundreds of publications and online genomics databases. By summarizing the current publicly available data, we present an insight for the design of useful data integrating tools and strategies for comparative omics studies in the future. Additionally, we discuss the future of omics applications in aquaculture, natural pharmaceutical biodiscovery and pest management, as well as to monitor the impact of environmental stressors.
Collapse
Affiliation(s)
- Anne H Klein
- Genecology Research Centre, University of the Sunshine Coast, Maroochydore DC, Queensland 4558, Australia
| | - Kaylene R Ballard
- Genecology Research Centre, University of the Sunshine Coast, Maroochydore DC, Queensland 4558, Australia
| | - Kenneth B Storey
- Institute of Biochemistry & Department of Biology, Carleton University, Ottawa, ON, Canada K1S 5B6
| | - Cherie A Motti
- Australian Institute of Marine Science (AIMS), Cape Ferguson, Townsville Queensland 4810, Australia
| | - Min Zhao
- Genecology Research Centre, University of the Sunshine Coast, Maroochydore DC, Queensland 4558, Australia
| | - Scott F Cummins
- Genecology Research Centre, University of the Sunshine Coast, Maroochydore DC, Queensland 4558, Australia
| |
Collapse
|
15
|
Freedman AH, Clamp M, Sackton TB. Error, noise and bias in de novo transcriptome assemblies. Mol Ecol Resour 2020; 21:18-29. [PMID: 32180366 DOI: 10.1111/1755-0998.13156] [Citation(s) in RCA: 30] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2019] [Revised: 01/25/2020] [Accepted: 03/10/2020] [Indexed: 12/21/2022]
Abstract
De novo transcriptome assembly is a powerful tool, and has been widely used over the last decade for making evolutionary inferences. However, it relies on two implicit assumptions: that the assembled transcriptome is an unbiased representation of the underlying expressed transcriptome, and that expression estimates from the assembly are good, if noisy approximations of the relative abundance of expressed transcripts. Using publicly available data for model organisms, we demonstrate that, across assembly algorithms and data sets, these assumptions are consistently violated. Bias exists at the nucleotide level, with genotyping error rates ranging from 30% to 83%. As a result, diversity is underestimated in transcriptome assemblies, with consistent underestimation of heterozygosity in all but the most inbred samples. Even at the gene level, expression estimates show wide deviations from map-to-reference estimates, and positive bias at lower expression levels. Standard filtering of transcriptome assemblies improves the robustness of gene expression estimates but leads to the loss of a meaningful number of protein-coding genes, including many that are highly expressed. We demonstrate a computational method, length-rescaled CPM, to partly alleviate noise and bias in expression estimates. Researchers should consider ways to minimize the impact of bias in transcriptome assemblies.
Collapse
Affiliation(s)
- Adam H Freedman
- Faculty of Arts and Sciences Informatics Group, Harvard University, Cambridge, MA, USA
| | - Michele Clamp
- Faculty of Arts and Sciences Informatics Group, Harvard University, Cambridge, MA, USA
| | - Timothy B Sackton
- Faculty of Arts and Sciences Informatics Group, Harvard University, Cambridge, MA, USA
| |
Collapse
|
16
|
Chandra Rajan K, Vengatesen T. Molecular adaptation of molluscan biomineralisation to high-CO 2 oceans - The known and the unknown. MARINE ENVIRONMENTAL RESEARCH 2020; 155:104883. [PMID: 32072987 DOI: 10.1016/j.marenvres.2020.104883] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/19/2019] [Revised: 01/11/2020] [Accepted: 01/19/2020] [Indexed: 06/10/2023]
Abstract
High-CO2 induced ocean acidification (OA) reduces the calcium carbonate (CaCO3) saturation level (Ω) and the pH of oceans. Consequently, OA is causing a serious threat to several ecologically and economically important biomineralising molluscs. Biomineralisation is a highly controlled biochemical process by which molluscs deposit their calcareous structures. In this process, shell matrix proteins aid the nucleation, growth and assemblage of the CaCO3 crystals in the shell. These molluscan shell proteins (MSPs) are, ultimately, responsible for determination of the diverse shell microstructures and mechanical strength. Recent studies have attempted to integrate gene and protein expression data of MSPs with shell structure and mechanical properties. These advances made in understanding the molecular mechanism of biomineralisation suggest that molluscs either succumb or adapt to OA stress. In this review, we discuss the fate of biomineralisation process in future high-CO2 oceans and its ultimate impact on the mineralised shell's structure and mechanical properties from the perspectives of limited substrate availability theory, proton flux limitation model and the omega myth theory. Furthermore, studying the interplay of energy availability and differential gene expression is an essential first step towards understanding adaptation of molluscan biomineralisation to OA, because if there is a need to change gene expression under stressors, any living system would require more energy than usual. To conclude, we have listed, four important future research directions for molecular adaptation of molluscan biomineralisation in high-CO2 oceans: 1) Including an energy budgeting factor while understanding differential gene expression of MSPs and ion transporters under OA. 2) Unraveling the genetic or epigenetic changes related to biomineralisation under stressors to help solving a bigger picture about future evolution of molluscs, and 3) Understanding Post Translational Modifications of MSPs with and without stressors. 4) Understanding carbon uptake mechanisms across taxa with and without OA to clarify the OA theories on Ω.
Collapse
Affiliation(s)
- Kanmani Chandra Rajan
- The Swire Institute of Marine Science and School of Biological Sciences, The University of Hong Kong, Hong Kong SAR, China; State Key Laboratory of Marine Pollution, Hong Kong SAR, China.
| | - Thiyagarajan Vengatesen
- The Swire Institute of Marine Science and School of Biological Sciences, The University of Hong Kong, Hong Kong SAR, China; State Key Laboratory of Marine Pollution, Hong Kong SAR, China.
| |
Collapse
|
17
|
Bashevkin SM, Dibble CD, Dunn RP, Hollarsmith JA, Ng G, Satterthwaite EV, Morgan SG. Larval dispersal in a changing ocean with an emphasis on upwelling regions. Ecosphere 2020. [DOI: 10.1002/ecs2.3015] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/03/2023] Open
Affiliation(s)
- Samuel M. Bashevkin
- Department of Environmental Science and Policy University of California, Davis One Shields Avenue Davis California 95616 USA
- Coastal and Marine Sciences Institute and Bodega Marine Laboratory University of California, Davis 2099 Westshore Road, P.O. Box 247 Bodega Bay California 94923 USA
- Delta Science Program Delta Stewardship Council 980 9th Street, Suite 1500 Sacramento California 95814 USA
| | - Connor D. Dibble
- Department of Environmental Science and Policy University of California, Davis One Shields Avenue Davis California 95616 USA
- Coastal and Marine Sciences Institute and Bodega Marine Laboratory University of California, Davis 2099 Westshore Road, P.O. Box 247 Bodega Bay California 94923 USA
| | - Robert P. Dunn
- Department of Environmental Science and Policy University of California, Davis One Shields Avenue Davis California 95616 USA
- Coastal and Marine Institute & Department of Biology San Diego State University 4165 Spruance Road San Diego California 92182 USA
| | - Jordan A. Hollarsmith
- Department of Environmental Science and Policy University of California, Davis One Shields Avenue Davis California 95616 USA
- Coastal and Marine Sciences Institute and Bodega Marine Laboratory University of California, Davis 2099 Westshore Road, P.O. Box 247 Bodega Bay California 94923 USA
- Department of Biological Sciences Simon Fraser University 8888 University Drive Burnaby British Columbia V5A 1S6 Canada
| | - Gabriel Ng
- Coastal and Marine Sciences Institute and Bodega Marine Laboratory University of California, Davis 2099 Westshore Road, P.O. Box 247 Bodega Bay California 94923 USA
- Department of Evolution and Ecology University of California, Davis One Shields Avenue Davis California 95616 USA
| | - Erin V. Satterthwaite
- Department of Environmental Science and Policy University of California, Davis One Shields Avenue Davis California 95616 USA
- Coastal and Marine Sciences Institute and Bodega Marine Laboratory University of California, Davis 2099 Westshore Road, P.O. Box 247 Bodega Bay California 94923 USA
- National Center for Ecological Analysis and Synthesis University of California Santa Barbara 735 State Street, Suite 300 Santa Barbara California USA
- Future Earth School of Global Environmental Sustainability Colorado State University 108 Johnson Drive Fort Collins Colorado 80523 USA
| | - Steven G. Morgan
- Department of Environmental Science and Policy University of California, Davis One Shields Avenue Davis California 95616 USA
- Coastal and Marine Sciences Institute and Bodega Marine Laboratory University of California, Davis 2099 Westshore Road, P.O. Box 247 Bodega Bay California 94923 USA
| |
Collapse
|
18
|
Rahman S, Schmidt D, Hughes JM. Genetic structure of Australian glass shrimp, Paratya australiensis, in relation to altitude. PeerJ 2020; 8:e8139. [PMID: 31942250 PMCID: PMC6955102 DOI: 10.7717/peerj.8139] [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: 10/26/2018] [Accepted: 11/01/2019] [Indexed: 11/29/2022] Open
Abstract
Paratya australiensis Kemp (Decapoda: Atyidae) is a widely distributed freshwater shrimp in eastern Australia. The species has been considered as an important stream organism for studying genetics, dispersal, biology, behaviour and evolution in atyids and is a major food source for stream dwelling fishes. Paratya australiensis is a cryptic species complex consisting of nine highly divergent mitochondrial DNA lineages. Previous studies in southeast Queensland showed that “lineage 4” favours upstream sites at higher altitudes, with cooler water temperatures. This study aims to identify putative selection and population structure between high elevation and low elevation populations of this lineage at relatively small spatial scales. Sample localities were selected from three streams: Booloumba Creek, Broken Bridge Creek and Obi Obi Creek in the Conondale Range, southeast Queensland. Six sample localities, consisting of 142 individuals in total were sequenced using double digest Restriction Site Associated DNA-sequencing (ddRAD-seq) technique. Among the 142 individuals, 131 individuals shared 213 loci. Outlier analysis on 213 loci showed that 27 loci were putatively under selection between high elevation and low elevation populations. Outlier analysis on individual streams was also done to test for parallel patterns of adaptation, but there was no evidence of a parallel pattern. Population structure was observed using both the 27 outliers and 186 neutral loci and revealed similar population structure in both cases. Therefore, we cannot differentiate between selection and drift here. The highest genetic differentiation was observed between high elevation and low elevation populations of Booloumba Creek, with small levels of differentiation in the other two streams.
Collapse
Affiliation(s)
- Sharmeen Rahman
- Griffith School of Environment and Australian Rivers Institute, Griffith University, Brisbane, QLD, Australia
| | - Daniel Schmidt
- Griffith School of Environment and Australian Rivers Institute, Griffith University, Brisbane, QLD, Australia
| | - Jane M Hughes
- Griffith School of Environment and Australian Rivers Institute, Griffith University, Brisbane, QLD, Australia
| |
Collapse
|
19
|
Sun Y, Zhang X, Wang Y, Day R, Yang H, Zhang Z. Immunity-related genes and signaling pathways under hypoxic stresses in Haliotis diversicolor: a transcriptome analysis. Sci Rep 2019; 9:19741. [PMID: 31874975 PMCID: PMC6930256 DOI: 10.1038/s41598-019-56150-2] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/14/2019] [Accepted: 12/03/2019] [Indexed: 12/30/2022] Open
Abstract
Due to increased temperatures and aquaculture density, thermal and hypoxia stresses have become serious problems for the aquaculture of abalone Haliotis diversicolor. Stresses lead to immunosuppression, which can cause severe negative impacts on aquaculture farms. To study the mechanism of immunosuppression after hypoxia stress and bacterial challenge, transcriptomes of H. diversicolor hemocytes involved in immunity were profiled. A total of 307,395,572 clean reads were generated and assembled into 99,774 unigenes. KEGG analysis indicated that 225 unigenes with immunologic function were mapped into immune-related pathways. Expression of 41 unigenes measured by quantitative real-time PCR (qRT-PCR) showed consistent results with that of transcriptome analysis. When exposure challenge of Vibrio parahaemolyticus, it is indicated that the PI3K-AKT, MAPK, NF-κB and P53 signal pathways were involved in the hypoxia-induced immunosuppression of H. diversicolor. Furthermore, when the AKT gene (HdAKT) was inhibited by double-stranded RNA (dsRNA), expression levels of HdAKT was lower than the blank and control group in hemocytes at 4 h, 12 h and 24 h (p < 0.05).
Collapse
Affiliation(s)
- Yulong Sun
- College of Animal Science, Fujian Agriculture and Forestry University, Fuzhou, 350002, P.R. China
- Key Laboratory of Marine Biotechnology of Fujian Province, Institute of Oceanology, Fujian Agriculture and Forestry University, Fuzhou, 350002, P.R. China
| | - Xin Zhang
- College of Animal Science, Fujian Agriculture and Forestry University, Fuzhou, 350002, P.R. China
- Key Laboratory of Marine Biotechnology of Fujian Province, Institute of Oceanology, Fujian Agriculture and Forestry University, Fuzhou, 350002, P.R. China
| | - Yilei Wang
- Fisheries College, Jimei University, Xiamen, 361021, China.
| | - Robert Day
- School of Biosciences, University of Melbourne, Parkville, Victoria, 3010, Australia
| | - Huiping Yang
- School of Forest Resources and Conservation, IFAS, University of Florida, 7922 NW 71st Street, Gainesville, FL, 32615, USA
| | - Ziping Zhang
- College of Animal Science, Fujian Agriculture and Forestry University, Fuzhou, 350002, P.R. China.
- Fisheries College, Jimei University, Xiamen, 361021, China.
| |
Collapse
|
20
|
Zimmerman SJ, Aldridge CL, Oh KP, Cornman RS, Oyler‐McCance SJ. Signatures of adaptive divergence among populations of an avian species of conservation concern. Evol Appl 2019; 12:1661-1677. [PMID: 31462921 PMCID: PMC6708427 DOI: 10.1111/eva.12825] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/20/2018] [Revised: 05/23/2019] [Accepted: 05/24/2019] [Indexed: 12/16/2022] Open
Abstract
Understanding the genetic underpinning of adaptive divergence among populations is a key goal of evolutionary biology and conservation. Gunnison sage-grouse (Centrocercus minimus) is a sagebrush obligate species with a constricted range consisting of seven discrete populations, each with distinctly different habitat and climatic conditions. Though geographically close, populations have low levels of natural gene flow resulting in relatively high levels of differentiation. Here, we use 15,033 SNP loci in genomic outlier analyses, genotype-environment association analyses, and gene ontology enrichment tests to examine patterns of putatively adaptive genetic differentiation in an avian species of conservation concern. We found 411 loci within 5 kbp of 289 putative genes associated with biological functions or pathways that were overrepresented in the assemblage of outlier SNPs. The identified gene set was enriched for cytochrome P450 gene family members (CYP4V2, CYP2R1, CYP2C23B, CYP4B1) and could impact metabolism of plant secondary metabolites, a critical challenge for sagebrush obligates. Additionally, the gene set was also enriched with members potentially involved in antiviral response (DEAD box helicase gene family and SETX). Our results provide a first look at local adaption for isolated populations of a single species and suggest adaptive divergence in multiple metabolic and biochemical pathways may be occurring. This information can be useful in managing this species of conservation concern, for example, to identify unique populations to conserve, avoid translocation or release of individuals that may swamp locally adapted genetic diversity, or guide habitat restoration efforts.
Collapse
Affiliation(s)
- Shawna J. Zimmerman
- Department of Ecosystem Science and Sustainability and Natural Resource Ecology Laboratory, Colorado State University in Cooperation with U.S. Geological SurveyFort Collins Science CenterFort CollinsColorado
| | - Cameron L. Aldridge
- Department of Ecosystem Science and Sustainability and Natural Resource Ecology Laboratory, Colorado State University in Cooperation with U.S. Geological SurveyFort Collins Science CenterFort CollinsColorado
| | - Kevin P. Oh
- U.S. Geological SurveyFort Collins Science CenterFort CollinsColorado
| | - Robert S. Cornman
- U.S. Geological SurveyFort Collins Science CenterFort CollinsColorado
| | | |
Collapse
|
21
|
Wang Y, Shahid MQ, Ghouri F, Ercişli S, Baloch FS. Development of EST-based SSR and SNP markers in Gastrodia elata (herbal medicine) by sequencing, de novo assembly and annotation of the transcriptome. 3 Biotech 2019; 9:292. [PMID: 31321198 DOI: 10.1007/s13205-019-1823-4] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2018] [Accepted: 06/23/2019] [Indexed: 01/28/2023] Open
Abstract
Tianma (Gastrodia elata Blume) has unique biological characteristics and high medicinal value. The wild resource of G. elata is being overutilized and should be conserved as it is already included in the list of endangered species in China. The population size of cultivated G. elata is small because of domestication bottleneck. Therefore, it is of utmost importance to evolve high-quality varieties and conserve wild resources of G. elata. In this study, we sequenced tuber transcriptomes of three major cultivated sub-species of Gastrodia elata, namely G. elata BI. f. elata, G. elata Bl. f. glauca S. Chow, and G. elata Bl. f. Viridis, and obtained about 7.8G clean data. The assembled high-quality reads of three sub-species were clustered into 56,884 unigenes. Of these, 31,224 (54.89%), 25,733 (45.24%), 22,629 (39.78%), and 11,856 (20.84%) unigenes were annotated by Nr, Swiss-Port, Eukaryotic Ortholog Groups (KOG), and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases, respectively. Here, a total of 3766 EST-SSRs and 128,921 SNPs were identified from the unigenes. The results not only offer huge number of genes that were responsible for the growth, development, and metabolism of bioactive components, but also a large number of molecular markers were detected for future studies on the conservation genetics and molecular breeding of G. elata.
Collapse
|
22
|
Wong JM, Gaitán-Espitia JD, Hofmann GE. Transcriptional profiles of early stage red sea urchins (Mesocentrotus franciscanus) reveal differential regulation of gene expression across development. Mar Genomics 2019; 48:100692. [PMID: 31227413 DOI: 10.1016/j.margen.2019.05.007] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2019] [Revised: 05/28/2019] [Accepted: 05/29/2019] [Indexed: 12/17/2022]
Abstract
The red sea urchin, Mesocentrotus franciscanus, is an ecologically important kelp forest species that also serves as a valuable fisheries resource. In this study, we have assembled and annotated a developmental transcriptome for M. franciscanus that represents eggs and six stages of early development (8- to 16-cell, morula, hatched blastula, early gastrula, prism and early pluteus). Characterization of the transcriptome revealed distinct patterns of gene expression that corresponded to major developmental and morphological processes. In addition, the period during which maternally-controlled transcription was terminated and the zygotic genome was activated, the maternal-to-zygotic transition (MZT), was found to begin during early cleavage and persist through the hatched blastula stage, an observation that is similar to the timing of the MZT in other sea urchin species. The presented developmental transcriptome will serve as a useful resource for investigating, in both an ecological and fisheries context, how the early developmental stages of this species respond to environmental stressors.
Collapse
Affiliation(s)
- Juliet M Wong
- Department of Ecology, Evolution and Marine Biology, University of California Santa Barbara, Santa Barbara, CA 93106, USA.
| | - Juan D Gaitán-Espitia
- The Swire Institute of Marine Science, School of Biological Sciences, The University of Hong Kong, Pokfulam Road, Hong Kong Special Administrative Region.
| | - Gretchen E Hofmann
- Department of Ecology, Evolution and Marine Biology, University of California Santa Barbara, Santa Barbara, CA 93106, USA.
| |
Collapse
|
23
|
Miller AD, Hoffmann AA, Tan MH, Young M, Ahrens C, Cocomazzo M, Rattray A, Ierodiaconou DA, Treml E, Sherman CDH. Local and regional scale habitat heterogeneity contribute to genetic adaptation in a commercially important marine mollusc (
Haliotis rubra
) from southeastern Australia. Mol Ecol 2019; 28:3053-3072. [DOI: 10.1111/mec.15128] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/16/2018] [Revised: 04/17/2019] [Accepted: 05/01/2019] [Indexed: 12/16/2022]
Affiliation(s)
- Adam D. Miller
- School of Life and Environmental Sciences Centre for Integrative Ecology, Deakin University Geelong Victoria Australia
- Deakin Genomics Centre Deakin University Geelong Victoria Australia
| | - Ary A. Hoffmann
- School of BioSciences Bio21 Institute, The University of Melbourne Parkville Victoria Australia
| | - Mun Hua Tan
- School of Life and Environmental Sciences Centre for Integrative Ecology, Deakin University Geelong Victoria Australia
- Deakin Genomics Centre Deakin University Geelong Victoria Australia
| | - Mary Young
- School of Life and Environmental Sciences Centre for Integrative Ecology, Deakin University Geelong Victoria Australia
| | - Collin Ahrens
- Hawkesbury Institute for the Environment, Western Sydney University Penrith New South Wales Australia
| | - Michael Cocomazzo
- School of Life and Environmental Sciences Centre for Integrative Ecology, Deakin University Geelong Victoria Australia
| | - Alex Rattray
- School of Life and Environmental Sciences Centre for Integrative Ecology, Deakin University Geelong Victoria Australia
| | - Daniel A. Ierodiaconou
- School of Life and Environmental Sciences Centre for Integrative Ecology, Deakin University Geelong Victoria Australia
| | - Eric Treml
- School of Life and Environmental Sciences Centre for Integrative Ecology, Deakin University Geelong Victoria Australia
| | - Craig D. H. Sherman
- School of Life and Environmental Sciences Centre for Integrative Ecology, Deakin University Geelong Victoria Australia
- Deakin Genomics Centre Deakin University Geelong Victoria Australia
| |
Collapse
|
24
|
Silliman K. Population structure, genetic connectivity, and adaptation in the Olympia oyster ( Ostrea lurida) along the west coast of North America. Evol Appl 2019; 12:923-939. [PMID: 31080505 PMCID: PMC6503834 DOI: 10.1111/eva.12766] [Citation(s) in RCA: 29] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2018] [Revised: 11/28/2018] [Accepted: 12/02/2018] [Indexed: 01/02/2023] Open
Abstract
Effective management of threatened and exploited species requires an understanding of both the genetic connectivity among populations and local adaptation. The Olympia oyster (Ostrea lurida), patchily distributed from Baja California to the central coast of Canada, has a long history of population declines due to anthropogenic stressors. For such coastal marine species, population structure could follow a continuous isolation-by-distance model, contain regional blocks of genetic similarity separated by barriers to gene flow, or be consistent with a null model of no population structure. To distinguish between these hypotheses in O. lurida, 13,424 single nucleotide polymorphisms (SNPs) were used to characterize rangewide population structure, genetic connectivity, and adaptive divergence. Samples were collected across the species range on the west coast of North America, from southern California to Vancouver Island. A conservative approach for detecting putative loci under selection identified 235 SNPs across 129 GBS loci, which were functionally annotated and analyzed separately from the remaining neutral loci. While strong population structure was observed on a regional scale in both neutral and outlier markers, neutral markers had greater power to detect fine-scale structure. Geographic regions of reduced gene flow aligned with known marine biogeographic barriers, such as Cape Mendocino, Monterey Bay, and the currents around Cape Flattery. The outlier loci identified as under putative selection included genes involved in developmental regulation, sensory information processing, energy metabolism, immune response, and muscle contraction. These loci are excellent candidates for future research and may provide targets for genetic monitoring programs. Beyond specific applications for restoration and management of the Olympia oyster, this study lends to the growing body of evidence for both population structure and adaptive differentiation across a range of marine species exhibiting the potential for panmixia. Computational notebooks are available to facilitate reproducibility and future open-sourced research on the population structure of O. lurida.
Collapse
|
25
|
Effects of Predator-Prey Interactions on Predator Traits: Differentiation of Diets and Venoms of a Marine Snail. Toxins (Basel) 2019; 11:toxins11050299. [PMID: 31130611 PMCID: PMC6563511 DOI: 10.3390/toxins11050299] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/24/2019] [Revised: 05/16/2019] [Accepted: 05/23/2019] [Indexed: 12/12/2022] Open
Abstract
Species interactions are fundamental ecological forces that can have significant impacts on the evolutionary trajectories of species. Nonetheless, the contribution of predator-prey interactions to genetic and phenotypic divergence remains largely unknown. Predatory marine snails of the family Conidae exhibit specializations for different prey items and intraspecific variation in prey utilization patterns at geographic scales. Because cone snails utilize venom to capture prey and venom peptides are direct gene products, it is feasible to examine the evolution of genes associated with changes in resource utilization. Here, we compared feeding ecologies and venom duct transcriptomes of individuals from three populations of Conus miliaris, a species that exhibits geographic variation in prey utilization and dietary breadth, in order to determine the extent to which dietary differences are correlated with differences in venom composition, and if expanded niche breadth is associated with increased variation in venom composition. While populations showed little to no overlap in resource utilization, taxonomic richness of prey was greatest at Easter Island. Changes in dietary breadth were associated with differences in expression patterns and increased genetic differentiation of toxin-related genes. The Easter Island population also exhibited greater diversity of toxin-related transcripts, but did not show increased variance in expression of these transcripts. These results imply that differences in dietary breadth contribute more to the structural and regulatory differentiation of venoms than differences in diet.
Collapse
|
26
|
Griffiths JS, Pan TCF, Kelly MW. Differential responses to ocean acidification between populations of Balanophyllia elegans corals from high and low upwelling environments. Mol Ecol 2019; 28:2715-2730. [PMID: 30770604 DOI: 10.1111/mec.15050] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/28/2018] [Revised: 02/07/2019] [Accepted: 02/08/2019] [Indexed: 12/22/2022]
Abstract
Ocean acidification (OA), the global decrease in surface water pH from absorption of anthropogenic CO2 , may put many marine taxa at risk. However, populations that experience extreme localized conditions, and are adapted to these conditions predicted in the global ocean in 2,100, may be more tolerant to future OA. By identifying locally adapted populations, researchers can examine the mechanisms used to cope with decreasing pH. One oceanographic process that influences pH is wind-driven upwelling. Here we compare two Californian populations of the coral Balanophyllia elegans from distinct upwelling regimes, and test their physiological and transcriptomic responses to experimental seawater acidification. We measured respiration rates, protein and lipid content, and gene expression in corals from both populations exposed to pH levels of 7.8 and 7.4 for 29 days. Corals from the population that experiences lower pH due to high upwelling maintained the same respiration rate throughout the exposure. In contrast, corals from the low upwelling site had reduced respiration rates, protein content and lipid-class content at low pH exposure, suggesting they have depleted their energy reserves. Using RNA-Seq, we found that corals from the high upwelling site upregulated genes involved in calcium ion binding and ion transport, most likely related to pH homeostasis and calcification. In contrast, corals from the low upwelling site downregulated stress response genes at low pH exposure. Divergent population responses to low pH observed in B. elegans highlight the importance of multi-population studies for predicting a species' response to future OA.
Collapse
Affiliation(s)
- Joanna S Griffiths
- Department of Biological Sciences, Louisiana State University, Baton Rouge, Louisiana
| | - Tien-Chien Francis Pan
- Department of Biological Sciences, University of Southern California, Los Angeles, California
| | - Morgan W Kelly
- Department of Biological Sciences, Louisiana State University, Baton Rouge, Louisiana
| |
Collapse
|
27
|
Rellstab C, Dauphin B, Zoller S, Brodbeck S, Gugerli F. Using transcriptome sequencing and pooled exome capture to study local adaptation in the giga‐genome of
Pinus cembra. Mol Ecol Resour 2019; 19:536-551. [DOI: 10.1111/1755-0998.12986] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/02/2018] [Revised: 11/29/2018] [Accepted: 12/05/2018] [Indexed: 12/21/2022]
Affiliation(s)
| | | | - Stefan Zoller
- ETH Zürich Genetic Diversity Centre Zürich Switzerland
| | - Sabine Brodbeck
- WSL Swiss Federal Research Institute Birmensdorf Switzerland
| | - Felix Gugerli
- WSL Swiss Federal Research Institute Birmensdorf Switzerland
| |
Collapse
|
28
|
Chen N, Huang Z, Lu C, Shen Y, Luo X, Ke C, You W. Different Transcriptomic Responses to Thermal Stress in Heat-Tolerant and Heat-Sensitive Pacific Abalones Indicated by Cardiac Performance. Front Physiol 2019; 9:1895. [PMID: 30687115 PMCID: PMC6334008 DOI: 10.3389/fphys.2018.01895] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2018] [Accepted: 12/17/2018] [Indexed: 12/12/2022] Open
Abstract
The Pacific abalone Haliotis discus hannai is one of the most economically important mollusks in China. Even though it has been farmed in southern China for almost 20 years, summer mortality remains the most challengeable problem for Pacific abalone aquaculture recently. Here, we determined the different heat tolerance ability for five selective lines of H. discus hannai by measuring the cardiac performance and Arrhenius breakpoint temperature (ABT). The Red line (RL) and Yangxia line (YL) were determined as the most heat-sensitive and most heat-tolerant line, respectively. Heart rates for RL were significantly lower than those of the YL at the same temperature (p < 0.05). The differentially expressed genes (DEGs), which were enriched in several pathways including cardiac muscle contraction, glutathione metabolism and oxidative phosphorylation, were identified between RL and YL at control temperature (20°C) and heat stress temperature (28.5°C, the ABT of the RL) by RNA-seq method. In the RL, 3370 DEGs were identified between the control and the heat-stress temperature, while only 1351 DEGs were identified in YL between these two temperature tests. Most of these DEGs were enriched in the pathways such as protein processing in endoplasmic reticulum, nucleotide binding and oligomerization domain (NOD) like receptor signaling, and ubiquitin mediated proteolysis. Notably, the most heat-tolerant line YL used an effective heat-protection strategy based on moderate transcriptional changes and regulation on the expression of key genes.
Collapse
Affiliation(s)
- Nan Chen
- State Key Laboratory of Marine Environmental Science, Xiamen University, Xiamen, China
- College of the Environment and Ecology, Xiamen University, Xiamen, China
| | - Zekun Huang
- State Key Laboratory of Marine Environmental Science, Xiamen University, Xiamen, China
- College of Ocean and Earth Sciences, Xiamen University, Xiamen, China
- State-Province Joint Engineering Laboratory of Marine Bioproducts and Technology, College of Ocean and Earth Sciences, Xiamen University, Xiamen, China
| | - Chengkuan Lu
- College of Ocean and Earth Sciences, Xiamen University, Xiamen, China
- State-Province Joint Engineering Laboratory of Marine Bioproducts and Technology, College of Ocean and Earth Sciences, Xiamen University, Xiamen, China
| | - Yawei Shen
- State Key Laboratory of Marine Environmental Science, Xiamen University, Xiamen, China
- College of Ocean and Earth Sciences, Xiamen University, Xiamen, China
- State-Province Joint Engineering Laboratory of Marine Bioproducts and Technology, College of Ocean and Earth Sciences, Xiamen University, Xiamen, China
| | - Xuan Luo
- State Key Laboratory of Marine Environmental Science, Xiamen University, Xiamen, China
- College of Ocean and Earth Sciences, Xiamen University, Xiamen, China
- State-Province Joint Engineering Laboratory of Marine Bioproducts and Technology, College of Ocean and Earth Sciences, Xiamen University, Xiamen, China
| | - Caihuan Ke
- State Key Laboratory of Marine Environmental Science, Xiamen University, Xiamen, China
- College of Ocean and Earth Sciences, Xiamen University, Xiamen, China
- State-Province Joint Engineering Laboratory of Marine Bioproducts and Technology, College of Ocean and Earth Sciences, Xiamen University, Xiamen, China
| | - Weiwei You
- State Key Laboratory of Marine Environmental Science, Xiamen University, Xiamen, China
- College of Ocean and Earth Sciences, Xiamen University, Xiamen, China
- State-Province Joint Engineering Laboratory of Marine Bioproducts and Technology, College of Ocean and Earth Sciences, Xiamen University, Xiamen, China
| |
Collapse
|
29
|
Xuereb A, D’Aloia CC, Daigle RM, Andrello M, Dalongeville A, Manel S, Mouillot D, Guichard F, Côté IM, Curtis JMR, Bernatchez L, Fortin MJ. Marine Conservation and Marine Protected Areas. POPULATION GENOMICS 2019. [DOI: 10.1007/13836_2018_63] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
|
30
|
Zhang X, Shi J, Sun Y, Habib YJ, Yang H, Zhang Z, Wang Y. Integrative transcriptome analysis and discovery of genes involving in immune response of hypoxia/thermal challenges in the small abalone Haliotis diversicolor. FISH & SHELLFISH IMMUNOLOGY 2019; 84:609-626. [PMID: 30366091 DOI: 10.1016/j.fsi.2018.10.044] [Citation(s) in RCA: 20] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/25/2018] [Revised: 10/18/2018] [Accepted: 10/20/2018] [Indexed: 06/08/2023]
Abstract
In recent years, the abalone aquaculture industry has been threatened by the deteriorating environmental conditions, such as hypoxia and thermal stress in the hot summers. It is necessary to investigate the molecular mechanism in response to these environmental challenges, and subsequently understand the immune defense system. In this study, the transcriptome profiles by RNA-seq of hemocytes from the small abalone Haliotis diversicolor after exposure to hypoxia, thermal stress, and hypoxia plus thermal stress were established. A total of 103,703,074 clean reads were obtained and 99,774 unigenes were assembled. Of the 99,774 unigenes, 47,154 and 20,455 had homologous sequences in the Nr and Swiss-Prot protein databases, while 16,944 and 10,840 unigenes could be classified by COG or KEGG databases, respectively. RNAseq analysis revealed that the differentially expressed genes (DEGs) after challenges of hypoxia, thermal stress, or hypoxia plus thermal stress were 24,189, 29,165 and 23,665, among which more than 3000 genes involved in at least 230 pathways, including several classical immune-related pathways. The genes and pathways that were involved in immune response to hypoxia/thermal challenges were identified by transcriptome analysis and further validated by quantitative real-time PCR and RNAi technology. The findings in this study can provide information on H. diversicolor innate immunity to improve the abalone aquaculture industry, and the analysis of the potential immune-related genes in innate immunity signaling pathways and the obtained transcriptome data can provide an invaluable genetic resource for the study of the genome and functional genes.
Collapse
Affiliation(s)
- Xin Zhang
- College of Animal Sciences, Fujian Agriculture and Forestry University, Fuzhou, Fujian Province, 350002, China
| | - Jialong Shi
- College of Animal Sciences, Fujian Agriculture and Forestry University, Fuzhou, Fujian Province, 350002, China
| | - Yulong Sun
- College of Animal Sciences, Fujian Agriculture and Forestry University, Fuzhou, Fujian Province, 350002, China
| | - Yusuf Jibril Habib
- College of Animal Sciences, Fujian Agriculture and Forestry University, Fuzhou, Fujian Province, 350002, China
| | - Huiping Yang
- School of Forest Resources and Conservation, Institute of Food and Agricultural Sciences, University of Florida, 7922 NW 71st Street, Gainesville, FL, 32653, USA
| | - Ziping Zhang
- College of Animal Sciences, Fujian Agriculture and Forestry University, Fuzhou, Fujian Province, 350002, China.
| | - Yilei Wang
- Fisheries College, Jimei University, Xiamen, 361021, China.
| |
Collapse
|
31
|
Xuereb A, Kimber CM, Curtis JMR, Bernatchez L, Fortin M. Putatively adaptive genetic variation in the giant California sea cucumber (
Parastichopus californicus
) as revealed by environmental association analysis of restriction‐site associated DNA sequencing data. Mol Ecol 2018; 27:5035-5048. [DOI: 10.1111/mec.14942] [Citation(s) in RCA: 29] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2018] [Revised: 10/17/2018] [Accepted: 10/26/2018] [Indexed: 12/22/2022]
Affiliation(s)
- Amanda Xuereb
- Department of Ecology and Evolutionary Biology University of Toronto Toronto Ontario Canada
| | - Christopher M. Kimber
- AVIAN Behavioural Genomics and Physiology Group, Department of Physics, Chemistry, and Biology (IFM) Linköping University Linköping Sweden
| | - Janelle M. R. Curtis
- Pacific Biological Station, Ecosystem Sciences Division Fisheries and Oceans Canada Nanaimo British Columbia Canada
| | - Louis Bernatchez
- Institut de Biologie Intégrative et des Systèmes Université Laval Québec Québec Canada
| | - Marie‐Josée Fortin
- Department of Ecology and Evolutionary Biology University of Toronto Toronto Ontario Canada
| |
Collapse
|
32
|
Kozakiewicz CP, Burridge CP, Funk WC, VandeWoude S, Craft ME, Crooks KR, Ernest HB, Fountain‐Jones NM, Carver S. Pathogens in space: Advancing understanding of pathogen dynamics and disease ecology through landscape genetics. Evol Appl 2018; 11:1763-1778. [PMID: 30459828 PMCID: PMC6231466 DOI: 10.1111/eva.12678] [Citation(s) in RCA: 26] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2018] [Revised: 06/24/2018] [Accepted: 06/28/2018] [Indexed: 12/30/2022] Open
Abstract
Landscape genetics has provided many insights into how heterogeneous landscape features drive processes influencing spatial genetic variation in free-living organisms. This rapidly developing field has focused heavily on vertebrates, and expansion of this scope to the study of infectious diseases holds great potential for landscape geneticists and disease ecologists alike. The potential application of landscape genetics to infectious agents has garnered attention at formative stages in the development of landscape genetics, but systematic examination is lacking. We comprehensively review how landscape genetics is being used to better understand pathogen dynamics. We characterize the field and evaluate the types of questions addressed, approaches used and systems studied. We also review the now established landscape genetic methods and their realized and potential applications to disease ecology. Lastly, we identify emerging frontiers in the landscape genetic study of infectious agents, including recent phylogeographic approaches and frameworks for studying complex multihost and host-vector systems. Our review emphasizes the expanding utility of landscape genetic methods available for elucidating key pathogen dynamics (particularly transmission and spread) and also how landscape genetic studies of pathogens can provide insight into host population dynamics. Through this review, we convey how increasing awareness of the complementarity of landscape genetics and disease ecology among practitioners of each field promises to drive important cross-disciplinary advances.
Collapse
Affiliation(s)
| | | | - W. Chris Funk
- Department of BiologyGraduate Degree Program in EcologyColorado State UniversityFort CollinsColorado
| | - Sue VandeWoude
- Department of Microbiology, Immunology, and PathologyColorado State UniversityFort CollinsColorado
| | - Meggan E. Craft
- Department of Veterinary Population MedicineUniversity of MinnesotaSt. PaulMinnesota
| | - Kevin R. Crooks
- Department of Fish, Wildlife, and Conservation BiologyColorado State UniversityFort CollinsColorado
| | - Holly B. Ernest
- Wildlife Genomics and Disease Ecology LaboratoryDepartment of Veterinary SciencesUniversity of WyomingLaramieWyoming
| | | | - Scott Carver
- School of Natural SciencesUniversity of TasmaniaHobartTasmaniaAustralia
| |
Collapse
|
33
|
Population genomic footprints of environmental pollution pressure in natural populations of the Mediterranean mussel. Mar Genomics 2018; 45:11-15. [PMID: 30447893 DOI: 10.1016/j.margen.2018.10.009] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2018] [Revised: 10/09/2018] [Accepted: 10/30/2018] [Indexed: 11/20/2022]
Abstract
Bivalve molluscs of the genus Mytilus are considered a model organism in ecotoxicology and are known to be well adapted to marine ecosystems affected by multiple anthropogenic factors, including pollution. In order to assess whether pollution interferes with the reproductive success of Mytilus and affects the diversity within and between populations, we sequenced the transcriptomes of 72 individuals from 9 populations of Mytilus galloprovincialis collected along a ca. 130-km north-south transect on the Western coast of the Iberian Peninsula. We found that polluted areas are acting as a barrier to gene flow, potentially because of the detrimental effect of anthropogenic chemicals on larvae carried from more pristine environments. Furthermore, we observed an increase in genetic diversity in populations from polluted site, which could be indicative of higher mutagenicity driven by the environment. We propose that a microevolutionary perspective is essential for a full characterization of human activities on the dispersal of M. galloprovincialis and that it should be incorporated into management, and conservation plans and policies in the context of the effects of pollution on populations.
Collapse
|
34
|
Brown AP, Arias-Rodriguez L, Yee MC, Tobler M, Kelley JL. Concordant Changes in Gene Expression and Nucleotides Underlie Independent Adaptation to Hydrogen-Sulfide-Rich Environments. Genome Biol Evol 2018; 10:2867-2881. [PMID: 30215710 PMCID: PMC6225894 DOI: 10.1093/gbe/evy198] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 09/11/2018] [Indexed: 12/23/2022] Open
Abstract
The colonization of novel environments often involves changes in gene expression, protein coding sequence, or both. Studies of how populations adapt to novel conditions, however, often focus on only one of these two processes, potentially missing out on the relative importance of different parts of the evolutionary process. In this study, our objectives were 1) to better understand the qualitative concordance between conclusions drawn from analyses of differential expression and changes in genic sequence and 2) to quantitatively test whether differentially expressed genes were enriched for sites putatively under positive selection within gene regions. To achieve this, we compared populations of fish (Poecilia mexicana) that have independently adapted to hydrogen-sulfide-rich environments in southern Mexico to adjacent populations residing in nonsulfidic waters. Specifically, we used RNA-sequencing data to compare both gene expression and DNA sequence differences between populations. Analyzing these two different data types led to similar conclusions about which biochemical pathways (sulfide detoxification and cellular respiration) were involved in adaptation to sulfidic environments. Additionally, we found a greater overlap between genes putatively under selection and differentially expressed genes than expected by chance. We conclude that considering both differential expression and changes in DNA sequence led to a more comprehensive understanding of how these populations adapted to extreme environmental conditions. Our results imply that changes in both gene expression and DNA sequence-sometimes at the same loci-may be involved in adaptation.
Collapse
Affiliation(s)
- Anthony P Brown
- School of Biological Sciences, Washington State University, 100 Dairy Road, Pullman, WA 99164
| | - Lenin Arias-Rodriguez
- División Académica de Ciencias Biológicas, Universidad Juárez Autónoma de Tabasco (UJAT), C.P. 86150, Villahermosa, Tabasco, México
| | - Muh-Ching Yee
- Stanford Functional Genomics Facility, CCSR 0120, Stanford, CA 94305
| | - Michael Tobler
- Division of Biology, Kansas State University, 116 Ackert Hall, Manhattan, KS 66506
| | - Joanna L Kelley
- School of Biological Sciences, Washington State University, 100 Dairy Road, Pullman, WA 99164
| |
Collapse
|
35
|
Hoey JA, Pinsky ML. Genomic signatures of environmental selection despite near-panmixia in summer flounder. Evol Appl 2018; 11:1732-1747. [PMID: 30344639 PMCID: PMC6183468 DOI: 10.1111/eva.12676] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2017] [Revised: 06/13/2018] [Accepted: 06/16/2018] [Indexed: 01/01/2023] Open
Abstract
Rapid environmental change is altering the selective pressures experienced by marine species. While adaptation to local environmental conditions depends on a balance between dispersal and natural selection across the seascape, the spatial scale of adaptation and the relative importance of mechanisms maintaining adaptation in the ocean are not well understood. Here, using population assignment tests, Approximate Bayesian Computation (ABC), and genome scans with double-digest restriction-site associated DNA sequencing data, we evaluated population structure and locus-environment associations in a commercially important species, summer flounder (Paralichthys dentatus), along the U.S. east coast. Based on 1,137 single nucleotide polymorphisms across 232 individuals spanning nearly 1,900 km, we found no indication of population structure across Cape Hatteras, North Carolina (F ST = 0.0014) or of isolation by distance along the coast using individual relatedness. ABC estimated the probability of dispersal across the biogeographic break at Cape Hatteras to be high (95% credible interval: 7%-50% migration). However, we found 15 loci whose allele frequencies were associated with at least one of four environmental variables. Of those, 11 were correlated with bottom temperature. For summer flounder, our results suggest continued fisheries management as a single population and identify likely response mechanisms to climate change. Broadly speaking, our findings suggest that spatial balancing selection can manifest in adaptive divergence on regional scales in marine fish despite high dispersal, and that these conditions likely result in the widespread distribution of adaptive alleles and a high potential for future genetic adaptation in response to changing environmental conditions. In the context of a rapidly changing world, a landscape genomics perspective offers a useful approach for understanding the causes and consequences of genetic differentiation.
Collapse
Affiliation(s)
- Jennifer A. Hoey
- Department of Ecology, Evolution, & Natural ResourcesRutgers UniversityNew BrunswickNew JerseyUSA
| | - Malin L. Pinsky
- Department of Ecology, Evolution, & Natural ResourcesRutgers UniversityNew BrunswickNew JerseyUSA
| |
Collapse
|
36
|
Nielsen ES, Henriques R, Toonen RJ, Knapp ISS, Guo B, von der Heyden S. Complex signatures of genomic variation of two non-model marine species in a homogeneous environment. BMC Genomics 2018; 19:347. [PMID: 29743012 PMCID: PMC5944137 DOI: 10.1186/s12864-018-4721-y] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2017] [Accepted: 04/23/2018] [Indexed: 12/26/2022] Open
Abstract
BACKGROUND Genomic tools are increasingly being used on non-model organisms to provide insights into population structure and variability, including signals of selection. However, most studies are carried out in regions with distinct environmental gradients or across large geographical areas, in which local adaptation is expected to occur. Therefore, the focus of this study is to characterize genomic variation and selective signals over short geographic areas within a largely homogeneous region. To assess adaptive signals between microhabitats within the rocky shore, we compared genomic variation between the Cape urchin (Parechinus angulosus), which is a low to mid-shore species, and the Granular limpet (Scutellastra granularis), a high shore specialist. RESULTS Using pooled restriction site associated DNA (RAD) sequencing, we described patterns of genomic variation and identified outlier loci in both species. We found relatively low numbers of outlier SNPs within each species, and identified outlier genes associated with different selective pressures than those previously identified in studies conducted over larger environmental gradients. The number of population-specific outlier loci differed between species, likely owing to differential selective pressures within the intertidal environment. Interestingly, the outlier loci were highly differentiated within the two northernmost populations for both species, suggesting that unique evolutionary forces are acting on marine invertebrates within this region. CONCLUSIONS Our study provides a background for comparative genomic studies focused on non-model species, as well as a baseline for the adaptive potential of marine invertebrates along the South African west coast. We also discuss the caveats associated with Pool-seq and potential biases of sequencing coverage on downstream genomic metrics. The findings provide evidence of species-specific selective pressures within a homogeneous environment, and suggest that selective forces acting on small scales are just as crucial to acknowledge as those acting on larger scales. As a whole, our findings imply that future population genomic studies should expand from focusing on model organisms and/or studying heterogeneous regions to better understand the evolutionary processes shaping current and future biodiversity patterns, particularly when used in a comparative phylogeographic context.
Collapse
Affiliation(s)
- Erica S Nielsen
- Evolutionary Genomics Group, Department of Botany and Zoology, University of Stellenbosch, Private Bag X1, Matieland,, 7602, South Africa
| | - Romina Henriques
- Evolutionary Genomics Group, Department of Botany and Zoology, University of Stellenbosch, Private Bag X1, Matieland,, 7602, South Africa
| | - Robert J Toonen
- Hawai'i Institute of Marine Biology, School of Ocean and Earth Science and Technology, University of Hawai'i at Mānoa, Kāne'ohe, HI, 96744, USA
| | - Ingrid S S Knapp
- Hawai'i Institute of Marine Biology, School of Ocean and Earth Science and Technology, University of Hawai'i at Mānoa, Kāne'ohe, HI, 96744, USA
| | - Baocheng Guo
- The Key Laboratory of Zoological Systematics and Evolution, Institute of Zoology Chinese Academy of Sciences, Beijing, 100101, China
| | - Sophie von der Heyden
- Evolutionary Genomics Group, Department of Botany and Zoology, University of Stellenbosch, Private Bag X1, Matieland,, 7602, South Africa.
| |
Collapse
|
37
|
Sandoval‐Castillo J, Robinson NA, Hart AM, Strain LWS, Beheregaray LB. Seascape genomics reveals adaptive divergence in a connected and commercially important mollusc, the greenlip abalone (
Haliotis laevigata
), along a longitudinal environmental gradient. Mol Ecol 2018; 27:1603-1620. [DOI: 10.1111/mec.14526] [Citation(s) in RCA: 83] [Impact Index Per Article: 13.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/29/2017] [Revised: 12/05/2017] [Accepted: 12/15/2017] [Indexed: 12/12/2022]
Affiliation(s)
- Jonathan Sandoval‐Castillo
- Molecular Ecology Laboratory College of Science and Engineering Flinders University Adelaide SA Australia
| | - Nick A. Robinson
- Nofima Ås Norway
- Sustainable Aquaculture Laboratory School of BioSciences University of Melbourne Parkville Vic Australia
| | - Anthony M. Hart
- Western Australian Fisheries and Marine Research Laboratories Department of Fisheries Western Australia Hillarys WA Australia
| | - Lachlan W. S. Strain
- Western Australian Fisheries and Marine Research Laboratories Department of Fisheries Western Australia Hillarys WA Australia
| | - Luciano B. Beheregaray
- Molecular Ecology Laboratory College of Science and Engineering Flinders University Adelaide SA Australia
| |
Collapse
|
38
|
Van Wyngaarden M, Snelgrove PVR, DiBacco C, Hamilton LC, Rodríguez‐Ezpeleta N, Zhan L, Beiko RG, Bradbury IR. Oceanographic variation influences spatial genomic structure in the sea scallop, Placopecten magellanicus. Ecol Evol 2018; 8:2824-2841. [PMID: 29531698 PMCID: PMC5838053 DOI: 10.1002/ece3.3846] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/22/2017] [Revised: 12/05/2017] [Accepted: 12/06/2017] [Indexed: 01/03/2023] Open
Abstract
Environmental factors can influence diversity and population structure in marine species and accurate understanding of this influence can both improve fisheries management and help predict responses to environmental change. We used 7163 SNPs derived from restriction site-associated DNA sequencing genotyped in 245 individuals of the economically important sea scallop, Placopecten magellanicus, to evaluate the correlations between oceanographic variation and a previously identified latitudinal genomic cline. Sea scallops span a broad latitudinal area (>10 degrees), and we hypothesized that climatic variation significantly drives clinal trends in allele frequency. Using a large environmental dataset, including temperature, salinity, chlorophyll a, and nutrient concentrations, we identified a suite of SNPs (285-621, depending on analysis and environmental dataset) potentially under selection through correlations with environmental variation. Principal components analysis of different outlier SNPs and environmental datasets revealed similar northern and southern clusters, with significant associations between the first axes of each (R2adj = .66-.79). Multivariate redundancy analysis of outlier SNPs and the environmental principal components indicated that environmental factors explained more than 32% of the variance. Similarly, multiple linear regressions and random-forest analysis identified winter average and minimum ocean temperatures as significant parameters in the link between genetic and environmental variation. This work indicates that oceanographic variation is associated with the observed genomic cline in this species and that seasonal periods of extreme cold may restrict gene flow along a latitudinal gradient in this marine benthic bivalve. Incorporating this finding into management may improve accuracy of management strategies and future predictions.
Collapse
Affiliation(s)
| | - Paul V. R. Snelgrove
- Department of BiologyMemorial University of NewfoundlandSt. John'sNLCanada
- Department of Ocean SciencesMemorial University of NewfoundlandSt. John'sNLCanada
| | | | | | | | - Luyao Zhan
- Faculty of Computer ScienceDalhousie UniversityHalifaxNSCanada
| | - Robert G. Beiko
- Faculty of Computer ScienceDalhousie UniversityHalifaxNSCanada
| | - Ian R. Bradbury
- Faculty of Computer ScienceDalhousie UniversityHalifaxNSCanada
- Fisheries and Oceans CanadaNorthwest Atlantic Fisheries CentreSt. John'sNLCanada
| |
Collapse
|
39
|
Graham AM, Lavretsky P, Muñoz-Fuentes V, Green AJ, Wilson RE, McCracken KG. Migration-Selection Balance Drives Genetic Differentiation in Genes Associated with High-Altitude Function in the Speckled Teal (Anas flavirostris) in the Andes. Genome Biol Evol 2018; 10:14-32. [PMID: 29211852 PMCID: PMC5757641 DOI: 10.1093/gbe/evx253] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 11/30/2017] [Indexed: 12/30/2022] Open
Abstract
Local adaptation frequently occurs across populations as a result of migration-selection balance between divergent selective pressures and gene flow associated with life in heterogeneous landscapes. Studying the effects of selection and gene flow on the adaptation process can be achieved in systems that have recently colonized extreme environments. This study utilizes an endemic South American duck species, the speckled teal (Anas flavirostris), which has both high- and low-altitude populations. High-altitude speckled teal (A. f. oxyptera) are locally adapted to the Andean environment and mostly allopatric from low-altitude birds (A. f. flavirostris); however, there is occasional gene flow across altitudinal gradients. In this study, we used next-generation sequencing to explore genetic patterns associated with high-altitude adaptation in speckled teal populations, as well as the extent to which the balance between selection and migration have affected genetic architecture. We identified a set of loci with allele frequencies strongly correlated with altitude, including those involved in the insulin-like signaling pathway, bone morphogenesis, oxidative phosphorylation, responders to hypoxia-induced DNA damage, and feedback loops to the hypoxia-inducible factor pathway. These same outlier loci were found to have depressed gene flow estimates, as well as being highly concentrated on the Z-chromosome. Our results suggest a multifactorial response to life at high altitudes through an array of interconnected pathways that are likely under positive selection and whose genetic components seem to be providing an effective genomic barrier to interbreeding, potentially functioning as an avenue for population divergence and speciation.
Collapse
Affiliation(s)
| | | | - Violeta Muñoz-Fuentes
- European Molecular Biology Laboratory, European Bioinformatics Institute, Wellcome Trust Genome Campus, Hinxton, Cambridge, United Kingdom
- Estación Biológica de Doñana, EBD-CSIC, Sevilla, Spain
| | - Andy J Green
- Estación Biológica de Doñana, EBD-CSIC, Sevilla, Spain
| | - Robert E Wilson
- Institute of Arctic Biology and University of Alaska Museum, University of Alaska, Fairbanks
| | - Kevin G McCracken
- Department of Biology, University of Miami
- Institute of Arctic Biology and University of Alaska Museum, University of Alaska, Fairbanks
- Rosenstiel School of Marine and Atmospheric Sciences, University of Miami
- John P. Hussman Institute for Human Genomics, University of Miami Miller School of Medicine
| |
Collapse
|
40
|
Zhao JS, Wang AY, Zhao HB, Chen YH. Transcriptome sequencing and differential gene expression analysis of the schistosome-transmitting snail Oncomelania hupensis inhabiting hilly and marshland regions. Sci Rep 2017; 7:15809. [PMID: 29150650 PMCID: PMC5693929 DOI: 10.1038/s41598-017-16084-z] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2017] [Accepted: 11/07/2017] [Indexed: 12/18/2022] Open
Abstract
The freshwater snail Oncomelania hupensis is the unique intermediate host of the blood fluke Schistosoma japonicum, which is the major cause of schistosomiasis. The snail inhabits two contrasting environments: the hilly and marshland regions. The hilly snails are smaller in size and have the typical smooth shell, whereas the marshland snails are larger and possess the ribbed shell. To reveal the differences in gene expression between the hilly and marshland snails, a total of six snails, three per environment, were individually examined by RNA sequencing technology. All paired-end reads were assembled into contigs from which 34,760 unigenes were predicted. Based on single nucleotide polymorphisms, principal component analysis and neighbor-joining clustering revealed two distinct clusters of hilly and marshland snails. Analysis of expression changes between environments showed that upregulated genes relating to immunity and development were enriched in hilly snails, while those associated with reproduction were over-represented in marshland snails. Eight differentially expressed genes between the two types of snails were validated by qRT-PCR. Our study identified candidate genes that could be targets for future functional studies, and provided a link between expression profiling and ecological adaptation of the snail that may have implications for schistosomiasis control.
Collapse
Affiliation(s)
- Jin-Song Zhao
- School of Basic Medicine, Wannan Medical College, Wuhu, 241002, China
| | - An-Yun Wang
- School of Public Health, Wannan Medical College, Wuhu, 241002, China
| | - Hua-Bin Zhao
- College of Life Sciences, Wuhan University, Wuhan, 430072, China
| | - Yan-Hong Chen
- College of Life Sciences, Wuhan University, Wuhan, 430072, China.
| |
Collapse
|
41
|
Bay RA, Rose NH, Logan CA, Palumbi SR. Genomic models predict successful coral adaptation if future ocean warming rates are reduced. SCIENCE ADVANCES 2017; 3:e1701413. [PMID: 29109975 PMCID: PMC5665595 DOI: 10.1126/sciadv.1701413] [Citation(s) in RCA: 86] [Impact Index Per Article: 12.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/01/2017] [Accepted: 10/10/2017] [Indexed: 05/02/2023]
Abstract
Population genomic surveys suggest that climate-associated genetic variation occurs widely across species, but whether it is sufficient to allow population persistence via evolutionary adaptation has seldom been quantified. To ask whether rapid adaptation in reef-building corals can keep pace with future ocean warming, we measured genetic variation at predicted warm-adapted loci and simulated future evolution and persistence in a high-latitude population of corals from Rarotonga, Cook Islands. Alleles associated with thermal tolerance were present but at low frequencies in this cooler, southerly locality. Simulations based on predicted ocean warming in Rarotonga showed rapid evolution of heat tolerance resulting in population persistence under mild warming scenarios consistent with low CO2 emission plans, RCP2.6 and RCP4.5. Under more severe scenarios, RCP6.0 and RCP8.5, adaptation was not rapid enough to prevent extinction. Population adaptation was faster for models based on smaller numbers of additive loci that determine thermal tolerance and for higher population growth rates. Finally, accelerated migration via transplantation of thermally tolerant individuals (1 to 5%/year) sped adaptation. These results show that cool-water corals can adapt to warmer oceans but only under mild scenarios resulting from international emissions controls. Incorporation of genomic data into models of species response to climate change offers a promising method for estimating future adaptive processes.
Collapse
Affiliation(s)
- Rachael A. Bay
- Department of Biology, Hopkins Marine Station of Stanford University, Pacific Grove, CA 93950, USA
- Corresponding author.
| | - Noah H. Rose
- Department of Biology, Hopkins Marine Station of Stanford University, Pacific Grove, CA 93950, USA
| | - Cheryl A. Logan
- School of Natural Sciences, California State University, Monterey Bay, Seaside, CA 93955, USA
| | - Stephen R. Palumbi
- Department of Biology, Hopkins Marine Station of Stanford University, Pacific Grove, CA 93950, USA
| |
Collapse
|
42
|
Harris SE, Munshi-South J. Signatures of positive selection and local adaptation to urbanization in white-footed mice (Peromyscus leucopus). Mol Ecol 2017; 26:6336-6350. [PMID: 28980357 DOI: 10.1111/mec.14369] [Citation(s) in RCA: 49] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2016] [Accepted: 09/25/2017] [Indexed: 02/06/2023]
Abstract
Urbanization significantly alters natural ecosystems and has accelerated globally. Urban wildlife populations are often highly fragmented by human infrastructure, and isolated populations may adapt in response to local urban pressures. However, relatively few studies have identified genomic signatures of adaptation in urban animals. We used a landscape genomic approach to examine signatures of selection in urban populations of white-footed mice (Peromyscus leucopus) in New York City. We analysed 154,770 SNPs identified from transcriptome data from 48 P. leucopus individuals from three urban and three rural populations and used outlier tests to identify evidence of urban adaptation. We accounted for demography by simulating a neutral SNP data set under an inferred demographic history as a null model for outlier analysis. We also tested whether candidate genes were associated with environmental variables related to urbanization. In total, we detected 381 outlier loci and after stringent filtering, identified and annotated 19 candidate loci. Many of the candidate genes were involved in metabolic processes and have well-established roles in metabolizing lipids and carbohydrates. Our results indicate that white-footed mice in New York City are adapting at the biomolecular level to local selective pressures in urban habitats. Annotation of outlier loci suggests selection is acting on metabolic pathways in urban populations, likely related to novel diets in cities that differ from diets in less disturbed areas.
Collapse
Affiliation(s)
- Stephen E Harris
- The Graduate Center, City University of New York (CUNY), New York, NY, USA
| | - Jason Munshi-South
- Louis Calder Center-Biological Field Station, Fordham University, Armonk, NY, USA
| |
Collapse
|
43
|
Song K, Li L, Zhang G. Bias and Correction in RNA-seq Data for Marine Species. MARINE BIOTECHNOLOGY (NEW YORK, N.Y.) 2017; 19:541-550. [PMID: 28884399 DOI: 10.1007/s10126-017-9773-5] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 06/12/2017] [Accepted: 08/18/2017] [Indexed: 06/07/2023]
Abstract
RNA-seq is a recently developed approach widely used for transcriptome profiling in biological analyses that use next-generation sequencing technologies. Accurate estimation of gene expression levels is critical for answering biological questions. Here, we show that the commonly used measure of gene expression levels, fragments per kilobase of transcript per million mapped reads (FPKM), is biased in transcript length, GC content, and dinucleotide frequencies in the RNA-seq analysis of marine species. We used a generalized linear model to correct the observed biases of FPKM. We used RNA-seq data sets from eight species obtained by different sequencing methods to evaluate the correction methods. Our work contributes to the understanding of potential technical artifacts in RNA-seq experiments for marine species, and presents a means by which more accurate gene expression measures can be obtained.
Collapse
Affiliation(s)
- Kai Song
- Key Laboratory of Experimental Marine Biology, Institute of Oceanology, Chinese Academy of Sciences, 7th Nanhai Rd, Qingdao, Shandong, 266071, China
- National & Local Joint Engineering Laboratory of Ecological Mariculture, Qingdao, Shandong, 266071, China
- Laboratory for Marine Fisheries and Aquaculture, Qingdao National Laboratory for Marine Science and Technology, Qingdao, Shandong, 266071, China
| | - Li Li
- Key Laboratory of Experimental Marine Biology, Institute of Oceanology, Chinese Academy of Sciences, 7th Nanhai Rd, Qingdao, Shandong, 266071, China.
- National & Local Joint Engineering Laboratory of Ecological Mariculture, Qingdao, Shandong, 266071, China.
- Laboratory for Marine Fisheries and Aquaculture, Qingdao National Laboratory for Marine Science and Technology, Qingdao, Shandong, 266071, China.
| | - Guofan Zhang
- Key Laboratory of Experimental Marine Biology, Institute of Oceanology, Chinese Academy of Sciences, 7th Nanhai Rd, Qingdao, Shandong, 266071, China.
- Laboratory for Marine Fisheries and Aquaculture, Qingdao National Laboratory for Marine Science and Technology, Qingdao, Shandong, 266071, China.
- Laboratory for Marine Biology and Biotechnology, Qingdao National Laboratory for Marine Science and Technology, Qingdao, Shandong, 266071, China.
| |
Collapse
|
44
|
Taylor ML, Roterman CN. Invertebrate population genetics across Earth's largest habitat: The deep-sea floor. Mol Ecol 2017; 26:4872-4896. [PMID: 28833857 DOI: 10.1111/mec.14237] [Citation(s) in RCA: 64] [Impact Index Per Article: 9.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2017] [Revised: 06/16/2017] [Accepted: 06/19/2017] [Indexed: 01/04/2023]
Abstract
Despite the deep sea being the largest habitat on Earth, there are just 77 population genetic studies of invertebrates (115 species) inhabiting non-chemosynthetic ecosystems on the deep-sea floor (below 200 m depth). We review and synthesize the results of these papers. Studies reveal levels of genetic diversity comparable to shallow-water species. Generally, populations at similar depths were well connected over 100s-1,000s km, but studies that sampled across depth ranges reveal population structure at much smaller scales (100s-1,000s m) consistent with isolation by adaptation across environmental gradients, or the existence of physical barriers to connectivity with depth. Few studies were ocean-wide (under 4%), and 48% were Atlantic-focused. There is strong emphasis on megafauna and commercial species with research into meiofauna, "ecosystem engineers" and other ecologically important species lacking. Only nine papers account for ~50% of the planet's surface (depths below 3,500 m). Just two species were studied below 5,000 m, a quarter of Earth's seafloor. Most studies used single-locus mitochondrial genes revealing a common pattern of non-neutrality, consistent with demographic instability or selective sweeps; similar to deep-sea hydrothermal vent fauna. The absence of a clear difference between vent and non-vent could signify that demographic instability is common in the deep sea, or that selective sweeps render single-locus mitochondrial studies demographically uninformative. The number of population genetics studies to date is miniscule in relation to the size of the deep sea. The paucity of studies constrains meta-analyses where broad inferences about deep-sea ecology could be made.
Collapse
Affiliation(s)
- M L Taylor
- Department of Zoology, University of Oxford, Oxford, UK
| | - C N Roterman
- Department of Zoology, University of Oxford, Oxford, UK
| |
Collapse
|
45
|
Benestan L, Moore JS, Sutherland BJG, Le Luyer J, Maaroufi H, Rougeux C, Normandeau E, Rycroft N, Atema J, Harris LN, Tallman RF, Greenwood SJ, Clark FK, Bernatchez L. Sex matters in massive parallel sequencing: Evidence for biases in genetic parameter estimation and investigation of sex determination systems. Mol Ecol 2017; 26:6767-6783. [DOI: 10.1111/mec.14217] [Citation(s) in RCA: 35] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/20/2016] [Revised: 03/23/2017] [Accepted: 03/29/2017] [Indexed: 12/26/2022]
Affiliation(s)
- Laura Benestan
- Département de Biologie; Université Laval; Québec QC Canada
- Institut de Biologie Intégrative et des Systèmes (IBIS); Université Laval; Québec QC Canada
| | - Jean-Sébastien Moore
- Département de Biologie; Université Laval; Québec QC Canada
- Institut de Biologie Intégrative et des Systèmes (IBIS); Université Laval; Québec QC Canada
| | - Ben J. G. Sutherland
- Département de Biologie; Université Laval; Québec QC Canada
- Institut de Biologie Intégrative et des Systèmes (IBIS); Université Laval; Québec QC Canada
| | - Jérémy Le Luyer
- Département de Biologie; Université Laval; Québec QC Canada
- Institut de Biologie Intégrative et des Systèmes (IBIS); Université Laval; Québec QC Canada
| | - Halim Maaroufi
- Institut de Biologie Intégrative et des Systèmes (IBIS); Université Laval; Québec QC Canada
| | - Clément Rougeux
- Département de Biologie; Université Laval; Québec QC Canada
- Institut de Biologie Intégrative et des Systèmes (IBIS); Université Laval; Québec QC Canada
| | - Eric Normandeau
- Département de Biologie; Université Laval; Québec QC Canada
- Institut de Biologie Intégrative et des Systèmes (IBIS); Université Laval; Québec QC Canada
| | | | - Jelle Atema
- Department of Biology; Boston University; Boston MA USA
| | - Les N. Harris
- Fisheries and Oceans Canada; Freshwater Institute; Winnipeg MB Canada
| | - Ross F. Tallman
- Fisheries and Oceans Canada; Freshwater Institute; Winnipeg MB Canada
| | - Spencer J. Greenwood
- Department of Biomedical Sciences & AVC Lobster Science Centre; Atlantic Veterinary College; University of Prince Edward Island; Charlottetown PE Canada
| | - Fraser K. Clark
- Department of Biomedical Sciences & AVC Lobster Science Centre; Atlantic Veterinary College; University of Prince Edward Island; Charlottetown PE Canada
| | - Louis Bernatchez
- Département de Biologie; Université Laval; Québec QC Canada
- Institut de Biologie Intégrative et des Systèmes (IBIS); Université Laval; Québec QC Canada
| |
Collapse
|
46
|
O'Leary JK, Barry JP, Gabrielson PW, Rogers-Bennett L, Potts DC, Palumbi SR, Micheli F. Calcifying algae maintain settlement cues to larval abalone following algal exposure to extreme ocean acidification. Sci Rep 2017; 7:5774. [PMID: 28720836 PMCID: PMC5515930 DOI: 10.1038/s41598-017-05502-x] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2017] [Accepted: 05/30/2017] [Indexed: 11/22/2022] Open
Abstract
Ocean acidification (OA) increasingly threatens marine systems, and is especially harmful to calcifying organisms. One important question is whether OA will alter species interactions. Crustose coralline algae (CCA) provide space and chemical cues for larval settlement. CCA have shown strongly negative responses to OA in previous studies, including disruption of settlement cues to corals. In California, CCA provide cues for seven species of harvested, threatened, and endangered abalone. We exposed four common CCA genera and a crustose calcifying red algae, Peyssonnelia (collectively CCRA) from California to three pCO2 levels ranging from 419–2,013 µatm for four months. We then evaluated abalone (Haliotis rufescens) settlement under ambient conditions among the CCRA and non-algal controls that had been previously exposed to the pCO2 treatments. Abalone settlement and metamorphosis increased from 11% in the absence of CCRA to 45–69% when CCRA were present, with minor variation among CCRA genera. Though all CCRA genera reduced growth during exposure to increased pCO2, abalone settlement was unaffected by prior CCRA exposure to increased pCO2. Thus, we find no impacts of OA exposure history on CCRA provision of settlement cues. Additionally, there appears to be functional redundancy in genera of CCRA providing cues to abalone, which may further buffer OA effects.
Collapse
Affiliation(s)
- Jennifer K O'Leary
- Hopkins Marine Station, Stanford University, Monterey, Pacific Grove, United States of America. .,California Sea Grant, Department of Biology, California Polytechnic State University, San Luis Obispo, United States of America.
| | - James P Barry
- Monterey Bay Aquarium Research Institute, Moss Landing, California, United States of America
| | - Paul W Gabrielson
- Biology Department, Herbarium, University of North Carolina, Chapel Hill, North Carolina, United States of America
| | - Laura Rogers-Bennett
- Bodega Marine Laboratory, University of California, Davis, California, United States of America.,California Department of Fish and Wildlife, Marine Region, Bodega Bay, California, United States of America
| | - Donald C Potts
- Department of Ecology and Evolutionary Biology, University of California, Santa Cruz, California, United States of America
| | - Stephen R Palumbi
- Hopkins Marine Station, Stanford University, Monterey, Pacific Grove, United States of America
| | - Fiorenza Micheli
- Hopkins Marine Station, Stanford University, Monterey, Pacific Grove, United States of America
| |
Collapse
|
47
|
Bay RA, Rose N, Barrett R, Bernatchez L, Ghalambor CK, Lasky JR, Brem RB, Palumbi SR, Ralph P. Predicting Responses to Contemporary Environmental Change Using Evolutionary Response Architectures. Am Nat 2017; 189:463-473. [DOI: 10.1086/691233] [Citation(s) in RCA: 95] [Impact Index Per Article: 13.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/03/2022]
|
48
|
Shiel BP, Hall NE, Cooke IR, Robinson NA, Strugnell JM. Epipodial Tentacle Gene Expression and Predetermined Resilience to Summer Mortality in the Commercially Important Greenlip Abalone, Haliotis laevigata. MARINE BIOTECHNOLOGY (NEW YORK, N.Y.) 2017; 19:191-205. [PMID: 28349286 PMCID: PMC5405107 DOI: 10.1007/s10126-017-9742-z] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/22/2016] [Accepted: 02/06/2017] [Indexed: 05/05/2023]
Abstract
"Summer mortality" is a phenomenon that occurs during warm water temperature spikes that results in the mass mortality of many ecologically and economically important mollusks such as abalone. This study aimed to determine whether the baseline gene expression of abalone before a laboratory-induced summer mortality event was associated with resilience to summer mortality. Tentacle transcriptomes of 35 greenlip abalone (Haliotis laevigata) were sequenced prior to the animals being exposed to an increase in water temperature-simulating conditions which have previously resulted in summer mortality. Abalone derived from three source locations with different environmental conditions were categorized as susceptible or resistant to summer mortality depending on whether they died or survived after the water temperature was increased. We detected two genes showing significantly higher expression in resilient abalone relative to susceptible abalone prior to the laboratory-induced summer mortality event. One of these genes was annotated through the NCBI non-redundant protein database using BLASTX to an anemone (Exaiptasia pallida) Transposon Ty3-G Gag Pol polyprotein. Distinct gene expression signatures were also found between resilient and susceptible abalone depending on the population origin, which may suggest divergence in local adaptation mechanisms for resilience. Many of these genes have been suggested to be involved in antioxidant and immune-related functions. The identification of these genes and their functional roles have enhanced our understanding of processes that may contribute to summer mortality in abalone. Our study supports the hypothesis that prestress gene expression signatures are indicative of the likelihood of summer mortality.
Collapse
Affiliation(s)
- Brett P Shiel
- Department of Ecology, Environment and Evolution, School of Life Sciences, La Trobe University, Kingsbury Drive, Melbourne, VIC, 3086, Australia.
| | - Nathan E Hall
- Department of Ecology, Environment and Evolution, School of Life Sciences, La Trobe University, Kingsbury Drive, Melbourne, VIC, 3086, Australia
- Life Sciences Computation Centre, VLSCI, Parkville, VIC, Australia
- Department of Molecular and Cell Biology, James Cook University, Townsville, Australia
| | - Ira R Cooke
- Life Sciences Computation Centre, VLSCI, Parkville, VIC, Australia
- Department of Molecular and Cell Biology, James Cook University, Townsville, Australia
- Department of Biochemistry, La Trobe Institute for Molecular Science, La Trobe University, Kingsbury Drive, Melbourne, VIC, 3086, Australia
| | - Nicholas A Robinson
- Nofima, P.O. Box 210, 1431, Ås, Norway
- Sustainable Aquaculture Laboratory-Temperate and Tropical (SALTT), School of BioSciences, The University of Melbourne, Parkville, VIC, 3010, Australia
| | - Jan M Strugnell
- Department of Ecology, Environment and Evolution, School of Life Sciences, La Trobe University, Kingsbury Drive, Melbourne, VIC, 3086, Australia
- Centre for Sustainable Tropical Fisheries and Aquaculture, College of Science and Engineering, James Cook University, Townsville, QLD, 4811, Australia
| |
Collapse
|
49
|
Aguilera F, McDougall C, Degnan BM. Co-Option and De Novo Gene Evolution Underlie Molluscan Shell Diversity. Mol Biol Evol 2017; 34:779-792. [PMID: 28053006 PMCID: PMC5400390 DOI: 10.1093/molbev/msw294] [Citation(s) in RCA: 42] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
Molluscs fabricate shells of incredible diversity and complexity by localized secretions from the dorsal epithelium of the mantle. Although distantly related molluscs express remarkably different secreted gene products, it remains unclear if the evolution of shell structure and pattern is underpinned by the differential co-option of conserved genes or the integration of lineage-specific genes into the mantle regulatory program. To address this, we compare the mantle transcriptomes of 11 bivalves and gastropods of varying relatedness. We find that each species, including four Pinctada (pearl oyster) species that diverged within the last 20 Ma, expresses a unique mantle secretome. Lineage- or species-specific genes comprise a large proportion of each species' mantle secretome. A majority of these secreted proteins have unique domain architectures that include repetitive, low complexity domains (RLCDs), which evolve rapidly, and have a proclivity to expand, contract and rearrange in the genome. There are also a large number of secretome genes expressed in the mantle that arose before the origin of gastropods and bivalves. Each species expresses a unique set of these more ancient genes consistent with their independent co-option into these mantle gene regulatory networks. From this analysis, we infer lineage-specific secretomes underlie shell diversity, and include both rapidly evolving RLCD-containing proteins, and the continual recruitment and loss of both ancient and recently evolved genes into the periphery of the regulatory network controlling gene expression in the mantle epithelium.
Collapse
Affiliation(s)
- Felipe Aguilera
- Centre for Marine Sciences, School of Biological Sciences, The University of Queensland, Brisbane, Australia
| | - Carmel McDougall
- Centre for Marine Sciences, School of Biological Sciences, The University of Queensland, Brisbane, Australia
| | - Bernard M. Degnan
- Centre for Marine Sciences, School of Biological Sciences, The University of Queensland, Brisbane, Australia
| |
Collapse
|
50
|
Barney BT, Munkholm C, Walt DR, Palumbi SR. Highly localized divergence within supergenes in Atlantic cod (Gadus morhua) within the Gulf of Maine. BMC Genomics 2017; 18:271. [PMID: 28359300 PMCID: PMC5374575 DOI: 10.1186/s12864-017-3660-3] [Citation(s) in RCA: 37] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2016] [Accepted: 03/25/2017] [Indexed: 01/01/2023] Open
Abstract
BACKGROUND Atlantic cod (Gadus morhua), is known to vary genetically across the North Atlantic, Greenland, and Newfoundland. This genetic variation occurs both spatially and temporally through decades of heavy fishing, and is concentrated in three linkage disequilibrium blocks, previously defined by pedigreed linkage mapping analysis. Variation within these genomic regions is correlated with both seawater temperature and behavioral ecotype. The full extent and nature of these linkage groups is important information for interpreting cod genetic structure as a tool for future fisheries management. RESULTS We conducted whole genome sequencing for 31 individual cod from three sub-populations in the Gulf of Maine. Across the genome, we found 3,390,654 intermediate to high frequency Single Nucleotide Polymorphisms (SNPs). We show that pairwise linkage analysis among these SNPs is a powerful tool to detect linkage disequilibrium clusters by recovering the three previously detected linkage groups and identifying the 1031 genes contained therein. Across these genes, we found significant population differentiation among spawning groups in the Gulf of Maine and between Georges Bank and Gulf of Maine. Coordinated divergence among these genes and their differentiation at both short and long spatial scales suggests that they are acting as linked supergenes in local adaptation of cod populations. CONCLUSIONS Differentiation between SNPs in linkage disequilibrium blocks is the major signal of genetic differentiation between all groups tested within the Gulf of Maine. Our data provide a map of genes contained in these blocks, allowing an enhanced search for neutral genetic structure for demographic inference and fisheries modeling. Patterns of selection and the history of populations may be possible to identify in cod using this description of linkage disequilibrium blocks and future data sets to robustly separate neutral and selected genetic markers.
Collapse
Affiliation(s)
- Bryan T. Barney
- Department of Biology, Hopkins Marine Station, Stanford University, 120 Ocean View Boulevard, Pacific Grove, CA 92950 USA
| | - Christiane Munkholm
- Gloucester Marine Genomics Institute, 55 Blackburn Center, Gloucester, MA 01930 USA
| | - David R. Walt
- Gloucester Marine Genomics Institute, 55 Blackburn Center, Gloucester, MA 01930 USA
| | - Stephen R. Palumbi
- Department of Biology, Hopkins Marine Station, Stanford University, 120 Ocean View Boulevard, Pacific Grove, CA 92950 USA
| |
Collapse
|